Nuclear energy density functional from chiral pion-nucleon dynamics

We calculate the nuclear energy density functional relevant for N=Z even-even nuclei in the systematic framework of chiral perturbation theory. The calculation includes the one-pion exchange Fock diagram and the iterated one-pion exchange Hartree and Fock diagrams. From these few leading order contributions in the small momentum expansion one obtains already a very good equation of state of isospin symmetric nuclear matter. We find that in the region below nuclear matter saturation density the effective nucleon mass $\widetilde M^*(\rho)$ deviates by at most 15% from its free space value $M$, with $0.89M<\widetilde M^*(\rho)<M$ for $\rho < 0.11 {\rm fm}^{-3}$ and $\widetilde M^*(\rho)>M$ for higher densities. The parameterfree strength of the $(\vec\nabla \rho)^2$-term, $F_\nabla(k_f)$, is at saturation density comparable to that of phenomenological Skyrme forces. The magnitude of $F_J(k_f)$ accompanying the squared spin-orbit density $\vec J ^2$ comes out somewhat larger. The strength of the nuclear spin-orbit interaction, $F_{so}(k_f)$, as given by iterated one-pion exchange is about half as large as the corresponding empirical value, however, with the wrong negative sign. The novel density dependencies of $\widetilde M^*(\rho)$ and $F_{\nabla,so,J}(k_f)$ as predicted by our parameterfree calculation should be examined in nuclear structure calculations (after introducing an additional short range spin-orbit contribution constant in density).Comment: 16 pages, 5 figure

Nuclear energy density functional from chiral pion-nucleon dynamics: Isovector spin-orbit terms

We extend a recent calculation of the nuclear energy density functional in the systematic framework of chiral perturbation theory by computing the isovector spin-orbit terms: $(\vec \nabla \rho_p- \vec \nabla \rho_n)\cdot(\vec J_p-\vec J_n) G_{so}(k_f)+ (\vec J_p-\vec J_n)^2 G_J(k_f)$. The calculation includes the one-pion exchange Fock diagram and the iterated one-pion exchange Hartree and Fock diagrams. From these few leading order contributions in the small momentum expansion one obtains already a good equation of state of isospin-symmetric nuclear matter. We find that the parameterfree results for the (density-dependent) strength functions $G_{so}(k_f)$ and $G_J(k_f)$ agree fairly well with that of phenomenological Skyrme forces for densities $\rho > \rho_0/10$. At very low densities a strong variation of the strength functions $G_{so}(k_f)$ and $G_J(k_f)$ with density sets in. This has to do with chiral singularities $m_\pi^{-1}$ and the presence of two competing small mass scales $k_f$ and $m_\pi$. The novel density dependencies of $G_{so}(k_f)$ and $G_J(k_f)$ as predicted by our parameterfree (leading order) calculation should be examined in nuclear structure calculations.Comment: 9 pages, 3 figure, published in: Physical Review C68, 014323 (2003

Chiral approach to nuclear matter: Role of two-pion exchange with virtual delta-isobar excitation

We extend a recent three-loop calculation of nuclear matter in chiral perturbation theory by including the effects from two-pion exchange with single and double virtual $\Delta(1232)$-isobar excitation. Regularization dependent short-range contributions from pion-loops are encoded in a few NN-contact coupling constants. The empirical saturation point of isospin-symmetric nuclear matter, $\bar E_0 = -16$MeV, $\rho_0 = 0.16$fm$^{-3}$, can be well reproduced by adjusting the strength of a two-body term linear in density (and weakening an emerging three-body term quadratic in density). The nuclear matter compressibility comes out as $K = 304$MeV. The real single-particle potential $U(p,k_{f0})$ is substantially improved by the inclusion of the chiral $\pi N\Delta$-dynamics: it grows now monotonically with the nucleon momentum $p$. The effective nucleon mass at the Fermi surface takes on a realistic value of $M^*(k_{f0})=0.88M$. As a consequence of these features, the critical temperature of the liquid-gas phase transition gets lowered to the value $T_c \simeq 15$MeV. In this work we continue the complex-valued single-particle potential $U(p,k_f)+i W(p,k_f)$ into the region above the Fermi surface $p>k_f$. The effects of $2\pi$-exchange with virtual $\Delta$-excitation on the nuclear energy density functional are also investigated. The effective nucleon mass associated with the kinetic energy density is $\widetilde M^*(\rho_0)= 0.64M$. Furthermore, we find that the isospin properties of nuclear matter get significantly improved by including the chiral $\pi N\Delta$-dynamics.Comment: 28 pages, 13 figure