Abstract:With a view toward addressing the explosive growth in the computational demands of nuclear structure and reactions modeling, we develop a novel quantum algorithm for neutron-nucleus simulations with general potentials, which provides acceptable bound-state energies even in the presence of noise, through the noise-resilient training method. In particular, the algorithm can now solve for any band-diagonal to full Hamiltonian matrices, as needed to accommodate a general central potential. This includes exponential Gaussian-like potentials and ab initio inter-cluster potentials (optical potentials). The approach can also accommodate the complete form of the chiral effective-field-theory nucleon-nucleon potentials used in ab initio nuclear calculations. We make this potential available for three different qubit encodings, including the one-hot (OHE), binary (BE), and Gray encodings (GE), and we provide a comprehensive analysis of the number of Pauli terms and commuting sets involved. We find that the GE allows for an efficient scaling of the model-space size $N$ (or number of basis states used) and is more resource efficient not only for tridiagonal Hamiltonians, but also for band-diagonal Hamiltonians having bandwidth up to $N$. We introduce a new commutativity scheme called distance-grouped commutativity (DGC) and compare its performance with the well-known qubit-commutativity (QC) scheme. We lay out the explicit grouping of Pauli strings and the diagonalizing unitary under the DGC scheme, and we find that it outperforms the QC scheme, at the cost of a more complex diagonalizing unitary. Lastly, we provide first solutions of the neutron-alpha dynamics from quantum simulations suitable for NISQ processors, using an optical potential rooted in first principles, and a study of the bound-state physics in neutron-Carbon systems, along with a comparison of the efficacy of the OHE and GE.

Ab Initio Simulations of Light Nuclear Systems Using Eigenvector Continuation and Auxiliary Field Monte Carlo

Ab Initionuclear Structure – the Large Sparse Matrix Eigenvalue Problem

Eigenvector continuation for emulating and extrapolating two-body resonances

Path Integral Quantum Monte Carlo Method for Light Nuclei

Floating block method for quantum Monte Carlo simulations

Eigenvector Continuation with Subspace Learning

Numerical Simulations of the Nonlinear Quantum Vacuum in the Heisenberg-Euler Weak-Field Expansion

Toward emulating nuclear reactions using eigenvector continuation

Wavefunction matching for solving quantum many-body problems

Auxiliary field Quantum Monte Carlo for dilute neutrons on the lattice

Solving nuclear structure problems with the adaptive variational quantum algorithm

Recent progress in the electroweak structure of light nuclei using quantum Monte Carlo methods

Large-scale Ab Initio Configuration Interaction Calculations for Light Nuclei

Hamiltonian Light-Front Field Theory: Recent Progress and Tantalizing Prospects

Neutron-nucleus dynamics simulations for quantum computers

Hydrodynamic and Non-hydrodynamic Excitations in Kinetic Theory -- A Numerical Analysis in Scalar Field Theory

Emulating ab initio computations of infinite nucleonic matter

Efficient Quantum Simulation of QCD Jets on the Light Front

Non-Perturbative Simulations of Quantum Field Theories using Complex Langevin Dynamics

Variational Monte Carlo Calculations of n+3H Scattering