{"id":"9638c090-8801-48e6-9301-ca9bbbc92270","arxiv_id":"2511.19191","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Rydberg atom absorption spectra in a BCS superfluid are predicted to shift by the pairing gap for dimer and mixed trimer states, while unshifted same-angular-momentum trimer peaks signal intact Cooper-pair trapping.","lead":"The paper proposes using a Rydberg atom embedded in an ultracold paired fermion gas as a local spectroscopic sensor. Its calculations show that the atom's absorption spectrum can reveal the superfluid pairing gap and whether Cooper pairs are broken or trapped intact.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gap-measurement claim rests on unverified 3D homogeneity of Δ after sudden Rydberg excitation; no time-dependent or self-consistent check is supplied.","rationale":"The paper is a proposal for a local spectroscopic probe of BCS pairing using a Rydberg impurity. The FDA calculation is a reasonable and technically careful approach for a quadratic mean-field Hamiltonian, with convergence checks (ℓ_max = 5, k_F R = 200) and a clear falsifiable prediction in Eq. (7). The reader's CONDITIONAL verdict is appropriate. The single most load-bearing assumption is that the order parameter Δ remains homogeneous after the sudden excitation. The paper explicitly acknowledges this assumption and cites 1D equilibrium back-reaction studies, but does not provide a 3D self-consistent or time-dependent check. Because the absorption spectrum is a long-time observable, the sudden-excitation timescale argument alone does not guarantee that Δ stays uniform during the Ramsey measurement; the order parameter can evolve on timescales of order 1/Δ, which are precisely the timescales probed in Fig. 3. This concern directly affects the central quantitative claim (Eq. 7) and the 2Δ trimer shifts. I agree with the reader's weakest-assumption identification. The Cooper-pair-trapping interpretation is also underived, but it is ancillary to the gap-measurement claim and is acknowledged by the authors as requiring future wavefunction-based studies; it does not move the verdict beyond CONDITIONAL. A concrete TDBdG simulation would settle whether the frozen-Δ model is adequate, so no change to the reader's verdict is needed.","tokens_in":14497,"tokens_out":16239,"duration_ms":164876,"concrete_test":"Perform a 3D time-dependent Bogoliubov–de Gennes (TDBdG) simulation with the same parameters (n_Ryd = 60, a_e = −15 a0, ρ0 = 5×10^11 cm^−3, k_F R = 200): initialize in the uniform-Δ BCS ground state, suddenly switch on VR(r), and evolve Δ(r,t) self-consistently. From the evolved state compute the Ramsey signal S(t) and absorption spectrum A(ω), and check whether the dimer peak shifts by Δ, mixed-ℓ trimer peaks by 2Δ, and the same-ℓ trimer peak remains unshifted to within a few percent of Δ over the full window t ≲ 100/ε_F. Report the maximum local deviation |Δ(r,t) − Δ0|/Δ0; if it exceeds ~10% for t ≳ 1/Δ, the frozen-Δ model is not quantitatively reliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the dimer-peak shift directly measures the superfluid gap via Eq. (7)—rests on the model Hamiltonian (4), in which the Rydberg potential VR(r) is added to the single-particle sector while the order parameter Δ is kept rigidly uniform. The paper states after Eq. (4) that the sudden excitation is fast compared to equilibration, 'ensuring that the order parameter Δ remains homogeneous in the vicinity of the impurity during the procedure.' This justification is incomplete. The Ramsey signal S(t) in Eq. (5) and the absorption spectrum A(ω) in Eq. (6) are computed for long times; Fig. 3 uses tε_F up to 100, and the quasiparticle regime is identified for t ≳ 1/Δ. Over these timescales, a local potential deep enough to bind a state at |E_D| ≈ 422 ε_F can locally deplete or reconstruct the pair condensate in 3D. The cited 1D self-consistent calculations [17,18] already show significant gap modification near a Rydberg impurity, and no 3D self-consistent or time-dependent check is provided. If Δ(r,t) develops spatial structure or changes magnitude on the measurement timescale, the simple energy bookkeeping behind Δ ≈ |E_D^(Δ)−E_D| and the 2Δ trimer shifts breaks down, so the 'direct measure' claim would require corrections. Additionally, Eq. (7) is presented as a numerical observation rather than derived analytically, which makes the assumed rigidity of Δ all the more load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes using a Rydberg impurity as a local spectroscopic probe of a 3D BCS superfluid. The impurity is excited suddenly, and the absorption spectrum A(ω) is computed from the Ramsey signal S(t) via the functional determinant approach adapted to a mean-field BCS Hamiltonian. The central claims are: (i) the deep dimer peak shifts by the superfluid gap, Eq. (7); (ii) mixed-angular-momentum trimer peaks shift by 2Δ while same-angular-momentum trimer peaks do not shift, corresponding to breaking versus intact trapping of Cooper pairs; and (iii) the gap suppresses the orthogonality catastrophe, producing a quasiparticle residue Z_B ~ (Δ/ε_F)^α with exponent α given by the scattering phase-shift sum, Eq. (9). The Supplement gives the FDA reduction and numerical details.","tokens_in":14863,"tokens_out":8424,"duration_ms":94245,"significance":"If the predictions hold, this is a valuable and experimentally plausible new probe of pairing in ultracold Fermi gases. The deep dimer shift, the 2Δ trimer shifts, and the no-shift same-angular-momentum trimer peaks are concrete, falsifiable signatures, and Eq. (9) is compared with an independently computed phase-shift expression rather than extracted from the same fit. The FDA derivation in Supplement S2 is standard and carefully laid out, and the numerical implementation uses realistic parameters for the 87Rb–40K system with convergence checks in angular momentum and system size. The main weakness is that all central spectral predictions are obtained with a rigidly uniform gap Δ; the justification for this approximation is only a timescale argument, and no 3D self-consistent or time-dependent check is supplied. This issue is load-bearing for the 'direct measure' claim and must be addressed before publication.","major_comments":[{"comment":"The assumption that Δ remains homogeneous after the sudden excitation is load-bearing but is justified only by a short timescale argument. The Rydberg potential binds a deep state at |E_D|≈422 ε_F, and the measurement uses times t≳1/Δ up to t ε_F=100 (Fig. 3); over such times a local potential of this strength can deplete or reconstruct the pair condensate. The cited 1D studies [17,18] already show significant gap back-reaction near a Rydberg impurity, and no 3D self-consistent BdG calculation or time-dependent estimate is provided. Since Eq. (7) and the 2Δ trimer shifts assume uniform Δ, a spatially dependent Δ(r,t) would modify the energy bookkeeping. Please add a 3D self-consistent calculation for the final static Hamiltonian, or at least a quantitative estimate of the healing time and of the corrections to Eq. (7), and state the regime in which the rigid-Δ result is valid.","section":"Model (Eq. (4))"},{"comment":"The claim in the main text and in the caption of Fig. 2(c) that all dimer peaks shift by Δ is contradicted by the Supplement. For the weakly bound n_b=3, ℓ=2 state, Eq. (S30) gives E_D^(Δ)≈E_D+Δ−δ_D(Δ), with δ_D visible for Δ/ε_F>0.2, and Eq. (S31) gives a separate anomalous δ_T(Δ) for the corresponding trimer. These deviations from Eq. (7) are not predicted theoretically. Since a real spectrum contains several bound states, an experiment using Eq. (7) needs either a reliable identification of the deep dimer peaks or a characterization of δ_D and δ_T. The main text should restrict the generality of its 'all peaks' statement and explicitly discuss this limitation.","section":"Supplement S4; Fig. 2(c)"},{"comment":"The exponent α in Eq. (9) is a central quantitative claim, but its numerical determination is not documented with systematic uncertainties. The quasiparticle weights Z_B and Z_γ are obtained by fitting the Ramsey signal to Eq. (8) over t≳1/Δ, after a Gaussian filter, within a window limited by t_max=100/ε_F; the power law in Δ is then read from the inset of Fig. 3. The inset appears to contain few points and no error bars or sensitivity to the fit interval. Please report the fit ranges, the number of independent Δ values, the systematic error from the finite time window, and a direct numerical-versus-analytic comparison of α with Eq. (9).","section":"Orthogonality catastrophe (OC) / Fig. 3 and Supplement S5"}],"minor_comments":[{"comment":"Typo: 'Strinkingly' should be 'Strikingly'.","section":"Introduction"},{"comment":"The angular-momentum cutoff is stated as ℓmax=5 in the main text and ℓmax≈5 in the Supplement; please make the notation consistent.","section":"Model and Supplement S3"},{"comment":"The display in Eq. (S15) is unclear: the two-component vector (u_n^2, v_n^2)^T is written as a fraction with a ± sign. Please rewrite it as an explicit two-component expression.","section":"Supplement S2, Eq. (S15)"},{"comment":"The phrase 'out-of-equilibrium dynamics of a three-dimensional configuration' could be misread as a time-dependent simulation of Δ. What is actually computed is the sudden-overlap spectrum with Δ held fixed. Please phrase this more precisely.","section":"After Eq. (4)"},{"comment":"The caption's statement 'All peaks shift proportional to the gap strength Δ' is too broad, given the anomalous weak-state shifts in Supplement S4. Please qualify it as applying to the deep dimer bound states considered in the main text.","section":"Fig. 2(c) caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is interesting and likely to attract wide attention. My recommendation of major revision is driven by the homogeneous-Δ approximation, which is central to Eq. (7) and the 2Δ trimer shifts. I would be convinced by a 3D self-consistent BdG check for the same parameters, even if it only shows that corrections are small in the sudden regime, together with a more careful statement of the validity range of Eq. (7) in view of the anomalous weak-state shifts."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here is my take. The paper does something genuinely useful: it works out, in a controlled mean-field model, what a Rydberg impurity's absorption spectrum looks like when immersed in a BCS superfluid, and it identifies clean signatures of pairing. The dimer peak shifts by Δ, the same-angular-momentum trimer peak doesn't shift, the mixed one shifts by 2Δ, and the orthogonality catastrophe is suppressed with an exponent given by scattering phase shifts. Those are concrete, falsifiable predictions for ultracold-gas experiments.\n\nThe technical core is sound. The FDA reduction to a single-particle problem is standard, and the implementation with the Bogoliubov transformation is careful. The OC exponent relation (Eq. 9) is a nice result: extracting Z_B from the Ramsey fits and comparing it to the independent phase-shift expression gives a genuine check, and the agreement is visible in the inset. The authors are also honest about the things they didn't do: the anomalous shifts of weakly bound states, the lack of a molecular wavefunction calculation for the intact-pair trapping, and the fact that Eq. (7) is a numerical observation rather than a derivation.\n\nThe soft spot that matters is the rigid Δ assumption. The authors state, after Eq. (4), that the sudden excitation 'ensures' Δ remains homogeneous, but they compute the Ramsey signal for times up to tε_F = 100. A potential that binds a state at 422 ε_F is strong, and the cited 1D back-reaction calculations show real local gap modification. No 3D self-consistent or time-dependent check is presented. If Δ(r,t) changes on the measurement timescale, the simple energy bookkeeping behind Eq. (7) fails, and the 'direct measure' claim needs qualification. This is not a fatal flaw in the model—within the rigid-Δ model the predictions are correct—but it is load-bearing enough that a referee should push on it.\n\nMinor quibble: the fitted quasiparticle weights come without error bars, and the intact-pair trapping interpretation is asserted more strongly than the calculation supports (the authors do flag this).\n\nWho is this for? Cold-atom theorists planning Rydberg spectroscopy of paired Fermi gases, and experimentalists looking for a local probe. It deserves a serious referee. I would send it out and ask for a quantitative estimate of back-reaction timescales or a 3D self-consistent calculation, and a derivation of Eq. (7) if possible. Otherwise the main claims are likely to stand.","headline":"A clean, promising theory paper that turns a Rydberg impurity into a local gap spectrometer; the main open question is whether the assumed rigid Δ survives the strong local potential over the long times used for the OC analysis.","tokens_in":15355,"tokens_out":3209,"would_cite":true,"duration_ms":34974,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.75.Ss","32.80.Ee","67.85.-d"],"model":"deepseek-v4-flash","headline":"Rydberg impurities immersed in a BCS superfluid give an absorption spectrum that encodes the pairing gap, distinguishes broken from intact Cooper pairs, and reveals how the gap stabilizes quasiparticles.","keywords":["BCS superfluid","Rydberg impurity","Cooper pairs","superfluid gap","absorption spectrum","orthogonality catastrophe","polaron quasiparticle","ultracold Fermi gases"],"falsifier":"A decisive check is a three-dimensional self-consistent calculation of the pairing gap around the Rydberg impurity: if the local gap at the Rydberg shell deviates from the bulk value by more than the linewidth, the dimer shift will depart from Δ and the nominally unshifted same-angular-momentum trimer peak will acquire a gap dependence. Experimentally, compare the dimer-peak shift against an independent bulk-gap measurement, such as rf spectroscopy, across interaction strengths.","tokens_in":14370,"feed_emoji":"⚛️","tokens_out":10610,"duration_ms":98632,"temperature":0.7,"pith_summary":"This paper proposes putting a Rydberg-excited atom inside a superfluid of ultracold paired fermions and using its optical absorption spectrum as a local, time-resolved sensor of the BCS state. The central claim is that the spectrum's molecular peaks are fingerprints of Cooper pairing: when one fermion is bound to the Rydberg atom, a Cooper pair must be broken and the peak shifts by the pairing gap, with Δ ≈ |E_D^(Δ) − E_D|; trimer peaks formed from two broken pairs shift by 2Δ, while trimer peaks that bind a whole Cooper pair stay fixed. This gives a direct local measurement of the gap and, by scanning the probe, a spatial map of pairing. The paper further claims that the same gap suppresses the orthogonality catastrophe that normally kills quasiparticles in an impurity in a Fermi sea, producing a sharp quasiparticle peak with weight scaling as (Δ/ε_F)^α and an exponent fixed by scattering phase shifts. A sympathetic reader would care because this is a concrete, experiment-ready route to look inside a strongly correlated superfluid at microscopic scales.","feed_headline":"A Rydberg atom reads out the BCS pairing gap","feed_subtitle":"Dimer peaks shift by the gap; unshifted trimer peaks signal intact Cooper pairs.","key_machinery":"The machinery is the Rydberg potential V_R(r) = 2π a_e/m_e |ψ_nRyd(r)|^2 together with the functional determinant approach applied to the BCS Hamiltonian. The impurity excitation is sudden, so the measured spectrum is determined by the Ramsey overlap S(t) = Tr{ρ e^{iH0 t} e^{-iH_R t}} between the gas with and without the Rydberg potential; the absorption spectrum is its Fourier transform. Adapting the determinant formula to Nambu fermions, the BCS Hamiltonian is diagonalized by a Bogoliubov transformation into two gapped branches, ω_{σn} = ±√(ξ_n^2 + Δ^2) in the balanced case, and S(t) becomes a determinant built from single-particle overlaps of the free and potential-perturbed states. The s","core_discovery":"The central discovery, on the paper's own terms, is that a suddenly excited Rydberg atom in a BCS superfluid acts as a local spectroscopic probe of pairing. Computing the absorption spectrum with the functional determinant approach adapted to the BCS state, the authors find that dimer peaks—one fermion occupying the Rydberg bound state—shift by the superfluid gap, Δ ≈ |E_D^(Δ) − E_D|, because forming the dimer breaks a Cooper pair and must pay its binding energy. Trimer peaks with mixed angular momenta shift by 2Δ, the cost of breaking two Cooper pairs, whereas trimer peaks with both fermions in the same angular-momentum channel do not shift, revealing that an intact Cooper pair can bind to","pith_inferences":["The paper leaves implicit that scanning the Rydberg probe through a trapped gas could resolve spatial structure of pairing such as vortex cores or the BCS–BEC crossover; the peak-shift rules provide a concrete observable for such maps.","A testable extension: because the exponent α is written purely in terms of scattering phase shifts at the Fermi energy, one could extract δ_ℓ(ε_F) from separate single-particle scattering measurements and compare the predicted (Δ/ε_F)^α scaling in the Ramsey signal.","The anomalous shifts for weakly bound states reported in the supplement mean the clean Δ and 2Δ rules are tied to deeply bound states; a robust measurement protocol should either avoid the weak states or model their corrections, which the paper leaves for future work."],"forward_implications":["A dimer-peak shift gives a local, quantitative readout of the superfluid gap; scanning the Rydberg beam spatially reconstructs the pairing gap, including its variation in inhomogeneous clouds.","Same-angular-momentum trimer peaks identify intact Cooper-pair trapping, a process that remains possible even when the coherence length is larger than the Rydberg potential, meaning pairs can be bound without losing their internal coherence.","The superfluid gap suppresses the orthogonality catastrophe; the atomic branch becomes a quasiparticle peak with weight scaling as (Δ/ε_F)^α, linking local pairing to polaron quasiparticle physics.","A Bogoliubov sideband at 2Δ appears on the atomic peak, providing a second, complementary route to measure the gap and to excite quasiparticles across the gapped dispersion.","The method generalizes to finite temperature and unbalanced spin mixtures, so it could probe pairing beyond the balanced mean-field BCS case."],"fun_headline_variants":["Rydberg probe measures superfluid gap","Dimer shift unlocks the BCS gap","Intact Cooper pairs read out by Rydberg","Spectroscopy exposes BCS pairing locally","Rydberg atom reads pair-breaking energy"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that during the sudden laser excitation the superconducting pairing field stays homogeneous around the Rydberg atom and the gas never equilibrates with it (stated after Eq. (4)); if the strong Rydberg potential locally deforms or depletes the pairs in three dimensions, the claimed peak-shift rules would change.","fun_headline_variants_meta":{"raw":{"variants":["Rydberg probe measures superfluid gap","Dimer shift unlocks the BCS gap","Intact Cooper pairs read out by Rydberg","Spectroscopy exposes BCS pairing locally","Rydberg atom reads pair-breaking energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000545,"raw_usage":{"total_tokens":2433,"prompt_tokens":723,"completion_tokens":1710,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":467,"completion_tokens_details":{"reasoning_tokens":1643}},"tokens_in":467,"tokens_out":1710,"duration_ms":14181,"temperature":1.0,"reasoning_tokens":1643,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T20:33:10.643604+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is a three-dimensional self-consistent calculation of the pairing gap around the Rydberg impurity: if the local gap at the Rydberg shell deviates from the bulk value by more than the linewidth, the dimer shift will depart from Δ and the nominally unshifted same-angular-momentum trimer peak will acquire a gap dependence. Experimentally, compare the dimer-peak shift against an independent bulk-gap measurement, such as rf spectroscopy, across interaction strengths.","supporting_citations":[],"review_version":1}