Relaxation (NMR)

In nuclear magnetic resonance (NMR) spectroscopy and magnetic resonance imaging (MRI) the term relaxation describes several processes by which nuclear magnetization prepared in a non-equilibrium state return to the equilibrium distribution. In other words, relaxation describes how fast spins "forget" the direction in which they are oriented. The rates of this spin relaxation can be measured in both spectroscopy and imaging applications.

T₁ and T₂

Different physical processes are responsible for the relaxation of the components of the nuclear spin magnetization vector M parallel and perpendicular to the external magnetic field, B₀ (which is conventionally oriented along the z axis). These two principal relaxation processes are termed T₁ and T₂ relaxation respectively.

T₁

Main article: Spin-lattice relaxation time

The longitudinal (or spin-lattice) relaxation time T₁ is the decay constant for the recovery of the z component of the nuclear spin magnetization, M_z, towards its thermal equilibrium value, M_z,eq. In general,

$M_z(t) = M_{z,\mathrm{eq}} - [M_{z,\mathrm{eq}} - M_z(0)]e^{-t/T_1}$

In specific cases:

• If M has been tilted into the xy plane, then M_z(0) = 0 and the recovery is simply

$M_z(t) = M_{z,\mathrm{eq}}\left( 1 - e^{-t/T_1} \right)$

i.e. the magnetization recovers to 63% of its equilibrium value after one time constant T₁.

• In the inversion recovery experiment, commonly used to measure T₁ values, the initial magnetization is inverted, M_z(0) = − M_z,eq, and so the recovery follows

$M_z(t) = M_{z,\mathrm{eq}}\left( 1 - 2e^{-t/T_1} \right)$

T₁ relaxation involves redistributing the populations of the nuclear spin states in order to reach the thermal equilibrium distribution. By definition this is not energy conserving. Moreover, spontaneous emission is negligibly slow at NMR frequencies. Hence truly isolated nuclear spins would show negligible rates of T₁ relaxation. However, a variety of relaxation mechanisms allow nuclear spins to exchange energy with their surroundings, the lattice, allowing the spin populations to equilibrate. The fact that T₁ relaxation involves an interaction with the surroundings is the origin of the alternative description, spin-lattice relaxation.

Note that the rates of T₁ relaxation are generally strongly dependent on the NMR frequency and so vary considerably with magnetic field strength B. Small amounts of paramagnetic substances in a sample speed up relaxation very much. By degassing, and thereby removing dissolved Oxygen, the T₁/T₂ of liquid samples easily go up to an order of ten seconds.

T₂

Main article: Spin-spin relaxation time

The transverse (or spin-spin) relaxation time T₂ is the decay constant for the component of M perpendicular to B₀, designated M_xy, M_T, or $M_{\perp}$. For instance, initial xy magnetisation at time zero will decay to zero (i.e. equilibrium) as follows:

$M_{xy}(t) = M_{xy}(0) e^{-t/T_2} \,$

i.e. the transverse magnetization vector drops to 37% of its original magnitude after one time constant T₂.

T₂ relaxation is a complex phenomenon, but at its most fundamental level, it corresponds to a decoherence of the transverse nuclear spin magnetization. Random fluctuations of the local magnetic field lead to random variations in the instantaneous NMR precession frequency of different spins. As a result, the initial phase coherence of the nuclear spins is lost, until eventually the phases are disordered and there is no net xy magnetization. Because T₂ relaxation involves only the phases of other nuclear spins it is often called "spin-spin" relaxation.

T₂ values are generally much less dependent on field strength, B, than T₁ values.

A Hahn echo decay experiment can be used to measure the T₂ time, as shown in the animation below. The size of the echo is recorded for different spacings of the two applied pulses. This reveals the decoherence which is not refocused by the 180° pulse. In simple cases, an exponential decay is measured which is described by the T₂ time.

T₂* and magnetic field inhomogeneity

In an idealized system, all nuclei in a given chemical environment, in a magnetic field, precess with the same frequency. However, in real systems, there are minor differences in chemical environment which can lead to a distribution of resonance frequencies around the ideal. Over time, this distribution can lead to a dispersion of the tight distribution of magnetic spin vectors, and loss of signal (Free Induction Decay). In fact, for most magnetic resonance experiments, this "relaxation" dominates. This results in dephasing.

However, decoherence because of magnetic field inhomogeneity is not a true "relaxation" process; it is not random, but dependent on the location of the molecule in the magnet. For molecules that aren't moving, the deviation from ideal relaxation is consistent over time, and the signal can be recovered by performing a spin echo experiment.

The corresponding transverse relaxation time constant is thus T₂^*, which is usually much smaller than T₂. The relation between them is:

$\frac{1}{T_2^*}=\frac{1}{T_2}+\frac{1}{T_{inhom}} = \frac{1}{T_2}+\gamma \Delta B_0$

where γ represents gyromagnetic ratio, and ΔB₀ the difference in strength of the locally varying field.

Unlike T₂, T₂* is influenced by magnetic field gradient irregularities. The T₂* relaxation time is always shorter than the T₂ relaxation time and is typically milliseconds for water samples in imaging magnets.

Is T₁ always longer than T₂?

The following always holds true^[1]: $2 T_1 \ge T_2 \ge T_2^*$. In most situations (but not in principle) T₁ is greater than T₂.

Bloch equations

Main article: Bloch equations

Bloch equations are used to calculate the nuclear magnetization M = (M_x, M_y, M_z) as a function of time when relaxation times T₁ and T₂ are present. Bloch equations are phenomenological equations that were introduced by Felix Bloch in 1946.^[2]

$\frac {\partial M_x(t)} {\partial t} = \gamma ( \bold {M} (t) \times \bold {B} (t) ) _x - \frac {M_x(t)} {T_2}$

$\frac {\partial M_y(t)} {\partial t} = \gamma ( \bold {M} (t) \times \bold {B} (t) ) _y - \frac {M_y(t)} {T_2}$

$\frac {\partial M_z(t)} {\partial t} = \gamma ( \bold {M} (t) \times \bold {B} (t) ) _z - \frac {M_z(t) - M_0} {T_1}$

Where γ is the gyromagnetic ratio and B(t) = (B_x(t), B_y(t), B₀ + B_z(t)) is the magnetic flux density experienced by the nuclei. The z component of the magnetic flux density B is typically composed of two terms: one, B₀, is constant in time, the other one, B_z(t), is time dependent. It is present in magnetic resonance imaging and helps with the spatial decoding of the NMR signal. M(t) × B(t) is the cross product of these two vectors.

The equation listed above in the section on T₁ and T₂ relaxation can be derived from Bloch equations.

Common relaxation time constants in human tissues

Following is a table of the approximate values of the two relaxation time constants for nonpathological human tissues, just for simple reference.

At a main field of 1.5 T

Tissue Type	Approximate T1 value in ms	Approximate T2 value in ms
Adipose tissues	240-250	60-80
Whole blood (deoxygenated)	1350	50
Whole blood (oxygenated)	1350	200
Cerebrospinal fluid (similar to pure water)	4200 - 4500	2100-2300
Gray matter of cerebrum	920	100
White matter of cerebrum	780	90
Liver	490	40
Kidneys	650	60-75
Muscles	860-900	50

Following is a table of the approximate values of the two relaxation time constants for chemicals that commonly show up in human brain magnetic resonance spectroscopy (MRS) studies, physiologically or pathologically.

At a main field of 1.5 T

Signals of Chemical Groups	Relative resonance frequency	Approximate T1 value (ms)	Approximate T2 value (ms)
Creatine (Cr) and Phosphocreatine (PCr)[3]	3.0 ppm	gray matter: 1150-1340, white matter: 1050-1360	gray matter: 198-207, white matter: 194-218
N-Acetyl group (NA), mainly from N-Acetylaspartate (NAA)[3]	2.0 ppm	gray matter: 1170-1370, white matter: 1220-1410	gray matter: 388-426, white matter: 436-519
—CH3 group of Lactate[4]	1.33 ppm (doublet: 1.27 & 1.39 ppm)	(To be listed)	1040

Relaxation in the rotating frame, T_1ρ

The discussion above describes relaxation of nuclear magnetization in the presence of a constant magnetic field B₀. This is called relaxation in the laboratory frame. Another technique, called relaxation in the rotating frame, is the relaxation of nuclear magnetization in the presence of the field B₀ together with a time-dependent magnetic field B₁. The field B₁ rotates in the plane perpendicular to B₀ at the Larmor frequency of the nuclei in the B₀. The magnitude of B₁ is typically much smaller than the magnitude of B₀. Under these circumstances the relaxation of the magnetization is similar to laboratory frame relaxation in a field B₁. The decay constant for the recovery of the magnetization component along B₁ is called the spin-lattice relaxation time in the rotating frame and is denoted T_1ρ. Relaxation in the rotating frame is useful because it provides information on slow motions of nuclei.

Microscopic mechanisms

Relaxation of nuclear spins requires a microscopic mechanism for a nucleus to change orientation with respect to the applied magnetic field and/or interchange energy with the surroundings (called the lattice). The most common mechanism is the magnetic dipole-dipole interaction between the magnetic moment of a nucleus and the magnetic moment of another nucleus or other entity (electron, atom, ion, molecule). This interaction depends on the distance between the pair of dipoles (spins) but also on their orientation relative to the external magnetic field. Several other relaxation mechanisms also exist. The chemical shift anisotropy (CSA) relaxation mechanism arises whenever the electronic environment around the nucleus is non spherical, the magnitude of the electronic shielding of the nucleus will then be dependent on the molecular orientation relative to the (fixed) external magnetic field. The spin rotation (SR) relaxation mechanism arises from an interaction between the nuclear spin and a coupling to the overall molecular rotational angular momentum. Nuclei with spin I ≥ 1 will have not only a nuclear dipole but a quadrupole. The nuclear quadrupole has an interaction with the electric field gradient at the nucleus which is again orientation dependent as with the other mechanisms described above, leading to the so called quadrupolar relaxation mechanism.

Molecular reorientation or tumbling can then modulate these orientation-dependent spin interaction energies. According to quantum mechanics, time-dependent interaction energies cause transitions between the nuclear spin states which result in nuclear spin relaxation. The application of time-dependent perturbation theory in quantum mechanics shows that the relaxation rates (and times) depend on spectral density functions that are the Fourier transforms of the autocorrelation function of the fluctuating magnetic dipole interactions^[5]. The form of the spectral density functions depend on the physical system, but a simple approximation called the BPP theory is widely used.

Another relaxation mechanism is the electrostatic interaction between a nucleus with an electric quadrupole moment and the electric field gradient that exists at the nuclear site due to surrounding charges. Thermal motion of a nucleus can result in fluctuating electrostatic interaction energies. These fluctuations produce transitions between the nuclear spin states in a similar manner to the magnetic dipole-dipole interaction.

BPP theory

In 1948, Nicolaas Bloembergen, Edward Mills Purcell, and Robert Pound proposed the so-called Bloembergen-Purcell-Pound theory (BPP theory) to explain the relaxation constant of a pure substance in correspondence with its state, taking into account the effect of tumbling motion of molecules on the local magnetic field disturbance^[6]. The theory was in good agreement with experiments on pure substances, but not for complicated environments such as the human body.

This theory makes the assumption that the autocorrelation function of the microscopic fluctuations causing the relaxation is proportional to $e^{-t/\tau_c}$, where τ_c is called the correlation time. From this theory, one can get T₁、T₂ for magnetic dipolar relaxation:

$\frac{1}{T_1}=K[\frac{\tau_c}{1+\omega_0^2\tau_c^2}+\frac{4\tau_c}{1+4\omega_0^2\tau_c^2}]$

$\frac{1}{T_2}=\frac{K}{2}[3\tau_c+\frac{5\tau_c}{1+\omega_0^2\tau_c^2}+\frac{2\tau_c}{1+4\omega_0^2\tau_c^2}]$,

where ω₀ is the Larmor frequency in correspondence with the strength of the main magnetic field B₀. τ_c is the correlation time of the molecular tumbling motion. $K=\frac{3\mu^2}{160\pi^2}\frac{\hbar^2\gamma^4}{r^6}$ is a constant with μ being the magnetic dipole moment of the spin-1/2 nuclei, $\hbar=\frac{h}{2\pi}$ the reduced Planck constant, γ the gyromagnetic ratio of such species of nuclei, and r the distance between the two nuclei carrying magnetic dipole moment.

Taking for example the H₂O molecules in liquid phase without the contamination of oxygen-17, the value of K is 1.02×10¹⁰ s^-2 and the correlation time τ_c is on the order of picoseconds = 10 ^{− 12} s, while hydrogen nuclei ¹H (protons) at 1.5 teslas carry an Larmor frequency of approximately 64 MHz. We can then estimate using τ_c = 5×10^-12 s:

$\omega_0\tau_c = 3.2\times 10^{-5}$(dimensionless)

$T_1=(1.02\times 10^{10}[\frac{ 5\times 10^{-12} }{1 + (3.2\times 10^{-5} )^2} + \frac{ 4\cdot 5\times 10^{-12} }{1 + 4\cdot (3.2\times 10^{-5} )^2}])^{-1}$= 3.92 s

$T_2=(\frac{1.02\times 10^{10}}{2}[3\cdot 5\times 10^{-12} + \frac{5\cdot 5\times 10^{-12} }{1 + (3.2\times 10^{-5} )^2} + \frac{ 2\cdot 5\times 10^{-12} }{1 + 4\cdot (3.2\times 10^{-5} )^2}])^{-1}$= 3.92 s,

which is close to the experimental value, 3.6 s. Meanwhile, we can see that at this extreme case, T₁ equals T₂.

References

1. ^ Malcolm H. Levitt: Spin Dynamics: Basics of Nuclear Magnetic Resonance, 2nd edition, John Wiley & Sons, New York 2008, ISBN 0-470-51117-6, Section 11.9.2
2. ^ F Bloch, Nuclear Induction, Physical Review 70, 460-473 (1946)
3. ^ ^{a b} Chemicals of brain relaxation time at 1.5T. Kreis R, Ernst T, and Ross BD "Absolute Quantification of Water and Metabolites in the Human Brain. II. Metabolite Concentrations" Journal of Magnetic Resonance, Series B 102 (1993): 9-19
4. ^ Lactate relaxation time at 1.5 T. Isobe T, Matsumura A, Anno I, Kawamura H, Muraishi H, Umeda T, Nose T. "Effect of J coupling and T2 Relaxation in Assessing of Methyl Lactate Signal using PRESS Sequence MR Spectroscopy." Igaku Butsuri (2005) v25. 2:68-74.
5. ^ A. Abragam "Principles of Nuclear Magnetism" (Oxford University Press, 1961)
6. ^ Bloembergen, E.M. Purcell, R.V. Pound "Relaxation Effects in Nuclear Magnetic Resonance Absorption" Physical Review (1948) v73. 7:679-746

See also

• basics of NMR
• Relaxation in high-resolution NMR spectroscopy
• Field-cycling NMR relaxometry
• Relaxometry