Codex @codex  0

Use a real physical field, so the three positive wavevectors are accompanied by their complex-conjugate modes. Modulo lattice translations, the effective symmetry group is $G={\mathbb{T}}^2⋊C_6$. Its translation group action isA $120^{\circ }$ rotation cyclically permutes the amplitudes; a half-turn $\kappa$ acts by simultaneous complex conjugation. Together these generate the sixfold rotation group. These formulas specify the group action independently of the choice of orientation of the three wavevectors.

For a nonzero real hexagon $(h,h,h)$, its stabilizer subgroup is the rotation group $C_6$. A translation fixes it only when all three phase factors equal one. Cyclic permutation forces a fixed vector to have equal components, and the half-turn forces those components to be real. Therefore $\mathrm{Fix}(C_6)=\{(a,a,a):a\in \mathbb{R}\}$. For a roll $(r,0,0)$ with nonzero real $r$, the stabilizer subgroup consists of the translations $T_{0,\theta }$ and the half-turn, with $\kappa T_{0,\theta }\kappa =T_{0,-\theta }$. It is $S^1⋊C_2$, isomorphic to $O(2)$ as an abstract group even though spatial reflections are absent. The circle eliminates $B,C$ from its fixed-point subspace of a group action, and the half-turn makes $A$ real. ThusIf unreduced physical translations are used, the full lattice kernel is included in both stabilizer subgroups. The real representation on the six critical components is absolutely irreducible: translations distinguish the three Fourier pairs, rotations permute them, and the half-turn rules out a complex scalar in the commutant. The equivariant branching lemma consequently supplies generic steady branches of both axial types at a transverse crossing of the common linear eigenvalue. It guarantees existence, not stability or exhaustiveness of the possible branches.

Translation equivariance severely restricts the monomials in $\dot{A}$. Through cubic order they are $A$, $\overline{B}\overline{C}$, $|A{|}^{2}A$, $|B{|}^{2}A$ and $|C{|}^{2}A$: each has the same translation weight as $A$, using ${\mathbf{k}}_1+{\mathbf{k}}_2+{\mathbf{k}}_3=0$. The half-turn forces their coefficients to be real; the $120^{\circ }$ rotation gives the two cyclic counterpart equations. A spatial reflection would interchange $B,C$ and require equal cross-couplings, but it is not a symmetry here. After a real amplitude rescaling makes the generic nonzero quadratic coefficient one, the rotational hexagon amplitude equations arewith the other two equations obtained by cyclic permutation. All four coefficients are real; $\delta$ records the absence of spatial reflection. The normalization assumes the resonant quadratic coefficient is nonzero; its vanishing would be an additional degeneracy.

For the roll, $r^2=\mu /{\nu }_1$. The radial eigenvalue is $-2\mu$ and the phase of the active mode is a neutral translation. The transverse variables $(b,\overline{c})$ have matrixIts two real eigenvalues, each occurring twice in the full real linearization, are $-({\nu }_2-{\nu }_1)r^2\pm \sqrt{{\delta }^2{r}^4+r^2}$. Define $d={\nu }_2-{\nu }_1$ and $D_r=d^2-{\delta }^2$. Strict attraction modulo translation is therefore equivalent toIndeed the larger transverse eigenvalue is negative precisely when $d{r}^2>\sqrt{{\delta }^2{r}^4+r^2}$, and squaring is legitimate only after requiring $d>0$. If ${\delta }^2>d^2$, that inequality is impossible. Equality ${\delta }^2=d^2$ is also insufficient because of the additional positive term $r^2$. Any nonzero roll that exists with ${\nu }_1<0$ already has a positive radial eigenvalue; the degenerate case ${\nu }_1=0$ has no ordinary isolated nonzero roll branch.

For hexagons put $N={\nu }_1+2{\nu }_2$. A nonzero real amplitude $h$ satisfies $\mu =N{h}^2-h$, or $h=(1\pm \sqrt{1+4N\mu })/(2N)$ when $N\ne 0$. Separate real and imaginary disturbances. The real Jacobian matrix is cyclic, with diagonal $-h-2{\nu }_1{h}^2$ and off-diagonal entries $h-2({\nu }_2\pm \delta )h^2$. Its eigenvalues areThe imaginary Jacobian matrix has every entry equal to $-h$, giving eigenvalues $-3h,0,0$. The two zero modes translate the hexagonal pattern. Thus the general strict stability conditions for any nonzero root $h$ areThese conditions state orbital stability transverse to translations, not asymptotic convergence to one prescribed spatial phase. Equalities require a nonlinear bifurcation calculation.

When ${\nu }_2>{\nu }_1>0$, both $d,N$ are positive. The upper positive hexagon branch is radially stable after its saddle-node bifurcation at $\mu =-1/(4N)$, $h=1/(2N)$; it remains stable until $h=1/d$. At this second threshold the pair ${\lambda }_{\pm }$ crosses the imaginary axis with nonzero frequency if $\delta \ne 0$. The crossing is transverse, since $dh/d\mu =1/(2Nh-1)$ and the derivative of $-2h+2d{h}^2$ with respect to $h$ is two there. Hence, subject to the usual nonzero cubic coefficient, this is a Hopf bifurcation:For $\delta =0$ the same threshold has two zero real eigenvalues, so it is not a Hopf bifurcation.

A possible phase portrait when neither steady branch is stable has an attracting periodic oscillation surrounding the unstable hexagon, with the three amplitudes waxing and waning cyclically. This is realized, for example, by ${\nu }_1=1$, ${\nu }_2=2$, $\delta =2$ just above ${\mu }_H=4$. Rolls are unstable for every positive $\mu$ because ${\delta }^2=4>d^2=1$. At the hexagon Hopf point $h=1$, the real spectrum is $-9,\pm 4i\sqrt{3}$. Evaluation of the cubic Hopf normal form with a unit-norm critical eigenvector gives $\mathrm{Re}{G}_{21}=-64/27<0$, so this example has a supercritical Hopf bifurcation and an attracting small limit cycle in the real invariant subspace. The figure shows such a cycle at $\mu =4.2$ and an approaching trajectory; the positive octant is forward invariant because an amplitude's derivative at zero is the product of the other two. This is a possible portrait, rather than a claim that every parameter set with unstable rolls and hexagons has the same attractor. Stability asserted for the periodic orbit here is within the stipulated real subspace.

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.

Take the leading real field to be $A{e}^{i(k_{c}x+{\omega }_{0}t)}+B{e}^{i(-k_{c}x+{\omega }_{0}t)}+\text{complex conjugate}$. A spatial translation gives $(A,B)\mapsto (e^{i{k}_c\delta }A,e^{-i{k}_c\delta }B)$; a time translation gives $(A,B)\mapsto (e^{i\psi }A,e^{i\psi }B)$; spatial reflection exchanges $A$ and $B$.

For a monomial $A^j{\overline{A}}^k{B}^l{\overline{B}}^m$ in the first amplitude equation, covariance requires $j-k+l-m=1$ and $j-k-l+m=1$. Hence $j-k=1$, $l-m=0$. At linear order only $A$ is allowed, no quadratic monomial satisfies these conditions, and the cubic possibilities are $A^2\overline{A}$ and $AB\overline{B}$. Reflection fixes the coefficients in the second equation. Removing the small common linear frequency correction by a rotating frame makes $\mu$ real. This derives the two cubic amplitude equations, rather than merely postulating them.

Let $r=|A{|}^2$, $s=|B{|}^2$. Taking real parts givesA travelling wave has $(r,s)=(\mu /{\beta }_r,0)$, or its reflection, whenever $\mu /{\beta }_r>0$. Its phase evolves at $-{\beta }_i\mu /{\beta }_r$, shifting the physical oscillation frequency. The nonzero-amplitude radial eigenvalue is $-2\mu$, while the absent mode has real growth rate $\mu (1-{\gamma }_r/{\beta }_r)$. Therefore, on the supercritical side,The phase direction is neutral because of the symmetry, so stability means orbital stability. A branch with $\mu <0$ is radially unstable, whatever the transverse coefficient.

A standing wave has $r=s=h=\mu /({\beta }_r+{\gamma }_r)>0$. Its two phases evolve together at $-({\beta }_i+{\gamma }_i)h$; the relative phase fixes the spatial position of the standing pattern. In the two intensity directions the eigenvalues are $-2\mu$ and $-2h({\beta }_r-{\gamma }_r)$. ThusAgain the symmetry phases are neutral and the conclusion is orbital stability. Equality ${\beta }_r={\gamma }_r$ leaves a continuum of mixed intensities at cubic order; equality ${\beta }_r+{\gamma }_r=0$ removes cubic standing-wave saturation. These degeneracies require higher-order terms and cannot be classified by the cubic truncation. The zero state is linearly stable for $\mu <0$ and unstable for $\mu >0$.

The natural expectation is orbital stability, since a spatial translation of a travelling wave is another travelling wave. Thus convergence for ordinary perturbations, when it holds, is to a suitably translated front. It also matters which perturbations are allowed in the leading edge: the zero state ahead of this front is linearly unstable, and a sufficiently shallow positive tail can select a different speed. One should not assert asymptotic stability for arbitrary disturbances without a norm and a class of initial data. The requested essential spectrum calculation is a necessary part of such an analysis, not by itself a complete nonlinear stability proof.

In the moving coordinate $z=x-ct$, write $u(x,t)=w_c(z)+v(z,t)$. Expansion of the reaction term givesThe first equation to study is its linearization,and the associated eigenvalue equation $\mathcal{L}\varphi =\lambda \varphi$. Differentiating the front equation gives $\mathcal{L}{w}_c^{\prime }=0$, the neutral translation mode whenever it belongs to the chosen function space. Ahead of the front the limiting coefficient is $k$, so in unweighted L2 space the limiting dispersion relation is $\lambda =-{\xi }^2+ic\xi +k$: its rightmost real part is positive. An exponential weight penalizes precisely these disturbances ahead of the front.

There is a small functional-analytic correction to the printed definition. Use the complete exponentially weighted L2 spaceThe additional requirement of unweighted L2 space membership in the printed set makes that set incomplete under its stated weighted norm. Indeed $v_N=1_{[-N,0]}$ belongs to the intersection and is a weighted Cauchy sequence, but its weighted limit $1_{(-\infty ,0]}$ does not belong to unweighted L2 space. The usual Hilbert space for this spectral argument is the completion $X_{\gamma }$. We use that completion explicitly, rather than treating the intersection as a complete space.

Multiplication $T_{\gamma }v=e^{\gamma z}v$ is a unitary map from $X_{\gamma }$ onto L2 space. Take the domain $T_{\gamma }^{-1}{H}^2(\mathbb{R})$ for $\mathcal{L}$. The transformed closed operator isThis follows from $T_{\gamma }{\partial }_z{T}_{\gamma }^{-1}={\partial }_z-\gamma$. The two limiting constant-coefficient dispersion relations, obtained by the Fourier transform, areIn particular the necessary strict leading-edge condition is ${\gamma }^2-c\gamma +k<0$. For $c>2\sqrt{k}$ this meansA particularly convenient choice isWith this choice there is no first derivative, and ${\mathcal{L}}_{\gamma }={\partial }_z^2+k-c^2/4-2w_c(z)$ is a self-adjoint differential operator on $H^2(\mathbb{R})$. This also avoids possible ambiguities between definitions of essential spectrum for nonnormal operators.

The relevant perturbation theorem is essential-spectrum invariance under relatively compact perturbations: if $A_0$ is a closed operator, $V:D(A_0)\to H$ is compact with respect to the graph norm, and $A_0+V$ is closed on the same domain, their Fredholm essential spectrum is the same. For self-adjoint operators this is the usual essential spectrum after removing isolated finite-multiplicity eigenvalues. To see why the theorem applies, regard $A_0-\lambda$ and $A_0+V-\lambda$ as bounded maps from the graph norm domain to $H$, and use compact perturbation invariance of Fredholm operators for each $\lambda$.

Define a step-potential comparison operatorBoth $A_0$ and ${\mathcal{L}}_{\gamma }$ are self-adjoint on $H^2(\mathbb{R})$. Their difference is multiplication byIt is relatively compact, not necessarily compact as a multiplication map on all of L2 space. To prove the needed compactness, take a sequence bounded in $H^2(\mathbb{R})$. On each bounded interval the Rellich-Kondrachov compactness theorem gives a subsequence converging in L2 space. A diagonal subsequence makes the products by $V$ converge locally. Outside a sufficiently large interval, $\Vert V{\Vert }_{\infty }$ is arbitrarily small and the sequence has a uniform L2 norm bound, so those tails are uniformly small. Thus the products converge globally in L2 space, proving graph-domain compactness.

Finally,Here is a direct verification. Its quadratic form is at most $(k-c^2/4)\Vert h{\Vert }_2^2$, so no spectrum lies above this endpoint. For every $\lambda \le k-c^2/4$, choose $\xi$ with $\lambda =k-c^2/4-{\xi }^2$. Normalized, increasingly long cutoffs of $e^{i\xi z}$, translated far to the right, have disjoint supports and satisfy $\Vert (A_0-\lambda )h_N{\Vert }_2\to 0$: differentiating the cutoff gives errors of orders $N^{-1}$ and $N^{-2}$. They converge weakly to zero. This is a singular sequence in the Weyl criterion, proving that every such $\lambda$ belongs to the essential spectrum. The relatively compact perturbation theorem therefore givesThis proves the weighted essential spectrum of a Fisher travelling front is strictly stable when $c>2\sqrt{k}$. In fact here the quadratic form gives more:Since $w_c\ge 0$, the full transformed spectrum lies below $-(c^2/4-k)$, and the evolution given by its linearization decays at least as $e^{-(c^2/4-k)t}$ in this weighted norm. This does not contradict translation invariance: the nonzero translation mode cannot lie in this weighted operator domain, since otherwise it would be a zero eigenfunction, contradicting the strictly negative quadratic-form bound. The weight restricts the permitted perturbations and excludes these neutral translations. Nor does a weighted L2 norm alone control the nonlinear remainder at $z\to -\infty$; a nonlinear theorem needs further boundedness or regularity assumptions.

At the minimal speed $c_{\ast }=2\sqrt{k}$, the best choice is $\gamma =\sqrt{k}$. Now ${\gamma }^2-c\gamma +k=(\gamma -\sqrt{k}{)}^2$, andThe essential spectrum has no positive part but touches zero; no positive exponential weight makes it strictly negative. Thus the critical front has no spectral gap and this argument gives no exponential decay rate there. Strict weighted spectral stability holds above the minimal speed; at the minimal speed it is marginal.

For the real-coefficient cubic-quintic amplitude equation, the trivial pattern $A=0$ exists for every $\mu$. A nonzero spatially uniform pattern has arbitrary constant phase and intensity $q=|A{|}^2>0$ satisfying $\mu +\alpha q-q^2=0$. ThusNonzero patterns exist for $\mu \ge -{\alpha }^2/4$. The upper branch $q_{+}$ exists throughout that range; a distinct positive lower branch $q_{-}$ exists only for $-{\alpha }^2/4<\mu <0$ and reaches zero at $\mu =0$.

For a homogeneous perturbation, write $A=R{e}^{i\varphi }$. Then $\dot{R}=R(\mu +\alpha R^2-R^4)$ and $\dot{\varphi }=0$. The radial eigenvalue at a nonzero pattern iswhile the phase eigenvalue is zero. Therefore the upper branch is stable modulo its arbitrary phase for $\mu >-{\alpha }^2/4$; the lower branch is unstable. The phase-neutral direction prevents asymptotic attraction to one particular complex amplitude, but permits orbital stability of the circle of equivalent patterns.

Conduction is asymptotically stable for $\mu <0$ and unstable for $\mu >0$. At $\mu =0$, the positive cubic term destabilizes arbitrarily small nonzero amplitudes, even though its linear eigenvalue vanishes. At $\mu =-{\alpha }^2/4$, the merged nonzero branch is semistable from above only: $\dot{R}=-R(R^2-\alpha /2{)}^2$, so amplitudes just below it eventually decrease toward zero. Thus the bistable interval is $-{\alpha }^2/4<\mu <0$, with the endpoints treated separately rather than declared stable merely because an eigenvalue is zero.

Eliminate the instantaneous Stokes flow velocity in favour of temperature. On a horizontal Fourier mode $k$, the Stokes temperature-slaving operator maps $\theta$ to $w={\mathcal{W}}_R\theta$, where $(D^2-k^2{)}^{2}w=R{k}^2\theta$ and $w=w^{\prime \prime }=0$. The temperature evolution has linear operator ${\mathcal{L}}_R=D^2-k^2+{\mathcal{W}}_R$ and bilinear map $\mathcal{N}(\varphi ,\chi )=-{\mathbf{u}}_R[\varphi ]\cdot \nabla \chi$. Under the homogeneous thermal Dirichlet boundary conditions, ${\mathcal{L}}_{R_c}$ is self-adjoint. Normalize its critical eigenfunction as $g(z)=\sin \pi z$ and set $s=k_c^2+{\pi }^2$; the critical vertical velocity is $f=sg$.

At order $ϵ$, the critical eigenfunction equation gives ${\mathcal{L}}_{R_c}(g{e}^{i{k}_{c}x})=0$. At order ${ϵ}^2$, the weakly nonlinear expansion contains the imposed second harmonic and the quadratic products of the critical mode: a horizontally uniform temperature correction proportional to $|A{|}^2$ and, in a general vertical-mode calculation, a second harmonic proportional to $A^2$. These corrections are found by solving the noncritical boundary value problems, with homogeneous thermal data except for the imposed forcing.

At order ${ϵ}^3$, the method of multiple scales produces the slow derivative $A_{\tau }$, the detuning term $R_2{\partial }_R{\mathcal{L}}_R$, and the two cross-advection terms involving first- and second-order fields. Project the $e^{i{k}_{c}x}$ component onto the adjoint eigenfunction using the vertical inner product. This is the solvability condition in the method of multiple scales: divide each resonant projection by $⟨g,g{⟩}_z$. The detuning supplies $\mu A$ with $\mu \propto R_2$; interactions of horizontal wavenumbers $2k_c$ and $-k_c$ permit $\beta A^{\ast }$ with $\beta \propto \gamma$; self-interaction through the slaved mean and second harmonic supplies $-\lambda |A{|}^{2}A$. Other products have the wrong horizontal wavenumber. Reflection $x\mapsto -x$ permits real coefficients with this cosine forcing. Thus the symmetry-allowed spatially forced convection amplitude equation isThere is a useful specialization that should not be silently missed. For the literal one-vertical-mode Stokes flow problem, the vanishing two-to-one forcing coefficient for Stokes convection makes $\beta =0$ at this order. To see this, write a positive second-harmonic forcing component as $(w,\theta )=(W(z),\Theta (z))e^{2i{k}_{c}x}$, incorporating the cosine's factor $1/2$. Its coupling to the negative critical harmonic has projected integrand, apart from sign and its factor $A^{\ast }$,The integral vanishes because $g=0$ at both plates, even though $\Theta (0)$ is nonzero. This proves the cancellation without solving the forced profiles. The permitted coefficient is therefore zero times $\gamma$; symmetry alone does not establish nonzero phase pinning for the equations actually supplied.

The same normalization makes the remaining coefficients explicit. Since ${\mathcal{W}}_{R}g=(R{k}_c^2/s^2)g$, $\mu =R_2{k}_c^2/s^2$. The quadratic second harmonic cancels for $f=sg$, while the uniform correction is ${\theta }_{20}=-s|A{|}^2\sin (2\pi z)/(2\pi )$. Projecting $-w_1{\partial }_z{\theta }_{20}$ gives $-s^2|A{|}^{2}A/2$. Thus for the literal model and this temperature normalization,A generic nonzero $\beta$ would require a nonvanishing projection in an amended physical model or a different forcing structure. It is still meaningful to classify the real-coefficient amplitude equation requested independently.

Write $A=x+iy$. Then $x_{\tau }=(\mu +\beta -\lambda (x^2+y^2))x$ and $y_{\tau }=(\mu -\beta -\lambda (x^2+y^2))y$. These are a gradient flow for $V=-\frac{1}{2}(\mu +\beta )x^2-\frac{1}{2}(\mu -\beta )y^2+\frac{\lambda }{4}(x^2+y^2{)}^2$, so local minima give stable equilibrium points. At the origin the two eigenvalues are $\mu +\beta$ and $\mu -\beta$. The origin has exponential asymptotic stability if $\mu <-|\beta |$, retains asymptotic stability with algebraic decay at $\mu =-|\beta |$, and is unstable if $\mu >-|\beta |$. At equality, ${\rho }^2=x^2+y^2$ obeys $({\rho }^2{)}_{\tau }\le -2\lambda {\rho }^4$, since both linear coefficients are nonpositive. Integrating this inequality proves attraction even in the zero-eigenvalue direction.

For $\beta >0$ the stable nonzero equilibrium points are real; for $\beta <0$ they are imaginary:The real branch has Jacobian matrix eigenvalues $-2(\mu +\beta ),-2\beta$; the imaginary branch has $-2(\mu -\beta ),2\beta$. The oppositely aligned branch, when it exists, is a saddle equilibrium. No mixed real-imaginary nonzero equilibrium is possible when $\beta \ne 0$.

For $\beta =0$, the origin is stable for $\mu \le 0$, with algebraic decay at zero. If $\mu >0$, the circle $|A{|}^2=\mu /\lambda$ is radially attracting. Each point has Lyapunov stability but has a neutral phase direction, so it does not have individual asymptotic stability; the circle has orbital stability. This is the literal model's unpinned family. The general nonzero-$\beta$ branches instead exhibit phase locking to one of two phases separated by $\pi$.

Define the requirement before listing planetary properties. For known life, a habitable environment must supply a suitable solvent, accessible chemical ingredients, usable free energy, and conditions allowing organized organic chemistry to persist for sufficient time. Habitability is the capacity to support life, not evidence that abiogenesis or colonization actually occurred. Here the main target is long-lived near-surface liquid water on a rocky planet; a local subsurface habitat can satisfy the biological requirements even if most of the surface is hostile.

Liquid water and compatible conditions are the central environmental requirements. The local temperature, pressure and composition must permit liquid water with adequate water activity. The pure-water triple point, approximately $273.16\mathrm{K}$ and $611\mathrm{Pa}$, is a useful reference, not a universal habitable threshold: dissolved salts alter the phase diagram, and subsurface pressure differs from atmospheric pressure. A brine being liquid does not establish biological usability if its water activity is too low. There is no reason to require a global mean surface temperature of exactly $288\mathrm{K}$, or to treat $0$ to $100{}^{\circ }\mathrm{C}$ at one atmosphere as a universal range for all habitats.

Composition, nutrients and energy must be accessible together. A planet needs reservoirs of water and biogenic elements, notably carbon, hydrogen, nitrogen, oxygen, phosphorus and sulfur, in forms available to a habitat. Bulk abundance is not enough if nutrients are locked away from reacting fluids. pH, salinity, toxic solutes and water-rock interaction influence accessibility. Photosynthesis can use sunlight, while chemolithotrophy can exploit redox disequilibrium generated by water-rock interaction, for example hydrogen production during serpentinization. A warm environment at chemical equilibrium supplies no sustained metabolic free energy merely because it is warm. Conversely, atmospheric oxygen and an ozone layer are not prerequisites for all life: they favor particular surface environments and metabolisms, whereas anaerobic life uses other pathways.

An appropriate atmosphere and surface energy balance favor surface habitability. The inventory and composition of the atmosphere control surface pressure, the greenhouse effect, heat transport, and protection from some harmful ultraviolet radiation. Too little warming can freeze the surface; too much can produce a runaway greenhouse and eventual water loss. The albedo, clouds, ocean coverage and atmospheric circulation matter alongside atmospheric carbon dioxide and water vapor. An illustrative globally averaged radiative balance iswhere $A$ is the Bond albedo, $a$ the orbital semimajor axis, $L_{\star }$ the specified stellar luminosity, and $\sigma$ the Stefan-Boltzmann constant. This determines an planetary effective radiating temperature, not the surface temperature. The latter requires a model of the greenhouse effect and heat transport. The stellar luminosity here is an imposed boundary condition; the planetary variables under discussion are the orbit, albedo and atmospheric properties, not characteristics of the host system.

Mass and interior properties govern retention and replenishment. Planetary mass, radius and surface gravity influence atmospheric escape, internal pressure, and cooling. In a simple thermal-retention diagnostic,where $v_{\mathrm{esc}}$ is the surface escape velocity, $G$ the gravitational constant, $M_p$ and $R_p$ the planetary mass and radius, and $m$, $T_{\mathrm{exo}}$ and $R_{\mathrm{exo}}$ refer to a gas particle and the exobase. A large Jeans escape parameter $\lambda$ suppresses ordinary thermal escape of that species; it does not exclude hydrodynamic escape or nonthermal losses. There is no single demonstrated minimum Earth-like mass valid for every atmospheric composition, age and habitat. Radiogenic heating, retained formation heat and possibly tidal heating sustain internal activity, volcanic outgassing and subsurface liquid water; excessive heating can be destructive. A massive primordial hydrogen envelope can also make a nominally rocky body unlike a temperate terrestrial surface.

The orbit must supply a tolerable irradiation history. With the host irradiation held specified, the semimajor axis must be compatible with the particular planet's radiation budget, not merely with a nominal distance label. The circumstellar habitable zone is a conditional surface-liquid water criterion; it neither guarantees habitability nor excludes internally heated subsurface habitats outside it. Orbital eccentricity $e$ produces an apoapsis-to-periapsis flux ratioFor fixed semimajor axis, the orbit-averaged incident flux is $⟨F⟩=L_{\star }/(4\pi a^2\sqrt{1-e^2})$. Both the average and the extremes matter; thermal inertia may buffer seasonal changes, so a perfectly circular orbit is favorable in some cases but not necessary. Sustained orbital stability matters more than requiring $e=0$.

Rotation and obliquity modify, rather than uniquely determine, the outcome. Obliquity, rotation period and their evolution control seasonality, the distribution of sunlight and the efficiency of atmospheric and oceanic heat transport. Very strong variations may repeatedly remove surface habitats; other configurations can retain refuges. Tidal locking does not automatically make a planet uninhabitable if an adequate atmosphere or ocean redistributes heat. Neither Earth's particular day length nor its particular obliquity is a necessary biological condition. Tidal heating should be assessed as a planetary energy budget, without assuming that more heating is always beneficial.

Long-lived recycling is favorable; its mechanism need not be unique. Exchange between interior and surface reservoirs can replace gases lost by atmospheric escape or burial, expose nutrients, and moderate climate through the carbonate-silicate cycle. Plate tectonics is one effective mechanism, but stagnant-lid planets can have volcanism, chemical weathering and crustal recycling sufficient for modeled intervals of surface habitability. These are conditional model results, not proof that such planets are inhabited.^[ref]

A planetary magnetic field can alter charged-particle access and atmospheric loss, but it is not a universal shield against all atmospheric escape and is not established as a necessary condition for habitability. An atmosphere, water column or rock overburden can also protect habitats. Likewise, exposed continents, a large moon, modern atmospheric oxygen and plate tectonics should not be promoted from plausible advantages to universal necessities. The defensible test is whether this particular body maintains accessible liquid water, chemical ingredients and usable free energy in sufficiently persistent environments. A planet can pass that test without any demonstration that life arose there.

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.

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.

In a supercritical Hopf bifurcation, a small periodic orbit attracts within the critical centre manifold on the side where the equilibrium point repels within that manifold. Orient the parameter so that its critical eigenvalue real part increases through zero. For the exact Hopf normal form $\dot{z}=(\mu +i\omega )z-|z{|}^{2}z$, with $\omega >0$, the nonzero constant radius solvesFor $0<r<r_{\ast }$ the radius grows; for $r>r_{\ast }$ it shrinks. A radial disturbance $\delta r$ therefore has leading equation $\dot{\delta r}=-2\mu \delta r$, since ${\partial }_r(\mu r-r^3){|}_{r_{\ast }}=-2\mu <0$. This is orbital stability: nearby orbits approach the circle, but their phase need not match a chosen representative. The exact angular frequency is $\omega$, so the period is $2\pi /\omega$ and the frequency is $\omega /(2\pi )$.

The Jacobian matrix at zero is $\left(\begin{array}{cc}\mu & -\omega \\ \omega & \mu \end{array}\right)$, with eigenvalues $\mu \pm i\omega$. Zero attracts for $\mu <0$ and repels for $\mu >0$. Exactly at $\mu =0$, $\dot{r}=-r^3$ givesThus zero still has asymptotic stability in this model, with algebraic decay, and there is no nonzero limit cycle at the threshold itself. As $\mu \downarrow 0$ on the cycle branch, both its radius $\sqrt{\mu }$ and its radial attraction rate $2\mu$ vanish. Small cycles consequently become slow to settle in radius even though their angular frequency stays fixed.

// Target: dynamical-systems.bigb