---
title: Pseudomode Elimination in Open Quantum Systems
url: https://www.emergentmind.com/topics/pseudomode-elimination
type: topic
---

# Pseudomode Elimination in Open Quantum Systems

Pseudomode elimination is the operation by which non-Markovian open-system dynamics is recovered from, or encoded through, an enlarged Markovian model containing a finite set of auxiliary damped modes called pseudomodes. In the canonical Garraway-type construction, a continuum bosonic environment is replaced by discrete harmonic oscillators associated with poles of the analytically continued spectral density, each coupled to a Markovian reservoir; tracing out those pseudomodes reproduces the reduced dynamics of the original system [2108.05755]. In later work, the same expression is used more broadly for closely related reductions: replacing any eliminated sector whose influence admits a rational self-energy by damped auxiliary modes, compressing exponential bath representations to near-minimal pseudomode sets, or replacing unphysical pseudomodes by physical ensembles together with extrapolation [2605.03946][2506.10308][2311.15240].

## 1. Canonical construction from spectral-density poles

In the standard bosonic setting, the system \(S\) is coupled linearly to an environment \(E\) through
\[
H = H_S + H_E + H_I,\qquad H_I=A\otimes B,\qquad B=\sum_k g_k(a_k+a_k^\dagger),
\]
with \(A=A^\dagger\). The environment is characterized by the two-time correlation function
\[
C(\tau)=\frac{1}{2\pi}\int_{-\infty}^{\infty} d\omega\,\gamma(\omega)e^{-i\omega\tau},\qquad \tau\ge 0,
\]
where \(\gamma(\omega)\) is the full spectral density associated with \(B\) [2108.05755].

If \(\gamma(\omega)\) is meromorphic in the complex plane and decays faster than \(O(1/|\omega|)\) as \(|\omega|\to\infty\), and if its poles in the lower half-plane are
\[
z_l=\xi_l-i\lambda_l,\qquad \lambda_l>0,
\]
with residues \(r_l=\operatorname{Res}_{\omega=z_l}[\gamma(\omega)]\), then contour integration yields
\[
C(\tau)=-i\sum_{l=1}^{N} r_l e^{-iz_l \tau}.
\]
Each exponential term defines a pseudomode: a discrete bosonic oscillator \(b_l\) with frequency \(\xi_l=\operatorname{Re}z_l\), damping rate \(2\lambda_l\), and system coupling
\[
g_l'=\sqrt{-i\,r_l}.
\]
Thus, in the canonical formulation, pseudomodes are in one-to-one correspondence with poles of the analytically continued spectral density [2108.05755].

This pole-based interpretation has a direct physical reading. The factor \(e^{-iz_l\tau}=e^{-i\xi_l\tau}e^{-\lambda_l\tau}\) represents an oscillatory memory contribution with decay time \(1/\lambda_l\). Slow pseudomodes, with small \(\lambda_l\), encode long-lived memory; fast pseudomodes encode short-time correlations and approach Markovian behavior. A closely related formulation uses the frequency-domain self-energy \(\Sigma(\omega)\): if \(\Sigma(\omega)\) is rational,
\[
\Sigma_{\alpha\beta}(\omega)=\sum_\ell \frac{r_{\alpha\beta}^{(\ell)}}{\omega-z_\ell},
\]
then the corresponding memory kernel is a finite sum of exponentials and can likewise be represented by a finite pseudomode set [2605.03946].

## 2. Markovian embedding and recovery of non-Markovian dynamics

Once the pseudomodes have been introduced, the enlarged system \(S+M\) is governed by a time-independent GKSL generator. In the formulation of Ref. [2108.05755], the auxiliary Hamiltonian is
\[
H_0=H_S+\sum_l \xi_l b_l^\dagger b_l + A\otimes \sum_l g_l'(b_l+b_l^\dagger),
\]
and the density matrix obeys
\[
\frac{d}{dt}\rho(t)=-i[H_0,\rho(t)] + 2\sum_l \lambda_l\left(b_l\rho b_l^\dagger-\frac12\{b_l^\dagger b_l,\rho\}\right).
\]
Because \(r_l\) may be complex, \(g_l'\) can be complex and \(H_0\) is generally non-Hermitian; the dissipator cancels this non-Hermiticity and yields a completely positive Markovian semigroup on the extended space [2108.05755].

Pseudomode elimination, in the strict sense, is then
\[
\rho_S(t)=\operatorname{Tr}_M[\rho(t)].
\]
For Gaussian environments, the reduced dynamics is completely determined by the two-time correlation function. The core equivalence argument is therefore correlation matching: in an extended auxiliary model \(S+M+R\), each pseudomode is coupled to an independent Markovian reservoir \(R\), and the auxiliary bath operator \(B'\) is constructed so that its correlation function satisfies
\[
C'(\tau)=-i\sum_{l=1}^{N} r_l e^{-iz_l\tau}=C(\tau).
\]
Under appropriate initial states, this implies
\[
\rho'_S(t)=\rho_S(t),
\]
so tracing out the pseudomodes reproduces exactly the physical non-Markovian dynamics of the original model [2108.05755].

The same logic appears in the older Friedrichs and Jaynes–Cummings context. There, a non-local integro-differential equation for an excited-state amplitude is converted into a finite system of ordinary differential equations by introducing pseudomode amplitudes, and pseudomode elimination means solving for those amplitudes and substituting them back into the system equation, or equivalently tracing them out from the enlarged GKSL dynamics [1904.01430]. In the single-Lorentzian two-level example, the pseudomode amplitude provides a direct measure of temporarily stored excitation, and the expected time the pseudomode remains excited is
\[
\langle t_P\rangle=\frac{1}{2\lambda},
\]
which was proposed as a simple physical meaning of the reservoir memory time [1701.07597].

## 3. Rational self-energies, generalized embeddings, and compressed mode sets

A major generalization replaces the pole expansion of a spectral density by a representability criterion in terms of response functions. In the Heisenberg-picture framework of Ref. [2605.03946], one partitions the full problem into retained operators \(S_\alpha\) and an eliminated sector, derives a generalized Langevin equation with memory kernel \(K_{\alpha\beta}(t)\), and defines the self-energy
\[
\Sigma_{\alpha\beta}(\omega)=\int_0^\infty dt\, e^{i\omega t}K_{\alpha\beta}(t).
\]
If \(\Sigma(\omega)\) is rational, the eliminated sector can be replaced by a finite set of damped auxiliary modes. The paper’s central claim is that pseudomode elimination is governed by rational self-energy, not by linearity: Kerr nonlinearities, cross-Kerr terms, bilinear exchange, and three-wave mixing can all be treated provided the relevant sector reduces to a rational propagator [2605.03946].

This perspective supports a wider class of auxiliary constructions. In quasi-Lindblad pseudomode theory, the bath correlation function is represented directly as a complex weighted sum of complex exponentials,
\[
C^A(\Delta t)= (\mathbf V-i\mathbf M)e^{-\mathbf Z\Delta t}(\mathbf V+i\mathbf M)^\dagger,
\]
so the global generator is trace preserving but not necessarily completely positive [2408.15529]. The representation is not unique: different gauge choices give different pseudomode couplings while leaving the reduced system dynamics unchanged when the global dynamics is simulated exactly. This non-uniqueness makes pseudomode elimination partly a problem of representation theory: pseudomodes can be rearranged, mixed, or effectively absorbed into one another without changing the reduced dynamics, although gauge choice can strongly affect numerical stability [2408.15529].

Coupled Lindblad pseudomode theory restores CPTP dynamics while retaining compact exponential representations. Its bath correlation function has the form
\[
C^{\mathrm c}(t)=g^\dagger e^{-iKt}g,\qquad K=H-i\Gamma,
\]
and the main theorem states that if a quasi-Lindblad and a coupled-Lindblad representation generate the same bath correlation function, then they generate the same reduced system dynamics. When a suitable feasibility condition is satisfied, the coupled-Lindblad model uses the same number of pseudomodes as the quasi-Lindblad one, and the number of coupled pseudomodes only needs to scale as \(\mathrm{polylog}(T/\varepsilon)\) in simulation time \(T\) and precision \(\varepsilon\) [2506.10308]. In the fermionic setting, a complementary construction based on analytic pole placement or AAA rational approximation followed by interpolative matrix decomposition yields a compressed pseudomode count
\[
N_{\mathrm{ID}}\sim \log(T)\log(1/\varepsilon),
\]
so pseudomode elimination becomes a controlled compression of redundant exponentials while preserving impurity observables up to the target error [2409.08816].

## 4. Representative models and benchmark systems

The spin-boson model is a standard benchmark because its bath correlation function combines a structured analytic part with thermal Matsubara tails. For the underdamped Brownian oscillator spectral density
\[
J(\omega)=\frac{\alpha \omega_0^2 \Gamma\, \omega}{(\omega_0^2-\omega^2)^2+\Gamma^2\omega^2},
\]
the correlation function splits as \(C(\tau)=C_0(\tau)+M(\tau)\), where \(C_0\) is already a finite sum of exponentials and \(M\) is an infinite Matsubara series. In Ref. [2108.05755], \(C_0\) is captured exactly with only two pseudomodes because complex couplings are allowed, while \(M\) is treated either by two additional real exponentials or by keeping only the first Matsubara term plus a local dephasing term. Benchmarking against numerically exact HEOM for \(\alpha=0.25\Delta\), \(\Gamma=0.05\Delta\), \(\beta\Delta=1\), \(\epsilon=0.5\Delta\), and \(\rho_S(0)=|e\rangle\langle e|\) shows excellent agreement over both short and long times [2108.05755].

In photosynthetic light harvesting, the pseudomode approach has been used to model vibronic non-Markovianity in simplified Fenna–Matthews–Olson dimers. With a single dominant vibronic mode at \(\omega_0=202\,\mathrm{cm}^{-1}\), Huang–Rhys factor \(S=0.02\), and effective coupling \(g\approx 29\,\mathrm{cm}^{-1}\), the model predicts coherence lifetimes on the \(\sim 10\,\mathrm{ps}\) scale, compared with \(\sim 250\,\mathrm{fs}\) in naive Markovian models, and oscillation frequencies of order \(29\,\mathrm{cm}^{-1}\) [1904.01430]. In a closely related two-level setting, the pseudomode amplitude was interpreted as the temporary storage of information and excitation in the reservoir, making pseudomode elimination a concrete picture of backflow rather than only an operator identity [1701.07597].

In dissipative cavity QED, a single pseudomode can encode Fano interference between direct emission and cavity-mediated emission. Ref. [2601.10087] shows that a previously known Lindblad master equation for a two-level system plus cavity with Fano interference can be rederived from a structured-reservoir model by introducing one pseudomode. The corresponding spectral function contains a constant part and a non-Lorentzian contribution forming the Fano profile, and the constant term is essential for obtaining a Lindblad master equation [2601.10087]. This example illustrates a recurrent theme: pseudomode elimination can convert what appears as non-Markovian reservoir structure into a compact explicit mode plus simple Markovian damping.

## 5. Purified, nonlinear, and thermodynamic extensions

Purified pseudomode theory reinterprets pseudomodes so that positive- and negative-time bath correlations are carried by separate auxiliary modes acting only on one side of the density matrix. In the purified input-output model, the original non-Markovian bath is viewed as already having been obtained by mathematically eliminating pseudomodes, and the construction effectively undoes that elimination in a way that preserves access to environmental observables and non-Gaussian initial bath states [2412.04264]. The enlarged dynamics is time-local, while tracing out the purified pseudomodes recovers the same reduced system dynamics; the same formalism also gives access to field observables such as \(\langle \phi^\dagger(t)\phi(t')\rangle\), cavity occupations, and emission spectra [2412.04264].

The nonlinear extension replaces linear couplings \(sX\) by
\[
H_{\text{int}}=sQ(X),\qquad Q(X)=\sum_{n=1}^\infty \alpha_n X^n.
\]
For such interactions, reproducing the reduced dynamics requires not only matching the two-time function \(C(t)\) but also the equal-time moment
\[
\sigma^2=C(0).
\]
Ref. [2502.18707] therefore augments the pseudomode construction with an extra zero-frequency pseudomode and then performs a Liouvillian Schrieffer–Wolff elimination of fast modes to derive a purified pseudomode master equation. The method is validated for spontaneous decay of a two-level atom in a lossy cavity and for the resonance fluorescence spectrum of a quantum dot with a phonon environment [2502.18707].

In strong-coupling quantum thermodynamics, pseudomodes are retained rather than immediately eliminated, because they provide a finite Markovian representation from which bath-dependent observables can be recovered. For bosonic baths linearly coupled to the system, exact expressions for heat, work, and average system–bath interaction energy can be written in terms of bath autocorrelation functions and two-time system correlators; after the pseudomode replacement, these quantities can be evaluated efficiently in terms of one-time expectation values of the system and the pseudomodes [2407.17886]. For example, the heat associated with bath \(j\) becomes
\[
Q_j(t_f)=-\int_0^{t_f} dt\, \operatorname{Tr}\!\left\{\widehat H_{IM,j}(t)\,\mathcal L_{M,j}[\rho_{SM}(t)]\right\},
\]
so the physical continuum bath has been eliminated in favor of a finite pseudomode network and Lindblad generator [2407.17886].

## 6. Alternative meanings, numerical reduction, and limitations

Outside the canonical open-system setting, the phrase “pseudomode elimination” acquires additional technical meanings. In the many-body recursion method, a dissipatively deformed Krylov-space Liouvillian yields an autocorrelation expansion
\[
C(t)=\sum_{l=0}^\infty (V_l^0)^2 e^{\lambda_l t},
\]
and elimination means truncating this expansion to the first few dominant pseudomodes. In the 2D XX and 2D Ising examples, only a few low-lying pseudomodes are needed to reproduce the correlation function accurately over the accessible time range [2407.12495].

A different reduction appears in the stochastic pseudomode model, which decomposes a Gaussian bosonic environment into a finite set of zero-temperature ancillary bosonic modes plus classical stochastic fields. For rational spectral densities, all parameters can be specified analytically, and the number of ancillary quantum degrees of freedom is reduced because Matsubara and other symmetric contributions are moved into the classical field sector [2301.07554]. The classical fields can even be imaginary-valued, which the authors show can decrease the entropy of the system, in contrast to real-valued fields [2301.07554]. Ref. [2311.15240] pushes this further by replacing explicitly unphysical pseudomodes with measurement results over ensembles of physical systems and an extrapolation procedure, thereby eliminating unphysical auxiliary modes as direct dynamical degrees of freedom while preserving the target reduced dynamics [2311.15240].

The main limitations are structural. Exactness typically requires that the bath correlation function be representable as a finite sum of exponentials, equivalently that the relevant spectral density or self-energy be rational or meromorphic with finitely many relevant poles [2108.05755][2605.03946]. Thermal Matsubara tails, low-frequency singular structure, or long-time algebraic behavior generally require truncation, fitting, or hybrid stochastic treatments [2108.05755][2301.07554]. In quasi-Lindblad formulations, the pseudomode representation is more compact but the global dynamics may lose positivity, and gauge choice becomes a stability problem rather than a merely formal one [2408.15529]. In coupled-Lindblad and fermionic compression schemes, the asymptotic mode count can be polylogarithmic in time and accuracy, but the construction still depends on the quality of rational approximation and on feasibility conditions or interpolation tolerances [2506.10308][2409.08816].

Taken together, these developments establish pseudomode elimination as both a structural theorem and a practical methodology. In its narrow sense it is the tracing out of auxiliary damped modes from an enlarged GKSL evolution to recover non-Markovian reduced dynamics. In its broader modern sense it also denotes the systematic reduction, compression, purification, or physical re-realization of auxiliary-mode descriptions so that environmental memory is represented with as few effective degrees of freedom as possible, while retaining control over reduced dynamics, observables, and numerical stability [2108.05755][2506.10308].

Source: https://www.emergentmind.com/topics/pseudomode-elimination