Open Quantum Systems/The Lindblad Form

Resource type: this resource contains a lecture or lecture notes.

Superoperators and dynamical maps

edit

The Liouville von Neumann equation is given by

 

We can define a superoperator   such that  . It is called a superoperator because it is an object that acts on an operator and results in a new operator. If the Hamiltonian is time-independent, we may formally integrate the Liouville von Neumann equation and obtain

 

where   is another superoperator that maps the density matrix from its initial form to its form at time   and therefore is called a dynamical map. It is related to the unitary evolution operator   according to

 

Quantum dynamical semigroups

edit

Now let us include also interaction between the system of interest and its environment. In the following, we will use   for the reduced density operator for the system and   for the partial trace over the environment. As the dynamics of the combination of system and environment is fully coherent, we have

 

which after taking the trace over the environment on both sides results in

 

In many typical situations, the initial state between the system and the environment is a product state of the form

 

Then, we may again think of the right hand side of the previous equation to define a superoperator representing a dynamical map  , but now for   alone! Furthermore, we may use the following decompositions [1]

 

where the operators   act only on the Hilbert space of the system and are given by

 

From the completeness of the states   on the Hilbert space of the environment, we may identify the relation

 

from which follows

 

i.e., the dynamical map is trace-preserving. Moreover, it is completely positive, mapping a positive density matrix onto another positive density matrix.[1]

In many important cases, we can make one further assumption on the dynamical map  . If correlations in the environment decay much faster than the timescale of the evolution in the system of interest, we may neglect memory effects describing how the system has previously interacted with the environment. This is also known as the Markov approximation. For example, consider a thermal state of the environment of the form

 

where   is the partition function,   is the inverse temperature and   is the energy of the state  . Then, if the environment is large and its dynamics is fast enough, any energy exchanged with the system will quickly dissipate away to form a new thermal state with almost exactly the same temperature. Then, from the viewpoint of the system, the state of the environment will appear to be almost constant all the time.

Formally, we can express the consequences of the Markov approximation on the dynamical map as [1]

 

Note that the constraint on the times being positive means that we can only piecewise propagate the system forward in time, i.e, the inverse of the dynamical map does usually not exist. This is in contrast to coherent dynamics, where there is an inverse operation corresponding to negative time arguments in the unitary evolution operator. Hence, while the dynamical maps of coherent systems form a group, the dynamical maps for open quantum systems only form a semigroup.

The generator of the semigroup is the Liouvillian  , which is a generalization of the superoperator appearing on the right hand side of the Liouville von Neumann equation. One important consequence of this generalization is that the von Neumann entropy is no longer a conserved quantity. However, the Liouvillian has to fulfill the property of being the generator of a completely positive and trace-preserving dynamical map. In the following, we will see how the most general form of the Markovian master equation will look like.

Most general form of the dynamics

edit

Similar to the case of a closed quantum system, we can write the dynamical map of an open quantum system as an exponential of the generator of the semigroup,

 

The superoperator   reduces to the one of the Liouville von Neumann equation in the case of purely coherent dynamics, but in general will have additional incoherent terms. Expanding the dynamical maps for short times  , we obtain

 

which in the limit   yields a first-order differential equation known as a quantum master equation,

 

Let us now derive an explicit form for the master equation. For this, we need to define an operator basis  . The inner product for operators is defined as

 

A complete orthonormal set consists of   operators, where   is the Hilbert space dimension. It is convenient to choose one of the operators as proportional to the identity, i.e.,   [1]. Then, all other operators are traceless. For example, in a two-level system, the remaining operators are proportional to the Pauli matrices.

We can now express the action of the dynamical map using this operator basis as

 

where the coefficients   is given by

 

with the operators   defined as above. The coefficients   form a postive matrix  , as for any  -dimensional vector  , we have [1]

 

Inserting this expansion into the quantum master equation, we obtain

 

where we separated off all terms containing  . As a next step, we define the following quantities

 

Note that the operator   is Hermitian, although   is not because the coefficients   are complex. Using these definitions, we find for the generator

 

Since the dynamical map is trace-preserving, the trace over quantum master equation has to vanish, i.e.,

 

from which we can read off that   has to be

 

Substituting this result back into the quantum master equation, we obtain

 

Finally, the matrix formed by the coefficients   is again Hermitian and positive, so we can diagonalize it to obtain positive eigenvalues  . Then, we find the most general form of a Markovian quantum master equation to be given by

 

where the operators   are appropriate linear combinations of the operators   obtained from the diagonalization procedure. This form of the quantum master equation is known as the Lindblad form, as Lindblad first showed that the generator of a Markovian master equation has to be of that form [2].

The system Hamiltonian is contained in the Hermitian operator  , but the latter can also include additional terms coming from the interaction with the environment. Furthermore, the eigenvalues   correspond to relaxation rates describing incoherent decay processes in the system. Typically, these decay processes will result in the system eventually reaching a stationary state characterized by  . However, such a stationary state does not necessarily mean the absence of any dynamics: for example, a single realization of two-level system in a maximally mixed state might still violently jump between both levels! Only when the ensemble average is taken, the dynamics will vanish.

References

edit
  1. 1.0 1.1 1.2 1.3 1.4 Breuer, Heinz-Peter; Petruccione, F (2002). The theory of open quantum systems. Oxford; New York: Oxford University Press. ISBN 9780198520634. 
  2. Lindblad, G. (1976-06-01). "On the generators of quantum dynamical semigroups". Communications in Mathematical Physics 48 (2): 119–130. doi:10.1007/BF01608499. ISSN 1432-0916 0010-3616, 1432-0916. http://link.springer.com/article/10.1007/BF01608499. Retrieved 2014-05-15.