The question of how to efficiently formulate Hamiltonian gauge theories is experiencing renewed interest due to advances in building quantum simulation platforms. We introduce a reformulation of an SU(2) Hamiltonian lattice gauge theory—a loop-string-hadron (LSH) formulation—that describes dynamics directly in terms of its loop, string, and hadron degrees of freedom, while alleviating several disadvantages of quantum-simulating the Kogut-Susskind formulation. This LSH formulation transcends the local loop formulation of d+1-dimensional lattice gauge theories by incorporating staggered quarks, furnishing the algebra of gauge-singlet operators, and being used to reconstruct dynamics between states that have Gauss’s law built in to them. LSH operators are then factored into products of “normalized” ladder operators and diagonal matrices, priming them for classical or quantum information processing. Self-contained expressions of the Hamiltonian are given up to d=3. The LSH formalism uses makes little use of structures specific to SU(2) and its conceptual clarity makes it an attractive approach to apply to other non-Abelian groups like SU(3).

Correlations and measures of entanglement in ground state wavefunctions of relativistic quantum field theories are spatially localized over length scales set by the mass of the lightest particle. We utilize this localization to design digital quantum circuits for preparing the ground states of free lattice scalar quantum field theories. Controlled rotations that are exponentially localized in their position-space extent are found to provide exponentially convergent wavefunction fidelity. These angles scale with the classical two-point correlation function, as opposed to the more localized mutual information or the hyper-localized negativity. We anticipate that further investigations will uncover quantum circuit designs with controlled rotations dictated by the measures of entanglement. This work is expected to impact quantum simulations of systems of importance to nuclear physics, high-energy physics, and basic energy sciences research.

Renormalization group ideas and effective operators are used to efficiently determine localized unitaries for preparing the ground states of non-interacting scalar field theories on digital quantum devices. With these methods, classically computed ground states in a small spatial volume can be used to determine operators for preparing the ground state in a beyond-classical quantum register, even for interacting scalar field theories. Due to the exponential decay of correlation functions and the double exponential suppression of digitization artifacts, the derived quantum circuits are expected to be relevant already for near-term quantum devices.

The great promise of quantum computers comes with the dual challenges of building them and finding their useful applications. We argue that these two challenges should be considered together, by co-designing full-stack quantum computer systems along with their applications in order to hasten their development and potential for scientific discovery. In this context, we identify scientific and community needs, opportunities, a sampling of a few use case studies, and significant challenges for the development of quantum computers for science over the next 2–10 years. This document is written by a community of university, national laboratory, and industrial researchers in the field of Quantum Information Science and Technology, and is based on a summary from a U.S. National Science Foundation workshop on Quantum Computing held on October 21–22, 2019 in Alexandria, VA.

Neutrino-nucleus cross section uncertainties are expected to be a dominant systematic in future accelerator neutrino experiments. The cross sections are determined by the linear response of the nucleus to the weak interactions of the neutrino, and are dominated by energy and distance scales of the order of the separation between nucleons in the nucleus. These response functions are potentially an important early physics application of quantum computers. Here we present an analysis of the resources required and their expected scaling for scattering cross section calculations. We also examine simple small-scale neutrino-nucleus models on modern quantum hardware. In this paper, we use variational methods to obtain the ground state of a three nucleon system (the triton) and then implement the relevant time evolution. In order to tame the errors in present-day NISQ devices, we explore the use of different error-mitigation techniques to increase the fidelity of the calculations.

An improved mapping of one-dimensional SU(2) non-Abelian gauge theory onto qubit degrees of freedom is presented. This new mapping allows for a reduced unphysical Hilbert space. Insensitivity to interactions within this unphysical space is exploited to design more efficient quantum circuits. Local gauge symmetry is used to analytically incorporate the angular momentum alignment, leading to qubit registers encoding the total angular momentum on each link. Results of multi-plaquette calculations on IBM’s quantum hardware are presented.

The evaluation of expectation values Tr[ρO] for some pure state ρ and Hermitian operator O is of central importance in a variety of quantum algorithms. Near optimal techniques developed in the past require a number of measurements N approaching the Heisenberg limit N=O(1/ϵ) as a function of target accuracy ϵ. The use of Quantum Phase Estimation requires however long circuit depths C=O(1/ϵ) making their implementation difficult on near term noisy devices. The more direct strategy of Operator Averaging is usually preferred as it can be performed using N=O(1/ϵ²) measurements and no additional gates besides those needed for the state preparation. In this work we use a simple but realistic model to describe the bound state of a neutron and a proton (the deuteron) and show that the latter strategy can require an overly large number of measurement in order to achieve a reasonably small relative target accuracy ϵr. We propose to overcome this problem using a single step of QPE and classical post-processing. This approach leads to a circuit depth C=O(ϵ^m) (with m≥0) and to a number of measurements N=O(1/ϵ²⁺ⁿ) for 0<n≤1. We provide detailed descriptions of two implementations of our strategy for n=1 and n≈0.5 and derive appropriate conditions that a particular problem instance has to satisfy in order for our method to provide an advantage.

Initializing a single site of a lattice scalar field theory into an arbitrary state requires O(2^nQ) entangling gates on a quantum computer with n_Q qubits per site. It is conceivable that, instead, initializing to functions that are good approximations to states may have utility in reducing the number of required entangling gates. In the case of a single site of a non-interacting scalar field theory, initializing to a symmetric exponential wavefunction requires n_Q − 1 entangling gates, compared with the 2^nQ−1 + n_Q − 3 + δ_nQ,1 required for a symmetric Gaussian wavefunction. In this work, we explore the initialization of 1-site (n_Q = 4), 2-site (n_Q = 3) and 3-site (n_Q = 3) non- interacting scalar field theories with symmetric exponential wavefunctions using IBM’s quantum simulators and quantum devices (Poughkeepsie and Tokyo). With the digitizations obtainable with n_Q = 3, 4, these tensor-product wavefunctions are found to have large overlap with a Gaussian wavefunction, and provide a suitable low-noise initialization for subsequent quantum simulations. In performing these simulations, we have employed a workflow that interleaves calibrations to mitigate systematic errors in production. The calibrations allow tolerance cuts on gate performance including the fidelity of the symmetrizing Hadamard gate, both in vacuum (|0⟩^⊗nQ ) and in medium (n_Q − 1 qubits initialized to an exponential function). The results obtained in this work are relevant to systems beyond scalar field theories, such as the deuteron radial wavefunction, 2- and 3-dimensional cartesian-space wavefunctions, and non-relativistic multi-nucleon systems built on a localized eigenbasis.

We show that using the loop-string-hadron (LSH) formulation of SU(2) lattice gauge theory (arXiv:1912.06133) as a basis for digital quantum computation easily solves an important problem of fundamental interest: implementing gauge invariance (or Gauss’s law) exactly. We first discuss the structure of the LSH Hilbert space in d spatial dimensions, its truncation, and its digitization with qubits. Error detection and mitigation in gauge theory simulations would benefit from physicality “oracles,” so we decompose circuits that flag gauge invariant wavefunctions. We then analyze the logical qubit costs and entangling gate counts involved with the protocols. The LSH basis could save or cost more qubits than a Kogut-Susskind-type representation basis, depending on how that is digitized as well as the spatial dimension. The numerous other clear benefits encourage future studies into applying this framework.

Entanglement suppression in the strong interaction S-matrix is shown to be correlated with approximate spin-flavor symmetries that are observed in low-energy baryon interactions, the Wigner SU(4) symmetry for two flavors and an SU(16) symmetry for three flavors. We conjecture that dynamical entanglement suppression is a property of the strong interactions in the infrared, giving rise to these emergent symmetries and providing powerful constraints on the nature of nuclear and hypernuclear forces in dense matter.