Refining Crystal Structures with Quadrupolar NMR ... - ACS Publications

Dec 20, 2017 - Department of Chemistry and Biochemistry, University of Windsor, Windsor, ON, Canada N9B 3P4. •S Supporting Information. ABSTRACT: Nu...
0 downloads 4 Views 1MB Size
Subscriber access provided by University of Florida | Smathers Libraries

Article

Refining Crystal Structures with Quadrupolar NMR and Dispersion-Corrected Density Functional Theory Sean T. Holmes, and Robert W. Schurko J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b12314 • Publication Date (Web): 20 Dec 2017 Downloaded from http://pubs.acs.org on January 3, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Refining Crystal Structures with Quadrupolar NMR and Dispersion-Corrected Density Functional Theory

Sean T. Holmes and Robert W. Schurko* Department of Chemistry and Biochemistry, University of Windsor, Windsor, ON, Canada N9B 3P4 *Author to whom correspondence should be addressed. E-mail: [email protected] Tel: (519) 253-3000 x3548 - Fax: (519) 973-7098

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Nuclear electric field gradient (EFG) tensors obtained from solid-state NMR spectroscopy are highly responsive to variations in structural features. The orientations and principal components of EFG tensors show great variation between different molecular structures; hence, extraction of EFG tensor parameters, either experimentally or computationally, provides a powerful means for structure determination and refinement. Here, dispersioncorrected plane-wave density functional theory (DFT) is used to refine atomic coordinates in organic crystals determined initially through single-crystal X-ray diffraction (XRD) or neutron diffraction methods. To accomplish this, an empirical parameterization of a two-body dispersion force field is illustrated, in which comparisons of experimental and calculated 14N, 17O, and 35Cl EFG tensor parameters are used to assess the quality of energy-minimized structures. The parameterization is based on a training set of 17 organic solids. The analysis is applied subsequently to the structural refinements of structural models from over 60 different materials. For the prediction of 35Cl EFG tensor parameters in particular, the optimization protocols described herein lead to a substantial improvement in agreement with experiment relative to structures obtained by XRD methods or by refinement with plane-wave DFT without the inclusion of the force field. The results further demonstrate that crystal structures with atomic coordinates refined with the present methods are able to pinpoint the positions of hydrogen atoms participating in H···Cl- hydrogen bonding with a higher degree of precision than is possible through neutron diffraction. This methodology, which is facile to implement within most DFT software packages, should prove to be very useful for future structural refinements using NMR crystallographic methods.

2 ACS Paragon Plus Environment

Page 2 of 39

Page 3 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

INTRODUCTION Solid-state nuclear magnetic resonance (SSNMR) spectroscopy is an extremely powerful technique for studying molecular-level structure and dynamics.1-3 The determination of threedimensional molecular structure is a central concern for crystallographers, solid-state chemists, and crystal engineers due to connections between molecular-level properties and bulk physicochemical properties.4 Relationships between molecular structures and NMR observables can be elucidated through first-principles calculations, typically performed in the framework of density functional theory (DFT).5-10 The interdisciplinary field of NMR crystallography combines SSNMR spectroscopy, X-ray diffraction (XRD) methods, and computational approaches, to provide unrivaled insight into molecular-level structures; these insights extend to the enhancement of structures solved by refinement of XRD powder patterns, or even to the solution of crystal structures independent of XRD data. Successes in this field include the ability to measure distances between spins,11-13 characterize disorder in solids,14-18 detect weak intermolecular interactions,19 study biologically-relevant materials that are often unsuited for other types of analyses,20-24 and reveal distortions in molecular symmetry that are not observed by diffraction methods alone.25 Several reports have illustrated that the precision of atomic positions obtained through SSNMR can rival that obtained through neutron diffraction.26-28 The majority of NMR crystallographic refinements have utilized measurements of isotropic chemical shifts or chemical shift tensors of spin-½ nuclides.29-34 There are numerous impediments associated with refining crystal structures using chemical shifts, including difficulties in measuring 1H chemical shifts (due to strong homonuclear dipolar coupling), complications in 13C or 15N SSNMR spectra due to overlapping peaks (especially in larger molecules), and the fact that much of the three-dimensional structural information revealed by

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 39

the chemical shift anisotropy is only obtainable from measurements of the principal elements of the chemical shift tensors (often, this requires lengthy two-dimensional NMR experiments). SSNMR of quadrupolar nuclides (I > ½) is particularly useful for probing molecular-level structure, since nuclear electric field gradient (EFG) tensors are responsive to even the smallest variation in structural features.35-40 For the study of organic solids, two of the most ubiquitous quadrupolar NMR nuclides are 14N41-50 and 17O,47, 51-60 which both have been used to study molecular-level structure and dynamics in solids.

14

N (I = 1) is a quadrupolar nucleus with a

high natural abundance (99.63%) and a low magnetogyric ratio [γ(14N) = 1.934 × 107 rad T-1 s-1]. Ultra-wideline 14N SSNMR patterns are inhomogeneously broadened by the first-order quadrupolar interaction, and can be up to several MHz in breadth, making their acquisition challenging; however, signal enhancing pulse sequences and the increasing availability of highfield NMR spectrometers have made such measurements tenable.61, 62 5

17

O SSNMR spectra (I =

/2) are also very challenging to obtain due to the low natural abundance (0.037%) and the low

magnetogyric ratio [γ(17O) = -3.628 × 107 rad T-1 s-1] of 17O the nucleus, and the possibility of substantial inhomogeneous line broadening that is complicated by the presence of both the second-order quadrupolar and chemical shift anisotropy interactions.63 Analysis of 17O sites in organic solids through SSNMR spectroscopy is becoming more tenable because of advances in sensitivity enhancing techniques (including dynamic nuclear polarization)13 and new synthetic pathways for isotopic labeling.64 There also exists a substantial body of literature dealing with nuclear quadrupole resonance (NQR) spectroscopy studies of 17O sites in organic solids.65-67 When solids contain cations or anions, which are frequently quadrupolar nuclides, analysis of the NMR quadrupolar powder patterns can reveal significant structural information because of the key role that these sites play in crystal packing.36, 37 Chloride anions in

4 ACS Paragon Plus Environment

Page 5 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

hydrochloride (HCl) and metal chloride salts act as hydrogen bond acceptors, forming both charge-assisted hydrogen bonds and weaker hydrogen bonds with neutral hydrogen-bond donating moieties.68 The variability in local chloride anion environments makes 35Cl (I = 3/2) SSNMR spectroscopy a powerful fingerprinting tool that is highly sensitive to solvation and polymorphism, as has been illustrated in 35Cl SSNMR analyses of chloride-anion-containing materials.69-81 These observations are particularly relevant for the characterization of active pharmaceutical ingredients (APIs), of which, greater than 50% of solid-state formulations are prepared as HCl salts.82 Prior work has demonstrated difficulties associated with first-principles predictions of 35

Cl EFG tensors of chloride anions,69-71, 77-80, 83 limiting the utility of this sensitive NMR probe

for resolving structural features of solids. It is not uncommon for calculated values of the individual principal elements of 35Cl EFG tensors to overestimate experimental values by several MHz when calculations are performed on structures obtained from diffraction methods, or DFT refinements of XRD-derived structures. Ab initio molecular dynamics study has illustrated that fast molecular dynamics (fs – ps time scales) do not have a large enough impact on calculated 35

Cl EFG tensor parameters to be the leading cause of this error.83, 84 Calculated 35Cl EFG

tensors differ substantially between structures refined with different model chemistries,85 suggesting that the source of the discrepancy between experiment and theory results from inadequate structures in the computational models. The positions of hydrogen atoms, which are essential components of the structures of HCl salts due to H···Cl- hydrogen bonding, are poorly resolved by single-crystal (sc)XRD analyses. It is therefore likely that the proper positioning of hydrogen atoms is the most significant factor for calculating 35Cl EFG tensors.

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

In ab initio calculations on models representing solids, dispersion, electrostatic, and exchange-correlation (xc) interactions control crystal packing.86 Calculations of physical quantities that are highly dependent upon structural parameters are affected strongly by the manner in which these interactions are treated.87-89 The dispersion interaction is often introduced into calculations as a parameterized force field that depends explicitly upon the set of atomic coordinates defining the positions of the atoms in the solid. Various quantities, including cohesive energies, lattice constants, bulk moduli, lattice vibrations, etc., have been employed to parameterize, or to benchmark, the performance of dispersion force fields to aid in their prediction in organic crystals.87-89 Calculations of nuclear EFG tensors are also suitable for parameterizing the dispersion force field, because these quantities are strongly responsive to structural variation, including both nuclear positions and electronic structure.85 One advantage to this approach is that anisotropic NMR interactions, such as nuclear EFG tensors, require proper modeling of three-dimensional structure to obtain agreement between experiment and theory. An empirically parameterized force field model based on quadrupolar NMR parameters could correct errors in both the treatments of dispersion and xc interactions, with wide-ranging applicability to the field of NMR crystallography. This study assesses the quality of crystal structure refinements through calculations of 14

N, 17O, and 35Cl EFG tensor parameters in neutral organic molecules and salts in the solid state.

Calculations are performed within the framework of plane-wave DFT, which inherently accounts for the long-range periodic order of crystal structures. Consistent with previous work on calculations of 35Cl EFG tensors in organic solids,69 we illustrate that conventional plane-wave DFT approaches alone are insufficient to refine atomic coordinates in crystals such that the resulting structures are able to accurately predict nuclear EFG tensor parameters.85 We

6 ACS Paragon Plus Environment

Page 6 of 39

Page 7 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

demonstrate that much of the discrepancy between experimental and calculated values arises from uncertainties or errors in the positions of atoms (especially hydrogen atoms) in crystal structures obtained by diffraction methods. Crystal structures are refined by combing DFT with a force field that is parameterized such that the refined structures yield accurate and reliable predictions of EFG tensors in organic solids. This analysis is applied subsequently to a test set consisting of 69 organic solids and 160 NMR-active nuclear sites. The resulting NMR-driven structures are able to model experimental nuclear EFG tensor parameters with a higher degree of accuracy than is obtained by previous calculations on XRD-based structures or structures refined by conventional plane-wave DFT. Finally, calculated 35Cl EFG tensor parameters obtained from calculation on the NMR-driven structures are demonstrated to be superior to those structures obtained by neutron diffraction methods.

COMPUTATIONAL DETAILS Overview of Computational Protocols. All calculations were performed on models of structures determined from previous scXRD or neutron diffraction studies, as indicated. Crystal structures were obtained from the Cambridge Structural Database (CSD) using CCDC ConQuest (version 1.19). For the XRD-derived structures of L-phenylalanine HCl, dibucaine HCl, and alprenolol HCl, several hydrogen atoms were absent from the experimentally-derived crystal structures, and were placed manually to satisfy the valences of the directly bonded atoms. In addition, a disordered oxygen atom in the crystal structure of 4-nitrobenzene was removed. All final sets of xyz atomic coordinates are provided in the Supporting Information (SI). Several strategies for refining the atomic coordinates in the crystal structures were explored, based on plane-wave DFT calculations,5 as implemented in the CASTEP module of

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 39

BIOVIA Materials Studio 2017 R2.90-92 All CASTEP calculations employed the generalized gradient approximation (GGA) revised Perdew-Burke-Ernzerhof (rPBE) xc functional,93 which performs well for the calculation of NMR parameters for 3rd- and 4th-period nuclides in solids.94, 95

Core-valence interactions were modeled with ultrasoft pseudopotentials generated on the

fly.96 Integrals over the Brillouin zone were evaluated using a Monkhorst-Pack grid with a kpoint spacing of 0.07 Å-1, and a plane-wave cutoff energy of 700 eV was used in all calculations. These parameters were chosen on the basis of a series of benchmark calculations on HCl salts (see Additional Methodological Information in the SI). Structural refinements employed the Broyden-Fletcher-Goldfarb-Shanno energy-minimizing scheme.97 Lattice constants remained fixed during the geometry optimizations; however, the computational methodologies outlined in this work could be applied to future analyses in which unit cell volumes vary during geometry optimizations. The thresholds for structural convergence include a maximum change in energy of 5 × 10-6 eV atom-1, a maximum displacement of 5 × 10-4 Å atom-1, and a maximum Cartesian force of 0.01 eV Å-1. Where indicated, dispersion was included in the calculations through the semi-empirical two-body force field model of Grimme,86, 98 or a modification thereof (vide infra). Practical implementation of the force field model in CASTEP is discussed in the SI. Statistical Analysis. The nuclear EFG is characterized by a symmetric, traceless, second-rank tensor with principal elements defined such that |V33| ≥ |V22| ≥ |V11|. The interaction between a nuclear electric quadrupole moment and the EFG with its origin at the nucleus is typically characterized by the quadrupolar coupling constant (CQ) and the asymmetry parameter (ηQ):

 =  /ℎ ,

η = ( −  )/ . 8 ACS Paragon Plus Environment

(1) (2)

Page 9 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In the above expressions, e is the elementary charge, h is Planck’s constant, and Q is the nuclear quadrupole moment [Q(14N) = 2.044 fm2, Q(17O) = -2.558 fm2, Q(35Cl) = -8.165 fm2].99 Because the sign of CQ cannot be determined from the NMR experiments in most cases, our analysis only compares the absolute values of the experimentally and theoretically derived values. We utilize several metrics to evaluate the agreement between calculated and experimental principal elements of EFG tensors. The first of these metrics is the root mean square relative error (rmsre), which is given by the following expression:   , − , 1  =   $ . !"# 3 , %





!"#





(3)

  Here, , and , are the calculated and experimental principal elements, respectively, of the !"#

EFG tensor of nuclear site m. The summations are over the M nuclear sites (m = 1, 2, . . . M) and the k principal elements of each EFG tensor (k = 1, 2, 3). Errors are reported relative to V33 (rather than Vkk) to avoid overweighting values of V11 near zero. The root mean square error (rmse) is given by the following expression:

 1 !"#      =  &, − , ' . ℎ 3 %







(4)

In addition, we use as a figure of merit the maximum positive deviation [max(+)] and maximum negative deviation [max(-)] between calculation and experiment in the magnitude of the principal elements, multiplied by eQ/h. Because the EFG tensor is traceless ( +  +  = 0), the

mean error for each dataset is essentially zero.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

RESULTS AND DISCUSSION Overview. Herein, we present a DFT-based computational protocol for refining the coordinates of atoms in crystal structures derived from diffraction methods, such that calculations on the resulting structures reproduce experimental nuclear EFG tensor parameters derived from NMR or NQR spectroscopy. Calculations employ a force field that is parameterized empirically to aid in the refinement of atomic coordinates in organic solids: specifically, the parameterization utilizes DFT calculations of 14N, 17O, and 35Cl EFG tensor parameters. The parameterization of the force field employs a training set of organic solids, including urea, formic acid, acetic acid, oxalic acid (two polymorphs), glycine (three polymorphs), glycine HCl, L-alanine HCl, L-valine HCl, L-cysteine HCl·H2O, L-threonine HCl, L-aspartic acid HCl, Lhistidine HCl·H2O, aminoguanidine HCl, and dopamine HCl. The training set contains 34 unique nuclear sites, with 8, 17, and 9 sites for 14N, 17O, and 35Cl nuclides, respectively. The crystal structures and quadrupolar NMR parameters of these solids are summarized in Table S5. Results are presented for structures in which the atomic coordinates determined from XRD are not refined, and for structures refined (i) with the reparameterized version of the force field, (ii) with a previous parameterization of the dispersion force field,86 and (iii) without any force field (i.e., structures are refined using DFT alone). Differences in structural parameters among the energy-minimized structures and XRD-derived structures are discussed. Calculations using the parameterized dispersion-corrected DFT structural refinement protocol are applied to the prediction of 14N, 17O, and 35Cl EFG tensors in a large test set of organic solids (Tables S7 – S10). Together with the materials included in the training set, the complete analysis consists of 86 solid materials and 160 nuclear sites, including, 58, 56, and 46

10 ACS Paragon Plus Environment

Page 10 of 39

Page 11 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

sites for 14N, 17O, and 35Cl nuclides, respectively. Finally, calculations of 35Cl EFG tensor parameters in a subset of materials that have also been characterized through neutron diffraction are examined. Parameterization of the Dispersion Force Field using the Training Set. In Grimme’s two-body semi-empirical dispersion model (DFT-D2),86 the isotropic two-body dispersion energy (*+,-# ) is treated as an atomic-pairwise correction to the Kohn-Sham energy (*./ ). The total energy in the DFT-D2 scheme is given by:

*01230 = *./ + *+,-# .

(5)

The dispersion energy is defined by the following unscaled expression: ;

*+,-# = − 

6,7 4.0. Figure 1 illustrates the rmsre between experimental and calculated principal elements of EFG tensors, obtained from calculations on structures with atomic coordinates refined for values of d between 2.5 and 10.0. The rmsre for the three sets of nuclides vary significantly over the range of values of 2.5 ≤ d ≤ 10, but vary little for values of d greater than 10. rmsre minima are observed at d = 4.00, 3.50, and 3.25 for 14N, 17O, and 35Cl sites, respectively. Considering all 34 sites together, regardless of the type of nucleus, a global minimum is observed at d = 3.25. Structures resulting from refinements employing the rPBE-D2 model chemistry, with the damping parameter set to a value of 3.25, are denoted rPBE-D2* structures. Structures obtained from refinements using the more common form of the dispersion function, where the damping parameter is set to a value of 20.0, are denoted rPBE-D2 structures. Table 1 summarizes statistical data associated with the prediction of EFG tensor parameters on training set using model structures obtained by XRD, and on structures with atomic coordinates refined at the rPBE, rPBE-D2, and rPBE-D2* levels. Calculations on structures refined at the rPBE-D2* level result in lower rmsre and rmse values than those obtained from the XRD structures and structures refined at the rPBE and rPBE-D2 levels. Calculations on rPBE and rPBE-D2 refined structures result in EFG tensor principal elements

12 ACS Paragon Plus Environment

Page 12 of 39

Page 13 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

that are similar in a statistical sense and frequently on a case-by-case basis: calculated values of the principal elements of 35Cl EFG tensors in the training set differ by only 160 kHz between the two classes of structures. In comparison, differences in the same values between rPBE-D2 and rPBE-D2* structures are an order of magnitude larger (1.25 MHz). Because nuclear EFG tensor parameters obtained for structures optimized at the rPBE and rPBE-D2 levels are similar, calculations at the rPBE-D2 level are not reported subsequently. The empirical structural corrections are based on experimental NMR parameters collected typically at room temperature, whereas the single crystal XRD structures were obtained between 77 K and room temperature. The effect of small lattice expansions/contractions associated with changes in temperature appear to have only small effects on calculated EFG tensor parameters. For instance, the calculated values of CQ(35Cl) for the chloride anion in Lhistidine HCl·H2O differ by only 120 kHz for structures obtained at room temperature and 100 K, following optimizations. The effects of fast molecular dynamics on calculated EFG tensor parameters are beyond the scope of the present study. As implemented in this study, the dispersion correction plays only an indirect role in calculations of EFG tensors: specifically, inclusion of dispersion does not affect the algorithm that provides the EFG tensor; rather, it affects the structural refinements, which in turn yield atomic coordinates and ground state electron densities that determine the nature of the EFG tensor. Calculations of the principal elements of the 35Cl EFG tensor on a single structural model of L-histidine HCl·H2O differ by only a few hundred hertz, whether or not the dispersion force field is included during the calculation of the EFG tensors (Table S4). This difference reflects the numerical convergence of the calculations.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analysis of Training Set Structures. Geometry optimization of the training set organic solids at any level leads to a reduction of the mean forces on the atoms by three to four orders of magnitude (Table S5). Mean atomic forces for the training set solids never differ by more than 0.006 eV Å-1 between structures obtained using different structural refinement protocols. However, variation in calculated EFG tensor parameters between structures refined at the rPBE and rPBE-D2* levels suggests that the energy-minimized structures differ substantially. Limiting the discussion of structural parameters to the nine hydrochloride salts in the training set, refinement of atomic coordinates at the rPBE-D2* level results in shorter covalent bond lengths than predicted at the rPBE level (Figure S1). Relative to bond lengths observed by XRD, covalent bonds between non-hydrogen atoms are predicted to be longer, following refinement at the rPBE level (by an average of 0.014 Å), and shorter, following refinement at the rPBE-D2* level (by an average of 0.012 Å). Covalent bonds involving hydrogen atoms are lengthened relative to XRD-derived structures by both structural refinement protocols (by 0.117 Å and 0.048 Å for rPBE and rPBE-D2* refinements, respectively). Root-mean-square deviations (RMSDs) in atomic positions highlight differences in structures prior to and following geometry optimization (Figure S2).26-28 Considering the positions of non-hydrogen atoms, the RMSD between XRD and rPBE structures is 0.057 Å on average, with a maximum deviation of 0.16 Å for any material. The equivalent deviations between XRD and rPBE-D2* structures are 0.044 Å and 0.11 Å. These data indicate that refinements at the rPBE-D2* level result in smaller deviations from the XRD structures than refinements at the rPBE level, on average. For hydrogen atoms, the RMSDs between XRDderived structures and refined structures are much larger (0.178 Å and 0.126 Å for the rPBE and rPBE-D2* refinements, respectively) than observed for the heavier atoms (Figure S3). The

14 ACS Paragon Plus Environment

Page 14 of 39

Page 15 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

hydrogen atom RMSDs are generally smaller between the two types of optimized structures (0.099 Å) than between the refined structures and the XRD-derived structures. RMSDs between scXRD-derived structures and dispersion-corrected DFT structures reported herein are well within the expected range reported by van de Streek.103 The large variation in calculated 35Cl EFG tensor parameters between the XRD, rPBE, and rPBE-D2* structures suggests that the chloride anion environments in the vicinity of these sites differs significantly between the various types of structures. Table 2 lists the H···Cl- contact types and contact distances for a subset of the training set organic solids, as obtained by XRD and refinement by the two model chemistries. Included are all H···Cl- contacts under 2.64 Å, as this is the sum of the van der Waals’ radii of the hydrogen and chlorine atoms used in the dispersion force field (Eq. 7).86 In nearly every case, rPBE-D2* H···Cl- contacts are longer than those obtained at the rPBE level. Differences in H···Cl- contacts are typically largest for the shortest contact, with variations over 0.1 Å, in some cases. This observation does not hold for Lhistidine HCl·H2O, where the first two contacts are of comparable length. Compared to H···Clcontacts determined by XRD, rPBE-D2* contacts differ by as much as 0.2 Å. Calculations of EFG tensor parameters for the test set. The performance of rPBED2* for materials outside of the initial training set is tested through calculations of 35Cl, 17O, and 14

N, EFG tensors from a wide array of organic solids. The results of the calculations are

illustrated in Figure 2, and statistical analysis is provided in Tables 3 and 4. The largest variation in calculated EFG tensor parameters between different structure types is observed for 35Cl sites. Considering all of the principal elements of the 35Cl EFG tensors together, rmse values of 1.51 MHz and 1.30 MHz were observed for the XRD-derived structures and rPBE refined structures, respectively. The errors for the rPBE structures are highly

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

systematic; the relationship between calculated and experimental values is modeled with a linear function with a slope of 1.32. This systematic overestimation of the magnitude of EFG tensor principal elements has been observed in previous work.83 We believe that the source of this error is related to the underestimation of H···Cl- contacts in the crystal structures when a proper treatment of dispersion is not included in the geometry optimization (vide supra). In extreme cases, calculated values of CQ for rPBE structures differ from experiment by over 3.5 MHz. In contrast, structural refinement at the rPBE-D2* level leads to a rmse of only 0.40 MHz, which is ca. 3 times smaller than the value obtained for rPBE-refined structures. The variation in calculated 35Cl EFG tensors is also observed in each of the three principal elements individually, with the best agreement with experiment obtained by calculations on rPBE-D2* structures (Table 4). Calculations of 17O EFG tensors are strongly influenced by the type of structure used in the calculation. Because all of the 17O SSNMR studies cited here involve oxygen atoms involved in covalent bonding, the variability in computed EFG tensor parameters is significantly smaller than observed for 35Cl sites. For XRD-derived structures, rPBE refined structures, and rPBE-D2* refined structures, the rmse values are 0.58 MHz, 0.49 MHz, and 0.36 MHz, respectively. Refinement at the rPBE level leads to notably worse predictions of eQV22/h than is predicted by either XRD structures or structures refined at the rPBE-D2* level. This difference is present throughout the majority of systems studies here, independent of the type of 17Ocontaining moiety, including oxygen sites that do not participate in hydrogen bonding. It is worth comparing the 17O EFG tensors reported herein to those from an earlier study.85 Previously, an rmse of 0.28 MHz was reported for 17O sites based on calculations using dispersion-corrected DFT, which is slightly lower than the rmse of 0.36 MHz reported herein.

16 ACS Paragon Plus Environment

Page 16 of 39

Page 17 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Three facts about the previous study should be noted: (i) the xc functional PW91 was used rather than rPBE; (ii) it employed a version of the DFT-D2 model parameterized using only 17O EFG tensor parameters, rather than the heteronuclear parameterization used here; (iii) it was based on a smaller number of materials, exhibiting a limited diversity of moieties. Based on these observations, it is possible that the introduction of a set of distinct values for the damping parameter, for all possible types of atomic pairs, will further improve upon the results obtained using a global damping parameter (i.e., a single value of d for all atomic-pairwise corrections). This is beyond the scope of the current study. Calculated 14N EFG tensor parameters appear to be the least sensitive to structural refinement (in terms of rmse) of the three nuclides discussed herein. The rmse for XRD-derived structures is 0.33 MHz. This quantity is reduced to 0.17 MHz and 0.15 MHz following structural refinement at the rPBE and rPBE-D2* levels, respectively. Differences in calculated 14N EFG tensor parameters between rPBE and rPBE-D2* refined structures are largest for the principal element eQV22/h. rPBE errors in 14N EFGs are most strongly pronounced for V22 of pseudotetrahedral nitrogen sites with 0.31 ≤ ηQ ≤ 0.70. There are nine 14N-containing materials with an experimentally-determined asymmetry parameter in this range (Table 5). In contrast, refinement at the rPBE-D2* level results in superior agreement with experimental values (Figure 3). Pseudotetrahedral 14N sites with ηQ in this range are indicative of species with NH3+ groups that act as hydrogen bond donors to nearby oxygen sites. Analysis of the rPBE-D2* refined structures reveals greater variation in N-H bond lengths in NH3+ groups than species with low values of ηQ (0.00 ≤ ηQ ≤ 0.30), in agreement with previous observations.43 The greater variation in N-H bond lengths results from the greater variation in O···H hydrogen bonding distances.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

When the N-H and O···H bond lengths are similar, the 14N site exhibits higher axial symmetry (lower ηQ). Comparison with Neutron Diffraction Structures. Structures derived from neutron diffraction are often considered superior to those derived from XRD due to the ability of the former to determine the positions of hydrogen atoms with a greater degree of precision.26-29 Proper positioning of hydrogen atoms is critical for accurate prediction of EFG tensor parameters, particularly those of 35Cl EFG tensors of chloride ion sites involved in hydrogen bonds. Comparison of calculation of 35Cl EFG tensor parameters between neutron-diffraction derived structures and rPBE or rPBE-D2* refinements of XRD-derived structures provides a robust method of assessing structural refinements. The crystal structures of glycine HCl, L-valine HCl, L-glutamic acid HCl, L-histidine HCl·H2O, L-lysine HCl·2H2O, and L-phenylalanine HCl have been determined from neutron diffraction,104-109 and are used here as a final figure of merit for assessing structural refinement by the rPBE-D2* model (Table 6). For each material, calculations of 35Cl EFG tensor parameters for rPBE-D2* refined structures results in superior agreement with experiment than calculations on models based on structures derived from neutron diffraction (Figure 4). The poorest agreement with experimental 35Cl EFG tensor parameters is obtained from rPBE-refined structures in each case. H···Cl- contacts determined by rPBE-D2* refinement are longer than those determined by neutron diffraction or rPBE refinement (Table 7), by an average of 0.046 Å and 0.078 Å, respectively. Neutron-diffraction-derived H···Cl- contacts are longer than rPBE contacts by an average of 0.032 Å. These observations points to the longer H···Cl- contacts determined by rPBE-D2* refinements being better models to the true intermolecular structure of the

18 ACS Paragon Plus Environment

Page 18 of 39

Page 19 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

hydrochloride salts since the best agreement with experimental 35Cl EFG tensor parameters is achieved for structures refined with this model chemistry.

CONCLUSIONS Calculated nuclear EFG tensors are sensitive probes of molecular-level structure. Not only are these tensor parameters useful in assessing the quality of models representing the structures of solids, but experimental determinations of nuclear EFG tensor parameters can be used to inform ab initio structural refinements. Plane-wave DFT calculations of nuclear EFG tensor parameters were used here as a basis for refining the positions of atoms in crystal structures determined initially through scXRD. The analysis was based on predictions of 14N, 17

O, and 35Cl EFG tensor parameters in 86 organic solids (17 and 69 molecules in the training set

and test set, respectively). Calculations on models of solid materials were carried out using plane-wave DFT with the inclusion of empirical force fields. The success of these calculations was predicated upon parameterization of the damping function in Grimme’s two-body dispersion force field86 to aid in structural refinements, by assessing the quality of the resulting structures through calculations of nuclear EFG tensor parameters. By adjusting the value of a single term in the force field (specifically, the damping parameter), one is able to achieve significantly improved agreement with experimental values for a large variety of organic solids. Calculations of nuclear EFG tensors on crystal structures optimized through this protocol are able to predict experimental values with a higher accuracy than those obtained through calculations on models representing XRD-derived structures, DFT refinements of XRD-derived structures, and even neutron diffraction structures (when such data are available for comparison).

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The structural refinement protocol described here produces structures that differ from XRD-derived structures, structures obtained with conventional DFT refinements of diffraction structures, and even neutron diffraction structures. The most significant structural differences are observed for H···Cl- hydrogen bonds, which can vary by more than 0.2 Å between different types of structural data/refinements. In contrast, bond lengths between non-hydrogen atoms are generally shortened by 0.012 Å, relative to XRD-derived structures. 35

Cl SSNMR has proven to be one of the most sensitive structural probes of chloride

anions in organic HCl or metal chloride salts. The variation in calculated quadrupolar coupling constants among different types of structural data can be a large as 4 MHz, which is ca. 50% as large as the largest values typically observed for these sites. The extraordinary sensitivity of 35Cl EFG tensors to local structure should prove to be a powerful tool in future NMR crystallography investigations, as this technique may be the only method able to determine the geometry H···Clbonding environments with a level of precision rivaling neutron diffraction. We note that the methodology developed herein has focused exclusively on organic solids; however, the introduction of an appropriate training set could allow for the extension of this work to include extended network solids. Additionally, these protocols may be expanded to include types of atoms or bond (e.g., covalently bound chlorine atoms) that were not considered in the present study. Finally, there are advanced computational methods that may be applied, such as the D3 correction and vdW-DFT;110 however, the former has not yet been implemented in CASTEP, and the latter requires the adjustment of a larger parameter set. Nonetheless, these may be good candidates for future explorations in crystal structure refinement utilizing EFG tensors. There are exciting prospects for this generalized method to be used for the refinement and prediction of new crystal structures utilizing EFG tensor parameters derived from this method, either

20 ACS Paragon Plus Environment

Page 20 of 39

Page 21 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

independently or alongside of additional computations of chemical shift tensors and/or experimentally derived X-ray and neutron diffraction data.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ASSOCIATED CONTENT Supporting Information Practical implementation of the force field model in Materials Studio; list of CCSD codes for relevant crystal structures; plane-wave DFT benchmark calculations; additional structural analysis of energy-minimized structures; summary of calculated nuclear EFG tensor parameters. Crystallographic Information Crystallographic information files for each of the refined structures discussed in this article.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]; Tel: (519) 253-3000 x3548; Fax: (519) 973-7098 Notes The authors declare no competing financial interests.

ACKNOWLEDGEMENTS We thank both Genentech and the Natural Sciences and Engineering Research Council of Canada (NSERC) for funding this research (NSERC in the form of a Discovery Grant). R.W.S. is also grateful for a 50th Anniversary Golden Jubilee Chair from the University of Windsor. This work was made possible, in part, by the facilities of the Shared Hierarchical Academic Research Computing Network (SHARCNET: www.sharcnet.ca). We are grateful to Prof. Maria Baias (New York University Abu Dhabi) and Prof. Oded Hod (Tel Aviv University) for very helpful discussions.

22 ACS Paragon Plus Environment

Page 22 of 39

Page 23 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

REFERENCES 1. Harris, R. K.; Wasylishen, R. E.; Duer, M. J., NMR Crystallography. John Wiley & Sons Ltd.: Singapore, 2009. 2. Ashbrook, S. E.; McKay, D., Combining solid-state NMR spectroscopy with first-principles calculations - a guide to NMR crystallography. J. Chem. Soc., Chem. Commun. 2016, 52, 7186-204. 3. Bryce, D. L.; Taulelle, F., NMR Crystallography. Acta Crystallogr., Sect. C 2017, 73, 126-127. 4. Desiraju, G. R.; Vittal, J. J.; Ramanan, A., Crystal Engineering: A Textbook. World Scientific Publishing Co. Pte. Ltd.: Singapore, 2011. 5. Bonhomme, C.; Gervais, C.; Babonneau, F.; Coelho, C.; Pourpoint, F.; Azais, T.; Ashbrook, S. E.; Griffin, J. M.; Yates, J. R.; Mauri, F.; Pickard, C. J., First-principles calculation of NMR parameters using the gauge including projector augmented wave method: a chemist's point of view. Chem. Rev. 2012, 112, 5733-5779. 6. Hofstetter, A.; Emsley, L., Positional variance in NMR crystallography. J. Am. Chem. Soc. 2017, 139, 2573-2576. 7. Dudenko, D. V.; Yates, J. R.; Harris, K. D. M.; Brown, S. P., An NMR crystallography DFT-D approach to analyse the role of intermolecular hydrogen bonding and p-p interactions in driving cocrystallisation of indomethacin and nicotinamide. CrystEngComm 2013, 15, 8797-8807. 8. Beran, G. J.; Hartman, J. D.; Heit, Y. N., Predicting molecular crystal properties from first principles: finite-temperature thermochemistry to NMR crystallography. Acc. Chem. Res. 2016, 49, 2501-2508. 9. Pawlak, T.; Jaworska, M.; Potrzebowski, M. J., NMR crystallography of a-poly(L-lactide). Phys. Chem. Chem. Phys. 2013, 15, 3137-3145. 10. Pawlak, T.; Potrzebowski, M. J., Fine refinement of solid-state molecular structures of leu- and metenkephalins by NMR crystallography. J. Phys. Chem. B 2014, 118, 3298-3309. 11. Laurencin, D.; Gervais, C.; Wong, A.; Coelho, C.; Mauri, F.; Massiot, D.; Smith, M. E.; Bonhomme, C., Implementation of high resolution 43Ca solid state NMR spectroscopy: toward the elucidation of calcium sites in biological materials. J. Am. Chem. Soc. 2009, 131, 13430-13440. 12. Fyfe, C. A.; Wong-Moon, K. C.; Huang, Y.; Grondey, H., INEPT experiments in solid-state NMR. J. Am. Chem. Soc. 1995, 117, 10397-10398. 13. Perras, F. A.; Wang, Z.; Naik, P.; Slowing, I. I.; Pruski, M., Natural abundance 17O DNP NMR provides Precise O−H distances and insights into the Brønsted acidity of heterogeneous catalysts. Angew. Chem. Int. Ed. 2017, 56, 9165-9169. 14. Moran, R. F.; Dawson, D. M.; Ashbrook, S. E., Exploiting NMR spectroscopy for the study of disorder in solids. Int. Rev. Phys. Chem. 2017, 36, 39-115. 15. Kerr, H. E.; Softley, L. K.; Suresh, K.; Nangia, A.; Hodgkinson, P.; Evans, I. R., A furosemideisonicotinamide cocrystal: an investigation of properties and extensive structural disorder. CrystEngComm 2015, 17, 6707-6715. 16. Kerr, H. E.; Mason, H. E.; Sparkes, H. A.; Hodgkinson, P., Testing the limits of NMR crystallography: the case of caffeine-citric acid hydrate. CrystEngComm 2016, 18, 6700-6707. 17. Harper, J. K.; Tishler, D.; Richardson, D.; Lokvam, J.; Pendrill, R.; Widmalm, G., Solid-state NMR characterization of the molecular conformation in disordered ethyl α-l-rhamnofuranoside. J. Phys. Chem. A 2013, 117, 5534-5541. 18. Powell, J.; Kalakewich, K.; Uribe-Romo, F. J.; Harper, J. K., Solid-state NMR and DFT predictions of differences in COOH hydrogen bonding in odd and even numbered n-alkyl fatty acids. Phys. Chem. Chem. Phys. 2016, 18, 12541-12549. 19. Plevin, M. J.; Bryce, D. L.; Boisbouvier, J., Direct detection of CH/p interactions in proteins. Nat. Chem. 2010, 2, 466-471. 20. Cavalli, A.; Salvatella, X.; Dobson, C. M.; Vendruscolo, M., Protein structure determination from NMR chemical shifts. Proc. Natl. Acad. Sci. USA 2007, 104, 9615-9620. 21. Mueller, L. J.; Dunn, M. F., NMR crystallography of enzyme active sites: probing chemically detailed, three-dimensional structure in tryptophan synthase. Acc. Chem. Res. 2013, 46, 2008-2017. 23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

22. Caulkins, B. G.; Bastin, B.; Yang, C.; Neubauer, T. J.; Young, R. P.; Hilario, E.; Huang, Y.-m. M.; Chang, C.-e. A.; Fan, L.; Dunn, M. F.; Marsella, M. J.; Mueller, L. J., Protonation states of the tryptophan synthase internal aldimine active site from solid-state NMR spectroscopy: direct observation of the protonated Schiff base linkage to pyridoxal-5′-phosphate. J. Am. Chem. Soc. 2014, 136, 12824-12827. 23. Caulkins, B. G.; Young, R. P.; Kudla, R. A.; Yang, C.; Bittbauer, T. J.; Bastin, B.; Hilario, E.; Fan, L.; Marsella, M. J.; Dunn, M. F.; Mueller, L. J., NMR crystallography of a carbanionic intermediate in tryptophan synthase: chemical structure, tautomerization, and reaction specificity. J. Am. Chem. Soc. 2016, 138, 15214-15226. 24. Gupta, R.; Hou, G.; Renirie, R.; Wever, R.; Polenova, T., 51V NMR crystallography of vanadium Chloroperoxidase and its directed evolution P395D/L241V/T343A mutant: protonation environments of the active site. J. Am. Chem. Soc. 2015, 137, 5618-5628. 25. Facelli, J. C.; Grant, D. M., Determination of molecular symmetry in crystalline naphthalene using solid-state NMR. Nature 1993, 365, 325-327. 26. Harper, J. K.; Iuliucci, R.; Gruber, M.; Kalakewich, K., Refining crystal structures with experimental 13 C NMR shift tensors and lattice-including electronic structure methods. CrystEngComm 2013, 15, 86938704. 27. Kalakewich, K.; Iuliucci, R.; Harper, J. K., Establishing accurate high-resolution crystal structures in the absence of diffraction data and single crystals—an NMR approach. Cryst. Growth Des. 2013, 13, 5391-5396. 28. Kalakewich, K.; Iuliucci, R.; Mueller, K. T.; Eloranta, H.; Harper, J. K., Monitoring the refinement of crystal structures with 15N solid-state NMR shift tensor data. J. Chem. Phys. 2015, 143, 194702. 29. Harper, J. K.; Grant, D. M., Enhancing crystal-structure prediction with NMR tensor data. Cryst. Growth Des. 2006, 6, 2315-2321. 30. Salager, E.; Day, G. M.; Stein, R. S.; Pickard, C. J.; Elena, B.; Emsley, L., Powder crystallography by combined crystal structure prediction and high-resolution 1H solid-state NMR spectroscopy. J. Am. Chem. Soc. 2010, 132, 2564-2566. 31. Hartman, J. D.; Day, G. M.; Beran, G. J. O., Enhanced NMR discrimination of pharmaceutically relevant molecular crystal forms through fragment-based ab initio chemical shift predictions. Cryst. Growth Des. 2016, 16, 6479-6493. 32. Dudenko, D. V.; Williams, P. A.; Hughes, C. E.; Antzutkin, O. N.; Velaga, S. P.; Brown, S. P.; Harris, K. D. M., Exploiting the synergy of powder X-ray diffraction and solid-state NMR spectroscopy in structure determination of organic molecular solids. J. Phys. Chem. C 2013, 117, 12258-12265. 33. Webber, A. L.; Emsley, L.; Claramunt, R. M.; Brown, S. P., NMR crystallography of campho[2,3c]pyrazole (Z′ = 6): combining high-resolution 1H-13C solid-state MAS NMR spectroscopy and GIPAW chemical-shift calculations. J. Phys. Chem. A 2010, 114, 10435-10442. 34. Brouwer, D. H.; Enright, G. D., Probing local structure in zeolite frameworks: ultrahigh-field NMR measurements and accurate first-principles calculations of zeolite 29Si magnetic shielding tensors. J. Am. Chem. Soc. 2008, 130, 3095-3105. 35. Wasylishen, R. E.; Ashbrook, S. E.; Wimperis, S., NMR of Quadrupolar Nuclei in Solid Materials. John Wiley and Sons: Chichester, 2012. 36. Schurko, R. W.; Jaroszewicz, M. J., Solid-state NMR of the light main group metals. In Encyclopedia of Inorganic and Bioinorganic Chemistry, John Wiley & Sons, Ltd: 2011. 37. Szell, P. M. J.; Bryce, D. L., Recent advances in chlorine, bromine, and iodine solid-state NMR spectroscopy. Annu. Rep. NMR Spectrosc. 2015, 84, 115-162. 38. Knop, O.; Palmer, E. M.; Robinson, R. W., Arrangements of point charges having zero electric-field gradient. Acta Crystallogr., Sect. A 1975, 31, 19-31. 39. Akitt, J. W.; McDonald, W. S., Arrangements of ligands giving low electric field gradients. J. Magn. Reson. 1984, 58, 401-412. 40. Autschbach, J.; Zheng, S.; Schurko, R. W., Analysis of electric field gradient tensors at quadrupolar nuclei in common structural motifs. Concepts Magn. Reson., Part A 2010, 36, 84-126.

24 ACS Paragon Plus Environment

Page 24 of 39

Page 25 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

41. O’Dell, L. A.; Schurko, R. W.; Harris, K. J.; Autschbach, J.; Ratcliffe, C. I., Interaction tensors and local dynamics in common structural motifs of nitrogen: a solid-state 14N NMR and DFT study. J. Am. Chem. Soc. 2011, 133, 527-546. 42. Veinberg, S. L.; Friedl, Z. W.; Harris, K. J.; O'Dell, L. A.; Schurko, R. W., Ultra-wideline 14N solidstate NMR as a method for differentiating polymorphs: glycine as a case study. CrystEngComm 2015, 17, 5225-5236. 43. Veinberg, S. L.; Friedl, Z. W.; Lindquist, A. W.; Kispal, B.; Harris, K. J.; O'Dell, L. A.; Schurko, R. W., 14N solid-state NMR spectroscopy of amino acids. ChemPhysChem 2016, 17, 4011-4027. 44. Veinberg, S. L.; Johnston, K. E.; Jaroszewicz, M. J.; Kispal, B. M.; Mireault, C. R.; Kobayashi, T.; Pruski, M.; Schurko, R. W., Natural abundance 14N and 15N solid-state NMR of pharmaceuticals and their polymorphs. Phys. Chem. Chem. Phys. 2016, 18, 17713-17730. 45. Garcia, M. L. S.; Smith, J. A. S., 14N and 2H quadrupole double resonance in salts of cytosine and adenine. J. Chem. Soc., Perkin Trans. 2 1983, 1401-1408. 46. Harbison, G. S.; Kye, Y.-S.; Penner, G. H.; Grandin, M.; Monette, M., 14N quadrupolar, 14N and 15N chemical shift, and 14N−1H dipolar tensors of sulfamic acid. J. Phys. Chem. B 2002, 106, 10285-10291. 47. O’Dell, L. A.; Ratcliffe, C. I.; Kong, X.; Wu, G., Multinuclear solid-state nuclear magnetic resonance and density functional theory characterization of interaction tensors in taurine. J. Phys. Chem. A 2012, 116, 1008-1014. 48. Naito, A.; Ganapathy, S.; Raghunathan, P.; McDowell, C. A., Determination of the nitrogen-14 quadrupole coupling tensor and the carbon-13 chemical shielding tensors in a single crystal of L-serine monohydrate. J. Chem. Phys. 1983, 79, 4173-4182. 49. Naito, A.; McDowell, C. A., Determination of the nitrogen-14 quadrupole coupling tensors and the carbon-13 chemical shielding tensors in a single crystal of L-asparagine monohydrate. J. Chem. Phys. 1984, 81, 4795-4803. 50. McDowell, C. A.; Naito, A.; Sastry, D. L.; Takegoshi, K., Determination of the 14N quadrupole coupling tensors in a single crystal of L-histidine hydrochloride monohydrate by NMR spectroscopy. J. Magn. Reson. 1986, 69, 283-292. 51. Kong, X.; Shan, M.; Terskikh, V.; Hung, I.; Gan, Z.; Wu, G., Solid-State 17O NMR of pharmaceutical compounds: salicylic acid and aspirin. J. Phys. Chem. B 2013, 117, 9643-9654. 52. Kong, X.; Terskikh, V.; Toubaei, A.; Wu, G., A solid-state 17O NMR study of platinum-carboxylate complexes: carboplatin and oxaliplatin. Can. J. Chem. 2015, 93, 945-953. 53. Nour, S.; Widdifield, C. M.; Kobera, L.; Burgess, K. M. N.; Errulat, D.; Terskikh, V. V.; Bryce, D. L., Oxygen-17 NMR spectroscopy of water molecules in solid hydrates. Can. J. Chem. 2016, 94, 189-197. 54. Wu, G.; Dong, S.; Ida, R.; Reen, N., A solid-state 17O nuclear magnetic resonance study of nucleic acid bases. J. Am. Chem. Soc. 2002, 124, 1768-1777. 55. Wu, G.; Mason, P.; Mo, X.; Terskikh, V., Experimental and computational characterization of the 17O quadrupole coupling and magnetic shielding tensors for p-nitrobenzaldehyde and formaldehyde. J. Phys. Chem. A 2008, 112, 1024-1032. 56. Pike, K. J.; Lemaitre, V.; Kukol, A.; Anupõld, T.; Samoson, A.; Howes, A. P.; Watts, A.; Smith, M. E.; Dupree, R., Solid-state 17O NMR of amino acids. J. Phys. Chem. B 2004, 108, 9256-9263. 57. Dong, S.; Ida, R.; Wu, G., A combined experimental and theoretical 17O NMR study of crystalline urea:  an example of large hydrogen-bonding effects. J. Phys. Chem. A 2000, 104, 11194-11202. 58. Dong, S.; Yamada, K.; Wu, G., Oxygen-17 nuclear magnetic resonance of organic solids. Z. Naturforsch. 2000, 55, 21-28. 59. Scheubel, W.; Zimmermann, H.; Haeberlen, U., 17O quadrupole coupling and nuclear magnetic shielding tensors in benzophenone. J. Magn. Reson. 1985, 63, 544-555. 60. Lemaitre, V.; Pike, K. J.; Watts, A.; Anupold, T.; Samoson, A.; Smith, M. E.; Dupree, R., New insights into the bonding arrangements of L- and D-glutamates from solid state 17O NMR. Chem. Phys. Lett. 2003, 371, 91-97. 61. Veinberg, S. L.; Lindquist, A. W.; Jaroszewicz, M. J.; Schurko, R. W., Practical considerations for the acquisition of ultra-wideline 14N NMR spectra. Solid State Nucl. Magn. Reson. 2017, 84, 45-58. 25 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

62. O’Dell, L. A., Direct detection of nitrogen-14 in solid-state NMR spectroscopy. Prog. Nucl. Magn. Reson. Spectrosc. 2011, 59, 295-318. 63. Wu, G., Solid-State 17O NMR studies of organic and biological molecules: recent advances and future directions. Solid State Nucl. Magn. Reson. 2016, 73, 1-14. 64. Métro, T.-X.; Gervais, C.; Martinez, A.; Bonhomme, C.; Laurencin, D., Unleashing the potential of 17 O NMR spectroscopy using mechanochemistry. Angew. Chem. Int. Ed. 2017, 56, 6803-6807. 65. Brosnan, S. G. P.; Edmonds, D. T.; Poplett, I. J. F., An 17O nuclear quadrupole resonance study of some carboxylic acids. J. Magn. Reson. 1981, 45, 451-460. 66. Poplett, I. J. F.; Smith, J. A. S., 17O and 2H quadrupole double resonance in some carboxylic acid dimers. J. Chem. Soc., Faraday Trans. 1981, 77, 1473-1485. 67. Lemaitre, V.; Smith, M. E.; Watts, A., A review of oxygen-17 solid-state NMR of organic materials towards biological applications. Solid State Nucl. Magn. Reson. 2004, 26, 215-235. 68. Desiraju, G. R.; Steiner, T., The Weak Hydrogen Bond in Structural Chemistry and Biology. Oxford University Press, Inc.: New York, 1999. 69. Hamaed, H.; Pawlowski, J. M.; Cooper, B. F. T.; Fu, R.; Eichhorn, S. H.; Schurko, R. W., Application of solid-state 35Cl NMR to the structural characterization of hydrochloride pharmaceuticals and their polymorphs. J. Am. Chem. Soc. 2008, 130, 11056-11065. 70. Hildebrand, M.; Hamaed, H.; Namespetra, A. M.; Donohue, J. M.; Fu, R.; Hung, I.; Gan, Z.; Schurko, R. W., 35Cl solid-state NMR of HCl salts of active pharmaceutical ingredients: structural prediction, spectral fingerprinting and polymorph recognition. CrystEngComm 2014, 16, 7334-7356. 71. Namespetra, A. M.; Hirsh, D. A.; Hildebrand, M. P.; Sandre, A. R.; Hamaed, H.; Rawson, J. M.; Schurko, R. W., 35Cl solid-state NMR spectroscopy of HCl pharmaceuticals and their polymorphs in bulk and dosage forms. CrystEngComm 2016, 18, 6213-6232. 72. Hirsh, D. A.; Rossini, A. J.; Emsley, L.; Schurko, R. W., 35Cl dynamic nuclear polarization solid-state NMR of active pharmaceutical ingredients. Phys. Chem. Chem. Phys. 2016, 18, 25893-25904. 73. Vogt, F. G.; Williams, G. R.; Copley, R. C. B., Solid-state NMR analysis of a boron-containing pharmaceutical hydrochloride salt. J. Pharm. Sci. 2013, 102, 3705-3716. 74. Vogt, F. G.; Williams, G. R.; Strohmeier, M.; Johnson, M. N.; Copley, R. C. B., Solid-state NMR analysis of a complex crystalline phase of ronacaleret hydrochloride. J. Phys. Chem. B 2014, 118, 1026610284. 75. Vogt, F. G.; Williams, G. R.; Johnson, M. N.; Copley, R. C. B., A spectroscopic and diffractometric study of polymorphism in ethyl 3-{3-[((2R)-3-{[2-(2,3-dihydro-1H-inden-2-yl)-1,1dimethylethyl]amino}-2-hydroxypropyl)oxy]-4,5-difluorophenyl}propanoate hydrochloride. Cryst. Growth Des. 2013, 13, 5353-5367. 76. Bryce, D. L.; Sward, G. D.; Adiga, S., Solid-state 35/37Cl NMR spectroscopy of hydrochloride salts of amino acids implicated in chloride ion transport channel selectivity: opportunities at 900 MHz. J. Am. Chem. Soc. 2006, 128, 2121-2134. 77. Bryce, D. L.; Gee, M.; Wasylishen, R. E., High-field chlorine NMR spectroscopy of solid organic hydrochloride salts: a sensitive probe of hydrogen bonding environments. J. Phys. Chem. A 2001, 105, 10413-10421. 78. Bryce, D. L.; Sward, G. D., Chlorine-35/37 NMR spectroscopy of solid amino acid hydrochlorides: refinement of hydrogen-bonded proton positions using experiment and theory. J. Phys. Chem. B 2006, 110, 26461-26470. 79. Chapman, R. P.; Bryce, D. L., A high-field solid-state 35/37Cl NMR and quantum chemical investigation of the chlorine quadrupolar and chemical shift tensors in amino acid hydrochlorides. Phys. Chem. Chem. Phys. 2007, 9, 6219-6230. 80. Chapman, R. P.; Hiscock, J. R.; Gale, P. A.; Bryce, D. L., A solid-state 35/37Cl NMR study of a chloride ion receptor and a GIPAW-DFT study of chlorine NMR interaction tensors in organic hydrochlorides. Can. J. Chem. 2011, 89, 822-834.

26 ACS Paragon Plus Environment

Page 26 of 39

Page 27 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

81. Attrell, R. J.; Widdifield, C. M.; Korobkov, I.; Bryce, D. L., Weak halogen bonding in solid haloanilinium halides probed directly via chlorine-35, bromine-81, and iodine-127 NMR spectroscopy. Cryst. Growth Des. 2012, 12, 1641-1653. 82. Bighley, L. D. B.; Berge, S. M.; Monkhouse, D. C., In Encyclopedia of Pharmaceutical Technology, Swarbrick, J. B.; Boylan, J. C., Eds. Marcel Dekker: 1995; Vol. 13, p 453. 83. Socha, O.; Hodgkinson, P.; Widdifield, C. M.; Yates, J. R.; Dracinsky, M., Exploring systematic discrepancies in DFT calculations of chlorine nuclear quadrupole couplings. J. Phys. Chem. A 2017. 84. Dracinsky, M.; Hodgkinson, P., A molecular dynamics study of the effects of fast molecular motions on solid-state NMR parameters. CrystEngComm 2013, 15, 8705-8712. 85. Holmes, S. T.; Iuliucci, R. J.; Mueller, K. T.; Dybowski, C., Semi-empirical refinements of crystal structures using 17O quadrupolar-coupling tensors. J. Chem. Phys. 2017, 146, 064201. 86. Grimme, S., Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787-1799. 87. Jurečka, P.; Černý, J.; Hobza, P.; Salahub, D. R., Density functional theory augmented with an empirical dispersion term. interaction energies and geometries of 80 noncovalent complexes compared with ab initio quantum mechanics calculations. J. Comput. Chem. 2007, 28, 555-569. 88. Civalleri, B.; Zicovich-Wilson, C. M.; Valenzano, L.; Ugliengo, P., B3LYP augmented with an empirical dispersion term (B3LYP-D*) as applied to molecular crystals. CrystEngComm 2008, 10, 405410. 89. Strutyński, K.; Melle-Franco, M.; Gomes, J. A. N. F., New parameterization scheme of DFT-D for graphitic materials. J. Phys. Chem. A 2013, 117, 2844-2853. 90. Clark, S. J.; Segall, M. D.; Pickard, C. J.; Hasnip, P. J.; Probert, M. J.; Refson, K.; Payne, M. C., First principles methods using CASTEP. Z. Kristallogr. 2005, 220, 567-570. 91. Profeta, M.; Mauri, F.; Pickard, C. J., Accurate first principles prediction of 17O NMR parameters in SiO2: assignment of the zeolite ferrierite spectrum. J. Am. Chem. Soc. 2003, 125, 541-548. 92. Pickard, C. J.; Mauri, F., All-electron magnetic response with pseudopotentials: NMR chemical shifts. Phys. Rev. B 2001, 63, 245101. 93. Hammer, B.; Hansen, L. B.; Nørskov, J. K., Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 1999, 59, 7413-7421. 94. Holmes, S. T.; Alkan, F.; Iuliucci, R. J.; Mueller, K. T.; Dybowski, C., Analysis of the bond-valence method for calculating 29Si and 31P magnetic shielding in covalent network solids. J. Comput. Chem. 2016, 37, 1704-1710. 95. Holmes, S. T.; Bai, S.; Iuliucci, R. J.; Mueller, K. T.; Dybowski, C., Calculations of solid-state 43Ca NMR parameters: a comparison of periodic and cluster approaches and an evaluation of DFT functionals. J. Comput. Chem. 2017, 38, 949-956. 96. Yates, J. R.; Pickard, C. J.; Mauri, F., Calculation of NMR chemical shifts for extended systems using ultrasoft pseudopotentials. Phys. Rev. B 2007, 76. 97. Pfrommer, B. G.; Côté, M.; Louie, S. G.; Cohen, M. L., Relaxation of crystals with the quasi-Newton method. J. Comput. Phys. 1997, 131, 233-240. 98. McNellis, E. R.; Meyer, J.; Reuter, K., Azobenzene at coinage metal surfaces: role of dispersive van der Waals interactions. Phys. Rev. B 2009, 80, 205414. 99. Pyykkö, P., Spectroscopic nuclear quadrupole moments. Mol. Phys. 2001, 99, 1617-1629. 100. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H., A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. 101. Chai, J.-D.; Head-Gordon, M., Long-range corrected hybrid density functionals with damped atomatom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615-6620. 102. Grimme, S., Density functional theory with London dispersion corrections. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 1, 211-228. 103. van de Streek, J.; Neumann, M. A., Validation of experimental molecular crystal structures with dispersion-corrected density functional theory calculations. Acta Crystallogr., Sect. B 2010, 66, 544-558. 27 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

104. Al-Karaghouli, A. R.; Cole, F. E.; Lehmann, M. S.; Miskell, C. R.; Verbist, J. J.; Koetzle, T. F., Precision neutron diffraction structure determination of protein and nucleic acid components. xvii. molecular and crystal structure of the amino acid glycine hydrochloride. J. Chem. Phys. 1975, 63, 13601366. 105. Koetzle, T. F.; Golic, L.; Lehmann, M. S.; Verbist, J. J.; Hamilton, W. C., Precision neutron diffraction structure determination of protein and nucleic acid components. xv. crystal and molecular structure of the amino acid L-valine hydrochloride. J. Chem. Phys. 1974, 60, 4690-4696. 106. Sequeira, A.; Rajagopal, H.; Chidambaram, R., A neutron diffraction study of the structure of Lglutamic acid·HCl. Acta Crystallogr., Sect. B 1972, 28, 2514-2519. 107. Koetzle, T. F.; Lehmann, M. S.; Verbist, J. J.; Hamilton, W. C., Precision neutron diffraction structure determination of protein and nucleic acid components. vii. the crystal and molecular structure of the amino acid L-lysine monohydrochloride dihydrate. Acta Crystallogr., Sect. B 1972, 28, 3207-3214. 108. Al-Karaghouli, A. R.; Koetzle, T. F., Neutron diffraction study of L-phenylalanine hydrochloride. Acta Crystallogr., Sect. B 1975, 31, 2461-2465. 109. Fuess, H.; Hohlwein, D.; Mason, S. A., Neutron diffraction study of L-histidine hydrochloride monohydrate. Acta Crystallogr., Sect. B 1977, 33, 654-659. 110. Kristian, B.; Valentino, R. C.; Kyuho, L.; Elsebeth, S.; Thonhauser, T.; Per, H.; Bengt, I. L., van der Waals forces in density functional theory: a review of the vdW-DF method. Rep. Prog. Phys. 2015, 78, 066501.

28 ACS Paragon Plus Environment

Page 28 of 39

Page 29 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 1. Statistical data associated with the prediction of EFG tensors for the training set organic solids, as calculated on XRD-derived structures and on rPBE, rPBE-D2, or rPBE-D2* refined structures. a Structure rmsre b rmse (MHz) c max(+) (MHz) d max(-) (MHz) e 14 N XRD 0.267 0.33 0.88 -0.21 rPBE 0.112 0.17 0.40 -0.08 rPBE-D2 0.089 0.15 0.44 -0.07 rPBE-D2* 0.077 0.10 0.21 -0.11 17 O XRD 0.071 0.54 1.34 -0.80 rPBE 0.055 0.44 0.79 -0.72 rPBE-D2 0.054 0.42 0.76 -0.74 rPBE-D2* 0.036 0.28 0.55 -0.39 35 Cl XRD 0.364 1.86 5.49 -3.66 rPBE 0.314 1.65 3.54 -0.43 rPBE-D2 0.307 1.63 3.56 -0.57 rPBE-D2* 0.055 0.23 0.61 -0.46 All XRD 0.233 rPBE 0.175 rPBE-D2 0.168 rPBE-D2* 0.053 a All calculations were performed using the rPBE plane wave DFT method, as implemented in the CASTEP module of Materials Studio. b Root mean square relative error. c Root mean square error. d Maximum positive error. e Maximum negative error.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 2. The shortest H⋯Cl- contacts for several amino acid hydrochlorides as determined via energy minimization and geometry optimization with DFT plane wave calculations. H···Cl- contacts (Å) b XRD rPBE rPBE-D2* Glycine HCl ROOH···Cl 1.974 1.981 2.071 + RNH3 ···Cl 2.159 2.114 2.186 RNH3+···Cl2.264 2.119 2.201 + RNH3 ···Cl 2.721 2.619 2.604 L-alanine HCl ROOH···Cl2.043 1.967 2.074 RNH3+···Cl2.205 2.126 2.197 RNH3+···Cl2.227 2.138 2.208 L-valine HCl ROOH···Cl 2.103 1.969 2.071 RNH3+···Cl2.277 2.153 2.238 + RNH3 ···Cl 2.325 2.220 2.261 RNH3+···Cl2.393 2.353 2.414 ROOH···Cl 2.206 1.986 2.093 L-cysteine HCl·H2O HOH···Cl2.394 2.157 2.224 HOH···Cl 2.395 2.237 2.304 RNH3+···Cl2.536 2.356 2.348 RNH3+···Cl2.165 2.142 2.215 L-histidine HCl·H2O HOH···Cl2.389 2.157 2.263 + RNH3 ···Cl 2.404 2.210 2.291 a Indicates the functional group contributing to the H⋯Cl contacts (i.e., ROOH⋯Cl-, RNH3+⋯Cl-, and HOH⋯Cl- denote carboxylic-, positively charged primary ammonia-, and water-type hydrogen contacts, respectively). b The shortest H⋯Cl- contacts (< 2.64 Å) as determined via XRD or geometry optimization at the rPBE and rPBE-D2* levels. Material

Contact type a

30 ACS Paragon Plus Environment

Page 30 of 39

Page 31 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 3. Statistical data associated with the prediction of EFG tensors for the test set organic solids, as calculated on XRD-derived structures, and on rPBE or rPBE-D2* refined structures. Structure rmse (MHz) b max(+) (MHz) c max(-) (MHz) d 14 N XRD 0.33 0.92 -1.06 rPBE 0.17 0.46 -0.61 rPBE-D2* 0.15 0.42 -0.44 17 O XRD 0.58 2.01 -2.50 rPBE 0.49 1.39 -1.15 rPBE-D2* 0.36 1.09 -1.33 35 Cl XRD 1.51 5.49 -3.66 rPBE 1.30 3.54 -0.99 rPBE-D2* 0.40 0.80 -1.21 a b c Root mean square relative error. Root mean square error. Maximum positive error. d Maximum negative error.

31 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 39

Table 4. Root mean square errors in calculated EFG tensor principal elements based on model structures obtained from XRD, or structures refined at the rPBE and rPBE-D2* levels. Structure

All principal values (MHz)

eQV₁₁/h (MHz)

eQV₂₂/h (MHz)

eQV₃₃/h (MHz)

ηQ

0.33 0.17 0.15

0.23 0.14 0.14

0.36 0.17 0.11

0.39 0.19 0.19

0.17 0.10 0.07

0.58 0.49 0.36

0.66 0.44 0.34

0.33 0.58 0.32

0.67 0.45 0.41

0.09 0.09 0.05

1.51 1.30 0.40

0.92 0.57 0.32

1.42 1.33 0.40

1.99 1.73 0.46

0.19 0.10 0.10

14

N XRD rPBE rPBE-D2* 17 O XRD rPBE rPBE-D2* 35 Cl XRD rPBE rPBE-D2*

32 ACS Paragon Plus Environment

Page 33 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 5. Experimental and calculated values of eQV22/h for all pseudotetrahedral 14N sites with 0.31 ≤ ηQ ≤ 0.70. Material L-isoleucine L-leucine L-lysine HCl·2H2O L-methionine L-threonine L-valine α-glycine β-glycine γ-glycine

Exp (MHz) -0.74 -0.74 -0.71 -0.87 -0.76 -0.74 -0.90 -0.90 -0.82

X-ray (MHz) -1.00 -1.23 -0.78 -1.68 -1.31 -1.11 -0.96 -1.17 -0.95

rPBE (MHz) -1.06 -1.10 -0.93 -1.16 -0.85 -1.03 -1.10 -1.13 -0.97

rmse (MHz) a 0.41 0.25 max(+) (MHz) b 0.82 0.36 max(-) (MHz) c a b c Root mean square error. Maximum positive error. Maximum negative error.

33 ACS Paragon Plus Environment

rPBE-D2* (MHz) -0.79 -0.83 -0.71 -0.86 -0.84 -0.78 -0.91 -0.87 -0.83 0.04 0.09 -0.03

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 39

Table 6. Experimental and calculated 35Cl EFG tensor parameters for six HCl salts, as determined for several types of structural data. Material

Glycine HCl L-valine HCl L-glutamic acid HCl L-histidine HCl·H2O L-lysine HCl·2H2O L-phenylalanine HCl

Exp CQ (MHz) 6.42 5.89 3.61 1.80 2.49 6.08

ηQ 0.61 0.51 0.65 0.72 0.42 0.52

rPBE CQ (MHz) -9.13 -8.22 5.47 2.77 4.23 -8.45

ηQ 0.82 0.53 0.69 0.60 0.50 0.43

Neutron CQ ηQ (MHz) -8.35 0.82 -8.15 0.54 3.89 0.30 2.71 0.64 3.45 0.41 -7.92 0.54

rmse (MHz) a 1.63 1.24 max(+) (MHz) b 3.14 2.43 max(-) (MHz) c -0.43 -0.50 a b c Root mean square error. Maximum positive error. Maximum negative error.

34 ACS Paragon Plus Environment

rPBE-D2* CQ ηQ (MHz) -6.57 0.76 -5.76 0.59 3.72 0.39 1.65 0.69 3.05 0.68 -6.07 0.48 0.35 0.80 -0.46

Page 35 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 7. The shortest H⋯Cl- contacts for several amino acid hydrochlorides as determined by neutron diffraction or via energy minimization and geometry optimization at the rPBE-D2* level. H···Cl contacts (Å) b rPBE Neutron rPBE-D2* Glycine HCl ROOH···Cl 1.981 2.008 2.071 RNH3+···Cl2.114 2.123 2.186 RNH3+···Cl2.119 2.160 2.201 RNH3+···Cl2.619 2.593 2.604 L-valine HCl ROOH···Cl 1.969 1.990 2.071 RNH3+···Cl2.153 2.161 2.238 RNH3+···Cl2.220 2.263 2.261 RNH3+···Cl2.353 2.355 2.414 L-glutamic acid HCl ROOH···Cl 2.008 2.073 2.118 RNH3+···Cl2.068 2.107 2.137 RNH3+···Cl2.154 2.137 2.190 + RNH3 ···Cl 2.142 2.215 2.165 L-histidine HCl·H₂O HOH···Cl 2.157 2.233 2.263 RNH3+···Cl2.210 2.265 2.291 L-lysine HCl·2H2O RNH3+···Cl2.103 2.150 2.198 RNH3+···Cl2.211 2.196 2.262 HOH···Cl2.211 2.273 2.278 HOH···Cl 2.267 2.329 2.384 L-phenylalanine HCl ROOH···Cl1.941 2.004 2.065 RNH3+···Cl2.090 2.142 2.201 RNH3+···Cl2.300 2.302 2.368 + RNH3 ···Cl 2.304 2.361 2.385 a Indicates the functional group contributing to the H⋯Cl contacts (i.e., ROOH⋯Cl-, RNH3+⋯Cl-, and HOH⋯Cl- denote carboxylic-, positively charged primary ammonia-, and water-type hydrogen contacts, respectively). b The shortest H⋯Cl- contacts (< 2.64 Å) as determined via neutron diffraction or geometry optimization at the rPBE and rPBE-D2* levels. Material

Contact type a

35 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ToC Graphic

36 ACS Paragon Plus Environment

Page 36 of 39

Page 37 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. Root mean square relative errors (rmsre) in the principal elements of computed EFG tensors of the training set organic solids as a function of the damping parameter, 2.5 ≤ d ≤ 10.0. By default, the value of d is set to 20 in CASTEP.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 39

Figure 2. Relationship between the principal elements of calculated and experimental nuclear EFG tensors. Panels illustrate the results for 35Cl sites (a - c), 17O sites (d – f), and 14N sites (g – i). Computed EFG tensor parameters are derived from calculations on XRD-derived structures (green), and structures with atomic coordinates refined at the rPBE (orange) or rPBE-D2* (blue) levels. The dashed lines represent perfect agreement between calculated and experimental values. The dotted red line in (b) represents the best fit between calculated and experimental values.

2 ACS Paragon Plus Environment

Page 39 of 39 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. Errors in calculated eQV22/h in pseudotetrahedral 14N sites with 0.31 ≤ ηQ ≤ 0.70. Calculations were performed on X-ray structures (green), rPBE structures (orange), and rPBE-D2* structures (blue). The materials are L-isoleucine (1), L-leucine (2), L-lysine HCl·2H2O (3), L-methionine (4), L-threonine (5), L-valine (6), -glycine (7), -glycine (8), and -glycine (9).

Figure 4. Root mean square error in the calculated principal elements of 35Cl EFG tensors based on structures obtained by rPBE refinement of X-ray diffraction structures (orange), neutron diffraction (grey), or rPBE-D2* refinement of X-ray diffraction structures. The materials are glycine HCl (1), L-valine HCl (2), L-glutamic acid HCl (3), L-histidine HCl·H2O (4), L-lysine HCl·2H2O (5), and L-phenylalanine HCl (6).

3 ACS Paragon Plus Environment