J. K. Freericks

We examine the detailed scenario for implementing n-control-qubit Toffoli gates and select gates on ion-trap quantum computers, especially those that shuttle ions into interaction zones. We determine expected performance of these gates with realistic parameters for an ion-trap quantum computer and taking into account the time variation of the exchange integrals. This allows us to estimate the errors due to spin-phonon entanglement as well. While there are challenges with implementing these gates, because their performance always has some degree of error, they should be feasible on current hardware, but they may be too slow to be used efficiently in quantum codes on noisy intermediate scale quantum computers.

Interacting quantum systems illustrate complex phenomena including phase transitions to novel ordered phases. The universal nature of critical phenomena reduces their description to determining only the transition temperature and the critical exponents. Numerically calculating these results for systems in new universality classes is complicated due to critical slowing down, requiring increasing resources near the critical point. An alternative approach analytically continues the calculation onto the complex plane and determines the partition function via its zeros. Here we show how to robustly perform this analysis on noisy intermediate scale trapped ion quantum computers in a scalable manner, using the XXZ model as a prototype. We illustrate the transition from XY-like behavior to Ising-like behavior as a function of the anisotropy. While quantum computers cannot yet scale to the thermodynamic limit, our work provides a pathway to do so as hardware improves, allowing the determination of critical phenomena for systems that cannot be solved otherwise.

We introduce a bang-bang shortcut to adiabaticity for the Dicke model, which we implement via a 2-D array of trapped ions in a Penning trap with a spin-dependent force detuned close to the center-of-mass drumhead mode. Our focus is on employing this shortcut to create highly entangled states that can be used in high-precision metrology. We highlight that the performance of the bang-bang approach is comparable to standard preparation methods, but can be applied over a much shorter time frame. We compare these theoretical ideas with experimental data which serve as a first step towards realizing this theoretical procedure for generating multi-partite entanglement.

We use a self-assembled two-dimensional Coulomb crystal of $\sim 70$ ions in the presence of an external transverse field to engineer a simulator of the Dicke Hamiltonian, an iconic model in quantum optics which features a quantum phase transition between a superradiant/ferromagnetic and a normal/paramagnetic phase. We experimentally implement slow quenches across the quantum critical point and benchmark the dynamics and the performance of the simulator through extensive theory-experiment comparisons which show excellent agreement. The implementation of the Dicke model in fully controllable trapped ion arrays can open a path for the generation of highly entangled states useful for enhanced metrology and the observation of scrambling and quantum chaos in a many-body system.

We use a self-assembled two-dimensional Coulomb crystal of $\sim 70$ ions in the presence of an external transverse field to engineer a quantum simulator of the Dicke Hamiltonian. This Hamiltonian has spin and bosonic degrees of freedom which are encoded by two hyperfine states in each ion and the center of mass motional mode of the crystal, respectively. The Dicke model features a quantum critical point separating two distinct phases: the superradiant (ferromagnetic) and normal (paramagnetic) phases. We experimentally explore protocols that aim to adiabatically prepare the superradiant ground state, a spin-boson cat state with macroscopic phonon occupation, which is well-suited for enhanced metrology and quantum information processing. We start in the normal phase, with all spins aligned along a large transverse field and ramp down the field across the critical point following various protocols. We measure the spin observables, both experimentally and in our simulations to characterize the state of the system at the end of the ramp. We find that under current operating conditions an optimally designed ramp is not sufficient to achieve significant fidelity with the superradiant ground state. However, our theoretical investigation shows that slight modifications of experimental parameters, together with modest reductions in decoherence rates and thermal noise can increase the cat-state fidelity to $\sim 75\%$ for $N \sim 20$ spins. Our results open a path for the use of large ensembles of trapped ions as powerful quantum sensors and quantum computers.

Recent work proposed a variant on Ramsey interferometry for coupled spin-$1/2$ systems that directly measures the retarded spin-spin Green's function. We expand on that work by investigating nonequilibrium retarded spin-spin Green's functions within the transverse-field Ising model. We derive the lowest four spectral moments to understand the short-time behavior and we employ a Lehmann-like representation to determine the spectral behavior. We simulate a Ramsey protocol for a nonequilibrium quantum spin system that consists of a coherent superposition of the ground state and diabatically excited higher-energy states via a temporally ramped transverse magnetic field. We then apply the Ramsey spectroscopy protocol to the final Hamiltonian, which has a constant transverse field. The short-time behavior directly relates to Lieb-Robinson bounds for the transport of many-body correlations, while the long-time behavior relates to the excitation spectra of the Hamiltonian. Compressive sensing is employed in the data analysis to efficiently extract that spectra.

We model the bang-bang optimization protocol as a shortcut to adiabaticity in the ground-state preparation of an ion-trap-based quantum simulator. Compared to a locally adiabatic evolution, the bang-bang protocol produces a somewhat lower ground-state probability, but its implementation is so much simpler than the locally adiabatic approach, that it remains an excellent choice to use for maximizing ground-state preparation in systems that cannot be solved with conventional computers. We describe how one can optimize the shortcut and provide specific details for how it can be implemented with current ion-trap-based quantum simulators.

The goal of adiabatic ground-state preparation is to start a simple quantum system in its ground state and adiabatically evolve the Hamiltonian to a complex one, maintaining the ground state throughout the evolution. In ion-trap-based quantum simulations, coherence times are too short to allow for adiabatic evolution for large chains, so the system evolves diabatically, creating excitations to higher energy states. Because the probability for diabatic excitation depends exponentially on the excitation energy and because the thermal distribution also depends exponentially on the excitation energy, we investigate whether the diabatic excitation can create a thermal distribution; as this could serve as an alternative for creating thermal states of complex quantum systems without requiring contact with a heat bath. In this work, we explore this relationship and determine situations where diabatic excitation can approximately create such a thermal state.

Theoretical models of spins coupled to bosons provide a simple setting for studying a broad range of important phenomena in many-body physics, from virtually mediated interactions to decoherence and thermalization. In many atomic, molecular, and optical systems, such models also underlie the most successful attempts to engineer strong, long-ranged interactions for the purpose of entanglement generation. Especially when the coupling between the spins and bosons is strong---such that it cannot be treated perturbatively---the properties of such models are extremely challenging to calculate theoretically. Here, exact analytical expressions for nonequilibrium spin-spin correlation functions are derived for a specific model of spins coupled to bosons. The spatial structure of the coupling between spins and bosons is completely arbitrary, and thus the solution can be applied to systems in any number of dimensions. The explicit and nonperturbative inclusion of the bosons enables the study of entanglement generation (in the form of spin squeezing) even when the bosons are driven strongly and near-resonantly, and thus provides a quantitative view of the breakdown of adiabatic elimination that inevitably occurs as one pushes towards the fastest entanglement generation possible. The solution also helps elucidate the effect of finite temperature on spin squeezing. The model considered is relevant to a variety of atomic, molecular, and optical systems, such as atoms in cavities or trapped ions. As an explicit example, the results are used to quantify phonon effects in trapped ion quantum simulators, which are expected to become increasingly important as these experiments push towards larger numbers of ions.

One of the goals in quantum simulation is to adiabatically generate the ground state of a complicated Hamiltonian by starting with the ground state of a simple Hamiltonian and slowly evolving the system to the complicated one. If the evolution is adiabatic and the initial and final ground states are connected due to having the same symmetry, then the simulation will be successful. But in most experiments, adiabatic simulation is not possible because it would take too long, and the system has some level of diabatic excitation. In this work, we quantify the extent of the diabatic excitation even if we do not know {\it a priori} what the complicated ground state is. Since many quantum simulator platforms, like trapped ions, can measure the probabilities to be in a product state, we describe techniques that can employ these measurements to estimate the probability of being in the ground state of the system after the diabatic evolution. These techniques do not require one to know any properties about the Hamiltonian itself, nor to calculate its eigenstate properties. All the information is derived by analyzing the product-state measurements as functions of time.

One of the challenges with quantum simulation in ion traps is that the effective spin-spin exchange couplings are not uniform across the lattice. This can be particularly important in Penning trap realizations where the presence of an ellipsoidal boundary at the edge of the trap leads to dislocations in the crystal. By adding an additional anharmonic potential to better control interion spacing, and a triangular shaped rotating wall potential to reduce the appearance of dislocations, one can achieve better uniformity of the ionic positions. In this work, we calculate the axial phonon frequencies and the spin-spin interactions driven by a spin-dependent optical dipole force, and discuss what effects the more uniform ion spacing has on the spin simulation properties of Penning trap quantum simulators. Indeed, we find the spin-spin interactions behave more like a power law for a wide range of parameters.

We develop the theory to describe the equilibrium ion positions and phonon modes for a trapped ion quantum simulator in an oblate Paul trap that creates two-dimensional Coulomb crystals in a triangular lattice. By coupling the internal states of the ions to laser beams propagating along the symmetry axis, we study the effective Ising spin-spin interactions that are mediated via the axial phonons and are less sensitive to ion micromotion. We find that the axial mode frequencies permit the programming of Ising interactions with inverse power law spin-spin couplings that can be tuned from uniform to $r^{-3}$ with DC voltages. Such a trap could allow for interesting new geometrical configurations for quantum simulations on moderately sized systems including frustrated magnetism on triangular lattices or Aharonov-Bohm effects on ion tunneling. The trap also incorporates periodic boundary conditions around loops which could be employed to examine time crystals.

The Coulomb repulsion between ions in a linear Paul trap give rise to anharmonic terms in the potential energy when expanded about the equilibrium positions. We examine the effect of these anharmonic terms on the accuracy of a quantum simulator made from trapped ions. To be concrete, we consider a linear chain of $\text{Yb}^{171+}$ ions stabilized close to the zigzag transition. We find that for typical experimental temperatures, frequencies change by no more than a factor of $0.01\%$ due to the anharmonic couplings. Furthermore, shifts in the effective spin-spin interactions (driven by a spin-dependent optical dipole force) also tend to be small for detunings to the blue of the transverse center-of-mass frequency. However, detuning the spin interactions near other frequencies can lead to nonnegligible anharmonic contributions to the effective spin-spin interactions. We also examine an odd behavior exhibited by the harmonic spin-spin interactions for a range of intermediate detunings, where nearest neighbor spins with a larger spatial separation on the ion chain interact more strongly than nearest neighbors with a smaller spatial separation.

Due to the experimental time constraints of state of the art quantum simulations with trapped ions, the direct preparation of the ground state by adiabatically ramping the field of a transverse field Ising model becomes more and more difficult as the number of particles increase. We propose a spectroscopy protocol that intentionally creates excitations through diabatic ramping of the transverse field and measures a low-noise observable as a function of time for a constant field to reveal the structure of the coherent dynamics of the resulting many-body states. To simulate the experimental data, noise from counting statistics and decoherence error are added. Compressive sensing is then applied to Fourier transform the simulated data into the frequency domain and extract the the low-lying energy excitation spectrum. By using compressive sensing, the amount of data in time needed to extract this energy spectrum is sharply reduced making such experiments feasible with current technology.

We derive the normal modes for a rotating Coulomb ion crystal in a Penning trap, quantize the motional degrees of freedom, and illustrate how they can by driven by a spin-dependent optical dipole force to create a quantum spin simulator on a triangular lattice with hundreds of spins. The analysis for the axial modes (oscillations perpendicular to the two-dimensional crystal plane) follow a standard normal-mode analysis, while the remaining planar modes are more complicated to analyze because they have velocity-dependent forces in the rotating frame. After quantizing the normal modes into phonons, we illustrate some of the different spin-spin interactions that can be generated by entangling the motional degrees of freedom with the spin degrees of freedom via a spin-dependent optical dipole force. In addition to the well-known power-law dependence of the spin-spin interactions when driving the axial modes blue of phonon band, we notice certain parameter regimes in which the level of frustration between the spins can be engineered by driving the axial or planar phonon modes at different energies. These systems may allow for the analog simulation of quantum spin glasses with large numbers of spins.

Frustration, or the competition between interacting components of a network, is often responsible for the complexity of many body systems, from social and neural networks to protein folding and magnetism. In quantum magnetic systems, frustration arises naturally from competing spin-spin interactions given by the geometry of the spin lattice or by the presence of long-range antiferromagnetic couplings. Frustrated magnetism is a hallmark of poorly understood systems such as quantum spin liquids, spin glasses and spin ices, whose ground states are massively degenerate and can carry high degrees of quantum entanglement. The controlled study of frustrated magnetism in materials is hampered by short dynamical time scales and the presence of impurities, while numerical modeling is generally intractable when dealing with dynamics beyond N~30 particles. Alternatively, a quantum simulator can be exploited to directly engineer prescribed frustrated interactions between controlled quantum systems, and several small-scale experiments have moved in this direction. In this article, we perform a quantum simulation of a long-range antiferromagnetic quantum Ising model with a transverse field, on a crystal of up to N = 16 trapped Yb+ atoms. We directly control the amount of frustration by continuously tuning the range of interaction and directly measure spin correlation functions and their dynamics through spatially-resolved spin detection. We find a pronounced dependence of the magnetic order on the amount of frustration, and extract signatures of quantum coherence in the resulting phases.

Linear Paul traps have been used recently to simulate the transverse field Ising model with long-range spin-spin couplings. We study the intrinsic effects of phonon creation (from the initial phonon ground state) on the spin-state probability and spin entanglement for such quantum spin simulators. While it has often been assumed that phonon effects are benign because they play no role in the pure Ising model, they can play a significant role when a transverse field is added to the model. We use a many-body factorization of the quantum time-evolution operator of the system, adiabatic perturbation theory and exact numerical integration of the Schrödinger equation in a truncated spin-phonon Hilbert space followed by a tracing out of the phonon degrees of freedom to study this problem. We find that moderate phonon creation often makes the probabilities of different spin states behave differently from the static spin Hamiltonian. In circumstances in which phonon creation is minor, the spin dynamics state probabilities converge to the static spin Hamiltonian prediction at the cost of reducing the spin entanglement. We show how phonon creation can severely impede the observation of kink transitions in frustrated spin systems when the number of ions increases. Many of our results also have implications for quantum simulation in a Penning trap.

A quantum simulator is a well controlled quantum system that can simulate the behavior of another quantum system which may require exponentially large classical computing resources to understand otherwise. In the 1980s, Feynman proposed the use of quantum logic gates on a standard controllable quantum system to efficiently simulate the behavior of a model Hamiltonian. Recent experiments using trapped ions and neutral atoms have realized quantum simulation of Ising model in presence of external magnetic fields, and showed almost arbitrary control in generating non-trivial Ising coupling patterns. Here we use laser-cooled trapped 171-Yb+ ions to simulate the emergence of magnetism in a system of interacting spins by implementing a fully-connected non-uniform ferromagnetic Ising model in a transverse magnetic field. To link this quantum simulation to condensed matter physics, we measure scalable correlation functions and order parameters appropriate for the description of larger systems, such as various moments of the magnetization. By increasing the Ising coupling strengths compared with the external field, the crossover from paramagnetism to ferromagnetic order sharpens as the system is scaled up from N = 2 to 9 trapped ion spins. This points toward the onset of a quantum phase transition that should become infinitely sharp as the system approaches the macroscopic scale. We compare the measured ground state order to theory, which may become intractable for non-uniform Ising couplings as the number of spins grows beyond 20- 30 and even NP complete for a fully-connected frustrated Ising model, making this experiment an important benchmark for large-scale quantum simulation.

We perform a quantum simulation of the Ising model with a transverse field using a collection of three trapped atomic ion spins. By adiabatically manipulating the Hamiltonian, we directly probe the ground state for a wide range of fields and form of the Ising couplings, leading to a phase diagram of magnetic order in this microscopic system. The technique is scalable to much larger numbers of trapped ion spins, where phase transitions approaching the thermodynamic limit can be studied in cases where theory becomes intractable.