Magnetic resonance (quantum mechanics) explained

In quantum mechanics, magnetic resonance is a resonant effect that can appear when a magnetic dipole is exposed to a static magnetic field and perturbed with another, oscillating electromagnetic field. Due to the static field, the dipole can assume a number of discrete energy eigenstates, depending on the value of its angular momentum (azimuthal) quantum number. The oscillating field can then make the dipole transit between its energy states with a certain probability and at a certain rate. The overall transition probability will depend on the field's frequency and the rate will depend on its amplitude. When the frequency of that field leads to the maximum possible transition probability between two states, a magnetic resonance has been achieved. In that case, the energy of the photons composing the oscillating field matches the energy difference between said states. If the dipole is tickled with a field oscillating far from resonance, it is unlikely to transition. That is analogous to other resonant effects, such as with the forced harmonic oscillator. The periodic transition between the different states is called Rabi cycle and the rate at which that happens is called Rabi frequency. The Rabi frequency should not be confused with the field's own frequency. Since many atomic nuclei species can behave as a magnetic dipole, this resonance technique is the basis of nuclear magnetic resonance, including nuclear magnetic resonance imaging and nuclear magnetic resonance spectroscopy.

Quantum mechanical explanation

\tfrac{1}{2}

system such as a proton; according to the quantum mechanical state of the system, denoted by

|\Psi(t)\rangle

, evolved by the action of a unitary operator

e-i{\hat/\hbar}

; the result obeys Schrödinger equation:

i\hbar

\partial
\partialt

\Psi=\hatH\Psi

States with definite energy evolve in time with phase

e-iEt/\hbar

, (

|\Psi(t)\rangle=|\Psi(0)\ranglee-iEt/\hbar

) where E is the energy of the state, since the probability of finding the system in state

|\langlex|\Psi(t)\rangle|2

=

|\langlex|\Psi(0)\rangle|2

is independent of time. Such states are termed stationary states, so if a system is prepared in a stationary state, (i.e. one of the eigenstates of the Hamiltonian operator), then P(t) = 1, i.e. it remains in that state indefinitely. This is the case only for isolated systems. When a system in a stationary state is perturbed, its state changes, so it is no longer an eigenstate of the system's complete Hamiltonian. This same phenomenon happens in magnetic resonance for a spin

\tfrac{1}{2}

system in a magnetic field.

The Hamiltonian for a magnetic dipole

m

(associated with a spin

\tfrac{1}{2}

particle) in a magnetic field
B0=B
0\hat{z}
is:

\hat{H}=-m

B0

=-\tfrac{\hbar}{2}\gamma\sigmazB0=-\tfrac{\hbar}{2}\omega0\begin{bmatrix}1&0\\0&-1\end{bmatrix}

Here

\omega0:=\gammaB0

is the Larmor precession frequency of the dipole for
B0
magnetic field and

\sigmaz

is z Pauli matrix. So the eigenvalues of

\hat{H}

are

-\tfrac{\hbar}{2}\omega0

and

\tfrac{\hbar}{2}\omega0

. If the system is perturbed by a weak magnetic field
B1
, rotating counterclockwise in x-y plane (normal to
B0
) with angular frequency

\omega

, so that
B1=\hat{i}B
1

\cos{\omegat}-\hat{j}B1\sin{\omegat}

, then

\begin{bmatrix}1\\0\end{bmatrix}

and

\begin{bmatrix}0\\1\end{bmatrix}

are not eigenstates of the Hamiltonian, which is modified into

\hat{H}=\gamma\begin{pmatrix}B0&

i\omegat
B
1e

\B1e-&-B0\end{pmatrix}.

It is inconvenient to deal with a time-dependent hamiltonian. To make

\hat{H}

time-independent requires a new reference frame rotating with
B1
, i.e. rotation operator

\hat{R}(t)

on

|\Psi(t)\rangle

, which amounts to basis change in Hilbert space. Using this on Schrödinger's equation, the Hamiltonian becomes:

\hat{H\prime}=\hat{R}(t)\hat{H}\hat{R}(t)\dagger+\tfrac{\hbar}{2}\omega\sigmaz

Writing

\hat{R}(t)

in the basis of

\sigmaz

as-

\hat{R}(t)=\begin{pmatrix}e-i{\omegat/2}&0\\0&ei{\omegat/2}\end{pmatrix}

Using this form of the Hamiltonian a new basis is found:
\prime}=\tfrac{\hbar}{2}\begin{pmatrix}\Delta\omega&-\omega
\hat{H
1\\-\omega

1&-\Delta\omega\end{pmatrix}

where

\Delta\omega=\omega-\omega0

and

\omega1=\gammaB1

This Hamiltonian is exactly similar to that of a two state system with unperturbed energies

\tfrac{\hbar}{2}\Delta\omega

and

-\tfrac{\hbar}{2}\Delta\omega

with a perturbation expressed by

\tfrac{\hbar}{2}\begin{pmatrix}0&-\omega1\\-\omega1&0\end{pmatrix}

; According to Rabi oscillation, starting with

\begin{bmatrix}1\\0\end{bmatrix}

state, a dipole in parallel to
B0
with energy

-\tfrac{\hbar}{2}\omega0

, the probability that it will transit to

\begin{bmatrix}0\\1\end{bmatrix}

state (i.e. it will flip) is

Now consider

\omega=\omega0

, i.e. the
B1
field oscillates at the same rate the dipole exposed to the
B0
field does. That is a case of resonance. Then at specific points in time, namely

t=

(2n+1)\pi
2+\Delta\omega
\sqrt{\omega2
1
}, the dipole will flip, going to the other energy eigenstate

\begin{bmatrix}0\\1\end{bmatrix}

with a 100% probability. When

\omega\not=\omega0

, the probability of change of energy state is small. Therefore, the resonance condition can be used, for instance, to measure the magnetic moment of a dipole or the magnetic field at a point in space.

A special case to show applications

A special case occurs where a system oscillates between two unstable levels that have the same life time

\tau

.[1] If atoms are excited at a constant, say n/time, to the first state, some decay and the rest have a probability

P12

to transition to the second state, so in the time interval between t and (t + dt) the number of atoms that jump to the second state from the first is

n(1-e-t/\tau)P12dt

, so at time t the number of atoms in the second state is

dN=n.e-t/\tau.(1-e-t/\tau)P12dt

=

n.e-t/\tau.P12dt

The rate of decay from state two depends on the number of atoms that were collected in that state from all previous intervals, so the number of atoms in state 2 is \int_^ ne^P_\ dt; The rate of decay of atoms from state two is proportional to the number of atoms present in that state, while the constant of proportionality is decay constant

λ

. Performing the integration rate of decay of atoms from state two is obtained as:

(n/2)\omega2/(\delta\omega

2+1/\tau
1

2)

From this expression many interesting points can be exploited, such

B0

so that

\omega0

in

\delta\omega

produces a Lorentz curve (see Cauchy–Lorentz distribution), detecting the peak of that curve, the abscissa of it gives

\omega0

, so now

\omega

(angular frequency of rotation of

B1

=

\gamma

(B0)max

, so from the known value of

\omega

and

(B0)max

, the gyromagnetic ratio

\gamma

of the dipole can be measured; by this method we can measure Nuclear spin where all electronic spins are balanced. Correct measurement of nuclear magnetic moment helps to understand the character of nuclear force.

\gamma

is known, by varying

\omega

, the value of

B0

can be obtained. This measurement technique is precise enough for use in sensitive magnetometers. Using this technique, the value of magnetic field acting at a particular lattice site by its environment inside a crystal can be obtained.
2+1/\tau
\sqrt{\omega
1

2}

, for several values of

\omega1

(i.e. of

B1

), we can plot d vs

\omega1

, and by extrapolating this line for

\omega1

, the lifetime of unstable states can be obtained from the intercept.

Rabi's method

See also: Stern–Gerlach experiment. The existence of spin angular momentum of electrons was discovered experimentally by the Stern–Gerlach experiment. In that study a beam of neutral atoms with one electron in the valence shell, carrying no orbital momentum (from the viewpoint of quantum mechanics) was passed through an inhomogeneous magnetic field. This process was not approximate due to the small deflection angle, resulting in considerable uncertainty in the measured value of the split beam.

Rabi's method was an improvement over Stern-Gerlach. As shown in the figure, the source emits a beam of neutral atoms, having spin angular momentum

\hbar/2

. The beam passes through a series of three aligned magnets. Magnet 1 produces an inhomogeneous magnetic field with a high gradient
\partialB
\partialz
(as in Stern–Gerlach), so the atoms having 'upward' spin (with

Sz=\hbar/2

) will deviate downward (path 1), i.e. to the region of less magnetic field B, to minimize energy. Atoms with 'downward' spin with

Sz=-\hbar/2

) will deviate upward similarly (path 2). Beams are passed through slit 1, to reduce any effects of source beyond. Magnet 2 produces only a uniform magnetic field in the vertical direction applying no force on the atomic beam, and magnet 3 is actually inverted magnet 1. In the region between the poles of magnet 3, atoms having 'upward' spin get upward push and atoms having 'downward' spin feel downward push, so their path remains 1 and 2 respectively. These beams pass through a second slit S2, and arrive at detector and get detected.

If a horizontal rotating field

B1

, angular frequency of rotation

\omega1

is applied in the region between poles of magnet 2, produced by oscillating current in circular coils then there is a probability for the atoms passing through there from one spin state to another (

Sz=+\hbar/2->-\hbar/2

and vice versa), when

\omega1

=

\omegap

, Larmor frequency of precession of magnetic moment in B. The atoms that transition from 'upward' to 'downward' spin will experience a downward force while passing through magnet 3, and will follow path 1'. Similarly, atoms that change from 'downward' to 'upward' spin will follow path 2', and these atoms will not reach the detector, causing a minimum in detector count. If angular frequency

\omega1

of

B1

is varied continuously, then a minimum in detector current will be obtained (when

\omega1

=

\omegap

). From this known value of

\omega1

(

=geB/{2\hbar}

, where g is 'Landé g factor'), 'Landé g-factor' is obtained which will enable one to have correct value of magnetic moment

\mu~(=gq\hbar/{4m})

. This experiment, performed by Isidor Isaac Rabi is more sensitive and accurate compared than Stern-Gerlach.

Correspondence between classical and quantum mechanical explanations

See also: Bloch equations. Spin angular momentum allows magnetic resonance phenomena to be explained via classical physics. When viewed from the reference frame attached to the rotating field, it seems that the magnetic dipole precesses around a net magnetic field

(\Delta\omega\hatz-\omega1\hatX)/\gamma

, where

\hatz

is the unit vector along uniform magnetic field

B0

and

\hatX

is the same in the direction of rotating field

B1

and

\delta\omega=\omega-\omega0

.

So when

\omega=\omega0

, a high precession amplitude allows the magnetic moment to be completely flipped. Classical and quantum mechanical predictions correspond well, which can be viewed as an example of the Bohr Correspondence principle, which states that quantum mechanical phenomena, when predicted in classical regime, should match the classical result. The origin of this correspondence is that the evolution of the expected value of magnetic moment is identical to that obtained by classical reasoning. The expectation value of the magnetic moment is

\langlem\rangle=\gamma\langleS\rangle

. The time evolution of

\langlem\rangle

is given by
i\hbard
dt

\langlem\rangle=\langle[m,\hatH]\rangle

\hatH=-mB(t)

so,

[mi,\hatH]=[mi,-mjBj]=[\gammaSi,-\gammaSj

2
B
j]=-\gamma

[Si,SjBj]=-\gamma2i\hbar[{SkBj-SjBk}],(ij,k)

So,

[mi,\hatH]=i\hbar\gamma[Bjmk-Bkmj]

andwhich looks exactly similar to the equation of motion of magnetic moment

m

in classical mechanics –
d
dt

m(t)=\gammam(t) x B(t)

This analogy in the mathematical equation for the evolution of magnetic moment and its expectation value facilitates to understand the phenomena without a background of quantum mechanics.

Magnetic resonance imaging

See also: Magnetic resonance imaging.

In magnetic resonance imaging (MRI) the spin angular momentum of the proton is used. The most available source for protons in the human body is represented by hydrogen atoms in water. A strong magnetic field

B

applied to water causes the appearance of two different energy levels for spin angular momentum,

+\gamma\hbarB/2

and

-\gamma\hbarB/2

, using

E=-\muB

.

According to the Boltzmann distribution, as the number of systems having energy

E

out of

N0

at temperature

T

is

N0e-E/kT

(where

k

is the Boltzmann constant), the lower energy level, associated with spin

\hbar/2

, is more populated than the other. In the presence of a rotating magnetic field more protons flip from

Sz=+\hbar/2

to

Sz=-\hbar/2

than the other way, causing absorption of microwave or radio-wave radiation (from the rotating field). When the field is withdrawn, protons tend to re-equilibrate along the Boltzmann distribution, so some of them transition from higher energy levels to lower ones, emitting microwave or radio-wave radiation at specific frequencies.

Instead of nuclear spin, spin angular momentum of unpaired electrons is used in EPR (electron paramagnetic resonance) in order to detect free radicals, etc.

Magnetic resonance as a quantum phenomenon

The phenomenon of magnetic resonance is rooted in the existence of spin angular momentum of a quantum system and its specific orientation with respect to an applied magnetic field. Both cases have no explanation in the classical approach and can be understood only by using quantum mechanics. Some people claim that purely quantum phenomena are those that cannot be explained by the classical approach. For example, phenomena in the microscopic domain that can to some extent be described by classical analogy are not really quantum phenomena. Since the basic elements of magnetic resonance have no classical origin, although analogy can be made with classical Larmor precession, MR should be treated as a quantum phenomenon.

See also

References

Notes and References

  1. Page-449, Quantum Mechanics, Vol.1, Claude Cohen-Tannoudji, Bernard Diu, Frank Laloe