Alexey V. Gorshkov

While quantum simulators promise to explore quantum many-body physics beyond classical computation, their capabilities are limited by the available native interactions in the hardware. On many platforms, accessible Hamiltonians are largely restricted to one- and two-body interactions, limiting access to multi-body Hamiltonians and to systems governed by simultaneous, non-commuting interaction terms that are central to condensed matter, quantum chemistry, and high-energy physics. We introduce and experimentally demonstrate a hybrid digital-analog protocol that overcomes these limitations by embedding analog evolution between shallow entangling-gate layers. This method produces effective Hamiltonians with simultaneous non-commuting three- and four-body interactions that are generated non-perturbatively and without Trotter error -- capabilities not practically attainable on near-term hardware using purely digital or purely analog schemes. We implement our scheme on a trapped-ion quantum processor and use it to realize a topological spin chain exhibiting prethermal strong zero modes persisting at high temperature, as well as models featuring three- and four-body interactions. Our hardware-agnostic and scalable method opens new routes to realizing complex many-body physics across quantum platforms.

Analog quantum simulators have provided key insights into quantum many-body dynamics. However, in such systems, both coherent and incoherent errors limit their scalability, hindering simulations in regimes that challenge classical simulations. In this work, we introduce an error mitigation technique that addresses and effectively suppresses a key source of error in leading simulator platforms: shot-to-shot fluctuations in the parameters for the Hamiltonian governing the system dynamics. We rigorously prove that amplifying this shot-to-shot noise and extrapolating to the zero-noise limit recovers noiseless results for realistic noise distributions. Experimentally, we demonstrate this technique on a 27-ion trapped-ion quantum simulator, extending the two-qubit exchange oscillation lifetime threefold. Numerically, we predict a significant enhancement in the effective many-body coherence time for Rydberg atom arrays under realistic conditions. Our scheme provides a possible route towards extending the effective coherence time in analog quantum experiments, enabling deeper explorations of quantum many-body dynamics.

We propose a cavity-mediated gate between two transmon qubits or other nonlinear superconducting elements. The gate is realized by driving both qubits at a frequency that is near-resonant with the frequency of the cavity. Since both qubits are subject to a cross-resonant drive, we call this gate a cross-cross-resonance gate. In close analogy with gates between trapped-ion qubits, in phase space, the state of the cavity makes a circle whose area depends on the state of the two qubits, realizing a controlled-phase gate. We propose two schemes for canceling the dominant error, which is the dispersive coupling. We also show that this cross-cross-resonance gate allows one to realize simultaneous gates between multiple pairs of qubits coupled via the same metamaterial composed of an array of coupled cavities or other linear mediators.

The liquid-vapor transition is a classic example of a discontinuous (first-order) phase transition. Such transitions underlie many phenomena in cosmology, nuclear and particle physics, and condensed-matter physics. They give rise to long-lived metastable states, whose decay can be driven by either thermal or quantum fluctuations. Yet, direct experimental observations of how these states collapse into a stable phase remain elusive in the quantum regime. Here, we use a trapped-ion quantum simulator to observe the real-time dynamics of ``bubble nucleation'' induced by quantum fluctuations. Bubbles are localized domains of the stable phase which spontaneously form, or nucleate, and expand as the system is driven across a discontinuous quantum phase transition. Implementing a mixed-field Ising spin model with tunable and time-dependent interactions, we track the microscopic evolution of the metastable state as the Hamiltonian parameters are varied in time with various speeds, bringing the system out of equilibrium. Site-resolved measurements reveal the emergence and evolution of finite-size quantum bubbles, providing direct insight into the mechanism by which the metastable phase decays. We also identify nonequilibrium scaling behavior near the transition, consistent with a generalized Kibble-Zurek mechanism. Our results demonstrate the power of quantum simulators to probe out-of-equilibrium many-body physics, including quantum bubble nucleation, a key feature of discontinuous quantum phase transitions, with application to studies of matter formation in the early universe.

Spontaneous particle-pair formation is a fundamental phenomenon in nature. It can, for example, appear when the potential energy between two particles increases with separation, as if they were connected by a tense string. Beyond a critical separation, new particle pairs can form, causing the string to break. String-breaking dynamics in quantum chromodynamics play a vital role in high-energy particle collisions and early universe evolution. Simulating string evolution and hadron formation is, therefore, a grand challenge in modern physics. Quantum simulators, well-suited for studying dynamics, are expected to outperform classical computing methods. However, the required experimental capabilities to simulate string-breaking dynamics have not yet been demonstrated, even for simpler models of the strong force. We experimentally probe, for the first time, the spatiotemporal dynamics of string-breaking in a (1+1)-dimensional $\mathbb{Z}_2$ lattice gauge theory using a fully programmable trapped-ion quantum simulator. We emulate external static charges and strings via site-dependent magnetic-field control enabled by a dual array of tightly focused laser beams targeting individual ions. First, we study how confinement affects isolated charges, finding that they freely spread without string tension but exhibit localized oscillations when tension is increased. Then, we observe and characterize string-breaking dynamics of a string stretched between two static charges after an abrupt increase in string tension. Charge pairs appear near the string edges and spread into the bulk, revealing a route to dynamical string-breaking distinct from the conventional Schwinger mechanism. Our work demonstrates that analog quantum simulators have achieved the necessary control to explore string-breaking dynamics, which may ultimately be relevant to nuclear and high-energy physics.

Studying high-energy collisions of composite particles, such as hadrons and nuclei, is an outstanding goal for quantum simulators. However, preparation of hadronic wave packets has posed a significant challenge, due to the complexity of hadrons and the precise structure of wave packets. This has limited demonstrations of hadron scattering on quantum simulators to date. Observations of confinement and composite excitations in quantum spin systems have opened up the possibility to explore scattering dynamics in spin models. In this article, we develop two methods to create entangled spin states corresponding to wave packets of composite particles in analog quantum simulators of Ising spin Hamiltonians. One wave-packet preparation method uses the blockade effect enabled by beyond-nearest-neighbor Ising spin interactions. The other method utilizes a quantum-bus-mediated exchange, such as the native spin-phonon coupling in trapped-ion arrays. With a focus on trapped-ion simulators, we numerically benchmark both methods and show that high-fidelity wave packets can be achieved in near-term experiments. We numerically study scattering of wave packets for experimentally realizable parameters in the Ising model and find inelastic-scattering regimes, corresponding to particle production in the scattering event, with prominent and distinct experimental signals. Our proposal, therefore, demonstrates the potential of observing inelastic scattering in near-term quantum simulators.

One of the most striking many-body phenomena in nature is the sudden change of macroscopic properties as the temperature or energy reaches a critical value. Such equilibrium transitions have been predicted and observed in two and three spatial dimensions, but have long been thought not to exist in one-dimensional (1D) systems. Fifty years ago, Dyson and Thouless pointed out that a phase transition in 1D can occur in the presence of long-range interactions, but an experimental realization has so far not been achieved due to the requirement to both prepare equilibrium states and realize sufficiently long-range interactions. Here we report on the first experimental demonstration of a finite-energy phase transition in 1D. We use the simple observation that finite-energy states can be prepared by time-evolving product initial states and letting them thermalize under the dynamics of a many-body Hamiltonian. By preparing initial states with different energies in a 1D trapped-ion quantum simulator, we study the finite-energy phase diagram of a long-range interacting quantum system. We observe a ferromagnetic equilibrium phase transition as well as a crossover from a low-energy polarized paramagnet to a high-energy unpolarized paramagnet in a system of up to $23$ spins, in excellent agreement with numerical simulations. Our work demonstrates the ability of quantum simulators to realize and study previously inaccessible phases at finite energy density.

Universality and scaling laws are hallmarks of equilibrium phase transitions and critical phenomena. However, extending these concepts to non-equilibrium systems is an outstanding challenge. Despite recent progress in the study of dynamical phases, the universality classes and scaling laws for non-equilibrium phenomena are far less understood than those in equilibrium. In this work, using a trapped-ion quantum simulator with single-spin resolution, we investigate the non-equilibrium nature of critical fluctuations following a quantum quench to the critical point. We probe the scaling of spin fluctuations after a series of quenches to the critical Hamiltonian of a long-range Ising model. With systems of up to 50 spins, we show that the amplitude and timescale of the post-quench fluctuations scale with system size with distinct universal critical exponents, depending on the quench protocol. While a generic quench can lead to thermal critical behavior, we find that a second quench from one critical state to another (i.e.~a double quench) results in a new universal non-equilibrium behavior, identified by a set of critical exponents distinct from their equilibrium counterparts. Our results demonstrate the ability of quantum simulators to explore universal scaling beyond equilibrium.

Simulating quantum systems is one of the most promising avenues to harness the computational power of quantum computers. However, hardware errors in noisy near-term devices remain a major obstacle for applications. Ideas based on the randomization of Suzuki-Trotter product formulas have been shown to be a powerful approach to reducing the errors of quantum simulation and lowering the gate count. In this paper, we study the performance of non-unitary simulation channels and consider the error structure of channels constructed from a weighted average of unitary circuits. We show that averaging over just a few simulation circuits can significantly reduce the Trotterization error for both single-step short-time and multi-step long-time simulations. We focus our analysis on two approaches for constructing circuit ensembles for averaging: (i) permuting the order of the terms in the Hamiltonian and (ii) applying a set of global symmetry transformations. We compare our analytical error bounds to empirical performance and show that empirical error reduction surpasses our analytical estimates in most cases. Finally, we test our method on an IonQ trapped-ion quantum computer accessed via the Amazon Braket cloud platform, and benchmark the performance of the averaging approach.

One-dimensional systems exhibiting a continuous symmetry can host quantum phases of matter with true long-range order only in the presence of sufficiently long-range interactions. In most physical systems, however, the interactions are short-ranged, hindering the emergence of such phases in one dimension. Here we use a one-dimensional trapped-ion quantum simulator to prepare states with long-range spin order that extends over the system size of up to $23$ spins and is characteristic of the continuous symmetry-breaking phase of matter. Our preparation relies on simultaneous control over an array of tightly focused individual-addressing laser beams, generating long-range spin-spin interactions. We also observe a disordered phase with frustrated correlations. We further study the phases at different ranges of interaction and the out-of-equilibrium response to symmetry-breaking perturbations. This work opens an avenue to study new quantum phases and out-of-equilibrium dynamics in low-dimensional systems.

Many-body open quantum systems balance internal dynamics against decoherence from interactions with an environment. Here, we explore this balance via random quantum circuits implemented on a trapped ion quantum computer, where the system evolution is represented by unitary gates with interspersed projective measurements. As the measurement rate is varied, a purification phase transition is predicted to emerge at a critical point akin to a fault-tolerent threshold. We probe the "pure" phase, where the system is rapidly projected to a deterministic state conditioned on the measurement outcomes, and the "mixed" or "coding" phase, where the initial state becomes partially encoded into a quantum error correcting codespace. We find convincing evidence of the two phases and show numerically that, with modest system scaling, critical properties of the transition clearly emerge.

String breaking is a central dynamical process in theories featuring confinement, where a string connecting two charges decays at the expense of the creation of new particle-antiparticle pairs. Here, we show that this process can also be observed in quantum Ising chains where domain walls get confined either by a symmetry-breaking field or by long-range interactions. We find that string breaking occurs, in general, as a two-stage process: First, the initial charges remain essentially static and stable. The connecting string, however, can become a dynamical object. We develop an effective description of this motion, which we find is strongly constrained. In the second stage, which can be severely delayed due to these dynamical constraints, the string finally breaks. We observe that the associated time scale can depend crucially on the initial separation between domain walls and can grow by orders of magnitude by changing the distance by just a few lattice sites. We discuss how our results generalize to one-dimensional confining gauge theories and how they can be made accessible in quantum simulator experiments such as Rydberg atoms or trapped ions.

We generalize past work on quantum sensor networks to show that, for $d$ input parameters, entanglement can yield a factor $\mathcal O(d)$ improvement in mean squared error when estimating an analytic function of these parameters. We show that the protocol is optimal for qubit sensors, and conjecture an optimal protocol for photons passing through interferometers. Our protocol is also applicable to continuous variable measurements, such as one quadrature of a field operator. We outline a few potential applications, including calibration of laser operations in trapped ion quantum computing.

We study the quasiparticle excitation and quench dynamics of the one-dimensional transverse-field Ising model with power-law ($1/r^α$) interactions. We find that long-range interactions give rise to a confining potential, which couples pairs of domain walls (kinks) into bound quasiparticles, analogous to mesonic bound states in high-energy physics. We show that these quasiparticles have signatures in the dynamics of order parameters following a global quench and the Fourier spectrum of these order parameters can be expolited as a direct probe of the masses of the confined quasiparticles. We introduce a two-kink model to qualitatively explain the phenomenon of long-range-interaction induced confinement, and to quantitatively predict the masses of the bound quasiparticles. Furthermore, we illustrate that these quasiparticle states can lead to slow thermalization of one-point observables for certain initial states. Our work is readily applicable to current trapped-ion experiments.

The study of quantum phase transitions requires the preparation of a many-body system near its ground state, a challenging task for many experimental systems. The measurement of quench dynamics, on the other hand, is now a routine practice in most cold atom platforms. Here we show that quintessential ingredients of quantum phase transitions can be probed directly with quench dynamics in integrable and nearly integrable systems. As a paradigmatic example, we study global quench dynamics in a transverse-field Ising model with either short-range or long-range interactions. When the model is integrable, we discover a new dynamical critical point with a non-analytic signature in the short-range correlators. The location of the dynamical critical point matches that of the quantum critical point and can be identified using a finite-time scaling method. We extend this scaling picture to systems near integrability and demonstrate the continued existence of a dynamical critical point detectable at prethermal time scales. Therefore, our method can be used to approximately locate the quantum critical point. The scaling method is also relevant to experiments with finite time and system size, and our predictions are testable in near-term experiments with trapped ions and Rydberg atoms.

Exactly solvable models have played an important role in establishing the sophisticated modern understanding of equilibrium many-body physics. And conversely, the relative scarcity of solutions for non-equilibrium models greatly limits our understanding of systems away from thermal equilibrium. We study a family of non-equilibrium models, some of which can be viewed as dissipative analogues of the transverse-field Ising model, in that an effectively classical Hamiltonian is frustrated by dissipative processes that drive the system toward states that do not commute with the Hamiltonian. Surprisingly, a broad and experimentally relevant subset of these models can be solved efficiently in any number of spatial dimensions. We leverage these solutions to prove a no-go theorem on steady-state phase transitions in a many-body model that can be realized naturally with Rydberg atoms or trapped ions, and to compute the effects of decoherence on a canonical trapped-ion-based quantum computation architecture.

Motivated by recent trapped-ion quantum simulation experiments, we carry out a comprehensive study of the phase diagram of a spin-1 chain with XXZ-type interactions that decay as $1/r^α$, using a combination of finite and infinite-size DMRG calculations, spin-wave analysis, and field theory. In the absence of long-range interactions, varying the spin-coupling anisotropy leads to four distinct phases: a ferromagnetic Ising phase, a disordered XY phase, a topological Haldane phase, and an antiferromagnetic Ising phase. If long-range interactions are antiferromagnetic and thus frustrated, we find primarily a quantitative change of the phase boundaries. On the other hand, ferromagnetic (non-frustrated) long-range interactions qualitatively impact the entire phase diagram. Importantly, for $α\lesssim3$, long-range interactions destroy the Haldane phase, break the conformal symmetry of the XY phase, give rise to a new phase that spontaneously breaks a $U(1)$ continuous symmetry, and introduce an exotic tricritical point with no direct parallel in short-range interacting spin chains. We show that the main signatures of all five phases found could be observed experimentally in the near future.

Continuous symmetry breaking (CSB) in low-dimensional systems, forbidden by the Mermin-Wagner theorem for short-range interactions, may take place in the presence of slowly decaying long-range interactions. Nevertheless, there is no stringent bound on how slowly interactions should decay to give rise to CSB in 1D quantum systems at zero temperature. Here, we study a long-range interacting spin chain with $U(1)$ symmetry and power-law interactions $V(r)\sim1/r^α$, directly relevant to ion-trap experiments. Using bosonization and renormalization group theory, we find CSB for $α$ smaller than a critical exponent $α_c(\le 3)$ depending on the microscopic parameters of the model. Furthermore, the transition from the gapless XY phase to the gapless CSB phase is mediated by the breaking of conformal symmetry due to long-range interactions, and is described by a new universality class akin to the Berezinskii-Kosterlitz-Thouless transition. Our analytical findings are in good agreement with a numerical calculation. Signatures of the CSB phase should be accessible in existing trapped-ion experiments.

Topological phases of matter are primarily studied in systems with short-range interactions. In nature, however, non-relativistic quantum systems often exhibit long-range interactions. Under what conditions topological phases survive such interactions, and how they are modified when they do, is largely unknown. By studying the symmetry-protected topological phase of an antiferromagnetic spin-1 chain with $1/r^α$ interactions, we show that two very different outcomes are possible, depending on whether or not the interactions are frustrated. While non-frustrated long-range interactions can destroy the topological phase for $α\lesssim3$, the topological phase survives frustrated interactions for all $α>0$. Our conclusions are based on strikingly consistent results from large-scale matrix-product-state simulations and effective-field-theory calculations, and we expect them to hold for more general interacting spin systems. The models we study can be naturally realized in trapped-ion quantum simulators, opening the prospect for experimental investigation of the issues confronted here.

The realization of a scalable quantum information processor has emerged over the past decade as one of the central challenges at the interface of fundamental science and engineering. Much progress has been made towards this goal. Indeed, quantum operations have been demonstrated on several trapped ion qubits, and other solid-state systems are approaching similar levels of control. Extending these techniques to achieve fault-tolerant operations in larger systems with more qubits remains an extremely challenging goal, in part, due to the substantial technical complexity of current implementations. Here, we propose and analyze an architecture for a scalable, solid-state quantum information processor capable of operating at or near room temperature. The architecture is applicable to realistic conditions, which include disorder and relevant decoherence mechanisms, and includes a hierarchy of control at successive length scales. Our approach is based upon recent experimental advances involving Nitrogen-Vacancy color centers in diamond and will provide fundamental insights into the physics of non-equilibrium many-body quantum systems. Additionally, the proposed architecture may greatly alleviate the stringent constraints, currently limiting the realization of scalable quantum processors.