Pith. sign in

REVIEW 3 major objections 5 minor 39 references

Oscillatory instability and stability of stationary solutions in the parametrically driven, damped nonlinear Schr\"odinger equation

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper claims that damping and parametric forcing can stabilize the otherwise blow-up-prone NLS soliton for nonlinearity exponent $\kappa\ge 2$, via a collision-induced eigenvalue mechanism that produces oscillatory stability.

desk verdict A careful analytic extension of the Barashenkov stability analysis to arbitrary κ, with a plausible but not fully verified numerical claim of oscillatory stability for κ≥2. read the letter →

arxiv 2411.17655 v2 pith:R4X4R5TW submitted 2024-11-26 nlin.PS math-phmath.MP

F. Carreño-Navas , R. Alvarez-Nodarse , N.R. Quintero This is my paper · ORCID
classification nlin.PSmath-phmath.MP MSC 35Q5535B3537K4034B24
keywords parametricallydrivendampednonlinearSchrödingerequationoscillatorystabilityinstabilitySturm-Liouvilleeigenvalueproblemsolitonnon-Kerrnonlinearityquadrupletblow-upstabilization
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper studies the parametrically driven, damped nonlinear Schrödinger equation $i\psi_t+\psi_{xx}-\psi+2|\psi|^{2\kappa}\psi=r\psi^*-i\rho\psi$ for arbitrary positive nonlinearity exponent $\kappa$. It constructs two exact stationary soliton solutions and shows that one of them, $\Psi_-$, is always unstable. The paper's central claim is that the other solution, $\Psi_+$, which is unstable for $\kappa\ge 2$ in the undriven problem, can be stabilized by the combined effect of damping and parametric forcing: once the drive crosses a critical amplitude, a collision of eigenvalues produces a complex quadruplet with $|\Lambda_i|>|\Lambda_r|$, and the soliton becomes oscillatory stable. For $\kappa<2$ the same collision produces an oscillatory instability, in line with the known $\kappa=1$ case. The results matter because they identify a mechanism by which non-Kerr solitons that would otherwise blow up can be made stable.

What carries the argument

The central object is the dissipationless eigenvalue problem obtained after the change of variables $X=\sqrt{\omega}\,x$, $T=\omega t$, and $\varepsilon=2(1-1/\omega)$, namely $L_0\tilde g_c=\Lambda f_c$ and $L_1 f_c=-\Lambda\tilde g_c$, with $L_0=-d^2/dX^2+1-\varepsilon-(2-\varepsilon)\psi^{2\kappa}$ and $L_1=-d^2/dX^2+1-(2-\varepsilon)(2\kappa+1)\psi^{2\kappa}$. The symmetries of this system force eigenvalues to appear in quadruplets $\{\Lambda,-\Lambda,\Lambda^*,-\Lambda^*\}$. Stability is decided by the inequality $\tilde\rho^2(\Lambda_i^2-\Lambda_r^2)\ge\Lambda_r^2\Lambda_i^2$, which requires $|\Lambda_i|>|\Lambda_r|$. The stability curve $r(\rho)$ separating stable and unstable regions is parametrized by $\varepsilon$ through $\rho(\varepsilon)=2\nu/(2-\varepsilon)$ and $r(\varepsilon)=\sqrt{\varepsilon^2+4\nu^2}/(2-\varepsilon)$, with $\nu=\Lambda_i\Lambda_r/\sqrt{\Lambda_i^2-\Lambda_r^2}$. For $\kappa\ge 2$, the mechanism is the collision of a real eigenvalue pair with a near-zero pair; the resulting complex quadruplet initially has $|\Lambda_i|<|\Lambda_r|$ and only acquires $|\Lambda_i|>|\Lambda_r|$ beyond a threshold $\tilde\varepsilon_c$.

What would settle it

Evolve Eq. (1) directly with parameters inside the claimed stable region for $\kappa=2.5$, for example $\rho\approx 2.5106$, $r\approx 2.5321$ (point E), starting from $\Psi_+$ plus a small perturbation; if the perturbation grows rather than decays, the oscillatory-stability claim fails. Alternatively, repeat the eigenvalue computation with $N=1001$ and larger $L$ and check whether any eigenvalue with $\operatorname{Re}(\lambda)>\tilde\rho$ appears outside the quadruplet before $\tilde\varepsilon_c$.

Watch

Extended reading notes

Core claim

The authors find two exact stationary solutions $\Psi_\pm(x,t)=B_\pm\operatorname{sech}^{1/\kappa}[\kappa\sqrt{\omega_\pm}x]\,e^{it-i\Theta_\pm/2}$, with $B_\pm=\big((\omega_\pm(\kappa+1))/2\big)^{1/(2\kappa)}$, $\omega_\pm=1\pm\sqrt{r^2-\rho^2}$, and $\Theta_\pm=\arcsin(\rho/r)$ or $\pi-\arcsin(\rho/r)$. Linearizing around these solutions and using the change of variables $X=\sqrt{\omega}\,x$, $T=\omega t$, $\varepsilon=2(1-1/\omega)$, they reduce the stability problem to the dissipationless pair $L_0\tilde g_c=\Lambda f_c$, $L_1 f_c=-\Lambda\tilde g_c$. For $\Psi_-$, a variational argument proves the spectrum contains a real positive eigenvalue, so this solution is always unstable. For $\Psi_+$, the spectrum is computed numerically. For $0.25\le\kappa<2$, a complex quadruplet emerges at $\varepsilon=\tilde\varepsilon_c$, and an analytic argument shows $|\Lambda_i|>|\Lambda_r|$, producing the familiar oscillatory instability. The principal numerical result is that for $2\le\kappa\le 3$, a stability region opens between the two branches of the curve $r(\rho)$ once $\varepsilon$ exceeds a threshold $\tilde\varepsilon_c$ (about $0.2$ for $\kappa=2.5$), where $|\Lambda_i|>|\Lambda_r|$; the otherwise unstable soliton is thereby stabilized by damping and parametric force.

Load-bearing premise

For the $\kappa\ge 2$ result, everything hinges on the numerical eigenvalue computation being complete and accurate: the paper tracks only the quadruplet and the continuum and does not perform direct time evolution, so a missed unstable eigenvalue or a discretization shift in the threshold would break the claimed stable window.

Editorial extensions

If this is right

  • For $\kappa\ge 2$, there is a concrete region in the $(\rho,r)$ parameter plane, namely the band between the two branches of $r(\rho)$, where the parametrically driven, damped NLS soliton is stable; outside this band the same soliton is unstable.
  • Stability can be controlled by tuning the nonlinearity exponent $\kappa$: smaller $\kappa$ enlarges the stable region, and for $\kappa<2$ the analytic proof that $|\Lambda_i|>|\Lambda_r|$ means the onset of the quadruplet coincides with the stability threshold.
  • Because the driven equation describes envelopes of small-amplitude breathers in parametrically driven, damped Klein-Gordon equations with higher-order nonlinearities, the stability diagram transfers to those breather solutions.
  • The $\Psi_-$ branch is unusable: it is unstable for every $\kappa>0$, so only the $\Psi_+$ branch can serve as a stable stationary state.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the oscillatory-stability mechanism is generic, the same eigenvalue-collision criterion could predict stabilization in other parametrically driven non-Kerr envelope equations, where the effective nonlinear exponent controls whether $|\Lambda_i|>|\Lambda_r|$ holds.
  • The two-branched shape of the stability curve for $\kappa\ge 2$ implies a re-entrant behavior in the drive amplitude: at fixed damping, increasing $r$ could destabilize and then restabilize the soliton, a prediction testable by direct time-domain simulation.
  • Since the driven equation breaks Galilean invariance, the stationary stability found here does not directly yield moving solitons; extending the curve to moving solutions would require a separate two-parameter continuation.
  • The quoted threshold $\tilde\varepsilon_c\approx 0.2$ for $\kappa=2.5$ comes from a fixed discretization; a convergence check with larger $N$ and $L$ could reveal how much the onset shifts, which would sharpen the stable-window prediction.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript studies the parametrically driven, damped nonlinear Schrödinger equation with power-law nonlinearity |ψ|^{2κ}ψ. It finds two exact stationary solutions Ψ±, linearizes around them, and reduces the linear stability problem to a dissipationless eigenvalue problem via the transformation Λ² = λ² − ρ̃² (Eq. (30)). For Ψ−, Appendix A proves the existence of a positive real eigenvalue, so Ψ− is always unstable. For Ψ+, the paper derives a stability criterion (Eq. (35)) and, for κ < 2, proves in Appendix B that any emerging complex quadruplet satisfies |Λ_i| > |Λ_r|, which yields oscillatory instability and allows computation of stability curves r(ρ). For κ ≥ 2, the paper relies on numerical eigenvalues to show that a real-pair collision produces a complex quadruplet and that for ε > ε̃_c the condition |Λ_i| > |Λ_r| holds, leading to a claimed oscillatory stability region between two branches of the stability curve. The abstract states this κ ≥ 2 result as the principal conclusion of the paper.

Significance. The paper contains several rigorous and useful contributions: the exact stationary solutions, the change of variables that reduces the three-parameter stability problem to a one-parameter eigenvalue problem, the analytic proof of instability of Ψ−, and the κ < 2 proof of |Λ_i| > |Λ_r| using a cited result from Weinstein. These parts are careful and internally consistent. If the κ ≥ 2 stabilization claim is confirmed by additional numerical evidence, it would be a significant extension of the known κ = 1 stability diagram to generalized nonlinearities, with potential implications for higher-order nonlinear Klein–Gordon breathers. However, the κ ≥ 2 claim is currently supported only by a numerical eigenvalue computation without a completeness check or convergence study, and the abstract's word "establishes" overstates the evidence presented.

major comments (3)
  1. [Abstract and Sec. 4.2] The claim that for κ ≥ 2 there exists a critical ε beyond which the unstable soliton becomes stable is not established by the reported evidence. The stability criterion in Eq. (35) is derived for a single eigenvalue; applying it to the tracked quadruplet is necessary but not sufficient unless all other discrete eigenvalues are shown to be harmless. The manuscript reports no count of the full discrete spectrum in the claimed stable window and no check that a second unstable eigenvalue with larger growth rate does not exist. Because for κ ≥ 2 the unforced soliton already has a real positive mode and the stabilization is a delicate real-axis collision, a missed eigenvalue or a discretization-induced shift of ε̃_c would invalidate the claimed region. The authors should add a spectral completeness check (e.g., monitoring all eigenvalues of the discretized problem for increasing N and L) and direct time integration of Eq. (1) in the claimed stable region.
  2. [Sec. 4.2, Figs. 4 and 5] The threshold ε̃_c ≈ 0.2 for κ = 2.5 and the entire two-branch structure of the stability curves are reported for a single discretization (L = 50, N = 501). The eigenvalue curves are near-tangent at the collision point, so the computed values of ρ(ε), r(ε), and the location of the claimed stability window are sensitive to discretization and to the numerical tolerance of the eigenvalue solver. A convergence table showing ε̃_c and representative curve points for (L, N) = (30, 501), (40, 801), and (50, 1001) should be included, together with an estimate of the numerical error in the reported values.
  3. [Sec. 4.2 and Sec. 5(iv)] The manuscript states in Sec. 5 that for 2 ≤ κ ≤ 3 it "numerically assess[es]" stability, yet the abstract and conclusions present the oscillatory stability as the principal established result. All statements for κ ≥ 2 are based solely on the linearized eigenvalue problem; no direct simulation of Eq. (1) is performed. Given the known nonlinear instability of the unforced NLS soliton for κ ≥ 2, a linear-stability analysis alone does not fully settle the dynamical stability question. The authors should either perform time evolution of the PDE in the claimed stable region or explicitly qualify the claim as a linear-stability prediction.
minor comments (5)
  1. [Sec. 1] The phrase "respectively" in Eq. (4) is ambiguous; the text should clarify which parameters correspond to which equation.
  2. [Sec. 3, item 1] The symmetry statement "if Λ is an eigenvalue, Λ* is also an eigenvalue" is correct, but the manuscript should also note that the complex conjugation symmetry combined with the real form of the operators implies the stated quadruplet structure; this is implicit but could be stated more cleanly.
  3. [Sec. 4.2] The sentence "Fixing κ, when ε is increased, we numerically identify in the spectrum four eigenvalues of order O(10⁻⁶)" is unclear: are these the zero modes that are numerically nonzero due to discretization, as described earlier in Sec. 4.1? Please clarify.
  4. [Appendix B, Eq. (61)] Proposition B1 is imported from Ref. [39] without proof. This is acceptable, but the precise hypotheses (on the potential and on κ) should be stated in the text, and the reference should be cited at the point of use.
  5. [Throughout] There are minor typographical issues, e.g., "overimposed" for "superimposed" in Sec. 4.1, and the phrase "r.h.s" in Sec. 2 should be spelled out.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the stability analysis is derived from the governing equation via substitution and algebraic manipulation; no fitted parameter is relabeled as a prediction.

full rationale

Walking the paper's derivation chain, the stationary solutions Ψ± follow from the phase-locking condition sinΘ=ρ/r (Eq. 12) and the radial equation (15), yielding the exact solutions (16)-(17). The linear stability problem is obtained by substituting the perturbed ansatz into Eq. (1), giving Eqs. (18)-(20), and the change of variables and ansatz borrowed from Ref. [1] convert the problem into the dissipationless eigenvalue system (32)-(33). The stability criterion (35) is an algebraic consequence of the definition Λ²=λ²−~ρ² and the condition ν≤~ρ; it is not imposed by fitting. The κ<2 result |Λ_i|>|Λ_r| is proved in Appendix B using an external standard result (Weinstein's Proposition 2.7), and the κ≥2 oscillatory-stability claim is a numerical eigenvalue computation for the same linearized operator, not a fit renamed as a prediction. The only apparent self-citation (Ref. [8], Mertens and Quintero) is background for empirical stability criteria and is not load-bearing; the cited change of variables is from Barashenkov et al., not from the present authors. The numerical-completeness concern (tracking one quadruplet, no direct time integration, no convergence table) is a correctness/evidence limitation, not a circularity: there is no exhibited equation in which the claimed result is identical to its input by construction.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claim rests on standard spectral theory, one external lemma (Weinstein), and a normal-mode ansatz for perturbations. No new entities are introduced and no parameters are fitted to data; the numerical stability curves are computed from the derived eigenvalue problem.

assumptions (4)
  • standard math Sturm-Liouville eigenvalue ordering: eigenvalues are ordered by the number of nodes of the eigenfunctions (used in Propositions A1 and A2).
    Invoked in Appendix A to conclude L1 has a negative nodeless eigenfunction from L1ψ_X=0 and that the lowest eigenvalue of L0 is -ε. Cited from Morse & Feshbach [37].
  • standard math Weinstein's Proposition B1: for κ<2, the operator L1 is nonnegative on the subspace orthogonal to ψ, i.e., inf_{⟨ξ|ψ⟩=0} ⟨ξ|L1ξ⟩ = 0.
    Cited from Weinstein (1985) [39]. Used to prove |Λ_i|>|Λ_r| for κ<2. The paper states it 'verifies' but does not prove the extension to the ε-dependent L1.
  • domain assumption Normal-mode ansatz for the perturbation (f,g)=e^{νT}Re(e^{iΩT}(fc,gc)), with λ=ν+iΩ, and the transformation (31) to a dissipationless eigenvalue problem.
    This assumes the linear stability problem separates into modes of this form and that stability of the original damped system follows from the condition ν≤ ilde ρ. The ansatz is standard and used in Ref [1], but it restricts the perturbation space to this family.
  • standard math Theorem 2.18 of Teschl [38] used to conclude L0 is strictly positive definite for ε<0.
    Cited in Proposition A2 to combine the positivity of the lowest eigenvalue and the continuous spectrum into strict positivity of the operator L0.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Oscillatory instability and stability of stationary solutions in the parametrically driven, damped nonlinear Schr\"odinger equation." pith.science (2026). https://pith.science/paper/R4X4R5TW

@misc{pith2026241117655,
  author       = {Pith},
  title        = {Pith review of: Oscillatory instability and stability of stationary solutions in the parametrically driven, damped nonlinear Schr\"odinger equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/R4X4R5TW}},
  note         = {Machine review of arXiv:2411.17655}
}
abstract

We found two stationary solutions of the parametrically driven, damped nonlinear Schr\"odinger equation with nonlinear term proportional to $|\psi(x,t)|^{2 \kappa} \psi(x,t)$ for positive values of $\kappa$. By linearizing the equation around these exact solutions, we derive the corresponding Sturm-Liouville problem. Our analysis reveals that one of the stationary solutions is unstable, while the stability of the other solution depends on the amplitude of the parametric force, damping coefficient, and nonlinearity parameter $\kappa$. An exceptional change of variables facilitates the computation of the stability diagram through numerical solutions of the eigenvalue problem as a specific parameter $\varepsilon$ varies within a bounded interval. For $\kappa <2$ , an {\it oscillatory instability} is predicted analytically and confirmed numerically. Our principal result establishes that for $\kappa \ge 2$, there exists a critical value of $\varepsilon$ beyond which the unstable soliton becomes stable, exhibiting {\it oscillatory stability}.

Figures

Figures reproduced from arXiv: 2411.17655 by the authors.

Figure 1
Figure 1. Critical value of ε˜c as a function of κ is depicted with a red dashed line. In the region above this curve, the necessary condition for stability |Λi | > |Λr | holds. Black solid line (overimposing with red dashed line for κ < 2) represents the values of ε˜ versus κ for which a quadruplet emerges. The dotted region is stable since Λr = 0. In the shadow region the condition |Λi | < |Λr | is satisfied. In the dot-das… view at source ↗
Figure 2
Figure 2. Crucial eigenvalues illustrating the emergence o [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Upper panel: The stability curves, traveled from t [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: For κ = 2.5, the real and imaginary parts of the eigenvalues which determine the soliton stability are represented as a function of ε. Parameters: L = 50 and N = 501. 5 Summary We have investigated the parametrically driven, damped nonlinear Schrödinger equation with a…
Figure 5
Figure 5. Figure 5: The stability curves are depicted with blue dashed [PITH_FULL_IMAGE:figures/full_fig_p012_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    I. V . Barashenkov, M. M. Bogdan, and V . I. Korobov. Stabil ity Diagram of the Phase-Locked Solitons in the Parametrically Driven, Damped Nonlinear Sc hrödinger Equation. Europhys. Lett., 15:113, 1991

    work page 1991

  2. [2]

    Bondila, I

    M. Bondila, I. V . Barashenkov, and M. M. Bogdan. Topograp hy of attractors of the parametrically driven nonlinear Schrödinger equation. Physica D, 87:314, 1995

    work page 1995

  3. [3]

    Wang and R

    X. Wang and R. Wei. Interactions and motions of double-so litons with opposite polarity in a parametrically driven system. Phys. Lett. A , 227:55, 1997

    work page 1997

  4. [4]

    W. Wang, X. Wang, J. Wang, and R. Wei. Dynamical behavior o f parametrically excited solitary waves in Faraday’s water trough experiment. Phys. Lett. A , 219:74, 1996

    work page 1996

  5. [5]

    Wang and R

    X. Wang and R. Wei. Oscillatory patterns composed of the p arametrically excited surface-wave solitons. Phys. Rev. E, 57:2405, 1998

    work page 1998

  6. [6]

    Gordillo and M

    L. Gordillo and M. A. García-Ñustes. Dissipation-drive n behavior of nonpropagating hydrody- namic solitons under confinement. Phys. Rev. Lett., 112:164101, 2014

    work page 2014

  7. [7]

    I. V . Barashenkov, E. V . Zemlyanaya and M. Bär. Traveling solitons in the parametrically driven nonlinear Schrödinger equation. Phys. Rev. E, 64:016603, 2001

    work page 2001

  8. [8]

    F. G. Mertens and N. R. Quintero. Empirical stability cri teria for parametrically driven solitons of the nonlinear Schrödinger equation. J. Phys. A: Math. Gen. , 53:315701, 2020. 16

    work page 2020

Show all 39 references
  1. [9]

    X. Wang. Parametrically excited nonlinear waves and the ir localizations. Phys. D: Nonlinear Phenom., 154:337, 2001

    2001

  2. [10]

    I. V . Barashenkov and E. V . Zemlyanaya. Stable Complexe s of Parametrically Driven, Damped Nonlinear Schrödinger Solitons. Phys. Rev. Lett., 83:2568, 1999

    1999

  3. [11]

    I. V . Barashenkov and E. V . Zemlyanaya. Soliton complex ity in the damped-driven nonlinear Schrödinger equation: Stationary to periodic to quasiperi odic complexes. Phys. Rev. E, 83:056610, 2011

    2011

  4. [12]

    M. M. Bogdan and O. V . Charkina. Structure of soliton bou nd states in the parametrically driven and damped nonlinear systems. Low Temp. Phys., 48:1062, 2022

    2022

  5. [13]

    J. M. Miles. Parametrically excited solitary waves. J. Fluid Mech, 148:451, 1984

    1984

  6. [14]

    A. Scott. Nonlinear Science: Emergence and Dynamics of Coherent Stru ctures, Oxford Texts in Applied and Engineering Mathematics, V ol. 8. Oxford University Press, 2003

    2003

  7. [15]

    Dauxois and M

    T. Dauxois and M. Peyrard. Physics of Solitons. Cambridge University Press, 2010

    2010

  8. [16]

    Sánchez and A

    A. Sánchez and A. R. Bishop. Collective coordinates and length-scale competition in spatially inhomogeneous soliton-bearing equations. SIAM Rev., 40:579, 1998

    1998

  9. [17]

    Sulem and C

    C. Sulem and C. L. Sulem. The Nonlinear Schrödinger Equation. New Y ork: Springer, 1999

    1999

  10. [18]

    A. Biswas. Perturbation of solitons with non-Kerr law n onlinearity. Chaos Solitons Fract., 13:815, 2002

    2002

  11. [19]

    V . S. Synakh and V . E. Zakharov. The nature of the self-fo cusing singularity. Sov. Phys. JETP , 41:465, 1975

    1975

  12. [20]

    E. W. Laedke and K. H. Spatschek. On localized solutions in nonlinear Faraday resonance. J. Fluid. Mech, 223:589, 1991

    1991

  13. [21]

    Denardo, B

    B. Denardo, B. Galvin, A. Greenfield, A. Larraza, S. Putt erman, and W. Wright. Observations of localized structures in nonlinear lattices: Domain walls a nd kinks. Phys. Rev. Lett., 68:1730, 1992

    1992

  14. [22]

    S. Longhi. Stable multipulse states in a nonlinear disp ersive cavity with parametric gain. Phys. Rev. E, 53:5520–5522, 1996

    1996

  15. [23]

    N. V . Alexeeva, I. V . Barashenkov, and G. P . Tsironis. Im purity-Induced Stabilization of Solitons in Arrays of Parametrically Driven Nonlinear Oscillators. Phys. Rev. Lett., 84:3053, 2000

    2000

  16. [24]

    M. G. Clerc, S. Coulibaly, and D. Laroze. Parametricall y driven instability in quasi-reversal sys- tems. Int. J. Bifurcation Chaos, 19:3525, 2009

    2009

  17. [25]

    G. J. de V alcárcel and K. Staliunas. Phase-bistable Ker r cavity solitons and patterns. Phys. Rev. A, 87:043802, 2013

    2013

  18. [26]

    Saxena, I

    A. Saxena, I. C. Christov, and A. Khare. Higher-Order Fi eld Theories: φ 6, φ 8 and Beyond. In Panayotis G. Kevrekidis and Jesús Cuevas-Maraver, editors , A Dynamical Perspective on the φ 4 Model: Past, Present and Future, pages 253–279. Springer International Publishing, Cham, 2019

    2019

  19. [27]

    Kirrmann, G

    P . Kirrmann, G. Schneider and A. Mielke. The validity of modulation equations for extended systems with cubic nonlinearities. Proc. Roy. Soc. Edinburgh, 122A:85, 1992

    1992

  20. [28]

    N. V . Alexeeva, I. V . Barashenkov, and D. E. Pelinovsky. Nonlinearity, 12:103, 1999. 17

    1999

  21. [29]

    Cooper, A

    F. Cooper, A. Khare, A. Comech, B. Mihaila, J. F. Dawson, and A. Saxena. Stability of exact solutions of the nonlinear Schrödinger equation in an exter nal potential having supersymmetry and parity-time symmetry. J. Phys. A: Math. Gen. , 50:015301, 2016

    2016

  22. [30]

    W. H. Press, S. A. Teukolsky, W. T. V etterling, and B. P . Flannery. Numerical Recipes in FORTRAN: the Art of Scientific Computing . Cambridge: Cambridge University Press, 1992

    1992

  23. [31]

    Carretero-González, D

    R. Carretero-González, D. J. Frantzeskakis, and P . G. K evrekidis. Nonlinear Waves & Hamiltonian Systems: From One To Many Degrees of Freedom, From Discrete To Continuum. Oxford University Press, 2024

    2024

  24. [32]

    J. M. Ortega and W. C. Rheinboldt. Chapter 7 - general ite rative methods. In J. M. Ortega and W. C. Rheinboldt, editors, Iterative Solution of Nonlinear Equations in Several V aria bles, pages 181–239. Academic Press, 1970

    1970

  25. [33]

    M. Y . Waziri, W. J. Leong, M. A. Hassan, and M. Monsi. An Ef ficient Solver for Systems of Nonlinear Equations with Singular Jacobian via Diagonal Up dating. Appl. Math. Sci. , 4:3403, 2010

    2010

  26. [34]

    Math Works, Inc., Natick, MA, 2020

    Matlab, The Language of Technical Computing . Math Works, Inc., Natick, MA, 2020

    2020

  27. [35]

    See Supplemental Material for movies showing the evolu tion of crucial eigenvalues for κ = 0.5, 1.5, 2.5 as ε increases

  28. [36]

    I. V . Barashenkov, T. Zhanlav, and M. M. Bogdan. in Nonlinear World. IV International Workshop on Nonlinear and Turbulent Processes in Physics. World Scientific, Kiev, 9-22 October 1989, 1990

    1989

  29. [37]

    P . M. Morse and H. Feshbach. Methods of Theoretical Physics , volume 1. McGraw-Hill, New Y ork, 1953

    1953

  30. [38]

    G. Teschl. Mathematical Methods in Quantum Mechanics: With Applicati ons to Schrödinger Op- erators. Graduate Studies in Mathematics, 99. American Mathematica l Society, Providence, R.I., 2009

    2009

  31. [39]

    M. I. Weinstein. Modulational Stability of Ground Stat es of Nonlinear Schrödinger Equations. SIAM J. Math. Anal. , 16:472, 1985. 18

    1985

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.