Paper detail

Nuclear energy density functional from chiral pion-nucleon dynamics revisited

We use a recently improved density-matrix expansion to calculate the nuclear energy density functional in the framework of in-medium chiral perturbation theory. Our calculation treats systematically the effects from $1π$-exchange, iterated $1π$-exchange, and irreducible $2π$-exchange with intermediate $Δ$-isobar excitations, including Pauli-blocking corrections up to three-loop order. We find that the effective nucleon mass $M^*(ρ)$ entering the energy density functional is identical to the one of Fermi-liquid theory when employing the improved density-matrix expansion. The strength $F_\nabla(ρ)$ of the $(\vec\nabla ρ)^2$ surface-term as provided by the pion-exchange dynamics is in good agreement with that of phenomenological Skyrme forces in the density region $ρ_0/2 <ρ<ρ_0$. The spin-orbit coupling strength $F_{so}(ρ)$ receives contributions from iterated $1π$-exchange (of the ``wrong sign'') and from three-nucleon interactions mediated by $2π$-exchange with virtual $Δ$-excitation (of the ``correct sign''). In the region around $ρ_0/2 \simeq 0.08 $fm$^{-3}$ where the spin-orbit interaction in nuclei gains most of its weight these two components tend to cancel, thus leaving all room for the short-range spin-orbit interaction. The strength function $F_J(ρ)$ multiplying the square of the spin-orbit density comes out much larger than in phenomenological Skyrme forces and it has a pronounced density dependence.