1
|
Ke Z, Weng J, Xu X. Calculating 13 C NMR chemical shifts of large molecules using the eXtended ONIOM method at high accuracy with a low cost. J Comput Chem 2023; 44:2347-2357. [PMID: 37572044 DOI: 10.1002/jcc.27201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Revised: 07/15/2023] [Accepted: 07/24/2023] [Indexed: 08/14/2023]
Abstract
Fragmentation-based methods for nuclear magnetic resonance (NMR) chemical shift calculations have become more and more popular in first-principles calculations of large molecules. However, there are many options for a fragmentation-based method to select, such as theoretical methods, fragmentation schemes, the number of levels of theory, etc. It is important to study the optimal combination of the options to achieve a good balance between accuracy and efficiency. Here we investigate different combinations of options used by a fragmentation-based method, the eXtended ONIOM (XO) method, for 13 C chemical shift calculations on a set of organic and biological molecules. We found that: (1) introducing Hartree-Fock exchange into density functional theory (DFT) could reduce the calculation error due to fragmentation in contrast to pure DFT functionals, while a hybrid functional, xOPBE, is generally recommended; (2) fragmentation schemes generated from the molecular tailoring approach (MTA) with small level parameter n, for example, n = 2 and the degree-based fragmentation method (DBFM) with n = 1, are sufficient to achieve satisfactory accuracy; (3) the two-level XO (XO2) NMR calculation is superior to the calculation with only one level of theory, as the second level (i.e., low level) of theory provides a way to well describe the long-range effect. These findings are beneficial to practical applications of fragmentation-based methods for NMR chemical shift calculations of large molecules.
Collapse
Affiliation(s)
- Zhipeng Ke
- Institute of Photochemistry and Photofunctional Materials, University of Shanghai for Science and Technology, Shanghai, China
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, China
| | - Jingwei Weng
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, China
- Hefei National Laboratory, Hefei, China
| |
Collapse
|
2
|
Zhang J, Kriebel CN, Wan Z, Shi M, Glaubitz C, He X. Automated Fragmentation Quantum Mechanical Calculation of 15N and 13C Chemical Shifts in a Membrane Protein. J Chem Theory Comput 2023; 19:7405-7422. [PMID: 37788419 DOI: 10.1021/acs.jctc.3c00621] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
In this work, we developed an accurate and cost-effective automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method to calculate the chemical shifts of 15N and 13C of membrane proteins. The convergence of the AF-QM/MM method was tested using Krokinobacter eikastus rhodopsin 2 as a test case. When the distance threshold of the QM region is equal to or larger than 4.0 Å, the results of the AF-QM/MM calculations are close to convergence. In addition, the effects of selected density functionals, basis sets, and local chemical environment of target atoms on the chemical shift calculations were systematically investigated. Our results demonstrate that the predicted chemical shifts are more accurate when important environmental factors including cross-protomer interactions, lipid molecules, and solvent molecules are taken into consideration, especially for the 15N chemical shift prediction. Furthermore, with the presence of sodium ions in the environment, the chemical shift of residues, retinal, and retinal Schiff base are affected, which is consistent with the results of the solid-state nuclear magnetic resonance (NMR) experiment. Upon comparing the performance of various density functionals (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB95, and OPBE), the results show that mPW1PW91 is a suitable functional for the 15N and 13C chemical shift prediction of the membrane proteins. Meanwhile, we find that the improved accuracy of the 13Cβ chemical shift calculations can be achieved by the employment of the triple-ζ basis set. However, the employment of the triple-ζ basis set does not improve the accuracy of the 15N and 13Cα chemical shift calculations nor does the addition of a diffuse function improve the overall prediction accuracy of the chemical shifts. Our study also underscores that the AF-QM/MM method has significant advantages in predicting the chemical shifts of key ligands and nonstandard residues in membrane proteins than most widely used empirical models; therefore, it could be an accurate computational tool for chemical shift calculations on various types of biological systems.
Collapse
Affiliation(s)
- Jinhuan Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Clara Nassrin Kriebel
- Institute of Biophysical Chemistry and Center for Biomolecular Magnetic Resonance, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
| | - Zheng Wan
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Man Shi
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Clemens Glaubitz
- Institute of Biophysical Chemistry and Center for Biomolecular Magnetic Resonance, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
3
|
Poidevin C, Stoychev GL, Riplinger C, Auer AA. High Level Electronic Structure Calculation of Molecular Solid-State NMR Shielding Constants. J Chem Theory Comput 2022; 18:2408-2417. [PMID: 35353527 PMCID: PMC9009078 DOI: 10.1021/acs.jctc.1c01095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
In this work, we
present a quantum mechanics/molecular mechanics
(QM/MM) approach for the computation of solid-state nuclear magnetic
resonance (SS-NMR) shielding constants (SCs) for molecular crystals.
Besides applying standard-DFT functionals like GGAs (PBE), meta-GGAs
(TPSS), and hybrids (B3LYP), we apply a double-hybrid (DSD-PBEP86)
functional as well as MP2, using the domain-based local pair natural
orbital (DLPNO) formalism, to calculate the NMR SCs of six amino acid
crystals. All the electronic structure methods used exhibit good correlation
of the NMR shieldings with respect to experimental chemical shifts
for both 1H and 13C. We also find that local
electronic structure is much more important than the long-range electrostatic
effects for these systems, implying that cluster approaches using
all-electron/Gaussian basis set methods might offer great potential
for predictive computations of solid-state NMR parameters for organic
solids.
Collapse
Affiliation(s)
- Corentin Poidevin
- Institut des Sciences Chimiques de Rennes, Av. Général Leclerc, 357000 Rennes, France
| | - Georgi L Stoychev
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | | | - Alexander A Auer
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
4
|
Li W, Dong H, Ma J, Li S. Structures and Spectroscopic Properties of Large Molecules and Condensed-Phase Systems Predicted by Generalized Energy-Based Fragmentation Approach. Acc Chem Res 2021; 54:169-181. [PMID: 33350806 DOI: 10.1021/acs.accounts.0c00580] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
ConspectusThe structures and spectroscopic properties of molecules and condensed-phase systems are usually experimentally characterized by X-ray, infrared (IR), Raman, nuclear magnetic resonance (NMR), and electronic absorption/emission spectra. Quantum mechanics (QM) calculations are critical in quantitatively understanding the relationship between the structure and physicochemical properties of various chemical systems. However, it is very challenging to apply traditional QM methods to large molecules and condensed-phase systems with large unit cells due to their steep computational scaling with the system size. To overcome this difficulty, theoretical chemists have developed various linear (or low) scaling QM methods, among which energy-based fragmentation methods have achieved great success for large molecules or clusters. One of the most popular energy-based fragmentation methods is the generalized energy-based fragmentation (GEBF) approach developed by us.In this approach, the ground-state energy of a large molecule can be evaluated from the ground-state energies of a series of embedded subsystems. In this Account, we focus on the recent developments and applicability of the GEBF approach for the structures and spectroscopic properties of complicated large molecules and condensed-phase systems. With new fragmentation schemes, the GEBF approach can now describe ionic liquid clusters and metal-containing supramolecular systems accurately and can provide accurate binding energies for host-guest complexes. In addition, the GEBF approach is now available for describing the localized excited states of large systems including a chromophore. More importantly, the GEBF approach under periodic boundary conditions (PBC-GEBF) has been developed to deal with periodic molecular crystals and liquids. Then, the ground-state energy (or property) per unit cell of a periodic condensed phase system can be predicted with QM calculations on nonperiodic embedded subsystems. This feature enables accurate electron correlation calculations on molecular crystals and liquids to be feasible on ordinary workstations. The PBC-GEBF approach has been applied to predict the crystal structures, lattice energies, and spectroscopic properties of some typical molecular crystals and solutions. By combining the GEBF method and machine learning (ML) method, a GEBF-ML force field has been developed for long normal alkanes, and the IR spectra of long alkanes can be obtained from the GEBF-ML molecular dynamics (MD) simulations. The GEBF and its periodic variant are expected to play increasingly important roles in investigating real-life chemical systems of broad interests at the ab initio levels.
Collapse
Affiliation(s)
- Wei Li
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Hao Dong
- Kuang Yaming Honors School and Institute for Brain Sciences, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Jing Ma
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Shuhua Li
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| |
Collapse
|
5
|
Unzueta PA, Beran GJO. Polarizable continuum models provide an effective electrostatic embedding model for fragment-based chemical shift prediction in challenging systems. J Comput Chem 2020; 41:2251-2265. [PMID: 32748418 DOI: 10.1002/jcc.26388] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 06/04/2020] [Accepted: 07/04/2020] [Indexed: 12/25/2022]
Abstract
Ab initio nuclear magnetic resonance chemical shift prediction provides an important tool for interpreting and assigning experimental spectra, but it becomes computationally prohibitive in large systems. The computational costs can be reduced considerably by fragmentation of the large system into a series of contributions from many smaller subsystems. However, the presence of charged functional groups and the need to partition the system across covalent bonds create complications in biomolecules that typically require the use of large fragments and careful descriptions of the electrostatic environment. The present work shows how a model that combines chemical shielding contributions from non-overlapping monomer and dimer fragments embedded in a polarizable continuum model provides a simple, easy-to-implement, and computationally inexpensive approach for predicting chemical shifts in complex systems. The model's performance proves rather insensitive to the continuum dielectric constant, making the selection of the optimal embedding dielectric less critical. The PCM-embedded fragment model is demonstrated to perform well across systems ranging from molecular crystals to proteins.
Collapse
Affiliation(s)
- Pablo A Unzueta
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| | - Gregory J O Beran
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| |
Collapse
|
6
|
Abstract
Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
7
|
Lino JBDR, Ramalho TC. Exploring Through-Space Spin-Spin Couplings for Quantum Information Processing: Facing the Challenge of Coherence Time and Control Quantum States. J Phys Chem A 2019; 123:1372-1379. [PMID: 30673241 DOI: 10.1021/acs.jpca.8b09425] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Nuclear magnetic resonance (NMR) is a powerful tool for studying quantum information processing (QIP). Recently quantum technologies have been proposed to overcome the challenges in large-scale NMR QIP. Furthermore, computational chemistry can promote its improvement. Nuclear spins-1/2 are natural qubits and have been used in most NMR quantum computation experiments. However, molecules that enable many qubits NMR QIP implementations should meet some requirements regarding their spectroscopic properties. Exceptionally large through-space (TS) P-P spin-spin coupling constants (SSCC or J) observed in 1,8-diphosphanaphthalenes (PPN) and in naphtho[1,8- cd]-1,2-dithiole phenylphosphines (NTP) were proposed and investigated to provide more accurate control within large-scale NMR QIP. Spectroscopic properties of PPN and NTP derivatives were explored by theoretical strategies using locally dense basis sets (LDBS). 31P chemical shifts (δ) calculated at the B3LYP/aug-cc-pVTZ-J level and TS P-P SSCCs at the PBE1PBE/pcJ-2 (LDBS-1) level are very close to the experimental data for the PPN molecule. Differently, for the NTP dimer, PBE1PBE/pcJ-2 (LDBS-2) predicts more accurate 31P δ, whereas PBE1PBE/Def2-TZVP (LDBS-1) forecasts more accurate TS P-P SSCCs. From our results, PPNo-F, PPNo-ethyl, and PPNo-NH2 were the best candidates for NMR QIP, in which the large TS SSCCS could face the need of long-time quantum gates implementations. Therefore, it could overcome natural limitations concerning the development of large-scale NMR.
Collapse
Affiliation(s)
| | - Teodorico Castro Ramalho
- Chemistry Department , Federal University of Lavras , 37200-000 Lavras , MG Brazil.,Center for Basic and Applied Research, Faculty of Informatics and Management , University Hradec Kralove , 50003 Hradec Kralove , Czech Republic
| |
Collapse
|
8
|
Kobayashi R, Amos RD, Reid DM, Collins MA. Application of the Systematic Molecular Fragmentation by Annihilation Method to ab Initio NMR Chemical Shift Calculations. J Phys Chem A 2018; 122:9135-9141. [DOI: 10.1021/acs.jpca.8b09565] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Rika Kobayashi
- ANU Supercomputer Facility, Leonard Huxley Building 56, Mills Road, Canberra, ACT 2601, Australia
| | - Roger D. Amos
- ANU Supercomputer Facility, Leonard Huxley Building 56, Mills Road, Canberra, ACT 2601, Australia
| | - David M. Reid
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Michael A. Collins
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
9
|
Zhao D, Song R, Li W, Ma J, Dong H, Li S. Accurate Prediction of NMR Chemical Shifts in Macromolecular and Condensed-Phase Systems with the Generalized Energy-Based Fragmentation Method. J Chem Theory Comput 2017; 13:5231-5239. [DOI: 10.1021/acs.jctc.7b00380] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dongbo Zhao
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Ruiheng Song
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Wei Li
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Jing Ma
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Hao Dong
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Shuhua Li
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| |
Collapse
|
10
|
Beran GJO, Hartman JD, Heit YN. Predicting Molecular Crystal Properties from First Principles: Finite-Temperature Thermochemistry to NMR Crystallography. Acc Chem Res 2016; 49:2501-2508. [PMID: 27754668 DOI: 10.1021/acs.accounts.6b00404] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Molecular crystals occur widely in pharmaceuticals, foods, explosives, organic semiconductors, and many other applications. Thanks to substantial progress in electronic structure modeling of molecular crystals, attention is now shifting from basic crystal structure prediction and lattice energy modeling toward the accurate prediction of experimentally observable properties at finite temperatures and pressures. This Account discusses how fragment-based electronic structure methods can be used to model a variety of experimentally relevant molecular crystal properties. First, it describes the coupling of fragment electronic structure models with quasi-harmonic techniques for modeling the thermal expansion of molecular crystals, and what effects this expansion has on thermochemical and mechanical properties. Excellent agreement with experiment is demonstrated for the molar volume, sublimation enthalpy, entropy, and free energy, and the bulk modulus of phase I carbon dioxide when large basis second-order Møller-Plesset perturbation theory (MP2) or coupled cluster theories (CCSD(T)) are used. In addition, physical insight is offered into how neglect of thermal expansion affects these properties. Zero-point vibrational motion leads to an appreciable expansion in the molar volume; in carbon dioxide, it accounts for around 30% of the overall volume expansion between the electronic structure energy minimum and the molar volume at the sublimation point. In addition, because thermal expansion typically weakens the intermolecular interactions, neglecting thermal expansion artificially stabilizes the solid and causes the sublimation enthalpy to be too large at higher temperatures. Thermal expansion also frequently weakens the lower-frequency lattice phonon modes; neglecting thermal expansion causes the entropy of sublimation to be overestimated. Interestingly, the sublimation free energy is less significantly affected by neglecting thermal expansion because the systematic errors in the enthalpy and entropy cancel somewhat. Second, because solid state nuclear magnetic resonance (NMR) plays an increasingly important role in molecular crystal studies, this Account discusses how fragment methods can be used to achieve higher-accuracy chemical shifts in molecular crystals. Whereas widely used plane wave density functional theory models are largely restricted to generalized gradient approximation (GGA) functionals like PBE in practice, fragment methods allow the routine use of hybrid density functionals with only modest increases in computational cost. In extensive molecular crystal benchmarks, hybrid functionals like PBE0 predict chemical shifts with 20-30% higher accuracy than GGAs, particularly for 1H, 13C, and 15N nuclei. Due to their higher sensitivity to polarization effects, 17O chemical shifts prove slightly harder to predict with fragment methods. Nevertheless, the fragment model results are still competitive with those from GIPAW. The improved accuracy achievable with fragment approaches and hybrid density functionals increases discrimination between different potential assignments of individual shifts or crystal structures, which is critical in NMR crystallography applications. This higher accuracy and greater discrimination are highlighted in application to the solid state NMR of different acetaminophen and testosterone crystal forms.
Collapse
Affiliation(s)
- Gregory J. O. Beran
- Department of Chemistry, University of California, Riverside, California 92521, United States
| | - Joshua D. Hartman
- Department of Chemistry, University of California, Riverside, California 92521, United States
| | - Yonaton N. Heit
- Department of Chemistry, University of California, Riverside, California 92521, United States
| |
Collapse
|
11
|
Collins MA. Can Systematic Molecular Fragmentation Be Applied to Direct Ab Initio Molecular Dynamics? J Phys Chem A 2016; 120:9281-9291. [DOI: 10.1021/acs.jpca.6b08739] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michael A. Collins
- Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia
| |
Collapse
|
12
|
Hartman J, Day GM, Beran GJO. Enhanced NMR Discrimination of Pharmaceutically Relevant Molecular Crystal Forms through Fragment-Based Ab Initio Chemical Shift Predictions. CRYSTAL GROWTH & DESIGN 2016; 16:6479-6493. [PMID: 27829821 PMCID: PMC5095663 DOI: 10.1021/acs.cgd.6b01157] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2016] [Revised: 09/09/2016] [Indexed: 05/10/2023]
Abstract
Chemical shift prediction plays an important role in the determination or validation of crystal structures with solid-state nuclear magnetic resonance (NMR) spectroscopy. One of the fundamental theoretical challenges lies in discriminating variations in chemical shifts resulting from different crystallographic environments. Fragment-based electronic structure methods provide an alternative to the widely used plane wave gauge-including projector augmented wave (GIPAW) density functional technique for chemical shift prediction. Fragment methods allow hybrid density functionals to be employed routinely in chemical shift prediction, and we have recently demonstrated appreciable improvements in the accuracy of the predicted shifts when using the hybrid PBE0 functional instead of generalized gradient approximation (GGA) functionals like PBE. Here, we investigate the solid-state 13C and 15N NMR spectra for multiple crystal forms of acetaminophen, phenobarbital, and testosterone. We demonstrate that the use of the hybrid density functional instead of a GGA provides both higher accuracy in the chemical shifts and increased discrimination among the different crystallographic environments. Finally, these results also provide compelling evidence for the transferability of the linear regression parameters mapping predicted chemical shieldings to chemical shifts that were derived in an earlier study.
Collapse
Affiliation(s)
- Joshua
D. Hartman
- Department
of Chemistry, University of California, Riverside, California 92521 United States
| | - Graeme M. Day
- School
of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gregory J. O. Beran
- Department
of Chemistry, University of California, Riverside, California 92521 United States
- E-mail:
| |
Collapse
|
13
|
Amos R, Kobayashi R. Ab Initio NMR Chemical Shift Calculations Using Fragment Molecular Orbitals and Locally Dense Basis Sets. J Phys Chem A 2016; 120:8907-8915. [DOI: 10.1021/acs.jpca.6b09190] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Roger Amos
- Australian National University, Leonard Huxley Bldg 56, Mills Rd, Canberra, ACT 2601, Australia
| | - Rika Kobayashi
- Australian National University, Leonard Huxley Bldg 56, Mills Rd, Canberra, ACT 2601, Australia
| |
Collapse
|
14
|
Hartman JD, Kudla RA, Day GM, Mueller LJ, Beran GJO. Benchmark fragment-based (1)H, (13)C, (15)N and (17)O chemical shift predictions in molecular crystals. Phys Chem Chem Phys 2016; 18:21686-709. [PMID: 27431490 DOI: 10.1039/c6cp01831a] [Citation(s) in RCA: 78] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The performance of fragment-based ab initio(1)H, (13)C, (15)N and (17)O chemical shift predictions is assessed against experimental NMR chemical shift data in four benchmark sets of molecular crystals. Employing a variety of commonly used density functionals (PBE0, B3LYP, TPSSh, OPBE, PBE, TPSS), we explore the relative performance of cluster, two-body fragment, and combined cluster/fragment models. The hybrid density functionals (PBE0, B3LYP and TPSSh) generally out-perform their generalized gradient approximation (GGA)-based counterparts. (1)H, (13)C, (15)N, and (17)O isotropic chemical shifts can be predicted with root-mean-square errors of 0.3, 1.5, 4.2, and 9.8 ppm, respectively, using a computationally inexpensive electrostatically embedded two-body PBE0 fragment model. Oxygen chemical shieldings prove particularly sensitive to local many-body effects, and using a combined cluster/fragment model instead of the simple two-body fragment model decreases the root-mean-square errors to 7.6 ppm. These fragment-based model errors compare favorably with GIPAW PBE ones of 0.4, 2.2, 5.4, and 7.2 ppm for the same (1)H, (13)C, (15)N, and (17)O test sets. Using these benchmark calculations, a set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided and their robustness assessed using statistical cross-validation. We demonstrate the utility of these approaches and the reported scaling parameters on applications to 9-tert-butyl anthracene, several histidine co-crystals, benzoic acid and the C-nitrosoarene SnCl2(CH3)2(NODMA)2.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California, Riverside, California 92521, USA.
| | | | | | | | | |
Collapse
|
15
|
Abstract
Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
16
|
Hartman JD, Monaco S, Schatschneider B, Beran GJO. Fragment-based (13)C nuclear magnetic resonance chemical shift predictions in molecular crystals: An alternative to planewave methods. J Chem Phys 2015; 143:102809. [PMID: 26374002 DOI: 10.1063/1.4922649] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We assess the quality of fragment-based ab initio isotropic (13)C chemical shift predictions for a collection of 25 molecular crystals with eight different density functionals. We explore the relative performance of cluster, two-body fragment, combined cluster/fragment, and the planewave gauge-including projector augmented wave (GIPAW) models relative to experiment. When electrostatic embedding is employed to capture many-body polarization effects, the simple and computationally inexpensive two-body fragment model predicts both isotropic (13)C chemical shifts and the chemical shielding tensors as well as both cluster models and the GIPAW approach. Unlike the GIPAW approach, hybrid density functionals can be used readily in a fragment model, and all four hybrid functionals tested here (PBE0, B3LYP, B3PW91, and B97-2) predict chemical shifts in noticeably better agreement with experiment than the four generalized gradient approximation (GGA) functionals considered (PBE, OPBE, BLYP, and BP86). A set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided based on these benchmark calculations. Statistical cross-validation procedures are used to demonstrate the robustness of these fits.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Stephen Monaco
- The Pennsylvania State University, The Eberly Campus, 2201 University Dr, Lemont Furnace, Pennsylvania 15456, USA
| | - Bohdan Schatschneider
- The Pennsylvania State University, The Eberly Campus, 2201 University Dr, Lemont Furnace, Pennsylvania 15456, USA
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
17
|
Reid DM, Collins MA. Approximating CCSD(T) Nuclear Magnetic Shielding Calculations Using Composite Methods. J Chem Theory Comput 2015; 11:5177-81. [DOI: 10.1021/acs.jctc.5b00546] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- David M. Reid
- Research School of Chemistry, The Australian National University, Canberra, ACT 0200, Australia
| | - Michael A. Collins
- Research School of Chemistry, The Australian National University, Canberra, ACT 0200, Australia
| |
Collapse
|