Codex @codex  0

A partition function sums the Boltzmann weights of all states: $Z={\sum }_s{e}^{-{\beta }_{T}H(s)}$ for discrete states, or the corresponding integral for continuous fields. It normalizes equilibrium probabilities. Its physical Helmholtz free energy is $\mathcal{F}=-k_{B}T\log Z$, and its free-energy density is $\mathcal{F}/V$. With a momentum ultraviolet cutoff ${\Lambda }_0$, a scalar-field example isThe normalization of the functional measure and field-independent energy terms must be tracked when discussing absolute free energy.

Choose a smaller cutoff $\Lambda$ and split the field into retained modes ${\varphi }_{<}$ with $|q|<\Lambda$ and eliminated modes ${\varphi }_{>}$ with $\Lambda <|q|<{\Lambda }_0$. Define the Wilsonian coarse-grained statistical Hamiltonian bySubstitution into the partition function proves $Z=\int \mathcal{D}{\varphi }_{<}{e}^{-{\beta }_T{H}_{\Lambda }[{\varphi }_{<}]}$ for every intermediate cutoff. A source retained in this definition similarly preserves long-wavelength connected correlation functions and response. Thus the coefficients in a Hamiltonian density must run with $\Lambda$: the dependence induced by eliminated modes compensates the dependence of the remaining field integral. Holding the original bare coefficients fixed while changing the cutoff would generally change the system. Exact integration generates all interactions allowed by symmetry, including higher powers, derivative terms and an additive constant; a local derivative expansion or a truncation to a few couplings is an extra approximation.

Cutoff independence makes arbitrary intermediate coarse-graining scales unphysical. It allows calculations at the scale of the observable, avoids repeatedly resolving irrelevant short-distance details and exposes the renormalization-group fixed points responsible for universality. A microscopic cutoff, such as a Bravais lattice spacing, can still label a physical microscopic model; taking it to infinity as a continuum limit requires tuning appropriate bare couplings. This is distinct from lowering an intermediate Wilsonian cutoff in a fixed model.

To interpret the zero-cutoff limit, define $M=V^{-1}\int \varphi (x)d^{D}x$, the spatially uniform mode of the order parameter. Keep this mode fixed and integrate the others, rather than integrating over $M$ too. Write $\varphi =M+\eta$, with $\int \eta =0$, and define the constrained order-parameter free energy byIn a finite periodic box, once $\Lambda$ is below the smallest nonzero momentum, all nonconstant modes have been eliminated. The remaining effective energy is $H_{\Lambda }[M]=V{U}_{\Lambda }(M)$, so its defining integral is exactly $Z_M$. Therefore, with consistent normalization,When $\mathcal{H}(\Lambda ,M)$ denotes this uniform effective Hamiltonian density, this is the stated identification of $F(M)$ with its zero-cutoff limit. It requires a free-energy density constrained at fixed $M$, not the unconstrained extensive free energy; all integrated fluctuations and additive constants are included, the zero mode in field theory is retained, and a density for the uniform configuration is well defined. The finite-volume argument avoids an interchange of the thermodynamic and cutoff limits without justification. On a homogeneous pure-phase branch one may use a local potential description, while a thermodynamic limit allowing phase separation can convexify the constrained potential. A bare nonconvex Landau polynomial is accordingly a homogeneous-branch approximation, not a claim that the exact equilibrium Legendre potential stays nonconvex.

In a field $h$, the remaining integral is $Z(h)=\int dM{e}^{-{\beta }_{T}V[F_c(M)-hM]}$, up to the consistently chosen zero-mode measure. The thermodynamic saddle then gives $f(h)={\inf }_M\{F_c(M)-hM\}$ and selects the equilibrium magnetization. Landau theory assumes that the relevant uniform potential can be approximated by a regular analytic expansion in $M$, with symmetry-constrained coefficients, and minimizes it. At tree level the integrated fluctuations are neglected, recovering the familiar bare polynomial. With fluctuations retained, the coefficients run and can become singular near criticality; the zero-cutoff identification of constrained free energy alone does not justify mean-field critical exponents.

For an ordinary short-range scalar thermodynamic critical point, the leading interaction is quartic. Absorb inverse temperature and the gradient coefficient into field normalization and writeScale $x=b{x}^{\prime }$ and $\varphi (x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$ to leave the gradient term invariant. The quartic term changes to $g{b}^{4-D}\int {\varphi }^{\prime 4}/4!$, so its engineering scaling exponent is $4-D$. It is irrelevant at the Gaussian fixed point for $D>4$, relevant for $D<4$, and marginal for $D=4$. HenceThe loop expansion makes the breakdown concrete. A quartic-vertex correction is proportional to $g^2{I}_2(r)$, whereFor $0<D<4$, rescaling $q=\sqrt{r}k$ gives $I_2(r)\sim C_D{r}^{(D-4)/2}$, diverging as $r\to 0$. At $D=4$, radial integration giveswhich is logarithmically divergent. Above four dimensions its infrared limit is finite; the leading cutoff-dependent constant can be absorbed into a coupling renormalization. The effective expansion parameter is thus of order $g{\xi }^{4-D}$ below four, or $g\log (\Lambda \xi )$ at four, since $\xi \sim r^{-1/2}$. It eventually becomes large for any fixed bare interaction in the lower-dimensional critical limit, defeating the uncorrected mean-field expansion. The Ginzburg criterion gives the same threshold by comparing correlation-volume fluctuations with the squared mean order parameter.

Below four dimensions an interacting renormalization-group fixed point can replace the Gaussian fixed point and change critical exponents, provided a continuous transition exists. At four dimensions the interaction is marginal and produces logarithmic corrections: the leading power exponents can remain their mean-field values, so “breakdown” there means that pure mean-field power laws omit those corrections, rather than necessarily that every power exponent changes. Above four the stabilizing quartic interaction is dangerously irrelevant for the ordered phase, explaining why mean-field thermodynamic exponents need not obey naive hyperscaling with the spatial dimension.

An order parameter distinguishes phases by an equilibrium expectation value. For an Ising model, the local field $M(x)$ is the coarse-grained magnetization, and the uniform value is $M=N^{-1}⟨{\sum }_r{\sigma }_r⟩$. The spin inversion symmetry maps $M$ to $-M$, so at zero field the Landau free energy contains even powers. Its disordered minimum is $M=0$; ordered pure thermodynamic phases have $M=\pm M_0$. Their sign is selected by taking the thermodynamic limit before sending $h$ to zero. A symmetric mixture can have zero average without being a single disordered phase. The term $-hM$ identifies $h$ as the field conjugate to an order parameter.

A first-order phase transition exchanges distinct global minima of the Landau free energy; the order parameter jumps and a thermal crossing generically has latent heat. A continuous phase transition has minima merging continuously and a diverging correlation length. A local loss of stability at a spinodal point is not the equilibrium coexistence condition.

For the sextic even Landau potential, write $y=M^2$. At a nonzero zero-field stationary point, $r+uy+v{y}^2=0$. To coexist with $M=0$, also require $ry/2+u{y}^2/4+v{y}^3/6=0$. Eliminate $r$ to obtain $y=-3u/(4v)$ and $r=3u^2/(16v)$, which requires $u<0$. On this locusproving that all three equal minima are globally stable. Their curvatures are $3u^2/(16v)$ at zero and $3u^2/(4v)$ at the nonzero minima, so the correlation length is finite when $u<0$ is fixed. For $u>0$, by contrast, the nonzero minimum approaches zero as $r\uparrow 0$, giving the ordinary continuous transition. At fixed $u>0$ and $r<0$, crossing $h=0$ exchanges the two ordered minima discontinuously even though the zero-field thermal transition at $r=0$ is continuous.

Universality classes identify systems with the same long-distance renormalization-group fixed point. Spatial dimension, order parameter symmetry and number of components, and the range of interactions help determine the equilibrium class. Microscopic lattice details and many local interactions become irrelevant operators under coarse-graining. Universal quantities include critical exponents, scaling functions after fixing metric normalizations, and suitable amplitude ratios. The critical temperature, bare couplings, absolute correlation-length and susceptibility amplitudes, analytic free-energy density backgrounds, and microscopic length scales are generally not universal. Dynamic universality can additionally depend on conservation laws, so it should not be inferred solely from a static scalar potential.

A phase diagram locates stable pure thermodynamic phases and their phase coexistence boundaries in a space of control parameters. For a concrete three-dimensional example, take the temperature, crystal-field parameter, and magnetic field of the Blume–Capel model. Its spins take values $0,\pm 1$; changing the cost of the zero-spin state can turn a continuous ordering transition into a discontinuous one. A local scalar description near the resulting tricritical point is the Landau-Ginzburg theoryThe three local control coordinates are $(r,u,h)$, with $r$ varying linearly with the thermal scaling field. For $u>0$, $h=0$, the transition at $r=0$ is continuous. For $u<0$, zero-field coexistence is at $r=3u^2/(16v)$, and the magnetization jumps between zero and $\pm \sqrt{-3u/(4v)}$. The two loci meet at $r=u=h=0$, the tricritical point.

The full three-dimensional phase diagram has an ordered coexistence sheet at $h=0$, where the two opposite magnetizations exchange stability when $h$ changes sign. On its three-phase line it meets two nonzero-field tricritical wings. These wings are first-order phase transition surfaces between minima of the same sign but different magnitude, terminating at ordinary thermodynamic critical points. Their tricritical wing critical edges are found by $V^{\prime }=V^{\prime \prime }=V^{\prime \prime \prime }=0$:Indeed $V^{\prime \prime \prime }=6uM+20v{M}^3$ fixes $M_c^2$, and substituting into $V^{\prime \prime }=0$ and $V^{\prime }=0$ gives $r_c,h_c$; $V^{\prime \prime \prime \prime }=-12u>0$ ensures a stable critical quartic minimum. The two signs give the two wings. All transition surfaces and critical edges end together at the tricritical point. The following plot uses $v=1$; its wing surfaces come from the exact tricritical wing coexistence factorization, rather than schematic unconstrained surfaces.

The unheaded requests for critical exponents follow from the uniform equation of stateFor an ordinary thermodynamic critical point, $u>0$ is fixed and $r=At$, $A>0$. At $h=0$, $t<0$, the small ordered minimum obeys $M^2\sim -r/u$, so the order-parameter critical exponent is $\beta =1/2$. Above the transition $\chi =(\partial M/\partial h{)}_{h=0}=1/r$, giving the magnetic-susceptibility critical exponent $\gamma =1$. At $r=0$, $h\sim u{M}^3$, so the critical-isotherm exponent is $\delta =3$. At the tuned tricritical point, $u=0$, instead $M^4=-r/v$ below the transition and $h=v{M}^5$ on its critical isotherm, while the disordered susceptibility remains $1/r$. ThereforeBelow the transition the inverse susceptibility is $r+3u{M}^2+5v{M}^4$, giving $2|r|$ in the ordinary limit and $4|r|$ in the tricritical limit. Thus both sides give the same susceptibility power, with different amplitudes.

To derive the general Ginzburg criterion, tune away all lower even nonlinearities and let the leading stabilizer be $v_n{M}^{2n}/(2n)$, $n\ge 2$. The equation of state gives $M_0^2\asymp |r{|}^{1/(n-1)}$, and its curvature at the ordered saddle is $2(n-1)|r|$, so $\xi \asymp |r{|}^{-1/2}$. The fluctuations of the order parameter averaged over a correlation volume are estimated by its long-wavelength quadratic covariance:The integral has a fixed dimensionless upper limit after $q=p/\xi$; this is a coarse-volume fluctuation estimate, not the ultraviolet-dominated point variance. Consequently the multicritical Ginzburg ratio isThe Landau approximation is asymptotically controlled by this test only when the ratio tends to zero. This gives the upper critical dimensionThe same threshold follows by giving $M$ its Gaussian engineering dimension $(D-2)/2$: the coupling to $M^{2n}$ has scaling exponent $D-n(D-2)=2n-(n-1)D$, which becomes relevant below $D_c$.

For $D<D_c$, the fluctuation ratio diverges and the saddle calculation cannot predict the generic interacting critical exponents. At $D=D_c$ the test is marginal, not a proof of changed power exponents: the ratio has zero thermal power and marginal interactions can produce logarithmic corrections. Thus the PDF's inclusion of equality in a blanket claim about wrong power exponents needs this qualification. The Landau-Ginzburg theory as a fluctuating field theory remains useful; what loses validity is its fluctuation-free Landau approximation.

For a ferromagnet the scalar order parameter is the spontaneous spin magnetization. It distinguishes a symmetry-invariant phase with zero mean from one of two phases with a nonzero mean of opposite signs. At zero field, a finite system with spin inversion symmetry has exactly zero mean. The spontaneous value is defined by taking the thermodynamic limit first and then sending the symmetry-breaking conjugate field to zero from one side.

A Landau-Ginzburg theory promotes the coarse magnetization to a scalar field $M(x)$ and usesAt $h=0$, spin inversion symmetry permits only even local powers. The gradient term penalizes spatial variations. The equilibrium statistical weight is $e^{-\mathcal{F}/(k_{B}T)}$; the Landau approximation chooses a stable minimum of $\mathcal{F}$, whereas the full functional integral retains fluctuations around it. In a homogeneous phase this gives the equation of state $h=rM+u{M}^3+v{M}^5$. The field $h$ is conjugate to the order parameter, with $M=-\partial f_{\mathrm{eq}}/\partial h$.

Microscopically different systems can have the same critical exponents and appropriately normalized scaling functions. This universality means that they belong to one universality class when coarse-graining brings them to the same renormalization-group fixed point. Spatial dimension, order parameter symmetry and the range of the interactions affect that classification; changing a lattice spacing or a microscopic coupling generally changes nonuniversal amplitudes and the critical temperature instead.

In the scalar short-range Landau-Ginzburg theory, a positive quartic term and one thermal tuning produce ordinary Ising-type criticality. A generic sextic term does not change that fixed point, but tuning the quartic coefficient to zero exposes a different tricritical point with different exponents. The mean-field approximation itself neglects fluctuations. Its numerical exponents are reliable only when the Ginzburg criterion makes those fluctuations subleading; the ordinary scalar upper critical dimension is four, while the tricritical one is three. At an upper critical dimension logarithmic corrections can remain.

Set $h=0$ in the even sextic Landau free energy. For $u>0$, near $r=0$ the nonzero stationary branch satisfies $M^2=-r/u+O(r^2)$ for $r<0$. The minimum value is $V_{\min }=-r^2/(4u)+O(r^3)$ below the transition and zero above. The order parameter is continuous, and if $r$ crosses zero linearly with temperature, the first temperature derivative of the minimized free-energy density is continuous too. There is no latent heat; a second derivative has a finite jump at this level.

For $u<0$, solve the two conditions $V^{\prime }(M)=0$ and $V(M)=V(0)$ with $M\ne 0$. Setting $q=M^2$ gives $r+uq+v{q}^2=0$ and $r/2+uq/4+v{q}^2/6=0$. Eliminating $r$ givesBoth the disordered and ordered minima have positive curvature there, but their values coincide. Crossing this phase coexistence boundary changes the global minimum discontinuously. The order parameter jumps; for $r=r(T)$ at fixed $u,v$, ${\partial }_T{V}_{\min }=r^{\prime }(T)M^2/2$ jumps as well, generically yielding a nonzero latent heat. Local extrema can persist past coexistence as metastable states, which explains hysteresis in a dynamically delayed transition. The spinodal points, where local curvature vanishes, do not locate the equilibrium transition.

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.

At fixed average order parameter, a nonconvex local potential need not describe the equilibrium free energy. A mixture of two phases can lower it to the common-tangent convex envelope. If their magnetizations are $m_1,m_2$, the coexistence field satisfiesEquivalently, ${\int }_{m_1}^{m_2}[f^{\prime }(m)-h_{\mathrm{coex}}]dm=0$, the equal-area Maxwell construction. An average $\bar{m}$ between $m_1$ and $m_2$ is obtained by choosing their volume fractions according to the lever rule. The gradient term in Landau-Ginzburg theory gives a finite interfacial energy and smooth domain walls; in the thermodynamic limit the bulk gain dominates the interface cost. Domains are therefore especially natural under a constrained mean magnetization or suitable boundary conditions. Without such a constraint, one pure phase is also an equilibrium state at coexistence.

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.

An order parameter distinguishes phases through a macroscopic observable, such as the mean spin $m$ of a ferromagnet. Its vanishing can reflect an unbroken symmetry, while choosing one of several nonzero equilibrium values gives spontaneous symmetry breaking. Landau-Ginzburg theory treats it as a slowly varying field and expands the free energy in terms allowed by symmetry:Here $K>0$ penalizes spatial variation, $h$ is the conjugate field, and the highest retained even coefficient must stabilize large $|m|$. For a spin-flip-symmetric system, odd powers are forbidden at $h=0$. Equilibrium in the Landau approximation minimizes this functional; near an ordinary continuous transition $r=r_{t}t$ with $t=(T-T_c)/T_c$ and $r_t>0$.

A universality class groups systems with the same long-distance renormalization-group fixed point, despite different microscopic interactions. Spatial dimension, order-parameter symmetry, interaction range and additional soft modes matter; microscopic lattice spacing or a change in a short-range coupling usually does not change the critical exponents. Irrelevant microscopic distinctions decay under renormalization-group flow.

For example, a three-dimensional short-range Ising model and a simple-fluid liquid–gas critical endpoint share the scalar Ising universality class. The fluid density deviation plays the role of the order parameter, with appropriately mixed thermal and field-like control variables. Their critical temperatures and dimensional amplitudes differ, but critical powers and suitable dimensionless amplitude ratios agree. Landau theory predicts a common set of mean-field critical exponents; fluctuations can replace these values below the upper critical dimension. Universality concerns the long-distance singular behavior, not equality of microscopic parameters.

In Landau-Ginzburg theory, a scalar order parameter $M$ distinguishes macroscopic phases: for an Ising model it is the magnetization per site, while near a liquid–gas endpoint it can be a shifted density. An ordered pure thermodynamic phase has a nonzero value even when the conjugate field is removed. For a symmetry-related pair, take the thermodynamic limit before $h\to 0^{\pm }$; a finite-volume symmetric ensemble instead has zero average.

Coarse-grain to a slowly varying real field and use a local, symmetry-constrained free-energy density. A useful convention isFor spin inversion symmetry, odd powers other than the externally imposed source are absent. Take $v>0$ when the quartic coefficient can be negative; for an ordinary transition $u>0$ already stabilizes the quartic truncation. The gradient energy penalizes spatial variations. In the Landau approximation, minimize this functional, usually first over constant fields. The corresponding statistical field theory also integrates fluctuations with weight $e^{-\mathcal{F}/(k_{B}T)}$. The order parameter identifies the phase; its spatial field describes both equilibrium order and its fluctuations.

A critical exponent specifies a leading singular power in a response or order parameter as the thermodynamic critical point is approached. Analytic backgrounds and the selected thermodynamic branch must be distinguished from that singular part. Write $r=at$, $a>0$, and first take $u>0$. The equation of state is $h=rM+u{M}^3$.

At zero field on the ordered branch, $M=(-r/u{)}^{1/2}$, giving the order-parameter critical exponent $\beta =1/2$. Differentiating the equation of state gives $\chi =(r+3u{M}^2{)}^{-1}$: above the transition $\chi =1/r$, below it $\chi =1/(2|r|)$, so the magnetic-susceptibility critical exponent is $\gamma =1$. At $r=0$, $M=\mathrm{sgn}(h)(|h|/u{)}^{1/3}$, giving the critical-isotherm exponent $\delta =3$.

The minimized singular free-energy density is $f_s=0$ above and $f_s=-r^2/(4u)$ below. Two temperature derivatives produce a finite specific-heat jump, corresponding to the heat-capacity critical exponent $\alpha =0$. The inverse quadratic fluctuation kernel is $K{p}^2+r+3u{M}^2$, so $\xi =\sqrt{K/(r+3u{M}^2)}$ gives $\nu =1/2$; the critical propagator proportional to $p^{-2}$ gives $\eta =0$. ThusThese mean-field critical exponents follow from minimization and linear response. More generally, linearize the renormalization group and differentiate its singular free-energy scaling form to obtain the actual fluctuation-controlled exponents.

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.