{"id":"59776aa7-9ffb-48ac-a8dc-e58a710aa3ba","arxiv_id":"1909.02649","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Measurements in cold 87Rb show the sigma-plus sigma-minus CPT dark-state pumping rate is an order of magnitude slower than the three-level formula, and a parameter-free 13-level Zeeman model matches the experiment to within 30%.","lead":"Cold-atom experiments show that coherent dark states in rubidium atoms are pumped about ten times slower than the standard three-level formula predicts, and a 13-level Zeeman model explains the gap. The corrected pumping rate also explains the observed light-shift dependence in Ramsey spectroscopy and shows that dark-state coherence loss is dominated by atoms moving out of the beam, not by internal decoherence.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 13-level model's no-fit agreement is the sole evidence for the Zeeman-manifold explanation, but its assumed uniform, coherence-free initial F=2 population after molasses is unmeasured and likely non-generic.","rationale":"The reader's weakest-assumption analysis correctly identifies the unsupported initial-state ansatz as the most load-bearing premise. The strongest claim is the factor-of-ten discrepancy and its explanation by the 13-level model; that explanation depends on the model's quantitative agreement, which in turn depends on the initial density matrix. The paper explicitly acknowledges the assumption but provides no justification, and an optical molasses is not expected to produce a strictly equal, coherence-free Zeeman population. Other concerns, such as absolute intensity calibration and the normalized light-shift comparison, are secondary because the model's agreement with the measured pumping rate across several intensities already constrains the intensity scale, and the light-shift comparison is presented only as a consistency check on the relation between pumping duration and light shift, not as the primary evidence. The appropriate response is to keep the conditional verdict and require the authors to supply the model equations, parameters, and either a measured initial state or a sensitivity analysis over plausible initial distributions.","tokens_in":8243,"tokens_out":16438,"duration_ms":199843,"concrete_test":"Re-run the 13-level simulation at I = 1.33 W/m2 with an initial density matrix given by a realistic D2 optical-molasses steady-state population distribution (and with small coherences included), instead of the uniform F=2 mixture; if the predicted 1/e pumping rate changes by more than 30% relative to the uniform-mixture prediction, the no-fit agreement is not robust and the attribution to the Zeeman manifold is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central attribution—that the roughly tenfold slowdown of the dark-state pumping rate is caused by the multilevel Zeeman structure—rests on the 13-level density-matrix model reproducing the measured rate to within ~30% with no fit parameters. The model's most fragile input is the initial density matrix: Section III states, 'Due to the cooling process, the initial conditions assume that all the Zeeman sub-levels of the F=2 level are equally populated without any coherence.' No in-situ measurement, calibration, or independent model of the optical-molasses steady state is provided. Optical molasses is known to produce non-thermal, polarization- and intensity-dependent Zeeman population distributions, and the pumping dynamics into the dark state depend on the initial populations of the coupled mF sublevels. If the true initial state deviates from the uniform mixture, the predicted pumping rate and the non-exponential shape of the build-up curve shift by an unquantified amount. Because the claimed no-fit agreement is the only evidence that the discrepancy is due to the Zeeman manifold rather than to an incorrect baseline or an uncalibrated input, the central claim is conditional on this assumption. The authors should either measure or bound the initial Zeeman distribution, or show that the model's predictions are insensitive to it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports measurements of dark-state optical pumping and decay in laser-cooled, free-falling 87Rb atoms interrogated in the σ+–σ− CPT configuration. The central result is that the measured pumping rate into the dark state is linear in the CPT beam intensity but approximately an order of magnitude slower than the three-level analytic formula of Eq. (1). The authors attribute this reduction to the multi-level Zeeman structure and support that attribution with a 13-level density-matrix model that, they state, reproduces the measured pumping rate to within about 30% with no fit parameters. They further show that using the pumping rate from the 13-level model in the Ramsey-CPT light-shift formula of Eq. (2) accounts for the measured dependence of the central-fringe light shift on preparation-pulse duration. Finally, the paper measures dark-state decay from Ramsey fringe amplitudes and reports that the decay is dominated by mechanical motion of the atoms out of the probe region, with a power- and Fourier-broadening-free CPT linewidth of 6±20 Hz.","tokens_in":8484,"tokens_out":3281,"duration_ms":35629,"significance":"If the central claim holds, the result is practically important for CPT-based atomic clocks and sensors: it shows that the commonly used three-level formula overestimates the dark-state pumping rate by roughly a factor of ten in a realistic multi-level atom, with direct consequences for light-shift evaluations in Ramsey-CPT spectroscopy. The measurement itself is direct and the comparison between the experimental pumping curve and the 13-level model is a meaningful test. The paper also provides a useful demonstration that in this cold, free-falling ensemble the dark-state coherence lifetime is limited by transit effects rather than by intrinsic decoherence. However, the force of the central claim depends on the 13-level model, whose equations and parameters are not given, and on an unmeasured assumption about the initial Zeeman populations; these gaps currently prevent the claimed no-fit agreement from being independently assessed.","major_comments":[{"comment":"The manuscript states that a 13-level density-matrix model is solved with Lindblad decay operators, but it does not provide the Hamiltonian, the Rabi frequencies for the σ+ and σ− fields, the Zeeman shifts, the spontaneous emission branching ratios, or the numerical integration method. Because the central claim that the multi-level Zeeman manifold explains the factor-of-ten discrepancy rests entirely on this model, the model equations and all parameter values must be included (or supplied as an openly available code/notebook) so that the no-fit agreement can be reproduced and checked.","section":"Section III, 13-level numerical model"},{"comment":"The model assumes that after optical molasses all Zeeman sublevels of F=2 are equally populated and that there are no coherences, justified only by 'Due to the cooling process.' No in-situ measurement, independent model, or sensitivity analysis is provided. Optical molasses is known to produce polarization- and intensity-dependent Zeeman population distributions, and the dark-state pumping dynamics depend on the populations of the coupled mF sublevels. The authors should either measure or bound the initial Zeeman distribution, or demonstrate quantitatively that the predicted pumping rate and build-up shape are insensitive to realistic deviations from the uniform, coherence-free assumption.","section":"Section III, initial conditions of the model"},{"comment":"The paper acknowledges that the pumping process is not purely exponential, yet the experimental 1/e rates in Fig. 3 are obtained by fitting exponentials to curves that visibly deviate from exponential form, and the claimed 'within 30%' agreement is quoted without quantifying this systematic uncertainty. The authors should define exactly how the 1/e rate is extracted from a non-exponential curve, report the fit uncertainties and goodness of fit, and show that the 30% agreement survives when the extraction procedure or the comparison metric is varied.","section":"Section III and Fig. 4, non-exponential pumping"},{"comment":"The correction of the light-shift formula is described only as 'replacing the pumping rate Ω² by the rate calculated from the 13-level model,' but it is not specified how this replacement is implemented in Eq. (2), where α and β both depend on Ω², γ, and δ. In addition, both theoretical curves in Fig. 5 are labeled 'normalized,' but the normalization procedure, the error bars on the experimental points, and the quantitative measure of 'good fit' are not given. The authors should state the modified formula, the normalization, and the fit residuals.","section":"Section IV and Fig. 5, light-shift comparison"},{"comment":"The natural CPT linewidth is reported as 6±20 Hz, which is consistent with a wide range of linewidths and is best interpreted as an upper bound rather than as direct evidence that atomic decoherence is negligible. The conclusion that transit motion dominates is based on the similarity of two decay curves in Fig. 6, but no quantitative model of the probe beam profile, the free-fall trajectory, and the expected transit-time decay is provided. Adding such a model would allow a quantitative comparison of the measured decay rates with the expected mechanical-loss rate, rather than relying on visual similarity.","section":"Section V, claimed negligible decoherence"}],"minor_comments":[{"comment":"The title contains a typographical artifact: 'dark states' appears as 'd ark states' in the manuscript text; this should be corrected.","section":"Title/abstract"},{"comment":"Eq. (2) is written as a proportionality with an unstated constant, and the quantities φ, α, and β are not fully defined in the text (for example, the relationship between φ and the frequency shift of the central fringe). The authors should give the full expression and define all symbols.","section":"Eq. (2)"},{"comment":"The inset legend and caption mention '3-level simulation' and '13-level simulation,' but the data points are identified only by shape; the caption should explicitly state which symbol corresponds to which simulation, and the figure should be legible when rendered at journal size.","section":"Fig. 3 inset"},{"comment":"The procedure for extracting the power-broadening-free and Fourier-broadening-free linewidth from the Voigt fits is only sketched. The authors should specify the Voigt parameters that were held fixed, how the Gaussian contribution was determined, and how the Fourier broadening was subtracted.","section":"Section V, CPT linewidth measurement"},{"comment":"The statement about slower pumping in lin⊥lin Cesium vapor cites a private communication [16] as support. A published reference or a more detailed description of that measurement would make the comparison verifiable.","section":"References"},{"comment":"There are several typographical errors, including 'appartus,' 'demostrated,' 'mechnical,' 'brodening,' and 'shutoff.' The manuscript should be carefully proofread.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is from an experienced NIST group and the experimental data appear to be taken carefully, but the manuscript currently does not provide enough information about the 13-level model or the initial-state assumption for the central claim to be independently verified. I would encourage the editor to require the authors to supply the model equations and parameters, or a public code repository, and to add a sensitivity analysis for the initial Zeeman distribution. The decay-rate conclusion would also be strengthened by a quantitative transit-loss model. These are fixable within the scope of the manuscript, so major revision rather than rejection seems appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the thing you should know: the central measurement is direct and clean. They prepare cold 87Rb atoms, apply sigma+ sigma- CPT light, track the Ramsey fringe amplitude as a function of the preparation pulse duration, and extract a pumping rate that is linear in intensity but about ten times slower than the standard three-level formula. That's a real effect, and it matters for CPT clock light-shift systematics. The paper does a good job presenting the data and making the case that the multi-level Zeeman structure is responsible, mainly via a 13-level density-matrix model that reproduces the measured rate within about 30% with no fitted parameters. The Ramsey light-shift connection is a nice payoff: putting the corrected rate into the theory makes the curve follow the measured shift. The citations look right—they credit the earlier Cs vapor measurement and their own previous light-shift work.\n\nWhere it gets softer: the 13-level model is described only verbally. No equations, no parameters, no sensitivity analysis. That makes the 30% agreement hard to evaluate. The initial-state assumption—equal, incoherent populations of the F=2 Zeeman sublevels after optical molasses—is exactly the kind of thing that can shift the predicted pumping rate. It might be fine, but they don't show it, and they don't show how robust the result is to realistic deviations. The light-shift comparison in Fig. 5 is normalized, so the 'quantitatively explain' claim in the abstract is overreach; they show the corrected shape matches, not that the absolute scale is predicted. And the natural linewidth result of 6 ± 20 Hz is really an upper bound; the text should be worded that way.\n\nNone of this sinks the paper. The core observation is solid and the explanation is physically sensible, with support from earlier Cs work. But it needs revision: give the model details, address the initial-state sensitivity, and either present an absolute light-shift comparison or soften the claim.\n\nWho is this for? Anyone doing CPT clocks, magnetometers, or dark-state work in alkali atoms—they'll want the corrected rate. I'd send it to peer review. Not a breakthrough, but a useful, honest experimental contribution that deserves expert evaluation. My own verdict would be conditional acceptance after the model details and robustness checks are added.","headline":"Solid cold-atom measurement of a tenfold-slower CPT pumping rate, with a plausible but under-specified 13-level model; warrants peer review with requests for model details and sensitivity analysis.","tokens_in":9007,"tokens_out":4624,"would_cite":true,"duration_ms":49845,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.65.-k","42.25.Bs"],"model":"deepseek-v4-flash","headline":"In laser-cooled 87Rb atoms probed with σ+–σ− light, the rate at which atoms are pumped into a coherent dark state is linear in intensity but about ten times slower than the standard three-level formula predicts; the paper shows the gap…","keywords":["coherent population trapping","dark state","optical pumping rate","Zeeman manifold","Ramsey-CPT spectroscopy","light shift","cold rubidium atoms","density-matrix model"],"falsifier":"Pump the atoms into a single Zeeman sublevel, such as $m_F=0$ or $m_F=2$, before applying the CPT pulses and measure the 1/e pumping rate versus intensity; the Zeeman-manifold explanation predicts a rate close to the three-level formula in that prepared state, whereas a tenfold slowdown would point to a different mechanism, such as imperfect polarization or a mis-modeled excited-state hyperfine structure.","tokens_in":8076,"feed_emoji":"⚛️","tokens_out":7650,"duration_ms":78187,"temperature":0.7,"pith_summary":"The paper measures how fast free-falling laser-cooled 87Rb atoms are driven into a coherent dark state by two resonant light fields in the $\\sigma^+-\\sigma^-$ configuration. The measured 1/e pumping rate grows linearly with light intensity, but it is about ten times smaller than the rate given by the standard three-level analytic formula. Using a 13-level density-matrix model that includes the Zeeman sublevels of the D1 transition and Lindblad decay, the authors reproduce the measured rate to within about 30% with no fit parameters. They then show that the slower pumping rate quantitatively explains the observed dependence of the Ramsey-CPT light shift on preparation-pulse duration, and that the dark-state coherence itself is extremely long, with ensemble decay dominated by atoms leaving the probe beam.","feed_headline":"Dark-state pumping measured ten times slower than predicted","feed_subtitle":"Rubidium's Zeeman structure explains the gap; a 13-level model matches the measured pump rate.","key_machinery":"The central object is a 13-level density-matrix simulation of the 87Rb D1 system (ground $F=1$ and $F=2$, excited $F'=2$, all Zeeman sublevels) with spontaneous emission included through Lindblad operators. The dark-state fraction is read from the magnitude of the coherence between the two clock states $|F=1,m_F=0\\rangle$ and $|F=2,m_F=0\\rangle$. This model supplies the replacement pumping rate used in the Ramsey-CPT light-shift formula, and it is tested against Ramsey fringe amplitudes, where the amplitude of the second CPT pulse tracks how much of the ensemble has been pumped into the dark state.","core_discovery":"For coherent population trapping, two resonant fields create a superposition of ground hyperfine states that no longer absorbs light. In the usual three-level picture the pumping rate into that dark state is $\\Gamma_p = (|\\Omega_1|^2 + |\\Omega_2|^2)/(4\\gamma_{\\mathrm{opt}})$, linear in intensity. The paper's central finding is that for 87Rb in the $\\sigma^+-\\sigma^-$ configuration the measured rate is still linear in total intensity but roughly an order of magnitude smaller than this formula. A numerical model of all 13 levels of the D1 transition ($F=1$, $F=2$, $F'=2$ with their Zeeman sublevels) reproduces the measured pumping rate within about 30%, while a numerical three-level model reproduces the analytic formula, showing that the discrepancy is caused by the multi-level Zeeman manifold. When the slower pumping rate is inserted into the Ramsey-CPT light-shift formula, it matches the measured central-fringe shift versus pumping duration. The paper also reports a natural CPT linewidth of $6 \\pm 20$ Hz after removing power and Fourier broadening, and attributes the observed decay of dark-state population to atoms falling out of the interrogation region rather than to atomic decoherence.","pith_inferences":["The roughly tenfold slowdown can be viewed as a geometric factor set by Clebsch–Gordan couplings among the Zeeman sublevels; a testable prediction is that pumping atoms prepared in the clock state $m_F=0$ would recover a rate close to the three-level formula, while atoms prepared in stretched states would pump even more slowly.","Since the model has no fitted parameters, the residual 30% discrepancy is a sensitive probe of the assumed initial state; deliberately preparing different Zeeman population distributions would map how the pumping rate depends on initial conditions and could tighten or falsify the model.","If ensemble decay is purely motional, then Ramsey-CPT contrast in free-fall clocks can be improved more directly by enlarging the probe volume or lowering the atom temperature than by further reducing atomic decoherence.","A direct measurement of the effective pumping coefficient in other alkalis or other polarization configurations could convert the numerical 13-level result into a simple analytic scaling law for CPT light-shift corrections."],"forward_implications":["Evaluations of the light shift in Ramsey-CPT clocks that rely on the three-level formula will overestimate the pumping speed by roughly tenfold in this configuration; using the corrected multi-level rate restores agreement with measured shifts.","The pumping rate remains proportional to intensity, so linear power scaling survives, but the proportionality coefficient must be computed from the full Zeeman manifold rather than from the three-level formula.","Three-level analytic and numerical descriptions are not adequate for $\\sigma^+-\\sigma^-$ pumping in alkali atoms; quantitative work on CPT clocks and sensors needs a model with all coupled sublevels.","Dark states in free-falling cold atoms can have very long coherence: the measured power- and Fourier-broadening-free CPT linewidth is $6 \\pm 20$ Hz.","Similar slower-than-three-level pumping should be expected in other multi-level systems, and the paper cites a cesium vapor measurement with lin $\\perp$ lin light showing the same qualitative behavior."],"supporting_citations":[{"why":"Supplies the three-level pumping-rate formula (Eq. 1) and the light-shift theory that the paper tests.","marker":"[13]"},{"why":"Introduces the Ramsey-CPT interrogation method and reports a cesium vapor measurement with a slower-than-formula pumping rate.","marker":"[12]"},{"why":"Gives the off-resonant light-shift analysis whose agreement with experiment improves when the measured slower pumping rate is used.","marker":"[15]"},{"why":"Provides the Ramsey-CPT central-fringe light-shift formula (Eq. 2) into which the corrected pumping rate is inserted.","marker":"[17]"},{"why":"Reports the previous measurement of the light shift in the same apparatus by locking a clock on the central fringe, consistent with the present results.","marker":"[10]"},{"why":"Identifies the $\\sigma^+-\\sigma^-$ configuration addressing $F'=2$ as the scheme that yields high-contrast CPT resonances.","marker":"[11]"}],"fun_headline_variants":["Dark-state pumping 10x slower than simple 3-level model","Zeeman manifold slows dark-state pumping tenfold","CPT pumping: 13-level model matches, 3-level off by 10x","Dark-state pumping rate: motion limits decay, Zeeman slows pumping"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 13-level model starts from the assumption that optical molasses leaves the five $F=2$ Zeeman sublevels equally populated with no coherences, so the measured 30% agreement rests on that assumed initial distribution.","fun_headline_variants_meta":{"raw":{"variants":["Dark-state pumping 10x slower than simple 3-level model","Zeeman manifold slows dark-state pumping tenfold","CPT pumping: 13-level model matches, 3-level off by 10x","Dark-state pumping rate: motion limits decay, Zeeman slows pumping"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000786,"raw_usage":{"total_tokens":3499,"prompt_tokens":1006,"completion_tokens":2493,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":2417}},"tokens_in":622,"tokens_out":2493,"duration_ms":20627,"temperature":1.0,"reasoning_tokens":2417,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:44:26.263132+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Pump the atoms into a single Zeeman sublevel, such as $m_F=0$ or $m_F=2$, before applying the CPT pulses and measure the 1/e pumping rate versus intensity; the Zeeman-manifold explanation predicts a rate close to the three-level formula in that prepared state, whereas a tenfold slowdown would point to a different mechanism, such as imperfect polarization or a mis-modeled excited-state hyperfine structure.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the three-level pumping-rate formula (Eq. 1) and the light-shift theory that the paper tests."},{"cited_title":"Kulin, B","cited_arxiv_id":null,"evidence_quote":"Introduces the Ramsey-CPT interrogation method and reports a cesium vapor measurement with a slower-than-formula pumping rate."},{"cited_title":"Zanon, S","cited_arxiv_id":null,"evidence_quote":"Gives the off-resonant light-shift analysis whose agreement with experiment improves when the measured slower pumping rate is used."},{"cited_title":"Guerandel, T","cited_arxiv_id":null,"evidence_quote":"Provides the Ramsey-CPT central-fringe light-shift formula (Eq. 2) into which the corrected pumping rate is inserted."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the $\\sigma^+-\\sigma^-$ configuration addressing $F'=2$ as the scheme that yields high-contrast CPT resonances."}],"review_version":1}