Resonances

Classical Resonances

Introduction

▼ Category:8.421

► Category:Atoms in electric fields

  Atoms

  Atoms in electric fields

  Atoms in Magnetic Fields

  Coherence

  Effects of the Nucleus on Atomic Structure

  Fine Structure and Lamb Shift

  Higher-Order Radiation Processes

  Interaction of an atom with an electromagnetic field

  Old Coherence

  Resonances

  Two-Photon Excitation

Classically, a resonance is a system where one or more variables change periodically such that when the system is no longer driven by an extended periodic process, the decay of the system's oscillation takes many (or at least several) oscillation periods. Equivalently, the response of the driven system as a function of drive frequency exhibits some form of peaked structure. The frequency corresponding to the maximum response of the system is called the resonance frequency ($f_{0}$)

<figure id="fig:example_resonances">

A resonance in the time (left) and frequency (right) domains. The time domain picture shows the (decaying) oscillation in some variable of the undriven system, while the frequency-domain pictures shows the steady-state response amplitude of the driven system. The time units are chosen such that $f_{0}=\omega _{0}/2\pi =1$ while the quality factor is $Q=10$. </figure>

The simplest resonance systems have a single resonance frequency, corresponding to the harmonic motion of the variable. Some form of dissipation (damping) results in a finite damping time ($\Delta t<\infty$) for the undriven system, and a finite width ($\Delta f>0$) of the resonance curve for the driven system. The quantities $\Delta {f}$ and $\Delta {t}$ are related by a Fourier transformation: any oscillation with time-varying amplitude must be made up of a superposition of different frequency components, giving the resonance curve a finite width in frequency space. Low dissipation results in a long decay time $\Delta t$, and a narrow frequency width $\Delta f$. The ratio of frequency width $\Delta f$ to resonance frequency $f_{0}$ is called the quality factor $Q$.

$Q={\frac {f_{0}}{\Delta f}}$

Atomic physics often deals with isolated atomic systems in a vacuum, providing resonances of high quality factor $Q$. For example, an optical transition, even in a room-temperature (ie Doppler-broadened) gas, will have a $Q$-factor of

Failed to parse (unknown function "\hertz"): {\displaystyle \frac{f_0}{\Delta f} \sim \frac{{10^{15}}{\hertz}}{{{10^9}\,{{\rm Hz}}}} \sim 10^6 }

In contrast, high $Q$ factors in solids typically require cryogenic temperatures: typical $Q$ factors of mechanical or electrical systems are $\sim 10^{3}$ at room temperature and $\sim 10^{6}$ at $\lesssim {{1}\,{\rm {K}}}$ temperatures.

High Q values are obtained in solids acting as optical resonators: $Q\gtrsim 10^{8}$ has been achieved in the whispering gallery modes of spheres or cylinders of high-purity glass .

One of the whispering-gallery resonators described in . This particular resonator had a quality factor $Q=10^{8}$. The very high surface quality required for such low losses is achieved by having surface tension shape the rim of the disk after it has been temporarily liquefied by a laser pulse; the inset shows the resonator before this treatment.

At the macroscopic scale, astronomical systems consisting of massive bodies travelling essentially undisturbed through the near-vacuum of space can also exhibit high $Q$ factors. The rotation of the earth has a $Q$ value of $10^{7}$, and pulsars of $10^{10}$.

A resonance is particularly useful for science and technology if it is reproducible by different realizations of the same system, and if it has a well-understood theoretical connection to fundamental constants. In atomic physics, quantum mechanics ensures both features: atoms of the same species are fundamentally identical (although the probing apparatus is not), and the quantum mechanical description of atomic structure is well understood and highly accurate, at least for atoms with not too many electrons. The resonance frequencies of hydrogen, and to a lesser degree helium, are directly tied via quantum mechanics to fundamental constants such as the Rydberg constant, $R_{\infty }={{1.0973731568525(73)\times 10^{7}}\,{{\rm {m}}^{-1}}}$ , or the fine structure constant $\alpha =1/137.03599911(46)$ . In fact, the Rydberg constant is the most accurately known constant in all of physics, and most fundamental constants have been determined by atomic physics or optical techniques. Additional motivation for measuring fundamental constants with ever-increasing accuracy has been provided by the speculation (and claims of evidence for it ), that the fundamental constants may not be so constant after all, but that their value is tied to the evolution (and therefore the age) of the universe. Furthermore unexpected resonances, or splittings of resonances, have in the past indicated the need for new theories or modifications of theories. For instance, "anomalous" effects in the Zeeman structure (magnetic-field dependance) of atomic lines led to the discovery of spin . while Lamb's measurement of the splitting between the $^{2}S_{1/2}$ and $^{2}P_{1/2}$ states of atomic hydrogen, on the order of 1000MHz instead of 0 as predicted by the Dirac theory, fostered the development of quantum electrodynamics, the quantum theory of light.

Classical Harmonic Oscillator

A classical harmonic oscillator (HO) is a system described by the second order differential equation

${\ddot {q}}+\gamma {\dot {q}}+\omega _{0}^{2}q=0$

such as a mass on a spring, or a charge, voltage, or current in an electric RLC resonant circuit. For weak damping ($\gamma ^{2}<4\omega _{0}^{2}$) the undriven system decays with an exponential envelope,

$e^{-{\frac {\gamma }{2}}t}e^{\pm {i\omega ^{\prime }t}}$

with $\omega ^{\prime }=\omega _{0}{\sqrt {1-\gamma ^{2}/4\omega _{0}^{2}}}$. In the cases we will consider, $\gamma ^{2}\ll {4\omega _{0}^{2}}$ so that $\omega ^{\prime }\approx \omega _{0}$. The decay rate constant for the amplitude is $\gamma /2$, while for the stored energy it is $\gamma$, i.e. the stored energy decays exponentially with a time constant $\tau =1/\gamma$. If the harmonic oscillator is driven at a frequency $\omega$ close to resonance, $\left|\omega _{0}-\omega \right|\ll \omega _{0}$, the amplitude response $q_{0}$ is a lorentzian:

$q_{0}={\frac {\text{const.}}{\omega _{0}-\omega +i{\frac {\gamma }{2}}}}={\frac {\text{const.}}{-\Delta +i{\frac {\gamma }{2}}}}$

where $\Delta =\omega -\omega _{0}$ is the detuning as shown on the right-hand side of <xr id="fig:example_resonances"/>. The full width at half maximum (FWHM) of a lorentzian curve is $\Delta \omega =\gamma$ leading to a quality factor

$Q={\frac {\omega _{0}}{\Delta \omega }}={\frac {\omega _{0}}{\gamma }}$

Note that $\omega _{0}$, $\Delta w$, and $\gamma$ are measured in angular frequency units \unit{\radian\per\second}, or \unit{\reciprocal\second} for short, and should not be confused with frequencies $f_{0}=\omega _{0}/2\pi$, $\Delta f=\Delta \omega /2\pi$ etc. When quoting values, we will often write explicitly $\omega _{0}=2\pi \times {{1}\,{\rm {MHz}}}$ rather than $\omega _{0}={{6.28\times 10^{6}}\,{{\rm {s}}^{-1}}}$ to remind ourselves that the quantity in question is an angular frequency. We will never 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 \omega_0=\unit{6.28\times10^6}{Hz}} : although formally dimensionally consistent, this invites confusion with real (laboratory or Hertzian) frequencies. For brevity, we will often write or say "frequency" when we mean angular frequency, but will quote real frequencies as $f_{0}={{1}\,{\rm {MHz}}}$. Note that the natural unit for decay rate constants is that of an angular frequency, as in $e^{-i\omega _{0}t-\gamma t}$. Therefore we will express decay rates as inverse times in angular frequency units: $\gamma ={{10^{4}}\,{{\rm {s}}^{-1}}}$ for instance, but not $\gamma ={{10^{4}}\,{\rm {Hz}}}$, nor $\gamma =2\pi \times {{1.6}\,{\rm {KHz}}}$. Damping time constants are the inverse of decay rate constants, $\tau =1/\gamma ={{100}\,{\mu {\rm {s}}}}$.

Resonance Widths and Uncertainty Relations

From $\Delta \omega =\gamma$ it follows that damping time $\tau$ and the resonance linewidth $\Delta \omega$ of the driven system obey: $\Delta \omega \cdot \tau =1$ or, assuming $E=\hbar \omega$, $\Delta E\tau =\hbar$: a finite number of frequency components or energies ($\Delta \omega >0$) is necessary to synthesize a pulse of finite ($\tau <\infty$) width in time.

<figure id="fig:gated_pulse">

Measurement in a limited time $\Delta t$, modelled as an ideal monochromatic sine wave of frequency $\omega _{0}/2\pi =1$ gated by an observation window of finite width $\Delta t$ </figure>

This is quite similar to an uncertainty relation familiar from Fourier transformation theory, according to which the measurement of a frequency in a finite time $\Delta t$ (see <xr id="fig:gated_pulse"/>) results in a frequency spread $\Delta \omega$ obeying $\Delta \omega \Delta t\geq 1/2$. By using $E=\hbar \omega$ we can further connect this with the usual Heisenberg uncertainty relation $\Delta E\Delta t\geq \hbar /2$

Question: <definition id="q:classical_line_splitting" noautocaption noblock>(Q%i)</definition> Does the Heisenberg uncertainty relation $\Delta E\Delta t\geq \hbar /2$, or rather the frequency uncertainty relation $\Delta \omega \Delta t\geq 1/2$ that follows from (classical) Fourier theory hold in classical systems? Can you measure the angular frequency of a classical harmonic oscillator in a time $\Delta t$ with an accuracy better than $\Delta \omega =1/2\Delta t$?

Answer: Yes - if you have a good model for the signal, for instance if you know that it is a pure sine wave, and if the signal to noise ratio (SNR) of the measurement is good enough, then you can fit the sine wave's frequency to much better accuracy than its spectral width $\Delta \omega$ (see Fig. <xr id="fig:line_splitting"/>). As a general rule, the line center can be found to within an uncertainty $\delta \omega \approx \Delta \omega /{\text{SNR}}$. This is known as splitting the line. Splitting the line by a factor of 100 is fairly straightforward, splitting it by a factor of a 1000 is a real challenge. As an extreme example of this, Caesium fountain clocks interrogate a 9 GHz transition for a typical measurement time of 1 s, yielding $\Delta \omega /\omega \approx 10^{-10}$, but the best ones can nonetheless achieve relative accuracies of $10^{-16}$ . Of course, this sort of performance requires a very good understanding of the line shape and of systematics.

<figure id="fig:line_splitting">

Splitting the line in time and frequency domains. Given a sufficiently good model of the signal or line shape, and adequate signal-to-noise, one can find the center frequency $\omega _{0}$ of the underlying signal to much better than $1{\text{ cycle}}/\Delta t$ or the line width. </figure>

Question: <definition id="q:quantum_line_splitting" noautocaption noblock>(Q%i)</definition> Can you measure the angular frequency of a quantum mechanical harmonic oscillator in a time $\Delta t$ to better than $\Delta \omega =1/2\Delta t$?

Question: <definition id="q:laser_line_splitting" noautocaption noblock>(Q%i)</definition> Can you measure the angular frequency of a laser pulse lasting a time $\Delta t$ to better than $\Delta \omega =1/2\Delta t$?

If you answered <xr id="q:quantum_line_splitting"/> and <xr id="q:classical_line_splitting"/> differently, you may have a problem with <xr id="q:laser_line_splitting"/>. Is a laser pulse quantum or classical? Clearly it is possible to beat $\Delta \omega =1/2\Delta t$ for a laser pulse. For instance, you can superimpose the laser pulse with a laser of fixed and known frequency on a photodiode, and fit the observed beat note. The stronger the probe pulse, the better the SNR of the beat note, and the more accurately you can extract the probe frequency. So how are we to reconcile quantum mechanics and classical mechanics, and what does the Heisenberg uncertainty really state?

The Heisenberg uncertainty relation makes a statement about predicting the outcome of a single measurement on a single system. For a single photon, the limit $1/2\Delta t$ does indeed hold. However, nothing in quantum mechanics says that you cannot improve your knowledge about the average values of the quantities characterizing the system by repeated measurements on identical copies of it. In the laser example above, the frequency resolution improves with probe pulse strength because we are measuring many photons - for the purposes of this discussion we could have measured the frequencies of the photons sequentially.

For independent measurement of $n$ uncorrelated photons the SNR is given by ${\sqrt {n}}$: the signal grows as $n$, while the noise is the shot noise of the photon number, which grows as ${\sqrt {n}}$. This SNR allows a frequency uncertainty of $1/2\Delta t{\sqrt {n}}$, known as the standard quantum limit. If we allow the use of entangled or correlated states of the photons, then the uncertainty can be as low as the Heisenberg limit $1/2\Delta tn$ .

Let us now revisit <xr id="q:quantum_line_splitting"/>. You may have heard that the quantum mechanical description of electromagnetism, quantum electrodynamics (QED), represents any single mode of the electromagnetic field by a harmonic oscillator of the same frequency. The ground state $n=0$ corresponds to no photons in the mode, the first excited state $n=1$ to one photon and so forth. This description is valid because photons to a very good approximation do not interact: the addition of a photon changes the system's energy always by the same increment. Since the many-photon pulse from <xr id="q:laser_line_splitting"/> maps onto a quantum harmonic oscillator, it seems to follow that <xr id="q:quantum_line_splitting"/> must also be answered in the affirmative.

On the other hand, a harmonic oscillator is a single quantum system to which the Heisenberg limit must apply. The resolution of this apparent contradiction is that <xr id="q:quantum_line_splitting"/> is asking about a measurement of the fundamental frequency of the harmonic oscillator, not the energy of a particular level. Therefore if you measure the energy of the n-th level to $\Delta E=\hbar /2\Delta t$, which can be done by exciting the $n$-th level from the ground state through a nonlinear process, then the uncertainty on the harmonic oscillator resonance frequency (that corresponds to the laser frequency in <xr id="q:quantum_line_splitting"/>) is given by $\Delta \omega _{0}=\Delta E/n=1/n2\Delta t$.

If the interaction with the oscillator is restricted to classical fields or forces, then there is an additional uncertainty of ${\sqrt {n}}$ on which particular level of the oscillator is being measured and so we recover the standard quantum limit $1/2\Delta t{\sqrt {n}}$ mentioned previously.

Magnetic Resonance: The Classical Magnetic Moment in a Spatially Uniform Field

Harmonic oscillator versus two-level systems

Now, unlike in classical mechanics, most resonances studied in atomic physics are not harmonic oscillators, but two-level systems. Unlike harmonic oscillators, two-level systems show saturation. When a harmonic oscillator is driven longer or faster, higher and higher excited states are populated: the oscillator amplitude can become arbitrarily large. In contrast, the amplitude of oscillation in a two-level system is limited to one half in the appropriate dimensionless units.(An analogous dimensionless amplitude for the harmonic oscillator would be the amplitude measured in units of the oscillator ground state size.)

Why, and under what conditions, do classical harmonic oscillator models of a two-level system work? A two-level system of energy difference $E$ can be approximated by a harmonic oscillator of frequency $\omega _{o}=E/\hbar$ when saturation effects in the two-level system are negligible, i.e. when the population of the second excited state in the harmonic oscillator problem is negligible, or equivalently when the population ratio of the first excited state to the the ground state is small: $P_{1}/P_{0}\ll 1$. This is the basis for the classical Lorentz model of the electron bound in the atom that describes many linear atomic properties (for instance the refractive index) very well. See "The origin of the refractive index" in chapter 31 of the Feynman Lectures on Physics .

When saturation comes into play, i.e. when the initial ground state is appreciably depleted, the harmonic oscillator ceases to be a good model. The limit on the oscillation amplitude in a two-level system suggests that a classical system with periodic evolution and a limit on the amplitude, namely rotation, could provide a better classical model of a two-level system. Indeed, the motion of a classical magnetic moment in a uniform field provides a model that captures almost all features of the quantum-mechanical two-level system, the exception being the projection onto one of the two possible outcomes in a measurement.

Magnetic Moment in a Static Field

The interaction energy of a magnetic moment ${\bf {\mu }}$ with a magnetic field ${\bf {B}}$ is given by

$W=-{\bf {\mu }}\cdot {\bf {B}}$

In a uniform field the force

${\bf {F}}=-{\bf {\nabla }}W=0$

vanishes, but the torque

${\bf {T}}={\bf {\mu }}\times {\bf {B}}$

does not. For an angular momentum ${\bf {L}}$ the equation of motion is then <equation id="eq:classical_precession_in_static_field" noautocaption> (%i) ${\dot {\bf {L}}}={\bf {T}}={\bf {\mu }}\times {\bf {B}}=\gamma {\bf {L}}\times {\bf {B}}$ </equation> where we have introduced the gyromagnetic ratio $\gamma$ as the proportionality constant between angular momentum and magnetic moment, as shown in Figure <xr id="fig:static_precession"/>.

<figure id="fig:static_precession">

Precession of a the magnetic moment and associated angular momentum about a static field ${\bf {B}}$.

</figure>

This describes the precession of the magnetic moment about the magnetic field with angular frequency

$\Omega _{L}=-\gamma B$

$\Omega _{L}$ is known as the Larmor frequency. For electrons we have Failed to parse (unknown function "\unit"): {\displaystyle \gamma_e=2\pi\times\unit{2.8}{\mega\hertz\per G}} , for protons Failed to parse (unknown function "\unit"): {\displaystyle \gamma_p=2\pi\times\unit{4.2}{\kilo\hertz\per G}} .

Transformation to a Rotating Coordinate System

A vector ${\bf {A}}$ rotating at constant angular velocity ${\bf {\Omega }}$ is described by

${\dot {\bf {A}}}={\bf {\Omega }}\times {\bf {A}}$

Then the rates of change of ${\bf {A}}$ measured in a coordinate system rotating at ${\bf {\Omega }}$ and in an inertial system are related by

${\dot {\bf {A}}}_{\text{in}}={\dot {\bf {A}}}_{\text{rot}}+{\bf {\Omega }}\times {\bf {A}}_{\text{in}}$

This follows immediately from the following facts:

• If ${\bf {A}}$ is constant in the rotating system then ${\dot {\bf {A}}}_{\text{in}}={\bf {\Omega }}\times {\bf {A}}_{\text{in}}$.
• If ${\bf {\Omega }}=0$ then ${\dot {\bf {A}}}_{\text{in}}={\dot {\bf {A}}}_{\text{rot}}$.
• Coordinate rotation is a linear transformation.

This transformation is sometimes abbreviated as the schematic rule <equation id="eq:rot_frame_transformation" noautocaption> (%i) $\left({\frac {d}{dt}}\right)_{\text{rot}}=\left({\frac {d}{dt}}\right)_{\text{in}}-{\bf {\Omega }}\times {\big (}{\big )}_{\text{in}}$ </equation> It follows that the angular momentum in a rotating frame obeys

${\dot {\bf {L}}}_{\text{rot}}={\dot {\bf {L}}}_{\text{in}}-{\bf {\Omega }}\times {\bf {L}}_{\text{in}}={\bf {L}}\times (\gamma {\bf {B}}+{\bf {\Omega }})$

If we choose ${\bf {\Omega }}={\bf {\Omega _{L}}}=-\gamma {\bf {B}}$, then ${\bf {L}}_{\text{rot}}$ is constant in the rotating frame. Often it is useful to think of a fictitious magnetic field ${\bf {B}}_{\text{fict}}={\bf {\Omega }}/\gamma$ that appears in a rotating frame.

Larmor's Theorem for a Charged Particle in a Magnetic Field

The vanishing of the torque on a magnetic moment when viewed in the correct rotating frame is reminiscent of Larmor's theorem for the motion of a charged particle in a magnetic field, which we now present.

According to answers.com, Lamor's theorem can be stated as For a system of charged particles, all having the same ratio of charge to mass, moving in a central field of force, the motion in a uniform magnetic induction B is, to first order in B, the same as a possible motion in the absence of B except for the superposition of a common precession of angular frequency equal to the Larmor frequency.

If the Lorentz force acts in an inertial frame,

${\bf {F}}_{\text{in}}=q{\bf {v}}_{\text{in}}\times {\bf {B}}$

then in the rotating frame, according to the rule <xr id="eq:rot_frame_transformation"/> we have

${\begin{aligned}{\dot {\bf {r}}}_{\text{rot}}&={\dot {\bf {r}}}_{\text{in}}-{\bf {\Omega }}\times {\bf {r}}_{\text{in}}\\{\ddot {\bf {r}}}_{\text{rot}}&={\ddot {\bf {r}}}_{\text{in}}-2{\bf {\Omega }}\times {\dot {\bf {r}}}_{\text{in}}+{\bf {\Omega }}\times \left({\bf {\Omega }}\times {\bf {r}}_{\text{in}}\right)\end{aligned}}$

resulting in a force ${\bf {F}}_{\text{rot}}=m{\ddot {\bf {r}}}_{\text{rot}}$ in the rotating frame given by

${\begin{aligned}{\bf {F}}_{\text{rot}}&=q{\bf {v}}_{\text{in}}\times {\bf {B}}-2m{\bf {\Omega }}\times {\bf {v}}_{\text{in}}+m{\bf {\Omega }}\times \left({\bf {\Omega }}\times {\bf {r}}_{\text{in}}\right)\\&=q{\bf {v}}_{\text{in}}\times \left({\bf {B}}+2{\frac {m}{q}}{\bf {\Omega }}\right)+m{\bf {\Omega }}\times \left({\bf {\Omega }}\times {\bf {r}}_{\text{in}}\right)\end{aligned}}$

where we have used $m{\ddot {r}}_{\text{in}}={\bf {F}}_{\text{in}}=q{\bf {v}}_{\text{in}}\times {\bf {B}}$. Choosing

$2{\frac {m}{q}}{\bf {\Omega }}=-{\bf {B}}=-B{\bf {\hat {e}}}_{z}$

yields

${\bf {F}}_{\text{rot}}={\frac {q^{2}}{4m}}B^{2}{\bf {\hat {e}}}_{z}\times \left({\bf {\hat {e}}}_{z}\times {\bf {r}}_{\text{in}}\right)\approx 0$

if the $B$ field is not too large. Thus the Lorentz force approximately disappears in the rotating frame.

Although Larmor's theorem is suggestive of the rotating co-ordinate transformation, it is important to realize that the two transformations, though identical in form, apply to fundamentally different systems. A magnetic moment is not necessarily charged- for example a neutral atom can have a net magnetic moment, and the neutron possesses a magnetic moment in spite of being neutral - and it experiences no net force in a uniform magnetic field. Furthermore, the rotating co-ordinate transformation is exact for a magnetic moment, whereas Larmor's theorem for the motion of a charged particle is only valid when the $B^{2}$ term is neglected.

Note: A breakdown of Larmor's theorem occurs for free electrons. Free electrons orbit in a magnetic field at the cyclotron frequency $\omega _{c}=qB/m$ which is twice Larmor's frequency. The reason is that the orbiting velocity is now given by $v=\omega _{c}r$. As can be seen in the equation above, in this case, the centrifugal force equals minus one half of the Coriolis force, and therefore twice the frequency of the rotating frame is needed to cancel the effects of the applied magnetic field.

Rotating Magnetic Field on Resonance

<figure id="fig:rotating_frame">

Field and moment vectors in the static and rotating frames for the case of resonant drive. </figure>

Consider a magnetic moment ${\bf {\mu }}$ precessing about a field ${\bf {B}}=B_{0}{\bf {\hat {e}}}_{z}$ with ${\bf {\mu }}=(\mu ,\theta ,\phi =-\omega _{0}t)$ in spherical coordinates, where $\omega _{0}=-\gamma B_{0}$ (as shown in fig:rotating_frame/>). Assume that we now apply an additional field ${\bf {B}}_{1}$, in the $xy$-plane rotating at $\omega _{0}$. To solve the resulting problem it is simplest to go into the rotating frame (fig:rotating_frame/>). Then ${\bf {B}}_{1}$ is stationary, say along ${\bf {\hat {e}}}_{x}$, and there is an additional fictitious field ${\bf {B}}_{\text{fict}}={\bf {\omega _{0}}}/\gamma =-{\bf {B}}_{0}$ which cancels the field ${\bf {B}}_{0}$. So in the rotating frame we are left just with the static field ${\bf {B}}_{1}$, and the magnetic moment precesses about ${\bf {B}}_{1}$ at the Rabi frequency

$\omega _{R}=\gamma B_{1}$

A magnetic moment initially along the $-{\bf {\hat {e}}}_{z}$ axis will point along the $+{\bf {\hat {e}}}_{z}$ axis at a time $T$ given by $\omega _{r}T=\pi$, while a magnetic moment parallel or antiparallel to the applied magnetic field ${\bf {B}}_{1}$ does not precess in the rotating frame.

Question: <definition id="q:rabi_freq_blue_detuned" noautocaption noblock>(Q%i)</definition> Assume the magnetic moment is initially pointing along the $-{\bf {\hat {e}}}_{z}$ axis. Assume that a rotating field $B_{1}\ll B_{0}$ is applied, but that it rotates at a frequency $\omega _{1}>\omega _{0}$, where $\omega _{0}$ is the Larmor frequency for the static field ${\bf {B}}=B_{0}{\bf {\hat {e}}}_{z}$. Compared to the on-resonant case, $\omega _{1}=\omega _{0}$, is the oscillation frequency of the magnetic moment.

• larger
• the same
• smaller

Question: Same question as <xr id="q:rabi_freq_blue_detuned"/> but for $\omega _{1}<\omega _{0}$.

Rotating Magnetic Field Off-Resonance

If the rotation frequency $\omega$ of ${\bf {B}}_{1}$ does not equal the Larmor frequency $\omega _{0}=\gamma B_{0}$ associated with the static field ${\bf {B}}_{0}$, then in the frame rotating with ${\bf {B}}_{1}$ at frequency $\omega$ the static field is not completely cancelled by the fictitious field 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/\gamma} , but a difference 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 {\bf{\hat e}} _z} remains, giving rise to a total effective field in the rotating frame

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 {\bf{B}} _\text{eff}= B_1 {\bf{\hat e}} _x + \left(B_0-\frac{\omega}{\gamma}\right) {\bf{\hat e}} _z }

The effective field is static, lies at an angle $\theta$ with the 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 \tan\theta=\frac{B_1}{B_0-\frac{\omega}{\gamma}} }

and is of magnitude

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 \left| B_\text{eff}\right| =\sqrt{B_1^2+\left(B_0-\frac{\omega}{\gamma}\right)^2} }

The magnetic moment precesses around it with an effective (sometimes called generalized) Rabi 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_R =\gamma B_\text{eff} =\sqrt{(\omega_0-\omega)^2+\omega_R^2}=\sqrt{\omega_R^2+\delta^2} }

where $\omega _{R}=\gamma B_{1}$ is the Rabi frequency associated 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 B_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 \delta=\omega-\omega_0} is the detuning from resonance with the Larmor 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=\gamma B_0} .

Geometrical Solution for the Classical Magnetic Moment in Static and Rotating Fields

<figure id="fig:rotating_coord_construction">

Geometrical relations for the spin in combined static and rotating magnetic fields, viewed in the frame co-rotating with the drive field ${\bf {B}}_{1}$. At lower right is a view looking straight down 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 {\bf{B}} _\text{eff}} axis.

</figure>

Referring to <xr id="fig:rotating_coord_construction"/>, we have

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} \phi &= \Omega_R t \\ C &= \mu \sin \alpha \\ A^2 &= \mu^2 ( 1 - \cos \alpha )^2 + C^2 \\ &= \mu^2 ( 1 - 2 \cos \alpha + \cos^2 \alpha + \sin^2 \alpha ) \\ &= 2 \mu^2 (1 - \cos \alpha ) \\ \Rightarrow \cos \alpha &= 1 - \frac{A^2}{2 \mu^2} \end{align}}

On the other hand

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} \frac{A}{2} &= \mu \sin \theta \sin \frac{\phi}{2} \\ A &= 2 \mu \sin \theta \sin \frac{\Omega_R t}{2} \end{align}}

so that <equation id="eq:classical_rabi_flopping" noautocaption> (%i) ${\begin{aligned}\cos \alpha &=1-{\frac {4\mu ^{2}\sin ^{2}\theta \sin ^{2}{\frac {\Omega _{R}t}{2}}}{2\mu ^{2}}}\\&=1-2\sin ^{2}\theta \sin ^{2}{\frac {\Omega _{R}t}{2}}\\\mu _{z}(t)&=\mu \cos \alpha \\&=\mu \left(1-2{\frac {\omega _{R}^{2}}{\delta ^{2}+\omega _{R}^{2}}}\sin ^{2}{\frac {\Omega _{R}t}{2}}\right)\\\mu _{z}(t)&=\mu \left(1-2{\frac {\omega _{R}^{2}}{\Omega _{R}^{2}}}\sin ^{2}{\frac {\Omega _{R}t}{2}}\right)\end{aligned}}$ </equation> 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 \omega_0 = \gamma B_0} the Larmor frequency of the static field, 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-\omega_0} the detuning, 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_R=\gamma B_1} the resonant and $\Omega _{R}={\sqrt {\omega _{R}^{2}+\delta ^{2}}}$ the generalized Rabi frequencies. Note that the precession is faster, but the amplitude smaller for an off-resonant field than for the resonant case. The above result is also the correct quantum-mechanical result.

If we define the quantity $P={\frac {\mu _{z}(t)-\mu _{z}(0)}{2\mu _{z}(0)}}$, we obtain

$P={\frac {\omega _{R}^{2}}{\Omega _{R}^{2}}}\sin ^{2}{\frac {\Omega _{R}t}{2}}$

In the quantum mechanical treatment (see below), $P$ is the probability to find the system in the second state, or the spin flip probability.

"Rapid" Adiabiatic Passage (Classical Treatment)

Rapid adiabatic passage is a technique for inverting a spin by (slowly) sweeping the detuning of a drive field through resonance. "Slowly" means slowly compared to the Larmor frequency $\gamma B_{\text{eff}}$ about the effective static field in the rotating frame for all times ("Rapid" means fast compared to relaxation processes which we neglect here).

The physical picture is as follows. Assume the detuning is initially negative ($\delta <0$, 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 \left| \delta\gg\omega_R\right| } ). Since

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 \tan\theta=\frac{B_1}{B_0-\frac{\omega}{\gamma}}=\frac{\omega_R}{\omega_0-\omega}=-\frac{\omega_R}{\delta} }

the effective magnetic field initially points of a small angle relative to 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 {\bf{\hat e}} _z} axis. If the detuning is increased slowly compared to the Larmor frequency, the spin will continue to precess tightly around 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 {\bf{B}} _\text{eff}} , which for $\delta =0$ points along the x axis, and 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 \delta\gg\omega_R} along 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 - {\bf{\hat e}} _z} axis (see <xr id="fig:rapid_adiabatic_passage"/>).

<figure id="fig:rapid_adiabatic_passage">

{Motion of the spin during rapid adiabatic passage, viewed in the frame rotating 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 {\bf{B}} _1} . The spin's rapid precession locks it to the direction of ${\bf {B}}_{\text{eff}}$ and thus it is dragged through an angle 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 \pi} as the frequency is swept through resonance. </figure>

Thus the magnetic moment, starting out 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 {\bf{B}} _0=B_0 {\bf{\hat e}} _z} , ends up pointing 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 - {\bf{B}} _0=B_0 {\bf{\hat e}} _z} . Note that in the rotating frame ${\bf {\mu }}$ remains always (almost) parallel to the effective field 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 {\bf{B}} _\text{eff}} . A similar precess is used in magnetic traps for atoms, but there 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 {\bf{B}} _\text{eff}} is a real, spatially dependent field constant in time. As the atom moves in this field, the fast precession of the magnetic moment about the local field keeps its direction locked to the local field, whose direction varies in the lab frame. Returning to rapid adiabatic passage, since the generalized Rabi frequency is smallest and equal to the resonant Rabi 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_R} at $\delta =0$, the adiabatic requirement is most severe there, i.e. near 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=\pi/2} . Near 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=\pi/2} we have, 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 B_z=B_0-\omega(t)/\gamma} ,

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 \left| \dot\theta\right| =\frac{ \left| \dot{B}_z\right| }{B_1}=\frac{ \left| \dot\omega\right| }{\gamma B_1} =\frac{ \left| \dot\omega\right| }{\omega_R}\overset{!}{\ll}\omega_R }

where the exclamation point in ${\overset {!}{\ll }}$ indicates a requirement which we impose. Consequently, if the evolution is to be adiabatic, we must have 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 \left| \dot{\omega}\right| \ll\omega_R^2} . This means that the change 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} of rotation 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} per Rabi period 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\pi/\omega_R} , $\left|\Delta \omega \right|=\left|{\dot {\omega }}\right|T=\left|{\dot {\omega }}\right|/\omega _{R}2\pi$, must be small compared to the Rabi 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_R} .

Note that the inversion of the spin is independent of whether 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} is swept up or down through the resonance.

Quantized Spin in a Magnetic Field

Equation of Motion for the Expectation Value

For the system we have been considering, the Hamiltonian 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 \hat{H} =- {\bf{\hat{\mu}}} \cdot {\bf{B}} _0 =-\gamma {\bf{\hat{L}}} \cdot {\bf{B}} _0 =-\gamma\hat{L}_z B_0 \text{ for } {\bf{B}} =B_0 {\bf{\hat{e}}} _z }

Recalling the Heisenberg equation of motion for any 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 \hat{O}} is

${\frac {d}{dt}}{\hat {O}}={\frac {i}{\hbar }}\left[{\hat {H}},{\hat {O}}\right]+{\frac {\partial {\hat {O}}}{\partial t}}$

where the last term refers to operators with an explicit time dependence, we have in this instance

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 \frac{d}{dt}\hat{\mu}_k =\gamma\frac{d}{dt}\hat{L}_k =\frac{i\gamma}{\hbar}\left[\hat{H},\hat{L}_k\right] =-\frac{i\gamma^2}{\hbar}B_0\left[\hat{L}_z,\hat{L}_k\right] }

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 \left[\hat{L}_k,\hat{L}_l\right]=i\hbar\epsilon_{klm}\hat{L}_m} with the Levi-Civita symbol 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_{klm}} , we have

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} \frac{d}{dt} \hat{\mu}_x &= - \frac{i \gamma^2}{\hbar} B_0 i \hbar \hat{L}_y = \gamma^2 \hat{L}_y B_0 = \hat{\mu}_y \gamma B_0 \\ \frac{d}{dt} \hat{\mu}_y &= - \mu_x \gamma B_0 \\ \frac{d}{dt} \hat{\mu}_z &= 0 \end{align}}

or in short

${\begin{aligned}{\frac {d}{dt}}{\bf {\hat {\mu }}}&=\gamma {\bf {\hat {\mu }}}\times {\bf {B}}\\{\frac {d}{dt}}{\bf {\hat {L}}}&=\gamma {\bf {\hat {L}}}\times {\bf {B}}\end{aligned}}$

These are just like the classical equations of motion <xr id="eq:classical_precession_in_static_field"/>, but here they describe the precession of the operator for the magnetic moment 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 {\bf{\hat\mu}} } or for the angular momentum 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 {\bf{\hat L}} } about the magnetic field at the (Larmor) angular 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_L=\gamma B_0} .

Note that:

• Just as in the classical model, these operator equations are exact; we have not neglected any higher order terms.

• Since the equations of motion hold for the operator, they must hold for the expectation value

$\langle {\dot {\bf {\hat {L}}}}\rangle =\gamma \langle {\bf {\hat {L}}}\rangle \times {\bf {B}}$

• We have not made use of any special relations for a spin-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/2} system, but just the general commutation relation for angular momentum. Therefore the result, precession about the magnetic field at the Larmor frequency, remains true for any value of angular momentum 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} .

• A spin-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/2} system has two energy levels, and the two-level problem with coupling between two levels can be mapped onto the problem for a spin in a magnetic field, for which we have developed a good classical intuition.

• If coupling between two or more angular momenta or spins within an atom results in an angular momentum 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 F>1/2} , the time evolution of this angular momentum in an external field is governed by the same physics as for the two-level system. This is true as long as the applied magnetic field is not large enough to break the coupling between the angular momenta; a situation known as the Zeeman regime. Note that if the coupled angular momenta have different gyromagnetic ratios, the gyromagnetic ratio for the composite angular momentum is different from those of the constituents.

• For large magnetic field the interaction of the individual constituents with the magnetic field dominates, and they precess separately about the magnetic field. This is the Paschen-Back regime.

• An even more interesting composite angular momentum arises when $N$ two-level atoms are coupled symmetrically to an external field. In this case we have an effective angular momentum $L=N/2$ for the symmetric coupling:

Level structure diagram for $n$ two-level atoms in a basis of symmetric states . The leftmost column corresponds to an effective spin-$n/2$ object. Other columns correspond to manifolds of symmetric states of 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 n} atoms with lower total effective angular momentum.

Again the equation of motion for the composite angular momentum 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 {\bf{\hat L}} } is a precession. This is the problem considered in Dicke's famous paper , in which he shows that this collective precession can give rise to massively enhanced couplings to external fields ("superradiance") due to constructive interference between the individual atoms.

The Two-Level System: Spin-1/2

Let us now specialize to the two-level system and calculate the time evolution of the occupation probabilities for the two levels.

<figure id="fig:two_level_spin_half">

Equivalence of two-level system with spin-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/2} . Note 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 \gamma>0} the spin-up state (spin aligned with field) is the ground state. For an electron, with $\gamma <0$, spin-up is the excited state. Be careful, as both conventions are used in the literature. </figure>

We have 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 \left\langle S_z\right\rangle =\frac{\hbar}{2}\left(P_\uparrow-P_\downarrow\right) =\frac{\hbar}{2}\left(P_g-P_e\right) =\frac{\hbar}{2}\left(1-2P_e\right) }

where in the last equation we have used the fact 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 P_e+P_g=1} . The signs are chosen for a spin 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 \gamma>0} , such as a proton (Figure <xr id="fig:two_level_spin_half"/>). For an electron, or any other spin with $\gamma <0$, the analysis would be the same but for the opposite sign 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 S_z} and the corresponding exchange 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 \uparrow} 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 \downarrow} . If the system is initially in the ground state, $P_{g}(t=0)=1$ (or the spin 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 + {\bf{\hat e}} _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 S_z(0)=+\hbar/2} ), the expectation value obeys the classical equation of motion <xr id="eq:classical_rabi_flopping"/>: <equation id="eq:rabi_transition_probability" noautocaption> (%i) 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} P_e(t) &= \frac{1}{2} - \frac{1}{\hbar} \left\langle S_z\right\rangle = \frac{1}{2} - \frac{ \left\langle S_z(0)\right\rangle }{\hbar} \left(1 - 2 \frac{\omega_R^2}{\Omega_R^2} \sin^2 \frac{\Omega_R t}{2}\right) \\ P_e(t) &= \frac{\omega_R^2}{\Omega_R^2} sin^2 \left(\frac{\Omega_R t}{2}\right) \end{align}} </equation>

Equation <xr id="eq:rabi_transition_probability"/> is the probability to find system in the excited state at time $t$ if it was in ground state at 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=0} . <xr id="fig:rabi_signal"/> shows a real-world example of such an oscillation.

<figure id="fig:rabi_signal">

Rabi oscillation signal taken in the Vuleti\'{c} lab shortly after this topic was covered in lecture in 2008. The amplitude of the oscillations decays with time due to spatial variations in the strength of the drive field (and hence of the Rabi frequency), so that the different atoms drift out of phase with each other. </figure>

Matrix form of Hamiltonian

With the matrix representation

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} \left|e\right\rangle &= \left|S_z = - \frac{1}{2}\right\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix} \\ \left|g\right\rangle &= \left|S_z = + \frac{1}{2}\right\rangle = \begin{pmatrix} 0 \\ 1 \end{pmatrix} \end{align}}

we can write 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} associated with the static field ${\bf {B}}_{0}=B_{0}{\bf {\hat {e}}}_{z}$ 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 \hat{H}_0 = - {\bf{\hat\mu}} \cdot {\bf{B}} _0 = - \gamma \hat{S}_z B_0 = \hbar \omega_0 \frac{\hat{S}_z}{\hbar} = \frac{1}{2} \hbar \omega_0 \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} = \frac{1}{2} \hbar \omega_0 \hat{\sigma}_z }

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 \omega_0=-\gamma B_0} is the Larmor frequency (recall 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 } is negative for an electron, making $\omega _{0}$ positive), 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 \hat{\sigma}_z= \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} }

is a Pauli spin matrix. The eigenstates 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 \left|e\right\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 \left|g\right\rangle } with eigenenergies $\pm \hbar \omega _{0}/2$. A spin initially 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 {\bf{\hat e}} _x} , corresponding 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 \left|\psi(t=0)\right\rangle =\frac{1}{\sqrt{2}}\left( \left|e\right\rangle + \left|g\right\rangle \right) =\frac{1}{\sqrt{2}}\left( \left|-\frac{1}{2}\right\rangle + \left|+\frac{1}{2}\right\rangle \right) }

evolves in time 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 \left|\psi(t)\right\rangle =\frac{1}{\sqrt{2}} \left(e^{-i\omega_0 t / 2} \left|e\right\rangle +e^{+i\omega_0 t / 2} \left|g\right\rangle \right) =\frac{e^{-i\omega_0 t/2}}{\sqrt{2}}\left( \left|e\right\rangle +e^{i\omega_0 t} \left|g\right\rangle \right) }

which describes a precession with angluar frequency $\omega _{0}$. The field 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 {\bf{B}} _1} , rotating 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} in 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 X-Y} plane corresponds to

${\begin{aligned}{\hat {H}}_{1}&=-{\bf {\hat {\mu }}}\cdot {\bf {B}}_{1}(t)=-{\bf {\hat {\mu }}}\cdot {\frac {\omega _{R}}{\gamma }}\left(-{\bf {\hat {e}}}_{x}\cos \omega t-{\bf {\hat {e}}}_{y}\sin \omega t\right)\\&=\omega _{R}\left({\hat {S}}_{x}\cos \omega t+{\hat {S}}_{y}\sin \omega t\right)\\&={\frac {\hbar \omega _{R}}{2}}\left({\hat {\sigma }}_{x}\cos \omega t+{\hat {\sigma }}_{y}\sin \omega t\right)\\&={\frac {\hbar \omega _{R}}{2}}\left({\begin{pmatrix}0&1\\1&0\end{pmatrix}}\cos \omega t+{\begin{pmatrix}0&-i\\i&0\end{pmatrix}}\sin \omega t\right)\\{\hat {H}}_{1}&={\frac {\hbar \omega _{R}}{2}}{\begin{pmatrix}0&e^{-i\omega t}\\e^{i\omega t}&0\end{pmatrix}}\end{aligned}}$

where have used the Pauli spin matrices 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_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 \sigma_y} . The full Hamiltonian is thus given by

<equation id="eq:dressed_atom_hamiltonian" noautocaption>

(%i) 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{H}=\frac{\hbar}{2} \begin{pmatrix} +\omega_0 & \omega_R e^{-i \omega t} \\ \omega_R e^{i \omega t} & - \omega_0 \end{pmatrix} } </equation> This is the famous "dressed atom" Hamiltonian in the so-called "rotating wave approximation". However, note that in the treatment here there are no approximations; it is exact. Its eigenstates and eigenvalues provide a very elegant, very intuitive solution to the two-state problem.

Solution of the Schrodinger Equation for Spin-1/2 in the Rotating Frame

The time-dependence of the above Hamiltonian can be eliminated by applying a unitary transformation to the system. Namely, a rotation about the z-axis accomplishes this. Call this rotation operator ${\bf {R}}$, which can be expressed

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} {\bf{R}} &= e^{i S_z \theta/\hbar} \\ &= \begin{pmatrix} e^{i \omega t/2} & 0 \\ 0 & e^{-i \omega t/2} \end{pmatrix}, \end{align}}

where 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 \theta = \omega t} . Then the transformed 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 \begin{align} H' = {\bf{R}}H\,{\bf{R}}^{\dagger}-i\hbar\,{\bf{R}}\frac{d\bf{R}^{\dagger}}{dt}, \end{align}}

is

${\begin{aligned}H'&={\frac {\hbar }{2}}{\begin{pmatrix}-\delta &\omega _{\rm {R}}\\\omega _{\rm {R}}&\delta \end{pmatrix}},\end{aligned}}$

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 \delta = \omega - \omega_0} . This time-independent Hamiltonian can be diagonalized exactly to obtain the eigenenergies and eigenstates. The general solution has the form,

${\begin{aligned}\psi (t)=a_{g}(t)|g\rangle +a_{e}(t)|e\rangle .\end{aligned}}$

In particular, the case $a_{g}(0)=1$ and $a_{e}(0)=0$ leads to the probability of finding the atom in the excited state to be

${\begin{aligned}P_{e}(t)=|a_{e}(t)|^{2}={\frac {\omega _{\rm {R}}^{2}}{\Omega _{\rm {R}}^{2}}}\sin ^{2}\left({\frac {\Omega _{\rm {R}}t}{2}}\right),\end{aligned}}$

where $\Omega _{\rm {R}}={\sqrt {\delta ^{2}+\omega _{\rm {R}}^{2}}}$ is the generalized Rabi frequency. This form is identical to the one obtained in the classical picture.

Solution of the Schrodinger Equation for Spin-1/2 in the Interaction Representation

The same system can be solved by a more brute-force method by isolating the "bare" time-dependence first and then using substitutions. The interaction representation consists of expanding 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 \left|\psi(t)\right\rangle } in terms of the eigenstates 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 \left|e\right\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 \left|g\right\rangle } of 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} , including their known time dependence $e^{iH_{0}t/\hbar }$ due 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 H_0} . That means we write 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 \left|\psi(t)\right\rangle =a_e(t) \left|e\right\rangle e^{-i\omega_0t/2}+a_g(t) \left|g\right\rangle e^{i\omega_0t/2} }

Substituting this into the Schrodinger 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 i\hbar\frac{d}{dt} \left|\psi(t)\right\rangle =\left(H_0+H_1\right) \left|\psi(t)\right\rangle }

then results in the equations of motion for the coefficients

${\begin{aligned}i{\dot {a}}_{g}&={\frac {\omega _{R}}{2}}e^{i\left(\omega -\omega _{0}\right)t}a_{e}\\i{\dot {a}}_{e}&={\frac {\omega _{R}}{2}}e^{-i\left(\omega -\omega _{0}\right)t}a_{g}\end{aligned}}$

Where we have used the matrix form of the Hamiltonian, <xr id="eq:dressed_atom_hamiltonian"/>. Introducing the detuning 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-\omega_0} , we have

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} i\dot{a}_g&= \frac{\omega_R}{2}e^{i\delta t}a_e \\ i\dot{a}_e&= \frac{\omega_R}{2}e^{-i\delta t}a_g \end{align}}

The explicit time dependence can be eliminated by the substitution

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} b_g&= e^{-i\delta t/2}a_g \\ b_e&= e^{i\delta t/2}a_e \end{align}}

As you will show (or have shown) in the problem set, this leads to solutions for $b_{g,e}$ 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 \begin{align} b_g(t)&= e^{i\Omega_R t/2}B_1+e^{-i\Omega_R t/2}B_2 \\ b_e(t)&= \frac{\omega_R}{\delta-\Omega_R}e^{i\Omega_R t/2}B_1 +\frac{\omega_R}{\delta+\Omega_R}e^{-i\Omega_R t/2}B_2 \end{align}}

with two constants that are determined by the initial conditions. 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 a_g(t=0)=1} 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 \left| a_e(t)\right| ^2=\frac{\omega_R^2}{\Omega_R^2}\sin^2\frac{\Omega_R t}{2} }

as already derived from the fact that the expectation value for the magnetic moment obeys the classical equation.

"Rapid" Adiabiatic Passage (Quantum Treatment)

Quantum mechanically, we have seen that the Hamiltonian of an interacting two-level system 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 \begin{align} H' &= \frac{\hbar}{2} \begin{pmatrix} -\delta & \omega_{\rm{R}} \\ \omega_{\rm{R}} & \delta \end{pmatrix}, \end{align}}

where $\delta \equiv \omega -\omega _{0}$. Let us call the "bare" i.e. unperturbed, eigenstates of the Hamiltonian 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 |\uparrow \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 |\downarrow \rangle} . If you examine the eigenstates of this Hamiltonian, you will realize that by changing the frequency of your coupling, the eigenstates will "evolve" from e.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 |\uparrow \rangle} to $|\downarrow \rangle$, or vice versa. This can be accomplished by sweeping 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} at a 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 \dot{\omega}} . This leads to the Landau-Zener problem of a sweep through an avoided crossing.

An avoided crossing. A system on one level can jump to the other if the parameter 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} that governs the energy levels is swept sufficiently rapidly.

Whether or not the system follows an energy level adiabatically depends on how rapidly the energy is changed, compared to the minimum energy separation. To cast the problem in quantum mechanical terms, imagine two non-interacting states whose energy separation, $\omega$, depends on some parameter 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} which varies linearly in time, and vanishes for some 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 x_0} . Now add a perturbation having an off-diagonal matrix element 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} which is independent of $x$, so that the energies 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 x_0} 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~V =\pm \hbar \omega_{\rm R}/2} . The probability that the system will "jump" from one adiabatic level to the other after passing through the "avoided crossing" (i.e., the probability of non-adiabatic behavior) 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 P_\text{na} =e^{-2\pi\Gamma} }

where

$\Gamma ={\frac {|\omega _{\rm {R}}|^{2}}{4{\dot {\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 \dot\omega \equiv d\omega /dt} . This result was originally obtained by Landau and Zener. The jumping of a system as it travels across an avoided crossing is called the Landau-Zener effect. Inserting the parameters for our magnetic field problem, we have

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_{\rm na} = \exp \left( -\frac{\pi}{2} \frac{\omega_R^2} {\left| \dot\omega\right|}\right) }

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 \omega^2_{\rm R} \gg \left| \dot\omega\right|} is satisfied, the probability of non-adiabatic behavior is exponentially small. This agrees with our classical analysis.

Let us now consider a diabatic sweep, where the probability to make a transition across the crossing is non-negligible, but small. This process can be described by the Schroedinger equation in a coherent way. In particular,

${\dot {a}}_{2}=i{\frac {H_{12}}{\hbar }}a_{1}\Rightarrow a_{2}={\frac {1}{2}}\omega _{\rm {R}}t\Rightarrow P\approx |a_{2}|^{2}\approx \omega _{\rm {R}}^{2}t_{\rm {eff}}^{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 t_{\rm eff}} is the effective sweep duration, or the dephasing 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 \Delta t} . This is determined 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 \Delta \omega \Delta t \approx \pi, }

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 \Delta \omega = \dot{\omega} \Delta t} . This then gives us $\Delta t\approx {\sqrt {1/{\dot {\omega }}}}$. Then we get

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 \approx \omega^2_{\rm R} t^2_{\rm eff} = \omega^2_{\rm R}/\dot{\omega}. }

Note that the sum of the probability to undergo a transition and the probability to not undergo a transition must equal one. This means 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 = 1 - P_{\rm na}} . Using our expression for $P_{\rm {na}}$ we obtained above, we find

$P=1-P_{\rm {na}}=2\pi \Gamma \approx {\frac {\omega _{\rm {R}}^{2}}{\dot {\omega }}}.$

Therefore we've determined the probability to undergo a diabatic transition two different ways and see that they agree.

Incidentally the "rapid" in adiabatic rapid passage is something of a misnomer. The technique was originally developed in nuclear magnetic resonance in which thermal relaxation effects destroy the spin polarization if one does not invert the population sufficiently rapidly. In the absence of such relaxation processes one can take as long as one pleases to traverse the anticrossings, and the slower the crossing the less the probability of jumping.

Decoherence Processes, Mixed States, Density Matrices

The purely Hamiltonian evolution described by the Schrodinger Equation leaves the system in a pure state. However, often we have to deal with incoherent processes such as the uncontrolled loss or addition of atoms, or other perturbations that change the system evolution in an uncontrollable way. If the incoherent process is merely a (state-dependent) loss of atoms to a third state then we can describe the system by a Hamiltonian with complex eigenstates

$E_{g,e}\rightarrow E_{g,e}+i{\frac {\Gamma _{g,e}t}{2}}$

The decay of the norm $\left\langle \psi \vert \psi \right\rangle$ corresponds to decay out of the $\left|g\right\rangle$, $\left|e\right\rangle$ system, as in <xr id="fig:decay_out"/>.

<figure id="fig:decay_out">

Decay out of 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 \left|e\right\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 \left|g\right\rangle} system, which can be modelled by a non-Hermitian Hamiltonian, the imaginary part of the eigenvalue corresponding to a decay rate. </figure>

However, other processes, such as spontaneous decay 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 \left|e\right\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 \left|g\right\rangle } or loss of coherence (well-defined phase relationship) between $\left|e\right\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 \left|g\right\rangle } due to uncontrolled level shifts (e.g. collisions, uncontrolled B-field fluctuations) cannot be dealt with within the Hamiltonian approach and require the density matrix formalism.

Density Operator

The density operator formalism is necessary when the preparation process does not prepare a single quantum 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 \left|\psi\right\rangle } , but a statistical mixture of quantum states 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 \left|\psi_i\right\rangle } with probabilities $P_{i}$. The density operator is the projection operator for that statistical mixture, or ensemble average.

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{\rho}=\Sigma P_i \left|\psi_i\right\rangle \left\langle \psi_i\right| =\overline{ \left|\psi\right\rangle \left\langle \psi\right| } }

From the Schrodinger equation it follows that the Hamiltonian evolution of the density operator is governed 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 i\hbar\dot{\hat{\rho}}=\left[\hat H,\hat\rho\right] }

In the presence of damping and decay processes there are additional terms, and the time evolution is governed by a so-called Master equation. The expectation value of any operator is 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 \left\langle \hat O\right\rangle = Tr \left(\hat O\hat \rho\right)= Tr \left(\hat \rho\hat O\right) }

The expectation value of the unity operator must be 1, so $Tr\left({\hat {\rho }}\right)=1$ The expectation value 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 \hat\rho} , 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 Tr \left(\hat \rho^2\right)} , is unity for a pure state , 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 <1} for a mixed state. In any basis $\left\{\left|n\right\rangle \right\}$, 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\rho} is specified in terms of its matrix 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_{n^\prime n}= \left\langle n^\prime\right| \hat \rho \left|n\right\rangle }

and the expectation value of any operator can be written as

${\begin{aligned}\left\langle {\hat {O}}\right\rangle &=Tr\left({\hat {O}}{\hat {\rho }}\right)=\Sigma _{n}\left\langle n\right|{\hat {O}}{\hat {\rho }}\left|n\right\rangle =\Sigma _{n^{\prime }n}\left\langle n\right|{\hat {O}}\left|n^{\prime }\right\rangle \left\langle n^{\prime }\right|{\hat {\rho }}\left|n\right\rangle \\\left\langle {\hat {O}}\right\rangle &=\Sigma _{nn^{\prime }}O_{nn^{\prime }}\rho _{n^{\prime }n}{\text{ with }}O_{nn^{\prime }}=\left\langle n\right|{\hat {O}}\left|n^{\prime }\right\rangle \end{aligned}}$

The diagonal elements 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_{n^\prime n}} ) are called the probabilities since 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_{nn}} is the probability to find the system in state $n$, while the off-diagonal elements are called coherences (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_{n^\prime n}} is the coherence between states 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^\prime} and $n$). If we 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 \rho_{n^\prime n}=Ae^{i\phi}} 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 A} , $\phi$ real 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 \phi_{n^\prime n}} is the phase between states 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^\prime} and $n$, while 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 A\leq\sqrt{\rho_{nn}\rho_{n^\prime n^\prime}}} . The coherence is maximum 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 A=\sqrt{\rho_{nn}\rho_{n^\prime n^\prime}}} . Finally note that the density operator is hermitian, so 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_{n^\prime n} = \rho_{nn^\prime}^\star} , as we would expect for a meaningfully-defined coherence between two states.

Density Matrix for a Two-Level System and Bloch Vector

We now introduce a slightly different parametrization of the two-level Hamiltonian <xr id="eq:dressed_atom_hamiltonian"/> and of the corresponding density matrix:

${\begin{aligned}H&={\frac {\hbar }{2}}{\begin{pmatrix}\omega _{0}&\omega _{R}e^{-i\omega t}\\\omega _{R}e^{i\omega t}&-\omega _{0}\end{pmatrix}}={\frac {\hbar }{2}}{\begin{pmatrix}\omega _{0}&V_{1}-iV_{2}\\V_{1}+iV_{2}&-\omega _{0}\end{pmatrix}}\\&={\frac {\hbar }{2}}\left(V_{1}{\hat {\sigma }}_{x}+V_{2}{\hat {\sigma }}_{y}+\omega _{0}{\hat {\sigma }}_{z}\right)\\\rho &={\begin{pmatrix}\rho _{11}&\rho _{12}\\\rho _{21}&\rho _{22}\end{pmatrix}}={\frac {1}{2}}{\begin{pmatrix}r_{0}+r_{3}&r_{1}-ir_{2}\\r_{1}+ir_{2}&r_{0}-r_{3}\end{pmatrix}}\\&={\frac {1}{2}}\left(r_{0}{\hat {1}}+r_{1}{\hat {\sigma }}_{x}+r_{2}{\hat {\sigma }}_{y}+r_{3}{\hat {\sigma }}_{z}\right)\end{aligned}}$

In the problem set you will show (or have shown) that the Hamiltonian 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 i\hbar\dot \rho=\left[H_1,\rho\right] }

implies the equation of motion

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 \frac{d {\bf{r}} }{dt} = {\bf{\omega}} \times {\bf{r}} }

for the vectors 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 {\bf{r}} =\left(r_1,r_2,r_3\right)} and ${\bf {\omega }}=\left(V_{1},V_{2},\omega _{0}\right)$ This slightly generalizes our previous result: for a pure state using the Heisenberg equations of motion we had shown that the spin-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/2} and the corresponding magnetic moment obey

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 \frac{d {\bf{S}} }{dt}= {\bf{\omega}} \times {\bf{S}} }

but we have now generalized this result to mixed states, i.e. to ensemble averages that have less than the maximum possible magnetic moment.

We have just derived the theorem by Feynman, Vernon and Hellwarth (1957) which states that the time evolution of the density matrix for the most general two-level system is isomorphic to the behavior of a classical moment in a suitable time-dependent magnetic field.

Note that Hamiltonian evolution does not change the purity of a state: a pure state always remains pure, no matter how violently it is rotated, and a mixed state retains its degree of mixedness (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 Tr \rho^2(t)=\text{const.}} ) Here we can check this explicitly by noting that

$Tr\rho ^{2}={\frac {1}{2}}\left(r_{0}^{2}+r_{1}^{2}+r_{2}^{2}+r_{3}^{2}\right)={\frac {1}{2}}\left(r_{0}^{2}+\left|{\bf {r}}\right|^{2}\right)$

Since ${\dot {\bf {r}}}\perp {\bf {r}}$, the length of ${\bf {r}}$ does not change, and hence $Tr\rho ^{2}$ is constant in time.

Phenomenological Treatment of Relaxation: Bloch Equations

Statistical mechanics tells us the form which the density matrix will ultimately take, but it does not tell us how the system will get there or how long it will take. All we know is that ultimately the density matrix will thermalize to

$\rho ^{T}={\frac {1}{Z}}e^{-H_{0}/kT},$

where Z is the partition function.

Since the interactions which ultimately bring thermal equilibrium are incoherent processes, the density matrix formulation seems like a natural way to treat them. Unfortunately in most cases these interactions are sufficiently complex that this is done phenomenologically. For example, the equation of motion for the density matrix might be modified by the addition of a damping term:

${\dot {\rho }}={\frac {1}{i\hbar }}[H,\rho ]-(\rho -\rho ^{T})/T_{e}$

which would (in the absence of a source of non-equilibrium interactions) drive the system to equilibrium with time 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 T_e} .

This equation is not sufficiently general to describe the behavior of most systems studied in resonance physics, which exhibit different decay times for the energy and phase coherence, called 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} 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} respectively.

• 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} - decay time for population differences between non-degenerate levels, e.g. for $r_{3}$ (also called the energy decay 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_2} - decay time for coherences (between either degenerate or non-degenerate states), i.e. 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 r_1} or $r_{2}$.

The reason is that, in general, it requires a weaker interaction to destroy coherence (the relative phase of the coefficients of different states) than to destroy the population difference, so some relaxation processes will relax only the phase, resulting 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 T_2 < T_1} . (caution: certain types of collisions violate this generality.)

The effects of thermal relaxation with the two decay times described above are easily incorporated into the vector model for the 2-level system since the z-component of the population vector 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 \mathbf{r}} represents the population difference 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} and $r_{y}$ represent coherences (i.e. off-diagonal matrix 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} ): The results 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 \dot{\mathbf{r}}_z = {\left( \mathbf{\omega}\times \mathbf{r} \right)_z} - {\left( r_z - r_z^T \right)}/T_1 }

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{\mathbf{r}}_{x,y} = {\left( \mathbf{\omega}\times \mathbf{r}\right)_{x,y}} - {\left(r_{x,y} - r_{x,y}^T\right)}/T_2 }

For a magnetic spin system $\mathbf {r}$ corresponds directly to the magnetic moment 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 \mu}} . The above equations were first introduced by Bloch in this context and are known as the Bloch equations.

The addition of phenomenological decay times does not generalize the density matrix enough to cover situations where atoms (possibly state-selected) are added or lost to a system. This situation can be covered by the addition of further terms 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 \dot{\rho}} . Thus a calculation on a resonance experiment in which state-selected atoms are added to a two-level system through a tube which also permits atoms to leave (e.g. a hydrogen maser) might look like

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{\rho} &= \frac{1}{i\hbar} [H, \rho] - \begin{pmatrix} (\rho_{11} - \rho_{11T})/T_1 & \rho_{12}/T_2 \\ \rho_{21}/T_2 & (\rho_{22}- \rho_{22T})/T_1 \end{pmatrix} +R \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix} - \rho /T_{\rm escape} - \rho /T_{\rm collision} \end{align} }

where R is the rate of addition of state-selected atoms. The last two terms express effects of atom escape from the system and of collisions (e.g. spin exchange) that can't easily be incorporated in $T_{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 T_2} .

The terms representing addition or loss of atoms will not have zero trace, and consequently will not maintain TrFailed 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 ) = 1} . Physically this is reasonable for systems which gain or lose atoms; the application of the density matrix to this case shows its power to deal with complicated situations. In most applications of the above equation, one looks for a steady state solution (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 \dot{\rho} = 0} ), so this does not cause problems.

References

• Armani2003: D. K. Armani, T. J. Kippenberg, S. M. Spillane, and K. J. Vahala. Ultra-high-Q toroid microcavity on a chip. Nature, 421:925–928, Feb 2003.
• codata2002: Peter J. Mohr, Barry N. Taylor, and David B. Newell. CODATA recommended values of the fundamental physical constants: 2002. Rev. Mod. Phys., 77, Jan 2005. Available online at http://physics.nist.gov/constants.
• Dicke1954: R. H. Dicke. Coherence in spontaneous radiation processes. Phys. Rev., 93(1):99, Jan 1954.
• Feynman1963: Richard Feynman, Robert B. Leighton, and Matthew L. Sands. The Feynman Lectures on Physics. Addison-Wesley, 1963. 3 volumes.
• Flambaum2004: V. V. Flambaum, D. B. Leinweber, A. W. Thomas, and R. D. Young. Limits on variations of the quark masses, QCD scale, and fine structure constant. Phys. Rev. D, 69(11):115006, 2004.
• Ramsey1949: Norman F. Ramsey. A new molecular beam resonance method. Phys. Rev., 76(7):996, Oct 1949.
• Ramsey1950: Norman F. Ramsey. A molecular beam resonance method with separated oscillating fields. Phys. Rev., 78(6):695–699, Jun 1950.
• Reinhold2006: E. Reinhold, R. Buning, U. Hollenstein, A. Ivanchik, P. Petitjean, and W. Ubachs. Indication of a cosmological variation of the proton-electron mass ratio based on laboratory measurement and reanalysis of h2 spectra. Phys. Rev. Lett., 96:151101, 2006.
• Santarelli1999: G. Santarelli, Ph. Laurent, P. Lemonde, A. Clairon, A. G. Mann, S. Chang, A. N. Luiten, and C. Salomon. Quantum projection noise in an atomic fountain: A high-stability cesium frequency standard. Phys. Rev. Lett., 82:4619, 1999.
• Uhlenbeck1926: G. E. Uhlenbeck and S. Goudsmit. Spinning electrons and the structure of spectra. Nature, 117(2938):264, Feb 1926.
• Wineland1992: D. J. Wineland, J. J. Bollinger, W. M. Itano, F. L. Moore, and D. J. Heinzen. Spin squeezing and reduced quantum noise in spectroscopy. Phys. Rev. A, 46(11):R6797–R6800, Dec 1992.
• Wineland1994: D. J. Wineland, J. J. Bollinger, W. M. Itano, and D. J. Heinzen. Squeezed atomic states and projection noise in spectroscopy. Phys. Rev. A, 50(1):67–88, Jul 1994.