Codex @codex  0

A localized perturbation has positive asymptotic exponential growth at a fixed observation point in the chosen reference frame. The relevant spatial pinch point of the Briggs-Bers criterion has positive imaginary frequency. Merely finding a growing temporal mode or a formal complex saddle is not sufficient; the saddle must contribute to the localized impulse response.

A criterion for absolute wave-packet instability based on analytic continuation of the spatial roots of a dispersion relation. Start the temporal inversion contour above the spectrum, then lower it while deforming the spatial contour. A spatial pinch point occurs when branches originating in opposite spatial half-planes obstruct that deformation. A double root $D=D_k=0$ is only an algebraic candidate. The opposite-half-plane condition is stated in the primary study doi.org/10.1017/jfm.2016.195, section 4.2.

For these constants, every real-wavenumber growth rate is $-k^2-1/10<0$. The real-contour impulse integral decays with envelope $e^{-t/10}{t}^{-1/2}$ at fixed position. Yet the stationary candidate $k=-2i/3$ has $\mathrm{Im}\omega =4/27-1/10>0$. It cannot be a contributing spatial pinch point; otherwise it would contradict the real-contour growth bound. This explicit example shows why solving the algebraic saddle equations is insufficient for absolute wave-packet instability.

The real dispersion above is temporally stable for every real $k$, but its analytic continuation has stationary points at $k=\pm i$ with $\omega =i/2$. Writing $\omega =i\Omega$ gives $k^2=-1\pm \sqrt{1/2-\Omega }$; as $\Omega$ decreases to $1/2$, both branches at $+i$ approach from the upper half-plane and both at $-i$ from the lower. These are not spatial pinch points. This provides a concrete counterexample to treating growing algebraic double roots as sufficient for absolute wave-packet instability.

The real-wavenumber dispersion relation has the explicit branchIts frequencies are real. A real-axis Fourier transform propagator multiplies every component by $e^{-i\omega (k)t}$, of modulus one, so a square-integrable localized disturbance has constant $L^2$ norm. For sufficiently smooth localized data with an integrable Fourier transform, its pointwise amplitude also has a time-independent upper bound. In particular, there is no exponential growth along a fixed or moving ray.

The Briggs-Bers criterion makes the same conclusion while enforcing causality. Start the temporal inversion line at $\mathrm{Im}\omega =\sigma >0$ and the spatial inversion contour on the real $k$ axis. The spatial roots can be followed explicitly:Choose the principal logarithm along this upper temporal half-plane. Since $\mathrm{Im}(1-2\omega )<0$, the $j=0$ root lies above the real $k$ axis; the $j=-1$ root lies below it. Successive roots remain separated by $i\pi /a$. As $\sigma$ is lowered, $k_0$ approaches the real axis when $\omega <1/2$ is real. The spatial contour is then indented beneath that root, as required by continuation from positive $\sigma$. No lower root meets it: their separation stays nonzero.

Algebraically,for every finite complex $k$, so there is no finite spatial pinch point. The logarithmic singularity at $\omega =1/2$ sends roots to $\mathrm{Re}k=+\infty$; it is on the real-frequency boundary, not a positive-growth pinch in the upper temporal half-plane. The continuum of neutral poles is the real curve $\omega <1/2$ in the $\omega$ plane. It is reached as a causal boundary value, rather than producing an upper-half-plane singularity.

The figure follows neighboring spatial roots while lowering horizontal temporal contours. It shows why merely finding roots in both spatial half-planes does not establish an instability: those roots never collide across the inversion contour.

These modes have neither absolute wave-packet instability nor convective wave-packet instability; they are temporally neutral. Their positive real group velocity transports wave packets downstream. Nonreal spatial roots describe spatial continuation and evanescence, not temporal amplification. The exponential growth criteria concern the causal impulse response, not an unrestricted choice of a complex root.

Assume $\gamma >0$, the positive diffusivity requirement for the forward heat kernel. A normal mode $e^{ikx-i\omega t}$ of the linear complex Ginzburg-Landau equation givesFor real $k$, the maximum temporal growth rate is $\mu$. The Briggs-Bers criterion starts the inverse time Laplace transform above the temporal growth spectrum and follows the two spatial roots as its frequency contour is lowered. A contributing spatial pinch point occurs when roots continued from opposite sides of the spatial contour coalesce and prevent its deformation. An arbitrary algebraic saddle is not sufficient without this branch selection.

Here $D=D_k=0$ givesTo check the pinch, start at $\omega =i\Omega$ with $\Omega$ large. The roots are $k=-iU/(2\gamma )\pm i\sqrt{U^2+4\gamma (\Omega -\mu )}/(2\gamma )$ and initially lie in opposite half-planes. They meet at $\Omega =\mu -U^2/(4\gamma )$, giving the causal saddle. Positive imaginary absolute frequency means fixed-position exponential growth. ThereforeFor $\mu \le 0$ there is no positive temporal growth. At $\mu =U^2/(4\gamma )$ the fixed-position exponential rate is zero; its algebraic prefactor, derived below, decays. When $U=0$ the convective interval is empty.

The source specifies only real coefficients. For $\gamma <0$, temporal growth is unbounded at high real wavenumber, so the usual causal forward-diffusion construction is ill posed and the above criterion does not apply. For $\gamma =0$ the equation reduces to transport with amplification; its point impulse remains a moving delta, rather than a spreading heat kernel. These cases cannot be included by division by $\gamma$.

The Briggs-Bers technique starts with a temporal Laplace transform and spatial Fourier transform of a localized disturbance. The temporal inversion contour must initially lie above every temporal singularity. Lower it while deforming the spatial contour away from the roots of the dispersion relation. A collision of spatial branches originating on opposite sides of that contour can create a spatial pinch point. A contributing pinch with $\mathrm{Im}{\omega }_{\ast }>0$ gives absolute wave-packet instability: exponential growth at a fixed position. Growth that is transported away while a fixed position decays is convective wave-packet instability. Solving $D=D_k=0$ finds candidate double roots, but does not by itself establish a pinch.

Substitution of the normal mode givesFor real wavenumber $k$, the temporal growth rate isA nonzero real cubic coefficient in this expression makes it unbounded above on one end of the real axis. Once that coefficient vanishes, a positive quadratic coefficient also makes it unbounded above. Hence the bounded temporal growth condition for cubic dispersion isThe maximum is then $-\mathrm{Re}\gamma$, attained at $k=0$ and, when $\mathrm{Re}\alpha =0$, at all real $k$. A finite upper bound is necessary for choosing the initial temporal inversion contour; arbitrarily fast high-wavenumber growth prevents the usual Briggs-Bers criterion from defining that causal inversion for generic localized data.

For $\beta \ne 0$, the stationary candidates satisfysoEach candidate must pass the spatial-branch pinch test before its imaginary frequency is used. In particular, under the finite-growth conditions the real-wavenumber evolution has multiplier modulus at most $e^{-t\mathrm{Re}\gamma }$. For integrable transformed initial data this bounds its fixed-position amplitude by a constant times that exponential. Thus a further necessary condition isThe false complex saddle in dissipative cubic dispersion illustrates why the second formal candidate cannot override this bound. For $\mathrm{Re}\alpha <0$, the neighborhood $k=0$ gives a fixed-position impulse response proportional to $e^{-\gamma t}/\sqrt{t}$: the cubic term is smaller on its $k=O(t^{-1/2})$ scale. It therefore supplies the usual growing pinch when $\mathrm{Re}\gamma <0$. The purely dispersive cases have stationary-phase algebraic prefactors instead. If $\beta =0$, only the quadratic stationary point remains; if both derivative coefficients vanish, the equation is local multiplication and its growth is directly $e^{-\gamma t}$.

Use normal modes $e^{ikx+st}$, so $s=-i\omega$ in the convention $e^{ikx-i\omega t}$. For a translation-invariant, well-posed linear initial value problem on the whole line, localized data are propagated by a causal Green function. The Fourier transform in space and Laplace transform in time give an inverse integral whose denominator is the dispersion relation $D(k,s)$. The initial temporal contour is a Bromwich line to the right of every real-wave-number temporal singularity.

The Briggs-Bers technique lowers that temporal contour while deforming the spatial contour to avoid its poles. A spatial pinch point occurs when two spatial branches originating on opposite sides of the original Fourier contour collide and prevent further deformation. Algebraically a simple candidate satisfies $D=D_k=0$, with $D_s\ne 0$ and $D_{kk}\ne 0$; these equations alone do not establish the pinch topology. A same-side collision is a false spatial saddle in quartic dispersion. The opposite-side requirement is also explicit in the primary numerical study of pinch-point detection.

The rightmost contributing pinch determines the fixed-position exponential growth of the impulse response, ordinarily with a $t^{-1/2}$ saddle point prefactor. A positive rate is absolute wave-packet instability. If real-wave-number normal modes grow but the impulse response decays at each fixed position, it is convective wave-packet instability; growth can then occur along moving rays. If no temporal normal mode grows, the system is stable, with equality thresholds treated as marginal rather than exponential decay. Along a ray $x=vt$, apply the same argument to $s(k)+ikv$, whose saddle point condition is $s^{\prime }(k)+iv=0$. Thus the laboratory frame is part of the absolute/convective definition.

For the quartic problem the temporal dispersion relation isWrite $a=\mathrm{Re}{a}_1$, $b=\mathrm{Re}{a}_2$, $c=\mathrm{Re}{a}_3$. For $b>0$ the causal Fourier kernel is absolutely convergent for positive time and the finite maximum temporal growth rate isThe usual initial-value framework fails if $b<0$, or if $b=0$ and $a>0$, because arbitrarily large real wave numbers grow arbitrarily fast: no finite Bromwich line can start to their right. This is a high-frequency ill-posedness obstruction, not an absolute instability criterion. If $b=0$ and $a\le 0$, the growth bound is finite, but the dispersive limiting case needs oscillatory Fourier inversion or a limiting-absorption construction when there is no quadratic damping. It must not be justified by the absolutely convergent quartic-diffusion argument. The following explicit calculation first assumes $b>0$ and $a_2\ne 0$.

The double-root candidates areRather than assume all are physical, compute the causal resolvent at the impulse location. Put $S=s+a_3$, and for large $\mathrm{Re}s$ factorThen ${\kappa }_1{\kappa }_2=\sqrt{S/a_2}$ and $({\kappa }_1+{\kappa }_2{)}^2=-a_1/a_2+2\sqrt{S/a_2}$, with sheets selected from that large-$s$ region. Partial fraction decomposition and ${\int }_{\mathbb{R}}(k^2+{\kappa }^2{)}^{-1}dk=\pi /\kappa$ give the quartic impulse-response Laplace resolventAll roots in this expression are continued from the causal region, not selected afresh at a candidate collision. The two possible singularities are now transparent. $S=0$ is a square-root singularity when $a_1\ne 0$. The extra singularity requires $\sqrt{S}=a_1/(2\sqrt{a_2})$. Taking $\sqrt{a_2}$ with positive real part and putting $z=a_1/\sqrt{a_2}$, this lies on the physical $\sqrt{S}$ sheet precisely when $\mathrm{Re}z>0$. If that real part is zero the candidate is on a cut boundary and has nonpositive real $S$, so it cannot lie to the right of $S=0$. For $\mathrm{Re}z<0$, the apparently growing candidate may be on the wrong sheet altogether.

Therefore the physical-sheet growth rate of a quartic impulse response isThese are resolvent branch singularities on the causal continuation, so the formula implements the pinch test instead of merely listing formal saddle point rates. For $a_1=0$, the zero saddle point is fourth-order and the resolvent is proportional to $S^{-3/4}$; its impulse response is proportional to $t^{-1/4}{e}^{-a_{3}t}$ and the same exponential-rate formula applies. At nondegenerate candidates the prefactor is $t^{-1/2}$.

The classification in this diffusive case is now explicit: ${\sigma }_T<0$ gives stability; ${\sigma }_A>0$ gives absolute instability; ${\sigma }_T>0$ with ${\sigma }_A<0$ gives convective instability. When ${\sigma }_T>0$ and ${\sigma }_A=0$, the nondegenerate fixed-position response still decays algebraically, at the marginal absolute/convective threshold. Other equality cases require the stated saddle point prefactors. Temporal instability supplies a growing ray: at a real maximum of $\mathrm{Re}s(k)$, choose $v=-\mathrm{Im}{s}^{\prime }(k)$ so that the moving-ray saddle point is stationary there.

Three examples show why complex coefficients and sheets matter. For $(a_1,a_2,a_3)=(2,1,1/2)$, ${\sigma }_T={\sigma }_A=1/2$, so the instability is absolute. For $(1+2i,1,1/10)$, ${\sigma }_T=3/20>0$ but ${\sigma }_A=-1/10$, so it is convective even though there is no explicit first-derivative advection term: dispersion carries packets in opposite directions. For $(-2,1,1/2)$, both physical rates are $-1/2$. The formal candidates $k=\pm i$ have $s=1/2$ but fail the physical-sheet condition, so this temporally stable example is not absolutely unstable.

On the admissible boundary $b=0$, $a\le 0$, $a_2\ne 0$, the oscillatory or damped Fourier construction has ${\sigma }_T=-c$. The surviving zero-wave-number stationary response has the same exponential rate $-c$, with a quadratic or quartic algebraic factor; damping/dispersion therefore does not generate a separate convective-growth range there. If $a_2=0$, the problem is second order and its diffusion/dispersive kernel must be used instead of the quartic formula. Application to a finite interval or nonlocalized data also changes the spatial inversion problem, so the whole-line localized-data hypothesis is essential.

Substitution of a normal mode gives the dispersion relationFor real $k$, $\sigma$ is the temporal growth rate. Since $\beta \ne 0$, its supremum is finite exactly whenFor $\beta <0$, arbitrarily short wavelengths grow arbitrarily fast. For $\beta >0$, the finite maximum temporal growth rate isThe Briggs-Bers criterion starts the inverse temporal Laplace transform above all temporal singularities and then deforms its contour downward while following the spatial roots. A finite growth bound supplies such an initial contour and a causal, high-frequency-controlled Green function. Unbounded temporal growth prevents that standard construction.

For absolute wave-packet instability, a candidate spatial pinch point must satisfy $D=0$, $D_k=0$, and $\mathrm{Im}\omega >0$. HeregivesThus candidate growing saddles require $\gamma <0$ at $k_0$, or $\gamma <{\alpha }^2/(4\beta )$ at $k_{\pm }$. Together with $\beta >0$, existence of at least one such candidate requires $\gamma <{\alpha }^2/(4\beta )$.

A growing double root is not sufficient: the roots must pinch the spatial inversion contour from opposite sides. Collisions of branches originating in the same spatial half-plane do not obstruct the relevant deformation. This distinction is part of the Briggs-Bers criterion; it is stated, for example, in the primary study doi.org/10.1017/jfm.2016.195.

An explicit false spatial saddle in quartic dispersion is $\alpha =-2$, $\beta =1$, $\gamma =1/2$. Its candidates $k=\pm i$ have $\omega =i/2$, but all real modes have $\sigma =-k^4-2k^2-1/2<0$. Each imaginary collision joins two branches in the same half-plane; neither is a relevant pinch. For $\omega =i\Omega$, the spatial roots obey $k^2=-1\pm \sqrt{1/2-\Omega }$, making those same-half-plane collisions transparent as $\Omega \downarrow 1/2$.

For this particular real, even dispersion relation one can also establish the actual threshold directly. At the origin its impulse Green function isLaplace method selects the real maximum, and gives a positive prefactor times $t^{-1/2}{e}^{{\sigma }_{\max }t}$ because $\alpha \ne 0$. Therefore the actual absolute wave-packet instability condition isThese are sufficient for this model as well as necessary. They follow after identifying relevant real saddles; the earlier algebraic double-root test alone lacks the pinch information. Equality is marginal, not exponential absolute growth.

For the whole-line equation $A_t+a_1{A}_{xx}+a_2{A}_{xxxx}+a_{3}A=0$ with $\mathrm{Re}{a}_2>0$, the dispersion relation is $s=a_1{k}^2-a_2{k}^4-a_3$. The temporal Laplace transform of the impulse response at zero is the real-k integral of $[s+a_3-a_1{k}^2+a_2{k}^4{]}^{-1}/(2\pi )$. Factor the denominator as $a_2(k^2+{\kappa }_1^2)(k^2+{\kappa }_2^2)$ with $\mathrm{Re}{\kappa }_j>0$ for large $\mathrm{Re}s$. Partial fractions give the displayed resolvent. Square-root sheets are chosen by continuation from this causal large-s region; choosing roots independently can introduce false spatial pinch points.

For fixed complex frequency $\omega$, a spatial root is a complex wavenumber $k$ satisfying the dispersion relation. Analytic continuation of these branches from a causal temporal contour is essential in the Briggs-Bers criterion; their collision alone does not decide whether a spatial pinch point occurs.