The effect of the applied RF field in an NMR experiment on the magnitude of the Spectral Density for a Dipolar Relaxation Mechanism is demonstrated theoretically. The effect is shown with Sin Cos Pulse as a concrete example. The order of magnitude of the magnetic moment where these effects will be significant for typical Rf amplitude values is derived. The effect may be of utility in providing an alternate method of control for MRI Tissue Contrast applications with further development.
In contemporary NMR methodologies, it is common to find experimental scenarios where the relaxation of the magnetization during a pulse train is important to be able to model and quantify [1-3]. In this note we suggest that for some molecular species the Rotational Diffusion can be affected and modified by the Magnetic Field Torque of the applied radio-frequency pulse. During the course of working on this concept, it has come to our attention that the Russian investigator Sitnitsky [4] has investigated this phenomenon.
This proposed influence may be important in some models for explaining experimental data, such as for Liquid Crystals [5]. We demonstrate the derivation of this effect on the spectral densities following the classic treatment of Abragam [6] and gives some ranges of parameters where this effect may be of importance.
We note that the proposed effects may be useful as another avenue to control the spin dynamics of an experimental system while the pulse is on. Also, the proposed effects have been dealt with rigorously in the Physics Literature [7].
The "Toy Model" we propose to explicate this effect is the following.
We envision a spin system, transformed to the so called Tilted Doubly Rotating Frame (TDRF [8]). In this frame there will be defined a so-called "effective field." We can write down an effective Hamiltonian for the applied RF of the following form:
$${H}^{RF}[t]={I}_{x}{\omega}_{1}[t]+{I}_{z}\Delta \omega [t]\text{(1)}$$
For the exposition here we consider the Sin/Cos pulse, defined as:
${\omega}_{1}[t]={\omega}_{1}^{M}\mathrm{Sin}[{\omega}_{1}^{M}t]\text{[2a]}$
$\Delta \omega [t]={\omega}_{1}^{M}\mathrm{Cos}[{\omega}_{1}^{M}t]\text{[2b]}$
Where ${\omega}_{1}^{M}$ is a constant (See for example the relevant papers of the Garwood Group [9-11].
In the TDRF the effective field can be seen from geometric arguments to be defined as:
${\omega}_{eff}[t]=\sqrt{{\omega}_{{}_{1}}{[t]}^{2}+\Delta \omega {[t]}^{2}}\text{[3]}$
Substituting Eq [2 a,b] into Eq[3], one easily appreciates that ${\omega}_{eff}[t]={\omega}_{1}^{M}$. So that as required for our argument the effective field defined as:
${B}_{eff}[t]=\frac{{\omega}_{1}^{M}}{{\gamma}_{particle}}\text{[4]}$
Where ${\gamma}_{particle}$ is the particle gyromagnetic ratio.
Now we consider a Molecular Species in solution with a defined dipole moment ${\mu}^{particle}$.
In a constant field ${B}_{eff}$, there is a potential energy of interaction [12] between the moment and the field defined as:
$U[\theta ]=-{\mu}^{particle}{B}_{eff}\mathrm{Cos}[\theta ]\text{[5]}$
Here the angle $$\theta $$ is defined as the angle between the vectoral directions of the dipole moment and the constant field.
We change to the convenient notation:
${K}_{0}={\mu}^{particle}{B}_{eff}\text{[6a]}$
So that:
$U[\vartheta ]=-{K}_{0}\mathrm{Cos}[\vartheta ]\text{[6b]}$
Knowing the geometry between the effective field and the magnetic moment in the TDRF, it is Seen that the angle $\vartheta $ is defined as:
$\vartheta [t]=Arc\mathrm{Tan}[\frac{{\omega}_{1}[t]}{\Delta \omega [t]}]\text{[7]}$
Using Eq [2a,b] in Eq [7] we see that:
$\vartheta [t]={\omega}_{1}^{M}t\text{[8]}$
Suppose we take the Nuclear Species of Interest to be in a molecule that we model and approximate as a sphere. We assume that the Rotational Brownian Motion can be represented as a series of small incremental rotations. We seek to find the Correlation Function which characterizes the rotational diffusion. As treated, in for example Abragam Chapter VIII or other places in the literature [13,14] we can define the Correlation Function in terms of the spatial part of the Dipolar Interaction Hamiltonian. If we adopt the notation of Abragam, we can define the Correlation Function as:
$$G[t,{K}_{u},{m}_{p}]={\displaystyle {\iint}_{{\Gamma}_{\vartheta}{\Gamma}_{{\vartheta}_{p}}}{F}^{\ast}[\theta ,}0,{m}_{p}]F[{\theta}_{p},0,{m}_{p}]W[\theta ,{\theta}_{p},Ku,t]\mathrm{Sin}[\theta ]d\theta \mathrm{Sin}[{\theta}_{p}]d{\theta}_{p}\text{[9]}$$
Where we set and consider the case where ${m}_{p}$ is zero.
So, to carry out this program we need to an expression for the Probability Density Function.
This PDF will be a solution of the so-called Smoluchowski Equation (SE), where the effects of the applied RF Torque will be included. As one can infer there are numerous assumptions one can apply to the formulation of the SE. The solution in general, (see for example the classic papers of Coffey's group [12] are not totally trivial, usually the derivation of series solutions which involve the solution of iterative expressions for the expansion coefficients, or continued fraction solutions.
We have chosen to present and use the solution of Sitnitsky [4] which is the most easily implemented solution we have found to program for demonstration of our methods.
Please see Appendix I for a detailed definition of the terms in the series expression for the PDF.
The PDF can be taken to be the approximate solution of the following partial differential Equation using our expression for the Potential Energy Term [4].
${\tau}_{R}\frac{\partial w[x,{x}_{i,}{K}_{u},t]}{\partial t}=\frac{{\partial}^{2}w[x,{x}_{i},{K}_{u},t]}{{\partial}^{2}x}-(2{K}_{u}+{\tau}_{R}C[x])\frac{\partial w[x,{x}_{i,}{K}_{u},t]}{\partial x}+w[x,{x}_{i},{K}_{u},t]({K}_{u}{}^{2}+{\tau}_{R}C[x]{K}_{u})\text{[10a]}$
$C[x]=\frac{b}{\mathrm{Cos}[{k}_{1}x]}{\displaystyle \underset{-1}{\overset{1}{\int}}d{x}^{\prime}}\mathrm{Cos}[{k}_{1}{x}^{\prime}](x-{x}^{\prime})Exp[-\lambda |x-{x}^{\prime}|]\text{[10b]}$
Here b is a constant defined in (4).
In Eq [10] x is defined to be as $\mathrm{Cos}[\theta ]$ where ${K}_{u}$ measure as defined previously the interaction between the Moment and the RF field with the definition:
${K}_{u}=\frac{{K}_{0}}{{k}_{B}T}\text{[11]}$
We can use the definition of the PDF to compute the Correlation Function as In Eq [9] and then compute the corresponding Spectral Density as:
$J[\omega ,{K}_{u},t,{m}_{p}]={\displaystyle \underset{0}{\overset{t}{\int}}G[{t}_{p},}{K}_{u},{m}_{p}]Exp[i\omega {t}_{p}]d{t}_{p}\text{[12]}$
In Figure 1 we show the dependence of the Spectral Density as given in Eq [12], for the case mp = 0.
Figure 1: Dependence of Spectral Density on the Ku parameter. (See text for definition of Ku). View Figure 1
Figure 2: Plot of Probability Density Function Theta. View Figure 2
As can be seen there is found to be an appreciable dependence of the Spectral Density on the parameter ${K}_{u}$ This is taken to indicate that the RF Field, with a range of values which will be discussed below, can affect the Spectral Density which is used to compute relaxation functions [6,13-15]. So that the RF field, through interaction on the Rotational Brownian Motion, can influence the values of the calculated relaxation functions during a pulse sequence.
To the knowledge of the author, this possibility has not been fully appreciated in the NMR literature.
The reader may wonder what is a lower bound on the magnetic moment of the particle of Interest for a typical value of the pulse amplitude.
In the Garwood papers [9-11], the pulse amplitude is typically on the order of Hz. Then we reason that the interaction energy of the magnetic moment with the field in the TDRF should be greater than the thermal energy of the surrounding liquid medium.
So, we propose:
${K}_{u}\gg 1$
${u}^{particle}{B}_{eff}\gg {k}_{B}T$
Or
${u}^{particle}\gg \frac{{k}_{B}T}{{B}_{eff}}$
We note that the units of a magnetic moment can be seen in CGS units to be $\frac{ergs}{Gauss}$.
So, at room temperature in CGS units, ${k}_{B}T$ is on the order of ${10}^{-14}ergs$.
${B}_{eff}$, is typically on the order of ${10}^{-4}$ Gauss.
${u}^{particle}\gg {10}^{-10}\frac{ergs}{Gauss}$
Practical lower bound on the magnetic moment of the particle for an effect of the RF field on the Rotational Brownian Motion of the particle and consequently on the Spectral Density for a dipolar relaxation mechanism.