{"id":"a5fb2245-ca83-4c9a-9bd4-f4bce3bfddec","arxiv_id":"2412.00838","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The algebraic double-soliton solutions of the massive Thirring model arise as a singular limit of Riemann-Hilbert double-pole solutions and correspond to a double embedded eigenvalue at zeta equals i.","lead":"This paper derives exponential and algebraic double-soliton solutions of the massive Thirring model from inverse scattering theory, and shows the algebraic ones correspond to double eigenvalues embedded in the continuous Lax spectrum. It matters because it gives the first spectral explanation of previously explicit rational solutions, resolving a conjecture about embedded eigenvalues in this integrable system.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exact 'double embedded eigenvalue' claim rests on an unquoted external lemma [15, Lemma 6.4]; Remark 2.4 verifies only two decay rates, not the lemma's full hypotheses, so algebraic multiplicity exactly 2 is not established in the paper.","rationale":"The paper's genuinely new contribution is the spectral interpretation of the algebraic double-solitons: they are claimed to realize double embedded eigenvalues, resolving a conjecture in [15]. The explicit formulas for ψ0 and ψ1 in Theorem 2.3 are a strong piece of evidence: they directly show the existence of an eigenvector and a generalized eigenvector, hence algebraic multiplicity at least 2. The paper also correctly notes (Remark 2.5) that the algebraic solution matches the independently obtained solution of [8], which supports the claim that (2.20)-(2.21) solves the MTM. The soft spot is exactly where the reader placed it: the passage from 'a Jordan chain of length two exists' to 'the eigenvalue is double, with exactly one eigenvector and one generalized eigenvector' requires an upper bound on the algebraic multiplicity. That upper bound is delegated to [15, Lemma 6.4], but the lemma is not stated and its hypotheses are not checked beyond the decay rates in Remark 2.4. This is a missing-support issue of the kind that should be flagged: the reader cannot verify the central spectral claim without consulting the external reference. The proposed concrete test — a direct asymptotic/ODE computation of the algebraic multiplicity for the explicit rational potential — would settle the matter without relying on the external lemma. Because the concern is narrow and fixable (either reproduce and verify the lemma, or supply a direct proof), the appropriate verdict is CONDITIONAL rather than REJECT or UNVERDICTED. The reader's ACCEPT is close, but the central claim is stated more strongly than the paper's internal argument supports.","tokens_in":28226,"tokens_out":8490,"duration_ms":77803,"concrete_test":"For the rational potential (2.20)-(2.21), compute the general solution of the Lax ODE ∂xψ = L(ualg,valg,i)ψ at ζ=i (e.g., by eliminating one component and using Frobenius analysis at x=±∞), and count the solutions lying in H^1(R). If the H^1 eigenspace is exactly one-dimensional, additionally solve the inhomogeneous equation (L(i)-i)ψ = ∂ζL|ζ=i ψ0 in H^1 and check that every solution is a scalar multiple of ψ1 modulo the kernel. If both checks pass, the algebraic multiplicity is exactly 2 and Remark 2.4 is vindicated; if not, the multiplicity statement in Theorem 2.3 must be weakened to 'at least double'.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Theorem 2.3 and the abstract claim that the algebraic double-soliton (2.20)-(2.21) has a double embedded eigenvalue ζ=i with 'only one eigenvector ψ0∈H^1 and one generalized eigenvector ψ1∈H^1'. What the paper actually proves is that ψ0 and ψ1 satisfy (2.18)-(2.19), i.e., that there is a Jordan chain of length two. This establishes algebraic multiplicity at least 2 and geometric multiplicity at least 1, but it does not by itself rule out additional independent H^1 eigenfunctions at ζ=i, nor does it rule out a longer Jordan chain. The only support for exactness is Remark 2.4, which says the pair satisfies 'the criterion for the spatial decay in [15, Lemma 6.4]' and then names only the rates ψ0=O(|x|^-2), ψ1=O(|x|^-1). The full hypotheses of that lemma are not reproduced, and no verification beyond the decay rates is supplied. If the lemma requires additional conditions — e.g., a non-resonance or normalization condition, a weighted decay condition on the potential, or a specific asymptotic normalization of the Jost solutions — the 'only one eigenvector' and 'double' parts of Theorem 2.3 are unsupported. The central conjecture (multiple embedded eigenvalues may exist) would still be resolved by the exhibited Jordan chain, so this is a correctness-of-characterization issue rather than a collapse of the paper's main application; nevertheless the precise spectral statement is load-bearing because the paper presents it as the resolution of the conjecture.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper constructs exponential and algebraic double-soliton solutions of the massive Thirring model in laboratory coordinates using the inverse scattering transform. For a double pole of the Riemann–Hilbert problem at λ0=e^{iγ}, the authors derive explicit exponential double-soliton formulas and show that the associated Lax pair has a Jordan chain of length two at ζ0=e^{iγ/2}. Taking a singular limit γ→π after a rescaling of the translation parameters, they obtain rational algebraic double-soliton solutions and exhibit an eigenvector and generalized eigenvector at the embedded value ζ=i. The paper interprets this as resolving the conjecture that multiple embedded eigenvalues can occur in the MTM spectral problem, and it also verifies consistency with earlier bilinear Hirota results.","tokens_in":28517,"tokens_out":6535,"duration_ms":56605,"significance":"If the spectral characterization is fully established, this is a substantial contribution: it provides the first explicit construction, via the Riemann–Hilbert method, of algebraic double-solitons of the massive Thirring model and connects them to a double embedded eigenvalue on the continuous spectrum. The derivation is systematic, with explicit residue-coefficient computations, closed-form potentials, and direct verification of the Jordan-chain equations (2.18)–(2.19). The consistency check with the bilinear results of [8] and the careful treatment of the singular limit are additional strengths. The paper is written with enough detail that the main algebraic steps can be reproduced, and the explicit formulas are a useful resource for further study of embedded eigenvalues in integrable systems.","major_comments":[{"comment":"The theorem states that ζ=i is a double embedded eigenvalue with only one eigenvector ψ0∈H^1 and one generalized eigenvector ψ1∈H^1. The proof in Section 4.2 constructs ψ0 and ψ1 satisfying (2.18)–(2.19), i.e., a Jordan chain of length two. This establishes algebraic multiplicity at least two and geometric multiplicity at least one, but it does not by itself exclude additional linearly independent H^1 eigenfunctions at ζ=i or a longer Jordan chain. The only support for the 'only one' and 'double' parts is Remark 2.4, which cites [15, Lemma 6.4] but verifies only the decay rates ψ0=O(|x|^{-2}) and ψ1=O(|x|^{-1}); the lemma's remaining hypotheses are not reproduced. Please either quote the lemma and verify all of its hypotheses for the explicit solution (2.20)–(2.21), or replace the 'double embedded eigenvalue with only one eigenvector' statement by the weaker (but sufficient for the conjecture) statement that the algebraic solution admits a Jordan chain of length two at ζ=i.","section":"Theorem 2.3 and Remark 2.4"},{"comment":"The passage from exponential to algebraic double-solitons is made by Taylor expanding D(M), N_u, and N_v in ε=π−γ and taking the limit after rescaling x̃0. This is a formal asymptotic argument; no uniform-convergence estimate or direct substitution into (1.1) is provided. Since (2.20)–(2.21) are claimed as exact solutions, either a direct verification of (1.1) or an explicit statement that they coincide with the bilinear solutions of [8] (under the transformation in Remark 2.5) should be part of the proof of Theorem 2.3. Remark 2.5 currently makes the connection only as a note after the theorem.","section":"Section 4"}],"minor_comments":[{"comment":"The displayed identity |u|^2+|v|^2 = (|N_u|^2+|N_v|^2)/|D(M)|^2 = 2i ∂_x log D(M)/D(M) is miswritten; the rightmost expression should presumably be 2i ∂_x log(D(M)/\\overline{D(M)}) or an equivalent bilinear identity. Please correct the formatting.","section":"Eq. (3.31)"},{"comment":"There is a typo 'folows' for 'follows', and the notation for the scattering coefficient is inconsistent: the residue term for λ0 is written with \\tilde α(λ) while the double-pole term uses α(λ), and similarly \\check α appears for λ̄0. Please unify the notation.","section":"Section 3.1, proof of Proposition 3.1"},{"comment":"The statement that the double-soliton solutions 'only have two non-trivial parameters' is confusing because Theorem 2.1 also contains the parameters c, x0, and t0 arising from the Lorentz and translation symmetries. Please clarify that these are gauge parameters and that γ and x̃0 are the shape parameters.","section":"Remark 2.2"},{"comment":"The factors exp(± i/4 ∫_x^∞ (|u|^2+|v|^2)dy) in (2.22)–(2.23) are evaluated at (u,v)=(ualg,valg), which decay only as O(|x|^{-1}); a sentence confirming that the integrals converge and that the resulting ψ0 and ψ1 belong to H^1(R,C^2) would be helpful.","section":"Theorem 2.3, definitions of n0 and n1"},{"comment":"The matching between the bilinear expressions (B.2)–(B.4) and the RH solution in Theorem 2.1 is stated after the substitutions, but the intermediate definitions of ξ and η in (B.2)–(B.4) are not given in the text of Appendix B; please add them for readability.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The paper is well suited to nlin.SI and contains a substantial amount of correct-looking, explicit computation. The main issue is the unsupported 'only one eigenvector' part of Theorem 2.3; the authors should either supply the full hypotheses of [15, Lemma 6.4] and verify them, or weaken the claim to the existence of a Jordan chain of length two, which already resolves the conjecture about multiple embedded eigenvalues. The formal limiting argument in Section 4 should also be backed by a direct verification or an explicit appeal to [8]. I would be happy to accept after these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real content here is the double-pole Riemann–Hilbert construction for the massive Thirring model in laboratory coordinates, plus the singular limit that produces the algebraic double-solitons. That is genuinely new: previous bilinear work gave the formulas, but not the IST derivation or the spectral interpretation. The exponential double-solitons match the bilinear construction in Appendix B, which is a good consistency check. The explicit eigenvector and generalized eigenvector for the algebraic solution, with the Jordan chain at ζ = i, is the strongest part of the paper, and it does resolve the conjecture that embedded eigenvalues of higher algebraic multiplicity can exist in this Lax problem.\n\nThe main soft spot is exactly what the stress-test note says. Theorem 2.3 states there is 'only one eigenvector' and one generalized eigenvector, i.e. algebraic multiplicity exactly two. What the paper actually proves is that the exhibited pair forms a Jordan chain of length two: that gives algebraic multiplicity at least two and geometric multiplicity at least one. To rule out additional independent H^1 eigenfunctions at ζ = i, Remark 2.4 defers to [15, Lemma 6.4] and only names the decay rates ψ0 = O(|x|^-2), ψ1 = O(|x|^-1). The lemma's full hypotheses are not quoted or checked. So the exact characterization is not established in the paper. This is not fatal to the main application, since the Jordan chain itself already demonstrates a multiple embedded eigenvalue, but it is a real gap in the stated theorem.\n\nA smaller caveat: the algebraic limit is taken by formal ε-expansions. No uniform convergence is shown, and the limiting functions are not directly substituted into (1.1). Given the explicit formulas and the spectral verification, I doubt this is wrong, but it is a loose end that a referee should ask about.\n\nOverall the paper is careful and the algebraic work is substantial. The computations are long and not machine-checked, but the structure is coherent and the comparison with known bilinear results is reassuring. This is a paper for integrable-systems people and for anyone working on embedded eigenvalues in Lax operators. It deserves a serious referee, not a desk rejection. I would send it out with a request to verify the hypotheses of [15, Lemma 6.4] and, if possible, a direct check that the limiting solution satisfies the PDE.","headline":"Solid RH-construction paper whose headline spectral claim is slightly ahead of what is proved: the Jordan chain is there, but the 'exactly one eigenvector' part rests on an external lemma the paper does not reproduce.","tokens_in":29073,"tokens_out":1443,"would_cite":true,"duration_ms":15332,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q51","37K15","35Q41"],"pacs":[],"model":"deepseek-v4-flash","headline":"The algebraic double-soliton solutions of the massive Thirring model carry a double embedded eigenvalue at $\\zeta=i$, confirming that multiple embedded eigenvalues occur in the model's Lax spectrum.","keywords":["massive Thirring model","double-soliton solutions","algebraic solitons","embedded eigenvalues","Riemann–Hilbert problem","inverse scattering transform","Lax spectrum","rational solutions"],"falsifier":"Compute the dimension of the $H^1$ solution space of $(\\partial_x-L(u_{\\rm alg},v_{\\rm alg},i))\\psi=0$ and the length of the Jordan chain at $\\zeta=i$. If the kernel contains a second independent $H^1$ eigenvector, or if the generalized eigenvector equation has no $H^1$ solution with $\\psi_0=O(|x|^{-2})$, then the double embedded eigenvalue claim fails; this can be checked by direct substitution of (2.20)–(2.21) into (2.18)–(2.19) or by numerically evaluating the Evans function on the imaginary axis.","tokens_in":27997,"feed_emoji":"🌊","tokens_out":5406,"duration_ms":49128,"temperature":0.7,"pith_summary":"The paper establishes that the massive Thirring model, a relativistic nonlinear Dirac equation in one spatial dimension, admits two families of double-soliton solutions, and it identifies what each family means for the Lax spectrum. Exponential double-solitons correspond to double isolated eigenvalues of the spectral problem, while algebraic double-solitons, which decay only as $O(|x|^{-1})$, correspond to a double embedded eigenvalue at $\\zeta=i$ sitting inside the continuous spectrum. If correct, this settles the conjecture that multiple embedded eigenvalues can occur in the MTM spectral problem, and it connects algebraic solitons to the inverse scattering transform through a singular limit of the Riemann–Hilbert problem.","feed_headline":"Algebraic double-solitons sit at a double embedded eigenvalue","feed_subtitle":"Exact rational solutions of the massive Thirring model realize the predicted Jordan-chain spectrum at ζ=i.","key_machinery":"The Riemann–Hilbert problem for the squared spectral parameter $\\lambda=\\zeta^2$ with reflectionless potential and a double pole at $\\lambda_0=e^{i\\gamma}\\in\\mathbb{C}^+$ is the central object. Residue and double-pole coefficients of the sectionally meromorphic matrix are computed from the scattering data and the symmetry $\\lambda_0\\to\\bar\\lambda_0$, producing a closed linear system whose Cramer's-rule solution gives explicit formulas for $u$ and $v$. The singular limit $\\gamma\\to\\pi$ with rescaled translation parameters $\\tilde x_0\\to\\tilde x_0\\epsilon^2$ is the mechanism that converts the double isolated eigenvalue into the double embedded eigenvalue at $\\zeta=i$, and the eigenvector and generalized eigenvector are obtained as suitably rescaled limits of $\\psi^{(+)}_1(\\zeta_0)$ and its derivative.","core_discovery":"The central discovery is the explicit algebraic double-soliton solution (2.20)–(2.21) and the proof that its Lax spectrum has a double embedded eigenvalue at $\\zeta=i$: one eigenvector $\\psi_0\\in H^1(\\mathbb{R},\\mathbb{C}^2)$ and one generalized eigenvector $\\psi_1\\in H^1(\\mathbb{R},\\mathbb{C}^2)$ satisfying (2.18)–(2.19) with $\\zeta_0=i$. The paper obtains this by first constructing exponential double-solitons from a reflectionless Riemann–Hilbert problem with a quadruplet of double poles at $\\lambda_0=e^{i\\gamma}$ and $\\bar\\lambda_0$, then taking $\\gamma\\to\\pi$ with a rescaling of the translation parameters so that the two double isolated eigenvalues coalesce into a symmetric pair of double embedded poles on the imaginary axis. The resulting rational functions are shown to reduce to the algebraic double-solitons found previously by the bilinear method, and the decay rates $\\psi_0=O(|x|^{-2})$, $\\psi_1=O(|x|^{-1})$ verify the criterion for embedded eigenvalues of higher algebraic multiplicity.","pith_inferences":["A natural extension is that double poles of higher order in the same limit should produce algebraic multi-solitons with longer Jordan chains at $\\zeta=i$, although this paper only constructs chains of length two.","One can test whether the double embedded eigenvalue at $\\zeta=i$ is structurally unstable under a generic perturbation, in analogy with the known instability of a pair of simple embedded eigenvalues that opens into a quadruplet of isolated eigenvalues.","Since the derivative NLS equation and the MTM share the same spectral problem in characteristic coordinates, the same singular-limit construction suggests that double embedded eigenvalues and algebraic double-solitons should also exist for derivative NLS on the zero background."],"forward_implications":["The algebraic double-soliton gives the first concrete rational solution of MTM on the zero background that realizes a double embedded eigenvalue in the Lax spectrum.","The exponential double-soliton describes two identical solitons whose separation grows as $\\log|t|/\\sin\\gamma$, the slow logarithmic dynamics associated with double eigenvalues.","The same Riemann–Hilbert construction, with $\\lambda_0=e^{i\\gamma}$ and the $\\gamma\\to\\pi$ limit, produces algebraic double-solitons without needing a direct RH formulation at an embedded eigenvalue.","Because the eigenvector and generalized eigenvector both lie in $H^1$, the algebraic double-soliton satisfies the decay criterion for higher-multiplicity embedded eigenvalues, so perturbative and stability analysis of embedded eigenvalues can now be tested on an exact solution."],"supporting_citations":[{"why":"states the decay criterion and the conjecture that higher-multiplicity embedded eigenvalues may exist in the MTM Lax spectrum.","marker":"[15]"},{"why":"constructs the algebraic double-soliton by the bilinear method, the solution whose spectral meaning is established here.","marker":"[8]"},{"why":"provides the tau-function bilinear formulation that independently reproduces the exponential double-soliton formula.","marker":"[4]"},{"why":"sets up the inverse scattering transform and potential-recovery formulas used to define the Riemann–Hilbert problem.","marker":"[9]"},{"why":"proves analyticity and the Riemann–Hilbert formulation for the MTM Jost functions in the squared spectral parameter.","marker":"[23]"},{"why":"introduces the double-eigenvalue mechanism that produces exponential double-solitons for the focusing NLS equation.","marker":"[34]"},{"why":"shows for the derivative NLS equation that algebraic double-solitons arise from exponential ones in the singular limit.","marker":"[33]"}],"fun_headline_variants":["Double embedded eigenvalues finally seen in Thirring model solitons","Algebraic solitons unlock double embedded eigenvalue puzzle","Massive Thirring model: double embedded eigenvalues solved","Conjecture resolved: algebraic double-solitons show Jordan-chain spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's conclusion that $\\zeta=i$ is exactly a double embedded eigenvalue with only one eigenvector stands on a previously proved criterion about how fast the generalized eigenvectors decay at infinity; the paper verifies those decay rates, but does not re-derive the full criterion or exclude other independent $H^1$ eigenvectors at the same spectral point.","fun_headline_variants_meta":{"raw":{"variants":["Double embedded eigenvalues finally seen in Thirring model solitons","Algebraic solitons unlock double embedded eigenvalue puzzle","Massive Thirring model: double embedded eigenvalues solved","Conjecture resolved: algebraic double-solitons show Jordan-chain spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000261,"raw_usage":{"total_tokens":1586,"prompt_tokens":933,"completion_tokens":653,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":549,"completion_tokens_details":{"reasoning_tokens":584}},"tokens_in":549,"tokens_out":653,"duration_ms":6289,"temperature":1.0,"reasoning_tokens":584,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:58:06.830733+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the dimension of the $H^1$ solution space of $(\\partial_x-L(u_{\\rm alg},v_{\\rm alg},i))\\psi=0$ and the length of the Jordan chain at $\\zeta=i$. If the kernel contains a second independent $H^1$ eigenvector, or if the generalized eigenvector equation has no $H^1$ solution with $\\psi_0=O(|x|^{-2})$, then the double embedded eigenvalue claim fails; this can be checked by direct substitution of (2.20)–(2.21) into (2.18)–(2.19) or by numerically evaluating the Evans function on the imaginary axis.","supporting_citations":[{"cited_title":"Klaus, D.E","cited_arxiv_id":null,"evidence_quote":"states the decay criterion and the conjecture that higher-multiplicity embedded eigenvalues may exist in the MTM Lax spectrum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"constructs the algebraic double-soliton by the bilinear method, the solution whose spectral meaning is established here."},{"cited_title":"Chen and B.-F","cited_arxiv_id":null,"evidence_quote":"provides the tau-function bilinear formulation that independently reproduces the exponential double-soliton formula."},{"cited_title":"Massive Thirring Model: Inverse Scattering and Soliton Resolution","cited_arxiv_id":"2307.15323","evidence_quote":"sets up the inverse scattering transform and potential-recovery formulas used to define the Riemann–Hilbert problem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"proves analyticity and the Riemann–Hilbert formulation for the MTM Jost functions in the squared spectral parameter."},{"cited_title":"Zakharov and A.B","cited_arxiv_id":null,"evidence_quote":"introduces the double-eigenvalue mechanism that produces exponential double-solitons for the focusing NLS equation."},{"cited_title":"Zhang, L","cited_arxiv_id":null,"evidence_quote":"shows for the derivative NLS equation that algebraic double-solitons arise from exponential ones in the singular limit."}],"review_version":1}