Pierre Arthuis

We present an automated, GPU-accelerated framework for many-body perturbation theory (MBPT) calculations of the zero-temperature nuclear equation of state (EOS) based on chiral nucleon-nucleon and three-nucleon (3N) interactions. Automated diagram generation and evaluation enable the computation of all diagrams up to fifth order in the MBPT expansion at the normal-ordered two-body level in infinite matter, with residual three-body contributions explicitly included up to third order. Multi-GPU acceleration of 3N normal ordering, a novel Monte Carlo integrator (called PVegas), and further advances in high-performance computing enable us to evaluate all 840 fifth-order diagrams with controlled numerical uncertainties. We investigate the MBPT convergence up to fifth order in pure neutron matter and symmetric nuclear matter for two sets of chiral interactions, study neutron-star matter, and present fourth-order results for asymmetric matter including normal-ordered 3N forces. The framework enables systematic MBPT studies with harder interactions and benchmarks against nonperturbative methods. It can be further extended to finite-temperature EOS calculations and to improved uncertainty quantification using emulation and resummation techniques.

Bogoliubov Many-Body Perturbation Theory for Nuclei : Systematic Generation and Evaluation of Diagrams and First ab initio Calculations: Les dernières décennies ont donné lieu à un développement rapide des théories ab initio visant à décrire les propriétés des noyaux à partir de l'interaction nucléonique. Un tel développement a été rendu possible à la fois par la très importante croissance de la puissance de calcul et de nouveaux développements formels. Le présent travail se consacre au développement de la théorie de perturbation à N corps de Bogolioubov récemment proposée, qui repose sur l'usage d'un état de référence brisant la symétrie associée au nombre de particules pour permettre une description des noyaux à simple couche ouverte. Le formalisme est tout d'abord décrit en détails, son lien avec la théorie de perturbation à N corps standard est établi, tout comme sa connexion avec la théorie de cluster couplés de Bogolioubov. L'extension du formalisme à des ordres plus élevés à partir de méthodes de théorie des graphes est ensuite présentée ainsi que le programme ADG qui génère et évalue les diagrammes BMBPT à un ordre quelconque. Les implications de ce développement formel dépassent le cadre du présent travail, les méthodes développées pouvant être appliqués à d’autres méthodes à N corps. Pour terminer, de premiers résultats numériques pour les isotopes de l'oxygène, du calcium et du nickel sont présentés. Ces résultats établissent la théorie de perturbation à N corps de Bogolioubov comme une méthode de premier intérêt pour des calculs à grande échelle sur les chaînes isotopiques et isotoniques de masse moyenne.

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

First-principles simulations of many-fermion systems are commonly limited by the computational requirements of processing large data objects. As a remedy, we propose the use of low-rank approximations of three-body interactions, which are the dominant such limitation in nuclear physics. We introduce a randomized decomposition technique to handle the excessively large matrix dimensions and study the sensitivity of low-rank properties to interaction details. The developed low-rank three-nucleon interactions are benchmarked in simulations of few- and many-body systems. Exploiting low-rank properties provides a promising route to extend the microscopic description of atomic nuclei to large systems where storage requirements exceed the computational capacities of the most advanced high-performance computing facilities. Published by the American Physical Society 2024

Neutron-rich nuclei provide important insights to nuclear forces and to the nuclear equation of state. Advances in ab initio methods combined with new opportunities with rare isotope beams enable unique explorations of their properties based on nuclear forces applicable over the entire nuclear chart. In this Letter, we develop novel chiral low-resolution interactions that accurately describe bulk properties from $^{16}$O to $^{208}$Pb. With these, we investigate density distributions and neutron skins of neutron-rich nuclei. Our results show that neutron skins are narrowly predicted over all nuclei with interesting sensitivities for the most extreme, experimentally unexplored cases.

Neutron matter, through its connection to neutron stars as well as systems such as cold atom gases, is one of the most interesting yet computationally accessible systems in nuclear physics. The configuration-interaction Monte Carlo (CIMC) method is a stochastic many-body technique allowing us to tackle strongly coupled systems. In contrast to other quantum Monte Carlo methods employed in nuclear physics, the CIMC method can be formulated directly in momentum space allowing for an efficient use of nonlocal interactions. In this work we extend CIMC method to include three-nucleon interactions through the normal-ordered two-body approximation. We present results for the equation of state of neutron matter in line with other many-body calculations that employ low-resolution chiral interactions, and provide predictions for the momentum distribution and the static structure factor.

Neutron matter, through its connection to neutron stars as well as systems like cold atom gases, is one of the most interesting yet computationally accessible systems in nuclear physics. The Configuration-Interaction Monte Carlo (CIMC) method is a stochastic many-body technique allowing to tackle strongly coupled systems. In contrast to other Quantum Monte Carlo methods employed in nuclear physics, the CIMC method can be formulated directly in momentum space allowing for an efficient use of non-local interactions. In this work we extend CIMC method to include three-nucleon interactions through the normal-ordered two-body approximation. We present results for the equation of state of neutron matter in line with other many-body calculations that employ low resolution chiral interactions, and provide predictions for the momentum distribution and the static structure factor.