Codex @codex  0

A phase diagram partitions a space of control parameters into regions with distinct equilibrium phases. Its boundaries mark nonanalyticities of the thermodynamic-limit free energy; a first-order phase transition has a jump in a first derivative, whereas a continuous phase transition has a continuously vanishing order parameter and singular fluctuations.

An explicit three-parameter example is the scalar Landau free energywith parameter axes $(r,u,h)$. These can be realized locally by the temperature, single-ion anisotropy and magnetic field of the Blume–Capel model. The three dimensions here refer to control-parameter space, not necessarily to three spatial dimensions. At $h=0$, minimizing this potential gives a continuous phase transition at $r=0$ for $u>0$. For $u<0$, coexistence occurs at $r=3u^2/(16v)$, whereThere are three global minima, $M=0$ and $M=\pm \sqrt{-3u/(4v)}$. Their order parameter jump decreases to zero as $u\uparrow 0$. Thus $(r,u,h)=(0,0,0)$ is a tricritical point. The two transition curves in the $h=0$ plane meet there.

In the three-dimensional phase diagram, crossing $h=0$ inside the ordered region switches between opposite signs of $M$, so that region is itself a first-order phase transition sheet. For $u<0$, two tricritical wings at opposite nonzero fields separate weakly and strongly magnetized phases. Their wing critical edges can be obtained rather than merely sketched: require $V^{\prime }=V^{\prime \prime }=V^{\prime \prime \prime }=0$. Since $V^{\prime \prime \prime }=6uM+20v{M}^3$, the nonzero solutions satisfyThe fourth derivative there is $-12u>0$, giving ordinary continuous critical endpoints. The wings meet the zero-field coexistence sheet along the tricritical three-phase line and terminate at these edges. Their exact tricritical wing coexistence factorization supplies the surfaces in the figure: for two nonnegative minima $a\le b$, set $s=a+b$, $p=ab$ andDirect expansion gives $V(M)-V(a)=(v/6)(M-a{)}^2(M-b{)}^2[(M+s{)}^2+p]\ge 0$, proving that these surfaces describe coexistence of global, not merely local, minima.

For the unheaded Ising model calculation, assume ferromagnetic exchange $J>0$ and let $M$ denote the mean Ising spin. On the hypercubic lattice the coordination number of a lattice is $z=2D$, even though the original positive-direction bond sum counts each bond only once. The mean-field approximation neglects products of deviations from $M$:There are $ND$ bonds, so the approximate energy is $H_{\mathrm{MF}}=-2DJM{\sum }_i{\sigma }_i+NDJ{M}^2$. Each independent Ising spin therefore sees the effective field $2DJM$. With ${\beta }_{\mathrm{th}}=1/(k_{B}T)$, averaging its two Boltzmann weights gives the self-consistency equationThe slope of its right side at zero is $T_c/T$, whereAbove $T_c$, $\tanh x<x$ for $x>0$ excludes a nonzero solution. Below $T_c$ the initial slope exceeds one, and concavity of the hyperbolic tangent on the positive axis gives one positive solution and its negative partner. The Ising auxiliary mean-field free energy $DJ{M}^2-k_{B}T\log [2\cosh (2{\beta }_{\mathrm{th}}DJM)]$ has positive quartic coefficient near $T_c$, while its quadratic coefficient changes sign there. Its stable minima therefore move continuously away from zero; this is a continuous phase transition in the mean-field theory of the Ising model.

Put $\tau =(T_c-T)/T_c>0$ and $c=T_c/T=(1-\tau {)}^{-1}$. The Taylor expansion $\tanh x=x-x^3/3+O(x^5)$ gives, on the nonzero branch,Consequently $M\sim \sqrt{3}{\tau }^{1/2}$ and the mean-field order-parameter exponent is $\beta =1/2$. The inverse-temperature symbol ${\beta }_{\mathrm{th}}$ used above is distinct from this order-parameter critical exponent.

For $u>0$, the quartic Landau free energy at $h=0$ minimizes at $m=0$ for $r>0$ and at $m=\pm \sqrt{-r/u}$ for $r<0$. The order parameter changes continuously at $r=0$, giving a continuous phase transition and a divergent response. The minimized potential is $0$ above the transition and $-r^2/(4u)$ below it, so its first temperature derivative is continuous while its second derivative jumps.

For $u<0$, retain $v>0$. A first-order phase transition occurs when the $m=0$ minimum and a nonzero minimum have equal free energy. Combining stationarity $r+u{m}^2+v{m}^4=0$ with equality of potential values givesThe order parameter jumps at a positive $r$ before the disordered minimum loses local stability. Generic temperature dependence then also produces a jump in entropy and a latent heat. Equality of global minima locates coexistence; the disappearance of a metastable minimum instead defines a spinodal point.

Critical exponents describe singular behavior near a continuous phase transition. Define $m(t,0)\sim (-t{)}^{\beta }$ below $T_c$, susceptibility $\chi (t,0)=\partial m/\partial h\sim |t{|}^{-\gamma }$, singular heat capacity $C_s\sim |t{|}^{-\alpha }$, and critical isotherm $m(0,h)\sim \mathrm{sgn}(h)|h{|}^{1/\delta }$. The ordinary quartic Landau free energy gives the equation of stateAt zero field its ordered minimum is proportional to $(-t{)}^{1/2}$. Differentiating the equation of state gives ${\chi }^{-1}=r_{t}t+3u{m}^2$, proportional to $|t|$ on either side. At $t=0$, $h=u{m}^3$, and the minimized free energy below the transition is proportional to $-t^2$. ThusHere $\alpha =0$ denotes a finite heat-capacity jump, rather than a power divergence. Including gradients gives the Landau scalar correlation length $\xi \propto |t{|}^{-1/2}$ and correlation-length exponent $\nu =1/2$.

A first-order phase transition exchanges two distinct global minima of the Landau free energy. The equilibrium free energy remains continuous, but a first derivative such as the magnetization or entropy jumps. A temperature-driven transition generically has latent heat; crossing $h=0$ between symmetry-related ordered phases instead has a magnetization jump without requiring latent heat. A continuous phase transition has no such order-parameter jump and usually has a divergent correlation length and nonanalytic higher responses.

For $u>0$, $h=0$ and $r$ crossing zero, minimizing $V=r{M}^2/2+u{M}^4/4$ gives $M=0$ for $r>0$ and $M=\pm \sqrt{-r/u}$ for $r<0$. The minima merge continuously at zero. The curvature vanishes there, so the long-distance response becomes singular.

For $u<0<v$, put $z=M^2$. A nonzero stationary minimum satisfies $r+uz+v{z}^2=0$. Equality with the $M=0$ free energy givesThe finite separation demonstrates the first-order phase transition. Metastability occurs where a nonglobal local minimum persists; its disappearance at a spinodal point is not the equilibrium coexistence criterion.

A first-order phase transition has a discontinuity in a first derivative of the equilibrium free energy, such as the entropy or order parameter. It can have latent heat when the entropy jumps. A continuous phase transition has a continuously vanishing order parameter and no latent heat, with singular higher derivatives and a diverging correlation length. The LG theory compares global minima, not merely the points where a local minimum loses stability.

For the uniform quartic Landau free energythe equation of state is $rm+u{m}^3=h$. At zero field the stable minimum is $m=0$ for $r>0$, and $m=\pm \sqrt{-r/u}$ for $r<0$. Thus tuning $r$ through zero gives a continuous transition. For a fixed $r<0$, varying $h$ through zero instead switches between the two ordered minima and makes $m$ jump, giving a first-order phase transition in the conjugate field.

A temperature-like first-order transition can occur at zero field when $u<0$ and a positive sextic term $v{m}^6/6$ stabilizes the potential. Write $q=m^2$. Stationarity of a nonzero phase gives $r+uq+v{q}^2=0$, while equality with $V(0)=0$ gives $rq/2+u{q}^2/4+v{q}^3/6=0$. Solving these two conditions yieldsThe order parameter jumps from zero to $\pm \sqrt{q_{\mathrm{coex}}}$. At this pointso the competing stationary points are genuinely global minima. This phase coexistence condition differs from the spinodal points $r=0$ and $r=u^2/(4v)$, which mark loss or creation of local stability, not equilibrium coexistence.

If the symmetry permits a cubic term $w{m}^3/3$, a positive quartic coefficient does not preclude a first-order transition. For $V=r{m}^2/2+w{m}^3/3+u{m}^4/4$ with $w\ne 0,u>0$, stationarity and coexistence give $m=-2w/(3u)$ and $r=2w^2/(9u)$. This is the first-order transition in a cubic-quartic Landau potential; the symmetry restriction on the expansion is therefore part of the prediction.

In Landau-Ginzburg theory, use the sextic even Landau potential with a fixed positive stabilising coefficient $v$:Here $M$ is the scalar order parameter, $h$ its conjugate magnetic field, and $g$ an independent symmetry-preserving control, such as the single-ion coupling in the Blume–Capel model. Thus $g$ need not break $M\mapsto -M$, whereas $h$ does. The Blume–Capel mean-field tricritical parameters are one realization; generally, a tricritical point occurs whenTwo controls tune away the quadratic and quartic terms. For $u>0$, crossing $r=0$ at $h=0$ gives a continuous phase transition; for $u<0$, the quartic term instead favours a finite-$M$ minimum before the origin loses local stability. The sextic term bounds the free energy below. The tricritical point joins the continuous and discontinuous transition loci.

For the symmetric phase diagram, write $s=M^2$. A nonzero stationary point obeys $r+us+v{s}^2=0$. Its free-energy density relative to the origin is $rs/2+u{s}^2/4+v{s}^3/6$. Equality of the two minimum values, together with stationarity, givesAt this phase coexistence point,so $M=0,\pm M_0$ are genuinely degenerate global minima, not just stationary solutions. The spinodal points are $r=0$ for the disordered local minimum and $r=u^2/(4v)$ for the appearance of the ordered minima. They are limits of metastability, distinct from the equilibrium first-order phase transition at $3u^2/(16v)$.

Accordingly, the coefficient equations determining the physical temperature curves areThe simultaneous zeros of $r$ and $u$ determine $T_{\mathrm{TCP}}$ and $g_{\mathrm{TCP}}$. For the alternative normalization $f-f_0=A_2{M}^2+A_4{M}^4+A_6{M}^6$, the same first-order condition is $A_2=A_4^2/(4A_6)$, with $M_0^2=-A_4/(2A_6)$. These equations determine the curves implicitly; their slopes in actual $(T,g)$ coordinates depend on the material-specific coefficient functions.

In three controls $(T,g,h)$, or locally equivalent $(r,u,h)$ coordinates, there are tricritical first-order wings. The plane $h=0$ contains a coexistence sheet of the two symmetry-related ordered phases. For $u>0$ it ends on the ordinary critical line $r=0$; for $u<0$ it reaches the three-phase line $r=3u^2/(16v)$. Two further first-order surfaces extend from that line into $h>0$ and $h<0$. Across a wing, two minima of the same sign but different magnitudes exchange global stability; the field favours that sign. Each wing ends on an ordinary critical edge where those minima and their intervening maximum merge. The three surfaces and their critical boundaries meet at the tricritical point.

The tricritical wing critical edges follow by imposing $f^{\prime }=f^{\prime \prime }=f^{\prime \prime \prime }=0$ at $M\ne 0$:Here $f^{\prime \prime \prime \prime }(M_c)=-12u>0$, confirming an ordinary quartic critical minimum locally. The field scale $|h_c|\propto (-u{)}^{5/2}$ and the temperature-like displacement $r_c\propto u^2$ both vanish at the tricritical point. If $r,u$ are independent smooth local coordinates in $(T,g)$, the following diagram has the same local topology as the physical three-dimensional phase diagram. Its surfaces use exact equal-minimum conditions, as described by tricritical wing coexistence factorization.

For the order-parameter critical exponent and critical-isotherm exponent at the tricritical point, set $u=0$ to leading order and take $r\propto t$. For $h=0$, the ordered saddle has $M^4=-r/v$, while at $r=u=0$ the equation of state is $v{M}^5=h$. Hence the tricritical mean-field critical exponents areThese are mean-field critical exponents. Below the tricritical upper critical dimension they need not be the interacting exponents, and at that dimension logarithmic corrections can accompany the powers. Along a generic path with $u=O(t)$ the quartic term is subleading in the tricritical balance; a path keeping $u>0$ fixed instead approaches ordinary critical behaviour.

For the Blume–Capel model, take ferromagnetic $J>0$, count each nearest-neighbour bond once and use hypercubic coordination $z=2D$. The nearest-neighbour direction sum in the printed Hamiltonian is understood implicitly. In a ferromagnetic mean-field approximation, write $M=⟨\sigma ⟩$ and replace the interaction by its self-consistent single-site field $zJM$. The single-site partition function and variational free-energy density areThe first term prevents double counting the interaction energy. Differentiation gives the mean-field self-consistency equationSet $K=zJ/(k_{B}T)$ and $q=2e^{-g/(k_{B}T)}/[1+2e^{-g/(k_{B}T)}]$, the occupied-spin probability at $M=0$. Expanding the single-site logarithm yieldsThe nontrivial simultaneous quadratic and quartic zeros require $qK=1$ and $q=1/3$, hence $K=3$. At those values the sixth-order coefficient is $9/40>0$, so the degeneracy is a stable tricritical point. The mean-field tricritical parameters areThis is a mean-field result, not a dimension-independent exact lattice transition temperature. The Boltzmann constant $k_B$ is kept explicit in this lattice calculation.

Use inverse temperature ${\beta }_T=(k_{B}T{)}^{-1}$, reserving $\beta$ without a subscript for the critical exponent. Write $c=qJ$; each undirected bond occurs once, so the number of bonds is $Nq/2$. In the mean-field theory of the Ising model, write each Ising spin as ${\sigma }_i=m+\delta {\sigma }_i$ and neglect the product of fluctuations:The resulting independent-spin statistical Hamiltonian isThe constant corrects the double counting of interaction energy. The Ising spin sums now factorize and can all be evaluated:This is the approximate partition function at an assumed mean field; equilibrium fixes $m$ self-consistently. The corresponding Ising auxiliary mean-field free energy isIts equilibrium value gives the mean-field Helmholtz free energy in the imposed field. Away from a stationary point its parameter $m$ is an assumed field variable, not necessarily the actual mean Ising spin of that independent-spin distribution.

For the full small-$m$ expansion at fixed $h$, set $t_h=\tanh ({\beta }_{T}h)$ and $s_h={\mathrm{sech}}^2({\beta }_{T}h)$. Differentiating $\log \cosh z$ gives successive derivatives $t_h$, $s_h$, $-2t_h{s}_h$, and $2s_h(2-3s_h)$ at $z={\beta }_{T}h$. ThereforeAt zero field, spin inversion symmetry eliminates odd powers and this simplifies toThe quadratic coefficient changes sign and the quartic coefficient is positive atThus the zero-field mean-field prediction is a continuous, continuous phase transition: the stable zero spin magnetization develops two symmetry-related nonzero minima continuously below $T_c$.

Both requested routes give the same mean-field self-consistency equation. First, the one-spin expectation in the effective field isSecond, differentiating the Ising auxiliary mean-field free energy giveswhose stationary condition is precisely that mean-field self-consistency equation. Choose its stable, lowest-free-energy branch rather than every algebraic solution.

An equally useful mean-field approximation parametrizes the trial distribution by its actual mean Ising spin. Its probabilities are $(1\pm m)/2$, giving the Bragg-Williams free energy of the Ising modelThis energy-minus-entropy function has expansionIts stationary equation is $h=k_{B}T\mathrm{arctanh}m-cm$, again equivalent to the mean-field self-consistency equation. The two functions differ away from equilibrium, but agree on stationary branches: for $z={\beta }_T(h+cm)$ and $m=\tanh z$, the entropy bracket equals $mz-\log (2\cosh z)$. Their small-$m$ coefficients therefore need not agree at arbitrary temperature; at $T_c$ they have the same leading critical quartic coefficient $c/12$. This distinction prevents confusing the auxiliary-field expansion with the physical-magnetization variational expansion.

For the order-parameter critical exponent, expand the equation of state at $h=0$:On a nonzero stable branch,Hence $\beta =1/2$. The spontaneous magnetization is understood by selecting a branch with an infinitesimal field after the thermodynamic limit; a finite symmetric sample has zero exact zero-field mean Ising spin.

For the magnetic susceptibility per site, define $\chi =\partial m/\partial h$ with $h$ in energy units. Implicit differentiation givesAbove $T_c$, $m=0$ and $\chi =1/[k_B(T-T_c)]$. Below $T_c$, evaluated on a selected ordered branch, the expansion of $m^2$ gives $\chi \sim 1/[2k_B(T_c-T)]$. Thus the magnetic-susceptibility critical exponent is $\gamma =1$ on both sides, with different amplitudes. At $T=T_c$, the equation of state becomesso the critical-isotherm exponent is $\delta =3$. These are mean-field critical exponents, not a claim that neglecting fluctuations gives the exact Ising transition in every dimension.

Discuss the zero-field phase boundaries, taking $h=0$ and writing $r={\mathcal{A}}_2$, $q={\mathcal{A}}_4$, $s={\mathcal{A}}_6>0$. This qualification matters because a conjugate field generally rounds the zero-field continuous transition. Nonzero stationary points satisfy $r+q{m}^2+s{m}^4=0$.

For $q>0$, the continuous phase transition occurs atOn the $r>0$ side the only stable phase is $m=0$; on the $r<0$ side the equilibrium global minima are $\pm m_0$ with $m_0^2\sim -r/q$. The magnetization therefore grows continuously from zero as $r$ changes sign.

For $q<0$, the first-order phase transition instead occurs while $r$ is positive. At phase coexistence, $m=0$ and $m=\pm m_0$ are equally deep global minimum points. With $a=m_0^2>0$, stationarity and equality of free energies giveEliminating $r$ gives $-q{a}^2/4-s{a}^3/3=0$, henceAt this coexistence value, the sextic even Landau potential factors as $f(m)=s{m}^2(m^2+3q/(4s){)}^2/6\ge 0$. This verifies that the three equal stationary values are global minima. Above this coexistence value of $r$ the disordered phase wins; below it the ordered phase wins, with a finite jump in the phase-selected magnetization. The values $r=0$ and $r=q^2/(4s)$ are the disordered and ordered spinodal points, respectively, not the coexistence boundary: the latter follows from the repeated root of $r+qa+s{a}^2=0$.

The two transition loci meet at the tricritical point,Two independent control parameters such as $T$ and $g$ can tune these two conditions. In a generic local phase diagram their coefficient map has a Jacobian matrix of rank two at the intersection; unspecified functions ${\mathcal{A}}_2(T,g)$ and ${\mathcal{A}}_4(T,g)$ do not by themselves determine the geometric orientation of the boundaries. The sextic term stabilizes this intersection and the negative-quartic first-order side.

At an ordinary zero-field continuous phase transition, the positive quartic coefficient remains nonzero while ${\mathcal{A}}_2\to 0$. Approaching from the ordered side therefore sends the variable $X$ in the tricritical crossover scaling form to $+\infty$, rather than to zero. With $Y=0$ the dimensionless minimizer obeys ${\psi }^4+X{\psi }^2-1=0$, soAt that stationary point, substitution gives ${\Phi }_{-}(X,0)=-{\psi }^2/4-{\psi }^6/12$. Its large-quartic asymptotic of tricritical scaling is consequently ${\Phi }_{-}(X,0)\sim -1/(4X)$, andIf ${\mathcal{A}}_2$ crosses zero linearly with temperature, this is quadratic in $T-T_c$, giving a finite heat capacity jump rather than a power divergence. ThereforeFor locally constant ${\mathcal{A}}_4>0$ and slope $a={\partial }_T{\mathcal{A}}_2{|}_{T_c}$, the leading singular heat capacity jump from this term is $\Delta C=T_c{a}^2/(2{\mathcal{A}}_4)$. Treating $\Phi$ as a nonzero constant while $X\to \infty$ would incorrectly assign the tricritical exponent to this ordinary transition.

Use the usual homogeneous, local, analytic Landau-Ginzburg theory; rotational symmetry and Z2 symmetry alone do not specify locality, analyticity, or the signs of the coefficients. The local derivative expansion contains every rotational scalar with an even total number of fields:The omitted terms include all allowed higher field powers and contracted spatial derivatives. Terms differing by a total derivative are equivalent for bulk behavior with suitable boundary conditions. The displayed gradient energy assumes spatial homogeneity. If only rotations about an origin are imposed, rotationally invariant position-dependent coefficients are also possible.

In the mean-field approximation, take a spatially uniform order parameter and minimize its Landau free energy. For $\gamma >0$, $u>0$ and $r=a(T-T_c)$ with $a>0$, keeping the quadratic and quartic terms givesThe single minimum splits into two minima as the temperature passes below the critical temperature. Choosing one minimum gives spontaneous symmetry breaking of the Z2 symmetry and a continuous phase transition. The order parameter is the expectation in a selected ordered state; averaging equally over both states gives zero. A Z2 symmetry does not guarantee a continuous phase transition: a negative quartic coefficient stabilized by a positive sextic coefficient can instead give a first-order phase transition.

Relaxing the Z2 symmetry permits odd terms. A nonzero linear source generally rounds this ordinary continuous phase transition. A cubic term can instead give a first-order phase transition: for $V(\varphi )=r{\varphi }^2/2+w{\varphi }^3/3+u{\varphi }^4/4$, $u>0$ and $w\ne 0$, equal minima occur atThe order parameter jumps from zero to ${\varphi }_{\ast }$. A continuous phase transition without an exact Z2 symmetry is still possible with additional tuning; odd terms do not logically exclude it.

For the multicritical even Landau potential, set $r={\mu }^2=a(T-T_c)$ and assume an integer $n\ge 2$, ${\lambda }_n>0$ and $\gamma >0$. For $n>2$, all lower even interactions must be absent or tuned to zero to obtain these mean-field critical exponents. The stationary equation isIts nonzero solutions are minima because $V^{\prime \prime }({\varphi }_0)=-(2n-2)r>0$. Thus the order-parameter critical exponent isSubstituting into the minimum free energy density givesThe singular heat capacity per unit volume is $c_s=-T{\partial }_T^2{f}_{\min }$. Smooth variation of the other coefficients does not change the leading power. ConsequentlyFor $n=2$, this means a finite heat capacity jump, with ${\alpha }_{\mathrm{MF}}=0$; for $n=3$, the tricritical point has ${\beta }_{\mathrm{MF}}=1/4$ and ${\alpha }_{\mathrm{MF}}=1/2$. The case $n=1$ would contain only a quadratic potential and would not stabilize an ordered phase.

The Gaussian fluctuation correction near a critical point is more singular than the mean-field approximation contribution whenThis is the upper critical dimension of an even scalar interaction. Equivalently, its coupling has engineering dimension $2n-(n-1)d$, positive below $d_c$. Dominant fluctuations signal breakdown of the mean-field approximation; the Gaussian expression is not then a calculation of the exact interacting heat-capacity critical exponent. At $d=d_c$, equality of the powers allows logarithmic corrections.

Write $r_i={\mu }_i^2$ and assume $g>0$, as required for quartic stability. In the mean-field approximation, the two-component radial quartic Landau potential isSet $x={\varphi }_1^2\ge 0$, $y={\varphi }_2^2\ge 0$. At fixed $\rho =x+y$, the quartic term is fixed and the quadratic term is minimized by assigning all of $\rho$ to the component with smaller $r_i$. Minimizing over $\rho$ then givesThe resulting phase diagram has these regions:

On $r_1=r_2=r<0$, the symmetry is enhanced to the orthogonal group $O(2)$, and the minima form the circleMixed-component minima occur on this line, rather than in an open coexistence region. Crossing the negative diagonal swaps the selected component discontinuously: this is a first-order phase transition, also called a spin-flop transition. For example, holding $r_2=r<0$ fixed and moving $r_1$ across $r$, the derivative $\partial f_{\min }/\partial r_1$ jumps from $-r/(8g)$ to zero. The angular degeneracy exactly on the line does not remove this discontinuity.

The positive half-axis $r_1=0$, $r_2>0$ is a continuous phase transition, as is $r_2=0$, $r_1>0$. Their endpoint at $(0,0)$ is a bicritical point, where the two continuous boundaries meet the first-order negative diagonal. Approaching the origin radially makes the ordered amplitude vanish continuously. The negative half-axes lie inside ordered phases and are not additional transition lines. There are two continuous half-axis boundaries and one first-order negative diagonal, meeting at the origin.

Take $b>0$ for a stable quartic Landau free energy. The cubic-term removal in a conserved quartic Landau free energy uses $\varphi =\psi -c/(3b)$ and gives, apart from a constant,For a closed binary fluid mixture, the compositional order parameter has fixed integral. Adding $d\int \varphi$ therefore adds the same constant to every allowed state. Equivalently, in the common-tangent construction for phase coexistence, a linear term changes the tangent slope, hence the chemical potential, by $d$, but leaves the coexistence compositions unchanged. This is why linear composition bias leaves coexistence compositions unchanged; it would not mean an unchanged composition at a fixed externally imposed chemical potential.

In the shifted variable, the binodal compositions are ${\psi }_{\pm }=\pm \sqrt{-a^{\prime }/b}$ for $a^{\prime }<0$. They merge atThe continuous phase transition therefore has the same scalar quartic Landau theory as the symmetric binary fluid mixture, with a shifted critical composition and control parameter. In the mean-field approximation, the composition difference vanishes with order-parameter critical exponent $1/2$. This argument presumes the critical composition lies in the physically accessible composition range; conservation affects the dynamics, not this static phase coexistence construction.

Within this analytic quartic Landau theory, equilibrium of a nonconserved scalar order parameter requires $f^{\prime }=0$. At a continuous phase transition, its restoring curvature vanishes, so $f^{\prime \prime }=0$. A nonzero $f^{\prime \prime \prime }$ would leave a cubic leading term in the Taylor expansion about that state and prevent it from being a local minimum. Thus a stable critical endpoint of a quartic Landau free energy requires all three derivatives to vanish, with $f^{\prime \prime \prime \prime }>0$. The condition $f^{\prime \prime }=0$ alone can instead identify a spinodal point.

For $f=\bar{b}{\lambda }^4+\bar{c}{\lambda }^3+\bar{a}{\lambda }^2+\bar{d}\lambda$, $\bar{b}>0$, the successive equations areThey giveAt these values, the free energy is exactly $f({\lambda }_c)+\bar{b}(\lambda -{\lambda }_c{)}^4$, proving stability. With no electric field, $\bar{d}=0$, so necessarily $\bar{c}=0$, followed by ${\lambda }_c={\bar{a}}_c=0$.

For $\bar{c}<0$ and $\chi >0$, the electric-field coupling to nematic order supplies the required negative ${\bar{d}}_c$ atTuning the quadratic coefficient to ${\bar{a}}_c$ simultaneously reaches the field-induced critical endpoint of the isotropic-nematic transition. The PDF's $\tilde{c}$ in the last clause is undefined; it is interpreted here as the previously defined $\bar{c}$. The endpoint is a prediction of the stable quartic, small-field mean-field approximation; higher-order terms or a field beyond the expansion's validity can shift it.

Average the trial field over an integer number of periods. The squared-field, squared-gradient and squared-Laplacian averages are $A^2/2$, $q_0^2{A}^2/2$ and $q_0^4{A}^2/2$, while $⟨{\varphi }^2|\nabla \varphi {|}^2⟩=A^4{q}_0^2/8$. Hence the single-mode smectic mean-field free energy isFor the physically stabilizing case $B>0$, minimization givesThe amplitude tends continuously to zero at $a_c$, so this mean-field approximation predicts a continuous phase transition. The ordered minimum has free-energy density $-(a_c-a{)}^2/(8B{q}_0^2)$.