--- title: Dynamical Hubbard Functional Overview url: https://www.emergentmind.com/topics/dynamical-hubbard-functional type: topic --- # Dynamical Hubbard Functional Overview Searching arXiv for recent papers on “dynamical Hubbard functional” and related formulations. “Dynamical Hubbard Functional” is not a single formalism but a family of closely related constructions in which Hubbard correlations are encoded in an explicitly time-, frequency-, or coordinate-dependent functional rather than in a purely static on-site parameter. In the literature, the term is used for a nonequilibrium grand-potential functional of the self-energy on the Keldysh–Matsubara contour, a Klein or Luttinger–Ward-type Green’s-function functional with a local screened interaction \(U(\omega)\), a density-functional exchange–correlation potential derived from DMFT, a Hubbard augmentation kernel in TDDFPT, and a structural interaction law \(U[Q]\) obtained from cRPA and embedded in DFT+DMFT [1603.07177][2302.12193][1004.2264][2409.19504][2510.26584]. The common element is the replacement of a static Hubbard correction by a dynamical object whose stationarity, functional derivative, or parametric dependence generates equations of motion, self-energies, response kernels, or lattice-controlled effective interactions. ## 1. Conceptual scope and principal meanings The term covers several technically distinct uses. In nonequilibrium self-energy functional theory, the functional variable is the contour-ordered self-energy \(\Sigma\), and the central object is a grand-potential functional \(\Omega[\Sigma]\). In Green’s-function approaches for correlated solids, the functional variable is the interacting one-particle Green’s function \(G(\omega)\), with a local dynamical exchange–correlation contribution \(\Phi_{\text{dynH}}[\mathbf{G}]\). In lattice DFT and TDDFT, the “dynamical Hubbard functional” is an exchange–correlation potential \(v_{xc}(n)\) derived from a dynamical many-body reference problem, typically DMFT. In TDDFPT, the relevant object is instead the linear variation of a Hubbard energy functional, which enters the dynamical susceptibility through a kernel \(f_U\). In lattice-coupled correlated materials, one also encounters a coordinate-dependent interaction \(U[Q]\), where a structural amplitude controls screened Coulomb interactions and Hubbard bands [1603.07177][2302.12193][1004.2264][2409.19504][2510.26584]. | Formulation | Central variable | Representative papers | |---|---|---| | Nonequilibrium SFT | \(\Omega[\Sigma]\) | [1603.07177], [1510.05866] | | Spectral Green’s-function dynH | \(E_{\text{dynH}}[G]\), \(\Phi_{\text{dynH}}[\mathbf{G}]\) | [2302.12193], [2508.18194], [2503.10893] | | DMFT-based lattice DFT/TDDFT | \(E_{xc}[n]\), \(v_{xc}(n)\) | [1004.2264], [1903.04984] | | TDDFPT Hubbard response | \(E_U\), \(V_U\), \(f_U\) | [2409.19504] | | Structural interaction functional | \(U[a]\), \(U[a(t)]\) | [2510.26584] | A plausible implication is that the phrase is best understood as a category label for dynamical functionalizations of Hubbard physics rather than as the name of a unique formal theory. That plurality is explicit across the arXiv literature. ## 2. Nonequilibrium self-energy functionals and real-time Mott physics In the nonequilibrium SFT formulation of the Fermi–Hubbard model, the basic real-time object is a contour-ordered Green’s function on the Keldysh–Matsubara contour \(\mathcal{C}\), \[ G(1,2) = -i \langle \mathcal{T}_{\mathcal C} c(1)c^\dagger(2)\rangle, \] with Dyson equation \[ G = \big(G_0^{-1}-\Sigma\big)^{-1}. \] The dynamical Hubbard functional is the grand-potential functional \[ \Omega[\Sigma'] = \Omega' + \frac{1}{\beta}\,\mathrm{Tr}\,\ln\big(G_0^{-1}-\Sigma'\big)^{-1} - \frac{1}{\beta}\,\mathrm{Tr}\,\ln G', \] evaluated on trial self-energies \(\Sigma'\) generated by a solvable reference system. Stationarity with respect to the time-dependent one-body parameters \(\boldsymbol{\theta}(t)\) of the reference system, \[ \frac{\delta \Omega[\Sigma_{\boldsymbol{\theta}}]}{\delta \boldsymbol{\theta}(t)}=0, \] defines the optimal approximation. In the two-site dynamical impurity approximation, the only variational parameter is the hybridization \(V'(t)\) between a correlated impurity and a single auxiliary bath site [1603.07177]. For the half-filled Hubbard Hamiltonian \[ H(t) = -T \sum_{\langle ij \rangle,\sigma} c^{\dagger}_{i\sigma} c_{j\sigma} + U(t)\sum_i \Big(n_i - \tfrac{1}{2}\Big)\Big(n_i - \tfrac{1}{2}\Big), \] the formalism was applied to sudden quenches and finite-time cosine ramps of \(U(t)\), with \(\beta=10\), \(U_{\rm ini}=0.01\), a one-dimensional chain of \(L=40\) sites, and time windows up to \(t_{\max}\le 25\). The optimized \(V'_{\rm opt}(t)\), the double occupancy \(D(t)=\langle n_\uparrow(t)n_\downarrow(t)\rangle\), the momentum distribution \(n_k(t)\), and \(E_{\rm tot}(t)\) diagnose the dynamical Mott transition. Two regimes emerge: weak coupling with prethermalization and a persistent Fermi-surface jump, and strong coupling with collapse-and-revival oscillations at frequency \(\sim U_{\rm fin}\). These are separated by a sharp critical interaction \(U_{\rm c}^{\rm dyn}\), at which the low-energy bath decouples dynamically, \[ V'_{\rm opt}(t)\to 0 \;\Rightarrow\; \Delta(t,t')\to 0. \] For a sudden quench, \(U_{\rm c}^{\rm dyn}\approx 4.61\) in the 1D-chain units used in the paper, corresponding to \(\approx 3.26\) after rescaling to DOS variance; for a ramp with \(\tau_{\rm ramp}=7\), \(U_{\rm c}^{\rm dyn}\approx 8.02\), and \(U_{\rm c}^{\rm dyn}(\tau_{\rm ramp})\) increases monotonically toward the equilibrium zero-temperature upper spinodal \(U_{\rm c2}(0)\approx 8.59\) [1603.07177]. The nonequilibrium variational-cluster approach provides the methodological precursor for stable real-time propagation within this framework. There, direct solution of the Euler equation was found to be numerically unstable, and the practical condition was replaced by \[ K^{(0)}[\lambda'](t_0)=0,\qquad \partial_t K^{(0)}[\lambda'](t)=0\quad (t>t_0), \] which preserves causal time stepping because updates at time \(t\) do not affect earlier times \(t'U_c^{\text{Mott}}\approx 14\,t\). In the time-dependent setting, the potential is used within an adiabatic local density approximation, \[ v_{xc,i}(t)=v_{xc}^{\text{DMFT}}\big(n_i(t)\big), \] which was benchmarked on cluster dynamics and applied to Bloch oscillations in a \(33\times5\times5\) trapped cluster [1004.2264]. For the Hubbard–Holstein model, the functional variables are the site occupations \(n_i\) and phonon coordinates \(x_i\). The Lang–Firsov transformation produces a screened local interaction \[ U_{\mathrm{eff}}=U-\frac{2g^2}{\omega}, \] and the DFT construction yields \(\eta_{xc}=0\) for the phonon xc field. For the isolated site, the analytic electronic xc potential depends on \(U'=U-2g^2/\omega\), and the location of the discontinuity changes with the sign of \(U_{\mathrm{eff}}\): for \(U_{\mathrm{eff}}>0\) it is at \(n=1\), whereas for \(U_{\mathrm{eff}}<0\) it shifts to \(n=0\) and \(n=2\). The corresponding ALDA was shown to describe linear conductance and real-time dynamics with good accuracy against exact benchmarks [1903.04984]. In TDDFPT for magnons, the functional is instead a nonempirical Hubbard augmentation of the response kernel. In the Dudarev form, \[ E_U=\frac{U_{\mathrm{eff}}}{2}\sum_{I,\sigma}\mathrm{Tr}\big[n^{I,\sigma}(1-n^{I,\sigma})\big], \] which gives the Hubbard potential \[ \big(V_U^{I,\sigma}\big)_{mm'}=U_{\mathrm{eff}}\Big(\frac{1}{2}\delta_{mm'}-n^{I,\sigma}_{mm'}\Big). \] Its linear variation defines the kernel contribution \(f_U\), and the full susceptibility obeys \[ \chi(\mathbf{q},\omega)=\big[1-\chi_0(\mathbf{q},\omega)K(\mathbf{q},\omega)\big]^{-1}\chi_0(\mathbf{q},\omega), \] with \(K=f_H+f_{xc}+f_U\). The method was formulated in a general noncollinear setting, implemented through the Liouville–Lanczos approach, and reported to satisfy the Goldstone condition without empirical rescaling of the exchange–correlation kernel [2409.19504]. In this context, the functional is “dynamical” because the Hubbard correction enters the full frequency-dependent response, even though \(U\) and \(J\) themselves are static parameters determined nonempirically by linear response. ## 5. Structural and mode-resolved interaction functionals A more recent usage makes the Hubbard interaction itself a functional of a lattice coordinate. For monolayer 1T-TaS\(_2\), the star-of-David charge-density-wave amplitude is parameterized by \[ a=\frac{d-d_{\rm eq}}{d_{\rm eq}}, \] where \(d\) is the central–outer Ta distance in the 13-Ta cluster. Constrained RPA gives the partially screened interaction \[ U(\omega)=\langle \phi\phi|[1-vP_r(\omega)]^{-1}v|\phi\phi\rangle, \] projected onto the narrow molecular Wannier orbital that forms the half-filled correlated band. Near equilibrium, the static screened interaction is approximately linear in the distortion, \[ U[a]\approx U_0+ga, \] with \(U_0=0.40\) eV at \(a=0\) and \(g\approx -0.10\) eV/\% for \(a\in[-1\%,+2\%]\) [2510.26584]. The structural dependence is large. The bare on-site interaction decreases from \(\approx 4.2\) eV at \(a=-4\%\) to \(\approx 1.4\) eV at \(a=+3\%\). Over the interval \(a=-1\%\) to \(a=+2\%\), the Wannier spread grows from \(24\) \(\text{\AA}^2\) to \(95\) \(\text{\AA}^2\), the effective bandwidth \(W_{\rm eff}\) grows from \(\approx 20\) meV to \(\approx 70\) meV, and the \(k\)-resolved Mott gap shrinks from \(\approx 0.5\) eV to \(\approx 0.1\) eV. At \(a\approx +3\%\), DFT+DMFT yields an insulator-to-correlated-metal transition with a quasiparticle peak at \(E_F\) [2510.26584]. Because the coherent star-of-David amplitude mode has \(\Omega\approx2.4\) THz, the time-dependent interaction is written as \[ a(t)=A\cos(\Omega t+\phi),\qquad U[a(t)]\approx U_0+ga(t). \] For \(A\approx0.5\%\), the expected modulation is \(\delta U\approx-50\) meV; for \(A\approx1\%\), \(\delta U\approx-0.1\) eV. The same work therefore proposes a minimal driven Hubbard Hamiltonian in which both hopping and interaction are functionals of the instantaneous lattice coordinate, \[ H(t)=\sum_{ij,\sigma} H^{(R_i-R_j)}_{a(t)} c^\dagger_{i\sigma}c_{j\sigma}+U[a(t)]\sum_i n_{i\uparrow}n_{i\downarrow}. \] This use of “dynamical Hubbard functional” is literal: the screened interaction becomes a calculable function of a collective structural mode. ## 6. Exact, field-theoretic, and renormalization-group constructions The term also appears in more formal many-body settings. In the Kotliar–Ruckenstein slave-boson representation of the two-site extended Hubbard model, an exact imaginary-time functional integral was constructed in the radial gauge, with fermions, a complex \(d\)-boson, radial slave-boson amplitudes, and time-dependent constraint fields. The discrete-time action includes a regularized time-shifted form of the KR \(z\)-factors and a slice-by-slice projection operator that enforces the constraints exactly. In the continuum limit, the projected functional integral reproduces the exact partition function and exact Matsubara Green’s function of the operator formulation, while resolving the normal-ordering and square-root objections associated with the KR representation [2402.05701]. In the two-dimensional Hubbard model, a “dynamical functional” viewpoint arises within functional renormalization group calculations that retain full Matsubara-frequency dependence of effective interactions and gap functions. The scale-dependent effective action \(\Gamma_\Lambda\) is evolved in the symmetric regime, then continued below \(\Lambda_c\) with dynamical magnetic and pairing gap functions \(\Delta_m(\nu)\) and \(\Delta_p(\nu)\). For \(t'/t=-0.16\) and \(U=3t\), this produced a regime with robust \(d\)-wave pairing coexisting with antiferromagnetism between \(p\approx 8\%\) and \(p\approx 20\%\), magnetic order up to \(p\approx 20\%\), and a Kosterlitz–Thouless superconducting dome centered around \(p\approx 15\%\) hole doping [2010.06267]. Here the “dynamical” content is the explicit retention of full frequency dependence in the functional flow. A geometric and algebraic variant appears in the supergroup approach to the Hubbard model. There, even Hubbard operators are related to generators of rotations in four-dimensional space, translations are obtained by deformation and contraction to an ISO(3)-type algebra, and a supercoherent-state construction yields an effective functional for spinless fermions in the large-\(U\) regime. The resulting effective Lagrangian, \[ \mathcal{L}_{\rm eff}=\mathcal{L}_0+\mathcal{L}_1, \] contains bosonic Cartan-form contributions and fermionic gradient terms on a curved group manifold [1102.2522]. This is the most geometrically literal interpretation of a Hubbard functional as an effective action defined on a nontrivial supergroup space. ## 7. Limitations, comparisons, and prospective directions Across these formulations, dynamical Hubbard functionals trade locality and controllability against completeness. In nonequilibrium SFT with a single bath site, particle number and spin are conserved, but energy conservation is not exact unless a continuum of variational degrees of freedom is included; the two-site reference also limits long-time thermalization and misses nonlocal correlations [1603.07177]. In Green’s-function dynH for solids, the self-energy is strictly local to the correlated manifold, the interaction is typically reduced to an orbital-averaged scalar \(U(\omega)\), and nonlocal self-energy effects are absent; the self-consistent SrVO\(_3\) studies therefore identify nonlocal extensions, intersite corrections, and combinations with DFT+\(V\) as natural continuations [2508.18194][2503.10893]. Density-functional realizations inherit the limitations of local and adiabatic approximations. In the 3D Hubbard TDDFT construction, the ALDA lacks memory and spatial nonlocality, and the derivative discontinuity near half filling can produce sizable deviations under fast driving or in strongly inhomogeneous clusters [1004.2264]. In the Hubbard–Holstein DFT, the lattice functional is local and memoryless, and the attractive regime was treated analytically only for the isolated site, while fuller DMFT-based attractive functionals were deferred [1903.04984]. In TDDFPT for magnons, the dynamical response is improved by the Hubbard kernel, but the ground-state correction still uses static \(U\) and \(J\), not an explicit \(U(\omega)\) [2409.19504]. The structural \(U[Q]\) formulation for 1T-TaS\(_2\) makes a different approximation: DFT+DMFT uses only the static \(U(\omega\to0)\), with no dynamic \(U(\omega)\), no nonlocal \(V\), and no EDMFT corrections. The authors explicitly identify EDMFT or GW+EDMFT as directions for refining both frequency-dependent screening and momentum-dependent self-energy effects [2510.26584]. A plausible synthesis is that future “dynamical Hubbard functionals” will increasingly couple three ingredients that are separated in current implementations: local dynamical self-energies, nonlocal screening or correlations, and explicit coupling to collective coordinates such as phonons, magnons, or structural modes. In that broader sense, the subject is less a closed formalism than an evolving research program. Its unifying objective is to cast Hubbard correlation physics into a functional language that remains variational, self-consistent, and computationally tractable while retaining genuinely dynamical information—whether in real time, in frequency space, or along a driven structural coordinate. Source: https://www.emergentmind.com/topics/dynamical-hubbard-functional