Ana Maria Rey

Scaling trapped-ion quantum sensors from single ions to large ensembles is a key challenge for next-generation precision measurements. At the same time, many ion species of interest for optical clocks and tests of fundamental physics lack closed cycling transitions required for direct laser cooling and state detection. Collective quantum logic spectroscopy addresses both limitations by coupling an ensemble of sensor, or spectroscopy, ions to one or more logic ions that provide sympathetic cooling and state readout. Here, we establish the fundamental performance limits and operating regimes of this protocol, identifying how the interaction strength, interrogation time, and logic-ensemble size govern sensitivity, dynamic range, and robustness to experimental imperfections. We show that quantum-limited sensitivity can be retained even with a single logic ion, while increasing the number of logic ions substantially improves readout efficiency and robustness. Beyond precision metrology, the same collective interface enables many-body measurements relevant to quantum information processing, including parity measurements and stabilizer-like syndrome extraction. Our results establish collective quantum logic spectroscopy as a scalable framework for optical clocks, quantum-enhanced sensing, and trapped-ion quantum information processing.

Quantum many-body systems driven far from equilibrium can exhibit chaos, entanglement, and non-classical correlations, yet directly observing these phenomena in large, closed quantum systems remains challenging. Here we realize the Dicke model -- a fundamental description of light-matter interactions -- in a two-dimensional crystal of approximately 100 trapped ions. The ions' internal state is optically coupled to the center of mass vibrational mode via an optical spin-dependent force, enabling unitary many-body dynamics beyond the mean-field and few-body limits. In the integrable regime, where the phonons can be adiabatically eliminated, we observe a dynamical phase transition between ferromagnetic to paramagnetic spin phases. In contrast, when the spins and phonons are strongly coupled, we observe clear signatures of non-integrable chaotic dynamics, including erratic phase-space trajectories and the exponential growth of excitations and entanglement quantified by the one-body Rényi entropy. By quenching from an unstable fixed point in the near-integrable regime, quantum noise can generate correlated spin-phonon excitations. Our numerical calculations, in clear agreement with experiment, reveal the generation of two-mode spin-phonon squeezing, 2.6 dB below the standard quantum limit (4.6 dB relative to the initial thermal state), followed by generalized vacuum Rabi collapses and revivals. Our results establish large ion crystals as scalable analog quantum simulators of non-equilibrium light-matter dynamics and provide a controlled platform for experimental studies of information scrambling and entanglement in closed many-body systems.

We show that quantum information scrambling can enable a generic SWAP gate between collective degrees of freedom in systems without universal local control. Our protocol combines the Hayden-Preskill recovery scheme, associated with the black hole information paradox, with quantum teleportation and runs them in parallel and in opposite directions, enabling bidirectional exchange of quantum states through global interactions alone. This approach cleanly distinguishes the roles of information spreading, entanglement, and chaos for enabling both coherent state transfer and recovery. We propose an experimental realization using the Dicke model, which can be realized in cavity-QED and trapped-ion platforms, highlighting the utility of holography in designing practical quantum gates.

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.

We propose the use of motional states of two interacting atoms trapped in a potential stroboscopically engineered by an optical tweezer as a means to implement a qubit-oscillator system, in analogy to those implemented in circuit quantum electrodynamics and trapped ions. In our setting, the center of mass degree of freedom of the atoms plays the role of a photon or phonon mode, while the interacting, relative mode acts as a qubit. No internal state is involved in our system, which makes this motional qubit robust to spin-dependent noise. We show that a universal set of bosonic operations, including displacement, rotation, squeezing, and the corresponding set of gates controlled by the qubit, can be implemented through precise temporal modulation of the optical tweezers. We numerically check that these gates can be generated with high fidelity, and discuss possible schemes for initial state preparation and final state readout. While we restrict the discussion to a single qubit-oscillator module, scalability can be achieved by coupling arrays of atoms via dipolar or Rydberg-dressed interactions.

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.

We propose the use of phonon$\unicode{x2013}$mediated interactions as an entanglement resource to engineer Einstein$\unicode{x2013}$Podolsky$\unicode{x2013}$Rosen (EPR) correlations and to perform teleportation of collective spin states in two$\unicode{x2013}$dimensional ion crystals. We emulate continuous variable quantum teleportation protocols between subsystems corresponding to different nuclear spin degrees of freedom. In each of them, a quantum state is encoded in an electronic spin degree of freedom that couples to the vibrational modes of the crystal. We show that high fidelity teleportation of spin-coherent states and their phase-displaced variant, entangled spin-squeezed states, and Dicke states, is possible for realistic experimental conditions in arrays from a few tens to a few hundred ions.

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.

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.

We study the non-equilibrium dynamics of dipoles confined in multiple stacked two-dimensional layers realising a long-range interacting quantum spin 1/2 XXZ model. We demonstrate that strong in-plane XXX interactions can protect a manifold of collective layer dynamics. This then allows us to map the many-body spin dynamics to bosonic models. In a bilayer configuration we show how to engineer the paradigmatic two-mode squeezing Hamiltonian known from quantum optics, resulting in exponential production of entangled pairs and generation of metrologically useful entanglement from initially prepared product states. In multi-layer configurations we engineer a bosonic variant of the Kitaev model displaying chiral propagation along the layer direction. Our study illustrates how the control over interactions, lattice geometry and state preparation in interacting dipolar systems uniquely afforded by AMO platforms such as Rydberg and magnetic atoms, polar molecules or trapped ions allow for the control over the temporal and spatial propagation of correlations for applications in quantum sensing and quantum simulation.

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.

Trapped ions boast long coherence times and excellent gate fidelities, making them a useful platform for quantum information processing. Scaling to larger numbers of ion qubits in RF Paul traps demands great effort. Another technique for trapping ions is via a Penning trap where a 2D crystal of hundreds of ions is formed by controlling the rotation of the ions in the presence of a strong magnetic field. However, the rotation of the ion crystal makes single ion addressability a significant challenge. We propose a protocol that takes advantage of a deformable mirror to introduce AC Stark shift patterns that are static in the rotating frame of the crystal. Through numerical simulations we validate the potential of this protocol to perform high-fidelity single-ion gates in crystalline arrays of hundreds of ions.

We propose protocols for the creation of useful entangled states in a system of spins collectively coupled to a bosonic mode, directly applicable to trapped-ion and cavity QED setups. The protocols use coherent manipulations of the spin-boson interactions naturally arising in these systems to prepare spin squeezed states exponentially fast in time. We demonstrate the robustness of the protocols by analyzing the effects of natural sources of decoherence in these systems and show their advantage compared to more standard slower approaches where entanglement is generated algebraically with time.

Developing the isolation and control of ultracold atomic systems to the level of single quanta has led to significant advances in quantum sensing, yet demonstrating a quantum advantage in real world applications by harnessing entanglement remains a core task. Here, we realize a many-body quantum-enhanced sensor to detect weak displacements and electric fields using a large crystal of $\sim 150$ trapped ions. The center of mass vibrational mode of the crystal serves as high-Q mechanical oscillator and the collective electronic spin as the measurement device. By entangling the oscillator and the collective spin before the displacement is applied and by controlling the coherent dynamics via a many-body echo we are able to utilize the delicate spin-motion entanglement to map the displacement into a spin rotation such that we avoid quantum back-action and cancel detrimental thermal noise. We report quantum enhanced sensitivity to displacements of $8.8 \pm 0.4~$dB below the standard quantum limit and a sensitivity for measuring electric fields of $240\pm10~\mathrm{nV}\mathrm{m}^{-1}$ in $1$ second ($240~\mathrm{nV}\mathrm{m}^{-1}/\sqrt{\mathrm{Hz}}$).

We study the far-from-equilibrium dynamical regimes of a many-body spin boson model with disordered couplings relevant for cavity QED and trapped ions experiments, using the discrete truncated Wigner approximation (DTWA). We focus on the dynamics of spin observables upon varying the disorder strength and the frequency of the photons, finding that the latter can considerably alter the structure of the system's dynamical responses. When the photons evolve at a similar rate as the spins, they can induce qualitatively distinct frustrated dynamics characterized by either logarithmic or algebraically slow relaxation. The latter illustrates resilience of glassy-like dynamics in the presence of active photonic degrees of freedom, suggesting that disordered quantum many body systems with resonant photons or phonons can display a rich diagram of non-equilibrium responses, with near future applications for quantum information science.

Macroscopic arrays of cold atoms trapped in optical cavities can reach the strong atom-light collective coupling regime thanks to the simultaneous interactions of the cavity mode with the atomic ensemble. In a recent work we reported a protocol that takes advantage of the strong and collective atom-light interactions in cavity QED systems for precise electric field sensing in the optical domain. We showed that it can provide between $10$-$20$~dB of metrological gain over the standard quantum limit in current cavity QED experiments operating with long-lived alkaline-earth atoms. Here, we give a more in depth discussion of the protocol using both exact analytical calculations and numerical simulations, and describe the precise conditions under which the predicted enhancement holds after thoroughly accounting for both photon loss and spontaneous emission, natural decoherence mechanisms in current experiments. The analysis presented here not only serves to benchmark the protocol and its utility in cavity QED arrays but also sets the conditions required for its applicability in other experimental platforms such as arrays of trapped ions.

We present a new method to compute short-time expectation values in large collective spin systems with generic Markovian decoherence. Our method is based on a Taylor expansion of a formal solution to the equations of motion for Heisenberg operators. This expansion can be truncated at finite order to obtain virtually exact results at short times that are relevant for metrological applications such as spin squeezing. In order to evaluate the expansion for Heisenberg operators, we compute the relevant structure constants of a collective spin operator algebra. We demonstrate the utility of our method by computing spin squeezing, two-time correlation functions, and out-of-time-ordered correlators for $10^4$ spins in strong-decoherence regimes that are otherwise inaccessible via existing numerical methods. Our method can be straightforwardly generalized to the case of a collective spin coupled to bosonic modes, relevant for trapped ion and cavity QED experiments, and may be used to investigate short-time signatures of quantum chaos and information scrambling.

The coupling of spins to long-wavelength bosonic modes is a prominent means to engineer long-range spin-spin interactions, and has been realized in a variety of platforms, such as atoms in optical cavities and trapped ions. To date, much of the experimental focus has been on the realization of long-range Ising models, but generalizations to other spin models are highly desirable. In this work, we explore a previously unappreciated connection between the realization of an XY model by off-resonant driving of single sideband of boson excitation (i.e.~a single-beam Mølmer-Sørensen scheme) and a boson-mediated Ising simulator in the presence of a transverse field. In particular, we show that these two schemes have the same effective Hamiltonian in suitably defined rotating frames, and analyze the emergent effective XY spin model through truncated Magnus series and numerical simulations. In addition to XY spin-spin interactions that can be non-perturbatively renormalized from the naive Ising spin-spin coupling constants, we find an effective transverse field that is dependent on the thermal energy of the bosons, as well as other spin-boson couplings that cause spin-boson entanglement not to vanish at any time. In the case of a boson-mediated Ising simulator with transverse field, we discuss the crossover from transverse-field Ising-like to XY-like spin behavior as a function of field strength.

Controllable arrays of ions and ultra-cold atoms can simulate complex many-body phenomena and may provide insights into unsolved problems in modern science. To this end, experimentally feasible protocols for quantifying the buildup of quantum correlations and coherence are needed, as performing full state tomography does not scale favorably with the number of particles. Here we develop and experimentally demonstrate such a protocol, which uses time reversal of the many-body dynamics to measure out-of-time-order correlation functions (OTOCs) in a long-range Ising spin quantum simulator with more than 100 ions in a Penning trap. By measuring a family of OTOCs as a function of a tunable parameter we obtain fine-grained information about the state of the system encoded in the multiple quantum coherence spectrum, extract the quantum state purity, and demonstrate the buildup of up to 8-body correlations. Future applications of this protocol could enable studies of many-body localization, quantum phase transitions, and tests of the holographic duality between quantum and gravitational systems.

The global coupling of few-level quantum systems ("spins") to a discrete set of bosonic modes is a key ingredient for many applications in quantum science, including large-scale entanglement generation, quantum simulation of the dynamics of long-range interacting spin models, and hybrid platforms for force and spin sensing. We present a general numerical framework for treating the out-of-equilibrium dynamics of such models based on matrix product states. Our approach applies for generic spin-boson systems: it treats any spatial and operator dependence of the two-body spin-boson coupling and places no restrictions on relative energy scales. We show that the full counting statistics of collective spin measurements and infidelity of quantum simulation due to spin-boson entanglement, both of which are difficult to obtain by other techniques, are readily calculable in our approach. We benchmark our method using a recently developed exact solution for a particular spin-boson coupling relevant to trapped ion quantum simulators. Finally, we show how decoherence can be incorporated within our framework using the method of quantum trajectories, and study the dynamics of an open-system spin-boson model with spatially non-uniform spin-boson coupling relevant for trapped atomic ion crystals in the presence of molecular ion impurities.

Quantum simulation of spin models can provide insight into complex problems that are difficult or impossible to study with classical computers. Trapped ions are an established platform for quantum simulation, but only systems with fewer than 20 ions have demonstrated quantum correlations. Here we study non-equilibrium, quantum spin dynamics arising from an engineered, homogeneous Ising interaction in a two-dimensional array of $^9$Be$^+$ ions in a Penning trap. We verify entanglement in the form of spin-squeezed states for up to 219 ions, directly observing 4.0$\pm$0.9 dB of spectroscopic enhancement. We also observe evidence of non-Gaussian, over-squeezed states in the full counting statistics. We find good agreement with ab-initio theory that includes competition between entanglement and decoherence, laying the groundwork for simulations of the transverse-field Ising model with variable-range interactions, for which numerical solutions are, in general, classically intractable.

Interacting quantum spin models are remarkably useful for describing different types of physical, chemical, and biological systems. Significant understanding of their equilibrium properties has been achieved to date, especially for the case of spin models with short-range couplings. However, progress towards the development of a comparable understanding in long-range interacting models, in particular out-of-equilibrium, remains limited. In a recent work, we proposed a semiclassical numerical method to study spin models, the discrete truncated Wigner approximation (DTWA), and demonstrated its capability to correctly capture the dynamics of one- and two-point correlations in one dimensional (1D) systems. Here we go one step forward and use the DTWA method to study the dynamics of correlations in 2D systems with many spins and different types of long-range couplings, in regimes where other numerical methods are generally unreliable. We compute spatial and time-dependent correlations for spin-couplings that decay with distance as a power-law and determine the velocity at which correlations propagate through the system. Sharp changes in the behavior of those velocities are found as a function of the power-law decay exponent. Our predictions are relevant for a broad range of systems including solid state materials, atom-photon systems and ultracold gases of polar molecules, trapped ions, Rydberg, and magnetic atoms. We validate the DTWA predictions for small 2D systems and 1D systems, but ultimately, in the spirt of quantum simulation, experiments will be needed to confirm our predictions for large 2D systems.

By applying complementary analytic and numerical methods, we investigate the dynamics of spin-$1/2$ XXZ models with variable-range interactions in arbitrary dimensions. The dynamics we consider is initiated from uncorrelated states that are easily prepared in experiments, and can be equivalently viewed as either Ramsey spectroscopy or a quantum quench. Our primary focus is the dynamical emergence of correlations and entanglement in these far-from-equilibrium interacting quantum systems: we characterize these correlations by the entanglement entropy, concurrence, and squeezing, which are inequivalent measures of entanglement corresponding to different quantum resources. In one spatial dimension, we show that the time evolution of correlation functions manifests a non-perturbative dynamic singularity. This singularity is characterized by a universal power-law exponent that is insensitive to small perturbations. Explicit realizations of these models in current experiments using polar molecules, trapped ions, Rydberg atoms, magnetic atoms, and alkaline-earth and alkali atoms in optical lattices, along with the relative merits and limitations of these different systems, are discussed.

Ising models, and the physical systems described by them, play a central role in generating entangled states for use in quantum metrology and quantum information. In particular, ultracold atomic gases, trapped ion systems, and Rydberg atoms realize long-ranged Ising models, which even in the absence of a transverse field can give rise to highly non-classical dynamics and long-range quantum correlations. In the first part of this paper, we present a detailed theoretical framework for studying the dynamics of such systems driven (at time t=0) into arbitrary unentangled non-equilibrium states, thus greatly extending and unifying the work of Ref. [1]. Specifically, we derive exact expressions for closed-time-path ordered correlation functions, and use these to study experimentally relevant observables, e.g. Bloch vector and spin-squeezing dynamics. In the second part, these correlation functions are then used to derive closed-form expressions for the dynamics of arbitrary spin-spin correlation functions in the presence of both T_1 (spontaneous spin relaxation/excitation) and T_2 (dephasing) type decoherence processes. Even though the decoherence is local, our solution reveals that the competition between Ising dynamics and T_1 decoherence gives rise to an emergent non-local dephasing effect, thereby drastically amplifying the degradation of quantum correlations. In addition to identifying the mechanism of this deleterious effect, our solution points toward a scheme to eliminate it via measurement-based coherent feedback.