{"id":"0639a6d8-007b-400a-8ac1-3966d09249ab","arxiv_id":"2602.12067","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a trial state with Huang–Yang-precision energy, the averaged momentum distribution of the dilute Fermi gas is rigorously shown to match Belyakov's 1961 formula, with subleading error.","lead":"This mathematics paper proves that a trial state of a dilute spin-1/2 Fermi gas, accurate enough to match the Huang–Yang ground-state energy, reproduces Belyakov's 1961 formula for the momentum distribution with explicit error bounds. It is the first rigorous check of that sixty-year-old perturbative prediction, and a step toward proving interactions do not destroy the Fermi surface.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 4.9's sketched number-operator bound is the linchpin of the error budget; an unjustified ρ^{1/2} or L factor in its proof would make the stated ρ^{5/3+1/9} error in Theorem 2.1(2) fail.","rationale":"The reader's conditional verdict already identifies Lemma 4.9 as a liability. My stress-test sharpens this: Lemma 4.9 is not merely 'sketched' but its displayed estimate contains an unexplained L^{3/2}ρ^{1/2} factor that is essential for the final error exponent. If that factor is wrong, the number-operator bound becomes too weak by a power of ρ, and the errors in Proposition 3.4 cease to be subleading to Belyakov's term. I verified the rest of the error budget at exponent level: with Lemma 4.9 as stated, the §5.5.2 terms scale like ρ^{31/18−2κ}, ρ^{7/3−κ}, ρ^{2−κ}, and ρ^{19/9−κ}, all subleading to ρ^{16/9}; and (2.6) follows from Lemma A.2. The energy import from preprints is a secondary concern because part (2) does not rely on it. Therefore the reader's CONDITIONAL verdict is appropriate; it should be UNCHANGED, with the condition being a complete proof of Lemma 4.9 and, ideally, independent verification of the energy imports.","tokens_in":49766,"tokens_out":15228,"duration_ms":146323,"concrete_test":"Complete the proof of Lemma 4.9 by deriving the displayed second inequality of §4.4 from Lemmas 4.2–4.8, explicitly tracking the L-dependence of ∫dxdy|χ<(x−y)| and all e^{±t k_F²} factors from ∥ût∥₂ and ∥v̂t∥₂. If the exact bound is not C L^{3/2}ρ^{1/2}... but, say, C L^{3/2}... or C L^{5/2}ρ^{1/2}..., recalculate the error budget in §5.5.2 and compare with ρ^{5/3+3α}; this settles whether (2.5) holds.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of the central claim (2.5) passes all error terms through Lemma 4.9, the bound ⟨N⟩ ≤ Cκ L³ ρ^{14/9−κ}. This lemma is only sketched in §4.4. The displayed second inequality in the sketch,\n\n|∂λ⟨T2;λΩ,NT2;λΩ⟩| ≤ C L^{3/2} ρ^{1/2} ||N^{1/2}T2;λΩ|| ∫ dt e^{−2tε} ||ût,↑||₂ ||v̂t,↑||₂,\n\nis not justified by the cited Lemmas 4.2–4.8 as written: Lemma 4.4 gives ∫dx ||a(v_t,x)ψ||² ≤ e^{2t k_F²}⟨N⟩, and the ∫dxdy |χ<(x−y)| factor carries volume dependence, so it is not immediate where L^{3/2}ρ^{1/2} comes from. If a complete derivation drops the ρ^{1/2} (or adds an extra L^{3/2}), the Gronwall argument yields a much weaker bound, e.g. ⟨N⟩ ≤ L³ρ^{5/9} instead of L³ρ^{14/9−κ}. Then the error term L^{-3}ρ^{7/9}⟨N⟩ appearing in §5.5.2 scales as ρ^{4/3}, which is not subleading to the Belyakov term ρ^{5/3+3α} for α<1/27. Thus Theorem 2.1(2) would not follow. The imported energy statement (2.4) from [GHNS24; GHNS25] is also unverified here, but it is not used in the derivation of (2.5), so Lemma 4.9 is the more load-bearing gap.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a trial state for the dilute spin-1/2 Fermi gas — the same state used in the Huang–Yang energy upper bound of [GHNS24] — and proves that its averaged momentum-space excitation density agrees with Belyakov's 1961 formula. The main result (Theorem 2.1) states that for the state Ψ = R T₁T₂Ω one has, for momenta near the Fermi surface, limsup_L |⟨Ψ, n^{exc}_{q,α}Ψ⟩ − n^{(Bel)}_{q,α}| ≤ C ρ^{5/3+1/9}, with the Belyakov term itself of order ρ^{5/3+3α}. The proof uses a second-order Duhamel expansion of ⟨T₁T₂Ω, n_g T₁T₂Ω⟩, estimates of the resulting double commutators with the quasi-bosonic generators B₁ and B₂, a number-operator bound, and a careful removal of the momentum cutoffs and ε-regularization. The paper also states an energy estimate for the same trial state, imported from [GHNS24; GHNS25].","tokens_in":50035,"tokens_out":31164,"duration_ms":331531,"significance":"If correct, the result converts a long-standing formal perturbative prediction — Belyakov's 1961 momentum distribution — into a rigorous statement about an explicit trial state in the low-density Fermi gas. The paper is transparent about the error budget and about the parameter window α < 1/27, and the leading term is pinned by the physical scattering length rather than by an adjustable constant. The comparison with Belyakov's formula is an external benchmark, not an internally imposed object, and the authors are explicit that simpler trial states are not expected to reproduce the Belyakov constant. The main caveat is that the proof relies on a number-operator estimate, Lemma 4.9, whose proof is only sketched and whose stated exponent is load-bearing for the main convergence claim.","major_comments":[{"comment":"Lemma 4.9 is the linchpin of the error budget: all error terms in §5.5.2 are absorbed using ⟨N⟩ ≤ C_κ L³ ρ^{14/9−κ}. The proof of the lemma is only sketched, and the displayed second inequality is not a direct consequence of the cited Lemmas 4.2–4.8. In particular, the factor ∫dxdy |χ_<(x−y)| in the first displayed line is of order L³, and Lemma 4.4 gives ∫dx ∥a(v_{t,x})ψ∥² ≤ e^{2tk_F²}⟨N⟩; the origin of the factor L^{3/2}ρ^{1/2} in the second line therefore requires a nontrivial argument that is not supplied. If the correct bound instead contained an additional L^{3/2} or a different power of ρ, the Gronwall argument would yield a weaker ⟨N⟩ bound, and the error term L^{−3}ρ^{7/9}⟨N⟩ in §5.5.2 would scale as ρ^{4/3} rather than ρ^{7/3−κ}, which would not be subleading to the claimed ρ^{5/3+1/9} error and would invalidate Theorem 2.1(2). The authors should provide the full Gronwall/commu","section":"§4.4, Lemma 4.9"},{"comment":"The energy statement (2.4) is imported from two preprints, [GHNS24] and [GHNS25]. While citing preprints is acceptable in mathematical physics, the paper should state clearly that this part of the theorem is conditional on the detailed content of those works, especially since the present paper uses a relaxed δ-range that is justified only by the same Lemma 4.9. The momentum claim (2.5) is independent of (2.4) in the proof, but the theorem as stated should be unambiguous about which parts are proved here and which are inherited.","section":"Theorem 2.1, Eq. (2.4)"}],"minor_comments":[{"comment":"The denominators appear to have a typo: inside the square, |r′| should be |r′|². The correct form appears later, e.g. in (5.12), but the two displayed Belyakov formulas are inconsistent as written.","section":"Eq. (2.3) and Eq. (3.19)"},{"comment":"Several minor typos: 'though of' should be 'thought of' (Introduction and Remark 2.6); in §5.5.2 the constant C_κ is typeset as 'C k' in several places.","section":"Throughout"},{"comment":"The proofs of Propositions 5.3 and 5.4 contain many 'similarly' and 'omitted details' passages, especially for terms I_{5;c}, I_{7;c}, and I_{9;b–d}. The pattern is clear, but a few more intermediate steps would improve verifiability.","section":"§5.3–§5.4"}],"recommendation":"major_revision","confidential_remarks":"The main obstacle to acceptance is the missing proof of Lemma 4.9. I do not see a fatal flaw in the rest of the error budget, and the overall strategy is convincing. If the authors provide a complete proof of Lemma 4.9 — or a precise reference to a fully worked argument for the stated exponent — I would support publication. The reliance on the authors' own preprints for the energy part of Theorem 2.1 should also be made explicit in the final version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this paper is the first rigorous derivation of Belyakov's 1961 momentum-distribution formula for a low-energy state of the dilute Fermi gas. That is a genuine step: earlier rigorous work stopped at energy bounds, and the physics calculations had known errors. Second, the proof architecture is explicit and the exponent budget is internally consistent; the soft spot is Lemma 4.9, and the stress-test note is right to focus there.\n\nWhat the paper does well. The trial state is the same Huang–Yang-precision state used in the GHNS upper bound, so the result connects energy and correlation structure. The three unitary transformations and the second-order Duhamel expansion are written out, and the error accounting in §5.5.2 does reproduce the claimed ρ^{5/3+1/9} error given Lemma 4.9. The coefficient a² is the squared scattering length entering through scattering theory, not a fitted constant. The two-sided bounds near the Fermi surface also look correct.\n\nSoft spots. The energy statement (2.4) is imported from two preprints with co-author overlap and is not proved here; that is acceptable only if those preprints are refereed and stable. More load-bearing, Lemma 4.9 is sketched in §4.4, and the displayed second inequality is not justified by the cited Lemmas 4.2–4.8 as written. The stress-test's worry is concrete: if a complete proof of that inequality loses a ρ^{1/2} or gains an L^{3/2}, the Gronwall argument gives a weaker bound and the error term in Theorem 2.1(2) need no longer be subleading. I don't see an actual counterexample, and the exponent budget works if the lemma is true, but this is exactly the place I would ask for a full proof before relying on the theorem. The authors also flag the genuine limitations: the result is for the averaged observable and for a trial state, not the true ground state; the unsmeared version would have an L³ error.\n\nBottom line. This deserves a serious referee. I would send it to review with a request to complete Lemma 4.9 and to tie the energy imports to refereed versions. The central claim is a real advance and the paper is honest about what it does not cover. I would cite it.","headline":"Genuinely new rigorous result on the momentum distribution of a dilute Fermi gas, but the proof's load-bearing number-operator lemma is only sketched and needs completion before I'd fully trust the error bar.","tokens_in":50818,"tokens_out":9781,"would_cite":true,"duration_ms":96620,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V70","82B10"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a dilute spin-1/2 Fermi gas, the trial state Ψ=RT1T2Ω has an averaged momentum distribution that matches Belyakov's 1961 formula at leading order for momenta near the Fermi surface.","keywords":["dilute Fermi gas","momentum distribution","Belyakov formula","Huang-Yang energy","excitation density","quasi-bosonic Bogoliubov transformation","thermodynamic limit","trial state"],"falsifier":"Compute, numerically or analytically, the expectation ⟨T2;λΩ, N T2;λΩ⟩ for a finite box at several values of the density; if the growth is faster than L³ρ^{14/9−κ} for any κ>0, the error terms in the Duhamel expansion cease to be subleading and the claimed agreement with Belyakov's formula fails.","tokens_in":49420,"feed_emoji":"⚛️","tokens_out":8018,"duration_ms":83156,"temperature":0.7,"pith_summary":"The paper studies a dilute gas of spin-1/2 fermions in the thermodynamic limit and computes the momentum distribution of a specific trial state, not the true ground state. That trial state, built as three unitary transformations of the free Fermi gas, is the same one used to prove the Huang–Yang upper bound for the ground-state energy. The main result is that its averaged excitation density matches, at leading order, Belyakov's 1961 perturbative formula for momenta close to the Fermi surface, with an error that is strictly smaller than the leading term. If correct, this turns a sixty-year-old formal calculation into a theorem about the momentum distribution of a state that is energetically almost indistinguishable from the ground state; it also explains why simpler trial states would give the wrong constant. The proof works by expanding the excitation-density observable to second order through the unitary transformations and showing that one double commutator reproduces Belyakov's integral while all others are error terms.","feed_headline":"Belyakov momentum formula holds for dilute Fermi trial state","feed_subtitle":"For momenta near the Fermi surface, the trial state's averaged momentum distribution matches the 1961 prediction at leading order.","key_machinery":"The machinery is the unitary conjugation of the density observable through the particle–hole transformation R and two quasi-bosonic Bogoliubov transformations T1 and T2. T1 and T2 are exponentials of quadratic fermion-pair operators with kernels derived from the zero-energy scattering equation (for high momenta) and a regularized scattering kernel (for low momenta). A second-order Duhamel expansion expresses the expectation value as three double commutators; the double commutator with B2 alone produces a constant term which, after removing the infrared cutoff and regularization, is exactly Belyakov's integral. The other double commutators are shown to be error terms using integral estimates","core_discovery":"The paper's central claim is Theorem 2.1: for the trial state Ψ=RT1T2Ω, the averaged excitation density lim sup_{L→∞} |⟨Ψ, n^exc_{q,α}Ψ⟩ − n^{(Bel)}_{q,α}| is bounded by Cρ^{5/3+1/9} in the dilute limit, and for momenta |q|≤Cρ^{1/3} the Belyakov term itself scales like ρ^{5/3+3α}. Thus the 1961 Belyakov formula, originally derived by formal perturbation theory, is shown to be the exact leading-order momentum-space correlation of this particular variational state, provided the imported energy precision and the sketched number-operator bound hold. The energy statement (2.4) is taken from other work; the momentum-distribution statement is proved here via a second-order Duhamel expansion of the","pith_inferences":["The paper proves the Belyakov formula only for the specific trial state; a natural extension is that every trial state accurate enough to resolve the Huang–Yang energy to third order satisfies the same momentum distribution, since the double commutator that produces the Belyakov term is state-independent up to controlled errors.","The ρ^{2−κ}∥ĝ∥₂ error suggests an optimal scale for the momentum smearing window: choosing α just below 1/27 balances the two error terms, and one might expect the window to be tunable to extract the Belyakov term with better precision.","Because the leading term is proportional to a² and independent of the interaction's detailed shape, cold-atom experiments measuring momentum distributions at very low density could in principle measure the s-wave scattering length from the magnitude of the Belyakov term.","The same conjugation-then-expand scheme could be applied to other one-body observables (e.g. static structure factor or occupation of momentum shells) with the same cost: an energy-precise state carries all correlations at the corresponding order."],"forward_implications":["For momenta |q|≤Cρ^{1/3}, the averaged excitation density of Ψ is, up to errors of order ρ^{5/3+1/9}, given by Belyakov's integral; this makes the 1961 formula a theorem for this state rather than a formal perturbative result.","The same state Ψ that provides the Huang–Yang energy upper bound also carries the exact leading-order momentum correlations, so energy precision and correlation precision come together.","The number-operator bound (Lemma 4.9) is the bottleneck: if it holds, all nine double-commutator error terms in Proposition 3.4 are subleading, and the final theorem follows.","Simpler variational states that only achieve the aρ² energy term produce the same order but a different constant, so matching Belyakov's formula is a stricter test of a trial state than matching the energy.","The averaging over momentum balls of radius ρ^{1/3+α} is indispensable; without it, the error would grow with the volume L³ and the theorem would be vacuous."],"fun_headline_variants":["Belyakov formula holds for Fermi trial state at leading order","Trial state recovers 1961 Belyakov momentum distribution","Dilute Fermi gas trial state reproduces Belyakov prediction","Mesoscopic momentum correlations match Belyakov in trial state"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that a certain number-operator bound (proved only in sketch) holds with an exponent around 14/9, and that the trial state's energy error is small enough (a power of the density slightly above 7/3); if either is off by even a small amount, the error in the momentum distribution is as large as the leading Belyakov term.","fun_headline_variants_meta":{"raw":{"variants":["Belyakov formula holds for Fermi trial state at leading order","Trial state recovers 1961 Belyakov momentum distribution","Dilute Fermi gas trial state reproduces Belyakov prediction","Mesoscopic momentum correlations match Belyakov in trial state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000649,"raw_usage":{"total_tokens":2748,"prompt_tokens":606,"completion_tokens":2142,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":350,"completion_tokens_details":{"reasoning_tokens":2070}},"tokens_in":350,"tokens_out":2142,"duration_ms":18658,"temperature":1.0,"reasoning_tokens":2070,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T06:04:15.814410+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, numerically or analytically, the expectation ⟨T2;λΩ, N T2;λΩ⟩ for a finite box at several values of the density; if the growth is faster than L³ρ^{14/9−κ} for any κ>0, the error terms in the Duhamel expansion cease to be subleading and the claimed agreement with Belyakov's formula fails.","supporting_citations":[],"review_version":1}