Chen Tang

Planar thermal equilibration is studied using direct numerical simulations of ultracold two-dimensional (2D) ion crystals in a Penning trap with a rotating wall. The large magnetic field of the trap splits the modes that describe in-plane motion of the ions into two branches: High frequency cyclotron modes dominated by kinetic energy and low frequency $\mathbf{E \times B}$ modes dominated by potential energy associated with thermal position displacements. Using an eigenmode analysis we extract the equilibration rate between these two branches as a function of the ratio of the frequencies that characterize the two branches and observe this equilibration rate to be exponentially suppressed as the ratio increases. Under experimental conditions relevant for current work at NIST, the predicted equilibration time is orders of magnitude longer than any relevant experimental timescales. We also study the coupling rate dependence on the thermal temperature and the number of ions. Besides, we show how increasing the rotating wall strength improves crystal stability. These details of in-plane mode dynamics help set the stage for developing strategies to efficiently cool the in-plane modes and improve the performance of single-plane ion crystals for quantum information processing.

Two-dimensional crystals of ions stored in Penning traps are a leading platform for quantum simulation and sensing experiments. For small amplitudes, the out-of-plane motion of such crystals can be described by a discrete set of normal modes called the drumhead modes, which can be used to implement a range of quantum information protocols. However, experimental observations of crystals with Doppler-cooled and even near-ground-state-cooled drumhead modes reveal an unresolved drumhead mode spectrum. In this work, we establish in-plane thermal fluctuations in ion positions as a major contributor to the broadening of the drumhead mode spectrum. In the process, we demonstrate how the confining magnetic field leads to unconventional in-plane normal modes, whose average potential and kinetic energies are not equal. This property, in turn, has implications for the sampling procedure required to choose the in-plane initial conditions for molecular dynamics simulations. For current operating conditions of the NIST Penning trap, our study suggests that the two dimensional crystals produced in this trap undergo in-plane potential energy fluctuations of the order of $10$ mK. Our study therefore motivates the need for designing improved techniques to cool the in-plane degrees of freedom.

A direct numerical simulation of many interacting ions in a Penning trap with a rotating wall is presented. The ion dynamics is modelled classically. Both axial and planar Doppler laser cooling are modeled using stochastic momentum impulses based on two-level atomic scattering rates. The plasmas being modeled are ultra-cold two-dimensional crystals made up of 100's of ions. We compare Doppler cooled results directly to a previous linear eigenmodes analysis. Agreement in both frequency and mode structure are obtained. Additionally, when Doppler laser cooling is applied, the laser cooled steady state plasma axial temperature agrees with the Doppler cooling limit. Numerical simulations using the approach described and benchmarked here will provide insights into the dynamics of large trapped-ion crystals, improving their performance as a platform for quantum simulation and sensing.