Time-Dependent Hartree Equations

• Time-dependent Hartree equations are coupled nonlinear Schrödinger equations derived via the Dirac–Frenkel projection, modeling quantum dynamics with self-consistent mean-field potentials.
• They are formulated within a rigorous variational framework that guarantees global existence, energy conservation, and precise error estimates across Coulomb and polynomial coupling regimes.
• Applications include semiclassical coherent-state propagation, dense fermionic mean-field limits, and multiconfigurational extensions, making these equations central to modern quantum dynamics studies.

Searching arXiv for recent and foundational papers on time-dependent Hartree equations to support the encyclopedia article. Time-dependent Hartree equations are coupled nonlinear Schrödinger equations that govern the evolution of product-state approximations in quantum dynamics. In the nonrelativistic setting they arise as the time-dependent Dirac–Frenkel projection of the many-body Schrödinger flow onto the manifold of uncorrelated product states, and they provide an effective mean-field description in which the nonlinear potential is generated self-consistently by the instantaneous density (Carles et al., 2022). In closely related contexts they also appear as limiting equations for large fermionic systems, as the single-configuration limit of multiconfigurational propagation schemes, and as the direct-field reduction of Hartree–Fock-type theories (Hoang et al., 16 Jul 2025, Bonfanti et al., 2018, Geng et al., 2024).

1. Variational origin and basic equations

The basic variational setting is the manifold of product states

$$\mathcal M=\{\varphi_1\otimes\cdots\otimes\varphi_N\mid\varphi_j\in L^2(\Bbb R^d)\},$$

together with the Dirac–Frenkel condition that the residual $i\partial_t u-Hu$ be orthogonal to the tangent space of $\mathcal M$ (Carles et al., 2022). In this formulation, the Hartree equations are the projected dynamics obtained by restricting the full many-body Schrödinger equation to uncorrelated product states.

For the simplest two-factor case, the resulting system is

$\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$

with

$$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$

This form makes explicit that each factor evolves under a one-body Hamiltonian plus a mean field generated by the density of the other factor (Carles et al., 2022).

In the $N$-orbital convolution form, the system reads

$$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$

with

$$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$

Two cases singled out in the literature are the convolution/Coulomb case $V(x)=\pm|x|^{-1}$ in $d=3$, and a more general coupling $i\partial_t u-Hu$0 with polynomial growth and no convolution structure (Carles et al., 2022).

A recurring point in the literature is that the nonlinear Hartree coupling need not be perturbative. The general Cauchy theory in (Carles et al., 2022) is formulated precisely for a setting in which “the nonlinear coupling cannot be considered as a perturbation.” This distinguishes the time-dependent Hartree system from analyses that rely on small-coupling arguments.

2. Functional framework, conservation laws, and the Cauchy problem

The well-posedness theory developed for the Hartree approximation uses one-body Hamiltonians

$i\partial_t u-Hu$1

with $i\partial_t u-Hu$2 and graph norms

$i\partial_t u-Hu$3

and similarly in the $i\partial_t u-Hu$4 variable (Carles et al., 2022). The working space is

$i\partial_t u-Hu$5

For general couplings, the analysis assumes a relative bound

$i\partial_t u-Hu$6

and commutator bounds of the form

$i\partial_t u-Hu$7

with analogous estimates in $i\partial_t u-Hu$8 (Carles et al., 2022). In the Coulomb/convolution setting, the theory covers $i\partial_t u-Hu$9 with $\mathcal M$0 and $\mathcal M$1 for some $\mathcal M$2, together with at-most-quadratic smooth parts of $\mathcal M$3 to ensure local Strichartz estimates (Carles et al., 2022).

The main results include a global existence and uniqueness theorem for the Coulomb/convolution case:

$\mathcal M$4

with

$\mathcal M$5

and conservation of the $\mathcal M$6 norms (Carles et al., 2022). For general polynomial coupling, if

$\mathcal M$7

then there exists a unique global solution

$\mathcal M$8

The total energy

$\mathcal M$9

is time-independent (Carles et al., 2022).

The proof strategy differs between the Coulomb and polynomial settings. For Coulomb/convolution couplings, the Duhamel formulation is combined with local Strichartz estimates and a fixed-point argument in a Banach space controlled by admissible $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$0 norms (Carles et al., 2022). For general polynomial couplings, the argument uses a recursive linearization inspired by the Cauchy problem for symmetric quasilinear hyperbolic equations. A linear step is solved with Kato–Rellich, high-order graph norms are estimated through commutator bounds, convergence is obtained in $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$1 by contraction for small $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$2, and globalization follows from energy coercivity (Carles et al., 2022).

3. Semiclassical propagation and coherent-state structure

A complementary line of analysis studies the time-dependent Hartree equation in $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$3 with external potential $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$4 and interaction kernel $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$5:

$\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$6

The assumptions are $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$7 and $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$8, with $\begin{cases} i\partial_t\phi^x = H_x\,\phi^x \;+\;\displaystyle\int w(x,y)\,|\phi^y(t,y)|^2\,dy\;\phi^x,\[1ex] i\partial_t\phi^y = H_y\,\phi^y \;+\;\displaystyle\int w(x,y)\,|\phi^x(t,x)|^2\,dx\;\phi^y, \end{cases}$9 real-valued, even, and all derivatives up to order three bounded in $$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$0 (Athanassoulis et al., 2010).

For coherent-state initial data centered at $$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$1,

$$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$2

the solution admits a semiclassical approximation by a time-dependent coherent state whose center follows a classical Hamiltonian flow (Athanassoulis et al., 2010). The effective classical Hamiltonian is

$$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$3

and the phase contains both the classical action and the Hartree contribution (Athanassoulis et al., 2010).

A central quantitative result is the $$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$4 error estimate

$$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$5

The derivation uses subtraction of the semiclassical ansatz from the exact solution, a source estimate of order $$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$6 arising from Taylor remainders of $$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$7 and $$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$8, unitarity of the linear propagator, and Grönwall’s inequality (Athanassoulis et al., 2010).

The same analysis yields a full formal asymptotic expansion

$$H_x=-\tfrac12\Delta_x+V_1(x),\qquad H_y=-\tfrac12\Delta_y+V_2(y).$$9

with recursively determined coefficients. The leading profile solves a “harmonic-Hartree” transport equation involving $N$0 and $N$1, while the next term solves a linear inhomogeneous Schrödinger equation driven by $N$2, $N$3, and the lower-order profile (Athanassoulis et al., 2010). This places the time-dependent Hartree equation within a semiclassical framework in which coherent states remain accurate up to an $N$4 error.

4. Mean-field derivation for dense fermions

In a distinct asymptotic regime, the time-dependent Hartree equations emerge as the effective large-$N$5 dynamics of dense interacting fermions. The microscopic system considered in (Hoang et al., 16 Jul 2025) consists of spinless fermions in $N$6 with antisymmetric wave function

$N$7

evolving under

$N$8

where $N$9 is real-valued, radial, and $$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$0-independent (Hoang et al., 16 Jul 2025).

Because the fermions are confined to a volume of order one, typical momenta are $$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$1, the kinetic energy is $$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$2, and the interaction energy is $$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$3 (Hoang et al., 16 Jul 2025). The nontrivial mean-field time scale is obtained by

$$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$4

On this rescaled time scale, the limiting dynamics is the $$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$5-orbital Hartree system

$$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$6

with orthonormal initial orbitals $$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$7 (Hoang et al., 16 Jul 2025).

A central technical step is a time-dependent gauge transformation

$$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$8

The gauged state satisfies

$$i\partial_t\psi_j(t,x) = -\tfrac12\Delta_x\psi_j(t,x) +\bigl(V*\rho(t,\cdot)\bigr)(x)\,\psi_j(t,x),\qquad j=1,\dots,N,$$9

with

$$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$0

The corresponding gauged Hartree orbitals

$$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$1

solve a one-body gauged equation with Hamiltonian $$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$2 expressed in terms of $$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$3 and related partial traces (Hoang et al., 16 Jul 2025). The paper’s interpretation is explicit: the gauge transformation eliminates the dominant interaction contribution and replaces it by more regular “magnetic-type” terms.

The proof introduces a counting functional for excitations outside the Slater space, an auxiliary dynamics generated by a quadratic truncation of the gauged Hamiltonian, propagation-of-regularity bounds for the gauged Hartree orbitals, and a weighted Duhamel estimate comparing the full and auxiliary evolutions (Hoang et al., 16 Jul 2025). The final estimate is a trace-norm approximation of the one-particle reduced density matrix:

$$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$4

This gives a rigorous derivation of the time-dependent Hartree equations in a strongly interacting dense fermion regime (Hoang et al., 16 Jul 2025).

A common misconception is that Hartree mean-field limits are available only in weak-coupling or semiclassical scalings. The dense-fermion result does not fit that template: the interaction is not scaled by $$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$5, and the proof instead relies on the gauge-transformed formulation (Hoang et al., 16 Jul 2025). This suggests that the effective Hartree description can persist in regimes where the raw two-body interaction is parametrically larger than in standard mean-field analyses.

5. Relation to Hartree–Fock, relativistic Hartree, and pairing theories

The time-dependent Hartree equation is closely related to, but not identical with, time-dependent Hartree–Fock. In the electronic TDHF formulation, the many-electron state is approximated by a single Slater determinant of time-dependent spin-orbitals,

$$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$6

and each orbital satisfies

$$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$7

with mean-field Hamiltonian $$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$8 (Gulania et al., 2023). In a finite basis,

$$\rho(t,x)=\sum_{k=1}^N|\psi_k(t,x)|^2,\qquad (V*\rho)(x)=\int_{\Bbb R^d}V(x-y)\,\rho(y)\,dy.$$9

and the Fock matrix is

$V(x)=\pm|x|^{-1}$0

The exchange contribution, encoded by the second term in brackets, is absent in the Hartree model (Gulania et al., 2023).

This difference appears especially clearly in the relativistic framework. In the time-dependent relativistic Hartree–Fock model with spherical symmetry, the single-particle Dirac equation is

$V(x)=\pm|x|^{-1}$1

with

$V(x)=\pm|x|^{-1}$2

where $V(x)=\pm|x|^{-1}$3 is the local direct Hartree mean field and $V(x)=\pm|x|^{-1}$4 is the nonlocal exchange term (Geng et al., 2024). If one omits all Fock terms, $V(x)=\pm|x|^{-1}$5, then

$V(x)=\pm|x|^{-1}$6

becomes “exactly the time-dependent relativistic Hartree (TD-RMF) equation” (Geng et al., 2024).

Pairing theories provide another reduction. In canonical-basis TDHFB, the generalized density matrix obeys

$V(x)=\pm|x|^{-1}$7

with normal density $V(x)=\pm|x|^{-1}$8 and pairing tensor $V(x)=\pm|x|^{-1}$9 (Ebata et al., 2010). Under the diagonal approximation for the pair potential in the canonical basis, the equations simplify to coupled evolution equations for canonical orbitals, occupations, and pair amplitudes. If $d=3$0 for all $d=3$1, then $d=3$2 freezes at $d=3$3 or $d=3$4, $d=3$5, and the single-particle states satisfy

$d=3$6

which is “precisely the usual TDHF (or in a local-density approximation, the time-dependent Hartree) equation”

$d=3$7

(Ebata et al., 2010).

The algorithmic literature reflects the same distinction. A hybrid quantum–classical TDHF implementation can propagate the one-body mean-field Hamiltonian

$d=3$8

with Trotterization, matchgate rotations, and Yang–Baxter-equation-based circuit compression to constant depth per time step (Gulania et al., 2023). The paper explicitly notes that TDHF is free-fermionic and classically simulatable in polynomial time, so the quantum formulation is positioned as a benchmark for real-time mean-field dynamics rather than as a reinterpretation of Hartree itself (Gulania et al., 2023).

6. Multiconfigurational generalizations and the Hartree limit

The time-dependent Hartree equations occupy a limiting position within multiconfigurational theories. In the standard MCTDH ansatz for an $d=3$9-mode wavefunction,

$i\partial_t u-Hu$00

the Dirac–Frenkel principle yields coupled equations for the coefficient tensor $i\partial_t u-Hu$01 and the single-particle functions (SPFs) (Bonfanti et al., 2018). The SPF equations involve the reduced single-particle density matrices and their inverses,

$i\partial_t u-Hu$02

which constitute a standard numerical bottleneck (Bonfanti et al., 2018).

If one chooses $i\partial_t u-Hu$03 for every mode, then there is only one SPF per mode and $i\partial_t u-Hu$04. The ansatz collapses to

$i\partial_t u-Hu$05

and the MCTDH equations reduce to the familiar time-dependent Hartree equations

$i\partial_t u-Hu$06

(Bonfanti et al., 2018). In this sense, MCTDH generalizes TDH by replacing a product-state ansatz with a time-dependent multiconfiguration expansion.

For indistinguishable particles, the multiconfigurational extension is MCTDH-X, comprising MCTDH-B for bosons and MCTDH-F for fermions (Lode et al., 2019). The wavefunction is expanded in time-dependent permanents or Slater determinants built from optimized orbitals:

$i\partial_t u-Hu$07

The coefficient equation is linear in $i\partial_t u-Hu$08,

$i\partial_t u-Hu$09

while the orbital equation is nonlinear because it depends on the one- and two-body reduced density matrices,

$i\partial_t u-Hu$10

(Lode et al., 2019). Restricted-active-space variants further truncate the CI space and introduce orbital-rotation terms between active subspaces (Lode et al., 2019).

Recent numerical analysis revisits MCTDH from the perspective of tangent-space projectors and projector splitting. The projector-splitting algorithm replaces explicit reduced-density inversion by an auxiliary set of non-orthogonal SPFs and QR or Cholesky factorizations of matricized coefficient tensors (Bonfanti et al., 2018). A closely related tensor-network extension, MPS-MCTDH, rewrites the coefficient tensor in matrix-product-state form and propagates local site tensors with effective Hamiltonians constructed by DMRG-like sweeps (Kurashige, 2018). The mean-field operators entering the SPF equations are then evaluated by contractions over the MPS representation (Kurashige, 2018).

The multiconfigurational literature therefore clarifies the position of TDH. It is not a rival to MCTDH-type methods but their lowest-rank product-state limit. A plausible implication is that many numerical innovations developed for MCTDH—projector splitting, tangent-space formulations, and tensor-network parametrizations—can be interpreted as systematic strategies for moving beyond the Hartree manifold while retaining variational time propagation (Bonfanti et al., 2018, Kurashige, 2018).

7. Applications and scope across physical settings

The application domain of time-dependent Hartree-type equations is broad, but the precise physical role depends on the variant under consideration. In the nonrelativistic product-state setting, the Hartree system is a model of dimension reduction for quantum dynamics and a leading-order mean-field approximation whose mathematical well-posedness supports subsequent approximation and error analyses (Carles et al., 2022). In the semiclassical regime, it describes the propagation of coherent states under self-consistent nonlinear potentials and admits explicit control of the departure from a classical-flow-based wave-packet ansatz (Athanassoulis et al., 2010).

In relativistic nuclear dynamics, the Hartree reduction appears when exchange is omitted from the time-dependent relativistic Hartree–Fock equation. The resulting TD-RMF system retains only the local Hartree self-energies $i\partial_t u-Hu$11 and $i\partial_t u-Hu$12 generated by the instantaneous scalar density and vector current (Geng et al., 2024). The spherical TD-RHF framework was then used to study the isoscalar giant monopole resonance of $i\partial_t u-Hu$13Pb, with numerical checks of conservation of total binding energy, particle number, and time-reversal invariance (Geng et al., 2024). These are properties of the RHF formulation, but the Hartree reduction is the direct-field limit of that framework.

In multiconfigurational many-body physics, the Hartree limit serves as a reference point against which correlation effects are measured. MCTDH-X benchmarks against the Harmonic Interaction Model and its time-dependent version show convergence to exact solutions as the number of orbitals $i\partial_t u-Hu$14 increases, and bosonic as well as fermionic applications include Bose-condensate dynamics, photoionization cross sections, time delays, high-harmonic generation, and transient absorption (Lode et al., 2019). The specific claim here is not that the Hartree equations alone capture those correlated phenomena, but that the time-dependent Hartree picture is embedded within a hierarchy of increasingly correlated variational theories (Lode et al., 2019).

Across these settings, the central structural feature remains unchanged: the equation is nonlinear because the potential depends on the evolving density, yet it remains one-body at the level of each orbital or factor. The literature reviewed here presents that structure in several mathematically distinct forms—product-state Dirac–Frenkel projection (Carles et al., 2022), semiclassical coherent-state propagation (Athanassoulis et al., 2010), dense-fermion mean-field derivation (Hoang et al., 16 Jul 2025), relativistic direct-field reduction (Geng et al., 2024), and multiconfigurational product-state limit (Bonfanti et al., 2018).