{"id":"2a641452-920f-41ca-bef9-8d73ee19a9d9","arxiv_id":"2505.00669","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The direct spectral problem for Paley-Wiener canonical systems is solved by recovering the spectral measure from step-function Hamiltonian data, with a convergence result for real Dirac systems.","lead":"This paper gives algorithms to compute the spectral measure of a canonical Hamiltonian system from the Hamiltonian, first for step-function Hamiltonians and then for smooth systems via step-function approximation. The convergence theorem is proved for real Dirac systems, giving a practical tool for direct spectral problems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 5.1's proof is non-rigorous: the energy estimate treats complex E−E_T as real, so the μ_T→μ convergence claim is unsupported as written.","rationale":"The reader's weakest_assumption focused on the restriction to real Dirac Hamiltonians, which is a genuine scope limitation: the abstract's convergence claim is broader than the theorem proved. However, I find a more immediate load-bearing problem inside the proof of Theorem 5.1 itself. The proof's differential inequality for D = E − E_T is invalid for complex D: writing (1/2)∂_t D² is not the derivative of |D|², and the subsequent inequalities therefore do not imply the claimed uniform convergence of |E_T(a,·)|² to |E(a,·)|². Because Theorem 5.1 is the only result establishing μ_T → μ for non-step Hamiltonians, the central claim of the paper is not rigorously supported as written. The flaw appears repairable by a standard energy estimate using |D|², so the appropriate disposition remains CONDITIONAL rather than REJECT. Even after such a repair, the real-Dirac restriction should be reflected in the abstract. Thus I do not change the reader's verdict, but I flag the proof step as the most load-bearing concern.","tokens_in":14420,"tokens_out":16272,"duration_ms":154400,"concrete_test":"Re-derive the key estimate rigorously: for D = E_T − E, use the correct real-axis relation E^# = \\overline{E} to compute (|D|²)' = 2 Re(f D^# \\bar{D}) + 2 Re((f_T−f)E_T^# \\bar{D}), bound this by 2|f||D|² + 2|f_T−f||E_T||D|, and apply Grönwall to obtain a uniform-in-x bound on |D(a,x)|² in terms of ∫_0^a |f_T−f| dt. If this bound holds with a constant independent of x, Theorem 5.1 is repairable; if the oscillatory factor e^{2itx} produces an unbounded term in x, the theorem as stated needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5, proof of Theorem 5.1. The paper aims to prove uniform convergence of |E_T(a,·)|² to |E(a,·)|² by estimating D = E − E_T. Two related problems appear. First, equation (17) states ∂_t E = f e^{2itx}E^# = f e^{2itx}E; on the real axis E^#(t,x) = \\overline{E(t,x)}, not E(t,x), and for the scattering function the correct real-axis identity is ∂_t E = f e^{2itx}E^#, not f e^{2itx}E. Second, the proof multiplies the complex differential equation for D by D and writes (1/2)∂_t D² = f e^{2itx}|D|² − [f_T−f]e^{2itx}E_T D. This is not the derivative of |D|², and no inequality of the form |∂_t D²| ≤ |f||D|² + C|f_T−f| follows for complex D. Consequently the integrated estimate |D(a,x)|² ≤ C ∫|f_T−f| dt, which is the justification for uniform convergence of |E_T|², is not established. Since Theorem 5.1 is the only argument that μ_T → μ for non-step Hamiltonians, the paper's central convergence claim lacks a rigorous proof as written. A standard energy estimate for Y=|D|² appears to repair the theorem, so the flaw is likely fixable rather than fatal, but the current text does not contain that argument. Separately, even if the proof is repaired, the abstract's unqualified statement about general non-step Hamiltonians is unsupported, since Theorem 5.1 is restricted to real Dirac systems.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This note studies direct spectral problems for Paley–Wiener canonical Hamiltonian systems on the half-line. For diagonal Hamiltonians that are step functions with uniform step size, the spectral measure is an even periodic Paley–Wiener measure, and the paper gives two algebraic algorithms to recover it from the sequence of step values h_11^n: one via Verblunsky coefficients (Section 3) and one via moments and Toeplitz determinants (Section 4). For non-step Hamiltonians, the paper defines step-function approximations H_T by averaging h_11 on intervals and claims in Theorem 5.1 that, for Hamiltonians arising from real Dirac systems, the associated spectral measures μ_T converge to μ in the sense that ∫|φ|^2 dμ_T → ∫|φ|^2 dμ for every φ in a Paley–Wiener space. The final section contains numerical illustrations of the algorithms and explicit limiting measures.","tokens_in":14756,"tokens_out":15787,"duration_ms":146423,"significance":"The step-function part is a clean application of OPUC theory that yields explicit, finite-dimensional formulas and complements the inverse spectral algorithm of Makarov–Poltoratski. If Theorem 5.1 is established, the paper would provide a practical periodization method for direct spectral problems in a nontrivial subclass of PW-systems, and the examples demonstrate that the method produces meaningful approximations. However, the main convergence theorem currently rests on an incorrect differential identity and an invalid energy estimate, so the contribution is not yet rigorously established.","major_comments":[{"comment":"The identity ∂_t E(t,x)=f(t)e^{2itx}E(t,x) is false on the real axis because E^#(t,x)=overline{E(t,x)}, not E(t,x). With the scattering function defined as E(t,z)=e^{itz}E(t,z), the correct real-axis equation is ∂_t E(t,x)=f(t)e^{2itx}overline{E(t,x)}. This error invalidates the differential equation for D=E−E_T that drives the proof of Theorem 5.1.","section":"Section 5, Eq. (17)"},{"comment":"The energy estimate treats the complex-valued D as real. The displayed identity (1/2)∂_t[D]^2=f e^{2itx}|D|^2−[f_T−f]e^{2itx}E_T D is not the derivative of |D|^2, and the subsequent inequality |∂_t D^2|≤... does not follow. The argument can likely be repaired by multiplying by \\bar{D} and using ∂_t|D|^2=2 Re(\\bar{D}∂_t D), but that argument is absent.","section":"Section 5, proof of Theorem 5.1"},{"comment":"The paper states that f_T is a discrete measure and asserts ∫_0^a |f_T(s)|ds → ∫_0^a |f(s)|ds and ∫_0^a |f_T−f|dt → 0 as T→0. Since f_T is defined implicitly through the logarithmic jumps of the averages of h_11, these convergence statements are not automatic and no proof is provided. The final estimate in Theorem 5.1 requires such an L^1-convergence statement.","section":"Section 5, paragraph after (14)"},{"comment":"The abstract claims convergence of the spectral measures for 'a non-step-function Hamiltonian' without the restriction to real Dirac systems. Theorem 5.1 is proved only for real Dirac Hamiltonians, and the introduction itself states that the class of PW-Hamiltonians is not characterized. The abstract and the opening paragraphs should be amended to reflect the actual scope.","section":"Abstract and Section 1"}],"minor_comments":[{"comment":"The definition of H^2(C+) is missing the dx in the integral; it should read sup_{y>0} ∫_R |f(x+iy)|^2 dx.","section":"Section 2.1"},{"comment":"The symbol ⊮ is used without definition; it appears to be an OCR artifact for a row vector of ones. Please define it or replace it with standard notation.","section":"Section 4.2 and 4.3"},{"comment":"The sentence contains a typo: 'as T → ∞' should be 'as T → 0' in the sentence about ∫_0^t |f_T(s)|ds.","section":"Section 5, sentence after (14)"},{"comment":"The formula for w_T(θ) should include absolute values around sin(Tθ), and the support interval should reflect the absolute value of the arcsine. Also, 'dTµ(x)' should be 'dμ_T(x)'.","section":"Example 6.5"},{"comment":"The derivation of (4) uses φ_n^*(1)=φ_n(1), which holds for even measures; this justification could be stated explicitly for the reader.","section":"Section 3.2, Eq. (4)"}],"recommendation":"major_revision","confidential_remarks":"The algebraic sections (3 and 4) appear correct and are valuable. The main concern is the rigor of Section 5; the proof of Theorem 5.1 needs to be rewritten. If the author can supply a correct energy estimate and clarify the f_T convergence, the paper would be suitable for publication. The abstract should be revised to avoid overclaiming. I recommend major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me give you a quick read of Zhang's paper. The genuinely useful part is the explicit finite-dimensional algorithm for the direct spectral problem for step-function Hamiltonians: from the h11 values you recover the Verblunsky coefficients by the Szegő recurrence, or the moments by Toeplitz determinant formulas. That part is coherent, standard, and checks out on the examples. The paper also tries to go beyond step functions by approximating a general Hamiltonian by step functions and proving convergence of the spectral measures. That theorem, 5.1, is the load-bearing new result, and here the paper has a real problem.\n\nThe proof of Theorem 5.1 is not rigorous as written. Equation (17) asserts that on the real axis E^#(t,x)=E(t,x). But for a scattering function, E^#(t,x) = \\overline{E(t,x)}, not E(t,x), and for a real Dirac system the scattering function is not real-valued on the real axis. Then the energy estimate multiplies the complex differential equation for D=E−E_T by D and treats (1/2)∂_t D^2 as a real energy. That step is invalid for complex D, so the estimate |D|^2 ≤ C∫|f_T−f| does not follow. The convergence theorem may well be true, and a standard estimate for |D|^2 would likely repair it, but the text does not contain that argument. Because the abstract claims convergence for general non-step Hamiltonians without the real Dirac restriction, it overstates the result. Also, Example 6.5 asserts the limiting measure (18) without a derivation, which is minor but annoying.\n\nWhat is good: the algebraic parts in Sections 3 and 4 are careful, the paper is clearly written, and it builds explicitly on standard OPUC theory and on [13], with no circularity. The step-function algorithm is a genuine complement to the inverse problem.\n\nWho should read it: anyone working on canonical systems or inverse spectral problems who wants an explicit direct algorithm. But treat Theorem 5.1 as unproven in the current version.\n\nRecommendation: send to peer review, but the referee should demand a corrected proof of Theorem 5.1 and a revised abstract. This is a fixable paper, not a reject.","headline":"Useful explicit algorithm for step-function Hamiltonians, but the convergence theorem for non-step Hamiltonians has a non-rigorous proof as written.","tokens_in":15254,"tokens_out":4216,"would_cite":true,"duration_ms":37987,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34A55","34L05","47B35"],"pacs":[],"model":"deepseek-v4-flash","headline":"For real Dirac systems, the direct spectral problem is solved by step-function approximation: the spectral measures of the approximating Hamiltonians converge to the true spectral measure on Paley-Wiener test functions.","keywords":["spectral measure","canonical Hamiltonian systems","Paley-Wiener spaces","de Branges spaces","real Dirac systems","step-function approximations","Toeplitz determinants","Verblunsky coefficients"],"falsifier":"For the real Dirac system with $h_{11}(t)=e^t$ (Example 6.5), compute numerically the norms $\\|\\varphi\\|_{L^2(\\mu_T)}$ for a fixed $\\varphi\\in PW_1$ as $T\\to0$ and compare them with the integral against the claimed limiting density $\\sqrt{4x^2-1}/(2|x|)\\,dx$; a systematic discrepancy would show the convergence in Theorem 5.1 fails.","tokens_in":14201,"feed_emoji":"📐","tokens_out":13260,"duration_ms":115135,"temperature":0.7,"pith_summary":"The paper solves the direct spectral problem for a class of Paley-Wiener canonical systems: given the Hamiltonian, find its spectral measure. The method is to approximate a general Hamiltonian by step functions, solve the direct problem for each step-function approximation using orthogonal polynomials and Toeplitz determinants, and then shrink the step size to zero. The central theorem proves that for Hamiltonians coming from real Dirac systems, the spectral measures of the step-function approximations converge to the spectral measure of the original system in the sense of $L^2$ norms on Paley-Wiener test functions. This gives a finite-dimensional algorithm for the direct problem and shows that the truncated-Toeplitz approach to inverse spectral problems can be run in reverse.","feed_headline":"Step-function Hamiltonians capture Dirac-system spectra as T shrinks","feed_subtitle":"For real Dirac systems, spectral measures of step approximations converge to the true measure on Paley–Wiener test functions.","key_machinery":"The argument moves through the Hermite-Biehler functions $E(t,z)=u(t,z)-iv(t,z)$ associated with the canonical system; the spectral measure satisfies the norm identity $\\|\\varphi\\|_{L^2(\\mu)}^2 = \\int_{\\mathbb{R}} |\\varphi(x)|^2\\,dx/|E(a,x)|^2$ for $\\varphi$ in the corresponding de Branges space. For real Dirac systems, $E$ obeys the scattering equation $\\partial_t E(t,x)=f(t)e^{2itx}\\overline{E(t,x)}$, which yields uniform bounds and, when $f$ is replaced by its step average $f_T$, a uniform convergence estimate $E_T(a,x)\\to E(a,x)$ as $T\\to0$. The step-function case is handled algebraically: the Szegő recurrence for orthogonal polynomials on the unit circle and finite Toeplitz determinants recover the moments and Verblunsky coefficients of the periodic spectral measure from the step heights $h_{11,n}=|\\varphi_n(1)|^2$.","core_discovery":"The paper's central result, Theorem 5.1, states that if $H$ is the det-normalized diagonal Hamiltonian of a real Dirac system, with $h_{11}(t)=\\exp(\\int_0^t f(s)ds)$ for a real-valued $f\\in L^1_{\\mathrm{loc}}$, and $H_T$ are the step-function approximations obtained by averaging $h_{11}$ over intervals of length $T$, then for every $a>0$ and every $\\varphi\\in PW_a$, the $L^2(\\mu_T)$ norm of $\\varphi$ converges to its $L^2(\\mu)$ norm as $T\\to0$. This means the spectral measures $\\mu_T$ converge to the true spectral measure $\\mu$ on the Paley-Wiener test functions that characterize PW-systems. The proof establishes uniform convergence of the scattering functions $E_T(a,x)$ to $E(a,x)$, which transfers to the reciprocals $1/|E_T(a,x)|^2$ because $|E(a,x)|$ stays bounded away from zero; Hölder's inequality then gives the norm convergence.","pith_inferences":["We conjecture that Theorem 5.1 extends to any Paley-Wiener Hamiltonian whose Hermite-Biehler functions satisfy a scattering-type equation with an L1-loc coefficient, not only the exact exponential form of real Dirac systems; the paper's proof only needs the uniform bound and the L1 approximation.","The closed-form limit for h11(t)=e^t suggests a more general correspondence between exponential growth rates of the Hamiltonian and the location of spectral support; testing h11(t)=e^{ct} for c different from 1 would probe this.","Because the paper treats only diagonal Hamiltonians and even spectral measures, a natural extension is to apply the same periodization to non-diagonal step-function Hamiltonians, where the Verblunsky coefficients are replaced by matrix-valued recurrences.","The connection to orthogonal polynomials may allow the direct spectral problem to be computed by fast Toeplitz algorithms, making the step-function approximation practical for numerical spectral computations."],"forward_implications":["For any real Dirac system, the spectral measure can be computed to arbitrary accuracy by averaging the Hamiltonian over small intervals and then running the finite-dimensional Toeplitz/orthogonal-polynomial algorithm for periodic measures.","The direct and inverse spectral problems for Paley-Wiener canonical systems become two directions of the same step-function approximation procedure, so a measure and its Hamiltonian can be recovered from each other consistently.","The convergence of the approximate spectral measures holds on every Paley-Wiener space PW_a, the exact class of test functions relevant to the sampling property of PW-measures.","For step-function Hamiltonians with uniform step size, the direct spectral problem reduces to finding the unique even measure on the circle with prescribed values |φ_n(1)|^2, and both the Verblunsky coefficients and the moments can be recovered algebraically.","The examples confirm the method on known cases and produce a closed-form limiting measure for the exponential Hamiltonian h11(t)=e^t."],"supporting_citations":[{"why":"Defines Paley-Wiener measures and systems, and gives the inverse spectral problem algorithm whose reversal yields the direct problem for step-function Hamiltonians.","marker":"[13]"},{"why":"Places real Dirac systems inside the Paley-Wiener Hamiltonian class, justifying the approximation setting of Theorem 5.1.","marker":"[2]"},{"why":"Introduces periodic approximations of spectral measures in inverse problems, the periodization strategy the paper transfers to the direct problem.","marker":"[18]"},{"why":"Supplies the Szegő recurrence, Verblunsky coefficients, and moment formulas used to recover the measure of a step-function Hamiltonian from |φ_n(1)|^2.","marker":"[21]"},{"why":"Provides the inversion formula for finite Toeplitz matrices that underpins the algebraic proof of the moment-recovery algorithm.","marker":"[22]"},{"why":"Gives the conversion of real Dirac systems to canonical systems, yielding the exponential Hamiltonian form used in the convergence proof.","marker":"[20]"},{"why":"Develops the de Branges space theory underlying the notion of spectral measures of canonical systems.","marker":"[6]"}],"fun_headline_variants":["Step-function limits unlock Paley–Wiener spectra","Direct spectral problem via step Hamiltonian limits","T→0 step Hamiltonians converge to true spectra","Canonical spectrum from step-function limits","PW spectra via step Hamiltonian approximation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Hamiltonian has the special exponential form produced by a real Dirac system, and the paper does not establish the same convergence for general Paley-Wiener Hamiltonians.","fun_headline_variants_meta":{"raw":{"variants":["Step-function limits unlock Paley–Wiener spectra","Direct spectral problem via step Hamiltonian limits","T→0 step Hamiltonians converge to true spectra","Canonical spectrum from step-function limits","PW spectra via step Hamiltonian approximation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001145,"raw_usage":{"total_tokens":4708,"prompt_tokens":858,"completion_tokens":3850,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":474,"completion_tokens_details":{"reasoning_tokens":3783}},"tokens_in":474,"tokens_out":3850,"duration_ms":31112,"temperature":1.0,"reasoning_tokens":3783,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:36:51.773333+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the real Dirac system with $h_{11}(t)=e^t$ (Example 6.5), compute numerically the norms $\\|\\varphi\\|_{L^2(\\mu_T)}$ for a fixed $\\varphi\\in PW_1$ as $T\\to0$ and compare them with the integral against the claimed limiting density $\\sqrt{4x^2-1}/(2|x|)\\,dx$; a systematic discrepancy would show the convergence in Theorem 5.1 fails.","supporting_citations":[{"cited_title":"Makarov and A","cited_arxiv_id":null,"evidence_quote":"Defines Paley-Wiener measures and systems, and gives the inverse spectral problem algorithm whose reversal yields the direct problem for step-function Hamiltonians."},{"cited_title":"Bessonov and R","cited_arxiv_id":null,"evidence_quote":"Places real Dirac systems inside the Paley-Wiener Hamiltonian class, justifying the approximation setting of Theorem 5.1."},{"cited_title":"Poltoratski and A","cited_arxiv_id":null,"evidence_quote":"Introduces periodic approximations of spectral measures in inverse problems, the periodization strategy the paper transfers to the direct problem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Szegő recurrence, Verblunsky coefficients, and moment formulas used to recover the measure of a step-function Hamiltonian from |φ_n(1)|^2."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the inversion formula for finite Toeplitz matrices that underpins the algebraic proof of the moment-recovery algorithm."},{"cited_title":"de Branges","cited_arxiv_id":null,"evidence_quote":"Develops the de Branges space theory underlying the notion of spectral measures of canonical systems."}],"review_version":1}