Michael L. Wall

We demonstrate the use of matrix product state (MPS) models for discriminating quantum data on quantum computers using holographic algorithms, focusing on classifying a translationally invariant quantum state based on $L$ qubits of quantum data extracted from it. We detail a process in which data from single-shot experimental measurements are used to optimize an isometric tensor network, the tensors are compiled into unitary quantum operations using greedy compilation heuristics, parameter optimization on the resulting quantum circuit model removes the post-selection requirements of the isometric tensor model, and the resulting quantum model is inferenced on either product state or entangled quantum data. We demonstrate our training and inference architecture on a synthetic dataset of six-site single-shot measurements from the bulk of a one-dimensional transverse field Ising model (TFIM) deep in its antiferromagnetic and paramagnetic phases. We experimentally evaluate models on Quantinuum's H1-2 trapped ion quantum computer using entangled input data modeled as translationally invariant, bond dimension 4 MPSs across the known quantum phase transition of the TFIM. Using linear regression on the experimental data near the transition point, we find predictions for the critical transverse field of $h=0.962$ and $0.994$ for tensor network discriminators of bond dimension $χ=2$ and $χ=4$, respectively. These predictions compare favorably with the known transition location of $h=1$ despite training on data far from the transition point. Our techniques identify families of short-depth variational quantum circuits in a data-driven and hardware-aware fashion and robust classical techniques to precondition the model parameters, and can be adapted beyond machine learning to myriad applications of tensor networks on quantum computers, such as quantum simulation and error correction.

The Fermi-Hubbard model describes ultracold fermions in an optical lattice and exhibits antiferromagnetic long-ranged order below the Néel temperature. However, reaching this temperature in the lab has remained an elusive goal. In other atomic systems, such as trapped ions, low temperatures have been successfully obtained by adiabatic demagnetization, in which a strong effective magnetic field is applied to a spin-polarized system, and the magnetic field is adiabatically reduced to zero. Unfortunately, applying this approach to the Fermi-Hubbard model encounters a fundamental obstacle: the $SU(2)$ symmetry introduces many level crossings that prevent the system from reaching the ground state, even in principle. However, by breaking the $SU(2)$ symmetry with a spin-dependent tunneling, we show that adiabatic demagnetization can achieve low temperature states. Using density matrix renormalization group (DMRG) calculations in one dimension, we numerically find that demagnetization protocols successfully reach low temperature states of a spin-anisotropic Hubbard model, and we discuss how to optimize this protocol for experimental viability. By subsequently ramping spin-dependent tunnelings to spin-independent tunnelings, we expect that our protocol can be employed to produce low-temperature states of the Fermi-Hubbard Model.

We investigate an extension of the quantum Ising model in one spatial dimension including long-range $1 / r^α$ interactions in its statics and dynamics with possible applications from heteronuclear polar molecules in optical lattices to trapped ions described by two-state spin systems. We introduce the statics of the system via both numerical techniques with finite size and infinite size matrix product states and a theoretical approaches using a truncated Jordan-Wigner transformation for the ferromagnetic and antiferromagnetic case and show that finite size effects have a crucial role shifting the quantum critical point of the external field by fifteen percent between thirty-two and around five-hundred spins. We numerically study the Kibble-Zurek hypothesis in the long-range quantum Ising model with Matrix Product States. A linear quench of the external field through the quantum critical point yields a power-law scaling of the defect density as a function of the total quench time. For example, the increase of the defect density is slower for longer-range models and the critical exponent changes by twenty-five per cent. Our study emphasizes the importance of such long-range interactions in statics and dynamics that could point to similar phenomena in a different setup of dynamical systems or for other models.

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.

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.