Local Mean-Field Density Functional Theory

• Local mean-field density functional theory is a framework where nonlocal many-body interactions are approximated by local densities and effective fields, enabling simplified descriptions of complex systems.
• It encompasses various formulations such as LMF–DFT for non-uniform liquids, SLDA for superfluid fermions, and DMFT-based lattice DFT for strongly correlated systems, each tailored to specific physical contexts.
• The approach offers computational efficiency and direct insights into structure and dynamics, but its precision relies on the careful division of potentials and controlled approximations like the local density theorem.

Searching arXiv for recent and foundational papers on local mean-field density functional theory and closely related formulations. Local mean-field density functional theory denotes a class of density-functional constructions in which inhomogeneous many-body systems are represented through local densities, local order parameters, or local effective one-body fields, while interaction effects are treated by mean-field, local-density, or auxiliary-field closures. In superfluid fermionic systems this viewpoint is realized by the Superfluid Local Density Approximation (SLDA), described as a “fully local, orbital-based density functional theory for superfluid fermions” (Bulgac, 2012). In classical non-uniform liquids, the Local Molecular Field (LMF) equation for an effective reference potential follows directly from the standard mean-field DFT treatment of attractive forces (Archer et al., 2012). In correlated lattice systems, local exchange-correlation potentials derived from Dynamical Mean Field Theory (DMFT), as well as density-and-pair-density functionals, provide local Kohn–Sham-like descriptions of Hubbard models (Karlsson et al., 2010, Lorenzana et al., 2012). A rigorous local-density approximation (LDA) theorem for short-range classical gases further establishes when the free energy with prescribed profile $\rho(x)$ is approximated by $\int f_T(\rho(x))\,dx$ (Jex et al., 2023).

1. Common formal structure

Across its major realizations, local mean-field DFT begins from a functional of one-body information and replaces either explicit many-body correlations or nonlocal interaction terms by local functionals, effective local fields, or local auxiliary equations. In classical DFT one writes the grand potential as

$$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$

with the Euler–Lagrange equation

$$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$

where $c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$ (Archer et al., 2012). In quantum formulations, one instead minimizes an energy functional with respect to orbitals, quasiparticle amplitudes, or generalized densities, obtaining Kohn–Sham, Bogoliubov–de Gennes, or collective Schrödinger equations (Bulgac, 2012, Lesinski, 2013).

The local variables differ substantially by domain, but the operational motif is recurrent: a uniform reference problem supplies either the free-energy density, the exchange-correlation energy, the pairing coupling, or a kinetic-renormalization factor; the inhomogeneous system is then treated through pointwise evaluation or local self-consistency. This suggests that “local” in this literature refers primarily to the representation of the functional and the resulting equations, not necessarily to the bare interaction itself.

Formulation	Local variables or fields	Systems emphasized
LMF–DFT	$\rho(r)$, $\phi_R(r)$	Non-uniform liquids (Archer et al., 2012)
Classical LDA	$\rho(x)$, $f_T(\rho)$	Short-range classical gases (Jex et al., 2023)
SLDA / ASLDA	$\rho(r)$, $\int f_T(\rho(x))\,dx$0, $\int f_T(\rho(x))\,dx$1, $\int f_T(\rho(x))\,dx$2	Unitary Fermi gas, nuclei (Bulgac, 2012, Bulgac et al., 2010)
DMFT-based lattice DFT	$\int f_T(\rho(x))\,dx$3, $\int f_T(\rho(x))\,dx$4	3D Hubbard model (Karlsson et al., 2010)
DPDFT	$\int f_T(\rho(x))\,dx$5, $\int f_T(\rho(x))\,dx$6, $\int f_T(\rho(x))\,dx$7	Inhomogeneous Hubbard model (Lorenzana et al., 2012)
Collective-coordinate DFT	$\int f_T(\rho(x))\,dx$8	Symmetry breaking and configuration mixing (Lesinski, 2013)
Local DFT-SIR-DH	$\int f_T(\rho(x))\,dx$9, $$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$0	Epidemic spreading with social interactions (Xu et al., 7 Sep 2025)

2. Classical non-uniform liquids and the LMF–DFT correspondence

The clearest classical local mean-field construction is the reformulation of Local Molecular Field Theory within DFT. The pair potential is split into a reference part and a remainder,

$$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$1

and a standard mean-field approximation is made for the excess free-energy difference,

$$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$2

Substitution into the Euler–Lagrange equation yields the effective reference field

$$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$3

or, equivalently,

$$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$4

In this form, LMF is not an external alternative to DFT but a restatement of the usual mean-field DFT treatment of attractions in terms of an effective one-body field (Archer et al., 2012).

Archer and Evans illustrated this correspondence for a hard-core Yukawa liquid adsorbed at a planar hard wall, using Rosenfeld’s Fundamental-Measure Theory as the hard-sphere reference functional. The self-consistent procedure is explicit: start from a trial $$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$5, solve the one-dimensional Euler–Lagrange equation for the hard-sphere profile $$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$6, update $$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$7 through the LMF equation involving the laterally integrated Yukawa tail, and iterate to self-consistency. The resulting density profile and thermodynamics are identical to those obtained by direct minimization of a DFT functional with

$$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$8

(Archer et al., 2012).

Near bulk liquid–gas coexistence, the same framework predicts a repulsive hump in $$\Omega[\rho]=F_{\rm id}[\rho]+F_{\rm ex}[\rho]+\int dr\,\rho(r)\bigl(\phi_{\rm ext}(r)-\mu\bigr),$$9 immediately outside the hard wall, of height $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$0, generating the low-density gas film characteristic of complete drying. For sufficiently long-range attraction the tail of $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$1 decays monotonically, whereas for shorter-range attractions it exhibits damped oscillations, reflecting the Fisher–Widom crossover in the asymptotic decay of $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$2. The one-dimensional exactly solvable model studied in the same work shows that the choice of splitting matters quantitatively: placing only a small fraction $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$3 of the attractive tail into the mean-field term, with $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$4, yields very accurate pressure and structure factor, whereas the standard RPA choice $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$5 does not (Archer et al., 2012). A persistent misconception is therefore that locality alone guarantees accuracy; in the LMF setting, the decisive ingredient is the “intelligent division” of the pair potential.

3. Local-density approximation as a controlled limit

A complementary route to local mean-field DFT is the LDA itself. For a short-range classical gas at temperature $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$6, the free energy with prescribed density profile can be approximated by the local functional

$$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$7

where $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$8 is the free energy per unit volume of the infinite homogeneous gas of density $$\ln[\Lambda^3\rho(r)]-c^{(1)}(r)+\phi_{\rm ext}(r)=\mu,$$9 (Jex et al., 2023). The rigorous result applies under assumptions of superstability, stability, and lower-regularity of the interaction, together with slow spatial variation of the density.

The theorem is quantitative. Introducing the local variation function

$c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$0

one obtains upper and lower bounds for $c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$1 in terms of $c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$2 and $c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$3, multiplied by an explicit decay factor $c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$4. The proof combines a Legendre-dual representation, coarse-graining into cubes, local thermodynamic-limit comparison in “main” cubes, and quantitative Ruelle-type bounds controlling local number fluctuations and inter-cube interactions (Jex et al., 2023).

This result sharpens the status of LDA in local mean-field theories. It is not merely a heuristic replacement of a nonlocal functional by a pointwise one; for slowly varying profiles it becomes an asymptotic statement with controlled errors. The same analysis indicates where the approximation should fail or need refinement: nonlocal corrections must be tied to density variation and are naturally expressed as gradient terms of the form

$c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$5

whose appearance is described in the discussion of refined cluster-expansion or gradient-expansion arguments (Jex et al., 2023).

4. Superfluid fermions and the Superfluid Local Density Approximation

In quantum many-body physics, the canonical local mean-field density functional is SLDA. Its basic local densities are

$c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$6

and the total energy is written as

$c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$7

The local superfluid energy density contains kinetic, particle-hole, and pairing parts. In the unitary Fermi gas (UFG), the unrenormalized SLDA functional takes the form

$c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$8

while in nuclear applications the particle-hole and pairing parts are generalized to isospin-dependent functionals involving $c^{(1)}(r)=-\delta F_{\rm ex}[\rho]/\delta\rho(r)$9 (Bulgac, 2012).

A central structural feature is ultraviolet renormalization. Both $\rho(r)$0 and $\rho(r)$1 diverge, but the combination

$\rho(r)$2

is finite, with $\rho(r)$3 the local pairing field. Because the pairing energy is proportional to $\rho(r)$4, SLDA yields the local gap equation

$\rho(r)$5

so the order parameter is purely local after renormalization (Bulgac, 2012).

Variation of $\rho(r)$6 gives the Kohn–Sham–Bogoliubov equations

$\rho(r)$7

$\rho(r)$8

with $\rho(r)$9. In the UFG, the three dimensionless constants $\phi_R(r)$0 are fixed by matching infinite-matter QMC data for the energy per particle, the pairing gap, and the quasiparticle spectrum. A related detailed SLDA fit reports $\phi_R(r)$1, $\phi_R(r)$2, $\phi_R(r)$3, $\phi_R(r)$4, and $\phi_R(r)$5 for the homogeneous system, while the 2012 review summarizes the homogeneous limit as recovering QMC energy with $\phi_R(r)$6 and gap $\phi_R(r)$7 (Bulgac et al., 2010, Bulgac, 2012).

Locality in SLDA is an approximation, not an exact microscopic statement. The pairing channel is taken to be zero-range, justified because typical Cooper-pair sizes in dilute or nuclear matter greatly exceed the microscopic interaction range, and gradient corrections are omitted at LDA level because numerics and QMC benchmarks show them to be small in the UFG (Bulgac, 2012). The asymmetric extension ASLDA promotes $\phi_R(r)$8 to polarization-dependent functions and is used when $\phi_R(r)$9, while the time-dependent extension TD-SLDA introduces a time-dependent action and a Galilean-invariant replacement of $\rho(x)$0 by $\rho(x)$1 to describe non-adiabatic dynamics, superfluid-to-normal transitions, and modes inaccessible in Landau–Ginzburg, Gross–Pitaevskii, or quantum hydrodynamics (Bulgac et al., 2010).

5. Strongly correlated lattice formulations

For the three-dimensional Hubbard model, a local mean-field density-functional scheme can be built by extracting the exchange-correlation potential from DMFT. The lattice DFT energy functional is

$\rho(x)$2

and in LDA one sets

$\rho(x)$3

Here $\rho(x)$4 is taken from the homogeneous 3D Hubbard model via DMFT,

$\rho(x)$5

When $\rho(x)$6 exceeds a critical $\rho(x)$7, $\rho(x)$8 develops a cusp at $\rho(x)$9, and the derivative $f_T(\rho)$0 jumps by a finite amount $f_T(\rho)$1, providing the DFT signature of the Mott–Hubbard metal–insulator transition at half-filling (Karlsson et al., 2010).

The time-dependent counterpart uses the adiabatic Local Density Approximation. The Kohn–Sham orbitals satisfy

$f_T(\rho)$2

with

$f_T(\rho)$3

This construction is explicitly memory-free: the static $f_T(\rho)$4 is reused at the instantaneous density. For a $f_T(\rho)$5-site cluster with a single interacting impurity at the center, the agreement with exact real-time densities is reported as excellent for moderate $f_T(\rho)$6 and smooth pulses; near half-filling and across the discontinuity at $f_T(\rho)$7, non-local and memory effects grow important. In Bloch-oscillation simulations on a $f_T(\rho)$8 lattice, the ALDA-DMFT dynamics reproduces $f_T(\rho)$9 and interaction-induced beat notes at $\rho(r)$0, but does not capture damping because genuine memory effects are absent (Karlsson et al., 2010).

A more elaborate local many-body embedding is the self-energy self-consistent DFT+DMFT scheme. There the full Green’s function is

$\rho(r)$1

and the DMFT self-energy is linearized near $\rho(r)$2 and inserted into the Kohn–Sham exchange-correlation potential for the correlated subspace. The scheme sets $\rho(r)$3, which removes the usual double-counting ambiguity. In SrVO$\rho(r)$4, the method yields $\rho(r)$5, keeps the $\rho(r)$6 quasiparticle bandwidth and Hubbard peaks at essentially the same positions as in one-shot calculations, and pushes the O $\rho(r)$7 bands down to $\rho(r)$8 (6 orbitals) and $\rho(r)$9 (3 orbitals), in excellent agreement with photoemission, without any ad hoc $\int f_T(\rho(x))\,dx$00 shift (Bhandary et al., 2019). Although this framework is more elaborate than an ordinary LDA, it remains local in the DMFT sense: the correlated physics enters through a local self-energy and its local feedback into the charge-density cycle.

6. Beyond ordinary density: pair density and collective coordinates

One route beyond ordinary local-density functionals is to enlarge the local variable set. In the adaptive pair-density functional theory of the inhomogeneous Hubbard model, the basic variables are the local density $\int f_T(\rho(x))\,dx$01 and the local double occupancy $\int f_T(\rho(x))\,dx$02. After a Legendre transform in the site-dependent interaction strengths $\int f_T(\rho(x))\,dx$03, one obtains the exact functional

$\int f_T(\rho(x))\,dx$04

The interacting kinetic energy is written through a reduction factor,

$\int f_T(\rho(x))\,dx$05

leading to Kohn–Sham-like equations with effective hoppings $\int f_T(\rho(x))\,dx$06. Under the site-factorized ansatz

$\int f_T(\rho(x))\,dx$07

the functional becomes formally identical to the Gutzwiller approximation, but the renormalization function $\int f_T(\rho(x))\,dx$08 is fixed from the exact uniform Bethe-ansatz solution through

$\int f_T(\rho(x))\,dx$09

(Lorenzana et al., 2012).

The numerical consequences are substantial. For a periodic one-dimensional chain of length $\int f_T(\rho(x))\,dx$10 with $\int f_T(\rho(x))\,dx$11 electrons and an alternating binary potential $\int f_T(\rho(x))\,dx$12, the per-site energy error remains below a few percent in DPDFT-FA even for $\int f_T(\rho(x))\,dx$13, whereas the LDA error can exceed $\int f_T(\rho(x))\,dx$14 in the strongly correlated regime $\int f_T(\rho(x))\,dx$15 and large $\int f_T(\rho(x))\,dx$16. In the single-electron impurity problem, exact physics requires $\int f_T(\rho(x))\,dx$17; LDA instead produces a spurious self-interaction energy saturating at $\int f_T(\rho(x))\,dx$18, while DPDFT-FA drives $\int f_T(\rho(x))\,dx$19 and hence $\int f_T(\rho(x))\,dx$20 as the impurity potential deepens (Lorenzana et al., 2012). The improvement is therefore tied to local adaptivity of pair density, not merely to a better scalar energy fit.

A second extension generalizes locality into a mixed single-particle/collective space. Defining the generalized density

$\int f_T(\rho(x))\,dx$21

one factorizes it as

$\int f_T(\rho(x))\,dx$22

where $\int f_T(\rho(x))\,dx$23 is a collective wave function and $\int f_T(\rho(x))\,dx$24 a normalized density slice satisfying $\int f_T(\rho(x))\,dx$25. The exact energy can then be written as

$\int f_T(\rho(x))\,dx$26

up to the rearrangement term. Minimization yields a local Schrödinger equation for $\int f_T(\rho(x))\,dx$27 and, at each $\int f_T(\rho(x))\,dx$28, Kohn–Sham equations

$\int f_T(\rho(x))\,dx$29

with $\int f_T(\rho(x))\,dx$30 (Lesinski, 2013). The translational example recovers $\int f_T(\rho(x))\,dx$31 and clean separation of center-of-mass and intrinsic dynamics, showing that locality can be formulated in an augmented collective manifold rather than only in ordinary real space.

7. Benchmarks, extensions, and recurring limitations

Local mean-field DFT methods are sustained by benchmark comparisons to controlled calculations. In the UFG, SLDA reproduces QMC energies in spherical traps with $\int f_T(\rho(x))\,dx$32 up to $\int f_T(\rho(x))\,dx$33 and errors $\int f_T(\rho(x))\,dx$34; sample values include $\int f_T(\rho(x))\,dx$35 with QMC $\int f_T(\rho(x))\,dx$36 and SLDA $\int f_T(\rho(x))\,dx$37 ($\int f_T(\rho(x))\,dx$38), $\int f_T(\rho(x))\,dx$39 with QMC $\int f_T(\rho(x))\,dx$40 and SLDA $\int f_T(\rho(x))\,dx$41 ($\int f_T(\rho(x))\,dx$42), and $\int f_T(\rho(x))\,dx$43 with QMC $\int f_T(\rho(x))\,dx$44 and SLDA $\int f_T(\rho(x))\,dx$45 ($\int f_T(\rho(x))\,dx$46). Time-dependent SLDA has been used for Anderson–Higgs oscillations of $\int f_T(\rho(x))\,dx$47, Bogoliubov sound modes, vortex crossing and reconnection, quantum turbulence, and quantum shock waves and domain walls in colliding clouds. In nuclear systems, the same framework describes odd–even mass staggering across 200+ nuclei with a single volume pairing strength, vortex cores in neutron matter, and giant dipole resonances in deformed, open-shell heavy nuclei in excellent agreement with experiment and with no adjustable parameters (Bulgac, 2012).

The same localization logic has recently been exported beyond traditional condensed-matter and nuclear settings. A local mean-field density functional theory model for epidemic spreading replaces convolution kernels for social interactions by a sum-of-exponentials approximation and auxiliary Debye–Hückel equations,

$\int f_T(\rho(x))\,dx$48

thereby converting a nonlocal DDFT system into a local DFT-SIR-DH PDE system. Its linear stability analysis yields both the classical reaction-driven instability condition $\int f_T(\rho(x))\,dx$49 and a cross-interaction condition for energy-driven instability; in a simplified self-repulsive case, the theory proves the existence of a unique global-in-time classical solution remaining nonnegative and bounded (Xu et al., 7 Sep 2025). This suggests that the local mean-field DFT idiom is portable wherever convolution-dominated interactions can be recast through local auxiliary fields.

Several limitations recur across the literature. In classical LMF, the split $\int f_T(\rho(x))\,dx$50 is not innocuous, and the one-dimensional solvable model shows that treating too much attraction in mean field leads to unphysical van der Waals loops or negative compressibility (Archer et al., 2012). In SLDA, omitted gradient corrections are small in the UFG at LDA level, but the approximation is expected to degrade in very tight traps or strong inhomogeneities (Bulgac, 2012). In ALDA-DMFT, the absence of memory precludes damping of Bloch oscillations (Karlsson et al., 2010). In DPDFT, the factorized-$\int f_T(\rho(x))\,dx$51 ansatz is approximate except in the infinite-coordination Gutzwiller limit (Lorenzana et al., 2012). The rigorous LDA theorem also clarifies that locality has a definite asymptotic regime—slowly varying densities and short-range interactions—rather than universal validity (Jex et al., 2023).

Taken together, these developments define local mean-field density functional theory not as a single formalism but as a family of mathematically related strategies: one starts from local densities or generalized local variables, imports uniform-system many-body information, and closes the inhomogeneous problem through local effective fields, local self-energies, local pair densities, or local auxiliary PDEs. The common payoff is computational tractability with direct access to structure, thermodynamics, and, in several cases, real-time dynamics; the common challenge is to control the information lost when genuinely nonlocal or history-dependent correlations are compressed into a local representation.