Ab Initio Gamow In-Medium Similarity Renormalization Group for Open-Shell Nuclei

Qi Yuan,Baishan Hu,Jianguo Li,Siqin Fan,Furong Xu
DOI: https://doi.org/10.1360/tb-2023-0163
2023-01-01
Abstract:Nuclei near the dripline are open quantum systems, which have strong coupling to the particle continuum. Many novel phenomena have been observed in these nuclei, such as halos, two proton/neutron emission, and the change of magic numbers. It is the frontier of current nuclear physics theory to study the properties of these weakly bound and unbound nuclei. Over the past few decades, ab initio calculations based on realistic nuclear forces derived from chiral effective field theory have made great progress and can now provide good descriptions of light- to heavy-mass nuclei. However, most ab initio methods are performed on the harmonic-oscillator (HO) basis or real-energy Hartree-Fock (HF) basis. The HO potential is an infinitely deep potential well. Hence, the HO basis is bound and localized, and not able to describe the coupling to the particle continuum, which plays an important role in the descriptions of weakly bound and unbound nuclei. The complex-energy Gamow-Berggren basis treats bound, resonant, and scattering continuum states on an equal footing, providing an efficient framework to calculate weakly bound and unbound nuclei. As one of the powerful ab initio methods, the in-medium similarity renormalization group method (IMSRG) diagonalizes the nuclear Hamiltonian by a term of continuous unitary transformations, providing an efficient tool to treat energy and other observables equally. Its extension to the Gamow-Berggren basis, named the Gamow in-medium similarity renormalization group (G-IMSRG), diagonalizes the continuum-coupled Hamiltonian by a term of orthogonal transformations and can well give the properties of weakly bound and unbound nuclei. However, it is currently limited to calculating closed (sub-) shell nuclei due to the convergence problem of the flow equation in complex-energy space. In the present paper, we developed the particle attached or removed equation of motion method to extend the G-IMSRG to the calculations of open-quantum systems with one particle attached to or removed from a closed (sub-) shell nucleus. As test grounds, we calculated the ground-state energies and the excitation energies of low-lying states of the neutron-rich O isotopes based on the chiral two-body and three-body nuclear force NNLOsat. Our calculations show that the three-body force plays a key role to reproduce the dripline of O isotopes and the coupling to the continuum also has important effects on the binding energies and spectra of neutron-rich O isotopes. Overall, the extended G-IMSRG can well describe the properties of weakly bound and unbound nuclei, and it can give the corresponding resonant width.
What problem does this paper attempt to address?