pith. sign in

arxiv: 2511.05142 · v2 · submitted 2025-11-07 · 🌊 nlin.PS

Stability of parametrically driven, damped nonlinear Dirac solitons

Bernardo S\'anchez-Rey , David Mellado-Alcedo , Niurka R. Quintero This is my paper

Pith reviewed 2026-05-18 00:16 UTC · model grok-4.3

classification 🌊 nlin.PS
keywords nonlinear Dirac equationparametric drivingdampingsoliton stabilitylinear stability analysisstationary solutionsnumerical simulations
0
0 comments X

The pith

One stationary soliton solution of the parametrically driven damped nonlinear Dirac equation is always unstable while the other is stabilized by large enough dissipation.

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

The paper examines the linear stability of two exact stationary solutions in the parametrically driven and damped nonlinear Dirac equation. Analysis of the eigenvalue problem from linearization shows that one solution is unstable for every choice of parameters. The second solution becomes stable once dissipation exceeds a threshold that traces a curve through parameter space, and this curve shifts with the driving frequency such that low-frequency cases remain stable everywhere. These analytic findings are checked against direct simulations that use a discretization method designed to keep errors small. The work clarifies which persistent wave structures can survive under combined driving and damping.

Core claim

Linearization of the parametrically driven damped nonlinear Dirac equation around its two exact stationary solutions produces an eigenvalue problem whose spectrum proves one solution unstable for all parameter values. For the second solution the spectrum crosses from unstable to stable when dissipation passes a critical value, yielding an explicit stability curve in parameter space whose shape depends on the driving frequency; low-frequency solitons lie entirely inside the stable region. Direct numerical evolution of the full nonlinear equation with a low-discretization-error algorithm reproduces the same stability diagram.

What carries the argument

the eigenvalue problem obtained by linearizing the driven damped nonlinear Dirac equation around each exact stationary solution

If this is right

  • The first solution cannot form a persistent structure under any combination of driving and damping.
  • Above the stability curve the second solution remains localized and stationary under continued parametric driving.
  • The location of the stability boundary moves with driving frequency, placing all low-frequency cases inside the stable region.
  • Numerical time-stepping that controls discretization error reproduces the analytic stability diagram without introducing spurious instabilities.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar stability thresholds may appear in other parametrically driven nonlinear wave equations once damping is introduced.
  • The stable regime could be used to select one soliton state over another in physical realizations of the model.
  • Extending the linear analysis to small perturbations that break exact stationarity would test how far the stability curve governs long-term behavior.

Load-bearing premise

Linearization around the exact stationary solutions captures the stability behavior of the full nonlinear system.

What would settle it

A numerical simulation or experiment in which the second solution grows or oscillates away from its stationary profile at dissipation values well above the predicted stability curve.

Figures

Figures reproduced from arXiv: 2511.05142 by Bernardo S\'anchez-Rey, David Mellado-Alcedo, Niurka R. Quintero.

Figure 1
Figure 1. Figure 1: FIG. 1. Top panel: Real part of the eigenvalues of the stability matrix (27) for the stationary solution [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Emergence of a complex quadruplet for [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Comparison between direct simulations of the NLDE (1), taking [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Left-hand panel: Eigenvalue spectrum for [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison between direct simulations of the NLDE (1), taking [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Left-hand panel: The solid blue circles represent the quadruplet for [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Relative error of the soliton charge at [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
read the original abstract

The linear stability of two exact stationary solutions of the parametrically driven, damped nonlinear Dirac equation is investigated. Stability is ascertained through the resolution of the eigenvalue problem, which stems from the linearization of this equation around the exact solutions. On the one hand, it is proven that one of these solutions is always unstable, which confirms previous analysis based on a variational method. On the other hand, it is shown that sufficiently large dissipation guarantees the stability of the second solution. Specifically, we determine the stability curve that separates stable and unstable regions in the parameter space. The dependence of the stability diagram on the driven frequency is also studied, and it is shown that low-frequency solitons are stable across the entire parameter space. These results have been corroborated with extensive simulations of the parametrically driven and damped nonlinear Dirac equation by employing a novel and recently proposed numerical algorithm that minimizes discretization errors.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The manuscript investigates the linear stability of two exact stationary solutions of the parametrically driven, damped nonlinear Dirac equation. Stability is determined by linearizing the PDE around these solutions and resolving the resulting eigenvalue problem. It is proven analytically that one solution is always unstable, confirming earlier variational analysis. For the second solution, sufficiently large dissipation is shown to guarantee stability, with the stability curve separating stable and unstable regions in parameter space determined numerically; the dependence on driven frequency is also examined, revealing that low-frequency solitons remain stable throughout the parameter space. These analytic and numerical results are corroborated by extensive time-stepping simulations that employ a novel algorithm asserted to minimize discretization errors.

Significance. If the central claims hold, the work supplies a concrete stability diagram for driven-damped nonlinear Dirac solitons, combining an analytic proof of unconditional instability for one branch with a dissipation threshold for the other. This is of interest for nonlinear wave systems in optics and condensed-matter analogs, where parametric driving and damping are common. The analytic confirmation of the variational prediction and the identification of a low-frequency stability regime constitute clear strengths; the numerical location of the stability boundary, however, rests on a recently proposed discretization whose error control is not yet fully documented.

major comments (1)
  1. The stability curve separating stable and unstable regions for the second solution is obtained by numerical resolution of the linearized eigenvalue problem. No convergence study, a priori error bound, or comparison against a standard spectral or finite-difference scheme is supplied for either the eigenvalue solver or the time-stepping code that corroborates the diagram. Because the dissipative Dirac operator is sensitive to high-frequency modes and because the analytic proof does not locate the curve, the precise position of the stability boundary remains the least secure element of the central claim.
minor comments (1)
  1. The abstract refers to 'the driven frequency' without specifying its symbol or range; a brief parenthetical definition would improve immediate readability.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive overall assessment of our manuscript and for the detailed comment on the numerical aspects. We respond to the major comment below and have incorporated revisions to strengthen the numerical validation.

read point-by-point responses

  1. Referee: The stability curve separating stable and unstable regions for the second solution is obtained by numerical resolution of the linearized eigenvalue problem. No convergence study, a priori error bound, or comparison against a standard spectral or finite-difference scheme is supplied for either the eigenvalue solver or the time-stepping code that corroborates the diagram. Because the dissipative Dirac operator is sensitive to high-frequency modes and because the analytic proof does not locate the curve, the precise position of the stability boundary remains the least secure element of the central claim.

    Authors: We acknowledge that the original manuscript did not include an explicit convergence study or direct comparisons with alternative schemes. The discretization employed for both the eigenvalue problem and the time-stepping simulations is the recently proposed scheme referenced in the manuscript, which is constructed to minimize truncation errors for high-frequency components in dissipative Dirac systems. To address the referee's concern, the revised version now contains a dedicated subsection reporting a grid-convergence study for the eigenvalue solver (showing stabilization of the critical dissipation threshold under successive refinements) together with selected benchmark comparisons of the time-stepping results against a standard Fourier spectral method. These additions document the robustness of the reported stability boundary without altering any of the central analytic or numerical conclusions. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper establishes stability via direct linearization of the parametrically driven damped nonlinear Dirac equation around its two exact stationary solutions, yielding an eigenvalue problem whose analytic resolution proves one solution always unstable (confirming prior variational results) and whose numerical resolution determines the stability curve for the second solution as a function of dissipation and driving frequency. These analytic and numerical findings are then corroborated by time-stepping simulations that employ a recently proposed algorithm for minimizing discretization errors. No load-bearing step reduces by construction to a fitted parameter, self-referential definition, or unverified self-citation chain; the eigenvalue problem and its resolution constitute independent content relative to the input stationary solutions, and the numerical corroboration is presented as external verification rather than a re-derivation of the same quantities.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The analysis relies on standard mathematical properties of linearization and eigenvalue problems for PDEs together with the existence of the two exact stationary solutions; no free parameters are fitted to data and no new physical entities are postulated.

axioms (1)
  • standard math Linearization of the nonlinear PDE around an exact stationary solution yields an eigenvalue problem whose spectrum determines linear stability.
    Invoked when the authors state that stability is ascertained through resolution of the eigenvalue problem stemming from linearization.

pith-pipeline@v0.9.0 · 5457 in / 1226 out tokens · 37169 ms · 2026-05-18T00:16:18.926489+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    R. Z. Sagdeev, D. A. Usikov, and G. M. Zaslavsky. Nonlinear Physics: From the Pendulum to Turbulence and Chaos . Contemporary concepts in physics. Harwood Academic Publishers, 1988

    work page 1988

  2. [2]

    F. G. Mertens, N. R. Quintero, F. Cooper, A. Khare, and A. Saxena. Nonlinear Dirac equation solitary waves in external fields . Phys. Rev. E , 86:046602, 2012

    work page 2012

  3. [3]

    E. N. Tsoy, I. M. Allayarov, and F. Kh. Abdullaev. Stable localized modes in asymmetric waveguides with gain and loss . Opt. Lett. , 39:4215, 2014

    work page 2014

  4. [4]

    V. V. Konotop and D. A. Zezyulin. Families of stationary modes in complex potentials. Opt. Lett. , 39:5535, 2014

    work page 2014

  5. [5]

    N. V. Alexeeva, I. V. Barashenkov, and A. Saxena. Spinor solitons and their PT -symmetric offspring . Ann. Phys. , 403:198, 2019

    work page 2019

  6. [6]

    D. J. Kaup. Perturbation theory for solitons in optical fibers. Phys. Rev. A , 42:5689, 1990

    work page 1990

  7. [7]

    J. Yang. Complete eigenfunctions of linearized integrable equations expanded around a soliton solution . J. Math. Phys. , 41:6614, 2000

    work page 2000

  8. [8]

    E. G. Charalampidis and S. I. Mistakidis. Two-component droplet phases and their spectral stability in one dimension . Phys. Rev. A , 111:013318, 2025

    work page 2025

  9. [9]

    Nogami, F

    Y. Nogami, F. M. Toyama, and Z. Zhao. Nonlinear Dirac soliton in an external field . J. Phys. A: Math. Gen. , 28:1413, 1995

    work page 1995

  10. [10]

    S\' a nchez and A

    A. S\' a nchez and A. R. Bishop. Collective Coordinates and Length-Scale Competition in Spatially Inhomogeneous Soliton-Bearing Equations . SIAM Rev , 40:579, 1998

    work page 1998

  11. [11]

    S\'anchez-Rey, J

    B. S\'anchez-Rey, J. Casado-Pascual, and N. R. Quintero. Kink ratchet induced by a time-dependent symmetric field potential. Phys. Rev. E , 94:012221, 2016

    work page 2016

  12. [12]

    K. V. Yershov, V. P. Kravchuk, D. D. Sheka, O. V. Pylypovskyi, D. Makarov, and Yu. Gaididei. Geometry-induced motion of magnetic domain walls in curved nanostripes . Phys. Rev. B , 98:060409, 2018

    work page 2018

  13. [13]

    B. Malomed. Variational methods in nonlinear fiber optics and related fields . Prog. in Optics, edited by E. Wolf , 43:71, North-Holland, Amsterdam, 2002

    work page 2002

  14. [14]

    A. I. Maimistov. Evolution of solitary waves which are approximately solitons of a nonlinear Schr\"odinger equation . Zh. Eksp. Teor. Fiz. , 104:3620, 1993

    work page 1993

  15. [15]

    N. R. Quintero, F. G. Mertens, and A. R. Bishop. Generalized traveling-wave method, variational approach, and modified conserved quantities for the perturbed nonlinear Schr\"odinger equation . Phys. Rev. E , 82:016606, 2010

    work page 2010

  16. [16]

    D. J. Gross and A. Neveu. Dynamical symmetry breaking in asymptotically free field theories . Phys. Rev. D , 10:3235, 1974

    work page 1974

  17. [17]

    S. Y. Lee, T. K. Kuo, and A. Gavrielides. Exact localized solutions of two-dimensional field theories of massive fermions with Fermi interactions . Phys. Rev. D , 12:2249, 1975

    work page 1975

  18. [18]

    Dreisow, M

    F. Dreisow, M. Heinrich, R. Keil, A. T \"u nnermann, S. Nolte, S. Longhi, and A. Szameit. Classical Simulation of Relativistic Zitterbewegung in Photonic Lattices . Phys. Rev. Lett. , 105:143902, 2010

    work page 2010

  19. [19]

    S. Longhi. Classical simulation of relativistic quantum mechanics in periodic optical structures . Appl. Phys. B , 104:453, 2011

    work page 2011

  20. [20]

    J. M. Zeuner, N. K. Efremidis, R. Keil, F. Dreisow, D. N. Christodoulides, A. T\"unnermann, S. Nolte, and A. Szameit. Optical analogues for massless Dirac particles and conical diffraction in one dimension . Phys. Rev. Lett. , 109:023602, 2012

    work page 2012

  21. [21]

    T. X. Tran, S. Longhi, and F. Biancalana. Optical analogue of relativistic Dirac solitons in binary waveguide arrays . Ann. Phys. , 340:179, 2014

    work page 2014

  22. [22]

    L. H. Haddad and L. D. Carr. The nonlinear Dirac equation in Bose-Einstein condensates: Foundation and symmetries . Phys. D: Nonlinear Phenom. , 238:1413, 2009

    work page 2009

  23. [23]

    L. H. Haddad, C. M. Weaver, and L. D. Carr. The nonlinear Dirac equation in Bose-Einstein condensates: I. Relativistic solitons in armchair nanoribbon optical lattices . New J. of Phys. , 17:063033, 2015

    work page 2015

  24. [24]

    A. U. Khawaja. Exact localized and oscillatory solutions of the nonlinear spin and pseudospin symmetric Dirac equations . Phys. Rev. A , 90:052105, 2014

    work page 2014

  25. [25]

    Ar\'evalo and L

    E. Ar\'evalo and L. Morales-Molina. Nonlinear excitations in the honeycomb lattice: Beyond the high-symmetry points of the band structure . Phys. Rev. A , 93:053816, 2016

    work page 2016

  26. [26]

    F. G. Mertens, F. Cooper, N. R. Quintero, S. Shao, A. Khare, and A. Saxena. Solitary waves in the nonlinear Dirac equation in the presence of external driving forces . J. Phys A: Math. and Theor. , 49:065402, 2016

    work page 2016

  27. [27]

    N. R. Quintero, S. Shao, R. Alvarez-Nodarse, and F. G. Mertens. Externally driven nonlinear Dirac equation revisited: theory and simulations . J. Phys. A: Math. Theor. , 52:155401, 2019

    work page 2019

  28. [28]

    N. R. Quintero and B. Sánchez-Rey. Exact stationary solutions of the parametrically driven and damped nonlinear Dirac equation . Chaos , 29:093129, 2019

    work page 2019

  29. [29]

    Alvarez and B

    A. Alvarez and B. Carreras. Interaction dynamics for the solitary waves of a nonlinear Dirac model . Phys. Lett. A , 86:327, 1981

    work page 1981

  30. [30]

    I. V. Barashenkov, M. M. Bogdan, and V. I. Korobov. Stability Diagram of the Phase-Locked Solitons in the Parametrically Driven, Damped Nonlinear Schr\"odinger Equation . Europhys. Lett. , 15:113, 1991

    work page 1991

  31. [31]

    I. V. Barashenkov, E. V. Zemlyanaya, and T. C. van Heerden. Time-periodic solitons in a damped-driven nonlinear Schr\"odinger equation . Phys. Rev. E , 83:056609, 2011

    work page 2011

  32. [32]

    N. V. Alexeeva, I. V. Barashenkov, A. A. Sukhorukov, and Yu. S. Kivshar. Optical solitons in PT -symmetric nonlinear couplers with gain and loss. Phys. Rev. A , 85:063837, 2012

    work page 2012

  33. [33]

    Garnier, F

    J. Garnier, F. Kh. Abdullaev, and M. Salerno. Solitons in strongly driven discrete nonlinear Schr\"odinger-type models . Phys. Rev. E , 75:016615, 2007

    work page 2007

  34. [34]

    N. V. Alexeeva, I. V. Barashenkov, and D. E. Pelinovsky. Dynamics of the parametrically driven NLS solitons beyond the onset of the oscillatory instability . Nonlinearity , 12:103, 1999

    work page 1999

  35. [35]

    Bondila, I

    M. Bondila, I. V. Barashenkov, and M. M. Bogdan. Topography of attractors of the parametrically driven nonlinear Schr \"o dinger equation . Physica D , 87:314, 1995

    work page 1995

  36. [36]

    Carreño-Navas, R

    F. Carreño-Navas, R. Alvarez-Nodarse, and N. R. Quintero. Oscillatory instability and stability of stationary solutions in the parametrically driven, damped nonlinear Schrödinger equation . Phys. D: Nonlinear Phenom. , 476:134611, 2025

    work page 2025

  37. [37]

    N. R. Quintero, B. S \'a nchez-Rey, F. Cooper, and F. G. Mertens. Length-scale competition in the parametrically driven nonlinear Dirac equation with a spatially periodic force . J. Phys. A: Math. Theor. , 52:285201, 2019

    work page 2019

  38. [38]

    Cooper, A

    F. Cooper, A. Khare, N. R. Quintero, B. Sánchez-Rey, F. G Mertens, and A. Saxena. Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity . J. Phys. A: Math. Theor. , 53:075203, 2020

    work page 2020

  39. [39]

    Berkolaiko and A

    G. Berkolaiko and A. Comech. On Spectral Stability of Solitary Waves of Nonlinear Dirac Equation in 1D . Math. Model. Nat. Phenom. , 7:13, 2012

    work page 2012

  40. [40]

    Cuevas-Maraver, P

    J. Cuevas-Maraver, P. Kevrekidis, and A. Saxena. Solitary waves in a discrete nonlinear Dirac equation . J. Phys. A: Math. Theor. , 48:055204, 2015

    work page 2015

  41. [41]

    Chugunova

    M. Chugunova. Spectral stability of nonlinear waves in dynamical systems . Doctoral Thesis, McMaster University, Hamilton, Ontario, Canada, 2007

    work page 2007

  42. [42]

    T. I. Lakoba. Numerical study of solitary wave stability in cubic nonlinear Dirac equations in 1D . Phys. Lett. A , 382:300, 2018

    work page 2018

  43. [43]

    Lakoba and J

    T. Lakoba and J. S. Jewell. Higher-order explicit schemes based on the method of characteristics for hyperbolic equations with crossing straight-line characteristics . Numer. Methods Partial Differ. Equ. , 37:2742, 2021

    work page 2021

  44. [44]

    Mellado-Alcedo and N

    D. Mellado-Alcedo and N. R. Quintero. Stability of nonlinear Dirac solitons under the action of external potential . Chaos , 34:013140, 2024

    work page 2024

  45. [45]

    Chugunova and D

    M. Chugunova and D. Pelinovsky. Block-Diagonalization of the Symmetric First-Order Coupled-Mode System . SIAM J. Appl. Dyn. Syst. , 5:66, 2006

    work page 2006

  46. [46]

    J. M. Ortega and W. C. Rheinboldt. Chapter 7 - General Iterative Methods . In J. M. Ortega and W. C. Rheinboldt, editors, Iterative Solution of Nonlinear Equations in Several Variables , pages 181--239 . Academic Press, 1970

    work page 1970

  47. [47]

    L. N. Trefethen. Spectral methods in MATLAB . SIAM, Philadelphia, 2000

    work page 2000

  48. [48]

    Cuevas-Maraver, P

    J. Cuevas-Maraver, P. Kevrekidis, A. Saxena, F. Cooper, and F. Mertens. Ordinary and Partial Differential Equations . Nova Sciences, New York, 2015

    work page 2015

  49. [49]

    Pelinovsky and Y

    D. Pelinovsky and Y. Shimabukuro. Transverse instability of line solitary waves in massive Dirac equations . J. Nonlinear Sci. , 26:365, 2016

    work page 2016

  50. [50]

    I. V. Barashenkov and E. V. Zemlyanaya. Oscillatory instabilities of gap solitons: a numerical study . Comput. Phys. Commun. , 126:22, 2000

    work page 2000

  51. [51]

    S. Shao, N. R. Quintero, F. G. Mertens, F. Cooper, A. Khare, and A. Saxena. Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity . Phys. Rev. E , 90:032915, 2014

    work page 2014

  52. [52]

    Alvarez and B

    A. Alvarez and B. Carreras. Energetic stability criterion for a nonlinear spinorial model . Phys. Rev. Lett. , 50:1230, 1983

    work page 1983

  53. [53]

    T. I. Lakoba. Study of instability of the Fourier split-step method for the massive Gross-Neveu model . J. Comput. Phys. , 402:109100, 2020

    work page 2020