Solutions of the optical Bloch equations

From amowiki
Jump to navigation Jump to search

The optical Bloch equations

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\rho} = -\frac{i}{\hbar}[H,\rho] -\frac{\Gamma}{2} \left[ { \sigma_+ \sigma_- \rho - 2 \sigma_- \rho \sigma_+ + \rho \sigma_+ \sigma_- } \right] \,, }

provide a time-dependent quantum description of a spontaneously emitting atom driven by a classical electromagnetic field. Considerable insight into the physical processes involved can be gained by studying these equations in the transient excitation limit, as well as the steady-state limit, as we see in this section. We begin by considering the coherent part of the evolution, then extend this to re-visit the Bloch sphere picture of the optical Bloch equations, which provides useful visualizations of transient responses and steady state solutions. Finally, we return to vacuum Rabi oscillations and investigate how cavity loss leads to damping of the oscillations, as an illustration of master equations which are more general than the optical Bloch equation.

Eigenstates of the Jaynes-Cummings Hamiltonian

The Hamiltonian involved in the optical Bloch equations,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle H = \frac{\hbar\omega_0}{2}(|e \rangle \langle e|-|g \rangle \langle g|) + \frac{\hbar\Omega}{2} (|g \rangle \langle e|+|e \rangle \langle g|) \,, }

describes the evolution of an atom in a classical field. We begin here by reviewing the coherent evolution under .

arises from the Jaynes-Cummings interaction we have previously considered in the context of cavity QED, describing a single two-level atom interacting with a single mode of the electromagnetic field:

The last term in this expression is , the dipole interaction between atom and field. By defining and , we may write this interaction as

In the frame of reference of the atom and field, recall that

When near resonance, , and because the and terms oscillate at nearly twice the frequency of , those terms can be dropped. Doing so is known as the rotating wave approximation, and it gives us a simplified interaction Hamiltonian

Under this approximation, it is useful to note that this interaction merely exchanges one quantum of excitation from atom to field, and back, so that the total number of excitations is a constant of the motion. We may thus write the total Hamiltonian, in the rotating wave approximation, as

where we have defined Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma_z = |e \rangle \langle e|-|g \rangle \langle g|} , and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta = \omega_0-\omega} . Below, we may use Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g = \hbar\Omega/2} to simplify writing.

What are the eigenstates of this Hamiltonian? It describes a two-level system coupled to a simple harmonic oscillator; when uncoupled, if Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta=0} , then the eigenstates are simply those of , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |g{\rangle}} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |n{\rangle}} , as shown here:

Transient and steady state solutions of the optical Bloch equations-obe-Hjc.png

When coupled, degenerate energy levels split (this is sometimes called dynamimc Stark splitting), with harmonic oscillator levels Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |n{\rangle}} and splitting into two energy levels separated by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 2g\sqrt{n+1}} . Since the coupling only pairs levels separated by one quantum of excitation, it is straightforward to show that the eigenstates of the Jaynes-Cummings Hamiltonian fall into well defined pairs of states, which we may label as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\pm,n{\rangle}} ; these are

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\pm,n \rangle = \frac{1}{\sqrt{2}} \left[ { |e,n \rangle \pm |g,n+1 \rangle } \right] \,, }

and they have energies

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\pm,n} = \hbar\omega(n+1) \pm g\sqrt{n+1} \,. }

When Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta\neq 0} , similar physics result, but with slightly more complicated expressions describing the eigenstates, as we shall see when we later return to the "dressed states" picture.

Strongly driven atom: Mollow triplet

An atom strongly coupled to a single mode electromagnetic field, or an atom driven strongly by a single mode field, will thus have an emission spectrum described by the coupled energy level diagram:

Transient and steady state solutions of the optical Bloch equations-obe-mollow-energy.png

where, to good approximation, the energy level differences are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega\pm\Omega} . These three lines which appear in the spectrum are known as the Mollow triplet:

Transient and steady state solutions of the optical Bloch equations-obe-mollow.png

The Mollow triplet is experimentally observed in a wide variety of systems. However, while our energy eigenstate analysis has predicted the number and frequencies of the emission lines, it fails to explain a key characteristic: the widths are not the same. If the central peak at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_0} has width HWHM Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma/2} , the two sidebands each have a HWHM of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 3\Gamma/4} . To explain this, we need the optical Bloch equations.

Optical Bloch equation evolution on the Bloch sphere

We have previously seen that an arbitrary qubit state can be represented as being a point on a unit sphere, located at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle (\theta,\phi)} in polar coordinates. Similarly, a density matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho} may be depicted as being a point inside or on the unit sphere, using

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho = \frac{I+\vec{r}\cdot{\sigma}}{2} \,, }

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{r}} is the Bloch vector representation of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho} .

Explicitly, if we let

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho = \left[ \begin{array}{cc}{\rho_{ee}}&{\rho_{eg}}\\{\rho_{ge}}&{\rho_{gg}}\end{array}\right] \,, }

then

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{array}{rcl} r_x &=& \langle X \rangle = \rho_{ge}+\rho_{eg} \\ r_y &=& \langle Y \rangle = (\rho_{ge}-\rho_{eg})/i \\ r_z &=& \langle Z \rangle = \rho_{gg}-\rho_{ee} \,. \end{array}}

Visualization: Bloch vector

Visulization of the evolution of a density matrix under the optical Bloch equations is thus helped by rewriting them in terms of a differential equation for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{r}} . A convenient starting point for this is the optical Bloch equation

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\rho} = -\frac{i}{\hbar}[H,\rho] -\frac{\Gamma}{2} \left[ { \sigma_+ \sigma_- \rho - 2 \sigma_- \rho \sigma_+ + \rho \sigma_+ \sigma_- } \right] \,, }

using the rotating frame Hamiltonian (suppressing Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hbar} )

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle H = \frac{\delta}{2} Z + g X \,. }

This gives us the equations of motion

Note how these equations of motion provide a simple set of flows on the Bloch sphere: the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta} terms correspond to a rotation in the plane, corresponds to a rotation in the plane, and drives a relaxation process which shrinks and components of the Bloch vector, while moving the component toward .


Physics: in-phase and quadrature components

Physically, what is the meaning of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_x} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_y} , and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_z} ? Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_z} is manifestly the population difference between the excited and ground states. The other two components may be interpreted by recognizing that the average dipole moment of the atom is

Thus, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_x} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_y} correspond to the phase components of the atomic dipole moment which are in-phase and in quadrature with the incident electromagnetic field.

Animation of full solutions

Here are plots of the full solutions of the optical Bloch equations on the Bloch sphere, for various cases of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta} , and (sample matlab files to solve the differential equations: obefun1.m plotobe2.m):

<jwplayer width="560" height="440" repeat="true" displayheight="420" image="http://feynman.mit.edu/8.422/plotobe6.png" autostart="false">http://feynman.mit.edu/8.422/plotobe6.flv</jwplayer>

Solutions of the optical Bloch equations

The optical Bloch equations, describing the evolution of a single atom coupled to the vacuum,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\rho} = -\frac{i}{\hbar}[H,\rho] -\frac{\Gamma}{2} \left[ { \sigma_+ \sigma_- \rho - 2 \sigma_- \rho \sigma_+ + \rho \sigma_+ \sigma_- } \right] \,, }

are a set of three coupled differential equations,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \dot{r}_x &= \delta\, r_y - \frac{\Gamma}{2} r_x \\ \dot{r}_y &= -\delta\, r_x - g\, r_z - \frac{\Gamma}{2} r_y \\ \dot{r}_z &= g\, r_y - \Gamma (r_z -1) \,, \end{align}}

written in the rotating frame of the atom's Hamiltonian. These equations can be solved analytically, but a great deal of intuition can be obtained from just studying the solutions in several limits. Here, we consider the steady state, transient, and weak excitation limits.

Transient response of the optical Bloch equations

The optical Bloch equations allow us to study the internal state of the atom as it changes due to the external driving field, and due to spontaneous emission.

Starting from the time-independent form of the equations,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \dot{r}_x &= \delta\, r_y - \frac{\Gamma}{2} r_x \\ \dot{r}_y &= -\delta\, r_x - g\, r_z - \frac{\Gamma}{2} r_y \\ \dot{r}_z &= g\, r_y - \Gamma (r_z -1) \,, \end{align}}

we may note that when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g\rightarrow 0} and at resonance, , the Bloch vector exhibits pure damping behavior, towards Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_z=1} , and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_x=r_y=0} . Note that a convenient way to write these differential equations is in matrix form,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\vec{r}} = \left[ \begin{array}{ccc} {-\Gamma/2}& {0} & {0} \\ {0} & -\Gamma/2 & -g \\ 0 & g & -\Gamma \end{array} \right] \vec{r} + \left[ \begin{array}{c} 0 \\ 0 \\ \Gamma \end{array} \right] }

When Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g\gg \Gamma} , Rabi oscillations occur, represnted by rapid rotations of the Bloch vector about . Since the relaxation along Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat{y}} occurs at rate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma/2} , and the relaxation about Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat{z}} occurs at rate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma} , we might expect that the average relaxation rate of the rotating components under such a strong driving field would be . The remaining component Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_x} does not rotate, because it sits along Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat{x}} , the axis of rotation. Thus, it relaxes with rate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma/2} . Computation of the eigenvalues of the equations of motion verify this qualitative picture, and show that for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g\gg \Gamma} , and , the eigenvalues of motion are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \pm ig + 3\Gamma/4} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma/2} .

These correspond to a main peak at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_0} with width Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma/2} , and two sidebands at , with widths Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 3\Gamma/4} , thus explaining the widths of the observed Mollow triplet lines.

Steady-state solution of the optical Bloch equations

The steady state solution of the optical Bloch equations are found by setting all the time derivatives to zero, giving a set of three simultaneous equations,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} 0 &= \delta\, r_y - \frac{\Gamma}{2} r_x \\ 0 &= -\delta\, r_x - g\, r_z - \frac{\Gamma}{2} r_y \\ 0 &= g\, r_y - \Gamma (r_z -1) \,. \end{align}}

The solutions are (up to overall minus signs which can be absorbed into definitions):

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} r_x &= (g\delta) \frac{1}{\delta^2 + g^2/2 + (\Gamma/2)^2} \\ r_y &= \left(\frac{g\Gamma}{2}\right) \frac{1}{\delta^2 + g^2/2 + (\Gamma/2)^2} \\ r_z &= \left(\delta^2+\frac{\Gamma^2}{4}\right) \frac{1}{\delta^2 + g^2/2 + (\Gamma/2)^2} \end{align}}

Physically, these are Lorentzians; the solution (the component in quadrature with the dipole) corresponds to an absorption curve with half-width

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sqrt{\frac{\Gamma^2}{4}+\frac{g^2}{2}} \,, }

and the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_x} solution (the component in-phase with the dipole) corresponds to a dispersion curve. And under a strong driving field, as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g\rightarrow\infty} , , indicating that the populations in the excited and ground states are equalizing. The steady-state population in the excited state is

an important result that will later be used in studying light forces.


These solutions can be re-expressed in a simplified manner by defining the saturation parameter

in terms of which we find

As , the atomic transitions become saturated, and the linewidth of the transition broadens from its natural value Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma} , becoming Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma' = \Gamma\sqrt{1+s}} on resonance, at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta=0} .

Here are plots of the steady state solution for varying Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g/\Gamma} :

<jwplayer width="560" height="440" repeat="true" displayheight="420" image="http://feynman.mit.edu/8.422/steady4.png" autostart="false">http://feynman.mit.edu/8.422/steady4.flv</jwplayer>

Weak excitation limit solution of the optical Bloch equations

Once again, the time-independent form of the optical Bloch equations is

where the remaining two components of the density matrix are given by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{gg} = 1-\rho_{ee}} , and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{eg} = \rho_{ge}^*} . It is insightful to study these equations in the limit of weak excitation, and for short evolution times.

It turns out that the solution of these equations to lowest order in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\Omega|} , and in the limit , with the initial conditions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{ee}=0} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{ge}=0} , gives

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{ee} = \frac{\frac{1}{4}|\Omega|^2}{(\omega_0-\omega_L)^2+\left(\frac{\Gamma}{2}\right)^2} \left[{ 1+e^{-\Gamma t}-2\cos[(\omega_0-\omega_L) t]e^{-\Gamma t/2}}\right] \,. }

Moreover the solution of these equations to lowest order in in the limit Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\Omega|t\ll 1} , with the initial conditions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{ee}=0} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{ge}=0} , gives

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{ee} = \frac{1}{4}|\Omega|^2 t^2 \,, }

irrespective of the values of and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma} .

Damped vacuum Rabi oscillations

The optical Bloch equations can be generalized to describe not just an atom interacting with the vacuum, but also an atom and a single cavity mode, each interacting with its own reservoir. This is the master equation for cavity QED, and using such a master equation we can revisit the phenomenon of vacuum Rabi oscillations and see what happens in the presence of damping.

Generalization of the optical Bloch equations

The starting point for generalizing the optical Bloch equations is the Lindblad form we previously saw in the full-derivation walkthrough,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\rho}_{A} = -\frac{i}{\hbar} [H,\rho_A] -\frac{\Gamma}{2} \left[ { \sigma_+ \sigma_- \rho_A - 2 \sigma_- \rho_A \sigma_+ + \rho_A \sigma_+ \sigma_- } \right] \,. }

In this expression, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma_+} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma_-} are "jump operators," and represent changes that occur to the atom when distinct dissipative events happen ("distinct" meaning that the environment changes between orthogonal states).

We may write the master equation for our more general scenario by replacing the atomic density matrix by a general density matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho} representing the atom and cavity field, and by replacing the atomic jump operators with general jump operators Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_k} ,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\rho} = -\frac{i}{\hbar} [H,\rho] +\sum_k \left[ { L_k \rho L^\dagger _k - \frac{1}{2}( L^\dagger _k L_k \rho + \rho L^\dagger _k L_k ) } \right] \,. }

Note that the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_k} include normalization factors which reflect their probabilities of occurrence. In other words, for the atom + vacuum model, .

Master equation for cavity QED system

For the cavity QED model, the atom and cavity field each have possible jump operators. In general, the atom and cavity may both couple to a thermal field with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle n_{th}} average photons. In such a case, the jump operators are Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_1 = \sqrt{\Gamma(n_{th}+1)}\sigma_-} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_3 = \sqrt{\Gamma n_{th}}\sigma_+} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_2 = \sqrt{\kappa(n_{th}+1)} a} , and , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma} parameterizes the spontaneous emission rate of the atom in free space, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \kappa} is parameterizes the cavity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q} factor.

Experimentally, typically the environment, the vacuum, is essentially at zero temperature, so Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle n_{th}=0} , in which case the only two relevant jump operators are and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_2 = \sqrt{\kappa}a} .

Damped vacuum Rabi oscillations

Vacuum Rabi oscillations, in the absence of damping, involve only two states of the atom and cavity: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |e,0{\rangle}} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |g,1{\rangle}} . When damping is added, the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |g,0{\rangle}} state must be included, since both the atom and cavity states can decay and loose their quanta of energy. Moreover, because only one quantum of excitation is involved in this system, we can observe the essential physics by considering the case when is zero (no spontaneous emission), but Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \kappa} is nonzero (the cavity is leaky). Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Omega_0} be the vacuum Rabi frequency, and denote this three state space by , , and . Written out explicitly in terms of the density matrix elements , the master equation is (for the case Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta=0} with the system starting in the state Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |1\rangle} )

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{array}{rcl} \dot{\rho}_{11} &=& i\frac{\Omega_0}{2}(\rho_{12} - \rho_{21}) \\ \dot{\rho}_{22} &=& -i\frac{\Omega_0}{2}(\rho_{12} - \rho_{21}) -\kappa \rho_{22} \\ \dot{\rho}_{12} &=& i\frac{\Omega_0}{2}(\rho_{11} - \rho_{22}) -\frac{\kappa}{2} \rho_{12} \\ \dot{\rho}_{21} &=& -i\frac{\Omega_0}{2}(\rho_{11} - \rho_{22}) -\frac{\kappa}{2} \rho_{21} \\ \dot{\rho}_{33} &=& \kappa \rho_{22} \,. \end{array}}

When the cavity damping rate is small, , then the vacuum Rabi oscillations are damped, with average damping rate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \kappa/2} since the atoms spend roughly half their time in the excited state.

When the cavity damping rate is large, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \kappa\gg\Omega_0} , then the atomic excitation is irreversibly damped, and no oscillations occur. Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle P_e} be the probability of being in the state. We can combine two of the above equations to write

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{array}{rcl} \dot{\rho}_{12} - \dot{\rho}_{21} &=& i\Omega_0(\rho_{11} - \rho_{22}) -\frac{\kappa}{2} (\rho_{12}-\rho_{21}) \,. \end{array}}

Since in this case Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\rho}_{12} + \dot{\rho}_{21}\approx 0} , and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{22}\approx 0} , it follows that

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle (\rho_{12}-\rho_{21}) \approx \frac{2i\Omega_0}{\kappa}(\rho_{11}-\rho_{22}) \,, }

Since in this regime, we slowly try to populate state while it quickly decays away, we can approximate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{22}\ll\rho_{11}} giving

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle (\rho_{12}-\rho_{21}) \approx \frac{2i\Omega_0}{\kappa}\rho_{11} \,, }

Finally plugging this back into the optical bloch equation for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{11}} yields

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\rho}_{11} = -\frac{\Omega_0^2}{\kappa}\rho_{11} \,, }

so decays exponentially, with rate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma_c = \Omega_0^2/\kappa} . Note that a large cavity damping rate slows down the decay of the atom's excited state. This is due to the fact that the stimulated transition (Rabi oscillation) to the lower state is slowed down by the cavity damping (which is a manifestation of the quantum Zeno effect).

Purcell factor: cavity enhanced spontaneous emission

How does this compare with the free space spontaneous emission rate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma} ? Recall that the vacuum Rabi frequency is (for the case that the dipole matrix element vector and cavity mode polarization are aligned)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Omega_0 = \frac{d E_0}{\hbar} \,, }

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle d} is the atomic dipole moment and is the electric field amplitude of a single photon at the atomic transition frequency Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_0} (Note that this is a standing wave in a cavity so it is a factor of two larger than the field for a running wave),

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_0 = \sqrt{\frac{2\hbar\omega_0}{\epsilon_0 V}} \,, }

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V} is the cavity volume. This gives

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Omega_0 = d \sqrt{\frac{2\omega_0}{\epsilon_0 \hbar V}} \,. }

Letting , we thus find that

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma_c = \frac{2d^2}{\epsilon_0\hbar} \frac{Q}{V} \, }

as the decay rate of the atom in the cavity.

Recall that the spontaneous emission rate of an atom in free space, as determined by Fermi's golden rule, is

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma = \frac{d^2 \omega_0^3}{3\pi\epsilon_0 \hbar c^3} }

The ratio of this rate to the decay rate in the cavity is

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta = \frac{\Gamma_c}{\Gamma} = \frac{3}{4\pi^2} \frac{Q \lambda^3}{V} \,, }

where we take Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \lambda = 2\pi c/\omega_0} as being the wavelength of the cavity field, which is assumed to be resonant with the atomic transition frequency . Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta} is independent of the atomic dipole strength, and determined solely by cavity parameters. Moreover, note that for small, high-Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q} cavities, the decay rate of the atom in the cavity can be much larger than the free space spontaneous emission rate. This "cavity enhanced" spontaneous emission rate was predicted by Purcell (1946), an observation credited as being the starting point of cavity QED.

References