No Arabic abstract
No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, identified by strong intraband E2 transitions and by rotational patterns for excitation energies, electromagnetic moments, and electromagnetic transitions. However, convergence rates differ significantly for different rotational observables and for different rotational bands. The choice of internucleon interaction may also substantially impact the convergence rates. Consequently, there is a substantial gap between simply observing the qualitative emergence of rotation in ab initio calculations and actually carrying out detailed quantitative comparisons. In this contribution, we illustrate the convergence properties of rotational band energy parameters extracted from NCCI calculations, and compare these predictions with experiment, for the isotopes 7-11Be, and for the JISP16 and Daejeon16 interactions.
The extension of ab initio quantum many-body theory to higher accuracy and larger systems is intrinsically limited by the handling of large data objects in form of wave-function expansions and/or many-body operators. In this work we present matrix factorization techniques as a systematically improvable and robust tool to significantly reduce the computational cost in many-body applications at the price of introducing a moderate decomposition error. We demonstrate the power of this approach for the nuclear two-body systems, for many-body perturbation theory calculations of symmetric nuclear matter, and for non-perturbative in-medium similarity renormalization group simulations of finite nuclei. Establishing low-rank expansions of chiral nuclear interactions offers possibilities to reformulate many-body methods in ways that take advantage of tensor factorization strategies.
Background: Solving nuclear many-body problems with an ab initio approach is widely recognized as a computationally challenging problem. Quantum computers offer a promising path to address this challenge. There are urgent needs to develop quantum algorithms for this purpose. Objective: In this work, we explore the application of the quantum algorithm of adiabatic state preparation with quantum phase estimation in ab initio nuclear structure theory. We focus on solving the low-lying spectra (including both the ground and excited states) of simple nuclear systems. Ideas: The efficiency of this algorithm is hindered by the emergence of small energy gaps (level crossings) during the adiabatic evolution. In order to improve the efficiency, we introduce techniques to avoid level crossings: 1) by suitable design of the reference Hamiltonian; 2) by insertions of perturbation terms to modify the adiabatic path. Results: We illustrate this algorithm by solving the deuteron ground state energy and the spectrum of the deuteron bounded in a harmonic oscillator trap implementing the IBM Qiskit quantum simulator. The quantum results agree well the classical results obtained by matrix diagonalization. Outlook: With our improvements to the efficiency, this algorithm provides a promising tool for investigating the low-lying spectra of complex nuclei on future quantum computers.
Nuclear structure models built from phenomenological mean fields, the effective nucleon-nucleon interactions (or Lagrangians), and the realistic bare nucleon-nucleon interactions are reviewed. The success of covariant density functional theory (CDFT) to describe nuclear properties and its influence on Brueckner theory within the relativistic framework are focused upon. The challenges and ambiguities of predictions for unstable nuclei without data or for high-density nuclear matter, arising from relativistic density functionals, are discussed. The basic ideas in building an ab initio relativistic density functional for nuclear structure from ab initio calculations with realistic nucleon-nucleon interactions for both nuclear matter and finite nuclei are presented. The current status of fully self-consistent relativistic Brueckner-Hartree-Fock (RBHF) calculations for finite nuclei or neutron drops (ideal systems composed of a finite number of neutrons and confined within an external field) is reviewed. The guidance and perspectives towards an ab initio covariant density functional theory for nuclear structure derived from the RBHF results are provided.
We propose a new Monte Carlo method called the pinhole trace algorithm for {it ab initio} calculations of the thermodynamics of nuclear systems. For typical simulations of interest, the computational speedup relative to conventional grand-canonical ensemble calculations can be as large as a factor of one thousand. Using a leading-order effective interaction that reproduces the properties of many atomic nuclei and neutron matter to a few percent accuracy, we determine the location of the critical point and the liquid-vapor coexistence line for symmetric nuclear matter with equal numbers of protons and neutrons. We also present the first {it ab initio} study of the density and temperature dependence of nuclear clustering.
The description of nuclei starting from the constituent nucleons and the realistic interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of the nuclear forces, with two-, three- and possibly higher many-nucleon components, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD-employing Hamiltonians constructed within chiral effective field theory. After a brief overview of the field, we focus on ab initio many-body approaches - built upon the No-Core Shell Model - that are capable of simultaneously describing both bound and scattering nuclear states, and present results for resonances in light nuclei, reactions important for astrophysics and fusion research. In particular, we review recent calculations of resonances in the $^6$He halo nucleus, of five- and six-nucleon scattering, and an investigation of the role of chiral three-nucleon interactions in the structure of $^9$Be. Further, we discuss applications to the $^7$Be$(p,gamma)^8$B radiative capture. Finally, we highlight our efforts to describe transfer reactions including the $^3$H$(d,n)^4$He fusion.