Universal Density-to-Current Map

• The topic establishes a bijective relationship between electron density and current, which is central to both static and time-dependent density functional theories.
• It highlights key challenges such as gauge-dependence and degeneracy in static systems that hinder a fully universal mapping using traditional approaches.
• Time-dependent frameworks, using nonlinear self-consistent Schrödinger equations and hydrodynamic methods, successfully yield explicit, one-to-one mappings for practical applications.

A universal density-to-current functional map is a fundamental construct in density functional theories that provides a one-to-one correspondence between an electronic density (or density history) and the quantum mechanical current density, possibly in the presence of electromagnetic fields or time dependence. This mapping is central in the development of both static and time-dependent (current-)density functional theories, determining whether the current, as a functional of the density, can be promoted to a fundamental variable alongside the density, and thus whether the full many-body system is determined by a sufficiently rich set of densities.

1. Foundational Concepts and Motivation

In conventional density functional theory (DFT), the ground-state properties of an interacting electron system are determined by the particle density $\rho(\mathbf r)$. However, in the presence of magnetic fields (vector potentials $\mathbf A(\mathbf r)$) or for fully time-dependent problems, the particle density alone is insufficient: either the gauge-invariant total current density $\mathbf j(\mathbf r)$ (static and time-dependent settings) or the paramagnetic current density $\mathbf j^{p}(\mathbf r)$ must also be considered as basic variables.

The motivation for establishing a universal density-to-current map is then twofold:

• To seek a rigorous Hohenberg–Kohn-like theorem involving both the density and current density.
• To construct universal functionals $F[\rho,\mathbf j]$ (or related maps) underpinning variational principles for the total energy or time evolution, suitable for practical density functional approximations.

2. Static Case: Obstacles and Negative Results

A major challenge in static ground-state CDFT is the proper and physically rigorous incorporation of the total current density as a basic variable. Classical approaches elevate the paramagnetic current density $\mathbf j^{p}$ (which is gauge-dependent) together with the density $\rho$; these variables are used in the universal Vignale–Rasolt functional $F_{\mathrm{VR}}[\rho, \mathbf j^{p}]$, which is rigorously convex, lower semicontinuous, and expectation-valued (Kvaal et al., 2020).

Attempts have been made to elevate the total (physical) current density, which is gauge-invariant, to the status of a basic variable in a universal functional. The principal proposal is Diener's variational principle, recast as a maximin (sup-inf) problem by Laestadius, Penz, and Tellgren (Laestadius et al., 2020):

• Define $F_{\rm D}(\rho, \mathbf k)$ as an infimum over paramagnetic currents penalized by their deviation from a trial total current field $\mathbf k$.
• Construct $$G_{\rm D}(\rho, \mathbf A) = \sup_{\mathbf k} \{ F_{\rm D}(\rho, \mathbf k) + \int \mathbf k \cdot \mathbf A \}$$.
• Hope that $G_{\rm D}(\rho, \mathbf A)$ reproduces the true Grayce–Harris functional $G_{\rm GH}(\rho, \mathbf A)$ and enables a universal mapping $(\rho, \mathbf j) \leftrightarrow (v, \mathbf A)$.

However, this framework fails due to convexity mismatch:

• $G_{\rm D}$ is always convex in $\mathbf A$, while $G_{\rm GH}$ is not; energy level crossings under magnetic fields (magnetic kinks) render $G_{\rm GH}$ non-convex.
• For some $(\rho, \mathbf A)$, $$G_{\rm GH}(\rho, \mathbf A) > G_{\rm D}(\rho, \mathbf A)$$; thus, minimization over $G_{\rm D}$ cannot yield correct ground-state energies.
• The associated Hohenberg–Kohn-type mapping also fails: there exist densities and vector potentials for which no minimizer of $F_{\rm D}(\rho, \mathbf k)$ satisfies the self-consistency condition relating $\mathbf k$ and the ground-state wave function's paramagnetic current.

Thus, within the Schrödinger variational framework, a universal total current density map is structurally precluded (Laestadius et al., 2020).

3. Degeneracies, Representability, and the Weak Hohenberg–Kohn Theorem

Even when considering the paramagnetic currents, the (density, current) pair $(\rho, \mathbf j^p)$ does not uniquely fix the set of degenerate ground states or associated external potentials. In the presence of degeneracy, it is possible to find different Hamiltonians (differing in, e.g., external magnetic fields) that yield the same pair $(\rho, \mathbf j^p)$ but have ground-state manifolds of different dimensions or structures.

Key findings (Laestadius et al., 2018):

• Degenerate ground states (with different angular momenta, for instance) can share observables $(\rho,\mathbf j^p)$.
• Explicit construction (e.g., 3 electrons in a 2D oscillator under varying field B) exhibits nontrivial degeneracy structures at fixed $(\rho,\mathbf j^p)$.
• Nevertheless, a weaker "ensemble" Hohenberg–Kohn theorem holds: for any ensemble-$(v, \mathbf A)$-representable $(\rho, \mathbf j^p)$, all Hamiltonians sharing that density pair share a nonempty set of ground states, at least as large as the rank of the representing density matrix.
• There is no universal one-to-one mapping $(\rho, \mathbf j^p) \mapsto (v, \mathbf A)$ when degeneracies are present.

This precludes a fully universal density-to-current mapping in the static case, except under nondegeneracy or stricter representability criteria (Laestadius et al., 2018).

4. Time-Dependent Frameworks: Existence and Construction of Universal Maps

The time-dependent case, especially within time-dependent current-density functional theory (TDCDFT), allows a much more robust framework for universal density-to-current mappings. Here the basic objects are the time-dependent density $n(\mathbf r, t)$ and physical current density $\mathbf j(\mathbf r, t)$.

The universal density-to-current (and current-to-density/potential) mapping is constructed nonperturbatively through a nonlinear self-consistent Schrödinger equation (NLSE) (Tokatly, 2011):

• For prescribed $\j(\r,t)$: solve the system

$\begin{aligned} i\hbar\partial_t\Psi(t) &= \hat H[\A]\Psi(t), \ \A(\r,t) &= \frac{m}{n(\r,t)}\langle\Psi(t)|\hat{\mathbf j}^p(\r)|\Psi(t)\rangle - \frac{m}{n(\r,t)}\j(\r,t), \end{aligned}$

under suitable initial conditions.

This produces a constructive map $$(\Psi_0, n, \mathbf j) \mapsto (\Psi[n, \mathbf j](t), \mathbf A[n,\mathbf j](\r,t))$$, yielding a one-to-one correspondence under mild regularity assumptions and ensuring $V$-representability and uniqueness. Analytic solutions and rigorous results are available for finite-dimensional systems and the lattice case, with the continuum still open to analysis (Tokatly, 2011).

5. Geometric and Hydrodynamic Approaches: Explicit Functional Forms

Recent advances provide direct nonlocal functional forms for the universal density-to-current map. In the geometric orbital-free formulation of TDDFT (Cancès et al., 12 Jan 2026):

• For a density evolution $\rho(\r, t)$, define the universal functional $$\mathcal{J}[\rho](\r, t) = j_{\Psi^{GP}[\rho]}(\r, t)$$, where $\Psi^{GP}[\rho](t)$ solves a constrained projected Schrödinger dynamics maintaining $\rho$.
• $\mathcal{J}[\rho]$ solves the continuity equation, with hydrodynamic velocity $v[\rho]$ extracted via an elliptic equation:

$\nabla\cdot[\rho v] = -\partial_t \rho,$

and thus

$\mathcal{J}[\rho] = \rho v[\rho].$

This hydrodynamic operator admits an explicit nonlocal Green’s function representation and, for prescribed initial states and interaction, yields a one-to-one and universal functional of the density (Cancès et al., 12 Jan 2026).

Moreover, this construction allows the definition of nonlocal mean-field operators in orbital-free or Kohn–Sham-type decompositions, providing viable alternatives to exchange–correlation potential-based methods.

6. Maxwell–Schrödinger DFT: Convex Analysis and Differentiable Functionals

By extending the variational degrees of freedom to include the induced Maxwell field, Maxwell–Schrödinger DFT (MDFT) formulates the problem in terms of the electron density $\rho(\mathbf r)$ and the total (internal plus external) magnetic field $\mathbf B(\mathbf r)$, equivalently its associated physical current. The universality emerges as follows (Tellgren, 2017):

• The ground-state energy is formulated as an infimal convolution (Moreau–Yosida regularization) of the standard Schrödinger energy and the magnetic self-energy, resulting in bounded curvature of the functional dependence on $\mathbf B$.
• Convex and concave analysis yields differentiable universal functionals $\bar f_{\rm MS}[\rho, \beta']$, with a one-to-one superdifferential/subdifferential structure relating density/current and potential/field variables—providing a rigorous Hohenberg–Kohn mapping even in the degenerate case.

The table below summarizes static versus time-dependent realizations:

Framework	Variables	Universal Map Status
Static (CDFT, $\mathbf j^p$)	$(\rho,\mathbf j^p)$	Well-defined functional, limited HK invertibility (Kvaal et al., 2020, Laestadius et al., 2018)
Static (total current)	$(\rho,\mathbf j)$	Attempts fail due to convexity/energy mismatch (Laestadius et al., 2020)
MDFT	$(\rho,\mathbf B)$ or $(\rho,\mathbf J)$	Full convexity, differentiability, HK-like mapping (Tellgren, 2017)
Time-dependent (TDCDFT, Geometric)	$n(\mathbf r, t)$, $\mathbf j(\mathbf r, t)$	Well-posed, one-to-one, explicit functionals (Tokatly, 2011, Cancès et al., 12 Jan 2026)

7. Implications, Practical Considerations, and Outlook

The existence or failure of a universal density-to-current map has both conceptual and practical consequences:

• In static ground-state theory, a universal map from $(\rho, \mathbf j)$ to external potentials is precluded by fundamental convexity mismatches and degeneracy structure; all functionals and practical calculations must be formulated in terms of $(\rho, \mathbf j^p)$ or employ extended frameworks (e.g., MDFT, QEDFT) (Laestadius et al., 2020, Tellgren, 2017).
• In time-dependent settings, geometric and NLSE-based constructions provide theoretically sound, uniquely defined, and practically applicable density-to-current functionals, with nonlocal but explicit kernels and operator structure (Cancès et al., 12 Jan 2026, Tokatly, 2011).
• Approaches such as MDFT, which regularize the magnetic field energy, restore full convexity and produce Fréchet-differentiable universal functionals, enabling robust variational principles and Hohenberg–Kohn-type theorems for the physical current (Tellgren, 2017).
• In practice, modern functionals, especially in the time-dependent context, may be usefully constructed using hydrodynamic, geometric, or Maxwell–Schrödinger-based principles, allowing greater representability and regularity.

A plausible implication is that future advances in density functional theory involving current densities—both static and time-dependent—will increasingly rely upon variational extensions (such as MDFT) and hydrodynamic/geometric mappings to construct reliable, nonlocal, and universally well-behaved functionals for practical use in electronic structure and quantum dynamics.