A theory is formulated to derive an effective interaction in nuclei in the framework of the unitarymodel-operator approach.A unitary operator is introduced to describe two-body correlations, and it is determined from the condition that the effective interaction should be decoupled between the lowmomentum model space and the high-momentum space.The unitary-transformed Hamiltonian is represented in a cluster expansion form.The one-and two-body cluster terms are mainly discussed.The theory is applied to the calculation of the energies of the ground-and one-body states with the Reid soft core potential.A large gain of the energy has been obtained for the ground state in comparison with the result of the usual G-matrix theory.As for the one-body energies obtained, agreement with the experimental values is fairly good.Although the discrepancies with the experiments still remain at the 1O~ 15 % level, the present approach shows to have some desirable features compared with the nuclear many-body theories proposed so far. § 1. IntroductionA fundamental problem in nuclear physics is to clarify the properties of nuclei in terms of the interaction among their constituent particles.Although there have been many theories or approaches, we have not yet reached the goal of reproducing successfully the experimental data.The investigation of nuclear properties, in the usual nonrelativistic approach, proceeds in two steps: The first one is to determine the interaction among nucleons, and the second is to solve the many-body SchrOdinger equation.Concerning the nucleon-nucleon force, it seems that the phenomenological potential, such as the Reid potential, would be incomplete and one should wait until quark dynamics, probably quantum chromodynamics (QCD), solves finally the problem o~ nuclear force.It is of no doubt that in order to fit the experiments we need more accurate knowledge on the nuclear force including the threebody force.However, if we want to reach the goal, we need also a good many-body theory.Regardless of the accuracy of the nuclear force, we shall not be able to solve the problem unless we have a way of constructing an accurate solution for a many-body system.We believe that many-body theory should share the responsibility for the discrepancy between theory and experiment.It will be sure that nuclear many-body theory is still to be improved.Among many theories for nuclear many-body problem, the most widely applied method is, needless to say, the linked-diagram perturbation theory.This approach may also be referred to as the folded-diagram theory based on the G-matrix interaction especially due 'to Brandow 1l and Kuo et al. 2l An outstanding problem in this approach is the convergence of the perturbation expansion.Since Schucan and Weidenmliller 3 ) pointed out the divergence of the perturbation expansion due to the presence of intruder states, this difficulty has long been a serious problem.Present status of this problem is
The solution of the nuclear A-body problem encounters severe limitations from the size of many-body operators. These limitations are typically related to both the (iterative) storing of the associated tensors and to the computational time related to their multiple contractions in the calculation of various quantities of interest. However, not all the degrees of freedom encapsulated into these tensors equally contribute to the description of many-body observables. Identifying systematic and dominating patterns, a relevant objective is to achieve an \emph{a priori} reduction to the most relevant degrees of freedom via a pre-processing of the A-body problem. The present paper is dedicated to the analysis of two different paradigms to do so. The factorization of tensors in terms of lower-rank ones, whose know-how has been recently transferred to the realm of nuclear structure, is compared to a reduction of the tensors' index size based on an importance truncation. While the objective is to eventually utilize these pre-processing tools in the context of non-perturbative many-body methods, benchmark calculations are presently performed within the frame of perturbation theory. More specifically, we employ the recently introduced Bogoliubov many-body perturbation theory that is systematically applicable to open-shell nuclei displaying strong correlations. This extended perturbation theory serves as a jumpstart for non-perturbative Bogoliubov coupled cluster and Gorkov self-consistent Green's function theories. Results obtained in "small" model spaces are equally encouraging for tensor factorization and importance truncation techniques. While the former requires significant numerical developments to be applied in large model spaces, the latter is presently applied in this context and demonstrates great potential to enable high-accuracy calculations at a much reduced computational cost.
Many-body perturbation theory provides a powerful framework to study the ground state and thermodynamic properties of nuclear matter as well as associated single-particle potentials and response functions within a systematic order-by-order expansion. However, computational challenges can emerge beyond the lowest orders of perturbation theory, especially when computing both single-particle potentials and response functions, which in general are complex-valued and require Cauchy principal value calculations of high-dimensional integrals. We demonstrate that normalizing flows are suitable for Monte Carlo importance sampling of both regular and irregular functions appearing in nuclear many-body calculations. Normalizing flows are a class of machine learning models that can be used to build and sample from complicated distributions through a bijective mapping from a simple base distribution. Furthermore, a well-trained model for a certain target integrand can be efficiently transferred to calculate related integrals with varying physical conditions. These features can enable more efficient tabulations of nuclear physics inputs to numerical simulations of supernovae and neutron star mergers across varying physical conditions and nuclear force models.