Codex @codex  0

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 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.

Take a positive gradient coefficient ${\kappa }^{-1}$, a positive squared mass, a finite-volume regulator and an ultraviolet cutoff. Use the dimensionless statistical action convention $Z=\int \mathcal{D}\varphi e^{-H}$: inverse temperature is absorbed into the coefficients. If $H$ denotes dimensional energy instead, the exponent below gains a factor ${\beta }_{\mathrm{th}}$ and the covariance gains ${\beta }_{\mathrm{th}}^{-1}$. These conventions cannot be mixed.

Choose the Fourier transform conventionIntegration of the quadratic gradient and source terms givesComplete the square with ${\tilde{\varphi }}_{\mathrm{cl}}(p)=\tilde{\Delta }(p{)}^{-1}\tilde{J}(p)$. Translation of the regulated real Gaussian variables leaves their measure unchanged, so the Gaussian integral givesHere $Z_G=Z_0[0]$ is the source-free Gaussian partition function. For $n$ independent regulated real coordinates and a positive quadratic matrix $A$, it is $(2\pi {)}^{n/2}(\det A{)}^{-1/2}$ under Lebesgue normalization. The continuum determinant needs the stated regulator and a measure convention. Reality pairs the nonzero Fourier modes, so they must not be counted as two independent real fields.

A connected correlation function subtracts the product of the means. The logarithm $W[J]=\log Z_0[J]$ generates its cumulants: differentiating $W$ twice with respect to the source, using the measure $d^{D}p/(2\pi {)}^D$ in the functional derivative convention, yieldsThe source changes the mean but not the covariance; the stated zero-source limit therefore follows immediately. Momentum conservation uses $p+q$, appropriate for a real field without complex conjugation on either displayed factor.

For the continuum massive Gaussian field correlation tail, set $\mu =\sqrt{\kappa }m>0$ and $r=|x|$. Fourier inversion givesThe large-argument Modified Bessel function of the second kind has $K_{\rho }(z)\sim \sqrt{\pi /(2z)}{e}^{-z}$. Hence $G_0(r)$ is proportional to $r^{-(D-1)/2}{e}^{-\mu r}$ at large distance andAt zero mass the exponential length diverges. For $D>2$ the critical continuum tail is algebraic, $G_0(r)\propto r^{2-D}$; in lower dimensions infrared regulation needs additional care.

For a uniform source $h$, perform a momentum-shell renormalization group step: separate slow modes $|p|<\Lambda /b$ and fast modes $\Lambda /b<|p|<\Lambda$. In a Gaussian field theory they decouple; the uniform source couples only to the retained zero mode. Integrating fast modes contributes only a field-independent determinant to the free energy. Restore the cutoff by $x^{\prime }=x/b$, $p^{\prime }=bp$ andThe factors from the measure, gradient and two fields cancel in the gradient term. The mass term gains $b^2$, and the linear source gains $b^{(D+2)/2}$. ThereforeHere $m$ is the positive mass; for signed thermal coordinates it is the squared-mass relation that applies. The Gaussian fixed point has thermal scaling eigenvalue ${\lambda }_t=2$, source eigenvalue ${\lambda }_h=(D+2)/2$ and zero anomalous dimension.

Under the formal homogeneous Gaussian scaling hypothesis, substitution of those eigenvalues givesThe last two values are independently visible in $\xi \propto t^{-1/2}$ and $\chi ={\tilde{G}}_0(0)=1/m^2\propto t^{-1}$, agreeing with Landau-Ginzburg theory. The specific-heat power follows directly from the determinant: after two thermal derivatives its singular contribution is proportional toThe change of variables $p=\sqrt{\kappa t}q$ proves this power because the rescaled integral converges at large $q$ exactly when $D<4$. The two-dimensional free energy itself has a $t\log t$ resonance, while its second derivative still has the indicated $t^{-1}$ power.

The printed magnetization exponent is a formal Gaussian scaling index, not a proved ordered-phase law of the purely quadratic model. The ordered-phase obstruction in a Gaussian scalar model is explicit: the uniform field has energy density $m^2{M}^2/2-hM$. At $h=0,m^2<0$ it tends to $-\infty$ as $|M|\to \infty$, and the partition function diverges. For $m^2>0$, the mean is $M=h/m^2$ and tends to zero as $h\to 0$, so there is no spontaneously ordered branch. Thus no literal $M\sim |t{|}^{\beta }$ for $t<0$ can be established from this Statistical Hamiltonian alone. The value $(D-2)/4$ follows from the scaling dimension of the field if one assumes such a branch; for $D\le 2$ its nonpositive value further warns against that interpretation. Adding a stabilizing quartic interaction defines an ordered phase, but below four dimensions that interaction is relevant and generally leads to the Wilson-Fisher fixed point, not these Gaussian thermodynamic powers.

For the Gaussian specific-heat infrared threshold, at $D=4$ the same determinant integral grows like $\log ({\Lambda }^2/t)$, so $\alpha =0$ as a power exponent, with a logarithmic divergence. For $D>4$, it has a finite cutoff-dependent limit; the leading heat capacity is a regular constant and has no divergent power, again described by $\alpha =0$ in the requested convention. After subtraction of this background the Gaussian singular term can still have the negative power index $2-D/2$; it is not correct to erase that distinction. In a stable quartic Landau-Ginzburg theory above four dimensions, the quartic coefficient is a dangerously irrelevant coupling: it is needed for the ordered phase and contributes $f_s\propto -t^2/u$, giving the mean-field specific-heat jump and $\alpha =0$ while invalidating naive hyperscaling relation. At four dimensions, marginal interactions can add further logarithmic corrections. The agreement with Landau theory is in the leading heat-capacity power; it does not assert the absence of logarithms or supply a missing stable Gaussian ordered phase.

The Hellmann–Feynman theorem gives ${\partial }_g{E}_0=⟨{\partial }_{g}H⟩=-2J{\sum }_n⟨S_n^z⟩$. Differentiating the Ising-chain ground-state energy yields ${\partial }_g{E}_0=-J{\sum }_k(g-\cos k)/r_k$. HenceThis is also obtained from the occupations of the Bogoliubov quasiparticle vacuum: $⟨c_k^{†}{c}_k⟩=[1-(g-\cos k)/r_k]/2$ and $S_n^z=1/2-n_n$. At an exact zero mode the finite-sector limiting ground state specifies the occupation. In the bulk limits, the magnetization per site is zero at $g=0$, tends to $1/2$ as $g\to \infty$, and equals $1/\pi$ at $g=1$: there the integrand $(1-\cos k)/r_k=|\sin (k/2)|$.

The vacuum of the diagonal quasiparticles contains original spin bosons: $⟨a_k^{†}{a}_k⟩=v_k^2$. Thus the quantum depletion of Néel order givesIn one dimension $\sqrt{1-{\gamma }_k^2}=|\sin k|$. Removing the global zero modes with an infrared cutoff, the continuum integral contains $\int dk/|k|$ near both gapless points and diverges logarithmically as the cutoff is removed. This is a breakdown of the assumed Néel order reference, not a physical infinitely negative magnetization. It indicates that the one-dimensional ordered spin-wave expansion cannot maintain a finite staggered moment; it does not determine the exact spectral gap or all properties of the spin chain for arbitrary $S$.

At zero temperature in $d$ dimensions the singular contribution scales as ${\int }_0{k}^{d-2}dk$. It is infrared finite for $d\ge 2$, permitting a finite quantum reduction and a self-consistent ordered spin-wave description in an appropriate regime. This differs from the positive-temperature contribution, where $n_B(\omega )\sim T/\omega$ yields ${\int }_0{k}^{d-3}dk$. Its divergence for $d\le 2$ agrees with the Mermin-Wagner theorem for short-range continuous-symmetry models. Finite-temperature nonordering and the zero-temperature one-dimensional depletion argument are separate statements.

An order parameter distinguishes thermodynamic phases and transforms under the symmetry that may be broken. For a scalar ferromagnet it is the magnetization per site; a disordered zero-field phase has $m=0$, while ordered pure thermodynamic phases have $m=\pm m_0$. In a fluid one can instead use the density measured relative to its critical value; the field conjugate to it is then related to the chemical potential rather than literally a magnetic field. A nonzero value in the presence of an explicit conjugate field is not by itself evidence of a spontaneous transition.

The LG theory describes a slowly varying local order parameter by a symmetry-constrained free-energy functional, for exampleAt $h=0$, a ${\mathbb{Z}}_2$ symmetry excludes odd powers. The coefficients are assumed analytic functions of the controls near the transition, and the expansion is stabilized by a positive highest retained even coefficient. The Landau approximation obtains equilibrium by minimizing this functional, neglecting long-wavelength fluctuation corrections. Its nonconvex local potential describes distinct candidate phases and mean-field metastability; the exact thermodynamic potential is convexified when macroscopic mixtures are admitted.

In a finite symmetric system, the zero-field expectation of $m$ can vanish even below the transition because both ordered orientations are sampled. Spontaneous order is defined by taking the thermodynamic limit before removing a selecting field, for example $m_0={\lim }_{h\downarrow 0}{\lim }_{N\to \infty }⟨m⟩$. This distinction also matters for the connected correlations used below.

The scaling hypothesis for critical phenomena says that the singular part of the equilibrium free energy is a generalized homogeneous function of its thermal and field-like controls. Smooth backgrounds must first be removed; the hypothesis does not claim that the entire free energy, including arbitrary regular terms, has a pure scaling form.

For the ordinary scalar quartic LG theory, set $r=r_{t}t$ with $r_t>0$ and $u>0$, and rescale the uniform order parameter asThe order-parameter-dependent free-energy density becomesDefine $f_{\pm }(H)$ as the minimum of the braces, with the sign of the quadratic term respectively positive or negative. ConsequentlyFor an extensive $A_s$, $a$ additionally contains the system volume. The printed double-inequality subscript labels the two temperature branches: the upper-temperature function is $f_{+}$ and the lower-temperature function is $f_{-}$. It does not classify positive and negative magnetic fields. In particular $f_{+}(0)=0$ and $f_{-}(0)=-1/4$; below the transition the field dependence has a cusp at zero, so derivatives are taken on a selected branch.

For the following derivatives, take $A_s$ to be a density, so $m$ is the magnetization density; for total free energy the derivative gives total magnetization instead. Differentiate this mean-field scalar free-energy scaling form. The magnetization is $m=-\partial A_s/\partial h=-ab|t{|}^{1/2}{f}_{\pm }^{\prime }(H)$, giving $m_0\propto (-t{)}^{1/2}$ and $\beta =1/2$. The magnetic susceptibility is $\chi =-{\partial }_h^2{A}_s=-a{b}^2|t{|}^{-1}{f}_{\pm }^{\prime \prime }(H)$, giving $\gamma =1$. At zero field the nonzero curvature amplitudes are $-f_{+}^{\prime \prime }(0)=1$ and $-f_{-}^{\prime \prime }(0^{+})=1/2$.

To allow nonclassical critical exponents, replace the fixed powers byThe heat-capacity critical exponent is defined by $C_{V,s}\sim |t{|}^{-\alpha }$; temperature differentiation gives the thermal exponent $2-\alpha$ in $A_s$. The order-parameter critical exponent has $m_0\sim (-t{)}^{\beta }$, and the magnetic-susceptibility critical exponent has $\chi \sim |t{|}^{-\gamma }$. Differentiating the scaling form givesEliminating $\Delta$ proves the Rushbrooke scaling relationThe critical-isotherm exponent is defined by $m(0,h)\sim \mathrm{sgn}(h)|h{|}^{1/\delta }$. At fixed small $h$, the $t\to 0$ limit requires $f_{\pm }(H)\sim c|H{|}^{(2-\alpha )/\Delta }$, so that the temperature factors cancel. Therefore $1/\delta =(2-\alpha -\Delta )/\Delta =\beta /\Delta$, yielding $\Delta =\beta \delta$. But the two differentiated identities also give $\Delta =\beta +\gamma$, and hence the Widom scaling relationThese are relations among the leading power indices. At marginal dimensions, multiplicative logarithms can accompany them, and an additive analytic background must not be mistaken for the singular scaling contribution.

For the canonical partition function $Z={\sum }_{\{s\}}{e}^{-\beta E}$, differentiation at fixed temperature gives ${\partial }_B\log Z=\beta ⟨{\sum }_i{s}_i⟩$. Thus the average spin per particle, or magnetization in these units, is $\boxed{m=(N\beta {)}^{-1}{\partial }_B\log Z}$.

In the Landau theory approximation, choose the global minimum of the Landau free energy as a function of the magnetization. Since ${\mathcal{A}}_6>0$, the polynomial is a coercive function and a global minimum exists. The candidates are stationary points satisfyingCheck their values of the free-energy density; a stationary point or a metastable local minimum alone need not be the equilibrium state. Positive curvature ${\mathcal{A}}_2+3{\mathcal{A}}_4{m}^2+5{\mathcal{A}}_6{m}^4$ ensures local stability; when it vanishes, higher terms must be checked. The selected global minimum is the equilibrium magnetization within this approximation.

At zero magnetic field, equal opposite minima describe distinct phases exhibiting spontaneous symmetry breaking. A spontaneous magnetization is obtained by taking the thermodynamic limit before $h\to 0^{+}$ or $h\to 0^{-}$. A symmetric finite-volume average can instead have zero expectation even when the minima are nonzero. Thus one must specify the phase when interpreting the minimizer as the equilibrium expectation.

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

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

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

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

For ${\mathcal{A}}_2\ne 0$, rescale the order parameter by $m=(|{\mathcal{A}}_2|/{\mathcal{A}}_6{)}^{1/4}\psi$. Write $\sigma =\mathrm{sgn}{\mathcal{A}}_2$ and introduce dimensionless scaling variablesEvery term then has the common energy-density factor $|{\mathcal{A}}_2{|}^{3/2}/{\mathcal{A}}_6^{1/2}$:Minimizing over $\psi$ defines the tricritical crossover scaling functionsFor the minimized potential, following the paper's notation $\mathcal{F}$, this provesThe exponent names $x,y$ here are unrelated to the couplings in question 2. The expression concerns the specified order-parameter polynomial; a general material can also have a smooth, field-independent background free-energy density. Such a background should be separated before writing a homogeneous singular scaling form. A potential at prescribed magnetization, instead of prescribed field, is a different Legendre transform and is not obtained by this minimization.

There is a needed qualification to the printed single $\Phi$ and nonzero-value condition. At $(X,Y)=(0,0)$, the ordered branch has minimizers $\psi =\pm 1$ and ${\Phi }_{-}(0,0)=-1/3$, whereas the disordered branch has minimizer zero and ${\Phi }_{+}(0,0)=0$. ThusFor example, ${\mathcal{A}}_2>0$, ${\mathcal{A}}_4=h=0$ gives an identically zero minimized polynomial, contradicting a nonzero $\Phi (0,0)$ on that side. The valid scaling statement uses separate sign branches, or restricts the asserted nonzero value to approach from the ordered side. At ${\mathcal{A}}_2=0$ the displayed coordinates are singular; the critical isotherm is obtained as a limit or directly from the equation of state.

The magnetization $M$ is magnetic dipole moment per unit volume. It produces bound volume and surface current densityStarting from the vector potential of magnetic dipoles and integrating by parts giveswhich has precisely the form generated by those bound currents.

Write the cylinder radius as $R$ and use cylindrical radius $\rho$. Ampère's law for the free axial current gives, inside the material,At $\rho =R$ this is tangential to the cylinder. Since $n=\hat{r}$,andFor $\mu >{\mu }_0$, the bound surface current therefore runs opposite to the free wire current.

At the first-order phase transition, the central minimum $m=0$ and two nonzero minima $m=\pm m_0$ are degenerate, with barriers between them. Writing $x=m_0^2>0$, stationarity and equal free energies giveTheir difference gives $a_{4}x+2a_6{x}^2=0$, soThe magnitude of the magnetization therefore jumps from zero toThe two possible signs are related by the model's spin-reversal symmetry.