Matter under extreme conditions behaves as a non-equilibrium, strongly-interacting quantum system. Such systems are in the non-perturbative regime of quantum chromodynamics (QCD), where their time evolution is classically intractable. Quantum computing offers a promising alternative, but present-day devices are limited in qubit count and circuit depth, so efficient formulations and algorithms are key to scaling toward experimentally relevant QCD simulations. A scalable framework is established for simulating energy loss and hadronization in 1+1-dimensional SU(2) lattice gauge theory, in which a heavy quark traverses a lattice of light quarks. Classical lattice QCD methods are translated into physics-informed circuits to prepare the interacting ground state. The evolution is run on IBM's quantum processor ibm_pittsburgh, where the charge distribution and energy loss are measured. A further study develops the quantum simulation of the fully gauge-fixed SU(2) Hamiltonian in a mixed basis, a formulation that is expected to remain efficient at all gauge couplings. Digitizing the continuous gauge fields with as few as three qubits per plaquette reproduces low-energy observables with
10−3 precision. Scaling to larger systems is limited by resource overhead, which motivates truncations that preserve the relevant observables. Sequency Hierarchy Truncation (SeqHT) exploits the hierarchy in basis-operator contributions. Demonstrated on the adiabatic preparation of the
λφ4 ground state on ibm_sherbrooke, SeqHT reduces circuit depth by ~30%, suppressing hardware noise in the extracted observables. All simulations in this dissertation are validated on quantum hardware against classical benchmarks. These studies advance quantum simulations of real-time QCD dynamics.
Hamiltonian quantum simulation of bosons on digital quantum computers requires truncating the Hilbert space to finite dimensions. The method of truncation and the choice of basis states can significantly impact the complexity of the quantum circuit required to simulate the system. For example, a truncation in the Fock basis where each boson is encoded with a register of
Q qubits, can result in an exponentially large number of Pauli strings required to decompose the truncated Hamiltonian. This, in turn, can lead to an exponential increase in
Q in the complexity of the quantum circuit. For lattice quantum field theories such as Yang-Mills theory and QCD, several Hamiltonian formulations and corresponding truncations have been put forward in recent years. There is no exponential increase in
Q when resorting to the orbifold lattice Hamiltonian, while we do not know how to remove the exponential complexity in
Q in the commonly used Kogut-Susskind Hamiltonian. Specifically, when using the orbifold lattice Hamiltonian, the continuum limit, or, in other words, the removal of the ultraviolet energy cutoff, is obtained with circuits whose resources scale like
Q, while they scale like
O(exp(Q)) for the Kogut-Susskind Hamiltonian: this can be seen as an exponential speed up in approaching the physical continuum limit for the orbifold lattice Hamiltonian formulation. We show that the universal framework, advocated by three of the authors (M.~H., S.~M., and E.~R.) and collaborators, provides a natural avenue to solve the exponential scaling of circuit complexity with
Q, and it is the reason why using the orbifold lattice Hamiltonian is advantageous.
Qubit, operator and gate resources required for the digitization of lattice
λϕ4 scalar field theories onto quantum computers in the NISQ era are considered, building upon the foundational work by Jordan, Lee and Preskill. The Nyquist-Shannon sampling theorem, introduced in this context by Macridin, Spetzouris, Amundson and Harnik building on the work of Somma, provides a guide with which to evaluate the efficacy of two field-space bases, the eigenstates of the field operator, as used by Jordan, Lee and Preskill, and eigenstates of a harmonic oscillator, to describe
0+1- and
1+1-dimensional scalar field theory. We show how techniques associated with improved actions, which are heavily utilized in Lattice QCD calculations to systematically reduce lattice-spacing artifacts, can be used to reduce the impact of the field digitization in
λϕ4, but are found to be inferior to a complete digitization-improvement of the Hamiltonian using a Quantum Fourier Transform. When the Nyquist-Shannon sampling theorem is satisfied, digitization errors scale as
∣log∣log∣ϵdig∣∣∣∼nQ (number of qubits describing the field at a given spatial site) for the low-lying states, leaving the familiar power-law lattice-spacing and finite-volume effects that scale as
∣log∣ϵlatt∣∣∼NQ (total number of qubits in the simulation). We find that fewer than
nQ=10 qubits per spatial lattice site are sufficient to reduce digitization errors below noise levels expected in NISQ-era quantum devices for both localized and delocalized field-space wavefunctions. For localized wavefunctions,
nQ=4 qubits are likely to be sufficient for calculations requiring modest precision.
Qubit, operator, and gate resources required for the digitization of lattice $\ensuremath{\lambda}{\ensuremath{\phi}}^{4}$ scalar field theories onto quantum computers are considered, building upon the foundational work by Jordan et al. [Quantum Inf. Comput. 14, 1014 (2014); Science 336, 1130 (2012)], with a focus towards noisy intermediate-scale quantum devices. The Nyquist-Shannon sampling theorem, introduced in this context by Macridin et al. [Phys. Rev. A 98, 042312 (2018)] building on the work of Somma [Quantum Inf. Comput. 16, 1125 (2016)], provides a guide with which to evaluate the efficacy of two field-space bases, the eigenstates of the field operator, as used by Jordan et al., and eigenstates of a harmonic oscillator, to describe (
0+1)- and (
d+1)-dimensional scalar field theory. We show how techniques associated with improved actions, which are heavily utilized in lattice QCD calculations to systematically reduce lattice-spacing artifacts, can be used to reduce the impact of the field digitization in $\ensuremath{\lambda}{\ensuremath{\phi}}^{4}$, but are found to be inferior to a complete digitization improvement of the Hamiltonian using a quantum Fourier transform. When the Nyquist-Shannon sampling theorem is satisfied, digitization errors scale as $|log|log|{\ensuremath{\epsilon}}_{\mathrm{dig}}|||\ensuremath{\sim}{n}_{Q}$ (number of qubits describing the field at a given spatial site) for the low-lying states, leaving the familiar power-law lattice-spacing and finite-volume effects that scale as $|log|{\ensuremath{\epsilon}}_{\mathrm{latt}}||\ensuremath{\sim}{N}_{Q}$ (total number of qubits in the simulation). For localized (delocalized) field-space wave functions, it is found that ${n}_{Q}\ensuremath{\sim}4(7)$ qubits per spatial lattice site are sufficient to reduce theoretical digitization errors below error contributions associated with approximation of the time-evolution operator and noisy implementation on near-term quantum devices. Only classical computing resources have been used to obtain the results presented in this work.
Digital quantum simulation of many-body dynamics faces a tension between algorithmic Trotter error and physical noise that accumulates with circuit depth. Typical higher-order product formulas mitigate the algorithmic error at the cost of deeper circuits. Our benchmark shows that, within the limited physical error budget, the feasible advantage is confined to absolute errors well below the
10−2 scale, which vanishes under noise levels of current quantum hardware. In sharp contrast to the previous Trotter error mitigation methods, we introduce a constant-depth multi-product formula (cd-MPF) that suppresses the Trotter error by combining multiple circuits at fixed circuit depth. We identify an auxiliary parameter
α, which reshapes the Trotter error terms while leaving the target evolution invariant. The classical linear combination of the measured expectation values cancels the leading
(Δt)2 algorithmic contribution and steepens the Trotter-error scaling with the two-qubit circuit depth
d from
d−2 to
d−4. Combined with physical-noise mitigation, our method can serve as a key ingredient for realizing long-time quantum dynamics simulation on near-term hardware.
We report a first-of-its-kind analysis on post-Trotter simulation of U(1), SU(2), and SU(3) lattice gauge theories including fermions in arbitrary spatial dimension. We provide explicit circuit constructions as well as T -gate counts and logical qubit counts for Hamiltonian simulation. We find a reduction of up to 25 orders of magnitude in space-time volume over Trotter methods for simulations of non-Abelian lattice gauge theories relevant to the standard model. This improvement results from our algorithm having polynomial scaling with the number of colors in the gauge theory, achieved by utilizing oracle constructions relying on the sparsity of physical operators, in contrast to the exponential scaling seen in state-of-the-art Trotter methods, which employ explicit mappings onto Pauli operators. Our work demonstrates that the use of advanced algorithmic techniques leads to dramatic reductions in the cost of simulating fundamental interactions, bringing it in step with resources required for first-principles quantum simulation of chemistry. Published by the American Physical Society 2024