---
title: Coherent Density Fluctuation Model (CDFM)
url: https://www.emergentmind.com/topics/coherent-density-fluctuation-model-cdfm
type: topic
---

# Coherent Density Fluctuation Model (CDFM)

Searching arXiv for recent and foundational papers on the Coherent Density Fluctuation Model (CDFM).
{"query":"Coherent Density Fluctuation Model CDFM nuclear symmetry energy finite nuclei arXiv", "max_results": 10, "sort_by": "relevance"}
{"query":"all:\"coherent density fluctuation model\" OR title:\"coherent density fluctuation model\"", "max_results": 10, "sort_by": "relevance"}
The Coherent Density Fluctuation Model (CDFM) is a correlation-based framework in which a finite nucleus is represented as a coherent superposition of uniform-density spherical pieces of nuclear matter, usually termed fluctons. In the $\delta$-function limit of the generator coordinate method, the model rewrites the one-body density matrix, the local density, and density-dependent observables of a finite nucleus as weighted integrals over flucton radii. The weight function is extracted directly from the radial density profile, which makes the formalism a practical bridge from nuclear-matter energy-density functionals to finite-nucleus quantities. Within this construction, CDFM has been used for nuclear symmetry energy and its volume and surface components, neutron pressure, incompressibility, giant monopole resonance energies, and superscaling analyses of inclusive electron and neutrino scattering [1607.03325] [2406.07190].

## 1. Formal definition and flucton representation

In CDFM, a nucleus of mass number $A$ is decomposed into fluctons of radius $x$, each carrying a constant density

$$
\rho_0(x)=\frac{3A}{4\pi x^3}.
$$

The one-body density matrix is written as a coherent superposition,

$$
\rho(\mathbf r,\mathbf r')=\int_0^\infty dx\,|F(x)|^2\,\rho_x(\mathbf r,\mathbf r'),
$$

with flucton density matrix

$$
\rho_x(\mathbf r,\mathbf r')=
3\rho_0(x)\,
\frac{j_1\!\big(k_F(x)|\mathbf r-\mathbf r'|\big)}{k_F(x)|\mathbf r-\mathbf r'|}\,
\Theta\!\left(x-\frac{|\mathbf r+\mathbf r'|}{2}\right).
$$

Here $j_1$ is the spherical Bessel function, $\Theta$ is the step function, and the flucton Fermi momentum is

$$
k_F(x)=\left(\frac{3\pi^2}{2}\rho_0(x)\right)^{1/3}\equiv\frac{\beta}{x},
\qquad
\beta=\left(\frac{9\pi A}{8}\right)^{1/3}\simeq1.52\,A^{1/3}.
$$

The diagonal density follows as

$$
\rho(r)=\int_0^\infty dx\,|F(x)|^2\,\rho_0(x)\,\Theta(x-r).
$$

For monotonic densities, $d\rho/dr\le 0$, the CDFM weight function is obtained from the radial derivative,

$$
|F(x)|^2=-\frac{1}{\rho_0(x)}\left.\frac{d\rho(r)}{dr}\right|_{r=x},
\qquad
\int_0^\infty |F(x)|^2 dx=1.
$$

This derivative form is central: it ties the model directly to the empirical or self-consistent density profile and makes the weight function peak in the surface region, where the density falls most rapidly. As a result, CDFM averages are intrinsically sensitive to subsaturation densities and to changes in surface diffuseness, shell structure, deformation, and neutron-skin development [2006.01471] [2506.08736].

## 2. Mapping nuclear-matter observables to finite nuclei

CDFM applies a universal folding prescription to any nuclear-matter quantity that depends on density. If $P^{NM}(\rho)$ denotes a nuclear-matter observable, its finite-nucleus counterpart is

$$
P(A)=\int_0^\infty dx\,|F(x)|^2\,P^{NM}\!\big[\rho_0(x)\big].
$$

This prescription underlies applications to symmetry energy, symmetry pressure, curvature coefficients, and incompressibility. In local-density implementations, the flucton density may be written either as $\rho_0(x)$ or as a local density $\rho(x)$ associated with the flucton; operationally both refer to the density sampled by the weight function.

For the finite-nucleus symmetry energy coefficient one writes

$$
S(A)=\int_0^\infty dx\,|F(x)|^2\,S^{NM}\!\big[\rho_0(x)\big].
$$

In many symmetry-energy studies a power-law density dependence is adopted,

$$
S[\rho(x)]=S(\rho_0)\left[\frac{\rho(x)}{\rho_0}\right]^\gamma,
$$

which yields

$$
S(A)=S(\rho_0)\int_0^\infty dx\,|F(x)|^2
\left[\frac{\rho(x)}{\rho_0}\right]^\gamma.
$$

The same structure is used for the symmetry-pressure and curvature integrals,

$$
P(A)=\int_0^\infty dx\,|F(x)|^2\,P^{NM}\!\big[\rho_0(x)\big],\qquad
K_{\rm sym}(A)=\int_0^\infty dx\,|F(x)|^2\,K_{\rm sym}^{NM}\!\big[\rho_0(x)\big].
$$

Several nuclear-matter inputs have been folded in this way. These include the Brueckner energy-density functional, standard Skyrme functionals, polynomially fitted Skyrme or RMF-based energy-density functionals, and Brueckner-Hartree-Fock calculations with Bonn B and Bonn CD potentials. Across these implementations, the surface-peaked structure of $|F(x)|^2$ makes the half-density region, typically around $\rho/\rho_0\approx 0.5$, especially important for finite-nucleus averages [1811.02897] [2211.05671].

## 3. Symmetry energy, volume–surface decomposition, and ratio schemes

A major CDFM application concerns the decomposition of the finite-nucleus symmetry energy into volume and surface parts. In droplet-model language, one writes

$$
S(A)=\frac{S^V}{1+\dfrac{S^S}{S^V}A^{-1/3}},
$$

or equivalently through a ratio parameter. Two conventions appear in the CDFM literature. In one convention,

$$
\kappa=\frac{S^V}{S^S},
$$

while in another,

$$
q(A)=\frac{S^S(A)}{S^V(A)}.
$$

This notational difference is substantive because some CDFM papers report $\kappa=S^V/S^S$, whereas the later “scheme II” work reformulates the problem directly in terms of $q=S^S/S^V$.

The earlier CDFM treatment, often called scheme I, adapts a Danielewicz-type expression to finite nuclei:

$$
\kappa(A)=
\frac{3}{r_0\rho_0\displaystyle\int_0^\infty dx\,|F(x)|^2\,x\,\rho_0(x)
\left\{
\frac{S^{NM}(\rho_0)}{S^{NM}[\rho_0(x)]}-1
\right\}}.
$$

This scheme was useful operationally, but later work identified conceptual and technical issues: it relied on replacing a half-infinite-matter density profile by a finite-nucleus CDFM density, and for some energy-density functionals the integrand could develop singular behavior.

The later scheme II derives the surface-to-volume ratio at the flucton level and then folds it with the same weight function. Its key flucton-level ingredient is

$$
q(x)=
\frac{x}{r_0}
\left(\frac{\rho_0(x)}{\rho_0}\right)^{1/3}
\left[
1-\frac{S^{NM}[\rho_0(x)]}{S^{NM}(\rho_0)}
\right],
$$

followed by

$$
q(A)=\int_0^\infty dx\,|F(x)|^2\,q(x).
$$

The finite-nucleus components are then reconstructed as

$$
S^V(A)=S(A)\left(1+\frac{q(A)}{A^{1/3}}\right),\qquad
S^S(A)=q(A)\,S(A)\left(1+\frac{q(A)}{A^{1/3}}\right).
$$

Scheme II was proposed as a more direct and physically motivated formulation that avoids the preliminary assumptions and mathematical ambiguities of scheme I. In applications to Ni, Sn, and Pb isotopic chains with SLy4 densities, it yields $S(A)$ in the interval $24$–$31$ MeV, $S^V(A)$ in the range $29$–$34$ MeV, and, for the Skyrme EDF, $S^S(A)$ in the range $14$–$18$ MeV. Its reported values of $1/q=S^V/S^S$ span about $1.0$–$2.4$ depending on chain and interaction, and the resulting systematics were described as closer to broad phenomenological ranges than the older scheme I values [2110.00456].

The same CDFM machinery also captures isotopic structure in the symmetry-energy observables. Pronounced kinks in $S(A)$, $S^V(A)$, $S^S(A)$, and the ratio are found near the doubly magic nuclei $^{78}$Ni and $^{132}$Sn, whereas no comparable kink is seen in the Pb chain. In finite-temperature calculations, these kinks blur and disappear as temperature rises from $T=0$ to $3$ MeV, consistent with thermal smearing of shell structure and increasing surface diffuseness [1811.02897].

## 4. Microscopic implementations and nuclear-structure applications

CDFM does not prescribe a unique microscopic source for $\rho(r)$; rather, it uses whatever self-consistent density is supplied by a chosen many-body framework. In the symmetry-energy applications of Antonov and collaborators, proton and neutron densities were generated with self-consistent HF+BCS calculations using Skyrme interactions such as SLy4, SGII, and Sk3, and then combined into a total density from which $|F(x)|^2$ was constructed [1607.03325]. In the finite-temperature extension, the temperature-dependent proton and neutron densities were obtained with the HFBTHO code, which solves the Skyrme Hartree-Fock-Bogoliubov problem in a cylindrical transformed deformed harmonic-oscillator basis, using SkM* and SLy4 and temperatures $T=0,1,2,3$ MeV [1811.02897].

Beyond the Ni, Sn, and Pb chains, CDFM has been applied to light and superheavy systems. In light deformed Ne, Na, Mg, Al, and Si isotopes, RMF+BCS densities with the NL3* parameter set were converted to spherical-equivalent densities through a two-Gaussian fit before the weight function was extracted. In these nuclei, the symmetry energy, neutron pressure, and symmetry-energy curvature coefficient show shell-related structure near $N=14$, $16$, $20$, and $28$, but the overall trends are irregular. The discussion attributes this behavior to strong deformation, shape coexistence, weak binding near the drip lines, and the sensitivity of $|F(x)|^2$ to the spherical-equivalent reconstruction. The same study connected enhanced symmetry energies in isotopes such as $^{29}$F, $^{28}$Ne, $^{29,30}$Na, and $^{31,35,36}$Mg with the region identified as the island of inversion [2106.03038].

In superheavy $Z=120$ isotopes, CDFM was used together with SHF densities from SLy4 and SkMP, and with E-RMF densities from G3 and IOPB-I. When a newly fitted Skyrme EDF was folded with SHF densities, the finite-nucleus symmetry energy, neutron pressure, and curvature exhibited a peak at $N=182$; earlier E-RMF-based CDFM results gave a corresponding peak at $N=184$. In contrast, Brueckner-EDF folding made these peaks absent or shifted, often toward $N\approx172$, and this was correlated with the Coester-band problem. In the same work, the volume contribution $S_V$ generally exceeded the surface contribution $S_S$ across the chain, with a visible $N=182$ peak in $S_V$ for SkMP [2211.05671].

A complementary RMF line of work constructed polynomial or simplified analytical expressions for the nuclear-matter energy per nucleon and then folded the resulting $K^{NM}$, $S^{NM}$, $L_{\rm sym}^{NM}$, and $K_{\rm sym}^{NM}$ with the CDFM weight. For closed and semi-closed-shell nuclei from $^{16}$O to $^{208}$Pb, these studies found that $P^A$, $S^A$, and $L_{\rm sym}^A$ lie in relatively narrow domains, whereas $K^A$ and $K_{\rm sym}^A$ vary more strongly with mass and with the underlying RMF parameter set, reflecting the stronger sensitivity of second derivatives to the density profile and the stiffness of the equation of state [2002.02135] [2403.02729].

## 5. Incompressibility and giant monopole resonance

CDFM has also been used to map nuclear-matter incompressibility into finite nuclei. The fundamental folding is

$$
K^A=\int_0^\infty dR\,|F(R)|^2\,K\!\big(\rho_0(R)\big),
$$

with

$$
K(\rho)=9\rho^2\frac{d^2}{d\rho^2}\left(\frac{E}{A}\right).
$$

This finite-nucleus incompressibility has then been related to the isoscalar giant monopole resonance (ISGMR) energy through two forms. One is the scaling relation

$$
E_{\rm ISGMR}=\hbar\sqrt{\frac{K^A}{m\langle r^2\rangle}},
$$

with

$$
\langle r^2\rangle=\frac{N}{A}\langle r_n^2\rangle+\frac{Z}{A}\langle r_p^2\rangle.
$$

The other is a Brueckner-like radius-parameter form,

$$
E_{\rm ISGMR}=\frac{\hbar}{r_0A^{1/3}}\sqrt{\frac{K^A}{m}},
$$

with the mass-dependent parameterization

$$
r_0(A)=1+\frac{2.40}{A^{2/3}}.
$$

Using Brueckner and BCPM nuclear-matter functionals together with SLy4 HF+BCS densities, CDFM calculations of $K^A$ and the corresponding ISGMR energies were reported for nuclei from Ca to Pb. With the scaling relation, Brueckner-based results were closer to experiment for lighter and mid-mass nuclei such as $^{40}$Ca and $^{58}$Ni, whereas BCPM gave higher energies and, with the mass-dependent $r_0(A)$ form, a more uniform overall agreement across isotopic chains such as Cd and Sn [2506.08736].

A distinct E-RMF-based CDFM analysis reported an anomalous trend: the calculated ISGMR energy $E_M$ was largest for IOPB-I, intermediate for G3, and smallest for NL3, even though the corresponding nuclear-matter incompressibilities satisfy $K_\infty({\rm NL3})>K_\infty({\rm G3})>K_\infty({\rm IOPB\mbox{-}I})$. The interpretation offered there is that vector-meson self- and cross-interactions modify the density dependence of $K(\rho)$ in the subsaturation surface region that dominates the CDFM integral, so the ordering of $K^A$ need not follow the ordering of $K_\infty$ at saturation [2508.05291].

These incompressibility applications underscore a recurrent CDFM theme: finite nuclei probe the density dependence of nuclear-matter quantities over the surface-dominated, subsaturation regime rather than only at the saturation point. This same logic motivates the use of neutron-skin information as an auxiliary input in ISGMR analyses, although the reported sensitivity of $E_{\rm ISGMR}$ to currently available neutron-skin uncertainties remains modest [2506.08736].

## 6. Superscaling, quasielastic scattering, and relativistic extensions

CDFM has also been reformulated as a finite-nucleus extension of the relativistic Fermi gas for inclusive electron and neutrino scattering. In this context, the scaling function is constructed as a coherent average of RFG kernels over the same density-fluctuation weight:

$$
f^{\rm QE}(\psi^*)=
\int_0^\infty dx\,|F(x)|^2\,f_{\rm RFG}^{\rm QE}[\psi^*(x)],
$$

with

$$
f_{\rm RFG}^{\rm QE}[\psi^*(x)] =
\frac{3}{4}
\left[
1-\left(\frac{k_Fx\psi^*}{\alpha}\right)^2
\right]
\theta\!\left(1-\left(\frac{k_Fx\psi^*}{\alpha}\right)^2\right),
$$

$$
\alpha=(9\pi A/8)^{1/3}\simeq1.52A^{1/3},
\qquad
\psi^*(x)=\frac{k_Fx\psi^*}{\alpha}.
$$

Here the global Fermi momentum is not a fit parameter but is obtained self-consistently from the density,

$$
k_F=\int_0^\infty dx\,|F(x)|^2\,\frac{\alpha}{x}.
$$

In the interacting-RFG version used for $^{12}$C, the scaling variable $\psi^*$ is built with a relativistic effective mass $m_N^*=0.8m_N$. The dimensionless variables are

$$
\eta_F^*=k_F/m_N^*,\qquad
\xi_F^*=\sqrt{1+\eta_F^{*2}}-1,\qquad
\kappa^*=q/(2m_N^*),\qquad
\lambda^*=\omega/(2m_N^*),\qquad
\tau^*=\kappa^{*2}-\lambda^{*2},
$$

and

$$
\psi^*=
\sqrt{\frac{\epsilon_0^*-1}{\epsilon_F^*-1}\,
{\rm sgn}(\lambda^*-\tau^*)}.
$$

The nuclear responses then factorize as

$$
R_K=
\frac{{\cal N}\,\xi_F^*}{m_N^*\eta_F^{*3}\kappa^*}
\,U_K\,f^{\rm QE}(\psi^*),
$$

so that the same CDFM scaling function enters both electron and charged-current neutrino cross sections.

With this formulation, CDFM has been applied to quasielastic $(e,e')$, $\nu_\mu$, and $\bar\nu_\mu$ scattering on $^{12}$C. The calculations use the empirical $^{12}$C density to obtain $|F(x)|^2$, include a relativistic effective mass $m_N^*=0.8m_N$, and add $2p$–$2h$ meson-exchange currents through external parametrizations. Reported comparisons with MiniBooNE, T2K, and MINER$\nu$A show good agreement over broad kinematic ranges. The quoted MEC contributions are roughly $20$–$25\%$ at the maximum for MiniBooNE-like kinematics, about $10\%$ in most T2K bins with values up to $\sim25\%$ at very forward angles, and roughly $30$–$50\%$ for MINER$\nu$A depending on kinematics and on the choice of the axial $N\to\Delta$ form factor parameter $C_5^A(0)$ [2406.07190] [2504.17141] [2604.15989].

This scattering program broadens the meaning of CDFM. In symmetry-energy work it functions as a density-to-observable map; in superscaling analyses it becomes a finite-nucleus extension of the relativistic Fermi gas in which density fluctuations, rather than a single sharp Fermi sphere, generate the nuclear response.

## 7. Limitations, conventions, and current methodological issues

Several limitations recur across CDFM applications. First, the model is usually combined with a local-density approximation. This is explicit in symmetry-energy and incompressibility calculations and implies that gradient corrections are neglected; such corrections may be non-negligible in the surface region and at finite temperature [1811.02897]. Second, the weight function requires a monotonic radial density. For spherical nuclei this is straightforward, but for deformed systems the density is often projected to a spherical-equivalent form, for example by two-Gaussian fitting, and this projection can alter the extracted $|F(x)|^2$ and the resulting surface observables [2106.03038].

Third, results are sensitive to the nuclear-matter input. In symmetry-energy decompositions, the choice of density dependence $S(\rho)$ is critical. Power-law forms with $\gamma\approx0.3$–$0.4$ and Tsang-type forms yield similar $S(A)$ and ratio systematics because they are close in the half-density region where $|F(x)|^2$ peaks, whereas M3Y-like or stiffer power-law forms can generate much larger ratio values and lower finite-nucleus symmetry energies [1811.02897]. This sensitivity is also visible in the superheavy $Z=120$ study, where Brueckner and fitted Skyrme inputs generate qualitatively different peak structures [2211.05671].

Fourth, there are nontrivial convention issues in the volume–surface decomposition itself. Some CDFM papers define $\kappa=S^V/S^S$, while scheme II introduces $q=S^S/S^V$ and reports $1/q=S^V/S^S$. A common misconception is that the different formulations are merely notational. In fact, they reflect different derivational routes, and scheme II was introduced precisely because the older scheme I involved a preliminary assumption and could suffer from integrand singularities for some EDFs [2110.00456].

Finally, in the scattering extensions, the quasielastic CDFM framework omits some channels or treats them approximately. Reported calculations neglect or defer explicit treatments of final-state interactions, long-range correlations such as RPA, and inelastic $\Delta$ or DIS sectors in selected applications. The low-$q$ region, especially at very forward angles, is identified as a regime where pure superscaling-based factorization is less reliable and where collective effects may need to be added [2406.07190] [2604.15989].

Taken together, these points place CDFM in a specific methodological niche. It is neither a stand-alone microscopic many-body theory nor a purely phenomenological fit. Its distinctive content lies in the coherent superposition of fluctons and the extraction of a surface-sensitive weight function from the density profile. That architecture has made it useful in domains as diverse as symmetry-energy decomposition, finite-temperature shell effects, superheavy shell indicators, incompressibility and ISGMR systematics, and finite-nucleus superscaling analyses.

Source: https://www.emergentmind.com/topics/coherent-density-fluctuation-model-cdfm