{"id":"d2b65785-e53f-441f-8c18-23a8eb2284b3","arxiv_id":"1908.03178","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A quantum-mechanical sum rule for femtoscopic correlation functions is made ultraviolet-finite by considering the difference of two correlation functions with the same asymptotic behavior, and is validated against exact Coulomb results and applied to neutron-proton correlations.","lead":"This paper repairs the 1995 sum rule for femtoscopic correlation functions, which had an ultraviolet divergent momentum integral, by subtracting a regulator correlation function. The improved sum rule is exact and is demonstrated on Coulomb systems and used to test the Lednicky-Lyuboshits neutron-proton model.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The regulator identity Eq. (17) is not justified for exact Coulomb wave functions: at r=0, |φ~_q(0)|^2-1 is the Gamow factor minus 1, which decays only as 1/q, so the momentum integral diverges and Eq. (18) is not derived without an explicit distributional regularization.","rationale":"The paper has a coherent and useful idea: subtract a regulator correlation function whose own momentum integral vanishes, thereby obtaining a convergent sum rule. The derivation is clean for short-range or approximate wave functions, and the numerical verification for repulsive Coulomb and for neutron-proton with r0 > 2 fm is genuine supporting evidence. The reader correctly identified Eq. (17) as the weakest assumption. My stress-test agrees but sharpens the point: Eq. (17) is not merely questionable for the approximate asymptotic wave function (33); it is already false as an ordinary pointwise identity for the exact Coulomb regulator, because at r=0 the integrand behaves as 1/q and the momentum integral diverges. This is not a pedantic objection, because the same long-range behavior is what forces the unphysical subtraction in the attractive Coulomb case, Eqs. (26)-(28). The paper's own sentence that 'the subtraction method applied here suggests that a regulator may not be a physical correlation function' is an admission that the general claim 'the ultraviolet divergence is canceled out' does not always hold. The central claim should therefore be stated conditionally: the improved sum rule is exact only when the regulator is chosen so that the relevant q-integral converges after the indicated cancellations, or when a distributional treatment is supplied. Because the numerical evidence supports the final relations in the tested cases, a conditional verdict rather than rejection is appropriate. The concrete test of the cutoff-regulated equality would settle whether the gap is merely formal or indicates a real limitation of the method.","tokens_in":9164,"tokens_out":14915,"duration_ms":162293,"concrete_test":"Compute the cutoff-regulated version of the identity used in Eq. (16) for exact repulsive Coulomb wave functions with a Gaussian source, e.g. r0=3 fm: evaluate I(Q) = ∫_{q<Q} d^3q [R(q)-R~(q)] and J(Q) = ∫ d^3r D_r(r) ∫_{q<Q} d^3q/(2π)^3 [|φ_q(r)|^2 - |φ~_q(r)|^2], using the same spherical cutoff Q for both. If I(Q) and J(Q) do not approach the same finite limit as Q→∞, or if either depends on the order of limits, then Eq. (18) is not a consequence of the stated assumptions and an explicit distributional regularization is required.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The pivotal step is Eq. (17): the regulator correlation function is dropped because ∫ d^3q/(2π)^3 (|φ~_q(r)|^2 - 1) is asserted to vanish for a nonidentical pair without bound states. As an ordinary pointwise identity, this fails for the exact repulsive Coulomb regulator used in Sec. V.A. For a repulsive Coulomb pair, the relative wave function at r=0 has modulus squared equal to the Gamow factor G_+(q)=x/(1-e^{-x}) with x=2π/(a_B q), and G_+(q)-1 ~ π/(a_B q) for large q. The integral 4π∫_0^Q q^2 dq (G_+(q)-1) therefore diverges linearly as Q→∞, so Eq. (17) cannot hold as a standard integral at r=0. The attractive case shows the same pathology: Eq. (27) gives q^2(G_+(q)+G_-(q)-2) → 2π^2/(3a_B^2), so the sum-rule integral (21) is linearly divergent even after using a regulator, and the paper must subtract a nonphysical 2π^2/(3a_B^2 q^2) term to obtain Eq. (28). The derivation of Eq. (18) requires an interchange of q and r integrations that is not valid for these long-range wave functions without a carefully specified regularization; the paper does not supply one. The numerical agreement in Figs. 2-4 shows that the final relations can be made to work, but it does not establish Eq. (17) as a legitimate identity, and the abstract's claim that the improved sum rule works for exact Coulomb correlation functions is too broad given the subtraction needed in the attractive case.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript revisits the femtoscopic sum rule of Mrówczyński (1995), which equates the momentum-space integral of a two-particle correlation function to a source-density term and bound-state formation rates via the quantum-mechanical completeness relation. The authors identify that for interacting pairs the original momentum integral is ultraviolet divergent, and propose an improved sum rule for a difference (or sum) of two correlation functions with the same large-momentum asymptotics, Eq. (18): ∫ d³q (R(q) − R̃(q)) = ±π³ D_r(0) − Σ_α A_α. The regulator R̃ is chosen to describe a pair of nonidentical particles without bound states, so that its own momentum integral is supposed to vanish. The rule is tested on exact Coulomb correlation functions for same-sign pions (repulsive case) and opposite-sign pions (attractive case), and is then used to assess the Lednicky–Lyuboshits neutron-proton correlation function, leading to the conclusion that the asymptotic np wave function is reliable for source radii r0 ≳ 2 fm.","tokens_in":9525,"tokens_out":6967,"duration_ms":75233,"significance":"If established, Eq. (18) provides a model-independent, parameter-free integral constraint that any femtoscopic correlation function must satisfy, and it would be a useful benchmark for approximate calculations; the paper's use of exact Coulomb wave functions and the neutron-proton example convincingly illustrate the intended application. The paper also deserves credit for clearly exposing the ultraviolet problem of the original sum rule and for testing the proposed rule against independently computed exact correlation functions rather than fitting parameters. However, the derivation contains a distributional gap at Eq. (17), and the attractive-Coulomb case is explicitly divergent even after the regulator, requiring an extra subtraction that is not derived from the completeness relation. These issues affect the central claim, so the result is promising but not yet fully established as stated.","major_comments":[{"comment":"The assertion ∫ d³q/(2π)³ (|φ̃_q(r)|² − 1) = 0 for the regulator is not an ordinary identity for exact Coulomb wave functions. For the repulsive Coulomb regulator used in Sec. V.A, |φ̃_q(0)|² equals the Gamow factor G_+(q), and G_+(q) − 1 ∼ π/(a_B q) at large q, so the integral 4π ∫ q² dq (G_+(q) − 1) diverges linearly. Eq. (17) can at best hold as a distributional identity away from r = 0, and the pointwise substitution into Eq. (16), together with the interchange of q- and r-integrations, requires a regularization that the paper does not specify.","section":"Section IV, Eq. (17)"},{"comment":"The attractive-Coulomb sum rule does not follow from the completeness-based derivation as written. Even after applying the regulator, Eq. (27) shows q²(G_+(q) + G_−(q) − 2) → 2π²/(3a_B²), so the integral in Eq. (24) is linearly divergent; the paper subtracts the nonphysical term 2π²/(3a_B² q²) to obtain Eq. (28). This subtraction is not derived from Eq. (17) or Eq. (18), and the finite-source case is explicitly left unresolved in the text ('We have not been able to convincingly show...'). The abstract's claim that the improved sum rule works for exact Coulomb correlation functions is therefore too broad without a clearly stated regularization prescription.","section":"Section V.B, Eqs. (27)-(28)"},{"comment":"The derivation of the main result relies on changing the order of q- and r-integrations after assuming that both integrals converge, as the paper itself notes. For the Coulomb regulators used in the tests, the q-integrals are not absolutely convergent, so Eq. (18) is not rigorously established by the argument given. The numerical saturation seen in Figs. 2-4 is encouraging evidence that the final relation can be made to work, but it does not supply the missing justification. I recommend reformulating the sum rule distributionally, or at least stating explicitly the cutoff or symmetric-integration prescription under which Eq. (17) and Eq. (18) hold.","section":"Section IV, Eqs. (16)-(18)"}],"minor_comments":[{"comment":"The regulator correlation function is written as R̃ in the text but sometimes appears as ~R in equations; the notation should be defined once after Eq. (16) and used consistently throughout.","section":"General notation"},{"comment":"The passage from the finite-source formation rate A_{nlm} defined by Eq. (13) to the point-like-source expression A_{nlm} = (2π)³|φ_{nlm}(0)|² in Eq. (23) should be displayed explicitly, including the use of D_r(r) = δ³(r).","section":"Section V.B, Eq. (21)"},{"comment":"The spelling 'Gamov factor' should be 'Gamow factor' in Eq. (22) and the surrounding text.","section":"Section V.B, Eq. (22)"},{"comment":"The numerical procedure for extracting the deuteron formation rate from the integrated correlation function should be described briefly in the text, since the figure compares four independent computations and no error estimates are given.","section":"Fig. 8"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a real problem and the proposed sum rule is likely correct in a suitably regularized sense, but the mathematical derivation needs tightening and the attractive-Coulomb case requires an honest statement of the extra subtraction. If the authors can supply a distributional formulation and revise the abstract's claim accordingly, the paper would be a useful contribution to the femtoscopy literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Radek,\n\nThe thing to know: this paper actually fixes the UV divergence that killed the 1995 sum rule, and the fix is simple enough to be useful. The improved rule, Eq. (18), says the momentum integral of R - R~ equals the bound-state term, provided the regulator pair has no bound states. That is a genuine step forward, and the repulsive-Coulomb and neutron-proton tests give it real support. I would credit the authors for being candid about the original flaw.\n\nWhat is new: the regulator idea itself, taking a difference of two correlation functions with the same asymptotic behavior instead of regulating one function ad hoc. It is not paradigm-shifting, but it is a clean contribution within an established program, and it gives femtoscopists a practical consistency check for approximate models like Lednicky-Lyuboshits. The neutron-proton application is the most valuable part: the sum rule gives a quantitative criterion for when the asymptotic wave function (33) breaks down, and the r0 > 2 fm conclusion looks sensible.\n\nNow the soft spots. The stress-test note is right: Eq. (17), the claim that the regulator's momentum integral vanishes identically, is not true as an ordinary integral for the exact repulsive Coulomb wave function. At r=0, |phi~_q(0)|^2 - 1 is the Gamow factor minus 1, which goes like 1/q, so the integral diverges linearly. The derivation of Eq. (18) needs a distributional or subtracted regularization that the paper never spells out. More concretely, the attractive Coulomb case in Sec. V.B still diverges after applying the regulator, and the authors subtract a nonphysical 2pi^2/(3a_B^2 q^2) term to get Eq. (28). That is fine as a mathematical trick, but the abstract says the improved sum rule \"works well for the exact Coulomb correlation functions\" without mentioning the subtraction. That is an overstatement and should be fixed.\n\nThe numerical agreement in Figs. 2-4 shows the final relations can be made to work, but it does not rescue Eq. (17) as a stated identity. The paper needs an explicit regularization scheme or a clear statement that the regulator is a limiting device rather than a physical correlation function. The authors almost say this at the end of Sec. V.B but don't carry it through to the abstract or the derivation.\n\nWho is this for: heavy-ion femtoscopy practitioners who compute correlation functions from models and need a quick sanity check. It deserves a serious referee, though the referee should push on the regulator justification and the abstract's scope. I would accept it with revision.\n\nCandidly,\n[Your name]","headline":"Useful, honest femtoscopy sum-rule paper with a real UV-divergence fix, but the attractive-Coulomb case leans on an unadvertised subtraction and the regulator identity is not pointwise justified for exact Coulomb wave functions.","tokens_in":10031,"tokens_out":942,"would_cite":true,"duration_ms":11727,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.Gz"],"model":"deepseek-v4-flash","headline":"The paper proves an exact integral identity for the difference of two femtoscopic correlation functions, repairing a 1995 sum rule whose momentum integral diverged for interacting pairs.","keywords":["femtoscopy","correlation function","sum rule","quantum completeness","ultraviolet divergence","neutron-proton correlation","Coulomb interaction","deuteron formation"],"falsifier":"Take an exactly solvable two-particle system with a known shallow bound state (e.g., a delta-shell or separable potential) and compute the correlation functions $R$ and $\\tilde{R}$ from the exact wave functions for several source radii; evaluate the left side of Eq. (18) with a momentum cutoff and compare with $\\pm\\pi^3D_r(0)-A$. If the saturated value differs from the right side by more than the numerical error, the sum rule fails. For the neutron-proton case, the same check can be made experimentally: measure the triplet and singlet correlation functions and ask whether $4\\pi\\int dq\\, q^2(R_t-R_s)$ equals $-A_D$ over a range of source sizes.","tokens_in":8961,"feed_emoji":"⚛️","tokens_out":5139,"duration_ms":52517,"temperature":0.7,"pith_summary":"The paper claims to repair a sum rule for femtoscopic correlation functions that was proposed in 1995 but turned out to be mathematically ill-defined: the momentum integral of a single correlation function diverges in exactly the cases of physical interest. The repair is to integrate not one correlation function but the difference (or sum) of two correlation functions chosen to have the same large-momentum behavior, so the divergent part cancels. The resulting identity, Eq. (18), relates the integrated difference to the two-particle source density at zero separation and to bound-state formation rates, and holds exactly when the second function describes a nonidentical pair with no bound states. Because the identity is exact, a model calculation of correlation functions can be tested by checking whether its integrated difference approaches the predicted constant. The paper demonstrates the rule on exact Coulomb correlation functions and uses it to determine when the standard neutron-proton correlation model breaks down.","feed_headline":"Difference of two correlation functions obeys an exact sum rule","feed_subtitle":"Repairing a 1995 sum rule, the paper gets finite momentum integrals and a test for model calculations.","key_machinery":"The load-bearing object is the completeness (closure) relation for the two-particle wave functions, Eq. (9), used in the form $\\int \\frac{d^3q}{(2\\pi)^3}(|\\varphi_q(r)|^2-1) = \\pm \\delta^{(3)}(2r) - \\sum_\\alpha |\\varphi_\\alpha(r)|^2$. For the regulator pair—nonidentical particles with no bound states—this integral is zero (Eq. (17)), which is exactly the piece that makes the original single-function integral divergent. The improved sum rule therefore works by cancellation: $R$ and $\\tilde{R}$ share the same large-$q$ tail, and their difference falls fast enough for the momentum integral to exist.","core_discovery":"The central claim is Eq. (18): for two correlation functions $R(q)$ and $\\tilde{R}(q)$ built from the same source, with $\\tilde{R}$ a regulator whose pair has no bound states and whose continuum completeness is that of free nonidentical particles, $\\int d^3 q\\,(R(q)-\\tilde{R}(q)) = \\pm \\pi^3 D_r(0) - \\sum_\\alpha A_\\alpha$, where $D_r$ is the relative source distribution and $A_\\alpha$ are bound-state formation rates. The proof runs through the quantum-mechanical closure relation: the integral over $q$ of $(|\\varphi_q(r)|^2-1)$ equals a delta term at $2r$ minus the bound-state densities, and for the nonidentical regulator that integral vanishes identically. Subtracting the two correlation functions cancels the common ultraviolet tail, so the momentum integral converges and the $r$- and $q$-integrations can be interchanged. The result restores the physical content of the original 1995 sum rule—integrated correlation measures source size and bound-state production—without its divergence.","pith_inferences":["The same regulator strategy could test models with open inelastic channels (e.g., $K^-p$ or $p\\bar{p}$), provided the regulator is chosen to carry the same channel content; the completeness relation would then include the inelastic states explicitly.","Because the sum rule holds for any source function, it could be used as a data-driven consistency check in experimental femtoscopy, comparing extracted sources and interaction parameters without committing to a specific model.","In the attractive case the extra subtraction in Eq. (28) shows that a regulator need not be a physical correlation function; one could construct purely formal regulators that cancel even the strongest divergent tails, at the price of adding known subtraction constants.","A natural extension is to apply the rule to correlation functions computed with realistic, energy-dependent and coupled-channel potentials; the failure of the rule would localize the regime where the approximate wave function or the completeness assumption breaks down."],"forward_implications":["Any femtoscopic correlation function computed in an approximate model must satisfy Eq. (18); checking the integrated difference against $\\pm\\pi^3D_r(0)-\\sum_\\alpha A_\\alpha$ gives a quantitative test of the model's accuracy and range of validity.","For exact Coulomb repulsion between identical pions, the integrated difference of symmetrized and unsymmetrized correlation functions reproduces $\\pi^3D_r(0)$ once $q_{\\rm max}r_0 \\gtrsim 1.5$, so the rule is verified for an exactly solvable case.","For the Coulomb-attractive case the same subtraction strategy works, but an extra analytic subtraction is needed (Eq. (28)) because the sum of Gamow factors still has a linearly divergent integrand; the rule then matches $-8\\pi^2\\zeta(3)/a_B^3$.","For neutron-proton pairs, the sum rule supplies the deuteron formation rate $A_D$ from the triplet-singlet correlation difference; the model is reliable only for source radius $r_0 > 2$ fm, and the commonly used correction factor (35) does not extend it much below that."],"supporting_citations":[{"why":"Original 1995 sum rule whose ultraviolet divergence this paper repairs.","marker":"[4]"},{"why":"Documented the divergence of the single-function integral and supplied the Coulomb correlation functions used in the exact tests.","marker":"[5]"},{"why":"Lednicky-Lyuboshits asymptotic wave function and correction factor used for the neutron-proton correlation functions.","marker":"[7]"},{"why":"Standard reference for the Coulomb wave functions and atomic bound-state densities used in the exact tests.","marker":"[8]"},{"why":"Provides the Hulthén wave function parameters used to compute the deuteron formation rate.","marker":"[12]"},{"why":"Supplies the neutron-proton scattering lengths and effective ranges used in the amplitude (34).","marker":"[13]"},{"why":"Experimental example where correlation functions determine scattering lengths, motivating the need for model tests.","marker":"[3]"}],"fun_headline_variants":["Subtraction trick fixes divergent femtoscopic sum rule","Repaired sum rule: difference of correlation functions converges","Finite sum rule from subtracting two correlation functions","Exact sum rule for difference of correlations tests models","1995 sum rule made convergent by subtracting correlations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the regulator pair's scattering states are complete in exactly the same way as free waves, so the momentum integral of $|\\tilde{\\varphi}_q(r)|^2-1$ vanishes; for neutron-proton this reduces to trusting an approximate asymptotic wave function that the sum rule itself shows is valid only for sources larger than about 2 fm.","fun_headline_variants_meta":{"raw":{"variants":["Subtraction trick fixes divergent femtoscopic sum rule","Repaired sum rule: difference of correlation functions converges","Finite sum rule from subtracting two correlation functions","Exact sum rule for difference of correlations tests models","1995 sum rule made convergent by subtracting correlations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000576,"raw_usage":{"total_tokens":2687,"prompt_tokens":882,"completion_tokens":1805,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":498,"completion_tokens_details":{"reasoning_tokens":1732}},"tokens_in":498,"tokens_out":1805,"duration_ms":13609,"temperature":1.0,"reasoning_tokens":1732,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:21:11.256997+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take an exactly solvable two-particle system with a known shallow bound state (e.g., a delta-shell or separable potential) and compute the correlation functions $R$ and $\\tilde{R}$ from the exact wave functions for several source radii; evaluate the left side of Eq. (18) with a momentum cutoff and compare with $\\pm\\pi^3D_r(0)-A$. If the saturated value differs from the right side by more than the numerical error, the sum rule fails. For the neutron-proton case, the same check can be made experimentally: measure the triplet and singlet correlation functions and ask whether $4\\pi\\int dq\\, q^2(R_t-R_s)$ equals $-A_D$ over a range of source sizes.","supporting_citations":[{"cited_title":"Mr´ owczy´ nski, Phys","cited_arxiv_id":null,"evidence_quote":"Original 1995 sum rule whose ultraviolet divergence this paper repairs."},{"cited_title":"Maj and St","cited_arxiv_id":null,"evidence_quote":"Documented the divergence of the single-function integral and supplied the Coulomb correlation functions used in the exact tests."},{"cited_title":"Lednicky and V","cited_arxiv_id":null,"evidence_quote":"Lednicky-Lyuboshits asymptotic wave function and correction factor used for the neutron-proton correlation functions."},{"cited_title":"Schiff, Quantum Mechanics (McGraw-Hill, New York, 1968)","cited_arxiv_id":null,"evidence_quote":"Standard reference for the Coulomb wave functions and atomic bound-state densities used in the exact tests."},{"cited_title":"Hodgson, Nuclear Reactions and Nuclear Structure (Clarendon Press, Oxford, 1971)","cited_arxiv_id":null,"evidence_quote":"Provides the Hulthén wave function parameters used to compute the deuteron formation rate."},{"cited_title":"McCarthy, Introduction to Nuclear Theory (John Wiley & Sons, New York, 1968)","cited_arxiv_id":null,"evidence_quote":"Supplies the neutron-proton scattering lengths and effective ranges used in the amplitude (34)."}],"review_version":1}