Adiabatic Cooling of Planar Motion in a Penning Trap Ion Crystal to Sub-Millikelvin Temperatures

Two-dimensional planar ion crystals in a Penning trap are a platform for quantum information science experiments. However, the low-frequency planar modes of these crystals are not efficiently cooled by laser cooling, which can limit the utility of the drumhead modes for quantum information processing. Recently, it has been shown that nonlinear mode coupling can enhance the cooling of the low-frequency planar modes. Here, we demonstrate in numerical simulations that this coupling can be dynamically tuned by adiabatically changing the rotation frequency of the ion crystal during experiments. Furthermore, we show that this technique can, in addition, produce lower temperatures for the low-frequency planar modes via an adiabatic cooling process. This result allows cooling of the planar modes to sub-millikelvin temperatures, resulting in improved spectral resolution of the drumhead modes at experimentally relevant rotation frequencies, which is crucial for quantum information processing applications.

💡 Research Summary

The paper investigates a method to cool the low‑frequency planar (E × B) modes of a two‑dimensional ion crystal confined in a Penning trap to sub‑millikelvin temperatures, a regime that is difficult to reach with standard Doppler laser cooling alone. The authors focus on 9Be⁺ crystals, which exhibit three families of normal modes: low‑frequency E × B modes (∼1–100 kHz), intermediate‑frequency drumhead modes (∼1 MHz), and high‑frequency cyclotron modes (∼10 MHz). While laser cooling efficiently damps the kinetic‑energy‑dominated cyclotron modes, the E × B modes are dominated by potential energy (characterized by a large ratio Rₙ = PE/KE ≫ 1) and thus remain relatively hot, broadening the drumhead spectrum and limiting quantum‑information‑processing (QIP) fidelity.

Previous work demonstrated that a near‑resonant nonlinear coupling between the E × B modes and drumhead modes can transfer energy from the former to the latter, enabling rapid cooling of the planar modes. However, this mechanism required operating the crystal at relatively high rotation frequencies ω_r, which reduces the Lamb‑Dicke parameter η for the drumhead modes and is undesirable for many QIP protocols that need small η. The present study therefore seeks a way to cool the planar modes while keeping ω_r low.

The central idea is to exploit adiabatic invariance of the action variable I for each harmonic normal mode. In the co‑rotating frame the Hamiltonian of a mode can be written as H = ω I, where ω is the mode frequency. If ω is varied slowly compared to the mode period, I remains constant and the mode energy scales as E_f = E_i · (ω_f/ω_i). The authors show analytically that the frequencies of the E × B modes are approximately proportional to the planar confinement parameter β, which itself depends on ω_r via β = ω_r(ω_c − ω_r)/ω_z² − ½. Consequently, a controlled reduction of ω_r leads to a proportional reduction of β and thus a reduction of the planar mode energies: E_f = E_i · (β_f/β_i).

To implement this adiabatic frequency sweep without destabilizing the crystal, three practical conditions are identified:

1. Fixed δ/β ratio – The rotating‑wall strength δ must be scaled together with β so that the aspect ratio of the crystal remains constant. Otherwise, a decreasing β with fixed δ would elongate the crystal, causing structural re‑configuration and heating.

2. Smooth ramp profile – Abrupt changes in the time derivative of ω_r generate impulsive Euler forces that excite a collective rocking mode. The authors adopt a half‑cosine ramp, ω_r(t) = ω_i + (ω_f − ω_i)/2 ·