{"id":"1f903e31-a0c0-4fdd-9790-0c2c745940ea","arxiv_id":"2412.08825","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"A discrete scattering overlap relation computes phase shifts from potential matrix elements in square-integrable bases, demonstrated on square-well and nucleon-nucleon potentials.","lead":"This paper derives an overlap formula for computing scattering phase shifts directly in square-integrable bases, the kind used in shell-model nuclear calculations. The formula avoids matching wave functions at large distances, a known bottleneck in many-body scattering calculations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (36) is exact only if V_{n,n'}=0 beyond Npot and u_n is free beyond; the paper concedes no such theorem exists, and its own Malfliet-Tjon results (Npot=50, residual oscillations) show the assumption is only approximate for realistic potentials.","rationale":"The paper's central contribution is Eq. (36), a discrete analog of the coordinate-space overlap relation. The derivation via telescoping sums is correct given the stated assumptions: tridiagonal T, truncated V, and free asymptotic form u_n=Af_n+Bg_n for n>Npot. The square-well tests (Npot=10, l=0-3) are convincingly accurate, and the use of the Lippmann-Schwinger equation rather than diagonalization is a reasonable numerical choice. The load-bearing weakness is not the algebra but the status of the truncation assumption. As the paper itself concedes in Sec. II, no theorem guarantees V_{n,n'}->0 in an oscillator basis; for the Malfliet-Tjon potential with its strong short-range repulsion, acceptable results require Npot=50 and visible oscillations remain, with smaller spaces making oscillations worse. Because the paper sets Nmax=Npot, the wavefunction coefficients themselves are generated in the truncated space, so the numerical examples do not establish convergence to the infinite-space limit that Eq. (36) assumes. The comparison of LS and diagonalization phase shifts only shows two finite-space methods agree with each other. Therefore the headline claim that Eq. (36) 'gives the scattering phase shift' is conditional on a cutoff assumption that is plausible for soft interactions but unproven and evidently strained for realistic hard-core potentials. A systematic Npot-convergence study, together with a direct check of the asymptotic form u_n≈Af_n+Bg_n, would settle whether the formula is exact in practice or only approximate. Until then, CONDITIONAL is the right verdict.","tokens_in":13217,"tokens_out":12411,"duration_ms":123968,"concrete_test":"For the Malfliet-Tjon potential, compute the SOR phase shifts via Eq. (36) for Npot = 20, 50, 100, 200 with b chosen by the ground-state-energy minimum at each Npot, and compare to coordinate-space Numerov phase shifts over E_lab = 0–300 MeV. Also compute the residual |u_n - (A f_n + B g_n)|/|u_n| for n just above Npot using the constants A and B extracted from Eq. (36). If the phase-shift error does not decrease monotonically or the residual stays O(1) as Npot grows, the truncation assumption fails and the central claim must be weakened to 'approximate for finite Npot'.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central derivation of Eq. (36) is algebraically sound, but its exactness rests on the J-matrix truncation assumption V_{n,n'}=0 for n,n'>Npot, together with the asymptotic form u_n=A f_n+B g_n for n>Npot (Eq. 25). Section II explicitly concedes there is no theorem that V_{n,n'} must eventually vanish as n grows. The numerical section sets Nmax=Npot, so the same finite space supplies both the potential matrix and the wavefunction; the LS solutions are therefore not demonstrated to be the exact infinite-space coefficients, and the comparison to matrix diagonalization checks only internal consistency. For the Malfliet-Tjon potential, which has a strong repulsive core, Npot=50 is needed and 'slight oscillations' persist in the phase shifts; smaller Npot increases the oscillations. This indicates Eq. (36) is at best an approximate overlap relation for realistic potentials, not the exact relation the abstract implies. The Daejeon-16 test additionally omits the coupled 3D1 channel, so it does not independently validate the formula in a realistic coupled-channel setting.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper derives a scattering overlap relation (SOR) for the J-matrix method. In Section III A the authors rederive the familiar coordinate-space integral relation for the phase shift from a Green's theorem argument. In Section III B they repeat the construction in an orthonormal L2 basis whose kinetic-energy matrix is tridiagonal. Assuming the potential matrix is truncated, V_{n,n'} = 0 for n,n' > Npot, and that the scattering coefficients have the free tail u_n = A f_n + B g_n for n > Npot, they obtain Eq. (36): tan δ = -[Σ f_n V_{n,m} u_m]/[α0 u0 + Σ g_n V_{n,m} u_m]. The derivation is a telescoping-sum argument analogous to the Wronskian/Casoratian manipulation in coordinate space. The paper then demonstrates the formula for square-well scattering (l=0,...,3) and for 3S1 nucleon-nucleon phase shifts using Malfliet-Tjon and Daejeon-16 interactions, with the oscillator length chosen by minimizing the ground-state energy.","tokens_in":13424,"tokens_out":9437,"duration_ms":91569,"significance":"The algebraic content of the paper is a useful and largely correct addition to the J-matrix literature. Eq. (36) is transparently derived, is a genuine discrete analog of the coordinate-space overlap relation, and has the attractive feature that the phase shift is obtained from potential matrix elements and coefficients only up to Npot, without an explicit match at large n. The authors are appropriately explicit about the main assumption: there is no theorem that V_{n,n'} eventually vanishes, so the relation is exact only for a model in which the potential is exactly truncated in the basis. The numerical demonstrations are consistent with this: the square-well results reproduce analytic phase shifts, and the LS and diagonalization routes agree internally. The paper does not claim to fit any phase-shift data; b and Npot are chosen by energy minimization and practical convergence, which is a strength. The main weakness is that the numerical tests set Nmax = Npot and therefore do not separately validate the exact infinite-space form of Eq. (36).","major_comments":[{"comment":"The exactness of Eq. (36) rests on the assumptions (i) V_{n,n'} = 0 for n,n' > Npot and (ii) u_n = A f_n + B g_n for n > Npot. The paper explicitly concedes in Section II that there is no theorem guaranteeing (i). The numerical section then sets Nmax = Npot, so the coefficients u_m entering Eq. (36) are the solutions of a truncated LS equation, not the exact coefficients of the infinite-space truncated-potential problem. The excellent square-well agreement therefore validates the combined truncation-plus-formula procedure, not Eq. (36) with exact input coefficients. To support the statement 'we no longer need Nmax > Npot', the authors should either construct the exact finite-dimensional tail-matched solution or report a convergence study of Eq. (36) with Npot fixed and Nmax increasing.","section":"III B, Eq. (36), and Sec. IV"},{"comment":"For the Malfliet-Tjon potential, the method requires Npot = 50 and even then 'slight oscillations' persist in the phase shifts, with oscillations growing for smaller Npot. This is acknowledged, but it means the practical claim that the SOR works in a small model space is demonstrated only for soft interactions. A quantitative convergence test (phase shift versus Npot at a few energies, with Nmax = Npot and with Nmax > Npot) would allow the reader to judge how much of the residual error is due to the potential truncation versus the finite-Nmax generation of u_m.","section":"IV, Fig. 2"}],"minor_comments":[{"comment":"The word 'reproducs' should be 'reproduces', and in Section II 'no long variational' should be 'no longer variational'.","section":"IV"},{"comment":"References [19] and [44] are the same article (Flores and Nollett) and should be merged or cross-referenced.","section":"References"},{"comment":"Because the 3D1 coupling is omitted, the resulting phase shifts should be described as the single-channel projection of the interaction, not as the physical isoscalar s-wave phase shifts of the full Daejeon-16 potential.","section":"IV, Daejeon-16 example"},{"comment":"The sentence 'we no longer need Nmax > Npot' should be qualified: Eq. (36) requires only coefficients up to Npot, but generating accurate coefficients may still require a larger space or an exact tail-matched treatment.","section":"III B"},{"comment":"The caption describes symbols and lines, but a legend would improve readability, especially since four angular momenta are shown.","section":"Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable contribution to the J-matrix scattering literature. The main gap is between the exact conditional statement of Eq. (36) and the numerical practice with Nmax = Npot; a revision that adds a convergence study and softens the exactness wording would settle the concern. I do not see grounds for rejection, but the load-bearing validation issue warrants major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThis is a clean, modest methods paper. The new thing is Eq. (36), a discrete-basis analog of the coordinate-space overlap relation that gives the phase shift directly from potential matrix elements and wave-function coefficients in a basis with tridiagonal kinetic energy. The derivation is transparent and self-contained. I believe this exact discrete formula is not in the cited J-matrix literature; the usual extraction uses boundary conditions or Green's functions. So the novelty is real, if modest.\n\nThe paper does well on its own terms. The square-well results for l=0 through 3 match analytic phase shifts, and the authors cross-checked Lippmann-Schwinger solutions against matrix diagonalization, which is good practice. The appendix's observation that phase-shift quality correlates with the ground-state energy for oscillator-length choice is useful, and the choice of b by minimizing the ground-state energy is a sensible heuristic.\n\nThe soft spots are about the numerical demonstration, not the algebra. Eq. (36) is exact only if the potential matrix elements vanish beyond Npot and the scattering coefficients take the free form there. The authors concede no such theorem exists. For the square well this is satisfied by construction. But for the Malfliet-Tjon potential, which has a strong repulsive core, they need Npot=50 and still see residual oscillations; smaller Npot makes them worse. That means the formula is an approximation for realistic potentials, and the paper would be stronger with convergence studies and error estimates, plus a direct comparison against the standard J-matrix phase-shift formula on the same examples. The Daejeon-16 test omits the coupled 3D1 channel, so it doesn't independently validate the method in a coupled-channel setting.\n\nThe stress-test note is essentially right: the truncation assumption is load-bearing and the paper is honest about it. But this is a standard working assumption in J-matrix theory, and the paper's transparency is a point in its favor.\n\nOverall, this deserves a serious referee. The central result is a useful addition to the toolkit. I'd recommend the authors add convergence/error analysis and a comparison with existing J-matrix methods before it is final.","headline":"Clean derivation of a new discrete overlap relation for J-matrix phase shifts, with honest but incomplete numerical support.","tokens_in":13978,"tokens_out":3322,"would_cite":false,"duration_ms":31467,"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 formula extracts scattering phase shifts directly from square-integrable basis calculations, bypassing large-distance wave-function matching.","keywords":["J-matrix method","scattering phase shifts","overlap relations","square-integrable basis","harmonic oscillator basis","configuration-interaction shell model","potential scattering","nucleon-nucleon interaction"],"falsifier":"Take a potential whose harmonic-oscillator matrix elements decay slowly or oscillate (e.g., a long-range $1/r$ tail or a hard-core potential) and compute $\\tan \\delta$ from Eq. (36) at a fixed energy for increasing $N_{\\rm pot}$; if the phase shift does not converge or changes by more than the target accuracy when $N_{\\rm pot}$ is doubled, the central exactness claim for that potential fails. Alternatively, construct a model where $V_{n,n'}$ is nonzero for $n,n' > N_{\\rm pot}$ and show the omitted terms change the result.","tokens_in":117,"feed_emoji":"⚛️","tokens_out":4390,"duration_ms":56985,"temperature":0.7,"pith_summary":"The paper derives a discrete analogue of the coordinate-space overlap relation for scattering phase shifts that works directly in any orthonormal square-integrable basis with tridiagonal kinetic energy. The central formula, Eq. (36), expresses $\\tan \\delta$ as a ratio of sums of potential matrix elements and wave-function coefficients only up to a cutoff $N_{\\rm pot}$, avoiding the need to resolve the wave function at large distances. This matters because many-body methods such as the configuration-interaction shell model work in harmonic-oscillator bases, where coordinate-space asymptotics are hard to extract. Demonstrations on square-well, Malfliet-Tjon, and Daejeon-16 potentials reproduce single-channel phase shifts, with the soft Daejeon-16 interaction needing only $N_{\\rm pot}=5$. The authors present the result as a first step toward scattering and reactions in square-integrable many-body frameworks.","feed_headline":"Phase shifts from overlap relations in square-integrable bases","feed_subtitle":"The J-matrix trick needs only potential matrix elements, opening scattering to CI shell-model frameworks.","key_machinery":"The central object is the scattering overlap relation (SOR) in a discrete $L^2$ basis. It is built from the Casoratian $f_N g_{N+1} - f_{N+1} g_N$, the discrete analog of the Wronskian, together with the free tridiagonal recursion for kinetic energy and the inhomogeneity $\\alpha_0 = (T_{0,0}-E)g_0 + T_{0,1}g_1$ that defines the irregular solution. These ingredients let the difference between the interacting and free equations telescope to a boundary term at $N$, producing Eq. (36) without ever evaluating the wave function at large radius. The paper also uses the Lippmann-Schwinger equation to generate the coefficients $u_n$ with fine control of energy.","core_discovery":"In the J-matrix setting, the scattering phase shift can be computed from an overlap relation that uses only matrix elements of the potential and the coefficients of the scattering wave function in the region where the potential is active. Concretely, if the potential matrix is truncated at $N_{\\rm pot}$ and the asymptotic free solutions $f_n$ and $g_n$ (regular and irregular coefficients of the tridiagonal free problem) are known, then\n$$\\tan \\delta = -\\frac{\\sum_{n=0}^{N_{\\rm pot}} f_n \\sum_m V_{n,m} u_m}{\\alpha_0 u_0 + \\sum_{n=0}^{N_{\\rm pot}} g_n \\sum_m V_{n,m} u_m},$$\nwhere $\\alpha_0 u_0$ is the discrete inhomogeneity at the first basis state. The paper argues this is the direct analogue of the coordinate-space Green's theorem relation $\\tan \\delta = -\\langle f|V|u\\rangle / (g(0)u'(0)-u(0)g'(0)+\\langle g|V|u\\rangle)$ and that it inherits the robustness of overlap relations to errors in the asymptotic wave function. Numerical tests show it reproduces analytic square-well phase shifts for $\\ell=0,\\dots,3$ and experimental $^3S_1$ phase shifts for two nucleon-nucleon interactions, with only a modest model space required for a soft interaction.","pith_inferences":["If the truncation assumption fails for realistic long-range interactions, Eq. (36) becomes an approximation; a convergence diagnostic would be to compare results as $N_{\\rm pot}$ increases or to estimate the omitted tail of $V_{n,m}$.","The same telescoping argument could likely be generalized to coupled channels, involving coupled Casoratians, which the authors list as future work.","The poor performance for the hard-core Malfliet-Tjon potential suggests the method's efficiency depends strongly on the softness of the interaction, so many-body applications may need softened potentials to get small-$N_{\\rm pot}$ phase shifts.","One could test the method on a potential with known slowly decaying harmonic-oscillator matrix elements to quantify the breakdown regime of the central truncation assumption."],"forward_implications":["In any square-integrable basis with tridiagonal kinetic energy (harmonic oscillator, Laguerre), phase shifts can be extracted from a small interior region of the Hamiltonian matrix.","Soft interactions such as Daejeon-16 allow accurate phase shifts with $N_{\\rm pot}=5$, suggesting that truncated many-body spaces may retain scattering information.","The method extends the coordinate-space overlap-integral technique to configuration-interaction shell-model frameworks without requiring a coordinate-space wave function at large $r$.","Since only potential matrix elements and local coefficients enter, the approach naturally adapts to ab initio interactions given as matrix elements in harmonic-oscillator space.","Ground-state energy minimization may serve as a practical guideline for choosing the oscillator length parameter $b$ in truncated calculations, as observed in the appendix."],"supporting_citations":[{"why":"Establishes the J-matrix method: kinetic energy is tridiagonal in an $L^2$ basis, the starting framework for the derivation.","marker":"[6]"},{"why":"Provides the regular and irregular free-solution coefficients and the inhomogeneous equation that defines $g_n$ and $\\alpha_0$.","marker":"[22]"},{"why":"Supplies the Casoratian and the J-matrix treatment of asymptotics used to identify the boundary term in the telescoping sum.","marker":"[24]"},{"why":"Gives existing J-matrix phase-shift extraction methods that the new overlap relation is compared against or extends.","marker":"[10]"},{"why":"Provides the coordinate-space integral relations that the paper generalizes to discrete basis representations.","marker":"[17]"},{"why":"Shows prior use of overlap integrals for nuclear reactions, motivating the discrete analogue.","marker":"[15]"},{"why":"Supplies the experimental nucleon-nucleon phase shifts used as the benchmark for the $^3S_1$ demonstration.","marker":"[47]"},{"why":"Defines the Malfliet-Tjon potential used as a challenging test case with a strong repulsive core.","marker":"[48]"},{"why":"Provides the Daejeon-16 interaction matrix elements used as the soft, modern interaction test in the paper.","marker":"[50]"}],"fun_headline_variants":["J-matrix overlap relation yields phase shifts directly","Phase shifts from overlap: a J-matrix shortcut","New overlap formula for phase shifts in J-matrix scattering","Scattering phase shifts via overlap: J-matrix trick","Overlap integrals give phase shifts in J-matrix method"],"cache_read_input_tokens":16128,"weakest_assumption_plain":"The formula is exact only under the assumption that the potential matrix elements vanish beyond a cutoff index $N_{\\rm pot}$, and the paper concedes there is no theorem guaranteeing this for realistic potentials; if that truncation is inaccurate, Eq. (36) is approximate rather than exact.","fun_headline_variants_meta":{"raw":{"variants":["J-matrix overlap relation yields phase shifts directly","Phase shifts from overlap: a J-matrix shortcut","New overlap formula for phase shifts in J-matrix scattering","Scattering phase shifts via overlap: J-matrix trick","Overlap integrals give phase shifts in J-matrix method"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000531,"raw_usage":{"total_tokens":2546,"prompt_tokens":920,"completion_tokens":1626,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":536,"completion_tokens_details":{"reasoning_tokens":1558}},"tokens_in":536,"tokens_out":1626,"duration_ms":10615,"temperature":1.0,"reasoning_tokens":1558,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:30:49.574778+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a potential whose harmonic-oscillator matrix elements decay slowly or oscillate (e.g., a long-range $1/r$ tail or a hard-core potential) and compute $\\tan \\delta$ from Eq. (36) at a fixed energy for increasing $N_{\\rm pot}$; if the phase shift does not converge or changes by more than the target accuracy when $N_{\\rm pot}$ is doubled, the central exactness claim for that potential fails. Alternatively, construct a model where $V_{n,n'}$ is nonzero for $n,n' > N_{\\rm pot}$ and show the omitted terms change the result.","supporting_citations":[{"cited_title":"Heller and Hashim A","cited_arxiv_id":null,"evidence_quote":"Establishes the J-matrix method: kinetic energy is tridiagonal in an $L^2$ basis, the starting framework for the derivation."},{"cited_title":"J- matrix method: Extensions to arbitrary angular momentum and to Coulomb scattering","cited_arxiv_id":null,"evidence_quote":"Provides the regular and irregular free-solution coefficients and the inhomogeneous equation that defines $g_n$ and $\\alpha_0$."},{"cited_title":"P-matrix and J-matrix approaches: Coulomb asymptotics in the harmonic oscillator representation of scattering theory","cited_arxiv_id":null,"evidence_quote":"Supplies the Casoratian and the J-matrix treatment of asymptotics used to identify the boundary term in the telescoping sum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives existing J-matrix phase-shift extraction methods that the new overlap relation is compared against or extends."},{"cited_title":"Kievsky, M","cited_arxiv_id":null,"evidence_quote":"Provides the coordinate-space integral relations that the paper generalizes to discrete basis representations."},{"cited_title":"One nucleon overlap integrals for light nuclei","cited_arxiv_id":null,"evidence_quote":"Shows prior use of overlap integrals for nuclear reactions, motivating the discrete analogue."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental nucleon-nucleon phase shifts used as the benchmark for the $^3S_1$ demonstration."},{"cited_title":"Three-nucleon calculations with realistic forces","cited_arxiv_id":null,"evidence_quote":"Defines the Malfliet-Tjon potential used as a challenging test case with a strong repulsive core."},{"cited_title":"N3lo nn interaction adjusted to light nuclei in ab exitu approach","cited_arxiv_id":null,"evidence_quote":"Provides the Daejeon-16 interaction matrix elements used as the soft, modern interaction test in the paper."}],"review_version":1}