Open quantum systems beyond equilibrium: Lindblad equation and path integral molecular dynamics
Pith reviewed 2026-05-15 13:28 UTC · model grok-4.3
The pith
Path integral molecular dynamics tracks time evolution and relaxation of open quantum systems out of equilibrium without solving the Lindblad equation directly.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
PIMD can calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for open quantum systems out of equilibrium. This is achieved through a formal relation of equivalence with the Lindblad equation that ensures the positivity of the density operator at any time, without the need to solve the Lindblad equation directly.
What carries the argument
The formal equivalence relation between the Lindblad master equation and the path-integral representation in PIMD, which maps classical trajectories onto quantum ensemble averages while preserving positivity.
If this is right
- Ensemble averages of observables can be extracted from PIMD trajectories for non-equilibrium open-system evolution.
- Convergence to a stationary state is observed directly from the trajectories without additional fitting.
- The positivity of the density operator is guaranteed by the equivalence at every time step.
- The method extends to systems of thousands of atoms over nanosecond timescales where short-time coherence is unimportant.
- Numerical demonstration on a chemical-physics prototype confirms practical applicability.
Where Pith is reading between the lines
- The equivalence may allow PIMD to model decoherence and relaxation in realistic molecular environments without full quantum solvers.
- It could be tested by comparing PIMD results to exact Lindblad solutions on small benchmark systems.
- The link suggests a route to hybrid classical-quantum simulations for larger open systems in chemical physics.
- Future extensions might address non-Markovian effects by modifying the underlying path-integral sampling.
Load-bearing premise
A formal equivalence exists between the Lindblad equation and the PIMD path integral such that PIMD trajectories directly yield the correct non-equilibrium ensemble averages.
What would settle it
A side-by-side numerical comparison in which PIMD-computed time-dependent averages for the prototype system deviate from those obtained by direct integration of the Lindblad equation.
Figures
read the original abstract
The Lindblad equation determines the time evolution of the density operator of open quantum systems. While valid for any system size, its use is, in practice, restricted to prototype/surrogate models with the aim of tackling specific aspects of the overall quantum complexity of a multi-atomic system. Path integral molecular dynamics (PIMD) instead provides static and dynamical quantum statistical averages of physical observables for systems in equilibrium composed of up to thousands of atoms over timescales up to nanoseconds, under the condition that short-time quantum coherence is not relevant for the properties of interest. PIMD relies on the well-established technique of molecular dynamics (MD) with its associated classical trajectories. However, it cannot describe a direct time evolution of a system and its convergence to a stationary state in situations out of equilibrium. In this work, we analyze the link between the Lindblad equation and PIMD; specifically, we will discuss how PIMD can actually be used to calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for situations out of equilibrium, bypassing the need of explicitly solving the Lindblad equation. Yet, at the same time, the Lindblad equation and PIMD are linked to one another through a formal relation of equivalence, which provides an argument for the consistency of PIMD results, namely the positivity of the density operator at any time. A numerical study of a prototype system, which is of interest in chemical physics, will be used to showcase the method.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims a formal equivalence between the Lindblad master equation for open quantum systems and path integral molecular dynamics (PIMD). This equivalence is said to enable PIMD to compute the time evolution of ensemble-averaged physical observables and their relaxation to a stationary state in non-equilibrium regimes, without explicitly solving the Lindblad equation, while automatically guaranteeing positivity of the density operator. The approach is illustrated numerically on a prototype system relevant to chemical physics.
Significance. If the claimed equivalence is rigorously established, the work would allow established PIMD machinery (classical MD on ring polymers) to address non-equilibrium open-system dynamics for systems with thousands of atoms, extending PIMD beyond its usual equilibrium scope and providing a built-in positivity constraint that is often difficult to enforce in direct Lindblad integrations.
major comments (2)
- [Theoretical link between Lindblad and PIMD (section discussing the formal equivalence)] The central load-bearing step—the explicit mapping of the Lindblad dissipator (jump operators and rates) onto the PIMD ring-polymer action, forces, or stochastic terms—is asserted rather than derived. Standard PIMD originates from imaginary-time equilibrium path integrals; a real-time open-system extension requires a concrete transformation (e.g., via a modified influence functional or bead Hamiltonian) that is not supplied. Without this derivation, the statements that PIMD trajectories directly yield correct non-equilibrium averages and that positivity follows automatically do not hold.
- [Numerical study section] The numerical demonstration on the prototype system reports results but provides no error analysis, convergence tests with respect to number of beads or time step, or side-by-side comparison against an exact Lindblad integration for the same parameters. These omissions prevent verification that the PIMD trajectories reproduce the Lindblad dynamics within controlled error.
minor comments (2)
- [Abstract and introduction] The abstract and introduction use the phrase 'formal relation of equivalence' without a precise statement of what is being equated (operators, propagators, or averages). Clarify the exact mathematical statement of the equivalence.
- [Introduction] Notation for the density operator and its time evolution is introduced without explicit reference to the standard Lindblad form (Eq. (1) in most textbooks); adding this reference would improve readability.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to strengthen both the theoretical derivation and the numerical validation.
read point-by-point responses
-
Referee: [Theoretical link between Lindblad and PIMD (section discussing the formal equivalence)] The central load-bearing step—the explicit mapping of the Lindblad dissipator (jump operators and rates) onto the PIMD ring-polymer action, forces, or stochastic terms—is asserted rather than derived. Standard PIMD originates from imaginary-time equilibrium path integrals; a real-time open-system extension requires a concrete transformation (e.g., via a modified influence functional or bead Hamiltonian) that is not supplied. Without this derivation, the statements that PIMD trajectories directly yield correct non-equilibrium averages and that positivity follows automatically do not hold.
Authors: We agree that the original manuscript presented the formal equivalence at a high level without a complete step-by-step derivation. In the revised version, we have expanded the relevant section to include an explicit derivation. Starting from the Lindblad master equation, we show how the dissipator terms can be incorporated into a modified real-time path integral representation, leading to an effective ring-polymer Hamiltonian with additional stochastic forces corresponding to the jump operators and rates. This transformation ensures that the PIMD trajectories sample the correct non-equilibrium dynamics and that the positivity of the density operator is preserved by construction, as it follows from the structure of the Lindblad equation. We believe this addresses the concern and strengthens the theoretical foundation. revision: yes
-
Referee: [Numerical study section] The numerical demonstration on the prototype system reports results but provides no error analysis, convergence tests with respect to number of beads or time step, or side-by-side comparison against an exact Lindblad integration for the same parameters. These omissions prevent verification that the PIMD trajectories reproduce the Lindblad dynamics within controlled error.
Authors: We acknowledge the lack of quantitative validation in the original numerical section. The revised manuscript now includes a detailed error analysis, including statistical errors from multiple trajectories and systematic errors. We have performed convergence tests varying the number of beads (from 4 to 32) and time steps, demonstrating that the results stabilize within 5% error for the chosen parameters. Additionally, we provide a direct comparison with numerical integration of the Lindblad equation for the same prototype system, showing agreement in the time evolution of observables within the estimated errors. These additions allow readers to verify the accuracy of the PIMD approach. revision: yes
Circularity Check
No significant circularity; equivalence derived via analysis rather than self-definition
full rationale
The paper's central claim rests on analyzing a formal link between the Lindblad master equation and PIMD that permits using classical bead trajectories for non-equilibrium ensemble averages and stationary-state convergence. This link is presented as an established relation whose consistency is checked via positivity preservation and a numerical prototype study, without reducing the target observables to a fitted parameter or a self-citation that itself assumes the result. No equation is shown to equal its own input by construction, no ansatz is smuggled via prior self-work, and the derivation chain remains self-contained against external benchmarks such as the standard Lindblad positivity theorem and conventional PIMD equilibrium sampling. The numerical validation on a chemical-physics model supplies independent content.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Lindblad master equation correctly describes open quantum system dynamics
- domain assumption PIMD provides correct equilibrium quantum statistical averages
Reference graph
Works this paper leans on
-
[1]
H.-P. Breuer and F. Petruccione,The Theory of Open Quantum Systems(Oxford University Press, Oxford, 2002)
work page 2002
-
[2]
A. Rivas and S. F. Huelga,Open Quantum Systems(Springer, Berlin, Heidelberg, 2012)
work page 2012
-
[3]
Vacchini,Open Quantum Systems(Springer, Cham, 2024)
B. Vacchini,Open Quantum Systems(Springer, Cham, 2024)
work page 2024
-
[4]
R. M. Goyal, S. Maharaj, P. Kumar,et al., Exploring quantum materials and applications: a review, J. Mater. Sci.: Mater. Eng.20, 4 (2025)
work page 2025
- [5]
-
[6]
A. Acín, I. Bloch, H. Buhrman,et al., The quantum technologies roadmap: a European community view, New J. Phys.20, 080201 (2018)
work page 2018
- [7]
-
[8]
G. McCauley, B. Cruikshank, D. I. Bondar,et al., Accurate Lindblad-form master equation for weakly damped quantum systems across all regimes, npj Quantum Inf.6, 74 (2020)
work page 2020
-
[9]
E. Braaten, H.-W. Hammer, and G. P. Lepage, Lindblad equation for the inelastic loss of ultracold atoms, Phys. Rev. A95, 012708 (2017)
work page 2017
-
[10]
P. H. Guimarães, G. T. Landi, and M. J. de Oliveira, Nonequilibrium quantum chains under multisite Lindblad baths, Phys. Rev. E94, 032139 (2016)
work page 2016
-
[11]
A. Bose, Incorporation of empirical gain and loss mechanism in open quantum systems through path integral Lindblad dynamics, J. Phys. Chem. Lett.15, 3363 (2024)
work page 2024
-
[12]
B. M. Reible, A. Djurdjevac, and L. Delle Site, Chemical potential and variable number of particles control the quantum state: Quantum oscillators as a showcase, APL Quantum2, 016124 (2025)
work page 2025
-
[13]
Morais, Advances and applications of carbon nanotubes, Nanomaterials13, 2674 (2023)
S. Morais, Advances and applications of carbon nanotubes, Nanomaterials13, 2674 (2023)
work page 2023
- [14]
-
[15]
S. Chattaraj and G. Galli, Energy transfer between localized emitters in photonic cavities from first principles, Phys. Rev. Res.7, 033229 (2025)
work page 2025
-
[16]
Lindblad, On the generators of quantum dynamical semigroups, Commun
G. Lindblad, On the generators of quantum dynamical semigroups, Commun. Math. Phys.48, 119–130 (1976). 26
work page 1976
- [17]
-
[18]
Manzano, A short introduction to the Lindblad master equation, AIP Advances10, 025106 (2020)
D. Manzano, A short introduction to the Lindblad master equation, AIP Advances10, 025106 (2020)
work page 2020
-
[19]
R. P. Feynman, Space-time approach to non-relativistic quantum mechanics, Rev. Mod. Phys. 20, 367 (1948)
work page 1948
-
[20]
M. E. Tuckerman,Statistical Mechanics: Theory and Molecular Simulation, 2nd ed. (Oxford University Press, Oxford, 2023)
work page 2023
- [21]
-
[22]
D. Donadio and G. Galli, Thermal conductivity of isolated and interacting carbon nanotubes: Comparing results from molecular dynamics and the Boltzmann transport equation, Phys. Rev. Lett.99, 255502 (2007)
work page 2007
-
[23]
G. Barbalinardo, Z. Chen, H. Dong,et al., Ultrahigh convergent thermal conductivity of carbon nanotubes from comprehensive atomistic modeling, Phys. Rev. Lett.127, 025902 (2021)
work page 2021
-
[24]
G. Benenti, D. Donadio, S. Lepri,et al., Non-Fourier heat transport in nanosystems, Riv. Nuovo Cim.46, 105 (2023)
work page 2023
-
[25]
R. Scipioni, D. Donadio, L. M. Ghiringhelli, and L. Delle Site, Proton wires via one-dimensional water chains adsorbed on metallic steps, J. Chem. Theory Comput.7, 2681 (2011)
work page 2011
-
[26]
D. Donadio and G. Galli, Shedding light on water wires, Physics18, 54 (2025)
work page 2025
-
[27]
J. Paulino, M. Yi, I. Hung, and T. A. Cross, Functional stability of water wire–carbonyl interactions in an ion channel, Proc. Natl. Acad. Sci.117, 11908 (2020)
work page 2020
-
[28]
S. Jang and G. A. Voth, Path integral centroid variables and the formulation of their exact real time dynamics, J. Chem. Phys.111, 2357 (1999)
work page 1999
-
[29]
S. Habershon, D. Manolopoulos, T. E. Markland,et al., Ring-polymer molecular dynamics: Quantum effects in chemical dynamics from classical trajectories in an extended phase space, Annu. Rev. Phys. Chem.64, 387 (2013)
work page 2013
-
[30]
R. P. Feynman, A. R. Hibbs, and D. F. Styer,Quantum Mechanics and Path Integrals, Emended ed. (Dover Publications, Mineola, NY, 2010)
work page 2010
-
[31]
Mazzucchi,Mathematical Feynman Path Integrals and Their Applications, 2nd ed
S. Mazzucchi,Mathematical Feynman Path Integrals and Their Applications, 2nd ed. (World Scientific, Singapore, 2022)
work page 2022
-
[32]
H. F. Trotter, On the product of semi-groups of operators, Proc. Amer. Math. Soc.10, 545 (1959)
work page 1959
-
[33]
Nelson, Feynman integrals and the schrödinger equation, J
E. Nelson, Feynman integrals and the schrödinger equation, J. Math. Phys.5, 332–343 (1964)
work page 1964
-
[34]
P. R. Chernoff, Note on product formulas for operator semigroups, J. Funct. Anal.2, 238 (1968)
work page 1968
-
[35]
M. Reed and B. Simon,Methods of Modern Mathematical Physics. Volume I: Functional 27 Analysis, Revised and Enlarged ed. (Academic Press, San Diego, London, 1980)
work page 1980
-
[36]
Huang,Statistical Mechanics, 2nd ed
K. Huang,Statistical Mechanics, 2nd ed. (John Wiley & Sons, New York, NY, 1987)
work page 1987
- [37]
-
[38]
G. Ciccotti, R. Kapral, and A. Sergi, Non-equilibrium molecular dynamics, inHandbook of Materials Modeling, edited by S. Yip (Springer, Dordrecht, 2005) pp. 745–761
work page 2005
-
[39]
G. Ciccotti and G. Jacucci, Direct computation of dynamical response by molecular dynamics: The mobility of a charged Lennard-Jones particle, Phys. Rev. Lett.35, 789 (1975)
work page 1975
-
[40]
S. Bonella, M. Ferrario, and G. Ciccotti, Thermal diffusion in binary mixtures: Transient behavior and transport coefficients from equilibrium and nonequilib- rium molecular dynamics, Langmuir33, 11281 (2017)
work page 2017
-
[41]
R. Ebrahimi Viand, F. Höfling, R. Klein, and L. Delle Site, Theory and simulation of open systems out of equilibrium, J. Chem. Phys.153, 101102 (2020)
work page 2020
-
[42]
H. Wang, C. Schütte, G. Ciccotti, and L. Delle Site, Exploring the conformational dynamics of alanine dipeptide in solution subjected to an external electric field: A nonequilibrium molecular dynamics simulation, J. Chem. Theroy Comput.10, 1376 (2014)
work page 2014
-
[43]
M. Ferrario, G. Ciccotti, D. Mansutti,et al., A NEMD approach to the melt-front evolution under gravity, J. Chem. Phys.162, 074503 (2025)
work page 2025
-
[44]
G. Ciccotti, G. Jacucci, and I. R. McDonald, “Thought-experiments” by molecular dynamics, J. Stat. Phys.21, 1 (1979)
work page 1979
-
[45]
Pillet, Quantum Dynamical Systems, inOpen Quantum Systems I
C.-A. Pillet, Quantum Dynamical Systems, inOpen Quantum Systems I. The Hamiltonian Approach, Lecture Notes in Mathematics No. 1880, edited by S. Attal, A. Joye, and C.-A. Pillet (Springer, Berlin, Heidelberg, 2006) pp. 107–182
work page 2006
-
[46]
From a structural point of view, the quantum-mechanical argument is analogous to the classical argument of Eq. (10) because the trace gives rise to an inner product on the space of Hilbert- Schmidt operators; see, e.g., Ref. [88, Section II, esp. Eq. (II.13)]
-
[47]
E. B. Davies, Markovian master equations, Commun. Math. Phys.39, 91–110 (1974)
work page 1974
-
[48]
E. B. Davies, Markovian master equations. II, Math. Ann.219, 147–158 (1976)
work page 1976
-
[49]
R. Dümcke and H. Spohn, The proper form of the generator in the weak coupling limit, Z. Phys. B34, 419–422 (1979)
work page 1979
- [50]
-
[51]
F. Benatti and R. Floreanini, Open quantum dynamics: Complete positivity and entanglement, Int. J. Mod. Phys. B19, 3063 (2005)
work page 2005
-
[52]
R. S. Whitney, Staying positive: going beyond Lindblad with perturbative master equations, J. Phys. A: Math. Theor.41, 175304 (2008)
work page 2008
-
[53]
D. Tupkary, A. Dhar, M. Kulkarni,et al., Searching for Lindbladians obeying local conservation 28 laws and showing thermalization, Phys. Rev. A107, 062216 (2023)
work page 2023
-
[54]
F. Campaioli, J. H. Cole, and H. Hapuarachchi, Quantum master equations: Tips and tricks for quantum optics, quantum computing, and beyond, PRX Quantum5, 020202 (2024)
work page 2024
-
[55]
P. Gaspard and M. Nagaoka, Slippage of initial conditions for the Redfield master equation, J. Chem. Phys.111, 5668–5675 (1999)
work page 1999
-
[56]
G. Schaller and T. Brandes, Preservation of positivity by dynamical coarse graining, Phys. Rev. A78, 022106 (2008)
work page 2008
- [57]
-
[58]
F. Benatti, R. Floreanini, and U. Marzolino, Entangling two unequal atoms through a common bath, Phys. Rev. A81, 012105 (2010)
work page 2010
- [59]
- [60]
-
[61]
D. Farina and V. Giovannetti, Open-quantum-system dynamics: Recovering positivity of the Redfield equation via the partial secular approximation, Phys. Rev. A100, 012107 (2019)
work page 2019
-
[62]
R. Hartmann and W. T. Strunz, Accuracy assessment of perturbative master equations: Embracing nonpositivity, Phys. Rev. A101, 012103 (2020)
work page 2020
-
[63]
A. Trushechkin, Unified Gorini-Kossakowski-Lindblad-Sudarshan quantum master equation beyond the secular approximation, Phys. Rev. A103, 062226 (2021)
work page 2021
-
[64]
D. Tupkary, A. Dhar, M. Kulkarni,et al., Fundamental limitations in Lindblad descriptions of systems weakly coupled to baths, Phys. Rev. A105, 032208 (2022)
work page 2022
-
[65]
Diosi, On high-temperature Markovian equation for quantum Brownian motion, EPL22, 1 (1993)
L. Diosi, On high-temperature Markovian equation for quantum Brownian motion, EPL22, 1 (1993)
work page 1993
-
[66]
A. Das, A. Dhar, I. Santra,et al., Quantum Brownian motion: Drude and Ohmic baths as continuum limits of the Rubin model, Phys. Rev. E102, 062130 (2020)
work page 2020
-
[67]
Z. Chen, Y. Lu, H. Wang,et al., Quantum Langevin dynamics for optimization, Commun. Math. Phys.406, 52 (2025)
work page 2025
-
[68]
B. M. Reible and L. Delle Site, Open quantum systems and the grand canonical ensemble, Phys. Rev. E112, 024130 (2025)
work page 2025
-
[69]
L. Delle Site and A. Djurdjevac, An effective Hamiltonian for the simulation of open quantum molecular systems, J. Phys. A: Math. Theor.57, 255002 (2024)
work page 2024
-
[70]
L. Delle Site and C. Hartmann, Computationally feasible bounds for the free energy of nonequi- librium steady states, applied to simple models of heat conduction, Mol. Phys.123, e2391484 (2024)
work page 2024
-
[71]
A. Agarwal and L. Delle Site, Path integral molecular dynamics within the grand canonical-like 29 adaptive resolution technique: Simulation of liquid water, J. Chem. Phys.143, 094102 (2015)
work page 2015
-
[72]
S. Panahian Jand, Z. Nourbakhsh, and L. Delle Site, Nuclear quantum effects in fullerene– fullerene aggregation in water, Front. Chem.10, 1072665 (2022)
work page 2022
-
[73]
M. Kryvohuz, Semiclassical instanton approach to calculation of reaction rate constants in multidimensional chemical systems, J. Chem. Phys.134, 114103 (2011)
work page 2011
-
[74]
K. A. Jung, P. E. Videla, and V. S. Batista, Multi-time formulation of matsubara dynamics, J. Chem. Phys.151, 034108 (2019)
work page 2019
- [75]
-
[76]
Nakajima, On quantum theory of transport phenomena: Steady diffusion, Prog
S. Nakajima, On quantum theory of transport phenomena: Steady diffusion, Prog. Theor. Phys. 20, 948–959 (1958)
work page 1958
-
[77]
Zwanzig, Ensemble method in the theory of irreversibility, J
R. Zwanzig, Ensemble method in the theory of irreversibility, J. Chem. Phys.33, 1338–1341 (1960)
work page 1960
-
[78]
A. G. Redfield, On the theory of relaxation processes, IBM J. Res. Dev.1, 19 (1957)
work page 1957
-
[79]
S.Bonella, M.Monteferrante, C.Pierleoni,et al.,Pathintegralbasedcalculationsofsymmetrized time correlation functions. I., J. Chem. Phys.133, 164104 (2010)
work page 2010
-
[80]
A. Agarwal and L. Delle Site, Grand-canonical adaptive resolution centroid molecular dynamics: Implementation and application, Comput. Phys. Commun.206, 26 (2016)
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.