---
title: Dissipaton Equations of Motion (DEOM)
url: https://www.emergentmind.com/topics/dissipaton-equations-of-motion-deom
type: topic
---

# Dissipaton Equations of Motion (DEOM)

Searching arXiv for recent and foundational DEOM papers to ground the article.
Dissipaton Equations of Motion (DEOM) is an exact, nonperturbative framework for open quantum systems in which a Gaussian environment is represented in terms of statistically independent quasi-particles called dissipatons. In this formulation, the bath is not treated merely as an eliminated influence or a hierarchy of auxiliary mathematical objects; rather, its fluctuations, hybrid bath observables, and system–bath entanglement are encoded in dissipaton density operators (DDOs), yielding a hierarchy that is formally equivalent to HEOM for reduced dynamics while providing a physically interpretable many-particle description of the environment [1912.00794]. Subsequent developments established phase-space dissipaton algebra for bath coordinates and conjugate momenta, imaginary-time and nonequilibrium thermodynamic constructions, dissipaton-embedded master equations, and extensions to quadratic and other nonlinear environment couplings, thereby broadening DEOM from Gaussian linear-coupling dynamics to phase-space, thermodynamic, transport, spectroscopy, impurity, and selected nonlinear settings [2211.08853].

## 1. Origins and conceptual basis

DEOM was introduced to address the standard difficulty of open quantum dynamics: the full system-plus-bath Hamiltonian is exact but usually too large to solve directly, whereas approximate reduced-density-matrix methods often rely on weak coupling, Markovianity, or comparable simplifications [1912.00794]. The framework retains the exactness of HEOM for Gaussian baths, but differs in physical interpretation: HEOM auxiliary operators are described as largely mathematical, while DEOM assigns explicit quasi-particle meaning to the bath degrees of freedom through dissipatons [2303.10666].

The standard setting is a Hamiltonian of the form
\[
H_{\rm tot}=H+h_{\rm bath}+H_{\rm sb},\qquad H_{\rm sb}=\hat Q \hat F,
\]
or, equivalently in simpler notation,
\[
H = H_S + h_B + \hat Q \hat x,
\]
with linear system–bath coupling and a Gaussian or Wick-like environment characterized by its two-point correlation function and spectral density [2303.10666; 1912.00794]. The equilibrium bath correlation is written as
\[
\langle \hat F(t)\hat F(0)\rangle
=
\frac{1}{\pi}\int_{-\infty}^{\infty} d\omega\,
\frac{e^{-i\omega t}J(\omega)}{1-e^{-\beta\omega}},
\]
and DEOM begins by expressing this correlation as a finite sum of exponentials,
\[
\langle \hat F(t)\hat F(0)\rangle=\sum_{k=1}^K \eta_k e^{-\gamma_k t},
\qquad
\langle \hat F(0)\hat F(t)\rangle=\sum_{k=1}^K \eta_{\bar k}^{*} e^{-\gamma_k t},
\]
with \(\gamma_{\bar k}=\gamma_k^*\) when complex-conjugate pairs occur [2303.10666].

Each exponential component is then assigned to a dissipaton, giving the decomposition
\[
\hat F=\sum_{k=1}^K \hat f_k,
\]
or, in solvation-coordinate notation,
\[
\hat x = \sum_{k=1}^K \hat f_k.
\]
The dissipatons satisfy pairwise correlations
\[
\hat f_k(t)\hat f_{k'}(0)=\delta_{kk'}\eta_k e^{-\gamma_k t},
\qquad
\hat f_{k'}(0)\hat f_k(t)=\delta_{kk'}\eta_{\bar k}^{*} e^{-\gamma_k t},
\]
so each dissipaton behaves as an independent stochastic or Brownian-like mode with decay rate \(\gamma_k\) and amplitude \(\eta_k\) [1912.00794; 2303.10666]. This quasi-particle picture is a defining conceptual feature of DEOM and underlies later interpretations of hybrid bath modes, bath moments, current operators, and thermodynamic observables [2211.08853].

## 2. Dissipaton algebra, DDOs, and the hierarchy

The fundamental dynamical variables of DEOM are the dissipaton density operators,
\[
\rho^{(n)}_{\mathbf n}(t)
=
{\rm tr}_B\!\left[
\big(\hat f_K^{n_K}\cdots \hat f_1^{n_1}\big)^\circ
\rho_{\rm tot}(t)
\right],
\]
with \(\mathbf n=(n_1,\dots,n_K)\) and \(n=\sum_k n_k\); the zeroth-tier member \(\rho_{\mathbf 0}^{(0)}(t)\) is the reduced density operator [1912.00794; 2303.10666]. In bosonic settings, \((\cdots)^\circ\) denotes an irreducible or symmetrized product; in fermionic settings it becomes antisymmetric [2303.10666; 2304.08259].

Closure of the hierarchy follows from the dissipaton algebra. One ingredient is the generalized diffusion relation,
\[
{\rm tr}\!\left[(\partial_t \hat f_k)\rho(t)\right]
=
-\gamma_k\,{\rm tr}\!\left[\hat f_k\rho(t)\right],
\]
which identifies dissipatons as generalized Brownian particles whose amplitudes relax diffusively [2303.10666]. The other ingredient is the generalized Wick theorem. For bosonic dissipatons, one representative relation is
\[
\rho_{\mathbf n}^{(n)}(t;\hat f_k^\times)
=
\rho_{\mathbf n_k^+}^{(n+1)}(t)
+n_k\,\eta_k\,\rho_{\mathbf n_k^-}^{(n-1)}(t),
\]
with a corresponding right-action rule involving \(\eta_{\bar k}^*\) [2303.10666]. These operator identities permit exact tier-raising and tier-lowering actions of bath operators on DDOs.

The resulting DEOM hierarchy for Gaussian linear coupling is
\[
\dot\rho^{(n)}_{\mathbf n}
=
-\left(i{\cal L}+\sum_k n_k\gamma_k\right)\rho^{(n)}_{\mathbf n}
-i\sum_k \hat Q^\times \rho^{(n+1)}_{\mathbf n_k^+}
-i\sum_k n_k\left(\eta_k \hat Q-\eta_{\bar k}^{*}\hat Q\right)\rho^{(n-1)}_{\mathbf n_k^-},
\]
or, in equivalent notations, the same tier structure with model-specific superoperators \({\cal A}\) and \({\cal C}\) [1912.00794; 2303.10666]. For fermionic reservoirs the hierarchy acquires Grassmannian sign structure, for example
\[
\dot\rho^{(n)}_{j_1\cdots j_n}
=
-\Big(i{\cal L}+\sum_{r=1}^n\gamma_{j_r}\Big)\rho^{(n)}_{j_1\cdots j_n}
-i\sum_j {\cal A}_{\bar j}\rho^{(n+1)}_{j_1\cdots j_n j}
-i\sum_{r=1}^n(-)^{n-r}{\cal C}_{j_r}\rho^{(n-1)}_{j_1\cdots j_{r-1}j_{r+1}\cdots j_n},
\]
as used for electronic impurity problems [1812.10007].

The formal equivalence between DEOM and HEOM for reduced dynamics is emphasized repeatedly, but DEOM differs by explicitly retaining bath quasi-particle content, thereby allowing direct access to hybrid system–bath observables, bath moments, and system–bath correlations [2303.10666; 1609.03685]. This also underlies the statement, made in multiple later works, that DEOM is not merely a solver for the reduced density matrix but a dissipaton-space mechanics of entangled system–environment dynamics [2211.08853].

## 3. Phase-space formulation and bath observables

A major extension of DEOM was the incorporation of the bath conjugate momentum, producing a phase-space dissipaton theory [1912.00794]. Earlier DEOM formulations treated the bath coordinate \(\hat x\) and hybrid bath solvation coordinate dynamics, but not the conjugate momentum \(\hat p\). The phase-space construction introduces
\[
\hat p = \sum_{k=1}^K \hat\varphi_k,
\]
where \(\hat\varphi_k\) is the momentum counterpart of dissipaton \(\hat f_k\) [1912.00794].

The central momentum-action formula is
\[
\rho^{(n)}_{\mathbf n}(t;\hat\varphi_k)
=
-\frac{\gamma_k}{\omega}
\Big[
\rho^{(n+1)}_{\mathbf n_k^+}(t)
-
n_k\eta_k \rho^{(n-1)}_{\mathbf n_k^-}(t)
\Big],
\]
with the conjugate form
\[
\rho^{(n)}_{\mathbf n}(t;\hat\varphi_k)
=
-\frac{\gamma_k}{\omega}
\Big[
\rho^{(n+1)}_{\mathbf n_k^+}(t)
-
n_k\eta_{\bar k}^* \rho^{(n-1)}_{\mathbf n_k^-}(t)
\Big],
\]
and \(\rho^{(n)}_{\mathbf n}(t;\hat p)=\sum_k \rho^{(n)}_{\mathbf n}(t;\hat\varphi_k)\) [1912.00794]. The construction was validated against two criteria. The necessary criterion is recovery of the canonical commutation relation,
\[
\rho^{(n)}_{\mathbf n}\!\left(t;[\hat x,\hat p]\right)= i\,\rho^{(n)}_{\mathbf n}(t),
\]
and the sufficient criterion is the operator-action identity
\[
\dot{\hat A}(t)\hat x(0) = -\omega\,\hat A(t)\hat p(0)
\]
for arbitrary \(\hat A\) [1912.00794].

The explicit commutator algebra confirms the canonical structure. For example,
\[
\rho^{(n)}_{\mathbf n}\!\left(t;[\hat f_k,\hat\varphi_{k'}]\right)
=
2\delta_{kk'}\frac{\gamma_k\eta_k}{\omega}\rho^{(n)}_{\mathbf n}(t),
\]
and summing over indices with
\[
\omega = -2i\sum_k \gamma_k \eta_k = 2i\sum_k \gamma_k \eta_{\bar k}^*
\]
recovers
\[
\rho^{(n)}_{\mathbf n}\!\left(t;[\hat x,\hat p]\right)= i\,\rho^{(n)}_{\mathbf n}(t)
\]
exactly in DEOM space [1912.00794].

This phase-space extension is practically important because current and transport observables often depend directly on bath momentum. In the heat-current example,
\[
\hat J \equiv -\frac{d h_B}{dt} = -i[H,h_B] = \hat Q\,\omega\,\hat p,
\]
so evaluation of the equilibrium current autocorrelation
\[
C(t)\equiv \delta \hat J(t)\,\delta \hat J(0)
\]
requires the momentum algebra [1912.00794]. More broadly, the Perspective on dissipaton theories explicitly describes DEOM as a phase-space, quasi-particle reformulation that supports Schrödinger-picture, Heisenberg-picture, and imaginary-time constructions in a unified setting [2211.08853].

A related development is the dissipaton-embedded quantum master equation (DQME), which maps each dissipaton to a continuous real variable \(x_k\) through
\[
\hat f_k/\zeta_k \longleftrightarrow x_k
\]
and encodes the hierarchy in a single distribution
\[
\hat\rho(\mathbf x,t)
=
\sum_{\mathbf n} \rho_{\mathbf n}^{(n)}(t)\prod_k \phi_{n_k}(x_k).
\]
The resulting DQME,
\[
\frac{\partial}{\partial t}\hat\rho(\mathbf x,t)
=
-i[H,\hat\rho]
+\sum_k \hat\Gamma_k \hat\rho
-i\sum_k \zeta_k[\hat Q,x_k\hat\rho]
-\sum_k \xi_k \left\{\hat Q,\frac{\partial \hat\rho}{\partial x_k}\right\},
\]
with
\[
\hat\Gamma_k = \gamma_k \frac{\partial}{\partial x_k} \left(\frac{\partial}{\partial x_k}+x_k\right),
\]
gives a Smoluchowski-type interpretation of dissipatons as generalized Brownian particles [2303.10666]. This suggests that DEOM can be regarded either as a hierarchy of DDOs or as a continuous dissipaton-coordinate dynamics, with exact equivalence at the level of the chosen bath decomposition.

## 4. Imaginary-time DEOM and thermodynamic formulations

DEOM was subsequently generalized from real-time dynamics to equilibrium and transient thermodynamics. One development introduced an imaginary-time DEOM via analytical continuation \(t\to -i\tau\), yielding equilibrium DDOs and the hybridization partition function [2008.04087]. In this construction, the imaginary-time density operator
\[
\hat\varrho(\tau)
=
\frac{e^{-\tau \hat H_{\rm tot}}e^{-(\beta-\tau)(\hat H+\hat h)}}{Z_0}
\]
is projected into dissipaton space,
\[
\varrho_{\mathbf n}^{(n)}(\tau)
=
\operatorname{tr}_{\rm bath}\!\left[
\big(f_K^{n_K}\cdots f_1^{n_1}\big)^{\circ}\hat\varrho(\tau)
\right],
\]
and propagation to \(\tau=\beta\) yields equilibrium DDOs and the hybridization partition function \(Z_{\rm hyb}\) [2008.04087].

This enables direct evaluation of the hybridization Helmholtz free energy,
\[
A_{\rm hyb}(T)=A(T)-A_0(T),
\qquad
Z_{\rm hyb}=e^{-\beta A_{\rm hyb}},
\]
for isothermal mixing of a local impurity and a Gaussian bath [2008.04087]. A thermodynamic integration formalism is also used, with \(\lambda\)-dependent coupling turned on from \(0\) to \(1\), so that equilibrium DDOs determine the reversible work and hence the free-energy change [2008.04087; 2211.08853].

The same framework was extended to nonequilibrium thermodynamics through \(\lambda(t)\)-DEOM, permitting work-distribution calculations and numerical verification of the Jarzynski equality and Crooks relation [2211.08853]. The Perspective states that dissipaton theories were developed for studying equilibrium and nonequilibrium thermodynamic mixing processes and that the Jarzynski equality and Crooks relation are accurately reproduced numerically [2211.08853].

Within the spin–boson illustration of the unified thermodynamic DEOM, the hybridization free energy is negative and decreases monotonically, consistent with spontaneity, while hybridization internal energy and entropy show turnover behavior as functions of temperature [2008.04087]. The same work distinguishes the thermodynamic entropy from the reduced-state von Neumann entropy and reports a transient negative hybridization entropy interpreted as an indication of “solvent-cage” formation [2008.04087]. These statements concern the specific thermodynamic analysis of that paper rather than a universal property of all DEOM applications.

A plausible implication is that dissipaton-space thermodynamics is not an auxiliary add-on but a direct consequence of retaining explicit system–bath hybrid variables. That implication is consistent with later presentations that place real-time DEOM, imaginary-time DEOM, and \(\lambda\)-based thermodynamic constructions within one formal family [2211.08853].

## 5. Extensions beyond Gaussian linear coupling

A central line of development concerns nonlinear environmental backactions, especially quadratic bath couplings. A 2016 generalization formulated DEOM for a non-Gaussian environment with Hamiltonian
\[
H_{\rm tot}=H+h+\hat Q(\alpha_1 \hat x+\alpha_2 \hat x^2),
\]
arguing that, unlike the linear bath-coupling case, the influence of a non-Gaussian environment cannot be completely characterized with the linear response theory [1608.07774]. The key algebraic advance is a generalized Wick theorem with dissipaton-pair contributions, needed because \(\hat x^2\) introduces operator products beyond single contractions [1608.07774].

The resulting extended hierarchy contains not only the usual \(n\leftrightarrow n\pm1\) couplings from linear terms but also \(n\leftrightarrow n\pm2\) channels and same-tier quadratic contributions. A representative form is
\[
\begin{aligned}
\dot\rho^{(n)}_{\mathbf n}
&= -\Big(i{\cal L}_{\rm eff}+\sum_k n_k\gamma_k\Big)\rho^{(n)}_{\mathbf n}
-i2\alpha_2\sum_{kj} n_k {\cal C}_k\,\rho^{(n)}_{\mathbf n,kj}^{-+} \\
&\quad -i\alpha_2\sum_{kj}\Big[ {\cal A}\rho^{(n+2)++}_{\mathbf n,kj}
+n_k(n_j-\delta_{jk}){\cal B}_{kj}\rho^{(n-2)--}_{\mathbf n,kj} \Big] \\
&\quad -i\alpha_1\sum_k\Big( {\cal A}\rho^{(n+1)+}_{\mathbf n,k}
+n_k{\cal C}_k\rho^{(n-1)-}_{\mathbf n,k} \Big),
\end{aligned}
\]
with superoperators defined by \({\cal A}\), \({\cal B}_{kj}\), and \({\cal C}_k\) [1608.07774]. A later comprehensive account emphasized that the extended DEOM is verified by an independent Brownian-solvation-mode embedded hierarchical quantum master equation and by thermodynamic checks reproducing Jarzynski and Crooks [2206.14375].

An alternative exact route to nonlinear coupling is the stochastic-fields-dressed DEOM (SFD-DEOM), in which the nonlinear bath term is decoupled by a Hubbard–Stratonovich transformation so that only linear bath coupling remains for each stochastic realization [2108.10013]. For the quadratic-coupling Hamiltonian
\[
H=H_{\rm S}+h_{\rm B}+\hat Q(\alpha_0+\alpha_1\hat x+\alpha_2\hat x^2),
\]
the transformed stochastic Hamiltonian becomes
\[
H(\xi_t)=H_0+h_{\rm B}+Q(\xi_t)\hat x,
\qquad
Q(\xi_t)=\alpha_1\hat Q+(1+i)\xi_t\sqrt{\alpha_2}\,\hat Q^{1/2},
\]
with an analogous bra-side expression involving \(\xi_t'\) [2108.10013]. The hierarchy is then of ordinary DEOM form but with stochastic system operators, and the exact reduced density is obtained after ensemble averaging [2108.10013].

More recent extended dissipaton theories treat quadratic fermionic couplings in electronic impurity systems. One such hierarchy for quadratic fermionic bath operators contains both explicit \(n+2\) and \(n-2\) couplings and was applied to the Kondo impurity model [2304.08259]. Another presentation unified linear and quadratic fermionic couplings in an extended DEOM and equivalent dissipaton-embedded master equation, with application to an adatom–graphene composite [2409.00669]. In both cases, the extension is described as exact, non-Markovian, and non-perturbative within the bath-correlation decomposition.

These developments clarify a recurring misconception. Standard DEOM is exact for Gaussian baths with linear coupling after exponential decomposition, but nonlinear or quadratic environmental couplings require additional algebraic structure, stochastic dressing, or an extended hierarchy. The later theories are extensions of DEOM, not automatic consequences of the original linear-coupling formalism [1608.07774; 2108.10013; 2206.14375].

## 6. Applications across spectroscopy, transport, impurity physics, and anharmonic environments

DEOM has been applied across a wide range of open-system problems in which exact non-Markovian dynamics and explicit bath observables are advantageous.

In excitonic and vibronic spectroscopy, DEOM was used to study Herzberg–Teller vibronic coupling in excitation energy transfer and two-dimensional coherent spectroscopy [1609.03685]. There the transition dipole depends linearly on the bath coordinate,
\[
\mu_a(\mathbf X_a)=\mu_a+\mu_a' \hat F_a/\lambda,
\]
so optical observables explicitly involve hybrid bath dynamics [1609.03685]. The simulations show that non-Condon vibronic coupling intensifies dynamical electronic-vibrational energy transfer, enhances total system-and-bath quantum coherence, and enriches the interpretation of nonlinear spectroscopy [1609.03685].

In photosynthetic reaction-center modeling, a mixed DEOM–Lindblad strategy treated the phonon bath exactly with DEOM and the photon bath by a Lindblad dissipator [2205.04353]. The mixed scheme showed transfer coherence up to a few hundreds femtoseconds and an environmental manipulation effect on current that differs qualitatively from fully Markovian or classical-bath descriptions [2205.04353]. This application uses DEOM primarily to preserve non-Markovian phonon memory and complex bath correlations.

In quantum impurity and Kondo physics, DEOM has been used as a numerically exact open-system method for triple quantum dots, double quantum dots, and impurity models [1812.10007; 2305.17686]. In the open triple quantum dot problem, DEOM was used to analyze long-range entanglement of Kondo clouds through spectral functions, spin–spin correlations, and occupancies, with a conduction-electron-like peak on the intermediate dot interpreted as a signature of overlap between two Kondo clouds [1812.10007]. In the double-impurity Anderson model, DEOM was employed to compute both impurity spectral functions and current-noise spectra in the Kondo regime, with the Kondo signatures reshaped by inter-dot Coulomb interaction, inter-dot hopping, and chemical-potential bias [2305.17686]. In the quadratic-coupling Kondo impurity model treated by extended DEOM, a sharp Kondo peak develops near the Fermi energy at low temperature, while perturbation theory fails and can violate the Friedel sum rule [2304.08259].

In transport and current-correlation problems, DEOM has also been adapted to many-level electron–phonon dynamics. For a one-dimensional Peierls-type model with Brownian-oscillator spectral density, DEOM was used to compute the finite-temperature real-time current autocorrelation function, optical conductivity, diffusion constant, and mobility [2508.14447]. That work exploited dissipaton algebra to handle the phonon-assisted current operator and used a momentum-space reformulation to reduce computational cost [2508.14447]. The reported conclusion is that with increasing damping, DEOM optical-conductivity profiles become qualitatively similar to transient-localization predictions, suggesting that low-frequency enhancements found for undamped vibrations are artifacts of the undamped limit [2508.14447].

In spin-lattice and spin-environment problems, DEOM has been generalized in two distinct directions. One work treated a continuous bath of independent spins by mapping the spin bath, in its linear-response limit, onto an effective Gaussian environment characterized by a temperature-dependent spectral density [2302.00215]. Another used extended DEOM for spin-lattice relaxation with linear and quadratic phonon couplings and showed that exact non-Markovian rate kernels can deviate substantially from Fermi’s golden rule, especially for quadratic coupling where the kernel exhibits a free-induction-decay-like feature and coupling-dependent damping [2404.04803].

In molecular spectroscopy with structured correlated environments, ext-DEOM was used to treat correlated vibration–solvent effects and Duschinsky rotations in electron transfer and optical spectra [2503.06876]. There the system–bath coupling takes the quadratic form
\[
\hat Q\Big({\bm\alpha}_1\cdot \hat{\bm q} + \hat{\bm q}^{T}{\bm\alpha}_2 \hat{\bm q}\Big),
\]
and the hierarchy contains \(\Delta n=\pm1\) channels from linear terms and \(\Delta n=+2,0,-2\) channels from quadratic terms [2503.06876]. The spectra show that Duschinsky rotation and solvent-induced correlations can substantially reshape peak positions, widths, and intensities [2503.06876].

Finally, a dissipaton-based approach to a locally probed one-dimensional chain formulated the dynamics directly in terms of c-number dissipaton moments rather than density operators, with cross-tier recursive couplings for higher-order chain–probe interactions [2603.29458]. This suggests a further methodological branch in which the dissipaton algebra is retained but the propagated objects are ordinary moments tailored to transport observables such as heat current.

Taken together, these applications indicate that DEOM functions as a general exact hierarchy formalism for Gaussian baths and as a platform for several controlled extensions. Its distinctive strength lies not only in exact reduced dynamics but in direct access to hybrid bath coordinates, momenta, current operators, thermodynamic quantities, and correlated system–bath observables across spectroscopy, transport, impurity physics, and selected nonlinear environments [2211.08853].

Source: https://www.emergentmind.com/topics/dissipaton-equations-of-motion-deom