{"id":"29882312-3043-4564-8b93-7433fba6b63e","arxiv_id":"2506.12631","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Exact solitary waves are found for the complementary generalized ABS nonlinear Dirac model, with full frequency range and a q-dependent bound-state threshold, stable for kappa < 2.","lead":"The authors derive exact wave shapes for a new family of nonlinear Dirac equations that mix two types of particle interactions. These solutions are stable over a wide parameter range, and their energy per unit charge does not depend on the coupling strength.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's VK stability claim 'stable for κ<2' is contradicted by its own Fig. 7: for q=3/2 and 1.5<κ<2, Q(ω) has both a minimum and a maximum, so dQ/dω>0 on an interval, violating Eq. (51).","rationale":"The exact construction of the solitary waves and the E/Q independence appear to be derived correctly and are not the weak point. The load-bearing weakness is in Section IV: the stability conclusion is not merely heuristic, it is contradicted by the paper's own description of Fig. 7. Because the VK criterion is a sign condition on dQ/dω, a Q(ω) curve with a minimum and a maximum necessarily has an interval of positive slope, hence VK instability. The text does not restrict 'stable for κ<2' to large q or to a specific ω window, so a reader cannot accept the stated stability region as written. The κ_c=5/2 threshold is also under-supported, being extrapolated from only two q values; this reinforces the conditional verdict, though the internal contradiction is more decisive. The proposed numerical check settles whether the stability statement can be repaired with qualifications or must be withdrawn, and the conditional verdict should remain until that is resolved.","tokens_in":11516,"tokens_out":7026,"duration_ms":83213,"concrete_test":"Recompute Q(ω) from Eqs. (28)-(30) for q=3/2 and κ=1.7, 1.8, 1.9 on a fine ω-grid, locate the reported ω_min and ω_max, and evaluate dQ/dω on that interval. If Q is non-monotone while E/Q<1, the statement that solutions are stable for κ<2 is false; if the positive-slope interval contains no bound states, the paper must state the precise ω-region for which the stability claim applies.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV uses the Vakhitov-Kolokolov criterion dQ/dω<0 and states that the solitary-wave bound states are stable for κ≤2. But the same section reports that for q=3/2, 'when 3/2<κ<2, Q has a minimum at ω_min ... Q also has a maximum ... denoted ω_max' (Fig. 7). On (ω_min, ω_max), dQ/dω>0, so condition (51) fails and the unqualified stability claim is internally inconsistent. This is not a disagreement with external consensus; it is the paper's own numerics contradicting its conclusion. The q=3/2 case is not marginal: it is singled out as far from the VV limit, so the claimed universal 'stable for κ<2' cannot hold without an explicit ω- and q-restriction. Separately, the bound-state cutoff κ_c=5/2 is inferred from q=10 and q=1.1 only, and the paper itself says κ_c depends on q, so the universal maximum is not demonstrated; but the VK contradiction is the more immediate obstacle to accepting the stability claims.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces the complementary gABS model, a 1+1-dimensional nonlinear Dirac theory with vector-vector minus scalar-scalar interactions parameterized by kappa>0 and q>1. It constructs exact stationary solitary-wave solutions by the standard reduction used in the authors' earlier work, computes the charge Q and energy E, shows that E/Q is independent of the coupling g, and investigates bound states (E/Q<1) and stability via the Vakhitov-Kolokolov criterion. It ends with a non-relativistic reduction to a modified nonlinear Schrödinger equation and a Derrick-type stability analysis. The main advertised results are: solutions exist for all kappa>0, q>1 and 0<omega<m; all such waves are single-humped; solitary-wave bound states exist only for kappa<=kappa_c(q), with a claimed universal maximum kappa_c=5/2; and the waves are stable for kappa<2.","tokens_in":11740,"tokens_out":4858,"duration_ms":59034,"significance":"If correct, the paper would provide an exact two-parameter family of solitary waves for a new nonlinear Dirac model, with the E/Q independence of g following as an exact scaling consequence rather than a numerical accident. Strengths of the manuscript include the explicit derivation leading to Eq. (27), the exact closed-form charge and energy at kappa=1 in Eqs. (46) and (48), and the clear identification of the q->infinity limit as the pure vector-vector model. However, the paper's most quantitatively new claims, namely the universal bound-state cutoff kappa_c=5/2 and the unqualified stability statement for kappa<2, are either extrapolated from a few numerical cases or contradicted by the paper's own numerics. These issues are load-bearing for the abstract and conclusions, so the paper needs substantial revision before the central claims can be accepted.","major_comments":[{"comment":"The conclusion that the solitary-wave bound states are stable for kappa<2 is internally inconsistent with the paper's own numerical results. For q=3/2 and 3/2<kappa<2, the text states that Q(omega) has a minimum at omega_min and a maximum at omega_max, so dQ/domega>0 on the interval (omega_min,omega_max). This directly violates the Vakhitov-Kolokolov condition dQ/domega<0 stated in Eq. (51). The stability claim must therefore be restricted either to specific omega-intervals on which dQ/domega<0 or to values of q for which monotonicity is proven; as written, the claim is not supported by the criterion the paper itself adopts.","section":"Sec. IV, Eq. (51), Fig. 7"},{"comment":"The universal bound-state cutoff kappa_c=5/2 is not demonstrated. The evidence consists of numerical results for q=10 and q=1.1 plus a remark that 'other simulations' lead to the same conclusion, while the text itself says kappa_c depends on q. In fact, for q=1.1 the paper reports that no bound states exist once kappa>=2, which is a different cutoff from 5/2. No analytical argument or systematic scan over q is given that would justify the claim that kappa_c(q)<=5/2 for all q>1. Since this is one of the abstract's headline results, it needs either a proof or a much more modest and precisely documented numerical statement.","section":"Sec. III E, Sec. VI"},{"comment":"The single-hump claim is asserted rather than proved. After writing the derivative dR^2/dx in Eq. (33), the paper simply states that x=0 is always a maximum and that R^2 always has a single hump. No sign analysis, uniqueness argument, or numerical tracking of the extrema is provided, even though the expression contains competing positive and negative contributions and the denominator can vary nontrivially with q. Because the contrast between single-humped and double-humped behavior is a central difference from the gABS model, this step needs a rigorous argument or at least a systematic numerical verification.","section":"Sec. III B, Eq. (33)"},{"comment":"The Derrick-type stability result in the non-relativistic reduction is stated more sharply than the displayed calculation supports. Equation (70) gives d^2H/dbeta^2 = 2(2-kappa)H1 - 2(2+kappa)H2, so positivity requires that H2 be sufficiently small relative to H1, not merely that kappa<2. The text asserts that H2 is positive and small, but no quantitative bound is provided, and the subsequent conclusion 'stable for kappa<2' is not rigorously implied by the equation. This should be rephrased as a leading-order statement or supplemented with bounds.","section":"Sec. V A, Eq. (70)"}],"minor_comments":[{"comment":"The argument of J is written as J(omega,kappa,p) but the model parameter is q; this is a typo that should be corrected.","section":"Eq. (41)"},{"comment":"Equation (5) is missing the equals sign and zero on the right-hand side; the displayed field equation ends with the term -g^2/q (bar-psi psi)^kappa and needs '=0'.","section":"Eq. (5)"},{"comment":"The symbol C is reused with a different meaning in Eq. (32): it previously denoted C(omega,kappa) in Eq. (29), but in Eq. (32) it absorbs additional x-independent constants. Using a different symbol, such as C_0 or C', would avoid confusion.","section":"Sec. III B, Eq. (32)"},{"comment":"The conclusions state that bound states exist only if kappa<2.5, whereas the body of the paper claims kappa<=kappa_c(q) with kappa_c depending on q. The concluding sentence should be worded to reflect the q-dependent cutoff and the fact that the universal statement is only an upper bound.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a direct continuation of the authors' own recent work on the gABS model, and the amount of self-citation is understandable. The novelty of the complementary model is clear. The main obstacle to acceptance is not novelty but internal consistency: the stability section contradicts its own numerical figure, and the universal cutoff claim is under-supported. If the authors restrict the stability statements to the parameter regions where their own numerics satisfy dQ/domega<0, and either prove or substantially weaken the kappa_c=5/2 universality claim, the paper could become publishable. I would not support rejection at this stage, since the exact solution construction appears sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful part of this paper is the exact solitary-wave construction for the complementary gABS model (positive VV, negative SS). The reduction is standard but correct, and the E/Q independence follows exactly from scaling. Unlike their earlier gABS paper, the frequency range is full 0<ω<m and the profiles appear to be single-humped; the new two-parameter family is worth having.\n\nThe soft spots are real. First, the VK stability claim \"stable for κ<2\" is contradicted by the paper's own Fig. 7: for q=3/2 and 3/2<κ<2, Q(ω) has a minimum and a maximum, so dQ/dω>0 on an interval, violating condition (51). The text acknowledges the min and max but still draws the unqualified conclusion. That needs an explicit ω-, κ-, and q-restriction, or a corrected analysis. Second, the universal bound-state cutoff κ_c=5/2 is inferred from q=10 and q=1.1 (plus unshown \"simulations\"); the paper itself says κ_c depends on q, so the claim that κ>5/2 never supports bound states is not demonstrated. Third, the single-hump proof is asserted after Eq. (33) rather than shown by analyzing the derivative; I suspect it is true, but it is not proven.\n\nThe derivation of the solutions themselves is sound, and the existence result is not in doubt. The stability section needs repair; the bound-state cutoff needs more evidence or a softened claim. This is a genuinely useful paper for the nonlinear Dirac solitary-wave community if those issues are fixed.\n\nI would send it to a serious referee, with the warning that the stability claims need scrutiny. I would cite it for the exact solutions and the E/Q identity.","headline":"Exact solutions for the complementary gABS model are solid, but the paper's own Fig. 7 undercuts its VK stability claim and the κ_c=5/2 cutoff is extrapolated from only two q values.","tokens_in":12280,"tokens_out":2485,"would_cite":true,"duration_ms":28165,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The complementary generalized ABS nonlinear Dirac model has exact single-humped solitary waves for every $\\kappa>0$ and every $0<\\omega<m$, with bound states only below a $q$-dependent threshold $\\kappa_c(q)\\le 5/2$.","keywords":["complementary gABS model","nonlinear Dirac equation","solitary waves","bound states","Vakhitov-Kolokolov criterion","modified nonlinear Schrödinger equation","charge-energy ratio","1+1 dimensions"],"falsifier":"Evaluate the integrals $I(\\omega,\\kappa,q)$ and $J(\\omega,\\kappa,q)$ at, say, $q=3$ and $\\kappa=2.6$ across the full range $0<\\omega<1$; if $E/Q$ dips below 1 at any $\\omega$, the claimed $\\kappa_c\\le5/2$ bound fails. For the stability claim, find any $q>1$ and $\\kappa<2$ for which $dQ/d\\omega\\ge0$ on an interior $\\omega$-interval, which would falsify the Vakhitov-Kolokolov stability conclusion for that region.","tokens_in":11319,"feed_emoji":"🌊","tokens_out":13181,"duration_ms":131665,"temperature":0.7,"pith_summary":"The paper introduces a complementary version of the generalized ABS (gABS) model: a nonlinear Dirac equation in 1+1 dimensions with vector-vector and scalar-scalar couplings of opposite signs, controlled by two parameters, $\\kappa>0$ and $q>1$. It establishes that this model admits exact solitary-wave solutions of the form $\\Phi(x)e^{-i\\omega t}$ for every $\\kappa>0$, every $q>1$, and every frequency $0<\\omega<m$, and that, unlike the gABS model, all such waves are single-humped. It further shows that although the energy $E$ and charge $Q$ each scale with the coupling constant $g$, their ratio $E/Q$ is independent of $g$, depending only on $\\omega$, $\\kappa$, and $q$. The most distinctive result is that solitary-wave bound states, defined by $E/Q<1$, exist only for $\\kappa\\le\\kappa_c(q)$, with numerical evidence placing the maximum possible threshold at $\\kappa_c=5/2$ and a stricter threshold $\\kappa_c=2$ at $q=1.1$. Stability is then examined with the Vakhitov-Kolokolov criterion and, in the non-relativistic limit, with a modified nonlinear Schr\\\"odinger equation, leading the authors to conclude stability for $\\kappa<2$.","feed_headline":"Exact solitons found for every κ in complementary Dirac model","feed_subtitle":"In the new two-parameter family, soliton energy-to-charge ratio is independent of the coupling constant g.","key_machinery":"The load-bearing object is the phase-angle reduction: writing the spinor as $R(x)(\\cos\\theta(x),\\sin\\theta(x))e^{-i\\omega t}$ turns the two-component nonlinear Dirac equation into a first-order equation for $\\theta$, $d\\theta/dx=\\kappa[m\\cos(2\\theta)-\\omega]$, which is solved exactly by $\\tan\\theta=\\alpha\\tanh(\\kappa\\beta x)$ with $\\beta=\\sqrt{m^2-\\omega^2}$. That solution, together with the algebraic relation $R^2=[(\\kappa+1)(m\\cos 2\\theta-\\omega)]/[g^2(1-q^{-1}\\cos^{\\kappa+1}(2\\theta))]^{1/\\kappa}$, converts the problem into two integrals, $I(\\omega,\\kappa,q)$ for $Q$ and $J(\\omega,\\kappa,q)$ for $H_2$. The key identity connecting them, $H_1=\\kappa H_2-\\kappa\\omega Q$, exposes the $g$-independence of $E/Q$, since each of $Q$ and $H_2$ scales as $g^{-2/\\kappa}$. Bound-state existence is read off from $E/Q<1$; stability is decided by the sign of $dQ/d\\omega$ according to the Vakhitov-Kolokolov criterion.","core_discovery":"The central claim is that the complementary gABS interaction, $L_I = \\frac{g^2}{\\kappa+1}[(\\bar{\\psi}\\gamma_\\mu\\psi)(\\bar{\\psi}\\gamma^\\mu\\psi)]^{(\\kappa+1)/2} - \\frac{g^2}{q(\\kappa+1)}(\\bar{\\psi}\\psi)^{\\kappa+1}$, is exactly solvable for solitary waves in the whole parameter range. Using the rest-frame ansatz $\\Psi(x,t)=R(x)(\\cos\\theta(x),\\sin\\theta(x))e^{-i\\omega t}$, the Dirac equation reduces to $d\\theta/dx=\\kappa[m\\cos(2\\theta)-\\omega]$, whose integral $\\tan\\theta=\\alpha\\tanh(\\kappa\\beta x)$ with $\\beta=\\sqrt{m^2-\\omega^2}$ yields a closed-form $R^2$ for every $0<\\omega<m$. From this profile the paper computes the charge and energy as convergent integrals, derives the relation $H_1=\\kappa H_2-\\kappa\\omega Q$, and obtains $E/Q=\\kappa H_2/Q+(1-\\kappa)\\omega$, which is manifestly independent of $g$ because both $Q$ and $H_2$ are proportional to $g^{-2/\\kappa}$. The paper then argues that bound states exist only for $\\kappa\\le\\kappa_c(q)$, with $\\kappa_c(q)$ numerically bounded by $5/2$, and uses the Vakhitov-Kolokolov condition together with a non-relativistic modified NLSE reduction to conclude stability for $\\kappa<2$.","pith_inferences":["Because $E/Q$ is $g$-independent, a system realizing this model could use the measured soliton energy-charge ratio as a direct probe of the bare parameters $\\kappa$ and $q$, bypassing calibration of the coupling constant; the paper does not draw this experimental corollary.","The paper's $\\kappa_c(q)$ threshold comes from numerical evaluation of the integrals at selected $q$ values; a systematic scan of $\\kappa_c(q)$ over $q\\in(1,\\infty)$, checking whether it interpolates from $2$ at $q\\to1$ toward $5/2$ at large $q$, would settle the claimed universality.","The Vakhitov-Kolokolov stability conclusion for $\\kappa<2$ is strained by the paper's own $q=3/2$ numerics, where $Q(\\omega)$ is non-monotonic for $1.5<\\kappa<2$; a direct spectral-stability computation of linearized modes would resolve whether instability windows exist there.","The integrable massive Thirring point sits at $q=\\infty$, and an open question the authors flag is whether integrability survives for large finite $q$; a perturbative check of conserved charges at order $1/q$ would be a natural next step."],"forward_implications":["In the whole $(\\kappa,q)$ plane with $\\kappa>0$ and $q>1$, exact single-humped solitary waves exist for every rest-frame frequency $0<\\omega<m$.","For any fixed $\\kappa$ and $q$, the ratio $E/Q$ is a function of $\\omega$, $\\kappa$, and $q$ alone, so measurements of the energy-to-charge ratio would not require knowledge of the coupling strength $g$.","No solitary-wave bound states exist for $\\kappa>\\kappa_c(q)$, and numerically $\\kappa_c(q)$ never exceeds $5/2$; at $q=1.1$ the threshold is already $\\kappa_c=2$.","Both the Vakhitov-Kolokolov analysis and the non-relativistic modified NLSE reduction point to stability for $\\kappa<2$, with a stable low-frequency region surviving for $\\kappa>2$ when $q$ is large.","In the non-relativistic limit, the complementary model reduces to the modified NLSE with the coupling rescaled by $(q-1)/q$, so the $\\mathrm{sech}^{1/\\kappa}(\\kappa\\beta x)$ soliton profile persists with an adjusted amplitude."],"supporting_citations":[{"why":"Introduces the ABS model with the vector-vector and scalar-scalar admixture that this paper generalizes to the complementary sign choice.","marker":"[1]"},{"why":"The authors' own gABS companion study; provides the model compared against and the ansatz and notation extended here.","marker":"[2]"},{"why":"Supplies the general-$\\kappa$ solitary-wave construction, gamma-matrix conventions, and the modified NLSE reduction used throughout.","marker":"[14]"},{"why":"The Vakhitov-Kolokolov criterion $dQ/d\\omega<0$ is used to infer linear stability of the solitary waves.","marker":"[15]"},{"why":"Derrick's scaling argument is used to assess stability in the non-relativistic modified NLSE regime.","marker":"[16]"}],"fun_headline_variants":["Solitary waves solved for all κ in complementary gABS model","Complementary gABS solitons: E/Q independent of coupling g","All (κ,q) admit solitary waves in complementary Dirac model","Exact solitons in complementary gABS: bound states for κ≤κc","Nonlinear Dirac solitons with g-free energy-to-charge ratio"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's universal conclusions rest on numerical evaluation of two integrals at only a few $q$ values, extrapolated to all $q>1$, even though at $q=3/2$ the charge $Q(\\omega)$ is non-monotonic for $1.5<\\kappa<2$, so neither the $\\kappa_c\\le5/2$ bound nor the 'stable for $\\kappa<2$' statement is proven uniformly.","fun_headline_variants_meta":{"raw":{"variants":["Solitary waves solved for all κ in complementary gABS model","Complementary gABS solitons: E/Q independent of coupling g","All (κ,q) admit solitary waves in complementary Dirac model","Exact solitons in complementary gABS: bound states for κ≤κc","Nonlinear Dirac solitons with g-free energy-to-charge ratio"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000756,"raw_usage":{"total_tokens":3514,"prompt_tokens":1249,"completion_tokens":2265,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":865,"completion_tokens_details":{"reasoning_tokens":2170}},"tokens_in":865,"tokens_out":2265,"duration_ms":18053,"temperature":1.0,"reasoning_tokens":2170,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:45:45.972894+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the integrals $I(\\omega,\\kappa,q)$ and $J(\\omega,\\kappa,q)$ at, say, $q=3$ and $\\kappa=2.6$ across the full range $0<\\omega<1$; if $E/Q$ dips below 1 at any $\\omega$, the claimed $\\kappa_c\\le5/2$ bound fails. For the stability claim, find any $q>1$ and $\\kappa<2$ for which $dQ/d\\omega\\ge0$ on an interior $\\omega$-interval, which would falsify the Vakhitov-Kolokolov stability conclusion for that region.","supporting_citations":[{"cited_title":"One would like to know the parameter range in the ( κ-q) plane for which the solitary wave bound state solutions are stable","cited_arxiv_id":null,"evidence_quote":"Introduces the ABS model with the vector-vector and scalar-scalar admixture that this paper generalizes to the complementary sign choice."},{"cited_title":"The interesting question is whether one can also construct the PT-invariant variant of the complementary gABS model, and if yes, can one also obtain its exact solutions","cited_arxiv_id":null,"evidence_quote":"The authors' own gABS companion study; provides the model compared against and the ansatz and notation extended here."},{"cited_title":"Cooper, A","cited_arxiv_id":null,"evidence_quote":"Supplies the general-$\\kappa$ solitary-wave construction, gamma-matrix conventions, and the modified NLSE reduction used throughout."},{"cited_title":"Ablowitz and Y","cited_arxiv_id":null,"evidence_quote":"The Vakhitov-Kolokolov criterion $dQ/d\\omega<0$ is used to infer linear stability of the solitary waves."},{"cited_title":"It is well known that for the NLSE, this method is a reliable tool in determining the regions of instability","cited_arxiv_id":null,"evidence_quote":"Derrick's scaling argument is used to assess stability in the non-relativistic modified NLSE regime."}],"review_version":1}