Codex @codex  0

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.

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.

For the second-harmonic feedback in a long-wave convection amplitude equation, at $\mu =2$ the linear steady operator is $L_0=-(1+{\partial }_x^2{)}^2$. Its eigenvalue on $e^{inx}$ is $-(1-n^2{)}^2$, so the critical Fourier modes are $n=\pm 1$. More generally the linear dispersion relation is $\lambda (k)=-1+\mu k^2-k^4$, giving first onset at $\mu =2,k=1$.

At order ${ϵ}^2$, the weakly nonlinear expansion givesWith ${\Theta }_1=A{e}^{ix}+\overline{A}{e}^{-ix}$,The constant part disappears under differentiation. Inverting $L_0$ on the second harmonic, whose eigenvalue is $-9$, givesup to a critical-harmonic correction absorbed into the definition of $A$. At order ${ϵ}^3$,The coefficient of $e^{ix}$ in the last three terms is respectivelyThe Fredholm solvability condition requires their sum to vanish, because $L_0$ annihilates the critical Fourier mode. For a nonzero amplitude,This requires ${\mu }_2/(3-8s^2/9)>0$. The cubic amplitude coefficient is positive for $s^2<27/8$, yielding a small-amplitude branch in a supercritical bifurcation; it is negative for $s^2>27/8$, yielding a subcritical bifurcation branch in this leading approximation. At $s^2=27/8$, the cubic coefficient vanishes and the quoted relation is singular: higher-order nonlinear terms and a different detuning balance are needed. One must not assert a finite $p$ there.

If a term $ϵ{\mu }_1$ were included in $\mu$, order ${ϵ}^2$ would also contain $-{\mu }_1{\Theta }_{1xx}={\mu }_1{\Theta }_1$. The quadratic nonlinearity produces only the zeroth and second harmonics, with the zeroth differentiated away, so it cannot balance a first-harmonic contribution at that order. Solvability would give ${\mu }_{1}A=0$. Hence a nontrivial critical-mode expansion forces ${\mu }_1=0$ and first balances detuning against cubic amplitude effects at order ${ϵ}^3$.

For the rotating-frame normalization of resonant amplitude forcing, use a rotating phase, $A=B\exp [i(\tilde{\Omega }T+\delta )]$. The real $\delta$ must occur inside the factor of $i$; the printed change of variables omits that $i$ on the phase constant. The correctly phased substitution givesA negative real forcing coefficient can be made positive by changing the forcing phase. Assuming $c>0$, setAll three terms then have the same coefficient scale, givingRenaming $\mathcal{T}$ as $T$ gives the requested canonical equation. This rescaling preserves forward time. If $c<0$, a forward-time normalization instead leaves a positive cubic term, and if $c=0$ a cubic normalization is impossible. If $e=0$, the unforced Landau amplitude equation must be treated separately. Thus the printed canonical form implicitly concerns saturation in a supercritical bifurcation with nonzero forcing.

Rotating Rayleigh-Bénard convection combines buoyancy-driven instability with the Coriolis force. Consider a plane layer of depth $d$ rotating uniformly about the vertical axis, heated from below. Adopt the Boussinesq approximation, fixed boundary temperatures and, for explicit formulas, impermeable stress-free boundary conditions. The conductive state of Rayleigh-Bénard convection is motionless with a linear temperature profile. The dimensionless controls are the Rayleigh number, Prandtl number and Taylor number:Here $\nu$ is kinematic viscosity and $\kappa$ thermal diffusivity. In thermal-diffusion time units, linear perturbations satisfyRotation does no direct mechanical work, since $\mathbf{u}\cdot ({\mathbf{e}}_z\times \mathbf{u})=0$, but couples vertical motion to vertical vorticity and changes the damping and oscillation balance.

For horizontal wavenumber $k$ and vertical mode $n$, put $m=n^2{\pi }^2$ and $a=k^2+m$. Use $w=W\sin (n\pi z)$, $\theta =\Theta \sin (n\pi z)$ and vertical vorticity $\zeta =Z\cos (n\pi z)$ with growth rate $\sigma$. Curling the momentum equation and eliminating pressure givesTheir determinant, without dividing by a possibly zero factor, is the rotating-convection growth-rate polynomialThis makes the linear stability analysis question precise: onset occurs when a root reaches zero real part and all other modes still decay.

A stationary neutral root has $\sigma =0$, givingRotation raises this stationary threshold. An oscillatory neutral root has $\sigma =i\varpi$ with $\varpi \ne 0$. Real-imaginary separation gives the oscillatory neutral curve of rotating convectionThis branch is admissible only when ${\varpi }^2>0$, requiring $P<1$ and sufficiently strong rotation. The restoring Coriolis force coupling permits an inertial/thermal oscillation whose phase-lagged buoyancy can overcome dissipation. At large $P$, temperature and momentum diffusion do not permit that overstability mechanism at primary onset, so the exchange of stabilities is stationary. The actual threshold is the minimum of the stationary and admissible oscillatory curves over all allowed modes, not an arbitrary formal value of $R_o$. Rigid plates require a different vertical eigenproblem and Ekman layers, so the explicit free-slip numbers are not universal.

At large Taylor number, the first vertical mode is selected in the ideal plane layer. Let $t=k^2,m={\pi }^2$. Minimizing the stationary curve givesThus the stationary neutral curve of rotating convection hasThe physical horizontal wavelength is $2\pi d/k_c$, hence decreases as ${\mathrm{Ta}}^{-1/6}$; its prefactor depends on the boundary convention. Thin nearly vertical cells reconcile the strong Coriolis force constraint with viscous and thermal diffusion. The oscillatory minimization replaces the right side of the wavenumber equation by $P^2\mathrm{Ta}m/(1+P{)}^2$, giving the same Taylor number exponent at fixed positive $P$. Where its frequency remains admissible,Equality is $8P^4-P-1=0$, whose positive root is approximately $P_{\ast }=0.6766$. Accordingly, for sufficiently rapid rotation in this free-slip problem, $P<P_{\ast }$ selects oscillatory onset and $P>P_{\ast }$ stationary onset. The weaker condition $P<1$ is only necessary for an oscillatory neutral mode; it does not by itself identify the first instability. Finite Taylor number, finite lateral geometry, allowed discrete wave numbers and plate conditions change the selection.

For the counterpropagating Hopf amplitudes in rotating convection near a simple oscillatory onset, the Hopf bifurcation produces slow complex amplitudes for counterpropagating roll waves. After separating the fast carrier oscillation, symmetry permits the cubic equationswith generally complex coefficients; an $i{\omega }_0{Z}_{\pm }$ term restores the fast frequency if desired. The real parts govern amplitude saturation and the imaginary parts give nonlinear frequency shifts. Write $a_s=\mathrm{Re}{g}_s$, $a_c=\mathrm{Re}{g}_c$. For a travelling wave from a supercritical bifurcation with only one amplitude nonzero, $|Z_{\pm }{|}^2=r/a_s$ with $r,a_s>0$, and the competing wave's growth rate is $r(1-a_c/a_s)$. It is amplitude-stable against that competitor when $a_c>a_s$. A standing wave has equal intensities $r/(a_s+a_c)$; provided this is positive, its intensity-difference mode is stable when $a_s>a_c$. Temporal and spatial phase symmetries leave neutral phase directions, so these are orbital/amplitude stability statements, not decay of every phase displacement.

These coefficients follow from nonlinear interactions and the Fredholm solvability condition obtained by projection onto an adjoint eigenfunction; symmetry alone cannot decide their signs. A negative saturating coefficient gives subcritical bifurcation behavior requiring higher-order terms. Spatial modulation leads to coupled complex Ginzburg–Landau equations with group velocities and diffusion; phase instabilities, mean-flow coupling and differently oriented rolls can destabilize a wave stable in the restricted two-amplitude system. A weakly nonlinear expansion of oscillations therefore predicts travelling waves or standing waves, frequency shifts, modulation and possible secondary mode competition, not a unique universal periodic state.

The Küppers–Lortz instability is a different route to time dependence: it destabilizes steady saturated rolls against oblique roll perturbations. For stationary-roll amplitudes of orientations ${\varphi }_j$, a leading competition system hasA pure roll has $|A_i{|}^2=r/g_0$ with $g_0>0$. An infinitesimal new roll at relative angle $\vartheta$ grows atFor sufficiently strong rotation in appropriate boundary and Prandtl number regimes, some finite oblique angle has $g(\vartheta )<g_0$, so a steady roll is unstable arbitrarily close above its stationary onset. Rotation is handed and allows $g(\vartheta )\ne g(-\vartheta )$, so replacement of one roll by another can favor a definite cyclic sense. Three or more competing orientations can form a heteroclinic cycle; whether it attracts depends on contraction/expansion rates and other modes. Noise, spatially varying domains and modulation can turn this competition into repeated orientation switching and irregular patterns.

The invading rolls are three-dimensional disturbances even when the original straight roll is described by a two-dimensional section. The finite-angle Küppers–Lortz instability mechanism should also be distinguished from the small-angle instability of rotating convection rolls mediated by large-scale mean flow at finite Prandtl number. Numerical thresholds and favored angles depend on mechanical boundaries and material parameters; the essential criterion is the cross-coupling relative to self-saturation. Rotation both changes primary onset and wavelength, and can prevent the resulting steady roll pattern from remaining a stable nonlinear state.

For the real logistic map, $f(x)=x$ factors asThere is always the fixed point $x_0=0$. If $\mu \ne 0$ there is also $x_{\ast }=1-1/\mu$, distinct from zero unless $\mu =1$. At $\mu =0$ the map is constant and zero is its only fixed point. Negative $\mu$ must be retained because the paper specifies $x\in \mathbb{R}$.

Since $f^{\prime }(x)=\mu (1-2x)$, the fixed-point multipliers areThe fixed point stability for an iteration criterion gives $x_0$ attracting for $-1<\mu <1$, and $x_{\ast }$ attracting for $1<\mu <3$. The respective fixed points are repelling when $|\mu |>1$ and $|2-\mu |>1$. At the unit-modulus boundaries the linear criterion needs a nonlinear supplement.

At $\mu =1$, the two fixed point branches meet and exchange stability: this is a transcritical bifurcation. Explicitly, with $\delta =\mu -1$,which has the required nonzero parameter and quadratic terms. At the critical value itself, $f(x)=x-x^2$ attracts small positive points and pushes negative points away, so the merged fixed point is not stable from both sides.

At $\mu =-1$ the origin’s fixed-point multiplier crosses $-1$, while at $\mu =3$ the nonzero fixed-point multiplier crosses $-1$. Both are period-doubling bifurcations. To check existence and stability of the new periodic orbit, factor the second iteration of a map:The two non-fixed roots areThey form a real period-two orbit for $\mu <-1$ or $\mu >3$. Its periodic-orbit multiplier isIt is less than one in absolute value just below $-1$ and just above $3$, respectively. Thus these are supercritical bifurcations in the direction where the parent fixed point loses stability. At the critical values the fixed points themselves are weakly attracting: at $\mu =-1$, $f^2(x)=x-2x^3+x^4$; at $\mu =3$, setting $u=x-2/3$ gives $u\mapsto u-18u^3-27u^4$ after two steps. Both shrink sufficiently small nonzero absolute values. The bifurcation values for finite fixed points are thereforeThe disappearance of the formula for $x_{\ast }$ at $\mu =0$ corresponds to that branch going to infinity, not a local bifurcation at a finite fixed point.

Fix zero mean and period $P$ in the Kuramoto-Sivashinsky equation. At $P=2\pi$ the first Fourier mode becomes neutral and all higher modes are damped. For $k=2\pi /P$, expand $u=a\sin (kx)+b\sin (2kx)+\cdots$. The first equations are $\dot{a}={\lambda }_{1}a+kab/2+\cdots$, $\dot{b}={\lambda }_{2}b-k{a}^2/2+\cdots$, with ${\lambda }_j=(jk{)}^2-(jk{)}^4$. Slaving the second harmonic gives $\dot{a}={\lambda }_{1}a+k^2{a}^3/(4{\lambda }_2)+\cdots$. Since ${\lambda }_2=-12$ at onset, this is a supercritical bifurcation. The small nonzero branch is stable in the sense of orbital stability against disturbances of the same fixed period, modulo translation; stability against disturbances of larger period does not follow from this local argument.