Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Non-Markovian dynamics of open quantum systems

Christen H. Fleming

Published Jan 1, 2011• 6 citations
OpenAlex

Abstract

An open quantum system is a quantum system that interacts with some environment whose degrees of freedom have been coarse grained away. This model describes non-equilibrium processes more general than scattering-matrix formulations. Furthermore, the microscopically-derived environment provides a model of noise, dissipation and decoherence far more general than Markovian (white noise) models. The latter are fully characterized by Lindblad equations and can be motivated phenomenologically. Non-Markovian processes consistently account for backreaction with the environment and can incorporate effects such as finite temperature and spatial correlations. We consider linear systems with bilinear coupling to the environment, or quantum Brownian motion, and nonlinear systems with weak coupling to the environment. For linear systems we provide exact solutions with analytical results for a variety of spectral densities. Furthermore, we point out an important mathematical subtlety which led to incorrect master-equation coefficients in earlier derivations, given nonlocal dissipation. For nonlinear systems we provide perturbative solutions by translating the formalism of canonical perturbation theory into the context of master equations. It is shown that unavoidable degeneracy causes an unfortunate reduction in accuracy between perturbative master equations and their solutions. We also extend the famous theorem of Lindblad, Gorini, Kossakowski and Sudarshan on completely positivity to non-Markovian master equations. Our application is primarily to model atoms interacting via a common electromagnetic field. The electromagnetic field contains correlations in both space and time, which are related to its relativistic (photon-mediated) nature. As such, atoms residing in the same field experience different environmental effects depending upon their relative position and orientation. Our more accurate solutions were necessary to assess sudden death of entanglement at zero temperature. In contrast to previous claims, we found that all initial states of two-level atoms undergo finite-time disentanglement. We were also able to access regimes which cannot be described by Lindblad equations and other simpler methods, such as near resonance. Finally we revisit the infamous Abraham-Lorentz force, wherein a single particle in motion experiences backreaction from the electromagnetic field. This leads to a number of well-known problems including pre-acceleration and runaway solutions. We found a more a more-suitable open-system treatment of the nonrelativistic particle to be perfectly causal and dissipative without any extraneous requirements for finite size of the particle, weak coupling to the field, etc..

Related papers

The theory of non-Markovian open quantum systems

Cesar Rodriguez

We study the role of correlations with the environment as the source of non-Markovian quantum evolutions. We first focus on the impact that correlations with the environment can have on the dynamical map that evolve the system. We expand the set of initial states of a system and its environment that are known to guarantee completely positive reduced dynamics for the system when the combined state evolves unitarily. We characterize the correlations in the initial state in terms of its quantum discord. The induced maps can be not completely positive when quantum correlations including, but not limited to, entanglement are present. We discuss the implications and limitations of the Markov approximation necessary to derive the Kossakowski-Lindblad master equation. A generalized non-Markovian master equation is derived from the dynamical map of systems correlated with their environment. The physical meaning of not completely positive maps is studied to obtain a consistent theory of non-Markovian quantum dynamics. These are associated to inverse maps necessary to establish correlations and they give rise to a canonical embedding map that is local in time. This master equation goes beyond the Kossakowski-Lindblad master equation. Non-equilibrium quantum thermodynamics can be be studied within this theory. Through out this discussion, the general dynamics of two interacting qubits is used as an example for illustrations.

•May 1, 2008•4 citations

Open Quantum Systems, Decoherence, Atom-Field Coupling

Reinhold A. Bertlmann, Nicolai Friis

Abstract Chapter 22 investigates the dynamics of open quantum systems that interact with their environment. We discuss the time-evolution of the joint system in terms of a time-ordered exponential of an integral over the Hamiltonian and determine a similar formulation for a Liouvillian super-operator. Based on these ideas we develop the notion of dynamical maps and apply the Born-Markov approximation to arrive at a dynamical semi-group and general Markovian master equations. We then turn to a specific example of such an equation of motion for the density operator: the Gorini-Kossakowski-Lindblad-Sudarshan (GKLS) equation which captures both dissipative dynamics and decoherence using Lindblad operators. We employ the GKLS equation to model the absorption, and spontaneous as well as stimulated emission of radiation by a two-level atom. We then discuss the Jaynes-Cummings model that describes the interaction of such an atom with a single mode of the electromagnetic field in a cavity

•Sep 7, 2023•0 citations

Non-Markovian dynamics and entanglement of two-level atoms in a common field

Christen H. Fleming, Nick Cummings, Charis Anastopoulos +1 more

We derive the stochastic equations and consider the non-Markovian dynamics of a system of multiple two-level atoms in a common quantum field. We make only the dipole approximation for the atoms and assume weak atom–field interactions. From these assumptions, we use a combination of non-secular open- and closed-system perturbation theory, and we abstain from any additional approximation schemes. These more accurate solutions are necessary to explore several regimes: in particular, near-resonance dynamics and low-temperature behavior. In detuned atomic systems, small variations in the system energy levels engender timescales which, in general, cannot be safely ignored, as would be the case in the rotating-wave approximation (RWA). More problematic are the second-order solutions, which, as has been recently pointed out (Fleming and Cummings 2011 Phys. Rev. E 83 031117), cannot be accurately calculated using any second-order perturbative master equation, such as RWA, Born–Markov, Redfield, etc. The latter problem, which applies to all perturbative open-system master equations, has a profound effect upon calculation of entanglement at low temperatures. We find that even at zero temperature all initial states will undergo finite-time disentanglement (sometimes termed ‘sudden death’), in contrast to the previous work. We also use our solution, without invoking the RWA, to characterize the necessary conditions for Dicke subradiance at finite temperature. We find that the sub-radiant states fall into two categories at finite temperature: one that is temperature independent and one that acquires temperature dependence. With the RWA, there is no temperature dependence in any case.

arXiv:1012.5067•Journal of Physics A Mathematical and Theoretical•Jan 19, 2012•30 citations

Theoretical and Numerical Studies of the Dynamics of Open Quantum Systems

Souvik Agasti

Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.

JYU dissertations•Jan 1, 2020•0 citations

Simulating Non-Markovian Dynamics in Open Quantum Systems

Meng Xu, Vasilii Vadimov, J. T. Stockburger +1 more

Recent advances in quantum technologies and related experiments have created a need for highly accurate, versatile, and computationally efficient simulation techniques for the dynamics of open quantum systems. Long-lived correlation effects (non-Markovianity), system-environment hybridization, and the necessity for accuracy beyond the Born-Markov approximation form particular challenges. Approaches to meet these challenges have been introduced, originating from different fields, such as hierarchical equations of motion, Lindblad-pseudomode formulas, chain-mapping approaches, quantum Brownian motion master equations, stochastic unravelings, and refined quantum master equations. This diversity, while indicative of the field's relevance, has inadvertently led to a fragmentation that hinders cohesive advances and their effective cross-community application to current problems for complex systems. How are different approaches related to each other? What are their strengths and limitations? Here we give a systematic overview and concise discussion addressing these questions. We make use of a unified framework which very conveniently allows to link different schemes and, this way, may also catalyze further progress. In line with the state of the art, this framework is formulated not in a fully reduced space of the system but in an extended state space which in a minimal fashion includes effective reservoir modes. This in turn offers a comprehensive understanding of existing methods, elucidating their physical interpretations, interconnections, and applicability.

quant-pharXiv:2601.02160•Reviews of Modern Physics•Jan 6, 2026•11 citations

Non-Markovian finite-temperature two-time correlation functions of system operators: Beyond the quantum regression theorem

Hsi‐Sheng Goan, Po‐Wen Chen, Chung-Chin Jian

An extremely useful evolution equation that allows systematically calculating the two-time correlation functions (CF's) of system operators for non-Markovian open (dissipative) quantum systems is derived. The derivation is based on perturbative quantum master equation approach, so non-Markovian open quantum system models that are not exactly solvable can use our derived evolution equation to easily obtain their two-time CF's of system operators, valid to second order in the system-environment interaction. Since the form and nature of the Hamiltonian are not specified in our derived evolution equation, our evolution equation is applicable for bosonic and/or fermionic environments and can be applied to a wide range of system-environment models with any factorized (separable) system-environment initial states (pure or mixed). When applied to a general model of a system coupled to a finite-temperature bosonic environment with a system coupling operator L in the system-environment interaction Hamiltonian, the resultant evolution equation is valid for both L = L(†) and L ≠ L(†) cases, in contrast to those evolution equations valid only for L = L(†) case in the literature. The derived equation that generalizes the quantum regression theorem (QRT) to the non-Markovian case will have broad applications in many different branches of physics. We then give conditions on which the QRT holds in the weak system-environment coupling case and apply the derived evolution equation to a problem of a two-level system (atom) coupled to the finite-temperature bosonic environment (electromagnetic fields) with L ≠ L(†).

arXiv:1103.1502•The Journal of Chemical Physics•Mar 28, 2011•41 citations
An open quantum system is a quantum system that interacts with some environment whose degrees of freedom have been coarse grained away. This model describes non-equilibrium processes more general than scattering-matrix formulations. Furthermore, the microscopically-derived environment provides a model of noise, dissipation and decoherence far more general than Markovian (white noise) models. The latter are fully characterized by Lindblad equations and can be motivated phenomenologically. Non-Markovian processes consistently account for backreaction with the environment and can incorporate effects such as finite temperature and spatial correlations. We consider linear systems with bilinear coupling to the environment, or quantum Brownian motion, and nonlinear systems with weak coupling to the environment. For linear systems we provide exact solutions with analytical results for a variety of spectral densities. Furthermore, we point out an important mathematical subtlety which led to incorrect master-equation coefficients in earlier derivations, given nonlocal dissipation. For nonlinear systems we provide perturbative solutions by translating the formalism of canonical perturbation theory into the context of master equations. It is shown that unavoidable degeneracy causes an unfortunate reduction in accuracy between perturbative master equations and their solutions. We also extend the famous theorem of Lindblad, Gorini, Kossakowski and Sudarshan on completely positivity to non-Markovian master equations. Our application is primarily to model atoms interacting via a common electromagnetic field. The electromagnetic field contains correlations in both space and time, which are related to its relativistic (photon-mediated) nature. As such, atoms residing in the same field experience different environmental effects depending upon their relative position and orientation. Our more accurate solutions were necessary to assess sudden death of entanglement at zero temperature. In contrast to previous claims, we found that all initial states of two-level atoms undergo finite-time disentanglement. We were also able to access regimes which cannot be described by Lindblad equations and other simpler methods, such as near resonance. Finally we revisit the infamous Abraham-Lorentz force, wherein a single particle in motion experiences backreaction from the electromagnetic field. This leads to a number of well-known problems including pre-acceleration and runaway solutions. We found a more a more-suitable open-system treatment of the nonrelativistic particle to be perfectly causal and dissipative without any extraneous requirements for finite size of the particle, weak coupling to the field, etc..
We study the role of correlations with the environment as the source of non-Markovian quantum evolutions. We first focus on the impact that correlations with the environment can have on the dynamical map that evolve the system. We expand the set of initial states of a system and its environment that are known to guarantee completely positive reduced dynamics for the system when the combined state evolves unitarily. We characterize the correlations in the initial state in terms of its quantum discord. The induced maps can be not completely positive when quantum correlations including, but not limited to, entanglement are present. We discuss the implications and limitations of the Markov approximation necessary to derive the Kossakowski-Lindblad master equation. A generalized non-Markovian master equation is derived from the dynamical map of systems correlated with their environment. The physical meaning of not completely positive maps is studied to obtain a consistent theory of non-Markovian quantum dynamics. These are associated to inverse maps necessary to establish correlations and they give rise to a canonical embedding map that is local in time. This master equation goes beyond the Kossakowski-Lindblad master equation. Non-equilibrium quantum thermodynamics can be be studied within this theory. Through out this discussion, the general dynamics of two interacting qubits is used as an example for illustrations.
Abstract Chapter 22 investigates the dynamics of open quantum systems that interact with their environment. We discuss the time-evolution of the joint system in terms of a time-ordered exponential of an integral over the Hamiltonian and determine a similar formulation for a Liouvillian super-operator. Based on these ideas we develop the notion of dynamical maps and apply the Born-Markov approximation to arrive at a dynamical semi-group and general Markovian master equations. We then turn to a specific example of such an equation of motion for the density operator: the Gorini-Kossakowski-Lindblad-Sudarshan (GKLS) equation which captures both dissipative dynamics and decoherence using Lindblad operators. We employ the GKLS equation to model the absorption, and spontaneous as well as stimulated emission of radiation by a two-level atom. We then discuss the Jaynes-Cummings model that describes the interaction of such an atom with a single mode of the electromagnetic field in a cavity
We derive the stochastic equations and consider the non-Markovian dynamics of a system of multiple two-level atoms in a common quantum field. We make only the dipole approximation for the atoms and assume weak atom–field interactions. From these assumptions, we use a combination of non-secular open- and closed-system perturbation theory, and we abstain from any additional approximation schemes. These more accurate solutions are necessary to explore several regimes: in particular, near-resonance dynamics and low-temperature behavior. In detuned atomic systems, small variations in the system energy levels engender timescales which, in general, cannot be safely ignored, as would be the case in the rotating-wave approximation (RWA). More problematic are the second-order solutions, which, as has been recently pointed out (Fleming and Cummings 2011 Phys. Rev. E 83 031117), cannot be accurately calculated using any second-order perturbative master equation, such as RWA, Born–Markov, Redfield, etc. The latter problem, which applies to all perturbative open-system master equations, has a profound effect upon calculation of entanglement at low temperatures. We find that even at zero temperature all initial states will undergo finite-time disentanglement (sometimes termed ‘sudden death’), in contrast to the previous work. We also use our solution, without invoking the RWA, to characterize the necessary conditions for Dicke subradiance at finite temperature. We find that the sub-radiant states fall into two categories at finite temperature: one that is temperature independent and one that acquires temperature dependence. With the RWA, there is no temperature dependence in any case.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.
Recent advances in quantum technologies and related experiments have created a need for highly accurate, versatile, and computationally efficient simulation techniques for the dynamics of open quantum systems. Long-lived correlation effects (non-Markovianity), system-environment hybridization, and the necessity for accuracy beyond the Born-Markov approximation form particular challenges. Approaches to meet these challenges have been introduced, originating from different fields, such as hierarchical equations of motion, Lindblad-pseudomode formulas, chain-mapping approaches, quantum Brownian motion master equations, stochastic unravelings, and refined quantum master equations. This diversity, while indicative of the field's relevance, has inadvertently led to a fragmentation that hinders cohesive advances and their effective cross-community application to current problems for complex systems. How are different approaches related to each other? What are their strengths and limitations? Here we give a systematic overview and concise discussion addressing these questions. We make use of a unified framework which very conveniently allows to link different schemes and, this way, may also catalyze further progress. In line with the state of the art, this framework is formulated not in a fully reduced space of the system but in an extended state space which in a minimal fashion includes effective reservoir modes. This in turn offers a comprehensive understanding of existing methods, elucidating their physical interpretations, interconnections, and applicability.
An extremely useful evolution equation that allows systematically calculating the two-time correlation functions (CF's) of system operators for non-Markovian open (dissipative) quantum systems is derived. The derivation is based on perturbative quantum master equation approach, so non-Markovian open quantum system models that are not exactly solvable can use our derived evolution equation to easily obtain their two-time CF's of system operators, valid to second order in the system-environment interaction. Since the form and nature of the Hamiltonian are not specified in our derived evolution equation, our evolution equation is applicable for bosonic and/or fermionic environments and can be applied to a wide range of system-environment models with any factorized (separable) system-environment initial states (pure or mixed). When applied to a general model of a system coupled to a finite-temperature bosonic environment with a system coupling operator L in the system-environment interaction Hamiltonian, the resultant evolution equation is valid for both L = L(†) and L ≠ L(†) cases, in contrast to those evolution equations valid only for L = L(†) case in the literature. The derived equation that generalizes the quantum regression theorem (QRT) to the non-Markovian case will have broad applications in many different branches of physics. We then give conditions on which the QRT holds in the weak system-environment coupling case and apply the derived evolution equation to a problem of a two-level system (atom) coupled to the finite-temperature bosonic environment (electromagnetic fields) with L ≠ L(†).