{"id":"d9640e8f-ca10-4996-9473-0b002fea3b91","arxiv_id":"2505.07570","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The authors show that the truncated Hamburger, Stieltjes, and Hausdorff moment problems reduce to generalized eigenvalue problems for Hankel matrices, and they reformulate classical determinacy criteria in this language.","lead":"This math paper connects classical moment problems, where you reconstruct a probability-like measure from its averages, to the boundary control method for discrete dynamical systems. It shows that solving a truncated moment problem is equivalent to solving a generalized eigenvalue problem built directly from Hankel matrices of the moments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Section 3 norming-constant step inverts ρ_k: W_N f_k = α_k φ_k and (C_N f_k, f_k) = 1 force ρ_k = α_k^{-2}, not α_k^2; as written the recovered measure has reciprocal masses.","rationale":"The reader's conditional verdict is appropriate, but the weakest assumption flagged in the reader's report — Remark 4's delegation of A_N recovery to prior work — is not the most load-bearing point. Section 3 explicitly bypasses the recovery of A_N: 'we recover the spectral measure ... without recovering the Jacobi matrix itself.' The load-bearing step is the claimed recovery of the norming constants from the generalized eigenvectors in the final procedure. There, a direct algebraic inconsistency exists: the normalization (C_N f_k, f_k) = 1 combined with W_N f_k = α_k φ_k forces α_k^2 ρ_k = 1, hence ρ_k = α_k^{-2}; the paper states ρ_k = α_k^2. Since (2.3) defines the measure as Σ 1/ρ_k δ_{λ_k}, the reconstructed measure has reciprocal masses. The N = 2 free Jacobi example makes this concrete: the paper's rule assigns masses 2, while the prescribed moments require masses 1/2. The error is localized and likely a typo — the inverse is forced by the surrounding equations — but as written, the central claim that the truncated moment problem is equivalent to the generalized spectral problem is not correct. A corrected formula ρ_k = α_k^{-2} would restore the argument. Hence the verdict remains CONDITIONAL: the construction is probably sound after a required correction, but the current text contains a concrete numerical error in the main procedure. This is a different concern from the reader's; I do not find the delegation in Remark 4 to be the primary risk, because the main path does not require recovering A_N.","tokens_in":12844,"tokens_out":41994,"duration_ms":321546,"concrete_test":"Implement the procedure for N = 2 with the free Jacobi matrix (a_1 = 1, b_1 = b_2 = 0). Given moments s_0 = 1, s_1 = 0, s_2 = 1, s_3 = 0, form S^0_2 = [[s_2, s_1], [s_1, s_0]] = I and S^1_2 = [[s_3, s_2], [s_2, s_1]] = [[0, 1], [1, 0]]. Solve S^1_2 g = λ S^0_2 g; for λ = 1 take g = (1, 1). Since \\tildeΛ_2 = I, f = g. Normalize via (C_2 f, f) = 1 with C_2 = S^0_2 = I, giving f = (1/√2, 1/√2). Compute α = (R f)_2 = f_1 = 1/√2. Reconstruct the measure via (2.3) once with ρ = α^2 = 1/2 and once with ρ = α^{-2} = 2. Check ∫ λ^k dρ for k = 0, 1, 2 against (1, 0, 1). Only ρ = 2 matches; ρ = 1/2 gives (4, 0, 4). This settles whether the norming-constant formula is inverted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The final procedure in Section 3 asserts that after solving the generalized spectral problem (3.10), normalizing controls by (C_N f_k, f_k) = 1 and setting α_k = (R f_k)_N, the norming constants are ρ_k = α_k^2. This is inconsistent with the paper's own definitions. From (2.12)–(2.13), (C_N f, g)_{F^N} = (W_N f, W_N g)_{H^N}. The normalization therefore gives 1 = (C_N f_k, f_k) = ||W_N f_k||^2. Step 2 states W_N f_k = α_k φ_k, where (φ_k, φ_l) = δ_{kl} ρ_k. Hence 1 = α_k^2 ρ_k, so ρ_k = α_k^{-2}, not α_k^2. Since the spectral function (2.3) is ρ_N = Σ_{λ_k < λ} 1/ρ_k, the paper's formula assigns masses α_k^{-2} instead of α_k^2 — the reciprocal masses. This is not a harmless scaling. For the free Jacobi matrix A_2 (a_1 = 1, b_1 = b_2 = 0), the moments are s_0 = 1, s_1 = 0, s_2 = 1, and the pencil (3.10) gives α = 1/√2 for each eigenvalue. The paper's rule yields ρ = 1/2 and hence the measure 2δ_1 + 2δ_{-1}, whose moments are s_0 = 4, s_2 = 4, not (1, 0, 1). The correct constant ρ = 2 gives masses 1/2, matching the prescribed moments. Thus the central claim that the truncated moment problem is equivalent to the generalized spectral problem fails as written; the formula must be ρ_k = α_k^{-2}.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies classical moment problems (Hamburger, Stieltjes, Hausdorff) via the boundary control method for discrete-time dynamical systems associated with Jacobi matrices. The authors introduce de Branges spaces of polynomials, express the reproducing kernel and Christoffel functions through the connecting operator, and reduce the truncated moment problem to the generalized spectral problem S^N_1 g = λ S^N_0 g for Hankel matrices of moments. They then propose a procedure to recover the spectral measure of the finite Jacobi matrix from the eigenvectors of this pencil, and use the spectral formulation to re-derive conditions for determinacy and indeterminacy of the full moment problems, including a new proof of the simultaneous determinacy of Hamburger and Stieltjes problems.","tokens_in":13293,"tokens_out":17812,"duration_ms":143429,"significance":"If the recovery procedure is corrected, the reduction of the truncated moment problem to a single Hankel pencil eigenvalue problem is a clean and potentially useful computational formulation. The paper also gives a dynamic-system interpretation of de Branges spaces and reproduces classical determinacy criteria as ratios of determinants of Hankel matrices, connecting them to the boundary control method. The main mathematical chain is consistent with the classical theory of Jacobi matrices and orthogonal polynomials. However, the norming-constant step in the recovery procedure contains an algebraic error that must be fixed before the central claim of the paper is valid. The paper relies substantially on prior work by the same authors, but that reliance is acknowledged; the main new idea is the direct reduction to Hankel pencils.","major_comments":[{"comment":"Steps 1-3 of the norming procedure are inconsistent with the definitions in Section 2. Since (C^N f,g)_{F^N} = (W^N f, W^N g)_{H^N} by (2.12)-(2.13), the normalization (C^N f_k, f_k)=1 gives ||W^N f_k||^2=1. Step 2 states W^N f_k = α_k φ_k, where by Definition 2 we have (φ_k,φ_k)=ρ_k. Hence 1 = α_k^2 ρ_k, so ρ_k = α_k^{-2}, not α_k^2 as written in Step 3. Because the spectral function (2.3) has jumps 1/ρ_k, the paper's formula assigns masses α_k^{-2} instead of α_k^2. For the free Jacobi matrix A_2 with a_1=1, b_1=b_2=0, the moments are s_0=1, s_1=0, s_2=1, and the pencil (3.10) gives α_k^2=1/2 for both eigenvalues; the paper's rule yields ρ_k=1/2 and hence the measure 2δ_1+2δ_{-1}, whose moments are s_0=4, s_2=4, not the prescribed (1,0,1). The correct constant ρ_k=2 gives the desired measure (1/2)δ_1+(1/2)δ_{-1}. The formula in Step 3 must be changed to ρ_k = α_k^{-2}.","section":"Section 3, procedure after Theorem 4"},{"comment":"Theorem 4 is a load-bearing assertion: it states that the dynamic matrices C^N and B^N are congruent via \\tilde Λ^N to the Hankel matrices S^N_0 and S^N_1, and this is precisely what converts the dynamic pencil (3.9) into the Hankel pencil (3.10). No proof or derivation is given in the manuscript. This step is essential to the central claim of the paper and should be proved in the text; it is likely an elementary computation from (2.7), (2.8), and (3.1)-(3.3), but the authors need to show it rather than leave it as an unproved statement.","section":"Section 3, Theorem 4"},{"comment":"The identification W^N f_k = α_k φ_k for the eigenvectors obtained from (3.9) is asserted without justification, and the constant α_k is defined only through the undefined symbol (R f_k)_N. This identification is what links the algebraic spectral problem to the spectral data of the Jacobi matrix A_N. The authors should either prove this fact directly or state precisely which result from [11,13] guarantees it, and give the definition of R appearing in α_k = (R f_k)_N.","section":"Section 3, Step 2 and Remark 4"}],"minor_comments":[{"comment":"The symbol R in α_k = (R f_k)_N is not defined; only R^T and R^T_N are introduced in Section 2. Please define R explicitly or replace this expression with a concrete formula.","section":"Section 3, Step 2"},{"comment":"The matrices H^{2T-1}_2 and H^{2T}_0 are used without definition. The tilded matrices \\tilde H^{2T}_0 and \\tilde H^{2T}_1 are defined in (4.11), but the plain H and the subscript 2 are never introduced; please state the definitions explicitly.","section":"Section 4, Proposition 7 and its proof"},{"comment":"The matrix S^{T+1}_{0,0} appears without definition. If it is the (T+1)-dimensional analogue of S^T_{0,0} defined after (4.9), please say so and give its explicit form.","section":"Section 4, Proposition 8 and proof of Proposition 9"},{"comment":"The bilinear form (1.2) and the scalar products (2.12)-(2.13) are written without complex conjugation. For complex polynomials this is not a positive-definite scalar product; please specify that the polynomials are real, or insert complex conjugates in the appropriate places.","section":"Section 2, equations (1.2), (2.12)-(2.13)"},{"comment":"There are numerous typos and slips, including 'Jacoi' in Section 3, 'Stiltjes' in Section 4, 'stetement' before Proposition 9, 'thethe' in the text, and 'week' for 'weak' in the Conclusion. A careful proofreading pass is needed.","section":"Throughout"},{"comment":"Theorems 3 and 5 are quoted from [14], which is an arXiv preprint. Please clarify whether [14] has been published in its current form and, if so, update the reference to the published version.","section":"Section 3, Theorems 3 and 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is in scope for math.SP and the underlying idea is appealing, but the norming-constant error in the Section 3 procedure is a serious flaw that directly breaks the central claim. The fix is straightforward (invert the formula), and after that the main reduction may be salvageable. The revised version should also provide a proof of Theorem 4, which is currently asserted without derivation, and clean up the undefined notation in Section 4."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis one is a reasonable reformulation paper with a concrete bug in the reconstruction step. The de Branges space perspective on the BC method is a nice unifying frame, and the reduction of the connecting operator and the pencil to Hankel matrices is clean. That part works. But the norming constants in Section 3 are inverted: with W_N f_k = α_k φ_k and (C_N f_k, f_k) = 1, the definitions give ρ_k = α_k^{-2}, not α_k^2. The simple N=2 example (A_2 with a_1=1, b_1=b_2=0) shows the paper's formula produces moments (4,0,4) instead of (1,0,1). So the central claim that the truncated problem reduces to the generalized spectral problem fails as written. It's a fixable slip—replace α_k^2 with α_k^{-2}—but it's load-bearing.\n\nThe paper does some things well. Theorem 4 is a useful algebraic identity, though it is stated without proof. The reformulations of the classical determinacy criteria in Hankel determinant language (Propositions 7-9) are correct and provide a convenient dictionary. The authors are honest about relying on their earlier recovery results, though that makes the self-citation chain heavy.\n\nSoft spots beyond the norming error: Section 4 has notational sloppiness (H^{2T-1}_2 and S^{T+1}_{0,0} are never defined), and several steps are dismissed as 'easy' when they carry actual work. None of those are fatal once the norming constant is fixed.\n\nWho is this for? People already inside the BC-method/de Branges community who want a unified statement of the moment-problem reformulation. It won't change how someone outside the area does moment problems.\n\nIf this lands on my desk, I'd send it to a referee with a note to check the norming constants carefully. The idea deserves a serious look, but the current version has a central error that must be corrected before it can be accepted.","headline":"A useful de Branges reformulation of the BC method for moment problems, but the Section 3 norming constants are inverted, so the reconstruction procedure fails as written on a simple example.","tokens_in":13774,"tokens_out":3312,"would_cite":false,"duration_ms":28970,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["47B36","44A60","15A22","46E22"],"pacs":[],"model":"deepseek-v4-flash","headline":"Truncated Hamburger, Stieltjes, and Hausdorff moment problems are equivalent to generalized spectral problems for Hankel moment matrices, whose eigenpairs supply the spectral data needed to reconstruct the measure.","keywords":["moment problem","Jacobi matrix","boundary control method","generalized spectral problem","Hankel matrix","de Branges space","truncated moment problem","determinacy"],"falsifier":"Solve $S^N_1 g = \\lambda S^N_0 g$ numerically for a measure whose moment sequence and Jacobi coefficients are known, and compare the resulting eigenvalues with the zeros of the polynomial $\\phi_{N+1}(\\lambda)$ generated by that Jacobi matrix; a mismatch for some $N$ would refute the claimed equivalence.","tokens_in":12676,"feed_emoji":"🧮","tokens_out":12631,"duration_ms":110705,"temperature":0.7,"pith_summary":"This paper establishes that the classical moment problems—recovering a measure on the real line, the positive half-line, or the unit interval from prescribed moments—can be recast as finite generalized spectral problems. The central vehicle is a discrete-time dynamical system associated with a Jacobi matrix, whose boundary-control inverse data are determined by the moments. The authors show that the generalized eigenvalue pencil $S^N_1 g = \\lambda S^N_0 g$, built from two Hankel matrices of moments, carries exactly the spectral data of the truncated Jacobi matrix, so the measure can be reconstructed from moment data without first recovering the Jacobi coefficients. They also rewrite known determinacy criteria for Hamburger and Stieltjes problems as limits of ratios of Hankel determinants and give a new proof that a Stieltjes moment sequence is determinate exactly when its interlaced Hamburger extension is determinate.","feed_headline":"One matrix pencil solves truncated moment problems","feed_subtitle":"Two Hankel moment matrices carry the spectral data of the Jacobi matrix, so a measure is recovered directly","key_machinery":"The central object is the generalized spectral pencil $S^N_1 g = \\lambda S^N_0 g$, with $S^N_0$ and $S^N_1$ Hankel matrices of moments, together with the factorization identities $C^N = \\tilde\\Lambda_N S^N_0 \\tilde\\Lambda_N^*$ and $B^N = \\tilde\\Lambda_N S^N_1 \\tilde\\Lambda_N^*$. Here $\\tilde\\Lambda_N$ is the Chebyshev-polynomial change-of-basis matrix conjugated by the reversal matrix, so the same linear transformation converts the boundary-control operators of the dynamical system into Hankel moment matrices. The pencil replaces the earlier spectral problem $B^N f = \\lambda C^N f$ for the connecting and spectral operators; its eigenpairs give the Dirichlet spectral data of the Jacobi matrix directly from moments. The identification of the polynomial space with the moment inner product as a de Branges space $B^N_J$ is what ties the moment problem to the boundary-control inverse problem and provides the reproducing kernel and Christoffel symbols.","core_discovery":"On the paper's own terms, the discovery is that solving a truncated moment problem of order $N$ is equivalent to solving the generalized spectral problem $S^N_1 g_k = \\lambda_k S^N_0 g_k$, where $S^N_0$ and $S^N_1$ are Hankel matrices formed from the moments $s_0,\\dots,s_{2N-1}$. The eigenvectors $g_k$ are related to the boundary controls $f_k$ by $g_k = (\\tilde\\Lambda_N)^* f_k$; the controls drive the discrete dynamical system into the eigenfunctions of the $N\\times N$ Jacobi matrix, and the eigenvalues $\\lambda_k$ are the Dirichlet eigenvalues of that matrix. From these eigenpairs the norming constants $\\rho_k$ are obtained through the connecting operator, and the measure $d\\rho_N$ is assembled as the spectral sum in (2.3). The same machinery factorizes the connecting operator and the spectral operator as $C^N = \\tilde\\Lambda_N S^N_0 \\tilde\\Lambda_N^*$ and $B^N = \\tilde\\Lambda_N S^N_1 \\tilde\\Lambda_N^*$, and identifies the polynomial space with the moment inner product as a de Branges space.","pith_inferences":["The explicit Hankel structure suggests that truncated moment problems could be tackled with structured generalized eigenvalue algorithms that exploit Toeplitz/Hankel low-rank structure, a practical direction the paper does not develop.","The $\\pm\\sqrt{\\mu}$ correspondence between the Stieltjes and Hamburger pencils hints at a broader duality: determinacy of an even Hamburger measure with vanishing odd moments may be equivalent to determinacy of its compressed Stieltjes moment sequence beyond the specific interlacing used in Proposition 9.","The determinant formula for the reproducing kernel in Remark 7 could be used to compute orthogonal-polynomial recurrence coefficients numerically from moment matrices alone, rather than from the Jacobi matrix."],"forward_implications":["Truncated moment problems can be solved as a single generalized eigenvalue computation on moment Hankel matrices, without a separate reconstruction of the Jacobi coefficients.","The finite spectral data of the truncated Jacobi matrix—eigenvalues and norming constants—are exactly the eigenpairs of the Hankel pencil, so any generalized eigenvalue solver yields the approximate measure.","The de Branges-space identification expresses reproducing kernels and Christoffel symbols in terms of moments, connecting orthogonal-polynomial quantities to boundary-control data.","The determinacy criteria for the Hamburger and Stieltjes problems are rewritten as limits of ratios of Hankel determinants, checkable directly from the moment sequence.","The interlacing construction shows that the Stieltjes problem and the associated Hamburger problem are determinate together, with eigenvalues paired as $\\pm\\sqrt{\\mu}$ between the two pencils."],"supporting_citations":[{"why":"Supplies the recovery of the Jacobi matrix from the response vector and the positivity criterion for the connecting operator, which underpins the claim that moment data determine the truncated Jacobi matrix.","marker":"[13]"},{"why":"Provides the de Branges-space construction for the discrete dynamical system and the reproducing kernel formulas used to identify the moment space.","marker":"[12]"},{"why":"Introduces the original relationship between moment problems and this boundary-control system, including the earlier spectral pencil that this paper rewrites in Hankel form.","marker":"[14]"},{"why":"Establishes the earlier dynamic inverse problem results, including the representation of response vectors, that the paper relies on for recovery from finite data.","marker":"[11]"},{"why":"Supplies the classical theory of moment problems as spectral problems of self-adjoint Jacobi operators, including the index-and-determinacy connection.","marker":"[18]"},{"why":"Provides the modern moment-problem apparatus of polynomial spaces, Christoffel functions, and determinacy criteria used in the reformulations.","marker":"[17]"},{"why":"Supplies the classical baseline for moment problems and orthogonal polynomials, including the construction of the measure from spectral data.","marker":"[1]"}],"fun_headline_variants":["One matrix pencil recovers the measure","Truncated moments as spectral problem","Boundary control cracks moment problems","Hankel pencil yields the measure","Moment problems via boundary control"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a finite list of moments up to order $2N-2$ determines the $N\\times N$ Jacobi matrix (equivalently, the response vector of the dynamical system), a recovery result the paper borrows from earlier work rather than proving here.","fun_headline_variants_meta":{"raw":{"variants":["One matrix pencil recovers the measure","Truncated moments as spectral problem","Boundary control cracks moment problems","Hankel pencil yields the measure","Moment problems via boundary control"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000553,"raw_usage":{"total_tokens":2591,"prompt_tokens":854,"completion_tokens":1737,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":470,"completion_tokens_details":{"reasoning_tokens":1678}},"tokens_in":470,"tokens_out":1737,"duration_ms":13306,"temperature":1.0,"reasoning_tokens":1678,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:13:43.990765+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve $S^N_1 g = \\lambda S^N_0 g$ numerically for a measure whose moment sequence and Jacobi coefficients are known, and compare the resulting eigenvalues with the zeros of the polynomial $\\phi_{N+1}(\\lambda)$ generated by that Jacobi matrix; a mismatch for some $N$ would refute the claimed equivalence.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the recovery of the Jacobi matrix from the response vector and the positivity criterion for the connecting operator, which underpins the claim that moment data determine the truncated Jacobi matrix."},{"cited_title":"Schr¨ odinger operator, Dirac system, discrete Schr¨ odinger operatorJournal of Mathemati- cal Analysis and Applications460927–53","cited_arxiv_id":null,"evidence_quote":"Provides the de Branges-space construction for the discrete dynamical system and the reproducing kernel formulas used to identify the moment space."},{"cited_title":"Dynamic inverse problem for special system associated with Jacobi matrices and classical moment problems","cited_arxiv_id":"1907.11153","evidence_quote":"Introduces the original relationship between moment problems and this boundary-control system, including the earlier spectral pencil that this paper rewrites in Hankel form."},{"cited_title":"2016Dynamical inverse problem for the discrete Schr¨ odinger operatorNanosystems: Physics, Chemistry, Mathematics7842–54","cited_arxiv_id":null,"evidence_quote":"Establishes the earlier dynamic inverse problem results, including the representation of response vectors, that the paper relies on for recovery from finite data."},{"cited_title":"1998The classical moment problem as a self-adjoint finite difference operator Advances in Math.13782–203 1St","cited_arxiv_id":null,"evidence_quote":"Supplies the classical theory of moment problems as spectral problems of self-adjoint Jacobi operators, including the index-and-determinacy connection."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the modern moment-problem apparatus of polynomial spaces, Christoffel functions, and determinacy criteria used in the reformulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classical baseline for moment problems and orthogonal polynomials, including the construction of the measure from spectral data."}],"review_version":1}