Boltzmann to Lindblad: Classical and Quantum Approaches to Out-of-Equilibrium Statistical Mechanics

Abstract: Open quantum systems play a central role in contemporary nanoscale technologies, including molecular electronics, quantum heat engines, quantum computation and information processing. A major theoretical challenge is to construct dynamical models that are simultaneously consistent with classical thermodynamics and complete positivity. In this work, we develop a framework that addresses this issue by extending classical stochastic dynamics to the quantum domain. We begin by formulating a generalized Langevin equation in which both friction and noise act symmetrically on the two Hamiltonian equations. From this, we derive a generalized Klein-Kramers equation expressed in terms of Poisson brackets, and we show that it admits the Boltzmann distribution as its stationary solution while satisfying the first and second laws of thermodynamics along individual trajectories. Applying canonical quantization to this classical framework yields two distinct quantum master equations, depending on whether the friction operators are taken to be Hermitian or non-Hermitian. By analyzing the dynamics of a harmonic oscillator, we determine the conditions under which these equations reduce to a Lindblad-type generator. Our results demonstrate that complete positivity is ensured only when friction and noise are included in both Hamiltonian equations, thus fully justifying the classical construction. Moreover, we find that the friction coefficients must satisfy the same positivity condition in both the Hermitian and non-Hermitian formulations, revealing a form of universality that transcends the specific operator representation. The formalism offers a versatile tool for deriving quantum versions of the thermodynamic laws and is directly applicable to a wide class of nonequilibrium nanoscale systems.

• The paper introduces a unified framework linking classical stochastic thermodynamics with quantum open system dynamics, ensuring Lindblad complete positivity.
• It generalizes Langevin and Fokker–Planck dynamics via symmetric friction and noise to meet thermodynamic laws and quantum coherence requirements.
• Numerical results confirm that only models enforcing symmetric friction and noise yield consistent quantum evolution with physically valid density matrices.

Classical and Quantum Approaches to Nonequilibrium Statistical Mechanics: From Boltzmann to Lindblad

Motivation and Scope

The paper "Boltzmann to Lindblad: Classical and Quantum Approaches to Out-of-Equilibrium Statistical Mechanics" (2512.11613) develops a systematic framework linking classical stochastic thermodynamics and quantum open system dynamics. The core objective is to construct dynamical models that are consistent with both classical thermodynamics and the requirement of complete positivity in quantum evolution, addressing key deficiencies in several well-used phenomenological quantum master equations.

Significant extensions are introduced to classical Langevin and Fokker–Planck–Klein–Kramers dynamics, leading to quantum master equations formulated via canonical quantization. The framework unifies treatment of friction and noise, justifies their symmetric inclusion in both equations of motion, identifies universal constraints for quantum complete positivity, and draws rigorous connections to the Lindblad form and thermodynamic structure.

Generalized Classical Stochastic Dynamics

The authors introduce a generalization of classical stochastic Hamiltonian dynamics, wherein both friction and noise act symmetrically on position and momentum equations. The resulting system of SDEs, for a Hamiltonian $\mathcal{H}_0 = K_0 + V_0$, possesses the structure:

• Momentum: $$\dot{\vec{p}}_i = -\frac{\partial \mathcal{H}_0}{\partial \vec{r}_i} - m_i\beta_p\frac{\partial \mathcal{H}_0}{\partial \vec{p}_i} + \sqrt{D_p m_i}\vec{n}_{p,i}(t) + \vec{f}_i(t)$$
• Position: $$\dot{\vec{r}}_i = \frac{\partial \mathcal{H}_0}{\partial \vec{p}_i} - m_i\beta_q\frac{\partial \mathcal{H}_0}{\partial \vec{r}_i} + \sqrt{D_q m_i}\vec{n}_{q,i}(t)$$

This symmetric construction generalizes the standard Langevin model, recovers it as a special case, and leads to a Fokker–Planck (Klein–Kramers) equation in phase space for the probability density $W(\vec{r}, \vec{p}, t)$, naturally written in Poisson bracket form for direct quantization.

Classically, the stationary state is guaranteed to be the Gibbs–Boltzmann distribution. The formalism produces explicit expressions for internal energy, heat, and work, and retrieves correct forms of the first and second laws of thermodynamics at the trajectory level. Importantly, the friction and diffusion coefficients must satisfy the fluctuation-dissipation relation $D_p = k_B T\beta_p$, $D_q = k_B T \beta_q$ for thermodynamic consistency.

Figure 1: The parametric constraints on friction coefficients ($\beta_p$, $\beta_q$) for Lindblad complete positivity delineate a region in the plane, shown for various temperatures via the parameter $\xi$.

Quantization and Construction of Quantum Master Equations

Canonical Quantization Procedure

Canonical quantization is performed by replacing Poisson brackets with commutators $[\cdot, \cdot]/(i\hbar)$, yielding stochastic dynamics for the quantum density matrix $\varrho$. The transition of classical dissipative and noisy operators to quantum analogs is nontrivial and admits two routes, depending on the Hermiticity imposed on friction operators.

Hermitian vs Non-Hermitian Friction Operators

• Hermitian Operators: By symmetrizing products (e.g., mapping $p W$ to $\frac{1}{2}(p \varrho + \varrho p)$), one constructs friction operators $\Theta^{p,q}$ which are themselves Hermitian. The explicit forms of these operators are derived such that the stationary solution $\varrho_{eq} \propto e^{-\mathcal{H}_0 / (k_B T)}$ is preserved.
• Non-Hermitian Operators: Alternatively, the friction operators may remain non-Hermitian ($\Xi^{p,q}$), simplifying certain analytical properties but decoupling direct association with quantum observables. The resulting equations are algebraically distinct, but for quadratic Hamiltonians, both approaches are shown to yield equivalent physical constraints.

Structure of Quantum Master Equations

Both quantization routes generate quantum master equations of the form:

$$\frac{d\varrho}{dt} = \frac{1}{i\hbar}[\mathcal{H}_0, \varrho] + \text{dissipation} + \text{noise}$$

The dissipation and noise terms involve double commutators and symmetrized commutators with the friction/noise operators, whose explicit forms depend on temperature and the system's energy spectrum.

Complete Positivity and Lindblad Structure

A central technical advance is establishing precise conditions under which the constructed quantum master equations assume Lindblad form, guaranteeing trace preservation and positivity for arbitrary initial states. This addresses well-known failures of the Redfield and phenomenological equations, which can give rise to negative density eigenvalues.

For the harmonic oscillator:

• The Lindblad structure requires (see (2512.11613), Eq. 98) that both $\beta_p > 0$ and $\beta_q > 0$, compelling inclusion of noise and friction in both Hamiltonian equations.
• The friction coefficients must satisfy a universal constraint, analytically derived, and independent of operator Hermiticity. This constraint is encapsulated by:

$$\left(\frac{ \cosh \xi - 1 }{ \sinh \xi } \right)^2 \leq \frac{ \beta_p }{ m^2 \omega^2 \beta_q } \leq \left(\frac{ \cosh \xi + 1 }{ \sinh \xi } \right)^2$$

where $\xi = \hbar \omega / (2k_B T)$.

• On this admissible domain, the quantum master equation matches the general Gorini–Kossakowski–Sudarshan–Lindblad (GKLS) generator, and reduces to the Quantum Optical Master Equation under specific parameter choices.

  Figure 2: Time-evolution of energy, entropy production, and density matrix positivity for the quantum harmonic oscillator under five models; only those respecting the symmetric friction/noise prescription ensure full Lindblad positivity.

Quantum Stochastic Thermodynamics

The quantum framework provides systematic, operator-level expressions for energy, work, heat, and both von Neumann and relative entropy production. The first law is retrieved as:

$$\frac{d\mathcal{E}}{dt} = \text{Power from external forces} + \text{Heat flux}$$

Heat flux, work, and entropy production are shown to have unique quantum-corrected forms tightly linked to the structure of the friction operators. The second law arises naturally as the monotonic decrease of quantum relative entropy under completely positive, trace-preserving evolution, i.e.,

$\frac{ d\mathcal{S}_p }{dt } \geq 0$

where entropy production vanishes only in the stationary (equilibrium) state.

Numerical Results: Harmonic Oscillator Example

The authors numerically solve the master equations for the quantum harmonic oscillator, starting from a variety of mixed and pure initial states, and using both Hermitian and non-Hermitian friction operator schemes along with established phenomenological models (including Caldeira–Leggett).

Key findings include:

• For models satisfying the friction/noise symmetry and the complete positivity constraint, energy converges correctly to the quantum Gibbs value, and the density matrix remains strictly positive for all times and initial states.
• Models lacking these properties (e.g., only friction/noise in momentum, standard Caldeira–Leggett) lead to negative density eigenvalues under certain pure initial conditions, violating physicality and indicating breakdown of complete positivity.

  Figure 3: Entropy production and negative eigenvalues of $\varrho(t)$ following a pure state initial condition. Only symmetric, constraint-satisfying models prevent the emergence of negative, unphysical states.

Theoretical and Practical Implications

• Universality of Symmetric Friction/Noise: The necessity to include friction and noise in both canonical variables emerges as a universal criterion for the construction of physically consistent quantum master equations, with direct classical-to-quantum correspondence.
• Thermodynamic Consistency: The approach justifies the definition of quantum heat, work, generalized equipartition, and entropy production, and their time evolution matches macroscopic thermodynamic laws.
• Operator Representation Independence: The universal constraint for complete positivity does not depend on the Hermitian structure of the friction operators, suggesting deep robustness.
• Immediate Applicability: The formalism provides a robust framework for consistently modeling nonequilibrium thermodynamics in open quantum systems—quantum heat engines, molecular electronics, dissipative quantum computing elements, quantum optomechanics, and strong-coupling nanodevices.

Conclusion

This study systematically derives and justifies a unified classical-quantum approach to out-of-equilibrium statistical mechanics, addressing longstanding issues concerning quantum dissipative dynamics and thermodynamic consistency. The necessity of symmetric inclusion of friction and noise is established both theoretically and numerically as the minimal requirement for generating Lindblad-type master equations that ensure completely positive quantum evolution. The generalized formalism is anticipated to offer a practical foundation for modeling a wide range of nonequilibrium quantum thermodynamic phenomena.

Future work includes extension of the method to arbitrary potentials and further numerical and analytic study of model-specific forms of friction operators to secure universality of complete positivity constraints in general quantum settings.