{"id":"d60aba78-1b26-4d12-89a5-8e2d13c1469e","arxiv_id":"2505.10622","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper presents a third-order nonlinear differential equation (Eq. 6) whose imaginary part, integrated over x, yields the correlation kernel K(E,E') for Schrodinger-type random matrix models, generalizing the Gel'fand-Dikii equation.","lead":"This paper derives a new differential equation for a quantity that, when integrated, gives the two-point correlation kernel of eigenvalues in many random-matrix-type systems. It generalizes the classic Gel'fand-Dikii equation for the local density of states to correlations between two different energies.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No uniqueness proof for the BVP defining S; numerical evidence is limited to near-diagonal Airy, so Eq. (6) is not established as a general definition of K(E,E′).","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing gap: no existence or uniqueness proof for the boundary-value problem defining S, and numerical support only for Airy with E′ close to E. The manuscript's own statement that the numerical method becomes unstable beyond about 10% energy separation is an explicit internal limitation, which reinforces rather than replaces the reader's concern. The derivation of Eq. (6) from the product of wavefunctions is algebraically sound, and the numerical checks for the diagonal case and for small separations are credible evidence that the equation is not vacuous. However, the paper's central claim is that Eq. (6) 'gives a definition' of ψ(E,x)ψ(E′,x) through the imaginary part of S; a definition requires that the specified differential equation plus boundary conditions have a unique solution with the claimed imaginary part in the intended generality. That well-posedness is unproven, and the tested regime is narrow. This does not invalidate the paper, but it supports the reader's conditional verdict: the equation is promising and correct where checked, yet the general claim requires additional proof or a much wider numerical demonstration. I therefore see no reason to change the reader's CONDITIONAL verdict.","tokens_in":7777,"tokens_out":9888,"duration_ms":105335,"concrete_test":"Solve the coupled system (3)+(6) for the Airy potential u(x)=-x on a large interval [-L,0] with E=1, ℏ=1, using a high-order spectral collocation method, and continue δ=(E′-E)/E from 0.01 to 0.5. Use two independent initial guesses, e.g., ψ(E,x)ψ(E′,x) and its negative, and compare the converged Im S pointwise with Ai(-(x+E))Ai(-(x+E′)) on the grid, as well as the integrated kernel via Eq. (5) against the Christoffel-Darboux form (21). If both guesses converge to the same Im S matching the product to 1e-6 for δ≥0.2, the concern is mitigated; if they converge to different branches or fail to match, the claimed definition is not supported for general energy separations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. (6), together with Eq. (5) and with R1,R2 solving Eq. (3), defines the correlation kernel through Im S for a wide class of potentials. The derivation only shows that S=ψ(E,x)ψ(E′,x) satisfies Eq. (6); it does not show that the boundary-value problem used numerically selects that solution, nor that its imaginary part is unique. The paper itself states that 'going much beyond a 10% difference becomes hard to control' in the numerical solutions section, and the only nontrivial test is the Airy potential with E′ within 1% of E. No existence or uniqueness theorem is supplied for the third-order nonlinear BVP, and the asymptotic boundary condition (19) with c2=0 is not shown to fix the correct branch. Because Eq. (6) is nonlinear and couples S to independently determined R1,R2, spurious branches with the same boundary data are not excluded. If such a branch exists, Im S need not equal ψ(E,x)ψ(E′,x), and Eq. (5) would not yield K(E,E′). This is the load-bearing gap behind the paper's central definitional claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter derives a third-order nonlinear differential equation, Eq. (6), for a two-energy object S-hat(E,E',x) associated with the one-dimensional Schrödinger operator H = -hbar^2 d^2/dx^2 + u(x). The derivation is direct algebra: starting from S = psi(E,x) psi(E',x), the Schrödinger equation is used to eliminate derivatives of the wavefunctions, and multiplication by 2S yields Eq. (6), which couples S to the diagonal resolvents R-hat(E,x) and R-hat(E',x) that separately satisfy the Gel'fand-Dikii equation (3). In the diagonal limit E'=E the equation reduces to Eq. (7), the differentiated Gel'fand-Dikii equation. The paper proposes (6), together with (5) and the large-|x| boundary condition (19), as a wavefunction-free definition of the correlation kernel K(E,E'), and it tests this on the Bessel (u=0) model, where an exact solution (11) is given, and on the Airy (u=-x) model, where a numerical boundary-value solution is reported to reproduce psi(E,x) psi(E',x) up to a normalization for energy differences up to about 10%. The closing remarks suggest applications to random-matrix two-point functions and to the W_{g,n} / Weil-Petersson volumes program of recent work.","tokens_in":8057,"tokens_out":31237,"duration_ms":276353,"significance":"If the defining property is established, Eq. (6) is a genuinely useful non-diagonal generalization of the Gel'fand-Dikii equation: it would permit perturbative hbar-expansions and non-perturbative construction of K(E,E') and of n-point correlation functions without explicit wavefunctions, with plausible applications to the volumes program of Refs. [15,16]. The algebraic derivation is clean and correct, and the diagonal reduction is a valid consistency check. The exact Bessel solution (11) and the Airy numerics (Figs. 1-3) give real supporting evidence for the product form of Im S-hat in the tested regimes. What is missing is the other half of the definitional claim: the paper shows that the product S=psi psi' satisfies (6), but it does not show that the boundary-value problem selects that branch, nor does it establish uniqueness of Im S-hat; the numerical evidence is confined to near-diagonal Airy energies and to the exactly solvable u=0 case, and the paper itself states that the numerical problem becomes hard to control beyond a 10% energy separation. The core idea is plausible and likely correct, but the advertised generality to a 'wide class of models' is not yet supported.","major_comments":[{"comment":"Eq. (6) is derived only as a necessary condition: the paper shows that if S = psi(E,x) psi(E',x), then (6) holds. The central claim, stated in the abstract and the closing remarks, that (6) together with (5) and the boundary condition (19) defines a unique object whose imaginary part is (up to normalization) psi(E,x) psi(E',x), requires existence and uniqueness for this third-order nonlinear boundary-value problem, and none is supplied. The numerical support is limited to the Airy potential with |E-E'|/E approximately less than 0.1, as the paper itself acknowledges ('going much beyond a 10% difference becomes hard to control'), plus the exactly solvable u=0 case; both test potentials are exactly solvable, so the generality claim is not exercised. I ask the authors either to prove that the boundary-value problem selects the product branch (for instance by showing uniqueness of the solution or of its imaginary part) or to reframe the definitional claim explicitly as a conjecture and add a test on a potential that is not exactly solvable.","section":"Derivation of the Equation; Numerical Solutions; Closing Remarks"},{"comment":"The integration constant c2 in the Bernoulli-type solution (17) is set to zero by fiat in Eq. (19). For E'=E this choice recovers the diagonal leading resolvent, but for E' != E there is no argument that the choice c2=0 (equivalently G(x)=1) is the unique boundary data whose continuation through (6) yields Im S-hat = +/- pi hbar psi(E,x) psi(E',x). Because (6) is nonlinear, different choices of c2 generically produce different members of the solution family, and the paper does not demonstrate that the physical branch is selected independently of this choice. This is the off-diagonal analogue of fixing the right-hand side of (3) to unity, and it needs a justification; as written, the branch selection is an additional input rather than a consequence of the equation.","section":"Leading Solutions, Eqs. (17)-(19)"},{"comment":"There is a sign-convention inconsistency among Eqs. (2), (4), (5), (10), and (11). With the stated normalization psi = A sin(sqrt(E) x / hbar), A = (pi hbar)^{-1/2} E^{-1/4}, the second equality in (10) gives Im R-hat = + pi hbar psi^2, so Eq. (2) produces a negative spectral density, contradicting Eq. (9) (which is positive). Similarly, Eq. (11) gives Im S-hat = + pi hbar psi psi', so Eq. (5) yields K(E,E') = - integral psi psi', which is the opposite sign of Eq. (4). The resolution is presumably a definite E +/- i0 convention for the resolvent (the branch sqrt(-E) = -i sqrt(E) would flip both signs), but as written the equations are mutually inconsistent, and the numerical statement 'the imaginary part ... coincides with psi psi'' is sign-ambiguous. Please fix one convention and propagate it consistently through Eqs. (10)-(12) and the numerical comparison.","section":"The Simplest Model, Eqs. (2), (5), (10)-(12)"}],"minor_comments":[{"comment":"The integration limits in Eq. (4) are written as integral from 0 to -infinity; the standard notation integral from -infinity to 0 would be clearer and less error-prone.","section":"Introduction, Eq. (4)"},{"comment":"The sentence 'The leading part of (7) shows that bR starts out as -1/2 (-E)^{1/2}' appears to contain two slips: for u=0 the hbar -> 0 limit of (3) gives R_0 = +/- 1/2 (-E)^{-1/2}, and the normalization is fixed by (3), not by (7), whose hbar -> 0 limit with u=0 merely gives R' = 0. The exponent and the reference to (7) should be corrected.","section":"The Simplest Model"},{"comment":"Please specify what is plotted in the comparisons. The exact products psi(E,x)^2 and psi(E,x) psi(E',x) for the Airy potential are oscillatory with infinitely many nodes in the classically allowed region, while the boundary condition (19) is described as capturing only the average fall-off; a statement of whether the agreement is pointwise or at the level of the envelope, and of the error metric used, is needed to interpret 'coincides with psi psi''.","section":"Numerical Solutions, Figs. 1-3"},{"comment":"The boundary-value problem is not fully specified. Please list the nine boundary conditions (which functions and which derivatives, at which endpoints) imposed in the bvp4c runs, in addition to the leading form (19); this is needed for reproducibility of the results in Figs. 1-3.","section":"Numerical Solutions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has a sound core: the derivation of Eq. (6) is correct, the diagonal reduction works, and the numerical checks, though narrow, support the product form in the tested regime. The main editorial question is whether a Letter-length paper can carry a definitional claim of this breadth without a uniqueness argument for the boundary-value problem and with a sign error in the central example section. I believe the gaps are fixable within scope: the authors should either supply a uniqueness/branch-selection argument or explicitly downgrade the claim, and they must repair the sign conventions in Eqs. (10)-(12). I would like the revised version to be seen by a referee again, with particular attention to whether the boundary-value problem for a non-solvable potential actually selects Im S-hat = +/- pi hbar psi psi' beyond the near-diagonal Airy case. The balance of the paper (introductory framing versus new content) is appropriate for a Letter, and the connections to the volumes program are potentially interesting but speculative as they stand."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper's main result is Eq. (6), a third-order nonlinear ODE for S(E,E',x) that reduces to the Gel'fand-Dikii equation when E=E'. The derivation is correct and simple: multiply the two Schrödinger equations for ψ(E) and ψ(E') and eliminate derivatives. What's new is that this gives an off-diagonal generalization of the resolvent equation, and the paper shows the exact solution for u=0 and numerical agreement for the Airy potential with E' near E. The numerics include a Christoffel-Darboux check of the integrated kernel, which is convincing in the tested range.\n\nThe soft spot is real: the claim that Eq. (6) plus boundary conditions 'defines' ψ(E)ψ(E') via Im S is not backed by an existence or uniqueness proof. The equation is nonlinear and third-order, and the asymptotic condition (19) with c2=0 is chosen because it gives the leading WKB product. Nothing rules out other solutions with the same boundary data whose imaginary part would not equal the product. The numerical evidence only covers E' within 1%-10% of E, and the paper itself notes that going beyond 10% is hard to control. So the paper's title 'for a class of correlation kernels' is fair, but the abstract's 'wide class of models' is currently supported by only two exactly solved cases (Bessel and Airy). That is not a fatal flaw—it is a solid contribution with limited verification.\n\nI would like to see the BVP question addressed, at least heuristically: is there an argument that the physical branch is selected by requiring S → ψψ as ℏ→0? The leading-order Bernoulli solution already hints at that, but it is not made rigorous. Also, more numerics for different potentials (harmonic oscillator, double well) would strengthen the claim.\n\nOverall, this is worth a serious referee. The equation is new, the derivation correct, and the applications to higher-genus corrections and volumes are plausible. I would accept it for review and ask for either a uniqueness argument or softer claims.","headline":"A correct and useful new off-diagonal generalization of the Gel'fand-Dikii equation, but the central definitional claim is only checked numerically near the diagonal and lacks a uniqueness argument.","tokens_in":8522,"tokens_out":2262,"would_cite":false,"duration_ms":22352,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A new differential equation defines correlation kernels without explicit wavefunctions.","keywords":["correlation kernel","Schrödinger equation","Gel'fand-Dikii equation","off-diagonal resolvent","random matrix theory","determinantal point process","Airy model","Bessel model"],"falsifier":"Take an exactly solvable potential with known wavefunction products, such as the harmonic oscillator, and solve Eq. (6) at an energy separation beyond the roughly 10% range tested for Airy; if the imaginary part of the selected solution fails to match the known product, or if the result depends on the integration constant $c_2$ in the boundary condition, the proposed definition is false.","tokens_in":7585,"feed_emoji":"⚛️","tokens_out":8658,"duration_ms":74871,"temperature":0.7,"pith_summary":"This paper proposes a new differential equation—Eq. (6)—for an object $\\widehat{S}(E,E',x)$ whose $x$-integral gives the correlation kernel $K(E,E')$ of a one-dimensional Schrödinger model $H=-\\hbar^2\\partial_x^2+u(x)$. The equation is an off-diagonal generalization of the classic Gel'fand–Dikii equation for the diagonal resolvent, and it reduces to that equation when $E=E'$. If it holds generally, it offers a way to compute $n$-point functions and gap probabilities without first solving for wavefunctions. The author verifies the claimed behavior for the Airy potential numerically and for the Bessel potentials analytically.","feed_headline":"Equation yields correlation kernels without wavefunctions","feed_subtitle":"Its imaginary part reproduces wavefunction products, enabling correlators without solving for states.","key_machinery":"The central object is $\\widehat{S}(E,E',x)$, an off-diagonal generalization of the diagonal resolvent $\\widehat{R}(E,x)$. The equation is obtained by starting from $S=\\psi(E,x)\\psi(E',x)$, using the Schrödinger equation to eliminate third derivatives of the wavefunctions, and multiplying by $2S$ so that the leftover products $\\psi_1\\psi_2'$ and $\\psi_2\\psi_1'$ become combinations of the known diagonal resolvents $R_1$ and $R_2$. What carries the argument is that all explicit wavefunction dependence disappears, leaving the closed system Eqs. (3) and (6). At leading order in $\\hbar$, the equation becomes a Bernoulli equation whose explicit solution supplies the large-$x$ boundary condition used to solve the full problem numerically.","core_discovery":"On the paper's own terms, the central discovery is that the product of two Schrödinger wavefunctions $\\psi(E,x)\\psi(E',x)$—not just the square $\\psi(E,x)^2$—obeys a closed third-order nonlinear differential equation, written in Eq. (6) for $\\widehat{S}(E,E',x)$, using only the potential $u(x)$ and the two diagonal resolvents $R_1$ and $R_2$. The physical solution is claimed to have imaginary part equal to $\\pi\\hbar\\,\\psi(E,x)\\psi(E',x)$ up to normalization, so integrating it through Eq. (5) yields the correlation kernel $K(E,E')$. The paper shows that Eq. (6) collapses to the Gel'fand–Dikii equation when $E=E'$, obtains an exact solution for the Bessel/trigonometric potentials, and reproduces wavefunction products numerically for the Airy potential at nearby energies.","pith_inferences":["If the boundary-value problem is unique beyond the tested cases, Eq. (6) should work as a numerical kernel generator for potentials where explicit wavefunctions are unknown; this is a testable extension the paper does not carry out.","The leading-order equation's link to complex Bernoulli equations and Lie–Hamilton systems hints at an integrable structure that could produce closed-form kernels for other potentials; the paper mentions but does not develop this.","Since the kernel also controls Fredholm-determinant gap probabilities, an off-diagonal version could enable two-parameter or joint gap statistics, a consequence the author leaves implicit."],"forward_implications":["For any potential $u(x)$ where the coupled system can be solved, the correlation kernel $K(E,E')$ follows from Eq. (5) without explicit wavefunctions; the imaginary part of $\\widehat{S}$ carries the required information.","A systematic small-$\\hbar$ expansion of $\\widehat{S}(E,E',x)$ becomes available, giving off-diagonal and $n$-point generalizations of the resolvent data that, in the diagonal case, produce the volumes $V_{g,1}$; the author notes the route to $V_{g,n}$.","In the exactly solvable Bessel case, the equation yields the closed-form kernel of Eq. (12), whose large-energy limit reproduces the universal connected two-point function.","Numerically solving the system as a boundary-value problem reproduces both the diagonal resolvent and off-diagonal wavefunction products for the Airy model, providing a non-perturbative computational route."],"supporting_citations":[{"why":"supplies the classic Gel'fand–Dikii equation for the diagonal resolvent that Eq. (6) generalizes","marker":"[7]"},{"why":"defines the Airy model used as the main numerical testbed for the proposed kernel equation","marker":"[9]"},{"why":"also uses the Airy model, alongside the Wigner semi-circle edge expansion","marker":"[10]"},{"why":"identifies the Bessel models whose exact wavefunctions and kernels are used to obtain the analytic solution","marker":"[11]"},{"why":"provides the Lie–Hamilton/Bernoulli solution used to write the leading-order complex solution for $\\widehat{S}_0$","marker":"[13]"},{"why":"shows the diagonal resolvent expansion yields volumes $V_{g,1}$, motivating the off-diagonal generalization to $V_{g,n}$","marker":"[15]"},{"why":"maps the relevant Schrödinger problem to Bessel's equation, supplying the wavefunction normalization used in the exact kernel","marker":"[23]"}],"fun_headline_variants":["Third-order equation for wavefunction products, no states needed","Two-energy Gel'fand–Dikii equation produces correlation kernels","Correlation kernels without solving for wavefunctions","New ODE reproduces wavefunction products, enabling kernels"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the boundary-value problem formed by Eq. (6), with $R_1$ and $R_2$ fixed by Eq. (3) and the asymptotic condition (19), has a unique solution whose imaginary part reproduces, up to normalization, the product wavefunction $\\psi(E,x)\\psi(E',x)$ for general potentials; the paper gives no existence or uniqueness proof and verifies this numerically only for small energy separations.","fun_headline_variants_meta":{"raw":{"variants":["Third-order equation for wavefunction products, no states needed","Two-energy Gel'fand–Dikii equation produces correlation kernels","Correlation kernels without solving for wavefunctions","New ODE reproduces wavefunction products, enabling kernels"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000923,"raw_usage":{"total_tokens":3920,"prompt_tokens":873,"completion_tokens":3047,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":489,"completion_tokens_details":{"reasoning_tokens":2980}},"tokens_in":489,"tokens_out":3047,"duration_ms":24047,"temperature":1.0,"reasoning_tokens":2980,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:07:06.732225+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take an exactly solvable potential with known wavefunction products, such as the harmonic oscillator, and solve Eq. (6) at an energy separation beyond the roughly 10% range tested for Airy; if the imaginary part of the selected solution fails to match the known product, or if the result depends on the integration constant $c_2$ in the boundary condition, the proposed definition is false.","supporting_citations":[{"cited_title":"Dyson, J","cited_arxiv_id":null,"evidence_quote":"supplies the classic Gel'fand–Dikii equation for the diagonal resolvent that Eq. (6) generalizes"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the Airy model used as the main numerical testbed for the proposed kernel equation"},{"cited_title":"In the form given in Eq","cited_arxiv_id":null,"evidence_quote":"also uses the Airy model, alongside the Wigner semi-circle edge expansion"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"identifies the Bessel models whose exact wavefunctions and kernels are used to obtain the analytic solution"},{"cited_title":"These models appear naturally in random matrix mod- els of Wishart form which have a “hard” edge [6, 24, 25]","cited_arxiv_id":null,"evidence_quote":"provides the Lie–Hamilton/Bernoulli solution used to write the leading-order complex solution for $\\widehat{S}_0$"},{"cited_title":"Solutions by quadratures of complex Bernoulli differential equations and their quantum deformation","cited_arxiv_id":"2312.16586","evidence_quote":"shows the diagonal resolvent expansion yields volumes $V_{g,1}$, motivating the off-diagonal generalization to $V_{g,n}$"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"maps the relevant Schrödinger problem to Bessel's equation, supplying the wavefunction normalization used in the exact kernel"}],"review_version":1}