Chuanyu Yanab,
François Kayserc and
Reiner Dieden*a
aLuxembourg Institute of Science and Technology, Department of “Materials Research and Technology”, Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg. E-mail: reiner.dieden@list.lu; Fax: (+352) 275 885; Tel: (+352) 275 888 4576
bUniversity of Luxembourg, Faculty of Science, Technology and Communication, Avenue des Hauts-Fourneaux, L-4365 Esch-sur-Alzette, Luxembourg
cGoodyear Innovation Center Luxembourg, Avenue Gordon Smith, L-7750 Colmar-Berg, Luxembourg
First published on 17th June 2020
{1H–29Si}–1H double cross polarization inverse detection (DCPi) solid-state NMR, has recently been shown to be a powerful tool for studying molecules adsorbed on the silica surface. In this contribution, we develop an improved version (MCPi) which incorporates a block of multiple contact pulses, and quantitatively compare the sensitivities of MCPi and DCPi over a typical range of experimental parameters. The MCPi pulse sequence aims at higher sensitivity and robustness for studying samples with various relaxation characteristics. In the case of dimethyl sulfoxide (DMSO) molecules adsorbed on the silica surface, MCPi performs equally well or up to 2.5 times better than DCPi over a wide range of parameters. The applicability to and performance of MCPi on composite materials was demonstrated using a sample of polymer–silica composite, where significantly higher sensitivity could be achieved at very long total mixing times. The results also showed that both techniques are surface specific in the sense that only the groups close to the surface can be detected.
For the characterization of silica and silicate surfaces, 29Si cross polarisation under magic angle spinning solid-state nuclear magnetic resonance spectroscopy (29Si CP/MAS ssNMR) has proven to be a valuable tool.15,20–24 Although direct polarisation (DP) can also provide valuable information its use is hampered by very long spin–lattice (T1) relaxation times.25,26
If the Hartmann–Hahn matching condition ωI = ωS is met, the polarisations of two nuclear species can equilibrate through dipolar coupling. If the polarisations are different, this cross-polarisation (CP) can lead to an enhancement of the lowest population by a factor of γI/γS (γ is the gyromagnetic ratio of the isotope).
The efficiency of CP depends essentially on two parameters, the rotating frame relaxation (TI1ρ) of the reservoir species (e.g. 1H), and the cross-polarisation constant, TIS. Since TIS can be very different for different I–S pair in the same sample, a single setting of CP parameters is not optimal for all S nuclei of the sample in most cases.
As improvements on the simple CP experiment, schemes involving multiple CP contacts have been proposed, mostly aiming at generating a more even signal enhancement and quantitative response.27–29
While the theoretical maximum enhancement (i.e. γ1H/γ29Si) is rarely achieved by one single CP step, it has been shown that multiple-step CP (multiCP) can significantly enhance the magnetisation for species with unequal cross-relaxation characteristics.28
In multiCP, the X nucleus magnetization can be pushed closer and closer to the maximum enhancement with increasing number of contacts, eventually achieving a homogeneous enhancement for all X-nuclei via multiple CP steps. The multiCP experiment was initially designed to obtain quantitative X nucleus response.28–30 Nevertheless, there are recent reports of using a multiCP block for the purpose of sensitivity enhancement.31,32
This is because magnetization transfer via CP relies on through-space dipolar couplings, the strength of which is proportional to 1/r3 (r is the distance between the dipolar coupled spin pair). Hence, the information derived from CP ssNMR is usually limited to a region of a few nanometers.38 The dipolar coupled protons and silicon atoms have been investigated by 2D 1H–29Si heteronuclear correlation (HETCOR) on the silicon39 and silica19,38,40–42 surfaces.
Similarly, X-filtered double-cross polarisation experiments have been reported to specifically record spectra of X-nuclei of small molecules in the indirect dimension of a 2D experiment. This technique consists of two distinct CP periods during the first of which magnetisation is transferred from 1H to the X nucleus, from where it is then back-transferred to proton via the second CP block.43–48
A 1D version of this experiment, {1H–29Si}–1H double cross polarization inverse detection (DCPi), has recently been applied to the study of molecules adsorbed on the silica and silicate surfaces.38,49–52 The {1H–29Si}–1H DCPi experiment38 consists of three steps. First, a 1H → 29Si CP period (tcp1) ensures the surface 29Si magnetization is established. Second, a train of 90° excitation pulses is used to saturate the 1H magnetization so that no transverse magnetization, due to direct excitation during the second spin lock period (tcp2), could be detected in the final spectrum. During the proton saturation period, the 29Si magnetization is stored on the longitudinal axis. In the third step, the 29Si magnetisation is rotated back onto the transverse plane by a 90° pulse and then, during the second cross-polarization period (tcp2), transferred back only to those 1H spins that are within reach, i.e. close to the surface. Because of its dependence on the 1H–29Si dipolar couplings, the intensity of the 1H signal depends on the inter-nuclear distance, although proton–proton spin diffusion can also interfere, which can be minimized by using short contact times (tcp2). By applying this approach one can tell whether a molecule/functional group is on (or close to) the silica surface.
This {1H–29Si}–1H DCPi technique could provide information complementary to other sensitivity-enhanced techniques (X-nucleus detected) such as the Carr–Purcell–Meiboom–Gill (CPMG)34,40,41,53,54 and dynamic nuclear polarization techniques (DNP).33,34,55–57 In terms of silica/silicate materials, the former can reveal the local proton environments on the surface, while the latter can provide information about the silicon groups of silica and silicates.
In addition to the applications reported so far, the {1H–29Si}–1H DCPi technique also has the potential to study competitive adsorption phenomena on the surface of silica and silicate, which is an important and challenging topic.58 These ubiquitous phenomena can occur, for example, among small molecules,59–62 among polymers63,64 and between polymer and surfactant.65,66
In order to further improve the sensitivity of this surface-specific technique for composite materials with a rather low silica content or nanoparticles with a relatively low surface loading,17,56,67,68 where even slight sensitivity improvement could reduce experiment time significantly, we report here a modified version of DCPi. This new version, multiple contact cross polarization inverse detection ssNMR (MCPi), uses a multiCP block to enhance the polarisation build-up during the initial CP period. Its sensitivity is quantitatively compared to that of DCPi over a typical range of conditions.
One of the reasons for replacing the first CP step of DCPi with a block of multiple CP steps is because it is expected to further push the 29Si magnetization towards the maximum level (γ1H/γ29Si).
Another reason for inserting the multiCP block is its robustness, i.e. its capability of achieving strong enhancements for spins with different cross-polarisation kinetics.28 Although it has always been assumed that TI1ρ > TIS, the situation can be the opposite in some cases, i.e. TI1ρ < TIS.69 To cover a wide range of relaxation characteristics, e.g. as for the different types of hydroxyl groups at the silica surface, without the need of putting too much effort on measuring the TI1ρ and TIS for every sample, multiCP can safely enhance the magnetization without losing it due to TI1ρ.
The solid samples were loaded into a 4 mm zirconium oxide (ZrO2) rotor with a Kel-F cap, while the liquid samples were loaded in a Kel-F HR-MAS insert, which was then inserted into the ZrO2 rotor. The rotors were spun at a MAS speed of 14 kHz (unless noted otherwise) at 303 K in a double resonance probe of a Bruker Avance III HD 600 MHz spectrometer, controlled by TopSpin 3.5pl7 software. The chemical shifts reported here were externally referenced vs. pure TMS by setting the methylene carbon of adamantane to 38.48 ppm. The 1H and 29Si 90° pulse lengths were 2.9 and 5.0 microseconds, respectively. In cross polarization experiments, the spin locking strength of 29Si was set to 50 kHz while that of 1H was a shaped ramp from 34 to 66 kHz. The number of scans (NS) for the {1H–29Si}–1H DCPi/MCPi and {1H}–29Si multiCP were 2048, the recycle delay (RD) was 4.0 seconds and acquisition time was 50 ms, unless noted otherwise. For the saturation train in DCPi and MCPi, we used 64 saturation pulses separated by a delay (τsat) of 1 ms, unless noted otherwise. In the MCPi experiments of the polymer–silica compound, the saturation scheme of Fig. 1B was replaced with two long low-power pulses (pulse length = 500 ms at a power level of 1 watts) with a 90°-phase shift, because even a train of 512 saturation pulses (inter-pulse delay τsat = 1 ms) was found to be insufficient for obtaining a clean background. For 1H–1H RFDR, the following acquisition parameters were used: RD = 4.0 seconds, mixing time = 18 ms, 128 increments on the indirect dimension and NS = 16. For the 29Si DP experiment at natural abundance, note that the T1 of 29Si is so long that our attempt to measure it in a practical timeframe was not successful. Thus, a RD of 120 s that is quite common in literature,16 was used here.
Fig. 1 The diagrams of two solid-state NMR pulse sequences. (A) {1H–29Si}–1H double cross polarization inverse detection (DCPi). (B) {1H–29Si}–1H multiple cross polarization inverse detection (MCPi). |
The signal-to-noise ratio (S/N) was calculated in TopSpin using the Bruker command “SINO” with the signal and noise ranges defined in Tables S1 and S2.† To make fair comparisons, the {1H –29Si}–1H DCPi and MCPi spectra were processed using the same processing parameters (512 time domain points, without apodization and fixed regions for calculating S/N ratio).
Fig. 2 (A) Schematic presentation of displacement of physically adsorbed water by dimethyl sulfoxide molecules on the silica nanoparticle surface. (B) 1H one-pulse spectrum of DMSO-adsorbed SiNP (0.20 g SiNP + 0.11 g DMSO-h6) obtained by 1H one-pulse experiment (recycle delay = 4 seconds). The assignment of chemical shifts is as follow: ● δ(HB-silanol) = 7.50 ppm; ◆ δ(HB-water) = 4.63 ppm; ★ δ(HB-DMSO) = 3.16 ppm. (C) The 1H–1H radio-frequency-driven-recoupling (RFDR) correlation spectrum of DMSO-adsorbed SiNP (same sample as shown in (B)), using the pulse sequence of Fig. S1D† (mixing time = 18 ms). Note: three diagonal peaks in the region of 2 to 0 ppm might be assigned to isolated silanols and a very weak peak at 0.47 ppm that cannot be assigned.15 |
Determination of the proton spin–lattice relaxation constant (TH1, Fig. S2†) showed that the slowest relaxing group is the DMSO protons (TH1 = 0.35 seconds). Based on the finding in Schmidt-Rohr's report,28 τrel (Fig. 1B) in the multiCP block was therefore set to 0.7 seconds (2 × TH1) in order to allow 1H to repolarize to 95% of the thermal equilibrium. For obtaining optimal sensitivity, a relaxation delay of RD = 0.44 s (i.e. 1.26 × TH1) was initially chosen. However, in order to ensure that all protons have fully relaxed in cases where the T1 (s) could not be clearly determined due to heavily overlapped 1H peaks, we chose RD = 4.0 (the second to last points of the T1 relaxation curves in Fig. S2†).
Another important parameter for running the {1H–29Si}–1H experiment is the number of CP contacts (m of Fig. 1B). To figure out the optimal m value, 29Si detected and 1H detected experiments (Fig. S1B and C†) were performed on the same sample used in Fig. 2. The optimization process is shown in Fig. 3. There are pros and cons for higher values of m. Theoretically, more contacts (higher m values) result in more magnetization transfer from 1H to 29Si and therefore more backward magnetization transfer, as was experimentally confirmed (see Fig. 3A and C). However, for qualitative analysis (as in this case), higher values of m result in longer experiment time. In addition, the closer the 29Si magnetization approaches the theoretical enhancement, the less efficient per time unit is each individual multiple CP step (Fig. 3B). Therefore, there is an optimal value for m, beyond which, there is no further marginal sensitivity gain per time unit (Fig. 3D). For the sample used in this case, the optimum was found to be m = 3, based on the maximum marginal gain of sensitivity of the silanol, water and DMSO signals (Fig. 3D).
With the key parameters being known (RD, τrel and m of Fig. 1), the sensitivities of the {1H–29Si}–1H DCPi and MCPi techniques were compared. The two experiments were performed on the same sample (Fig. S4 and S5†) and the S/N ratios were calculated (Tables S1 and S2†). As expected, one immediately noticeable difference with the 1H one-pulse is the relative intensities of the silanol, water and DMSO signals. In the one-pulse experiment (Fig. 2B), the relative peak intensities from highest to lowest is DMSO, water and silanol, whereas it is the opposite in the inverse detection experiments (Fig. 4A), which makes sense because the peak intensity of the one-pulse experiment are solely determined by concentration whereas it is modulated by inter-nuclear distance.
Fig. 4 The sensitivity comparison between {1H–29Si}–1H MCPi and DCPi. Same sample as used in Fig. 2B. (A) 1H spectra obtained using the following parameters: RD = 0.44 s, tcp1 = 1 ms, tcp2 = 5 ms for both MCPi and DCPi, and m = 3 for MCPi. (B) Sensitivity ratio of MCPi over DCPi (kMCPi/kDCPi) when RD = 0.44 s (i.e. 1.26 × TH1), tcp1 = variable and tcp2 = 5 ms, and m = 3 for MCPi. (C) Same experiments as performed for (B), except RD = 4.0 s. |
The overall sensitivity of each technique depends on the signal intensity, which was gained from one or multiple CP steps. Under three different conditions, the signals in the {1H–29Si}–1H MCPi spectra were all stronger than in DCPi spectra. Overall, {1H–29Si}–1H MCPi yields higher signal to noise ratio (S/N) than DCPi (Table S1 and Fig. S4†).
By normalizing the S/N ratios to experiment time per scan (i.e. the recycle delay + tcp1 × m + τsat × n + τrel × m + tcp2 + acquisition time), we could compare the intrinsic sensitivities of {1H–29Si}–1H MCPi and DCPi (Table S1†) using eqn (1)–(6). In a fully relaxed spin system, the S/N is proportional to the square root of the number of scans (NS) (eqn (1)), which is proportional to the total experimental time (Expt, eqn (2)). Therefore, the S/N is proportional to the square root of the total experimental time (eqn (3)). Thus, the sensitivity of a given pulse sequence can be expressed as a sensitivity factor (k) times the square root of total experimental time (eqn (4)). Under the same experimental conditions (i.e. identical sample, temperature, receiver gain, number of scans, processing parameters, etc.), the k value can be determined by normalizing the S/N to total experimental time (eqn (5)). Hence, the time-efficiency of two sequences can be compared by comparing their k values (eqn (6)).
(1) |
NS ∝ Expt | (2) |
(3) |
(4) |
(5) |
(6) |
The tcp1 dependence of the k ratios between {1H–29Si}–1H MCPi and DCPi under variable tcp1 are shown in Fig. 4B and C. In NMR experiments, the sensitivity per time unit can be maximized by setting the recycle delay to 1.26 times the longest T1 in a sample.72–74 As shown in Fig. 4B, when a short relaxation delay of RD = 0.44 s was used, the overall sensitivity of MCPi for the molecule of interest, DMSO, was better or equivalent than that of DCPi over the tcp1 range. For silanol and water peaks, the sensitivity of MCPi is similar to that of DCPi, which makes sense because these groups are located closely to the surface silicon atoms so that magnetization transfer via CP is fast and a single CP contact is probably sufficient for transferring most of the proton magnetization to silicon. However, the optimum of 1.26 × T1 is derived from an idealized situation where all the 1H T1 values in a sample can be clearly determined. Although this choice of recycle delay was applicable to the DMSO-adsorbed SiNP, that may be different in the case of composite materials where the 1H peaks are heavily overlapped and some of the peaks might be so overwhelming that the T1 of 1H of other surface groups cannot be clearly determined. In that case, one needs to resort to using a long recycle delay in order to avoid saturating the surface groups (see the next section of a polymer–silica composite).
As shown in Fig. 4C, the results of RD = 4.0 s showed that the sensitivity improvement of {1H–29Si}–1H MCPi compared to DCPi, was nearly three-fold at tcp1 = 1 ms (meaning that it would take about 9 times longer time for DCPi to achieve the same S/N level as MCPi). At longer tcp1, the improvement factor is around 1.5. Theoretically speaking, {1H–29Si}–1H MCPi would be even more efficient over the DCPi in situations where the proton spin–lattice relaxation in the rotating frame is much faster (shorter T1ρ) (i.e. long spin locking pulse is inapplicable). In fact, a contact pulse of 1 ms is quite common.28,38,50
Moreover, we further explored the dependence of signal intensity on tcp2 in the experiments of {1H–29Si}–1H MCPi (Fig. S6†). The rates of S/N evolution with tcp2 differ between DMSO, silanol and water due to the difference in the distance relative to the silica surface. This behaviour is analogous to the one reported previously for the regular DCPi.38
Fig. 5 Schematic representation of the silica surface in a polymer–silica composite, where the silica has been silanized by trimethylchlorosilane and the polymer is polyisoprene (see Experimental). |
The 29Si spectra in Fig. 6A and B show that after the reaction with an excessive amount of silane, there is a significant amount of trimethylsilyl group grafted onto the silica surface and a noticeable decrease of the silanol peak (at −102 ppm). The 1H spectra in Fig. 6C and D show that once the silanized silica was mixed with polymer, the 1H peaks of polymer are so dominant that not all of the silanized silica 1H peaks (3 to 10 ppm) can be distinguished. Thus, instead of performing a regular 1H T1 measurement, the recycle delay needed for reaching the thermal equilibrium was estimated to be >3.5 s in a set of 1H one-pulse spectra with variable recycle delays (Fig. S7†). As a result, a recycle delay of 4.0 s and therefore an estimated inter-CP block delay τrel = 0.7 s were used for the following MCPi experiments. These two key parameters were confirmed to be efficient for building up the 29Si magnetization rapidly (Fig. S8†).
The 1H spectra in Fig. 7A were obtained by MCPi at various m values. As mentioned earlier, the MCPi pulse sequence of Fig. 1B actually can also achieve the function of DCPi simply by setting the multiCP loop number, m = 1 and inter-CP delay τrel = 0, because the contribution of the first 90° pulse on 29Si is negligible due to the long T1 relaxation of 29Si (Fig. 3C).
Fig. 7 (A) 1H MCPi spectra of polymer–silica composite obtained by a slightly modified pulse sequence (Fig. S9†): RD = 4.0 s, tcp1 = tcp2 = 5 ms, τrel = 0.7 s and m = 1 (black), 3 (blue) and 6 (green). For m = 1, τrel = 0 s, this is effectively a DCPi experiment. Grey boxes are the regions where the polyisoprene signals would be expected. (B) The ratio of sensitivity factors k, between MCPi and DCPi calculated by eqn (6) using the S/N and experiment duration (Expt) values from Table S3.† |
The sensitivity factor, k (see eqn (5)), for each m value was calculated plugging in the S/N ratios and the experiment duration (Expt) of Table S3.† As can be seen in Fig. 7B, the MCPi pulse sequence achieves a higher sensitivity by increasing the number of multiCP loops, m, with total mixing times that could not normally be achieved with a single contact. This result was overall in agreement with the results of the DMSO-adsorbed silica case. As mentioned in the Introduction, the observable inter-nuclear distance has an upper limit of a few nanometers. Even though a wide range of tcp2 values have been used (i.e. tcp2 = 5 ms in Fig. 7A and tcp2 = 1, 10 & 15 ms in Fig. S10†), none of the polyisoprene 1H peaks was observed, which is probably due to the weak dipolar coupling between the polyisoprene and silica surface.
Overall, MCPi is expected to perform better than DCPi in three situations: (1) the sample (e.g. composite materials) has multiple relatively fast 1H spin–lattice relaxations in the rotating frame TH1ρ, (2) the silica surface is not (or only slightly) protonated, and (3) the overall silica content in the sample is very low. In the first situation, the built-up 29Si magnetization will be constantly drained away due to TH1ρ. This will be attenuated using multiple, but short, cross polarisation steps, which will also be better suited for a variety of TH1ρs. In the second scenario, i.e. when there are only few protons near the silica surface (e.g. in pyrogenic silica or modified silicas with few grafted or adsorbed groups), the proton–silicon dipolar coupling is weaker and the CP kinetics is slower. Therefore, the total contact time has to be very long, which can be safely achieved by a multiCP block. In the third case, when the total amount of silica in sample, e.g. a composite material, is low, therefore requiring long experimental time in order to achieve signal-to-noise ratios suitable for analysis, even a small sensitivity enhancement means significant saving of experiment time.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ra04995f |
This journal is © The Royal Society of Chemistry 2020 |