Open quantum system
Template:Short description Template:More footnotes needed Template:Use American English In physics, an open quantum system is a quantum mechanical system that interacts with an external quantum system, known as the environment or a bath. In general, these interactions significantly change the dynamics of the system, such that the information contained in the system is lost to its environment. Because no quantum system is completely isolated from its surroundings,<ref>Template:Cite book</ref> it is important to develop a theoretical framework for treating these interactions 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,<ref>Template:Cite journal</ref> quantum measurement theory, quantum statistical mechanics, quantum information science, quantum thermodynamics, quantum cosmology, quantum biology, and semi-classical approximations.<ref>Template:Cite journal</ref>
Quantum system and environment
Open quantum systems are sometimes described by a composite system. The interior system of interest, or principal system, is surrounded by an environment and together they create a closed quantum system.<ref>Template:Cite book</ref> The term "open quantum system" may also be used to refer to the interior system and it may be called the "reduced system" in connection with its reduced density matrix.<ref name=Breuer-2007>Template:Cite book</ref> When the environment has an infinite number of degrees of freedom and thus the modes of the environment form a continuum, it is called a reservoir; this generally results in irreversible dynamics for the interior system. A heat bath or simply bath is a reservoir in thermal equilibrium.<ref name=Breuer-2007/>
Labeling the Hilbert space of the interior system <math>H_S</math> and the environment by <math>H_B</math>, the Hilbert space of the combined system is given by the tensor product <math display="block">H = H_S \otimes H_B.</math> The density matrix for the combined system can be written in terms of the state vectors in the combined Hilbert space, <math>\psi_\alpha</math> as: <math display="block">\rho = \Sigma_\alpha w_\alpha |\psi_\alpha\rangle\langle\psi_\alpha|</math> The normalized positive weights <math>w_\alpha</math> describe a particular configuration of the system or one of a large ensemble of identical systems. The expectation value of a observable Template:Mvar in subspace <math>H_S</math> can be written as <math display="block">\langle A \rangle = \mathrm{tr}_S\{ A \rho_S \}</math> where <math display="block"> \rho_S = \mathrm{tr}_B \rho</math> is the reduced density matrix for system Template:Mvar. This reduced density matrix is primary focus of study for open quantum systems. The symbol <math>\mathrm{tr}_B</math> represents the partial trace over system Template:Mvar and similarly <math>\mathrm{tr}_S</math> represents the partial trace over system Template:Mvar.<ref name=Breuer-2007/>Template:Rp
Models
Recent work has shown that thermodynamic consistency imposes commutativity constraints on the dissipative part of quantum dynamical maps, leading to a general structure for master equations in open systems.<ref>Template:Cite journal</ref> The time evolution of quantum systems can be determined by solving the effective equations of motion, also known as 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.<ref>Template:Citation</ref> 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. 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">Template:Cite journal File:CC-BY icon.svg This article contains quotations from this source, which is available under the Creative Commons Attribution 4.0 International (CC BY 4.0) license.</ref>
It is also common to assume that the time evolution of the system is Markovian, meaning that the way the system's state changes depends only on its current state, with no "memory" of what the state was in the past. 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. This is especially the case of systems where the bath is subjected to a time-dependent external driving force.<ref name="Gamba"/>
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 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 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>Template:Cite book</ref> An alternative approach to fix the Redfield has been proposed by J. Thingna, J.-S. Wang, and Peter Hänggi.<ref>Template:Cite journal</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>Template:Cite journal</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, 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 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>Template:Cite journal</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 nitrogen-vacancy (N-V) center in diamonds. In this example, the color center is the system and the bath consists of carbon-13 (13C) impurities which interact with the system via the magnetic dipole-dipole interaction.<ref>Template:Cite journal</ref>
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 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.
Driven systems
In the presence of time-dependent external fields acting on the system but not on the bath, only the projected dynamics of the system is affected. It contains an additional external force term. The total Hamiltonian of the system is given by:
- <math> H_{T} = H_{\rm S} + H_{\rm B} + H_{\rm SB} + H_{ext,\rm S} </math>
where <math> H_{ext,S}</math> denotes the external field contribution acting on the system only, but not on the bath. If, instead, the time-dependent external field acts on both the system and the bath (environment) as well, the total Hamiltonian is:
- <math> H_{T} = H_{\rm S} + H_{\rm B} + H_{\rm SB} + H_{ext,\rm SB} </math>
where <math> H_{ext,SB}</math> denotes the external field contribution acting on both the system and the bath. Then the fluctuation-dissipation theorem has an extra term with respect to the case where the bath is undriven, which depends on the time-correlation of the external field. As a result, the dynamics is intrinsically non-Markovian even if the system and the bath are fully Markovian when the external field is switched off.<ref name="Gamba">Template:Cite journal</ref>
History
Template:Primary sources In 1927, John von Neumann introduced the concept of density matrices to describe the statistical state of a quantum system, providing a mathematical framework for handling mixed states due to environmental interactions.<ref>Template:Cite journal</ref><ref name="Schlosshauer2007">Template:Cite journal</ref> This formalism became a cornerstone for studying open quantum systems, allowing for the description of mixed states arising from system-environment coupling.<ref>Template:Cite journal</ref>
In the 1940s, Felix Bloch developed an approach to understand nuclear spin systems in a magnetic field, as in nuclear magnetic resonance (NMR).<ref>Template:Cite journal</ref><ref name="Chruscinski2022">Template:Cite journal</ref> His work laid the groundwork for studying relaxation and decoherence processes, introducing the first quantum master equation.<ref>Template:Cite journal</ref> This was later generalized to treat generic two-state quantum systems.<ref>Template:Cite journal</ref>
In the 1950s and 1960s, the quantum master equation was developed to describe the time evolution of open quantum systems. Alfred Redfield introduced the Redfield equation in 1957, a Markovian master equation for weakly coupled systems.<ref>Template:Cite journal</ref><ref name="Chruscinski2022"/> However, it does not always guarantee density matrix positivity.<ref>Template:Cite book</ref>
In 1970, H. Dieter Zeh formalized decoherence in his seminal paper.<ref>Template:Cite journal</ref><ref name="Schlosshauer2007"/> In the 1970s, the Lindblad equation, independently proposed by Göran Lindblad, Andrzej Kossakowski, and George Sudarshan (GKSL equation), provided a general Markovian framework ensuring complete positivity and trace preservation.<ref>Template:Cite journal</ref><ref>Template:Cite journal</ref><ref name="Bassi2017">Template:Cite journal</ref>
In 1981, Amir Caldeira and Anthony J. Leggett introduced the Caldeira–Leggett model, modeling the environment as harmonic oscillators linearly coupled to the system.<ref>Template:Cite journal</ref> This model is widely used for studying quantum dissipation and decoherence.<ref>Template:Cite journal</ref>
Theoretical methods for non-Markovian dynamics, accounting for environmental memory effects, were developed in the 1950s and 1960s through the Nakajima–Zwanzig formalism.<ref>Template:Cite journal</ref><ref>Template:Cite journal</ref><ref name="Chruscinski2022"/> These approaches gained broader application in the 1980s and 1990s, advancing the understanding of quantum decoherence and entanglement in strongly coupled systems.<ref>Template:Cite journal</ref>
See also
References
<references></references>
Further reading
- Template:Cite book
- Template:Cite book
- Template:Cite book
- Template:Cite book
- Template:Cite book
- Template:Cite book
- Template:Cite journal
- Template:Cite book
- Template:Cite journal
- Template:Cite book
- Template:Cite book
- Template:Cite book