pith. sign in

arxiv: 2602.22541 · v2 · pith:RB5US35Wnew · submitted 2026-02-26 · 🪐 quant-ph · physics.plasm-ph

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

classification 🪐 quant-ph physics.plasm-ph
keywords laser coolingPenning trapion crystalsquantum sensingDoppler coolingnumerical simulationE x B modes3D crystals
0
0 comments X

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.

The paper develops a numerical framework to simulate laser cooling of up to 10^5 ions stored in a Penning trap. This framework is applied to characterize and optimize the cooling of ellipsoidal 3D crystals. New pathways to enhanced cooling are found by adding an axial component to the potential energy-dominated E×B modes. Greatly enhanced cooling of the perpendicular kinetic energy to below 1 mK is observed in prolate ion crystals, which enables a simplified cooling beam setup. Specific values of trap and laser beam parameters are proposed for optimal cooling in a variety of examples, illustrating the feasibility of preparing large 3D crystals for high-sensitivity quantum sensing protocols.

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

Figures reproduced from arXiv: 2602.22541 by Allison L. Carter, Athreya Shankar, Daniel H.E. Dubin, John J. Bollinger, John Zaris, Scott E. Parker, Wes Johnson.

Figure 1
Figure 1. Figure 1: FIG. 1: The setup to cool the crystal relies on two [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: (a) Two examples of 3D crystal normal modes [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: (a) The shape of [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Theory estimates for the maximum [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The rms ion displacements during the last 1 ms [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: The kinetic and potential energy cooling of the [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: The perpendicular kinetic energy and total po [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: The simulated perpendicular kinetic energy af [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: (a) Mode spectra and mode [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: (a) Local equilibrium ion configuration of an [PITH_FULL_IMAGE:figures/full_fig_p014_11.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [Introduction / Model] Notation for the axial perturbation strength and laser detuning should be defined explicitly in the first appearance rather than assumed from context.
  2. [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

2 responses · 0 unresolved

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

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

0 steps flagged

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

1 free parameters · 1 axioms · 0 invented entities

Ledger inferred from abstract only. The simulation rests on standard ion-trap physics; free parameters are the optimized trap and laser settings reported as outcomes of the numerical search.

free parameters (1)
  • axial component strength and laser detuning values
    Specific numerical values proposed for optimal cooling; these are outputs of the simulation rather than external inputs.
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.
    Invoked throughout the characterization of ellipsoidal crystals and the identification of enhanced cooling pathways.

pith-pipeline@v0.9.0 · 5749 in / 1412 out tokens · 80427 ms · 2026-05-21T13:06:49.113379+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

48 extracted references · 48 canonical work pages · 1 internal anchor

  1. [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)

  2. [2]

    Safavi-Naini, R

    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)

  3. [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)

  4. [4]

    Shankar, E

    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)

  5. [5]

    M. Qiao, Z. Cai, Y. Wang, B. Du, N. Jin, W. Chen, P. Wang, C. Luan, E. Gao, X. Sun, H. Tian, J. Zhang, and K. Kim, arXiv (2022), arXiv:2204.07283

  6. [6]

    J. H. Pham, J. Y. Z. Jee, A. Rischka, M. J. Biercuk, and R. N. Wolf, In situ-tunable spin–spin interactions in a Penning trap with in-bore optomechanics, Adv. Quantum Technol.7, 2400086 (2024), arXiv:2401.17742 [quant-ph]

  7. [7]

    Swingle, G

    B. Swingle, G. Bentsen, M. Schleier-Smith, and P. Hay- den, Measuring the scrambling of quantum information, Phys. Rev. A94, 040302 (2016)

  8. [8]

    Measuring out-of-time-order correlations and multiple quantum spectra in a trapped ion quantum magnet

    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]

  9. [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)

  10. [10]

    Pezz` e, A

    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)

  11. [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)

  12. [12]

    Affolter, K

    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)

  13. [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)

  14. [14]

    Toscano, D

    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)

  15. [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)

  16. [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)

  17. [17]

    Kiesenhofer, H

    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

  18. [18]

    S. A. Guo et al., A site-resolved two-dimensional quan- tum simulator with hundreds of trapped ions, Nature 630, 613 (2024), arXiv:2311.17163 [quant-ph]

  19. [19]

    Hawaldar, P

    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)

  20. [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...

  21. [21]

    D. R. Leibrandt, S. G. Porsev, C. Cheung, and M. S. Safronova, Prospects of a thousand-ion sn 2+ Coulomb- crystal clock with sub-10 −19 inaccuracy, Nature Com- mun.15, 5663 (2024), arXiv:2205.15484 [physics.atom- ph]

  22. [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)

  23. [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)

  24. [24]

    Drewsen, C

    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)

  25. [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)

  26. [26]

    Jordan, K

    E. Jordan, K. A. Gilmore, A. Shankar, A. Safavi- Naini, J. G. Bohnet, M. J. Holland, and J. J. Bollinger, Near ground-state cooling of two-dimensional trapped- ion crystals with more than 100 ions, Phys. Rev. Lett. 122, 053603 (2019)

  27. [27]

    Shankar, C

    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)

  28. [28]

    Johnson, A

    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)

  29. [29]

    Johnson, B

    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)

  30. [30]

    Mortensen, E

    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)

  31. [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)

  32. [32]

    Poindron, J

    A. Poindron, J. Pedregosa-Gutierrez, and C. Champ- enois, Thermal bistability in laser-cooled trapped ions, Phys. Rev. A108, 013109 (2023)

  33. [33]

    Zaris, W

    J. Zaris, W. Johnson, A. Shankar, J. J. Bollinger, and S. E. Parker, Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method, Journal of Plasma Physics91, E53 (2025)

  34. [34]

    Greengard and V

    L. Greengard and V. Rokhlin, A fast algorithm for par- ticle simulations, Journal of Computational Physics73, 325 (1987)

  35. [35]

    Greengard and V

    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)

  36. [36]

    Cheng, L

    H. Cheng, L. Greengard, and V. Rokhlin, A fast adap- tive multipole algorithm in three dimensions, Journal of Computational Physics155, 468 (1999)

  37. [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)

  38. [38]

    D. H. E. Dubin, Normal modes, rotational inertia, and thermal fluctuations of trapped ion crystals, Physics of Plasmas27, 102107 (2020)

  39. [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)

  40. [40]

    D. H. E. Dubin, Theory of structural phase transitions in a trapped Coulomb crystal, Phys. Rev. Lett.71, 2753 (1993)

  41. [41]

    D. H. E. Dubin, Theory of electrostatic fluid modes in a cold spheroidal non-neutral plasma, Phys. Rev. Lett.66, 2076 (1991)

  42. [42]

    A. W. Trivelpiece and R. W. Gould, Space charge waves in cylindrical plasma columns, Journal of Applied Physics 30, 1784 (1959)

  43. [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)

  44. [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)

  45. [45]

    Patacchini and I

    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)

  46. [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)

  47. [47]

    Totsuji, T

    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)

  48. [48]

    Huang, F

    X.-P. Huang, F. Anderegg, E. M. Hollmann, C. F. Driscoll, and T. M. O’Neil, Steady-state confinement of non-neutral plasmas by rotating electric fields, Phys. Rev. Lett.78, 875 (1997)