View source for Quantized spin in a magnetic field

=== Equation of Motion for the Expectation Value ===
For the system we have been considering, the Hamiltonian is
:<math>
\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
</math>
Recalling the Heisenberg equation of motion for any operator
<math>\hat{O}</math> is
:<math>
\frac{d}{dt}\hat{O}
=\frac{i}{\hbar}\left[\hat{H},\hat{O}\right]+\frac{\partial\hat{O}}{\partial t}
</math>
where the last term refers to operators with an explicit time dependence, we
have in this instance
:<math>
\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]
</math>
Using <math>\left[\hat{L}_k,\hat{L}_l\right]=i\hbar\epsilon_{klm}\hat{L}_m</math> with the
Levi-Civita symbol <math>\epsilon_{klm}</math>, we have
:<math>\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}</math>
or in short
:<math>\begin{align} 
  \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{align}</math>
These are just like the classical equations of motion
\ref{eq:classical_precession_in_static_field}, but here they describe the
precession of the operator for the magnetic moment <math> {\bf{\hat\mu}} </math> or for the
angular momentum <math> {\bf{\hat L}} </math> about the magnetic field at the (Larmor)
angular frequency <math>\omega_L=\gamma B_0</math>.

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
:<math>
    \langle{\dot{ {\bf{\hat L}} }}\rangle =\gamma\langle{ {\bf{\hat L}} }\rangle\times {\bf{B}} 
</math>

* We have not made use of any special relations for a spin-<math>1/2</math> 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 <math>l</math>.

* A spin-<math>1/2</math> 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 <math>F>1/2</math>, 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 <math>N</math> two-level atoms are coupled symmetrically to an external field.  In this case we have an effective angular momentum <math>L=N/2</math> for the symmetric coupling:

<blockquote>
::[[Image:resonances-Dicke-states.png|thumb|600px|none|]]
Level structure diagram for <math>n</math> two-level atoms in a basis of symmetric states \cite{Dicke1954}.  The leftmost column corresponds to an effective spin-<math>n/2</math> object.  Other columns correspond to manifolds of symmetric states of the <math>n</math> atoms with lower total effective angular momentum.
</blockquote>
: Again the equation of motion for the composite angular momentum <math> {\bf{\hat L}} </math> is a precession.  This is the problem considered in Dicke's famous paper \cite{Dicke1954}, 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.

<blockquote>
::[[Image:resonances-two-level.png|thumb|500px|none|]]
{fig:two_level_spin_half} Equivalence of two-level system with spin-<math>1/2</math>.  Note that for
<math>\gamma>0</math> the spin-up state (spin aligned with field) is the ground state.  For an electron, with <math>\gamma<0</math>, spin-up is the excited state.  Be careful, as both conventions are used in the literature.
</blockquote>

We have that
:<math>
 \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)
</math>
where in the last equation we have used the fact that <math>P_e+P_g=1</math>.  The signs
are chosen for a spin with <math>\gamma>0</math>, such as a proton (Figure
\ref{fig:two_level_spin_half}).  For an electron, or any other spin with
<math>\gamma<0</math>, the analysis would be the same but for the opposite sign of <math>S_z</math>
and the corresponding exchange of <math>\uparrow</math> and <math>\downarrow</math>.  If the system is
initially in the ground state, <math>P_g(t=0)=1</math> (or the spin along
<math>+ {\bf{\hat e}} _z</math>, <math>S_z(0)=+\hbar/2</math>), the expectation value obeys the classical
equation of motion \ref{eq:classical_rabi_flopping}:
:{eqn:rabi_transition_probability}
:<math>\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}</math>

Equation \ref{eqn:rabi_transition_probability} is the probability to find system
in the excited state at time <math>t</math> if it was in ground state at time <math>t=0</math>.
Figure \ref{fig:rabi_signal} shows a real-world example of such an oscillation.

<blockquote>
::[[Image:resonances-Rabi.png|thumb|400px|none|]]
{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.
</blockquote>

=== Matrix form of Hamiltonian ===
With the matrix representation
:<math>\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}</math>
we can write the Hamiltonian <math>H_0</math> associated with the static
field <math> {\bf{B}} _0 = B_0  {\bf{\hat e}} _z</math> as
:<math>
  \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
</math>
where <math>\omega_0=\gamma B_0</math> is the Larmor frequency, and
:<math>
\hat{\sigma}_z=
  \begin{pmatrix}
      1 & 0 \\
      0 & -1
  \end{pmatrix}
</math>
is a Pauli spin matrix.  The eigenstates are <math> \left|e\right\rangle </math>, <math> \left|g\right\rangle </math> with
eigenenergies <math>\pm\hbar\omega_0/2</math>.  A spin initially
along <math> {\bf{\hat e}} _x</math>, corresponding to
:<math>
 \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)
</math>
evolves in time as
:<math>
 \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)
</math>
which describes a precession with angluar frequency <math>\omega_0</math>.
The field <math> {\bf{B}} _1</math>, rotating at <math>\omega</math> in the <math>X-Y</math> plane corresponds to
:<math>\begin{align} 
  \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{align}</math>
where have used the Pauli spin matrices <math>\sigma_x</math>, <math>\sigma_y</math>. The full
Hamiltonian is thus given by

:{eq:dressed_atom_hamiltonian}
:<math>
\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}
</math>
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 Interaction Representation ===

The interaction representation consists of expanding the state <math> \left|\psi(t)\right\rangle </math>
in terms of the eigenstates <math> \left|e\right\rangle </math>, <math> \left|g\right\rangle </math> of the Hamiltonian <math>H_0</math>,
including their known time dependence <math>e^{iH_0 t/\hbar}</math> due to <math>H_0</math>.
That means we write here
:<math>
 \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}
</math>
Substituting this into the Schrodinger equation
:<math>
i\hbar\frac{d}{dt}  \left|\psi(t)\right\rangle =\left(H_0+H_1\right) \left|\psi(t)\right\rangle 
</math>
then results in the equations of motion for the coefficients
:<math>\begin{align} 
i\dot{a}_e&= \frac{\omega_R}{2}e^{i\left(\omega_0-\omega\right)t}a_g \\
i\dot{a}_g&= \frac{\omega_R}{2}e^{-i\left(\omega_0-\omega\right)t}a_e
\end{align}</math>
Where we have used the matrix form of the Hamiltonian,
\ref{eq:dressed_atom_hamiltonian}.  Introducing the detuning
<math>\delta=\omega-\omega_0</math>, we have
:<math>\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}</math>
The explicit time dependence can be eliminated by the sustitution
:<math>\begin{align} 
b_g&= e^{-i\delta t/2}a_g \\
b_e&= e^{i\delta t/2}a_e
\end{align}</math>
As you will show (or have shown) in the problem set, this leads to solutions for
<math>b_{g,e}</math> given by
:<math>\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 t/2}B_1
        +\frac{\omega_R}{\delta+\Omega_R}e^{-i\Omega_R t/2}B_2
\end{align}</math>
with two constants that are determined by the initial conditions.
For <math>a_g(t=0)=1</math> we find
:<math>
 \left| a_e(t)\right| ^2=\frac{\omega_R^2}{\Omega_R^2}\sin^2\frac{\Omega_R t}{2}
</math>
as already derived from the fact that the expectation value for
the magnetic moment obeys the classical equation.


=== "Rapid" Adiabiatic Passage (quantum treatment) ===
(in the works)

The quantum mechanical treatment yields a probability for non-adiabatic transition (probability for the
magnetic moment not following the magnetic field) given by
:<math>
P_\text{na}=exp\left(-\frac{\pi}{2}\frac{\omega^2_R}{ \left| \dot\omega\right| }\right)
</math>
in agreement with the qualitative classical criterion derived above.