Alexander Schuckert

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.

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.

We study the ground-state phase diagram of a one-dimensional $\mathbb{Z}_2$ lattice gauge theory coupled to soft-core bosonic matter at unit filling, inspired by the Higgs sector of the standard model. Through a combination of analytical perturbative approaches, exact diagonalization, and density-matrix-renormalization-group simulations, we uncover a rich phase diagram driven by gauge-field-mediated resonant pair hopping and the confinement of single particles. The pair hopping results in a bunching state with superextensive energy and macroscopic particle number fluctuations at strong electric field strengths and weak on-site interactions. The bunching state crosses over into a pair superfluid phase as the on-site interaction increases, characterized by a finite superfluid density and powerlaw-decaying pair correlations. At large on-site interaction strengths and driven by effective interactions induced by the gauge constraint, the superfluid transitions into an incompressible pair Mott insulator phase. At weak field strengths and on-site interactions, we find a plasma-like region, where single bosons exhibit large short-range correlations and the ground state is composed almost equally of states with even and odd local boson occupation. The presence of a bunching state with large number fluctuations, which is difficult to study using classical numerics, motivates experimental realizations in hybrid boson-qubit quantum simulation platforms such as circuit QED, neutral atoms, and trapped ions. Our findings highlight the rich interplay between gauge fields and soft-core bosonic matter.

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.

We develop a hybrid oscillator-qubit processor framework for quantum simulation of strongly correlated fermions and bosons that avoids the boson-to-qubit mapping overhead encountered in qubit hardware. This framework gives exact decompositions of particle interactions such as density-density terms and gauge-invariant hopping, as well as approximate methods based on the Baker-Campbell Hausdorff formulas including the magnetic field term for the $U(1)$ quantum link model in $(2+1)$D. We use this framework to show how to simulate dynamics using Trotterisation, perform ancilla-free partial error detection using Gauss's law, measure non-local observables, estimate ground state energies using a oscillator-qubit variational quantum eigensolver as well as quantum signal processing, and we numerically study the influence of hardware errors in circuit QED experiments. To show the advantages over all-qubit hardware, we perform an end-to-end comparison of the gate complexity for the gauge-invariant hopping term and find an improvement of the asymptotic scaling with the boson number cutoff $S$ from $\mathcal{O}(\log(S)^2)$ to $\mathcal{O}(1)$ in our framework as well as, for bosonic matter, a constant factor improvement of better than $10^4$. We also find an improvement from $\mathcal{O}(\log(S))$ to $\mathcal{O}(1)$ for the $U(1)$ magnetic field term. While our work focusses on an implementation in superconducting hardware, our framework can also be used in trapped ion, and neutral atom hardware. This work establishes digital quantum simulation with hybrid oscillator-qubit hardware as a viable and advantageous method for the study of qubit-boson models in materials science, chemistry, and high-energy physics.

Quantum computing with discrete variable (DV, qubit) hardware is approaching the large scales necessary for computations beyond the reach of classical computers. However, important use cases such as quantum simulations of physical models containing bosonic modes, and quantum error correction are challenging for DV-only systems. Separately, hardware containing native continuous-variable (CV, oscillator) systems has received attention as an alternative approach, yet the universal control of such systems is non-trivial. In this work, we show that hybrid CV-DV hardware offers a great advantage in meeting these challenges, offering a powerful computational paradigm that inherits the strengths of both DV and CV processors. We provide a pedagogical introduction to CV-DV systems and the multiple abstraction layers needed to produce a full software stack connecting applications to hardware. We present a variety of new hybrid CV-DV compilation techniques, algorithms, and applications, including the extension of quantum signal processing concepts to CV-DV systems and strategies to simulate systems of interacting spins, fermions, and bosons. To facilitate the development of hybrid CV-DV processor systems, we introduce formal Abstract Machine Models and Instruction Set Architectures -- essential abstractions that enable developers to formulate applications, compile algorithms, and explore the potential of current and future hardware for realizing fault-tolerant circuits, modules, and processors. Hybrid CV-DV quantum computations are beginning to be performed in superconducting, trapped ion, and neutral atom platforms, and large-scale experiments are set to be demonstrated in the near future. We present a timely and comprehensive guide to this relatively unexplored yet promising approach to quantum computation and providing an architectural backbone to guide future development.

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.

Calculating the equilibrium properties of condensed matter systems is one of the promising applications of near-term quantum computing. Recently, hybrid quantum-classical time-series algorithms have been proposed to efficiently extract these properties from a measurement of the Loschmidt amplitude $\langle ψ| e^{-i \hat H t}|ψ\rangle$ from initial states $|ψ\rangle$ and a time evolution under the Hamiltonian $\hat H$ up to short times $t$. In this work, we study the operation of this algorithm on a present-day quantum computer. Specifically, we measure the Loschmidt amplitude for the Fermi-Hubbard model on a $16$-site ladder geometry (32 orbitals) on the Quantinuum H2-1 trapped-ion device. We assess the effect of noise on the Loschmidt amplitude and implement algorithm-specific error mitigation techniques. By using a thus-motivated error model, we numerically analyze the influence of noise on the full operation of the quantum-classical algorithm by measuring expectation values of local observables at finite energies. Finally, we estimate the resources needed for scaling up the algorithm.

Preparing finite temperature states in quantum simulators of spin systems, such as trapped ions or Rydberg atoms in optical tweezers, is challenging due to their almost perfect isolation from the environment. Here, we show how finite-temperature observables can be obtained with an algorithm motivated from the Jarzynski equality and equivalent to the one in Lu, Banuls and Cirac, PRX Quantum 2, 020321 (2021). It consists of classical importance sampling of initial states and a measurement of the Loschmidt echo with a quantum simulator. We use the method as a quantum-inspired classical algorithm and simulate the protocol with matrix product states to analyze the requirements on a quantum simulator. This way, we show that a finite temperature phase transition in the long-range transverse field Ising model can be characterized in trapped ion quantum simulators. We propose a concrete measurement protocol for the Loschmidt echo and discuss the influence of measurement noise, dephasing, as well as state preparation and measurement errors. We argue that the algorithm is robust against those imperfections under realistic conditions.

The eigenstate thermalization hypothesis (ETH) offers a universal mechanism for the approach to equilibrium of closed quantum many-body systems. So far, however, experimental studies have focused on the relaxation dynamics of observables as described by the diagonal part of ETH, whose verification requires substantial numerical input. This leaves many of the general assumptions of ETH untested. Here, we propose a theory-independent route to probe the full ETH in quantum simulators by observing the emergence of fluctuation-dissipation relations, which directly probe the off-diagonal part of ETH. We discuss and propose protocols to independently measure fluctuations and dissipations as well as higher-order time ordered correlation functions. We first show how the emergence of fluctuation dissipation relations from a nonequilibrium initial state can be observed for the 2D Bose-Hubbard model in superconducting qubits or quantum gas microscopes. Then we focus on the long-range transverse field Ising model (LTFI), which can be realized with trapped ions. The LTFI exhibits rich thermalization phenomena: For strong transverse fields, we observe prethermalization to an effective magnetization-conserving Hamiltonian in the fluctuation dissipation relations. For weak transverse fields, confined excitations lead to non-thermal features resulting in a violation of the fluctuation-dissipation relations up to long times. Moreover, in an integrable region of the LTFI, thermalization to a generalized Gibbs ensemble occurs and the fluctuation-dissipation relations enable an experimental diagonalization of the Hamiltonian. Our work presents a theory-independent way to characterize thermalization in quantum simulators and paves the way to quantum simulate condensed matter pump-probe experiments.

Generic short-range interacting quantum systems with a conserved quantity exhibit universal diffusive transport at late times. We employ non-equilibrium quantum field theory and semi-classical phase-space simulations to show how this universality is replaced by a more general transport process in a long-range XY spin chain at infinite temperature with couplings decaying algebraically with distance as $r^{-α}$. While diffusion is recovered for $α>1.5$, longer-ranged couplings with $0.5<α\leq 1.5 $ give rise to effective classical Lévy flights; a random walk with step sizes drawn from a distribution with algebraic tails. We find that the space-time dependent spin density profiles are self-similar, with scaling functions given by the stable symmetric distributions. As a consequence, for $0.5<α\leq1.5$ autocorrelations show hydrodynamic tails decaying in time as $t^{-1/(2α-1)}$ and linear-response theory breaks down. Our findings can be readily verified with current trapped ion experiments.