Codex @codex  0

A codimension-two bifurcation requires two independent parameter conditions to hold simultaneously. Simultaneous stationary and Hopf bifurcation thresholds require both kinds of critical mode in the reduced amplitude equation. A crossing of modes at different wavenumbers is different from a zero-frequency double-zero degeneration of one fixed-wave-number cubic.

Near an oscillatory onset of rotating Rayleigh-Bénard convection, opposite travelling roll waves have coupled Hopf bifurcation amplitudes. Symmetry gives the displayed cubic amplitude equations, with generally complex self-coupling and cross-coupling coefficients. Write $a_s=\mathrm{Re}{g}_s$ and $a_c=\mathrm{Re}{g}_c$. For $r>0,a_s>0$, a single travelling wave has intensity $r/a_s$ and is stable against its counterpropagating competitor when $a_c>a_s$. Equal standing-wave intensities are $r/(a_s+a_c)$; with $r>0$ and positive denominator they are stable to intensity differences when $a_s>a_c$. Phase symmetries give orbital, rather than strict phase-decay, stability. Spatial modulation, mean flows and oblique convection rolls require additional stability tests.

A forcing at the carrier wavenumber breaks continuous spatial translation symmetry and adds a constant to its complex amplitude equation. In a frame following a moving forcing pattern, its phase drift gives $i\Lambda A$. A nonzero forcing coefficient is normalized to one by scaling amplitude by its cube root and time by its two-thirds power. A steady intensity $s=|A{|}^2$ satisfies $s[(\mu -s{)}^2+{\Lambda }^2]=1$, and its phase is uniquely determined by the steady complex equation.

For a self-adjoint operator $L$, solving $Lu=f$ requires $f$ to be orthogonal to every $v\in \ker L$, because $⟨v,Lu⟩=⟨Lv,u⟩=0$. When the range is closed, this condition is also sufficient: $(\mathrm{ran}L{)}^{\perp }=\ker L$ and closedness gives $\mathrm{ran}L=(\ker L{)}^{\perp }$. In finite dimensions the range is automatically closed. This Fredholm solvability condition determines amplitude equations by projecting perturbative forcing onto critical eigenfunctions.

A finite-wavenumber convection instability with a resonant three-wave triad can yield three coupled amplitude equations. For the gradient nonlinearity $\nabla \cdot (|\nabla \Theta {|}^2\nabla \Theta )$ and weak quadratic asymmetry $-ϵ\gamma \nabla \cdot (\Theta \nabla \Theta )$, projection onto a triad of length squared $K$ gives the stated equation and its cyclic copies. The quadratic coefficient is $+\gamma K$, while both self and cross cubic coefficients are $-3K^2$ for angles of $120$ degrees. The finite critical mode requires positive linear damping in the underlying long-wave equation.

The Kuramoto-Sivashinsky equation $u_t+u{u}_x+u_{xx}+u_{xxxx}=0$ combines destabilizing long-wave growth, stabilizing fourth-order diffusion, and nonlinear transport. Its potential form is ${\Phi }_t=-{\Phi }_{xx}-{\Phi }_{xxxx}+({\Phi }_x{)}^2$, with $u=-2{\Phi }_x$. It arises as a long-wave phase amplitude equation near loss of the Benjamin-Feir stability condition. On a domain of period $P$, the zero solution's Fourier modes grow at rate $k^2-k^4$, $k=2\pi n/P$.

For a real steady amplitude and zero-mean field, twice integrating the field equation yields $B=\kappa (⟨A^2⟩-A^2)$, $\kappa =\mu /\sigma$. The amplitude equation becomes $A^{\prime \prime }+(1-\kappa ⟨A^2⟩)A+(\kappa -1)A^3=0$. Using $(R\mathrm{sech}cX{)}^{\prime \prime }=c^{2}A-2c^2{A}^3/R^2$ gives $\kappa ⟨A^2⟩=1+c^2$ and $R^2=2c^2/(\kappa -1)$, requiring $\kappa >1$. Its finite-window mean is $2R^2\tanh (cL/2)/(cL)$. For $cL\gg 1$, consistency gives $c+c^{-1}=4\kappa /[L(\kappa -1)]$. The leading roots coalesce at $c=1$; for $L>2$ the leading upper parameter limit is $\kappa \le L/(L-2)$, while retaining the finite-window factor makes the necessary inequality strict. The isolated sech profile has an exponentially small derivative mismatch on a periodic interval, so exact periodic tails need correction. For $L\le 2$ the algebra supplies no finite upper bound of this form.

Magnetoconvection is thermal convection coupled to the evolution and forces of a magnetic field. Conserved magnetic flux can supply a slow field that must be retained along with the convective complex amplitude, rather than being eliminated into a single local amplitude equation.

For ${\psi }_1=\cos (\pi x)\sin (\pi z)$, the first temperature mode is ${\theta }_1=(2\pi {)}^{-1}\sin (\pi x)\sin (\pi z)$. The streamfunction advection bracket is $J({\psi }_1,{\theta }_1)=(\pi /4)\sin (2\pi z)$, independent of $x$. The second-order equations therefore have ${\psi }_2=0$ andThis modifies the vertical temperature gradient. Its interaction with the critical roll produces the cubic term in the amplitude equation.

The Newell–Whitehead–Segel equation is an amplitude equation for slowly modulated stationary rolls in an isotropic pattern-forming medium. With the displayed sign convention, a roll $A=R{e}^{iqX}$ has $R^2=\mu -q^2$. A transverse phase modulation has linear growth rate $q{k}^2-k^4/4$, so long-wave transverse instability occurs for $q>0$. Reversing the carrier-wave convention reverses the sign attached to its detuning $q$; one must derive the stability sign from the operator actually used.

The spatial mean is conserved: averaging every term in the partial differential equation over the periodic cell gives ${\partial }_t⟨T⟩=0$. Adding a constant to $T$ also preserves the equation. We therefore fix the mean, conveniently at zero. This qualification is essential: without it the constant Fourier mode contributes an additional neutral direction, so the centre subspace would have five real dimensions rather than four.

On a Fourier mode $e^{i(mx+ny)}$, the linearized linear operator has eigenvalueOn the fixed-mean subspace, the first eigenvalues to reach zero are those with $m^2+n^2=1$. At $\mu =1$ these are the modes $e^{\pm ix},e^{\pm iy}$; every other nonconstant Fourier mode is strictly damped. The reality condition leaves two complex amplitudes, hence a four-real-dimensional centre manifold on the fixed-mean subspace. The trivial state loses linear stability as $\mu$ increases through one.

Use slow time $\tau ={\epsilon }^{2}t$, and write $T=\epsilon W+{\epsilon }^3{T}_3+\cdots$, whereThere is no quadratic nonlinearity, so no forced second-order correction is needed. At order ${\epsilon }^3$, with $\mu =1+{\epsilon }^2\nu$, projection onto the critical Fourier modes givesThe correction $T_3$ does not contribute to the critical projection because the linear operator vanishes on those modes. To find the nonlinear coefficient of $e^{ix}$, split the divergence term intoThe first term contributes $-3|A{|}^{2}A$. In the second, only the zero-$y$ Fourier coefficient of $W_y^2$ contributes; it is $2|B{|}^2$, and differentiating $W_x$ once more gives $-2|B{|}^{2}A$. The third term has zero $y$-average, and the fourth has no $x$-dependent critical component. Interchanging $x,y$ gives the other amplitude equation. Thus the conserved-mean convection amplitude equations areTo establish the selected planform, put $A=r{e}^{i\theta }$, $B=s{e}^{i\psi }$. The phases are constant and the magnitudes obeyFor $\nu >0$, the square state has $r=s=\sqrt{\nu /5}$. Its magnitude Jacobian matrix is $-(2\nu /5)\left(\begin{array}{cc}3& 2\\ 2& 3\end{array}\right)$, with eigenvalues $-2\nu$ and $-2\nu /5$. Its two neutral phases correspond to translation symmetry. The pure roll has $r^2=\nu /3,s=0$; an infinitesimal $B$ disturbance grows at rate $\nu -2\nu /3=\nu /3>0$. The other pure roll is equivalent, while the trivial state is unstable for $\nu >0$. Consequently squares, with $|A{|}^2=|B{|}^2=\nu /5$, are the stable small-amplitude planform, modulo translations. For $\nu <0$ the trivial state is stable.

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.

Near a spatially extended Hopf bifurcation, a critical oscillatory normal mode has a slowly varying complex amplitude. For a selected travelling-wave branch, a method of multiple scales expansion takes the formHere the distance from onset is of order ${\epsilon }^2$, $v_0$ is the critical eigenvector, and $v_g$ is the group velocity. The comoving coordinate removes the first spatial derivative of the envelope. Projection at the first resonant order gives $A_T=\rho A+D_0{A}_{XX}-G_0|A{|}^{2}A$. Spatial translation and the temporal phase of the Hopf oscillation give a constant-phase symmetry $A\mapsto e^{i\theta }A$, which permits the cubic term $|A{|}^{2}A$ but forbids a generic $A^2$ term. The coefficients are complex because the envelope has both growth and frequency detuning. For a supercritical branch with positive real diffusion, real rescalings and a uniform phase rotation give the complex Ginzburg–Landau equationA subcritical Hopf branch requires higher saturation terms instead. If spatial reflection makes opposite travelling waves critical together, two coupled amplitude equations are required initially; the scalar equation applies to a selected single-wave branch when its competing amplitude is stable.

Distant boundaries can still determine which travelling pattern is observed. They can select its phase and wavenumber, reflect a travelling disturbance, inject a competing wave, or let a growing wave packet leave the domain. In particular, amplification in a comoving frame can be convective wave-packet instability rather than absolute wave-packet instability at a fixed laboratory point. A finite-domain outcome then need not reflect intrinsic bulk modulation dynamics. Periodic boundary conditions are the appropriate idealization when studying bulk instabilities without end selection or reflected-wave forcing: they eliminate physical end layers and make the allowed envelope Fourier modes discrete. They do not assert that real distant boundaries are always negligible. Choose a periodic length compatible with the carrier, fix the phase winding, and compare perturbations within that same periodic problem. On an open physical domain the group velocity and actual boundary conditions must be restored.

The optimum carrier wavenumber corresponds to the spatially uniform envelope $A_f=e^{-icT}$, often called a flat state. Its phase is arbitrary. Write $A=e^{-icT}(1+r)e^{i\varphi }$ with real small amplitude and phase disturbances. To first order,A Fourier mode $e^{\lambda T+ikX}$ therefore has matrix and dispersion relationThe trace is negative. The amplitude mode at $k=0$ has eigenvalue $-2$, while the phase mode is neutral because constant phase shifts preserve the solution. Put $D=1+bc$. For every nonzero real $k$, the determinant is positive if $D\ge 0$, whereas $D<0$ makes it negative for $0<k^2<-2D/(1+b^2)$. ThusThis is the Benjamin-Feir stability condition. At $D=0$, nonzero modes still decay, but the leading long-wave decay is fourth order rather than diffusive. On a periodic envelope interval of length $L$, the smallest nonzero wavenumber is $2\pi /L$, so strict stability of all allowed nonzero modes requiresA sufficiently short periodic domain can therefore exclude the unstable band even when the infinite-domain Benjamin-Feir stability condition fails.

The amplitude is damped while the phase is slow. For a slowly varying phase, the leading slaved amplitude is $r=-\frac{1}{2}\{b{\varphi }_{XX}+({\varphi }_X{)}^2\}+\cdots$. Substitution in the phase equation gives the diffusion term $D{\varphi }_{XX}$ and the nonlinear frequency shift $(c-b)({\varphi }_X{)}^2$. The latter also follows directly from a constant phase gradient: an envelope with gradient $q$ has amplitude squared $1-q^2$ and frequency $c+(b-c)q^2$. Constant phase symmetry rules out undifferentiated $\varphi$; translation invariance gives constant coefficients. The leading comoving Ginzburg-Landau equation also has spatial parity, so odd linear derivatives are absent at this order. Its neutral phase eigenvalue expands asThis coefficient follows by inserting a power series for $\lambda$ in the quadratic dispersion relation. Close to $D=0$, it becomes $E=(1+b^2)/2>0$ at leading order. The long-wave phase reduction of the complex Ginzburg-Landau equation is consequentlyThe fourth derivative regularizes negative phase diffusion. Terms omitted here have higher order in the long-wave expansion. First-order advection has been removed by the group-velocity frame; higher odd derivatives in a more general travelling-wave problem are beyond this leading amplitude approximation.

Take $D=-d$ with $d>0$ small. At the instability boundary $c=-1/b$, so $g=-(1+b^2)/b$ is nonzero. Set $\xi =\sqrt{d/E}X$, $\tau =d^{2}T/E$, and $\varphi =(d/g)\Phi$. At leading order the phase dynamics become the potential Kuramoto-Sivashinsky equation,With $u=-2{\Phi }_{\xi }$ this isThe variable $u$ is a phase-gradient disturbance, not the original physical field. Periodicity of the phase fixes ${\int }_0^{P}ud\xi =0$ for the zero-winding branch; the Kuramoto-Sivashinsky equation conserves this mean. The scaled period is $P=L\sqrt{d/E}$.

At fixed period $P$, the flat phase has growth rates ${\lambda }_n=k_n^2-k_n^4$, $k_n=2\pi n/P$. Its first loss of stability occurs at $P=2\pi$. To find what replaces it, put $k=2\pi /P$ and consider $P$ just above $2\pi$. After choosing an origin by translation, expand a small steady modulation as $u=a\sin (k\xi )+b_2\sin (2k\xi )+\cdots$. Projection of $-u{u}_{\xi }$ onto the first two harmonics yieldsSince ${\lambda }_2=-12$ at onset, this harmonic is slaved to $b_2=k{a}^2/(2{\lambda }_2)+\cdots$. The fundamental amplitude obeysThe primary periodic bifurcation of the Kuramoto-Sivashinsky equation is therefore supercritical. For $0<P-2\pi \ll 1$,The radial eigenvalue of this nonzero branch is $-2{\lambda }_1+O({\lambda }_1^2)<0$; all higher harmonics remain damped and the modulation's spatial phase is neutral. Thus the small finite-amplitude modulation is stable in the sense of orbital stability, at fixed mean and fixed period, despite the flat state's instability. The two signs of $a$ represent translations of the same modulation. In the potential formulation the spatial mean of $\Phi$ advances at $⟨({\Phi }_{\xi }{)}^2⟩$: the original wave acquires a nonlinear frequency correction rather than requiring a time-independent absolute phase.

For larger $P$, further linear modes become active and the small-amplitude argument no longer determines stability. A periodic modulation $\bar{u}$ must then be tested with the operator $v_{\tau }=-v_{\xi \xi }-v_{\xi \xi \xi \xi }-{\partial }_{\xi }(\bar{u}v)$ on mean-zero periodic disturbances. Its translation mode is always neutral, and secondary oscillatory or more complicated modulation can occur. Moreover, stability against perturbations of the same period does not prove stability against disturbances of a larger period; those require a separate sideband or Floquet theory calculation. The controlled conclusion near the first threshold is the existence of a stable small modulation family, not universal stability of all periodic Kuramoto-Sivashinsky solutions.

The third planform is fixed by a diagonal reflection and changes sign under a quarter-turn. Extend the preceding group representation byThe other diagonal reflection changes only the sign of $B$. Consequently the first two components of an equivariant dynamical system have factors $A$ and $B$, whereas its third component is even in each of $A,B$. The quarter-turn exchanges $A^2$ and $B^2$ while negating $C$. Its admissible quadratic coupling is therefore $A^2-B^2$ in the third equation and opposite $AC,BC$ couplings in the first two. At cubic order the only new terms are $A{C}^2,B{C}^2,C^3,C(A^2+B^2)$. The linear action on the two-dimensional irreducible mode has one coefficient ${\mu }_1$, while the distinct scalar representation has its own coefficient ${\mu }_2$. A small splitting of the two marginal linear modes supplies these two unfolding parameters. With real coefficients the resulting amplitude equations areThese equations follow by listing all equivariant monomials through degree three; the displayed minus signs are coefficient conventions, not consequences of symmetry.

In the specified reduced system put $A=B=s$, $C=0$. A nonzero branch exists for ${\mu }_1>0$, with $q=s^2={\mu }_1/(1+\lambda )$. Write perturbations as $u=(\delta A+\delta B)/2$, $v=(\delta A-\delta B)/2$, and $c=\delta C$. Linearization givesFor the square-pattern interaction with a sign-changing scalar mode, the remaining eigenvalues have sum $\tau =p+{\mu }_2$ and product $D=p{\mu }_2+4{\alpha }_1{\gamma }_{1}q=p({\mu }_2+K)$, where $K=2{\alpha }_1{\gamma }_1/(\lambda -1)>0$. A real two-dimensional stability matrix has both eigenvalues in the left half-plane exactly when its trace is negative and its determinant is positive. ThereforeIf ${\mu }_2<-K$, the block has negative determinant and one growing direction, explaining the instability at sufficiently small ${\mu }_2$. At the lower endpoint ${\mu }_2=-K$, one eigenvalue is zero and the other is $p-K<0$. The symmetry $(A,B,C)\mapsto (B,A,-C)$ fixes the diagonal state and negates its critical $(v,c)$ direction. Thus a generic nonlinear unfolding has a pitchfork bifurcation to unequal-amplitude states with nonzero $C$.

At the upper endpoint ${\mu }_2=-p$, the trace vanishes andThe symmetric eigenvalue remains negative. The critical conjugate pair crosses the imaginary axis with derivative $d\mathrm{Re}\sigma /d{\mu }_2=1/2$. If its cubic Hopf normal form coefficient is nonzero, this is a Hopf bifurcation and nearby periodic solutions exist on one side. Their stability and the side of existence depend on that coefficient; the linear calculation does not determine them. It does establish the generic possibility of oscillating square patterns. The equality $p=K$ is a double-zero degeneracy, outside this nondegenerate Hopf bifurcation argument.

Represent the leading real pattern by $A(T)e^{i{k}_{c}x}+\overline{A}(T)e^{-i{k}_{c}x}$. Before forcing, continuous translation symmetry acts as $A\mapsto e^{i\theta }A$, and spatial reflection acts as $A\mapsto \overline{A}$. Equivariance excludes quadratic terms and permits a linear term and $|A{|}^{2}A$ at cubic order. Symmetry under reflection makes their coefficients real. A supercritical bifurcation has a saturating cubic coefficient; rescaling amplitude and time makes it $-1$. The imperfection has a Fourier component at the same wavenumber, so its resonant projection contributes an additive real coefficient. On the spatially uniform wave envelope subspace the amplitude equation is therefore $A_T=\mu A-|A{|}^{2}A+\alpha$. Allowing modulations would additionally produce the spatial terms of the real Ginzburg–Landau equation.

For a forcing pattern moving in the positive $x$ direction, its laboratory-frame coefficient is $\alpha e^{-i\Omega T}$, where $\Omega$ is proportional to $k_{c}c$ after converting to slow time. The comoving amplitude is $B=e^{i\Omega T}A$, since the fast carrier becomes $e^{i{k}_c(x-ct)}$. Differentiation givesA constant phase choice makes $\alpha >0$, and complex conjugation reverses the sign of $\Omega$, so take $\Lambda =\Omega >0$. Put $B={\alpha }^{1/3}a$ and $\tau ={\alpha }^{2/3}T$. The rescaled linear coefficients are $\mu /{\alpha }^{2/3}$ and $\Lambda /{\alpha }^{2/3}$, and the forcing becomes one. Below the rescaled variables are renamed $A,T,\mu ,\Lambda$, giving the directly forced cubic amplitude equation.

At a periodic steady state, two spatial integrations of $0=(B+\delta A^2{)}_{XX}$ make $B+\delta A^2$ constant. The zero-mean constraint fixes that constant to $\delta ⟨A^2⟩$. Define $\beta =1-\delta >0$ and $\nu =\mu -\delta ⟨A^2⟩$; the stationary amplitude equation isFor a broad localized plateau, the two transition layers occupy $O(\ell )$ length while the plateau and exterior occupy $2hL$ and $2(1-h)L$. Provided both lengths are large compared with the front width $\ell$, spatial averaging gives $⟨A^2⟩=h{R}^2+O(R^2\ell /L)$.

Multiplying the stationary equation by $A_X$ gives the first integralIn the well-separated-front approximation, $A\to 0$ with $A_X\to 0$ outside, so $E=0$. On the plateau $A\to R\ne 0$, giving both $\nu +\alpha R-\beta R^2=0$ and $W(R)=0$. Substitution of the first into the second gives $R^3(\beta R/4-\alpha /6)=0$, henceThe equal-potential condition is therefore $2{\alpha }^2=-9\nu \beta$. It is a front balance, not an integration of the stationary equation without its $A_X$ multiplier.

These plateau values also produce actual isolated fronts. At this balance $W(A)=-\beta A^2(A-R{)}^2/4$, so for $\alpha >0$ the increasing front satisfies $A_X=\sqrt{\beta /2}A(R-A)$. It isFor $\alpha <0$ the sign-reversed construction gives the negative plateau. Two oppositely oriented well-separated fronts give the leading periodic mesa solution approximation; their weak interaction and tails supply finite-$L$ corrections.

Finally $\nu =\mu -\delta h{R}^2$ gives the coexistence fraction of a conserved-field amplitude mesa:The leading existence range is the interval obtained from $0<h<1$:This assumes $\alpha \ne 0$ and $L\gg \ell$, with $hL,(1-h)L\gg \ell$. At the two displayed ends a transition pair can no longer be treated as well separated, so they are asymptotic bounds rather than exact finite-period bifurcation values. For $\alpha =0$ this particular nonzero plateau balance collapses to $R=0$ and does not yield the mesa solution family.

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.