Three-body Momentum Representation SRG Evolution and Operator Factorization

E. R. Anderson, R. J. Furnstahl, Robert J. Perry +2 more

Brian Dainton

We first examine how T-matrix equivalence drives the flow of similarity renormalization group (SRG) evolved potential matrix elements to a universal form, with the ultimate goal of gaining insight into universality for three-nucleon forces. In agreement with observations made previously for Lee-Suzuki transformations, regions of universal potential matrix elements are restricted to where half-on-shell T-matrix equivalence holds, but the potentials must also reproduce binding energies. We find universality in local energy regions, reflecting a local decoupling by the SRG. To continue the study in the 3-body sector, we create a simple 1-D spinless boson "theoretical laboratory" for a dramatic improvement in computational efficiency. We introduce a basis-transformation, harmonic oscillator (HO) basis, which is used for current many-body calculations and discuss the imposed truncations. When SRG evolving in a HO-basis, we show that the evolved matrix elements, once transformed back into momentum-representation, differ from those when evolving with momentum representation. This is because the generator in each basis is not exactly the same due to the truncation, which creates evolution artifacts in the 3-body potential matrix elements. In the 2- body sector, this can be avoided by increasing the basis size, but it remains unclear whether this is possible in the 3-body sector, as truncation errors in the 3-body sector are more difficult to avoid, and the computational power required is greatly increased for three-body evolution.

Bin Zhu, R. Wirth, H. Hergert

One of the main challenges for ab initio nuclear many-body theory in the coming decade is the growth of computational and storage costs as calculations are extended to increasingly heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform singular value decompositions (SVD) of current nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization schemes and resolution scales. Next, we develop and implement the similarity renormalization group (SRG) evolution of the interaction in terms of the relevant singular vectors, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions naturally allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Finally, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the in-medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.

Micah D. Schuster, Sofia Quaglioni, Calvin M. Johnson +2 more

Kyle Wendt, R. J. Furnstahl, Robert J. Perry

B. D. Carlsson

The Similarity Renormalization Group (SRG) flow equation is explored for systems of two and three spinless bosons in a momentum-space partial-wave basis. The two- and three-body binding energies as well as the phaseshifts are used to gauge that the transformation is unitary and to study how well the SRG decouples high- and low-energy physics. I consider four different potentials with different characteristics: Two simplied nucleon potentials and two inter-atomic helium potentials (a soft-core potential and the state-of-the-art LM2M2 potential that is fitted to a wealth of experimental data). An initial three-body force is included for two of these potentials. Even with only two-body terms in the initial hamiltonian, SRG induced many-body forces are shown to arise during the transformation. These induced forces are computed for the three-body system and their evolution is studied as a function of the flow parameter. In all cases the SRG transformed potentials display greatly improved decoupling. This is achieved with a three-body binding-energy deviation of less than 0.1% in all cases except for the soft-core helium potential.\n