Microscopic nucleus-nucleus optical potentials from nuclear matter
with uncertainty analysis from chiral forces
T. R. Whitehead1
1Facility for Rare Isotope Beams, Michigan State University, East Lansing, Michigan 48824, USA
Nucleus-nucleus optical potentials are constructed from an energy density functional approach
first outlined by Brueckner et al. The interaction term of the energy density functional comes
from the complex nucleon self-energy computed in nuclear matter with two- and three-body chiral
nuclear forces. Nuclear density distributions are calculated from Skyrme functionals constrained to
the equations of state calculated from the same chiral forces used for the self-energy. Predictions for
elastic scattering cross sections and fusion cross sections are compared to experimental data. Very
good agreement is found with experiment for elastic scattering of heavier nucleus-nucleus systems at
energies in the range of 20 < E < 90 MeV/N, while accurate descriptions of lighter and lower-energy
systems may require the inclusion of collective excitations.
Introduction - As the field of low-energy nuclear theory
increasingly turns its focus to exotic nuclei, theoretical
predictions with uncertainty quantification will be essen-
tial for experimental efforts at radioactive beam facilities.
In past studies of nuclei near stability, phenomenological
reaction models tuned to experimental data were consid-
ered sufficient. However, as modern experimental capa-
bilities allow for studies of rare isotopes and as advances
in nuclear many-body theory and our understanding of
the nuclear force enable predictive microscopic calcula-
tions of exotic nuclei, it is imperative that nuclear re-
action models be further developed to accommodate the
needs of experimental efforts.
Optical potentials play a central role in the theoretical
modeling of nuclear reactions and the interpretation of
experimental data. There has been much recent develop-
ment in the area of microscopic nucleon-nucleus optical
potentials [1–5]. However, most reaction experiments in
the rare isotope era will involve nucleus-nucleus interac-
tions and therefore require nucleus-nucleus optical po-
tentials. The study of reaction channels such as knock-
out, transfer, charge-exchange, and fusion could bene-
fit substantially from the implementation of microscopic
nucleus-nucleus optical potentials with quantified uncer-
tainties. In contrast to nucleon-nucleus optical poten-
tials, less attention has been paid to the development
of microscopic nucleus-nucleus optical potentials, with a
few notable exceptions that focus on light systems. These
works follow the double-folding approach [11] where ei-
ther a free-space NN interaction [6–8] or a g-matrix inter-
action [9, 10] is folded with nuclear density distributions.
When a free-space NN interaction is employed, the dou-
ble folding approach yields a real interaction which must
be supplemented by an imaginary term to account for
absorption. Typically this is done by assuming the real
and imaginary radial dependences to be equal and ad-
justing the strength of the imaginary term to reproduce
experimental data. Recently in Ref. [7], a more theo-
retically rigorous approach is taken by utilizing the dis-
persion relation to derive the strength of the imaginary
term from the real part. Another important distinction
of the work by Durant et al. [6–8] amongst the double
folding approaches is the estimation of the theoretical un-
certainty by varying the radial cutoff of a local chiral NN
interaction. Aside from these works, the quantification
of theoretical uncertainties in the calculation of nucleus-
nucleus optical potentials has been largely neglected. In
Refs. [9, 10] a double folding approach is carried out us-
ing the Melbourne g-matrix interaction modified to ap-
proximately account for effects from three-body forces.
By utilizing the g-matrix instead of a free-space interac-
tion, in-medium effects are accounted for. In Ref. [10] the
authors show the importance of including collective ex-
citations in a coupled-channels framework. They achieve
substantially better agreement with data in the coupled-
channels calculations, however discrepancies still persist
at large scattering angles.
An alternative method for constructing the nucleus-
nucleus optical potential is through an energy density
functional based on the nuclear matter single-nucleon po-
tential. Brueckner et al. first outlined how energy density
functionals may be used to construct real nucleus-nucleus
interactions in Refs. [12, 13]. This energy density func-
tional approach may be generalized to include the imag-
inary term of the nucleus-nucleus optical potential on
equal footing to the real term [14–22] in contrast to the
double folding approach when carried out with an NN
interaction. Another advantage of utilizing many-body
calculations in nuclear matter is that Pauli blocking ef-
fects and correlations are included.
The present work constructs a nucleus-nucleus optical
potential starting from many-body perturbation theory
calculations of the self-energy in nuclear matter, which
are also the basis of the WLH global nucleon-nucleus op-
tical potential [1]. The results of the current work may
be combined with the WLH model for a consistent micro-
scopic treatment of the effective interactions in few-body
reaction models. Two-body nuclear forces from chiral ef-
fective field theory (EFT) calculated to N3LO with three-
body forces at N2LO [23–27] are used in the calculation
of the nuclear matter self-energy. To estimate the theo-
retical uncertainty from the nuclear force in predictions
of nuclear reaction observables, three chiral interactions
with momentum space cutoffs in the range of Λ= 414-
arXiv:2210.03031v1 [nucl-th] 6 Oct 2022