Zohreh Davoudi

Simulating out-of-equilibrium dynamics of quantum field theories in nature is challenging with classical methods, but is a promising application for quantum computers. Unfortunately, simulating interacting bosonic fields involves a high boson-to-qubit encoding overhead. Furthermore, when mapping to qubits, the infinite-dimensional Hilbert space of bosons is necessarily truncated, with truncation errors that grow with energy and time. A qubit-based quantum computer, augmented with an active bosonic register, and with qubit, bosonic, and mixed qubit-boson quantum gates, offers a more powerful platform for simulating bosonic theories. We demonstrate this capability experimentally in a hybrid analog-digital trapped-ion quantum computer, where qubits are encoded in the internal states of the ions, and the bosons in the ions' motional states. Specifically, we simulate nonequilibrium dynamics of a (1+1)-dimensional Yukawa model, a simplified model of interacting nucleons and pions, and measure fermion- and boson-occupation-state probabilities. These dynamics populate high bosonic-field excitations starting from an empty state, and the experimental results capture well such high-occupation states. This simulation approaches the regime where classical methods become challenging, bypasses the need for a large qubit overhead, and removes truncation errors. Our results, therefore, open the way to achieving demonstrable quantum advantage in qubit-boson quantum computing.

Recent advances in quantum technologies have enabled quantum simulation of gauge theories -- some of the most fundamental frameworks of nature -- in regimes far from equilibrium, where classical computation is severely limited. These simulators, primarily based on neutral atoms, trapped ions, and superconducting circuits, hold the potential to address long-standing questions in nuclear, high-energy, and condensed-matter physics, and may ultimately allow first-principles studies of matter evolution in settings ranging from the early universe to high-energy collisions. Research in this rapidly growing field is also driving the convergence of concepts across disciplines and uncovering new phenomena. In this Review, we highlight recent experimental and theoretical developments, focusing on phenomena accessible in current and near-term quantum simulators, including particle production and string breaking, collision dynamics, thermalization, ergodicity breaking, and dynamical quantum phase transitions. We conclude by outlining promising directions for future research and opportunities enabled by available quantum hardware.

We present a digital quantum computation of two-hadron scattering in a $Z_2$ lattice gauge theory in 1+1 dimensions. We prepare well-separated single-particle wave packets with desired momentum-space wavefunctions, and simulate their collision through digitized time evolution. Multiple hadronic wave packets can be produced using the efficient, systematically improvable algorithm of this work, achieving high fidelity with the target initial state. Specifically, employing a trapped-ion quantum computer (IonQ Forte), we prepare up to three meson wave packets using 11 and 27 system qubits, and simulate collision dynamics of two meson wave packets for the smaller system. Results for local observables are consistent with numerical simulations at early times, but decoherence effects limit evolution into long times. We demonstrate the critical role of high-fidelity initial states for precision measurements of state-sensitive observables, such as $S$-matrix elements. Our work establishes the potential of quantum computers in simulating hadron-scattering processes in strongly interacting gauge theories.

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.

Simulating non-equilibrium phenomena in strongly-interacting quantum many-body systems, including thermalization, is a promising application of near-term and future quantum computation. By performing experiments on a digital quantum computer consisting of fully-connected optically-controlled trapped ions, we study the role of entanglement in the thermalization dynamics of a $Z_2$ lattice gauge theory in 2+1 spacetime dimensions. Using randomized-measurement protocols, we efficiently learn a classical approximation of non-equilibrium states that yields the gap-ratio distribution and the spectral form factor of the entanglement Hamiltonian. These observables exhibit universal early-time signals for quantum chaos, a prerequisite for thermalization. Our work, therefore, establishes quantum computers as robust tools for studying universal features of thermalization in complex many-body systems, including in gauge theories.

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.

Quantum simulation holds promise of enabling a complete description of high-energy scattering processes rooted in gauge theories of the Standard Model. A first step in such simulations is preparation of interacting hadronic wave packets. To create the wave packets, one typically resorts to adiabatic evolution to bridge between wave packets in the free theory and those in the interacting theory, rendering the simulation resource intensive. In this work, we construct a wave-packet creation operator directly in the interacting theory to circumvent adiabatic evolution, taking advantage of resource-efficient schemes for ground-state preparation, such as variational quantum eigensolvers. By means of an ansatz for bound mesonic excitations in confining gauge theories, which is subsequently optimized using classical or quantum methods, we show that interacting mesonic wave packets can be created efficiently and accurately using digital quantum algorithms that we develop. Specifically, we obtain high-fidelity mesonic wave packets in the $Z_2$ and $U(1)$ lattice gauge theories coupled to fermionic matter in 1+1 dimensions. Our method is applicable to both perturbative and non-perturbative regimes of couplings. The wave-packet creation circuit for the case of the $Z_2$ lattice gauge theory is built and implemented on the Quantinuum H1-1 trapped-ion quantum computer using 13 qubits and up to 308 entangling gates. The fidelities agree well with classical benchmark calculations after employing a simple symmetry-based noise-mitigation technique. This work serves as a step toward quantum computing scattering processes in quantum chromodynamics.

Strongly-coupled gauge theories far from equilibrium may exhibit unique features that could illuminate the physics of the early universe and of hadron and ion colliders. Studying real-time phenomena has proven challenging with classical-simulation methods, but is a natural application of quantum simulation. To demonstrate this prospect, we quantum compute non-equal time correlation functions and perform entanglement tomography of non-equilibrium states of a simple lattice gauge theory, the Schwinger model, using a trapped-ion quantum computer by IonQ Inc. As an ideal target for near-term devices, a recently-predicted [Zache et al., Phys. Rev. Lett. 122, 050403 (2019)] dynamical quantum phase transition in this model is studied by preparing, quenching, and tracking the subsequent non-equilibrium dynamics in three ways: i) overlap echos signaling dynamical transitions, ii) non-equal time correlation functions with an underlying topological nature, and iii) the entanglement structure of non-equilibrium states, including entanglement Hamiltonians. These results constitute the first observation of a dynamical quantum phase transition in a lattice gauge theory on a quantum computer, and are a first step toward investigating topological phenomena in nuclear and high-energy physics using quantum technologies.

Tracking the dynamics of physical systems in real time is a prime application of digital quantum computers. Using a trapped-ion system with up to six qubits, we simulate the real-time dynamics of a lattice gauge theory in 1+1 dimensions, i.e., the lattice Schwinger model, and demonstrate non-perturbative effects such as pair creation for times much longer than previously accessible. We study the gate requirement of two formulations of the model using the Suzuki-Trotter product formula, as well as the trade-off between errors from the ordering of the Hamiltonian terms, the Trotter step size, and experimental imperfections. To mitigate experimental errors, a recent symmetry-protection protocol for suppressing coherent errors and a symmetry-inspired post-selection scheme are applied. This work demonstrates the integrated theoretical, algorithmic, and experimental approach that is essential for efficient simulation of lattice gauge theories and other complex physical systems.

Trapped-ion quantum simulators, in analog and digital modes, are considered a primary candidate to achieve quantum advantage in quantum simulation and quantum computation. The underlying controlled ion-laser interactions induce all-to-all two-spin interactions via the collective modes of motion through Cirac-Zoller or Molmer-Sorensen schemes, leading to effective two-spin Hamiltonians, as well as two-qubit entangling gates. In this work, the Molmer-Sorensen scheme is extended to induce three-spin interactions via tailored first- and second-order spin-motion couplings. The scheme enables engineering single-, two-, and three-spin interactions, and can be tuned via an enhanced protocol to simulate purely three-spin dynamics. Analytical results for the effective evolution are presented, along with detailed numerical simulations of the full dynamics to support the accuracy and feasibility of the proposed scheme for near-term applications. With a focus on quantum simulation, the advantage of a direct analog implementation of three-spin dynamics is demonstrated via the example of matter-gauge interactions in the U(1) lattice gauge theory within the quantum link model. The mapping of degrees of freedom and strategies for scaling the three-spin scheme to larger systems, are detailed, along with a discussion of the expected outcome of the simulation of the quantum link model given realistic fidelities in the upcoming experiments. The applications of the three-spin scheme go beyond the lattice gauge theory example studied here and include studies of static and dynamical phase diagrams of strongly interacting condensed-matter systems modeled by two- and three-spin Hamiltonians.

Quantum field theories are the cornerstones of modern physics, providing relativistic and quantum mechanical descriptions of physical systems at the most fundamental level. Simulating real-time dynamics within these theories remains elusive in classical computing. This provides a unique opportunity for quantum simulators, which hold the promise of revolutionizing our simulation capabilities. Trapped-ion systems are successful quantum-simulator platforms for quantum many-body physics and can operate in digital, or gate-based, and analog modes. Inspired by the progress in proposing and realizing quantum simulations of a number of relativistic quantum field theories using trapped-ion systems, and by the hybrid analog-digital proposals for simulating interacting boson-fermion models, we propose hybrid analog-digital quantum simulations of selected quantum field theories, taking recent developments to the next level. On one hand, the semi-digital nature of this proposal offers more flexibility in engineering generic model interactions compared with a fully-analog approach. On the other hand, encoding the bosonic fields onto the phonon degrees of freedom of the trapped-ion system allows a more efficient usage of simulator resources, and a more natural implementation of intrinsic quantum operations in such platforms. This opens up new ways for simulating complex dynamics of e.g., Abelian and non-Abelian gauge theories, by combining the benefits of digital and analog schemes.

Gauge field theories play a central role in modern physics and are at the heart of the Standard Model of elementary particles and interactions. Despite significant progress in applying classical computational techniques to simulate gauge theories, it has remained a challenging task to compute the real-time dynamics of systems described by gauge theories. An exciting possibility that has been explored in recent years is the use of highly-controlled quantum systems to simulate, in an analog fashion, properties of a target system whose dynamics are difficult to compute. Engineered atom-laser interactions in a linear crystal of trapped ions offer a wide range of possibilities for quantum simulations of complex physical systems. Here, we devise practical proposals for analog simulation of simple lattice gauge theories whose dynamics can be mapped onto spin-spin interactions in any dimension. These include 1+1D quantum electrodynamics, 2+1D Abelian Chern-Simons theory coupled to fermions, and 2+1D pure Z2 gauge theory. The scheme proposed, along with the optimization protocol applied, will have applications beyond the examples presented in this work, and will enable scalable analog quantum simulation of Heisenberg spin models in any number of dimensions and with arbitrary interaction strengths.