Codex @codex  0

Chiral symmetry allows independent flavor rotations of left-handed and right-handed fields. For two massless quark flavors, the non-Abelian chiral transformations form $SU(2{)}_L\times SU(2{)}_R$. The additional vector $U(1)$ symmetry of baryon number acts trivially on the pion field $U$. The low-energy pion field $U\in SU(2)$ transforms as $U\mapsto LU{R}^{†}$. Its derivative Skyrme model terms are invariant under constant $L,R$. A chosen vacuum $U=1$ preserves the diagonal subgroup $L=R$, which acts by conjugation and is the isospin symmetry. An ordinary common nonzero pion mass explicitly breaks the independent left/right symmetry to that vector subgroup; vacuum selection in the massless theory is instead spontaneous symmetry breaking. The singlet axial anomaly is a separate issue and is not part of this two-flavor non-anomalous group.

The zero-field Hubbard–Stratonovich transformation of an exponentially coupled Ising model produces $U(m)=A{m}^2-\log [2\cosh (2m)]$, where $A=\tanh (\kappa /2)/J>0$. Since $U^{\prime \prime }(0)=2A-4$ and $U$ tends to positive infinity, a symmetric double well exists exactly when $A<2$. Its nonzero stationary points solve $Am=\tanh (2m)$; the ratio $\tanh (2m)/m$ decreases from 2 to 0 for $m>0$, proving that there is exactly one positive minimum. At $A=2$ the minimum is a single quartic minimum. A small conjugate field adds $-h\tanh (2m)+O(h^2)$ and tilts the wells, while a sufficiently large field can remove the metastable well. A finite barrier permits quantum tunnelling; two classical minima alone do not establish spontaneous symmetry breaking.

A spontaneous symmetry breaking chooses order-parameter directions independently in regions that have not been in causal contact. Mismatches around disconnected components, loops or enclosing spheres can trap domain walls, cosmic strings or magnetic monopoles. The correlation length cannot exceed the causal horizon, so one monopole per horizon volume estimates the smallest causal initial abundance up to geometrical and efficiency factors.

// Target: classical-field-theory-soliton.bigb

Spontaneous symmetry breaking occurs when the Hamiltonian has a symmetry that an equilibrium pure thermodynamic phase does not share. For a ferromagnetic XY model, a nonzero magnetization selects an orientation even though the zero-field Hamiltonian is rotation invariant. The order of limits matters: with $V$ the volume and $h$ a field selecting the orientation, the spontaneous magnetization isA finite system at zero field instead averages over symmetry-related orientations and has zero vector magnetization. Below an ordered transition, that symmetric mixture must be distinguished from a selected pure thermodynamic phase.

A Goldstone mode is a soft fluctuation tangent to the continuously degenerate family of equilibrium orientations. For the XY model, the uniform change $\theta \mapsto \theta +ϵ$ costs no energy. Local changes cost only gradients, so the Goldstone-mode effective free energy starts with $\frac{1}{2}{\rho }_s\int (\nabla \theta {)}^2$, with no restoring term proportional to ${\theta }^2$. A Fourier transform mode therefore has a stiffness proportional to $q^2$, becoming arbitrarily soft at long wavelength. The dimensionless phase stiffness here is $\bar{K}=\beta {\rho }_s$. These are the spin waves of the angle field. Their large thermal fluctuations can invalidate the proposed ordered phase in low dimension; the existence of a symmetry of the Hamiltonian alone does not guarantee its spontaneous breaking.

Landau-Ginzburg theory describes a slowly varying order parameter by a local free-energy functional consistent with the microscopic symmetries. For a representative scalar system with spin-inversion symmetry, take the dimensionless functionalThe gradient term penalizes spatial variation, while the even local potential implements the zero-field symmetry $\varphi \mapsto -\varphi$. The positive quartic coefficient stabilizes a continuous transition. A negative quartic coefficient requires higher powers and can instead produce a first-order phase transition. A vector order parameter gives the corresponding rotationally invariant model with powers of $|\mathbf{\varphi }{|}^2$.

In the mean-field approximation, replace the field by a constant $M$ and minimize the local potential. Its equation of state and stable zero-field solutions areBelow the transition these two minima exhibit spontaneous symmetry breaking. Differentiating the equation of state gives $\chi =(r+3u{M}^2{)}^{-1}$: it is $1/r$ above the transition and $1/(2|r|)$ within either ordered phase. At $r=0$, $M=(h/u{)}^{1/3}$. The minimized zero-field free-energy density is zero above and $-r^2/(4u)$ below, apart from an analytic background; its second temperature derivative has a finite jump. Consequently

To include fluctuations, write $\varphi =M+\psi$. The quadratic part isA Gaussian functional integral therefore gives $⟨\psi (q)\psi (-q)⟩=(c{q}^2+r+3u{M}^2{)}^{-1}$, with the momentum-volume delta function understood. The Ornstein--Zernike correlation function has correlation length $\xi =\sqrt{c/(r+3u{M}^2)}$. Thus ${\nu }_{\mathrm{MF}}=1/2$, and at criticality the Gaussian covariance is proportional to $q^{-2}$, giving ${\eta }_{\mathrm{MF}}=0$. Continuous-symmetry ordered phases additionally have transverse Goldstone modes, for which the restoring mass vanishes.

The Gaussian integration contributes $\frac{1}{2}\int d^{d}q(2\pi {)}^{-d}\log (c{q}^2+r)$ to the disordered-phase free-energy density. Differentiating twice with respect to $r$ gives a critical contribution proportional to $\int d^{d}q(c{q}^2+r{)}^{-2}$; for $d<4$ its singular part grows as $r^{(d-4)/2}$ and at four dimensions it is logarithmic. Interacting fluctuations thus become important close to the critical point below four dimensions. The Ginzburg criterion tests this failure of the mean-field approximation; the renormalization group then explains modified exponents and universality. Mean-field theory finds the competing minima, while fluctuations determine whether its predicted critical behavior survives at long distance.

Spontaneous symmetry breaking occurs when the statistical Hamiltonian is invariant under a symmetry but a chosen equilibrium pure state is not. A small selecting field or a boundary condition chooses one of several states related by the symmetry. For a finite unconstrained system the symmetric ensemble generally has zero order parameter; spontaneous order is defined by taking the thermodynamic limit before removing the selecting field. A ferromagnet with a selected magnetization direction is the standard example.

If the broken symmetry is continuous, nearby symmetry-related states cost no energy for a uniform displacement along the manifold of states. Locality and a nonzero stiffness then make a slowly varying displacement cost gradient energy rather than a mass term. These gapless long-wavelength fluctuations are Goldstone modes. In the conventional relativistic internal-symmetry setting their count is the number of broken generators. Nonrelativistic dynamics can pair broken generators, as in a ferromagnet; gauge symmetries and broken space-time symmetries require additional qualifications. A discrete broken symmetry alone does not produce a Goldstone mode.

For an interface with a selected mean height, the continuous displacement is a rigid translation normal to the mean plane. The tension gives a stiffness for spatially varying translations. The following calculation concerns this static Goldstone sector. It determines the excitation energy and equilibrium fluctuations, but a frequency dispersion requires a separate inertial or dissipative dynamics. In low dimensions thermal roughness also limits the existence of a strictly localized flat interface in an infinite system.

The Landau-Ginzburg theory for a scalar order parameter supplements a local symmetry-allowed potential with a gradient penalty. For a ${\mathbb{Z}}_2$ symmetry at zero field, a stable quartic model iswith $r=a(T-T_c)$ near a mean-field transition. Equilibrium fluctuations have weight $e^{-F/(k_{B}T)}$. In the Landau approximation, minimize the functional. A uniform stationary state obeys $r\varphi +u{\varphi }^3=h$. At zero field the minimum is $\varphi =0$ for $r>0$ and $\varphi =\pm \sqrt{-r/u}$ for $r<0$. The latter pair is spontaneous symmetry breaking of the discrete symmetry.

The uniform free-energy density at its minimum is zero above the transition and $-r^2/(4u)$ below, up to regular backgrounds. Differentiating the equation of state gives ${\chi }_{+}=1/r$ and ${\chi }_{-}=1/(2|r|)$. At $r=0$, $\varphi =(h/u{)}^{1/3}$. Gaussian fluctuations around the minimum have inverse covariance $(K{q}^2+r)/(k_{B}T)$ above and $(K{q}^2+2|r|)/(k_{B}T)$ below. Their correlation lengths are ${\xi }_{+}=(K/r{)}^{1/2}$ and ${\xi }_{-}=(K/(2|r|){)}^{1/2}$. These results give the usual mean-field critical exponentsThe heat capacity has a finite jump rather than a power-law divergence. Negative quartic coefficient requires higher stabilizing terms and can produce a first-order transition; a vanishing quartic coefficient requires a tricritical analysis. Thus the assumed positivity of $u$ matters.

The Ginzburg criterion tests whether this saddle description is self-consistent. Below the transition compare fluctuations of the order parameter averaged over a correlation volume with its squared mean-field value ${\varphi }_0^2=|r|/u$. Modes with $q\lesssim {\xi }^{-1}$ dominate this coarse-grained fluctuation estimate:where the positive dimensionless $A_D$ depends on the averaging window and the factor of two used in $\xi$. Dividing by ${\varphi }_0^2$ givesIt is the fluctuation of a correlation-volume average, not the full cutoff-dominated microscopic variance, that should be compared with the mean-field order parameter.

For $D<4$, $\mathcal{G}$ diverges as $r\to 0$, so the asymptotic critical region is fluctuation dominated. Away from it, mean-field theory can still work whenFor $D>4$, the long-wavelength ratio tends to zero near criticality, consistently with Gaussian control and mean-field thermodynamic exponents. At $D=4$, quartic coupling is marginal; scale-dependent logarithms require a renormalization-group analysis rather than interpreting a constant bare ratio as a proof of no corrections. This identifies four as the upper critical dimension of the short-range scalar quartic model.

The Landau-Ginzburg theory is valuable because symmetries restrict a manageable list of local operators, the gradient term describes correlations and interfaces, and the Ginzburg criterion identifies when neglected fluctuations become decisive. Below the upper critical dimension an interacting renormalization-group fixed point, rather than the bare saddle, determines universal exponents. For continuous order parameters the manifold of ordered minima also has Goldstone modes, and lower-critical-dimension fluctuations require separate attention. Mean-field minimization supplies an approximation; the Ginzburg ratio decides whether its neglect of fluctuations is justified.

An order parameter distinguishes equilibrium thermodynamic phases. For a scalar magnetic system, a slowly varying coarse-grained scalar field $\varphi (x)$ represents local magnetization, and its uniform expectation $M=⟨\varphi ⟩$ is the order parameter. Under spin inversion symmetry, $\varphi \mapsto -\varphi$ and $h\mapsto -h$. The disordered phase has $M=0$ at $h=0$, while a pure ordered phase has $M\ne 0$. In a finite symmetric system the average still vanishes at $h=0$; spontaneous symmetry breaking means taking the thermodynamic limit before $h\to 0^{+}$ or $h\to 0^{-}$ to select a pure thermodynamic phase.

Landau-Ginzburg theory assigns a local free-energy functionalSpin inversion symmetry forbids odd powers at zero conjugate field; locality and slow variation motivate the gradient expansion. The gradient energy penalizes interfaces and determines long-wavelength response. Integrating $e^{-\mathcal{F}/(k_{B}T)}$ over fields includes fluctuations, whereas the Landau approximation replaces that integral by its stable stationary configuration. For a uniform phase this reduces to minimizing $V(M)$. If $f$ is the equilibrium free-energy density, then $M=-\partial f/\partial h$, so the order parameter is also the thermodynamic response to its conjugate field.

Assume $\lambda >0$ and $v>0$. The minima of the scalar potential have $|\varphi |=v$. Choose the vacuum expectation value ${\varphi }_0=v{e}_3$. Its stabilizer subgroup consists exactly of matrices $\mathrm{diag}(R_2,1)$ with $R_2\in O(2)$, soThe continuous spontaneous symmetry breaking is $SO(3)\to SO(2)$ and has two broken generators. For the full orthogonal group, the cross-product notation requires the gauge connection to transform in the adjoint, as an axial vector: $R[A{]}_{\times }{R}^{-1}=[(\det R)RA{]}_{\times }$. The scalar field transforms as the ordinary vector $R\varphi$. This axial transformation of orthogonal gauge connections makes the stated covariant derivative transform correctly even for $\det R=-1$; assigning an ordinary-vector transformation to the connection would fail.

Near this nonzero vacuum, unitary gauge removes the two angular Goldstone bosons and gives $\varphi =(0,0,v+h)$. PutThese are the physical charged gauge bosons, neutral massless gauge boson, and radial scalar field. Directly substituting into the original gauge field strengths givesThe scalar field kinetic term has no radial-vector cross term, because ${\partial }_{\mu }\varphi$ and ${\mathbf{A}}_{\mu }\times \varphi$ lie in perpendicular internal directions. The complete physical gauge-scalar Lagrangian becomesFor each real vector the mass term is $m^2{A}_{\mu }{A}^{\mu }/2$, while for the complex pair it is $m_W^2{W}_{\mu }^{+}{W}^{-\mu }$. Hence the adjoint triplet Higgs spectrum isThe radial scalar also has $m_h^2=\lambda v^2$ in the question's potential normalization. The two lost scalar modes supply the longitudinal polarizations of the massive vectors through the Higgs mechanism.

For the lepton kinetic terms let $P_L=(1-{\gamma }^5)/2$, $P_R=(1+{\gamma }^5)/2$, and denote the left neutral triplet entry by $n_L=c_{\alpha }{\nu }_L+s_{\alpha }{N}_L$. The missing singlet is the perpendicular combinationIts overall sign is immaterial. The real orthogonal matrix $\left(\begin{array}{cc}{c}_{\alpha }& s_{\alpha }\\ -s_{\alpha }& c_{\alpha }\end{array}\right)$ is orthogonal, so the sum of the kinetic terms of $n_L$ and $L^S$ is ${\bar{\nu }}_{L}i/\partial {\nu }_L+{\bar{N}}_{L}i/\partial N_L$ with no mixed terms. Combining with the right triplet gives ordinary Dirac field kinetic terms for $E^{+},N,e$, and only a left-handed kinetic term for the massless neutrino. This is orthogonal completion of mixed neutral lepton kinetic terms. The heavy charged field is $E^{+}$, as printed in the PDF, not the TeX's $N^{+}$.

To compute the charge commutator, the canonical anticommutation relations implyThe quartic terms cancel on reordering the fermions. Thus well-defined normal-ordered integrated bilinears obey fermionic bilinear charge algebra, $[Q_C,Q_D]=Q_{[C,D]}$. Left/right mixed commutators vanish since $P_L{P}_R=0$.

In either triplet, writeThe displayed current has a factor two relative to the projected bilinears, since $1\mp {\gamma }^5=2P_{L,R}$. Consequently its charges are $T^{+}=2(Q_K^L+Q_K^R)$ and $T^{-}=(T^{+}{)}^{†}$. The neutral combinations have canonical anticommutators because $c_{\alpha }^2+s_{\alpha }^2=1$. Applying the matrix commutator giveswhere $Q_{\ell }=e\int :(E^{+†}{E}^{+}-e^{†}e):$ is the electromagnetic charge in the lepton sector. The neutrino, heavy neutral lepton and singlet contribute zero. This establishes triplet charged-current closure on electromagnetic charge, with the requested proportionality independent of normalization choices for the weak Noether charges. The PDF correctly writes $J^{-}=(J^{+}{)}^{†}$; the TeX repeats $J^{+}$ incorrectly. Equal-time algebra does not require conservation of the separate fermionic weak currents after symmetry breaking. For full gauge-theory charges, the corresponding charged-boson contributions must also be included.

An order parameter distinguishes phases macroscopically. In an Ising model it is the local coarse-grained magnetization; the uniform mean is $M=V^{-1}\int \varphi (x)d^{D}x$. In a selected ordered phase $M\ne 0$, whereas the disordered phase has $M=0$ at zero conjugate field. The ordered signs are related by spin inversion symmetry. A finite zero-field system may average both signs to zero, so spontaneous symmetry breaking is defined using the thermodynamic limit before removing a symmetry-selecting field. For other systems the order parameter can be a vector, a density difference, or a tensor rather than a single magnetization.

Landau-Ginzburg theory makes the order parameter a spatial field and expands a local free energy consistent with the symmetries. For a short-range scalar system with spin inversion symmetry at $h=0$, a useful expression isThe linear source is the field conjugate to an order parameter. Odd terms at zero field are forbidden by spin inversion symmetry. The coefficients depend smoothly on control parameters away from singular fluctuation effects. With $u>0$ the quartic term suffices locally; if $u<0$, take $v>0$ for stability. The Landau approximation minimizes this functional; statistical field theory also integrates fluctuations about those minima.

For fractal geometry, let $C_n$ be the $n$th middle-third construction of the Cantor set. It contains $2^n$ intervals of length $3^{-n}$, with total length $(2/3{)}^n$. Each finite $C_n$ contains intervals, so its Hausdorff dimension is one. The limiting set $C={\bigcap }_n{C}_n$ has zero length andIndeed, $2^n$ intervals of scale $3^{-n}$ cover it, giving the upper dimension bound from $2^n{3}^{-ns}\to 0$ whenever $s>\log 2/\log 3$. For the lower bound, assign mass $2^{-n}$ to each level-$n$ interval. Any sufficiently small interval of length $r$ meets only a bounded number of construction intervals at the comparable scale, giving a bound $\mu (I)\le K{r}^{\log 2/\log 3}$. Covering $C$ by such intervals then forces their dimension-weighted sizes to sum to at least $1/K$, yielding the matching lower bound. The noninteger dimension is deducible from a simple iterative rule even though it is absent from every finite-stage set at arbitrarily fine resolution. A finite construction nevertheless exhibits the same scaling over intermediate resolutions. Thus the order of idealization and probing ever smaller scales matters: an actual finite-resolution pattern can have useful fractal behaviour without possessing the exact limiting dimension at all scales.

For a phase transition, finite-spin partition functions are finite sums of positive exponentials at real finite temperature and field. Their free energies are therefore analytic there. Genuine thermodynamic nonanalyticity can appear only after a suitable infinite-system limit. The Curie–Weiss model gives an explicit demonstration. With $s_i=\pm 1$, takeGrouping configurations by $m=N^{-1}{\sum }_i{s}_i$ and using their binomial multiplicities gives the limiting free energy as the minimum ofThe number of possible magnetizations grows only linearly with $N$, so their largest exponential contribution determines the limiting free energy. Differentiating the variational expression givesAt $h=0$ and $\beta J>1$, its two stable minima have magnetizations $\pm m_{\ast }\ne 0$. Every finite zero-field system instead has exactly zero mean magnetization by spin-inversion symmetry. ConsequentlyThis noncommutation of limits explains how spontaneous symmetry breaking emerges. Below the critical temperature, the limiting free energy has a cusp in $h$ at zero, while finite systems have a sharp but smooth crossover. Away from the critical point, the relative weights of the two phases behave approximately as $e^{2\beta Nh{m}_{\ast }}$; the rounded magnetization is correspondingly about $m_{\ast }\tanh (\beta Nh{m}_{\ast })$. Exact singularity belongs to the idealized limit, whereas the narrow crossover and long-lived phases can be physically useful at finite size. The theoretical reduction is the derivation and its approximation control, not a claim that a finite analytic partition function is already nonanalytic.

For superselection, consider an infinite spin lattice with quasi-local observable algebra $\mathcal{A}$. These are zero-temperature equilibrium phases of a ferromagnetic Ising spin Hamiltonian with nearest-neighbour interactions $-J{\sigma }_z^{(i)}{\sigma }_z^{(j)}$. At finite $N$, the all-up and all-down vectors belong to one Hilbert space, and the coherent stateis distinguishable from the corresponding mixture by a suitably chosen global observable. But if $A_R$ acts on only $R<N$ spins, its off-diagonal matrix element between these vectors vanishes: the untouched spins contribute $⟨\uparrow |\downarrow ⟩=0$. Hence its expectations already equal those of the equal-weight mixture, independent of $\theta$. Increasingly global observables can recover the phase at finite $N$, but the infinite product of spin flips is not a norm limit of finite-support observables and is absent from $\mathcal{A}$.

In the infinite-volume limit the up and down phases have disjoint state representations. The averaged magnetization $M_N=N^{-1}{\sum }_i{\sigma }_z^{(i)}$ commutes asymptotically with any fixed local observable: $\Vert [M_N,A_R]\Vert \le 2R\Vert A_R\Vert /N$. In the phase representations its limiting values are $+1$ and $-1$, distinguishing sectors. Finite spin-flip vectors form a dense set on which these averages tend to their respective scalar values, so boundedness extends the strong limits to the phase Hilbert spaces. An intertwiner between the representations would obey $(-I)T=T(+I)$ and must therefore vanish, explaining their disjointness. In the combined representation the observable algebra is block-diagonal between them, so no quasilocal measurement detects their relative phase. This yields an emergent superselection rule relative to the chosen observable algebra. It is not a derivation of physical wave-function collapse. Fractal dimension, a thermodynamic singularity and sector separation each exhibit new limiting structure while retaining an explicit microscopic construction.

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

Use spin one-half states $|\uparrow ⟩,|\downarrow ⟩$ with $S^z=\pm 1/2$, so ${\sigma }^i=2S^i$ and $H=-J{\sum }_n({\sigma }_n^x{\sigma }_{n+1}^x+g{\sigma }_n^z)$. At $g=0$, every ferromagnetic bond has its lowest energy $-J$ when adjacent $x$ spins agree. On the connected periodic chain the normalized ground states arewhere $|\pm x⟩=(|\uparrow ⟩\pm |\downarrow ⟩)/\sqrt{2}$. Their span is the two-dimensional ground-state subspace; any normalized superposition is another ground state. They are exchanged by the global spin-flip discrete symmetry generated by ${\prod }_n{\sigma }_n^z$.

As $g\to +\infty$, the field dominates and selectsThe limiting polarized state is unique. At large finite $g$, exchange admixes virtual spin flips, so the product state is the limiting wavefunction rather than an exact finite-field eigenstate. The bulk expansion below gives $E_0/N=-J[g+1/(4g)+O(g^{-3})]$. The exact twofold degeneracy stated at $g=0$ should be distinguished from thermodynamic spontaneous symmetry breaking in the ordered phase: a finite chain at nonzero field can have a split pair of symmetry eigenstates.

A homogeneous static saddle minimizes $f(\rho )=-\mu \rho +g{\rho }^2/2$ for $\rho =|\psi {|}^2\ge 0$. With $g>0$ and $\mu >0$,The mean-field Bose-Einstein condensate has a macroscopic coherent occupation of the uniform mode. Its continuous particle-number phase symmetry is $\psi \mapsto e^{i\chi }\psi$, the circle group $U(1)$ generated by particle number. Choosing one phase breaks it, and the saddle manifold is a circle. For $\mu \le 0$, the constrained minimum is the vacuum ${\rho }_0=0$, so the condensed saddle requires the stated positive-density regime.

Spontaneous symmetry breaking is understood through a phase-selected thermodynamic description: an exact finite-volume number eigenstate has vanishing field expectation. The resulting Goldstone boson is a gapless phase/sound mode, while phase gradients carry the superfluid velocity $\mathbf{v}=\nabla \varphi /m$. This is the saddle-point conclusion requested, not a proof of true condensate order in every dimension or temperature. Long-wavelength phase fluctuations can invalidate the assumed order, as quantified below.

Let $V(\varphi )=m^2{\varphi }^2/2+\lambda ({\varphi }^2{)}^2/4$, where ${\varphi }^2={\sum }_i{\varphi }_i^2$. Its stationary-point equation is $(m^2+\lambda {\varphi }^2){\varphi }_i=0$.

For $m^2>0$, the unique minimum is $\varphi =0$. The quadratic mass term is $m^2{\delta }_{ij}$, so there are $n$ real scalar fields, all of mass $\sqrt{m^2}$, and the global special orthogonal group $SO(n)$ remains unbroken.

For $m^2<0$, write $v^2=-m^2/\lambda$. The vacuum manifold is the sphere ${\varphi }^2=v^2$. Choose $⟨\varphi ⟩=v{e}_n$. The transformations preserving this vacuum expectation value rotate the first $n-1$ components, giving $SO(n)\to SO(n-1)$. At this vacuum the potential's Hessian matrix isWriting $\varphi =({\pi }_1,\dots ,{\pi }_{n-1},v+\eta )$, the radial mode $\eta$ therefore hasThese are one massive radial scalar and $n-1$ physical Goldstone bosons. The Goldstone theorem applies to the broken global generators; the flat angular directions of the vacuum manifold explain their zero masses. The spontaneous symmetry breaking is a choice of vacuum, not a change in the invariant Lagrangian. For $n=1$ there are no continuous angular directions: only the radial scalar remains, and the potential's discrete $\varphi \mapsto -\varphi$ symmetry is broken.

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.

Take the couplings in the Hamiltonian to include inverse temperature, so the Boltzmann factor is $e^{-H}$. If physical energies are used instead, first replace $J,L,D,K$ by their products with $1/(k_{B}T)$. Split the on-site term equally between its two bonds. For spin states $s,t\in \{1,0,-1\}$ the spin-chain transfer matrix has entriesUsing the specified coupling coordinates, in the order $(1,0,-1)$ this isFor example, $e^{D/2}=x/z$ and $e^{-J+L+D}=y/z$. This is the nearest-neighbour Blume–Emery–Griffiths model with an additive constant and the paper's sign convention for $D$.

Summing the periodic spin chain gives the partition function $Z_N=\mathrm{tr}{W}^N={\sum }_{a=1}^3{\lambda }_a^N$. For finite real couplings all entries of $W$ are strictly positive; the Perron–Frobenius theorem gives a unique positive dominant eigenvalue with ${\lambda }_1>|{\lambda }_a|$ for $a\ne 1$. Since $W$ is a symmetric matrix, all its eigenvalues are real. Thus in the thermodynamic limit,The positivity condition is stronger and more useful than mere ordering by signed value: subdominant eigenvalues can be negative. Also, the printed strict ordering between the other two is not guaranteed for all couplings. For instance, $J=L=D=0$ makes $W$ the positive constant matrix $e^{K}11^T$ with two equal zero eigenvalues. Degeneracy there does not affect the largest-eigenvalue limit; no explicit generic eigenvalues are needed.

For spin magnetization, introduce a dimensionless field $h$ through $-h{\sum }_i{\sigma }_i$, or insert $S=\mathrm{diag}(1,0,-1)$ into the trace. Thenwhere $v_1$ is a normalized eigenvector of the largest eigenvalue and $W_{st}(h)=e^{h(s+t)/2}{W}_{st}(0)$. Spin inversion symmetry gives $v_{1,+}=v_{1,-}$, since the positive eigenvector of the largest eigenvalue is unique. Thereforeat zero field, both at finite $N$ and in the finite-coupling thermodynamic limit. The finite-$N$ result follows directly by pairing each configuration with its spin-reversed partner. Eigenvalues at a single fixed field do not determine a general observable: a field derivative of the largest eigenvalue, or its eigenvector, is required. Positivity and the real analytic dependence on the couplings exclude a finite-temperature spontaneous symmetry breaking transition in this one-dimensional finite-range chain; singular zero-temperature coupling limits require separate treatment.

For even $N$, spin decimation on alternate sites sums the middle spin of each two-bond segment, hence the coarse spin-chain transfer matrix is $W^{\prime }=W^2$. DefineDirect multiplication gives $W^2=e^{2K}{z}^{-2}\left(\begin{array}{ccc}A& B& C\\ B& E& B\\ C& B& A\end{array}\right)$. Matching its entry ratios to the original parameterization yields the spin-1 chain decimation recursionThe remaining overall positive factor is absorbed into $K^{\prime }$. Keeping that factor preserves the free energy as well as normalized spin probabilities. Indeed $\mathrm{tr}(W^2{)}^{N/2}=\mathrm{tr}{W}^N$, exactly; on an odd ring an unmatched boundary segment needs separate handling rather than assuming a uniform two-site block decomposition.

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

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

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

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

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

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

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

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