Open Quantum Systems/The Quantum Optical Master Equation

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

Weak-coupling approximation

edit

In the following, we will consider a system   that is weakly coupled to a bath  . The Hamiltonian of the combined system is given by

 

where   denotes the part of Hamiltonian only acting on  ,   only acts on  , and   accounts for the interaction between the two. We first perform a unitary transformation   into the interaction picture, i.e. (assuming  ),

 

Inserting the transformed density matrix into the Liouville von Neumann equation, we obtain

 

which can be brought into the convenient form

 

The interaction Hamiltonian is now time-dependent according to

 

We can formally integrate the Liouville--Von Neumann equation and obtain

 

This expression can be inserted back into the Liouville von Neumann equation and after taking the trace over the bath, we arrive at [1]

 

where we have assumed that the initial state is such that the interaction does not generate any (first-order) dynamics in the bath, i.e.,

 

So far, we have made little progress, as the right-hand side of the equation of motion still contains the density operator of the full system. To obtain a closed equation of motion for   only, we assume that the interaction is weak such that the influence on the bath is small. This is also known as the Born approximation. Then, we may treat the bath as approximately constant and write for the total density operator

 

Then, we obtain a closed integro-differential equation for the density operator  ,

 

Such an integro-differential equation is very difficult to handle as the dynamics at time   depends on the state of the system in all previous times. The equation of motion can be brought into a time-local form by replacing   by  . As we will see later on, this is not really an approximation, but already implicit in the weak-coupling assumption. This step brings us to an equation known as the Redfield master equation [1],

 

This equation is still a non-Markovian master equation and does not guarantee to conserve positivity of the density matrix due to approximations we have made.

Markov approximation

edit

To obtain a Markovian master equation, we first substitute   by   in the integrand, which does not change the bounds of the integration. Then, we can understand the parameter   as indicating how far we go backwards in time to account for memory effects, which can be characteristic timescale  , over which correlations in the bath decay. Under the Markov approximation, these memory effects are short-lived and therefore the integrand decays very quickly for  . Then, we may replace the upper bound of the integration by infinity, and obtain a Markovian master equation,

 

The typical timescale of the dynamics generated by this master equation is characterized by  , which is the relaxation time of the system due to the interaction with the bath. For the Markov approximation to be valid, this relaxation time has to be long compared to the bath correlation time, i.e,  . For quantum optical systems,   is the inverse of the optical frequency, i.e., several inverse THz, while the lifetime of an optical excitation is in the inverse MHz range. Therefore, the Markov approximation is well justified in quantum optical systems.

The two approximation we have made so far are often grouped together as the Born-Markov approximation. However, they still do not guarantee that the resulting master equation generates a quantum dynamical semigroup and hence cannot be cast into a Lindblad form. For this, another approximation is necessary.

Secular approximation

edit

The secular approximation involves discarding fast oscillating terms in the Markovian master equation. It is therefore similar to the rotating-wave approximation (RWA) used in NMR and quantum optics. However, applying the RWA directly on the level on the interaction Hamiltonian can cause problems such as an incorrect renormalization of the system Hamiltonian [2]. The secular approximation (which is sometimes also called RWA), however, is carried out on the level of the quantum master equation.

To be explicit, let us write the interaction Hamiltonian in the form

 

where the Hermitian operators   and   only act on system and bath, respectively. Assume, we have already diagonalized the system Hamiltonian  , so we know its eigenvalues   and its projectors onto eigenstates  . Then, we can project the operators   on subspaces with a fixed energy difference   [1],

 

As the eigenvectors of   form a complete set, we can recover   by summing over all frequencies,

 

Then, we can write the interaction Hamiltonian as

 

Using the relation

 

and its Hermitian conjugate, we can write the interaction Hamiltonian in the interaction picture as

 

where we have introduced the interaction picture operators of the bath,

 

Inserting this expression into the Markovian master equation leads us to [1]

 

where we have used the one-sided Fourier transform of the bath correlation functions,

 

In the case where the state of the bath   is an eigenstate of the bath Hamiltonian  , the bath correlations do not depend on time.

The secular approximation is then referred to as the omission of all terms with  , as these terms oscillate fast and average out. Again, in quantum optical systems, this comes from the fact that optical transition frequencies are much larger than the decay rates of excited states, i.e.,  . Hence, we obtain

 

This master equation can now be cast into a Lindblad form. For this, we split real and imaginary parts of the coefficients   according to

 

where the real part can be written as

 

and forms a positive matrix [1]. Then, diagonalizing the coefficient matrix yields the Lindblad form for the dynamics

 

where the Lamb shift Hamiltonian   commutes with the system Hamiltonian and hence results in a renormalization of the energy levels. This master equation is still expressed in the interaction picture, it can be transformed back to the Schrödinger picture by adding the system Hamiltonian   to the coherent part of the dynamics.

Interaction with the radiation field

edit

In the case of an atom interacting with the quantized radiation field, the Hamiltonian for the latter is given by a sum of harmonic oscillators,

 

where   is the polarization index [3]. If the wavelength of the radiation field is much larger than the spatial extent of the atomic wavefunction, the dominant interaction term is given by the electric dipole term,

 

where   is the dipole operator of the atom and   is the quantized electric field,

 

We now assume the reference state of the bath,   is the vacuum without any photons. Then, we can use the following relations for the reservoir correlations [1]

 

Using this, we find for the spectral correlation tensor

 

Going to the continuum limit, we can use

 

Fortunately, this means that the spectral correlation tensor is already in a diagonal form and we can directly obtain a Lindblad form for the master equation. For the  -integration, we make use of the relation

 

where   denotes the Cauchy principal value. We can thus split the spectral correlation tensor into its real and imaginary part to obtain the decay rate and the Lamb shift, respectively. Note. however, that the expression for the Lamb shift is divergent as we have neglected relativistic effects that become relevant for large values of  . Since it only gives a slight renormalization of the energy levels of the atom, we neglect it in the following. For simplicity, let us consider the case where we are interested in the transition between only two atomic levels differing by the frequency  . Then, we finally arrive at the master equation in Lindblad form,

 

where we have used the spontaneous emission rate

 

and the spin flip operators

 

The master equation can be solved by expanding it into Pauli matrices. The vector   is known as the Bloch vector and the equations of motion for its component decouple,

 

From this, we see that the off-diagonal elements of the density matrix decay exponentially with the rate  , a phenomenon which is called decoherence. The  -component of the Bloch vector decays exponentially with the rate   to a steady state of  , i.e., the atom ends up in the state with lower energy.

Resonance fluorescence

edit

Let us assume now that in addition to the radiation field being in the vacuum state, there is a single mode that is being driven by an external laser. As the laser only affects a single mode, we can ignore its effect on the dissipative dynamics, i.e., the only consequence results from the Hamiltonian of the form

 

where we have assumed that the laser is on resonance with the atom. Denoting the two levels of the atom by  , and  , respectively, we can go into the rotating frame of the laser driving by making the unitary transformation

 

Inserting back into the master equation, we see that the rotation exactly compensates the energy difference between   and  , and the laser term becomes

 

In the rotating wave approximation we neglect the fast oscillating term. Then, we can write the laser Hamiltonian as

 

where we have introduced the Rabi frequency  . The total quantum master equation then reads

 

The interesting aspect about this master equation is that it exhibits a competition between the coherent dynamics generated by the laser and the dissipative dynamics arising from the decay into the vacuum of the radiation field. As before, we can write the master equation as an equation of motion for the Bloch vector [1],

 

using the matrix   given by

 

and the vector  ,

 

This equation of motion is called the optical Bloch equation. Its stationary state   can be found from the condition   and is given by

 

Note that the population of the excited state,

 

is always less than   even in the limit of strong driving, i.e.,  . Thus, it is not possible to create population inversion in a two level system in the stationary state by coherent driving. The population will merely saturate at  .

It is also interesting to look at the relaxation dynamics towards the stationary state. For this, we introduce another vector expressing the difference between the Bloch vector and the stationary solution [1], i.e.,

 

The dynamics of this vector is described by a homogeneous differential equation,

 

We find the eigenvalues of   to be

 

where   is given by

 

Since all eigenvalues have negative real parts, all coefficients of the   vector will eventually decay and the stationary state   is reached. If the atom is initially in the state  , the resulting dynamics is given by

 

References

edit
  1. 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 Breuer, Heinz-Peter; Petruccione, F (2002). The theory of open quantum systems. Oxford; New York: Oxford University Press. ISBN 9780198520634. 
  2. Fleming, Chris; N. I. Cummings, Charis Anastopoulos, B. L. Hu (2010-10-08). "The rotating-wave approximation: consistency and applicability from an open quantum system analysis". Journal of Physics A: Mathematical and Theoretical 43 (40): 405304. doi:10.1088/1751-8113/43/40/405304. ISSN 1751-8121. http://iopscience.iop.org/1751-8121/43/40/405304. Retrieved 2014-05-22. 
  3. Walls, D. F.; Milburn, Gerard J. (1994). Quantum Optics. Springer-Verlag. ISBN 9780387588315.