{"id":"c5bc6e96-be77-4278-b63d-b3c1f9791719","arxiv_id":"2504.13299","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Exact solitons are constructed for a generalized two-parameter nonlinear Dirac model, with E/Q independent of the coupling g and stability predicted by the Vakhitov-Kolokolov criterion for kappa < 2.","lead":"This paper derives exact stationary solitary wave shapes for a two-parameter family of 1+1 dimensional nonlinear Dirac equations mixing scalar-scalar and vector-vector interactions. The explicit profiles and their stability map give analytic benchmarks for nonlinear Dirac models used in photonic lattices, Bose-Einstein condensates, and conducting polymers.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Solitary-wave construction appears sound, but the stability boundaries are not established: the Vakhitov-Kolokolov condition dQ/dω<0 is only necessary, not sufficient, for linear stability of nonlinear Dirac waves.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the VK criterion is treated as sufficient for nonlinear Dirac stability without proof. My independent check of the algebraic derivation confirms that the exact-solution construction itself is sound: Eq. (19) and Eq. (20) are mutually consistent, Eq. (24) gives a positive localized R for all ω in the stated open interval, and the scaling argument for E/Q independent of g is valid. The paper's principal contribution therefore survives. What does not survive is the stability overreach: the conclusions claim stable regimes based solely on a necessary condition, and the paper even concedes that the κ>2 instability region awaits numerical verification. A revised version should either add a spectral stability computation or explicitly soften the stability claims to 'VK-consistent' rather than 'stable'. The reader's secondary point about Eq. (32) also lands: the sign/coefficient in Eq. (32) does not reproduce the claimed threshold in the p→∞ limit, although the final threshold quoted in Eq. (33) is correct when re-derived from Eq. (24). Because the central existence result is solid and the stability issue is a missing justification rather than a demonstrated contradiction, the appropriate outcome remains conditional acceptance, matching the reader's verdict.","tokens_in":13400,"tokens_out":25753,"duration_ms":227067,"concrete_test":"Implement a spectral linear-stability solver for the linearization of Eq. (2) about the exact profile (24), and scan a small grid in the predicted stable region, including at least κ=1, p=2, ω/m=0.5 and κ=1.5, p=5, ω/m=0.8, where dQ/dω<0. If any eigenvalue with positive real part is found despite dQ/dω<0, the Vakhitov-Kolokolov-based stability boundary is invalid; if none is found in a well-resolved scan, the concern is substantially weakened. This directly supplies the missing spectral check from Sec. IV.F.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The existence construction is internally consistent: Eq. (19) follows from H=0, Eq. (20) integrates to θ=arctan(α tanh(κβx)), and R from Eq. (24) is real and localized whenever ω/m>1/p^{1/(κ+1)} because cos^{κ+1}(2θ) stays above (ω/m)^{κ+1}>1/p. So the central claim of exact solitary waves and of E/Q independent of g is not where the argument fails. The load-bearing weakness is in Sec. IV.F-G: the statements that solutions are stable for κ<2 and in a small window beyond κ>2 rely entirely on the Vakhitov-Kolokolov condition (Eq. (56)). For nonlinear Dirac equations this condition is known to be necessary but not sufficient for spectral stability; one must also analyze the spectrum of the linearized operator, excluding eigenvalues with positive real part and checking the relevant Krein-signature/coercivity conditions. None of that analysis is supplied. Moreover, Sec. IV.G explicitly says the predicted instability for κ>2 'will have to be verified by numerics' and phrases the conclusion as 'We believe', so the instability region is also not established. The concluding sentence that VK analysis 'has shown' stability is therefore an overstatement. A secondary internal error, less central but worth correcting, is that Eq. (32) contradicts Eq. (33) and the p→∞ limit; the correct α_c² from expanding Eq. (24) is (p−1)/(p(2κ+1)+1), not (p−1)κ/[2p+(p+1)κ].","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs exact stationary solitary wave solutions of a two-parameter family of nonlinear Dirac equations in 1+1 dimensions, with Lagrangian (7): a scalar-scalar interaction of power κ+1 and a vector-vector interaction of power (κ+1)/2 with relative weight set by p>1. For the ansatz Ψ(x,t)=Φ(x)e^{-iωt}, the authors derive the profile from Eqs. (19)-(24), obtaining localized real solutions under the frequency condition 1/p^{1/(κ+1)} < ω/m < 1. They then compute charge and energy, show E/Q is independent of the coupling g, discuss single vs double hump transitions, use the Vakhitov-Kolokolov criterion to claim stability for κ<2 and a small window above κ=2, and finally give a nonrelativistic reduction to a modified nonlinear Schrödinger equation with a Derrick-theorem stability argument.","tokens_in":13686,"tokens_out":15318,"duration_ms":131571,"significance":"If the existence construction is correct, this is a valuable extension of the ABS model: it gives explicit soliton profiles over the full (κ,p) parameter region and identifies a sharp frequency threshold, with E/Q independent of g as a concrete and falsifiable prediction. The derivation from the stated Lagrangian is largely analytic and self-contained, and the κ=1, p=2 ABS limit and the p→∞ scalar-scalar limit are natural consistency checks. The claimed stability results, however, go beyond what is proved in the paper, and several of the formulas underlying the charge, energy, and hump-transition calculations are internally inconsistent. The core existence result may survive, but the quantitative and stability claims require substantial reworking.","major_comments":[{"comment":"The integral I(ω,κ,p) is inconsistent with Eq. (27) and the change of variables y=tanh(κβx). Equation (27) gives R² ∝ (1−y²)^{1/κ}, and since dx = dy/[κβ(1−y²)], the Jacobian leads to an integrand proportional to (1−y²)^{1/κ−1}. The manuscript instead writes (1−y²)^{1−1/κ}, which differs from the correct exponent at every κ except κ=1. This error propagates into Q(ω,κ,p), E(ω,κ,p), dQ/dω, and therefore into the Vakhitov-Kolokolov plots and stability diagrams in Secs. IV.C-IV.F. The formulas and the numerical results derived from them need to be rederived and recomputed.","section":"Sec. IV.B, Eq. (36b)"},{"comment":"The double/single-hump threshold quoted in Eqs. (33)-(34) does not follow from Eq. (32). Setting the bracket in Eq. (32) to zero gives α² = κ(p−1)/[2p+(p+1)κ], which is equivalent to ω/m = (p+κ)/[p(κ+1)], not the displayed (pκ+1)/[p(κ+1)]. The two expressions agree only for κ=1. In addition, Eq. (29) uses (1−y²)^{2/κ} whereas Eq. (27) gives (1−y²)^{1/κ}; this inconsistency appears to be the source of the incorrect threshold. The claimed transition curve in the (κ,p) plane is therefore not established by the derivation as written.","section":"Sec. IV.A, Eqs. (32)-(34)"},{"comment":"The central stability claim is supported only by the Vakhitov-Kolokolov condition dQ/dω<0 (Eq. (56)). For nonlinear Dirac equations this condition is necessary but not sufficient for spectral stability; one also needs an analysis of the spectrum of the linearized operator, including Krein-signature or coercivity conditions. No such analysis is supplied. Furthermore, Sec. IV.G explicitly says the instability for κ>2 \"will have to be verified by numerics\" and phrases the prediction as \"We believe\", yet the final section states \"we have further shown that these solutions are not only stable for κ<2 but even for a small window beyond κ>2\". This is an overstatement. The stability statements should be either downgraded to conjectures or supported by a full linear-stability analysis.","section":"Secs. IV.F, IV.G, and Conclusions"}],"minor_comments":[{"comment":"Equation (15) contains an unbalanced parenthesis and the first g_s appears without its square; the surrounding notation alternates between g_s, g_v and g_s^2, g_v^2.","section":"Sec. III, Eq. (15)"},{"comment":"The statement in Eq. (74) that \"H2 is positive and small relative to H1\" is used to conclude stability from Derrick's theorem, but no estimate of the relative size of H1 and H2 in the regime of validity of the nonrelativistic reduction is provided.","section":"Sec. V, Eq. (74)"},{"comment":"The paper presents many numerical plots of Q, dQ/dω, and E/Q, but does not state the numerical accuracy of the integrals used; a brief description of the quadrature method and error estimates would improve reproducibility.","section":"Sec. IV, Figs. 4-10"},{"comment":"References [26] and [27] are listed but not discussed in the text; please either cite them in the stability discussion or remove them.","section":"Bibliography"}],"recommendation":"major_revision","confidential_remarks":"The paper is an honest analytic construction but the internal algebraic inconsistencies in Sec. IV are load-bearing: they affect the quantitative charge and energy formulas and the double/single-hump threshold. The stability section overreaches relative to what the VK criterion can establish for nonlinear Dirac equations. The authors should be asked to correct the algebra, recompute the affected plots, and resubmit with either a complete linear-stability analysis or carefully hedged stability claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the central derivation is in good shape. The authors construct exact stationary solutions for the (kappa,p) family, get a clean existence domain, and prove E/Q independence of g. The stress-test concern about the existence formulas did not land: Eqs. (19)-(24) check out and the localization condition is as stated. Where the paper is shaky is stability.\n\nWhat is actually new: the two-parameter family itself, the exact profiles, the hump-to-single-hump transition, and the E/Q independence. For a specialist in nonlinear Dirac equations these are useful closed-form benchmarks. The paper is a legitimate extension of the ABS model and the earlier S-S/V-V work, and the self-citations to [16] and [17] are appropriate, not circular.\n\nThere is a genuine internal inconsistency: Eq. (32) does not match Eq. (33) or the p->infinity limit. Expanding Eq. (24) gives the threshold alpha_c^2 = (p-1)/(p(2kappa+1)+1), not the coefficient in Eq. (32). This is minor but should be corrected.\n\nThe larger problem is stability. Section IV.F invokes the Vakhitov-Kolokolov condition dQ/domega < 0 as if it were sufficient for linear stability of nonlinear Dirac solitary waves. It is not. For NLDEs, VK is necessary but not sufficient; one must also analyze the spectrum of the linearized operator. None of that is supplied. Section IV.G then explicitly says the predicted instability for kappa > 2 \"will have to be verified by numerics\" and phrases the conclusion as \"We believe.\" So the concluding statement that VK analysis \"has shown\" stability is an overstatement. The kappa < 2 stability claim may be true, but it is not established by the analysis in this paper.\n\nThe non-relativistic reduction is plausible and the Derrick-type argument is fine as a heuristic, but it uses an approximate Hamiltonian and should be labeled as such.\n\nBottom line: the exact-soliton part is valuable and likely correct; the stability part needs either real spectral analysis or substantially softened claims. The paper deserves a serious referee, but the referee should ask for those changes before publication.\n\nRecommendation: send to peer review.","headline":"Sound exact-soliton construction for a two-parameter NLDE family; the stability claims outrun the evidence.","tokens_in":14279,"tokens_out":1838,"would_cite":true,"duration_ms":19287,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q51","35Q55","35Q41","37K40"],"pacs":["03.65.Pm","05.45.Yv"],"model":"deepseek-v4-flash","headline":"This paper constructs exact stationary solitary-wave solutions for every κ>0, p>1 in a generalized ABS nonlinear Dirac model, with E/Q independent of the coupling constant.","keywords":["exact solitary waves","nonlinear Dirac equation","generalized ABS model","scalar-scalar interaction","vector-vector interaction","Vakhitov-Kolokolov criterion","nonrelativistic reduction","modified nonlinear Schrödinger equation"],"falsifier":"Run direct numerical time evolution of the exact solitons for $\\kappa>2$ at frequencies where $dQ/d\\omega>0$ (the Vakhitov--Kolokolov unstable region) and where $E/Q>1$; linear stability would be settled by computing the spectrum of the linearized Dirac operator around $\\Phi(x)$ and looking for eigenvalues with positive real part. If the solitons remain stable where $dQ/d\\omega>0$, the Vakhitov--Kolokolov criterion is not sufficient for this model.","tokens_in":13155,"feed_emoji":"🌊","tokens_out":10851,"duration_ms":88655,"temperature":0.7,"pith_summary":"This paper constructs exact stationary solitary-wave solutions of a 1+1-dimensional nonlinear Dirac equation whose self-interaction mixes a scalar-scalar term of power $\\kappa+1$ with a vector-vector term of power $(\\kappa+1)/2$, with an adjustable relative coefficient set by $p>1$. The authors show that for every $\\kappa>0$ and $p>1$ a localized rest-frame solution $\\Psi(x,t)=\\Phi(x)e^{-i\\omega t}$ exists whenever $1/p^{1/(\\kappa+1)}<\\omega/m<1$, filling the entire allowed parameter plane. The exact profile is given in closed form, the charge and energy are finite, and the ratio $E/Q$ is independent of the coupling constant $g$. Stability is addressed through the Vakhitov--Kolokolov criterion, giving stability for $\\kappa<2$ and a small window beyond, with an unstable region predicted for larger $\\kappa$. A non-relativistic reduction yields a modified nonlinear Schrödinger equation whose solitons are stable for $\\kappa<2$.","feed_headline":"Exact solitary waves fill the κ–p plane of a generalized Dirac model","feed_subtitle":"A tunable vector interaction yields explicit solitons whose energy-to-charge ratio ignores the coupling strength","key_machinery":"The central object is the radial decomposition $\\Phi(x)=R(x)(\\cos\\theta(x), \\sin\\theta(x))^T$ in the gamma-matrix representation $\\gamma^0=\\sigma_3$, $\\gamma^1=i\\sigma_1$. Substituting it into the equation of motion turns the stationary problem into a first-order equation $d\\theta/dx=\\kappa[m\\cos(2\\theta)-\\omega]$, solved by $\\tan\\theta=\\alpha\\tanh(\\kappa\\beta x)$, together with an algebraic equation for $R^{2\\kappa}$. This factorization reduces the entire two-parameter family to quadratures: $Q$ and $H_2$ become one-dimensional integrals over $y=\\tanh(\\kappa\\beta x)$, from which $E/Q$ and the charge-frequency relation are obtained, and the Vakhitov--Kolokolov derivative $dQ/d\\omega$ can be evaluated. The same radial machinery drives the non-relativistic reduction, replacing the interaction by the effective coupling $g^2(p-1)/p$ in the modified nonlinear Schrödinger equation.","core_discovery":"On the paper's own terms, the central discovery is an exact construction covering the whole admissible parameter plane: for the generalized ABS model with interaction Lagrangian $L_I = \\frac{g^2}{(\\kappa+1)}(\\bar{\\psi}\\psi)^{\\kappa+1} - \\frac{g^2}{p(\\kappa+1)}[\\bar{\\psi}\\gamma_\\mu\\psi\\,\\bar{\\psi}\\gamma^\\mu\\psi]^{(\\kappa+1)/2}$, every pair $\\kappa>0$, $p>1$ admits localized solitary waves $\\Psi(x,t)=\\Phi(x)e^{-i\\omega t}$ with rest-frame frequency in the window $1/p^{1/(\\kappa+1)}<\\omega/m<1$. The spinor profile is expressed through $\\theta(x)=\\arctan[\\alpha\\tanh(\\kappa\\beta x)]$ and $R^{2\\kappa}(x)=\\frac{(\\kappa+1)p}{g^2}\\frac{m\\cos(2\\theta)-\\omega}{p\\cos^{\\kappa+1}(2\\theta)-1}$, with $\\alpha=\\sqrt{(m-\\omega)/(m+\\omega)}$ and $\\beta=\\sqrt{m^2-\\omega^2}$. Because both $Q$ and $H_2$ scale as $g^{-2/\\kappa}$, the ratio $E/Q$ depends only on $\\omega,\\kappa,p$ and not on $g$. The bound state exists throughout the $(\\kappa,p)$ plane, the transition from double-humped to single-humped density occurs at $\\omega/m=(p\\kappa+1)/(p(\\kappa+1))$, and the stability boundary is governed by $dQ/d\\omega$.","pith_inferences":["Editorial inference: because $E/Q$ at $\\kappa=1$ is available in closed form, this family provides a clean analytic test bed for whether the binding threshold $E/Q<m$ and the stability criterion coincide in nonlinear Dirac models.","Editorial inference: the sharp double-to-single hump boundary could serve as an observable marker in honeycomb-lattice photonics or spin-orbit-coupled Bose-Einstein condensates, where nonlinear Dirac equations appear as envelope equations.","Editorial inference: the same radial factorization may extend to the PT-symmetric and pseudoscalar variants listed as open problems, since the construction relies only on the ratio of scalar and vector densities.","Editorial inference: a numerical spectrum of the linearized operator around $\\Phi(x)$ would settle whether the Vakhitov--Kolokolov criterion is sufficient, and if it fails the stability boundary would need a full spectral analysis."],"forward_implications":["For every $\\kappa>0$ and $p>1$, exact closed-form solitary waves exist on the whole allowed frequency interval $1/p^{1/(\\kappa+1)}<\\omega/m<1$, so the generalized ABS model has a complete stationary-soliton family.","$E/Q$ is independent of $g$, so fixing the charge determines the coupling as a function of $\\omega,\\kappa,p$, enabling direct comparison of soliton families without tuning the coupling.","The charge density changes from double-humped to single-humped at $\\omega/m=(p\\kappa+1)/(p(\\kappa+1))$, giving a sharp geometric signature of the transition.","The Vakhitov--Kolokolov analysis predicts linear stability for $\\kappa<2$ and a small window beyond, with instability for larger $\\kappa$, a prediction that direct numerical evolution can test.","In the non-relativistic limit the model reduces to a modified nonlinear Schrödinger equation with effective coupling $g^2(p-1)/p$, and its solitons are stable for $\\kappa<2$."],"supporting_citations":[{"why":"Supplies the gamma-matrix representation, the radial reduction method, and the known S-S and V-V solitary wave results that the present construction extends.","marker":"[16]"},{"why":"Defines the ABS model at $\\kappa=1, p=2$ whose two-parameter generalization is the subject of this paper.","marker":"[17]"},{"why":"Provides the Vakhitov--Kolokolov criterion $dQ/d\\omega<0$ used to determine the stability regions.","marker":"[18]"},{"why":"Supplies Derrick's theorem used to assess scale-invariance stability of the non-relativistic modified NLSE solitons.","marker":"[19]"}],"fun_headline_variants":["Exact Dirac solitons for every admissible (κ,p)","Dirac solitons' E/Q ignores coupling constant","All admissible (κ,p) admit exact Dirac solitons","Exact solitons span κ–p plane in Dirac model","Exact Dirac solitons shift from double to single hump"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Vakhitov--Kolokolov condition $dQ/d\\omega<0$ is sufficient, not merely necessary, for linear stability of these nonlinear Dirac solitary waves; the paper applies it without proving the spectral and coercivity conditions of the linearized operator, and for $\\kappa>2$ it explicitly defers numerical verification of the predicted unstable region.","fun_headline_variants_meta":{"raw":{"variants":["Exact Dirac solitons for every admissible (κ,p)","Dirac solitons' E/Q ignores coupling constant","All admissible (κ,p) admit exact Dirac solitons","Exact solitons span κ–p plane in Dirac model","Exact Dirac solitons shift from double to single hump"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001771,"raw_usage":{"total_tokens":7150,"prompt_tokens":1271,"completion_tokens":5879,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":887,"completion_tokens_details":{"reasoning_tokens":5792}},"tokens_in":887,"tokens_out":5879,"duration_ms":39164,"temperature":1.0,"reasoning_tokens":5792,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:13:38.390819+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run direct numerical time evolution of the exact solitons for $\\kappa>2$ at frequencies where $dQ/d\\omega>0$ (the Vakhitov--Kolokolov unstable region) and where $E/Q>1$; linear stability would be settled by computing the spectrum of the linearized Dirac operator around $\\Phi(x)$ and looking for eigenvalues with positive real part. If the solitons remain stable where $dQ/d\\omega>0$, the Vakhitov--Kolokolov criterion is not sufficient for this model.","supporting_citations":[{"cited_title":"Mikhailov, JETP Lett 23, 320 (1976), E.A","cited_arxiv_id":null,"evidence_quote":"Provides the Vakhitov--Kolokolov criterion $dQ/d\\omega<0$ used to determine the stability regions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies Derrick's theorem used to assess scale-invariance stability of the non-relativistic modified NLSE solitons."}],"review_version":1}