Simulating the real-time dynamics of quantum field theories (QFTs) is one of the most promising applications of quantum simulators. Regularizing a bosonic QFT for quantum simulation purposes typically involves a truncation in Hilbert space in addition to a discretization of space. Here, we discuss how to perform such a regularization of scalar QFTs by explicitly constructing suitable many-body lattice Hamiltonians using multilevel or qudit systems and show that this enables quantitative predictions in the continuum limit by extrapolating results obtained for large-spin models. With extensive matrix-product-state simulations, we numerically demonstrate the sequence of extrapolations that leads to quantitative agreement of observables for the integrable sine-Gordon (sG) QFT. We further show how to prepare static and moving-soliton excitations and we analyze their scattering dynamics in the continuum limit, in agreement with a semiclassical model and with quantitative analytical predictions. Finally, we illustrate how a nonintegrable perturbation of the sG model gives rise to dynamics reminiscent of string breaking and plasma oscillations in gauge theories. Our methods are directly applicable in state-of-the-art analog quantum simulators, opening the door to quantitatively investigating a wide variety of scalar-field theories and tackling long-standing questions in nonequilibrium QFT such as the fate of the false vacuum.
Simulating the real-time dynamics of quantum field theories (QFTs) is one of the most promising applications of quantum simulators. Regularizing a bosonic QFT for quantum simulation purposes typically involves a truncation in Hilbert space in addition to a discretization of space. Here, we discuss how to perform such a regularization of scalar QFTs by explicitly constructing suitable many-body lattice Hamiltonians using multilevel or qudit systems and show that this enables quantitative predictions in the continuum limit by extrapolating results obtained for large-spin models. With extensive matrix-product-state simulations, we numerically demonstrate the sequence of extrapolations that leads to quantitative agreement of observables for the integrable sine-Gordon (sG) QFT. We further show how to prepare static and moving-soliton excitations and we analyze their scattering dynamics in the continuum limit, in agreement with a semiclassical model and with quantitative analytical predictions. Finally, we illustrate how a nonintegrable perturbation of the sG model gives rise to dynamics reminiscent of string breaking and plasma oscillations in gauge theories. Our methods are directly applicable in state-of-the-art analog quantum simulators, opening the door to quantitatively investigating a wide variety of scalar-field theories and tackling long-standing questions in nonequilibrium QFT such as the fate of the false vacuum.
The solution of gauge theories is one of the most promising applications of quantum technologies. Here, we discuss the approach to the continuum limit for U(1) gauge theories regularized via finite-dimensional Hilbert spaces of quantum spin-S operators, known as quantum link models. For quantum electrodynamics (QED) in one spatial dimension, we numerically demonstrate the continuum limit by extrapolating the ground state energy, the scalar, and the vector meson masses to large spin lengths S, large volume N, and vanishing lattice spacing a. By analytically solving Gauss' law for arbitrary S, we obtain a generalized PXP spin model and count the physical Hilbert space dimension analytically. This allows us to quantify the required resources for reliable extrapolations to the continuum limit on quantum devices. We use a functional integral approach to relate the model with large values of half-integer spins to the physics at topological angle Θ=π. Our findings indicate that quantum devices will in the foreseeable future be able to quantitatively probe the QED regime with quantum link models.
Quantum simulation of lattice gauge theories offers a path to first-principles calculations of real-time dynamics and finite-density phenomena that are inaccessible to classical Euclidean Monte Carlo methods. Realizing this potential requires a pipeline of controlled approximations---lattice discretization, truncation, digitization, state preparation, time evolution, measurement, and continuum extrapolation---each contributing errors and computational costs. General algorithmic asymptotic complexity arguments capture only part of the cost picture: the practical resources required for end-to-end simulation depend on numerical pre-factors, the cost of the subroutines that various algorithms call, and structural features of the lattice theory, all of which the asymptotic scaling leaves unspecified. This thesis develops contributions that aim to make such dependencies explicit and quantitative. We begin with a rigorous analysis of the asymptotic gate complexity of product formulas and Quantum Signal Processing (QSP) for the time evolution of general classes of lattice Hamiltonians, identifying the parameter regimes in which each method is favorable and estimating the associated pre-factors numerically. Since the cost of QSP-based simulation is dominated (in regimes of practical interest) by the block-encoding subroutine it calls, we also develop several novel block-encoding constructions tailored to bosonic lattice field theories and show that signal-processing techniques can be repurposed as highly efficient block-encoding subroutines. We then turn to the problem of truncating the Kogut--Susskind electric-field Hilbert space. Existing error bounds rely on energy-conservation arguments and require the truncation to grow with simulation time. By exploiting Hilbert space fragmentation we derive the tightest error bounds known to date: they decay factorially in the truncation parameter and are independent of simulation time. We validate them numerically for U(1) pure gauge theory and the Schwinger model. Finally, we address the continuum limit in the presence of approximate time evolution. For product formulas, we show that Trotter errors correspond to irrelevant operators that vanish in the continuum limit, which yields a renormalization trajectory independent of the Trotter step size. More generally, we introduce the Statistically-Bounded Time Evolution (SBTE) protocol, an algorithm-independent framework that treats approximation errors as a systematic uncertainty to be driven below the working statistical precision. This decouples the continuum-limit prescription from the choice of time-evolution algorithm and provides a basis for end-to-end cost comparisons across simulation methods.
We use bosonic tensor networks to connect the lattice sine-Gordon model in the Hamiltonian formulation quantitatively to its continuum theory. Matching the lattice vertex operator to its conformal normalization at the free-boson ultraviolet fixed point yields the exact relation between the bare lattice coupling and the renormalized continuum mass parameter. The soliton mass then approaches Zamolodchikov's exact continuum prediction throughout the studied parameter range, without adjustable parameters. Using uniform matrix product states and a quasiparticle ansatz, we also recover the relativistic soliton dispersion and the two lightest breather masses at the percent level. We simulate real-time collisions of Gaussian soliton-antisoliton wave packets near a reflectionless point and extract the Wigner spatial displacement. We compare this displacement with the exact continuum prediction obtained from the momentum derivative of the transmission phase, recovering its characteristic rapidity dependence. Our bosonic simulations provide a foundation for nonintegrable extensions, offer lessons for renormalization in other Hamiltonian lattice field theories, including gauge theories, and provide benchmarks for continuous-variable quantum simulations.
Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use 2+1 dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.
Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use 2 + 1 dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.
In the last decade the non-equilibrium dynamics of isolated quantum systems has raised a huge interest for both its conceptual and practical implications. Simulation of quantum dynamics represents by now a very active research effort, with highly controllable quantum devices composed of hundreds of qubits already being realized in numerous setups. These offer the unprecedented possibility to access real-time quantum dynamics and strongly correlated quantum many-body states, opening numerous avenues not only for future quantum technologies, but also for long-standing theoretical problems. One of the most sought-after applications is the computation of real-time evolution in theories of fundamental interactions. Since Wilson's proposal in 1974, the equilibrium phase diagrams of lattice gauge theories have been successfully investigated with Monte Carlo simulations in a wide range of parameters. Nonequilibrium properties, on the other hand, cannot be accessed with conventional techniques. Thanks to the recent advances in numerics (e.g. tensor networks) and experiments (quantum simulation), new promising tools are now available for tackling these questions. Remarkable works have already demonstrated that applying these techniques to simple theories, such as one-dimensional quantum electrodynamics, can help to detect the signatures of confinement and pair creation. Despite these impressive achievements, investigating the non-equilibrium dynamics of gauge theories with quantum simulators turned out to be extremely challenging. This is mostly due to the complex constrained dynamics of these theories, which make their simulation with realistic atomic platforms particularly difficult. At the same time, the dynamical constraints can play a key role in determining novel interesting phenomena and peculiar non-equilibrium properties. These phenomena -- which include, for example, quantum scars, fractons, Hilbert space fragmentation -- can lead to anomalous dynamics, eventually evading the paradigm of thermalization, which predicts the relaxation of physical observable to their thermodynamic values. The aim of this thesis is twofold: first, we want to address through a combination of numerical and analytical tools some of the open questions regarding the non-equilibrium dynamics of quantum many-body systems with constrained Hilbert spaces, including lattice gauge theories. In addition, the work reported here is part of a large effort of the community towards the quantum simulation of complex phenomena of condensed matter and high-energy physics.