The Lipari-Szabo Model-Free Analysis as a Method for Study of Molecular Motion in Solid State Heteronuclear Systems Using NMR Off-Resonance.

In the present work, a new method for measuring motional parameters using the off-resonance technique was described. The Lipari-Szabo model-free formalism was used to analyze molecular dynamics in a heteronuclear system [1, 2]. Cross-relaxation solid state nuclear magnetic resonance off-resonance experiments were performed on a homebuilt pulse spectrometer operating at the frequency of 30.2 MHz for protons at temperature 173 K. The proton spins were spin-locked in the effective field [Formula: see text] while 19F spins were continuously saturated for a long time. It was possible to carry out these experiments because a uniquely designed probe was able to produce three slightly differing frequencies on and off-resonance for protons and the frequency of 28.411 MHz for fluorine [3-6]. The off-resonance frequencies can be changed within the range of 30.2-30.6 MHz.

Introduction

Solid state nuclear magnetic resonance (NMR) off-resonance measurements provide valuable information about molecular dynamics in a heteronuclear system [7-10]. So far only a few authors have presented the results of study of the cross-relaxation in the heteronuclear systems in the solid state [7, 10–12]. Complex analysis of these systems and apparatus requirements pose problems which are difficult to overcome. The construction of an efficient probe with high sensitivity, working at three different but close to each other frequencies, gives a possibility to obtain very interesting data on molecular dynamics in the heteronuclear systems [3–6, 10]. A uniquely designed and constructed probe of NMR spectrometer contains two independently tuned and electrically isolated coils installed in the coaxial position on the dewar. Three slightly differing frequencies on and off-resonance for protons and the frequency of 28.411 MHz for fluorine nuclei with no electrical interaction with one another can be produced.

This work puts forward an innovative and useful NMR off-resonance technique, in which detected proton 1H spins are spin-locked in the effective field while fluorine 19F spins are continuously saturated for a long time. The saturation of fluorine spins for a long time simplifies the procedure for the analysis of experimental data, on account of the suppression of the heteronuclear 1H-19F dipolar coupling. For the analysis of NMR off-resonance data a model-free approach is applied. It is based on a rigorous theory advanced by Lipari and Szabo [1, 2] and is advantageous because it allows us to directly receive valuable parameters with well-defined physical significance [13-17]. According to this approach, molecular motion can be specified by set molecular parameters: a generalized order parameter S2and an effective correlation time τ. The order parameter S2 measures the spatial restriction of the internuclear vector motions and describes the amplitude of molecular motion. Its values are in the range from 0 to 1. When its value is equal to 0, it corresponds to isotropic motion; the value 1 of parameter S means that the motion is completely restricted. The effective correlation time τ is a measure of the rate of motion and it is linked with a single correlation time τ, describing the overall motion, by the equation:where τ is the internal correlation time.

The object of this study is a solid state sample of guanidinium tetrafluoroborate [C(NH2)3]+BF4- (GHBF4) [18]. The work focuses on research methodology. The sample is used as a model in order to demonstrate that the test results agree with the previous results [10]. Experimental data was obtained using cross-relaxation solid state NMR off-resonance technique and were analyzed using the Lipari-Szabo model-free formalism. From the dependence of the ratio of the cross-relaxation rate σ to the proton relaxation rate ρ on the effective frequency ω, three motion parameters: the order parameter S2, the correlation time τM and the internal correlation time τi were determined. It is a new application of the Lipari-Szabo model-free analysis to the interpretation of NMR off-resonance data on molecular dynamics in heteronuclear systems.

Theory

Let us consider a system consisting of two spins I = 1/2 and S = 1/2 with magnetogyric ratios γI and γS, respectively. In our case spin I refers to protons and S refers to fluorines. These two Zeeman reservoirs are connected to the lattice and to each other. Heteronuclear dipolar interaction mediates cross-relaxation between these Zeeman reservoirs and is characterized by the cross-relaxation rate σ. Both heteronuclear and homonuclear dipolar interactions mediate the spin–lattice relaxation process of each spin system separately with the proton and the fluorine spin–lattice relaxation rates ρI and ρS, respectively. The equations describing the time evolution of the longitudinal spin polarizations: 〈I〉(t) and 〈S〉(t) in the laboratory frame were led out by Solomon [19]. The macroscopic differential equation for the longitudinal rotating frame relaxation of spin I modulated by dipolar interaction with spin S was derived by Peng and co authors [8]. The decay of the polarizations 〈I’〉(t) along the effective field is given by the following equation [8, 10]:where ρ is the off-resonance spin lattice relaxation rate in the rotating frame, σ is the cross-relaxation rate, I’(∞) is the equilibrium polarization for the I spins in the rotating frame, S(∞) is the polarization of the S spins at thermal equilibrium in the laboratory frame projected on the direction of the effective field ; N and N describe the number of nuclei I and S, respectively.

When the S spins are saturated for a long time, the equation 〈S〉(t) = 0 is satisfied and Eq. (2) can be rewritten as:The analytical solution to Eq. (3) is monoexponential [9, 10] and is presented by the following equation:where 〈I’〉(0) is the polarization of spins I at time t = 0 dependent on the initial conditions of the experiment.

Taking into account the following dependences:

Equation (4) can be rewritten to form:In the NMR off-resonance experiment in the rotating frame, the effective magnetic field is given by the sum of the vector , parallel to the static field , and perpendicular vector , which represents the amplitude of the radio-frequency (rf) field. The component has the amplitude proportional to the difference between the resonance angular frequency ω0 and the angular frequency ω of field B. The effective field forms an angle β with the static magnetic field and is defined by the equation:In the NMR off-resonance experiment in the rotating frame with the saturating of fluorines (Fig. 1) after the completion of the off-resonance rf pulse, the parallel and the perpendicular components of magnetization M’(t), recover to the thermodynamic equilibrium in the static magnetic field B0. After the time t (about 200 μs for solid state), much longer than spin–spin relaxation time T, the perpendicular component of magnetization disappears and then the second 90° rf on-resonance pulse or a “solid echo” sequence on resonance is applied.Taking into account the following equations [10]:and introducing these expressions into Eq. (6) we get:In Eq. (9) the magnetization M’(0) of spins I at time t = 0 depends on the initial conditions of the experiment. According to the initial conditions of the experiment for t = 0 the following equation is satisfied:and Eq. (9) assumes the form:The monoexponential Eq. (11) was used to fit the experimental data, where: N and N are the numbers of proton and fluorine spins, ρ is the proton spin–lattice relaxation rate in the rotating frame, σ is the cross-relaxation rate in the rotating frame.

Fig. 1

Pulse sequence for NMR off-resonance experiment in the rotating frame with the saturating of fluorines

By measuring the amplitude of NMR signal M(t) as a function of duration t with the saturating of fluorines we are able to obtain the proton relaxation rate ρI and the ratio of the cross-relaxation rate σ to the proton relaxation rate by fitting formula (11) to the experimental data.

The ratio contains information about nuclear Overhauser effect in the off-resonance rotating frame NOE and can be expressed by a suitable spectral density function:The NMR off-resonance data were analyzed using the Lipari-Szabo model-free approach, which assumes that the dipolar spectral density can be written as:whereS2 is the order parameter, which is a measure of the restriction of the motion,

τ is the long rotational correlation time describing the global rotation of the molecule,

τ is the effective correlation time,

τ is the short internal correlation time, describing the fast reorientation of the local group.

Introducing Eq. (13) into Eq. (12) and assuming that:after a little algebra one gets:whereEquation (15) takes three independent parameters:From the ratio:the correlation time τ can be determined.

And then from Eq. (19) we are able to estimated the value of the order parameter S.

Experiment and results

The NMR off-resonance experiment in the rotating frame with the saturating of fluorines was performed on a solid state sample of guanidinium tetrafluoroborate [C(NH2)3]+BF4− (GHBF4), synthesized in the Division of Molecular Crystals [18] at temperature 173 K on a homebuilt pulse spectrometer with a 10 % error [6]. A specially designed and constructed probe of NMR spectrometer contains two independently tuned and electrically isolated coils installed in the coaxial position on the dewar [3–6, 10]. Figure 1 presents a pulse sequence used in NMR off-resonance experiment. Two rf pulses are applied at the same time: one at the frequency off-resonance for the proton spins I, which produces the effective field and the second rf pulse at the frequency of 28.411 MHz, which saturates the fluorine spin S. The fluorines are saturated throughout the whole duration of the off-resonance rf pulse. After the pulse sequence at the resonance frequency for protons of 30.2 MHz, the solid echo is obtained. The echo amplitude is proportional to the proton magnetization M(t) and depends on the duration t of the the off-resonance rf pulse. By fitting Eq (11) to the experimental data, the proton relaxation rate ρ and the ratio of the cross-relaxation rate σ to the proton relaxation rate ρ were estimated. Measurements were performed as a function of the intensity of the effective field at the constant angle β equal to 10° by simultaneous changes in both the amplitude of rf field B and the angular frequency ω of field B. The dependence of σ/ρ ratio as a function of the effective frequency ω in the rotating frame is presented in Fig. 2. The model-free approach was used to analyze experimental results. Equation (15) was used to fit experimental data. Parameters a, b and c, described by Eq. (19) were obtained by the method of least squares. Then the correlation time τ was estimated based on Eq. (20) and subsequently the order parameter S.

Fig. 2

Dependence of σ /ρ ratio as a function of the effective frequency ω in the rotating frame

The calculated internal correlation time τ is equal to 0.44 × 10−8 s, the long rotational correlation time describing the global rotation of the molecule τ = 1.04 × 10−6 s and S is equal to 0.4.

Conclusion

The cross-relaxation NMR off-resonance experiment and the application of the Lipari-Szabo spectral density function make it possible to extract numerical values of the dimensionless order parameter S2 and of two correlation times, which describe two molecular processes: one with the overall correlation times τ and the other with the internal correlation time τ.

Using the model-free approach allows us to study the molecular dynamics without the assumption of a specific model for internal motion.

In general, the measurements can be performed (1) at a constant angle β and variable effective frequency ωe, which enables us to separate relaxation mechanisms, (2) at a constant effective frequency ωe and variable angle β.