---
title: 'VDAT: Variational Discrete Action Theory'
url: https://www.emergentmind.com/topics/variational-discrete-action-theory-vdat
type: topic
---

# VDAT: Variational Discrete Action Theory

Variational Discrete Action Theory (VDAT) denotes, in the cited literature, a set of discrete variational frameworks that place a discrete action functional at the center of analysis and derive dynamics, control laws, integrable lattice equations, or ground-state approximations from stationarity, constrained minimization, or related closure conditions. In one major line of work, VDAT is the higher-order, constrained, structure-preserving framework obtained by discretizing Hamilton’s principle directly; in another, it is the pluri-Lagrangian multi-time theory of surface-dependent discrete actions; in another, it is an integer-time quantum many-body formalism built from the sequential product density matrix (SPD) ansatz and a discrete action in compound space; and in another, it is a variational-inference formulation of control on discrete-event decision processes [1306.0298] [1307.0523] [2011.14510] [1905.02606].

## 1. Terminological scope and shared variational structure

A common source of confusion is that the same acronym is used for non-equivalent constructions. The cited works collectively show that VDAT is not restricted to one mathematical object or application domain. This suggests that the unifying feature is methodological rather than disciplinary: the primitive quantity is discrete action, and the derived objects are update maps, corner equations, free-energy objectives, or integer-time Green’s functions, depending on context [1306.0298] [1312.1440] [2011.14509] [1905.02606].

| Usage in the literature | Primitive discrete object | Characteristic result |
|---|---|---|
| Geometric mechanics | \(S_d=\sum L_d\) | symplectic variational integrator |
| Integrable lattice systems | action on discrete surfaces | closure \(dL=0\) and surface independence |
| Discrete-event control | variational free-energy / ELBO | EP-style policy improvement |
| Quantum many-body theory | SPD with integer-time discrete action | SCDA exact in \(d=\infty\) |

In the geometric-mechanics literature, VDAT is explicitly described as “the viewpoint that discrete dynamical schemes should be derived by applying Hamilton’s principle directly at the discrete level,” with a discrete Lagrangian \(L_d\) approximating the action integral of a continuous Lagrangian \(L\). In the multiform literature, the analogous primitive object is a discrete \(d\)-form whose action is evaluated on arbitrary \(d\)-surfaces in a higher-dimensional multi-time. In the quantum many-body literature, the primitive object is the SPD together with an integer-time discrete action and corresponding Green’s-function formalism. In the discrete-event control literature, control is recast as variational inference on trajectory distributions, with rewards entering as log-potentials [1306.0298] [1307.0523] [2011.14510] [1905.02606].

## 2. Classical discrete mechanics, constraints, and structure preservation

In the mechanical formulation, a first-order discrete Lagrangian is \(L_d:Q\times Q\to\mathbb{R}\), and for a uniform time-step \(h\) the discrete action of a path \(q=(q_0,\dots,q_N)\) is
$$
S_d[q]=\sum_{k=0}^{N-1}L_d(q_k,q_{k+1};h).
$$
Stationarity with fixed endpoints yields the discrete Euler–Lagrange equations
$$
D_2L_d(q_{k-1},q_k)+D_1L_d(q_k,q_{k+1})=0,\qquad 1\le k\le N-1.
$$
Under the standard regularity assumption that \(D_{12}L_d(q_k,q_{k+1})\) is invertible, the update \((q_{k-1},q_k)\mapsto(q_k,q_{k+1})\) is well posed. The associated left and right discrete Legendre transforms,
$$
F^-L_d(q_k,q_{k+1})=(q_k,p_k^-),\qquad p_k^-=-D_1L_d(q_k,q_{k+1}),
$$
$$
F^+L_d(q_k,q_{k+1})=(q_{k+1},p_k^+),\qquad p_k^+=D_2L_d(q_k,q_{k+1}),
$$
pull back the canonical symplectic form on \(T^*Q\) to the same discrete two-form \(\omega_d\), and the discrete flow preserves \(\omega_d\) exactly. If a Lie group \(G\) acts on \(Q\) and \(L_d\) is \(G\)-invariant, the discrete momentum map \(J_d\) is conserved along the discrete flow [1306.0298].

The higher-order extension replaces \(T^{(r)}Q\) by \(Q^{r+1}\). A higher-order discrete Lagrangian is \(L_d:Q^{r+1}\to\mathbb{R}\), the discrete action is
$$
S_d[q]=\sum_{k=0}^{N-r}L_d(q_{(k:k+r)}),
$$
and stationarity with fixed boundary blocks \(q_{(0:r-1)}\) and \(q_{(N-r+1:N)}\) yields the higher-order discrete Euler–Lagrange equations
$$
D_{r+1}L_d(q_{(k:k+r)})+D_rL_d(q_{(k+1:k+1+r)})+\cdots+D_1L_d(q_{(k+r:k+2r)})=0.
$$
For higher-order discrete constraints \(\Phi_d:Q^{s+1}\to\mathbb{R}^m\), one introduces multipliers \(\lambda_k\in\mathbb{R}^m\) and the augmented action
$$
S_d^c[q,\lambda]=\sum_{k=0}^{N-r}\left[L_d(q_{(k:k+r)})+\lambda_k^T\Phi_d(q_{(k:k+s)})\right].
$$
The resulting constrained update is implicit, but under the nonsingularity condition of Proposition 3.2 it defines a unique local map \(T_d\). The discrete Poincaré–Cartan one-forms \(\theta_d^\pm\) satisfy \(d\theta_d^-=d\theta_d^+\), so \(\Omega_d:=-d\theta_d^-=-d\theta_d^+\); after restriction to the constraint manifold, \(\Omega_M\) is symplectic and \((T_d)^*\Omega_M=\Omega_M\). If both \(L_d\) and \(\Phi_d\) are \(G\)-invariant, the higher-order discrete momentum map is preserved [1306.0298].

Time dependence can be incorporated by treating time as a discrete variable:
$$
S_d[t,q]=\sum_{i=0}^{N-r}(t_{i+r}-t_i)L_d(t_{(i:i+r)},q_{(i:i+r)}).
$$
Variation with respect to \(t\) yields a discrete energy balance. If \(L_d\) is autonomous, \(D_1L_d=0\) and the discrete energy \(E_d\) is exactly preserved, yielding a symplectic and energy-momentum preserving method. The same line of development extends to stochastic Hamiltonian systems: stochastic discrete Hamiltonian variational integrators are derived from a stochastic action functional based on a type-II stochastic generating function, they are symplectic almost surely, preserve integrals of motion related to Lie group symmetries, include stochastic symplectic Runge–Kutta methods as a special case, and achieve global mean-square order \(1.0\) for scalar noise and for commutative multidimensional noise when only \(\Delta W\)-type quadrature is used [1306.0298] [1609.00463].

A parallel canonical treatment emphasizes that discrete variational systems have both pre- and post-constraint surfaces. For a nearest-neighbor action \(S=\sum_kL_k(x_k,x_{k+1})\), the pre- and post-momenta are
$$
p_k^-:=-\partial_{x_k}L_{k-1,k}(x_{k-1},x_k),\qquad
p_k^+:=\partial_{x_k}L_{k,k+1}(x_k,x_{k+1}),
$$
with momentum matching \(p_k^-=p_k^+\) on solutions. Singular Legendre maps yield pre-constraints \(C_k^-\) and post-constraints \(C_k^+\). A notable discrete peculiarity is that first-class constraints need not generate symmetries unless they coincide as both pre and post; on evolving phase spaces, the number of constraints at a fixed step depends on the initial and final step of evolution. This formulation is explicitly designed to handle constant and evolving phase spaces, lattice field theory, and discrete gravity [1303.4294].

## 3. Pluri-Lagrangian systems and discrete integrability

In the integrable-systems literature, the relevant VDAT object is a Lagrangian multiform. A \(d\)-dimensional pluri-Lagrangian problem on \(\mathbb{Z}^m\), with \(m>d\), asks for fields \(x\) such that for any oriented \(d\)-dimensional manifold \(\Sigma\), the action
$$
S_\Sigma=\int_\Sigma L
$$
is stationary with respect to all interior vertex variations. For \(d=2\), one uses a discrete Lagrangian \(2\)-form on oriented elementary squares \(\sigma_{ij}\), and the central local equations are the corner equations supported on \(3\)D corners. The discrete exterior derivative on an elementary cube is
$$
S^{ijk}=dL(\sigma_{ijk})=\Delta_kL(\sigma_{ij})+\Delta_iL(\sigma_{jk})+\Delta_jL(\sigma_{ki}),
$$
and the multi-time Euler–Lagrange equations are the eight equations \(\partial S^{ijk}/\partial(\text{corner variable})=0\). For three-point \(2\)-forms of ABS type,
$$
L(\sigma_{ij})=\mathcal{L}(x,x_i,x_j;\alpha_i,\alpha_j)
= L(x,x_i;\alpha_i)-L(x,x_j;\alpha_j)-\Lambda(x_i,x_j;\alpha_i,\alpha_j),
$$
the cube action reduces to an octahedral expression independent of \(x\) and \(x_{ijk}\), so the eight equations reduce to six nontrivial corner equations [1307.0523].

The closure relation
$$
\Delta_iL(\sigma_{jk})+\Delta_jL(\sigma_{ki})+\Delta_kL(\sigma_{ij})=0
$$
on solutions is the variational integrability criterion. It implies that the action depends only on the boundary of the surface, hence is invariant under local flips. The paper “What is integrability of discrete variational systems?” proves that for the three-point \(2\)-forms encoding the ABS list, the corner equations are consistent in the sense of minimal rank \(2\) per cube, and that the corresponding \(2\)-forms are closed not only on solutions of the underlying quad-equations but also on general solutions of the corner equations. The same work also exhibits a pluri-Lagrangian system not coming from a multidimensionally consistent system of quad-equations, showing that the pluri-Lagrangian theory goes beyond the ABS setup [1307.0523].

A closely related formulation treats the action as depending simultaneously on the field and on the choice of discrete surface. For \(2\)-dimensional systems embedded in higher-dimensional lattices, the generalized Euler–Lagrange equations include both ordinary field variations and local surface variations, and the closure relation around an elementary cube enforces surface independence. For \(3\)-dimensional systems, a Lagrangian \(3\)-form on oriented cubes satisfies an analogous hypercube closure condition. This is the multiform version of multidimensional consistency, and the paper explicitly argues that the variational principle can be considered as the defining equations for the Lagrangians themselves [1312.1440].

A further development concerns convex variational principles for discrete Laplace-type equations induced by integrable quad-equations. On the white graph of a bipartite quad-graph, the generalized discrete action
$$
S=\sum_{e=(x,y)\in E(\Lambda_W)}L(X,Y;\theta_e)-\sum_{x\in V(\Lambda_W)}\Theta_xX
$$
has Euler–Lagrange equations
$$
\sum_{e=(x,y)}\phi(X,Y;\theta_e)=\Theta_x.
$$
Under the reality conditions derived in the paper, \(S\) becomes real; for Q1–Q3 the sign or unit-circle conditions on the labels are necessary and sufficient for strict convexity or concavity, while for Q4 they are sufficient. Convexity then gives existence and uniqueness for Dirichlet boundary value problems. In the Q3 case, the functional is variationally equivalent to the Bobenko–Springborn circle-pattern functionals, relating discrete integrability, geometry, and convex action principles [1111.6273].

## 4. Optimal control, control-as-inference, and data-driven variational learning

A control-theoretic usage of VDAT appears in discrete variational optimal control. Here the forced discrete Lagrange–d’Alembert principle uses discrete forces
$$
f_k^-:Q\times Q\times U\to T^*_{q_k}Q,\qquad
f_k^+:Q\times Q\times U\to T^*_{q_{k+1}}Q,
$$
and yields the forced discrete Euler–Lagrange equations
$$
D_2L_d(q_{k-1},q_k)+D_1L_d(q_k,q_{k+1})
+f_{k-1}^+(q_{k-1},q_k,u_{k-1}^+)+f_k^-(q_k,q_{k+1},u_k^-)=0.
$$
The optimal control problem is then reformulated as a variational integrator of an augmented, higher-dimensional discrete Lagrangian system on spaces such as \(T^*Q\times T^*Q\), \(\mathfrak{g}^*\times G\times\mathfrak{g}^*\), or nonholonomic analogues. The resulting discrete necessary conditions are again Euler–Lagrange equations, and therefore inherit the preservation properties of variational integrators on configuration manifolds, Lie groups, underactuated systems, and nonholonomic systems with symmetries [1203.0580].

A distinct usage appears in discrete-event control of complex systems. There, VDAT is defined as a variational framework for optimal control in systems whose interventions are discrete actions or events. The dynamical model is a Discrete Event Decision Process (DEDP) with event intensities
$$
h_v(s_t,c_v)=c_v\prod_{m=1}^Mg_v^{(m)}(s_t^{(m)}),
$$
trajectory distribution
$$
p(\tau)=p(s_0)\prod_{t=0}^{T-1}p(a_t\mid s_t;\theta)\,p(v_t\mid s_t,a_t)\,\mathbf{1}[s_{t+1}=s_t+\Delta_{v_t}],
$$
and control-as-inference target
$$
p(\tau\mid R)\propto p(\tau)\exp\!\left(\sum_t r(s_t,a_t,v_t)\right).
$$
Approximate inference is carried out with a Bethe free-energy approximation and an EP-style forward–backward message-passing algorithm on projected singleton kernels, with overall per policy-evaluation sweep \(O(H\cdot M\cdot S^2\cdot V)\). In the reported transportation benchmarks, VDAT on SynthTown achieved TRPE \(24.21\) and EC \(75\), versus GPS \(14.79\) and EC \(100\), AC \(16.76\) and EC \(100\), and PG \(11.77\) and EC \(150\); on Berlin, VDAT achieved TRPE \(11.52\) and EC \(200\), while GPS, AC, and PG had negative rewards and failed to converge within a reasonable number of epochs [1905.02606].

A data-driven extension appears in “Variational Learning of Euler-Lagrange Dynamics from Data.” Although the authors do not use the term VDAT, their approach is explicitly characterized as relying on discrete Euler–Lagrange constraints derived from a discrete action, on variational integrators for prediction, and on variational backward error analysis. The discrete action is
$$
S_d(\{q_k\})=\sum_{k=0}^{N-1}L_d(q_k,q_{k+1};h),
$$
with midpoint or trapezoidal \(L_d\), and the learning target is an inverse modified Lagrangian \(L_{\mathrm{invmod}}\) satisfying the chosen discrete variational integrator. The method uses only position data, compensates discretization errors by inverse variational backward error analysis, and yields improved long-horizon energy behavior on the pendulum and Hénon–Heiles systems. In the pendulum experiment, the modified energy \(H^{[[2]]}\) along LSI trajectories remained within an oscillation band of width \(\sim 10^{-6}\), versus \(\sim 10^{-4}\) for \(H^{[[0]]}\) [2112.12619].

## 5. Quantum many-body VDAT: SPD, integer time, and SCDA

In the quantum many-body literature, VDAT is a variational theory for ground-state properties of quantum Hamiltonians built from the sequential product density matrix
$$
\hat{\rho}_{\mathcal N}
=\exp(\boldsymbol{\gamma}_1\cdot\spom)\,\hat{P}_1\cdots
\exp(\boldsymbol{\gamma}_{\mathcal N}\cdot\spom)\,\hat{P}_{\mathcal N},
$$
where \([\spom]_{ij}=\hat{a}_i^\dagger \hat{a}_j\), \(\boldsymbol{\gamma}_\tau\) are bilinear variational parameters, and \(\hat{P}_\tau\) are interacting projectors. The integer \(\mathcal N\) labels “integer time,” meaning the ordered sequence index of the SPD rather than physical time. The variational manifold is nested, so increasing \(\mathcal N\) monotonically improves the approximation, and \(\mathcal N\to\infty\) recovers the exact solution [2011.14510].

The hierarchy reproduces familiar ansätze at small \(\mathcal N\). The cited works state that \(\mathcal N=1\) recovers Hartree–Fock; \(\mathcal N=2\) recovers Gutzwiller, Baeriswyl, Jastrow, and unitary/variational coupled cluster; and \(\mathcal N=3\) recovers Gutzwiller–Baeriswyl and Baeriswyl–Gutzwiller. The quantum VDAT formalism then introduces an integer-time interaction representation, a compound-space discrete action, a generating function \(Z(\boldsymbol{g}_0)\), an integer-time Dyson equation,
$$
(\boldsymbol{g}^{-1}-\mathbf{1})=(\boldsymbol{g}_0^{-1}-\mathbf{1})\exp(-\boldsymbol{\Sigma})^T,
$$
and an integer-time Bethe–Salpeter structure for two-particle correlators. This construction generalizes path-integral and Green’s-function techniques to integer time and makes SPD evaluation tractable [2011.14509].

Two exact limits are central. First, for impurity-like local projectors, SPD-l can be exactly evaluated by summing a finite number of integer-time diagrams. For the one-bath-site Anderson impurity model with
$$
\hat{\rho}=e^{g\Delta\hat d}\,e^{\gamma \hat K}\,e^{g\Delta\hat d},
$$
the variational minimum yields
$$
\mathcal E=-\frac14\sqrt{64t^2+U^2},
$$
which is the exact ground-state energy for that problem. Second, for the Hubbard model in \(d=\infty\), the self-consistent canonical discrete action approximation (SCDA) is the integer-time analogue of DMFT and exactly evaluates SPD-d by assuming locality of the integer-time self-energy. In this sense, VDAT offers a variational hierarchy that is exact in \(d=\infty\) while retaining a cost profile far below continuous-time impurity solvers [2011.14509].

The earlier proposal “Variational Discrete Action Theory” emphasizes the practical consequence of this hierarchy: \(\mathcal N=2\) exactly reproduces the Gutzwiller approximation in \(d=\infty\), whereas \(\mathcal N=3\), which exactly evaluates generalized Gutzwiller–Baeriswyl states, already provides a “truly minimal yet precise description of Mott physics with a cost similar to the GA.” The same work explicitly presents SCDA update equations for \(\boldsymbol{\mathcal G}\), \(\boldsymbol{g}_{\mathrm{loc}}\), and \(\boldsymbol{S}_{\mathrm{loc}}\), making the analogy with DMFT structural rather than metaphorical [2011.14510].

## 6. Multi-orbital algorithms, qubit parametrization, and the 1RDM functional

Subsequent many-body developments concentrate on multi-orbital Hubbard models in \(d=\infty\). “Precise ground state of multi-orbital Mott systems via the variational discrete action theory” introduces a decoupled minimization algorithm for \(\mathcal N=2\)–\(4\). The paper states that \(\mathcal N=2\) rigorously recovers the multi-orbital Gutzwiller approximation, while \(\mathcal N=3\) “precisely captures the competition between the Hubbard \(U\), Hund \(J\), and crystal field \(\Delta\) in the two orbital Hubbard model over all parameter space, with a negligible computational cost.” For sufficiently large \(U/t\) and \(J/U\), \(\Delta\) drives a first-order transition within the Mott insulating regime, and in the large orbital-polarization limit with finite \(J/U\), interactions remain nontrivial even for small \(U/t\). The same work reports that \(\mathcal N=3\) can be run in seconds on a single CPU core for two orbitals [2206.08433].

The gauge-constrained \(\mathcal N=3\) algorithm sharpens SCDA by using a gauge constraint that automatically satisfies the self-consistency condition on integer-time Green’s functions. Closed-form expressions are derived for general density-density interactions, allowing straightforward application up to the seven-orbital Hubbard model. Reported single-core timings for one SPD evaluation at \(\mathcal N=3\) range from \(10^{-4}\,\mathrm{s}\) for \(N_{\mathrm{orb}}\le 3\) to \(3\,\mathrm{s}\) at \(N_{\mathrm{orb}}=7\). In the SU\((2N_{\mathrm{orb}})\) benchmarks, the \(\mathcal N=3\) critical interaction \(U_c\) closely follows the DMFT fit, whereas \(\mathcal N=2\) gives a larger \(U_c\) [2304.14616].

A later reformulation maps the local Hilbert space to a \(2N_{\mathrm{orb}}\)-qubit system with
$$
\hat n_\ell=\frac12(1-\hat\sigma_\ell^z),
$$
and rewrites the \(\mathcal N=3\) VDAT trial energy in terms of the momentum density distribution, the shape of a reference Fermi surface, and a pure qubit state. For the SU\((2N_{\mathrm{orb}})\) Hubbard model, this qubit parametrization yields an explicit procedure for computing \(U_c\) and gives the large-\(N_{\mathrm{orb}}\) asymptotic expression
$$
U_c/t\approx \frac{32}{3\pi}N_{\mathrm{orb}}+\frac{20}{3\pi}.
$$
The same paper shows that the qubit parametrization applies also to \(\mathcal N=2\), where the G-type variant yields an identical expression to slave-spin mean-field theory [2508.20111].

The same year, the constrained-search program was pushed to the level of one-body reduced density-matrix functional theory. Using the \(\mathcal N=3\) VDAT ansatz, the authors explicitly construct a 1RDM functional for the multi-orbital Hubbard model in the thermodynamic limit and report that “non-analytic behavior emerges in our 1RDMF at fixed integer filling, which gives rise to the Mott transition.” They further explain this by separating the constrained search into multiple stages and show how a nonzero Hund exchange drives the continuous Mott transition to become first-order. This use of VDAT is noteworthy because it translates the discrete-action variational hierarchy into an explicit orbital functional capable of representing Mott and Hund physics up to seven orbitals [2508.14110].

Taken together, these developments indicate a mature quantum-many-body branch of VDAT: \(\mathcal N=2\) reproduces GA; \(\mathcal N=3\) is the practical high-accuracy regime; gauge fixing and qubit parametrization reduce SCDA to explicit local optimization problems; and the same structure can be repackaged as a 1RDM functional. A plausible implication is that, within this literature, VDAT serves both as a variational hierarchy and as a constructive route from wave-function ansätze to explicit functionals.

Source: https://www.emergentmind.com/topics/variational-discrete-action-theory-vdat