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

arXiv:quant-ph/0311077 · Physical Review A

Is the dynamics of open quantum systems always linear?

K. M. Fonseca-Romero, Peter Talkner A, Peter Hänggi A

Published May 17, 2004• 59 citations
arXivHTMLPDF

Abstract

We study the influence of the preparation of an open quantum system on its reduced time evolution. In contrast to the frequently considered case of an initial preparation where the total density matrix factorizes into a product of a system density matrix and a bath density matrix the time evolution generally is no longer governed by a linear map nor is this map affine. Put differently, the evolution is truly nonlinear and cannot be cast into the form of a linear map plus a term that is independent of the initial density matrix of the open quantum system. As a consequence, the inhomogeneity that emerges in formally exact generalized master equations is in fact a nonlinear term that vanishes for a factorizing initial state. The general results are elucidated with the example of two interacting spins prepared at thermal equilibrium with one spin subjected to an external field. The second spin represents the environment. The field allows the preparation of mixed density matrices of the first spin that can be represented as a convex combination of two limiting pure states, i.e., the preparable reduced density matrices make up a convex set. Moreover, the map from these reduced density matrices onto the corresponding density matrices of the total system is affine only for vanishing coupling between the spins. In general, the set of the accessible total density matrices is nonconvex.

DOI
10.1103/physreva.69.052109

Related papers

Deviation from linear reduced dynamics always occurs for each non-factorisable system-environment unitary evolution

Iman Sargolzahi

In the simplest approximation, the reduced dynamics of a quantum system $S$ interacting with its environment $E$ is considered to be given by a completely positive map. But, in general, this is not the case. In fact, the reduced dynamics of the system in not even linear, in general. Whether the reduced dynamics is linear or not is determined by two factors: the set of possible initial states of the system-environment $\mathcal{S}=\left\lbrace ρ_{SE} \right\rbrace $, and the joint system-environment unitary evolution $U$. When $U$ is factorisable as $U=U_S\otimes U_E$, then we can choose $\mathcal{S}=\mathcal{D}$, where $\mathcal{D}$ is the set of all system-environment density operators. In other words, when $U$ is factorisable, the reduced dynamics of the system $S$ is linear (in fact unitary) for arbitrary initial state of the system-environment $ρ_{SE}$. We show that this result cannot be generalized to any non-factorisable $U$: For any non-factorisable unitary evolution of the whole system-environment $U$, the set $\mathcal{S}$ must be chosen as a proper subset of $\mathcal{D}$ to achieve linear reduced dynamics. As a byproduct, considering a convex set of possible initial states of the system-environment $\mathcal{S}$ such that $\mathrm{Tr}_{E} \ \mathcal{S}=\mathrm{Tr}_{E} \ \mathcal{D}$, we show that when the reduced dynamics of the system, for one system-environment unitary evolution $U_1$, is positive, but not completely positive, this implies that reduced dynamics is not linear for another $U_2$.

quant-pharXiv:2608.18977•Aug 20, 2026•0 citations

Mapping the Schrödinger picture of open quantum dynamics

Thomas F. Jordan, Anil Shaji, E. C. G. Sudarshan

For systems described by finite matrices, an affine form is developed for the maps that describe evolution of density matrices for a quantum system that interacts with another. This is established directly from the Heisenberg picture. It separates elements that depend only on the dynamics from those that depend on the state of the two systems. While the equivalent linear map is generally not completely positive, the homogeneous part of the affine maps is, and is shown to be composed of multiplication operations that come simply from the Hamiltonian for the larger system. The inhomogeneous part is shown to be zero if and only if the map does not increase the trace of the square of any density matrix. Properties are worked out in detail for two-qubit examples.

arXiv:quant-ph/0505123•Physical Review A•Jan 18, 2006•20 citations

Decoherence in open quantum systems: influence of the intrinsic bath dynamics

V. V. Ignatyuk, V. G. Morozov

The non-Markovian master equation for open quantum systems is obtained by generalization of the standard Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations is written for the reduced density matrices of the bath ϱB(t) and of the system ϱS(t). A formal solution of the equation for ϱB(t) in the 2nd approximation in interaction yields a specific extra term related to the intrinsic bath dynamics. This term is nonlinear in the reduced density matrix ϱS(t), and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we apply the generalized ZN projection scheme to a simple dephasing model. We study the obtained kinetic equation both in the Markovian approximation and beyond it (for the term related to the intrinsic bath dynamics) and compare the results with the exact ones.

arXiv:2112.05595•Condensed Matter Physics•Jan 1, 2022•4 citations

Influence of initial correlations on evolution over time of an open quantum system

Victor F. Los

A novel approach to accounting for the influence of initial system-bath correlations on the dynamics of an open quantum system, based on the conventional projection operator technique, is suggested. To avoid the difficulties of treating the initial correlations, the conventional Nakajima-Zwanzig inhomogeneous generalized master equations (GMEs) for a system's reduced statistical operator and correlation function are exactly converted into the homogeneous GMEs (HGMEs), which take into account the initial correlations in the kernel governing the evolution of these HGMEs. In the second order (Born) approximation in the system-bath interaction, the obtained HGMEs are local in time and valid at all timescales. They are further specialized for a realistic equilibrium Gibbs initial (at $t=t_0$) system+bath state (for a system reduced statistical operator an external force at $t>t_0$ is applied) and then for a bath of oscillators (Boson field). As an example, the evolution of a selected quantum oscillator (a localized mode) interacting with a Boson field (Fano-like model) is considered at different timescales. It is shown explicitly how the initial correlations influence the oscillator evolution process. In particular, it is shown that the equilibrium system's correlation function acquires at the large timescale the additional constant phase factor conditioned by survived initial system-bath correlations.

arXiv:2405.07299•May 12, 2024•0 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

Master equation for open quantum systems: Zwanzig-Nakajima projection technique and the intrinsic bath dynamics

V. V. Ignatyuk, V. G. Morozov

The non-Markovian master equation for open quantum systems is obtained by generalization of the ordinary Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations for the reduced density matrices of the bath $\varrho_{B}(t)$ and the system $\varrho_{S}(t)$ are written down. A formal solution of the equation for $\varrho_{B}(t)$, having been inserted in the equation for the reduced density matrix of the system, in the 2-nd approximation in interaction yields a very specific extra term in the generalized master equation. This term, being nonlinear in $\varrho_{S}(t)$, is related to the intrinsic bath dynamics and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we applied the generalized ZN projection scheme to a simple dephasing model. It is shown that consideration of the lowest order in interaction is insufficient to describe time evolution of the system coherence adequately. We explain this fact by analyzing the exact and approximate forms of $\varrho_{B}(t)$ and give some hints how to take the dynamic correlations (which originate from the spin-bath coupling) into account.

arXiv:1911.12312•Nov 27, 2019•0 citations
We study the influence of the preparation of an open quantum system on its reduced time evolution. In contrast to the frequently considered case of an initial preparation where the total density matrix factorizes into a product of a system density matrix and a bath density matrix the time evolution generally is no longer governed by a linear map nor is this map affine. Put differently, the evolution is truly nonlinear and cannot be cast into the form of a linear map plus a term that is independent of the initial density matrix of the open quantum system. As a consequence, the inhomogeneity that emerges in formally exact generalized master equations is in fact a nonlinear term that vanishes for a factorizing initial state. The general results are elucidated with the example of two interacting spins prepared at thermal equilibrium with one spin subjected to an external field. The second spin represents the environment. The field allows the preparation of mixed density matrices of the first spin that can be represented as a convex combination of two limiting pure states, i.e., the preparable reduced density matrices make up a convex set. Moreover, the map from these reduced density matrices onto the corresponding density matrices of the total system is affine only for vanishing coupling between the spins. In general, the set of the accessible total density matrices is nonconvex.
In the simplest approximation, the reduced dynamics of a quantum system SSS interacting with its environment EEE is considered to be given by a completely positive map. But, in general, this is not the case. In fact, the reduced dynamics of the system in not even linear, in general. Whether the reduced dynamics is linear or not is determined by two factors: the set of possible initial states of the system-environment $\mathcal{S}=\left\lbrace ρ_{SE} \right\rbrace $, and the joint system-environment unitary evolution UUU. When UUU is factorisable as U=US⊗UEU=U_S\otimes U_EU=US​⊗UE​, then we can choose S=D\mathcal{S}=\mathcal{D}S=D, where D\mathcal{D}D is the set of all system-environment density operators. In other words, when UUU is factorisable, the reduced dynamics of the system SSS is linear (in fact unitary) for arbitrary initial state of the system-environment ρSEρ_{SE}ρSE​. We show that this result cannot be generalized to any non-factorisable UUU: For any non-factorisable unitary evolution of the whole system-environment UUU, the set S\mathcal{S}S must be chosen as a proper subset of D\mathcal{D}D to achieve linear reduced dynamics. As a byproduct, considering a convex set of possible initial states of the system-environment S\mathcal{S}S such that TrE S=TrE D\mathrm{Tr}_{E} \ \mathcal{S}=\mathrm{Tr}_{E} \ \mathcal{D}TrE​ S=TrE​ D, we show that when the reduced dynamics of the system, for one system-environment unitary evolution U1U_1U1​, is positive, but not completely positive, this implies that reduced dynamics is not linear for another U2U_2U2​.
For systems described by finite matrices, an affine form is developed for the maps that describe evolution of density matrices for a quantum system that interacts with another. This is established directly from the Heisenberg picture. It separates elements that depend only on the dynamics from those that depend on the state of the two systems. While the equivalent linear map is generally not completely positive, the homogeneous part of the affine maps is, and is shown to be composed of multiplication operations that come simply from the Hamiltonian for the larger system. The inhomogeneous part is shown to be zero if and only if the map does not increase the trace of the square of any density matrix. Properties are worked out in detail for two-qubit examples.
The non-Markovian master equation for open quantum systems is obtained by generalization of the standard Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations is written for the reduced density matrices of the bath ϱB(t) and of the system ϱS(t). A formal solution of the equation for ϱB(t) in the 2nd approximation in interaction yields a specific extra term related to the intrinsic bath dynamics. This term is nonlinear in the reduced density matrix ϱS(t), and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we apply the generalized ZN projection scheme to a simple dephasing model. We study the obtained kinetic equation both in the Markovian approximation and beyond it (for the term related to the intrinsic bath dynamics) and compare the results with the exact ones.
A novel approach to accounting for the influence of initial system-bath correlations on the dynamics of an open quantum system, based on the conventional projection operator technique, is suggested. To avoid the difficulties of treating the initial correlations, the conventional Nakajima-Zwanzig inhomogeneous generalized master equations (GMEs) for a system's reduced statistical operator and correlation function are exactly converted into the homogeneous GMEs (HGMEs), which take into account the initial correlations in the kernel governing the evolution of these HGMEs. In the second order (Born) approximation in the system-bath interaction, the obtained HGMEs are local in time and valid at all timescales. They are further specialized for a realistic equilibrium Gibbs initial (at t=t0t=t_0t=t0​) system+bath state (for a system reduced statistical operator an external force at t>t0t>t_0t>t0​ is applied) and then for a bath of oscillators (Boson field). As an example, the evolution of a selected quantum oscillator (a localized mode) interacting with a Boson field (Fano-like model) is considered at different timescales. It is shown explicitly how the initial correlations influence the oscillator evolution process. In particular, it is shown that the equilibrium system's correlation function acquires at the large timescale the additional constant phase factor conditioned by survived initial system-bath correlations.
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.
The non-Markovian master equation for open quantum systems is obtained by generalization of the ordinary Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations for the reduced density matrices of the bath ϱB(t)\varrho_{B}(t)ϱB​(t) and the system ϱS(t)\varrho_{S}(t)ϱS​(t) are written down. A formal solution of the equation for ϱB(t)\varrho_{B}(t)ϱB​(t), having been inserted in the equation for the reduced density matrix of the system, in the 2-nd approximation in interaction yields a very specific extra term in the generalized master equation. This term, being nonlinear in ϱS(t)\varrho_{S}(t)ϱS​(t), is related to the intrinsic bath dynamics and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we applied the generalized ZN projection scheme to a simple dephasing model. It is shown that consideration of the lowest order in interaction is insufficient to describe time evolution of the system coherence adequately. We explain this fact by analyzing the exact and approximate forms of ϱB(t)\varrho_{B}(t)ϱB​(t) and give some hints how to take the dynamic correlations (which originate from the spin-bath coupling) into account.