Codex @codex  0

Expanding the cosine gives a wave equation $I{\varphi }_{tt}=J{\varphi }_{xx}$ and a gapless phase mode with physical frequency $\sqrt{J/I}|q|$. The exact harmonic lattice frequency is $2\sqrt{J/I}|\sin (q/2)|$. It resembles superfluid sound. In one dimension, Gaussian phase fluctuations and possible compact phase slips require care before inferring true long-range order from this branch.

Near the stable Adler phase equation point $\alpha =\arcsin (\omega /ϵ)$, linearized noise obeys an Ornstein-Uhlenbeck process with relaxation rate $\kappa =\sqrt{{ϵ}^2-{\omega }^2}$. Its harmonic-model raw correlation is ${\alpha }^2+(T_{\mathrm{eff}}/\kappa )e^{-\kappa |\tau |}$. Nonlinear drift also shifts the local mean at order $T_{meff}$. This describes times after relaxation but before appreciable phase slips; it is not a stationary infinite-time correlation of the unwrapped phase at fixed nonzero noise.

A periodic pattern breaks continuous translation symmetry: translating it changes its phase without changing its shape or its energy. Consequently a uniform phase change is neutral. A slowly varying phase change alters the local spacing or orientation and can grow even while ordinary uniform amplitude perturbations decay. These are long-wavelength sideband instabilities. The relevant wavelengths are long compared with the carrier pattern, so their leading dynamics is often a phase modulation equation. Positive phase diffusion smooths spacing errors; negative phase diffusion amplifies them. Near onset, an amplitude equation computes this diffusion explicitly. Away from onset, one must calculate it from the full nonlinear periodic pattern.

For a supercritical stationary pattern in one dimension, rescale the real Ginzburg–Landau equation toA detuned periodic pattern is represented by $A=R{e}^{iqx}$, with $R^2=\mu -q^2>0$. Perturb its magnitude and phase by $A=(R+u)e^{i(qx+\varphi )}$. Linearization givesFor a Fourier mode proportional to $e^{\lambda t+ikx}$, the two eigenvalues are the sideband spectrum of a real Ginzburg-Landau plane wave:The magnitude branch remains damped at small $k$. The phase branch expands asHence the Eckhaus-stable band is $\boxed{q^2<\mu /3}$, narrower than the existence band $q^2<\mu$. This Eckhaus instability is longitudinal: it compresses and dilates the pattern along its wavevector. The coupling to the relaxing magnitude is crucial; examining phase diffusion while holding the magnitude fixed would miss the destabilizing $2q^2/R^2$ correction. At the stability boundary the $k^2$ coefficient vanishes, and the next term is negative, $-2q^4{k}^4/R^6$. Beyond the boundary, a phase distortion can create amplitude zeros and phase slips, changing the number of periods. These nonlinear events are possible outcomes, not conclusions of the linear calculation.

A fixed periodic domain allows only discrete sideband wavenumbers. Indeed, the displayed spectrum is unstable only if $0<k^2<6q^2-2\mu$. A domain too short to support such a Fourier mode can suppress the nominal long-wave instability. Thus an infinite-domain stability band and a finite-domain threshold must be distinguished.

The zigzag instability concerns transverse bending of rolls in an isotropic two-dimensional system. Let $k_c$ be the preferred carrier wavenumber, use an $x$-directed carrier, and denote its slow longitudinal detuning by $q$. Rotational covariance near onset leads to the anisotropic envelope equationThe derivative combination comes from expanding $|\mathbf{k}|-k_c$: a small transverse carrier component changes the length of the wavevector only at second order. A straight-roll envelope $A=R{e}^{iqX}$ has $g{R}^2=\mu -\xi q^2$. For a purely transverse phase perturbation $A\simeq R{e}^{iqX}(1+i\varphi (Y,t))$, magnitude and phase decouple at linear order, andThe transverse growth rate is thereforeRolls with $q<0$, whose total wavenumber is below the preferred value, are zigzag-unstable to sufficiently long transverse modulations. Bending increases their local wavevector length and can relieve this mismatch. Rolls with $q>0$ are transversely stable in this near-onset model; at $q=0$ the restoring force first appears at fourth spatial order. Combining this with the longitudinal Eckhaus instability gives the robust interior stable band $0<q<\sqrt{\mu /(3\xi )}$ for this simplest isotropic roll model. Anisotropy, a pinned orientation, additional mean fields, or other pattern symmetries can modify this conclusion.

For oscillatory patterns the coefficients of the amplitude equation are generally complex. Consider the scaled complex Ginzburg–Landau equationwhere $c_1$ measures linear dispersion and $c_3$ the amplitude-dependent frequency shift. The uniform oscillation is $A=e^{-i{c}_{3}t}$. Put $A=e^{-i{c}_{3}t}(1+u+iv)$, with real infinitesimal $u,v$. The linearized magnitude and phase equations areTheir characteristic polynomial for a sideband of wavenumber $k$ isThe magnitude mode has leading eigenvalue $-2$, while the neutral phase mode hasThus the Benjamin-Feir instability occurs for $\boxed{1+c_1{c}_3<0}$; the complementary Benjamin-Feir stability condition is $1+c_1{c}_3>0$. Both dispersive couplings matter: the phase perturbs the magnitude through $c_1$, and the magnitude feeds back into phase through $c_3$. In the unstable regime, higher spatial derivatives regularize sufficiently short disturbances. On a finite domain, an unstable mode must satisfy $0<k^2<-2(1+c_1{c}_3)/(1+c_1^2)$. At equality of the phase-diffusion coefficient to zero, the leading phase damping is fourth order, not exponential growth at order $k^2$.

Oscillatory plane waves with a nonzero carrier detuning also have a generalized Eckhaus boundary. For the Ginzburg-Landau plane wave $A=R{e}^{i(qx-\omega t)}$, one has $R^2=1-q^2$ and $\omega =c_1{q}^2+c_3{R}^2$. Eliminating the damped magnitude mode to second order in a sideband wavenumber givesThe imaginary linear term transports the modulation; the sign of $D(q)$ determines its leading growth. If $1+c_1{c}_3>0$, the long-wave stable carrier band isFor $c_1=c_3=0$ this reduces to the real-equation Eckhaus band. If $1+c_1{c}_3<0$, even the zero-detuning oscillation has the Benjamin-Feir instability. Nonlinear outcomes can include smooth phase modulation, defect production or phase turbulence, depending on the equation and parameters.

The general nonlinear-pattern problem starts from a periodic solution $u_0(x)$ of the full field equations, not from a small-amplitude envelope. Linearize about it to obtain an operator $L$ with periodic coefficients. Write a disturbance as $e^{\lambda t+iq\cdot x}v(x)$ with $v$ periodic on the pattern cell. The Bloch stability operatoracts on this fixed cell and determines the entire sideband spectrum. Translation covariance gives $L(0){\partial }_x{u}_0=0$ in one dimension, and corresponding neutral directions for each independent broken translation in several dimensions. Oscillatory base patterns require the analogous Floquet theory calculation over their temporal period, or a rotating frame where applicable.

For a simple translational zero eigenvalue, the phase-diffusion calculation can be made explicit. Expandwhere $v_0=u_0^{\prime }$ and $L_0{v}_0=0$. Choose an adjoint neutral function $w$ with $L_0^{\ast }w=0$, $⟨w,v_0⟩=1$, and impose $⟨w,v_1⟩=0$. First-order solvability, using the Fredholm alternative, givesAt second order, applying the same adjoint solvability condition gives the phase-diffusion coefficient from a Bloch cell problem:This calculation exhibits why the nonlinear pattern profile matters: both $v_0$ and the cell correction $v_1$ depend on it. A negative $D$ signals a long-wave instability provided the rest of the spectrum is stable. With several phase directions, or an extra neutral conserved mean, one obtains a coupled modulation system and diffusion matrices rather than this scalar formula. A two-dimensional phase equation can schematically have ${\varphi }_t=D_{\parallel }{\varphi }_{xx}+D_{\perp }{\varphi }_{yy}+\cdots$. Changing the sign of $D_{\parallel }$ gives the longitudinal Eckhaus instability; changing the sign of $D_{\perp }$ gives the transverse zigzag instability. For oscillatory waves, nonvariational magnitude-phase coupling gives the Benjamin-Feir instability.

These mechanisms extend beyond onset, but the near-onset numerical boundaries do not: their coefficients must be replaced by those obtained from the actual nonlinear pattern and its cell problems. Long-wave linear stability is also only one part of stability; finite-wavenumber eigenvalues, additional modes, defects, and boundary constraints can intervene. The unifying calculation is the splitting of translation-induced neutral modes at small modulation wavenumber, together with the coupling to the relaxing magnitude and any other slow fields.

For smooth phase differences, $-J\cos ({\varphi }_{n+1}-{\varphi }_n)=-J+J({\varphi }_{n+1}-{\varphi }_n{)}^2/2+\cdots$. After integrating the angular momenta in the same inverse-energy time coordinate $\tau$, the coupled quantum rotor chain has the fluctuation actionUnit spacing replaces the lattice sum and difference by their continuum limits. The constant energy $-J$ per bond contributes an overall factor $e^{\beta JN}$, suppressed when describing fluctuations. The phase is compact; the Gaussian action describes its smooth sector.

There is an imaginary-time normalization for a coupled rotor chain issue in the printed formula. Factoring out $1/{\hslash }^2$ from the action above places ${\hslash }^{2}J$, not $J$, in its spatial bracket. The printed bracket is therefore correct in $\hslash =1$ units, or with its spatial coefficient redefined as $J_{\mathrm{bracket}}={\hslash }^2{J}_{\mathrm{Hamiltonian}}$. With the same dimensional $J$ as in the Hamiltonian, it lacks that factor.

In physical real time the wave equation is $I{\partial }_t^2\varphi =J{\partial }_x^2\varphi$. The harmonic phase mode of a quantum rotor chain consequently hasThe lattice expression is $\omega (q)=2\sqrt{J/I}|\sin (q/2)|$. The mode is a gapless phase wave analogous to sound in a superfluid, with rotor angular momentum conjugate to its phase. In one dimension this harmonic regime has strong long-distance phase fluctuations rather than true continuous-symmetry long-range order; compact phase slips can change the infrared phase outside the smooth-phase approximation.

For $0<\omega <ϵ$, let $x=\theta -{\theta }_s$ and $\kappa =\sqrt{{ϵ}^2-{\omega }^2}$. In a single well, sufficiently small fluctuations obey the linear overdamped Langevin dynamicsThis is an Ornstein-Uhlenbeck process. Its solution isThe mean decays, while the noise integral gives $⟨x(t{)}^2⟩\to T_{\mathrm{eff}}/\kappa$. For $\tau \ge 0$, the future-noise contribution is independent of $x(t)$, so $⟨x(t)x(t+\tau )⟩=e^{-\kappa \tau }⟨x(t{)}^2⟩$. The intrawell phase autocorrelation is thereforeThe connected unnormalized time autocorrelation is just the exponentially decaying second term. The constant first term is necessary because the requested raw correlation is not centered at the equilibrium point. This is the harmonic-well approximation. If a full first-order expansion in $T_{\mathrm{eff}}$ is wanted, the quadratic term in the drift is $\omega x^2/2$. Local stationarity gives $⟨x⟩=\omega T_{\mathrm{eff}}/(2{\kappa }^2)+O(T_{\mathrm{eff}}^2)$, so the raw correlation gains a constant $\alpha \omega T_{\mathrm{eff}}/{\kappa }^2$ at that order. The connected correlation is unchanged to first order; the boxed formula is the conventional linearized result rather than a complete nonlinear first-order raw-correlation expansion.

The required time window is long compared with ${\kappa }^{-1}$ but short compared with the escape time, with the lag also short compared with escape. Fluctuations of size $\sqrt{T_{\mathrm{eff}}/\kappa }$ must be small compared with the distance to a neighboring maximum, and the relevant barriers must greatly exceed $T_{\mathrm{eff}}$. For a fixed nonzero noise strength, the unwrapped phase eventually makes phase slips; its raw correlation does not have this stationary infinite-time limit. Thus “large time” here refers to local relaxation within the occupied well, not to the limit after arbitrarily many barrier crossings.

For thermally activated phase slips, the forward saddle is at $\pi -\alpha$ and the backward saddle at $-\pi -\alpha$. Their barriers above the same minimum areIn the Kramers escape rate approximation, both saddles have curvature $-\kappa$ and the minimum has curvature $+\kappa$, so the prefactors agree:This is the forward-backward bias of phase slips. It favors forward motion because each forward step lowers the tilted potential by $2\pi \omega$. For $\omega \ll ϵ$, $\Delta V_{\pm }=2ϵ\mp \pi \omega +O({\omega }^2/ϵ)$, which gives the same ratio; the exact barrier difference is independent of $ϵ$. These probabilities may be read as rates in a common short observation interval, or as normalized competing exit probabilities. The $T_{\mathrm{eff}}$ in the exponent is the noise scale as defined in the Langevin equation; no additional $k_B$ factor is inserted.

The coupled quantum rotor chain has a compact phase and conjugate integer number $n_j={\ell }_j/\hslash$. Expanding the cosine givesIn the phase-space Euclidean path integral, using inverse-energy time, each rotor has the action term ${\ell }_j^2/(2I)-i{\ell }_j{\dot{\varphi }}_j/\hslash$. Completing its Gaussian square produces $I{\dot{\varphi }}_j^2/(2{\hslash }^2)$. Equivalently, apply the short-time free-rotor kernel locally and retain winding sectors. The smooth-sector continuum limit at unit spacing is consequentlyThe overall constant is omitted when the reference bond energy is subtracted. With the same physical J as the Hamiltonian, the printed spatial coefficient has an extra factor $1/{\hslash }^2$. The imaginary-time normalization for a coupled rotor chain requires ${\hslash }^{2}J/2$ inside a bracket with common prefactor $1/{\hslash }^2$. Setting $\hslash =1$ hides the discrepancy; it cannot establish the printed formula with explicit physical units.

Returning to physical time gives $I{\partial }_t^2\varphi =J{\partial }_x^2\varphi$. The harmonic phase mode of a quantum rotor chain therefore hasThis is a gapless sound-like collective phase mode, analogous to neutral superfluid sound. In rotor-chain phase-number hydrodynamics, number fluctuations supply the compressibility and the intersite coupling supplies phase stiffness. The mode is a phase-density oscillation, not a free single-rotor transition. The global phase zero mode of a finite rotor chain is a compact rotor with inertia NI, which should not be counted as an independent zero-frequency oscillator.

The compactness qualification of a harmonic rotor chain matters here. The smooth harmonic theory requires small neighbouring phase differences and cannot settle the entire strong-fluctuation phase diagram. In one dimension its zero-temperature phase variance grows logarithmically, so the gapless regime has algebraic phase correlations rather than a true condensate. Thermal fluctuations and compact phase slips require additional care; phase-slip proliferation can produce a gapped phase outside the Gaussian regime. ^[ref]

A model retaining an oscillation's phase while neglecting amplitude dynamics. Interactions can synchronize relative phases through phase locking; noise can cause phase slips. The Adler phase equation is a simple example for a coupled relative phase.