{"id":"874836a2-c8a7-403e-a919-2ce86cb1f47c","arxiv_id":"2411.17655","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For nonlinearity power κ ≥ 2, parametric driving combined with damping can stabilize an otherwise unstable soliton, a mechanism the authors call oscillatory stability.","lead":"This paper studies solitons in a parametrically driven, damped nonlinear Schrödinger equation and finds two exact stationary solutions. It shows that with the right amount of forcing and damping, an otherwise unstable soliton can become stable, and it maps the stable and unstable regions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The κ≥2 stabilisation claim rests on an unverified completeness of the discrete spectrum; Eq. (35) is only sufficient if no other unstable eigenvalues exist.","rationale":"The analytic portions are solid: the exact stationary solutions are correctly derived, the instability of Ψ− follows from a clean variational argument, and the proof of |Λ_i|>|Λ_r| for κ<2 is a legitimate application of Weinstein's lemma. The central novel claim, however, is numerical. Eq. (35) is a necessary condition for the tracked quadruplet, but stability of the full linear problem requires that every eigenvalue satisfies ν≤ρ̃. The paper only tracks the quadruplet and the continuum and does not demonstrate that no additional discrete eigenvalue crosses into the unstable half-plane inside the claimed window. Direct numerical time evolution is the standard check that would settle both linear completeness and the physical relevance of the stability region. The reader's weakest assumption identified exactly this issue; I agree with that assessment. The reader's conditional verdict already expresses the appropriate recommendation, so no adjustment is needed.","tokens_in":15694,"tokens_out":6378,"duration_ms":58629,"concrete_test":"Directly time-integrate Eq. (1) from Ψ+ plus a small generic perturbation at an interior point of the claimed stable region, e.g., for κ=2.5 at the midpoint between the two branches at a fixed ρ in the lower panel of Fig. 5, for times t ≳ 10^4; repeat with a perturbation aligned to the previously unstable real mode. If |ψ| stays bounded and perturbations decay for both choices, the central claim is supported; if either run grows, the claimed stability region is not real. A convergence check of the same run at L=100, N=2001 would separate physical drift from discretization error.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The principal claim in Sec. 5 (and the abstract) that for κ≥2 the unstable soliton becomes stable for ε>ε̃_c depends on more than the analytic derivations: the stability criterion (35) is applied as a necessary-and-sufficient condition for the single tracked quadruplet, but nothing in the paper rules out a second discrete eigenvalue whose growth rate ν exceeds ρ̃ in the claimed stable window. Section 4.2 reports only the collision of a real pair with the zero modes and the subsequent quadruplet, and Fig. 5 draws the stable region between the two branches from that quadruplet alone. The numerical basis is L=50, N=501, with no reported convergence table for ε̃_c, no error bounds, and no direct integration of Eq. (1). Because for κ≥2 the unforced soliton is already unstable through a real positive mode, the stabilisation mechanism is a delicate real-axis collision; a missed eigenvalue or a discretization-induced shift of ε̃_c would leave the claimed region either empty or incorrect. Thus the abstract's 'establishes' overstates the evidence unless this completeness and accuracy gap is closed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":15817,"tokens_out":3831,"duration_ms":36645,"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":[{"comment":"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.","section":"Abstract and Sec. 4.2"},{"comment":"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.","section":"Sec. 4.2, Figs. 4 and 5"},{"comment":"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.","section":"Sec. 4.2 and Sec. 5(iv)"}],"minor_comments":[{"comment":"The phrase \"respectively\" in Eq. (4) is ambiguous; the text should clarify which parameters correspond to which equation.","section":"Sec. 1"},{"comment":"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.","section":"Sec. 3, item 1"},{"comment":"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.","section":"Sec. 4.2"},{"comment":"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.","section":"Appendix B, Eq. (61)"},{"comment":"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.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the analytic parts are solid. My main concern is the overstatement of the κ ≥ 2 result in the abstract and conclusions relative to the numerical evidence. If the authors add the requested spectral completeness check, convergence study, and ideally direct time evolution, the claim could become publishable. I recommend major revision rather than rejection because the underlying mechanism is plausible and the analytic framework is sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Main take: this is a careful, mostly analytic extension of the 1991 Barashenkov–Bogdan–Korobov stability analysis from κ=1 to arbitrary κ>0. The genuinely new piece is the κ≥2 claim that damping plus parametric drive can stabilize the otherwise unstable NLS soliton via a real-axis collision followed by oscillatory stability. That claim is interesting and plausible, but it rests on numerical eigenvalue computations whose completeness and accuracy are not yet demonstrated. The abstract's 'establishes' overstates the evidence; with a convergence study and ideally a direct time integration, the result would be convincing.\n\nWhat is good: the stationary solutions (16)–(17) are exact and their existence conditions are independent of κ; the dissipationless eigenvalue reduction in Sec. 3 is clean; the proof that Ψ− is always unstable (Appendix A) is a solid variational argument; and Appendix B uses Weinstein's lemma to prove |Λi|>|Λr| for all κ<2 when a quadruplet exists. That last bit is a real analytical result, and it is new as far as I can tell. The κ=1 case is correctly recovered. The stability-curve computation via the ε parametrization is elegant, and the paper lays out the logic carefully. Figures 1–5 support the qualitative picture.\n\nThe soft spot is the κ≥2 section. The stabilization claim hinges on criterion (35) being necessary and sufficient for the full discrete spectrum, but the paper tracks one quadruplet and the continuum, and does not demonstrate that no other discrete mode can become unstable inside the claimed stable window between the two branches in Fig. 5. The numerics use L=50, N=501, with larger N mentioned but no convergence table for ε̃c; there are no error bounds, and no direct integration of Eq. (1) to confirm dynamical stability. Given that at ε=0 the soliton has a real positive eigenvalue, the mechanism is a delicate collision; a missed eigenvalue or a discretization shift of ε̃c could change the stable region. The stress-test concern lands here. This is not a fatal flaw—the numerics plausibly capture the right eigenvalues—but the abstract's strong wording should be softened to 'numerical evidence indicates' or similar.\n\nThe citation pattern is fair: Ref. [1] is credited for κ=1, and the relevant NLS stability literature is present. No sign of self-citation padding.\n\nWho this is for: people working on parametrically driven/damped NLS, non-Kerr nonlinearities, and stability of envelope solitons. It deserves a serious referee; a referee should ask for a convergence check of the κ≥2 eigenvalues, a statement about completeness of the tracked spectrum, and ideally one direct time evolution. My recommendation is conditional accept with moderate revision.","headline":"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.","tokens_in":16449,"tokens_out":2163,"would_cite":true,"duration_ms":19699,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35B35","37K40","34B24"],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["parametrically driven damped nonlinear Schrödinger equation","oscillatory stability","oscillatory instability","Sturm-Liouville eigenvalue problem","soliton stability","non-Kerr nonlinearity","eigenvalue quadruplet","blow-up stabilization"],"falsifier":"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$.","tokens_in":15427,"feed_emoji":"🌊","tokens_out":6999,"duration_ms":57483,"temperature":0.7,"pith_summary":"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.","feed_headline":"Driven damping stabilizes solitons that normally blow up","feed_subtitle":"For nonlinearity exponent $\\kappa\\ge 2$, a stable window opens beyond a critical drive amplitude.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the $\\kappa=1$ stability diagram and the exceptional change of variables that this paper extends to arbitrary $\\kappa$.","marker":"[1]"},{"why":"Establishes the connection between the parametrically driven, damped NLS equation and breathers of damped nonlinear Klein-Gordon systems, which motivates the stability transfer.","marker":"[2]"},{"why":"Documents the well-known blow-up instability of the NLS soliton for $\\kappa\\ge 2$, the instability that this paper shows can be overcome.","marker":"[17]"},{"why":"Identifies oscillatory instability as the collision of two internal modes, the mechanism generalized here to oscillatory stability.","marker":"[28]"},{"why":"Provides the $\\kappa=1$ proof of $\\Psi_-$ instability that Appendix A extends to all $\\kappa>0$.","marker":"[36]"},{"why":"Delivers the operator-theoretic result used to prove that $L_0$ is strictly positive definite for $\\varepsilon<0$, a step in the $\\Psi_-$ instability proof.","marker":"[38]"},{"why":"Supplies the positivity result (Proposition B1) used to prove $|\\Lambda_i|>|\\Lambda_r|$ for $\\kappa<2$.","marker":"[39]"},{"why":"Provides the fixed Newton method and discretization approach used for the numerical eigenvalue computations that underpin the $\\kappa\\ge 2$ result.","marker":"[31]"}],"fun_headline_variants":["Damping and drive tame soliton blow-up for strong nonlinearity","Stable window opens for nonlinear solitons at high exponent","Oscillatory stability transition in driven damped NLS","Parametric drive flips soliton instability to stability","κ≥2 unlocks stable solitons under parametric drive"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Damping and drive tame soliton blow-up for strong nonlinearity","Stable window opens for nonlinear solitons at high exponent","Oscillatory stability transition in driven damped NLS","Parametric drive flips soliton instability to stability","κ≥2 unlocks stable solitons under parametric drive"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000637,"raw_usage":{"total_tokens":3001,"prompt_tokens":1078,"completion_tokens":1923,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":694,"completion_tokens_details":{"reasoning_tokens":1849}},"tokens_in":694,"tokens_out":1923,"duration_ms":12554,"temperature":1.0,"reasoning_tokens":1849,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:54:37.887234+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the $\\kappa=1$ stability diagram and the exceptional change of variables that this paper extends to arbitrary $\\kappa$."},{"cited_title":"Bondila, I","cited_arxiv_id":null,"evidence_quote":"Establishes the connection between the parametrically driven, damped NLS equation and breathers of damped nonlinear Klein-Gordon systems, which motivates the stability transfer."},{"cited_title":"Sulem and C","cited_arxiv_id":null,"evidence_quote":"Documents the well-known blow-up instability of the NLS soliton for $\\kappa\\ge 2$, the instability that this paper shows can be overcome."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies oscillatory instability as the collision of two internal modes, the mechanism generalized here to oscillatory stability."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the $\\kappa=1$ proof of $\\Psi_-$ instability that Appendix A extends to all $\\kappa>0$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Delivers the operator-theoretic result used to prove that $L_0$ is strictly positive definite for $\\varepsilon<0$, a step in the $\\Psi_-$ instability proof."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the positivity result (Proposition B1) used to prove $|\\Lambda_i|>|\\Lambda_r|$ for $\\kappa<2$."},{"cited_title":"Carretero-González, D","cited_arxiv_id":null,"evidence_quote":"Provides the fixed Newton method and discretization approach used for the numerical eigenvalue computations that underpin the $\\kappa\\ge 2$ result."}],"review_version":1}