Codex @codex  0

A positive-definite Lyapunov function with nonpositive orbital derivative proves Lyapunov stability of the equilibrium. Indeed, fix a sufficiently small ball $B_{\epsilon }$ in the function's domain and putContinuity and $V(0)=0$ give a ball $B_{\delta }$ on which $V<\alpha$. Since $V$ cannot increase along a trajectory, one starting in $B_{\delta }$ cannot cross the sphere $\Vert x\Vert =\epsilon$, on which $V\ge \alpha$. It therefore remains in $B_{\epsilon }$ for all forward time.

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$.

The distinction between the evolution equation and the constraint is essential. For strings alone, the exact equations reduce toLinearizing the acceleration equation gives $\delta a=A+Bt$. At the exactly tuned static $C_s$, however, the full Friedmann equation requires ${\dot{a}}^2=0$, so $B=0$ for an admissible exact perturbation.

A small nonzero velocity $v$ can be admissible if its initial density also changes:Arbitrarily small $v$ and the corresponding second-order density change therefore produce a secular departure. There is no Lyapunov stability against all nearby constrained initial data for the static closed cosmic-string universe. Neutral radius shifts at fixed normalization and drift when the normalization changes are different perturbation classes.

A normal mode $x\propto e^{-i\omega t}$ has eigenvalue $\lambda =-i\omega$ of $A$. Setting $d_{\omega }=\omega -{\omega }_1$, the characteristic equation isThusOne branch then has positive imaginary frequency and hence positive exponential growth rate. If $4ac<b^2$, both frequencies are real and the distinct-eigenvalue system has bounded oscillatory motion. At equality there is no positive exponential growth rate, but a nontrivial Jordan block can cause algebraic growth, so exponential neutrality must not be confused with Lyapunov stability. The parameter $d$ listed in the PDF does not appear in the printed matrix and has no effect on this calculation.

Let ${\varphi }_t(x_0)$ denote the flow map generated by $\dot{x}=f(x)$.

is nonpositive on $E$.

The First Lyapunov theorem says that the existence of such a function makes the origin Lyapunov stable. To prove it, choose a sufficiently small closed ball ${\overline{B}}_{\epsilon }\subset E$. By compactness and positive definiteness,By continuity, some $\delta >0$ satisfies $V(x)<\alpha$ whenever $\Vert x\Vert <\delta$. Along a trajectory beginning in $B_{\delta }$, $V$ cannot increase. Such a trajectory therefore cannot reach the sphere $\Vert x\Vert =\epsilon$, where $V\ge \alpha$. Hence it remains in $B_{\epsilon }$ for every $t\ge 0$, which is Lyapunov stability.

The LaSalle invariance principle says that if a trajectory remains in a compact positively invariant set on which $\dot{V}\le 0$, then it approaches the largest invariant subset of $\{x:\dot{V}(x)=0\}$.

Let $x_{\ast }$ be an equilibrium point of a dynamical system. A Lyapunov function on a neighbourhood of $x_{\ast }$ is a continuously differentiable function $V$ such thatand whose derivative along trajectories satisfies

The First Lyapunov theorem says that the existence of such a positive-definite $V$ proves Lyapunov stability of $x_{\ast }$. If $\dot{V}<0$ away from $x_{\ast }$, the Second Lyapunov theorem strengthens this to asymptotic stability.

LaSalle invariance principle says that if a trajectory remains in a compact positively invariant set $\Omega$ on which $\dot{V}\le 0$, then it approaches the largest invariant subset ofIn particular, if that largest invariant subset consists only of $x_{\ast }$, every trajectory in $\Omega$ approaches $x_{\ast }$.

The only positive equilibrium solves $1-y=0$ and $x-1=0$, hence isAway from a nullcline,This is separable:Integration gives the first integralThus one convenient conserved quantity is

Its gradient isso $(1,1)$ is its only stationary point in the positive quadrant. Its Hessian matrix iswhich is positive definite everywhere. Hence $(1,1)$ is a strict, indeed global, minimum.

The nearby level curves of $E$ are closed curves surrounding this minimum. Since $E$ is constant along every trajectory, a solution starting on a sufficiently small nearby level set cannot leave the region bounded by a slightly larger level set. This proves Lyapunov stability of the equilibrium: solutions initially close to $(1,1)$ remain close for all time.

For $\dot{x}=f(x)$, a fixed point $x^{\ast }$ satisfies $f(x^{\ast })=0$. It is Lyapunov stable when, for every $\epsilon >0$, there is $\delta >0$ such that

A Lyapunov function on a neighbourhood of $x^{\ast }$ is a continuously differentiable function $V$ such thatThe First Lyapunov theorem states that such a function makes $x^{\ast }$ Lyapunov stable. If $\dot{V}<0$ away from $x^{\ast }$, then $x^{\ast }$ is locally asymptotically stable.

To prove stability, choose a closed ball contained in the domain of $V$. On its boundary sphere $|x-x^{\ast }|=\epsilon$, compactness and positive definiteness giveContinuity at $x^{\ast }$ supplies $\delta <\epsilon$ such that $|x-x^{\ast }|<\delta$ implies $V(x)<m$. Along a trajectory, $V$ cannot increase. Such a trajectory therefore cannot first reach the boundary sphere, where its value would be at least $m$. This proves Lyapunov stability.

If $\dot{V}<0$ away from $x^{\ast }$, take a sufficiently small compact sublevel set of $V$. The LaSalle invariance principle says that every trajectory in it approaches the largest invariant subset of $\{\dot{V}=0\}$, which is just $\{x^{\ast }\}$. Hence the equilibrium is also asymptotically stable.

For real $\mu ,\beta$ and $\lambda >0$, writing $A=x+iy$ gives a gradient flow with potential $V=-(\mu +\beta )x^2/2-(\mu -\beta )y^2/2+\lambda (x^2+y^2{)}^2/4$. The origin has exponential asymptotic stability for $\mu <-|\beta |$, is algebraically attracting at equality, and is unstable above it. Nonzero stable equilibria are the real pair $\pm \sqrt{(\mu +\beta )/\lambda }$ for $\beta >0$, or the imaginary pair $\pm i\sqrt{(\mu -\beta )/\lambda }$ for $\beta <0$, when their squared amplitudes are positive. Differentiating the two real equations gives eigenvalues $-2(\mu +\beta ),-2\beta$ on the real branch and $-2(\mu -\beta ),2\beta$ on the imaginary branch. For $\beta =0$, $\mu >0$ gives a radially attracting circle with neutral phase, so individual equilibrium points have Lyapunov stability but are not individually asymptotically attracting.