Athreya Shankar

We propose and characterize a protocol for rapid near-ground state cooling of the center-of-mass (c.m.) mode of a large trapped ion crystal. When the initial mean thermal occupation of the mode $\bar{n}_i$ is small compared to the number of ions $N$, a red sideband drive implements a beam-splitter type SWAP operation between the mode and the collective spin of the $N$ ions, with the latter effectively serving as a quantum harmonic oscillator. Subsequently, a reset of the spins removes the entropy, leading to near-ground state cooling of the c.m. mode. We term this protocol as quantum beam-splitter cooling (QBSC). We analyze the impact of several practical imperfections on the final temperature achievable under QBSC, including finite ion number, off-resonant carrier and blue-sideband contributions, and the impact of the sideband drives arising from spectator modes. In addition, we outline practical strategies to eliminate the carrier drive. Furthermore, we show that measuring the population statistics of the ions at the end of the SWAP operation can enable near-optimal quantum beam-splitter thermometry (QBST), with the classical Fisher information approaching the quantum Fisher information of a thermal state. We discuss the connection of QBSC with continuous sideband cooling and compare QBST with a recently proposed rapid adiabatic passage-based thermometry scheme. Our work constitutes an example of harnessing many-body effects to open new routes to laser cooling and thermometry in large trapped ion crystals.

Large ion crystals in Penning traps provide a platform for quantum simulation and sensing with hundreds of spins, but their continuous rigid-body rotation has so far limited flexible local qubit control. Here we demonstrate programmable site-selective spin control across large rotating ${}^{9}\mathrm{Be}^{+}$ crystals in a Penning trap. A tightly focused off-resonant laser beam drives local $R_z$ phase rotations via differential AC Stark shifts. Beam steering synchronised with crystal rotation enables addressing of arbitrary ions throughout the crystal. Ramsey-based characterisation shows $R_z(π)$ gate fidelity of 94.6% and nearest-neighbour crosstalk of 1.2%. We use this capability to prepare spatially structured spin patterns, generating a biskyrmion spin texture in a single-layer crystal, then extending the method to bilayer crystals we perform layer-selective addressing operations. We further demonstrate dual-quadrature Ramsey sensing by imprinting a relative $π/2$ phase shift between spatial sub-ensembles, enabling simultaneous measurement of orthogonal spin components within a single experimental realisation. These results establish programmable local control in large rotating ion crystals, opening new routes for engineering spatially structured quantum states in multidimensional trapped-ion systems.

Quantum simulation with controllable many-body platforms offers a powerful route to exploring complex phases and dynamics that are difficult to access in natural materials. Among these, topological spin textures such as skyrmions are central to modern condensed-matter physics and play a key role in chiral quantum many-body systems. Their controlled realisation in large, programmable quantum platforms, however, remains an outstanding challenge. Here, we report deterministic generation and site-resolved reconstruction of topological spin textures in a two-dimensional crystal of more than 150 trapped ions. Using globally applied spin-dependent forces, we generate skyrmion configurations and reconstruct the full vector spin field with single-ion resolution, obtaining a winding number of 0.99$\pm$0.02 and a mean local fidelity of 0.87$\pm$0.04. In addition, we implement single-ion-resolved control to deterministically prepare domain-wall states, extending our approach to a broader class of non-uniform spin textures. These results establish trapped-ion crystals as a platform for engineering complex spin textures and open the door to exploring topology-dependent nonequilibrium dynamics in long-range interacting quantum systems.

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.

We develop a unified many-body theory of two-photon dark-state laser cooling, the workhorse for preparing trapped ions close to their motional quantum ground state. For ions with a $Λ$ level structure, driven by Raman lasers, we identify an ion-number-dependent crossover between weak and strong coupling where both the cooling rate and final temperature are simultaneously optimized. We obtain simple analytic results in both extremes: In the weak coupling limit, we show a Lorentzian spin-absorption spectrum determines the cooling rate and final occupation of the motional state, which are both independent of the number of ions. We also highlight the benefit of including an additional spin dependent force in this case. In the strong coupling regime, our theory reveals the role of collective dynamics arising from phonon exchange between dark and bright states, allowing us to explain the enhancement of the cooling rate with increasing ion number. Our analytic results agree closely with exact numerical simulations and provide experimentally accessible guidelines for optimizing cooling in large ion crystals, a key step toward scalable, high-fidelity trapped-ion quantum technologies.

Displacement sensing is a fundamental task in metrology. However, the development of quantum-enhanced sensors that fully utilize the available degrees of freedom in many-body quantum systems remains an outstanding challenge. We propose novel many-body displacement sensing schemes that use spin-dependent squeezed (SDS) states -- hybrid spin-boson states whose bosonic squeezed quadrature is conditioned on an auxiliary spin. We prove that SDS states are \emph{optimal}, i.e. their quantum Cramér-Rao bound saturates the Heisenberg limit. We propose explicit measurement sequences that can be readily implemented in systems such as trapped ions. We also introduce a scalable state-preparation protocol and numerically demonstrate the preparation of $8.7$~dB of spin-dependent squeezing $15$ times faster than the standard approach using second-order sidebands in trapped ions. The potential applications of our sensing protocols range from measuring single-photon scattering to searches for dark matter.

Three-dimensional (3D) crystals offer a route to scale up trapped ion systems for quantum sensing and quantum simulation applications. However, engineering coherent spin-motion couplings and effective spin-spin interactions in large crystals poses technical challenges associated with decoherence and prolonged timescales to generate appreciable entanglement. Here, we explore the possibility to speed up these interactions in 3D crystals via parametric amplification. We derive a general Hamiltonian for the parametric amplification of spin-motion coupling that is applicable to crystals of any dimension in both rf Paul traps and Penning traps. Unlike in lower dimensional crystals, we find that the ability to faithfully (uniformly) amplify the spin-spin interactions in 3D crystals depends on the physical implementation of the spin-motion coupling. We consider the light-shift (LS) gate, and the so-called phase-insensitive and phase-sensitive Mølmer-Sørensen (MS) gates, and find that only the latter gate can be faithfully amplified in general 3D crystals. We discuss a situation where non-uniform amplification can be advantageous. We also reconsider the impact of counter-rotating terms on parametric amplification and find that they are not as detrimental as previous studies suggest.

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.

Electromagnetically Induced Transparency (EIT) cooling is a well-established method for preparing trapped ion systems in their motional ground state. However, isolating a three-level system, as required for EIT cooling, is often challenging or impractical. Nonetheless, multilevel systems can inherently host dark states. In this work, we extend the EIT cooling framework to such multilevel systems. We develop a formalism to accurately determine the cooling rate in the weak sideband coupling regime and provide an approximate estimate for cooling rates beyond this regime, without the need for explicit simulation of the motional degree of freedom. We clarify the connection between the cooling rate and the absorption spectrum, offering a pathway for efficient near-ground-state cooling of ions with complex electronic structures.

We propose a protocol for the preparation of generalized Greenberger-Horne-Zeilinger (GHZ) states of $N$ atoms each with $d=3$ or $4$ internal levels. We generalize the celebrated one-axis twisting (OAT) Hamiltonian for $N$ qubits to qudits by including OAT interactions of equal strengths between every pair of qudit levels, a protocol we call as balanced OAT (BOAT). Analogous to OAT for qubits, we find that starting from a product state of an arbitrary number of atoms $N$, dynamics under BOAT leads to the formation of GHZ states for qutrits ($d=3$) and ququarts ($d=4$). While BOAT could potentially be realized on several platforms where all-to-all coupling is possible, here we propose specific implementations using trapped ion systems. We show that preparing these states with a fidelity above a threshold value rules out lower dimensional entanglement than that of the generalized GHZ states. For qutrits, we also propose a protocol to bound the fidelity that requires only global addressing of the ion crystal and single-shot readout of one of the levels. Our results open a path for the scalable generation and certification of high-dimensional multipartite entanglement on current atom-based quantum hardware.

We simulate the dynamics, including laser cooling, of 3D ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions' equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions' absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1,000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3D ion crystals could be well-suited as platforms for future quantum science experiments.

We study optimal quantum sensing of multiple physical parameters using repeated measurements. In this scenario, the Fisher information framework sets the fundamental limits on sensing performance, yet the optimal states and corresponding measurements that attain these limits remain to be discovered. To address this, we extend the Fisher information approach with a second optimality requirement for a sensor to provide unambiguous estimation of unknown parameters. We propose a systematic method integrating Fisher information and Bayesian approaches to quantum metrology to identify the combination of input states and measurements that satisfies both optimality criteria. Specifically, we frame the optimal sensing problem as an optimization of an asymptotic Bayesian cost function that can be efficiently solved numerically and, in many cases, analytically. We refer to the resulting optimal sensor as a `quantum compass' solution, which serves as a direct multiparameter counterpart to the Greenberger-Horne-Zeilinger state-based interferometer, renowned for achieving the Heisenberg limit in single-parameter metrology. We provide exact quantum compass solutions for paradigmatic multiparameter problem of sensing two and three parameters using an SU(2) sensor. Our metrological cost function opens avenues for quantum variational techniques to design low-depth quantum circuits approaching the optimal sensing performance in the many-repetition scenario. We demonstrate this by constructing simple quantum circuits that achieve the Heisenberg limit for vector field and 3D rotations estimation using a limited set of gates available on a trapped-ion platform. Our work introduces and optimizes sensors for a practical notion of optimality, keeping in mind the ultimate goal of quantum sensors to precisely estimate unknown parameters.

Trapped ion systems are a leading platform for quantum information processing, but they are currently limited to 1D and 2D arrays, which imposes restrictions on both their scalability and their range of applications. Here, we propose a path to overcome this limitation by demonstrating that Penning traps can be used to realize remarkably clean bilayer crystals, wherein hundreds of ions self-organize into two well-defined layers. These bilayer crystals are made possible by the inclusion of an anharmonic trapping potential, which is readily implementable with current technology. We study the normal modes of this system and discover salient differences compared to the modes of single-plane crystals. The bilayer geometry and the unique properties of the normal modes open new opportunities, in particular in quantum sensing and quantum simulation, that are not straightforward in single-plane crystals. Furthermore, we illustrate that it may be possible to extend the ideas presented here to realize multilayer crystals with more than two layers. Our work increases the dimensionality of trapped ion systems by efficiently utilizing all three spatial dimensions and lays the foundation for a new generation of quantum information processing experiments with multilayer 3D crystals of trapped ions.

We propose a highly feasible technique with no experimental overhead to rapidly cool the in-plane degrees of freedom of large two-dimensional ion crystals in Penning traps. Through simulations, we demonstrate that our approach enables the in-plane modes to cool down to a temperature of around 1 mK in less than 10 ms. Our technique relies on near-resonant coupling of the poorly cooled in-plane motions and the efficiently cooled out-of-plane motions, and is achieved without introducing additional potentials. The rapid cooling enabled by our approach is in contrast to typical operating conditions, where our simulations of the laser cooling dynamics suggest that the ion crystal's in-plane motion cools very slowly on a timescale of several hundreds of milliseconds, a rate likely slower than experimental heating rates. Our work sets the stage for sub-Doppler laser cooling of the planar motion, and more robust and versatile quantum simulation and quantum sensing experiments with two-dimensional crystals in Penning traps.

Penning traps have been used for performing quantum simulations and sensing with hundreds of ions and provide a promising route toward scaling up trapped ion quantum platforms because of the ability to trap and control up to thousands of ions in 2D and 3D crystals. A leading source of decoherence in laser-based multiqubit operations on trapped ions is off-resonant spontaneous emission. While many trapped ion quantum computers or simulators utilize clock qubits, other systems rely on Zeeman qubits, which require a more complex calculation of this decoherence. We examine theoretically the impacts of spontaneous emission on quantum gates performed with trapped ions in a high magnetic field. We consider two types of gates -- light-shift and Molmer-Sorensen gates -- and compare the decoherence errors in each. We also compare different detunings, polarizations, and required intensities of the laser beams used to drive the gates. We show that both gates can have similar performance at their optimal operating conditions and examine the experimental feasibility of various operating points. By examining the magnetic field dependence of each gate, we demonstrate that when the $P$ state fine structure splitting is large compared to the Zeeman splittings, the theoretical performance of the Molmer-Sorensen gate is significantly better than that of the light-shift gate. Additionally, for the light-shift gate, we make an approximate comparison between the fidelities that can be achieved at high fields with the fidelities of state-of-the-art two-qubit trapped ion quantum gates. We show that, with regard to spontaneous emission, the achievable infidelity of our current configuration is about an order of magnitude larger than that of the best low-field gates, but we also discuss alternative configurations with potential error rates that are comparable with state-of-the-art trapped ion gates.

Two-dimensional $p+ i p$ superconductors and superfluids are systems that feature chiral behavior emerging from the Cooper pairing of electrons or neutral fermionic atoms with non-zero angular momentum. Their realization has been a longstanding goal because they offer great potential utility for quantum computation and memory. However, they have so far eluded experimental observation both in solid state systems as well as in ultracold quantum gases. Here, we propose to leverage the tremendous control offered by rotating two-dimensional trapped-ion crystals in a Penning trap to simulate the dynamical phases of two-dimensional $p+ip$ superfluids. This is accomplished by mapping the presence or absence of a Cooper pair into an effective spin-1/2 system encoded in the ions' electronic levels. We show how to infer the topological properties of the dynamical phases, and discuss the role of beyond mean-field corrections. More broadly, our work opens the door to use trapped ion systems to explore exotic models of topological superconductivity and also paves the way to generate and manipulate skyrmionic spin textures in these platforms.

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.

We study, both experimentally and theoretically, electromagnetically induced transparency cooling of the drumhead modes of planar 2-dimensional arrays with up to $N\approx 190$ Be${}^+$ ions stored in a Penning trap. Substantial sub-Doppler cooling is observed for all $N$ drumhead modes. Quantitative measurements for the center-of-mass mode show near ground state cooling with motional quantum numbers of $\bar{n} = 0.3\pm0.2$ obtained within $200~μs$. The measured cooling rate is faster than that predicted by single particle theory, consistent with a quantum many-body calculation. For the lower frequency drumhead modes, quantitative temperature measurements are limited by apparent damping and frequency instabilities, but near ground state cooling of the full bandwidth is strongly suggested. This advancement will greatly improve the performance of large trapped ion crystals in quantum information and quantum metrology applications.

Penning traps, with their ability to control planar crystals of tens to hundreds of ions, are versatile quantum simulators. Thermal occupations of the motional drumhead modes, transverse to the plane of the ion crystal, degrade the quality of quantum simulations. Laser cooling using electromagnetically induced transparency (EIT cooling) is attractive as an efficient way to quickly initialize the drumhead modes to near ground-state occupations. We numerically investigate the efficiency of EIT cooling of planar ion crystals in a Penning trap, accounting for complications arising from the nature of the trap and from the simultaneous cooling of multiple ions. We show that, in spite of challenges, the large bandwidth of drumhead modes (hundreds of kilohertz) can be rapidly cooled to near ground-state occupations within a few hundred microseconds. Our predictions for the center-of-mass mode include a cooling time constant of tens of microseconds and an enhancement of the cooling rate with increasing number of ions. Successful experimental demonstrations of EIT cooling in the NIST Penning trap [E. Jordan, K. A. Gilmore, A. Shankar, A. Safavi-Naini, M. J. Holland, and J. J. Bollinger, "Near ground-state cooling of two-dimensional trapped-ion crystals with more than 100 ions", (2018), submitted.] validate our predictions.

Ultranarrow-linewidth atoms coupled to a lossy optical cavity mode synchronize, i.e. develop correlations, and exhibit steady-state superradiance when continuously repumped. This type of system displays rich collective physics and promises metrological applications. These features inspire us to investigate if analogous spin synchronization is possible in a different platform that is one of the most robust and controllable experimental testbeds currently available: ion-trap systems. We design a system with a primary and secondary species of ions that share a common set of normal modes of vibration. In analogy to the lossy optical mode, we propose to use a lossy normal mode, obtained by sympathetic cooling with the secondary species of ions, to mediate spin synchronization in the primary species of ions. Our numerical study shows that spin-spin correlations develop, leading to a macroscopic collective spin in steady-state. We propose an experimental method based on Ramsey interferometry to detect signatures of this collective spin; we predict that correlations prolong the visibility of Ramsey fringes, and that population statistics at the end of the Ramsey sequence can be used to directly infer spin-spin correlations.