{"id":"52c5fad1-633d-414a-924f-e40f10f26ef0","arxiv_id":"2607.17446","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In the fractional parametrically driven damped NLS equation, quiescent ψ+ solitons are stable below a Lévy-index-dependent threshold, while lossless moving solitons are stable only in intermediate or high velocity windows.","lead":"This paper maps which light pulses — solitons — can survive in a laser-cavity model that combines parametric pumping, loss, and an unusual 'fractional' form of diffraction. It reports that stationary solitons are stable only below a drive threshold that depends on the fractional order, while moving solitons exist and become stable only in specific speed bands.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table I’s α=2 entries for γ≤0.3 contradict Eq. (7); the numerical stability criterion misclassifies near-marginal solitons, so the quantitative thresholds in Tables I–III are unreliable.","rationale":"The reader’s weakest assumption—that the eigenvalue solver misclassifies near-marginal solitons in the weak-loss regime—is precisely the load-bearing concern. The paper itself supplies a decisive benchmark: the exact α=2 result, Eq. (7), directly contradicts Table I at small γ. This is an internal inconsistency, not a matter of external consensus. The qualitative stability structure (ψ+ stable, ψ− unstable at rest; high-velocity stabilization of ψ−) is plausible and partly anchored by the exact limits, but the quantitative α-dependence of the thresholds—the paper’s main new content—cannot be trusted until the numerical criterion is re-verified. The reader’s conditional verdict is appropriate; no change is needed, but the concern is real and should be resolved by the proposed test.","tokens_in":14673,"tokens_out":4514,"duration_ms":44111,"concrete_test":"Using the exact α=2 solution in Eqs. (4)–(6) at γ=0.1, h=0.5 (and h=1.0), compute the eigenvalue problem (22) with the same 512-mode Fourier collocation and a perturbation subspace that excludes the translational mode. Exact theory (Eq. (7)) predicts Re(δ)≤0; any positive eigenvalue is spurious. Then directly integrate Eq. (11) for these parameters; if the soliton persists for t ≥ 10^3 with no growth, Table I’s h_c=0.12 is falsified. This isolates the numerical criterion from the fractional-α results.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s novel quantitative content is the α-dependence of the stability threshold h_c (Table I) and of the traveling-soliton velocity windows (Tables II–III). All are generated by the eigenvalue criterion of Eq. (22) solved with 512-mode Fourier collocation. That criterion is demonstrably unreliable in the weak-loss limit. The paper itself quotes the exact α=2 result, Eq. (7): ψ+ is stable for h < √(1+γ²). Yet Table I reports h_c=0.06 at γ=0 and 0.12 at γ=0.1, whereas the exact values are 1.0 and ≈1.005. Entries become consistent only for γ≥0.4. This pattern indicates that near the conservative limit eigenvalues lie close to the imaginary axis, so numerical truncation/noise produces spurious Re(δ)>0 and truly stable solitons are flagged unstable. No convergence analysis, grid-resolution dependence, or eigenvalue tolerance is reported, so there is no way to judge whether the α<2 thresholds suffer the same artifact. The qualitative structure—ψ+ stable at rest, ψ− stabilized at high velocity—may survive, but the α-dependent numbers, which are the advertised contribution, are unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies solitons of the one-dimensional fractional parametrically driven damped nonlinear Schrödinger equation (Eq. (1)). For quiescent solitons it derives the zero-background stability condition h ≤ sqrt(1+γ²) (Eq. (14)), the norm condition h ≥ γ (Eq. (20)), and numerically constructs two soliton families, ψ+ and ψ−, reporting the stability threshold h_c of ψ+ in Table I. For the lossless case γ=0 it computes traveling solitons and reports stable velocity intervals in Tables II and III. The paper also presents collision simulations. The central quantitative claims are the α-dependence of h_c and of the velocity windows.","tokens_in":14884,"tokens_out":3599,"duration_ms":40242,"significance":"The analytical parts—background dispersion relations and the norm condition—are correct and provide useful constraints for this model. The paper is the first systematic study of the effect of fractional diffraction in the parametrically driven damped NLS setting, and the qualitative picture (ψ+ stable at rest, ψ− stabilized at high velocity) is plausible. However, the numerical stability results are not validated against the exact α=2 limit and are demonstrably wrong there; since Tables I–III are produced by the same eigenvalue criterion, the quantitative content of the paper is currently unsupported. With a corrected and benchmarked stability solver, the work could become a valuable contribution.","major_comments":[{"comment":"Table I contradicts the paper's own exact α=2 result. For γ=0.1, α=2.0, Table I reports h_c=0.12, whereas Eq. (7) gives h_c=√(1+γ²)=1.005. For γ=0, Table I gives h_c=0.06, while the exact statement in §II says ψ+ is stable for h<1. The discrepancy indicates that the eigenvalue solver (Eq. (22), 512-mode Fourier collocation) misclassifies near-marginal solitons as unstable at small γ, likely because eigenvalues lie close to the imaginary axis and no tolerances or convergence checks are reported. Since Table I is the central quantitative result for the α-dependence of the stability threshold, this is a load-bearing error.","section":"§III, Eq. (7) and Table I"},{"comment":"The traveling-soliton stability windows in Tables II and III are computed with the same eigenvalue criterion, and specifically in the conservative limit γ=0 where the solver has just been shown to fail (Table I at γ=0 gives h_c=0.06 instead of 1 for α=2). No benchmark against the known α=2 traveling-soliton results of Ref. [42] is provided, and no grid-resolution or tolerance study is reported. The velocity intervals are therefore unreliable as quantitative predictions, and the advertised α-dependence of the stability windows is unsupported.","section":"§IV, Tables II–III"},{"comment":"The manuscript does not report any convergence analysis, dependence on the number of Fourier modes, or eigenvalue acceptance thresholds for the stability criterion. Given the demonstrated failure at the exactly solvable α=2 limit, the α<2 entries in Table I and Tables II–III cannot be trusted without such validation. At minimum, the authors should rerun the stability calculations with an improved eigensolver and demonstrate that the exact α=2 thresholds are reproduced to a stated accuracy.","section":"§III, numerical method"}],"minor_comments":[{"comment":"Grammatical error: 'Collision between moving solitons are considered too' should be 'Collisions between moving solitons are considered too'.","section":"Abstract"},{"comment":"Typo: 'collision' is misspelled as 'ollision' in panel (a) description.","section":"§V, caption of Fig. 7"},{"comment":"The introduction cites many relevant works, but the transition from conservative fractional NLS to the driven-damped model would benefit from a brief explicit statement of how the parametric drive is realized experimentally in the fractional cavity context.","section":"§I, references"},{"comment":"The phrase 'down branches' should be 'lower branches'.","section":"§III, Fig. 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The numerical stability issue is severe enough that the current quantitative claims cannot be accepted. The authors need to replace the stability computations with a validated solver that reproduces the exact α=2 results, and then re-derive Tables I–III. If the corrected results still show the same qualitative structure, the paper could be suitable for publication; as it stands, the central advertised contribution is unreliable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper introduces the fractional version of the parametrically driven damped nonlinear Schrödinger equation and maps out its quiescent and traveling solitons. The model is new and physically motivated—fractional GVD in laser cavities is an active experimental direction—and the analytic background-stability analysis (Eqs. 12–15) and the norm bound h≥γ are correct and clearly presented. The qualitative picture—ψ+ stable in a window, ψ− always unstable at rest, ψ− stabilized at high velocity, α shrinking the windows—is plausible and, at α=2 and larger γ, the stability thresholds in Table I match the known exact condition h_c = sqrt(1+γ^2). The collision studies are a nice extra.\n\nThe problem is the quantitative core. Table I is supposed to give the critical h_c above which ψ+ loses stability. For α=2, the exact condition is h_c = sqrt(1+γ^2), which the paper itself quotes as Eq. (7). Yet Table I lists h_c = 0.06 at γ=0 and 0.12 at γ=0.1, where the exact values are 1.00 and 1.005. The entries do converge to the exact curve only for γ≥0.4. This is not a minor rounding issue; it is off by an order of magnitude in the weak-loss regime. The likely cause is the eigenvalue solver: with 512 Fourier modes and no reported convergence checks or tolerances, eigenvalues near the imaginary axis in the conservative limit acquire spurious positive real parts, so stable solitons are flagged unstable. The same pipeline produces Tables II and III for traveling solitons at γ=0, where the issue is even more acute.\n\nSo the paper's advertised contribution—the α- and γ-dependence of the stability thresholds and the velocity windows—is unsupported as it stands. The qualitative structure may survive re-analysis, and the large-γ part of Table I is consistent, but the tables need to be recomputed with a properly validated solver and checked against the α=2 exact results across the full γ range. The authors should also report discretization and tolerance studies, and probably make the code available.\n\nThis is a serious paper from a credible group, and I'd send it to peer review—a good referee will catch this and the authors can fix it. But I wouldn't cite the numerical maps in their current form.\n\nBest,","headline":"The fractional PDDNLSE model and its qualitative soliton picture are worth knowing about, but Table I's α=2 entries contradict the exact stability condition for small γ, so the quantitative stability maps are not trustworthy until fixed.","tokens_in":15518,"tokens_out":6396,"would_cite":false,"duration_ms":57578,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35R11","37K40"],"pacs":[],"model":"deepseek-v4-flash","headline":"In a fractional driven-damped Schrödinger equation, only one standing soliton branch is stable, and motion can stabilize the other.","keywords":["fractional nonlinear Schrödinger equation","Riesz derivative","Lévy index","parametric drive","dissipative solitons","traveling solitons","linear stability","laser cavity"],"falsifier":"Compare the numerically computed stability threshold h_c at α=2 with the exact formula sqrt(1+γ²). For γ=0 and γ=0.1 the paper's Table I lists h_c=0.06 and 0.12, whereas the exact values are 1 and about 1.005; this large discrepancy is directly checkable and would settle whether the numerical stability criterion is trustworthy.","tokens_in":14461,"feed_emoji":"🌊","tokens_out":4678,"duration_ms":52736,"temperature":0.7,"pith_summary":"This paper studies a one-dimensional Schrödinger model in which ordinary diffraction is replaced by a fractional Riesz derivative, while cubic self-focusing is balanced by a parametric drive and linear loss. It claims that this equation supports two families of standing solitons: a large-amplitude branch that is stable only for a finite window of driving strength, and a small-amplitude branch that is always unstable when at rest. In the lossless version, moving solitons exist below a maximum velocity, and the unstable resting branch becomes stable at high speeds. The paper argues that the fractional order, the Lévy index, strongly controls both the existence conditions and the stability thresholds, effectively expanding the variety of nonlinear modes available in fractional media.","feed_headline":"Motion stabilizes unstable solitons in fractional driven waves","feed_subtitle":"A fractional laser-cavity model predicts stable solitons whose existence windows shrink as the Lévy index drops.","key_machinery":"The central object is Eq. (1), the fractional parametrically driven damped nonlinear Schrödinger equation: iψ_t − iV ψ_ξ + (iγ+ω)ψ − (−∂²_ξ)^(α/2)ψ + 2|ψ|²ψ = hψ*. The Riesz fractional derivative of order α, with 1<α≤2, replaces the usual Laplacian and introduces nonlocal diffraction. Because the parametric drive and the fractional derivative break Galilean invariance, traveling solitons must be sought in a co-moving frame. Stability is decided by linearizing around a soliton and solving the eigenvalue problem of Eq. (22): if any eigenvalue has positive real part, the soliton is unstable. In the non-fractional limit α=2, exact soliton solutions and an exact stability bound are known, and tho","core_discovery":"The central claim is that the fractional parametrically driven damped nonlinear Schrödinger equation admits quiescent solitons of two species, ψ+ and ψ−, which exist above the parametric-gain threshold h=γ. According to the authors' numerical analysis, ψ− is always unstable, while ψ+ is stable for h between γ and a critical value h_c(α,γ) that shrinks as the Lévy index α decreases and vanishes before α reaches 1. In the conservative limit γ=0, the same equation supports traveling solitons with velocities below a maximum; ψ+ solitons are stable in intermediate velocity bands, whereas ψ− solitons, unstable when stationary, become stable at high velocities. The fractional diffraction operator i","pith_inferences":["The reported α-dependence of h_c may be overstated if the numerical stability criterion is unreliable in the weakly damped limit; a testable next step is to recompute the thresholds with an independent eigensolver.","Because exact results are known at α=2, one can directly validate the numerical tables by comparing h_c with sqrt(1+γ²); the paper's Table I appears inconsistent with that benchmark at small γ.","The pattern that the unstable branch becomes stable only when moving suggests a general mechanism in driven-damped systems: translation can push growth rates across the imaginary axis, which may also apply to bound states of driven solitons.","In an experimental cavity, the predicted velocity-dependent stability could be probed by launching pulses with controlled relative velocities and observing whether fast small-amplitude pulses survive collisions while slow ones diffract."],"forward_implications":["A laser cavity with emulated fractional group-velocity dispersion should produce standing solitons that are robust only for a finite pumping interval, with that interval closing as the Lévy index approaches 1.","Traveling solitons exist despite the breaking of Galilean invariance, but only below a maximum velocity that decreases with smaller fractional order.","Motion itself can act as a stabilizer: the otherwise unstable small-amplitude branch becomes stable at high velocities.","Collisions between high-velocity stable solitons of the small-amplitude branch are quasi-elastic, while large-amplitude soliton collisions display attraction caused by nonlocal fractional coupling.","Lowering the fractional order shrinks the stability windows in both driving strength and velocity, offering a control knob for soliton existence."],"fun_headline_variants":["Moving solitons gain stability in fractional driven waves","Fractional drive stabilizes solitons in motion","Lower Levy index narrows soliton existence windows","High-velocity solitons stabilize in fractional drive"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The paper's stability thresholds rest on a numerical eigenvalue calculation, solved by Fourier collocation with 512 modes, that must correctly separate stable from unstable solitons; in the weakly damped, non-fractional limit that calculation disagrees with the exact analytical result, so the reported thresholds may not be reliable.","fun_headline_variants_meta":{"raw":{"variants":["Moving solitons gain stability in fractional driven waves","Fractional drive stabilizes solitons in motion","Lower Levy index narrows soliton existence windows","High-velocity solitons stabilize in fractional drive"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000662,"raw_usage":{"total_tokens":2825,"prompt_tokens":673,"completion_tokens":2152,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":417,"completion_tokens_details":{"reasoning_tokens":2090}},"tokens_in":417,"tokens_out":2152,"duration_ms":18558,"temperature":1.0,"reasoning_tokens":2090,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T18:03:50.118220+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the numerically computed stability threshold h_c at α=2 with the exact formula sqrt(1+γ²). For γ=0 and γ=0.1 the paper's Table I lists h_c=0.06 and 0.12, whereas the exact values are 1 and about 1.005; this large discrepancy is directly checkable and would settle whether the numerical stability criterion is trustworthy.","supporting_citations":[],"review_version":1}