Pith. sign in

REVIEW 3 major objections 3 minor 40 references

Laser cooling trapped-ion crystal modes beyond the Lamb-Dicke regime

T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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+.

desk verdict 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. read the letter →

arxiv 2411.18818 v1 pith:Z6QDYVPR submitted 2024-11-27 physics.atom-ph quant-ph

John P. Bartolotta , Brian Estey , Michael Foss-Feig , David Hayes , Christopher N. Gilbreth This is my paper · ORCID
classification physics.atom-phquant-ph PACS 32.80.Pj37.10.Ty
keywords lasercoolingtrappedionsLamb-DickeregimesemiclassicalapproximationEITcapturerangebarium-138rate
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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)$.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 3 minor

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.

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 (3)
  1. [Section IV and Appendix C] 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.
  2. [Section IV, Fig. 4(a)] 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.
  3. [Section III and Appendix C] 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.
minor comments (3)
  1. [Main text, Eq. (25) and Appendix B] 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.
  2. [Appendix A] There is a typo in the paragraph after Eq. (A15): "Hamiltnoian" should be "Hamiltonian." Also, in Section VII, "interemediate-temperature" should be "intermediate-temperature."
  3. [Section V and Appendix D] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: PACMAN cooling rates are computed from the semiclassical dynamics and benchmarked against independent fully-quantum and experimental data; the acknowledged coherent-state limitation is a correctness caveat, not a circular reduction.

full rationale

The paper's central claim is that its semiclassical PACMAN method predicts energy-dependent cooling rates beyond the Lamb-Dicke regime. The derivation starts from Ehrenfest's theorem and two explicitly stated approximations (Eqs. A10-A11), then defines the cooling rate as WSC,mu(n) = -R_mu(n)/n_mu (Eq. 15), where R_mu is the time-averaged semiclassical energy-change rate computed from the laser force and coherent-state parametrized motion. This is a model construction, not a fit to the quantities it later predicts. The benchmark in Section III compares WSC to a fully-quantum master equation initialized independently in a phase-averaged coherent state; agreement at epsilon > 0.1 is a nontrivial dynamical result, not an identity, since the FQ model includes spontaneous-recoil dissipation to second order while PACMAN omits it. The experimental comparison in Section IV is also out-of-sample: laser intensities, detunings, and linewidths are calibrated by fitting the fluorescence lineshape (Fig. 4b), a separate observable from the re-cooling time tau plotted in Fig. 4a, and no cooling-rate or tau data are used to adjust the model. The paper's own Appendix C explicitly concedes that the coherent-state parametrization is not universally valid and that the experimental high-energy case may thermalize; this is an honest correctness limitation and an argument against overinterpreting Fig. 4a, but it is not circularity, because the prediction is not defined in terms of the experimental outcome. The only self-citation appears in the speculative concluding remark about entropy removal via spontaneous emission and is not load-bearing. No equation reduces a predicted quantity to a fitted parameter, no uniqueness result is imported from the authors' prior work, and no known result is merely renamed. The derivation chain is therefore self-contained, with external quantum simulations and experiment serving as genuine benchmarks.

Assumptions & free parameters 1 free parameters · 6 assumptions · 0 invented entities

No new physical entities are introduced. The ledger contains no per-mode fitted scale; the only data-fitted numbers are the six experimental laser parameters calibrated to the EIT fluorescence lineshape. The model's burden sits in the stated semiclassical and fixed-energy approximations, which are transparent and partially validated by the fully-quantum comparison.

free parameters (1)
  • Ba+ experiment laser calibration set (I493, I650, Delta493, Delta650, gamma493, gamma650) = I493 = 0.784 mW/mm^2, I650 = 0.906 mW/mm^2, Delta493/2pi = 25.35 MHz, Delta650/2pi = 23.94 MHz, gamma493/2pi = 0 kHz…
    These six parameters are fitted to the experimental EIT fluorescence lineshape in Fig. 4(b) before computing the cooling rate. The cooling-time comparison is therefore a calibrated prediction, not a fully parameter-free one.
assumptions (6)
  • domain assumption Semiclassical factorization and classical trajectory approximation (Eqs. A10-A11)
    Internal and motional operator expectations are factorized and motional operators are replaced by their classical expectation values. This is the central load-bearing approximation and is stated explicitly in Appendix A.
  • domain assumption Rotating-wave approximation and neglect of recoil, inter-ion internal coherences, and reabsorption of spontaneously emitted photons (Section II)
    These omissions simplify the equations of motion. The paper acknowledges that spontaneous-recoil diffusion is not tracked, which is acceptable for rate predictions but limits the model near steady state.
  • domain assumption PACMAN fixed-energy averaging assumes the laser cooling timescale is long compared to the secular motion and that energy loss per cycle is small (Eq. 13)
    The rate is obtained by averaging the instantaneous power over a long time at fixed total mechanical energy. This requires the cooling dynamics to evolve much more slowly than the trap period.
  • domain assumption Harmonic normal-mode decomposition and coherent-state energy parameterization E = n*omega (Eqs. 8-11)
    The method as presented relies on a harmonic approximation to the trap plus Coulomb potential, so trap and Coulomb nonlinearities are assumed negligible at the energy scales studied.
  • domain assumption Adiabatic elimination of the excited state, dynamic steady-state truncation at finite qmax, and large-detuning Taylor expansion in Appendix B (Eqs. B4-B16, B19)
    The analytic capture-range and cooling-rate formulas are derived with these additional approximations, which are standard but not rigorously controlled at all energies.
  • domain assumption P-representation and Fokker-Planck evolution for phase-averaged distributions with diffusive heating neglected (Appendix D, Eq. D7)
    The evolution of arbitrary phase-averaged motional distributions is obtained by inserting the semiclassical rates into a Fokker-Planck equation, omitting diffusion terms.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Laser cooling trapped-ion crystal modes beyond the Lamb-Dicke regime." pith.science (2026). https://pith.science/paper/Z6QDYVPR

@misc{pith2026241118818,
  author       = {Pith},
  title        = {Pith review of: Laser cooling trapped-ion crystal modes beyond the Lamb-Dicke regime},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Z6QDYVPR}},
  note         = {Machine review of arXiv:2411.18818}
}
abstract

Laser cooling methods for trapped ions are most commonly studied at low energies, i.e., in the Lamb-Dicke regime. However, ions in experiments are often excited to higher energies for which the Lamb-Dicke approximation breaks down. Here we construct a non-perturbative, semiclassical method for predicting the energy-dependent cooling dynamics of trapped-ion crystals with potentially many internal levels and motional modes beyond the Lamb-Dicke regime. This method allows accurate and efficient modeling of a variety of interesting phenomena, such as the breakdown of EIT cooling at high energies and the simultaneous cooling of multiple high-temperature modes. We compare its predictions both to fully-quantum simulations and to experimental data for a broadband EIT cooling method on a Raman $S$-$D$ transition in $^{138}$Ba${^+}$. We find the method can accurately predict cooling rates over a wide range of energies relevant to trapped ion experiments. Our method complements fully quantum models by allowing for fast and accurate predictions of laser-cooling dynamics at much higher energy scales.

Figures

Figures reproduced from arXiv: 2411.18818 by the authors.

Figure 1
Figure 1. FIG. 1. Left [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Cooling rates [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Left [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: FIG. 5. Cooling rates [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. EIT cooling performance as a function of energy. Common parameters are: [PITH_FULL_IMAGE:figures/full_fig_p013_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Characterizing the quantum motional distribution [PITH_FULL_IMAGE:figures/full_fig_p014_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references · 30 canonical work pages

  1. [1]

    J. I. Cirac and P. Zoller, Phys. Rev. Lett.74, 4091 (1995)

    1995

  2. [2]

    S. M. Brewer, J.-S. Chen, A. M. Hankin, E. R. Clements, C. W. Chou, D. J. Wineland, D. B. Hume, and D. R. Leibrandt, Phys. Rev. Lett. 123, 033201 (2019)

    2019

  3. [3]

    Eschner, G

    J. Eschner, G. Morigi, F. Schmidt-Kaler, and R. Blatt, J. Opt. Soc. Am. B 20, 1003 (2003)

    work page 2003

  4. [4]

    R. D. Delaney, L. R. Sletten, M. J. Cich, B. Estey, M. I. Fabrikant, D. Hayes, I. M. Hoffman, J. Hostet- ter, C. Langer, S. A. Moses, A. R. Perry, T. A. Peterson, A. Schaffer, C. Volin, G. Vittorini, and W. C. Burton, Phys. Rev. X 14, 041028 (2024)

    work page 2024

  5. [5]

    A. M. Hankin, E. R. Clements, Y. Huang, S. M. Brewer, J.-S. Chen, C. W. Chou, D. B. Hume, and D. R. Leibrandt, Phys. Rev. A 100, 033419 (2019)

    work page 2019

  6. [6]

    M. W. van Mourik, P. Hrmo, L. Gerster, B. Wilhelm, R. Blatt, P. Schindler, and T. Monz, Phys. Rev. A 105, 033101 (2022)

    work page 2022

  7. [7]

    Stenholm, Rev

    S. Stenholm, Rev. Mod. Phys. 58, 699 (1986)

    1986

  8. [8]

    H. J. Metcalf and P. v. d. Straten, Laser cooling and trapping (Springer, 2002)

    work page 2002

Show all 40 references
  1. [9]

    J. H. Wesenberg, R. J. Epstein, D. Leibfried, R. B. Blakestad, J. Britton, J. P. Home, W. M. Itano, J. D. Jost, E. Knill, C. Langer, et al., Phys. Rev. A 76, 053416 (2007)

    2007

  2. [10]

    Dalibard and C

    J. Dalibard and C. Cohen-Tannoudji, Journal of Physics B: Atomic and Molecular Physics 18, 1661 (1985)

    1985

  3. [11]

    Molmer, Journal of Physics B: Atomic, Molecular and Optical Physics 27, 1889 (1994)

    K. Molmer, Journal of Physics B: Atomic, Molecular and Optical Physics 27, 1889 (1994). 9

    1994

  4. [12]

    Morigi, J

    G. Morigi, J. Eschner, J. I. Cirac, and P. Zoller, Phys. Rev. A 59, 3797 (1999)

    1999

  5. [13]

    Morigi and J

    G. Morigi and J. Eschner, Phys. Rev. A 64, 063407 (2001)

    2001

  6. [14]

    Rabl, Phys

    P. Rabl, Phys. Rev. B 82, 165320 (2010)

    2010

  7. [15]

    Lau and M

    H.-K. Lau and M. B. Plenio, Phys. Rev. B 94, 054305 (2016)

    2016

  8. [16]

    J. I. Cirac, R. Blatt, P. Zoller, and W. D. Phillips, Phys. Rev. A 46, 2668 (1992)

    1992

  9. [17]

    Morigi, J

    G. Morigi, J. Eschner, and C. H. Keitel, Phys. Rev. Lett. 85, 4458 (2000)

    2000

  10. [18]

    A. J. Rasmusson, M. D’Onofrio, Y. Xie, J. Cui, and P. Richerme, Phys. Rev. A 104, 043108 (2021)

    2021

  11. [19]

    R. J. Cook, Phys. Rev. A 20, 224 (1979)

    1979

  12. [20]

    Allevi, M

    A. Allevi, M. Bondani, P. Marian, T. A. Marian, and S. Olivares, J. Opt. Soc. Am. B 30, 2621 (2013)

    2013

  13. [21]

    D. J. Wineland and W. M. Itano, Phys. Rev. A 20, 1521 (1979)

    1979

  14. [22]

    Y. M. Hu, W. L. Yang, Y. Y. Xu, F. Zhou, L. Chen, K. L. Gao, M. Feng, and C. Lee, New Journal of Physics 13, 053037 (2011)

    2011

  15. [23]

    E. Peik, J. Abel, T. Becker, J. von Zanthier, and H. Walther, Phys. Rev. A 60, 439 (1999)

    1999

  16. [24]

    D. M. Meekhof, C. Monroe, B. E. King, W. M. Itano, and D. J. Wineland, Phys. Rev. Lett. 76, 1796 (1996)

    1996

  17. [25]

    K. V. Kepesidis, S. D. Bennett, S. Portolan, M. D. Lukin, and P. Rabl, Phys. Rev. B 88, 064105 (2013)

    2013

  18. [26]

    Carmichael, Statistical Methods in Quantum Optics 1: Master Equations and Fokker-Planck Equations , Physics and astronomy online library (Springer, 1999)

    H. Carmichael, Statistical Methods in Quantum Optics 1: Master Equations and Fokker-Planck Equations , Physics and astronomy online library (Springer, 1999)

    1999

  19. [27]

    Home, Entanglement of Two Trapped-Ion Spin Qubits, Ph.D

    J. Home, Entanglement of Two Trapped-Ion Spin Qubits, Ph.D. thesis, University of Oxford (2006)

    2006

  20. [28]

    Leibfried, R

    D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. Mod. Phys. 75, 281 (2003)

    2003

  21. [29]

    S. J. van Enk and G. Nienhuis, Phys. Rev. A 46, 1438 (1992)

    1992

  22. [30]

    J. P. Bartolotta, S. B. J¨ ager, J. T. Reilly, M. A. Norcia, J. K. Thompson, G. Smith, and M. J. Holland, Phys. Rev. Res. 4, 013218 (2022)

    2022

  23. [31]

    Reiter and A

    F. Reiter and A. S. Sørensen, Phys. Rev. A 85, 032111 (2012)

    2012

  24. [32]

    Marquet, F

    C. Marquet, F. Schmidt-Kaler, and D. F. V. James, Ap- plied Physics B 76, 199 (2003)

    2003

  25. [33]

    P. He, P. M. Tengdin, D. Z. Anderson, A. M. Rey, and M. Holland, Phys. Rev. A 95, 053403 (2017)

    2017

  26. [34]

    Mølmer, Y

    K. Mølmer, Y. Castin, and J. Dalibard, J. Opt. Soc. Am. B 10, 524 (1993)

    1993

  27. [35]

    M. R. N. Surajit Sen, Tushar Kanti Dey and G. Gan- gopadhyay, Journal of Modern Optics 62, 166 (2015), https://doi.org/10.1080/09500340.2014.960019

    2015

  28. [36]

    Morigi, Phys

    G. Morigi, Phys. Rev. A 67, 033402 (2003)

    2003

  29. [37]

    semiclassical approximation

    A. M. J. Oz-Vogt and M. Revzen, Jour- nal of Modern Optics 38, 2339 (1991), https://doi.org/10.1080/09500349114552501. Appendix A: Derivation of the semiclassical model Here we derive the semiclassical model in full detail. We begin with a quantum master equation that describe...

    1991 doi

  30. [38]

    This procedure results in the effective Hamiltonian ˆHeff = X j=1,2 ∆j 1 + |Ωj|2 4∆2 j + Γ2 ! |j⟩ ⟨j| + ΩR 2 |1⟩ ⟨2| + H.c

    Adiabatic elimination To reduce the complexity of the resulting solutions, we adiabatically eliminate the excited state by use of an effective operator formalism [31]. This procedure results in the effective Hamiltonian ˆHeff = X j=1,2 ∆j 1 + |Ωj|2 4∆2 j + Γ2 ! |j⟩ ⟨j| + ΩR 2 ...

  31. [39]

    (B14) In our case, S(t) = ( ρ11, Re[ρ12], Im[ρ12])T , and A(t) is a matrix which couples the elements of S(t)

    Dynamic steady state The dynamics of an open quantum system can often be written in terms of optical Bloch equations (OBEs) [35] d dt S(t) = A(t)S(t). (B14) In our case, S(t) = ( ρ11, Re[ρ12], Im[ρ12])T , and A(t) is a matrix which couples the elements of S(t). If the mo- tion...

  32. [40]

    (B13) analytically, we make a few simplifications

    Analytic approximation to the cooling dynamics Before attempting to solve Eq. (B13) analytically, we make a few simplifications. First, we symmetrize the pa- rameters so that ∆ 1 = ∆2 = ∆ (operating near the EIT dark resonance), Ω 1 = Ω2 = Ω (maximizing the cooling rate [36]),...

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.