---
title: Time-Dependent Variational Principle (TDVP)
url: https://www.emergentmind.com/topics/time-dependent-variational-principle-tdvp
type: topic
---

# Time-Dependent Variational Principle (TDVP)

The time-dependent variational principle (TDVP) is a foundational method for simulating quantum dynamics in high-dimensional and strongly correlated systems by restricting time evolution to a chosen variational manifold. Originally formulated for wavefunction dynamics, the TDVP provides a mathematically rigorous path to optimal, symmetry-preserving approximations for time evolution and yields diverse algorithmic frameworks, most notably for matrix product states (MPS), tensor networks, and Gaussian state spaces [1103.0936][1408.5056][1908.03090][1912.01125][1907.04837][2010.15725]. The TDVP has been extended to dissipative dynamics, mixed states, and open quantum systems, with implementations ranging from analytic Gaussian ansätze to neural-network-based representations of quantum states [2104.00013][1206.3102]. This article details its variational architecture, tangent-space projection, algorithmic forms, error structure, application breadth, and recent methodological advances.

## 1. Variational Principle and Projected Dynamics

The TDVP is a variational projection of exact quantum dynamics onto a tractable submanifold of states, defined by parameters $\theta = (\theta_1, \ldots, \theta_N)$. For Schrödinger evolution,
\[
i\,\partial_t\,|\Psi(\theta)\rangle = \hat H\,|\Psi(\theta)\rangle,
\]
TDVP demands that $|\Psi(\theta(t))\rangle$ evolve such that its time derivative is the orthogonal projection of $\hat H\,|\Psi(\theta)\rangle$ onto the manifold’s tangent space:
\[
i\,P_{T_\Psi}\,\partial_t|\Psi(\theta)\rangle = P_{T_\Psi}\,\hat H\,|\Psi(\theta)\rangle,
\]
where $P_{T_\Psi}$ is the projector onto the tangent space at $|\Psi(\theta)\rangle$ [1103.0936].

This is variationally justified via Dirac-Frenkel stationarity of the action
\[
S[\Psi(\theta)] = \int dt\,\langle\Psi(\theta)|i \partial_t - \hat H|\Psi(\theta)\rangle,
\]
leading to Euler–Lagrange-type equations
\[
\langle\partial_j\Psi|\partial_i\Psi\rangle\,\dot{\theta}^i = -\langle\partial_j\Psi|\hat H| \Psi\rangle,
\]
with $\partial_i = \frac{\partial}{\partial \theta^i}$ and the Gram matrix $G_{ji} = \langle\partial_j\Psi|\partial_i\Psi\rangle$. After gauge fixing, $G$ becomes invertible and the flow reduces to a system of explicit ODEs [1103.0936][1408.5056].

In the McLachlan variant, one minimizes the norm $\|i\,\partial_t|\Psi\rangle - \hat H|\Psi\rangle\|$ over tangent vectors, resulting in equivalent projective equations [1912.01125].

## 2. Tangent Space, Gauge, and Projector Construction

For a variational manifold parameterized by tensors or wavefunction coefficients, the tangent space $T_\Psi$ consists of all infinitesimal variations generated by the parameters. The orthogonal projector can be written as
\[
P_{T} = \sum_{i,j} |\partial_i \Psi\rangle\, (G^{-1})^{ij} \langle\partial_j \Psi|,
\]
where $G$ is the metric/Grothendieck (Gram) matrix. This formalism ensures that only physical, gauge-fixed tangent vectors contribute to the dynamics. When redundancies due to parameterization (such as MPS gauge freedom) exist, gauge fixing—e.g., by canonical forms—renders $G$ invertible and the equations numerically stable [1103.0936][1408.5056][1908.03090].

For MPS, the canonical site decomposition and tangent-space construction permit an explicit splitting of the projector into sums of one-site and bond projectors [1103.0936][1408.5056], providing the algebraic backbone for practical local-update algorithms.

## 3. Algorithmic Realizations: MPS, Tensor Networks, Gaussians

TDVP admits algorithmic implementations across various variational classes:

- **Matrix Product States (MPS):** For 1D quantum chains, the TDVP is most often implemented on the MPS manifold, either in uniform (infinite) or finite-size settings. The dynamics is realized by sweeping local (one-site or two-site) updates, in which each tensor (or pair) is evolved under an effective Hamiltonian constructed from the network’s environment. Imaginary-time TDVP delivers efficient ground-state optimization resembling DMRG, while real-time TDVP is symplectic, energy-conserving, and numerically stable [1103.0936][1408.5056].

- **Tree Tensor Networks (TTNs):** For arbitrary loop-free tensor networks, the TDVP constructs tangent projectors respecting the network’s gauge and locality. Updates take the form of single-site or multi-site evolutions interleaved with canonicalization sweeps. The flexibility of TTN geometry allows application to impurity problems and hybridization with Matrix Product Operators (MPOs) [1908.03090].

- **Gaussian States and Bosonic/Fermionic Systems:** TDVP on Gaussian (coherent/squeezed) manifolds generates coupled ODEs for displacements and covariances, capturing both mean-field and quantum fluctuation effects. The formalism provides accurate dynamics for field theories, bosonic lattice systems, and optimal approximations for ground states and linear response [1912.01125][1907.04837][2010.15725].

## 4. Error Structure, Conservation Laws, and Complexity

TDVP evolution at every step is optimal within the chosen manifold: the residual $i\partial_t|\Psi\rangle - \hat H|\Psi\rangle$ is minimized in norm or orthogonally projected out. The only source of deviation from exact dynamics is the geometric (variational) error—how much $\hat H|\Psi\rangle$ fails to lie in $T_\Psi$—which can be monitored at each time step. In contrast to Trotter-based algorithms, there is no Trotter error; time step errors enter only through the ODE solution [1103.0936][1408.5056].

Crucially, TDVP exactly conserves norm and energy (for time-independent Hamiltonians) and any symmetry generator commuting with $\hat H$, a property absent in TEBD or split-step integrators [1103.0936][1809.01400]. The computational complexity for each time step scales as $O(D^3)$ for MPS of bond dimension $D$, matching TEBD but often with substantially reduced truncation-induced error, especially in real-time evolution [1103.0936][1408.5056][2508.10096]. For Gaussian manifolds, the ODEs are of size $O(N^2)$ for $N$-mode systems [1907.04837].

## 5. Extensions: Dissipative, Mixed-State, and Open Quantum Dynamics

TDVP generalizes to mixed-state and dissipative dynamics governed by Lindblad equations. In these settings, the equation of motion for the density matrix $\rho$,
\[
\dot{\rho} = \mathcal{L}[\rho],
\]
is projected onto the tangent space of a variational manifold of mixed states. For pure-state ansatzes (e.g., MPS with statistical sampling), a Fokker–Planck or stochastic differential equation is derived for the variational parameters, with the drift and diffusion determined by projected Lindbladian action [1411.5546].

In the general mixed-state setting, the projection depends on a choice of monotone Riemannian metric (e.g., quantum Fisher, Bures–Helstrom), inducing a family of TDVP flows. For fermionic Gaussian states, all such metrics coincide and the TDVP reduces to “Gaussification”—preserving Gaussianity under dissipative evolution [1206.3102]. Algorithms extend to neural-network ansatzes for open quantum systems by recasting the density operator as a POVM-based probability distribution parameterized by deep autoregressive networks and using TDVP to project the Lindblad flow, with stochastic sampling for force and metric evaluation [2104.00013].

In open system molecular dynamics, constrained TDVP variants using the square-root NOSSE formulation enforce trace conservation and recover energy conservation in the closed-system limit, eliminating basis-induced artificial dissipation [1501.02025].

## 6. Numerical Strategies and Modern Extensions

Recent advances have focused on improving numerical stability, flexibility, and applicability:

- **Controlled Bond Expansion (CBE):** One-site TDVP is supplemented with local bond-dimension increases whenever the projection error or variance bound exceeds a prescribed threshold, improving accuracy in entangling dynamics but retaining O($D^3$) scaling [2208.10972].

- **Stochastic Adaptive Bond Growth:** The stochastic adaptive 1-TDVP (SA-1TDVP) algorithm grows MPS bond dimension by sampling new singular values according to level-spacing statistics of the entanglement Hamiltonian, achieving near-2TDVP accuracy at a fraction of the computational cost [2110.12703].

- **Krylov-Enriched TDVP:** Ancilla-Krylov-enriched TDVP augments the tangent space dynamically with global Krylov vectors, dramatically reducing projection error and allowing larger time steps without sacrificing unitarity [2005.06104].

- **Quantum Circuit Simulation:** In circuit contexts, TDVP has been tailored for quantum circuits by “local” projector splitting, efficiently handling long-range gates and entanglement diffusion, and outperforming TEBD for large-scale circuits by spreading bond growth globally [2508.10096].

- **Clifford and Neural Network Augmentations:** Clifford-dressed and Clifford-circuit-augmented TDVP interleaves periodically applied entanglement-cooling Clifford gates with TDVP sweeps, transferring stabilizer entanglement to the Clifford layer and reducing the necessary MPS bond dimension for accurate long-time simulation. Neural network-based TDVP attaches variational learning to density-matrix propagation using neural architectures [2407.01692][2407.03202][2104.00013].

## 7. Applications, Performance, and Scope

TDVP underpins state-of-the-art simulations across many areas:

- **Quantum Lattice Dynamics:** TDVP-MPS is optimal and stable for global quantum quenches, transport, and real-time correlation functions in 1D and quasi-1D systems, including those with long-range or time-dependent Hamiltonians [1103.0936][1408.5056][2208.10972].
  
- **Ground State and Imaginary-Time Optimization:** Imaginary-time TDVP unifies time evolution and ground-state search, effectively recovering the DMRG algorithm in the infinite-time-step limit [1408.5056].

- **Non-equilibrium and Finite-Temperature Dynamics:** Gaussian TDVP frameworks provide rigorous access to spectral, linear response, and thermal correlations in Bose-Hubbard, anharmonic lattice, and related bosonic systems [1907.04837][2010.15725].

- **Open Quantum Systems:** Lindblad-TDVP variants capture dissipative Many-Body dynamics far beyond exact diagonalization sizes [2104.00013]. For Gaussian states, TDVP provides unique, geometry-independent evolution [1206.3102].

- **Tensor Network Generalization:** TDVP extends seamlessly to tree tensor networks (TTN), tree-like tensor structures for impurity problems, and quantum chemistry applications [1908.03090].

- **Quantum Circuits and Benchmarking:** Quantum circuit simulation using local-TDVP achieves new limits on 49-qubit circuit depth and complexity, outperforming TEBD especially for hardware-efficient and long-range circuits [2508.10096].

Overall, the time-dependent variational principle offers a robust, systematically optimal, and symmetry-respecting framework for quantum dynamics and optimization, adaptable to a wide range of variational manifolds and system classes, with broad impact on condensed matter, quantum information, and quantum chemistry [1103.0936][1408.5056][2104.00013][2208.10972].

Source: https://www.emergentmind.com/topics/time-dependent-variational-principle-tdvp