View source for Open quantum system

{{Use American English|date=January 2019}}{{Short description|Quantum mechanical system that interacts with a quantum-mechanical environment}}
{{In line citation|date=January 2019}}
In [[physics]], an '''open quantum system''' is a [[Quantum mechanics|quantum]]-mechanical system that interacts with an external [[quantum system]], which is known as the ''environment'' or a ''bath''. In general, these interactions significantly change the dynamics of the system and result in [[quantum dissipation]], such that the information contained in the system is lost to its environment. Because no quantum system is completely isolated from its surroundings,<ref>{{cite book |first1=H.-P. |last1=Breuer |first2=F. |last2=Petruccione |title=The Theory of Open Quantum Systems |publisher=Oxford University Press |date=2007 |page=vii |quote="Quantum mechanical systems must be considered as open systems"}}</ref> it is important to develop a theoretical framework for treating these interactions in order to obtain an accurate understanding of quantum systems.

Techniques developed in the context of open quantum systems have proven powerful in fields such as [[quantum optics]], [[Measurement in quantum mechanics|quantum measurement theory]], quantum [[statistical mechanics]], [[quantum information]] science, [[quantum thermodynamics]], [[quantum cosmology]], [[quantum biology]], and semi-classical approximations.

==Quantum system and environment==

A complete description of a quantum system requires the inclusion of the environment. Completely describing the resulting combined system then requires the inclusion of its environment, which results in a new system that can only be completely described if its environment is included and so on. The eventual outcome of this process of embedding is the state of the whole universe described by a [[wavefunction]] <math>\Psi</math>. The fact that every quantum system has some degree of openness also means that no quantum system can ever be in a [[pure state]].
[[File:System-bath.jpg|thumb|System bath partition]]
Even if the combined system is in a pure state and can be described by a wavefunction <math> \Psi </math>, a subsystem in general cannot be described by a wavefunction. This observation motivated the formalism of [[density matrices]], or density operators, introduced by [[John von Neumann]]<ref>{{Citation|last=von Neumann|first=John|title=Wahrscheinlichkeitstheoretischer Aufbau der Quantenmechanik|journal=Göttinger Nachrichten|volume=1|pages=245–272|year=1927|authorlink=John von Neumann}}</ref> in 1927 and independently, but less systematically by [[Lev Landau]] in 1927 and [[Felix Bloch]] in 1946. In general, the state of a subsystem is described by the density operator <math> \rho </math> and the expectation value of an observable <math> A </math> by the scalar product <math> (\rho \cdot A) = \rm{tr}\{ \rho A \} </math>. There is no way to know if the combined system is pure from the knowledge of observables of the subsystem alone. In particular, if the combined system has [[quantum entanglement]], the state of the subsystem is not pure.

== Dynamics ==
In general, the time evolution of closed quantum systems is described by unitary operators acting on the system.  For open systems, however, the interactions between the system and its environment make it so that the dynamics of the system cannot be accurately described using unitary operators alone.

The time evolution of quantum systems can be determined by solving the effective equations of motion, also known as [[Master equation|master equations]], that govern how the density matrix describing the system changes over time and the dynamics of the observables that are associated with the system.  In general, however, the environment that we want to model as being a part of our system is very large and complicated, which makes finding exact solutions to the master equations difficult, if not impossible.  As such, the theory of open quantum systems seeks an economical treatment of the dynamics of the system and its observables.  Typical observables of interest include things like energy and the robustness of [[quantum coherence]] (i.e. a measure of a state's coherence).  Loss of energy to the environment is termed [[quantum dissipation]], while loss of coherence is termed [[quantum decoherence]].

Due to the difficulty of determining the solutions to the master equations for a particular system and environment, a variety of techniques and approaches have been developed. A common objective is to derive a reduced description wherein the system's dynamics are considered explicitly and the bath's dynamics are described implicitly. The main assumption is that the entire system-environment combination is a large closed system.  Therefore, its time evolution is governed by a [[unitary transformation]] generated by a global [[Hamiltonian (quantum mechanics)|Hamiltonian]]. For the combined system bath scenario the global Hamiltonian can be decomposed into:

:<math> H=H_{\rm S}+H_{\rm B}+H_{\rm SB} </math>

where <math>H_{\rm S}</math> is the system's Hamiltonian, <math>H_{\rm B} </math> is the bath Hamiltonian and <math>H_{\rm SB}</math> is the system-bath interaction. The state of the system can then be obtained from a partial trace over the combined system and bath: <math>\rho_{\rm S} (t) =\rm{tr}_{\rm B} \{\rho_{SB} (t)\} </math>.<ref name="Kosloff20132">{{cite journal|last=Kosloff|first=Ronnie|year=2013|title=Quantum Thermodynamics: A Dynamical Viewpoint|journal=Entropy|volume=15|issue=6|pages=2100–2128|arxiv=1305.2268|bibcode=2013Entrp..15.2100K|doi=10.3390/e15062100|issn=1099-4300|doi-access=free}}  [[File:CC-BY_icon.svg|50x50px]] This article contains quotations from this source, which is available under the [[creativecommons:by/4.0/|Creative Commons Attribution 4.0 International (CC BY 4.0)]] license.</ref>

Another common assumption that is used to make systems easier to solve is the assumption that the state of the system at the next moment depends only on the current state of the system.  in other words, the system doesn't have a memory of its previous states. Systems that have this property are known as [[Markovian property|Markovian]] systems. This approximation is justified when the system in question has enough time for the system to relax to equilibrium before being perturbed again by interactions with its environment.  For systems that have very fast or very frequent perturbations from their coupling to their environment, this approximation becomes much less accurate.

=== Markovian equations ===
When the interaction between the system and the environment is weak, a time-dependent [[perturbation theory]] seems appropriate for treating the evolution of the system. In other words, if the interaction between the system and its environment is weak, then any changes to the combined system over time can be approximated as originating from only the system in question. Another typical assumption is that the system and bath are initially uncorrelated <math> \rho(0)=\rho_{\rm S} \otimes \rho_{\rm B} </math>. This idea originated with [[Felix Bloch]] and was expanded upon by Alfred Redfield in his derivation of the [[Redfield equation]]. The Redfield equation is a Markovian master equation that describes the time evolution of the density matrix of the combined system. The drawback of the Redfield equation is that it does not conserve the [[Positive element|positivity]] of the density operator.

A formal construction of a local equation of motion with a [[Markovian property]] is an alternative to a reduced derivation. The theory is based on an axiomatic approach. The basic starting point is a [[completely positive map]]. The assumption is that the initial system-environment state is uncorrelated <math> \rho(0)=\rho_{\rm S} \otimes \rho_{\rm B} </math> and the combined dynamics is generated by a [[unitary operator]]. Such a map falls under the category of [[Kraus operator]]. The most general type of a time-homogeneous master equation with the [[Markovian property]] describing non-unitary evolution of the density matrix ρ that is trace-preserving and completely positive for any initial condition is the Gorini–Kossakowski–Sudarshan–[[Lindblad equation]] or GKSL equation:
:<math>\dot\rho_{\rm S}=-{i\over\hbar}[H_{\rm S},\rho_{\rm S}]+{\cal L}_{\rm D}(\rho_{\rm S}) </math>

<math> H_{\rm S}</math> is a ([[Hermitian]]) [[Hamiltonian (quantum mechanics)|Hamiltonian]] part and <math>{\cal L}_{\rm D}</math>:
:<math>{\cal L}_{\rm D}(\rho_{\rm S})=\sum_n \left(V_n\rho_{\rm S} V_n^\dagger-\frac{1}{2}\left(\rho_{\rm S} V_n^\dagger V_n + V_n^\dagger V_n\rho_{\rm S}\right)\right)</math>
is the dissipative part describing implicitly through system operators <math> V_n </math> the influence of the bath on the system.
The [[Markov property]] imposes that the system and bath are uncorrelated at all times <math> \rho_{\rm SB}=\rho_{\rm S} \otimes \rho_{\rm B} </math>.
The GKSL equation is unidirectional and leads any initial state <math> \rho_{\rm S}</math> to a steady state solution which is an invariant of the equation of motion <math> \dot \rho_{\rm S}(t \rightarrow \infty ) = 0 </math>.
The family of maps generated by the GKSL equation forms a [[Quantum dynamical semigroup]]. In some fields, such as [[quantum optics]], the term [[Lindblad superoperator]] is often used to express the quantum master equation for a dissipative system. E.B. Davis derived the GKSL with [[Markovian property]] master equations using [[perturbation theory]] and additional approximations, such as the rotating wave or secular, thus fixing the flaws of the [[Redfield equation]]. Davis construction is consistent with the Kubo-Martin-Schwinger stability criterion for thermal equilibrium i.e. the [[KMS state]].<ref>{{cite book|first=Heinz-Peter|last=Breuer|author2=F. Petruccione|title=The Theory of Open Quantum Systems|publisher=Oxford University Press|year=2007|isbn=978-0-19-921390-0}}</ref> An alternative approach to fix the Redfield has been proposed by J. Thingna, J.-S. Wang, and P. Hänggi<ref>{{Cite journal|last1=Thingna|first1=Juzar|last2=Wang|first2=Jian-Sheng|last3=Hänggi|first3=Peter|date=2012-05-21|title=Generalized Gibbs state with modified Redfield solution: Exact agreement up to second order|journal=The Journal of Chemical Physics|language=en|volume=136|issue=19|pages=194110|doi=10.1063/1.4718706|pmid=22612083|issn=0021-9606|arxiv=1203.6207|bibcode=2012JChPh.136s4110T|s2cid=7014354}}</ref> that allows for system-bath interaction to play a role in equilibrium differing from the KMS state.

In 1981, [[Amir Caldeira]] and [[Anthony J. Leggett]] proposed a simplifying assumption in which the bath is decomposed to normal modes represented as harmonic oscillators linearly coupled to the system.<ref>A. Caldeira and A. J. Leggett, ''Influence of dissipation on quantum tunneling in macroscopic systems'', Physical Review Letters, vol. 46, p. 211, 1981.</ref> As a result, the influence of the bath can be summarized by the bath spectral function. This method is known as the [[Caldeira-Leggett model|Caldeira–Leggett model]], or harmonic bath model. To proceed and obtain explicit solutions, the [[path integral formulation]] description of [[quantum mechanics]] is typically employed.  A large part of the power behind this method is the fact that harmonic oscillators are relatively well-understood compared to the true coupling that exists between the system and the bath.  Unfortunately, while the Caldeira-Leggett model is one that leads to a physically consistent picture of quantum dissipation, its [[Ergodicity|ergodic]] properties are too weak and so the dynamics of the model do not generate wide-scale [[quantum entanglement]] between the bath modes.

An alternative bath model is a spin bath.<ref>{{Cite journal|last1=Prokof'ev|first1=N. V.|last2=Stamp|first2=P. C. E.|date=2000|title=Theory of the spin bath|journal=Reports on Progress in Physics|language=en|volume=63|issue=4|pages=669|doi=10.1088/0034-4885/63/4/204 | arxiv=cond-mat/0001080|bibcode=2000RPPh...63..669P |s2cid=55075035|issn=0034-4885}}</ref> At low temperatures and weak system-bath coupling, the Caldeira-Leggett and spin bath models are equivalent. But for higher temperatures or strong system-bath coupling, the spin bath model has strong ergodic properties. Once the system is coupled, significant entanglement is generated between all modes. In other words, the spin bath model can simulate the Caldeira-Leggett model, but the opposite is not true.

An example of natural system being coupled to a spin bath is a [[N-V center|nitrogen-vacancy (N-V) center]] in diamonds. In this example, the color center is the system and the bath consists of [[carbon-13]] (<sup>13</sup>C) impurities which interact with the system via the magnetic dipole-dipole [http://iopscience.iop.org/article/10.1088/1367-2630/18/9/093001/meta interaction].

For open quantum systems where the bath has oscillations that are particularly fast, it is possible to average them out by looking at sufficiently large changes in time. This is possible because the average amplitude of fast oscillations over a large time scale is equal to the central value, which can always be chosen to be zero with a minor shift along the vertical axis. This method of simplifying problems is known as the secular approximation.

=== Non-Markovian equations ===
Open quantum systems that do not have the Markovian property are generally much more difficult to solve. This is largely due to the fact that the next state of a non-Markovian system is determined by each of its previous states, which rapidly increases the memory requirements to compute the evolution of the system. Currently, the methods of treating these systems employ what are known as [[Projection (linear algebra)|projection operator]] techniques. These techniques employ a projection operator <math>\mathcal{P}</math>, which effectively applies the trace over the environment as described previously.  The result of applying <math>\mathcal{P}</math> to <math>\rho</math>(i.e. calculating <math>\mathcal{P}\rho</math>) is called the ''relevant part'' of <math>\rho</math>.  For completeness, another operator <math>\mathcal{Q}</math> is defined so that <math>\mathcal{P}+\mathcal{Q}=\mathcal{I}</math> where <math>\mathcal{I}</math> is the identity matrix.  The result of applying <math>\mathcal{Q}</math> to <math>\rho</math>(i.e. calculating <math>\mathcal{Q}\rho</math>) is called the ''irrelevant part'' of <math>\rho</math>.  The primary goal of these methods is to then derive a master equation that defines the evolution of <math>\mathcal{P}\rho</math>.

One such derivation using the projection operator technique results in what is known as the [[Nakajima–Zwanzig equation]]. This derivation highlights the problem of the reduced dynamics being non-local in time:

: <math>\partial_t{\rho }_\mathrm{S}=\mathcal{P}{\cal L}{{\rho}_\mathrm{S}}+\int_{0}^{t}{dt'\mathcal{K}({t}'){{\rho }_\mathrm{S}}(t-{t}')}.</math>

Here the effect of the bath throughout the time evolution of the system is hidden in the memory kernel <math> \kappa (\tau)</math>. While the Nakajima-Zwanzig equation is an exact equation that holds for almost all open quantum systems and environments, it can be very difficult to solve.  This means that approximations generally need to be introduced to reduce the complexity of the problem into something more manageable.  As an example, the assumption of a fast bath is required to lead to a time local equation: <math> \partial_t \rho_S = {\cal L } \rho_S </math>. Other examples of valid approximations include the weak-coupling approximation and the single-coupling approximation.

In some cases, the projection operator technique can be used to reduce the dependence of the system's next state on all of its previous states.  This method of approaching open quantum systems is known as the time-convolutionless projection operator technique, and it is used to generate master equations that are inherently local in time. Because these equations can neglect more of the history of the system, they are often easier to solve than things like the Nakajima-Zwanzig equation.

Another approach emerges as an analogue of classical dissipation theory developed by [[Ryogo Kubo]] and Y. Tanimura. This approach is connected to [[hierarchical equations of motion]] which embed the density operator in a larger space of auxiliary operators such that a time local equation is obtained for the whole set and their memory is contained in the auxiliary operators.

==See also==

* [[Lindblad equation]]
* [[Markov property]]
* [[Master equation]]
* [[Quantum thermodynamics]]

==References==
<references group=""></references>

=== Unclassified references ===
* {{cite book|title=Quantum Theory and Its Stochastic Limit|last1=Accardi|first1=Luigi|last2=Lu|first2=Yun Gang|last3=Volovich|first3=I.V.|publisher=Springer Verlag|year=2002|isbn=978-3-540-41928-0|location=New York}}
* {{cite book|title=Quantum Dynamical Semigroups and Applications|last1=Alicki|first1=Robert|last2=Lendi|first2=Karl|publisher=Springer Verlag|year=1987|isbn=978-0-387-18276-6|location=Berlin}}
* {{cite book|title=Open Quantum Systems II: The Markovian Approach|last1=Attal|first1=Stéphane|last2=Joye|first2=Alain|last3=Pillet|first3=Claude-Alain|publisher=Springer|year=2006|isbn=978-3-540-30992-5}}
* {{cite book|title=Quantum Theory of Open Systems|last=Davies|first=Edward Brian|publisher=Academic Press|year=1976|isbn=978-0-12-206150-9|location=London}}
* {{cite book|title=Information Dynamics and Open Systems: Classical and Quantum Approach|last1=Ingarden|first1=Roman S.|last2=Kossakowski|first2=A.|last3=Ohya|first3=M.|publisher=Springer Verlag|year=1997|isbn=978-0-7923-4473-5|location=New York}}
* {{cite book|url=https://books.google.com/books?id=tvMcxdPNmRAC&q=.+Non-Equilibrium+Entropy+and+Irreversibility|title=Non-Equilibrium Entropy and Irreversibility|last=Lindblad|first=G.|publisher=Delta Reidel|year=1983|isbn=978-1-4020-0320-2|location=Dordrecht}}
* {{cite journal|last1=Okolowicz|first1=J.|last2=Płoszajczak|first2=M.|last3=Nazarewicz|first3=W.|year=2012|title=On the Origin of Nuclear Clustering|journal=Progress of Theoretical Physics Supplement|volume=196|pages=230–243|arxiv=1202.6290|bibcode=2012PThPS.196..230O|doi=10.1143/PTPS.196.230|s2cid=119109268}}
* {{cite book|title=Quantum Mechanics of Non-Hamiltonian and Dissipative Systems|last=Tarasov|first=Vasily E.|publisher=Elsevier Science|year=2008|isbn=978-0-08-055971-1|location=Amsterdam, Boston, London, New York}}
* {{cite book|title=Quantum Dissipative Systems|last=Weiss|first=Ulrich|publisher=World Scientific|year=2012|isbn=978-981-4374-91-0|edition=4th}}
* {{cite book|title=Quantum Measurement and Control|last1=Wiseman|first1=Howard M.|last2=Milburn|first2=Gerard J.|date=2010|publisher=Cambridge University Press|isbn=978-0-521-80442-4}}

== External links ==

* {{Wikiversity-inline|Open Quantum Systems}}
{{DEFAULTSORT:Open Quantum System}}
[[Category:Quantum mechanics]] 
[[Category:Physical systems]]