Showing 1–50 of 295 results for author: Savage, M

The Magic in Nuclear and Hypernuclear Forces

Authors: Caroline E. P. Robin, Martin J. Savage

Abstract: Toward an improved understanding of the role of quantum information in nuclei and exotic matter, we examine the magic (non-stabilizerness) in low-energy strong interaction processes. As stabilizer states can be prepared efficiently using classical computers, and include classes of entangled states, it is magic and fluctuations in magic, along with entanglement, that determine resource requirements… ▽ More Toward an improved understanding of the role of quantum information in nuclei and exotic matter, we examine the magic (non-stabilizerness) in low-energy strong interaction processes. As stabilizer states can be prepared efficiently using classical computers, and include classes of entangled states, it is magic and fluctuations in magic, along with entanglement, that determine resource requirements for quantum simulations. As a measure of fluctuations in magic induced by scattering, the "magic power" of the S-matrix is introduced. Using experimentally-determined scattering phase shifts and mixing parameters, the magic power in nucleon-nucleon and hyperon-nucleon scattering, along with the magic in the deuteron, are found to exhibit interesting features. The $Σ^-$-baryon is identified as a potential candidate catalyst for enhanced spreading of magic and entanglement in dense matter, depending on in-medium decoherence. △ Less

Submitted 20 May, 2024; v1 submitted 16 May, 2024; originally announced May 2024.

Comments: 14 pages, 6 figures. v2: additional references, modifications to uncertainty bands

Report number: IQuS@UW-21-077

Steps Toward Quantum Simulations of Hadronization and Energy-Loss in Dense Matter

Authors: Roland C. Farrell, Marc Illa, Martin J. Savage

Abstract: A framework for simulating the real-time dynamics of particles in dense matter using quantum computers is developed. As a demonstration, we perform classical simulations of heavy-hadrons propagating through a dense medium in the Schwinger model. Measurements of the time-dependent energy and charge density are used to identify mechanisms responsible for energy loss and hadron production (hadronizat… ▽ More A framework for simulating the real-time dynamics of particles in dense matter using quantum computers is developed. As a demonstration, we perform classical simulations of heavy-hadrons propagating through a dense medium in the Schwinger model. Measurements of the time-dependent energy and charge density are used to identify mechanisms responsible for energy loss and hadron production (hadronization). A study of entanglement dynamics highlights the importance of quantum coherence between the particles that make up the dense medium. Throughout this work, care is taken to isolate, and remove, phenomena that arise solely from a finite lattice spacing. It is found that signatures of entanglement are more sensitive to lattice artifacts than other observables. Toward quantum simulations, we present an efficient method and the corresponding quantum circuits for preparing ground states in the presence of heavy mesons. These circuits are used to estimate the resources required to simulate in-medium energy loss and hadronization in the Schwinger model using quantum computers. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 22 pages, 10 figures, 3 tables

Report number: IQuS@UW-21-076, NT@UW-24-06

Quantum Simulating Nature's Fundamental Fields

Authors: Christian W. Bauer, Zohreh Davoudi, Natalie Klco, Martin J. Savage

Abstract: Simulating key static and dynamic properties of matter -- from creation in the Big Bang to evolution into sub-atomic and astrophysical environments -- arising from the underlying fundamental quantum fields of the Standard Model and their effective descriptions, lies beyond the capabilities of classical computation alone. Advances in quantum technologies have improved control over quantum entanglem… ▽ More Simulating key static and dynamic properties of matter -- from creation in the Big Bang to evolution into sub-atomic and astrophysical environments -- arising from the underlying fundamental quantum fields of the Standard Model and their effective descriptions, lies beyond the capabilities of classical computation alone. Advances in quantum technologies have improved control over quantum entanglement and coherence to the point where robust simulations are anticipated to be possible in the foreseeable future. We discuss the emerging area of quantum simulations of Standard-Model physics, challenges that lie ahead, and opportunities for progress in the context of nuclear and high-energy physics. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 22 pages, 4 figures, Nature Reviews Physics, Perspective article, version from February 2023

Report number: IQuS@UW-21-052

Journal ref: Nature Rev.Phys. 5 (2023) 7, 420-432

Qu8its for Quantum Simulations of Lattice Quantum Chromodynamics

Authors: Marc Illa, Caroline E. P. Robin, Martin J. Savage

Abstract: We explore the utility of $d=8$ qudits, qu8its, for quantum simulations of the dynamics of 1+1D SU(3) lattice quantum chromodynamics, including a mapping for arbitrary numbers of flavors and lattice size and a re-organization of the Hamiltonian for efficient time-evolution. Recent advances in parallel gate applications, along with the shorter application times of single-qudit operations compared w… ▽ More We explore the utility of $d=8$ qudits, qu8its, for quantum simulations of the dynamics of 1+1D SU(3) lattice quantum chromodynamics, including a mapping for arbitrary numbers of flavors and lattice size and a re-organization of the Hamiltonian for efficient time-evolution. Recent advances in parallel gate applications, along with the shorter application times of single-qudit operations compared with two-qudit operations, lead to significant projected advantages in quantum simulation fidelities and circuit depths using qu8its rather than qubits. The number of two-qudit entangling gates required for time evolution using qu8its is found to be more than a factor of five fewer than for qubits. We anticipate that the developments presented in this work will enable improved quantum simulations to be performed using emerging quantum hardware. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 20 pages, 4 figures, 1 table

Report number: IQuS@UW-21-074

Quantum Simulations of Hadron Dynamics in the Schwinger Model using 112 Qubits

Authors: Roland C. Farrell, Marc Illa, Anthony N. Ciavarella, Martin J. Savage

Abstract: Hadron wavepackets are prepared and time evolved in the Schwinger model using 112 qubits of IBM's 133-qubit Heron quantum computer ibm_torino. The initialization of the hadron wavepacket is performed in two steps. First, the vacuum is prepared across the whole lattice using the recently developed SC-ADAPT-VQE algorithm and workflow. SC-ADAPT-VQE is then extended to the preparation of localized sta… ▽ More Hadron wavepackets are prepared and time evolved in the Schwinger model using 112 qubits of IBM's 133-qubit Heron quantum computer ibm_torino. The initialization of the hadron wavepacket is performed in two steps. First, the vacuum is prepared across the whole lattice using the recently developed SC-ADAPT-VQE algorithm and workflow. SC-ADAPT-VQE is then extended to the preparation of localized states, and used to establish a hadron wavepacket on top of the vacuum. This is done by adaptively constructing low-depth circuits that maximize the overlap with an adiabatically prepared hadron wavepacket. Due to the localized nature of the wavepacket, these circuits can be determined on a sequence of small lattices using classical computers, and then robustly scaled to prepare wavepackets on large lattices for simulations using quantum computers. Time evolution is implemented with a second-order Trotterization. To reduce both the required qubit connectivity and circuit depth, an approximate quasi-local interaction is introduced. This approximation is made possible by the emergence of confinement at long distances, and converges exponentially with increasing distance of the interactions. Using multiple error-mitigation strategies, up to 14 Trotter steps of time evolution are performed, employing 13,858 two-qubit gates (with a CNOT depth of 370). The propagation of hadrons is clearly identified, with results that compare favorably with Matrix Product State simulations. Prospects for a near-term quantum advantage in simulations of hadron scattering are discussed. △ Less

Submitted 11 June, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: 54 pages, 19 figures, 20 tables

Report number: IQuS@UW-21-069, NT@UW-24-1

Journal ref: Phys.Rev.D 109 (2024) 11, 114510

Time Scales in Many-Body Fast Neutrino Flavor Conversion

Authors: Ramya Bhaskar, Alessandro Roggero, Martin J. Savage

Abstract: Time scales associated with many-body fast neutrino flavor conversions in core-collapse supernova are explored in the context of an effective two-flavor model with axial symmetry. We present a preliminary study of time scales obtained from a linear stability analysis and from the distributions of Loschmidt echo crossing times (intimately connected to dynamical phase transitions in non-equilibrium… ▽ More Time scales associated with many-body fast neutrino flavor conversions in core-collapse supernova are explored in the context of an effective two-flavor model with axial symmetry. We present a preliminary study of time scales obtained from a linear stability analysis and from the distributions of Loschmidt echo crossing times (intimately connected to dynamical phase transitions in non-equilibrium systems) determined by time evolution with the exact many-body Hamiltonian. Starting from a tensor-product initial state describing systems of $N$ neutrinos, with $N/2$ electron-type and $N/2$ heavy-type, with uniform angular distributions, the Loschmidt echo crossing times, $t_{\mathcal{L}_{\times}}$, are found to exhibit two distinct time scales that are exponentially separated. The second peak structure at longer times, effectively absent for $N=4$, develops with increasing $N$. When re-scaled in terms of $\log t_{\mathcal{L}_{\times}}$, the distributions are found to become increasingly well described by the sum of two stable distributions. The distribution of Loschmidt echo crossing times differs somewhat from the results of the (numerical) linear stability analysis, which exhibits a peak at finite frequency and a second peak consistent with zero frequency. The exact analysis suggests that the zero-frequency instability manifests itself as a modest flavor-conversion time scale. △ Less

Submitted 4 March, 2024; v1 submitted 23 December, 2023; originally announced December 2023.

Comments: 13 pages, 5 figures; Updated references, and acknowledgments section

Report number: IQuS@UW-21-069

Quantum Computing for Nuclear Physics

Authors: Martin J. Savage

Abstract: Future quantum computers are anticipated to be able to perform simulations of quantum many-body systems and quantum field theories that lie beyond the capabilities of classical computation. This will lead to new insights and predictions for systems ranging from dense non-equilibrium matter, to low-energy nuclear structure and reactions, to high-energy collisions. I present an overview of digital q… ▽ More Future quantum computers are anticipated to be able to perform simulations of quantum many-body systems and quantum field theories that lie beyond the capabilities of classical computation. This will lead to new insights and predictions for systems ranging from dense non-equilibrium matter, to low-energy nuclear structure and reactions, to high-energy collisions. I present an overview of digital quantum simulations in nuclear physics, with select examples relevant for studies of quark matter. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Conference proceedings for Quark Matter 2023, Sept 2023, Houston, Texas, USA. 6 pages and 2 figures

Report number: IQuS@UW-21-066

Scalable Circuits for Preparing Ground States on Digital Quantum Computers: The Schwinger Model Vacuum on 100 Qubits

Authors: Roland C. Farrell, Marc Illa, Anthony N. Ciavarella, Martin J. Savage

Abstract: The vacuum of the lattice Schwinger model is prepared on up to 100 qubits of IBM's Eagle-processor quantum computers. A new algorithm to prepare the ground state of a gapped translationally-invariant system on a quantum computer is presented, which we call Scalable Circuits ADAPT-VQE (SC-ADAPT-VQE). This algorithm uses the exponential decay of correlations between distant regions of the ground sta… ▽ More The vacuum of the lattice Schwinger model is prepared on up to 100 qubits of IBM's Eagle-processor quantum computers. A new algorithm to prepare the ground state of a gapped translationally-invariant system on a quantum computer is presented, which we call Scalable Circuits ADAPT-VQE (SC-ADAPT-VQE). This algorithm uses the exponential decay of correlations between distant regions of the ground state, together with ADAPT-VQE, to construct quantum circuits for state preparation that can be scaled to arbitrarily large systems. These scalable circuits can be determined using classical computers, avoiding the challenging task of optimizing parameterized circuits on a quantum computer. SC-ADAPT-VQE is applied to the Schwinger model, and shown to be systematically improvable, with an accuracy that converges exponentially with circuit depth. Both the structure of the circuits and the deviations of prepared wavefunctions are found to become independent of the number of spatial sites, $L$. This allows for a controlled extrapolation of the circuits, determined using small or modest-sized systems, to arbitrarily large $L$. The circuits for the Schwinger model are determined on lattices up to $L=14$ (28 qubits) with the qiskit classical simulator, and subsequently scaled up to prepare the $L=50$ (100 qubits) vacuum on IBM's 127 superconducting-qubit quantum computers ibm_brisbane and ibm_cusco. After introducing an improved error-mitigation technique, which we call Operator Decoherence Renormalization, the chiral condensate and charge-charge correlators obtained from the quantum computers are found to be in good agreement with classical Matrix Product State simulations. △ Less

Submitted 22 April, 2024; v1 submitted 8 August, 2023; originally announced August 2023.

Comments: 15 pages + appendices. 17 figures, 13 tables

Report number: IQuS@UW-21-060, NT@UW-23-13

Journal ref: PRX Quantum 5 (2024) 2, 020315

Multi-Body Entanglement and Information Rearrangement in Nuclear Many-Body Systems

Authors: S. Momme Hengstenberg, Caroline E. P. Robin, Martin J. Savage

Abstract: We examine how effective-model-space (EMS) calculations of nuclear many-body systems rearrange and converge multi-particle entanglement. The generalized Lipkin-Meshkov-Glick (LMG) model is used to motivate and provide insight for future developments of entanglement-driven descriptions of nuclei. The effective approach is based on a truncation of the Hilbert space together with a variational rotati… ▽ More We examine how effective-model-space (EMS) calculations of nuclear many-body systems rearrange and converge multi-particle entanglement. The generalized Lipkin-Meshkov-Glick (LMG) model is used to motivate and provide insight for future developments of entanglement-driven descriptions of nuclei. The effective approach is based on a truncation of the Hilbert space together with a variational rotation of the qubits (spins), which constitute the relevant elementary degrees of freedom. The non-commutivity of the rotation and truncation allows for an exponential improvement of the energy convergence throughout much of the model space. Our analysis examines measures of correlations and entanglement, and quantifies their convergence with increasing cut-off. We focus on one- and two-spin entanglement entropies, mutual information, and $n$-tangles for $n=2,4$ to estimate multi-body entanglement. The effective description strongly suppresses entropies and mutual information of the rotated spins, while being able to recover the exact results to a large extent with low cut-offs. Naive truncations of the bare Hamiltonian, on the other hand, artificially underestimate these measures. The $n$-tangles in the present model provide a basis-independent measures of $n$-particle entanglement. While these are more difficult to capture with the EMS description, the improvement in convergence, compared to truncations of the bare Hamiltonian, is significantly more dramatic. We conclude that the low-energy EMS techniques, that successfully provide predictive capabilities for low-lying observables in many-body systems, exhibit analogous efficacy for quantum correlations and multi-body entanglement in the LMG model, motivating future studies in nuclear many-body systems and effective field theories relevant to high-energy physics and nuclear physics. △ Less

Submitted 29 December, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 24 pages, 14 figures

Report number: IQuS@UW-21-053

Journal ref: Eur. Phys. J. A 59, 231 (2023)

Quantum Simulations of SO(5) Many-Fermion Systems using Qudits

Authors: Marc Illa, Caroline E. P. Robin, Martin J. Savage

Abstract: The structure and dynamics of quantum many-body systems are the result of a delicate interplay between underlying interactions, which leads to intricate entanglement structures. Despite this apparent complexity, symmetries emerge and have long been used to determine the relevant degrees of freedom and simplify classical descriptions of these systems. In this work, we explore the potential utility… ▽ More The structure and dynamics of quantum many-body systems are the result of a delicate interplay between underlying interactions, which leads to intricate entanglement structures. Despite this apparent complexity, symmetries emerge and have long been used to determine the relevant degrees of freedom and simplify classical descriptions of these systems. In this work, we explore the potential utility of quantum computers with arrays of qudits in simulating interacting fermionic systems, when the qudits can naturally map these relevant degrees of freedom. The Agassi model of fermions is based on an underlying $so(5)$ algebra, and the systems it describes can be partitioned into pairs of modes with five basis states, which naturally embed in arrays of $d=5$ qudits (qu5its). Classical noiseless simulations of the time evolution of systems of fermions embedded in up to twelve qu5its are performed using Google's cirq software. The resource requirements of the qu5it circuits are analyzed and compared with two different mappings to qubit systems, a physics-aware Jordan-Wigner mapping and a state-to-state mapping. We find advantages in using qudits, specifically in lowering the required quantum resources and reducing anticipated errors that take the simulation out of the physical space. A previously unrecognized sign problem has been identified from Trotterization errors in time evolving high-energy excitations. This has implications for quantum simulations in high-energy and nuclear physics, specifically of fragmentation and highly inelastic, multi-channel processes. △ Less

Submitted 29 December, 2023; v1 submitted 19 May, 2023; originally announced May 2023.

Comments: 34 pages, 21 figures

Report number: IQuS@UW-21-049

Journal ref: Phys. Rev. C 108, 064306 (2023)

Quantum Information Science and Technology for Nuclear Physics. Input into U.S. Long-Range Planning, 2023

Authors: Douglas Beck, Joseph Carlson, Zohreh Davoudi, Joseph Formaggio, Sofia Quaglioni, Martin Savage, Joao Barata, Tanmoy Bhattacharya, Michael Bishof, Ian Cloet, Andrea Delgado, Michael DeMarco, Caleb Fink, Adrien Florio, Marianne Francois, Dorota Grabowska, Shannon Hoogerheide, Mengyao Huang, Kazuki Ikeda, Marc Illa, Kyungseon Joo, Dmitri Kharzeev, Karol Kowalski, Wai Kin Lai, Kyle Leach , et al. (76 additional authors not shown)

Abstract: In preparation for the 2023 NSAC Long Range Plan (LRP), members of the Nuclear Science community gathered to discuss the current state of, and plans for further leveraging opportunities in, QIST in NP research at the Quantum Information Science for U.S. Nuclear Physics Long Range Planning workshop, held in Santa Fe, New Mexico on January 31 - February 1, 2023. The workshop included 45 in-person pa… ▽ More In preparation for the 2023 NSAC Long Range Plan (LRP), members of the Nuclear Science community gathered to discuss the current state of, and plans for further leveraging opportunities in, QIST in NP research at the Quantum Information Science for U.S. Nuclear Physics Long Range Planning workshop, held in Santa Fe, New Mexico on January 31 - February 1, 2023. The workshop included 45 in-person participants and 53 remote attendees. The outcome of the workshop identified strategic plans and requirements for the next 5-10 years to advance quantum sensing and quantum simulations within NP, and to develop a diverse quantum-ready workforce. The plans include resolutions endorsed by the participants to address the compelling scientific opportunities at the intersections of NP and QIST. These endorsements are aligned with similar affirmations by the LRP Computational Nuclear Physics and AI/ML Workshop, the Nuclear Structure, Reactions, and Astrophysics LRP Town Hall, and the Fundamental Symmetries, Neutrons, and Neutrinos LRP Town Hall communities. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: A white paper for the 2023 nuclear physics long-range planning activity, emerging from the workshop "Quantum Information Science for U.S. Nuclear Physics Long Range Planning'', held in Santa Fe, New Mexico on January 31 - February 1, 2023. 26 pages with 7 figures

Quantum Simulations in Effective Model Spaces (I): Hamiltonian Learning-VQE using Digital Quantum Computers and Application to the Lipkin-Meshkov-Glick Model

Authors: Caroline E. P. Robin, Martin J. Savage

Abstract: The utility of effective model spaces in quantum simulations of non-relativistic quantum many-body systems is explored in the context of the Lipkin-Meshkov-Glick model of interacting fermions. We introduce an iterative hybrid-classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entang… ▽ More The utility of effective model spaces in quantum simulations of non-relativistic quantum many-body systems is explored in the context of the Lipkin-Meshkov-Glick model of interacting fermions. We introduce an iterative hybrid-classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entanglement into the effective model space, and the associated ground-state wavefunction. HL-VQE is found to provide an exponential improvement in Lipkin-Meshkov-Glick model calculations, compared to a naive truncation without Hamiltonian learning, throughout a significant fraction of the Hilbert space. Quantum simulations are performed to demonstrate the HL-VQE algorithm, using an efficient mapping where the number of qubits scales with the $\log$ of the size of the effective model space, rather than the particle number, allowing for the description of large systems with small quantum circuits. Implementations on IBM's QExperience quantum computers and simulators for 1- and 2-qubit effective model spaces are shown to provide accurate and precise results, reproducing classical predictions. This work constitutes a step in the development of entanglement-driven quantum algorithms for the description of nuclear systems, that leverages the potential of noisy intermediate-scale quantum (NISQ) devices. △ Less

Submitted 23 August, 2023; v1 submitted 14 January, 2023; originally announced January 2023.

Comments: 30 pages, 17 figures, v3/v4: minor modifications

Report number: IQuS@UW-21-041

Journal ref: Phys. Rev. C 108, 024313 (2023)

Preparation for Quantum Simulation of the 1+1D O(3) Non-linear σ-Model using Cold Atoms

Authors: Anthony N. Ciavarella, Stephan Caspar, Hersh Singh, Martin J. Savage

Abstract: The 1+1D O(3) non-linear σ-model is a model system for future quantum lattice simulations of other asymptotically-free theories, such as non-Abelian gauge theories. We find that utilizing dimensional reduction can make efficient use of two-dimensional layouts presently available on cold atom quantum simulators. A new definition of the renormalized coupling is introduced, which is applicable to sys… ▽ More The 1+1D O(3) non-linear σ-model is a model system for future quantum lattice simulations of other asymptotically-free theories, such as non-Abelian gauge theories. We find that utilizing dimensional reduction can make efficient use of two-dimensional layouts presently available on cold atom quantum simulators. A new definition of the renormalized coupling is introduced, which is applicable to systems with open boundary conditions and can be measured using analog quantum simulators. Monte Carlo and tensor network calculations are performed to determine the quantum resources required to reproduce perturbative short-distance observables. In particular, we show that a rectangular array of 48 Rydberg atoms with existing quantum hardware capabilities should be able to adiabatically prepare low-energy states of the perturbatively-matched theory. These states can then be used to simulate non-perturbative observables in the continuum limit that lie beyond the reach of classical computers. △ Less

Submitted 4 April, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 12 pages, 5 figures, 2 tables, published version

Report number: IQuS@UW-21-038, INT-PUB-22-051

Accelerating Progress Towards Practical Quantum Advantage: The Quantum Technology Demonstration Project Roadmap

Authors: Paul Alsing, Phil Battle, Joshua C. Bienfang, Tammie Borders, Tina Brower-Thomas, Lincoln D. Carr, Fred Chong, Siamak Dadras, Brian DeMarco, Ivan Deutsch, Eden Figueroa, Danna Freedman, Henry Everitt, Daniel Gauthier, Ezekiel Johnston-Halperin, Jungsang Kim, Mackillo Kira, Prem Kumar, Paul Kwiat, John Lekki, Anjul Loiacono, Marko Loncar, John R. Lowell, Mikhail Lukin, Celia Merzbacher , et al. (10 additional authors not shown)

Abstract: Quantum information science and technology (QIST) is a critical and emerging technology with the potential for enormous world impact and is currently invested in by over 40 nations. To bring these large-scale investments to fruition and bridge the lower technology readiness levels (TRLs) of fundamental research at universities to the high TRLs necessary to realize the promise of practical quantum… ▽ More Quantum information science and technology (QIST) is a critical and emerging technology with the potential for enormous world impact and is currently invested in by over 40 nations. To bring these large-scale investments to fruition and bridge the lower technology readiness levels (TRLs) of fundamental research at universities to the high TRLs necessary to realize the promise of practical quantum advantage accessible to industry and the public, we present a roadmap for Quantum Technology Demonstration Projects (QTDPs). Such QTDPs, focused on intermediate TRLs, are large-scale public-private partnerships with a high probability of translation from laboratory to practice. They create technology demonstrating a clear 'quantum advantage' for science breakthroughs that are user-motivated and will provide access to a broad and diverse community of scientific users. Successful implementation of a program of QTDPs will have large positive economic impacts. △ Less

Submitted 20 March, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: Updated title and abstract

Multi-Neutrino Entanglement and Correlations in Dense Neutrino Systems

Authors: Marc Illa, Martin J. Savage

Abstract: The time-evolution of multi-neutrino entanglement and correlations are studied in two-flavor collective neutrino oscillations, relevant for dense neutrino environments, building upon previous works. Specifically, simulations performed of systems with up to 12 neutrinos using Quantinuum's H1-1 20 qubit trapped-ion quantum computer are used to compute n-tangles, and two- and three-body correlations,… ▽ More The time-evolution of multi-neutrino entanglement and correlations are studied in two-flavor collective neutrino oscillations, relevant for dense neutrino environments, building upon previous works. Specifically, simulations performed of systems with up to 12 neutrinos using Quantinuum's H1-1 20 qubit trapped-ion quantum computer are used to compute n-tangles, and two- and three-body correlations, probing beyond mean-field descriptions. n-tangle re-scalings are found to converge for large system sizes, signaling the presence of genuine multi-neutrino entanglement. △ Less

Submitted 31 May, 2023; v1 submitted 16 October, 2022; originally announced October 2022.

Comments: 7 pages + 10 pages of Supplementary Material, 9 figures and 11 tables; modified Sec 3; Published version

Report number: IQuS@UW-21-034

Journal ref: Phys. Rev. Lett. 130, 221003 (2023)

The Future of High Energy Physics Software and Computing

Authors: V. Daniel Elvira, Steven Gottlieb, Oliver Gutsche, Benjamin Nachman, S. Bailey, W. Bhimji, P. Boyle, G. Cerati, M. Carrasco Kind, K. Cranmer, G. Davies, V. D. Elvira, R. Gardner, K. Heitmann, M. Hildreth, W. Hopkins, T. Humble, M. Lin, P. Onyisi, J. Qiang, K. Pedro, G. Perdue, A. Roberts, M. Savage, P. Shanahan , et al. (3 additional authors not shown)

Abstract: Software and Computing (S&C) are essential to all High Energy Physics (HEP) experiments and many theoretical studies. The size and complexity of S&C are now commensurate with that of experimental instruments, playing a critical role in experimental design, data acquisition/instrumental control, reconstruction, and analysis. Furthermore, S&C often plays a leading role in driving the precision of th… ▽ More Software and Computing (S&C) are essential to all High Energy Physics (HEP) experiments and many theoretical studies. The size and complexity of S&C are now commensurate with that of experimental instruments, playing a critical role in experimental design, data acquisition/instrumental control, reconstruction, and analysis. Furthermore, S&C often plays a leading role in driving the precision of theoretical calculations and simulations. Within this central role in HEP, S&C has been immensely successful over the last decade. This report looks forward to the next decade and beyond, in the context of the 2021 Particle Physics Community Planning Exercise ("Snowmass") organized by the Division of Particles and Fields (DPF) of the American Physical Society. △ Less

Submitted 8 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: Computational Frontier Report Contribution to Snowmass 2021; 41 pages, 1 figure. v2: missing ref and added missing topical group conveners. v3: fixed typos

State Preparation in the Heisenberg Model through Adiabatic Spiraling

Authors: Anthony N. Ciavarella, Stephan Caspar, Marc Illa, Martin J. Savage

Abstract: An adiabatic state preparation technique, called the adiabatic spiral, is proposed for the Heisenberg model. This technique is suitable for implementation on a number of quantum simulation platforms such as Rydberg atoms, trapped ions, or superconducting qubits. Classical simulations of small systems suggest that it can be successfully implemented in the near future. A comparison to Trotterized ti… ▽ More An adiabatic state preparation technique, called the adiabatic spiral, is proposed for the Heisenberg model. This technique is suitable for implementation on a number of quantum simulation platforms such as Rydberg atoms, trapped ions, or superconducting qubits. Classical simulations of small systems suggest that it can be successfully implemented in the near future. A comparison to Trotterized time evolution is performed and it is shown that the adiabatic spiral is able to outperform Trotterized adiabatics. △ Less

Submitted 1 June, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 22 pages, 8 figures, published version, fixed missing acknowledgment

Report number: IQuS@UW-21-033

Journal ref: Quantum 7, 970 (2023)

Preparations for Quantum Simulations of Quantum Chromodynamics in 1+1 Dimensions: (II) Single-Baryon $β$-Decay in Real Time

Authors: Roland C. Farrell, Ivan A. Chernyshev, Sarah J. M. Powell, Nikita A. Zemlevskiy, Marc Illa, Martin J. Savage

Abstract: A framework for quantum simulations of real-time weak decays of hadrons and nuclei in a 2-flavor lattice theory in one spatial dimension is presented. A single generation of the Standard Model is found to require 16 qubits per spatial lattice site after mapping to spin operators via the Jordan-Wigner transformation. Both quantum chromodynamics and flavor-changing weak interactions are included in… ▽ More A framework for quantum simulations of real-time weak decays of hadrons and nuclei in a 2-flavor lattice theory in one spatial dimension is presented. A single generation of the Standard Model is found to require 16 qubits per spatial lattice site after mapping to spin operators via the Jordan-Wigner transformation. Both quantum chromodynamics and flavor-changing weak interactions are included in the dynamics, the latter through four-Fermi effective operators. Quantum circuits which implement time evolution in this lattice theory are developed and run on Quantinuum's H1-1 20-qubit trapped ion system to simulate the $β$-decay of a single baryon on one lattice site. These simulations include the initial state preparation and are performed for both one and two Trotter time steps. The potential intrinsic error-correction properties of this type of lattice theory are discussed and the leading lattice Hamiltonian required to simulate $0νββ$-decay of nuclei induced by a neutrino Majorana mass term is provided. △ Less

Submitted 10 April, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 27 pages, 11 figures. Typo fixed in v3

Report number: IQuS@UW-21-030, NT@UW-22-12

Journal ref: Phys. Rev. D 107, 054513 (2023)

Report of the Snowmass 2021 Topical Group on Lattice Gauge Theory

Authors: Zohreh Davoudi, Ethan T. Neil, Christian W. Bauer, Tanmoy Bhattacharya, Thomas Blum, Peter Boyle, Richard C. Brower, Simon Catterall, Norman H. Christ, Vincenzo Cirigliano, Gilberto Colangelo, Carleton DeTar, William Detmold, Robert G. Edwards, Aida X. El-Khadra, Steven Gottlieb, Rajan Gupta, Daniel C. Hackett, Anna Hasenfratz, Taku Izubuchi, William I. Jay, Luchang Jin, Christopher Kelly, Andreas S. Kronfeld, Christoph Lehner , et al. (13 additional authors not shown)

Abstract: Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure… ▽ More Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure and spectrum, to serve as a numerical laboratory to reach beyond the Standard Model, or to invent and improve state-of-the-art computational paradigms, the lattice-gauge-theory program is in a prime position to impact the course of developments and enhance discovery potential of a vibrant experimental program in High-Energy Physics over the coming decade. This projection is based on abundant successful results that have emerged using lattice gauge theory over the years: on continued improvement in theoretical frameworks and algorithmic suits; on the forthcoming transition into the exascale era of high-performance computing; and on a skillful, dedicated, and organized community of lattice gauge theorists in the U.S. and worldwide. The prospects of this effort in pushing the frontiers of research in High-Energy Physics have recently been studied within the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021), and the conclusions are summarized in this Topical Report. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 57 pages, 1 figure. Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Topical Group Report for TF05 - Lattice Gauge Theory

Report number: UMD-PP-022-08, LA-UR-22-29361, FERMILAB-CONF-22-703-T

Snowmass Computational Frontier: Topical Group Report on Quantum Computing

Authors: Travis S. Humble, Gabriel N. Perdue, Martin J. Savage

Abstract: Quantum computing will play a pivotal role in the High Energy Physics (HEP) science program over the early parts of the 21$^{st}$ Century, both as a major expansion of our capabilities across the Computational Frontier, and in synthesis with quantum sensing and quantum networks. This report outlines how Quantum Information Science (QIS) and HEP are deeply intertwined endeavors that benefit enormou… ▽ More Quantum computing will play a pivotal role in the High Energy Physics (HEP) science program over the early parts of the 21$^{st}$ Century, both as a major expansion of our capabilities across the Computational Frontier, and in synthesis with quantum sensing and quantum networks. This report outlines how Quantum Information Science (QIS) and HEP are deeply intertwined endeavors that benefit enormously from a strong engagement together. Quantum computers do not represent a detour for HEP, rather they are set to become an integral part of our discovery toolkit. Problems ranging from simulating quantum field theories, to fully leveraging the most sensitive sensor suites for new particle searches, and even data analysis will run into limiting bottlenecks if constrained to our current computing paradigms. Easy access to quantum computers is needed to build a deeper understanding of these opportunities. In turn, HEP brings crucial expertise to the national quantum ecosystem in quantum domain knowledge, superconducting technology, cryogenic and fast microelectronics, and massive-scale project management. The role of quantum technologies across the entire economy is expected to grow rapidly over the next decade, so it is important to establish the role of HEP in the efforts surrounding QIS. Fully delivering on the promise of quantum technologies in the HEP science program requires robust support. It is important to both invest in the co-design opportunities afforded by the broader quantum computing ecosystem and leverage HEP strengths with the goal of designing quantum computers tailored to HEP science. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 50 pages

Report number: FERMILAB-PUB-22-437-QIS,IQuS@UW-21-028

Simulating Heisenberg Interactions in the Ising Model with Strong Drive Fields

Authors: Anthony N. Ciavarella, Stephan Caspar, Hersh Singh, Martin J. Savage, Pavel Lougovski

Abstract: The time-evolution of an Ising model with large driving fields over discrete time intervals is shown to be reproduced by an effective XXZ-Heisenberg model at leading order in the inverse field strength. For specific orientations of the drive field, the dynamics of the XXX-Heisenberg model is reproduced. These approximate equivalences, valid above a critical driving field strength set by dynamical… ▽ More The time-evolution of an Ising model with large driving fields over discrete time intervals is shown to be reproduced by an effective XXZ-Heisenberg model at leading order in the inverse field strength. For specific orientations of the drive field, the dynamics of the XXX-Heisenberg model is reproduced. These approximate equivalences, valid above a critical driving field strength set by dynamical phase transitions in the Ising model, are expected to enable quantum devices that natively evolve qubits according to the Ising model to simulate more complex systems. △ Less

Submitted 23 October, 2023; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 10 pages, 5 figures, accepted version

Report number: IQuS@UW-21-028, INT-PUB-22-020

Lattice QCD and Particle Physics

Authors: Andreas S. Kronfeld, Tanmoy Bhattacharya, Thomas Blum, Norman H. Christ, Carleton DeTar, William Detmold, Robert Edwards, Anna Hasenfratz, Huey-Wen Lin, Swagato Mukherjee, Konstantinos Orginos, Richard Brower, Vincenzo Cirigliano, Zohreh Davoudi, Bálint Jóo, Chulwoo Jung, Christoph Lehner, Stefan Meinel, Ethan T. Neil, Peter Petreczky, David G. Richards, Alexei Bazavov, Simon Catterall, Jozef J. Dudek, Aida X. El-Khadra , et al. (57 additional authors not shown)

Abstract: Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). △ Less

Submitted 2 October, 2022; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

Report number: FERMILAB-CONF-22-531-T

Preparations for Quantum Simulations of Quantum Chromodynamics in 1+1 Dimensions: (I) Axial Gauge

Authors: Roland C. Farrell, Ivan A. Chernyshev, Sarah J. M. Powell, Nikita A. Zemlevskiy, Marc Illa, Martin J. Savage

Abstract: Tools necessary for quantum simulations of $1+1$ dimensional quantum chromodynamics are developed. When formulated in axial gauge and with two flavors of quarks, this system requires 12 qubits per spatial site with the gauge fields included via non-local interactions. Classical computations and D-Wave's quantum annealer Advantage are used to determine the hadronic spectrum, enabling a decompositio… ▽ More Tools necessary for quantum simulations of $1+1$ dimensional quantum chromodynamics are developed. When formulated in axial gauge and with two flavors of quarks, this system requires 12 qubits per spatial site with the gauge fields included via non-local interactions. Classical computations and D-Wave's quantum annealer Advantage are used to determine the hadronic spectrum, enabling a decomposition of the masses and a study of quark entanglement. Color edge states confined within a screening length of the end of the lattice are found. IBM's 7-qubit quantum computers, ibmq_jakarta and ibm_perth, are used to compute dynamics from the trivial vacuum in one-flavor QCD with one spatial site. More generally, the Hamiltonian and quantum circuits for time evolution of $1+1$ dimensional $SU(N_c)$ gauge theory with $N_f$ flavors of quarks are developed, and the resource requirements for large-scale quantum simulations are estimated. △ Less

Submitted 10 April, 2023; v1 submitted 4 July, 2022; originally announced July 2022.

Comments: 53 pages, 29 figures. Typo fixed in v3

Report number: IQuS@UW-21-027, NT@UW-22-05

Journal ref: Phys. Rev. D 107, 054512 (2023)

Quantum Simulation for High Energy Physics

Authors: Christian W. Bauer, Zohreh Davoudi, A. Baha Balantekin, Tanmoy Bhattacharya, Marcela Carena, Wibe A. de Jong, Patrick Draper, Aida El-Khadra, Nate Gemelke, Masanori Hanada, Dmitri Kharzeev, Henry Lamm, Ying-Ying Li, Junyu Liu, Mikhail Lukin, Yannick Meurice, Christopher Monroe, Benjamin Nachman, Guido Pagano, John Preskill, Enrico Rinaldi, Alessandro Roggero, David I. Santiago, Martin J. Savage, Irfan Siddiqi , et al. (6 additional authors not shown)

Abstract: It is for the first time that Quantum Simulation for High Energy Physics (HEP) is studied in the U.S. decadal particle-physics community planning, and in fact until recently, this was not considered a mainstream topic in the community. This fact speaks of a remarkable rate of growth of this subfield over the past few years, stimulated by the impressive advancements in Quantum Information Sciences… ▽ More It is for the first time that Quantum Simulation for High Energy Physics (HEP) is studied in the U.S. decadal particle-physics community planning, and in fact until recently, this was not considered a mainstream topic in the community. This fact speaks of a remarkable rate of growth of this subfield over the past few years, stimulated by the impressive advancements in Quantum Information Sciences (QIS) and associated technologies over the past decade, and the significant investment in this area by the government and private sectors in the U.S. and other countries. High-energy physicists have quickly identified problems of importance to our understanding of nature at the most fundamental level, from tiniest distances to cosmological extents, that are intractable with classical computers but may benefit from quantum advantage. They have initiated, and continue to carry out, a vigorous program in theory, algorithm, and hardware co-design for simulations of relevance to the HEP mission. This community whitepaper is an attempt to bring this exciting and yet challenging area of research to the spotlight, and to elaborate on what the promises, requirements, challenges, and potential solutions are over the next decade and beyond. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: This is a whitepaper prepared for the topical groups CompF6 (Quantum computing), TF05 (Lattice Gauge Theory), and TF10 (Quantum Information Science) within the Computational Frontier and Theory Frontier of the U.S. Community Study on the Future of Particle Physics (Snowmass 2021). 103 pages and 1 figure

Report number: UMD-PP-022-04, LA-UR-22-22100, RIKEN-iTHEMS-Report-22, FERMILAB-PUB-22-249-SQMS-T, IQuS@UW-21-027, MITRE-21-03848-2

Journal ref: PRX Quantum 4, 027001, 2023

Some Conceptual Aspects of Operator Design for Quantum Simulations of Non-Abelian Lattice Gauge Theories

Authors: Anthony Ciavarella, Natalie Klco, Martin J. Savage

Abstract: In the Kogut-Susskind formulation of lattice gauge theories, a set of quantum numbers resides at the ends of each link to characterize the vertex-local gauge field. We discuss the role of these quantum numbers in propagating correlations and supporting entanglement that ensures each vertex remains gauge invariant, despite time evolution induced by operators with (only) partial access to each verte… ▽ More In the Kogut-Susskind formulation of lattice gauge theories, a set of quantum numbers resides at the ends of each link to characterize the vertex-local gauge field. We discuss the role of these quantum numbers in propagating correlations and supporting entanglement that ensures each vertex remains gauge invariant, despite time evolution induced by operators with (only) partial access to each vertex Hilbert space. Applied to recent proposals for eliminating vertex-local Hilbert spaces in quantum simulation, we describe how the required entanglement is generated via delocalization of the time evolution operator with nearest-neighbor controls. These hybridizations, organized with qudits or qubits, exchange classical operator preprocessing for reductions in quantum resource requirements that extend throughout the lattice volume. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 4 pages, 1 figure, Contribution to proceedings of the 2021 Quantum Simulation for Strong Interactions (QuaSi) Workshops at the InQubator for Quantum Simulation (IQuS)

Report number: IQuS@UW-21-024

Basic Elements for Simulations of Standard Model Physics with Quantum Annealers: Multigrid and Clock States

Authors: Marc Illa, Martin J. Savage

Abstract: We explore the potential of D-Wave's quantum annealers for computing some of the basic components required for quantum simulations of Standard Model physics. By implementing a basic multigrid (including "zooming") and specializing Feynman-clock algorithms, D-Wave's Advantage is used to study harmonic and anharmonic oscillators relevant for lattice scalar field theories and effective field theories… ▽ More We explore the potential of D-Wave's quantum annealers for computing some of the basic components required for quantum simulations of Standard Model physics. By implementing a basic multigrid (including "zooming") and specializing Feynman-clock algorithms, D-Wave's Advantage is used to study harmonic and anharmonic oscillators relevant for lattice scalar field theories and effective field theories, the time evolution of a single plaquette of SU(3) Yang-Mills lattice gauge field theory, and the dynamics of flavor entanglement in four neutrino systems. △ Less

Submitted 14 November, 2022; v1 submitted 24 February, 2022; originally announced February 2022.

Comments: 27 pages + Supplemental Material, 13 figures, 11 tables; improved Section III B; Published version

Report number: IQuS@UW-21-020

Journal ref: Phys. Rev. A 106, 052605 (2022)

Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives

Authors: Ingo Tews, Zohreh Davoudi, Andreas Ekström, Jason D. Holt, Kevin Becker, Raúl Briceño, David J. Dean, William Detmold, Christian Drischler, Thomas Duguet, Evgeny Epelbaum, Ashot Gasparyan, Jambul Gegelia, Jeremy R. Green, Harald W. Grießhammer, Andrew D. Hanlon, Matthias Heinz, Heiko Hergert, Martin Hoferichter, Marc Illa, David Kekejian, Alejandro Kievsky, Sebastian König, Hermann Krebs, Kristina D. Launey , et al. (20 additional authors not shown)

Abstract: This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges… ▽ More This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: Perspective pieces of the virtual INT program 21-1b "Nuclear Forces for Precision Nuclear Physics", 107 pages

Report number: INT-PUB-22-002, LA-UR-22-20419, UMD-PP-022-02, FERMILAB-PUB-22-090-T

Journal ref: Few-Body Systems 63, 67 (2022)

Entanglement Structures in Quantum Field Theories: Negativity Cores and Bound Entanglement in the Vacuum

Authors: Natalie Klco, D. H. Beck, Martin J. Savage

Abstract: The many-body entanglement between two finite (size-$d$) disjoint vacuum regions of non-interacting lattice scalar field theory in one spatial dimension -- a $(d_A \times d_B)_{\rm mixed}$ Gaussian continuous variable system -- is locally transformed into a tensor-product "core" of $(1_A \times 1_B)_{\rm mixed}$ entangled pairs. Accessible entanglement within these core pairs exhibits an exponenti… ▽ More The many-body entanglement between two finite (size-$d$) disjoint vacuum regions of non-interacting lattice scalar field theory in one spatial dimension -- a $(d_A \times d_B)_{\rm mixed}$ Gaussian continuous variable system -- is locally transformed into a tensor-product "core" of $(1_A \times 1_B)_{\rm mixed}$ entangled pairs. Accessible entanglement within these core pairs exhibits an exponential hierarchy, and as such identifies the structure of dominant region modes from which vacuum entanglement could be extracted into a spatially separated pair of quantum detectors. Beyond the core, remaining modes of the "halo" are determined to be AB-separable in isolation, as well as separable from the core. However, state preparation protocols that distribute entanglement in the form of $(1_A \times 1_B)_{\rm mixed}$ core pairs are found to require additional entanglement in the halo that is obscured by classical correlations. This inaccessible (bound) halo entanglement is found to mirror the accessible entanglement, but with a step behavior as the continuum is approached. It remains possible that alternate initialization protocols that do not utilize the exponential hierarchy of core-pair entanglement may require less inaccessible entanglement. Entanglement consolidation is expected to persist in higher dimensions and may aid classical and quantum simulations of asymptotically free gauge field theories, such as quantum chromodynamics. △ Less

Submitted 16 October, 2022; v1 submitted 20 October, 2021; originally announced October 2021.

Comments: 11 pages, 3 figures, 10 pages supplemental material. The manuscript has been re-organized and elaborated to improve clarity

Report number: IQuS@UW-21-012

Hierarchical Qubit Maps and Hierarchical Quantum Error Correction

Authors: Natalie Klco, Martin J. Savage

Abstract: We consider hierarchically implemented quantum error correction (HI-QEC), in which the fidelities of logical qubits are differentially optimized to enhance the capabilities of quantum devices in scientific applications. By employing qubit representations that propagate hierarchies in simulated systems to those in logical qubit noise sensitivities, heterogeneity in the distribution of physical-to-l… ▽ More We consider hierarchically implemented quantum error correction (HI-QEC), in which the fidelities of logical qubits are differentially optimized to enhance the capabilities of quantum devices in scientific applications. By employing qubit representations that propagate hierarchies in simulated systems to those in logical qubit noise sensitivities, heterogeneity in the distribution of physical-to-logical qubits can be systematically structured. For concreteness, we estimate HI-QEC's impact on surface code resources in computing low-energy observables to fixed precision, finding up to $\sim 60\%$ reductions in qubit requirements plausible in early error corrected simulations. Hierarchical qubit maps are also possible without error correction in qubit and qudit systems where fidelities are non-uniform, either unintentionally or by design. Hierarchical optimizations are another element in the co-design process of quantum simulations for nuclear and particle physics. △ Less

Submitted 4 September, 2021; originally announced September 2021.

Comments: 11 pages, 5 figures

Report number: IQuS@UW-21-010

Standard Model Physics and the Digital Quantum Revolution: Thoughts about the Interface

Authors: Natalie Klco, Alessandro Roggero, Martin J. Savage

Abstract: Advances in isolating, controlling and entangling quantum systems are transforming what was once a curious feature of quantum mechanics into a vehicle for disruptive scientific and technological progress. Pursuing the vision articulated by Feynman, a concerted effort across many areas of research and development is introducing prototypical digital quantum devices into the computing ecosystem avail… ▽ More Advances in isolating, controlling and entangling quantum systems are transforming what was once a curious feature of quantum mechanics into a vehicle for disruptive scientific and technological progress. Pursuing the vision articulated by Feynman, a concerted effort across many areas of research and development is introducing prototypical digital quantum devices into the computing ecosystem available to domain scientists. Through interactions with these early quantum devices, the abstract vision of exploring classically-intractable quantum systems is evolving toward becoming a tangible reality. Beyond catalyzing these technological advances, entanglement is enabling parallel progress as a diagnostic for quantum correlations and as an organizational tool, both guiding improved understanding of quantum many-body systems and quantum field theories defining and emerging from the Standard Model. From the perspective of three domain science theorists, this article compiles thoughts about the interface on entanglement, complexity, and quantum simulation in an effort to contextualize recent NISQ-era progress with the scientific objectives of nuclear and high-energy physics. △ Less

Submitted 10 July, 2021; originally announced July 2021.

Comments: 63 pages, 5 figures

Report number: IQuS@UW-21-007

Entanglement Spheres and a UV-IR connection in Effective Field Theories

Authors: Natalie Klco, Martin J. Savage

Abstract: We show that long-distance quantum correlations probe short-distance physics. Two disjoint regions of the latticized, massless scalar field vacuum are numerically demonstrated to become separable at distances beyond the negativity sphere, which extends to infinity in the continuum limit. The size of this quantum coherent volume is determined by the highest momentum mode supported in the identical… ▽ More We show that long-distance quantum correlations probe short-distance physics. Two disjoint regions of the latticized, massless scalar field vacuum are numerically demonstrated to become separable at distances beyond the negativity sphere, which extends to infinity in the continuum limit. The size of this quantum coherent volume is determined by the highest momentum mode supported in the identical regions, each of diameter $d$. More generally, effective field theories (EFTs), describing a system up to a given momentum scale $Λ$, are expected to share this feature -- entanglement between regions of the vacuum depends upon the UV-completion beyond a separation proportional to $Λ$. Through calculations extended to three-dimensions, the magnitude of the negativity at which entanglement becomes sensitive to UV physics in an EFT (lattice or otherwise) is conjectured to scale as $\sim e^{-Λd}$, independent of the number of spatial dimensions. It is concluded that two-region vacuum entanglement at increasing separations depends upon the structure of the theory at increasing momentum scales. This phenomenon may be manifest in perturbative QCD processes. △ Less

Submitted 28 August, 2021; v1 submitted 27 March, 2021; originally announced March 2021.

Comments: 9 pages, 3 figures, 4 appendices (4 figures and 2 tables)

Report number: IQuS@UW-21-005

A Trailhead for Quantum Simulation of SU(3) Yang-Mills Lattice Gauge Theory in the Local Multiplet Basis

Authors: Anthony Ciavarella, Natalie Klco, Martin J. Savage

Abstract: Maintaining local interactions in the quantum simulation of gauge field theories relegates most states in the Hilbert space to be unphysical -- theoretically benign, but experimentally difficult to avoid. Reformulations of the gauge fields can modify the ratio of physical to gauge-variant states often through classically preprocessing the Hilbert space and modifying the representation of the field… ▽ More Maintaining local interactions in the quantum simulation of gauge field theories relegates most states in the Hilbert space to be unphysical -- theoretically benign, but experimentally difficult to avoid. Reformulations of the gauge fields can modify the ratio of physical to gauge-variant states often through classically preprocessing the Hilbert space and modifying the representation of the field on qubit degrees of freedom. This paper considers the implications of representing SU(3) Yang-Mills gauge theory on a lattice of irreducible representations in both a global basis of projected global quantum numbers and a local basis in which controlled-plaquette operators support efficient time evolution. Classically integrating over the internal gauge space at each vertex (e.g., color isospin and color hypercharge) significantly reduces both the qubit requirements and the dimensionality of the unphysical Hilbert space. Initiating tuning procedures that may inform future calculations at scale, the time evolution of one- and two-plaquettes are implemented on one of IBM's superconducting quantum devices, and early benchmark quantities are identified. The potential advantages of qudit environments, with either constrained 2D hexagonal or 1D nearest-neighbor internal state connectivity, are discussed for future large-scale calculations. △ Less

Submitted 21 February, 2021; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: 51 pages, 17 figures (references added, single qubit results re-analyzed, typos fixed)

Report number: IQuS@UW-21-001

Journal ref: Phys. Rev. D 103, 094501 (2021)

Low-energy Scattering and Effective Interactions of Two Baryons at $m_π\sim 450$ MeV from Lattice Quantum Chromodynamics

Authors: Marc Illa, Silas R. Beane, Emmanuel Chang, Zohreh Davoudi, William Detmold, David J. Murphy, Kostas Orginos, Assumpta Parreño, Martin J. Savage, Phiala E. Shanahan, Michael L. Wagman, Frank Winter

Abstract: The interactions between two octet baryons are studied at low energies using lattice QCD (LQCD) with larger-than-physical quark masses corresponding to a pion mass of $m_π\sim 450$ MeV and a kaon mass of $m_{K}\sim 596$ MeV. The two-baryon systems that are analyzed range from strangeness $S=0$ to $S=-4$ and include the spin-singlet and triplet $NN$, $ΣN$ ($I=3/2$), and $ΞΞ$ states, the spin-single… ▽ More The interactions between two octet baryons are studied at low energies using lattice QCD (LQCD) with larger-than-physical quark masses corresponding to a pion mass of $m_π\sim 450$ MeV and a kaon mass of $m_{K}\sim 596$ MeV. The two-baryon systems that are analyzed range from strangeness $S=0$ to $S=-4$ and include the spin-singlet and triplet $NN$, $ΣN$ ($I=3/2$), and $ΞΞ$ states, the spin-singlet $ΣΣ$ ($I=2$) and $ΞΣ$ ($I=3/2$) states, and the spin-triplet $ΞN$ ($I=0$) state. The $s$-wave scattering phase shifts, low-energy scattering parameters, and binding energies when applicable, are extracted using Lüscher's formalism. While the results are consistent with most of the systems being bound at this pion mass, the interactions in the spin-triplet $ΣN$ and $ΞΞ$ channels are found to be repulsive and do not support bound states. Using results from previous studies at a larger pion mass, an extrapolation of the binding energies to the physical point is performed and is compared with experimental values and phenomenological predictions. The low-energy coefficients in pionless EFT relevant for two-baryon interactions, including those responsible for $SU(3)$ flavor-symmetry breaking, are constrained. The $SU(3)$ symmetry is observed to hold approximately at the chosen values of the quark masses, as well as the $SU(6)$ spin-flavor symmetry, predicted at large $N_c$. A remnant of an accidental $SU(16)$ symmetry found previously at a larger pion mass is further observed. The $SU(6)$-symmetric EFT constrained by these LQCD calculations is used to make predictions for two-baryon systems for which the low-energy scattering parameters could not be determined with LQCD directly in this study, and to constrain the coefficients of all leading $SU(3)$ flavor-symmetric interactions, demonstrating the predictive power of two-baryon EFTs matched to LQCD. △ Less

Submitted 23 March, 2021; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: 69 pages, 31 figures and 25 tables; published version

Report number: ICCUB-20-020, UMD-PP-020-7, MIT-CTP/5238, INT-PUB-20-038, FERMILAB-PUB-20-498-T

Journal ref: Phys. Rev. D 103, 054508 (2021)

A Survey of 3 to 5.4 Micron Emission From Planetary Nebulae using SOFIA/FLITECAM

Authors: Erin C. Smith, Sarah E. Logsdon, Ian S. McLean, Elizabeth Fletcher, William D. Vacca, E. E. Becklin, 5, Sachindev Shenoy, Maureen Savage, Ryan T. Hamilton

Abstract: Here we present the results of an airborne 3-5.4$μ$m spectroscopic study of three young, Carbon-rich planetary nebulae IC 5117, PNG 093.9-00.1, and BD $+$30 3639. These observations were made using the grism spectroscopy mode of the FLITECAM instrument during airborne science operations onboard NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA). The goal of this study is to characteri… ▽ More Here we present the results of an airborne 3-5.4$μ$m spectroscopic study of three young, Carbon-rich planetary nebulae IC 5117, PNG 093.9-00.1, and BD $+$30 3639. These observations were made using the grism spectroscopy mode of the FLITECAM instrument during airborne science operations onboard NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA). The goal of this study is to characterize the 3.3 $μ$m and 5.25 $μ$m PAH dust emission in planetary nebulae and study the evolution of PAH features within evolved stars before their incorporation into new stellar systems in star-forming regions. Targets were selected from IRAS, KAO and ISO source lists, and were previously observed with FLITECAM on the 3-meter Shane telescope at Lick Observatory to allow direct comparison between the ground and airborne observations. We measure PAH emission equivalent width and central wavelength, classify the shape of the PAH emission, and determine the PAH/Aliphatic ratio for each target. The 3.3 $μ$m PAH emission feature is observed in all three objects. PNG 093.9-00.1 exhibits NGC 7027-like aliphatic emission in the 3.4 to 3.6 $μ$m region while IC 5117 and BD +30 3639 exhibit less aliphatic structure. All three PNs additionally exhibit PAH emission at 5.25 $μ$m. △ Less

Submitted 31 August, 2020; originally announced August 2020.

Comments: 17 Pages, 3 figures

Nuclear matrix elements from lattice QCD for electroweak and beyond-Standard-Model processes

Authors: Zohreh Davoudi, William Detmold, Kostas Orginos, Assumpta Parreño, Martin J. Savage, Phiala Shanahan, Michael L. Wagman

Abstract: Over the last decade, numerical solutions of Quantum Chromodynamics (QCD) using the technique of lattice QCD have developed to a point where they are beginning to connect fundamental aspects of nuclear physics to the underlying degrees of freedom of the Standard Model. In this review, the progress of lattice QCD studies of nuclear matrix elements of electroweak currents and beyond-Standard-Model o… ▽ More Over the last decade, numerical solutions of Quantum Chromodynamics (QCD) using the technique of lattice QCD have developed to a point where they are beginning to connect fundamental aspects of nuclear physics to the underlying degrees of freedom of the Standard Model. In this review, the progress of lattice QCD studies of nuclear matrix elements of electroweak currents and beyond-Standard-Model operators is summarized, and connections with effective field theories and nuclear models are outlined. Lattice QCD calculations of nuclear matrix elements can provide guidance for low-energy nuclear reactions in astrophysics, dark matter direct detection experiments, and experimental searches for violations of the symmetries of the Standard Model, including searches for additional CP violation in the hadronic and leptonic sectors, baryon-number violation, and lepton-number or flavor violation. Similarly, important inputs to neutrino experiments seeking to determine the neutrino-mass hierarchy and oscillation parameters, as well as other electroweak and beyond-Standard-Model processes can be determined. The phenomenological implications of existing studies of electroweak and beyond-Standard-Model matrix elements in light nuclear systems are discussed, and future prospects for the field toward precision studies of these matrix elements are outlined. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Report number: MIT-CTP/5229, UMD-PP-020-4, ICCUB-20-018, FERMILAB-PUB-20-445-T

Geometric Quantum Information Structure in Quantum Fields and their Lattice Simulation

Authors: Natalie Klco, Martin J. Savage

Abstract: An upper limit to distillable entanglement between two disconnected regions of massless non-interacting scalar field theory has an exponential decay defined by a geometric decay constant. When regulated at short distances with a spatial lattice, this entanglement abruptly vanishes beyond a dimensionless separation, defining a negativity sphere. In two spatial dimensions, we determine this geometri… ▽ More An upper limit to distillable entanglement between two disconnected regions of massless non-interacting scalar field theory has an exponential decay defined by a geometric decay constant. When regulated at short distances with a spatial lattice, this entanglement abruptly vanishes beyond a dimensionless separation, defining a negativity sphere. In two spatial dimensions, we determine this geometric decay constant between a pair of disks and the growth of the negativity sphere toward the continuum through a series of lattice calculations. Making the connection to quantum field theories in three-spatial dimensions, assuming such quantum information scales appear also in quantum chromodynamics (QCD), a new relative scale may be present in effective field theories describing the low-energy dynamics of nucleons and nuclei. We highlight potential impacts of the distillable entanglement structure on effective field theories, lattice QCD calculations and future quantum simulations. △ Less

Submitted 9 August, 2020; originally announced August 2020.

Comments: 9 pages, 3 figures

Report number: INT-PUB-20-031

Journal ref: Phys. Rev. D 103, 065007 (2021)

Entanglement Rearrangement in Self-Consistent Nuclear Structure Calculations

Authors: Caroline Robin, Martin J. Savage, Nathalie Pillet

Abstract: Entanglement properties of $^4$He and $^6$He are investigated using nuclear many-body calculations, specifically the single-nucleon entanglement entropy, and the two-nucleon mutual information and negativity. Nuclear wavefunctions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field t… ▽ More Entanglement properties of $^4$He and $^6$He are investigated using nuclear many-body calculations, specifically the single-nucleon entanglement entropy, and the two-nucleon mutual information and negativity. Nuclear wavefunctions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field theory. Entanglement measures within single-particle bases, the harmonic oscillator (HO), Hartree-Fock (HF), natural (NAT) and variational natural (VNAT) bases, are found to exhibit different degrees of complexity. Entanglement in both nuclei is found to be more localized within NAT and VNAT bases than within a HO basis for the optimal HO parameters, and, as anticipated, a core-valence (tensor product) structure emerges from the full six-body calculation of $^{6}$He. The two-nucleon mutual information shows that the VNAT basis, which typically exhibits good convergence properties, effectively decouples the active and inactive spaces. We conclude that measures of one- and two-nucleon entanglement are useful in analyzing the structure of nuclear wave functions, in particular the efficacy of basis states, and may provide useful metrics toward developing more efficient schemes for ab initio computations of the structure and reactions of nuclei, and quantum many-body systems more generally. △ Less

Submitted 12 May, 2021; v1 submitted 17 July, 2020; originally announced July 2020.

Comments: 18 pages, 10 figures, v2: accepted version

Report number: INT-PUB-20-029

Journal ref: Phys. Rev. C 103, 034325 (2021)

Fixed-Point Quantum Circuits for Quantum Field Theories

Authors: Natalie Klco, Martin J. Savage

Abstract: 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 f… ▽ More 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. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: 9 pages, 3 figures

Report number: INT-PUB-20-003

Quantum Computer Systems for Scientific Discovery

Authors: Yuri Alexeev, Dave Bacon, Kenneth R. Brown, Robert Calderbank, Lincoln D. Carr, Frederic T. Chong, Brian DeMarco, Dirk Englund, Edward Farhi, Bill Fefferman, Alexey V. Gorshkov, Andrew Houck, Jungsang Kim, Shelby Kimmel, Michael Lange, Seth Lloyd, Mikhail D. Lukin, Dmitri Maslov, Peter Maunz, Christopher Monroe, John Preskill, Martin Roetteler, Martin Savage, Jeff Thompson

Abstract: 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 commu… ▽ More 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. △ Less

Submitted 29 July, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Journal ref: PRX Quantum 2, 017001 (2021)

Systematically Localizable Operators for Quantum Simulations of Quantum Field Theories

Authors: Natalie Klco, Martin J. Savage

Abstract: 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 lattice scalar quantum field theories. Controlled rotations that are exponentially localized in their pos… ▽ More 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 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 correlation between sites and 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. △ Less

Submitted 4 April, 2020; v1 submitted 7 December, 2019; originally announced December 2019.

Comments: 22 pages, 7 figures, 6 appendices

Report number: INT-PUB-19-058

Journal ref: Phys. Rev. A 102, 012619 (2020)

Sparsening Algorithm for Multi-Hadron Lattice QCD Correlation Functions

Authors: W. Detmold, D. J. Murphy, A. V. Pochinsky, M. J. Savage, P. E. Shanahan, M. L. Wagman

Abstract: Modern advances in algorithms for lattice QCD calculations have steadily driven down the resources required to generate gauge field ensembles and calculate quark propagators, such that, in cases relevant to nuclear physics, performing quark contractions to assemble correlation functions from propagators has become the dominant cost. This work explores a propagator sparsening algorithm for forming… ▽ More Modern advances in algorithms for lattice QCD calculations have steadily driven down the resources required to generate gauge field ensembles and calculate quark propagators, such that, in cases relevant to nuclear physics, performing quark contractions to assemble correlation functions from propagators has become the dominant cost. This work explores a propagator sparsening algorithm for forming correlation functions describing multi-hadron systems, such as light nuclei, with reduced computational cost. The algorithm constructs correlation functions from sparsened propagators defined on a coarsened lattice geometry, where the sparsened propagators are obtained from propagators computed on the full lattice. This algorithm is used to study the low-energy QCD ground-state spectrum using a single Wilson-clover lattice ensemble with $m_π \approx 800$ MeV. It is found that the extracted ground state masses and binding energies, as well as their statistical uncertainties, are consistent when determined from correlation functions constructed from sparsened and full propagators. In addition, while evidence of modified couplings to excited states is observed in sparsened correlation functions, it is demonstrated that these effects can be removed, if desired, with an inexpensive modification to the sparsened estimator. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 18 pages, 5 figures, 4 tables

Report number: MIT-CTP/5127, INT-PUB-19-026

Journal ref: Phys. Rev. D 104, 034502 (2021)

SU(2) non-Abelian gauge field theory in one dimension on digital quantum computers

Authors: Natalie Klco, Jesse R. Stryker, Martin J. Savage

Abstract: 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 r… ▽ More 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. The results of a multi-plaquette calculation on IBM's quantum hardware are presented. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 8 pages, 4 figures, 3 tables

Report number: INT-PUB-19-033

Journal ref: Phys. Rev. D 101, 074512 (2020)

FOBOS: A Next-Generation Spectroscopic Facility at the W. M. Keck Observatory

Authors: K. Bundy, K. Westfall, N. MacDonald, R. Kupke, M. Savage, C. Poppett, A. Alabi, G. Becker, J. Burchett, P. Capak, A. Coil, M. Cooper, D. Cowley, W. Deich, D. Dillon, J. Edelstein, P. Guhathakurta, J. Hennawi, M. Kassis, K. -G. Lee, D. Masters, T. Miller, J. Newman, J. O'Meara, J. X. Prochaska , et al. (16 additional authors not shown)

Abstract: High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term… ▽ More High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: Submitted to the Astro2020 Decadal Survey as a facilities white paper

Minimally-Entangled State Preparation of Localized Wavefunctions on Quantum Computers

Authors: Natalie Klco, Martin J. Savage

Abstract: Initializing a single site of a lattice scalar field theory into an arbitrary state with support throughout the quantum register requires ${\cal O}(2^{n_Q})$ 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… ▽ More Initializing a single site of a lattice scalar field theory into an arbitrary state with support throughout the quantum register requires ${\cal O}(2^{n_Q})$ 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^{n_Q-1} + n_Q-3 + δ_{n_Q,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 attainable 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 improvement and \emph{Somma Inflation}. 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 ($|{\bf 0}\rangle^{\otimes n_Q}$) 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. △ Less

Submitted 28 April, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: 47 pages, 17 figures

Report number: INT-PUB-19-013

Journal ref: Phys. Rev. A 102, 012612 (2020)

Status and Future Perspectives for Lattice Gauge Theory Calculations to the Exascale and Beyond

Authors: Bálint Joó, Chulwoo Jung, Norman H. Christ, William Detmold, Robert G. Edwards, Martin Savage, Phiala Shanahan

Abstract: In this and a set of companion whitepapers, the USQCD Collaboration lays out a program of science and computing for lattice gauge theory. These whitepapers describe how calculation using lattice QCD (and other gauge theories) can aid the interpretation of ongoing and upcoming experiments in particle and nuclear physics, as well as inspire new ones. In this and a set of companion whitepapers, the USQCD Collaboration lays out a program of science and computing for lattice gauge theory. These whitepapers describe how calculation using lattice QCD (and other gauge theories) can aid the interpretation of ongoing and upcoming experiments in particle and nuclear physics, as well as inspire new ones. △ Less

Submitted 22 November, 2019; v1 submitted 22 April, 2019; originally announced April 2019.

Comments: 44 pages. 1 of USQCD whitepapers,

Journal ref: Eur. Phys. J. A (2019) 55: 199

Branching principles of animal and plant networks identified by combining extensive data, machine learning, and modeling

Authors: Alexander B Brummer, Panagiotis Lymperopoulos, Jocelyn Shen, Elif Tekin, Lisa P. Bentley, Vanessa Buzzard, Andrew Gray, Imma Oliveras, Brian J. Enquist, Van M. Savage

Abstract: Branching in vascular networks and in overall organismic form is one of the most common and ancient features of multicellular plants, fungi, and animals. By combining machine-learning techniques with new theory that relates vascular form to metabolic function, we enable novel classification of diverse branching networks--mouse lung, human head and torso, angiosperm and gymnosperm plants. We find t… ▽ More Branching in vascular networks and in overall organismic form is one of the most common and ancient features of multicellular plants, fungi, and animals. By combining machine-learning techniques with new theory that relates vascular form to metabolic function, we enable novel classification of diverse branching networks--mouse lung, human head and torso, angiosperm and gymnosperm plants. We find that ratios of limb radii--which dictate essential biologic functions related to resource transport and supply--are best at distinguishing branching networks. We also show how variation in vascular and branching geometry persists despite observing a convergent relationship across organisms for how metabolic rate depends on body mass. △ Less

Submitted 6 August, 2020; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 55 pages, 8 figures, 8 tables

Entanglement Suppression and Emergent Symmetries of Strong Interactions

Authors: Silas R. Beane, David B. Kaplan, Natalie Klco, Martin J. Savage

Abstract: 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 thes… ▽ More 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. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: 6 pages, 5 pages supplemental, 4 figures

Report number: INT-PUB-18-056, NT@UW-18-19

Journal ref: Phys. Rev. Lett. 122, 102001 (2019)

Theoretical aspects of quantum electrodynamics in a finite volume with periodic boundary conditions

Authors: Zohreh Davoudi, James Harrison, Andreas Jüttner, Antonin Portelli, Martin J. Savage

Abstract: First-principles studies of strongly-interacting hadronic systems using lattice quantum chromodynamics (QCD) have been complemented in recent years with the inclusion of quantum electrodynamics (QED). The aim is to confront experimental results with more precise theoretical determinations, e.g. for the anomalous magnetic moment of the muon and the CP-violating parameters in the decay of mesons. Qu… ▽ More First-principles studies of strongly-interacting hadronic systems using lattice quantum chromodynamics (QCD) have been complemented in recent years with the inclusion of quantum electrodynamics (QED). The aim is to confront experimental results with more precise theoretical determinations, e.g. for the anomalous magnetic moment of the muon and the CP-violating parameters in the decay of mesons. Quantifying the effects arising from enclosing QED in a finite volume remains a primary target of investigations. To this end, finite-volume corrections to hadron masses in the presence of QED have been carefully studied in recent years. This paper extends such studies to the self-energy of moving charged hadrons, both on and away from their mass shell. In particular, we present analytical results for leading finite-volume corrections to the self-energy of spin-0 and spin-$\frac{1}{2}$ particles in the presence of QED on a periodic hypercubic lattice, once the spatial zero mode of the photon is removed, a framework that is called $\mathrm{QED}_{\mathrm{L}}$. By altering modes beyond the zero mode, an improvement scheme is introduced to eliminate the leading finite-volume corrections to masses, with potential applications to other hadronic quantities. Our analytical results are verified by a dedicated numerical study of a lattice scalar field theory coupled to $\mathrm{QED}_{\mathrm{L}}$. Further, this paper offers new perspectives on the subtleties involved in applying low-energy effective field theories in the presence of $\mathrm{QED}_{\mathrm{L}}$, a theory that is rendered non-local with the exclusion of the spatial zero mode of the photon, clarifying recent discussions on this matter. △ Less

Submitted 19 March, 2019; v1 submitted 13 October, 2018; originally announced October 2018.

Comments: 57 pages, 10 figures, version accepted for publication in Phys. Rev. D

Report number: INT-PUB-18-051, UMD-PP-018-07

Journal ref: Phys. Rev. D 99, 034510 (2019)

Simulations of Subatomic Many-Body Physics on a Quantum Frequency Processor

Authors: Hsuan-Hao Lu, Natalie Klco, Joseph M. Lukens, Titus D. Morris, Aaina Bansal, Andreas Ekström, Gaute Hagen, Thomas Papenbrock, Andrew M. Weiner, Martin J. Savage, Pavel Lougovski

Abstract: Simulating complex many-body quantum phenomena is a major scientific impetus behind the development of quantum computing, and a range of technologies are being explored to address such systems. We present the results of the largest photonics-based simulation to date, applied in the context of subatomic physics. Using an all-optical quantum frequency processor, the ground-state energies of light nu… ▽ More Simulating complex many-body quantum phenomena is a major scientific impetus behind the development of quantum computing, and a range of technologies are being explored to address such systems. We present the results of the largest photonics-based simulation to date, applied in the context of subatomic physics. Using an all-optical quantum frequency processor, the ground-state energies of light nuclei including the triton ($^3$H), $^{3}$He, and the alpha particle ($^{4}$He) are computed. Complementing these calculations and utilizing a 68-dimensional Hilbert space, our photonic simulator is used to perform sub-nucleon calculations of the two-body and three-body forces between heavy mesons in the Schwinger model. This work is a first step in simulating subatomic many-body physics on quantum frequency processors---augmenting classical computations that bridge scales from quarks to nuclei. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 23 pages, 5 figures, 20 pages supplemental material

Journal ref: Phys. Rev. A 100, 012320 (2019)

Digitization of Scalar Fields for Quantum Computing

Authors: Natalie Klco, Martin J. Savage

Abstract: Qubit, operator and gate resources required for the digitization of lattice $λφ^4$ scalar field theories onto quantum computers are considered, building upon the foundational work by Jordan, Lee and Preskill, with a focus towards noisy intermediate-scale quantum (NISQ) devices. The Nyquist-Shannon sampling theorem, introduced in this context by Macridin, Spentzouris, Amundson and Harnik building o… ▽ More Qubit, operator and gate resources required for the digitization of lattice $λφ^4$ scalar field theories onto quantum computers are considered, building upon the foundational work by Jordan, Lee and Preskill, with a focus towards noisy intermediate-scale quantum (NISQ) devices. The Nyquist-Shannon sampling theorem, introduced in this context by Macridin, Spentzouris, 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 $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 $λφ^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 |ε_{\rm dig}|||\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 |ε_{\rm latt}||\sim N_Q$ (total number of qubits in the simulation). For localized(delocalized) field-space wavefunctions, it is found that $n_Q\sim4(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. △ Less

Submitted 24 May, 2019; v1 submitted 30 August, 2018; originally announced August 2018.

Comments: 48 pages, 19 figures, 4 tables

Report number: INT-PUB-18-044

Journal ref: Phys. Rev. A 99, 052335 (2019)