{"id":"4648fc2e-76da-4534-8ee9-3119a1373871","arxiv_id":"2504.16424","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Singular values of complex symmetric tridiagonal Hamiltonians are computed as eigenvalues of a Hermitian block-tridiagonal partner via matrix continued fractions, with a fixed-point convergence analysis.","lead":"The paper presents a method to compute the real singular values of non-Hermitian tridiagonal quantum Hamiltonians by turning them into eigenvalues of a Hermitian block-tridiagonal matrix and evaluating the resolvent with matrix continued fractions. It illustrates convergence on a simple constant-coefficient example and gives a partial fixed-point analysis.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The MCF convergence claim is demonstrated only for constant-coefficient iterations, not for the variable-coefficient Hamiltonians (7), so the infinite-dimensional part of Lemma 1 remains unsupported.","rationale":"The reader's weakest-assumption identification is correct and aligns with the main structural gap. The paper's central mathematical idea is plausible and the finite-dimensional block-tridiagonal descendant construction is valid. What is missing is a rigorous limiting argument connecting the constant-coefficient fixed-point analysis of Section 4.2 to the actual variable-coefficient recurrence (16) for infinite N. The illustrative example only iterates the asymptotic scalar maps, not the matrix continued fraction for a genuine discretized Hamiltonian with, say, growing β_k. The paper's own β=2 example is not a contradiction but a demonstration that no real fixed point exists in that parameter region, which reinforces the need for a criterion rather than a proof. The concrete test proposed here would directly probe the method on a variable-coefficient Hamiltonian and compare MCF roots against exact finite-N singular values, thereby settling whether the fixed-point analysis supports the central claim. Since the reader already assigned a conditional verdict based on the same concern, no verdict adjustment is needed.","tokens_in":10505,"tokens_out":6148,"duration_ms":67922,"concrete_test":"Take α_k=1, γ_k=1/2, β_k=k and form finite truncations H_N of (7) for N=20, 40, 80, 160; compute their singular values by dense SVD. Run recurrence (16) with F_{N+1}=0 to obtain F_1(σ), and locate the values of σ where det F_1^{-1}(σ)=0. If, as N grows, these roots converge to the SVD singular values, the variable-coefficient MCF claim is confirmed. If the roots fail to converge, or converge to values different from the SVD results, even in a region where the asymptotic scalar map (20) with these β_k, γ_k has a real attracting fixed point, then the fixed-point analysis is insufficient and the infinite-dimensional reading of Lemma 1 is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The finite-dimensional construction in Section 4 is sound: for fixed N, the block-tridiagonal matrix (13) is the Pushnitski-Stampach descendant of (7), and the recurrence (16) terminates at F_{N+1}=0. The load-bearing gap is the N=∞ limit. Section 4.2 replaces the k-dependent recurrence by the constant-coefficient scalar maps (20) and studies their real fixed points, but this establishes neither necessity nor sufficiency for convergence of the original MCF. No contraction argument or Jacobian stability check is supplied; the stability of a candidate fixed point is only asserted. Moreover, the fixed-point analysis treats a forward iteration of a single map, while the actual recurrence runs backwards over changing β_k and γ_k, and the error from setting F_{M+1}=0 is not controlled. The numerical illustration in Table 1 iterates constant maps (20), not the MCF for a discretized Hamiltonian (7) with growing β_k, so it does not demonstrate the claimed quick convergence for the physical application. The paper itself flags the β=2, γ=1/2 choice as producing useless oscillatory results and acknowledges discrete-spectrum and boundedness assumptions, but no general criterion separates convergent from divergent regimes.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies complex symmetric tridiagonal Hamiltonians H of the form (7), interpreted as discretizations of one-dimensional Schrödinger operators with complex local potentials. It proposes to compute the singular values σ_n of H as eigenvalues of an auxiliary Hermitian block-tridiagonal operator H defined in (13)–(14), and to evaluate the resolvent of H by a matrix continued fraction (MCF) recurrence (16). The main claims are Lemma 1 (spectral equivalence for N≤∞), a fixed-point-based convergence analysis in §4.2 culminating in Lemma 2, and a numerical illustration in Table 1 purporting to show quick convergence. The finite-dimensional block construction is a reordering of the Hermitian dilation [[0,H],[H†,0]] and is correct; the burden of the paper lies in the infinite-dimensional and convergence claims.","tokens_in":10753,"tokens_out":20229,"duration_ms":194630,"significance":"If established, the method would give a practical route from complex resonance spectra to real singular values through a Hermitian block-tridiagonal eigenproblem, and the explicit 2×2 parametrization (17)–(19) is a genuinely useful computational reduction. Credit should be given for the clean finite-dimensional construction and for the clear scalar continued-fraction background in §3.3. However, the advertised contributions—the N=∞ spectral equivalence and the fixed-point proof of MCF convergence—are not established as written, and the numerical illustration exercises a simplified constant-coefficient map rather than the MCF for the variable-coefficient Hamiltonian (7).","major_comments":[{"comment":"The statement that Lemma 1 holds at any finite or infinite Hilbert-space dimension N≤∞ is not supported by the proof. The proof defines singular values as eigenvalues of the dilation (15) and then cites the Pushnitski–Stampach renumbering, but for N=∞ the equivalence requires operator-theoretic hypotheses on H (for example, a precise definition of singular values for non-compact operators and a self-adjoint realization of the dilation). The manuscript explicitly steps back from such hypotheses ('we will proceed in a more pragmatic manner') and merely assumes a discrete non-degenerate spectrum. Thus the infinite-dimensional part of Lemma 1 remains an assertion, and the abstract's claim about the discretized Schrödinger operator (7) is not justified. If the intended scope is finite N or a rigorously controlled truncation, that restriction should be stated and the claims adjusted accordingly.","section":"4, Lemma 1"},{"comment":"The convergence analysis replaces the k-dependent MCF recurrence (16) by the constant-coefficient scalar maps (20) and studies their real fixed points, but it does not prove convergence of the original backward recurrence for the variable-coefficient Hamiltonian (7). No relation between the asymptotic constant map and the actual k-dependent iteration is established, and no truncation-error bound for the boundary condition F_{M+1}=0 is supplied. Moreover, the stability of the selected fixed point is asserted rather than proved: unlike the scalar case in (11)–(12), no Jacobian or spectral-radius computation is given for the three-dimensional map (20). The β=2, γ=1/2 example is informative, but it does not provide a general criterion separating convergent from divergent regimes, and the paper itself acknowledges this limitation in §4.3.","section":"4.2, MCF convergence"},{"comment":"Table 1 iterates the constant-coefficient maps (20) with α=σ=1, β=4, γ=1/2; it does not compute the MCF (16) for a discretized Hamiltonian (7) with k-dependent diagonal elements β_k, γ_k. The caption should state this explicitly, because as presented the table invites the reading that the MCF for the physical application converges quickly. Consequently the abstract's statement that 'numerical MCF convergence is found quick' is not demonstrated by the numerical evidence in the paper; the example shows only convergence of the simplified constant-map iteration.","section":"4.3, Table 1"}],"minor_comments":[{"comment":"Lemma 2's proof is described only as 'standard elimination ... Gröbner basis'; since the quartic (23) and the relation (24) are load-bearing algebraic inputs, the derivation should be supplied in an appendix or made reproducible.","section":"4.2, Lemma 2"},{"comment":"The definition of singular values is cited to a Wikipedia article; for the operator-theoretic claims in Lemma 1, a mathematical reference with precise hypotheses on singular values of non-compact operators should be used.","section":"4, preliminaries"},{"comment":"The notation in (17) would be clearer if it explicitly stated that u_k, v_k, x_k, y_k are the parameters of F_k^{-1}, not of F_k; the subsequent map (18)–(19) is easier to follow when this convention is highlighted.","section":"4.1, Eq. (17)"}],"recommendation":"major_revision","confidential_remarks":"To the editor: I found no circularity in the central derivation; the recurrence and the finite-dimensional block construction are assembled from standard material. The main risk is that the paper's central convergence claim is supported only by a heuristic constant-coefficient analysis and by a toy numerical iteration, while the N=∞ part of Lemma 1 is asserted rather than proved. A revised version that restricts the rigorous statements to finite N or to a controlled truncation, and that clearly separates the heuristic convergence analysis from the proven results, would substantially strengthen the paper. The finite-dimensional part appears sound and worth publishing after such revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: the finite-dimensional identity in this paper is real and worth borrowing — singular values of a complex symmetric tridiagonal matrix come out of a Hermitian block-tridiagonal matrix via a terminating matrix continued fraction. The advertised N = ∞ convergence, however, rests on a constant-coefficient placeholder, not on the actual variable-coefficient recurrence. Don't rely on the convergence claim without checking it.\n\nWhat's genuinely new: the explicit MCF recurrence (16), the reduction of the 2x2 block map to the scalar maps (20), and the quartic fixed-point equation (23). For fixed finite N, the construction is sound: the recurrence terminates at F_{N+1}=0, and the spectral equivalence with the block matrix [[0,H],[H†,0]] after renumbering is exactly the Pushnitski–Stampach idea. That's a clear, legitimate extension.\n\nThe soft spots are in the infinite-dimensional part. Lemma 1 states N ≤ ∞, but the proof only handles finite matrices. For N = ∞ the author appeals to a pragmatic assumption of discrete non-degenerate spectrum, which is a framing, not a proof. The convergence section replaces the k-dependent β_k and γ_k by constants and studies fixed points of the asymptotic map. That tells you about the constant map, not why the original recurrence converges. There is no error bound on truncating at F_{M+1}=0 and no stability calculation — the stable fixed point in the example is simply asserted, and Table 1 iterates the constant map rather than an MCF for a discretized Hamiltonian (7) with growing diagonal. The β=2, γ=1/2 example shows a failure but gives no general criterion. These are real gaps.\n\nWho should read it: people computing resonances in tridiagonal models who want a Hermitian reformulation; the finite-N machinery could be useful. It deserves a serious referee, but the N=∞ claims need heavy revision — either restrict Lemma 1 to finite N or provide a genuine contraction argument. I'd bring it to a reading group with a caveat.","headline":"Useful finite-dimensional singular-value trick; the N=∞ convergence claim is asserted, not demonstrated.","tokens_in":11236,"tokens_out":3512,"would_cite":false,"duration_ms":34822,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["15A18","47B36","40A15","81Q12"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that singular values of a complex symmetric tridiagonal Hamiltonian can be computed as eigenvalues of a Hermitian block-tridiagonal partner via a matrix continued fraction, whose convergence is governed by real fixed…","keywords":["non-Hermitian Hamiltonians","complex symmetric matrices","tridiagonal matrices","singular values","matrix continued fractions","block-tridiagonal Hermitian operator","fixed-point convergence","quantum resonances"],"falsifier":"Run the iteration (20) with a parameter triple for which the quartic polynomial $P(u)$ has no real root—for example $\\sigma=1$, $\\gamma=1/2$, $\\beta=2$—and check whether the sequences $u,x,y$ settle to a fixed point. If they nevertheless converge to the correct singular values, the fixed-point criterion would not be necessary; if they diverge, the stated limitation is confirmed. A broader settlement would compare the singular values from the MCF computation with a direct full SVD of finite truncations of $H$.","tokens_in":10299,"feed_emoji":"⚛️","tokens_out":8086,"duration_ms":70526,"temperature":0.7,"pith_summary":"The paper is about unstable quantum states described by non-Hermitian tridiagonal Hamiltonians whose energy eigenvalues are complex. It proposes replacing the difficult task of finding those complex energies by the easier task of finding the real, non-negative singular values of the Hamiltonian. Its central claim is Lemma 1: for a complex symmetric tridiagonal Hamiltonian matrix, the singular values are exactly the eigenvalues of an auxiliary Hermitian block-tridiagonal matrix, so they can be obtained from a matrix continued fraction expansion of the auxiliary resolvent. A fixed-point analysis of the recurrence then tells when the continued fraction converges, and an illustrative discrete Schrödinger equation with a complex local potential shows fast convergence in a diagonal-dominated case and failure when the fixed points become complex. A sympathetic reader would care because the method gives real spectral data for open quantum systems without solving a complex eigenvalue problem.","feed_headline":"A Hermitian partner computes singular values of complex Hamiltonians","feed_subtitle":"Matrix continued fractions make the singular values computable, with a fixed-point rule deciding when convergence holds.","key_machinery":"The central object is the complex symmetric tridiagonal matrix Hamiltonian $H$ of Eq. (7), together with its Hermitian partner $\\mathbb{H}$ of Eq. (13), a block-tridiagonal matrix with $2\\times2$ blocks. The argument runs on the matrix continued fraction recurrence (16), obtained by replacing scalar entries in the standard analytic continued fraction factorization of $H-E$ by Hermitian two-by-two blocks, and then on its reduction to the three coupled scalar iteration maps (20) with real parameters $u,x,y$. The convergence of the whole scheme is controlled by fixed points of these maps, and Lemma 2 packages the fixed-point condition as the quartic equation $P(u)=0$; this real-root condition is what makes the method either converge quickly or fail, as the two illustrative examples demonstrate.","core_discovery":"For every complex symmetric tridiagonal Hamiltonian $H$ of the form (7), with real $\\alpha_k, \\beta_k, \\gamma_k$ and $N \\leq \\infty$, the singular values $\\sigma_n$ of $H$ coincide with the eigenvalues of the Hermitian block-tridiagonal matrix $\\mathbb{H}$ in (13), whose $2\\times2$ blocks are $A_k = \\begin{pmatrix} 0 & \\beta_k+i\\gamma_k \\\\ \\beta_k-i\\gamma_k & 0 \\end{pmatrix}$ and $B_k = \\begin{pmatrix} 0 & \\alpha_k \\\\ \\alpha_k & 0 \\end{pmatrix}$. The proof goes through the doubled matrix $\\widetilde{H} = \\begin{pmatrix} 0 & H \\\\ H^\\dagger & 0 \\end{pmatrix}$ and a renumbering of basis elements. Since $\\mathbb{H}$ is Hermitian, its resolvent factorizes into a matrix continued fraction through the recurrence (16) for two-by-two matrices $F_k$, and the paper reduces this recurrence to the coupled scalar maps (20) for $u,x,y$. Lemma 2 gives a quartic polynomial $P(u)$ whose real roots are the candidate fixed points of these maps, so convergence of the continued fraction is decided by the existence of a real attracting fixed point. The paper's examples show quick numerical convergence for one parameter set and non-convergence for a nearby set.","pith_inferences":["If the fixed-point condition is sharpened, the transition between real and complex roots of $P(u)$ could provide a boundary in parameter space separating computable from non-computable regimes; mapping this boundary numerically would be a direct test of the method's domain.","The block-tridiagonal partner has the structure of a tight-binding model on a doubled chain, so the convergence of the continued fraction may correspond to a localization transition of the associated transfer-matrix map, connecting the result to transport phenomenology.","A natural extension is to $M\\times M$ blocks: the quartic fixed-point condition would become a higher-degree algebraic condition, and one could test whether real-root obstructions of the same kind control convergence there.","Since the iteration maps (20) are low-dimensional, the derivative at the stable fixed point could yield a quantitative bound on the contraction rate, giving an a priori stopping criterion for practical computations."],"forward_implications":["Singular values of complex-symmetric tridiagonal resonance Hamiltonians can be obtained as real eigenvalues of a Hermitian block-tridiagonal operator, bypassing direct complex spectral computations.","The matrix continued fraction recurrence gives a practical numerical scheme whose convergence rate is governed by the fixed-point map (20), with the quartic polynomial (23) as a closed-form diagnostic.","The method inherits a sharp boundary: when the fixed-point equation has no real root, as for $\\sigma=1$, $\\gamma=1/2$, $\\beta=2$, the continued fraction fails, so not every complex local potential can be treated in this way.","The Hermitian partner $\\mathbb{H}$ has all resolvent poles on the real half-axis, so the singular values are directly readable from the secular equation $\\det F_1^{-1}(\\sigma)=0$ in the MCF formalism."],"supporting_citations":[{"why":"Supplies the spectral equivalence between singular values of $H$ and eigenvalues of the doubled Hermitian operator, the basis of Lemma 1.","marker":"[11]"},{"why":"Introduces the matrix continued fraction solution for anharmonic-oscillator eigenvalues, the origin of the MCF recurrence (16).","marker":"[21]"},{"why":"Shows the analytic continued fraction resolvent of the anharmonic oscillator with poles at energy eigenvalues, the scalar template the MCF generalizes.","marker":"[17]"},{"why":"Provides the fixed-point perturbation approach to continued fraction convergence that section 3.3 transplants to the matrix case.","marker":"[19]"},{"why":"Reviews complex symmetric operators and supports the relevance of the class of matrices in Eq. (7).","marker":"[12]"},{"why":"Foundational treatment of complex symmetric operators, underlying the transposition-symmetric Hamiltonian ansatz.","marker":"[13]"},{"why":"Supplies a family of anharmonic oscillators where matrix continued fractions converge quickly, used as a comparative success case.","marker":"[28]"}],"fun_headline_variants":["Hermitian partner reveals singular values of complex Hamiltonians","Matrix continued fractions compute singular values with fixed-point proof","Complex Hamiltonians meet Hermitian trick for singular values","Singular values via Hermitian block operator and continued fractions","Fixed-point rule guarantees convergence in matrix continued fractions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof of convergence assumes that for large indices the recurrences can be replaced by constant-parameter maps with real parameters that possess a real attracting fixed point, and this fails for some complex potentials (the $\\beta=2$, $\\gamma=1/2$ example), so the method is not guaranteed to converge in general.","fun_headline_variants_meta":{"raw":{"variants":["Hermitian partner reveals singular values of complex Hamiltonians","Matrix continued fractions compute singular values with fixed-point proof","Complex Hamiltonians meet Hermitian trick for singular values","Singular values via Hermitian block operator and continued fractions","Fixed-point rule guarantees convergence in matrix continued fractions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000265,"raw_usage":{"total_tokens":1620,"prompt_tokens":972,"completion_tokens":648,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":572}},"tokens_in":588,"tokens_out":648,"duration_ms":6869,"temperature":1.0,"reasoning_tokens":572,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:04:03.751654+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the iteration (20) with a parameter triple for which the quartic polynomial $P(u)$ has no real root—for example $\\sigma=1$, $\\gamma=1/2$, $\\beta=2$—and check whether the sequences $u,x,y$ settle to a fixed point. If they nevertheless converge to the correct singular values, the fixed-point criterion would not be necessary; if they diverge, the stated limitation is confirmed. A broader settlement would compare the singular values from the MCF computation with a direct full SVD of finite truncations of $H$.","supporting_citations":[{"cited_title":"Singh, S","cited_arxiv_id":null,"evidence_quote":"Shows the analytic continued fraction resolvent of the anharmonic oscillator with poles at energy eigenvalues, the scalar template the MCF generalizes."},{"cited_title":"Znojil, Fixed point perturbation theory and the potential r2 + λr 2/ (1 + gr 2)","cited_arxiv_id":null,"evidence_quote":"Provides the fixed-point perturbation approach to continued fraction convergence that section 3.3 transplants to the matrix case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reviews complex symmetric operators and supports the relevance of the class of matrices in Eq. (7)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundational treatment of complex symmetric operators, underlying the transposition-symmetric Hamiltonian ansatz."},{"cited_title":"Znojil, Symmetrically anharmonic oscillators","cited_arxiv_id":null,"evidence_quote":"Supplies a family of anharmonic oscillators where matrix continued fractions converge quickly, used as a comparative success case."}],"review_version":1}