Optimizing Doppler laser cooling protocols for quantum sensing with 3D ion crystals in a Penning trap
Pith reviewed 2026-05-21 13:06 UTC · model grok-4.3
The pith
Numerical simulations reveal enhanced laser cooling pathways for large 3D ion crystals in Penning traps by adding axial components to E×B modes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We develop a powerful numerical framework to simulate laser cooling of up to 10^5 ions stored in a Penning trap. We apply this framework to characterize and optimize the cooling of ellipsoidal 3D crystals. We document new pathways to enhanced cooling based on the addition of an axial component to the potential energy-dominated E×B modes. Furthermore, we observe greatly enhanced cooling of the perpendicular kinetic energy to below 1 mK in prolate ion crystals, enabling a simplified cooling beam setup for such crystals. We propose specific values of trap and laser beam parameters which lead to optimal cooling in a variety of examples.
What carries the argument
Numerical simulation framework for laser cooling dynamics of up to 10^5 ions under combined trap potentials, E×B drifts, and Doppler forces, with added axial perturbations to the modes.
Load-bearing premise
The numerical model correctly captures the coupling between the E×B modes, axial perturbations, and laser cooling forces for large ion numbers without introducing artifacts that alter the predicted temperature reductions or optimal parameters.
What would settle it
An experiment with a prolate 3D ion crystal using the proposed optimal trap and laser parameters that measures perpendicular kinetic energy well above 1 mK would show the predicted cooling enhancement does not occur.
Figures
read the original abstract
Large, 3D trapped ion crystals offer improved sensitivity in quantum sensing protocols, and are expected to be implemented as platforms in near-future experiments. However, numerical techniques used to study the laser cooling of such crystals are inefficient as the number of ions, $N$, in the crystal increases. Here we develop a powerful numerical framework to simulate laser cooling of up to $10^5$ ions stored in a Penning trap. We apply this framework to characterize and optimize the cooling of ellipsoidal 3D crystals. We document new pathways to enhanced cooling based on the addition of an axial component to the potential energy-dominated $\boldsymbol{E}\times\boldsymbol{B}$ modes. Furthermore, we observe greatly enhanced cooling of the perpendicular kinetic energy to below 1 mK in prolate ion crystals, enabling a simplified cooling beam setup for such crystals. We propose specific values of trap and laser beam parameters which lead to optimal cooling in a variety of examples. This work illustrates the feasibility of preparing large 3D crystals for high-sensitivity quantum science protocols, motivating their use in future experiments.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a numerical framework to simulate Doppler laser cooling of up to 10^5 ions in ellipsoidal 3D crystals confined in a Penning trap. It characterizes cooling dynamics, introduces axial components to the potential-energy-dominated E×B modes as a new pathway for enhanced cooling, and reports that prolate crystals achieve perpendicular kinetic energies below 1 mK, enabling simplified beam geometries. Specific trap and laser parameters are proposed for optimal performance across examples.
Significance. If the large-N simulations are accurate, the work would be significant for quantum sensing: it shows that 3D ion crystals can be cooled to temperatures suitable for high-sensitivity protocols with reduced experimental complexity. The framework's ability to reach N=10^5 is a technical advance over prior methods, and the prolate-crystal result offers a concrete, falsifiable prediction for future experiments.
major comments (2)
- [Numerical Methods / Simulation Framework] The manuscript provides no benchmarks of the large-N integrator against small-N exact solutions or established literature results for Penning-trap Doppler cooling (e.g., comparisons of equilibrium temperatures or mode spectra for N<100). This validation is load-bearing for the central claim of sub-1 mK perpendicular cooling, because any under-resolution of the Doppler shift, cyclotron frequency, or axial-radial coupling could produce spurious temperature reductions.
- [Results for Prolate Crystals] The reported temperature reductions in prolate crystals (below 1 mK) are presented without accompanying convergence tests in particle number, timestep, or mean-field approximations. Without these, it is unclear whether the enhanced cooling arises from the added axial E×B perturbation or from numerical artifacts at N=10^5.
minor comments (2)
- [Introduction / Model] Notation for the axial perturbation strength and laser detuning should be defined explicitly in the first appearance rather than assumed from context.
- [Figures] Figure captions for the cooling trajectories should state the exact N, trap parameters, and averaging procedure used to extract the final temperatures.
Simulated Author's Rebuttal
We thank the referee for the careful review and constructive comments on our manuscript. We address the major concerns point by point below and will revise the manuscript to incorporate additional validation and convergence tests.
read point-by-point responses
-
Referee: [Numerical Methods / Simulation Framework] The manuscript provides no benchmarks of the large-N integrator against small-N exact solutions or established literature results for Penning-trap Doppler cooling (e.g., comparisons of equilibrium temperatures or mode spectra for N<100). This validation is load-bearing for the central claim of sub-1 mK perpendicular cooling, because any under-resolution of the Doppler shift, cyclotron frequency, or axial-radial coupling could produce spurious temperature reductions.
Authors: We agree that explicit benchmarks against small-N results are important for validating the integrator. Our framework extends established molecular-dynamics methods for Penning traps, but we did not present direct comparisons in the original submission. In the revised manuscript we will add a dedicated validation subsection showing equilibrium temperatures and mode spectra for N=20–80 ions, compared to literature values and smaller-scale exact integrations. We will also report timestep and frequency-resolution checks to confirm that Doppler shifts and axial-radial couplings are adequately resolved. revision: yes
-
Referee: [Results for Prolate Crystals] The reported temperature reductions in prolate crystals (below 1 mK) are presented without accompanying convergence tests in particle number, timestep, or mean-field approximations. Without these, it is unclear whether the enhanced cooling arises from the added axial E×B perturbation or from numerical artifacts at N=10^5.
Authors: We accept that convergence tests are needed to substantiate the prolate-crystal results. In the revision we will add explicit checks demonstrating that the sub-1 mK perpendicular temperatures remain stable when the timestep is halved, when N is varied between 5×10^4 and 10^5, and under changes to any mean-field approximations. These tests will be placed in the Results section and will show that the cooling improvement is attributable to the axial component of the E×B modes rather than numerical artifacts. revision: yes
Circularity Check
Numerical simulation framework is self-contained with no circular derivation steps
full rationale
The paper develops and applies a numerical framework for forward simulation of laser cooling dynamics in large Penning-trap ion crystals, using established models of E×B modes, axial perturbations, and velocity-dependent laser forces. No load-bearing results are shown to reduce by construction to fitted parameters, self-defined quantities, or self-citation chains within the paper; the reported temperature reductions and optimal parameters emerge from explicit integration of the physical equations rather than internal renaming or fitting loops. While the skeptic correctly flags the need for small-N validation to confirm absence of integrator artifacts, this is a question of model fidelity and convergence, not circularity in the derivation chain itself.
Axiom & Free-Parameter Ledger
free parameters (1)
- axial component strength and laser detuning values
axioms (1)
- domain assumption The E×B mode structure and its coupling to Doppler cooling forces remain accurately represented by the numerical model at large N.
Lean theorems connected to this paper
-
IndisputableMonolith.Foundation.RealityFromDistinctionreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We develop a powerful numerical framework to simulate laser cooling of up to 10^5 ions stored in a Penning trap... We observe greatly enhanced cooling of the perpendicular kinetic energy to below 1 mK in prolate ion crystals.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
J. W. Britton, B. C. Sawyer, A. C. Keith, C.-C. J. Wang, J. K. Freericks, H. Uys, M. J. Biercuk, and J. J. Bollinger, Engineered two-dimensional Ising interactions in a trapped-ion quantum simulator with hundreds of spins, Nature484, 489 (2012)
work page 2012
-
[2]
A. Safavi-Naini, R. J. Lewis-Swan, J. G. Bohnet, M. G¨ arttner, K. A. Gilmore, J. E. Jordan, J. Cohn, J. K. Freericks, A. M. Rey, and J. J. Bollinger, Verification of a many-ion simulator of the Dicke model through slow quenches across a phase transition, Phys. Rev. Lett.121, 040503 (2018)
work page 2018
-
[3]
J. Cohn, A. Safavi-Naini, R. J. Lewis-Swan, J. G. Bohnet, M. G¨ arttner, K. A. Gilmore, J. E. Jordan, A. M. Rey, J. J. Bollinger, and J. K. Freericks, Bang-bang shortcut to adiabaticity in the Dicke model as realized in a Pen- ning trap experiment, New Journal of Physics20, 055013 (2018)
work page 2018
-
[4]
A. Shankar, E. A. Yuzbashyan, V. Gurarie, P. Zoller, J. J. Bollinger, and A. M. Rey, Simulating dynamical phases of chiralp+ipsuperconductors with a trapped ion magnet, PRX Quantum3, 040324 (2022)
work page 2022
- [5]
- [6]
-
[7]
B. Swingle, G. Bentsen, M. Schleier-Smith, and P. Hay- den, Measuring the scrambling of quantum information, Phys. Rev. A94, 040302 (2016)
work page 2016
-
[8]
M. G¨ arttner, J. G. Bohnet, A. Safavi-Naini, M. L. Wall, J. J. Bollinger, and A. M. Rey, Measuring out-of-time- order correlations and multiple quantum spectra in a trapped ion quantum magnet, Nature Phys.13, 781 (2017), arXiv:1608.08938 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[9]
J. G. Bohnet, B. C. Sawyer, J. W. Britton, M. L. Wall, A. M. Rey, M. Foss-Feig, and J. J. Bollinger, Quantum spin dynamics and entanglement generation with hun- dreds of trapped ions, Science352, 1297 (2016)
work page 2016
-
[10]
L. Pezz` e, A. Smerzi, M. K. Oberthaler, R. Schmied, and P. Treutlein, Quantum metrology with nonclassical states of atomic ensembles, Rev. Mod. Phys.90, 035005 (2018)
work page 2018
-
[11]
K. A. Gilmore, J. G. Bohnet, B. C. Sawyer, J. W. Brit- ton, and J. J. Bollinger, Amplitude sensing below the zero-point fluctuations with a two-dimensional trapped- ion mechanical oscillator, Phys. Rev. Lett.118, 263602 (2017)
work page 2017
-
[12]
M. Affolter, K. A. Gilmore, J. E. Jordan, and J. J. Bollinger, Phase-coherent sensing of the center-of-mass motion of trapped-ion crystals, Phys. Rev. A102, 052609 (2020)
work page 2020
-
[13]
K. A. Gilmore, M. Affolter, R. J. Lewis-Swan, D. Barber- ena, E. Jordan, A. M. Rey, and J. J. Bollinger, Quantum- enhanced sensing of displacements and electric fields with two-dimensional trapped-ion crystals, Science373, 673 (2021)
work page 2021
-
[14]
F. Toscano, D. A. R. Dalvit, L. Davidovich, and W. H. Zurek, Sub-Planck phase-space structures and Heisenberg-limited measurements, Phys. Rev. A73, 023803 (2006)
work page 2006
-
[15]
R. N. Wolf, J. H. Pham, J. Y. Z. Jee, A. Rischka, and M. J. Biercuk, Efficient site-resolved imaging and spin- state detection in dynamic two-dimensional ion crystals, Phys. Rev. Appl.21, 054067 (2024)
work page 2024
-
[16]
B. J. McMahon, K. R. Brown, C. D. Herold, and B. C. Sawyer, Individual-ion addressing and readout in a Pen- ning trap, Phys. Rev. Lett.133, 173201 (2024)
work page 2024
-
[17]
D. Kiesenhofer, H. Hainzer, A. Zhdanov, P. C. Holz, M. Bock, T. Ollikainen, and C. F. Roos, Controlling two-dimensional Coulomb crystals of more than 100 ions in a monolithic radio-frequency trap, PRX Quantum4, 020317 (2023). 16
work page 2023
- [18]
-
[19]
S. Hawaldar, P. Shahi, A. L. Carter, A. M. Rey, J. J. Bollinger, and A. Shankar, Bilayer crystals of trapped ions for quantum information processing, Phys. Rev. X 14, 031030 (2024)
work page 2024
-
[20]
H. N. Hausser, J. Keller, T. Nordmann, N. M. Bhatt, J. Kiethe, H. Liu, I. M. Richter, M. von Boehn, J. Rahm, S. Weyers, E. Benkler, B. Lipphardt, S. D¨ orscher, K. Stahl, J. Klose, C. Lisdat, M. Filzinger, N. Hunte- mann, E. Peik, and T. E. Mehlst¨ aubler, 115in+−172yb+ Coulomb crystal clock with 2.5×10 −18 systematic un- certainty, Phys. Rev. Lett.134, 0...
work page 2025
- [21]
-
[22]
W. M. Itano, J. J. Bollinger, J. N. Tan, B. Jelenkovi´ c, X.-P. Huang, and D. J. Wineland, Bragg diffraction from crystallized ion plasmas, Science279, 686 (1998)
work page 1998
-
[23]
W. M. Itano, L. R. Brewer, D. J. Larson, and D. J. Wineland, Perpendicular laser cooling of a rotating ion plasma in a Penning trap, Phys. Rev. A38, 5698 (1988)
work page 1988
-
[24]
M. Drewsen, C. Brodersen, L. Hornekær, J. S. Hangst, and J. P. Schifffer, Large ion crystals in a linear Paul trap, Phys. Rev. Lett.81, 2878 (1998)
work page 1998
-
[25]
D. H. E. Dubin and T. M. O’Neil, Trapped nonneutral plasmas, liquids, and crystals (the thermal equilibrium states), Rev. Mod. Phys.71, 87 (1999)
work page 1999
- [26]
-
[27]
A. Shankar, C. Tang, M. Affolter, K. Gilmore, D. H. E. Dubin, S. Parker, M. J. Holland, and J. J. Bollinger, Broadening of the drumhead-mode spectrum due to in- plane thermal fluctuations of two-dimensional trapped ion crystals in a Penning trap, Phys. Rev. A102, 053106 (2020)
work page 2020
-
[28]
W. Johnson, A. Shankar, J. Zaris, J. J. Bollinger, and S. E. Parker, Rapid cooling of the in-plane motion of two- dimensional ion crystals in a Penning trap to millikelvin temperatures, Phys. Rev. A109, L021102 (2024)
work page 2024
-
[29]
W. Johnson, B. Bullock, A. Shankar, J. Zaris, J. J. Bollinger, and S. E. Parker, Adiabatic cooling of planar motion in a Penning-trap ion crystal to sub-millikelvin temperatures, Phys. Rev. A112, 063110 (2025)
work page 2025
-
[30]
A. Mortensen, E. Nielsen, T. Matthey, and M. Drewsen, Observation of three-dimensional long-range order in small ion Coulomb crystals in an rf trap, Phys. Rev. Lett. 96, 103001 (2006)
work page 2006
-
[31]
C. Tang, D. Meiser, J. J. Bollinger, and S. E. Parker, First principles simulation of ultracold ion crystals in a Penning trap with Doppler cooling and a rotating wall potential, Physics of Plasmas26, 073504 (2019)
work page 2019
-
[32]
A. Poindron, J. Pedregosa-Gutierrez, and C. Champ- enois, Thermal bistability in laser-cooled trapped ions, Phys. Rev. A108, 013109 (2023)
work page 2023
- [33]
-
[34]
L. Greengard and V. Rokhlin, A fast algorithm for par- ticle simulations, Journal of Computational Physics73, 325 (1987)
work page 1987
-
[35]
L. Greengard and V. Rokhlin, A new version of the fast multipole method for the Laplace equation in three di- mensions, Acta Numerica6, 229–269 (1997)
work page 1997
- [36]
-
[37]
C.-C. J. Wang, A. C. Keith, and J. K. Freericks, Phonon- mediated quantum spin simulator employing a planar ionic crystal in a Penning trap, Phys. Rev. A87, 013422 (2013)
work page 2013
-
[38]
D. H. E. Dubin, Normal modes, rotational inertia, and thermal fluctuations of trapped ion crystals, Physics of Plasmas27, 102107 (2020)
work page 2020
-
[39]
S. Jain, J. Alonso, M. Grau, and J. P. Home, Scalable arrays of micro-Penning traps for quantum computing and simulation, Phys. Rev. X10, 031027 (2020)
work page 2020
-
[40]
D. H. E. Dubin, Theory of structural phase transitions in a trapped Coulomb crystal, Phys. Rev. Lett.71, 2753 (1993)
work page 1993
-
[41]
D. H. E. Dubin, Theory of electrostatic fluid modes in a cold spheroidal non-neutral plasma, Phys. Rev. Lett.66, 2076 (1991)
work page 2076
-
[42]
A. W. Trivelpiece and R. W. Gould, Space charge waves in cylindrical plasma columns, Journal of Applied Physics 30, 1784 (1959)
work page 1959
-
[43]
S. B. Torrisi, J. W. Britton, J. G. Bohnet, and J. J. Bollinger, Perpendicular laser cooling with a rotating- wall potential in a Penning trap, Phys. Rev. A93, 043421 (2016)
work page 2016
-
[44]
T. B. Mitchell, J. J. Bollinger, W. M. Itano, and D. H. E. Dubin, Stick-slip dynamics of a stressed ion crystal, Phys. Rev. Lett.87, 183001 (2001)
work page 2001
-
[45]
L. Patacchini and I. Hutchinson, Explicit time-reversible orbit integration in particle in cell codes with static homogeneous magnetic field, Journal of Computational Physics228, 2604 (2009)
work page 2009
-
[46]
D. H. E. Dubin and T. M. O’Neil, Computer simulation of ion clouds in a Penning trap, Phys. Rev. Lett.60, 511 (1988)
work page 1988
-
[47]
H. Totsuji, T. Kishimoto, C. Totsuji, and K. Tsuruta, Competition between two forms of ordering in finite Coulomb clusters, Phys. Rev. Lett.88, 125002 (2002)
work page 2002
- [48]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.