Quantum Fokker–Planck Master Equation

Updated 6 July 2026

Quantum Fokker–Planck master equation is a family of quantum evolution equations combining drift, diffusion, and dissipation across operator, phase-space, and noncommutative frameworks.
It unifies canonical quantization, Wigner-transform methods, and noncommutative techniques to ensure trace preservation, detailed balance, and controlled entropy production.
Applications span continuous measurement, quantum Brownian motion, and even cosmological models, linking microscopic dynamics with coarse-grained diffusion phenomena.

In the literature surveyed here, the quantum Fokker–Planck master equation denotes a family of irreversible quantum evolution equations that combine drift, diffusion, and dissipation in operator space, phase space, or enlarged system–detector spaces. Depending on the framework, the dynamical object may be a density operator, a Wigner function, a reduced diagonal density functional, or an element of a noncommutative $L^p$ space; correspondingly, the generator may appear as a GKLS dissipator, a Wigner–Poisson/Fokker–Planck PDE, or an operator-valued Fokker–Planck equation for classical filter variables coupled to quantum backaction (Oliveira, 2023, Alejo et al., 2018, Labuschagne et al., 2021, Annby-Andersson et al., 2021).

1. Scope and competing definitions

A central feature of the subject is that there is no single universally fixed “quantum Fokker–Planck master equation.” Rather, the term covers several mathematically distinct but structurally related constructions. Some are obtained by canonical quantization of classical Fokker–Planck generators; some arise from Wigner transforms of Lindblad or quantum Brownian-motion equations; some are exact or approximate reduced equations for continuously monitored systems; and some are formulated directly on noncommutative $L^p$ spaces over von Neumann algebras (Oliveira, 2023, Escalante, 2023, Labuschagne et al., 2021).

Formulation	State variable	Representative structure
Operator-space FP/GKLS	$\rho$	current–divergence form, transformable to Lindblad form
Wigner/phase-space FP	$W(x,p)$ or $w(x,\xi)$	drift, friction, diffusion, pseudo-differential potential
Noncommutative FP	$f_t \in L^2(M)$	$H_q=C+V$ built from derivations and positive maps
Measurement-feedback QFPME	$\hat\rho(D)$ or $\hat\rho(G)$	quantum term plus FP drift–diffusion in detector/filter variables

In the operator-space approach, the quantum Fokker–Planck structure is a continuity equation in operator space. In phase-space approaches, the equation acts on a Wigner function and resembles a Kramers or Ornstein–Uhlenbeck equation supplemented by quantum potential terms. In measurement-and-feedback theory, the Fokker–Planck variable is neither position nor momentum but a classical detector output or filtered signal, and the equation propagates a joint system–signal state (Oliveira, 2023, Alejo et al., 2018, Annby-Andersson et al., 2021).

A further distinction concerns what is meant by “quantum.” In some works the evolution remains fully quantum at the level of $\rho$ or $L^p$ 0. In others, a classical Fokker–Planck equation is derived from an underlying Schrödinger dynamics after coarse-graining away coherences in a preferred basis. That latter use retains a quantum origin but not a fully quantum reduced state description (Lone, 28 May 2025).

2. Operator formulations, GKLS structure, and noncommutative semigroups

One operator-theoretic construction begins from a canonical quantization of the classical Fokker–Planck equation and writes the dissipative dynamics in a current–divergence form,

$L^p$ 1

with

$L^p$ 2

This formulation makes the quantum analogues of probability current and detailed balance explicit, and the same paper shows that the resulting dissipator can be transformed into Lindblad form with a positive Kossakowski matrix, thereby ensuring trace preservation and positivity of the density operator (Oliveira, 2023).

In that same framework, thermal contact is encoded by choosing

$L^p$ 3

so that detailed balance is expressed as

$L^p$ 4

The Gibbs state is then a fixed point, and the Fokker–Planck structure yields explicit expressions for entropy production and entropy flux. The same work emphasizes that the entropy production rate need not vanish for a closed system except in equilibrium (Oliveira, 2023).

A different, more abstract version is developed on $L^p$ 5-finite von Neumann algebras equipped with Haagerup $L^p$ 6 spaces. There the quantum Fokker–Planck equation is

$L^p$ 7

where the “quantized Laplacian”

$L^p$ 8

is built from weak $L^p$ 9-closed derivations, while the potential term is $\rho$ 0 for a unital Schwarz map $\rho$ 1. Under Detailed Balance II, $\rho$ 2 generates conservative $\rho$ 3-Markov $\rho$ 4-semigroups, induced by integrable Markov semigroups on $\rho$ 5; if $\rho$ 6 commutes with $\rho$ 7, the semigroups satisfy DBII as well (Labuschagne et al., 2021).

That noncommutative theory also gives an asymptotic analysis by spectral methods. If $\rho$ 8 has a spectral gap $\rho$ 9, then

$W(x,p)$ 0

and even without a gap one has strong convergence to $W(x,p)$ 1 as $W(x,p)$ 2. Convergence in $W(x,p)$ 3 is supplemented by a noncommutative Csiszár–Kullback inequality,

$W(x,p)$ 4

which controls approach to equilibrium in terms of relative entropy (Labuschagne et al., 2021).

By contrast, the high-temperature Markovian Caldeira–Leggett equation,

$W(x,p)$ 5

is explicitly identified as not being of Lindblad form; the low-temperature and positivity issues that accompany this structure motivate several later reformulations (Bolivar, 2010).

3. Wigner-space and phase-space equations

In kinetic and open-systems theory, the quantum Fokker–Planck master equation often appears in Wigner representation. For a reduced density matrix $W(x,p)$ 6, the Wigner function is

$W(x,p)$ 7

and the nonlinearity enters through the Wigner pseudo-differential operator $W(x,p)$ 8. In the non-Markovian Unruh–Zurek and Hu–Paz–Zhang models with self-consistent Poisson coupling, the equations take the form

$W(x,p)$ 9

and

$w(x,\xi)$ 0

with time-dependent integral coefficients encoding memory of the bath and a nonlinear Hartree potential satisfying $w(x,\xi)$ 1, $w(x,\xi)$ 2. The analysis constructs explicit Gaussian propagators, proves local existence by Banach fixed point, and obtains global mild solutions by kinetic energy estimates and Lieb–Thirring control of the density (Alejo et al., 2018).

For an Ohmic quantum Langevin dynamics, a simpler reduced formulation arises by integrating out momentum and studying the momentum-independent Wigner function $w(x,\xi)$ 3. The resulting exact reduced equation is

$w(x,\xi)$ 4

with $w(x,\xi)$ 5 and $w(x,\xi)$ 6. The adjoint density-matrix equation becomes

$w(x,\xi)$ 7

which is considerably simpler than the full phase-space Schramm–Jung–Grabert equation while preserving the exact $w(x,\xi)$ 8-statistics (Colmenares, 2017).

A related non-Markovian Caldeira–Leggett construction starts instead from a classical non-Markovian Fokker–Planck equation and quantizes it by a “dynamical quantization” procedure. In position representation the result is

$w(x,\xi)$ 9

with $f_t \in L^2(M)$ 0, $f_t \in L^2(M)$ 1, and a concrete non-Markovian choice $f_t \in L^2(M)$ 2. In Wigner representation for a free particle, this becomes a quantum Rayleigh equation with time-dependent momentum diffusion (Bolivar, 2010).

Another phase-space version follows from a Markovian quantum Langevin equation for a trapped particle in a quantum-thermal Ohmic bath. In the underdamped harmonic regime one obtains the semiclassical Kramers-type equation

$f_t \in L^2(M)$ 3

with

$f_t \in L^2(M)$ 4

For generic nonlinear $f_t \in L^2(M)$ 5, the quantum Fokker–Planck equation acquires the higher-derivative correction

$f_t \in L^2(M)$ 6

where $f_t \in L^2(M)$ 7 is given by an infinite series in $f_t \in L^2(M)$ 8. For harmonic confinement, the series truncates exactly; for nonlinear potentials, truncations must be treated cautiously because of Pawula’s theorem (Furutani et al., 2023).

4. Memory, regularity, and Markovian limits

Non-Markovianity enters these equations in different ways. In the UZ and HPZ Wigner–Poisson models, it appears through integral coefficients $f_t \in L^2(M)$ 9, $H_q=C+V$ 0, $H_q=C+V$ 1, and $H_q=C+V$ 2 built from bath spectral densities, thermal factors $H_q=C+V$ 3, and time integrals over the full past of the system. This modifies both drift and dissipation and produces propagators with singular short-time scaling. In particular, the paper derives

$H_q=C+V$ 4

and

$H_q=C+V$ 5

with an associated significant lack of Sobolev regularity relative to Markovian counterparts. The global theory therefore relies on mild formulations, time-dependent Gaussian bounds, continuity equations, and kinetic-energy control rather than on uniform parabolic smoothing (Alejo et al., 2018).

A gravitational variant yields a non-Markovian quantum master equation and associated Wigner Fokker–Planck equation for two quantum masses interacting with a graviton bath at all temperatures. The graviton spectral density is

$H_q=C+V$ 6

so it is Ohmic at the field level, but the system coupling is quadratic in the separation coordinates, which generates quartic derivative structures and energy-basis decoherence terms. In the high-temperature short-memory regime the master equation reduces to the ABH structure $H_q=C+V$ 7, whereas in the low-temperature non-Markovian regime the off-diagonal momentum-space elements decay as

$H_q=C+V$ 8

The zero-temperature part is logarithmic in time, and the temperature-dependent part is quadratic in time, which differs sharply from the exponential Markovian behavior (Cho et al., 16 Apr 2025).

A different limit appears in derivations from reversible Schrödinger dynamics. Expanding the state in the eigenbasis of a time-independent $H_q=C+V$ 9, one may write the populations $\hat\rho(D)$ 0 as

$\hat\rho(D)$ 1

and then neglect the coherence correction $\hat\rho(D)$ 2 under coarse-graining. This yields a classical master equation

$\hat\rho(D)$ 3

and, in the continuum limit, a Fokker–Planck equation

$\hat\rho(D)$ 4

A plausible implication is that some occurrences of “quantum Fokker–Planck” in the literature refer not to a fully quantum reduced state equation, but to a classical diffusion equation derived from an underlying quantum dynamics after basis-dependent decoherence and coarse-graining (Lone, 28 May 2025).

5. Continuous measurement, feedback, and signal filtering

In continuous-measurement theory, the quantum Fokker–Planck master equation describes the joint dynamics of a quantum system and a classical detector output. For continuous monitoring of a Hermitian observable $\hat\rho(D)$ 5 with measurement strength $\hat\rho(D)$ 6, and for a detector with finite bandwidth $\hat\rho(D)$ 7 modeled by an Ornstein–Uhlenbeck filter, the operator-valued density $\hat\rho(D)$ 8 obeys

$\hat\rho(D)$ 9

where

$\hat\rho(G)$ 0

Here $\hat\rho(G)$ 1 may depend nonlinearly on the filtered outcome $\hat\rho(G)$ 2, so the formalism covers both linear and nonlinear feedback. In the fast-measurement limit $\hat\rho(G)$ 3, detector elimination yields a reduced Markovian master equation

$\hat\rho(G)$ 4

which generalizes the Wiseman–Milburn equation and supports analytic treatments of threshold engines and driven qubits (Annby-Andersson et al., 2021).

That construction has been extended to arbitrary linear filtering. If the filter state $\hat\rho(G)$ 5 satisfies

$\hat\rho(G)$ 6

with $\hat\rho(G)$ 7, then the joint operator-valued density $\hat\rho(G)$ 8 obeys

$\hat\rho(G)$ 9

This encompasses single-pole low-pass filters, cascaded low-pass filters, band-pass filters, and finite-order ODE kernels. In harmonic-oscillator cooling with one filtered layer, the ensemble-averaged steady-state energy is

$\rho$ 0

with optimum $\rho$ 1, for which $\rho$ 2. Two- and three-layer cascades and band-pass filters admit explicit analytical formulas and reveal bandwidth-dependent stability boundaries (Sousa, 20 Sep 2025).

These feedback QFPMEs differ conceptually from phase-space quantum Brownian-motion equations. The Fokker–Planck variable is now the filtered signal rather than a particle coordinate or momentum, the diffusion tensor is $\rho$ 3, and the Lindblad term $\rho$ 4 represents measurement backaction. The structure is therefore a hybrid quantum–classical transport equation rather than a Wigner kinetic equation (Annby-Andersson et al., 2021, Sousa, 20 Sep 2025).

6. Numerical reformulations and domain-specific extensions

For numerical work on Markovian open systems, the Wigner–Fokker–Planck equation can be transformed by Fourier transform in momentum to symmetrized position coordinates $\rho$ 5. Starting from

$\rho$ 6

the transform removes the nonlocal pseudo-differential potential term and yields the local convection–diffusion equation

$\rho$ 7

In one dimension, decomposing $\rho$ 8 gives a coupled convection–diffusion system suitable for NIPG-DG discretization. The paper reports that for piecewise linear DG functions the method exhibits $\rho$ 9 convergence in $L^p$ 00 for the harmonic steady state and $L^p$ 01 projection error, and it identifies a cost crossover in favor of the transformed equation for meshes with $L^p$ 02 intervals in momentum or $L^p$ 03 (Escalante, 2023).

Quantum Fokker–Planck methods also appear in cosmology. For a massless scalar field with quartic interaction in de Sitter space, coarse-graining at

$L^p$ 04

and integrating out short modes yields a functional equation for the reduced diagonal density $L^p$ 05,

$L^p$ 06

For the coarse-grained homogeneous field this reduces to

$L^p$ 07

with stationary distribution $L^p$ 08. Here again the phrase “quantum Fokker–Planck” refers to a reduced equation derived directly from the quantum theory, but acting on a coarse-grained probability density rather than on a full density matrix (Collins et al., 2017).

Applications to condensed-matter and atomic systems follow the same pattern of drift–diffusion plus quantum correction. In Josephson circuits and Bose Josephson junctions, the trapped-particle Ohmic-bath formalism produces stationary semiclassical distributions of Boltzmann form with $L^p$ 09 replacing $L^p$ 10, and the corresponding quantum-corrected Fokker–Planck equation carries higher-order velocity derivatives for nonlinear potentials. In superconducting circuits the cutoff may be chosen as $L^p$ 11, while in one-dimensional atomic gases the cutoff is $L^p$ 12; both cases are treated as experimentally accessible applications of the same quantum-thermal Fokker–Planck framework (Furutani et al., 2023).

Across these domains, the common content is not a single canonical equation but a recurrent structural motif: deterministic Hamiltonian or Liouvillian transport, irreversible friction or drift, diffusion fixed by environmental fluctuations or measurement imprecision, and a quantum correction term whose form depends on the chosen representation. The differences among formulations—GKLS versus Wigner, Markovian versus non-Markovian, microscopic versus coarse-grained, operator-valued versus classical-valued—are therefore constitutive rather than peripheral to the subject (Oliveira, 2023, Alejo et al., 2018, Annby-Andersson et al., 2021).