# Solomon Equations

- Page ID
- 61683

## Introduction

The Solomon equations describe the relaxation between two coupled spins, I and S. This page will be dedicated to developing a theoretical treatment beginning from the Bloch equations. We will investigate the ramifications of the relaxation and the consequences on the overall relaxation rates.

## Theoretical Treatment

We begin by describing a system of two spins, I and S. Each spin has two states; low energy \(\alpha\) and high energy \(\beta\) and we can then come up with a 4 level energy system for these two spins. This is shown in the figure below.

**Figure 1.** Energy level diagram for a 2 spin system.

A spin can have allowed relaxation induced transitions between the \(\alpha\) and \(\beta\) states which occur at a rate (W^{1}_{IS}) as well as forbidden relaxation induced transitions between like states (W^{2}), which are shown in the figure below. Forbidden transitions W^{2} corresponds to both spins flipping, known as a flip-flip transition, while W^{0} is where each spin flips between states, known as a flip-flop transition. The superscripts denote the order (zero, single, double) type of transition.

**Figure 2.** Transition rates of relaxation induced transitions.

We can now define the change in the population of each level in with 1=\(\alpha \alpha\) 2=\(\beta \alpha\), 3= \(\alpha \beta\) and 4=\(\beta \beta\) during a given amount of time. For example, for energy level 1, the population is initially n_{1}^{0}. As time progress, the some of this population will go to energy levels 2, 3, and 4 at rates, W^{1}_{I}, W^{1}_{S}, W^{2}, respectively. Meanwhile the state will gain population from states 2, 3, and 4 at at rates, W^{1}_{I}, W^{1}_{S}, W^{2}, respectively. It may seem that this would give a net result of 0, but the populations of each state is given by the Boltzmann distribution and the argument may be easily made that the rate will depend on the number of spins in the given energy level. Therefore we may now write:

\[\dfrac{dn_1}{dt}=-W_S^1n_1-W_I^1n_1-W^2n_1+W_S^1n_3+W_I^1n_2+W^2n_4\] (eq. 1)

the rates for the other energy levels may then be written as:

\[\dfrac{dn_2}{dt}=+W_I^1n_1+W_S^2n_4+W^0n_3-W_S^2n_2-W_I^1n_2-W^0n_2\] (eq. 2)

\[\dfrac{dn_3}{dt}=+W_S^1n_1+W_I^2n_4+W^0n_2-W_I^2n_3-W_S^1n_3-W^0n_3\](eq. 3)

\[\dfrac{dn_4}{dt}=-W_S^2n_4-W_I^2n_4-W^2n_4+W_I^2n_3+W_S^2n_2+W^2n_1\](eq. 4)

We can then calculate the spin magnetization for the I and S spins, along the Z axis. The reason for the Z axis is because we are concerned with the recovery of the magnetization to it's equilibrium value. Therefore, the net magnetization for spin I will be the difference in populations between the I transitions and is expressed as

\[I_z=n_1-n_2+n_3-n_4\] (eq. 5)

and spin S will be

\[S_z=n_1-n_3+n_2-n_4.\] (eq. 6)

Further we can define the difference in populations between the 2 I-spin transitions, n_{1}-n_{2}-n_{3}+n_{4, }which is the same for the 2 S-spin transitions and will be denoted 2I_{z}S_{z}.

\[2I_zS_z=n_1-n_2-n_3+n_4.\] (eq. 7)

The total population of the system is then given by

\(T=n_1+n_2+n_3+n_4\) (eq. 8)

And n_{i} can then be re-written as a combination of I_{z}, S_{z}, 2I_{z}S_{z}, and T such that

\(n_1=\dfrac{1}{4}(T+I_z+S_z+2I_zS_z)\) (eq. 9)

\(n_2=\dfrac{1}{4}(T-I_z+S_z-2I_zS_z)\) (eq. 10)

\(n_3=\dfrac{1}{4}(T+I_z-S_z-2I_zS_z)\) (eq. 11)

\(n_4=\dfrac{1}{4}(T-I_z-S_z+2I_zS_z)\) (eq. 12)

Which can be verified by plugging equations 5-8 into equations 9-12.

### Derivation of The Solomon Equations

The Solomon equations are the time-dependent operators I_{z}, S_{z}, and 2I_{z}S_{z} expressed in terms of W and the operators. We will take the following approach to derive each operator:

1) Take the time derivative of each operator.

2) Express the derivative in terms of W by substituting equations 1-4 for \(\dfrac{dn_i}{dt}\).

3) Group the n_{i} terms.

4) Substitute the equations 9-12 for n_{i} and group the like W terms. This will give the final equations.

Each step will be numbered according to this list.

#### Derivation of I_{z}

1) \(\dfrac{dI_z}{dt}=\dfrac{dn_1}{dt}-\dfrac{dn_2}{dt}+\dfrac{dn_3}{dt}-\dfrac{dn_4}{dt}\)

2) \(=-W_S^1n_1-W_I^1n_1-W^2n_1+W_S^1n_3+W_I^1n_2+W^2n_4\)

\(-(+W_I^1n_1+W_S^2n_4+W^0n_3-W_S^2n_2-W_I^1n_2-W^0n_2)\)

\(+(+W_S^1n_1+W_I^2n_4+W^0n_2-W_I^2n_3-W_S^1n_3-W^0n_3)\)

\(-(-W_S^2n_4-W_I^2n_4-W^2n_4+W_I^2n_3+W_S^2n_2+W^2n_1)\)

3) \(\dfrac{dI_z}{dt}=-2n_1[W_I^1+W^2]+2n_2[W_I^1+W^0]-2n_3[W_I^2+W^0]+2n_4[W_I^2+W^2]\)

4) \(W^2[-2\dfrac{1}{4}*(T+I_z+S_z+2I_zS_z)+2\frac{1}{4}(T-I_z-S_z+2I_zS_z)]=W^2(-I_z-S_z)\)

\(W_I^1[-2\dfrac{1}{4}(T+I_z+S_z+2I_zS_z)+2\dfrac{1}{4}(T-I_z+S_z-2I_zS_z)=W_1^1(-I_z-2I_zS_z)\)

\(W^0[2\dfrac{1}{4}(T-I_z+S_z-2I_zS_z)+-2\dfrac{1}{4}(T+I_z-S_z-2I_zS_z)]=W^0(-I_z+S_z)\)

\(W_I^2[-2\dfrac{1}{4}(T+I_z-S_z-2I_zS_z)+2\dfrac{1}{4}(T-I_z-S_z+2I_zS_z)]=W_I^2(-I_z+2I_zS_z)\)

**\(\dfrac{dI_z}{dt}=-(W_I^1+W_I^2+W^2+W^0)I_z-(W^2-W^0)S_z-(W_I^1-W_I^2)2I_zS_z\)**

#### Derivation of S_{z}

\[S_z=n_1-n_3+n_2-n_4.\]

1) \(\dfrac{dS_z}{dt}=\dfrac{n_1}{dt}-\dfrac{n_3}{dt}+\dfrac{n_2}{dt}-\dfrac{n_4}{dt}\)

2) \(\dfrac{dS_z}{dt}=[-W_S^1n_1-W_I^1n_1-W^2n_1+W_S^1n_3+W_I^1n_2+W^2n_4]\)

\(-[+W_S^1n_1+W_I^2n_4+W^0n_2-W_I^2n_3-W_S^1n_3-W^0n_3]\)

\(+[W_I^1n_1+W_S^2n_4+W^0n_3-W_S^2n_2-W_I^1n_2-W^0n_2]\)

\(-[-W_S^2n_4-W_I^2n_4-W^2n_4+W_I^2n_3+W_S^2n_2+W^2n_1]\)

3) \(\dfrac{dS_z}{dt}=\)\(-2n_1[W_S^1+W^2]\)\(-2n_2[W^0+W_S^2]\)\(+2n_3[W_S^1+W^0]\)\(+2n_4[W^2+W_S^2]\)

4) \(2\dfrac{1}{4}W_S^1[(-T-I_z-S_z-2S_zI_z)+(T+I_z-S_z-2I_zS_z)]=W_S^1(-S_z-2I_zS_z)\)

\(2\dfrac{1}{4}W^2[(-T-I_z-S_z-2S_zI_z)+(T-I_z-S_z+2I_zS_z)]=W^2(I_z+S_z)\)

\(2\dfrac{1}{4}W^0[(-T+I_z-S_z+2S_zI_z)+(T+I_z-S_z-2I_zS_z)]=W^0(I_z-S_z)\)

\(2\dfrac{1}{4}W_S^21[(-T-I_z+S_z+2S_zI_z)+(T-I_z-S_z+2I_zS_z)]=W_S^2(-S_z-2I_zS_z)\)

**\(\dfrac{dS_z}{dt}=-S_z(W_S^1+W^2+W^0+W_S^2)+2I_zS_z(W_S^2-W_S^1)+I_z(W^0-W^2)\)**

#### Deriving 2I_{z}S_{z}

\[2I_zS_z=n_1-n_2-n_3+n_4.\]

1)\(\dfrac{d2I_zS_z}{dt}=\dfrac{n_1}{dt}-\dfrac{n_2}{dt}-\dfrac{n_3}{dt}+\dfrac{n_4}{dt}\)

2)\(\dfrac{d2I_zS_z}{dt}=[-W_S^1n_1-W_I^1n_1-W^2n_1+W_S^1n_3+W_I^1n_2+W^2n_4]\)

\(-[W_I^1n_1+W_S^2n_4+W^0n_3-W_S^2n_2-W_I^1n_2-W^0n_2]\)

\(-[W_S^1n_1+W_I^2n_4+W^0n_2-W_I^2n_3-W_S^1n_3-W^0n_3]\)

\(+[-W_S^2n_4-W_I^2n_4-W^2n_4+W_I^2n_3+W_S^2n_2+W^2n_1]\)

3)\(\dfrac{d2I_zS_z}{dt}=-2n_1(W_S^1+W_I^1)+2n_2(W_I^1+W_S^2)+2n_3(W_S^1+W_I^2)-2n_4(W_S^2+W_I^2)\).

4)\(W_S^1[-2\dfrac{1}{4}(T+I_z+S_z+2I_zS_z)+2\dfrac{1}{4}(T+I_z-S_z-2I_zS_z)]=-W_S^1[S_z+2I_zS_z]\)

\(W_I^1[-2\dfrac{1}{4}(T+I_z+S_z+2I_zS_z)+2\dfrac{1}{4}(T-I_z+S_z-2I_zS_z)]=W_I^1[I_z+2I_zS_z]\)

\(W_S^2[-2\dfrac{1}{4}(T-I_z+S_z-2I_zS_z)-2\dfrac{1}{4}(T-I_z-S_z+2I_zS_z)]=-W_S^2[S_z-2I_zS_z]\)

\(W_I^2[-2\dfrac{1}{4}(T+I_z-S_z-2I_zS_z)-2\dfrac{1}{4}(T-I_z-S_z+2I_zS_z)]=W_I^2[I_z-2I_zS_z]\)

**\(\dfrac{d2I_zS_z}{dt}=S_z(W_S^2-W_S^1)+I_z(W_I^2-W_I^1)-2I_zS_z(W_S^2+W_S^1+W_I^2+W_I^1)\)**

## Relaxation Mechanisms of the Solomon equations

The Solomon equations for operators I_{z}, S_{z}, and 2I_{z}S_{z} are given below (in case the reader skips the derivation)

\(\dfrac{I_z}{dt}=-(W_I^1+W_I^2+W^2+W^0)I_z-(W^2-W^0)S_z-(W_I^1-W_I^2)2I_zS_z\)

\(\dfrac{S_z}{dt}=-S_z(W_S^1+W^2+W^0+W_S^2)+2I_zS_z(W_S^2-W_S^1)+I_z(W^0-W^2)\)

\(\dfrac{2I_zS_z}{dt}=S_z(W_S^2-W_S^1)+I_z(W_I^2-W_I^1)-2I_zS_z(W_S^2+W_S^1+W_I^2+W_I^1)\)

Now, these equations are not quite right, as the treatment that we gave actually considers only the perturbed spins. Therefore we may now re-write these equations to account for the the spins that are perturbed by defining the perturbed spins as the perturbed population+equilibrium population=I_{z,} while the equilibrium population is I^{0}_{z}. The equations may now be re-written as

\(\dfrac{I_z-I_z^0}{dt}=-(W_I^1+W_I^2+W^2+W^0)(I_z-I_z^0)-(W^2-W^0)(S_z-S_z^0)-(W_I^1-W_I^2)\)

\(\dfrac{S_z-S_z^0}{dt}=-(S_z-S_z^0)(W_S^1+W^2+W^0+W_S^2)+2I_zS_z(W_S^2-W_S^1)+(I_z-I_z^0)(W^0-W^2)\)

\(\dfrac{2I_zS_z}{dt}=(S_z-S_z^0)(W_S^2-W_S^1)+(I_z-I_z^0)(W_I^2-W_I^1)2I_zS_z(W_S^2+W_S^1+W_I^2+W_I^1)\)

### Self Relaxation

The Solomon equations for a given spin, here described for I_{z} but applicable for S_{z} and 2I_{z}S_{z}, show that the rate of relaxation is tied to I_{z}, S_{z}, and 2I_{z}S_{z}. The rate constants associated with I_{z} corresponds to the self relaxation of the spin magnetization on its own. This is given by the overall rate

\(W_I^1+W_I^2+W^2+W^0\)

Note that the rates W^{2} and W^{0} correspond to processes cross-relaxation processes that need a spin S to occur. Therefore, the relaxation rate for I_{z} is not dependent on any other spins in the system. The summation of the rate constants is denoted by R_{I}.

### Cross-Relaxation

The term (W^{0}-W^{2})S_{z} is present in the Solomon equations for I_{z} and corresponds to the the rate at which S_{z} transfers magnetization to I_{z }spin magnetization in a phenomenon is known as cross-relaxation. To properly illustrate the effect, take \(W^2=W^0=0\). In this case, relaxation of a spin I is independent of spin S. On the other hand, if \(W^2 and/or W^0\neq0\) but spin S is not perturbed, (i.e, (\(S_z-S_z^0)=0\)) cross relaxation will not occur either. The cross relaxation term is denoted as \(\sigma_{IS}\).

### 2I_{z}S_{z} Relaxation

The relaxation associated with \(W_I^1-W_I^2\) is the transfer of magnetization from I_{z}S_{z} to I, denoted as \(\Delta_{I}\), which will only occur if the rate constants are different. Meanwhile the operatore 2IzSz also shows self-relaxation as RIS, denoted by

\(R_{IS}=W_I^1+W_I^2+W_S^1+W_S^2\)

## References

1) I.Solomon. "Relaxation Processing in a System of Two Spins" *Phys. Rev. *99, **1955**, 559.

2) J. Keeler. "Relaxation" Chapter 8 **2004**

3) M. Levitt. "Spin Dyanmics"