Multinuclear Solid-State Nuclear Magnetic ... - ACS Publications

Jan 6, 2012 - Herein, we present multinuclear solid-state NMR measurements for taurine, which features three challenging quadrupolar nuclei: 14N (nucl...
1 downloads 0 Views 1MB Size
Article pubs.acs.org/JPCA

Multinuclear Solid-State Nuclear Magnetic Resonance and Density Functional Theory Characterization of Interaction Tensors in Taurine Luke A. O’Dell,*,† Christopher I. Ratcliffe,† Xianqi Kong,‡ and Gang Wu‡ †

Steacie Institute for Molecular Sciences, National Research Council, 100 Sussex Drive, Ottawa, Ontario, K1A 0R6, Canada Department of Chemistry, Queen’s University, 90 Bader Lane, Kingston, Ontario, K7L 3N6, Canada



S Supporting Information *

ABSTRACT: A variety of experimental solid-state nuclear magnetic resonance (NMR) techniques has been used to characterize each of the elements in 2-aminoethane sulfonic acid (taurine). A combination of 15N cross-polarization magic angle spinning (CPMAS), 14N ultrawideline, and 14N overtone experiments enabled a determination of the relative orientation of the nitrogen electric field gradient and chemical shift tensors. 17O spectra recorded from an isotopically enriched taurine sample at multiple magnetic fields allowed the three nonequivalent oxygen sites to be distinguished, and NMR parameters calculated from a neutron diffraction structure using density functional theory allowed the assignment of the 17O parameters to the correct crystallographic sites. This is the first time that a complete set of 17O NMR tensors are reported for a sulfonate group. In combination with 1H and 13C MAS spectra, as well as a previously reported 33S NMR study, this provides a very broad set of NMR data for this relatively simple organic molecule, making it a potentially useful structure on which to test DFT calculation methods (particularly for the quadrupolar nuclei 14N, 17O, and 33S) or NMR crystallography approaches. of a combination of 15N and 14N methods, we obtain a full set of NMR parameters for the nitrogen site, including a determination of the relative orientation of the EFG and chemical shift (CS) tensors. Parameters for the three nonequivalent oxygen sites are also obtained from an isotopically enriched sample and assigned to the correct crystallographic sites with the aid of DFT calculations. Solid-state 1H and 13C NMR spectra are also presented, which, in combination with a previous 33S study,11 constitutes a near-complete set of NMR results for this model system.

1. INTRODUCTION Solid-state nuclear magnetic resonance (NMR) spectroscopy is one of the premier techniques for the structural characterization of materials, capable of providing detailed information on the local structure and dynamics in a broad range of systems. For the study of crystalline samples, a combination of solid-state NMR experiments and density functional theory (DFT) calculations can be particularly powerful. Gauge-including projector augmented wave (GIPAW) methods1 that take into account the infinite periodic lattice can be used to calculate NMR parameters to a high degree of accuracy, thereby allowing the exquisite sensitivity of the measured NMR parameters to the crystal or molecular structure to be fully exploited, for example to refine an X-ray diffraction structure. Such approaches are now often referred to as NMR crystallography.2 Spin-half nuclei such as 1H, 13C, or 29Si have been used to refine structures via their isotropic chemical shifts,3 chemical shift anisotropies,4 or spin diffusion properties.5 Parameters measured from quadrupolar nuclei have also been used to refine crystal structures,6−9 but due to various associated difficulties such as low sensitivities or large anisotropic interactions, such nuclei are far more seldom used in NMR crystallography. This is unfortunate, since quadrupolar nuclei are capable of providing a wealth of local structural information.10 In particular, the electric field gradient (EFG) offers a probe of the local electronic environment that is highly sensitive to the surrounding structure, particularly hydrogen bonding arrangements. Herein, we present multinuclear solid-state NMR measurements for taurine, which features three challenging quadrupolar nuclei: 14N (nuclear spin I = 1), 17O (I = 5/2), and 33S (I = 3/2). This is therefore a potentially useful system on which to test out new experimental or computational methodologies. By use © 2012 American Chemical Society

2. EXPERIMENTAL DETAILS a. Sample Details. For the solid-state 1H, 13C, 14/15N NMR experiments, a natural abundance sample of taurine (>99% pure, purchased from Sigma Aldrich) was used. For the 17O NMR experiments, 17O-enriched taurine was synthesized in the following manner. Sodium sulfite (300 mg, 2.38 mmol) was dissolved in [17O]-water (1.16 g, 41% 17O-enriched, purchased from CortecNet). The solution was kept at room temperature for 23 h. 2-Bromoethylamine hydrochloride (585 mg, 2.86 mmol) was then added to the solution. The mixture was stirred at 65 ± 3 °C in an oil bath for 13 h. The [17O]-water was recovered. The residual material was refluxed in methanol (30 mL) for 30 min, cooled to room temperature, collected (filtration), washed with ethanol (2 × 10 mL), and dried at 80 °C. The crude product was dissolved in 1 mL of 40% hydrobromic acid and filtered through sintered glass. Ethanol (35 mL) was added to the filtrate to precipitate out the product. The solid material was then washed Received: November 10, 2011 Revised: January 6, 2012 Published: January 6, 2012 1008

dx.doi.org/10.1021/jp210844t | J. Phys. Chem. A 2012, 116, 1008−1014

The Journal of Physical Chemistry A

Article

with ethanol (4 × 10 mL) and dried at 80 °C, giving 220 mg of white solid (74% yield). The 17O enrichment of the product was determined to be 23% by solution 17O NMR. b. Solid-State NMR Experiments. The 1H, 13C, and 15N NMR spectra were recorded on an 11.7 T Bruker Avance III spectrometer using Varian double-resonance probes with a 3.2 mm rotor diameter and 20 kHz MAS rate for the 1H experiment and a 4.0 mm rotor diameter and 10 kHz MAS rate for the 13C and 15N experiments. The 1H spectrum was acquired using a windowed DUMBO pulse sequence12 for homonuclear decoupling, with a 2.8 μs π/2 pulse, a 24 μs DUMBO pulse length (ν1 = 71 kHz), a recycle delay of 5 and 4 scans acquired. For the 13C experiment, a 2 μs excitation pulse was applied with continuous-wave 1 H decoupling of ∼40 kHz power, a recycle delay of 2 min and 32 scans acquired. The 15N spectrum was acquired with crosspolarization from the 1H nuclei using a 1H π/2 pulse of 4 μs, 1.5 ms constant-amplitude contact pulses (ν1 = 55 kHz), a recycle delay of 2.5, and 1000 scans acquired. The 1H and 13C spectra are referenced to tetramethylsilane and the 15N spectrum to solid NH4Cl. The 14N overtone NMR spectrum was recorded on a 21.1 T Bruker Avance II spectrometer (Ultrahigh-field NMR Facility for Solids, Ottawa, Canada) using a home-built double-resonance static probe with a 7.0 mm inner-diameter coil size. A 25-μs on-resonance excitation pulse (ν1 = 80 kHz as measured on the

2

H nucleus, which is close in frequency to the 14N overtone transition) was used with continuous-wave 1H decoupling of ∼80 kHz power, a recycle delay of 1 s, and 50,000 scans acquired. The spectrum is referenced such that 0 kHz corresponds to exactly twice the 14N fundamental (single-quantum) frequency of solid NH4Cl. The 17O NMR spectra were recorded at 11.7 and 21.1 T on Bruker Avance and Bruker Avance II spectrometers, respectively. The static and 1D MAS spectra were acquired using a spin−echo pulse sequence (π/2−τ−π/2−acquire). At 11.75 T, both nonspinning and MAS experiments were performed on a Bruker 4.0 mm double resonance MAS probe, the latter with a MAS rate of 12.0 kHz and 512 scans acquired with a recycle delay of 1 s. To obtain a high-quality static spectrum at this field, 8854 scans were acquired with a recycle delay of 5 s. At 21.1 T, nonspinning experiments were performed on a home-built static probe with the sample packed in a 5.0 mm Teflon tube, while MAS and 3QMAS experiments were performed on a Bruker 4.0 mm MAS probe with the sample packed in a Si3N4 rotor spinning at 12.5 kHz. For both MAS and static spectra, a spin−echo sequence was used with a pulse delay of 50 μs, 256 scans acquired and a recycle delay of 20 s. In the 3QMAS experiment, the t1 increment was synchronized with the sample spinning and a total of 48 t1 increments were collected with 192 scans each and a recycle delay of 15 s. Spectral processing was carried out using the NUTS software (Acorn NMR). The 17O spectra were simulated using Dmfit.13 The 14N overtone powder pattern was simulated using RMNSim. 14 c. DFT Calculations. DFT calculations were carried out on a neutron diffraction structure15 (CCSD code TAURIN03) using the CASTEP software16 and the Materials Studio 4.3 program (Accelrys) running on a Linux server with 8 processing cores and 32 GB of RAM. Perdew, Burke, and Ernzerhof (PBE) functionals were used with a plane wave basis set cutoff of 610 eV and a 3 × 1 × 2 Monkhorst-Pack k-space grid. Electric field gradient and chemical shielding tensors were

Figure 1. The molecular structure of taurine, taken directly from the crystal structure15 with individual sites labeled.

Figure 2. (a) 1H NMR spectrum obtained from taurine at 11.7 T and 20 kHz MAS, using a windowed DUMBO homonuclear decoupling pulse sequence. The peak at 7.6 ppm corresponds to the three distinct NH3 protons, while the peak at 3.3 ppm arises from the four CH2 protons. (b) 13C NMR spectrum obtained from the same sample at 11.7 T and 10 kHz MAS. The two carbon sites are resolved with the peak at 46.5 ppm corresponding to the carbon bonded directly to the sulfur and the peak at 35.8 ppm arising from the carbon bonded to the nitrogen. 1009

dx.doi.org/10.1021/jp210844t | J. Phys. Chem. A 2012, 116, 1008−1014

The Journal of Physical Chemistry A

Article

calculated for all nuclei. Because of the relatively high degree of accuracy for this crystal structure, no geometry optimization was performed. DFT optimizations of the proton positions in this structure have previously been shown to yield calculated 14N and 33 S NMR parameters that agree less well with experimental values than those obtained from the unoptimized structure.11 Calculated isotropic shielding values σcalc for 1H, 13C, and 17O were converted to isotropic chemical shifts δcalc using the expression δcalc = σref − σcalc, with values for σref previously reported

as σref(1H) = 30.8 ppm, σref(13C) = 170.0 ppm, and σref(17O) = 265 ppm.9 For 14N, the relation δcalc = 201.4 − 1.05σcalc was used, having been determined in a previous study.17 To attempt to account for thermal vibrations in the DFT calculations of the NMR parameters, a molecular dynamics simulation was also carried out using the CASTEP software16 following a procedure recently reported.18 This used a cutoff energy of 550 eV and a 3 × 1 × 2 Monkhorst Pack grid. The temperature was set at 290 K using a Nose thermostat in the NVT ensemble. The time step was set to 1 fs with a total time of 6 ps.

Table 1. 1H and 13C Isotropic Chemical Shift Values for the Seven Crystallographically Distinct Hydrogen Sites and Two Distinct Carbon Sites in Taurine, Calculated from the Crystal Structure Using DFTa

3. RESULTS AND DISCUSSION 3.1. 1H and 13C NMR Spectra. The zwitterionic molecular structure of 2-aminoethane sulfonic acid (taurine) in its crystalline form is shown in Figure 1. The crystal structure is monoclinic with four molecules per unit cell and a cell volume of 486 Å3.15 The asymmetric unit contains a single molecule, so all atoms in the molecule are crystallographically distinct, resulting in seven hydrogen sites, two carbon sites, three oxygen sites, and a single nitrogen and sulfur site. A 1H MAS NMR spectrum obtained from taurine at 11.7 T is shown in Figure 2a. Two peaks are observed and are assigned to the NH3 protons (7.6 ppm) and CH2 protons (3.3 ppm) on the basis of the isotropic chemical shifts calculated from the crystal structure (Table 1). The calculated shifts for the four crystallographically distinct CH2 protons are clustered within 0.5 ppm of each other, and the experimental 1H NMR peaks are not individually resolved even with the use of the DUMBO

site

calculated δiso/ppm

experimental δiso/ppm

H1 H2 H3 H4 H5 H6 H7 C1 C2

8.1 5.8 6.4 3.2 3.7 3.4 3.4 43.9 32.9

7.6(1)

3.3(1)

46.5(1) 35.8(1)

a

Because of peak overlap, experimental 1H isotropic chemical shift values correspond to the peak maxima for the two groups of sites indicated.

Figure 3. (a) 15N CPMAS NMR spectrum of taurine at 11.7 T, (b) 14N overtone spectrum obtained from a static powder sample at 21.1 T, (c) simulation of the overtone powder pattern made using the parameters given in Table 2 (including the quadrupolar interaction and the chemical shift anisotropy), and (d) simulation neglecting the effects of the chemical shift anisotropy. 1010

dx.doi.org/10.1021/jp210844t | J. Phys. Chem. A 2012, 116, 1008−1014

The Journal of Physical Chemistry A

Article

Table 2. Experimental (Bold) and Calculated NMR Parameters for the Nitrogen and Three Nonequivalent Oxygen Sites in Taurinea 14

N exptl 14 N calcd 17 O (O1) exptl 17 O (O1) calcd 17 O (O2) exptl 17 O (O2) calcd 17 O (O3) exptl 17 O (O3) calcd

δiso/ppm

Ω/ppm

κ

CQ/MHz

ηQ

α/deg

β/deg

γ/deg

−6.0(2) −9.1 170.4(2) 178.2 179.0(2) 187.3 187.4(2) 190.8

18(5) 12 88(10) 115 88(10) 88 78(10) 65

−0.16(30) 0.05 0.75(10) 0.36 0.68(10) 0.53 0.69(10) 0.67

1.19(1) 1.29 6.70(2) −7.17 6.65(2) −6.97 6.80(2) −6.76

0.18(1) 0.36 0.14(4) 0.08 0.16(4) 0.16 0.05(4) 0.04

270(50) 251 90(10) 266 0(10) 123 0(10) 10

77(3) 86 −20(2) −15 0(2) 12 −5(2) −11

90(15) 352 90(10) 250 0(10) 12 0(10) 17

Calculated isotropic shielding values were converted to chemical shifts δiso based on previously reported empirical relations for nitrogen17 and oxygen.9 Also shown are the span (Ω) and skew (κ) of the CS tensor, the quadrupolar coupling constant (CQ) and asymmetry parameter (ηQ), and the Euler angles (α, β, and γ) describing the relative orientations of the CS and EFG tensors. Definitions of these parameters are standard and can be found elsewhere.17Parameters for the 33S nucleus are reported in reference 11. a

homonuclear decoupling pulse sequence. The three NH3 protons are also unresolved in the experimental spectrum, despite the calculated isotropic shift values being spread over a range of 2.3 ppm. The convergence of the NH3 proton chemical shifts to a single, relatively sharp peak is due to the rotational dynamics of this amine group which are unaccounted for in the DFT calculations. There have been a large number of NMR studies of NH3 group rotations in crystalline amino acids, particularly 1H relaxation measurements19 and 2 H line shape studies,20 and typically, jump rates for the NH3 groups in these systems are on the order of 106 s−1 at room temperature.20−22 The molecular dynamics simulation carried out on this structure showed considerable librational motion of the NH3 group, but interchange of the hydrogen sites was not observed due to the extremely short time scale of the simulation (6 ps). To account for these fast thermal motions in our calculations, a new set of NMR parameters for all nuclei were obtained from the MD simulation by averaging the tensor elements calculated from 46 individual frames spaced 50 fs apart. This approach was successfully used previously to improve calculated 13C chemical shift anisotropies;18 however, in this case it did not result in an improved agreement with our experimental values for this system, and these parameters are therefore not shown. The two carbon sites are well-resolved in the 13C NMR spectrum with isotropic chemical shifts of 46.5 and 35.8 ppm (Figure 3), and on the basis of the calculated values (Table 1) these can be assigned to sites C1 and C2, respectively. This assignment is consistent with the observation of a much shorter spin−lattice relaxation time for the peak at 35.8 ppm. The C2 carbon is directly bonded to a nitrogen site; therefore all C2 carbons experience 13C−14/15N dipolar couplings, which will contribute to relaxation (though we note that no residual dipolar splittings were observed for the C2 peak, being too small to resolve at this field). In contrast, C1 is bonded to a sulfur site with a corresponding low abundance of the NMRactive 33S isotope (0.76%), thus contributions to spin−lattice relaxation from 13C−33S couplings will be insignificant. We note that the calculated separation of these peaks (11.0 ppm), a quantity insensitive to errors in the value of σref(13C) used, matches remarkably well with the experimental value (10.7 ppm). 3.2. 14N and 15N NMR Spectra. An accurate isotropic chemical shift for the nitrogen site (−6.0 ppm) was obtained from the 15N CPMAS spectrum (Figure 3a). Along with the quadrupolar coupling constant (CQ) and asymmetry parameter (ηQ) that were recently determined11 using ultrawideline 14N

Figure 4. Calculated orientations of the principal components of the electric field gradient (purple, |V33| ≥ |V22| ≥ |V11|) and chemical shift tensors (green, δ11 ≥ δ22 ≥ δ33) of (a) nitrogen, (b) sulfur, and (c) the three oxygen sites in the molecular frame of taurine.

techniques,23 this provides constraints for the fitting of the 14N overtone powder pattern obtained at 21.1 T (Figure 3b). 14N overtone NMR spectra, in which the nuclei are irradiated at approximately twice the Larmor frequency in order to directly excite and observe the Δm = 2 “overtone” transition,24 are insensitive to the first-order quadrupolar interaction, which for 14 N is typically on the order of several MHz.25 The 14N overtone powder pattern widths are therefore typically on the order of kHz, and perturbations due to both the second-order quadrupolar interaction and chemical shift anisotropy can be resolved, particularly at high field strengths where the effects of the CSA are more pronounced. Such experiments therefore provide a potentially straightforward way of determining the relative orientations of the nitrogen EFG and CS tensors from a powder sample. The simulated 14N overtone powder pattern 1011

dx.doi.org/10.1021/jp210844t | J. Phys. Chem. A 2012, 116, 1008−1014

The Journal of Physical Chemistry A

Article

Figure 5. Experimental MAS (a and b) and static (c and d) 17O NMR spectra (black) obtained from an isotopically enriched sample of taurine at the magnetic field strengths shown. Spectra (a) and (b) were recorded at 12.5 and 12.0 kHz MAS, respectively. Asterisks denote spinning sidebands. Red traces show fitted simulations, with the individual simulated lineshapes for sites O1, O2, and O3 shown below in green, purple, and blue, respectively.

between the quadrupolar and shielding interactions, but to the best of our knowledge, such an approach has so far only been applied to cases where the two tensors are approximately coincident (e.g., nitrate or tetraalkylammonium ions).26 The fitting of the static overtone powder pattern using constraints from 15N and ultrawideline 14N experiments represents a more straightforward approach that is analogous to the way Euler angles are commonly determined for half-integer quadrupolar nuclei, where static central-transition powder patterns are fitted using constraints from MAS spectra (vide infra). DFT calculations allow the (calculated) orientations of the EFG and CS tensors to be visualized in the molecular frame, and this is illustrated for the nitrogen site in Figure 4a. The largest principal component of the EFG tensor (V33) and the least shielded component of the CS tensor (δ11) are predicted to be aligned approximately parallel to the N−C bond. These orientations are consistent with previous DFT studies of amino groups.27,28 The correct CS tensor orientation for the 33S nucleus, which was mislabeled in a previous publication,11 is also shown in Figure 4b for completeness. The least shielded component of

shown in Figure 3c was fitted to the experimental spectrum by manually adjusting the CS tensor parameters (the span Ω and skew κ) and Euler angles (α, β, and γ) while keeping δiso, CQ, and ηQ fixed to their experimental values. The resultant, complete set of 14N interaction parameters are presented in Table 2. Experimental uncertainties are larger for certain parameters than others, indicating that these parameters had a relatively minor effect on the shape of the simulated overtone powder pattern (e.g., the pattern shape was extremely sensitive to the value of the Euler angle β, but far less so to α). Overall, the agreement between the experimental and calculated 14N parameters is reasonably good. The overestimation of the calculated CQ by 100 kHz is of the same magnitude observed in a previous 14 N study.17 Some additional example simulations with altered interaction parameters are provided as Supporting Information to illustrate the extent to which the shape of the static overtone powder pattern changes as the parameters are adjusted. The relative orientation of nitrogen EFG and CS tensors can also in principle be ascertained by fitting 14N MAS spinning sideband manifolds taking into account second-order cross-terms 1012

dx.doi.org/10.1021/jp210844t | J. Phys. Chem. A 2012, 116, 1008−1014

The Journal of Physical Chemistry A

Article

shown in Table 2) peak positions confirms the accuracy of the reported parameters. For the sulfonate group in taurine, the observed values of CQ (6.65−6.80 MHz) and δiso (170−187 ppm) for 17O are significantly larger than those reported for p-toluenesulfonic acid monohydrate (4.79 MHz and 140 ppm).31 Clearly, this discrepancy is due to the difference in hydrogen bonding between the two compounds. In p-toluenesulfonic acid monohydrate, each of the three sulfonate oxygen atoms is involved in a very strong O···H−O hydrogen bond with the O···O distance being 2.520, 2.525, and 2.538 Å.32 In taurine, however, the sulfonate oxygen atoms are involved in weaker O···H−N hydrogen bonding with the O···N distances ranging between 2.789 and 3.017 Å.15 Interestingly, among the three oxygen atoms in taurine, the isotropic chemical shift appears to exhibit stronger sensitivity toward the strength of hydrogen bonding than the quadrupole coupling constant.

4. SUMMARY Experimental NMR parameters have been obtained for 1H, 13C, 14/15 N, and 17O in crystalline taurine. The use of a combination of 15N CPMAS, 14N static ultrawideline, and 14N overtone NMR spectroscopy provided a relatively straightforward method for determining the relative orientation of the nitrogen EFG and CS tensors. 17O NMR parameters were also measured for the three nonequivalent oxygen sites using multiple-field static, MAS, and 3QMAS experiments. Along with a previous 33S study,11 this constitutes a large and varied set of data on this relatively simple crystalline structure, making it a potentially useful system on which to test DFT calculation methods or NMR crystallography approaches, particularly for the quadrupolar nuclei. NMR parameters calculated from an unoptimized neutron diffraction structure showed a good overall agreement with experimental values, though there appears to be some room for improvement. This may come from an improved crystal structure, more accurate computational methods (e.g., improved pseudopotentials or methods that account for dynamics), or a combination of the two. Ultimately, it is anticipated that more accurate DFT calculations for quadrupolar nuclei will allow a refinement of crystal structures such as this using the broad range of experimental NMR parameters available from such nuclei.

Figure 6. (a) Experimental 17O 3QMAS spectrum obtained from an isotopically enriched sample of taurine at 21.1 T. (b) Comparison between calculated and observed peak positions along the isotropic dimension. The uncertainty in the experimental values is estimated to be ±0.5 ppm.

the 33S CS tensor (δ11) lies approximately parallel to the S−C bond. 3.3. 17O NMR Spectra. Taurine has three inequivalent oxygen sites in similar structural environments, but despite considerable overlap of the 17O powder patterns, these sites could be distinguished on the basis of multiple-field MAS and static NMR spectra (Figure 5). EFG parameters and isotropic chemical shifts were measured from the MAS spectra and then fixed during fits of the static spectra, allowing the CS tensors and Euler angles to be quantified. The calculated 17O parameters for these sites show reasonably good agreement with the experimental values, allowing the assignment of the NMR parameters to the correct crystallographic sites. The EFG and CS tensor orientations for all three oxygen sites are shown in Figure 4c. In all three cases, the largest principal component of the EFG tensor (V33) and the most shielded component of the CS tensor (δ33) are aligned approximately parallel to the O−S bond. This is the first time that a complete set of 17O EFG and CS tensors have been determined for a SO3− group.29 To further confirm the accuracy of the 17O NMR parameters, we obtained a 3QMAS spectrum from the enriched taurine. As seen from Figure 6, three well-resolved peaks are observed in the isotropic dimension. For 17O, the peak position along the isotropic dimension of a sheared 3QMAS spectrum is determined by Fiso (ppm) = δiso (ppm) + (3/850)(CQ/ν0)2(1 + ηQ2/3) × 10−6 (ppm).30 As seen from Figure 6, the excellent agreement between the observed and predicted (using the 17O NMR parameters



ASSOCIATED CONTENT

S Supporting Information *

Further examples of simulated 14N static overtone powder patterns. This material is available free of charge via the Internet at http://pubs.acs.org.

■ ■

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected].

ACKNOWLEDGMENTS Prof. Philip Grandinetti is thanked for providing the RMNSim software used to simulate the 14N overtone spectra. Dr. Andreas Brinkmann is acknowledged for assistance in setting up the 1H DUMBO experiment and for helpful discussions. Access to the 21.1 T NMR instrument was provided by the National UltrahighField NMR Facility for Solids, Ottawa (www.nmr900.ca). We are grateful to Drs. Eric Ye and Victor Terskikh for assistance in recording the 17O NMR spectra at 21.1 T. G.W. thanks the Natural Sciences and Engineering Research Council (NSERC) of Canada for support. 1013

dx.doi.org/10.1021/jp210844t | J. Phys. Chem. A 2012, 116, 1008−1014

The Journal of Physical Chemistry A



Article

REFERENCES

(1) Pickard, C. J.; Mauri, F. Phys. Rev. B 2001, 63, 245101. (2) Harris, R. K.; Wasylishen, R. E.; Duer, M. J. NMR Crystallography; Wiley: New York, 2009. (3) Harris, R. K.; Hodgkinson, P.; Zorin, V.; Dumez, J.-N.; ElenaHerrmann, B.; Emsley, L.; Salager, E.; Stein, R. S. Magn. Reson. Chem. 2010, 48, 103−112. (4) Brouwer, D. H.; Moudrakovski, I. L.; Darton, R. J.; Morris, R. E. Magn. Reson. Chem. 2010, 48, 113−121. (5) Salager, E.; Stein, R. S.; Pickard, C. J.; Elena, B.; Emsley, L. Phys. Chem. Chem. Phys. 2009, 11, 2610−2621. (6) Hansen, M. R.; Madsen, G. K. H.; Jakobsen, H. J.; Skibsted, J. J. Phys. Chem. A 2005, 109, 1989−1997. (7) Zhou, B.; Giavani, T.; Bildsøe, H.; Skibsted, J.; Jakobsen, H. J. Chem. Phys. Lett. 2005, 402, 133−137. (8) Widdifield, C. M.; Bryce, D. L. Phys. Chem. Chem. Phys. 2009, 11, 7120−7122. (9) Wong, A.; Smith, M. E.; Terskikh, V.; Wu, G. Can. J. Chem. 2011, 89, 1087−1094. (10) MacKenzie, K. J. D.; Smith, M. E. Multinuclear Solid-state NMR of Inorganic Materials: Pergamon, 2002. (11) O’Dell, L. A.; Ratcliffe, C. I. J. Phys. Chem. A 2011, 115, 747− 752. (12) Sakellariou, D.; Lesage, A.; Hodgkinson, P.; Emsley, L. Chem. Phys. Lett. 2000, 319, 253−260. (13) Massiot, D.; Fayon, F.; Capron, M.; King, I.; Le Calve, S.; Alonso, B.; Durand, J.-O.; Bujoli, B.; Gan, Z.; Hoatson, G. Magn. Reson. Chem. 2002, 40, 70−76. (14) Trease, N. M.; Grandinetti, P. J. J. Chem. Phys. 2008, 128, 052318. (15) Briant, C. E.; Jones, D. W. J. Chem. Crystallogr. 1997, 27, 481− 483. (16) Clark, S. J.; Segall, M. D.; Pickard, C. J.; Hasnip, P. J.; Probert, M. J.; Refson, K.; Payne, M. C. Z. Kristallogr. 2005, 220, 567−570. (17) O’Dell, L. A.; Schurko, R. W.; Harris, K. J.; Autschbach, J.; Ratcliffe, C. I. J. Am. Chem. Soc. 2011, 133, 527−546. (18) Dumez, J.-N.; Pickard, C. J. J. Chem. Phys. 2009, 130, 104701. (19) Andrew, E. R.; Hinshaw, W. S.; Hutchins, M. G.; Canepa, P. C. Chem. Phys. Lett. 1974, 26, 50−52. (20) Gu, Z.; Ebisawa, K.; McDermott, A. Solid State Nucl. Magn. Reson. 1996, 7, 161−172. (21) Kristensen, J. H.; Hoatson, G. L.; Vold, R. L. J. Chem. Phys. 1999, 110, 4533−4553. (22) Kitchin, S. J.; Ahn, S.; Harris, K. D. M. J. Phys. Chem. A 2002, 106, 7228−7234. (23) O’Dell, L. A.; Schurko, R. W. J. Am. Chem. Soc. 2009, 131, 6658−6659. (24) Tycko, R.; Opella, S. J. J. Chem. Phys. 1987, 86, 1761−1774. (25) O’Dell, L. A. Prog. Nucl. Magn. Reson. Spec. 2011, 59, 295−318. (26) Giavani, T.; Johannsen, K.; Jacobsen, C. J. H.; Blom, N.; Bildsøe, H.; Skibsted, J.; Jakobsen, H. J. Solid State Nucl. Magn. Reson. 2003, 24, 218−235. (27) Gervais, C.; Dupree, R.; Pike, K. J.; Bonhomme, C.; Profeta, M.; Pickard, C. J.; Mauri, F. J. Phys. Chem. A 2005, 109, 6960−6969. (28) O’Dell, L. A.; Schurko, R. W. Phys. Chem. Chem. Phys. 2009, 11, 7069−7077. (29) Wu, G. Prog. Nucl. Magn. Reson. Spec. 2008, 52, 118−169. (30) Wu, G.; Dong, S. J. Am. Chem. Soc. 2001, 123, 9119−9125. (31) Dong, S.; Yamada, K.; Wu, G. Z. Naturforsch. 2000, 55a, 21−28. (32) Lundgren, J.-O.; Williams, J. M. J. Chem. Phys. 1973, 58, 788− 796.

1014

dx.doi.org/10.1021/jp210844t | J. Phys. Chem. A 2012, 116, 1008−1014