Use coordinates $x^0=t$, $x^1=x$ and first compute the mixed canonical stress-energy tensor supplied by spacetime translations:Define the energy density $\mathcal{E}={\varphi }_{1,t}{\varphi }_{2,t}+c^2{\varphi }_{1,x}{\varphi }_{2,x}$ and $\mathcal{P}={\varphi }_{1,t}{\varphi }_{2,x}+{\varphi }_{2,t}{\varphi }_{1,x}$. Direct substitution givesThe two wave equations implyFor example, differentiating $\mathcal{E}$ and replacing each ${\varphi }_{j,tt}$ by $c^2{\varphi }_{j,xx}$ gives exactly $c^2{\partial }_x\mathcal{P}$ by the product rule. These are respectively the two component equations $\boxed{{\partial }_{\mu }{T}^{\mu }{}_{\nu }=0}$, establishing stress-energy conservation explicitly.

If the given metric $g^{\mu \nu }=\mathrm{diag}(1,-1)$ is used to raise the second index in these unrescaled coordinates, the result is $T^{\mu \nu }=\left(\begin{array}{cc}\mathcal{E}& -\mathcal{P}\\ -c^2\mathcal{P}& \mathcal{E}\end{array}\right)$, still conserved. To display the symmetric relativistic form with light speed $c$, choose $X^0=ct$, $X^1=x$ instead. With ${\partial }_{X^0}=c^{-1}{\partial }_t$, the corresponding tensor isThis states the coordinate convention explicitly rather than setting $c=1$ in some components but not others.

The infinitesimal connected internal symmetry of a classical field theory is $\delta {\varphi }_1=ϵ{\varphi }_1$, $\delta {\varphi }_2=-ϵ{\varphi }_2$. For this exact invariance of the Lagrangian density, the Noether current is $J^{\mu }={\sum }_j{\mathcal{L}}_{{\partial }_{\mu }{\varphi }_j}\delta {\varphi }_j/ϵ$. In coordinates $(t,x)$ it isVarying the two fields gives the two wave equations ${\varphi }_{1,tt}-c^2{\varphi }_{1,xx}=0$ and ${\varphi }_{2,tt}-c^2{\varphi }_{2,xx}=0$. Direct differentiation cancels the first-derivative cross terms:The Noether charge $Q=\int J^{0}dx$ is constant when the spatial boundary flux vanishes. The disconnected transformations from part (a) do not have independent infinitesimal generators of this type.

Apply the linear transformation of the scalar fields to either time- or space-derivative product. Its coefficients become $\alpha \gamma$ multiplying the square of the first derivative, $\beta \delta$ multiplying the square of the second, and $\alpha \delta +\beta \gamma$ multiplying their cross product. Invariance for arbitrary fields therefore requiresThe complete real solution isThere are four connected components, distinguished by the diagonal/antidiagonal choice and the sign of $r$. Only the positive diagonal family is in the identity component; write $r=e^s$ to obtain ${\varphi }_1\mapsto e^s{\varphi }_1$, ${\varphi }_2\mapsto e^{-s}{\varphi }_2$. This is the reciprocal rescaling symmetry of two cross-coupled scalar fields. The antidiagonal transformations exchange the fields as well as rescale them.

For a homogeneous saddle-point approximation define $r=a(T-T_c)>0$. With the cubic source coupling, minimizing the Landau free energy givesIts curvature $2r+6b{m}^2$ is positive, so there is a unique local minimum and global minimum. The solution near $B=0$ hasbecause substituting this leading term makes the cubic self-interaction of order $B^9$. Therefore the cubic source coupling and third-order susceptibility giveThus the first nonzero response diverges like $(T-T_c{)}^{-1}$ from above. The factor $3!$ is essential: the coefficient of $B^3$ is not itself the third derivative. This is susceptibility to the specified cubic source variable, rather than the usual third-order response to a linearly coupled field.

The coarse-grained order parameter is a real scalar field because the microscopic spin has one component. At $B=0$, spin inversion symmetry takes $m\mapsto -m$, so an analytic local Landau free energy contains only even powers. Near a continuous phase transition, the quadratic coefficient changes sign linearly, represented by $a(T-T_c)$ with $a>0$. A positive quartic coefficient stabilizes the density when the quadratic term becomes negative, giving the two minima $m=\pm \sqrt{a(T_c-T)/b}$ below $T_c$.

The ferromagnetic interaction favours neighbouring spins pointing alike, so spatial variation costs energy; its leading long-wavelength form is the positive gradient term $c({\partial }_{x}m{)}^2$. The source term $-Bm$ favours alignment with the applied magnetic field and explicitly breaks spin inversion symmetry at fixed nonzero $B$. Finally $f_0(T)$ is the spin-independent background free-energy density.

The quartic sign can also be checked from part (b): its expansion containswhose quartic coefficient is positive. This is consistent with the proposed Landau-Ginzburg theory, not a tricritical truncation requiring a stabilizing sixth-order term. The microscopic square lattice is two-dimensional; the written gradient describes a profile depending on one coordinate. For general profiles, replace it by $c[({\partial }_{x}m{)}^2+({\partial }_{y}m{)}^2]$.

The small-argument Taylor series of the single-site response isLet $K=4\beta J$. For a small nonzero order parameter at zero field, the self-consistency equation becomesSince $K_c=3/2$ and $K=K_c/(1+t)$, where $t=T/T_c-1$, this givesConsequentlyThe positive branch is the expression printed in the question, while its negative partner is required by spin inversion symmetry. The square-root behavior gives the mean-field critical exponent $1/2$ for the spontaneous magnetization.

At zero field, write the right-hand side of the self-consistency equation as $G_T(S)=g(4\beta JS)$, with $g(x)=2\sinh x/[1+2\cosh x]$. The origin is always a solution, and the slope there isThe critical crossing occurs when this equals one, givingThis is a continuous phase transition in the mean-field approximation. For $x>0$,so the positive branch is a strictly concave function, while $g$ is odd and bounded by one. For $T\ge T_c$, the slope at zero is at most one and strict concavity puts $G_T(S)$ below $S$ for every $S>0$, leaving only the unmagnetized solution. For $T<T_c$, it initially lies above $S$ but has $G_T(1)<1$; hence there is one positive intersection and, by spin inversion symmetry, one negative intersection.

The curvature of the Helmholtz free energy in the mean-field approximation is $4J[1-G_T^{\prime }(S)]$ per site. It is negative at $S=0$ below $T_c$ and positive at the nonzero intersections, identifying the ordered minima and the unstable zero solution.

The critical temperature and branches here belong to the mean-field approximation. In the exact finite system, spin inversion symmetry forces zero average magnetization at zero field; selecting a broken-symmetry phase requires the thermodynamic limit with a vanishing field of specified sign.

The self-consistency equation requires that the mean spin computed using the single-site Boltzmann distribution equals the spin $S$ used to form the effective field. With $h=4JS+B$,ThusEquivalently, differentiating part (b) gives $N^{-1}\partial F_{\mathrm{MF}}/\partial S=4J[S-⟨s{⟩}_h]$, so stationarity of the Helmholtz free energy in the mean-field approximation is exactly this condition. Among multiple stationary points one must select the stable minimum.

The square lattice has coordination number of a lattice $z=4$. The mean-field decoupling isAccounting for the factor $1/2$ in the directed-neighbour sum gives the Hamiltonian in the mean-field approximation:The constant term is necessary to correct the double counting of the interaction. Let $h=4JS+B$. The one-site partition function isso $Z_{\mathrm{MF}}=e^{-2\beta JN{S}^2}{Z}_1^N$. The Helmholtz free energy is thereforeThe equilibrium mean-field approximation is obtained by minimizing this expression over the allowed $S$. The $s=0$ state supplies the extra $1$ in $Z_1$; omitting it would instead describe the two-state Ising model. This is the Blume-Capel mean-field transition at zero single-ion anisotropy.

The mean-field approximation replaces the fluctuating environment of each microscopic variable by a uniform average field. For spins, write $S_i=S+\delta S_i$ and drop products $\delta S_i\delta S_j$ between different sites. The resulting single-site problem can be solved independently, and its average spin must reproduce the assumed order parameter $S$: this is the self-consistency equation.

It is reliable when correlated thermal fluctuations about the mean are small, for example with many interacting neighbours, long-range interactions or high spatial dimension. Near a thermodynamic critical point the quantitative test is the Ginzburg criterion. For short-range scalar ordering the upper critical dimension is four, with logarithmic qualifications at four; a two-dimensional nearest-neighbour lattice is below this and its asymptotic critical region of a phase transition is not described accurately by simple mean-field theory.^[ref]

The dispersion relation shows that an infinite rod has unstable long-wave modes. A finite rod excludes sufficiently long wavelengths, so stability depends on its length and the precise end conditions. The fourth spatial derivative requires two conditions at each end; fixing transverse positions alone does not specify whether slopes are constrained.

First take displacement-fixed ends with free rotation. Setting $\varphi (0)=\varphi (l)=0$ leaves slope variations free; the bending boundary term in the first variation then gives the corresponding natural boundary conditions on the slope, ${\varphi }_{xx}(0)={\varphi }_{xx}(l)=0$. The spatial normal modes are $\sin (n\pi x/l)$, andThe first mode destabilizes first. ConsequentlyAt equality, the first mode has zero restoring force and is marginal, with a static component and a possible linear-in-time component in the linear theory. Every integer $n<l\sqrt{\alpha /\beta }/\pi$ gives an exponentially growing mode.

If "fixed" instead means clamped, the clamped boundary condition is $\varphi ={\varphi }_x=0$ at both ends. The critical compression is thenFor example, the static critical shape $1-\cos (2\pi x/l)$ obeys the four clamp conditions and the zero-frequency equation at this load. To establish that it is the first threshold, write $v={\varphi }_x$: it has zero endpoints and zero mean, and periodic extension gives ${\int }_0^l{v}_x^{2}dx\ge (2\pi /l{)}^2{\int }_0^l{v}^{2}dx$ by the zero-mean periodic Poincaré inequality, with equality at $v\propto \sin (2\pi x/l)$. This proves the critical Rayleigh quotient. Both cases are instances of the compressed Euler-Bernoulli beam and its endpoint-dependent buckling threshold. The source does not specify the endpoint rotation, so both interpretations are stated explicitly.

For this density the Euler-Lagrange field equation with second derivatives givesInsert a plane wave $\varphi =A{e}^{i(kx-\omega t)}$. The second and fourth spatial derivatives give $-k^2$ and $k^4$, respectively, soFor propagating modes, $k^2>\alpha /\beta$ and the phase velocity has magnitude $|\omega /k|=\sqrt{(\beta k^2-\alpha )/\rho }$, which depends on the wave number. Thus the waves are dispersive. For $0<k^2<\alpha /\beta$, the squared frequency is negative: the corresponding solutions grow or decay exponentially rather than propagate with a real frequency. The uniform $k=0$ displacement is a zero-frequency shift mode.

A constant shift $\varphi \mapsto \varphi +ϵ$ leaves every derivative unchanged, and the density has no undifferentiated $\varphi$. Thus it is exactly invariant under this scalar-field shift symmetry. Localize the shift parameter temporarily to $ϵ(x,t)$:Integrating the last term by parts once identifies the Noether current for a constant shift with second derivatives:The Euler-Lagrange field equation is $\rho {\varphi }_{tt}+\alpha {\varphi }_{xx}+\beta {\varphi }_{xxxx}=0$, precisely the continuity equation ${\partial }_t{j}^0+{\partial }_x{j}^x=0$. Therefore the conserved quantity, when the endpoint current vanishes or the field has suitable periodic/decaying conditions, isthe total transverse momentum. More explicitly $dQ/dt=-[\alpha {\varphi }_x+\beta {\varphi }_{xxx}{]}_{\mathrm{ends}}$. The fixed endpoints introduced later can exert transverse forces and break the global shift symmetry of the boundary problem, so vanishing endpoint current must not be inferred merely from fixed displacement.

Interpret $\varphi$ as the small transverse displacement of the elastic beam. The first term is the positive kinetic energy density $\rho {\varphi }_t^2/2$, with $\rho$ the mass per unit length. The Lagrangian density is kinetic minus potential energy, so the potential density encoded by the other two terms isFor a compressed rod, transverse deflection releases compression energy: the small-slope geometric change in length is proportional to ${\varphi }_x^2/2$, so the applied compressive load lowers the potential by $\alpha {\varphi }_x^2/2$. This explains the positive sign of the $\alpha$ term in the Lagrangian. Tension would give the opposite sign and a restoring gradient term.

Bending costs elastic energy. At small slopes the curvature is ${\varphi }_{xx}$, so positive bending stiffness $\beta$ gives potential $\beta {\varphi }_{xx}^2/2$, and hence a negative term in the Lagrangian. For a uniform slender beam $\beta$ is its flexural rigidity. The competition between destabilizing compression and stabilizing bending is responsible for Euler buckling of an elastic filament.

Vary the action $S=\int dtdx\mathcal{L}$ with compactly supported $\delta \varphi$, or impose vanishing $\delta \varphi$ and its first derivatives on the spacetime boundary. The first variation isOne integration by parts on the first-derivative term and two on the second-derivative term remove all derivatives from $\delta \varphi$. The boundary terms vanish, leavingThe fundamental lemma of the calculus of variations gives the Euler-Lagrange field equation with second derivatives:Repeated indices run over time and space. The coefficient of the symmetric second derivative is defined by the displayed variation; if ${\varphi }_{tx}$ is instead treated as a single independent mixed derivative variable, it is counted once rather than twice.

A common rotation changes neither the separation nor the spring energy. Its infinitesimal direction is $(1,1)$, exactly the zero-frequency normal mode. Thus the continuum of rotated equilibrium points requires that zero mode; a nonzero restoring frequency in that direction would contradict rotational symmetry.

The second mode has $I_1{\eta }_1+I_2{\eta }_2=0$, so it changes the relative angle while leaving the inertia-weighted common angle unchanged. It therefore carries zero perturbation in the conserved angular momentum associated with common motion. Indeed the coordinatesseparate the Lagrangian exactly:The cyclic coordinate $q$ is free, while $\delta$ has restoring or inverted curvature. This also explains why all members of either equilibrium family have the same frequencies. Reflection symmetry makes the potential even in $\delta$, and the stable oscillatory mode is consistent with time reversal in classical mechanics through its two opposite frequency signs.

Let ${\delta }_0={\theta }_0-{\phi }_0$ be $0$ or $\pi$, and write the angular displacements as ${\eta }_1,{\eta }_2$. To quadratic order the restoring potential energy is $\kappa \cos {\delta }_0({\eta }_1-{\eta }_2{)}^2/2$. The linearized normal mode equation isThe determinant gives ${\omega }^2[{\omega }^2{I}_1{I}_2-\kappa \cos {\delta }_0(I_1+I_2)]=0$. Thus the normal modes of two spring-coupled coaxial rotors areFor the aligned family, the relative mode oscillates with positive frequencyFor the antipodal family, ${\omega }_{\mathrm{rel}}^2=-{\Omega }^2$, so that displacement has solutions $e^{\pm \Omega t}$ and the equilibrium is unstable.

The common zero mode has ${\eta }_1={\eta }_2=C+Dt$. It describes a neutral displacement or uniform common rotation, not a restoring oscillation. Thus the aligned family is stable with respect to relative motion and has orbital stability modulo common rotations. A particular static point is not strictly Lyapunov stable against unrestricted small velocity perturbations: arbitrarily small nonzero $J$ eventually moves the common angle far from that point. With $J=0$, the common angle is fixed and the aligned relative equilibrium point is stable. This distinguishes the symmetry-induced drift from the exponential instability of the antipodal family.

The angular Euler-Lagrange equations areA static equilibrium point has both angular velocities zero and $\sin (\theta -\phi )=0$. Hence the complete equilibrium families, with angles taken modulo $2\pi$, areThe first family is aligned in azimuth and minimizes the spring length; the second is antipodal and maximizes it. Their potentials are $k[(a-A{)}^2+h^2]/2$ and $k[(a+A{)}^2+h^2]/2$, respectively. Common rotations generate a whole family rather than isolated equilibria.

The two continuous, time-independent symmetries are time translation and a common rotation of both angles. Their Noether charges are respectively the total energy and vertical angular momentum:For the second symmetry, $\theta \mapsto \theta +ϵ$ and $\phi \mapsto \phi +ϵ$, so the Noether theorem gives the sum of their conjugate momenta. Independent rotations are not symmetries of the spring coupling.

There are also discrete symmetries: simultaneous reflection $(\theta ,\phi )\mapsto (-\theta ,-\phi )$, and time reversal in classical mechanics, which reverses both angular velocities. These leave the dynamics invariant but provide no additional infinitesimal Noether charges. Combined with common rotations, the reflections give the full planar orthogonal symmetry of the angular geometry. If $I_1=I_2$, interchange of the two angular coordinates is an additional discrete symmetry of the reduced dynamics.

There is a further useful time-dependent symmetry, up to a total derivative, on a continuous lift of the angles: adding the same constant angular velocity, $(\theta ,\phi )\mapsto (\theta +vt,\phi +vt)$. It changes $L$ byIts conserved quantity is $\boxed{I_1\theta +I_2\phi -tJ}$. It is the integration constant for free common rotation, rather than a new globally single-valued angle charge on the configuration torus. The globally defined continuous charges are $E$ and $J$ above.

Let $\theta$ and $\phi$ be the lower and upper azimuthal angles, measured from a common fixed direction. The squared spring length isBy Hooke's law, a spring with zero natural length has potential energy $k{d}^2/2$. ThereforeConstant gravitational energies can be omitted. It is useful to define the moments of inertia $I_1=m{a}^2$, $I_2=M{A}^2$ and the angular coupling $\kappa =kaA$. The part of the potential that affects motion is simply $-\kappa \cos (\theta -\phi )$. In particular, $h$ changes the constant spring energy but not the angular equations, because the spring has zero natural length.