{"id":"5d960665-b0d0-4bed-9749-afab7bee84c3","arxiv_id":"2505.00485","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Half-line Dirac operators with L2 potentials admit an explicit, two-sided weighted stability estimate between potentials and their Schur spectral functions, proved via an exact Kronig-Penney model and Schur's algorithm.","lead":"This paper proves the first explicit two-sided stability estimate for the spectral correspondence of half-line Dirac operators with square-integrable potentials. The proof relies on solving an exact Kronig-Penney model with delta interactions using Schur's algorithm, then approximating general potentials by such models.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Proof of Theorem 1.3 as printed does not establish (1.9): Section 3.4 uses r=e^{-Aℓ} giving weights A/2 instead of A, and Lemma 3.10's proof claims an ℓ-uniform Lipschitz bound that is false by a factor 1/ℓ.","rationale":"The reader's verdict of CONDITIONAL is appropriate: the theorem is very likely true, but the printed proof does not derive it. The reader identified the scaling mismatch in Section 3.4, which is correct and easily fixable by setting r=e^{-2Aℓ}. However, the reader's stated weakest assumption (reliance on external facts from Denisov) is not the load-bearing concern; those are established theorems. The more serious internal gap is the invalid Lipschitz estimate in Lemma 3.10, which is used to justify the Riemann-sum limit connecting the discrete Schur coefficients to the continuous Fourier norms. Even after fixing the weight, Lemma 3.10's proof remains incomplete. Since both gaps are repairable and the surrounding structure (Theorem 3.1, the Kronig-Penney solution, and the non-uniform continuity results) is sound, the verdict should remain CONDITIONAL rather than moving to REJECT or UNVERDICTED.","tokens_in":48444,"tokens_out":27331,"duration_ms":264498,"concrete_test":"For g_ℓ = χ_{[-π/(2ℓ), π/(2ℓ)]} and h_ℓ = √(2π) hat g_ℓ, compute |h_ℓ(2ℓk+δ)-h_ℓ(2ℓk)| for small δ; show it is (π/(2ℓ))δ||g_ℓ||_1 + o(δ), so no ℓ-independent C can satisfy the Claim in Lemma 3.10. Independently, re-run Section 3.4 with r=e^{-2Aℓ} and verify that the limits reproduce exactly (1.9) with c1=√(π/2), c2=2√(2π).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is Theorem 1.3, whose proof in Section 3.4 has two independent gaps. First, the exponential weight mismatch: with r=e^{-Aℓ}, the left-hand sum tends to the L1_A norm and the right-hand sum (as stated in Lemma 3.10) tends to the W1_{A/2} norm, neither of which appears in (1.9). Replacing r by e^{-2Aℓ} repairs the weights and yields the constants c1=√(π/2), c2=2√(2π) after taking limits. Second, Lemma 3.10 contains a false estimate: h=√(2π)(hat g_{qℓ}-hat g_{q̃ℓ}) satisfies |h'(ξ)| ≤ ||x(g_{qℓ}-g_{q̃ℓ})||_1 ≤ (π/(2ℓ))||g_{qℓ}-g_{q̃ℓ}||_1, since the supports have length π/ℓ. The claimed bound |h(ξ)-h(2ℓk)| ≤ C|ξ-2ℓk| ||g_{qℓ}-g_{q̃ℓ}||_1 with an absolute constant C is therefore false; the true bound carries a factor 1/ℓ. Consequently the O(ℓ^{1/2}) estimate for (3.37) does not follow, and the limit passage from the discrete sums to the continuous Fourier norms is not justified. The external theorems of Denisov cited in the paper are not the weakest point; the internal proof of the approximation step is.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the spectral correspondence q ↦ f_q between L2 potentials on the half-line and Schur functions in S2(C+) for Dirac operators. The main new result, Theorem 1.3, claims an explicit two-sided uniform estimate with constants c1 = √(π/2), c2 = 2√(2π) between the weighted L1 norm of the potential difference and the Wiener-type W1_A norm of the difference of Schur functions, for A ≥ 12 max(||q||², ||q̃||²). The proof strategy is to approximate general L2 potentials by discrete measures supported on ℓZ+, to solve the direct and inverse spectral problems for the resulting Kronig-Penney models exactly via Schur's algorithm, to establish quantitative continuity estimates for Schur iterates (Theorem 3.1), and then to pass to the limit ℓ → 0 using approximation lemmas. The paper also proves Theorem 1.2 on the failure of uniform continuity in both directions of the Sylvester-Winebrenner homeomorphism, and derives several corollaries for the relativistic Kronig-Penney model. The central announced result therefore depends on the discrete approximation machinery in Section 3 and on the limiting passage in Lemma 3.10 and Section 3.4.","tokens_in":48694,"tokens_out":7523,"duration_ms":70057,"significance":"If Theorem 1.3 is correct, it would be the first explicit two-sided uniform quantitative stability estimate for the Sylvester-Winebrenner spectral correspondence for half-line Dirac operators with general L2 potentials. The paper combines a substantial exact-solvable model (Dirac operators with δ-interactions on a half-lattice) with Schur's algorithm and gives explicit constants. It also proves a non-uniformity theorem in both directions, which is a valuable cautionary complement to the positive estimate. The authors provide detailed proofs of many auxiliary statements, and the use of the external bijection theorem and sum rule from Denisov's work is legitimate rather than circular. However, the proof of the main theorem as printed contains two distinct technical gaps: a weight mismatch in the limiting step and a false uniform Lipschitz estimate in Lemma 3.10. These gaps occur precisely at the load-bearing passage from the discrete ℓ-dependent estimates to the continuous Fourier norms in (1.9).","major_comments":[{"comment":"The choice r = e^{-Aℓ} is incompatible with Lemma 3.10. Lemma 3.10 asserts a limit with weight e^{-Aξ/2} on the right-hand side, whereas Theorem 1.3 and (3.39) require the weight e^{-Aξ}. With r = e^{-Aℓ}, equation (3.36) shows that the sum converges to (1/√(2π))∫_{R+} e^{-Aξ/2}|f̂_q - f̂_q̃| dξ, not to the W1_A norm. Consequently the limiting inequality obtained from (3.43) bounds W1_{A/2}, not W1_A, and (1.9) is not established as printed. Replacing r by e^{-2Aℓ} would align the weights: the Riemann sum then tends to (1/√(2π))∫ e^{-Aξ}|f̂_q - f̂_q̃| dξ, and combining (3.43) with Lemma 3.7 with weight 2A yields exactly the constants c1 = √(π/2) and c2 = 2√(2π) after ε → 0. This repair is local, but as it stands the printed proof does not prove the stated theorem.","section":"Section 3.4, Eq. (3.39)-(3.43)"},{"comment":"The claimed uniform Lipschitz bound |h(ξ) - h(2ℓk)| ≤ C|ξ - 2ℓk| · ||g_{qℓ} - g_{q̃ℓ}||_{L1} is false. Since h is the Fourier transform of g_{qℓ} - g_{q̃ℓ}, one has |h'(ξ)| ≤ (1/√(2π)) ∫ |x| |g_{qℓ}(x) - g_{q̃ℓ}(x)| dx ≤ (π/(2ℓ)) ||g_{qℓ} - g_{q̃ℓ}||_{L1}, because the supports of g_{qℓ} and g_{q̃ℓ} are contained in an interval of length π/ℓ. Using this correct bound, each term in the sum in (3.38) is of size O(ℓ^{1/2} e^{-Aℓk}) rather than O(ℓ^{3/2} e^{-Aℓk}), and the total error becomes O(ℓ^{-1/2}) instead of O(ℓ^{1/2}). The convergence claimed in (3.37) is therefore not justified by the argument given. This is a load-bearing step in the discrete-to-continuous limit, so the proof of Lemma 3.10 needs a substantially different argument or an additional structural property of the functions g_{qℓ}.","section":"Lemma 3.10, proof of (3.37)-(3.38)"}],"minor_comments":[{"comment":"In the displayed line preceding (3.37), the expression \"2r|ĝ_{qℓ}(2ℓk) - ...|\" should read \"2ℓ|...|\"; the variable r is not defined at that point and the intended Riemann sum step is 2ℓ.","section":"Section 3.3, Eq. (3.36)"},{"comment":"In the sentence \"where the factor 2π appears in the last inequality\", the reference should be to the last equality; the normalization of m_T accounts for a factor 2π, not an inequality.","section":"Lemma 3.8, proof"},{"comment":"The theorem statement and (3.39) assume A ≥ 12 max(||q||², ||q̃||²), which forces A > 0, but (3.39) is introduced with \"A ∈ R\". This should be corrected to A > 0 for clarity.","section":"Section 3.4, opening paragraph"},{"comment":"There are several typographical issues in the displayed formulas, e.g., stray OCR artifacts such as \"/greaterorequalslant\" in place of ≥ and \"OPERA TORS\" in the title header. These do not affect the mathematics but should be cleaned up in the final version.","section":"Throughout"},{"comment":"The proof of Proposition 4.14 relies on Example 6.1 from [9] for the existence of non-unique reflection coefficients; this is an acceptable external input, but the dependence should be stated explicitly in the main text so that the reader understands that Theorem 1.2 inherits the non-injectivity result from [9].","section":"Section 4.4, Proposition 4.12"}],"recommendation":"major_revision","confidential_remarks":"The two gaps identified in Section 3.4 and Lemma 3.10 are technical but load-bearing; the weight mismatch is straightforward to repair, while the false Lipschitz bound suggests that the approximation step needs a genuinely new argument. I do not see grounds for rejection, because the overall strategy is coherent and the remaining parts of the paper appear carefully executed, but the current arXiv version does not yet contain a correct proof of the advertised main theorem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Main take: this paper has two genuinely new pieces—an explicit two-sided stability estimate for the half-line Dirac spectral correspondence in the L2 class, and an exact description of the massless Dirac Kronig-Penney model via Schur's algorithm. The second is solid and elegant. The first is plausible and important, but the printed proof has a gap that is not just a typo.\n\nThe good parts: Theorem 1.4 and its corollaries are well proved. The reduction of measure-valued Dirac operators on a lattice to Schur's algorithm works, and the applications (Rakhmanov, Baxter, Szegő-Golinskii-Ibragimov analogues) are nice. Section 4 on non-uniform continuity of the Sylvester-Winebrenner map contains explicit constructions and reads correctly. The historical remarks are accurate and useful.\n\nThe problem is in the proof of Theorem 1.3, specifically Lemma 3.10. Two issues. First, the scaling: with r=e^{-Aℓ}, the conclusion has e^{-Aξ/2}, not e^{-Aξ}, so it does not match (1.9). That can be fixed by taking r=e^{-2Aℓ}. Second and more serious, the proof of Lemma 3.10 claims an ℓ-uniform Lipschitz bound for h=√(2π)(ĝ_qℓ−ĝ_q̃ℓ). But h is the Fourier transform of a function supported on an interval of length π/ℓ, so its derivative is bounded by (π/ℓ) times the L1 norm, not by an absolute constant times that L1 norm. The factor 1/ℓ changes the Riemann-sum error from O(ℓ^{1/2}) to O(1/√ℓ), which diverges. As printed, the passage from discrete Fourier coefficients to the continuous Fourier norm is unjustified. This is the load-bearing step for Theorem 1.3. I do not see a way to save it with a small edit; it needs a different argument, possibly using the special structure of Schur functions or a different approximation scheme.\n\nThe theorem may well be true—the strategy is natural and the estimate has the right scaling—but the current manuscript does not prove it. The weakest point is not the reliance on Denisov's bijection and sum rule; those are external and legitimate. It is the internal limiting argument in Lemma 3.10.\n\nWho is this for? Spectral theorists working on Dirac operators, inverse scattering, and OPUC. The Kronig-Penney part alone is worth a serious referee. I would send it to referees, with instructions that the proof of Theorem 1.3 needs significant repair before publication.","headline":"Solid Kronig-Penney result and a plausible main claim, but the proof of Theorem 1.3 has a gap in Lemma 3.10 that is more than a typo; the paper deserves serious refereeing but needs repair.","tokens_in":49313,"tokens_out":12970,"would_cite":true,"duration_ms":125593,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34L40"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves the first explicit two-sided uniform stability estimate for the Dirac spectral correspondence with $L^2$ potentials.","keywords":["Dirac operators","spectral continuity","Schur's algorithm","Kronig–Penney model","inverse spectral problem","nonlinear Fourier transform","Wiener algebra","delta-interactions"],"falsifier":"Take two compactly supported piecewise-constant potentials $q,\\tilde q$ with known matrix-exponential solutions, compute $f_q,f_{\\tilde q}$ by solving the Dirac system, choose $A=12\\max(\\|q\\|_{L^2}^2,\\|\\tilde q\\|_{L^2}^2)$, and check whether the ratio $\\|f_q-f_{\\tilde q}\\|_{W^1_A}/\\|q-\\tilde q\\|_{L^1_{2A}}$ lies between $\\sqrt{\\pi/2}$ and $2\\sqrt{2\\pi}$. A single pair with ratio outside this interval, computed to numerical tolerance, would refute Theorem 1.3's constants; a search over random piecewise-constant pairs would also test whether the stated threshold on $A$ is needed.","tokens_in":48150,"feed_emoji":"📐","tokens_out":7536,"duration_ms":73202,"temperature":0.7,"pith_summary":"The paper establishes that the half-line Dirac spectral correspondence—mapping a square-integrable potential $q$ to its Schur function $f_q$—has explicit, uniform, two-sided continuity estimates. Concretely, for any two $L^2$ potentials $q,\\tilde q$, the weighted $L^1$ distance between the potentials and the weighted Wiener-algebra distance between their spectral data are comparable up to universal constants, provided the weight parameter is large enough relative to the two norms. Previous knowledge of this Sylvester–Winebrenner correspondence was qualitative: it is a homeomorphism, but not uniformly continuous. The paper supplies the first uniform bound with explicit constants, and it is new even when one potential is zero. The proof works by solving the inverse spectral problem exactly for Dirac operators with $\\delta$-interactions on a lattice via Schur's algorithm, then approximating arbitrary $L^2$ potentials by such discrete measures.","feed_headline":"Dirac spectral map obeys explicit two-sided stability bounds","feed_subtitle":"For L2 potentials, weighted distances to spectral Schur functions match up to universal constants.","key_machinery":"The engine is Schur's algorithm. For a Schur function $F$ in the unit disk, the iteration $(zF_{k+1})(z)=(F_k(z)-F_k(0))/(1-\\overline{F_k(0)}F_k(z))$ produces recurrence coefficients $F_k(0)$; for a periodic Schur function in the upper half-plane the same iteration reduces to a sequence of disk recurrence coefficients. The paper proves (Theorem 1.4) that for a Dirac potential supported on the lattice $\\ell\\mathbb Z_+$, the point mass at $\\ell k$ is exactly $\\kappa(F_k(0))$, where $\\kappa(w)=\\frac{w}{2|w|}\\log\\frac{1+|w|}{1-|w|}$. This exact dictionary turns the spectral correspondence for Kronig–Penney-type potentials into a two-sided estimate for Schur's algorithm in weighted Wiener algebras (Theorem 3.1). General $L^2$ potentials are then approximated by lattice measures with $\\ell\\to 0$, with convergence of potentials and spectral data controlled by lemmas proved in Sections 3.2–3.3.","core_discovery":"Theorem 1.3 states that for $q,\\tilde q\\in L^2(\\mathbb R_+)$ with Schur functions $f_q,f_{\\tilde q}$, one has $$ c_1\\|q-\\tilde q\\|_{$L^{1}$_{2A}(\\mathbb R_+)} \\le \\|f_q-f_{\\tilde q}\\|_{$W^{1}$_A(\\mathbb R_+)} \\le c_2\\|q-\\tilde q\\|_{$L^{1}$_{2A}(\\mathbb R_+)}, $$ with $c_1=\\sqrt{\\pi/2}$, $c_2=2\\sqrt{2\\pi}$, and $A\\ge 12\\max(\\|q\\|^2_{L^2},\\|\\tilde q\\|^2_{L^2})$. Here the weighted $L^1$ norm is $\\int_{\\mathbb R_+}|q(\\xi)|e^{-A\\xi}d\\xi$ and the spectral norm is the $L^1$ norm of the Fourier transform of $f_q-f_{\\tilde q}$ with weight $e^{-A\\xi}$. The paper also proves that on bounded sets this homeomorphism is not uniformly continuous in either direction (Theorem 1.2), so the two-sided weighted estimate is the sharp type of uniform control that holds.","pith_inferences":["My inference: the same Schur-algorithm route could yield explicit stability estimates for other exactly solvable one-dimensional models, such as Schr\\\"odinger operators with point interactions, by following the same two-step approximation; the paper does not claim this.","My inference: the constants $\\sqrt{\\pi/2}$ and $2\\sqrt{2\\pi}$ are unlikely to be sharp; the proof passes through an $\\varepsilon$-approximation and a general Schur-theorem estimate with lossy factors, so I would expect the optimal ratio to be smaller and possibly expressible in terms of the $\\kappa$ map alone.","My inference: a natural testable extension is to replace the exponential weights by polynomial weights and check whether a two-sided estimate of the same form holds with constants depending on the polynomial degree; the structure of the proof does not immediately preclude or imply this."],"forward_implications":["The direct problem is stable with rate: if one potential changes slightly in weighted $L^1$, the spectral Schur function changes proportionally in weighted Wiener norm, independent of the pair.","The inverse problem inherits the same bound: spectral data close in the weighted Wiener metric force potentials close in the weighted $L^1$ metric.","Approximation by lattice potentials converges quantitatively: the discrete measures $q_\\ell$ in (3.17) give spectral data whose loss is controlled as $\\ell\\to0$, so the exactly solvable model is a usable numerical and discretization scheme.","The result supplies explicit constants for the nonlinear Fourier transform associated with the massless Dirac operator, so it bears on quantitative stability of inverse scattering for nonlinear Schr\\\"odinger-type problems.","Uniform continuity in the original unweighted metrics fails (Theorem 1.2); the weighted norms in Theorem 1.3 identify the correct quantitative framework."],"supporting_citations":[{"why":"Supplies the bijection $L^2(\\mathbb R_+)\\to S_2(\\mathbb C_+)$ and the sum rule (1.5) that the proof uses in Lemma 3.8 and Section 3.4.","marker":"[15]"},{"why":"Introduces the Sylvester–Winebrenner spectral correspondence whose continuity the paper makes quantitative.","marker":"[54]"},{"why":"Provides the Schur-algorithm and orthogonal-polynomial facts, including recurrence coefficients determining $F$ and Fourier-coefficient convergence, used in the proof of Theorem 3.1 and in the lattice dictionary.","marker":"[51]"},{"why":"Sets up Dirac operators with measures and their reduction to canonical Hamiltonian systems, used to prove Theorem 1.4 for lattice potentials.","marker":"[61]"},{"why":"Gives the Krein–de Branges theorem for canonical systems invoked in the approximation step of Lemma 3.9.","marker":"[47]"},{"why":"Supplies the solution theory for measure differential equations that defines Dirac operators with $\\delta$-interactions.","marker":"[44]"}],"fun_headline_variants":["Explicit two-sided stability for Dirac spectral map","Quantitative spectral continuity for Dirac operators","Dirac spectra: sharp weighted stability bounds","Two-sided weighted estimate for Dirac spectral correspondence","Uniform control of Dirac spectral data with explicit constants"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument relies on the previously established theorem that the map $q\\mapsto f_q$ is a bijection from $L^2(\\mathbb R_+)$ onto the class of Schur functions with integrable $\\log(1-|f|^2)$, together with the sum rule $\\|q\\|_{L^2}^2=\\frac1\\pi\\int -\\log(1-|f_q|^2)$; if either fails, the limiting step that transfers the lattice estimates to general $L^2$ potentials would not go through.","fun_headline_variants_meta":{"raw":{"variants":["Explicit two-sided stability for Dirac spectral map","Quantitative spectral continuity for Dirac operators","Dirac spectra: sharp weighted stability bounds","Two-sided weighted estimate for Dirac spectral correspondence","Uniform control of Dirac spectral data with explicit constants"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1468,"prompt_tokens":873,"completion_tokens":595,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":489,"completion_tokens_details":{"reasoning_tokens":528}},"tokens_in":489,"tokens_out":595,"duration_ms":6717,"temperature":1.0,"reasoning_tokens":528,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:43:02.876604+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take two compactly supported piecewise-constant potentials $q,\\tilde q$ with known matrix-exponential solutions, compute $f_q,f_{\\tilde q}$ by solving the Dirac system, choose $A=12\\max(\\|q\\|_{L^2}^2,\\|\\tilde q\\|_{L^2}^2)$, and check whether the ratio $\\|f_q-f_{\\tilde q}\\|_{W^1_A}/\\|q-\\tilde q\\|_{L^1_{2A}}$ lies between $\\sqrt{\\pi/2}$ and $2\\sqrt{2\\pi}$. A single pair with ratio outside this interval, computed to numerical tolerance, would refute Theorem 1.3's constants; a search over random piecewise-constant pairs would also test whether the stated threshold on $A$ is needed.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the bijection $L^2(\\mathbb R_+)\\to S_2(\\mathbb C_+)$ and the sum rule (1.5) that the proof uses in Lemma 3.8 and Section 3.4."},{"cited_title":"Sylvester and D","cited_arxiv_id":null,"evidence_quote":"Introduces the Sylvester–Winebrenner spectral correspondence whose continuity the paper makes quantitative."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Schur-algorithm and orthogonal-polynomial facts, including recurrence coefficients determining $F$ and Fourier-coefficient convergence, used in the proof of Theorem 3.1 and in the lattice dictionary."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets up Dirac operators with measures and their reduction to canonical Hamiltonian systems, used to prove Theorem 1.4 for lattice potentials."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Krein–de Branges theorem for canonical systems invoked in the approximation step of Lemma 3.9."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the solution theory for measure differential equations that defines Dirac operators with $\\delta$-interactions."}],"review_version":1}