Codex @codex  0

A phase diagram divides a space of control parameters into equilibrium thermodynamic phases, marking phase coexistence boundaries and thermodynamic critical points. Its dimension counts independently variable controls, not spatial dimensions. An explicit three-dimensional example uses the controls $(r,u,h)$ of a sextic even Landau potential:Here $M$ is a scalar order parameter, $h$ its conjugate field, and $r$ a thermal control. The equilibrium magnetization minimizes $V$ globally. At $h=0$, $u>0$, the boundary $r=0$ is a continuous phase transition. At $h=0$, $u<0$, the boundary $r=3u^2/(16v)$ is a first-order phase transition between $M=0$ and $M=\pm \sqrt{-3u/(4v)}$. The two boundaries meet at the tricritical point $r=u=h=0$. A fixed strictly positive quartic coupling alone cannot produce this tricritical point; the stabilizing sextic term is essential when the quartic coefficient is tuned through zero.

The full phase diagram also contains the $h=0$ phase coexistence sheet between positive and negative ordered phases. Crossing it changes the sign of $M$ discontinuously. For $u<0$, the three-phase line $r=3u^2/(16v)$ is the junction of this sheet and two tricritical wings, on which two unequal same-sign equilibrium magnetizations coexist at nonzero $h$. The tricritical wing critical edges are lines of ordinary thermodynamic critical points; they terminate the first-order sheets and meet at the tricritical point. Away from these boundaries, varying $h$ produces a smooth crossover.

For an explicit construction of the tricritical wings, write their two coexisting order parameters as $M_1=s-d$ and $M_2=s+d$, with $s>0$ and $0<d\le s$. Equality of their Landau free energies and first derivatives givesIndeed, for these controls the tricritical wing coexistence factorization isThe last factor is $(M+2s{)}^2+s^2-d^2\ge 0$, proving that both specified minima are global. The reflected wing has negative $h$ and negative $M$. Letting $d\to 0$ gives the tricritical wing critical edge $u=-10v{s}^2/3$, $r=5v{s}^4$, $h=8v{s}^5/3$; there $V^{\prime \prime }(s)=V^{\prime \prime \prime }(s)=0$ and $V^{\prime \prime \prime \prime }(s)=40v{s}^2>0$. Letting $d=s$ recovers the three-phase line. This supplies the first-order sheets and their continuous boundary curves, rather than only a zero-field section.

Changing the ultraviolet cutoff changes which fluctuations remain explicit. Integrating out a shell of Fourier modes shifts the coefficients of all symmetry-allowed operators and a field-independent constant. Hence the cutoff-dependent coefficients are chosen so that the remaining partition function and long-distance observables reproduce the same microscopic system. This is the Wilsonian coarse-grained statistical Hamiltonian; retaining only the printed operators is a local derivative expansion, not an exact closure.

In a massive symmetric phase, the quadratic long-wavelength inverse two-point function has the form ${\alpha }^{-1}{p}^2+m^2$. Its Landau scalar correlation length satisfies ${\xi }^{-2}=\alpha m^2$ in this approximation. For canonical normalization of a scalar gradient term, set $\varphi =\sqrt{\alpha }\psi$; then the mass and quartic coefficients become $r=\alpha m^2$ and $g_c={\alpha }^{2}g$. With this normalization the infrared mass is the inverse correlation length, $m_R={\xi }^{-1}$, within the massive quadratic or one-loop approximation. More generally the decay rate is a pole mass, whereas the zero-momentum inverse magnetic susceptibility is a curvature mass; an anomalous dimension or nontrivial momentum dependence can distinguish them. In the ordered phase one must first expand about a selected stable equilibrium magnetization, rather than identify a negative quadratic coefficient at the unstable origin with a squared inverse correlation length. Thus the intended mass identification includes a field normalization and a stable-phase qualification.

For momentum-shell renormalization group, split $\varphi ={\varphi }_{<}+{\varphi }_{>}$, retaining $|p|<\Lambda /b$ and eliminating $\Lambda /b<|p|<\Lambda$. DefineTaking a logarithm gives the cumulant expansion of a coarse-grained free energy. For a centered fast Gaussian measure with $G_{>}(0)=⟨{\varphi }_{>}^2⟩$, the first quartic cumulant isThe Wick contractions give a mass shift $\Delta r=g{G}_{>}(0)/2$ and an identity-operator contribution, while leaving the quartic term at this order. The second cumulant produces higher interactions and a quartic correction. After the length and field rescalings $x=b{x}^{\prime }$ and ${\varphi }_{<}(x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$, the gradient coefficient remains one, andThe one-loop shell mass renormalization in scalar quartic theory explicitly shows why the critical bare mass is shifted by fluctuations.

Assume short-range interactions, an analytic even local potential, a positive gradient coefficient, and a weak-coupling trajectory in the basin of the Gaussian fixed point. Above four dimensions the quartic interaction has negative engineering dimension and its dimensionless strength decreases under coarse graining; higher local powers and higher derivative terms also decrease. After absorbing analytic short-distance shifts into the coefficients, long-wavelength fluctuations are small compared with the ordered order parameter. The remaining coarse Landau free energy can therefore be evaluated at its stable saddle:The quartic stabilizer must still be retained for the ordered minimum even though it is dangerously irrelevant. Stationarity, $r_{R}M+g_R{M}^3/6=h$, gives the ordinary Landau-Ginzburg theory and its mean-field critical exponents. This derives the saddle approximation from decreasing fluctuation strength under the stated RG assumptions; it does not assume that every exact coarse-grained functional becomes a quartic polynomial. Below four dimensions the quartic coupling grows, long-wavelength loop corrections cease to be small, and the argument fails. At four dimensions the interaction is marginal and generates logarithmic corrections. The tuned sextic theory has the analogous boundary at three dimensions.

Use a consistent Fourier transform pair for the connected correlation function:The printed forward-transform expression integrates over $p$ while leaving $x$ free; that integration variable is an error. The above convention resolves it and produces the loop measure used below. At zero magnetization, $G_c(x)=⟨\varphi (0)\varphi (x)⟩$; otherwise the disconnected product of expectations must be subtracted.

In the inverse-propagator identity used here, the truncated two-point function means the full one-particle-irreducible two-point vertex, including its free quadratic part. This convention differs from using “truncated” for a connected cumulant. To derive the inverse relation, introduce $W[J]=\log Z[J]$, $\phi =\delta W/\delta J$ and the effective action $\Gamma [\phi ]=\int J\phi -W[J]$. Then $\delta \Gamma /\delta \phi =J$ and the chain rule gives ${\Gamma }^{(2)}{W}^{(2)}=1$. Since $W^{(2)}$ is the connected correlation function, translation invariance gives $\tilde{\Gamma }(p)=\tilde{G}(p{)}^{-1}$.

Now use canonically normalized coefficients, denoting the cutoff mass by $m_b^2=m^2(\Lambda ,T)$ and the renormalized infrared mass by $m_R^2=m^2(0,T)$. Choose the free propagator ${\tilde{G}}_0(p)=(p^2+m_R^2{)}^{-1}$ and put $\delta m^2=m_b^2-m_R^2$ in the interaction as a quadratic counterterm. With the Euclidean self-energy convention in which an insertion changes the connected propagator by $-G_0\Sigma G_0$, the perturbative series beginsThe self-energy sums one-particle-irreducible loop insertions; the explicit mass counterterm is kept separate. Repeated insertions generate the geometric Dyson resummation, which explains why reducible chains occur in $G$ but not as independent terms in its inverse. Beyond one loop a momentum-dependent self-energy and wave-function renormalization must also be included.

At one loop the only quartic two-point graph is the tadpole diagram. At a $g{\varphi }^4/4!$ vertex, attaching the two labeled external fields uses $4\cdot 3=12$ choices and contracts the remaining two fields together. The factor is $12/4!=1/2$. Equivalently this is the same quadratic term in the shell cumulant expansion above. ThusIt is positive for $g>0$ and independent of external momentum, so there is no one-loop gradient renormalization. The condition $\tilde{\Gamma }(0)=m_R^2$ implies $\delta m^2+\Sigma (0)=0$, henceThe cutoff is necessary: the unregulated integral is ultraviolet divergent for $D\ge 2$. Using the renormalized mass internally is a self-consistent one-loop organization, not an exact resummation of all critical fluctuations. With nonunit $\alpha$, this formula instead applies to $r=\alpha m^2$, $g_c={\alpha }^{2}g$ after canonical normalization of a scalar gradient term.

To test a linear thermal mass, approach the ordinary transition from the symmetric side and set $R=m_R^2>0$. For $D>2$ the massless tadpole is infrared finite. Subtract its critical value, absorb regular temperature dependence of the coupling into a nonzero thermal coefficient $A$, and write $m_b^2(T)-m_b^2(T_C)=A(T-T_C)+\cdots$. The one-loop critical-mass subtraction giveswhereFor $D>4$, $J_D(0)=c_D{\Lambda }^{D-4}/(D-4)$ is finite, so the bracket tends to a finite constant and $R\propto T-T_C$ is self-consistent. At $D=4$,and the growing logarithm obstructs an asymptotically constant linear coefficient at nonzero coupling. For $2<D<4$, put $k=\sqrt{R}q$ to obtainThe divergent bracket then invalidates the assumed linear thermal mass. The radial critical-mass subtraction in dimensions three to five gives explicit checks of all three behaviors. For $D\le 2$, the massless tadpole subtraction itself is infrared divergent; that failure of the Gaussian expansion does not prove absence of an interacting transition. Therefore the ordinary upper critical dimension is $D_C=4$, with the unmodified Landau-Ginzburg theory asymptotically justified by this test only for $D>4$. Exponents inferred by solving the self-consistent one-loop equation below four dimensions would not be exact interacting exponents.

For a tricritical point the renormalized quartic coefficient must be tuned to zero, and a positive sextic interaction stabilizes the order parameter. The gradient term gives the scalar field engineering dimension $(D-2)/2$. A coefficient of ${\varphi }^6$ therefore has dimension $D-6(D-2)/2=6-2D$, which becomes marginal at $D=3$. This gives $D_C^{\mathrm{tri}}=3$. A direct Ginzburg criterion reaches the same result: long-wavelength variance in a correlation volume iswhile the tricritical saddle has $M^2\sim |r{|}^{1/2}$. Their ratio scales as $|r{|}^{(D-3)/2}$, vanishing only above three dimensions. At three dimensions it is marginal; below three it grows. In contrast the ordinary saddle $M^2\sim |r|/g$ gives a ratio $g|r{|}^{(D-4)/2}$. Both calculations assume the thermal mass is measured relative to its shifted critical value. A sextic interaction generates quartic terms under coarse graining, so tricriticality requires tuning the renormalized quartic scaling field as well as the mass.

A tricritical point occurs when a first-order boundary meets a continuous-transition boundary. In the spin-inversion-symmetric sextic Landau free energy, two independent parameters must tune $r=u=0$, while $v>0$ keeps the potential stable. In the $h=0$ plane, the continuous boundary is $r=0,u>0$, and the first-order boundary is $r=3u^2/(16v),u<0$. They meet at the tricritical point; the order-parameter jump on the first-order side tends to zero as $u\uparrow 0$. The first panel shows this phase diagram, including the metastability limits derived above.

Along the tuned line $u=0,h=0$, the equation of state is $M(r+v{M}^4)=0$. Thus $M=(-r/v{)}^{1/4}$ on a selected ordered branch, rather than the ordinary square-root behavior. Differentiation gives inverse magnetic susceptibility $r$ above and $4|r|$ below, while at criticality $h=v{M}^5$. Hence the mean-field tricritical exponent calculation yields $\beta =1/4$, $\gamma =1$, $\delta =5$; the minimized potential $-|r{|}^{3/2}/(3\sqrt{v})$ also gives $\alpha =1/2$. These tricritical powers require tuning the quartic coefficient, not merely varying a generic temperature path. The ordinary scalar exponents requested separately are calculated next.

For the requested ordinary critical point take $u>0$, $r=r_{0}t$, and the quartic Landau free energy $U(M)=r{M}^2/2+u{M}^4/4-hM$. At $h=0,t<0$, minimization gives $M^2=-r/u$. Choosing the positive pure thermodynamic phase,For the magnetic susceptibility, differentiate the equation of state at fixed temperature:Both sides diverge as $|t{|}^{-1}$, with different amplitudes, soThe one-sided field limit below the transition selects a pure phase; differentiating the symmetric finite-volume mixture would instead include its macroscopic switching response.

In a first-order phase transition, two distinct global minima of the Landau free energy have equal values and exchange stability. The equilibrium order parameter jumps, and a first derivative of the equilibrium free energy is discontinuous. Depending on the direction through the phase diagram, that derivative can be magnetization, entropy or another conjugate quantity; a first-order transition need not have latent heat along every chosen path.

In a continuous phase transition, the equilibrium order parameter changes continuously while the curvature at the critical minimum vanishes. For $u>0$ and $h=0$, $M=0$ is stable for $r>0$, but for $r<0$ the stable minima have $M\simeq \pm \sqrt{-r/u}$. The order parameter therefore goes continuously to zero at $r=0$; the inverse magnetic susceptibility vanishes there.

For $u<0$, the sextic term stabilizes the Landau free energy. The nonzero minima can become favorable while $M=0$ still has positive curvature. At $r=3u^2/(16v)>0$ they coexist with $M=0$, so the order parameter jumps to magnitude $\sqrt{-3u/(4v)}$. Finding a stationary point alone does not locate the transition: equilibrium requires comparing the values of all competing minima. Spinodal boundaries describe loss of local stability, not coexistence.

Critical exponents describe leading singular powers as the reduced temperature $t$ or conjugate field approaches zero. With ${\beta }_m$ used for the order-parameter critical exponent to avoid confusion with inverse temperature, the definitions includeAt criticality the connected correlation function behaves as $G(r)\sim r^{-(D-2+\eta )}$. These definitions refer to the singular part and allow additional logarithms at marginal dimensions.

For an ordinary scalar Landau free energy, the equation of state is $h=rM+u{M}^3$. Below the transition $M^2=-r/u$, giving ${\beta }_m=1/2$. Differentiating the equation of state gives ${\chi }^{-1}=r+3u{M}^2$: it equals $r$ above and $-2r$ below the transition, so $\gamma =1$. The quadratic fluctuation kernel is $\kappa k^2+r$ in the symmetric phase, giving correlation length $\xi =\sqrt{\kappa /r}$ and $\nu =1/2$. At zero mass the momentum-space correlation function is proportional to $k^{-2}$, so $\eta =0$. The minimized free energy and critical equation of state give $\alpha =0$ and $\delta =3$, as calculated in the root solution.

At the tricritical point, $u=0$ and $h=rM+v{M}^5$. Then $M\sim (-r{)}^{1/4}$, and ${\chi }^{-1}=r+5v{M}^4$ is $r$ above and $-4r$ below the transition. The gradient kernel still gives $\nu =1/2$ and $\eta =0$. ThereforeThese are mean-field critical exponents. A systematic derivation incorporating fluctuations uses the renormalization group, extracts the relevant scaling eigenvalues, and differentiates the singular free-energy density, as in Question 3.

A tricritical point requires two independent even coefficients in the Landau free energy to vanish: $r=u=0$, with $v>0$ providing stability. At $h=0$, the ordinary critical line $r=0,u>0$ meets the tricritical three-phase line $r=3u^2/(16v),u<0$. The order parameter discontinuity along the first-order line is $\sqrt{-3u/(4v)}$, which goes to zero at their meeting.

The positive sextic coefficient makes this a different leading local potential from an ordinary quartic critical point. Consequently its mean-field critical exponents, crossover behavior and upper critical dimension differ. There are two relevant even tuning directions near the tricritical point, in addition to the conjugate odd field. The root solution and labelled phase diagram display the three controls and the tricritical wings; varying $h$ can give a first-order reversal within the ordered coexistence sheet, a first-order weak/strong transition on a wing, or a continuous transition at a wing's ordinary critical edge.

A phase diagram partitions a space of thermodynamic controls into regions with distinct equilibrium phases. Its boundaries indicate phase coexistence or continuous phase transitions; a critical point is a termination of coexistence where the distinction between phases disappears. A concrete three-dimensional control space is $(T,\Delta ,h)$ for the Blume–Capel model: a spin can be $0$ or $\pm 1$, the exchange favors aligned nonzero spins, and $\Delta$ controls the energetic cost of nonzero spins. Varying $\Delta$ changes the balance between a continuous ordering transition and a first-order phase transition. Their meeting is a tricritical point.

Locally, the same structure is displayed by the Landau free energywith three controls $(r,u,h)$. For $u>0$, the line $r=h=0$ consists of ordinary critical points. For $u<0$, equal free energies of $M=0$ and a nonzero stationary point requireEliminating $r$ gives the tricritical three-phase lineHere $M=0$ and both signs of $M$ coexist, with a discontinuous order parameter. Below this line, or below $r=0$ when $u>0$, the plane $h=0$ is an ordered phase coexistence sheet: crossing it reverses the sign of $M$. Two tricritical wings also extend to nonzero $h$, separating weakly and strongly magnetized states of the same sign. Their ordinary critical points satisfy $f^{\prime }=f^{\prime \prime }=f^{\prime \prime \prime }=0$. Since $f^{\prime \prime \prime }=6uM+20v{M}^3$, their tricritical wing critical edges areAt these edges $f^{\prime \prime \prime \prime }=-12u>0$, giving an ordinary quartic critical expansion around $M_c$. All these structures end at the tricritical point $r=u=h=0$. The diagrams show both the zero-field section and the genuinely three-dimensional coexistence structure.

For the final critical exponent calculation, take $r=at$ with $a>0$, $t=(T-T_c)/T_c$, and hold the other controls fixed. At an ordinary critical point, the quartic term dominates the sextic term. For $h=0$ and $r<0$,whereas on the disordered side the singular contribution is zero. Two temperature derivatives give a finite jump in heat capacity, so $\alpha =0$. At $r=0$, the equation of state is $h=u{M}^3+O(M^5)$, so $\delta =3$.

Along the tricritical trajectory $u=0$, the nonzero minimum obeys $M^4=-r/v$ andThus two temperature derivatives produce $|t{|}^{-1/2}$. At $r=u=0$, $h=v{M}^5$. Therefore the mean-field critical exponents areWith $A$ the thermodynamic free energy, the physical sign is $C=-T{\partial }_T^{2}A$, rather than the positive sign printed in the paper. This correction affects the sign of the heat capacity, not its critical exponent. The tricritical result assumes that the independent quartic control is tuned to zero; a generic trajectory through nearby ordinary critical points need not have these exponents.

A tricritical point occurs when the continuous and first-order transition lines meet. In the symmetric sextic Landau free energy this requires tuning two coefficients, $r=u=0$, while $v>0$ stabilizes the potential. For $u>0$, the zero-field transition lies on $r=0$ and is continuous. For $u<0$, the coexistence curve is $r=3u^2/(16v)$. The order-parameter jump satisfies $M^2=-3u/(4v)$ and shrinks to zero as the tricritical point is approached. The labelled phase diagram below uses $(u,r)$ as local coordinates; two independent physical controls, such as temperature and composition, can provide these coordinates.

Along the tricritical tuning $u=0$, the equation of state is $h=rM+v{M}^5$. At zero field the ordered solution is $M=(-r/v{)}^{1/4}$, and the critical isotherm has $M=\mathrm{sgn}(h)(|h|/v{)}^{1/5}$. The inverse magnetic susceptibility is $r$ above and $4|r|$ below. The minimized singular potential is $-|r{|}^{3/2}/(3\sqrt{v})$ below, so differentiating twice thermally gives a $|t{|}^{-1/2}$ singularity. These calculations giveThe correlation-length and correlation-function exponents follow from the same quadratic gradient term. A generic path with a nonzero positive quartic coefficient instead has ordinary critical behaviour; observing tricritical powers requires the extra tuning. The sextic coupling is marginal in three dimensions, giving tricritical upper critical dimension three and logarithmic corrections there.

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.

At a tricritical point, both $r$ and $u$ vanish, while $v>0$ stabilizes the Landau free energy. It separates ordinary continuous transitions for $u>0$ from first-order coexistence for $u<0$. The zero-field order-parameter jump shrinks to zero as $u\uparrow 0$, and the tricritical wings terminate there in the full $(r,u,h)$ diagram. For $r<0$, $u=0$, minimizing gives $M^4=-r/v$ andTwo thermal derivatives give the heat-capacity critical exponent $\alpha =1/2$, compared with the finite mean-field heat-capacity jump, $\alpha =0$, at an ordinary quartic transition. The remaining tricritical mean-field critical exponents are $(\beta ,\gamma ,\delta )=(1/4,1,5)$, and the mean-field correlation-length critical exponent is $\nu =1/2$. To observe tricritical scaling one must tune the quartic scaling field as well as the thermal field: a generic path with fixed $u>0$ instead has ordinary critical behavior. Its upper critical dimension is three, as the multicritical Ginzburg ratio shows.

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 critical exponent is the leading power of a singular observable as a control variable approaches a critical value. Use $r\propto t=(T-T_c)/T_c$ and distinguish the exponents $\alpha ,\beta ,\gamma ,\delta ,\nu ,\eta$ from coefficient names in a Hamiltonian. Their standard definitions areThe Landau approximation predicts these powers by minimizing the Landau free energy and differentiating the equation of state. At an ordinary quartic point, $M^2=-r/u$, ${\chi }^{-1}=r$ above and $2|r|$ below, and $h=u{M}^3$ on the critical isotherm. The equilibrium singular free energy is proportional to $-r^2$ below, giving a heat-capacity jump. The quadratic fluctuation denominator is $\kappa k^2+V^{\prime \prime }$, so $\nu =1/2$ and $\eta =0$. Thus the ordinary mean-field critical exponents areAt a tuned tricritical point, the preceding solution gives $M\propto (-r{)}^{1/4}$, $V_{\min }\propto -|r{|}^{3/2}$ and ${\chi }^{-1}=4|r|$ below. At $r=0$ the equation of state is $h=v{M}^5$. The curvature still vanishes linearly, and the Gaussian gradient term still gives zero anomalous dimension. ThereforeThese mean-field tricritical exponent calculations require approach along the tuned quartic-zero direction. Fluctuations modify the powers below the appropriate upper critical dimension, and marginal effects can add logarithms at it.

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.

In a smooth scalar Landau free energy, a line of coexistence between two distinct minima can end by those minima and the intervening maximum coalescing. At a generic endpoint with no competing third phase,The first condition is stationarity, the second removes the inverse susceptibility, and the third is needed because nonzero cubic curvature at a zero-quadratic point would prevent it being a stable minimum. After shifting $M=M_c+\psi$ and mixing two control parameters, the local potential is $V_c+R{\psi }^2/2+U{\psi }^4/4-H\psi +\cdots$, $U>0$. Its two-phase coexistence line $H=0,R<0$ reaches $R=H=0$, with $\Delta M=2\sqrt{-R/U}\to 0$ and diverging magnetic susceptibility. This is the familiar critical endpoint of a quartic Landau free energy.

The word “must” needs the two-phase endpoint assumptions. A first-order line can instead end at a triple point or a boundary of the available control parameters. There is an explicit scalar counterexample in the sextic potential: fix $u<0<v$ and vary $(r,h)$. The $+M/-M$ coexistence line at $h=0$, $r<3u^2/(16v)$ ends when the $M=0$ phase also becomes globally stable. At that triple point the nonzero minima remain $\pm \sqrt{-3u/(4v)}$; their jump does not vanish. Thus the literal unrestricted assertion is false. The continuous-endpoint argument above establishes the intended case and distinguishes it from this counterexample.