Effects of radial radio-frequency field inhomogeneity on MAS solid-state NMR experiments
Radio-frequency field inhomogeneity is one of the most common imperfections in NMR experiments. They can lead to imperfect flip angles of applied radio-frequency (rf) pulses or to a mismatch of resonance conditions, resulting in artefacts or degraded performance of experiments. In solid-state NMR under magic angle spinning (MAS), the radial component becomes time-dependent because the rf irradiation amplitude and phase is modulated with integer multiples of the spinning frequency. We analyse the influence of such time-dependent MAS-modulated rf fields on the performance of some commonly used building blocks of solid-state NMR experiments. This analysis is based on analytical Floquet calculations and numerical simulations, taking into account the time dependence of the rf field. We find that, compared to the static part of the rf field inhomogeneity, such time-dependent modulations play a very minor role in the performance degradation of the investigated typical solid-state NMR experiments.
Radio-frequency (rf) field inhomogeneity describes the spatial inhomogeneity of the rf field inside the coil or sample volume and is one of the major experimental imperfections that leads to artefacts or reduced efficiency in NMR experiments. The magnitude of the rf field amplitude distribution over the sample space can be estimated with a nutation experiment (Torrey, 1949; Barnaal and Lowe, 1963). Measuring such nutation spectra of thin sample slices placed along the rotor axis allows the characterization of the spatial rf field distribution along the coil axis (Nishimura et al., 2001; Paulson et al., 2004). The full spatial distribution, however, is only accessible using gradient methods (Guenneugues et al., 1999; Odedra and Wimperis, 2013) that are typically not available in solid-state NMR probes. Alternative approaches include the measurement of the rf field amplitude using the ball shift experiment (Maier and Slater, 1952), numerical simulations based on finite elements, or approximative analytical solutions of the Maxwell equations (Engelke, 2002; Tošner et al., 2017, 2018). The design of the coil geometry has a major influence on the magnitude and distribution of the rf field amplitude over the active sample volume, and different geometries have been proposed to improve the rf homogeneity (Idziak and Haeberlen, 1982; Privalov et al., 1996; Li et al., 2006). However, in solid-state NMR probes, solenoid coils along the sample spinning axis are most commonly used due to the high, achievable rf field amplitudes. The gap between the rotor and the coil is minimized in order to optimize the filling factor. This design choice typically leads to large rf field inhomogeneity that can manifest itself in reduced efficiency in experiments such as cross-polarization (Hartmann and Hahn, 1962; Stejskal et al., 1977), homonuclear decoupling (Bielecki et al., 1989, 1990; Mote et al., 2016), heteronuclear decoupling (Purusottam et al., 2015; Frantsuzov et al., 2017), symmetry-based recoupling sequences (Levitt, 2007), or even pulsed recoupling experiments like rotational-echo double-resonance (REDOR; Nishimura et al., 2001).
Reducing the magnitude of the rf field inhomogeneity can be achieved experimentally by physically restricting the sample along the rotor axis or even to a sphere in the centre of the rotor (Lindon et al., 2009). Alternatively, gradients can be used for a virtual sample restriction (Charmont et al., 2000), but since gradients are not very common in solid-state NMR probes, this approach is rarely used. Another possibility is radio or nutation-frequency-selective pulses that can be used for the same purpose (Charmont et al., 2002; Aebischer et al., 2020). All these methods, however, are accompanied by a reduction in signal due to the restriction of the measured sample volume to a smaller part of the coil volume.
In solid-state NMR under magic angle spinning (MAS) conditions, the radial component of the rf field is modulated by time (Levitt et al., 1988; Tekely and Goldman, 2001; Goldman and Tekely, 2001; Tošner et al., 2017), leading to further potential complications in the experiments. Such MAS-induced time-dependent radio-frequency fields could give rise to additional or modified resonance conditions or to other changes in the effective Hamiltonian generated by the pulse sequence. The importance of such time-dependent terms was first described in rotary resonance recoupling (Levitt et al., 1988), where it leads to changes in the observed line shape. Besides the appearance of additional sidebands in cross-polarization experiments (Tekely and Goldman, 2001; Goldman and Tekely, 2001), nutation spectra (Elbayed et al., 2005), and reported phase distortions and loss of magnetization in MLEV-16 sequences under MAS (Piotto et al., 2001), there have been very few studies of the effects of such modulations of the amplitude and phase of the rf field caused by MAS rotation. These modulations have been included in the design of a heteronuclear polarization transfer scheme based on optimal control strategies (Tošner et al., 2018), where impressive gains have been shown. These improvements prompted us to investigate potential effects of such MAS-modulated radio-frequency field amplitudes and phases on basic building blocks in common solid-state NMR pulse sequences in more detail. The approach we have chosen is rather simple. We use analytical approaches based on Floquet theory (Leskes et al., 2010; Scholz et al., 2010) and numerical simulations based on computed rf field distributions in a typical MAS rotor to characterize the time evolution of the density operator under MAS rotation with and without time-dependent rf field amplitudes and phases. The computational approach allows us to investigate the time evolution of the density operator in different spatial parts of the rotor and to magnify the amplitude or phase modulations of the rf field to obtain a better picture of their importance.
To illustrate the magnitude and distribution of the rf field amplitude and phase over the active sample volume in some typical MAS NMR probes, rf field distributions were calculated based on Engelke (2002) and Tošner et al. (2017). Figure 1 shows the relative amplitude and phase of the rf field in a cylindrical coordinate system as a function of z (the axis along the rotor axis) and r (the radial direction) for three common probe designs with MAS rotors of 3.2, 1.9, and 1.3 mm outer diameter. In these plots, the angle was chosen. The maximum intensity of the rf amplitude distribution was used as the reference point for the relative amplitude ωrel(r); hence, a value of 1 means that the amplitude experienced at this position corresponds to the nominal rf amplitude. The rf phase was computed relative to the centre point of the rotor with . One can clearly see the decay in the rf amplitude towards the edges of the rotor along the rotor axis (large z values), while phase errors mostly occur for large r and z values.
The radial dependence of the relative rf field amplitude and phase as a function of the angle ϑ is shown in Fig. 2 for different values of z and r for the 3.2 mm MAS probe. Under sample rotation, the angle ϑ varies as a function of time, and the rf field amplitude and phase are periodically modulated with the rotor frequency. The trajectories in Fig. 2 clearly show that the magnitude of these amplitude and phase modulations increases towards the edges of the rotor. We, therefore, expect crystallites located at large r and z values to experience the strongest modulations of the rf field amplitude and phase.
For numerical simulations of spin dynamics, the numerical values for ω1,rel(r) and ϕrel(r), obtained from simulations of the rf field distribution shown in Figs. 1 and 2, were used directly as input. For analytical calculations based on Floquet theory, a parameterization of the values using a Fourier series with the MAS frequency, ωr, is more convenient and was obtained by fitting the following expressions:
to the amplitude and phase changes of the rf field (see Fig. 2). Typically, terms of the Fourier series up to n=4 were used in the fits to characterize the time-dependent amplitude ωrel(t) and phase ϕrel(t).
3.1 Floquet description
In the high field approximation, the total Hamiltonian in the rotating frame under MAS for a homonuclear spin system comprised of NI spins is given by the following:
Fourier components of the spatial tensors ω(n) of the chemical shift of a spin Ip and couplings between two spins Ip and Iq are given by the following:
for the isotropic and anisotropic chemical shifts, the scalar J, and the anisotropic dipolar coupling. The sets of Euler angles describe the orientation of the tensors in the rotor-fixed frame, and denote the reduced Wigner matrix elements. The elements of the tensor in its principal axis system are denoted by ρℓ,m. The most general rf Hamiltonian for Eq. (3) is given by the following:
where both the rf amplitude ω1(t) and the rf phase ϕ(t) can be time-dependent due to the irradiation scheme. Under MAS, the radial part of the rf inhomogeneity will lead to additional modulations of ω1(t) and ϕ(t) that are periodic with the rotor frequency.
The spin-system Hamiltonian can be transformed into an interaction frame with respect to by the following:
The propagator characterizing the interaction frame transformation is given by the following:
where is the Dyson time-ordering operator (Dyson, 1949) that ensures proper time-ordering of non-commuting operators in products. If the rf irradiation is periodic with τm, the spin operators in Eq. (3) will transform according to the following:
where the transformation behaviour of the Cartesian spin operator is characterized by the Fourier coefficients (Scholz et al., 2009, 2010). These coefficients only depend on the rf irradiation scheme and are independent of the spin system details. Generally, the interaction frame trajectories of the spin operators are characterized by two basic frequencies, i.e. the modulation frequency of the pulse sequence and an additional effective nutation frequency that can be determined from the overall flip angle over one period of the pulse scheme . This additional effective field is zero if the propagator over a full cycle of the pulse sequence is unity (Tan et al., 2016).
The interaction frame Hamiltonian can thus be expanded as a Fourier series with three basic frequencies, as follows:
with Fourier components as follows:
where denotes the Kronecker delta. For convenience, the general two-spin Fourier coefficients were defined. They can be computed as the convolution of single-spin coefficients as follows:
The scalar product of the I spin vector operators remains time invariant under rf irradiation and can be incorporated into the Fourier coefficients. In triple-mode Floquet theory, the first-order effective Hamiltonian is given by all contributions that satisfy the resonance condition, as follows:
and thus, the sum of non-resonant () and resonant terms is as follows:
Analogously, in the following the second-order effective Hamiltonian is given by:
The summation is restricted to values of ν, κ, and λ for which is satisfied in order to avoid singularities.
3.1.1 Theoretical description including rf inhomogeneity
For spatial rf field distributions that do not have cylindrical rotation symmetry, MAS will lead to a periodic modulation of the rf field amplitude and phase experienced by a spin packet. At a given position, the general rf Hamiltonian, including these additional modulations, can be expressed as follows:
where ω1,nom(t), and ϕnom(t) correspond to the nominal rf amplitude and phase, i.e. to the values corresponding to a perfectly homogeneous rf field. They are determined by the pulse scheme under investigation and are periodic with . Deviations from these nominal values are introduced by the relative rf amplitude and phase ω1,rel(t) and ϕrel(t) that are periodic with . The overall period of the rf Hamiltonian will only be of finite length if the modulation frequency of the pulse scheme ωm and the rotor frequency ωr are commensurate:
where c corresponds to the number of rotor cycles required for the synchronization of the rf irradiation and the MAS rotation. Such a synchronization condition is generally fulfilled for rotor-synchronized pulse schemes such as most recoupling sequences (Nielsen et al., 2012). For irradiation schemes that are typically applied asynchronously to avoid resonance conditions, a careful selection of the synchronization condition is required. In principle, the treatment can be generalized to all cases where the largest common divisor of ωr and ωm is not too small. Such a synchronization is required to make the modulations of the rf field Hamiltonian by MAS cyclic over the basic repetition time of the sequence. The interaction frame Hamiltonian can then be written as follows:
where the substitution was used. The summation over the index k runs from −∞ to +∞; the sum over n′ in Eq. (23) is, thus, also unrestricted. The resulting interaction frame Hamiltonian is modulated with due to the time dependence of the spatial part of the Hamiltonian during MAS (Fourier number n). The radial rf inhomogeneity leads to an additional modulation of the spin part with the rotor frequency. As the modulation frequency of the pulse sequence is commensurate with ωr, these two modulations can be combined and a single Fourier number n′ can be used. Therefore, the Hamiltonian is given by a Fourier series with only two basic frequencies (ωm and ωeff), and the triple-mode Floquet analysis is reduced to a bimodal treatment. Alternatively, one can continue with the triple-mode Floquet description and assume resonance conditions between ωm and ωr. For resonant phenomena, the triple-mode Floquet approach is more suited since the additional effective field can lead to changes in the resonance conditions, leading to asynchronous sequences (Hellwagner et al., 2017; Tan et al., 2015). For non-resonant phenomena, both descriptions will give equivalent results. The periodic rf amplitude and phase modulations due to the radial inhomogeneity merely affect the interaction frame transformation, which is fully characterized by ωm and ωeff, as long as these rf field modulations are periodic on the length of the interaction frame trajectory.
Possible resonance conditions for the bimodal interaction frame Hamiltonian of Eq. (23) are given by the following:
As effective fields due to the rf inhomogeneity will typically be small compared to ωr, and ℓ is limited to a maximum value of two, it is reasonable to assume that only non-resonant terms will contribute to the effective Hamiltonian. The second-order approximation of can, thus, be written as follows:
where the summation over ν and λ is restricted to values satisfying .
Each of the Fourier components in Eq. (23) is composed of a sum of several terms since there are multiple combinations of n and k resulting in the same n′. As the index n is limited to values between ±2 (limited by the rank of the spatial tensor), the are given by the following:
In the first-order approximation, the non-resonant contribution to the effective Hamiltonian is simply given by the Fourier component as follows:
The second-order scaling factors , , and for χ=x are given by the following:
Similar expressions result for χ=y and z and can be found in the Supplement (Sect. S1).
All first- and second-order scaling factors can be obtained from the Fourier coefficients characterizing the rf interaction frame trajectory of the Cartesian spin operator and, thus, do not depend on the details of the spin system. The effects of the additional rf field modulations due to the radial contribution to the rf inhomogeneity will lead to changes in the interaction frame trajectories and, therefore, changes in the scaling factors for the effective Hamiltonians compared to those calculated assuming a perfectly homogeneous rf field.
4.1 Numerical simulations
The effect of the rf field inhomogeneity on common solid-state NMR pulse sequences was investigated by numerical simulations in the usual rotating frame using the GAMMA spin-simulation environment (Smith et al., 1994). Unless otherwise noted, spin dynamics were simulated at a B0 field of 14.1 T, corresponding to a proton resonance frequency of 600 MHz. Powder averaging was implemented according to the Zaremba–Conroy–Wolfsberg (ZCW; Cheng et al., 1973) scheme, using between 100 and 10 000 crystallite orientations. A summary of the simulation parameters, such as MAS frequencies and nominal rf field strengths for the individual pulse sequences that were investigated, is given in the Supplement (Table S1).
For all experimental schemes treated here, the general form of the time-dependent rf field Hamiltonian is given by Eq. (19). Deviations from the nominal rf amplitude and phase due to the rf inhomogeneity are introduced by the relative amplitude and phase, denoted by , and ϕrel(r,t), respectively. These parameters depend on the position of the crystallite in the sample space r and the rotor orientation. Under MAS, they are modulated by the rotor frequency. The time-dependent trajectories were computed numerically (see Sect. 2) and are given as an input to the numerical simulations.
Simulations were performed for volume elements of the rz plane indicated in yellow in Fig. 1a, with an initial orientation of . For nutation experiments (see Sect. 5.1), several ϑ0 values were considered, as the spin dynamics are, in principle, dependent on ϑ0 due to the non-commuting Hamiltonians at different time points during a rotor period. The sample space was restricted to the length of the coil along the rotor axis (indicated by dashed lines in Fig. 1b). Potential effects of the radial inhomogeneity should be similar for samples exceeding the length of the coil but might be more pronounced since the magnitude of rf amplitude and phase modulations increases towards the rotor edges (see Fig. 2). Simulation results of the individual volume elements were summed up during data processing and weighted with r to account for the increase in volume with radial distance. The coil sensitivity (reciprocity theorem; Hoult and Richards, 1976; Tošner et al., 2017) was taken into account by additional weighting of each rz element with the average relative rf amplitude over a rotor cycle .
In order to separate the effect of the static rf field inhomogeneity from time-dependent effects due to amplitude and phase modulation arising from sample rotations, the spin dynamics were simulated under different conditions. Amplitude and phase modulations were considered separately and either treated as time dependent or as the static average over a rotor period. The four following cases considered in this work are denoted as C1–C4, where:
C1 – time-averaged constant amplitude and zero phase;
C2 – time-dependent amplitude and zero phase;
C3 – time-averaged constant amplitude and time-dependent phase;
C4 – time-dependent amplitude and time-dependent phase;
Simulations with the time-averaged constant phase were not performed as constant phase offsets are small (maximum of less than 5∘), and the absolute phase of the rf irradiation has no influence on the outcome of an experiment. For reference, Table 1 summarizes the treatment of amplitude and phase modulations for these four cases.
Experiments were performed on a 500 MHz Bruker Avance III HD NMR spectrometer equipped with a Bruker 1.9 mm triple-resonance MAS probe in double-resonance configuration at a temperature of 285 K. All powdered samples (natural-abundance glycine and natural-abundance L-histidine ) were purchased from commercial sources and used without further purification. Nutation spectra of glycine were recorded as two-dimensional experiments without sign discrimination in t1 (simple sine amplitude modulation) at MAS frequencies of 15 and 30 kHz. The nominal rf field amplitude was calibrated to 100 kHz using a nutation spectrum. Spectra were recorded with 512 t1 increments with 12 scans each and a time increment for the nutation pulse of 2.5 µs. The spectral width in the direct dimension was set to 100 kHz, and 1024 complex data points were acquired. MATLAB (The MathWorks Inc., Natick, MA, USA) was used for data processing using a cosine-squared window function. The two-dimensional proton–proton correlation spectra of L-histidine with frequency-switched Lee–Goldburg (FSLG) decoupling (Bielecki et al., 1989, 1990; Mote et al., 2016) in the indirect dimension were recorded at MAS frequencies of 14 and 28 kHz. Spectra were acquired with 512 t1 increments with eight scans each and time increments between 43.2 and 48 µs. States-type (States et al., 1982) data acquisition was used for phase-sensitive detection and sign discrimination in t1. The spectral width in the direct dimension was set to 200 kHz, and 1024 complex data points were recorded. Nutation-frequency-selective I-BURP-2 (Geen and Freeman, 1991) pulses in the spin lock frame were used for sample restriction (Aebischer et al., 2020). The rf field amplitudes were calibrated using a nutation spectrum and set to 100 kHz during hard pulses and spin lock. The FSLG decoupling was implemented using shaped pulses with a time resolution of 100 ns for the phase ramp. Shape files with 80, 108, and 160 points were used, corresponding to nutation frequencies about the effective field of 250, 185.2, and 125 kHz. The carrier was placed outside the spectral region of interest, and its position is indicated by an arrow. Spectra were processed in MATLAB with zero-filling to 4096×4096 data points and the application of a cosine-squared apodization. The 1D spectra shown were obtained by summation over the relevant spectral region in ω2. Frequency axes in parts per million (ppm) were determined by comparison of the peak positions observed for the α and δ2 proton resonances in histidine with those found in the literature (Mithu et al., 2013).
In this section, we discuss how a number of common solid-state NMR experiments are affected by the MAS time-modulated radio-frequency fields. This is done by analytical calculations based on the Floquet description presented in Sect. 3.1, numerical simulations, and, in some cases, using experimental data.
5.1 Nutation spectroscopy
5.1.1 Numerical simulations and experimental results
Nutation spectra represent a simple method for characterizing the rf field distribution in the sample. Such spectra were simulated for one-spin systems, and the rf inhomogeneity was included in the rf Hamiltonian of Eq. (19), setting the nominal phase of the rf field to zero corresponding to rf irradiation along the x axis. The nominal rf field amplitude ν1,nom was set to 100 kHz, and the four cases C1–C4 (see Table 1) were studied. As only isotropic spin interactions were considered, simulations were performed for a single crystallite orientation.
Simulated nutation spectra, using the rf field profiles of the 3.2 and 1.3 mm MAS probes at different spinning frequencies (15 and 30 kHz, respectively), are shown in Fig. 3. The overall nutation profile of the 1.3 mm probe is narrower, indicating a more homogeneous rf field distribution inside the coil and, thus, a less pronounced drop-off of the static rf field amplitude along the rotor axis. Phase modulation of the rf field (C3 and C4) leads to sidebands at ( are visible). The intensity of these sidebands increases with increasing MAS frequency. Amplitude modulation of the applied rf field, on the other hand (C2 and C4), leads to sidebands at (only m=1 visible). These sidebands are significantly weaker than those arising due to phase modulations, and their intensity increases with decreasing spinning frequency. The reduced intensity of these amplitude modulation sidebands can be explained by the fact that their position depends on the magnitude of the static rf field amplitude and, thus, varies depending on the position within the sample space. Both types of sidebands are weaker in the 1.3 mm rf profile, indicating that rf amplitude and phase modulations are less pronounced in comparison to the 3.2 mm rf profile. The phases of the sidebands depend on the initial position ϑ0 of the simulated rz plane (see Fig. S1 in the Supplement). However, the obtained spectra are very similar for all initial orientations, and no significant influence of ϑ0 on the effects of the radial rf inhomogeneity has been observed.
Experimental 1H nutation spectra of natural-abundance glycine measured at a proton resonance frequency of 500 MHz, using a Bruker 1.9 mm MAS probe, are shown in Fig. 4 for MAS frequencies of 30 kHz (Fig. 4a) and 15 kHz (Fig. 4b). As was observed in the simulated nutation spectra (see Fig. 3), sidebands at due to rf phase modulations are visible in the experimental spectra. Moreover, sidebands at are visible that replicate the shape of the main nutation profile at 100 kHz. At a lower MAS frequency, these sidebands increase in intensity, whereas those at multiples of the rotor frequency are attenuated. In the simulated spectra shown in Fig. 3, the sidebands at 115 and 130 kHz were significantly weaker and did not have the same shape as the overall nutation profile. As is shown in Fig. S2, no such sidebands are observed in the experimental nutation spectra of natural-abundance adamantane (see Fig. S2a), indicating that they arise from the MAS modulation of anisotropic interactions. This is confirmed by the simulated nutation spectra obtained for a dipolar-coupled two-spin system (see Fig. S2b), where strong sidebands at ν1±νr are obtained that nicely replicate the shape of the main nutation profile for all four cases (C1–C4).
5.1.2 Floquet analysis
Phase modulation of the rf field leads to non-commuting terms in the rf Hamiltonian at different points in time, thus prohibiting an analytical determination of the time evolution of the magnetization during the nutation experiment. However, insight can be gained from the interaction frame trajectory of spin operators that can be computed numerically. For the spin operator, this trajectory can be expanded as follows:
Fourier analysis of the time-dependent aχ(t) coefficients then yields the frequency components present in the nutation spectrum. Such interaction frame trajectories of the spin operator with rf irradiation along the x axis were computed numerically in MATLAB with a time resolution of 50 ns. A nominal rf field strength of 100 kHz was chosen and a MAS frequency of 30 kHz assumed. The effects of MAS time-dependent rf amplitude and phase modulations were studied separately. Modulations were modelled as Fourier series (see also Eqs. 1 and 2), and magnitude and phase coefficients were given as input.
Absolute values and phases of the aχ(t) coefficients are shown in Figs. 5a and b as a function of the magnitude of amplitude modulations with ωr and 2⋅ωr ( and ). The static amplitude offset was set to 1, and all with n≥3 were set to zero. No phase modulation was taken into account. Amplitude modulation with the rotor frequency (Fig. 5a) leads to sidebands at , with m being any integer, whereas amplitude modulation with 2⋅ωr (Fig. 5b) leads to sidebands at . The intensity of these sidebands increases with the magnitude of the modulation in both cases. However, sidebands arising from rf amplitude modulation with the base frequency ωr are significantly stronger. The phase of the amplitude modulation for the spectra in Fig. 5 was set to zero as it only influences the phases of the centre band and sidebands. Static amplitude offsets will simply shift the entire spectrum. Figure 5c and d show the absolute values and phases of aχ(t) for phase modulations with ωr and 2⋅ωr ( and ). The static phase offset and all with n≥3 were set to zero. No rf amplitude modulations were taken into account. The frequency range shown in the figure was limited to 0–90 kHz since rf phase modulations lead to sidebands at that are significantly less intense than the main band at the nominal rf amplitude (ν1= 100 kHz). In contrast to the sidebands observed for rf amplitude modulation, phase modulation with n⋅ωr exclusively leads to sidebands at . Compared to the sidebands arising from amplitude modulation, the intensities of these phase modulation sidebands are considerably lower. Their intensity increases with the Fourier number n of the modulation (sidebands arising from are more intense than the ones from ). The phase of the modulation () was set to zero again, as it only affects the phase of the aχ(t) coefficients (see also Fig. S1).
These results are in good agreement with the simulated and experimental nutation spectra shown in Figs. 3 and 4, where two separate families of sidebands arose for rf amplitude and phase modulation. As described earlier, the higher intensity observed for phase modulation sidebands in these spectra can be explained by the overall larger magnitude of the phase modulations and the independence of the sideband position from the average rf field amplitude. The position of the amplitude modulation sidebands, on the other hand, shifts with the average static rf field amplitude, leading to a broadening of the sidebands.
We have also looked at the effects of the radial rf field inhomogeneity in the context of spin lock experiments which is closely related to the nutation experiment. This connects to the first experimental observation of such effects in rotary resonance recoupling experiments, where additional peaks in the centre of the expected dipolar doublet have been observed. This was attributed to time-dependent phase modulations (Levitt et al., 1988). A more detailed study of this experiment is described in Sect. S4.
Hartmann–Hahn cross-polarization (Hartmann and Hahn, 1962; Stejskal et al., 1977) is probably the most ubiquitous pulse sequence element in solid-state NMR. Under MAS, the sum, or difference, of the two rf field amplitudes has to be matched to an integer multiple of the spinning frequency as follows:
Due to the rf field inhomogeneity across the sample, this condition cannot be fulfilled simultaneously in the entire sample volume, and only certain parts of the sample will participate in the polarization transfer, thus decreasing the resulting signal intensity. One popular strategy to overcome this volume selectivity is ramped-amplitude cross-polarization (Metz et al., 1994) or the adiabatic modulation of the rf field amplitude during the contact time (Hediger et al., 1995, see Sect. S5 for more details). In this work magnetization transfers at the n=1 zero-quantum matching conditions in heteronuclear CN, HN, and HC two-spin systems, were simulated for standard, ramped-amplitude, and adiabatic-passage CP experiments for the 3.2 mm MAS probe at a proton resonance frequency of 600 MHz. A MAS frequency of 20 kHz was assumed, and the nominal rf fields and contact times were set for CN as follows: 85 kHz on C, 65 kHz on N, and 5 ms. The nominal rf fields and contact times were set for HN as follows: 70 kHz on H, 50 kHz on N, and 1 ms. The nominal rf fields and contact times were set for HC as follows: 90 kHz on H, 70 kHz on C, and 1 ms. The anisotropy of the dipolar-coupling tensor was estimated from average bond lengths, and values of 1.9, 25, and −46 kHz were used for the CN, HN, and HC simulations, respectively. Chemical shifts as well as J coupling constants were set to zero. The simulations were performed with a time resolution of 250 ns, and the x magnetization of both the source and the destination spin was detected every 5 µs. Powder averaging was performed over 1154 crystallite orientations.
The simulated time evolution of the spin-locked x magnetization on both spins for the cases C1–C4 (see Table 1) is shown in Fig. 6 for all three spin pairs. The rf amplitude on the source spin was either kept constant (left-hand column), modulated with a linear ramp (middle column), or used a tangential modulation (right-hand column; see Fig. S4 for more details). In CN spin pairs, a tangential modulation of the rf amplitude on one of the spins leads to a significant improvement of the transfer efficiency (up to 35 %) compared to both the standard and the ramped amplitude CP experiment (less than 20 %). For spin pairs with stronger dipolar couplings, such as HN and HC, both the ramped amplitude CP and the adiabatic passage CP lead to similar transfer efficiencies of up to 70 %. In all simulated experiments, only marginal differences between the obtained transfer efficiencies for the four cases C1–C4 are observed for all spin pairs. Time-dependent modulations of the rf amplitude and phase due to the radial rf inhomogeneity, therefore, do not seem to have a significant effect on the magnetization build up on the destination spin. On the source spin, some magnetization is lost when rf phase modulations are present (C3 and C4), which is also observed for one-spin spin lock simulations. As the radial contributions to the rf inhomogeneity are weaker in the 1.9 and 1.3 mm probes, similar results would be expected for these probes. Overall, these simulation results suggest that the effect of the radial inhomogeneity on CP polarization transfers is negligible. Only the static rf amplitude offset over the relevant sample space is important due to the volume selectivity it causes.
Moreover, polarization transfers in NCA and NCO two-spin systems using the tm-SPICE sequences (Tošner et al., 2018) were simulated. These pulse schemes were developed using optimal control (OC) strategies, taking into account the MAS modulations of the rf field due to the radial rf inhomogeneity. The resulting magnetization transfers are shown in Fig. 7 for the 3.2 mm MAS probe at a proton resonance frequency of 400 MHz. Nominal rf amplitudes on both channels were set to 40 kHz, and a spinning speed of 20 kHz was assumed. The shape files for the pulse sequences contain 1750 points, with a time resolution of 2 µs, corresponding to a contact time of 3.5 ms. The time resolution for the propagation was set to 250 ns, and the expectation value of the and operators was detected every 5 µs. Powder averaging was performed over 1154 crystallite orientations. Spin system parameters (J couplings, chemical shift anisotropy, CSA, and dipolar coupling tensors) were taken from Tošner et al. (2018) and can be found in Tables S2 and S3. Impressive transfer efficiencies of around 60 % are obtained for both NCA (Fig. 7a) and NCO (Fig. 7b) spin pairs. However, only minor differences between the four cases C1–C4 are observed in these simulations. This reflects the fact that the optimization of this sequence took rf amplitude and phase modulations of different magnitude, as well as varying initial phases of these modulations, into account. Therefore, the sequence performs well under all possible conditions encountered in the rotor, leading to an increase in the robustness of the resulting sequences towards static and time-dependent rf inhomogeneity. The broad range of considered conditions stabilizes the optimization towards a broader minimum that gives a better transfer over the complete rotor.
5.3 Rotational-echo double resonance
In REDOR recoupling (Gullion and Schaefer, 1989a, b), the heteronuclear dipolar coupling is reintroduced by trains of two rotor-synchronized π pulses per rotor cycle. This technique allows the quantitative measurement of dipolar couplings in heteronuclear spin pairs and has become a valuable tool in the characterization of structure (Hong, 2006; Rienstra et al., 2002; Michal and Jelinski, 1997; Jia et al., 2015) and dynamics (Schanda et al., 2010).
Numerical simulations of REDOR recoupling were performed for CN and HN spin pairs at a proton resonance frequency of 600 MHz, with an XY-4 phase cycling scheme (Gullion et al., 1990) as it is commonly implemented to generate a pure Ising-type Hamiltonian. Dipolar couplings were estimated from average bond lengths, and the anisotropy of the coupling was set to 2 and 24 kHz in CN and HN, respectively. Chemical shift tensors and scalar J couplings were set to zero. Resulting simulated REDOR curves for CN spin pairs are shown in Fig. 8 for a 1.3 mm (Fig. 8a) and a 3.2 mm (Fig. 8b) probe, assuming a spinning frequency of 20 kHz. The nominal rf field strengths were set to 100 kHz (62.5 kHz) on C and 65 kHz (50 kHz) on N in the 1.3 mm (3.2 mm) probe. A time resolution of 250 ns was chosen for the propagation, and 538 crystallite orientations were used for the powder averaging. Compared to the theoretical REDOR curve (dashed line; analytical expression including finite-pulse effects can be found in Jaroniec et al., 2000), considerably lower recoupling efficiencies are obtained for C1–C4 in both probes (Nishimura et al., 2001). However, only minor differences between the four cases are observed with amplitude modulations (C2 and C4), leading to a slight deterioration of the recoupling performance. These effects are very similar for both probes. The relative timing of the rotor-synchronized π pulses in the REDOR sequence with respect to the time-dependent rf field amplitude and phase modulations only has a marginal effect on the recoupling performance (see Fig. S6 for further details). Simulated REDOR curves for the HN spin system in the 1.9 mm probe, assuming a spinning frequency of 40 kHz, are shown in Fig. 8c and d. As the dipolar coupling in the HN spin pair is too large to allow sufficient sampling of the REDOR curve, modified REDOR implementations were simulated in which one (Gullion and Schaefer, 1989b, Fig. 8c) or both (Jain et al., 2019, Fig. 8d) of the pulses in the basic building block are shifted. The corresponding pulse sequences are shown in Fig. S5. These schemes lead to a scaling of the effective dipolar coupling and, thus, allow sufficient sampling of the REDOR curve even for strongly dipolar-coupled spin pairs. A time resolution of 125 ns was chosen, and the nominal rf fields were set to 125 kHz on H and 50 kHz on N. Powder averaging was performed for 10 000 crystallite orientations. Theoretical REDOR curves are again indicated by the dashed lines (analytical expressions, including finite-pulse effects, can be found in Schanda et al., 2011, for Fig. 8c and Jain et al., 2019, for Fig. 8d). For both shifting regimes, only slight deviations from the theoretical curves are observed. Moreover, resulting REDOR curves for C1–C4 are identical, indicating that time-dependent amplitude and phase modulations have no effect on the recoupling performance in strongly dipolar-coupled spin pairs. Overall, the REDOR sequence seems to be predominantly affected by the static rf inhomogeneity which causes deviations in the pulse flip angles from the desired 180∘ due to average rf field amplitude deviations.
5.4 Symmetry-based CN recoupling – C7 and POST-C7
Symmetry-based C sequences represent an important class of homonuclear recoupling sequences. Since the first introduction of the original C sequence (Lee et al., 1995), many other symmetry-based sequences have been proposed and characterized; however, only the C and the POST-C7 sequence (Hohwy et al., 1998), where the basic C is replaced by the cyclically permuted C POST element, were considered in this work.
Numerical simulations of polarization transfers during C7 and POST-C7 recoupling were performed for CC two-spin systems under conditions typical for a 3.2 mm MAS probe at a carbon resonance frequency of 150 MHz. The nominal rf field amplitude was set to 70 kHz, and a spinning frequency of 10 kHz assumed. A time resolution of approximately 714 ns was chosen, and 538 crystallite orientations were used for the powder averaging. The time evolution of the expectation values of and in a CC spin pair with isotropic chemical shifts that are symmetric around zero () is shown in Fig. 9 for C7 (Fig. 9a) and POST-C7 (Fig. 9b) for the cases C1–C4. The anisotropy of the dipolar coupling tensor was estimated from average bond lengths and set to 4.5 kHz. For both sequences, transfer efficiencies of approximately 70 % are achieved in the first transient for a mixing time of approximately 10 ms. Amplitude modulations due to the radial rf field inhomogeneity (C2 and C4) lead to a slight deterioration in the recoupling performance. This effect is more pronounced for C7, indicating the improved robustness of the POST-C7 sequence. At longer times, a loss of magnetization on both spins is observed when time-dependent amplitude modulations are taken into account. Interestingly, the overall order of the cases is different for the two sequences. For C7 recoupling, higher transfer efficiencies are observed for C4 in comparison to C2, whereas this order is reversed for POST-C7. Thus POST-C7 seems to be more sensitive to combined amplitude and phase modulations (C4). Figure 9c (C7) and d (POST-C7) show simulation results of the same spin system for a spatially restricted sample (central third along the rotor axis). This restriction of the sample space mitigates the effects of amplitude modulations, and only very marginal differences between the four cases are observed for both sequences. In order to further investigate the robustness of the two sequences, simulations in a second model system with a large CSA tensor were performed at a lower external magnetic field (75 MHz carbon resonance frequency). The parameters of this spin system were based on phthalic acid (Hellwagner et al., 2017, see Table S4). Significantly lower overall transfer efficiencies (approximately 50 %) are observed (see Fig. 9e and f), and for both sequences, rf field amplitude modulations (C2 and C4) further deteriorate the recoupling efficiency. This decrease is less pronounced for POST-C7 (see Fig. 9f), again indicating its improved robustness. We have not carried out a complete Floquet analysis of the effects caused by the time-dependent rf field amplitudes. We believe that the decreased efficiency is due to the appearance of effective fields that shift the resonance condition slightly as is the case for pulse transients (Hellwagner et al., 2017).
5.5 Frequency-switched Lee–Goldburg decoupling
5.5.1 Numerical simulation and experimental results
Frequency-switched Lee–Goldburg (FSLG) decoupling is a homonuclear dipolar decoupling technique that can be used in combination with MAS to improve resolution of spectra for dipolar-coupled homonuclear spin systems (Lee and Goldburg, 1965; Goldburg and Lee, 1963; Bielecki et al., 1989, 1990). The experiment is based on off-resonance rf irradiation, leading to a truncation of the second-rank spin tensor of the homonuclear dipolar coupling by an effective radio-frequency field inclined at an angle with respect to the static magnetic field. Experimentally, FSLG can also be implemented using on-resonance irradiation with a constant rf field amplitude and a continuous phase ramp to generate the frequency offset. The total cycle time is divided into two intervals of equal length during which the phase is rotated in opposite directions (inverting the offset) and with a phase jump of 180∘ in between.
The effects of the radial part of the rf field inhomogeneity on the residual linewidth under FSLG decoupling were simulated for a homonuclear dipolar-coupled three-spin system in a 3.2 mm MAS probe at a proton resonance frequency of 600 MHz, assuming a MAS frequency of 12.5 kHz. The nominal rf field amplitude was set to 102.06 kHz, corresponding to a tilting along the magic angle of an effective field with a strength of 125 kHz (FSLG cycle time of 16 µs). Using such a synchronization between the FSLG sequence and the sample spinning makes the simulations much more efficient than an asynchronous implementation, while at the same time avoiding all resonance conditions up to and including the second order. The FSLG decoupling was implemented using a phase ramp with a time resolution of 50 ns. The same time resolution was chosen for the propagation of the Hamiltonian. The initial density operator was set to and transverse magnetization components detected every 48 µs (three FSLG cycles). A total of 8192 data points were acquired, and the free induction decay (FID) processed in MATLAB. Powder averaging was performed over 1154 orientations. The parameters characterizing the chemical shift and dipolar coupling tensors were chosen to mimic a CH2 group, with couplings to an additional remote spin, and can be found in Tables S5 and S6. Scalar J couplings were neglected and set to zero.
Simulated spectra of the three-spin system are shown in Fig. 10a for C1–C4. In all four cases, all resonances show strong asymmetric features on the left-hand side of the spectral line due to the distribution of the chemical shift scaling factors (Hellwagner et al., 2020). Considerable additional line broadening is observed when amplitude modulations are taken into account (C2 and C4). As the same linewidths were obtained in simulations of an asynchronous implementation of FSLG decoupling (MAS frequency of approximately 14.1 kHz; see Fig. S7), the broadening is not caused by resonance effects. This effect is observed for all three resonances but is most pronounced for the CH2 resonance around 1.25 kHz. The additional time dependence of the rf phase in C4 results in no additional broadening, and the two remaining cases (C1 and C3) are indistinguishable. Phase modulation, therefore, does not seem to have an influence on the obtained linewidth. Figure 10 also shows simulated FSLG spectra for radial slices of the simulated rz plane at r=0.65 (Fig. 10b) and 1.3 mm (Fig. 10c) and for a spatially restricted sample space (central third along the rotor axis; –1.35 mm; all r values; Fig. 10d). As the magnitude of the rf amplitude modulations increases towards the rotor edges (see Fig. 2), significantly stronger broadening is observed for radial slices closer to the coil windings. Spatial restriction of the sample space to the central third leads to a reduction in the linewidth. This line narrowing is significantly more pronounced for the resonance at −2.75 kHz, where the foot on the left-hand side of the resonance is eliminated for all four cases. For C2 and C4, broadening is observed even in this spatially restricted sample, indicating that time-dependent amplitude modulations without a static rf amplitude offset still result in contributions to the residual linewidth. Simulations were also performed for a six-spin system, and qualitatively similar results were obtained.
In order to observe this broadening experimentally, the sample space has to be restricted to areas close to the coil windings where strong rf field amplitude modulations occur. This could, in principle, be achieved by physically restricting the sample using cylindrical spacers. However, homogeneous packing in such a sample is difficult to achieve. Alternatively, nutation-frequency-selective pulses, as described in Aebischer et al. (2020), can be used to select the desired areas which also correspond to high average rf field amplitudes. Figure 11 shows the FSLG decoupled proton spectra of L-histidine measured at a proton resonance frequency of 500 MHz in a Bruker 1.9 mm MAS probe, using a 2 ms I-BURP-2 pulse in the spin lock frame for the B1 field selection of areas where the rf field amplitude corresponds to the nominal value (see Fig. S9 for a simulated inversion profile). Spectra were recorded with different B1 field strengths for the FSLG decoupling at spinning frequencies of 14 and 28 kHz. No significant improvement in the obtained linewidth is observed for higher MAS frequencies and stronger decoupling field strengths. This indicates that the residual linewidth in these spectra is not decoupling limited which, therefore, prohibits the experimental characterization of the additional broadening caused by rf field amplitude modulations due to the radial rf field inhomogeneity.
5.5.2 Floquet analysis
In order to gain physical insight into the origin of the observed line broadening in FSLG-decoupled spectra due to rf field amplitude modulations, scaling factors for the first- and second-order contributions to the effective Hamiltonian were computed (see Sect. 3.1.1 for more details). As a simple measure for the magnitude of contributing first-order terms the norms of one-spin coefficients,
and two-spin coefficients,
Interaction frame trajectories using the rf field distribution in a 3.2 mm MAS probe during FSLG decoupling were computed numerically in MATLAB with a time resolution of 50 ns and the Fourier coefficients extracted. The spinning frequency was chosen to be 12.5 kHz and the effective field strength along the magic angle set to 125 kHz, corresponding to a modulation frequency of the rf Hamiltonian of 62.5 kHz. This leads to the synchronization of the MAS rotation and the rf irradiation after a single rotor cycle or five FSLG cycles. This choice of frequencies should avoid all resonance conditions up to and including the second order. Relative rf field amplitude and phase modulations were modelled as Fourier series, and fitted Fourier coefficients up to the fourth order were used as input (see Sect. 2). In analogy to the treatment of the rf field inhomogeneity in numerical simulations, amplitude and phase modulations were considered separately, and the four cases C1–C4 summarized in Table 1 studied.
A full FSLG cycle assuming a time-independent rf Hamiltonian and an ideal phase ramp with a 180∘ phase shift in the middle consists of two β rotations with an opposite direction about the effective field. The overall propagator would, thus, be the unity operator and ωeff=0 in Eq. (12). However, the time-dependent modulations of the rf amplitude and phase due to the radial rf inhomogeneity can give rise to an additional effective field. The magnitude of this field as a function of the position within the sample space in the 3.2 mm probe is shown in Fig. 12. As the static rf field inhomogeneity does not lead to additional effective fields, only cases in which either the rf field amplitude, the phase, or both are time dependent (C2, C3 and C4) are shown. Effective fields arise mainly at the edges of the rotor (large r and z) where modulations are strongest. Amplitude and phase modulations alone (C2 and C3 respectively) lead to very small effective fields (max. 50 Hz), whereas larger effective fields (up to 400 Hz) result for combined amplitude and phase modulations (C4). In comparison to the rotor frequency and the basic modulation frequency of the FSLG sequence, these additional fields are small and will most likely not have any significant effects, except for a small change in the effective field direction and magnitude.
In the first-order approximation, the relevant scaling factors are those of the chemical shift ( with ) and those of the dipolar coupling ( with ). These can contribute to the first-order effective Hamiltonian (see Eq. 27), since the modulation by the rf field amplitude can be compensated by the time dependence due to MAS. The resulting norm of the coefficients (Eq. 38) as a function of the position within the sample space is shown in Fig. 13a for C1–C4. As coefficients are symmetric (), only those corresponding to and 2 are shown. The scaling of the isotropic chemical shift (k=0) is close to the ideal value of cos (θm)≈0.577 in regions of the rotor where the rf field amplitude is comparable to the nominal value. Towards the edges of the rotor, the rf field amplitude decreases, leading to a smaller tilt angle of the effective field during FSLG and, thus, an increase in the scaling factor. The time-modulated part of the rf field inhomogeneity does not appear to have any influence on the isotropic chemical shift, as no significant differences between the four cases are observed. Time-dependent amplitude modulations (C2 and C4) lead to additional non-zero coefficients for chemical shift contributions with k≠0 that can also contribute to the first-order effective Hamiltonian for (see Eq. 27) where parts of the CSA tensor become time independent. These contributions will be strongest at the very edges of the sample space (large r and z), but non-zero coefficients are also obtained in the central third of the rotor close to the coil windings.
Under ideal conditions, the FSLG decoupling scheme leads to the averaging of the anisotropic dipolar coupling in the first-order approximation, and the corresponding scaling factors would be zero. However, dipolar coupling terms are reintroduced when rf modulations are taken into account. The norm of the relevant coefficients (Eq. 39) is shown in Fig. 13b for the 3.2 mm MAS probe. Again, only the terms can contribute to the first-order effective Hamiltonian and partially reintroduce Fourier components of the dipolar coupling. As was the case for the chemical shift scaling factors, the coefficients are symmetric, and thus, only those for k=1 and 2 are shown. Amplitude modulations (C2 and C4) lead to significant k=1 scaling factors, and non-zero coefficients are not only obtained at the very edges of the rotor but also in the central third close to the coil windings. The additional phase modulation in C4 does not have an influence, and amplitude modulations alone thus seem to be responsible for the reintroduction of the first-order coupling terms. The contribution of individual coefficients to the norm are shown in Fig. S10 for C4. Significant are obtained for , , and terms. These first-order time-independent homonuclear coupling terms contribute to the residual linewidth under FSLG decoupling and, thus, lead to an additional line broadening. Numerical simulations taking only the first-order effective Hamiltonian into account confirmed that the observed line broadening for C4 in simulated spectra (see Fig. 10) can indeed be attributed to the first-order contributions to the effective Hamiltonian (see Fig. S11).
In principle, the second-order effective Hamiltonian during FSLG decoupling contains three types of commutator cross-terms. However, contributions from chemical shift cross-terms () only contain one-spin operators and will, thus, lead to an additional effective field and will only weakly influence the residual linewidth under FSLG by changing the direction or magnitude of the effective field. The same is true for the one-spin component of the dipolar–dipolar cross-terms (). This leaves only two sources of coupling terms in the second-order effective Hamiltonian, namely the three-spin contribution of dipolar–dipolar cross-terms and commutators between chemical shift and dipolar terms. Out of these two, the former will be most relevant for the residual linewidth as dipolar couplings are generally much larger than typical chemical shifts. The corresponding scaling factors for the three-spin contribution in were computed according to Eq. (35). Because the effective fields generated by the modulations of the rf field amplitude and phase are small (see Fig. 12), ν=0 was excluded from the summation in order to avoid near-resonance conditions and the norm,
was computed to characterize the strength of the three-spin coupling terms. Logarithmic contour plots for the resulting for C1–C4 in the 3.2 mm MAS probe are shown in Fig. 14. Scaling factors are shown for n1=1 (Fig. 14a) and n1=2 (Fig. 14b) and all possible values of the index n2. It can be seen that time-dependent rf modulations (C2, C3, and C4) increase the magnitude of the second-order cross-terms significantly compared to static rf amplitude and phase offsets alone (C1). However, the observed scaling factors are still negligible compared to the magnitude of the first-order terms, and no significant effect on the linewidth would be expected. In all four cases, substantially higher scaling factors are obtained for pairs of indices where . These increase strongly towards the edges of the rotor along the rotor axis, but since no difference between the cases C1–C4 is observed, they do not seem to be influenced by the radial part of the rf inhomogeneity.
The analysis of the scaling factors of the terms contributing to the effective Hamiltonian up to the second order suggests that the static part of the rf inhomogeneity has a significant influence on the isotropic chemical shift scaling and also leads to stronger second-order contributions. However, the overall magnitude of these second-order terms remains small compared to first-order contributions. Time-dependent rf amplitude modulations have pronounced first-order effects and lead to the reintroduction of anisotropic chemical shift and dipolar coupling terms that cause line broadening (see Fig. S11). No such effects were observed for phase modulations.
Magic angle spinning in combination with inhomogeneous radial rf fields leads to a time-dependent modulation of the rf field amplitude and phase. We have investigated the effect of these time-dependent rf fields on some common solid-state NMR pulse sequences using numerical simulations and an analytical approach based on Floquet theory. In none of the investigated building blocks used in solid-state NMR experiments could we find significant effects from such time-dependent rf fields. In nutation spectra, two distinct families of sidebands, arising due to rf field amplitude and rf field phase modulations, respectively, were observed in simulated and experimental spectra. The intensity of these sidebands can help to characterize the strength of the modulations and, thus, to give insights into the radial contribution to the rf field inhomogeneity for a given MAS probe. In the polarization transfer sequences, like Hartmann–Hahn cross-polarization, REDOR, and C7, only minor effects were observed that will most likely be of no consequence for experimental implementations. In all these sequences, the static rf field inhomogeneity over the sample volume played a much larger role and leads to significant performance degradation.
In simulations of homonuclear FSLG decoupling, considerable line broadening was observed for rf field amplitude modulations. Floquet analysis of the effective Hamiltonian up to the second order revealed that this broadening is most likely due to the reintroduction of homonuclear coupling terms to the first order caused by the MAS modulation of the rf field amplitude. However, no experimental characterization of this effect was possible as the experimentally obtained linewidths were not limited by the homonuclear decoupling. Overall, the results presented in this work suggest that the influence of the MAS modulation of the rf field amplitude and phase in many pulse sequences is small and, thus, negligible for typical experimental implementations. Moreover, they manifest themselves in areas of the sample space close to the rotor edges and can, thus, be reduced by physical or virtual sample restriction. Nevertheless, these modulations can become relevant in the development of new pulse sequences based on optimal control strategies and should be taken into account in their development in order to increase their robustness towards rf inhomogeneity and enlarge the NMR-responsive sample volume.
The supplement related to this article is available online at: https://doi.org/10.5194/mr-2-523-2021-supplement.
ME designed the research. ZT provided the data about rf field inhomogeneity in probes. KA carried out all measurements and simulations with some help from ME. All the authors discussed and interpreted the results and were involved in writing the paper.
Matthias Ernst is an executive editor of MR.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We would like to thank Perunthiruthy K. Madhu, Kaustubh Mote, and Johannes Hellwagner for the insightful discussions about theory and the experimental implementation of homonuclear decoupling. Beat H. Meier and Alexander Barnes are acknowledged for providing measurement time for the project.
This research has been supported by the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. 200020_188988) and the Czech Science Foundation (GACR; grant no. 20-00166J).
This paper was edited by Bernd Reif and reviewed by Malcolm Levitt and one anonymous referee.
Aebischer, K., Wili, N., Tošner, Z., and Ernst, M.: Using nutation-frequency-selective pulses to reduce radio-frequency field inhomogeneity in solid-state NMR, Magn. Reson., 1, 187–195, https://doi.org/10.5194/mr-1-187-2020, 2020. a, b, c, d
Aebischer, K., Tošner, Z., and Ernst, M.: Effects of radial radio-frequency field inhomogeneity on MAS solid-state NMR experiments, ETH Zurich [data set], https://doi.org/10.3929/ethz-b-000488476, 2021. a
Bielecki, A., Kolbert, A. C., De Groot, H. J. M., Griffin, R. G., and Levitt, M. H.: Frequency-switched Lee–Goldburg sequences in solids, Advances in Magnetic and Optical Resonance, 14, 111–124, 1990. a, b, c
Charmont, P., Lesage, A., Steuernagel, S., Engelke, F., and Emsley, L.: Sample Restriction Using Magnetic Field Gradients in High-Resolution Solid-State NMR, J. Magn. Reson., 145, 334–339, https://doi.org/10.1006/jmre.2000.2117, 2000. a
Charmont, P., Sakellariou, D., and Emsley, L.: Sample Restriction Using Radiofrequency Field Selective Pulses in High-Resolution Solid-State NMR, J. Magn. Reson., 154, 136–141, https://doi.org/10.1006/jmre.2001.2467, 2002. a
Cheng, V., Suzukawa, H., and Wolfsberg, M.: Investigations of a nonrandom numerical method for multidimensional integration, J. Chem. Phys., 59, 3992–3999, 1973. a
Elbayed, K., Dillmann, B., Raya, J., Piotto, M., and Engelke, F.: Field modulation effects induced by sample spinning: application to high-resolution magic angle spinning NMR, J. Magn. Reson., 174, 2–26, https://doi.org/10.1016/j.jmr.2004.11.017, 2005. a
Frantsuzov, I., Vasa, S. K., Ernst, M., Brown, S. P., Zorin, V., Kentgens, A. P. M., and Hodgkinson, P.: Rationalising Heteronuclear Decoupling in Refocussing Applications of Solid-State NMR Spectroscopy, ChemPhysChem, 18, 394–405, https://doi.org/10.1002/cphc.201601003, 2017. a
Geen, H. and Freeman, R.: Band-selective radiofrequency pulses, J. Magn. Reson., 93, 93–141, 1991. a
Guenneugues, M., Berthault, P., and Desvaux, H.: A Method for Determining B1 Field Inhomogeneity. Are the Biases Assumed in Heteronuclear Relaxation Experiments Usually Underestimated?, J. Magn. Reson., 136, 118–126, https://doi.org/10.1006/jmre.1998.1590, 1999. a
Gullion, T. and Schaefer, J.: Detection of Weak Heteronuclear Dipolar Coupling by Rotational-Echo Double-Resonance Nuclear Magnetic Resonance, Advances in Magnetic and Optical Resonance, 13, 57–83, https://doi.org/10.1016/B978-0-12-025513-9.50009-4, 1989b. a, b
Hediger, S., Meier, B. H., and Ernst, R. R.: Adiabatic passage Hartmann–Hahn cross polarization in NMR under magic angle sample spinning, Chem. Phys. Lett., 240, 449–456, 1995. a
Hellwagner, J., Sharma, K., Tan, K. O., Wittmann, J. J., Meier, B. H., Madhu, P. K., and Ernst, M.: Optimizing symmetry-based recoupling sequences in solid-state NMR by pulse-transient compensation and asynchronous implementation, J. Chem. Phys., 146, 244202, https://doi.org/10.1063/1.4989542, 2017. a, b, c
Hellwagner, J., Grunwald, L., Ochsner, M., Zindel, D., Meier, B. H., and Ernst, M.: Origin of the residual line width under frequency-switched Lee–Goldburg decoupling in MAS solid-state NMR, Magn. Reson., 1, 13–25, https://doi.org/10.5194/mr-1-13-2020, 2020. a
Hohwy, M., Jakobsen, H. J., Eden, M., Levitt, M. H., and Nielsen, N. C.: Broadband dipolar recoupling in the nuclear magnetic resonance of rotating solids: A compensated C7 pulse sequence, J. Chem. Phys., 108, 2686–2694, 1998. a
Idziak, S. and Haeberlen, U.: Design and construction of a high homogeneity rf coil for solid-state multiple-pulse NMR, J. Magn. Reson., 50, 281–288, https://doi.org/10.1016/0022-2364(82)90058-0, 1982. a
Jain, M. G., Mote, K. R., Hellwagner, J., Rajalakshmi, G., Ernst, M., Madhu, P. K., and Agarwal, V.: Measuring strong one-bond dipolar couplings using REDOR in magic-angle spinning solid-state NMR, J. Chem. Phys., 150, 134201, https://doi.org/10.1063/1.5088100, 2019. a, b
Jaroniec, C. P., Tounge, B. A., Rienstra, C. M., Herzfeld, J., and Griffin, R. G.: Recoupling of heteronuclear dipolar interactions with rotational-echo double-resonance at high magic-angle spinning frequencies, J. Magn. Reson., 146, 132–139, 2000. a
Jia, L., Liang, S., Sackett, K., Xie, L., Ghosh, U., and Weliky, D. P.: REDOR solid-state NMR as a probe of the membrane locations of membrane-associated peptides and proteins, J. Magn. Reson., 253, 154–165, https://doi.org/10.1016/j.jmr.2014.12.020, 2015. a
Lee, Y. K., Kurur, N. D., Helmle, M., Johannessen, O. G., Nielsen, N. C., and Levitt, M. H.: Efficient dipolar recoupling in the NMR of rotating solids. A sevenfold symmetric radiofrequency pulse sequence, Chem. Phys. Lett., 242, 304–309, 1995. a
Levitt, M. H.: Symmetry-Based Pulse Sequences in Magic-Angle Spinning Solid-State NMR, in: eMagRes, edited by: Harris, R. K. and Wasylishen, R. L., John Wiley, Chichester, https://doi.org/10.1002/9780470034590.emrstm0551, 2007. a
Li, C., Mo, Y., Hu, J., Chekmenev, E., Tian, C., Gao, F. P., Fu, R., Gor'kov, P., and Brey, W. Cross, T. A.: Analysis of RF heating and sample stability in aligned static solid-state NMR spectroscopy, J. Magn. Reson., 180, 51–57, https://doi.org/10.1016/j.jmr.2006.01.013, 2006. a
Lindon, J. C., Beckonert, O. P., Holmes, E., and Nicholson, J. K.: High-resolution magic angle spinning NMR spectroscopy: Application to biomedical studies, Prog. Nucl. Mag. Res. Sp., 2, 79–100, 2009. a
Michal, C. A. and Jelinski, L. W.: REDOR 3D: Heteronuclear Distance Measurements in Uniformly Labeled and Natural Abundance Solids, J. Am. Chem. Soc., 119, 9059–9060, https://doi.org/10.1021/ja9711730, 1997. a
Mithu, V. S., Tan, K. O., and Madhu, P. K.: Selective inversion of 1H resonances in solid-state nuclear magnetic resonance: Use of double-DANTE pulse sequence, J. Magn. Reson., 237, 11–16, 2013. a
Nishimura, K., Fu, R., and Cross, T. A.: The Effect of RF Inhomogeneity on Heteronuclear Dipolar Recoupling in Solid State NMR: Practical Performance of SFAM and REDOR, J. Magn. Reson., 152, 227–233, https://doi.org/10.1006/jmre.2001.2410, 2001. a, b, c
Odedra, S. and Wimperis, S.: Imaging of the B1 distribution and background signal in a MAS NMR probehead using inhomogeneous B0 and B1 fields, J. Magn. Reson., 231, 95–99, https://doi.org/10.1016/j.jmr.2013.04.002, 2013. a
Paulson, E. K., Martin, R. W., and Zilm, K. W.: Cross polarization, radio frequency field homogeneity, and circuit balancing in high field solid state NMR probes, J. Magn. Reson., 171, 314–323, https://doi.org/10.1016/j.jmr.2004.09.009, 2004. a
Piotto, M., Bourdonneau, M., Furrer, J., Bianco, A., Raya, J., and Elbayed, K.: Destruction of magnetization during TOCSY experiments performed under magic angle spinning: effect of radial B1 inhomogeneities, J. Magn. Reson., 149, 114–118, https://doi.org/10.1006/jmre.2001.2287, 2001. a
Privalov, A. F., Dvinskikh, S. V., and Vieth, H.-M.: Coil Design for Large-Volume High-B1 Homogeneity for Solid-State NMR Applications, J. Magn. Reson. Ser. A, 123, 157–160, https://doi.org/10.1006/jmra.1996.0229, 1996. a
Purusottam, R. N., Bodenhausen, G., and Tekely, P.: Effects of inherent rf field inhomogeneity on heteronuclear decoupling in solid-state NMR, Chem. Phys. Lett., 635, 157–162, https://doi.org/10.1016/j.cplett.2015.06.051, 2015. a
Rienstra, C. M., Tucker-Kellogg, L., Jaroniec, C. P., Hohwy, M., Reif, B., McMahon, M. T., Tidor, B., Lozano-Pérez, T., and Griffin, R. G.: De novo determination of peptide structure with solid-state magic-angle spinning NMR spectroscopy, P. Natl. Acad. Sci. USA, 99, 10260–10265, https://doi.org/10.1073/pnas.152346599, 2002. a
Schanda, P., Meier, B. H., and Ernst, M.: Quantitative analysis of protein backbone dynamics in microcrystalline ubiquitin by solid-state NMR spectroscopy, J. Am. Chem. Soc., 132, 15957–15967, 2010. a
Schanda, P., Meier, B. H., and Ernst, M.: Accurate measurement of one-bond H–X heteronuclear dipolar couplings in MAS solid-state NMR, J. Magn. Reson., 210, 246–259, 2011. a
Smith, S., Levante, T., Meier, B., and Ernst, R.: Computer simulations in magnetic resonance. An object-oriented programming approach, J. Magn. Reson., 106, 75–105, 1994. a
States, D. J., Haberkorn, R. A., and Ruben, D. J.: A two-dimensional nuclear Overhauser experiment with pure absorption phase in four quadrants, J. Magn. Reson., 48, 286–292, 1982. a
Tan, K. O., Rajeswari, M., Madhu, P. K., and Ernst, M.: Asynchronous symmetry-based sequences for homonuclear dipolar recoupling in solid-state nuclear magnetic resonance, J. Chem. Phys., 142, 065101, https://doi.org/10.1063/1.4907275, 2015. a
Tan, K. O., Agarwal, V., Meier, B. H., and Ernst, M.: A generalized theoretical framework for the description of spin decoupling in solid-state MAS NMR: Offset effect on decoupling performance, J. Chem. Phys., 145, 094201, https://doi.org/10.1063/1.4961909, 2016. a, b
Tošner, Z., Purea, A., Struppe, J. O., Wegner, W., Engelke, F., Glaser, S. J., and Reif, B.: Radiofrequency fields in MAS solid state NMR probes, J. Magn. Reson., 284, 20–32, https://doi.org/10.1016/j.jmr.2017.09.002, 2017. a, b, c, d
Tošner, Z., Sarkar, R., Becker-Baldus, J., Glaubitz, C., Wegner, S., Engelke, F., Glaser, S. J., and Reif, B.: Overcoming Volume Selectivity of Dipolar Recoupling in Biological Solid-State NMR Spectroscopy, Angew. Chem. Int. Edit., 57, 14514–14518, https://doi.org/10.1002/anie.201805002, 2018. a, b, c, d, e