{"id":"e3b494cf-a945-474b-8b51-7f3445c6626b","arxiv_id":"1909.00963","paper_version":6,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A 4x4 Riemann-Hilbert formalism gives large-n asymptotics for Toeplitz+Hankel determinants with independent symbols, conditional on a non-degeneracy bound and verified for an explicit family.","lead":"This paper develops a new 4x4 matrix Riemann-Hilbert framework for studying sums of Toeplitz and Hankel determinants, and uses it to derive the large-size asymptotics for a solvable class of such determinants. The main result is a step toward spectral asymptotics of Hankel matrices, important in random matrix theory and statistical mechanics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 1.1 rests on non-degeneracy (1.9), verified only for an explicit family; if E(n) can vanish or decay too fast for other admissible symbols, formula (1.10) is not established.","rationale":"I read the paper in good faith and traced the main line of argument. The 4x4 RHP construction, the reduction to the model problem, the explicit solution for d\\tilde d=1, and the small-norm Riemann-Hilbert estimates form a coherent and carefully executed derivation. The proof of Theorem 1.1 is internally consistent: the non-degeneracy condition (1.9) is exactly what is needed to make the error term in (4.21) relatively small, and the final algebra leading to h_{n-1} = -α(0)E(n)/E(n-1) checks out, including the shift relating R1,jk(0;n) to R^{(1)}_{jk}(n+1). The explicit example (4.38)-(4.39) demonstrates that the class satisfying (1.9) is nonempty, with the coefficient κ claimed to be nonzero away from a small exceptional set. I find no fatal flaw or internal contradiction. The strongest caveat, already identified by the reader, is that condition (1.9) is not verified for the full class of symbols satisfying the other hypotheses; if E(n) has zeros or decays faster than any allowed r^n, the theorem is vacuous for those symbols. The paper itself acknowledges this in Section 5.5, listing the characterization of C0 as an open problem. There is also a mild overstatement in the abstract concerning the interval-supported Hankel symbol, since Section 3 only reduces that case to the same model problem without proving an asymptotic theorem; however, this does not affect the correctness of Theorem 1.1. The concern is about scope and genericity, not about the validity of the conditional theorem, so the reader's ACCEPT verdict remains appropriate.","tokens_in":42309,"tokens_out":20699,"duration_ms":197801,"concrete_test":"For a numerically generated family of admissible smooth symbols, e.g. d(e^{iθ}) = exp(i∑_{k=1}^{K}(c_k e^{ikθ} - \\bar c_k e^{-ikθ})) so that d\\tilde d=1 on T, and φ of Szegő type with inner singularities at known radius, compute E(n) in (1.6) by direct contour quadrature for n = 50, 100, ..., 500. Check whether |E(n)| remains bounded below by C r^n for some r∈[r0,1) and C>0. If E(n) is found to vanish or to decay faster than r^n for a statistically significant fraction of random admissible symbols, the non-degeneracy condition is nongeneric and Theorem 1.1 applies only to a restricted, currently uncharacterized subclass. A complementary check is to re-evaluate the Watson-lemma coefficient κ in (4.41) numerically on a dense parameter grid with α1 real and α1∉{±1/2}, confirming κ≠0 and that the explicit family indeed satisfies (1.9).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central asymptotic formula (1.10) is derived by reconstructing the Y-RHP from the X-RHP, which requires one of the six conditions of Lemma 2.7; the proof specifically relies on condition (2.36). The key estimate (4.19)-(4.21) shows that this condition holds only when |E(n)| is bounded below by Cr^n as in (1.9). Thus the theorem is logically sound but entirely conditional: if E(n) has zeros for arbitrarily large n, or decays faster than r^n for every admissible r in [r0,1), then the lower bound fails, step (4.21) breaks down, and the ratio formula is not obtained for those symbol pairs. The paper verifies (1.9) only for the explicit family (4.38)-(4.39) via the Watson-lemma estimate (4.40), with the coefficient κ asserted to be generically nonzero. For the general class of smooth symbols satisfying d\\tilde d=1, no verification is given, and Section 5.5 explicitly leaves the characterization of the subclass C0 (where E(n) obeys (1.11)) as an open problem. This is an acknowledged limitation of scope rather than an internal inconsistency, but it means the advertised general asymptotic result is not yet established beyond the explicit family and whatever further cases satisfy (1.9).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a 4x4 Riemann-Hilbert formalism for Toeplitz+Hankel determinants in which the Toeplitz and Hankel symbols are not assumed to be related. For the case where the Hankel symbol is supported on the unit circle, the authors carry out a steepest descent analysis and reduce the problem to a model Riemann-Hilbert problem. For the special class w = dφ with d satisfying d(z)d(z^{-1})=1, they construct the model solution explicitly and derive a conditional asymptotic formula for the ratio h_{n-1}=D_n/D_{n-1}, namely (1.10), under the non-degeneracy hypothesis (1.9). They also show that the same model problem emerges when the Hankel symbol is supported on an interval [a,b], and they provide a concrete one-parameter family (4.38)-(4.39) for which the non-degeneracy condition is verified. Section 5 lists open questions, including the characterization of the classes C and C0.","tokens_in":42522,"tokens_out":15586,"duration_ms":158519,"significance":"If the proof is completed, this is a substantial contribution to the asymptotic theory of Toeplitz+Hankel determinants. The 4x4 framework is new in allowing unrelated Toeplitz and Hankel symbols, it connects the unit-circle and interval-supported Hankel cases through the same model problem, and it yields large-n asymptotics for the associated orthogonal polynomials as well as for the norm parameter h_n. The paper is honest about the conditional nature of the main result: Theorem 1.1 explicitly depends on (1.9), and Section 5.5 acknowledges that the class of symbols satisfying this condition is not characterized. The worked example (4.38)-(4.41) providing a nonempty subclass C0 is a genuine strength, as is the comparison with the operator-theoretic results of Basor and Ehrhardt. No parameter fitting or circular reasoning is involved: the asymptotics are expressed through the explicit functional E(n) in (1.6). However, the proof of the final ratio formula contains an index discrepancy that must be repaired before the theorem is established as stated.","major_comments":[{"comment":"There is a load-bearing index error in the passage from (4.36) to (4.37). With the definitions as printed, (4.12) gives R_{1,jk}(0;n) = (1/2πi)∫_{Γ'_i} μ^n g_{jk}(μ) dμ for jk=12,14,23,43, while (4.31) defines R^{(1)}_{jk}(n) = (1/2πi)∫_{Γ'_i} μ^{n-2} g_{jk}(μ) dμ for the same index set. Hence the correct relation is R_{1,jk}(0;n) = R^{(1)}_{jk}(n+2), not R_{1,jk}(0;n) = R^{(1)}_{jk}(n+1) as stated in the text. Consequently the denominator in (4.36) equals E(n-2) under the displayed definitions, not E(n-1), and the identity (4.37), which is used to prove formula (1.10) of Theorem 1.1, does not follow. If the intended definition in (4.31) is μ^{n-1} on Γ'_i, the exponent should be corrected consistently; otherwise the derivation leading to (4.37) needs to be repaired. This issue also propagates to the explicit asymptotics stated after (4.41), which rely on the ratio E(n)/E(n-1).","section":"Section 4.2, Eqs. (4.31) and (4.37)"},{"comment":"The non-degeneracy hypothesis (1.9) is essential to the proof: it is used in (4.19)-(4.21) to guarantee condition (4.18) of Lemma 2.7, which in turn is needed to reconstruct the Y-RHP from the X-RHP. The hypothesis is verified only for the explicit family (4.38)-(4.39), and Section 5.5 explicitly leaves open the characterization of the classes C and C0. This is not an internal inconsistency, since Theorem 1.1 is honestly conditional, but the advertised scope of the asymptotic result is accordingly narrower than a first reading of the abstract suggests. I would ask the authors to state prominently that Theorem 1.1 applies subject to (1.9), and that beyond the explicit family the condition is verified only in particular cases; this is a scope clarification rather than a request for new theorems.","section":"Section 1, Theorem 1.1; Section 5.5"}],"minor_comments":[{"comment":"The sentence \"This fact is confirmed numerically\" is not accompanied by any numerical data, figure, or description of the computation; either provide the numerical verification or remove the sentence.","section":"Section 4.2, after Eq. (4.41)"},{"comment":"The abstract's phrase \"This in turn will allow us to find the asymptotics\" could be read as unconditional; it would be helpful to mention explicitly that the asymptotic theorem requires the non-degeneracy condition (1.9), whose validity for the general class is presently open.","section":"Abstract and Introduction"},{"comment":"The sufficiency condition (4.43) for the interval-supported case is stated without derivation; a brief indication of where this condition comes from, or a reference to a future work, would improve readability.","section":"Remark 4.2, Eq. (4.43)"},{"comment":"There are occasional spelling and formatting slips (for example, \"offset\" for \"offset\" and some inconsistent spacing in displayed equations) that should be corrected in copy-editing.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a real advance, and the main theorem is honest about being conditional. The 4x4 Riemann–Hilbert formulation for Toeplitz+Hankel determinants with independent Toeplitz and Hankel symbols is new; prior work needed identical symbols or a special class with d(±1)=1. The paper removes the d(±1)=1 condition, solves the model problem for the class d\\tilde d=1, and derives the asymptotics h_{n-1} = -α(0) E(n)/E(n-1)(1+O(e^{-c_1 n})). The derivation is self-contained, uses standard small-norm Riemann–Hilbert theory, and E(n) is expressed through explicit contour integrals of the symbol data rather than fitted to the target determinant. The explicit family (4.38)–(4.39), with the Watson-lemma estimate (4.40), shows the non-degeneracy condition is not vacuous. That is real progress.\n\nThe soft spots are proportionate. The non-degeneracy condition (1.9) is load-bearing: reconstruction of the Y-RHP from the X-RHP uses condition (2.36), which is exactly the lower bound |E(n)| ≥ C r^n. The paper verifies this only for an explicit family; Section 5.5 leaves the characterization of the relevant class C0 open. So the advertised general asymptotics for all smooth d\\tilde d=1 symbols is not established; the theorem is precisely conditional. That is a limitation of scope, not an internal inconsistency, and the authors say so explicitly. The interval-supported Hankel case is also slightly oversold in the abstract: Section 3 reduces to the same model problem for the pair (φ, -\\tilde u), but no asymptotic theorem is proved for that case. Minor issue: some proofs in Section 3 are omitted as \"similar\" to earlier ones; given the paper's length, this is acceptable but worth noting. The citation pattern is fine: [2], [13], and [14] are used accurately, and the connection to the operator-theoretic approach of Basor–Ehrhardt is explained rather than claimed loosely.\n\nWho is this for? Specialists in Riemann–Hilbert methods and asymptotic analysis of structured determinants, and people working on Hankel spectral asymptotics or the zig-zag Ising model. It deserves serious refereeing. The conditional theorem is stated cleanly, the proof is detailed, and the explicit solvable family gives the paper concrete content beyond a general formalism.","headline":"A genuinely new 4x4 Riemann–Hilbert framework for Toeplitz+Hankel determinants with independent symbols, with a main theorem that is honestly conditional on a non-degeneracy hypothesis; deserves a serious referee.","tokens_in":43135,"tokens_out":1943,"would_cite":true,"duration_ms":20977,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["15B05","30E15","35Q15"],"pacs":[],"model":"deepseek-v4-flash","headline":"A 4×4 Riemann–Hilbert problem computes large-$n$ asymptotics of Toeplitz+Hankel determinants whose Toeplitz and Hankel symbols are not assumed related.","keywords":["Toeplitz+Hankel determinants","Riemann-Hilbert problem","nonlinear steepest descent","orthogonal polynomials","asymptotic analysis","Szegő-type symbols","unit circle","Hankel matrices"],"falsifier":"Take the explicit family (4.38)–(4.39) with a parameter choice where the constant $\\kappa$ in (4.41) vanishes, so that $E(n)$ violates condition (1.9), and compute $h_{n-1}=D_n/D_{n-1}$ numerically to high precision for large $n$; if the asymptotic $-\\alpha(0)E(n)/E(n-1)$ still holds, the non-degeneracy condition is not necessary, while a breakdown would confirm it is load-bearing.","tokens_in":42054,"feed_emoji":"📐","tokens_out":10028,"duration_ms":97673,"temperature":0.7,"pith_summary":"The paper develops a 4×4 Riemann–Hilbert method for the large-$n$ asymptotics of determinants of sums of Toeplitz and Hankel matrices, with no a priori relation between the Toeplitz symbol $\\varphi$ and the Hankel symbol $w$. The main theorem gives the ratio of consecutive determinants, $h_{n-1}=D_n/D_{n-1}$, as $-\\alpha(0)E(n)/E(n-1)$ with an exponentially small relative error, where $E(n)$ is an explicit combination of contour integrals and $\\alpha$ is the Szegő function of $\\varphi$. The same model problem emerges when the Hankel weight is supported on an interval inside $(0,1)$, so the analysis covers both geometries. This matters because such determinants control norms of associated orthogonal polynomials and arise, for instance, in the characteristic polynomials of Hankel matrices.","feed_headline":"4x4 Riemann-Hilbert scheme cracks Toeplitz+Hankel asymptotics","feed_subtitle":"The ratio of consecutive determinants is pinned down by an explicit functional up to exponentially small error.","key_machinery":"The load-bearing object is a 4×4 Riemann–Hilbert problem (the $X$-problem) whose jump matrix is an ordinary multiplicative jump; it is obtained by doubling a 2×2 problem with the Carleman shift $z\\mapsto z^{-1}$. Nonlinear steepest descent reduces it to a model Riemann–Hilbert problem on the unit circle. For the solvable class, the condition $d\\tilde d=1$ makes one entry of the model jump vanish, and the jump is factorized using the Szegő functions $\\alpha,\\beta$ of $\\varphi$ and $d$, yielding an explicit model solution $\\Lambda$. The error is controlled by a small-norm problem for $R=S\\Lambda^{-1}$, whose jump is exponentially close to the identity; the leading terms $R_{1,23}$ and $R_{1,43}$ combine into the functional $E(n)$ appearing in the determinant-ratio formula.","core_discovery":"The paper's central claim, Theorem 1.1, is that for smooth nonvanishing Szegő-type symbols $\\varphi$ and $w=d\\varphi$ on the unit circle, with $d\\tilde d=1$ on the circle, the determinant ratio obeys $$h_{n-1}=-\\$\\alpha$(0)\\frac{E(n)}{E(n-1)}\\left(1+O($e^{{-c_1 n}}$)\\right),\\qquad n\\to\\infty,$$ where $E(n)=2\\alpha(0)R_{1,43}(0;n)-C_\\rho(0)R_{1,23}(0;n)$ is built from two exponentially small contour integrals and $c_1>0$. The proof reconstructs the original 2×2 orthogonal-polynomial problem from the 4×4 problem, solves the model problem explicitly for this symbol class, and controls the correction by a small-norm Riemann–Hilbert analysis. As a consequence the paper also obtains large-$n$ asymptotics of the orthogonal polynomials themselves, and for an explicit family shows $D_n\\sim C(-b_1)^n n^{\\alpha_1-1}$, an oscillatory behavior the authors report as confirmed numerically.","pith_inferences":["If the non-degeneracy condition (1.9) is generic, as the explicit example suggests, the functional $E(n)$ may be the universal object controlling smooth Toeplitz+Hankel asymptotics; testing other symbol families would reveal how far that universality reaches.","A direct numerical check of the explicit family should reproduce $-b_1(n/(n-1))^{\\alpha_1-1}$ up to $O(n^{-2})$ at moderate $n$; this would independently verify the error terms in (4.37).","Applying the same 4×4 problem to $D_n(-\\lambda,w;0,0)$ would give Hankel eigenvalue asymptotics, a case the paper identifies as simpler in symbol but harder because the model jump no longer simplifies.","Because the method drops the condition $d(\\pm1)=1$ required by the operator-theoretic approach, overlapping cases provide a built-in consistency test: the two routes must agree wherever both apply."],"forward_implications":["For every symbol pair satisfying Theorem 1.1, $D_n\\neq0$ for all sufficiently large $n$, and the norm $h_{n-1}$ is known up to exponentially small relative error.","Because the same model problem governs the interval-supported Hankel case, the method transfers to $D_n(\\varphi,w;1,s)$ whenever the induced function $d=-\\varphi^{-1}\\tilde u$ satisfies $d\\tilde d=1$ on the unit circle.","For the explicit family (4.38)–(4.39), the determinant has the oscillatory asymptotic $D_n\\sim C(-b_1)^n n^{\\alpha_1-1}$, with the constant $C$ still to be determined.","The asymptotics in Remark 4.1 give explicit expressions for the orthogonal polynomials $P_n$ inside, on, and outside the unit circle, uniformly in $z$."],"supporting_citations":[{"why":"Supplies the class of symbol pairs (1.4) whose model problem the paper solves, and the operator-theoretic formulas to which the result is compared.","marker":"[2]"},{"why":"Introduced the Riemann–Hilbert approach to Toeplitz+Hankel determinants with coinciding symbols and provides the evenness and singularity framework that this work extends to independent symbols.","marker":"[13]"},{"why":"Supplies the small-norm Riemann–Hilbert theory used to prove existence and exponentially small error estimates for the correction R.","marker":"[16, 17]"},{"why":"Provides the standard Riemann–Hilbert formulation of Toeplitz determinants used in the discussion of general offsets r and s.","marker":"[1]"}],"fun_headline_variants":["Exact Toeplitz+Hankel asymptotics via 4x4 RH","4x4 RH solves Toeplitz+Hankel determinant asymptotics","Explicit asymptotics for Toeplitz+Hankel via 4x4 RH","RH steepest descent: exact Toeplitz+Hankel asymptotics","4x4 Riemann-Hilbert pins Toeplitz+Hankel asymptotics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof requires the exponentially small functional $E(n)$ to stay bounded below by $C r^n$ for some $r<1$ as $n\\to\\infty$; if $E(n)$ decays faster than every such bound or has zeros for large $n$, the reconstruction of the 2×2 problem from the 4×4 problem is not justified.","fun_headline_variants_meta":{"raw":{"variants":["Exact Toeplitz+Hankel asymptotics via 4x4 RH","4x4 RH solves Toeplitz+Hankel determinant asymptotics","Explicit asymptotics for Toeplitz+Hankel via 4x4 RH","RH steepest descent: exact Toeplitz+Hankel asymptotics","4x4 Riemann-Hilbert pins Toeplitz+Hankel asymptotics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001114,"raw_usage":{"total_tokens":5051,"prompt_tokens":1141,"completion_tokens":3910,"prompt_tokens_details":{"cached_tokens":1024},"prompt_cache_hit_tokens":1024,"prompt_cache_miss_tokens":117,"completion_tokens_details":{"reasoning_tokens":3802}},"tokens_in":117,"tokens_out":3910,"duration_ms":283841,"temperature":1.0,"reasoning_tokens":3802,"cache_read_input_tokens":1024,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:31:00.772730+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the explicit family (4.38)–(4.39) with a parameter choice where the constant $\\kappa$ in (4.41) vanishes, so that $E(n)$ violates condition (1.9), and compute $h_{n-1}=D_n/D_{n-1}$ numerically to high precision for large $n$; if the asymptotic $-\\alpha(0)E(n)/E(n-1)$ still holds, the non-degeneracy condition is not necessary, while a breakdown would confirm it is load-bearing.","supporting_citations":[{"cited_title":"Asymptotic formulas for determinants of a special class of Toeplitz + Hankel matrices","cited_arxiv_id":"1603.00506","evidence_quote":"Supplies the class of symbol pairs (1.4) whose model problem the paper solves, and the operator-theoretic formulas to which the result is compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard Riemann–Hilbert formulation of Toeplitz determinants used in the discussion of general offsets r and s."}],"review_version":1}