Unraveling quantum open system dynamics

Open quantum system dynamics are traditionally studied with the master equation, a differential equation for the density matrix describing the system. An alternative, and equivalent approach, utilizing pure-state wavefunctions and stochastic evolution, is presented in this section. This approach, known variously as the "quantum monte-carlo wavefunction" technique, or the method of quantum jumps or quantum trajectories, is based on the fact that density matrices can be represented as stochastic combinations of pure states. We apply this technique to modeling an atom coupled to a vacuum, and show how the monte-carlo approach generates average evolution equivalent to the optical Bloch equations. The approach offers new ways to physically interpret open quantum system evolution, but must be used with caution, since there are an infinite number of equivalent interpretations, many of which may appear to conflict with each other.

▼ Category:8.422

► Category:Quantum Light

► Category:Optical Bloch Equations

► Category:Bose-Einstein condensates

► Category:Ion traps and quantum information

► Category:Light forces

► Category:Photon-atom interactions

  Introduction to Atomic Physics

  Quantum light: states and dynamics

  Optical Bloch equations

  Photon-atom interactions

  Frontiers in Atomic Physics

  Ion traps and quantum information

  Quantum gases

Quantum jumps

Much of traditional atomic physics has focused on ensembles of many atoms; correspondingly, the density matrix formalism we have employed in studying an atom coupled to the vacuum provides an ensemble description of quantum dynamics. In fact, the concept of ensembles is vital to the foundations of quantum mechanics; only ensemble averages can be measured, and one might even say that in a sense, only ensemble averages are "real." Schr\"dinger, in 1952, wrote (and we quote here from Gerry and Knight, "Introductory Quantum Optics," Cambridge University Press, 2005):

... we never experiment with just one electron or atom or (small) molecule. In thought experiments, we sometimes assume that we do; this invariably entails ridiculous consequences\ldots In the first place, it is fair to state that we are not experimenting with single particles, any more than we can raise ichthysauria in the zoo.

Despite this rather extreme opinion, modern atomic physics experiments are now broadly moving in the direction of coherently addressing and manipulating single atoms and single molecules in the laboratory. Experiments have been successful with both (hot and cold) neutral atoms and ions, and with a variety of single (hot) molecules. These date back to Dehmelt's 1975 proposal to trap single ions to provide an atomic time standard.

Such experiments with single quantum systems motivate a description of open quantum system dynamics in terms of statistical averages over single realizations, as illustrated by the following experiment. Consider a single three-level ion, such as strontium, with a $8$ ns lifetime 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 422} nm S-P transition, and a 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 \sim 400} ms lifetime 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 674} nm S-D transition:

Experimental observation of the blue fluorescence from this atom, when strongly excited 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 422} nm, and weakly excited 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 674} nm, shows a very interesting signature: quantum jumps. When the atom is in the metastable D state, 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 422} nm fluorescence ceases, because no photons can be scattered from the strong transition until it transitions back to the ground S state. A typical observation record might look like this:

Naturally, the average fluorescence rate will be consistent with that expected by the relative strengths of the red and blue excitations, but this would be just a constant value. Such quantum jumps, however, are ubiquitous in experiments involving single atoms and other single quanta, and it would be nice to have a theoretical prescription for describing such observations. In addition, it would be nice to have available possible interpretations of what how the individual systems might be evolving to give some observed ensemble behavior, such as an exponential decay.

Quantum Monte-Carlo wavefunction technique for atom + vacuum

Consider a single two-level atom evolving under the Hamiltonian 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_0} , while also undergoing spontaneous emission. This is the same scenario we studied to formulate the optical Bloch equations. It can be modeled using the quantum Monte-Carlo wavefunction (QMCWF) technique, using the following steps.

Let us assume the atom begins in some 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 |\psi{\rangle}} , and let the atomic levels be $|e{\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{\rangle}} . 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 \Gamma} be the decay rate of the excited state, as determined, for example, from Fermi's golden rule. The QMCWF technique introduces a new non-Hermitian "Hamiltonian"

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 = H_0 - i\frac{\Gamma}{2} |e \rangle \langle e| \,, }

which governs evolution of the state, according to the following. For each differential time step 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 dt} , 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 |\psi{\rangle}} evolves as:

1. Compute 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 dp = \Gamma dt | \langle e| \psi{\rangle}|^2} 2. 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 0\leq \epsilon \leq 1} be a uniformly distributed random number 3. 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 \epsilon<dp} 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 |\psi \rangle \leftarrow |g{\rangle}} 4. 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 \epsilon\ge dp} 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 |\psi \rangle \leftarrow \left[ { e^{-i H dt}/\sqrt{1-dp}} \right] |\psi{\rangle}} 5. Go to 1

Physically, 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 dp} is the probability of the atom jumping from 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{\rangle}} to 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}} , that is, spontaneously emitting a photon, during the time interval 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 dt} . Note that typically 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 dp\ll 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 dp = 1-{\langle}\psi| e^{iH^\dagger dt}e^{-iH dt}|\psi{\rangle} \,, }

so $1-dp$ is the normalization of the state after 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 dt} of evolution under the non-Hermitian Hamiltonian $H$. Essentially, the QMCWF scheme models an imaginary observer watching for a spontaneously emitted photon coming from the atom. If a photon is emitted, then the atom transitions immediately into 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}} ; this is the rule given in step 3, and such a transition is known as a quantum jump. Otherwise, if no photon is emitted, the state changes nevertheless, albeit by only a small amount. Because the observer saw no photon, and thus the probability of the atom being in $|e{\rangle }$ is slightly diminished. Specifically, 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 H_0=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 |\psi{\rangle}=\alpha |g{\rangle}+\beta|e{\rangle}} , then according to the rule in step 4,

$|\psi \rangle \leftarrow {\frac {\alpha |g\rangle +\beta e^{-\Gamma dt/2}|e{\rangle }}{\sqrt {|\alpha |^{2}+|\beta |^{2}e^{-\Gamma dt}}}}\,.$

This is exactly the same physics we have seen earlier, in the beamspliter model of the optical Bloch equations.

Equivalence of QMCWF and the optical Bloch equations

The equivalence of QMCWF and the optical Bloch equations (OBE) is shown by demonstrating that the evolution of the 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(t) = \overline{|\psi(t){\rangle}{\langle}\psi(t)|} }

satisfies the OBE, where the average (denoted by the overline) is taken over instances of running the QMCWF procedure. This follows from computing the density matrix for the state after one QMCWF procedure step:

${\begin{array}{rcl}\rho (t+dt)&=&dp\,|g\rangle \langle g|+(1-dp){\frac {e^{-iHdt}|\psi {\rangle }{\langle }\psi |e^{+iH^{\dagger }dt}}{1-dp}}\\&\approx &dp\,|g\rangle \langle g|+(1-iHdt)\rho (t)(1+iH^{\dagger }dt)\\&\approx &dp\,|g\rangle \langle g|+\rho (t)-i(H\rho (t)-\rho (t)H^{\dagger })dt\\&=&\Gamma dt\langle e|\rho (t)|e\rangle \,|g\rangle \langle g|+\rho (t)-i(H\rho (t)-\rho (t)H^{\dagger })dt\\&=&\Gamma dt|g\rangle \langle e|\rho (t)|e\rangle \langle g|+\rho (t)-i\left[H_{0},\rho (t)\right]-{\frac {\Gamma }{2}}\left(|e\rangle \langle e|\rho (t)+\rho (t)|e\rangle \langle e|\right)\,.\end{array}}$

Writing this as a coarse-grained differential equation, taking the limit of small 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 dt} , we find

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} \frac{d}{dt} \rho(t) &\approx& \frac{\rho(t+dt)-\rho(t)}{dt} \\ &=& - i \left [H_0,\rho(t) \right] -\frac{\Gamma}{2} \left( |e \rangle \langle e| \rho(t) + \rho(t)|e \rangle \langle e| \right) + \Gamma dt |g \rangle \langle e| \rho(t) |e \rangle \langle g| \,. \end{array}}

This is the optical Bloch equation.

Quantum Monte-Carlo wavefunction technique: general case

The method of the quantum Monte-Carlo wavefunction technique can be applied not just to the atom + vacuum scenario, but also to model any open quantum system dynamics. The relaxation part of a master equation may be described by a Lindblad operator

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 { \mathcal L}(\rho) = -\frac{1}{2} \sum_k C^\dagger _k C_k \rho + \rho C^\dagger _k C_k + \sum_k C_k \rho C^\dagger _k \,, }

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 C_k} are known as "quantum jump" operators. The corresponding QMCWF procedure employs the non-Hermitian Hamiltonian

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 = H_0 - \frac{i\hbar}{2} \sum_k C^\dagger _k C_k \,, }

and the steps

1. Compute 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 dp = \sum_k dp_k} , 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 dp_k = dt {\langle}\psi| C^\dagger _k C_k|\psi{\rangle}} 2. 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 0\leq \epsilon \leq 1} be a uniformly distributed random number 3. 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 \epsilon<dp} 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 |\psi \rangle \leftarrow C_k|\psi \rangle / \sqrt{dp_k/dt}} 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 k} randomly chosen with probability 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 dp_k/dp} . 4. 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 \epsilon\ge dp} 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 |\psi \rangle \leftarrow \left[ { e^{-i H dt}/\sqrt{1-dp}} \right] |\psi{\rangle}} 5. Go to 1

There are some computational advantages to the QMCWF approach, over numerical solution of the differential equations normally obtained with master equations. In particular, an 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} dimensional Hilbert space described by a density matrix requires 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 O(N^2)} variables, whereas a pure-state wavefunction requires only 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 O(N)} variables. Of course, calculating expectation values means that the stochastic evolution of the QMCWF technique must be repeated 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 \sim N} times, so the overall effect is a tradeoff between storage space and computational time. However, the QMCWF technique is immediately parallelizable, and often (but not always) desired observables converge quickly.

For more about this, see, for example, the nice article <refbase>5239</refbase>

Example: spontaneous emission

The dynamics of an atom interacting with the vacuum can be desribed not just by the usual optical Bloch equations, but also, by an infinite variety of alternate pictures. This can be seen thorough a simple example using the quantum monte-carlo wavefunction technique.

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 X} , $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 Z} denote the Pauli matrices, as usual. Recall that 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] + \left[ { L \rho L^\dagger - \frac{1}{2}( L^\dagger L \rho + \rho L^\dagger L ) } \right] }

can be written 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 \dot{\vec{r}} = \left[ \begin{array}{ccc} {-\Gamma/2}& {0} & {0} \\ {0} & -\Gamma/2 & 0 \\ 0 & 0 & -\Gamma \end{array} \right] \vec{r} - \left[ \begin{array}{c} 0 \\ 0 \\ \Gamma \end{array} \right] \,, }

using

$\rho ={\frac {I+{\vec {r}}\cdot {\vec {\sigma }}}{2}}\,.$

In the "standard" unraveling of spontaneous emission into quantum jumps, in the absence of a classical driving field, there is no Hamiltonian evolution (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=I} ) and one jump operator, 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=\sqrt{\Gamma}\sigma_- = \sqrt{\Gamma} ( X-i Y)/2} . This gives exactly the standard optical Bloch equations, as was seen above. Unraveling the dynamics with the quantum monte carlo wavefunction technique, we obtain evolution trajectories which look like this:

<jwplayer width="700" height="820" repeat="true" displayheight="800" image="http://feynman.mit.edu/8.422/simse5a.png" autostart="false">http://feynman.mit.edu/8.422/simse5a.flv</jwplayer>

The average over many such evolutions leads to the exponential decay expected from the optical Bloch equations:

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

An alternate unraveling of spontaneous emission

However, it turns out that identical dynamics are obtained, on average, with the same form of master equation, but 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 H=-\Gamma Y/4} and $L={\sqrt {\Gamma }}(I+X-iY)/2$, instead.

<jwplayer width="700" height="820" repeat="true" displayheight="800" image="http://feynman.mit.edu/8.422/simse5.png" autostart="false">http://feynman.mit.edu/8.422/simse5.flv</jwplayer>

The average over many such evolutions also leads to the exponential decay expected from the optical Bloch equations:

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

Analytically, the exact agreement can be proven by writing the master equation as a set of differential equations 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}} , with the alternate choice 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 H} 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} .

The origin of this alternate unraveling of spontaneous emission arises in the lack of knowledge about what basis the environment is observing the system in, as can be seen in the discussion about another representation of open quantum system dynamics, the operator sum representation.

Example: phase and amplitude damping

Spontaneous emission, which is mainly what we have studied so far, with the OBE and QMCWF, is described by the density matrix evolution

$\rho (t=0)=\left[{\begin{array}{cc}{1-c}&{b^{*}}\\{b}&{c}\end{array}}\right]\longrightarrow \rho (t)=\left[{\begin{array}{cc}{1-ce^{-\Gamma t}}&{b^{*}e^{-\Gamma t/2}}\\{be^{-\Gamma t/2}}&{ce^{-\Gamma t}}\end{array}}\right]\,.$

Here, 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} is the spontaneous emission rate of the atom, parameterizing the loss of energy from the atom (the decay of the diagonal elements 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} to 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 \langle g|} ), and the simultaneous loss of phase coherence of the atom (the decay of the off-diagonal elements 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} to zero). The quantum jump operator for spontaneous emission, in a two-level atom, 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 C = \sqrt{\Gamma} \sigma^-} . More generally, however, we can have two parameters which describe separate decay of the diagonal and off-diagonal elements,

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(t) = \left[ \begin{array}{cc}{1-c e^{-t/T_1}}&{b^* e^{-t/T_2}}\\ {b e^{-t/T_2}}&{c e^{-t/T_1}}\end{array}\right] \,, }

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 T_1} is the energy loss rate, 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 T_2} is the rate of loss of quantum coherence, also known as the phase damping rate. Both of these effects are known as decoherence.

In particular, we define decoherence as any process which can turn pure states into statistical mixtures. 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} is pure if either 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 {\rm Tr}(\rho^2)=1} , or 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 = |\psi{\rangle}{\langle}\psi|} for some 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 |\psi{\rangle}} , or the entropy 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 S(\rho) = -{\rm Tr}(\rho \log \rho)=0} . Otherwise, it is mixed.

Phase damping, which gives rise to 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 T_2} , is in a sense the most "quantum" kind of noise; it comes along with spontaneous emission, but can also exist by itself. Historically, the term "decoherence" has sometimes been identified exclusively with phase damping, because of its important role in the emergence of classical behavior from quantum systems, but today, decoherence is used as a much more general term, because there are so many ways to loose coherence from a quantum system other than just phase damping.

Open system dynamics have infinitely many equivalent unravelings

The various physical origins of phase damping are interesting to consider, because they can teach us an important fact about models of decoherence for single quantum systems: just 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 \rho} can be unraveled in infinitely many ways as $\sum _{k}p_{k}|\psi _{k}{\rangle }{\langle }\psi _{k}|$, decoherence processes may also be unraveled in an infinite number of ways, each equally equivalent and equally physically meaningful.

We consider three different physical models of phase damping, on a two-level atom.

Random phase noise

A two-level atom of 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} excited by far off-resonance light, 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\gg\omega_0} , experiences an AC stark shift of amount proportional to the light intensity. If this intensity fluctuates, then the atom's phase is randomly modulated, causing the evolution

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 \rangle \rightarrow e^{i\theta}|e{\rangle} \,, }

(6.2.1)

where $\theta$ is the phase imparted by the AC stark shift. Suppose 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} is modeled as a Gaussian distributed random variable, with mean zero and variance 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 2\lambda t} , such 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 {\rm prob}(\theta) = \frac{1}{\sqrt{4\pi\lambda t}} e^{-\theta^2/4\lambda t} \,. }

If the initial state of the atom is $\psi \rangle =a|g\rangle +b|e{\rangle }$, then after time 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 t} , it evolves into the average state described by the 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(t) = \int_{-\infty}^{+\infty} \left[ \begin{array}{cc}{|a|^2}&{ab^* e^{i\theta}}\\{a^*b e^{-i\theta}}&{|b|^2}\end{array}\right] {\rm prob}(\theta) \,d\theta \,, }

which is found to be

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(t) = \left[ \begin{array}{cc} {|a|^2}& {ab^* e^{-\lambda t}} \\ {a^*b e^{-\lambda t}}& {|b|^2} \end{array}\right] \,. }

If the atomic Hamiltonian is $H_{0}=\Omega |e\rangle \langle e|$, and the atom's initial state 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 |\psi \rangle = (|g{\rangle}+|e \rangle )/\sqrt{2}} , then the dipole moment of the atom shows a simple Rabi oscillation:

When random phase noise is imposed on the atom, then its dipole moment decays with time. The evolution of a single atom, according to the "trajectory" described by the random walk of Eq.(6.2.1), is a noisy Rabi oscillation:

Elastic collisions

Another physical origin for phase damping is elastic collisions. Assume the two-level atom bounces along a waveguide, interacting with the walls without loosing kinetic energy, but changing its trajectory slightly at each bounce, in a manner depending on the state of the atom. This can be modeled by a Hamiltonian interacting the atom with a single mode environment,

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_{SE} = |e \rangle \langle e| \otimes \left[ { \gamma |0 \rangle \langle 1| + \gamma^* |1 \rangle \langle 0| } \right] \,, }

with coupling constant 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} . During a small differential timestep $dt$, an initial atomic state $a|g\rangle +b|e{\rangle }$ coupled to an environment $|0{\rangle }$ evolves to become

${\begin{array}{rcl}(a|g\rangle +b|e\rangle )\otimes |0{\rangle }&\rightarrow &a|g{\rangle }|0\rangle +b|e\rangle (\cos \theta |0\rangle +\sin \theta |1\rangle )\\&=&\left[{a|g\rangle +b\cos \theta |e\rangle }\right]|0{\rangle }+\left[{b\sin \theta |e\rangle }\right]|1{\rangle }\,,\end{array}}$

where $e^{-\lambda dt}=\cos \theta$. This expression is very similar to that obtained for the gedankenexperiment used in the QMCWF model of spontaneous emission; the difference is that when a photon is observed in the environment, the atom does not collapse into $|g{\rangle }$, but rather, into 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{\rangle}} . In other words, it does not loose energy; it only looses information about what state it was in, prior to the collapse.

Just as in the proof of the equivalence of QMCWF to the OBE, we can compute the density matrix evolution which this model gives rise to, by writing down an expression 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(t)} , based on Eq.(\ref{eq:ecmodel}),

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} \rho(t) &=& \left[ { a|g \rangle + b\cos\theta |e \rangle } \right] \left[ { \langle g|a^* + \langle e| b^*\cos\theta } \right] + |b|^2 \sin^2 \theta |e \rangle \langle e| \\ &=& \left[ \begin{array}{cc}{|a|^2}&{ab^* e^{-\lambda t}}\\{a^*b e^{-\lambda t}}&{|b|^2}\end{array}\right] \,. \end{array}}

Note that this is exactly the same evolution as we obtained for the random phase noise model, Eq.(\ref{eq:pdamp}).

Despite the density matrix evolution being identical to that of the random phase model, the elastic collision model implies a different single-particle evolution trajectory. In contrast to the noisy Rabi oscillations previously seen, for the elastic collisions, the atomic state initially decays, then jumps into the $|e{\rangle }$ state at some random time; this is illustrated by this sample trajectory:

Phase flips

A third physical model for phase damping is based on quantum jumps. The Lindblad operator for phase damping is evidently

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 { \mathcal L}(\rho) = -\frac{1}{T_2} \left[ { |e \rangle \langle e| \rho |g \rangle \langle g| + |g \rangle \langle g| \rho |e \rangle \langle e| } \right] \,. }

Equivalently, it may be rewritten in standard form 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 { \mathcal L}(\rho) = -\frac{1}{2} \left[ { C^\dagger C \rho + \rho C^\dagger C } \right] + C \rho C^\dagger \,, }

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 C = \frac{1}{\sqrt{2 T_2}} \left[ { |g \rangle \langle g| - |e \rangle \langle e| } \right] \,. }

Note that $C^{\dagger }C$ is proportional to the identity, and 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 H=H_0} , such that no evolution occurs due to this relaxation process except when a quantum jump occurs. Moreover, the effect of a quantum jump is to flip the phase of the atom 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 -1} , changing 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{\rangle} \rightarrow -|e{\rangle}} . If such a flip happens with probability $(1-e^{-\lambda t})/2$ at each moment in time, then the density matrix for this evolution 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 \begin{array}{rcl} \rho(t) &=& \frac{1+e^{-\lambda t}}{2} \left[ { a|g \rangle + b|e \rangle } \right] \left[ { \langle g|a^* + \langle e| b^* } \right] + \frac{1-e^{-\lambda t}}{2} \left[ { a|g \rangle - b|e \rangle } \right] \left[ { \langle g|a^* - \langle e| b^* } \right] \\ &=& \left[ \begin{array}{cc}{|a|^2}&{ab^* e^{-\lambda t}}\\{a^*b e^{-\lambda t}}&{|b|^2}\end{array}\right] \,. \end{array}}

This is again the same density matrix dynamics as previously obtained for the random phase noise and elastic collision models. However, the trajectories of individual evolutions is different; at each moment in time, the two-level atom either remains completely unchanged, or its excited state flips sign, inverting its dipole moment:

Discussion

We have seen three models of phase damping, all of which produce the same density matrix evolution, but each of which has very different microscopic trajectories for individually evolving systems. Which is correct?

The answer is that all of them are correct, and yet none are. Any of the three can be used for physical intuition and interpretation, but only as long as the only conclusions drawn depend on statistical averages. In fact, in the absence of control over the environment, no experiment can distinguish between phase damping processes described by these three models, even in principle. This strong statement arises from the fact that there are an infinite number of ways a (mixed) density matrix can be written as a statistical mixture of pure states; correspondingly, there are an infinite number of "unravelings" of density matrix time evolutions, into statistical evolutions of pure state wavefunctions.

The freedom of interpretation which arises in studying decoherence processes arises from a unitary degree of freedom. For example, in the QMCWF model, a gedankenexperiment is performed, in which the atom is allowed to decay and the emitted photon is captured. The evolution of the atom cannot depend on what measurement basis is used for the photodetection. This basis choice is a unitary transform which can be chosen arbitrarily in the gedankenexperiment, and different choices lead to the different unravelings of the master equation into trajectories.

References

• <refbase>5239</refbase>
• <refbase>5249</refbase>
• <refbase>5250</refbase>
• <refbase>5251</refbase>