Relativistic Mean Field Approximation

• Relativistic Mean Field (RMF) approximation is a Lorentz-invariant framework that uses self-consistent Dirac equations to incorporate energy-, density-, and isospin-dependent nucleon self-energies.
• The approach maps DBHF-derived self-energy information onto local mean fields with ILDA corrections to accurately predict nuclear binding energies, radii, and spin–orbit splittings.
• Despite its precision, the RMF model underpredicts radii and spin–orbit splittings, highlighting the need for extensions to incorporate three-body forces and long-range correlations.

The relativistic mean field (RMF) approximation is a non-perturbative, covariant framework for modeling the structure and dynamics of nuclear matter and finite nuclei. Grounded in Lorentz-invariant quantum field theory, RMF bridges the gap between free-space nucleon–nucleon (NN) interactions and nuclear observables by embedding dynamical nucleon self-energies, extracted from ab initio treatments such as Dirac–Brueckner–Hartree–Fock (DBHF), into local mean fields. These self-energies are energy-, density-, and isospin-dependent and encode the essential medium modifications of nucleonic degrees of freedom. The approximation yields a set of self-consistent Dirac equations for nucleons coupled to meson mean fields (or, in DBHF-rooted approaches, to parameterized scalar and vector potentials), facilitating an accurate description of nuclear binding, radii, single-particle spectra, and spin–orbit splittings. The RMF scheme is systematized via local density approximations (LDA) and surface corrections (ILDA), with microscopic input from realistic NN interactions and calibration to empirical nuclear data (Sun et al., 2020).

1. Fundamental Structure: DBHF Self-Energy and Dirac Equation

The RMF approach, as implemented in DBHF-inspired models, originates from the in-medium two-body scattering problem using one-boson-exchange NN potentials, e.g., Bonn A. The key step is the extraction of the nucleon self-energy $\Sigma^\tau(E,k,\rho,\beta)$ via a subtracted T-matrix technique. Lorentz covariance yields:

$$\Sigma^\tau(E, k, \rho, \beta) = \Sigma_S^\tau(k, E, \rho, \beta) - \gamma_0\,\Sigma_0^\tau(k, E, \rho, \beta) + \gamma\cdot k\,\Sigma_V^\tau(k, E, \rho, \beta),$$

where $\tau$ denotes isospin (p, n), $E$ is the single-particle energy, $k$ the momentum, $\rho$ the total density, and $\beta$ the asymmetry parameter.

From these, the Schrödinger-equivalent scalar and vector potentials for nucleons emerge as: $$U_S^\tau(E, \rho, \beta) = \frac{\Sigma_S^\tau - M\,\Sigma_V^\tau}{1 + \Sigma_V^\tau}, \quad U_0^\tau(E, \rho, \beta) = \frac{-\Sigma_0^\tau + \epsilon\,\Sigma_V^\tau}{1 + \Sigma_V^\tau},$$ with $\epsilon = E + M$.

The single-particle dynamics are governed by the Dirac equation: $$[\alpha\cdot p + \gamma_0(M + U_S^\tau) + U_0^\tau]\,\Psi^\tau = \epsilon\,\Psi^\tau,$$ which incorporates medium-induced modifications to nucleon masses and energies (Sun et al., 2020).

2. Parameterization of DBHF Self-Energies and Density Functional Mapping

To make the DBHF self-energy accessible for finite systems, its dependence on density, energy, and asymmetry is mapped onto explicit functional forms, e.g., for symmetric matter: $$U_{S(0)}(E, \rho, 0) = (b_{11}\,E + b_{12})\,\rho^{2/3} + (b_{21}\,E + b_{22})\,\rho, \quad U_{0(0)}(E, \rho, 0) = ...$$ Here, $b_{ij}$ are fitted coefficients, with additional high-density corrections for $\rho > 0.16$ fm$^{-3}$.

Isospin asymmetry modifies these forms by an almost linear term in $\beta$: $$U_S^\tau(E, \rho, \beta) = U_S(E, \rho, 0)\,\left[1 + (c_1^\tau E + c_2^\tau)\rho + (c_3^\tau E + c_4^\tau)\right]\beta,$$ with $c_i^\tau$ obtained from DBHF fits (Sun et al., 2020). These parameterizations reproduce scalar and vector potentials to better than a few percent over the relevant density and energy domains.

3. Finite Nucleus Implementation: Local Density and Surface Corrections

In finite nuclei, the nucleon density distributions $\rho_\tau(r)$ and local asymmetry $\beta(r)$ are used to construct spatially dependent mean fields for each single-particle state $i$: $$U_{i,S(0)}(r) = U_{S(0)}^{\tau_i}[E_i, \rho(r), \beta(r)].$$ Self-consistent iteration of the Dirac equation updates the nucleon densities and energies until convergence.

The standard LDA scheme neglects finite-range effects associated with NN interactions, leading to systematic errors in surface properties. To correct this, the improved local density approximation (ILDA) applies Gaussian smearing to the mean fields: $$U_{S(0)}^{\text{ILDA}}(r) = \frac{1}{(t\sqrt{\pi})^3} \int d^3 r'\,U_{S(0)}(r')\,\exp\left[-|r - r'|^2/t^2\right],$$ where $t \sim 0.5$–1.0 fm is a single surface width parameter adjusted to reproduce empirical surface energies. This ILDA approach substantially improves bulk prediction of binding energies, maintaining a microscopic connection to the underlying DBHF effective interaction (Sun et al., 2020).

4. Self-Consistent Observables: Binding Energy, Radii, Spin–Orbit Splitting

With ILDA mean fields, the Dirac equation is solved including the Coulomb potential. The total binding energy is computed as:

$$E = \sum_{i<F} \epsilon_i - \frac{1}{2} \sum_{i<F} \int\left[U_{i,S}(r) \rho_{S,i}(r) + U_{i,0}(r) \rho_i(r)\right] d^3r - \frac{1}{2} \int A_0(r) \rho_C(r) d^3r + E_\text{CM}$$

where $\rho_{S,i} = \bar\Psi_i \Psi_i$ (scalar density), $\rho_i = \Psi_i^\dagger \Psi_i$ (vector density), and $E_\text{CM} \simeq -(3/4)\cdot41\,A^{-1/3}$ MeV corrects for center-of-mass motion.

Charge radii are obtained from the proton density and folded with the proton form factor; spin–orbit splittings arise naturally from the relativistic enhancement of the small Dirac component (Sun et al., 2020).

5. Numerical Systematics: Accuracy and Surface Effects

ILDA-RMF calculations across closed-shell nuclei (from $^{16}$O to $^{208}$Pb) show that:

• LDA (with $t \to 0$) overbinds by several MeV and underestimates rms charge radii by 10–20%, especially in light nuclei.
• ILDA (with fitted $t$) achieves binding energies in excellent agreement with experiment ($\Delta(E/A) \lesssim 0.1$ MeV) and reduces radii errors to 5–10%.
• The difference $(E^\text{ILDA} - E^\text{LDA})/A$ scales as $17.2\,A^{-1/3} - 2.5$ MeV, manifesting standard liquid-drop surface systematics.
• Spin–orbit splittings are reproduced within $\sim$15% but systematically smaller than empirical values, reflecting residual deficiencies in the scalar–vector balance and the absence of core-polarization corrections.

The single surface width parameter $t$ captures the leading effect of NN finite-range physics while maintaining direct connection to DBHF mean-field input (Sun et al., 2020).

6. Limitations and Prospective Extensions

Despite its successes, the DBHF-based RMF approach displays systematic $\sim$10% underprediction of radii and spin–orbit splittings, indicating missing long-range correlations and three-body force effects in the underlying self-energy. Prospective improvements suggested include:

• Extending the subtracted T-matrix framework to account for three-nucleon correlations.
• Incorporating density-functional corrections to mimic beyond-mean-field fluctuations.
• Generalizing the scheme for open-shell and deformed nuclei, including pairing and configuration mixing.

Such developments would enhance the predictive power of the approach and further establish the DBHF-rooted RMF as a comprehensive, microscopic theory of nuclear structure (Sun et al., 2020).

7. Summary and Impact

The relativistic mean field approximation, especially in its DBHF-based incarnation combined with ILDA surface corrections, provides a quantitatively reliable, microscopically motivated framework for the calculation of the bulk properties of finite nuclei. It achieves accuracy in binding energies across the nuclear chart with minimal phenomenological input (a single surface parameter), reproduces key structural observables, and systematically connects empirical nuclear data to realistic NN interactions. Its current limitations reveal the need for higher-order correlations and three-body effects, a frontier for future microscopic nuclear modeling (Sun et al., 2020).