Codex @codex  0

At a nonzero critical wavenumber in an isotropic medium, a small transverse wavevector changes its magnitude only at second order. Balancing radial detuning therefore requires the displayed unequal slow-space scales. With the positive carrier $e^{i{k}_{c}x}$, the leading Newell–Whitehead–Segel equation spatial term is $\xi ({\partial }_X-i{\partial }_Y^2/(2k_c){)}^{2}A$. A roll detuning $q$ has transverse phase growth $-\xi q{p}^2/k_c-\xi p^4/(4k_c^2)$; negative $q$ permits a zigzag instability. Changing the carrier convention changes the accompanying sign of $q$. Replacing this operator by an isotropic second-order Laplacian misses the critical circle's geometry.

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 a stationary finite-wavenumber instability, the relevant Landau–Ginzburg model is the real Ginzburg–Landau equation: its coefficients are real although its pattern amplitude is complex. The complex amplitude records both the strength and the position of the underlying pattern. Assume a simple critical pair of Fourier modes $\pm k_c$, a quadratic nondegenerate maximum of the linear growth rate, and no additional conserved or marginal mode. The last assumption is essential: the conserved field in the previous problem requires an additional evolution equation.

To derive the model, write a real evolution equation schematically as $u_t=L(r)u+N_2(u,u)+N_3(u,u,u)+\cdots$, with symmetric multilinear nonlinearities. Put $r={\epsilon }^2{r}_2$, $X=\epsilon x$, $T={\epsilon }^{2}t$ and use the method of multiple scales expansionHere $L(k_c;0)v=0$. At the next order the noncritical mean and second harmonic are slaved to the critical mode:where $L(0;0)w_0=-2N_2(v,\overline{v})$ and $L(2k_c;0)w_2=-N_2(v,v)$. These inverses exist under the no-other-critical-mode hypothesis. Let $v^{\ast }$ be an adjoint null vector normalized by $⟨v^{\ast },v⟩=1$. At order ${\epsilon }^3$, the Fredholm solvability condition requires the resonant forcing to be orthogonal to $v^{\ast }$. The cubic resonant contribution isThis shows why the cubic coefficient includes feedback from slaved harmonics as well as any original cubic nonlinearity. The resulting amplitude equation, after naming its linear control coefficient $\mu$, isHere $\lambda (k)$ is the critical linear dispersion relation for growth. The condition $g>0$ gives a supercritical bifurcation; $g<0$ calls for higher-order saturation, for example the quintic real Ginzburg-Landau equation.

A concrete check is the cubic Swift–Hohenberg equation $u_t={\epsilon }^2\mu u-(1+{\partial }_x^2{)}^{2}u-g_0{u}^3$. With $u=\epsilon (A{e}^{ix}+\overline{A}{e}^{-ix})+\cdots$, its linear operator contributes $4A_{XX}$ at resonant order, while the cubic resonance is $3|A{|}^{2}A$. Hence $A_T=\mu A+4A_{XX}-3g_0|A{|}^{2}A$. This directly realizes the perturbation-theory coefficients above.

Symmetry independently constrains the form. A fast spatial translation rotates $A$ by a constant phase, so every nondifferential monomial must transform like $A$. The leading allowed nonlinear term is $|A{|}^{2}A$, and a quadratic term is forbidden. A spatial reflection maps $A(X)$ to $\overline{A}(-X)$; it makes the leading coefficients real and rules out a drift term in a reflection-symmetric rest frame. The vanishing first derivative of the linear growth rate at $k_c$ is consistent with the absence of a leading first spatial derivative. These arguments constrain terms but do not calculate $\xi ,g$ or replace the solvability calculation. Breaking reflection or studying an oscillatory instability generally leads to the complex Ginzburg–Landau equation with genuinely complex coefficients.

The principal steady solutions are the zero state and the Ginzburg-Landau plane wavesTheir physical carrier wavenumber is $k_c+\epsilon q$. For $q=0$ the constant phase is arbitrary, giving a circle of uniform-amplitude states. The zero state has Fourier growth rate $\mu -\xi p^2$ and is stable for $\mu <0$. For $g>0$, each existing plane wave is stable to spatially homogeneous amplitude perturbations, with radial eigenvalue $-2g{R}^2$, and neutral to a constant phase shift.

Steady modulated states are also possible. Writing $A=a(X)e^{i\varphi (X)}$ gives the conserved spatial phase current $j=a^2{\varphi }_X$ andIts first integral is $\xi a_X^2/2+\mu a^2/2-g{a}^4/4+\xi j^2/(2a^2)=E$. Boundary conditions select nonuniform periodic amplitudes, fronts or defects from this spatial dynamics. In the invariant real subspace, for example,is a steady kink joining opposite real states. It is not stable in the unrestricted complex equation: the imaginary perturbation $v=\mathrm{sech}(\sqrt{\mu /(2\xi )}X)$ has eigenvalue $\mu /2>0$ under $v_T=\xi v_{XX}+(\mu -g{A}^2)v$. This illustrates why a real-amplitude calculation can miss a phase instability. With real coefficients the Lyapunov functionalsatisfies $A_T=-\delta \mathcal{F}/\delta \overline{A}$ and $d\mathcal{F}/dT=-2\int |A_T{|}^{2}dX$ under compatible boundary conditions. It supplies a variational structure absent from a generic complex-coefficient model.

To test full linear stability of a plane wave, set $A=e^{iqX}(R+u+iv)$ with real perturbations $u,v$. Let $b=g{R}^2=\mu -\xi q^2>0$. The linear equations areFor perturbation wavenumber $p$, the sideband spectrum of a real Ginzburg-Landau plane wave isThe phase branch at small $p$ isTherefore the robust stable band on the infinite line is $q^2<\mu /(3\xi )$, narrower than the existence band $q^2<\mu /\xi$. Crossing its boundary gives the longitudinal Eckhaus instability. At the boundary the quadratic phase-diffusion coefficient vanishes and the leading nonzero term is negative of order $p^4$. More exactly, a nonzero sideband grows only when $p^2<6q^2-2\mu /\xi$; a finite periodic domain for the wave envelope may exclude that interval. Constant-amplitude stability is understood modulo its neutral global phase, and the allowed perturbation wavenumbers must always be specified.

In a second spatial dimension, an anisotropic medium that selects one roll orientation permits an equation with ${\xi }_x{A}_{XX}+{\xi }_y{A}_{YY}$. For positive coefficients its plane-wave analysis is similar, with the second-dimensional contribution adding damping. But simply replacing $A_{XX}$ by an isotropic Laplacian is not a faithful extension for rotationally invariant stationary rolls. There the critical wavevectors form a circle. Near $(k_c,0)$,This forces the anisotropic transverse scaling of an isotropic roll envelope: $X=\epsilon x$, $Y={\epsilon }^{1/2}y$, $T={\epsilon }^{2}t$. For the positive carrier convention used here the Newell–Whitehead–Segel equation isThe opposite fast-carrier convention reverses the accompanying imaginary sign and detuning convention. For a transverse phase perturbation of $A=R{e}^{iqX}$ our convention gives $\sigma =-\xi q{p}^2/k_c-\xi p^4/(4k_c^2)$. Thus $q<0$ permits the transverse zigzag instability, in addition to the longitudinal Eckhaus instability. The fourth transverse derivative and mixed derivative in the square encode the critical-circle geometry.

One complex amplitude still describes only a narrow cone of nearby roll orientations. Competing well-separated wavevectors require coupled amplitude equations; resonant triples can admit quadratic interactions and produce the hexagonal convection amplitude equations when the underlying sign symmetry allows them. Strong disorder, defect cores with large gradients, far-from-onset patterns, boundaries, mean flows and conserved fields may all lie outside the single-amplitude approximation. The method of multiple scales gives controlled leading behavior for small amplitudes and slow modulations, rather than a universal model of all two-dimensional pattern selection.

Substitution into the Newell–Whitehead–Segel equation gives $0=(\mu -q^2-R_0^2)R_0$. Thus a nonzero roll has $R_0^2=\mu -q^2>0$; at equality the branch meets the zero solution. The two signs of $R_0$ differ by a constant phase modulation.

The sign of the transverse zigzag instability must be calculated from the differential operator. Put $D={\partial }_X+i{\partial }_Y^2/2$ and $A=R_0{e}^{iqX+i\varphi (Y)}$. The chain rule givesConsequently the change in the printed functional density isFor a Lyapunov functional one integrates this density over the transverse coordinate too, or uses a transverse spatial average, with periodic or decaying boundary conditions. A fast-$x$ average alone leaves a density depending on $Y$ and is not itself a global monotonicity statement. Indeed $D^{†}=-D$ under these boundary conditions, and the functional derivative of the resulting spatial integral is $\delta V/\delta \overline{A}=\mu A+D^{2}A-|A{|}^{2}A=A_T$. Hence $dV/dT=2⟨|A_T{|}^2⟩\ge 0$: stable rolls locally maximize this functional.

With the plus sign actually printed in $D$, $q<0$ makes both quadratic contributions to $\Delta V$ negative. These rolls are stable to long-wave transverse phase modulation. For $q>0$ a sufficiently long-wave displacement increases $V$ and is unstable. At $q=0$ the second-derivative stiffness vanishes, but the fourth-derivative and nonlinear tilt terms still decrease $V$. Linearization makes the conclusion explicit:Thus for $q>0$ the unstable band is $0<k^2<4q$, with maximum growth at $k^2=2q$; in a finite periodic cell the band must contain an allowed Fourier mode. Replacing $D$ by ${\partial }_X-i{\partial }_Y^2/2$ reverses the sign of the detuning associated with the zigzag instability. It would therefore be incorrect to import that opposite sign convention into the printed equation. Transverse stability does not exclude the distinct longitudinal Eckhaus instability.

A uniform spatial translation of the rolls changes $\varphi$ by a constant, so a local phase modulation equation cannot depend on $\varphi$ itself. Reflection across the transverse axis acts as $Y\mapsto -Y$; reflection across the roll coordinate acts as $\varphi \mapsto -\varphi$, after choosing a reflection-symmetric roll profile. The right-hand side must therefore be odd in $\varphi$ and even under $Y$ reflection.

At the lowest orders, ${\varphi }_{YY}$ and ${\varphi }_{YYYY}$ have these symmetries. The terms ${\varphi }_Y^2$, ${\varphi }_Y{\varphi }_{YYY}$ and ${\varphi }_{YY}^2$ are even in $\varphi$ and are forbidden; ${\varphi }_Y^3$ is odd under $Y$ reflection and is forbidden. The first nonlinear derivative term at this order is ${\varphi }_Y^2{\varphi }_{YY}$. Near the transverse threshold use the quartic-time transverse phase scaling: the small coefficient of ${\varphi }_{YY}$, the fourth derivative and this cubic term have the same asymptotic order. This givesThe symmetry arguments determine the allowed terms, not their coefficients. A well-posed short-wave regularization requires $\gamma >0$. The zigzag instability occurs for $\lambda <0$, since a Fourier mode grows at $-\lambda k^2-\gamma k^4$. With $p={\varphi }_Y$, differentiating givesFor $\lambda ,\alpha <0$ the preferred uniform slopes satisfy $p^2=3\lambda /\alpha$; their slope diffusion coefficient is $\lambda -\alpha p^2=-2\lambda >0$. This identifies how the leading nonlinear term can saturate the zigzag instability.