{"id":"af659dfa-ec18-45c0-84f5-a3d411cdf528","arxiv_id":"2506.12531","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive exact renormalization and binding-energy relations for the Efimov effect in three identical bosons from the Skorniakov-Ter-Martirosian equation.","lead":"This paper derives an exact relation between the three-body coupling, the ultraviolet cutoff, and the Efimov three-body parameter for three identical bosons, plus an exact formula for the binding energies of Efimov states. If correct, it replaces previously approximate or purely numerical constants in Efimov physics with closed-form integral expressions.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exact renormalization relation (7) follows from the proof only if ξ(t) and ξ′(t) are bounded on t>0; this is assumed in SM §II.C, never proven, and is the step that forces f to be linear via Liouville. If the bound fails, the derived H(φ) does not follow.","rationale":"The reader identified the same weakest assumption: boundedness of ξ and ξ′ in §II.C of the supplementary material. I agree, and I would elevate it from a 'minor caveat' to the hinge of the proof. The entire Wiener-Hopf construction, including the definition of f, the application of Liouville's theorem, and the resulting linear form f=κ−iω, depends on the decay estimate (45). Without that estimate, the analyticity arguments only show that f is entire; they do not force f to be linear, so the exact functional dependence of H on φ is not established. The numerical agreement (claimed discrepancies below 1e−7 for the renormalization relation and 1e−6 for the energies) makes it very plausible that the boundedness assumption is satisfied on the physical branch, so this is a proof-completeness concern rather than an indication of numerical error. I recommend conditional acceptance: the authors should supply a direct proof of boundedness of ξ and ξ′ from Eq. (11), or verify it for the explicit solution (51) as described in the concrete test. A related but minor issue is that the main text states F is analytic in a strip with a<2, whereas the nearest singularities of F are at ±i, so the correct bound is a<1; the proof only needs existence of some a>0, so this does not affect the argument but should be corrected.","tokens_in":16447,"tokens_out":30584,"duration_ms":377142,"concrete_test":"Take the explicit solution (51) with κ related to H by Eq. (20) (equivalently Eq. (7)), compute F+ via the contour integral (17), and numerically Fourier-invert ξtilde to obtain ξ(t) for t∈(0,∞) for a dense set of H values spanning the full range shown in Fig. 1, including H near the pole where 1−s0 tan(φ−φ0) vanishes. Check that sup_{t>0}(|ξ(t)|+|ξ′(t)|) is finite and that Eq. (45) holds with a finite C1. If the sup is finite for all such H, the a priori assumption is justified and the proof can be made rigorous; if it diverges for any H, the Liouville step and therefore Eq. (7) must be revisited.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central step is the proof that f(ω), defined in Eq. (18), is linear. This relies on Liouville's theorem, and the polynomial-growth bound is obtained by using, first, |ξtilde(ω)| ≤ C/(1+|ω|) for Imω ≥ ε. That bound is derived in Supplementary §II.C from the explicit assumption that ξ(t) and ξ′(t) are bounded on t>0. The assumption is not derived from the integral equation (11), and it is not verified a posteriori from the explicit solution (51). If the true solution had a non-integrable singularity in ξ′ near t=0, or an unbounded derivative at large t, the decay estimate in Eq. (45) would fail, the Liouville argument would be invalid, and the exact form of H(φ), Eqs. (7)-(8), would not follow from the presented argument. Because this is the only step that fixes the functional form of the renormalization relation, it is load-bearing. The numerical validation gives strong evidence that the final formula is nevertheless correct, but an exact proof requires the regularity of ξ to be established rather than assumed; otherwise the claim of exactness is not fully supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the Skorniakov-Ter-Martirosian integral equation for three identical bosons at two-body resonance with a finite momentum cutoff. The authors claim exact closed-form results: the renormalization relation H = α(1 + s0 tan(φ − φ0))/(1 − s0 tan(φ − φ0)) with α ≈ 0.87866 and φ0 ≈ 0.05281, and the asymptotic Efimov binding-energy relation s0 ln(Λ*/κ*) = −s0 ln(√3 e^{−π/(2s0)}) ≈ 1.01807 mod π. The derivation combines a Wiener-Hopf-style analytic factorization of the zero-energy kernel (Regime I) with an exact sine solution of the finite-energy equation after a change of variables (Regime II). High-precision numerical solution of the discretized STM equation is used to validate the formulas, with reported discrepancies below 10^{-7} for the renormalization relation and about 10^{-6} for higher bound-state energies.","tokens_in":16682,"tokens_out":10776,"duration_ms":129554,"significance":"If correct, these results settle two long-standing questions in universal three-body physics: the exact cutoff dependence of the three-body coupling and the universal binding-energy spectrum. The constants α and φ0 are computed from parameter-free integrals over the known kernel rather than fitted to the target results, and the numerical validation is independent and high-precision. The analytic method, a generalization of the Wiener-Hopf technique, is likely transferable to other Efimov-type problems. The main caveat is that the proof of the renormalization relation assumes a regularity property that is stated but not established, which makes the exactness claim conditional in its present form.","major_comments":[{"comment":"The proof that f(ω) is linear, which fixes the functional form of the renormalization relation H(φ) and hence Eqs. (7)-(8), relies on the estimate |ξ~(ω)| ≤ C1/(1+|ω|) for Im ω ≥ ε. This estimate is derived from the explicit assumption that ξ(t) and ξ′(t) are bounded on t>0, stated just before Eq. (45). The assumption is not proved from the integral equation (11) and is not verified from the explicit solution (51); it is therefore load-bearing. If the true solution had an unbounded derivative near t=0 or at large t, the Liouville argument would fail and the claimed exact form of H(φ) would not follow from the presented argument. Please prove the boundedness from the equation or otherwise justify it, or state the main theorem as conditional on this regularity assumption.","section":"Supplementary §II.C (Eqs. (44)-(50))"},{"comment":"The paper cites the previous numerical result s0 ln(Λ*/κ*) ≈ 0.971 mod π and then claims the exact value 1.01807. The two numbers differ by about 0.047, which is not attributable to the reported numerical uncertainty. Since Eq. (11) is a central result, the manuscript should explain the discrepancy, for example by showing that the earlier value used a different phase convention or by demonstrating that the earlier numerical calculation was less accurate. Without such an explanation, readers cannot assess whether the new value supersedes or conflicts with the earlier literature.","section":"Eq. (6) versus Eq. (11) and Fig. 2"}],"minor_comments":[{"comment":"The statement that F(ω) is analytic and zero-free in the strip |Im ω| ≤ a for any a < 2 is false for a > 1, because G(ω) has poles at ω = ±i. Since the proof only needs the existence of some a > 0 (for example, any a < 1), the statement should be corrected.","section":"Main text, after Eq. (16)"},{"comment":"Please report the numerical method used to evaluate the integrals in Eq. (8) for α and φ0, and provide the values with more digits or a reproducibility statement. This would strengthen the parameter-free claim and allow readers to verify the quoted constants.","section":"Main text, Numerical Validation"},{"comment":"In Eq. (24), the approximation τ ≈ ln(k/√(|E|/3)) should explicitly state the validity condition k ≫ √|E|; the condition currently appears only in the surrounding text.","section":"Eq. (24)"}],"recommendation":"major_revision","confidential_remarks":"The regularity gap in Supplementary §II.C is the main substantive obstacle to acceptance. If the authors can supply a rigorous proof of the boundedness of ξ and ξ′ (or a non-circular verification using the explicit solution), and if they clarify the discrepancy with the literature value 0.971, the paper would be a strong contribution. The numerical evidence is already very convincing."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it likely solves a long-standing problem: the exact renormalization relation for three identical bosons at unitarity, giving H = α(1 + s0 tan(φ − φ0))/(1 − s0 tan(φ − φ0)) with α ≈ 0.87866 and φ0 ≈ 0.05281, and a companion exact formula for the Efimov binding energies. Second, the derivation has one genuinely load-bearing assumption that the authors state but do not prove, so the claim of exactness is slightly stronger than the proof currently supports.\n\nWhat is genuinely new: the Wiener-Hopf treatment of the finite-cutoff STM equation, the explicit factorization of the kernel, and the resulting closed-form integrals for α and φ0. The binding-energy result comes from an elegant change of variables that turns the nonzero-energy equation into one solved exactly by sin(s0 τ). The previous literature has only approximate analytic (α = 1) or numerical (α ≈ 0.879, 0.971) values, so the improvement is real and well-motivated.\n\nWhat the paper does well: the numerical validation is strong — discrepancies below 10^-7 for the renormalization relation and about 10^-6 for the binding energies, with no parameter fitting to the target results. The physical reasoning is clear, and the supplementary material is detailed.\n\nThe soft spot: the proof that f(ω) is linear, which fixes the functional form of H(φ), relies on the bound |ξ(ω)| ≤ C/(1 + |ω|), obtained from the assumption that ξ(t) and ξ′(t) are bounded on t > 0 (Supplementary §II.C, Eq. 45). This boundedness is not derived from the integral equation, and the authors do not verify it a posteriori from the explicit solution (51). If the true solution had an unbounded derivative near t = 0, the Liouville argument would fail, and the exact form of H(φ) would not follow. This is a real gap, but it is the only serious one, and the numerical evidence strongly suggests the final formula is nonetheless correct. The authors should either prove the boundedness or present it as a condition under which the result is established.\n\nMinor notes: the integrals defining α and φ0 are evaluated numerically, but they are exact analytic expressions; that is not a flaw. The paper is concise to the point of being terse, but the supplement fills in the steps.\n\nWho this is for: anyone working on Efimov physics, few-body universality, or three-body effective field theory, especially in cold atoms or nuclear physics. It deserves a serious referee — the gap in the proof should be raised in review, but the result is important enough that the paper should not be desk-rejected.","headline":"Strong candidate for the exact Efimov renormalization relation, with a real but possibly fixable gap in the proof of the key Liouville step.","tokens_in":17219,"tokens_out":1589,"would_cite":true,"duration_ms":22030,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives the exact renormalization relation and shallow Efimov binding energies for three identical bosons at unitarity, fixing previously approximate constants.","keywords":["Efimov effect","three-body parameter","renormalization","Skorniakov-Ter-Martirosian equation","unitary Bose gases","Wiener-Hopf method","Efimov binding energies","universal few-body physics"],"falsifier":"Directly compute the zero-energy solution of Eq. (12) at high precision near $t=0$ and form its Fourier transform $\\tilde\\xi(\\omega)$: if $|\\tilde\\xi(\\omega)|$ grows faster than $C/(1+|\\omega|)$ for $\\mathrm{Im}\\,\\omega \\ge \\epsilon$, the Liouville step collapses and the exact relation (7) would need modification. A physics-level cross-check is the ratio of successive shallow Efimov binding energies, which the paper predicts to be governed by $s_0\\ln(\\Lambda_*/\\kappa_*) = 1.01807 \\pmod{\\pi}$ rather than the previously reported $0.971$.","tokens_in":16243,"feed_emoji":"⚛️","tokens_out":10392,"duration_ms":105835,"temperature":0.7,"pith_summary":"This paper claims to settle, in the low-energy limit, the two open constants of three-body physics for identical bosons at infinite scattering length: the renormalization relation connecting the three-body coupling to the cutoff, and the precise binding energies of the infinite tower of Efimov states. The authors solve the Skorniakov–Ter-Martirosian equation exactly at zero energy with a generalization of the Wiener-Hopf method, obtaining $H = \\alpha\\,(1 + s_0\\tan(\\varphi-\\varphi_0))/(1 - s_0\\tan(\\varphi-\\varphi_0))$ with $\\alpha \\approx 0.87866$ and $\\varphi_0 \\approx 0.05281$, and they derive the shallow-state binding phase $s_0\\ln(\\Lambda_*/\\kappa_*) \\approx 1.01807 \\pmod{\\pi}$. Both formulas are checked numerically, with discrepancies below $10^{-7}$ for the renormalization relation and about $10^{-6}$ for the third bound state. If correct, these results replace previously approximate or purely numerical constants with exact analytic ones, giving a complete low-energy characterization of three identical bosons at unitarity.","feed_headline":"Three bosons at unitarity get exact binding-energy law","feed_subtitle":"A Wiener-Hopf analysis fixes the three-body constant to 0.87866 and the Efimov phase to 1.01807.","key_machinery":"The load-bearing object is the Skorniakov–Ter-Martirosian integral equation with a finite cutoff, analysed by a Wiener-Hopf-style extension of the wavefunction $\\xi(t)$ to the whole real line with an auxiliary error term. The Fourier-domain equation is factorized through $F(\\omega) = F_+(\\omega)F_-(\\omega)$, and the key function $f(\\omega) = (s_0^2-\\omega^2)(1-i\\omega)F_+(\\omega)\\tilde\\xi(\\omega)$ is shown to be analytic in the whole complex plane. A bound $|\\tilde\\xi(\\omega)| \\le C/(1+|\\omega|)$ on the upper half-plane makes $|f(\\omega)|/(1+|\\omega|)$ bounded, so Liouville's theorem forces $f(\\omega) = \\kappa - i\\omega$ to be linear; this linearity fixes $\\kappa = (\\alpha-H)/(\\alpha+H)$ and hence the renormalization relation (7). For the binding energies, a second exact factorization carries the argument: in the variable $\\tau = \\mathrm{arcsinh}(k/\\sqrt{-4E/3})$, the kernel splits as $K(\\tau,\\sigma) = G(\\tau-\\sigma) - G(\\tau+\\sigma)$, so that $\\sin(\\omega\\tau)$ is an exact eigenfunction, giving the phase $\\eta = \\pi/2$ without approximation.","core_discovery":"At two-body resonance, with a finite cutoff $\\Lambda$, the three-body coupling $h = H/\\Lambda^2$ is related to the three-body parameter $\\Lambda_*$ by the exact renormalization relation $$H = \\$\\alpha$\\,\\frac{1 + s_0\\tan(\\varphi-\\varphi_0)}{1 - s_0\\tan(\\varphi-\\varphi_0)},$$ where $\\varphi = s_0\\ln(\\Lambda_*/\\Lambda)$, with $\\alpha$ and $\\varphi_0$ given by the integrals in Eq. (8), numerically $\\alpha \\approx 0.87866$ and $\\varphi_0 \\approx 0.05281$. The binding momenta $\\kappa_*$ of the shallow Efimov states satisfy $s_0\\ln(\\Lambda_*/\\kappa_*) = -s_0\\ln\\!\\big(\\sqrt{3}\\,e^{-\\pi/(2s_0)}\\big) \\approx 1.01807 \\pmod{\\pi}$, which fixes the ratio of consecutive binding energies in the limit $n\\to\\infty$. The paper also gives the shallow bound-state wavefunction for all $k \\in (0,\\Lambda)$, taking the sine form $\\sin\\!\\big(s_0\\,\\mathrm{arcsinh}(k/\\sqrt{-4E/3})\\big)$ at small $k$. Together these results determine the entire low-energy behavior of the system once $\\Lambda_*$ is known.","pith_inferences":["Because the exact phase $1.01807$ differs from the widely used $0.971$, experiments that extract $\\Lambda_*$ from a measured Efimov energy with the old phase systematically shift $\\Lambda_*$ by a factor of roughly $e^{(1.01807-0.971)/s_0} \\approx 1.05$; re-analyzing published Efimov spectra with Eq. (11) is an immediate, low-cost test of the claim.","The mechanism that makes $f(\\omega)$ linear is generic: for any kernel whose Fourier transform has the same positivity, decay, and analyticity properties, the same Wiener-Hopf construction fixes the renormalization relation, so the integral formulas for $\\alpha$ and $\\varphi_0$ are a template for other three-body problems, not a one-off calculation.","A natural next step is to include corrections in $1/(\\Lambda a_s)$ at large but finite scattering length; the exact zero-$\\Lambda_*$ relation should control the leading-order shift of the Efimov spectrum, which cold-atom experiments varying $a_s$ could test."],"forward_implications":["The three-body parameter requires no numerical fitting: for any cutoff $\\Lambda$, the coupling $H$ that keeps $\\Lambda_*$ fixed is exactly $H = \\alpha\\,(1+s_0\\tan(\\varphi-\\varphi_0))/(1-s_0\\tan(\\varphi-\\varphi_0))$.","Consecutive Efimov binding energies in the shallow limit obey a ratio fixed by the exact phase $1.01807 \\pmod{\\pi}$, replacing the earlier approximate $0.971$.","The shallow bound-state wavefunction is now known on the whole momentum interval $0