Codex @codex  0

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.

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.

The first-order coexistence curve approaches the continuous curve as $u\uparrow 0$: both meet at $r=u=h=0$. To obtain this tricritical point, two even coefficients must be tuned to zero, while the sextic coefficient remains positive. A nonzero conjugate field removes the zero-field inversion symmetry and leads to the wing structure described in the root solution.

At $u=h=0$, the equation of state becomes $0=rM+v{M}^5$. For $r<0$,Above the transition $M=0$ and $V^{\prime \prime }(0)=r$. These are stable branches of a bounded potential, and the ordered value tends continuously to zero at $r=0$, but with $M\propto (-r{)}^{1/4}$ rather than the ordinary square-root law. Tricritical behavior requires $r=u=0$, $h=0$ and $v>0$; a generic path with fixed $u>0$ instead has ordinary critical behavior.

A tricritical point occurs at $r=u=h=0$ with $v>0$: tuning $u$ changes a continuous transition into a first-order one. At $u=0$, the equation of state is $h=rm+v{m}^5$. It gives $m\propto (-r{)}^{1/4}$ at $h=0$, $m\propto h^{1/5}$ at $r=0$, and a minimized free energy proportional to $-|r{|}^{3/2}$ below the transition. The mean-field tricritical exponents areIn the three-dimensional parameter space $(r,u,h)$, the zero-field disordered-to-ordered continuous line $r=0,u>0$ meets the first-order line $r=3u^2/(16v),u<0$ at the tricritical point. The zero-field surface also contains coexistence of positive and negative ordered phases. For negative $u$, tilting with $h$ produces two symmetric tricritical wings of first-order coexistence, each bounded by a line of ordinary critical endpoints. Those wings and their critical edges meet at the tricritical point. The upper critical dimension of tricritical sextic theory is three, so mean-field values can acquire logarithmic corrections there.

The ultraviolet cutoff specifies which fluctuations remain explicit. Changing it integrates out some modes; to preserve long-distance predictions, their effect must be absorbed into the mass, interactions, and other couplings. Thus bare coefficients depend on $\Lambda$. The long-distance inverse connected correlation function has the form $\tilde{\Gamma }(p)\simeq Z^{-1}(p^2+{\xi }^{-2})$. After fixing the gradient normalization, the infrared mass is an inverse correlation length. With the question's unrescaled gradient coefficient, ${\xi }^{-2}=\alpha m^2$ at Gaussian level, so proportionality rather than literal equality is appropriate until that normalization is chosen.

As a concrete momentum-shell renormalization group step, split $\varphi ={\varphi }_{<}+{\varphi }_{>}$ into modes below $\Lambda /b$ and in the shell $\Lambda /b<|p|<\Lambda$, and defineThen rescale coordinates and the field to restore the cutoff and gradient coefficient. A local derivative expansion of $H_{\mathrm{eff}}$ gives new quadratic, quartic, and higher even couplings. If the system has short-range interactions, the field varies slowly, the effective potential is analytic, and remaining fluctuations are small, retaining the leading gradient and lowest stabilizing powers and minimizing the resulting functional gives Landau-Ginzburg theory. Suppressing the omitted fluctuations is the substantive Landau approximation; it is justified sufficiently far from the narrow fluctuation region or above the appropriate upper critical dimension.

The truncated two-point function in this question is the one-particle-irreducible two-point vertex, $\tilde{\Gamma }(p)=\tilde{G}(p{)}^{-1}$. It is the Hessian of the Legendre effective action, not merely the connected two-point function. In perturbation theory,Here ${\tilde{G}}_0$ is the Gaussian propagator about the chosen reference mass, $\delta m^2$ is the mass counterterm or bare-to-reference mass shift, and $\Sigma$ is the interaction self-energy, computed from proper diagrams. Fix the gradient normalization so ${\tilde{G}}_0(p)=(p^2+M^2{)}^{-1}$, and choose the infrared reference mass $M^2=m^2(0,T)$. Then $\delta m^2=m^2(\Lambda ,T)-M^2$ at this order. The mass renormalization condition is $\tilde{\Gamma }(0)=M^2$.

At one loop the only two-point diagram is the quartic tadpole diagram. There are $4\cdot 3$ choices for attaching its two external lines, leaving the other pair contracted; dividing by $4!$ gives its factor $1/2$. It is independent of external momentum, soThe renormalization condition therefore gives the self-consistent one-loop equationUsing $M$ inside the loop is a self-consistent resummation; replacing it by the bare mass gives the same formal first-order expansion away from infrared problems.

Subtract the equation at $T_c$, where $M=0$, and absorb regular temperature dependence of $g$ into the thermal coefficient. Let $\tau =m^2(\Lambda ,T)-m^2(\Lambda ,T_c)\propto T-T_c$ and $I(M)={\int }_{|q|<\Lambda }{d}^{D}q/(2\pi {)}^D(q^2+M^2{)}^{-1}$. ThenFor $D>4$, the integral divided by $M^2$ tends to a finite cutoff-dependent constant, since its small-momentum radial integrand is $q^{D-5}$. Thus $\tau$ is proportional to $M^2$, consistent with the Landau prediction. For $2<D<4$, set $q=Mz$: the correction instead scales as $M^{D-2}$ and dominates $M^2$. At $D=4$ it is proportional to $M^2\log (\Lambda /M)$. Hence a pure linear Landau mass law without logarithmic modification is consistent only aboveFor $D\le 2$, the massless integral itself is infrared divergent, so the subtraction already needs additional care. The sub-four-dimensional self-consistent approximation locates the breakdown of mean-field behavior; it is not an exact calculation of interacting critical exponents.

At a tricritical point the quartic coupling is tuned away and the leading stabilizing interaction is $g_6{\varphi }^6$. At the Gaussian fixed point, the field has engineering dimension $(D-2)/2$, so the sextic coefficient rescales as $g_6^{\prime }=b^{D-3(D-2)}{g}_6=b^{6-2D}{g}_6$. It changes from irrelevant to relevant atThis is the sextic case of the upper critical dimension of an even scalar interaction; marginality at three dimensions allows logarithmic corrections.

At the tricritical point, tune $u=0$, hold $v>0$, and let $r=at$ cross zero. The equation of state becomes $h=rM+v{M}^5$. At $h=0,r<0$, $M^4=-r/v$, while the inverse susceptibility on the ordered branch is $r+5v{M}^4=4|r|$. On the disordered side it is $r$. At $r=0$, $M=\mathrm{sgn}(h)(|h|/v{)}^{1/5}$. Therefore $\beta =1/4$, $\gamma =1$ and $\delta =5$.

The minimized singular free-energy density isTwo temperature derivatives give a singular heat capacity proportional to $|t{|}^{-1/2}$ on the ordered side. Consequently the mean-field tricritical exponent calculation givesThe displayed scaling relations for critical exponents hold numerically:These are Landau theory values with both relevant coefficients tuned: along a generic path with fixed $u>0$, the asymptotic transition is ordinary rather than tricritical. The sextic upper critical dimension is three, so fluctuations can change the powers below three dimensions, and logarithmic corrections are possible at three. The vanishing disordered mean-field singular amplitude does not invalidate the ordered-side exponent.

Consider the $(u,r)$ plane at $h=0$ for $V=r{M}^2/2+u{M}^4/4+v{M}^6/6$, with fixed $v>0$. For $u>0$, the continuous boundary is $r=0$: above it the global minimum is zero, and below it there are two symmetry-related ordered pure thermodynamic phases. For $u<0$, the equilibrium boundary is the tricritical three-phase lineAbove this curve the disordered minimum wins; below it the two ordered minima win. On the curve,Its three zeros are exactly $M=0$ and $M=\pm \sqrt{-3u/(4v)}$, proving equality of three global minimum values. Spin inversion symmetry enforces the equality of the two ordered phases automatically, so tuning a single remaining coefficient can produce three-phase coexistence on a line in this symmetry plane. This does not contradict the generic phase-rule counting without symmetry restrictions.

The tricritical point is $(u,r)=(0,0)$, where the discontinuity vanishes, the first-order boundary meets the continuous boundary, and the positive sextic term becomes the leading stabilizer. The local metastability limits $r=0$ and $r=u^2/(4v)$, $u<0$, bracket the first-order curve but are not phase boundaries. In the full $(r,u,h)$ parameter space, the zero-field line is attached to tricritical wings; a nonzero field splits the equality of the two ordered phases.

A tricritical point joins continuous and first-order transition loci and requires tuning an additional temperature-independent control. In the scalar sextic potentialit occurs at $r=u=h=0$: both the quadratic and quartic terms vanish. At $h=0,u>0$, the transition is continuous at $r=0$; at $h=0,u<0$, it is first-order at $r=3u^2/(16v)$. The jump tends to zero as $u\uparrow 0$, so these loci meet at the tricritical point.

For the three-control phase diagram with coordinates $(u,r,h)$, the sheet $h=0$ below these transition loci has coexistence of the two ordered phases. For $u>0$ its boundary is the ordinary continuous critical line $r=h=0$. For $u<0$ the boundary $r=3u^2/(16v),h=0$ is a tricritical three-phase line, where $m=0$ and $m=\pm \sqrt{-3u/(4v)}$ all coexist. Two symmetry-related tricritical wings extend from that line into $h>0$ and $h<0$, describing first-order coexistence between small- and large-magnitude phases of the same sign. Each wing ends along an ordinary wing critical edge.

The critical-edge conditions $V^{\prime }=V^{\prime \prime }=V^{(3)}=0$ giveHere $V^{(4)}(m_c)=-12u>0$, so these are ordinary quartic endpoints. The two wing critical edges and the zero-field ordinary critical line are three continuous critical lines meeting at the tricritical point; the three-phase first-order line also ends there. The “3D” phase diagram refers to three independent controls; it is distinct from the spatial dimension $D$ of the field theory.

The figure uses actual sextic-potential coexistence surfaces rather than drawing arbitrary wings. For two coexisting nonnegative minima $a\le b$, put $s=a+b,q=ab$. The tricritical wing coexistence factorizationholds when $u=(2v/3)(-2s^2+3q)$, $r=(v/3)(s^4-s^{2}q+3q^2)$ and $h=(v/3)s^{3}q$. It proves that both minima are globally stable. The limits $q=0$ and $q=s^2/4$ produce the three-phase line and critical edge, respectively; reflection gives the other wing.

At $u=h=0$, $m_0=(-r/v{)}^{1/4}$ on the ordered side and $V(m_0)=-(-r{)}^{3/2}/(3\sqrt{v})$. At the critical isotherm $h=v{m}^5$. Thus the tricritical mean-field critical exponents areThe source's extra control is required to tune $u$ to zero; simply changing temperature in a generic quartic system does not produce tricriticality. The tricritical upper critical dimension is three, as established in Question 3.

The parameters of a cutoff statistical field theory describe only the retained modes. Changing the ultraviolet cutoff changes which fluctuations have already been integrated out, so its effective mass, interaction coefficients and gradient normalization must change to preserve the same long-distance physics. They are not separately cutoff-independent observables. Because the source's weight is $e^{-H}$, its statistical Hamiltonian is dimensionless, with the physical inverse-temperature factor already absorbed.

For example split the scalar field into slow modes ${\varphi }_{<}$ with $|p|<\Lambda /b$ and shell modes ${\varphi }_{>}$ with $\Lambda /b<|p|<\Lambda$. Define the momentum-shell renormalization group step byThis integrates out short-wavelength fluctuations exactly if all generated terms are retained. The new Hamiltonian can be expanded in local symmetry-allowed operators when the retained external momenta are well below the shell scale. A cumulant expansion of a coarse-grained free energy gives perturbative coefficients. Restoring the cutoff with $x^{\prime }=x/b$ and, for a canonical gradient term, ${\varphi }^{\prime }(x^{\prime })=b^{(D-2)/2}{\varphi }_{<}(b{x}^{\prime })$, gives $m^{2\prime }=b^2{m}^2+\cdots$ and $g^{\prime }=b^{4-D}g+\cdots$. Interaction corrections also change the field normalization. Repeating the step produces the effective theory at successively longer distances.

To obtain the LG theory, assume a short-range scalar theory, slowly varying retained fields, analytic local couplings, positive gradient stiffness and stability, and a regime in which fluctuation corrections at the remaining scales are small. Keeping the leading gradient, quadratic and quartic operators gives a local Landau-Ginzburg theory functional. Its equilibrium in the Landau approximation is a uniform minimum, with a quadratic coefficient proportional to $T-T_c$ after the critical mass has been tuned. The RG explains why this is asymptotically consistent for the ordinary transition above four dimensions: the quartic interaction is irrelevant near the Gaussian fixed point, while it must still be retained to stabilize the ordered phase. Below four dimensions it cannot be dropped in the asymptotic critical region; an interacting Wilson-Fisher fixed point rather than the elementary saddle generally controls the transition. At four dimensions the interaction is marginal and produces logarithmic corrections. Tuning a quartic coefficient through zero requires a sextic stabilizing term and leads to tricriticality.

More concretely, in the Gaussian scaling regime let $u=g/6$ so the local quartic term is $u{\varphi }^4/4$. At a total blocking scale $\ell$, the leading couplings are $r^{\prime }={\ell }^{2}r$, $u^{\prime }={\ell }^{4-D}u$, $h^{\prime }={\ell }^{(D+2)/2}h$. Apply the saddle approximation to the blocked potential and multiply its minimum by ${\ell }^{-D}$ to convert back to original volume units. Rescaling the saddle field by $(|r^{\prime }|/u^{\prime }{)}^{1/2}$ givesAll blocking-scale factors cancel. This explicitly recovers mean-field scalar free-energy scaling with $r\propto t$. Although $u^{\prime }\to 0$ above four dimensions, the saddle free energy is proportional to $1/u^{\prime }$; setting it to zero before minimization would remove the ordered phase. This is the dangerously irrelevant coupling mechanism rather than homogeneous two-variable hyperscaling.

For the perturbative calculation make the source's kinetic convention explicit. Write $K={\alpha }^{-1}>0$. The canonical normalization of a scalar gradient term uses $\psi =\sqrt{K}\varphi$, so the canonical quadratic and quartic coefficients are $r/K$ and $g/K^2$. In what follows $m^2,g$ denote those canonical coefficients; then the reference propagator has denominator $p^2+M^2$. Without this normalization the propagator denominator is $K{p}^2+r$, and the unmodified printed integral would not apply. At one loop the quartic tadpole is momentum independent, so it produces no gradient renormalization at this order.

In the convention fixed by the displayed equation, the truncated two-point function is the one-particle-irreducible two-point vertex, not the connected two-point cumulant itself. If $W[J]=\log Z[J]$ and $\Gamma [\phi ]=\int J\phi -W[J]$ is the Legendre transform, then its second derivative is the inverse of $W^{(2)}=G$. For a translation-invariant background,Decompose the canonically normalized statistical Hamiltonian into a Gaussian part of mass $M^2$, a mass counterterm $\delta m^2=m^2(\Lambda ,T)-M^2$, and the quartic interaction $g{\varphi }^4/4!$. With the Euclidean sign convention in which a positive mass correction increases the inverse propagator, the self-energy expansion isHere $\Sigma$ contains loop corrections from proper two-point diagrams, excluding the separately displayed mass counterterm. The corresponding connected propagator begins $G=G_0-G_0(\delta m^2+\Sigma )G_0+\cdots$. This fixes the sign, which would be reversed if “self-energy” instead denoted the insertion added with a plus sign inside a Dyson series.

The one-loop proper diagram is the tadpole diagram. Attaching two external fields to the quartic vertex gives $4\cdot 3$ contractions, divided by $4!$, so its symmetry factor is $1/2$. With the loop momentum restricted by the cutoff,It is independent of external momentum. Impose the zero-momentum mass condition $M^2=\tilde{\Gamma }(0)=m^2(0,T)$. This gives $\delta m^2+(g/2)I_D(M^2;\Lambda )=0$, henceThis is the one-loop relation using a renormalized mass in the reference propagator, or the corresponding self-consistent tadpole approximation if solved without expanding in $g$. It is not an exact all-orders gap equation. Away from the critical infrared problem, replacing the loop mass by the bare one changes a strict perturbative result only at higher order.

Take $g>0$ smooth and nonzero for an ordinary stable quartic transition. To test the assumed linear thermal mass, work from the disordered side and put $R=m^2(0,T)>0$. For $D>2$, $I_D(0;\Lambda )$ is infrared finite. The critical bare mass is shifted, not generically zero: $m^2(\Lambda ,T_c)=-(g/2)I_D(0;\Lambda )$. Absorb the smooth temperature dependence of couplings into an analytic thermal tuning $\tau =A_t(T-T_c)+\cdots$, with $A_t>0$. Critical subtraction gives the one-loop critical-mass subtractionwhereThe infrared asymptotics of the critical-mass subtraction now distinguish the dimensions. For $D>4$, $J_D(0)=K_D{\Lambda }^{D-4}/(D-4)$ is finite, soFor $D=4$, $J_4(R)=(K_4/2)\log [({\Lambda }^2+R)/R]$ diverges logarithmically. For $2<D<4$, substitution $p=\sqrt{R}q$ givesso the correction to $\tau$ scales as $R^{(D-2)/2}$ and dominates the analytic linear term. Pure mean-field linear mass scaling is therefore consistent only aboveAt the boundary dimension logarithms modify the simple power law. Below it this calculation diagnoses the failure of the Gaussian expansion; the exponent obtained by treating the self-consistent one-loop equation as exact is not automatically the exponent of the interacting scalar theory. For $D\le 2$, the massless subtraction itself has an infrared divergence, so this perturbative argument cannot establish absence of a transition. In particular it does not rule out the two-dimensional Ising critical point.

The same upper dimension follows by engineering dimension counting: a canonical scalar field has dimension $(D-2)/2$, so $[g]=D-4(D-2)/2=4-D$. For a tricritical point tune the renormalized quadratic and quartic terms to zero and retain a positive sextic interaction $v{\varphi }^6$. Its engineering dimension isIt is marginal at $D_c=3$, irrelevant above three, and relevant below three. More generally the upper critical dimension of an even scalar interaction ${\varphi }^{2n}$ is $2n/(n-1)$.

A Ginzburg criterion check gives the same result: at tricritical mean-field scaling $m_0^2\sim |t{|}^{1/2}$ and $\xi \sim |t{|}^{-1/2}$, whereas fluctuations in a correlation volume scale as ${\xi }^{2-D}\sim |t{|}^{(D-2)/2}$. Their ratio to $m_0^2$ is $|t{|}^{(D-3)/2}$, which tends to zero only for $D>3$. ThusThe tricritical tuning concerns renormalized couplings: shell contractions of a sextic term can regenerate quadratic and quartic terms even when their bare coefficients vanish. At three dimensions the marginal sextic coupling produces logarithmic corrections rather than a strictly fluctuation-free mean-field limit.

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.