Abstract
Three-nucleon (3N) interactions are key for an accurate solution of the nuclear many-body problem. However, fully taking into account 3N forces constitutes a computational challenge and hence approximate treatments are commonly employed. The method of normal ordering has proven to be a powerful tool that allows to systematically include 3N interactions in an efficient way, but traditional normal-ordering frameworks require the representation of 3N interactions in a large single-particle basis, typically necessitating a truncation of 3N matrix elements. While this truncation has only a minor impact for light and medium-mass nuclei, its effects become sizable for heavier systems and hence limit the scope of ab initio calculations. In this work, we present a novel normal-ordering framework that allows to circumvent this limitation by performing the normal ordering directly in a Jacobi basis. We discuss in detail the new framework, benchmark it against established results, and present calculations for ground-state energies and charge radii of heavy nuclei, such as Sn132 and Pb208.
| Original language | English |
|---|---|
| Article number | 024310 |
| Journal | Physical Review C |
| Volume | 107 |
| Issue number | 2 |
| DOIs | |
| State | Published - Feb 2023 |
| Externally published | Yes |
Funding
We thank P. Arthuis and T. Miyagi for helpful discussions. We also thank T. Miyagi for providing matrix elements for our benchmark calculations. This work was supported in part by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant Agreement No. 101020842), the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)–Projektnummer 279384907–SFB 1245 and Germany's Excellence Strategy–EXC 2118 –390831469, and by the BMBF Contract No. 05P21RDFNB. Computations were in part performed with an allocation of computing resources at the Jülich Supercomputing Center. This work was completed in part at the Helmholtz GPU Hackathon 2021, part of the Open Hackathons program. The authors would like to acknowledge OpenACC-Standard.org for their support.