REVIEW 4 major objections 4 minor 25 references
Solitary waves in the complementary generalized ABS model
T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read 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$.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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$.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [Sec. IV, Eq. (51), Fig. 7] 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.
- [Sec. III E, Sec. VI] 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.
- [Sec. III B, Eq. (33)] 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.
- [Sec. V A, Eq. (70)] 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.
minor comments (4)
- [Eq. (41)] 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.
- [Eq. (5)] 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'.
- [Sec. III B, Eq. (32)] 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.
- [Sec. VI] 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.
Circularity Check
No significant circularity: all central claims follow from the stated Lagrangian by explicit derivation; the numerical extrapolations and the VK stability inconsistency are correctness risks, not circular reductions.
full rationale
The paper derives the solitary-wave solutions of the complementary gABS Lagrangian (Eq. 4) by direct manipulation of the field equation (Eq. 5) and the stationary ansatz (Eq. 7). The ODE for the chiral angle (Eq. 19) and the density profile (Eq. 27) are obtained algebraically from the model, with no fitted parameter entering the solution. The claimed g-independence of E/Q follows from the explicit expressions (28)-(42): both Q and H2 carry the same factor (g^2)^(-1/kappa), so the ratio cancels; this is a derived identity, not a definition. The single-hump property and the full range 0<omega<m follow from positivity of R^2 in Eq. (27). The bound-state region kappa<=kappa_c(q) is obtained by numerically evaluating the derived integrals at selected q, and the universal bound kappa_c<=5/2 is an extrapolation from those samples; this is an evidentiary limitation, not a circular reduction. The paper's stability claim 'stable for kappa<2' is not supported by its own Fig. 7 for q=3/2 and 1.5<kappa<2, where Q(omega) is non-monotonic, but that is an internal correctness problem, not circularity. Citations to the authors' prior work [2,14] supply notation, the general solution method, and standard nonrelativistic reduction steps; the new results are not imported from those references but are re-derived here. No step in the derivation chain is equivalent by construction to its own input.
Assumptions & free parameters
assumptions (6)
- domain assumption Stationary rest-frame ansatz psi(x,t)=R(x)(cos theta(x), sin theta(x)) e^{-i omega t} with 0<omega<m.
- domain assumption Vanishing boundary conditions imply T11=0 (integration constant of motion is zero).
- domain assumption Bound state criterion E/Q < 1 (with m=1).
- domain assumption Vakhitov-Kolokolov criterion dQ/domega < 0 for linear stability.
- domain assumption Derrick's theorem for mNLSE stability analysis.
- standard math Gamma matrix representation gamma0=sigma3, gamma1=i sigma1.
Cite this review
Pith. "Pith review of Solitary waves in the complementary generalized ABS model." pith.science (2026). https://pith.science/paper/EYZSIP3T
@misc{pith2026250612631,
author = {Pith},
title = {Pith review of: Solitary waves in the complementary generalized ABS model},
year = {2026},
howpublished = {\url{https://pith.science/paper/EYZSIP3T}},
note = {Machine review of arXiv:2506.12631}
}
abstract
We obtain exact solutions of the nonlinear Dirac equation in 1+1 dimension of the form $ \Psi(x,t) =\Phi(x) \rme^{-\rmi \omega t}$ where the nonlinear interactions are a combination of vector-vector and scalar-scalar interactions with the interaction Lagrangian given by $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}$, where $\kappa>0$ and $q>1$. This is the complement of the generalization of the ABS model \cite{abs} that we recently studied \cite{ak} and denoted as the gABS model. We show that like the gABS model, in the complementary gABS models the solitary wave solutions also exist in the entire $(\kappa, q)$ plane and further in both models energy of the solitary wave divided by its charge is {\it independent} of the coupling constant $g$. However, unlike the gABS model here all the solitary waves are single humped, any value of $0 < \omega < m$ is allowed and further unlike the gABS model, for this complementary gABS model the solitary wave bound states exist only in case $\kappa \le \kappa_c$, where $\kappa_c$ depends on the value of $q$. Here $\omega$ and $m$ denote frequency and mass, respectively. We discuss the regions of stability of these solutions as a function of $\omega,q,\kappa$ using the Vakhitov-Kolokolov criterion. Finally we discuss the non-relativistic reduction of the two-parameter family of this complementary generalized ABS model to a modified nonlinear Schr\"odinger equation (NLSE) and discuss the stability of the solitary waves in the domain of validity of the modified NLSE.
Figures
Figures from the paper (7 more)
Reference graph
Works this paper leans on
-
[1]
The most important and perhaps the obvious question is about the stability of these solutions. One would like to know the parameter range in the ( κ-q) plane for which the solitary wave bound state solutions are stable
-
[2]
ABS [1] have also considered the PT-invariant generalization of their model. 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
-
[3]
There has been a lot of discussion in the literature about the behavior of the SS as well as VV solitary waves in various external fields [17, 18]. It is then natural to examine the behavior of solitary waves of the complementary gABS model under various external fields
-
[4]
Another interesting model is VV plus (1/q) SS interaction ( q > 0)
In this paper we have considered two-parameter families of models with VV minus (1/q) SS interaction and have shown that the solitary wave bound states exist only if κ <2.5 irrespective of the value of q >1, while in case q = 1.1 the condition gets stricter and solitary wave bound states exist only if κ <2. Another interesting model is VV plus (1/q) SS in...
-
[5]
Another possible interaction is the PS-PS (pseudoscalar-pseudoscalar) interaction
In this paper we have considered two-parameter families of models with a novel ad- mixture of the VV and SS interactions. Another possible interaction is the PS-PS (pseudoscalar-pseudoscalar) interaction. Can one construct a NLD model with the admixture of the PS-PS with the VV and/or with the SS models? Going further, can one construct a NLD model with a...
-
[6]
We hope to address some of these questions in the near future
Finally, perhaps the most interesting question is if one can find relevance of the com- plementary gABS model (at least for one of the allowed q values) in the context of some physical phenomena, for example in the Bose-Einstein condensate or photonics or condensed matter or some other physical system. We hope to address some of these questions in the nea...
-
[7]
A. B. Aceves and S. Wabnitz, Phys. Lett. A 141, 37 (1989); C. M. de Sterke and J. E. Sipe, in Progress in Optics, edited by E. Wolf (Elsevier, Amsterdam, 1994), Vol. XXXIII; C. M. de Sterke, D. G. Salinas and J. E. Sipe, Phys. Rev. E 54, 1969 (1996)
work page 1989
-
[8]
N.V Alexeeva, I. V. Barashenkov and A. Saxena, Annals of Physics 403, 198 (2019)
work page 2019
Show all 25 references
-
[9]
Khare, F
A. Khare, F. Cooper, J. Dawson and A. Saxena. ArXiv:2504.13299v1 [nlin.PS] 17 Apr 2025. 17
2025 arXiv
-
[10]
Ivanenko, ZhETF (JETP) 8, 260 (1938)
D. Ivanenko, ZhETF (JETP) 8, 260 (1938)
1938
-
[11]
Thirring, Ann
W.E. Thirring, Ann. Phys. 3, 91 (1958)
1958
-
[12]
Soler, Phys
M. Soler, Phys. Rev. D 1, 2766 (1970)
1970
-
[13]
Gross and A
D.J. Gross and A. Neveu, Phys. Rev. D 10, 3235 (1974)
1974
-
[14]
Cooper, A
F. Cooper, A. Khare, B. Mihaila and A. Saxena, Phy. Rev. E 82, 036604 (2010). We follow this paer’s notations as well as several of its results
2010
-
[15]
Ablowitz and Y
M.J. Ablowitz and Y. Zhu, Phys. Rev. A 82, 013840 (2010)
2010
-
[16]
It is well known that for the NLSE, this method is a reliable tool in determining the regions of instability
to study the stability with respect to the scale transformations keeping the mass M = Z ∞ −∞ dx |ψ|2 , (66) of the soliton fixed. It is well known that for the NLSE, this method is a reliable tool in determining the regions of instability. Derrick’s theorem states that if we m...
-
[17]
O Peleg, G Bartal, B Freedman, O Manela, M Segev and D N Christodoulides, Phys. Rev. Lett. 98, 103901 (2007)
2007
-
[18]
J Dalibard, F Gerbier, G Juzeli¯ unas, and P ¨Ohberg, Rev. Mod. Phys. 83, 1523 (2011)
2011
-
[19]
D. K. Campbell and A. R. Bishop, Phys. Rev. B 24, 4859 (1981); Nucl. Phys. B 200, 297 (1982)
1981
-
[20]
Aubry, Physica D 103, 201 (1997)
S. Aubry, Physica D 103, 201 (1997)
1997
-
[21]
Lee, T.K
S.Y. Lee, T.K. Kuo and A. Gavrielides, Phys. Rev. D 12, 2249 (1975)
1975
-
[22]
Vakhitov and A.A
N.G. Vakhitov and A.A. Kolokolov, Radiophys. Quantum Electron 16, 783 (1973)
1973
-
[23]
Derrick, J
G.H. Derrick, J. Math. Phys. 5, 1252 (1964)
1964
-
[24]
Mertens, N.R
F.G. Mertens, N.R. Quintero and A.R. Bishop, Phys. Rev. E81, 016608 (2010); F. Mertens, N.R. Quintero, I.V. Barashenkov and A.R. Bishop, Phys. Rev. E84, 026614 (2011); F.G. Mertens, F. Cooper, N.R. Quintero, A. Khare and A. Saxena, J. Phys. A 49, 065402 (2016)
2010
-
[25]
Cooper, A
F. Cooper, A. Khare, N.R. Quintero, F.G. Mertens and A. Saxena, Phys. Rev. E 85, 046607 (2012); F. Cooper, N.R. Quintero, B. Sanchez-Ray, F.G. Mertens and A. Saxena, J. Phys. A 53 455702 (2020) 18
2012
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.