{"id":"302eb061-901d-419b-a68d-28f586174297","arxiv_id":"1909.00086","paper_version":3,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"Feshbach molecules of 40K87Rb thermally equilibrate during magnetoassociation, and sub-Poissonian density fluctuations in both Feshbach and ground-state molecules independently confirm the molecular Fermi-Dirac occupation.","lead":"Ultracold potassium-rubidium molecules reach thermal equilibrium during formation as a degenerate Fermi gas, and their density fluctuations are directly measured to be sub-Poissonian. The results provide a new, independent thermometer for molecular quantum gases and show how elastic collisions enable efficient cooling.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ~6-collision thermalization estimate is the linchpin; the zero-range cross-section extrapolation is approximate, but ramp-rate and fluctuation-thermometry evidence leave the ACCEPT verdict intact.","rationale":"The paper's strongest claim has two logical parts: (1) Feshbach molecules reach thermal equilibrium during magnetoassociation, and (2) that equilibrium is maintained through STIRAP and is independently confirmed by sub-Poissonian fluctuations. Part (1) depends on the atom-dimer cross-section measurement and the derived collision count, which is the least secure quantitative step. The reader's weakest assumption correctly identifies this. The concern is real but not, in my judgment, fatal: the large measured |aad| is consistent between two independent techniques (oscillation damping and cross-species thermalization), the overlap-density variation in Fig. 1 shows internal consistency, the ramp-rate dependence in Fig. 2 has the expected optimum, and the fluctuation-thermometry agreement in Fig. 4(b) independently establishes that the final molecular sample is in thermal equilibrium. None of the approximations in the collision-count estimate is obviously wrong by an order of magnitude, and the text is transparent about the neglected effects. The central result is therefore acceptable, though the causal 'rapid thermalization' phrasing would be strengthened by a more realistic kinetic simulation. I recommend no change to the reader's ACCEPT verdict.","tokens_in":15895,"tokens_out":17340,"duration_ms":191533,"concrete_test":"Recompute the number of elastic collisions during a 5 ms, 555-to-545.5 G Feshbach ramp with a kinetic Monte Carlo or Boltzmann simulation using the measured aad(B) from the c=0.74(5) fit, time-dependent K and KRb* densities and temperatures, Fermi-Dirac final-state suppression of K collisions, and a velocity-dependent T-matrix cross section; compare the cumulative collision number at the end of the ramp with the ~3-collision thermalization threshold. If the simulated Nel falls below threshold, the abstract's causal 'thermalization ... enables' statement should be downgraded to 'the final molecular sample is in thermal equilibrium', and the verdict should be CONDITIONAL on additional collision-rate data; if Nel remains above threshold, the ACCEPT verdict is fully supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central causal assertion that Feshbach molecules rapidly thermalize during the ramp rests on the estimated number of elastic collisions per molecule, computed in the Supplementary Material ('Number of elastic collisions'). That estimate multiplies an overlap density n=1.6e13 cm^-3, a field-averaged cross section built from the zero-range formula sigma = 4 pi a_ad^2/(1 + k_th^2 a_ad^2), and a 0.6 ms near-resonance window, using Fermi-temperature energy scales (T_K^F=680 nK, T_KRb*^F=220 nK) in vrel and k_th. Each factor is an approximation: the measured |aad| fit c=0.74(5) is anchored at B-B0 < -0.3 G, outside the near-resonance field range that dominates the integral; the overlap density is Boltzmann and time-independent during association; and the text itself notes that Pauli blocking of K collisions and the state distributions of K and KRb* are not included. If the true collision number falls below the roughly 2.7 collisions conventionally required for s-wave thermalization, then the claim that collisions cause the observed degeneracy is not quantitatively established, because the low T/T_F and thermal-looking momentum distributions could be adiabatically inherited from the deeply degenerate K-Rb starting mixture. The ramp-rate optimum and the agreement between expansion and fluctuation thermometry support but do not fully isolate the causal role of elastic collisions. This is the same soft spot the reader identified; it does not invalidate the final-state equilibrium finding, but the 'rapid thermalization' claim is only as secure as the collision-count estimate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports measurements on a degenerate Fermi gas of 40K87Rb molecules produced by magnetoassociation from a Bose–Fermi mixture. The authors measure the K–KRb* atom–dimer elastic scattering length from damping of center-of-mass oscillations and from cross-species thermalization, obtaining a scattering length magnitude well described by a scaling factor c = 0.74(5) relative to the K–Rb scattering length. They observe that the molecular T/T_F is minimized at intermediate Feshbach ramp rates, indicating a competition between thermalization and inelastic losses. They further measure sub-Poissonian number fluctuations in K, Feshbach molecules, and ground-state molecules after time-of-flight, and find that T/T_F inferred from fluctuation suppression agrees with expansion thermometry. The central claims are that Feshbach molecules rapidly thermalize with the parent atoms during the association ramp, that equilibrium is essentially maintained through coherent transfer to the ground state, and that fluctuation measurements directly probe the molecular Fermi–Dirac distribution.","tokens_in":16201,"tokens_out":9453,"duration_ms":90690,"significance":"These results constitute the first measurement of atom–dimer elastic scattering for a heteronuclear molecule in this regime and provide a quantitative explanation for the production of degenerate molecular samples at T/T_F = 0.3. The fluctuation thermometry is a new tool for molecules, benchmarked on non-degenerate K atoms and validated by both simulation and agreement with expansion thermometry. If the claims hold, the work strengthens the path to degenerate polar gases and offers a direct probe of quantum statistics in molecules. The experimental analysis is careful: the damping model is checked against the full two-species equations, the two independent measurements of the scattering length agree, and the STIRAP efficiency and imaging corrections are explicitly characterized.","major_comments":[],"minor_comments":[{"comment":"The estimate of Nel = 5.9 relies on a field-averaged cross section that extrapolates the measured c = 0.74 to the near-resonance region, a Boltzmann overlap density, and Fermi-temperature energy scales, while explicitly neglecting Pauli blocking of K collisions and the state distributions of K and KRb*. A sensitivity analysis (for example, varying the near-resonance contribution or the overlap density) or a more cautious wording such as 'on the order of several collisions' would strengthen the quantitative support for the 'more than 6 elastic collisions' statement, although this does not affect the central conclusion given the independent ramp-rate and cross-species thermalization evidence.","section":"Supplementary Material, 'Number of elastic collisions'"},{"comment":"The rise in T/T_F at fast ramp rates is attributed to hindered thermalization, but non-adiabatic effects during molecule association could also contribute; a brief discussion of this alternative would aid the interpretation.","section":"Main text, Fig. 2"},{"comment":"The variance suppression formula is presented without derivation; a short explanation of the local density approximation and the treatment of the expanded potential V(x,z) would improve accessibility.","section":"Main text, Eq. (1)"},{"comment":"The placeholder 'xxxxx' for the Supplementary Materials link should be replaced with the actual DOI or journal URL.","section":"Main text, references"},{"comment":"The simulation-extracted slope of 0.4 is compared to the experimental value 0.45(4) before the 2.2 scaling factor is applied; stating this explicitly would avoid confusion about whether the simulation captures the resolution correction.","section":"Supplementary, 'Finite bin size effects'"},{"comment":"The chemical formula 40K87Rb appears with inconsistent superscript formatting; please use a consistent notation (e.g., 40K87Rb with superscripts or a plain text equivalent).","section":"Throughout"}],"recommendation":"accept","confidential_remarks":"The manuscript is a strong experimental letter suitable for the journal. The collision-number estimate in the supplement is approximate, but the authors acknowledge its limitations, and the central thermalization claim is independently supported by the ramp-rate dependence and the cross-species thermalization measurement. No concerns about novelty or scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this is a solid experimental letter that reports two genuinely new measurements—the first field-dependent K–KRb* atom–dimer scattering length, and the first observation of sub-Poissonian number fluctuations in a degenerate molecular gas. Both claims survive close reading. The thermalization narrative—Feshbach molecules equilibrate with the atomic bath before coherent transfer—is supported by multiple independent observables: the large |a_ad| from collisional damping, the cross-species thermalization agreement, the ramp-rate optimum, and consistency between expansion and fluctuation thermometry.\n\nThe strongest part is the fluctuation experiment. They calibrate imaging carefully (α, saturation, readout), correct for finite bin size with both a measured scaling factor and a simulation, benchmark on non-degenerate K, and account for STIRAP-induced hole burning. The polylogarithm fit to the variance profile giving T/T_F in agreement with expansion thermometry is convincing.\n\nThe soft spot is the quantitative claim of 'more than 6 elastic collisions per molecule'. That number is an estimate built from the zero-range cross-section formula, a Boltzmann overlap density, Fermi-temperature energy scales, and a field-average over a range where the fitted |a_ad| is extrapolated. The authors themselves flag the omissions (Pauli blocking, state distributions). If the true collision number falls below ~3, the causal claim that collisions produce the observed degeneracy loses quantitative grounding. But the ramp-rate optimum and the agreement between thermometry methods argue that something is thermalizing the sample; the final-state equilibrium result does not depend on the exact collision count. So this is a caveat, not a flaw.\n\nMinor: the Rb–KRb* scattering length data are noisy, but they're clearly labeled and supplementary.\n\nThis paper deserves serious peer review and is publishable after minor revision—ideally with a sensitivity analysis of the collision-number estimate. It will be cited by anyone working on degenerate molecular gas production.","headline":"Solid experimental letter: first atom–dimer scattering length for K–KRb* and first sub-Poissonian molecular gas fluctuations; thermalization claim is strong, with the collision-count estimate as the main quantitative soft spot.","tokens_in":16734,"tokens_out":1858,"would_cite":true,"duration_ms":18364,"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":"Feshbach molecules produced from a degenerate Bose–Fermi mixture rapidly reach thermal equilibrium with both parent atomic species, and the equilibrium is essentially preserved through coherent transfer to the ground state.","keywords":["ultracold polar molecules","Feshbach molecules","degenerate Fermi gas","atom-dimer scattering length","thermalization","sub-Poissonian number fluctuations","Fermi-Dirac distribution","STIRAP transfer"],"falsifier":"Measure the same K–KRb* elastic cross section by an independent method, e.g., from the loss rate of a trapped KRb* sample in a controlled K bath due to three-body recombination at the same fields and temperatures, and compare the extracted |a_ad| with the c = 0.74(5) scaling; a significant discrepancy would invalidate the thermalization-rate claim. Alternatively, prepare a molecular sample at T/T_F < 0.1 and check whether the variance profile still matches the predicted form with the separately measured T/T_F; if it does not, the claim that fluctuations directly probe the Fermi–Dirac distribution fails.","tokens_in":15728,"feed_emoji":"⚛️","tokens_out":7896,"duration_ms":64810,"temperature":0.7,"pith_summary":"The paper shows that during magnetoassociation, KRb Feshbach molecules scatter elastically off the surrounding potassium atoms often enough to fully thermalize before the ramp ends, and that subsequent coherent transfer to the rovibrational ground state leaves the gas essentially in equilibrium. The authors extract a large atom–dimer scattering length for K–KRb* collisions, |a_ad| ≈ 0.74 times the K–Rb scattering length, implying more than six elastic collisions per molecule in a typical ramp. They also measure sub-Poissonian density fluctuations in both Feshbach and ground-state molecules, a first for molecules, and show that the fluctuation profile gives a T/T_F consistent with expansion thermometry. Why this matters: reaching high phase-space density in polar molecules has been limited by weak ground-state interactions, and this work shows the thermalization bottleneck can be passed before the molecules ever reach the ground state.","feed_headline":"Feshbach molecules thermalize as they are formed, preserving degeneracy","feed_subtitle":"Atom-dimer collisions equilibrate KRb pairs; fluctuation measurements confirm the Fermi-Dirac distribution.","key_machinery":"The two load-bearing tools are (i) collisional damping of center-of-mass oscillations, where the decay rate of intentionally excited KRb* oscillations in a K cloud gives the elastic collision rate $\\Gamma = n \\sigma v_{\\rm rel}$ and hence the atom–dimer scattering length through the s-wave formula $\\sigma = 4\\pi a_{\\rm ad}^2/(1 + k_{\\rm th}^2 a_{\\rm ad}^2)$; and (ii) local number-fluctuation thermometry, where the variance-to-mean ratio $\\sigma^2_N/N = \\mathrm{Li}_1(-\\zeta e^{-V/kT})/\\mathrm{Li}_2(-\\zeta e^{-V/kT})$ is fit to binned image data to extract the fugacity and T/T_F. A third element is the Feshbach ramp-rate sweep, which maps the competition between thermalizing elastic collisions and inelastic losses and confirms that thermalization saturates for intermediate ramp rates.","core_discovery":"The central discovery is that the production of a degenerate Fermi gas of polar 40K87Rb molecules is self-thermalizing: elastic collisions between Feshbach molecules and free potassium atoms during the magnetic-field ramp equilibrate the molecular sample, and the equilibrium survives the STIRAP transfer with only a small perturbation. This is established by measuring the damping of molecular center-of-mass oscillations in a potassium bath, which yields a K–KRb* s-wave atom–dimer scattering length |a_ad| that scales with the K–Rb scattering length a(B) as |a_ad| = (0.74 ± 0.05) a(B) across the tested detunings. The measured cross section implies at least six elastic collisions per molecule during a standard ramp, and the ramp-rate dependence of the final T/T_F shows an optimum around 0.5–3 G/ms, consistent with elastic thermalization winning over inelastic loss. As an independent probe of degeneracy, the paper measures the variance-to-mean ratio of molecule number in small image bins and observes sub-Poissonian fluctuations whose spatial profile matches the Fermi–Dirac prediction, giving a T/T_F in agreement with expansion measurements. The STIRAP transfer, at 85% efficiency, is shown to increase the apparent T/T_F only modestly (e.g., from 0.44 to 0.49 for one set), so the ground-state gas remains essentially in equilibrium.","pith_inferences":["If the variance-suppression profile is sensitive to two-body correlations, the same imaging technique could detect the onset of pairing or a BCS-type crossover in a molecular Fermi gas, where standard expansion thermometry would need to be modified.","The measured c ≈ 0.74, below the universal prediction of 1.09, suggests that finite-range effects matter at these temperatures; a systematic study of |a_ad| versus collision energy could map where universality sets in.","Because STIRAP acts as a binomial state-depletion process, the fluctuation method can be inverted to measure the STIRAP conversion efficiency per momentum state, providing a diagnostic that is independent of total-number calibrations."],"forward_implications":["Degenerate polar ground-state molecules can be produced without post-transfer evaporative cooling, because the thermalization bottleneck is passed during the Feshbach ramp.","The measured K–KRb* atom–dimer scattering length, |a_ad| ≈ 0.74 a(B), provides a quantitative benchmark for three-body and atom–dimer scattering theory in heteronuclear systems at finite collision energy.","Sub-Poissonian number fluctuations in molecules can serve as a local, model-free probe of the molecular Fermi–Dirac distribution, complementing expansion thermometry in regimes where expansion is distorted.","The optimal magnetoassociation ramp rate for KRb is set by the balance between elastic thermalization and inelastic loss, with T/T_F minimized at 0.3 for rates near 1 G/ms.","STIRAP efficiency sets a floor on the achievable ground-state degeneracy unless ground-state thermalization occurs; improving transfer efficiency directly improves T/T_F."],"supporting_citations":[{"why":"Supplies the method and initial conditions for producing a degenerate Fermi gas of KRb molecules, the system under study.","marker":"[18]"},{"why":"Provides the universal prediction a_ad = 1.09a for this mass ratio, the baseline against which the measured c = 0.74(5) is compared.","marker":"[46]"},{"why":"Established sub-Poissonian number-fluctuation thermometry in Fermi gases, the technique adapted here to molecules.","marker":"[34]"},{"why":"Supplies the polylog variance-profile model (Eq. 1) used to extract T/T_F from fluctuation data.","marker":"[35]"},{"why":"Demonstrated that bosonic Feshbach molecules thermalize through elastic collisions, the mechanism this paper extends to the heteronuclear fermionic case.","marker":"[22]"},{"why":"Characterized inelastic K–KRb* and Rb–KRb* collisions, used to choose the initial Rb number and interpret ramp-rate losses.","marker":"[29]"},{"why":"Pioneered the center-of-mass oscillation-damping method for measuring elastic cross sections, the technique used to extract a_ad.","marker":"[41]"},{"why":"Provides the Feshbach-resonance and scattering-length formalism, including the s-wave cross-section formula.","marker":"[45]"}],"fun_headline_variants":["Self-thermalizing KRb Fermi gas holds its degeneracy","Atom-dimer scattering equilibrates molecular gas, preserves quantum state","Sub-Poissonian fluctuations probe Fermi-Dirac distribution directly","Degenerate polar molecules thermalize via collisions with atoms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the measured damping of center-of-mass oscillations is accurately converted to an atom–dimer scattering length using the s-wave zero-range formula and a Boltzmann overlap density, approximations that may be inadequate at the high collision energies and large scattering lengths probed.","fun_headline_variants_meta":{"raw":{"variants":["Self-thermalizing KRb Fermi gas holds its degeneracy","Atom-dimer scattering equilibrates molecular gas, preserves quantum state","Sub-Poissonian fluctuations probe Fermi-Dirac distribution directly","Degenerate polar molecules thermalize via collisions with atoms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000207,"raw_usage":{"total_tokens":1395,"prompt_tokens":932,"completion_tokens":463,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":392}},"tokens_in":548,"tokens_out":463,"duration_ms":4474,"temperature":1.0,"reasoning_tokens":392,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T06:01:40.390487+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the same K–KRb* elastic cross section by an independent method, e.g., from the loss rate of a trapped KRb* sample in a controlled K bath due to three-body recombination at the same fields and temperatures, and compare the extracted |a_ad| with the c = 0.74(5) scaling; a significant discrepancy would invalidate the thermalization-rate claim. Alternatively, prepare a molecular sample at T/T_F < 0.1 and check whether the variance profile still matches the predicted form with the separately measured T/T_F; if it does not, the claim that fluctuations directly probe the Fermi–Dirac distribution fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the universal prediction a_ad = 1.09a for this mass ratio, the baseline against which the measured c = 0.74(5) is compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Characterized inelastic K–KRb* and Rb–KRb* collisions, used to choose the initial Rb number and interpret ramp-rate losses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Pioneered the center-of-mass oscillation-damping method for measuring elastic cross sections, the technique used to extract a_ad."}],"review_version":1}