{"id":"4ef8abc7-d210-45d4-8726-a1406e3623f3","arxiv_id":"2411.18818","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A power-averaged semiclassical model predicts trapped-ion EIT cooling rates and capture ranges beyond the Lamb-Dicke regime, with experimental validation in 138Ba+.","lead":"A group at Quantinuum built a semiclassical simulation method, PACMAN, that predicts laser cooling rates for trapped-ion crystals when the ions are too hot for the usual Lamb-Dicke approximation. They tested it against full quantum simulations and a barium-ion EIT cooling experiment, finding good agreement over a wide energy range.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Experimental validation relies on an unverified coherent-state parametrization of a thermalizing distribution; the paper's own Appendix C admits this can fail.","rationale":"The reader's weakest-assumption identification, the semiclassical factorization in Eqs. (A10)-(A11), is foundational, but the sharper weakness for the central claim is the step from WSC(n) to the experimentally measured thermalized dynamics. The paper's own Appendix C explicitly limits the coherent-state parametrization and notes the experimental regime may thermalize; Sec. IV nevertheless uses a coherent-state trajectory to predict tau. This is a specific, textually grounded gap that can be tested by applying the paper's own thermal-averaging formula, Eq. (32), to the experimental protocol. The fully-quantum benchmark at lower n is real supporting evidence, and the experimental calibration via the fluorescence lineshape is a strength, but neither addresses this distribution-mismatch issue. Because the reader's verdict is already CONDITIONAL and this concern reinforces rather than overturns the need for conditionality, I would keep the verdict unchanged. The internal inconsistency between Eq. (25) and Eq. (B22) is a related secondary issue but does not by itself undermine the core method as directly as the coherent-to-thermal extrapolation does.","tokens_in":19803,"tokens_out":7989,"duration_ms":77624,"concrete_test":"Recompute the predicted re-cooling time tau for the Sec. IV experiment using thermal-state averaging: at each mean occupancy nbar, use W_th(nbar) = (1/nbar) * integral of n WSC(n) P_th(nbar, n) dn from Eq. (32), and integrate dnbar/dt = -W_th(nbar) nbar with the same fitted laser parameters, then compare with the experimental points in Fig. 4(a). If the thermal-averaged tau disagrees with the coherent-state prediction by more than the experimental scatter, or shifts the predicted capture range away from about 3000, the coherent-state assumption is load-bearing; if it reproduces Fig. 4(a) equally well, the concern is resolved. A complementary check is to run the fully-quantum master equation with an initial thermal or TPAC state at the largest feasible n0 (e.g., n0 = 20-50) and compare WFQ with the thermal-averaged WSC as well as with WSC(n0).","verdict_should_be":"UNCHANGED","load_bearing_attack":"PACMAN's cooling rate WSC(n) is defined and benchmarked for coherent (PAC) states: Eqs. (13)-(15) and Fig. 2. The experiment in Sec. IV, however, initializes the axial COM mode at n0 about 1000, and the measured occupancy is extracted by fitting a thermal distribution (Sec. IV). To predict the re-cooling time tau, the paper integrates the coherent-state rate dn/dt = -WSC(n) n, implicitly assuming the motional distribution remains a coherent state throughout cooling. Appendix C checks this only for n0 = 100 in a simplified Lambda-system and explicitly states that the experimental case—much higher energy and richer internal structure—may have resulted in additional heating and that the distribution can be substantially altered. If the distribution thermalizes, Eq. (32) or the Fokker-Planck equation (D7), not the coherent-state trajectory, is the correct semiclassical prediction, and these can differ materially near the EIT capture range where WSC(n) crosses zero. The agreement in Fig. 4(a) is therefore not yet a clean test of the claimed energy-dependent cooling rates; without experimental error bars, it cannot rule out a fortuitous match. This is the load-bearing gap in the central claim that PACMAN accurately predicts cooling rates in the high-energy regime.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces PACMAN, a semiclassical method for computing energy-dependent laser-cooling rates of trapped-ion crystal modes beyond the Lamb-Dicke regime. The method evolves the internal quantum master equation along classical trajectories, uses the resulting laser force to compute the average power removed from each harmonic mode, and defines a cooling rate W_SC(n) as the negative of that power divided by the coherent-state occupation n. The authors benchmark W_SC(n) against a fully quantum master-equation simulation for a Lambda-system EIT-cooling scheme, obtaining agreement for epsilon > 0.1, and they apply the method to a broadband EIT-cooling experiment on a 138Ba+–171Yb+ crystal. They further derive an analytic approximation for the high-energy EIT capture range, extend the formalism to phase-averaged (e.g., thermal) motional distributions, and demonstrate multi-mode cooling-rate calculations for a twelve-mode YBBY crystal.","tokens_in":19986,"tokens_out":6405,"duration_ms":61875,"significance":"If the method is valid, it addresses a practical need: efficient prediction of laser-cooling dynamics at energies where the Lamb-Dicke expansion fails and where fully quantum simulations are intractable. The derivation in Appendix A is careful and explicit about the semiclassical approximations, and the comparison with the fully quantum model in Fig. 2 is a genuine test of the cooling rates in that simplified system. The extension to arbitrary phase-averaged distributions via Eq. (28) and the Fokker-Planck formulation in Appendix D are valuable and give the method broader applicability than a single coherent-state trajectory. The analytic capture-range estimate in Appendix B is also a useful design tool. However, the experimental validation in Section IV is not yet a clean test of the energy-dependent cooling rates because the prediction uses a coherent-state trajectory while the data are analyzed with a thermal distribution, and the paper itself concedes in Appendix C that the experimental case may thermalize and require a Fokker-Planck treatment. The experimental comparison also lacks error bars, so the strength of the claimed agreement is difficult to assess.","major_comments":[{"comment":"The central experimental comparison conflates two different motional distributions. PACMAN predicts the time evolution of a coherent-state occupation n(t) via dn/dt = -W_SC(n)n [Eqs. (15) and (30)], while the experimental data are analyzed by fitting a thermal distribution to extract nbar_th(t) and then fitting those thermal occupancies to the exponential form of Eq. (24). If the motional distribution thermalizes during cooling, the appropriate semiclassical prediction is not the coherent trajectory but the thermal rate of Eq. (32), or the Fokker-Planck evolution of Eq. (D7). Appendix C explicitly states that \"the experimental results in Section IV suggest that the modes thermalize before reaching steady state\" and that the high-energy, multi-level experimental case \"may have resulted in some additional heating\" requiring a Fokker-Planck approach. To make Fig. 4(a) a valid test of the energy-dependent W_SC(n), the authors should either compute tau using the thermal-distribution expression of Eq. (32) and compare it with the experimental thermal-fit results, or demonstrate quantitatively that the coherent-state and thermal-state predictions are indistinguishable over the n0 range shown in Fig. 4(a). As presented, the agreement in Fig. 4(a) does not cleanly validate the claimed cooling rates in the high-energy regime.","section":"Section IV and Appendix C"},{"comment":"The experimental comparison lacks uncertainty quantification. The re-cooling time tau is obtained from an exponential fit to only five cooling times satisfying 0.5 < nbar_th(t) < 20, and no error bars or run-to-run variations are reported in Fig. 4(a). Without an estimate of the statistical or systematic uncertainty in tau, the claim of \"strong agreement\" over a broad energy range is not yet supported. The authors should provide error bars on the experimental points, state the number of repeated measurements, and give the fitted exponential parameters or an equivalent measure of the fit quality.","section":"Section IV, Fig. 4(a)"},{"comment":"The validation against the fully quantum model is limited to the simplified Lambda system and to energies up to n = 100 in the continuous simulation; for n > 100, the FQ cooling rates are obtained from a truncated Fock-state window rather than from a full long-time simulation. This is a reasonable practical choice, but it means the agreement in Fig. 2 at the highest energies rests on the assumption that the truncated window faithfully represents the cooling dynamics. Given that the experimental system operates at n0 ~ 1000 with a richer internal level structure, the authors should state more explicitly how the truncated-window FQ calculations were validated at intermediate energies and whether the same window approximation was used for any of the experimental parameters. This would clarify the range of energies for which W_SC(n) is directly benchmarked.","section":"Section III and Appendix C"}],"minor_comments":[{"comment":"The capture-range formula in the main text, ncap approx (Omega^2 - 2 omega^2)/(2 eta^2 omega^2), differs from Eq. (B22), which gives ncap approx (Omega^2 - omega^2)/(2 eta^2 omega^2). Please reconcile the factor of 2 in the numerator and ensure the formula is stated consistently in both places.","section":"Main text, Eq. (25) and Appendix B"},{"comment":"There is a typo in the paragraph after Eq. (A15): \"Hamiltnoian\" should be \"Hamiltonian.\" Also, in Section VII, \"interemediate-temperature\" should be \"intermediate-temperature.\"","section":"Appendix A"},{"comment":"The derivation of Eq. (28) assumes that the P-function evolves by the Fokker-Planck equation (D7) with diffusive heating neglected. This is stated, but the sentence should also note explicitly that the resulting thermal-state prediction of Eq. (32) is therefore not a full quantum prediction and may differ from the coherent-state trajectory even when the mean occupancies coincide.","section":"Section V and Appendix D"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid methods contribution with a clean derivation and a convincing fully quantum benchmark in the simplified Lambda system. The main issue is that the experimental validation, which is the headline evidence for the method's practical accuracy at high energies, is undermined by the coherent-versus-thermal distribution mismatch and by the absence of error bars. Both are fixable: the authors can compute the thermal-distribution prediction using Eq. (32) and compare it with the measured tau values, and they can report experimental uncertainties. If those changes are made, the paper would likely be acceptable. The authors are honest about the limitations in Appendix C, and I do not see a circularity problem in the parameter fitting to the fluorescence lineshape, since the cooling-rate comparison is out-of-sample."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper has a useful core and an honest experimental section, but the headline claim that it accurately predicts cooling rates over a wide energy range is only partly backed by the evidence. The new piece is the fixed-energy power-averaged (PACMAN) formulation of semiclassical cooling, which extends earlier work by Rabl and Lau & Plenio to multi-ion, multi-level crystals and comes with an analytic estimate of the EIT capture range. The derivation in Appendix A is clean, and the comparison to the fully-quantum master equation for the single-mode Lambda system agrees well for epsilon > 0.1. That part is solid.\n\nThe experimental comparison is the most valuable part. They calibrate laser parameters on the fluorescence lineshape and then predict re-cooling times for a range of initial excitations, plus the capture range near n ~ 3000. That is a genuine out-of-sample check, and the qualitative agreement is encouraging. But there are two soft spots. First, the experiment extracts thermal occupancies from Rabi flops, while the theory integrates the coherent-state rate dn/dt = -WSC(n) n. If the distribution thermalizes during cooling, as Appendix C says it might for the experimental conditions, the correct prediction uses the Fokker-Planck equation (D7), not the coherent trajectory. Without experimental error bars on the re-cooling times, I can't tell how tight the agreement really is. The paper should address this directly by applying the thermal-state formalism from Section V to the experimental parameters, or at least discussing how much the distribution changes. Second, the capture-range formula in Eq. (25) differs from Eq. (B22) by a factor of two in the numerator (Omega^2 - 2omega^2 vs Omega^2 - omega^2). That's a minor error, but it should be fixed.\n\nOverall, the method is sound, the limitations are honestly stated, and the work deserves peer review. I'd accept it after revision. The Fokker-Planck check and error bars are necessary before the central claim is fully supported.","headline":"PACMAN is a useful extension of semiclassical cooling theory; the FQ validation is solid, but the experimental agreement is not yet clean without error bars and a distribution-model check.","tokens_in":20571,"tokens_out":4954,"would_cite":true,"duration_ms":40645,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["32.80.Pj","37.10.Ty"],"model":"deepseek-v4-flash","headline":"A semiclassical model called PACMAN predicts trapped-ion laser-cooling rates beyond the Lamb-Dicke regime and reproduces a broadband EIT cooling experiment on 138Ba+.","keywords":["laser cooling","trapped ions","Lamb-Dicke regime","semiclassical approximation","EIT cooling","capture range","barium-138","cooling rate"],"falsifier":"Measure the capture threshold directly: prepare the axial center-of-mass mode of a barium-ytterbium crystal in coherent states at occupancies spanning the predicted zero of $W_{\\mathrm{SC}}(n)$ (near $n \\approx 3000$ for the paper's parameters), cool for a fixed short time, and record whether the energy decreases or increases. With the Raman Rabi frequency held fixed, the model predicts the threshold should scale like $n_{\\mathrm{cap}} \\propto \\Delta$; varying the 493 and 650 nm laser detunings over a factor of a few and checking that the heating boundary follows that scaling would settle the claim. The existing data already constrain the boundary to lie between the experimental failure at $n_0 \\geq 2000$ and the predicted $n_{\\mathrm{cap}} \\approx 3000$, so a scan with fine $n_0$ steps would resolve the remaining gap.","tokens_in":19561,"feed_emoji":"⚛️","tokens_out":19177,"duration_ms":149157,"temperature":0.7,"pith_summary":"The paper aims to establish that a semiclassical method, called PACMAN, can predict laser-cooling rates for trapped-ion crystals at energies where the standard Lamb-Dicke expansion breaks down, and that these predictions are quantitatively reliable. The authors derive an energy-dependent cooling rate $W_{\\mathrm{SC}}(n)$ by treating each vibrational mode as a classical oscillator at coherent-state energy $n\\omega$, evolving the ions' internal quantum states under the Doppler-shifted laser fields, and time-averaging the resulting power. They show this rate matches a fully quantum master equation wherever that simulation is tractable (for $\\epsilon = \\eta\\sqrt{2n+1} \\gtrsim 0.1$), and matches a broadband EIT cooling experiment on $^{138}\\mathrm{Ba}^{+}$ across a wide range of initial energies: predicted re-cooling times agree with measured ones, and the predicted capture range, the energy above which cooling becomes heating, sits near $n \\approx 3000$, consistent with the experiment's failure to re-cool excitations of $n_0 \\geq 2000$. This matters because ions in real experiments are routinely excited to high energies by transport and collisions, and at those energies the standard expansion fails while fully quantum simulation becomes numerically intractable, leaving experiments without a practical prediction of whether their cooling protocol will work.","feed_headline":"Semiclassical model predicts ion cooling beyond Lamb-Dicke regime","feed_subtitle":"Matches full quantum simulation and a 138Ba+ EIT experiment, predicting where cooling flips to heating.","key_machinery":"The load-bearing object is PACMAN, the power-averaged cooling method for analyzing $\\bar{n}$: a prescription for an energy-dependent, mode-resolved cooling rate. Its two semiclassical replacements are the factorization of internal-motional correlations and the substitution of positions and momenta by their expectation values, which turn Ehrenfest's theorem into Newton's equations $d\\mathbf{r}_j/dt = \\mathbf{p}_j/m_j$ and $d\\mathbf{p}_j/dt = \\mathbf{F}^{\\mathrm{laser}}_j + \\mathbf{F}^{\\mathrm{trap}}_j + \\mathbf{F}^{\\mathrm{Coulomb}}_j$, coupled to quantum master equations for each ion's internal state through the laser phase $\\phi^{(l,j)}_{\\alpha\\beta}(\\mathbf{r}_j,t) = \\int_0^t [\\Delta^{(l)}_{\\alpha\\beta} - \\mathbf{k}_l \\cdot \\mathbf{v}_j(t')]\\,dt'$. The motion of each normal mode is fixed to a classical coherent-state trajectory at energy $n_\\mu\\omega_\\mu$, and the instantaneous power $\\mathbf{F}^{\\mathrm{laser}}_\\mu(t)\\cdot\\mathbf{v}_\\mu(t)$ is averaged over many trap periods to define $W_{\\mathrm{SC},\\mu}(\\mathbf{n})$. Two supporting pieces carry the benchmarking and the generalization: the fully quantum and Lamb-Dicke master-equation models whose regime split is indexed by $\\epsilon = \\eta\\sqrt{2n+1}$, and a P-function formalism whose moment equation $\\frac{d}{dt}\\langle \\hat{n}_\\mu^k\\rangle = -k\\pi^N \\int d^N n\\, n_\\mu^k W_{\\mathrm{SC},\\mu}(\\mathbf{n}) P(\\mathbf{n})$ extends the coherent-state rates to thermal, PAC, and other phase-averaged distributions. An analytic reduction in Appendix B identifies the cooling as a competition between coherent Raman forces and dissipative Doppler forces, yielding the capture-range estimate $n_{\\mathrm{cap}} \\approx (\\Omega^2 - \\omega^2)/(2\\eta^2\\omega^2)$.","core_discovery":"The central claim, stated on the paper's own terms, is that a non-perturbative semiclassical prescription computes the average laser-cooling rate of every crystal mode as a function of the mode energies, and that this function is accurate deep into the regime where the Lamb-Dicke expansion fails. The prescription factorizes internal and motional expectation values ($\\langle \\hat{\\sigma} e^{i \\mathbf{k}\\cdot\\hat{\\mathbf{r}}}\\rangle \\approx \\langle\\hat{\\sigma}\\rangle \\langle e^{i\\mathbf{k}\\cdot\\hat{\\mathbf{r}}}\\rangle$), replaces the motional operators by their classical averages, evolves the internal density matrix under the resulting Doppler-shifted laser phases, and defines $W_{\\mathrm{SC},\\mu}(\\mathbf{n}) = -R_\\mu(\\mathbf{n})/n_\\mu$ from the time-averaged power. The paper's evidence for the claim has two anchors. Numerically, the semiclassical cooling rate reproduces the fully quantum cooling rate for $\\epsilon > 0.1$ while the Lamb-Dicke model reproduces it only for $\\epsilon < 0.1$. Experimentally, PACMAN reproduces the measured re-cooling time of a coherently excited axial mode of a barium-ytterbium crystal over more than three orders of magnitude in $n_0$, and its zero crossing predicts the EIT capture range near $n \\approx 3000$, matching the observed inability to cool excitations of $n_0 \\geq 2000$ back toward the ground state. A corollary the paper draws is that cooling rates beyond the Lamb-Dicke regime are generally lower than the exponential Lamb-Dicke rate and can become negative (runaway heating), so pre-cooling toward Doppler temperatures is required for EIT cooling to succeed.","pith_inferences":["The paper computes the capture range for a pure coherent state; feeding the same $W_{\\mathrm{SC}}$ through Eq. (28) for a thermal distribution should shift the effective heating threshold to lower average energy, because thermal states place some probability above the zero crossing. Experiments starting from Doppler-cooled thermal states rather than coherent kicks could test this directly.","Treating $W_{\\mathrm{SC}}(n)$ as a design objective suggests an optimization rule the paper does not spell out: choose detunings and Rabi frequencies to maximize the energy removed along the actual cooling trajectory, $\\int W_{\\mathrm{SC}}(n(t))\\,dt$, rather than the zero-temperature rate, which would likely favor different parameters for intermediate-temperature operation.","The semiclassical factorization is checked against a fully quantum simulation in only one case (Appendix C); repeating that fidelity analysis across ion masses, trap frequencies, and laser geometries could map where neglected recoil diffusion breaks the method, notably for heavy coolant ions or long cooling times where the distribution visibly thermalizes, as the experiment hints.","Because the method handles arbitrary internal level structures, capture-range engineering could be applied to sympathetic cooling of heteronuclear crystals, choosing the coolant transition and detuning so that $n_{\\mathrm{cap}}$ exceeds the energy left in poorly cooled modes that have small Lamb-Dicke parameters."],"forward_implications":["EIT cooling has a finite energy capture range: a mode below $n_{\\mathrm{cap}}$ is cooled, but a mode above it is heated, so sufficient pre-cooling, toward Doppler temperature, is a precondition for EIT cooling to reach the ground state.","Cooling rates degrade from the exponential Lamb-Dicke prediction already for $n \\gtrsim 10$; the semiclassical rate surface, not the zero-temperature rate, is the right input for designing and optimizing cooling protocols at intermediate temperatures.","The cooling rate of any one mode depends on the energies of all the other modes through the shared laser forces, so a mode cools much faster when it alone carries the energy than when the same energy is spread across all modes.","For thermal or other phase-averaged initial states, the moment equation (28) lets one track the moments of the motional distribution through the semiclassical rates, giving a direct route to cooling dynamics from a Doppler-temperature initial state.","The capture range can be extended to arbitrarily high energies by increasing the laser detuning while holding the Raman Rabi frequency fixed, at the cost of more laser power and slower intermediate-energy cooling."],"supporting_citations":[{"why":"Supplies the Lamb-Dicke-regime rate-equation model whose exponential rate W is the n-to-0 limit that PACMAN must reproduce and one of the two benchmarks.","marker":"[16]"},{"why":"The EIT-cooling scheme and dark state on which both modeled and experimental protocols are built.","marker":"[17]"},{"why":"The semiclassical factorization of internal and motional expectation values and the mode-energy dependence that PACMAN adopts.","marker":"[14]"},{"why":"The high-energy EIT heating picture and the Fokker-Planck and dynamic-steady-state tools used in the appendices and moment formalism.","marker":"[15]"},{"why":"Couples internal and motional dynamics through the laser phases witnessed by the ion, the backbone of PACMAN's coupled equations.","marker":"[9]"},{"why":"Defines the phase-averaged coherent (PAC) state used as the initial motional state for the fully quantum comparison and in the P-function formalism.","marker":"[20]"},{"why":"Provides the Glauber-Sudarshan P representation and Fokker-Planck framework used to extend cooling rates to thermal and other phase-averaged distributions.","marker":"[26]"},{"why":"The carrier, red, and blue sideband Rabi-flop thermometry used to extract average occupancies in the cooling experiment.","marker":"[24]"}],"fun_headline_variants":["Semiclassical shortcut predicts ion cooling beyond Lamb-Dicke","Ion-cooling model works beyond Lamb-Dicke limit","Beyond Lamb-Dicke: semiclassical cooling rates hit the mark","Semiclassical method predicts ion cooling at high energies","Cooling rates beyond Lamb-Dicke? Semiclassical model says yes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method replaces each ion's quantum motion with its average classical position and velocity, and neglects the random recoil kicks of spontaneous emission and any entanglement between motion and internal state; if the motional wavepacket ever spreads to a significant fraction of the cooling wavelength, the predicted rates, including the capture range, could drift from the true quantum result, and the paper demonstrates that the packet stays narrow in only one tested case.","fun_headline_variants_meta":{"raw":{"variants":["Semiclassical shortcut predicts ion cooling beyond Lamb-Dicke","Ion-cooling model works beyond Lamb-Dicke limit","Beyond Lamb-Dicke: semiclassical cooling rates hit the mark","Semiclassical method predicts ion cooling at high energies","Cooling rates beyond Lamb-Dicke? Semiclassical model says yes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00077,"raw_usage":{"total_tokens":3484,"prompt_tokens":1092,"completion_tokens":2392,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":708,"completion_tokens_details":{"reasoning_tokens":2318}},"tokens_in":708,"tokens_out":2392,"duration_ms":15933,"temperature":1.0,"reasoning_tokens":2318,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:51:24.248613+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the capture threshold directly: prepare the axial center-of-mass mode of a barium-ytterbium crystal in coherent states at occupancies spanning the predicted zero of $W_{\\mathrm{SC}}(n)$ (near $n \\approx 3000$ for the paper's parameters), cool for a fixed short time, and record whether the energy decreases or increases. With the Raman Rabi frequency held fixed, the model predicts the threshold should scale like $n_{\\mathrm{cap}} \\propto \\Delta$; varying the 493 and 650 nm laser detunings over a factor of a few and checking that the heating boundary follows that scaling would settle the claim. The existing data already constrain the boundary to lie between the experimental failure at $n_0 \\geq 2000$ and the predicted $n_{\\mathrm{cap}} \\approx 3000$, so a scan with fine $n_0$ steps would resolve the remaining gap.","supporting_citations":[{"cited_title":"Morigi, J","cited_arxiv_id":null,"evidence_quote":"The EIT-cooling scheme and dark state on which both modeled and experimental protocols are built."},{"cited_title":"Rabl, Phys","cited_arxiv_id":null,"evidence_quote":"The semiclassical factorization of internal and motional expectation values and the mode-energy dependence that PACMAN adopts."},{"cited_title":"Lau and M","cited_arxiv_id":null,"evidence_quote":"The high-energy EIT heating picture and the Fokker-Planck and dynamic-steady-state tools used in the appendices and moment formalism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Couples internal and motional dynamics through the laser phases witnessed by the ion, the backbone of PACMAN's coupled equations."},{"cited_title":"Allevi, M","cited_arxiv_id":null,"evidence_quote":"Defines the phase-averaged coherent (PAC) state used as the initial motional state for the fully quantum comparison and in the P-function formalism."}],"review_version":1}