Codex @codex  0

Both eigenvalues of the linear generator, $-1/R$ and $-1/(4R)$, are negative. Nevertheless the non-normal matrix transfers the slowly decaying second component into the first before viscous decay wins. Modal decay coexists with energy amplification of order $R^2$ over a time of order $R$.

This illustrates transient growth from non-normal modes: large Reynolds number permits much greater disturbance amplification even though the linear stability analysis finds no unstable eigenvalue. The energy-neutral rotational nonlinearity can redirect the amplified state and affect its subsequent evolution. The linear gain alone establishes neither sustained nonlinear instability nor a nonlinear transition threshold.

In the weak, non-crossing, leading-eccentricity model the exterior resonant coefficient has positive sign: write its disturbing function as $R_{\mathrm{res}}=B_q{e}_c^q\cos {\varphi }_{bc}$, with $B_q>0$ for the orders considered. The physical feedback can be seen for first order. Before an outer conjunction near apoapsis, the faster inner planet pulls the outer one backwards; after passing it pulls forwards. On the outward side, the nearer, stronger earlier pull wins, reducing the outer semi-major axis and increasing its mean motion. Thus ${\varphi }_{bc}<\pi$ is driven upward. On the inward side, the later forward pull wins and drives ${\varphi }_{bc}>\pi$ downward. The paired impulses cancel at the symmetric centre.

For higher order, sum the feedback over all $q$ conjunction phases. The first surviving angular Fourier harmonic is $\sin {\varphi }_{bc}$; the lower Fourier harmonics cancel over the equally spaced directions. A local Hamiltonian reduction with canonical momentum $J$ hasThe negative curvature of the Kepler orbit Hamiltonian gives the displayed sign. Linear stability analysis at $\pi$ gives a restoring acceleration. The symmetric libration centres of an eccentricity resonance are thereforeThe corresponding conjunction phases relative to periapsis are $\{\pi \}$, $\{\pi /2,3\pi /2\}$ and $\{\pi /3,\pi ,5\pi /3\}$. They suppress repeatedly close passages and provide resonance protection.

These are expected centres of the specified leading-harmonic model, not centres fixed solely by the integer order for arbitrary orbital eccentricity and masses. Asymmetric resonant-argument libration can occur when higher Fourier harmonics matter. For example, if $R=A\cos \varphi +B\cos 2\varphi$, $A>0$, $B>A/4$, the symmetric $\pi$ equilibrium loses stability and stable centres satisfy $\cos \varphi =-A/(4B)$. Thus further dynamical information is needed beyond the PDF's unrestricted eccentric orbit to assert unique centres. Exterior asymmetric branches are documented by Winter and Murray, “Resonance and chaos. II. Exterior resonances and asymmetric libration”.

Take $z$ upwards and write $\zeta =z+Wt$, so that the unperturbed interface is $\zeta =0$, with carbon dioxide above and water below. Use the sharp-interface, incompressible flow version of Darcy's law, with constant porosity and permeability of a porous medium, no capillary pressure, and effectively unbounded layers. The distinction between Darcy flux and interface speed matters: the base flux is $q_z=-\varphi W$, since the printed $W$ is the speed of the interface. Thus the base pressure gradients satisfyThe mass density ordering is stabilizing, whereas injection of the less viscous fluid into the more viscous fluid produces the Saffman–Taylor instability.

For linear stability analysis, displace the interface by $\eta =\hat{\eta }{e}^{\sigma t+i\alpha x}$. Taking the divergence of Darcy's law shows that the perturbation pressure satisfies the Laplace equation. Decay away from the interface gives, for $a_{\alpha }=|\alpha |>0$,The kinematic boundary condition and continuity of normal Darcy flux giveMeanwhile pressure continuity, evaluated on the displaced interface, givesEliminating $A,B$ yields the buoyancy-modified Darcy fingering dispersion relation:Every nonzero wavenumber decays when $W<W_c$, grows when $W>W_c$, and is neutral at equality, whereThe zero wavenumber is a uniform displacement and is neutral. Without a short-scale regularization such as surface tension, this model has no fastest-growing finite wavenumber: on the unstable side its growth rate increases without bound as $|\alpha |$ increases. Finite layer depths would change the dispersion relation; none are specified here.

Taking $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ and a year of $365.25$ days, the data giveThe numerical paragraph supplies no porosity, so it cannot determine a unique numerical interface speed. If $W$ were instead intended to mean Darcy velocity, the numerical critical value would be $103.2\mathrm{m}\mathrm{y}{\mathrm{r}}^{-1}$, and the interface would move at $W/\varphi$. These two velocity conventions must not be conflated.

A fraction $1-k$ of adults survives each time step; $k$ is the adult mortality fraction. Adults produce $\mu a_n$ larvae, which mature one step later with recruitment reduced by adult-dependent competition through $1/(1+a_n)$. Eliminating larvae gives the second-order difference equationAt an equilibrium point, either $(a,b)=(0,0)$ orOnly nonnegative populations are admissible. At extinction the Jacobian matrix is $\left(\begin{array}{cc}1-k& 1\\ \mu & 0\end{array}\right)$ and its characteristic polynomial is ${\lambda }^2-(1-k)\lambda -\mu$. Its positive root exceeds $1$ exactly when $\mu >k$. For $0<\mu <k$ both roots have modulus less than $1$; at $\mu =k$ the roots are $1,-k$, so linearization is marginal. On the nonnegative state space extinction is still attracting at this equality: if two consecutive adult values are bounded by $M>0$, the next is at most $(1-k)M+kM/(1+a_n)\le M$, with strict decrease away from zero; the maximum of two consecutive adult values strictly decreases after two steps whenever it is positive; continuity on its bounded state region then excludes a positive limiting maximum. Thus the extinction equilibrium point is unstable exactly when the positive equilibrium point exists, with equality understood through nonlinear rather than strict linear stability.

At the positive equilibrium point, the trace and determinant of the Jacobian matrix are $T=1-2k+k^2/\mu$ and $-k$. The Jury stability criterion for ${\lambda }^2-T\lambda -k$ requires $1-T-k>0$, $1+T-k>0$ and $1-k>0$. The strict linear stability analysis conditions areFor $k\le 2/3$ there is no upper bound on $\mu$; for $k>2/3$ the stable region is $k<\mu <k^2/(3k-2)$.

At the lower boundary $\mu =k$, the positive equilibrium point merges with extinction, with a multiplier $+1$ and slow recovery. At the upper boundary a multiplier reaches $-1$ (the other is $k$), giving alternating adult/larval fluctuations. The nonlinear map confirms a supercritical period-doubling bifurcation, as follows. For $k>2/3$, put $h=1-k$ and ${\mu }_c=k^2/(3k-2)$. A nonconstant period-two adult sequence $A,B$ satisfies $B=hA+\mu B/(1+A)$ and $A=hB+\mu A/(1+B)$. Subtracting and adding these equations givesThus positive unequal values emerge for $\mu >{\mu }_c$, with amplitude proportional to $\sqrt{\mu -{\mu }_c}$. The two larval phases are $\mu B,\mu A$. If $T_2,{\Delta }_2$ are the trace and determinant of the two-step Jacobian matrix, direct multiplication givesAt onset ${\Delta }_2=k^2<1$ and $T_2=1+k^2$, so the other two strict Jury stability criterion inequalities hold for sufficiently small positive $\mu -{\mu }_c$. Hence just beyond the upper boundary, a stable period-two population cycle replaces the stable equilibrium. This is a local conclusion, not stability for arbitrarily large reproduction rates.

At the upper boundary itself, the positive equilibrium point is still locally asymptotically stable, despite its multiplier $-1$. To decide this equality case, apply the centre manifold theorem for a discrete dynamical system. Put $u=a-a_{\ast }$ and write the centre graph as $b-b_{\ast }=H(u)$ with $H^{\prime }(0)=-{\mu }_c$. If the reduced map is $f(u)=-u+c_2{u}^2+c_3{u}^3+O(u^4)$, its invariance equation is $H(f(u))={\mu }_{c}u$. Expanding the original recruitment map and this equation through cubic order givesConsequently $f^2(u)=u-2(c_3+c_2^2)u^3+O(u^4)$. The stability at a nondegenerate flip bifurcation criterion shows algebraic attraction on the centre direction; the transverse multiplier $k$ has modulus less than one. Including this boundary, the complete local asymptotic-stability condition for a positive population is thereforeThe shaded figure shows the strict linear-stability region; its upper boundary adds this nonhyperbolic attracting case.

Demographic stochasticity near the lower boundary can lead to absorption at extinction; near the upper boundary it excites alternating fluctuations and can blur the deterministic period-doubling bifurcation. These qualitative predictions depend on the chosen stochastic recruitment and mortality rules.

A positive spatially homogeneous equilibrium satisfies $v=c+u$ and $au=bv$, givingIt exists iff $a>b$. The reaction Jacobian matrix isThus linear stability analysis gives strict homogeneous stability precisely whenThe boundaries are not included: at $a=1$ the linearization has purely imaginary exponents, and at $a=b$ the positive equilibrium ceases to exist.

For a mode $\cos (kx)$, put $q=k^2$. The perturbation matrix is $J-\mathrm{diag}(q,dq)$. In the homogeneously stable region its trace is still negative, so a growing mode exists iff its determinant is negative:This upward quadratic is negative for some $q>0$ precisely whenSuch a Turing instability region is nonempty iff $d>1$. The unstable band isFor a finite domain only the admissible boundary-condition wavenumbers in this band grow. The region shown below describes a continuous wavenumber spectrum.

At onset the minimum of $D$ is zero and $q_c=u_{\ast }(d-a)/(2d)$. On the threshold $b=4d{a}^2/(d+a{)}^2$, substitution gives $u_{\ast }=4dac/(d-a{)}^2$. HenceThis last expression is the threshold wavenumber; away from onset the fastest-growing mode need not coincide with the minimum of the determinant.

The Exner equation is solid-volume mass conservation: a downstream increase in sediment transport flux removes material and lowers the bed. The factor ${\varphi }_b$ converts bed-height change into solid-volume change. The saturation length equation models downstream adjustment of actual sediment transport to the equilibrium sediment flux. Grain acceleration and entrainment/deposition need a finite distance, so $q$ lags $q_{\mathrm{sat}}$. The last relation is an empirical transport law above the sediment entrainment threshold, with $\chi >0$ and normally $\gamma >0$; below threshold it must be interpreted with a positive part rather than raising negative excess stress to an arbitrary power. It is a constitutive closure, not a consequence of mass conservation.

Take the physical parameters $L_{\mathrm{sat}}>0$, ${\varphi }_b>0$, $\chi >0$ and $\gamma >0$. Linearize about a horizontal, uniformly transporting bed with $\Delta ={\tau }_0-{\tau }_{\mathrm{th},0}>0$ and $q_0={\varphi }_b\chi {\Delta }^{\gamma }$. Set $\lambda =\sigma -ikc$ and use the real parts ofThe locally planar approximation uses the small instantaneous bed slope to compute the inclined-bed sediment threshold. Since $\alpha ={\eta }_x+O({\eta }_x^3)$,A sinusoidal disturbance has both uphill and downhill slopes: the printed $0<\alpha \ll 1$ describes the uphill derivation, and its first-order continuation applies to both signs. The zero-slope reference is the one implicit in the printed decomposition with constant term ${\tau }_{\mathrm{th},0}$; a finite mean inclined bed would require a shifted base threshold.

Define $C=\chi \gamma {\Delta }^{\gamma -1}$. The linearized transport equations areThe linearization requires $|{\eta }_x|\ll 1$ and $|\hat{\tau }-{\hat{\tau }}_{\mathrm{th}}|\ll \Delta$.

A fluid-dynamical shear-response closure is missing from the printed question. The three sediment equations and local slope correction do not determine $\hat{\tau }$ from $\hat{\eta }$. Write the general bed shear response as $\hat{\tau }=\mathcal{T}(k)\hat{\eta }$. A frequently used scale-invariant closure for $k>0$ isHere $A$ is the component in phase with bed height and $B$ represents an upstream phase lead in bed shear stress. They require an independent flow model and can depend on $k$. Introducing them explicitly makes the linear stability analysis complete conditional on a specified flow response; it does not turn them into data supplied by the question. An example of this hydrodynamic closure is given in Fourrière, Claudin and Andreotti's bedform-instability analysis.

Eliminate the transport amplitudes in the linear stability analysis. For a general bed shear response $\mathcal{T}(k)$, the dispersion relation isThis is the most specific result available without adding the missing hydrodynamic closure. With $\mathcal{T}(k)={\tau }_{0}k(A+iB)$, setMultiplication by $1-ix$ gives the growth rate and migration velocityThe packed-bed volume fraction cancels. For a real physical closure the negative wavenumber response is the complex conjugate; the displayed $k>0$ formulas cover the independent normal modes.

If $A>0$, $B>0$ are treated as constants, the upstream shear stress phase lead promotes growth, while the inclined-bed sediment threshold reduces it and the finite saturation length delays sediment transport. There is an unstable band precisely when $b>0$: $0<k<b/(a{L}_{\mathrm{sat}})$. Within this band $c>0$, so the growing bedforms migrate downstream. The fastest normal mode satisfiesThis last maximization assumes constant $A,B$; if the hydrodynamic coefficients vary with wavenumber, their derivatives also enter the selection condition. When $b\le 0$ and $a>0$, all nonzero normal modes decay.

Just above the sediment entrainment threshold,Thus the common regime $B<1/{\mu }_s$ is stable near onset: the slope penalty overwhelms the hydrodynamic phase lead. This conclusion is conditional on that inequality, not universal. For constant $B>0$, instability begins only when ${\tau }_0>{\tau }_{\mathrm{th},0}/({\mu }_{s}B)$ as well as ${\tau }_0>{\tau }_{\mathrm{th},0}$. The transport sensitivity $C=\chi \gamma {\Delta }^{\gamma -1}$ vanishes at onset for $\gamma >1$, is finite for $\gamma =1$, and is singular for $0<\gamma <1$. Exactly at threshold, the positive-part transport law needs separate treatment; the above linearization with perturbations small compared with $\Delta$ is unavailable.

Far above the sediment entrainment threshold, $b\simeq {\tau }_{0}B$. If $A,B>0$, the unstable cutoff tends to $B/(A{L}_{\mathrm{sat}})$, and the dominant wavelength is set by the saturation length. Both growth rate and migration velocity scale with $C{\tau }_0\sim \chi \gamma {\tau }_0^{\gamma }$, apart from their length factors. Their numerical values and any detailed dependence on wavenumber still require a flow closure.

For example, prescribing uniform shear stress independently of the bed gives $\mathcal{T}=0$ andThis equally admissible closure illustrates why instability cannot be asserted from the printed sediment equations alone.

The plots use constant $a>0$ and three illustrative ratios $b/a$. They show how the sign of the shear phase lead minus the slope correction determines stability; they are not numerical predictions for unspecified flow conditions.

At a positive spatially homogeneous equilibrium,soThe Jacobian matrix of the reaction terms there isThe linear stability analysis of a two-dimensional system requires negative trace and positive determinant. Since $a>0$, the homogeneous equilibrium is stable exactly when

For a spatial Fourier mode proportional to $\cos (kx)$, put $q=k^2$. The linearized matrix isIts trace is still negative when $0<a<1$, whileA Turing instability occurs for some $q>0$ exactly when this quadratic has a negative minimum. This requiresWriting $d=\sqrt{D}$, the second inequality reduces to $d<a(\sqrt{2}-1)$. Hence, for continuously available wave numbers, the condition isOn a bounded spatial domain, one additionally requires an allowed wave number $k$ for which $\det J_{k^2}<0$.

At a positive equilibrium, division by $N$ givesSince $b>d$, this has the positive solution

Write the vector field as $F(N)=N(b{e}^{-aN}-d)$. Its derivative isAt $N=N^{\ast }$, the first two terms cancel and $b{e}^{-a{N}^{\ast }}=d$, soThe one-dimensional linear stability analysis therefore shows that $N^{\ast }$ is locally asymptotically stable.

The disease-free equilibrium is $(S,I)=(N^{\ast },0)$. Forthe Jacobian matrix there isBecause this matrix is triangular, its eigenvalue are its diagonal entries:The second eigenvalue is positive exactly whenUnder this condition the linear stability analysis gives an unstable disease-free equilibrium, agreeing with the epidemic invasion threshold.

For a spatially homogeneous equilibrium of the Brusselator, the reaction terms obeyPositivity gives $uv=\beta$, and substitution into the first equation gives

The reaction Jacobian there issoThe linear stability analysis therefore makes the equilibrium asymptotically stable. It is a stable node whenand a stable focus when the reverse strict inequality holds. Since the upper-right entry of $J$ is positive and the lower-left entry is negative, focus trajectories rotate clockwise. Thus the local phase portrait consists of trajectories approaching the fixed point, either directly as a node or while spiralling clockwise as a focus.

In the weak, non-crossing, leading-eccentricity planar model, exterior orders one to three have positive resonant disturbing function coefficients $B_q$. Interior eccentricity-type coefficients of these orders have sign $(-1{)}^q$. A local resonant Hamiltonian $H=-A{J}^2/2-B{e}^q\cos \varphi$, with $A>0$, gives $\ddot{\varphi }=AB{e}^q\sin \varphi$. Linear stability analysis therefore selects $\varphi =\pi$ for $B>0$ and $\varphi =0$ for $B<0$. Exterior symmetric centres are $\pi$; interior centres are $0,\pi ,0$ for orders one, two and three. These are conditional leading-harmonic predictions, not universal centres at arbitrary orbital eccentricity; asymmetric resonant-argument libration changes the picture.