\texttt{NuHamil}: A numerical code to generate nuclear two- and three-body matrix elements from chiral effective field theory

Takayuki Miyagi
DOI: https://doi.org/10.1140/epja/s10050-023-01039-y
2023-07-09
Abstract:The applicability of nuclear {\it ab initio} calculations has rapidly extended over the past decades. However, starting research projects is still challenging due to the required numerical expertise in the generation of underlying nuclear interaction matrix elements and many-body calculations. To ease the first issue, in this paper we introduce the numerical code \texttt{NuHamil} to generate the nucleon-nucleon (NN) and three-nucleon (3N) matrix elements expressed in a spherical harmonic-oscillator basis, inputs of many-body calculations. The ground-state energies for the selected doubly closed shell nuclei are calculated with the no-core shell-model (NCSM) and in-medium similarity renormalization group (IMSRG). The code is written in modern Fortran, and OpenMP+MPI hybrid parallelization is available for the 3N matrix-element calculations.
Nuclear Theory
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the development of numerical codes for generating nucleon - nucleon (NN) and three - nucleon (3N) matrix elements in nuclear physics. Specifically, these matrix elements are the basic inputs for nuclear ab initio calculations, but their generation process is complex and requires specialized knowledge of numerical techniques, which has become an obstacle to new research projects. To solve this problem, the author introduced a numerical code named NuHamil, which can generate NN and 3N matrix elements expressed in the single - particle harmonic - oscillator (HO) basis, and these matrix elements can be used as inputs for most ab initio calculation methods. By providing such a tool, NuHamil aims to simplify the threshold for researchers to enter this field and reduce the technical challenges they face in preparing the necessary matrix element inputs. The paper also discusses how to use the free - space similarity renormalization group (SRG) to soften nuclear interactions in order to accelerate convergence in many - body calculations. In addition, the paper presents some benchmark many - body calculation results using the no - core shell model (NCSM) and the in - medium similarity renormalization group (IMSRG), demonstrating the effectiveness and accuracy of the NuHamil code. Finally, the paper provides instructions for using the code and conclusions, emphasizing the importance and potential applications of NuHamil in nuclear physics research.