Covariant Energy Density Functionals

• CEDFs are covariant energy density functionals that implement relativistic density functional theory by modeling nucleons as Dirac particles in self-consistent meson or point-coupling fields.
• They are calibrated against nuclear observables such as binding energies, charge radii, and deformations, with formulations including NL, DD-ME, and point-coupling models.
• Recent developments focus on improving calibration methods, incorporating beyond-mean-field correlations, and optimizing numerical techniques to reduce uncertainties.

Searching arXiv for the cited CEDF papers to ground the article in the specified literature. Covariant energy density functionals (CEDFs) are the functional realizations of covariant density functional theory (CDFT), a relativistically consistent framework for the microscopic description of bulk and spectral nuclear properties. In this approach, nuclei are described by Dirac nucleons moving in self-consistent scalar and vector mean fields generated either by effective meson exchange or by point-coupling interactions. CEDFs are used for the global description of binding energies, charge radii, deformations, neutron skins, drip lines, fission barriers, and dense-matter equations of state, and they have been benchmarked across large portions of the nuclear chart with explicit assessments of systematic and statistical uncertainties (Afanasjev, 2015, Afanasjev et al., 2021).

1. Covariant formulation and mean-field structure

The defining feature of a CEDF is Lorentz covariance. In the meson-exchange formulation, the basic Lagrangian density couples the nucleon Dirac field $\psi$ to the isoscalar-scalar $\sigma$, isoscalar-vector $\omega_\mu$, isovector-vector $\rho_\mu$, and electromagnetic $A_\mu$ fields. A representative form is

$$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$

with the corresponding meson and photon kinetic and mass terms. In nonlinear meson models such as NL3*, the density dependence is encoded through the $\sigma$ self-interaction

$$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$

whereas in density-dependent meson-exchange models and in density-dependent point-coupling models the couplings are explicit functions of the baryon or vector density, so that no explicit $U(\sigma)$ is required (Afanasjev, 2015, Afanasjev et al., 2016).

Variation of the functional yields the Dirac equation for single-particle spinors,

$$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$

with scalar and vector self-energies generated by the covariant fields. In the Hartree approximation, the scalar field is typically $\sigma$0, while the time-like vector field contains $\sigma$1, $\sigma$2, and Coulomb contributions. The total energy functional in the static mean-field approximation is built from the occupied Dirac orbitals plus meson-field, Coulomb, and, when the couplings are density-dependent, rearrangement contributions (Afanasjev, 2015, Afanasjev et al., 2021).

Point-coupling CEDFs eliminate explicit mesons in favor of local contact interactions and gradient terms. A representative interaction part is

$\sigma$3

This representation is central to DD-PC1, PC-PK1, and later neural-network or reduced-basis developments built on point-coupling RMF structure (Afanasjev et al., 2021, Li et al., 29 Jun 2026, Anderson et al., 2024).

Pairing correlations are usually treated within the relativistic Hartree-Bogoliubov (RHB) framework. In the global surveys and many later developments, the particle-hole channel is defined by the CEDF and the particle-particle channel by a finite-range Gogny-type force or a finite-range separable pairing force tuned to reproduce pairing gaps (Afanasjev et al., 2021, Agbemava et al., 2014).

2. Major families and representative parametrizations

Three major families of CEDFs are in widespread use: nonlinear meson models (NL), density-dependent meson-exchange models (DD-ME), and point-coupling models (PC). These families differ mainly in how they encode density dependence and finite-range effects, but all retain the covariant scalar-vector structure of the mean field (Afanasjev et al., 2021).

Family	Representative functionals	Characteristic features
NL	NL3*, NL5(C,D,E)	Nonlinear $\sigma$4 self-interactions
DD-ME	DD-ME2, DD-ME$\sigma$5	Density-dependent meson-nucleon couplings
PC	DD-PC1, PC-PK1	Zero-range contact terms plus gradients

NL3* is a nonlinear meson-coupling functional fitted to bulk properties of 12 spherical nuclei and selected nuclear-matter properties; in another formulation it is described as fitted to binding energies, charge radii, neutron skins of selected nuclei, and nuclear incompressibility. DD-ME2 and DD-ME$\sigma$6 are density-dependent meson-exchange functionals calibrated to binding energies and charge radii, with DD-ME$\sigma$7 also incorporating pairing information and an isovector-scalar $\sigma$8 channel. DD-PC1 and PC-PK1 are point-coupling functionals fitted, respectively, to deformed and spherical nuclear data sets, with DD-PC1 emphasizing binding energies and deformations of 64 nuclei and PC-PK1 including charge radii and spin-orbit splittings (Afanasjev, 2015, Afanasjev et al., 2021, Afanasjev et al., 2016).

The density dependence of DD-ME2 is commonly written as

$\sigma$9

with typical values $\omega_\mu$0, $\omega_\mu$1, and $\omega_\mu$2 in one summary and $\omega_\mu$3 in another parameter listing tied to a different source presentation. DD-PC1 is characterized by density-dependent four-fermion couplings in the scalar, vector, and isovector-vector channels and associated gradient terms (Afanasjev et al., 2021, Afanasjev et al., 2016, Taninah et al., 2020).

A standard point of comparison among major CEDFs is the set of nuclear-matter saturation properties. One frequently quoted compilation gives, for NL3*, DD-ME2, DD-ME$\omega_\mu$4, DD-PC1, and PC-PK1, respectively:

• $\omega_\mu$5,
• $\omega_\mu$6,
• $\omega_\mu$7,
• $\omega_\mu$8,
• $\omega_\mu$9,
• $\rho_\mu$0. In that assessment, NL3* and PC-PK1 violate the empirical “SET2b” constraints in $\rho_\mu$1 and $\rho_\mu$2 (Afanasjev, 2015).

3. Calibration targets, observables, and uncertainty measures

CEDFs are calibrated to selected combinations of finite-nucleus observables and nuclear-matter properties. The data entering the fit may include binding energies, charge radii, neutron skins, deformations, spin-orbit splittings, pairing information, and pseudo-data for symmetric matter, pure neutron matter, or other infinite-matter constraints. This heterogeneous calibration practice is central to later discussions of why similarly performing functionals can have different nuclear-matter characteristics, and conversely why functionals close to empirical nuclear-matter windows need not perform best for finite nuclei (Afanasjev, 2015, Afanasjev et al., 2016).

For ground-state studies, the standard observables are the binding energy $\rho_\mu$3, the charge radius

$\rho_\mu$4

and the quadrupole deformation

$\rho_\mu$5

For charge-radii systematics one often starts from the point-proton radius and includes finite-size and Darwin-Foldy corrections,

$\rho_\mu$6

with the compact approximation

$\rho_\mu$7

Differential radii along isotopic chains are defined by

$\rho_\mu$8

and odd-even staggering in charge radii by the three-point indicator

$\rho_\mu$9

These definitions are central in the global charge-radius literature within CDFT (Afanasjev et al., 2021).

For mass systematics, the global rms deviation for a set of $A_\mu$0 even-even nuclei is

$A_\mu$1

Theoretical systematic uncertainty in the “global-spread” approach is quantified at each $A_\mu$2 by

$A_\mu$3

that is, the spread among representative functionals. Statistical uncertainty within a single parametrization is propagated from the covariance matrix,

$A_\mu$4

The same distinction between systematic spreads and covariance-based statistical errors is used in charge-radius assessments, where the spread over NL3*, DD-ME2, DD-ME$A_\mu$5, and DD-PC1 is taken as a systematic error estimate (Afanasjev, 2015, Afanasjev et al., 2021).

4. Global performance across the nuclear chart

Large-scale axial RHB calculations with modern CEDFs have established a consistent but not uniform level of accuracy across the nuclear landscape. In one broad survey of even-even nuclei with $A_\mu$6 between the two-proton and two-neutron drip lines, the rms mass deviations across 835 nuclei were reported as $A_\mu$7 MeV for NL3*, $A_\mu$8 MeV for DD-ME2, $A_\mu$9 MeV for DD-ME$$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$0, and $$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$1 MeV for DD-PC1 (Agbemava et al., 2014). In a closely related later assessment over approximately 640 even-even nuclei, the ordering was DD-PC1 ($$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$2 MeV), DD-ME$$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$3 ($$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$4 MeV), DD-ME2 ($$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$5 MeV), NL3* ($$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$6 MeV), with PC-PK1 quoted at $$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$7 MeV in the saturation-property comparison (Afanasjev, 2015, Afanasjev et al., 2016).

Systematic spreads in binding energies are modest near stability and large at the neutron-rich edge. In the 2015 assessment, $$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$8 is about $$\mathcal{L} = \bar\psi \Bigl[ \gamma^\mu \bigl( i\partial_\mu - g_\omega \omega_\mu - g_\rho \vec\tau\!\cdot\!\vec\rho_\mu - e A_\mu (1-\tau_3)/2 \bigr) - (m+g_\sigma \sigma) \Bigr]\psi +\cdots$$9 MeV in the valley of stability and can grow up to about $\sigma$0 MeV toward the neutron drip line (Afanasjev, 2015). In the 2014 global survey, mass spreads increase from roughly $\sigma$1–$\sigma$2 MeV near stability to more than $\sigma$3 MeV near the two-neutron drip line (Agbemava et al., 2014). These results establish that extrapolation uncertainty is not a peripheral issue but a structural feature of present-day CEDFs.

Charge radii are described more accurately and with smaller functional spreads than masses. Across the nuclear chart, NL3*, DD-ME2, DD-ME$\sigma$4, DD-PC1, PC-PK1, and NL5(C–E) achieve

$\sigma$5

with a mean of about $\sigma$6 fm, corresponding to about $\sigma$7 on a typical $\sigma$8 fm charge radius. Systematic spreads in predicted absolute radii are $\sigma$9 fm for medium-heavy and heavy nuclei and grow to about $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$0 fm for $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$1 (Afanasjev et al., 2021). In earlier global surveys, rms deviations for 351 measured even-even charge radii were $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$2 fm for NL3*, $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$3 fm for DD-ME2, $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$4 fm for DD-ME$$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$5, and $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$6 fm for DD-PC1, with improved values of $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$7–$$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$8 fm if light He and anomalous Cm isotopes were excluded (Agbemava et al., 2014).

Two-particle separation energies, drip lines, and deformation maps reveal where the functional dependence is most consequential. Reported rms deviations are about $$U(\sigma)=\frac12 m_\sigma^2\sigma^2+\frac13 g_2\sigma^3+\frac14 g_3\sigma^4,$$9–$U(\sigma)$0 MeV for $U(\sigma)$1 and about $U(\sigma)$2–$U(\sigma)$3 MeV for $U(\sigma)$4, with the proton side generally better reproduced than the neutron side (Agbemava et al., 2014). The two-proton drip line predicted by DD-ME2 and DD-ME$U(\sigma)$5 agrees with experiment within two neutrons for $U(\sigma)$6, whereas the two-neutron drip line exhibits much larger model spreads, often by tens of neutrons between major shell closures. Shell gaps at $U(\sigma)$7, $U(\sigma)$8, and $U(\sigma)$9 remain robust in all EDFs examined in that survey (Agbemava et al., 2014).

For deformations, the spreads $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$0 are $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$1 in spherical and well-deformed regions but can exceed $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$2 at shape-coexistence boundaries (Agbemava et al., 2014). In actinides and superheavy nuclei, systematic studies with DD-PC1, DD-ME2, NL3*, and PC-PK1 identify spherical shell closures at $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$3, $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$4, and $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$5 as major sources of uncertainty in ground-state deformations and fission barriers. Theoretical uncertainties in barrier heights are moderate in well-deformed actinides but can peak around $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$6 MeV at $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$7, $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$8, and along $$\Bigl[ -\,i\boldsymbol\alpha\!\cdot\!\nabla + \beta\,(m+S(\mathbf r)) + V(\mathbf r) \Bigr]\psi_k(\mathbf r) = \varepsilon_k\psi_k(\mathbf r),$$9; near $\sigma$00, NL3* and PC-PK1 give barriers below $\sigma$01 MeV while DD-ME2 and DD-PC1 stay above $\sigma$02 MeV (Taninah et al., 2020).

5. Correlations, shell structure, and beyond-mean-field physics

A central result of the CEDF literature is that there is no strong one-to-one correlation between global finite-nucleus performance and any single nuclear-matter parameter such as $\sigma$03, $\sigma$04, or $\sigma$05. A representative example compares DD-ME2 and DD-PC1: despite similar $\sigma$06 values ($\sigma$07 versus $\sigma$08 MeV), their mass predictions differ by up to about $\sigma$09 MeV over large portions of the chart. Conversely, DD-ME2 and DD-ME$\sigma$10 have nearly identical $\sigma$11 values and display $\sigma$12 MeV for about half of all nuclei up to $\sigma$13. The conclusion drawn in the global assessments is that differences in shell structure and in the finite-nucleus fit play a role as large as bulk nuclear-matter parameters (Afanasjev, 2015, Afanasjev et al., 2016).

This conclusion directly addresses a recurrent misconception: strict enforcement of empirical nuclear-matter windows is neither necessary nor sufficient for optimal finite-nucleus performance. The 2016 analysis states that functionals that come close to satisfying current nuclear-matter constraints can have problems in the description of existing nuclear data, whereas functionals carefully fitted to finite nuclei but violating some nuclear-matter constraints may perform better for binding energies and other ground- and excited-state properties. The stated reason is that finite nuclei depend not only on bulk nuclear matter but also on shell effects and surface properties (Afanasjev et al., 2016).

Charge-radius systematics sharpen this point further. CEDF calculations reproduce absolute and differential charge radii globally, but outliers appear in regions of shape coexistence and in light nuclei with soft surfaces, where static mean field neglects beyond-mean-field correlations (Afanasjev et al., 2021). In deformed actinides and light superheavy nuclei, inclusion of octupole deformation increases $\sigma$14 by about $\sigma$15 in low-$\sigma$16 Th and U isotopes, improving agreement with laser-spectroscopy data (Afanasjev et al., 2021). For odd-even staggering, recent CEDF studies identify two additional mechanisms beyond the traditional DFT picture: level-ordering changes induced by particle-vibration coupling (PVC), and fragmentation of the blocked state. In the cited schematic model, replacing

$\sigma$17

boosts $\sigma$18 by roughly $\sigma$19, bringing the predicted staggering close to experiment for $\sigma$20–$\sigma$21 (Afanasjev et al., 2021).

Self-consistency is also essential for differential radii. A specific example from the Pb isotopes shows that filling the $\sigma$22 subshell from $\sigma$23Pb to $\sigma$24Pb increases the rms radius of the $\sigma$25 proton orbitals by about $\sigma$26 fm, but only by about $\sigma$27 fm for $\sigma$28 orbitals. This orbital-dependent polarization cannot be captured with a frozen proton core or a uniform radius formula $\sigma$29, and it is identified as crucial for reproducing the magnitude of radius kinks at shell closures (Afanasjev et al., 2021).

Beyond-mean-field quadrupole dynamics substantially improves mass systematics. A global 5DCH study with PC-PK1 for 575 even-even nuclei found dynamic correlation energies associated with rotational and quadrupole vibrational motion ranging from about $\sigma$30 to $\sigma$31 MeV, mostly between $\sigma$32 and $\sigma$33 MeV. After including these correlations, the rms mass deviation was reduced to $\sigma$34 MeV, smaller than the mean-field values quoted for NL3* ($\sigma$35 MeV), DD-ME2 ($\sigma$36 MeV), DD-ME$\sigma$37 ($\sigma$38 MeV), and DD-PC1 ($\sigma$39 MeV); the rms deviation for two-nucleon separation energies was reduced by about $\sigma$40 relative to a cranking prescription (Lu et al., 2015).

Functional extensions at the Hartree-Bogoliubov level have followed the same logic. A 2022 study including tensor terms in the vector-isoscalar channel reported improvements in RMS binding energies, spin-orbit splittings, and shell gaps across the chart, including deformed nuclei. In infinite matter, the Dirac mass increased from $\sigma$41 for DD-ME2 to $\sigma$42 with tensor terms, while the binding-energy RMS over 828 nuclei improved from about $\sigma$43 MeV for DD-ME2 to $\sigma$44 MeV without tensor and $\sigma$45 MeV with tensor; the spin-orbit RMS changed from $\sigma$46 MeV to $\sigma$47 MeV and then $\sigma$48 MeV, and the shell-gap RMS from $\sigma$49 MeV to $\sigma$50 MeV and then $\sigma$51 MeV (Mercier et al., 2022).

6. Numerical optimization, emulation, astrophysical refinement, and emerging directions

The post-2020 CEDF literature has focused heavily on calibration methodology and numerical control. A major advance is the anchor-based optimization approach for global fits to approximately 2000 even-even masses, in which the objective function augments the usual least-squares term with quadratic anchor penalties on selected parameters: $\sigma$52 With a derivative-based minimizer and anchors taken from established functionals or microscopic expectations, the number of full mass-table evaluations is reduced from about $\sigma$53 to about $\sigma$54, so the total cost becomes roughly $\sigma$55 a single global RHB run (Afanasjev et al., 1 Nov 2025).

These optimization protocols are closely tied to the control of basis-truncation and atomic-to-nuclear mass-conversion errors. An empirical infinite-basis extrapolation,

$\sigma$56

using energies computed at three values of $\sigma$57, reduces the fermionic basis-truncation error from about $\sigma$58 MeV at $\sigma$59 to $\sigma$60 keV across the known chart (Afanasjev et al., 1 Nov 2025). A subsequent study incorporated both fermionic and bosonic infinite-basis corrections and total electron binding energies in the fit protocol. It defined pseudodata as

$\sigma$61

and found that neglect of these corrections had induced a global calculation error of order $\sigma$62 MeV or higher for the three major classes of CEDFs. For the new Z-type functionals, the reported rms deviations for 882 even-even nuclei are $\sigma$63 MeV for NL5(Z), $\sigma$64 MeV for PC-Z, and $\sigma$65 MeV for DD-MEZ, with negligible numerical errors of order $\sigma$66–$\sigma$67 keV; DD-MEZ reaches $\sigma$68 MeV in a complementary summary (Osei et al., 22 Jul 2025, Afanasjev et al., 1 Nov 2025).

Reduced-order and optimized-basis methods address the same bottleneck from a different angle. A universal reduced basis for the single-particle Dirac equation, built by proper orthogonal decomposition of high-fidelity snapshots, reproduces bound levels with energy errors $\sigma$69 in dimensionless units using reduced bases of size about $\sigma$70–$\sigma$71 per $\sigma$72 channel. Direct reduced-basis diagonalization yields a speed-up of order $\sigma$73 relative to full Runge-Kutta integration, and linearization of Woods-Saxon-like potentials pushes the speed-up to order $\sigma$74 (Anderson et al., 2024). In point-coupling calculations with harmonic-oscillator bases, a globally optimized frequency scaling

$\sigma$75

combined with empirical formulas for the minimal shell cutoff, reduces the necessary fermionic space by about $\sigma$76–$\sigma$77 shells for $\sigma$78 keV accuracy and cuts CPU times by about $\sigma$79–$\sigma$80 (Dalbah et al., 29 May 2026).

Machine-learning extensions have begun to treat the energy density itself as the learning target. A 2026 point-coupling study used a physics-guaranteed neural network whose inputs at each radial point are

$\sigma$81

and whose output is the local energy density. When trained on a correction to an existing CEDF, with

$\sigma$82

the approach improves the binding-energy accuracy from $\sigma$83 keV to $\sigma$84 keV in the known region and retains extrapolation accuracy of about $\sigma$85 MeV up to $\sigma$86 steps beyond known nuclei. The study reports that shell effects in two-neutron separation energies are effectively captured (Li et al., 29 Jun 2026).

CEDFs also continue to be refined from the dense-matter side. Bayesian refinement of FSUGold2 and FSUGarnet with pulsar masses, NICER radii, tidal deformabilities, and $\sigma$87EFT neutron-matter constraints leads to posterior values such as $\sigma$88 MeV for refined FSUGold2 and $\sigma$89 MeV for refined FSUGarnet, with corresponding maximum neutron-star masses $\sigma$90 and $\sigma$91 (Salinas et al., 2023). At the same time, the PREX/CREX tension remains unresolved within standard isovector sectors; the cited work states that CEDFs with only two isovector couplings lie on a nearly linear $\sigma$92–$\sigma$93 correlation and cannot simultaneously reproduce the reported skins of both nuclei (Salinas et al., 2023). A later Bayesian study of density-dependent coupling parametrizations finds that current modeling requires freedom in the isoscalar channel at least up to the skewness coefficient $\sigma$94 and in the isovector channel at least up to the curvature coefficient $\sigma$95 in order to capture supra-saturation variations in the equation of state and symmetry energy (Wei et al., 1 Dec 2025).

Taken together, these developments suggest a convergent picture rather than a settled one. Present CEDFs provide global mean-field mass accuracies at the $\sigma$96–$\sigma$97 MeV level, charge-radius accuracies at the few-$\sigma$98 fm level, and controlled uncertainty estimates across the chart, but the isovector channel, shell-structure deficiencies, and the consistent inclusion of beyond-mean-field correlations remain the principal open problems. The explicit research directions identified in recent work are broader calibration datasets, quantified covariance analyses, explicit beyond-mean-field treatments such as 5DCH or GCM during the fit, improved spectroscopic quality through tensor or $\sigma$99-meson channels, and more flexible density dependencies or higher-order terms that better decouple isoscalar, isovector, and shell-structure aspects of the functional (Afanasjev, 2015, Afanasjev et al., 1 Nov 2025)