1
|
Pokorný V, Touš P, Štejfa V, Růžička K, Rohlíček J, Czernek J, Brus J, Červinka C. Anisotropy, segmental dynamics and polymorphism of crystalline biogenic carboxylic acids. Phys Chem Chem Phys 2022; 24:25904-25917. [PMID: 36260017 DOI: 10.1039/d2cp03698c] [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: 06/16/2023]
Abstract
Carboxylic acids of the Krebs cycle possess invaluable biochemical significance. Still, there are severe gaps in the availability of thermodynamic and crystallographic data, as well as ambiguities prevailing in the literature on the thermodynamic characterization and polymorph ranking. Providing an unambiguous description of the structure, thermodynamics and polymorphism of their neat crystalline phases requires a complex multidisciplinary approach. This work presents results of an extensive investigation of the structural anisotropy of the thermal expansion and local dynamics within these crystals, obtained from a beneficial cooperation of NMR crystallography and ab initio calculations of non-covalent interactions. The observed structural anisotropy and spin-lattice relaxation times are traced to large spatial variations in the strength of molecular interactions in the crystal lattice, especially in the orientation of the hydrogen bonds. A completely resolved crystal structure for oxaloacetic acid is reported for the first time. Thanks to multi-instrumental calorimetric effort, this work clarifies phase behavior, determines third-law entropies of the crystals, and states definitive polymorph ranking for succinic and fumaric acids. These thermodynamic observations are then interpreted in terms of first-principles quasi-harmonic calculations of cohesive properties. A sophisticated model capturing electronic, thermal, and configurational-entropic effects on the crystal structure approaches captures the subtle Gibbs energy differences governing polymorph ranking for succinic and fumaric acids, representing another success story of computational chemistry.
Collapse
Affiliation(s)
- Václav Pokorný
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovského nám. 2, CZ-162 06 Prague 6, Czech Republic
| | - Petr Touš
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
| | - Vojtěch Štejfa
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
| | - Květoslav Růžička
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
| | - Jan Rohlíček
- Institute of Physics, Czech Academy of Sciences, Na Slovance 2, CZ-182 21 Prague 8, Czech Republic
| | - Jiří Czernek
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovského nám. 2, CZ-162 06 Prague 6, Czech Republic
| | - Jiří Brus
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovského nám. 2, CZ-162 06 Prague 6, Czech Republic
| | - Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
| |
Collapse
|
2
|
Schriber JB, Nascimento DR, Koutsoukas A, Spronk SA, Cheney DL, Sherrill CD. CLIFF: A component-based, machine-learned, intermolecular force field. J Chem Phys 2021; 154:184110. [PMID: 34241025 DOI: 10.1063/5.0042989] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Computation of intermolecular interactions is a challenge in drug discovery because accurate ab initio techniques are too computationally expensive to be routinely applied to drug-protein models. Classical force fields are more computationally feasible, and force fields designed to match symmetry adapted perturbation theory (SAPT) interaction energies can remain accurate in this context. Unfortunately, the application of such force fields is complicated by the laborious parameterization required for computations on new molecules. Here, we introduce the component-based machine-learned intermolecular force field (CLIFF), which combines accurate, physics-based equations for intermolecular interaction energies with machine-learning models to enable automatic parameterization. The CLIFF uses functional forms corresponding to electrostatic, exchange-repulsion, induction/polarization, and London dispersion components in SAPT. Molecule-independent parameters are fit with respect to SAPT2+(3)δMP2/aug-cc-pVTZ, and molecule-dependent atomic parameters (atomic widths, atomic multipoles, and Hirshfeld ratios) are obtained from machine learning models developed for C, N, O, H, S, F, Cl, and Br. The CLIFF achieves mean absolute errors (MAEs) no worse than 0.70 kcal mol-1 in both total and component energies across a diverse dimer test set. For the side chain-side chain interaction database derived from protein fragments, the CLIFF produces total interaction energies with an MAE of 0.27 kcal mol-1 with respect to reference data, outperforming similar and even more expensive methods. In applications to a set of model drug-protein interactions, the CLIFF is able to accurately rank-order ligand binding strengths and achieves less than 10% error with respect to SAPT reference values for most complexes.
Collapse
Affiliation(s)
- Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Daniel R Nascimento
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Steven A Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| |
Collapse
|
3
|
Varandas AJC. Canonical and explicitly-correlated coupled cluster correlation energies of sub-kJ mol -1 accuracy via cost-effective hybrid-post-CBS extrapolation. Phys Chem Chem Phys 2021; 23:9571-9584. [PMID: 33885095 DOI: 10.1039/d1cp00357g] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cost-effectiveness and accuracy are two basic pillars in electronic structure calculations. While cost-effectiveness enhances applicability, high accuracy is sustained when employing advanced computational tools. With the gold standard method of ab initio quantum chemistry at the focal point, canonical CCSD(T) and modern explicitly correlated CCSD(T)-F12 calculations are employed hand in hand to develop accurate hybrid post-CBS extrapolation schemes, which are validated using popular training sets involving a total of 130 molecules. By using raw valence-only calculations at CCSD(T)/VDZ and CCSD(T)/VQZ-F12 levels of theory, the novel scheme leads to the prediction of absolute energies that differ on average (-0.170 ± 0.224) kcal mol-1 from the highest affordable CCSD(T)-F12b/V(Q,5)Z-F12 extrapolations, but only (-0.048 ± 0.228) kcal mol-1 from the post-CBS extrapolated values based on CBS(D,T), CBS(D,Q) and CBS(T,Q) energies. From the cost-effectiveness standpoint, the approach is a kind of pseudo one-point extrapolation scheme since its cost is basically that of the highest-rung raw energy where it is based. Variants that imply no additional cost are also discussed, emerging h-pCBS(dt,dq)ab as the most effective. The approach can also be used with PNO-based local correlation methods that gained popularity due to allowing coupled-cluster calculations even for large molecules at reduced computational cost, namely local PNO-CCSD(T) and PNO-CCSD(T)-F12b. To gauge the approach performance, both the hydrogen molecule and the O-C2H5 torsion path of ethyl-methyl-ether, an extra molecule here considered with presupposed existence in astrophysical objects, are also studied. Additionally, the nonbonding interactions in the A24 test set are revisited per se. The results show that the title approach may be useful in high-accuracy quantum chemistry, with further improvements requiring the inclusion of contributions beyond the theory here employed such as the ones due to relativistic and nonadiabatic effects.
Collapse
Affiliation(s)
- A J C Varandas
- School of Physics and Physical Engineering, Qufu Normal University, 273165 Qufu, China
| |
Collapse
|
4
|
Kodrycka M, Patkowski K. Platinum, gold, and silver standards of intermolecular interaction energy calculations. J Chem Phys 2019; 151:070901. [DOI: 10.1063/1.5116151] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Monika Kodrycka
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| |
Collapse
|
5
|
Červinka C, Beran GJO. Towards reliable ab initio sublimation pressures for organic molecular crystals - are we there yet? Phys Chem Chem Phys 2019; 21:14799-14810. [PMID: 31225538 DOI: 10.1039/c9cp01572h] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Knowledge of molecular crystal sublimation equilibrium data is vital in many industrial processes, but this data can be difficult to measure experimentally for low-volatility species. Theoretical prediction of sublimation pressures could provide a useful supplement to experiment, but the exponential temperature dependence of sublimation (or any saturated vapor) pressure curve makes this challenging. An uncertainty of only a few percent in the sublimation enthalpy or entropy can propagate to an error in the sublimation pressure exceeding several orders of magnitude for a given temperature interval. Despite this fundamental difficulty, this paper performs some of the first ab initio predictions of sublimation pressure curves. Four simple molecular crystals (ethane, methanol, benzene, and imidazole) have been selected for a case study showing the currently achievable accuracy of quantum chemistry calculations. Fragment-based ab initio techniques and the quasi-harmonic approximation are used for calculations of cohesive and phonon properties of the crystals, while the vapor phase is treated by the ideal gas model. Ab initio sublimation pressure curves for model compounds are compared against their experimental counterparts. The computational uncertainties are estimated, weak points of the computational methodology are identified, and further improvements are proposed.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
6
|
Cardozo TM, Galliez AP, Borges I, Plasser F, Aquino AJA, Barbatti M, Lischka H. Dynamics of benzene excimer formation from the parallel-displaced dimer. Phys Chem Chem Phys 2018; 21:13916-13924. [PMID: 30570626 DOI: 10.1039/c8cp06354k] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Excimers play a key role in a variety of excited-state processes, such as exciton trapping, fluorescence quenching, and singlet-fission. The dynamics of benzene excimer formation in the first 2 ps after S1 excitation from the parallel-displaced geometry of the benzene dimer is reported here. It was simulated via nonadiabatic surface-hopping dynamics using the second-order algebraic diagrammatic construction (ADC(2)). After excitation, the benzene rings take ∼0.5-1.0 ps to approach each other in a parallel-stacked structure of the S1 minimum and stay in the excimer region for ∼0.1-0.4 ps before leaving due to excess vibrational energy. The S1-S2 gap widens considerably while the rings visit the excimer region in the potential energy surface. Our work provides detailed insight into correlations between nuclear and electronic structure in the excimer and shows that decreased ring distance goes along with enhanced charge transfer and that fast exciton transfer happens between the rings, leading to the equal probability of finding the exciton in each ring after around 1.0 ps.
Collapse
Affiliation(s)
- Thiago Messias Cardozo
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ 21941-909, Brazil.
| | - Andre Pessoa Galliez
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ 21941-909, Brazil.
| | - Itamar Borges
- Departamento de Química, Instituto Militar de Engenharia, Praça General Tibúrcio, 80, 22290-270 Rio de Janeiro, Brazil
| | - Felix Plasser
- Department of Chemistry, Loughborough University, LE11 3TU, UK
| | - Adelia J A Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P. R. China and Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA and Institute of Theoretical Chemistry, University of Vienna, Währingerstraße 17, 1090 Vienna, Austria.
| | | | - Hans Lischka
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P. R. China and Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA and Institute of Theoretical Chemistry, University of Vienna, Währingerstraße 17, 1090 Vienna, Austria.
| |
Collapse
|
7
|
Przybytek M. Dispersion Energy of Symmetry-Adapted Perturbation Theory from the Explicitly Correlated F12 Approach. J Chem Theory Comput 2018; 14:5105-5117. [DOI: 10.1021/acs.jctc.8b00470] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michał Przybytek
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
8
|
Faber FA, Christensen AS, Huang B, von Lilienfeld OA. Alchemical and structural distribution based representation for universal quantum machine learning. J Chem Phys 2018; 148:241717. [PMID: 29960351 DOI: 10.1063/1.5020710] [Citation(s) in RCA: 204] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce a representation of any atom in any chemical environment for the automatized generation of universal kernel ridge regression-based quantum machine learning (QML) models of electronic properties, trained throughout chemical compound space. The representation is based on Gaussian distribution functions, scaled by power laws and explicitly accounting for structural as well as elemental degrees of freedom. The elemental components help us to lower the QML model's learning curve, and, through interpolation across the periodic table, even enable "alchemical extrapolation" to covalent bonding between elements not part of training. This point is demonstrated for the prediction of covalent binding in single, double, and triple bonds among main-group elements as well as for atomization energies in organic molecules. We present numerical evidence that resulting QML energy models, after training on a few thousand random training instances, reach chemical accuracy for out-of-sample compounds. Compound datasets studied include thousands of structurally and compositionally diverse organic molecules, non-covalently bonded protein side-chains, (H2O)40-clusters, and crystalline solids. Learning curves for QML models also indicate competitive predictive power for various other electronic ground state properties of organic molecules, calculated with hybrid density functional theory, including polarizability, heat-capacity, HOMO-LUMO eigenvalues and gap, zero point vibrational energy, dipole moment, and highest vibrational fundamental frequency.
Collapse
Affiliation(s)
- Felix A Faber
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials, Department of Chemistry, University of Basel, Basel, Switzerland
| | - Anders S Christensen
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials, Department of Chemistry, University of Basel, Basel, Switzerland
| | - Bing Huang
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials, Department of Chemistry, University of Basel, Basel, Switzerland
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials, Department of Chemistry, University of Basel, Basel, Switzerland
| |
Collapse
|
9
|
Dutta NN, Patkowski K. Improving “Silver-Standard” Benchmark Interaction Energies with Bond Functions. J Chem Theory Comput 2018; 14:3053-3070. [DOI: 10.1021/acs.jctc.8b00204] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Narendra Nath Dutta
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
10
|
Burns LA, Faver JC, Zheng Z, Marshall MS, Smith DGA, Vanommeslaeghe K, MacKerell AD, Merz KM, Sherrill CD. The BioFragment Database (BFDb): An open-data platform for computational chemistry analysis of noncovalent interactions. J Chem Phys 2017; 147:161727. [PMID: 29096505 PMCID: PMC5656042 DOI: 10.1063/1.5001028] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2017] [Accepted: 08/20/2017] [Indexed: 11/14/2022] Open
Abstract
Accurate potential energy models are necessary for reliable atomistic simulations of chemical phenomena. In the realm of biomolecular modeling, large systems like proteins comprise very many noncovalent interactions (NCIs) that can contribute to the protein's stability and structure. This work presents two high-quality chemical databases of common fragment interactions in biomolecular systems as extracted from high-resolution Protein DataBank crystal structures: 3380 sidechain-sidechain interactions and 100 backbone-backbone interactions that inaugurate the BioFragment Database (BFDb). Absolute interaction energies are generated with a computationally tractable explicitly correlated coupled cluster with perturbative triples [CCSD(T)-F12] "silver standard" (0.05 kcal/mol average error) for NCI that demands only a fraction of the cost of the conventional "gold standard," CCSD(T) at the complete basis set limit. By sampling extensively from biological environments, BFDb spans the natural diversity of protein NCI motifs and orientations. In addition to supplying a thorough assessment for lower scaling force-field (2), semi-empirical (3), density functional (244), and wavefunction (45) methods (comprising >1M interaction energies), BFDb provides interactive tools for running and manipulating the resulting large datasets and offers a valuable resource for potential energy model development and validation.
Collapse
Affiliation(s)
- Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - John C Faver
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - Zheng Zheng
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - Michael S Marshall
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel G A Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Kenno Vanommeslaeghe
- Department of Analytical Chemistry and Pharmaceutical Technology (FABI), Center for Pharmaceutical Research (CePhaR), Vrije Universiteit Brussel (VUB), Laarbeeklaan 103, B-1090 Brussels, Belgium
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| | - Kenneth M Merz
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
11
|
Manna D, Kesharwani MK, Sylvetsky N, Martin JML. Conventional and Explicitly Correlated ab Initio Benchmark Study on Water Clusters: Revision of the BEGDB and WATER27 Data Sets. J Chem Theory Comput 2017; 13:3136-3152. [DOI: 10.1021/acs.jctc.6b01046] [Citation(s) in RCA: 69] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Debashree Manna
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Manoj K. Kesharwani
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| |
Collapse
|
12
|
Sirianni DA, Burns LA, Sherrill CD. Comparison of Explicitly Correlated Methods for Computing High-Accuracy Benchmark Energies for Noncovalent Interactions. J Chem Theory Comput 2016; 13:86-99. [DOI: 10.1021/acs.jctc.6b00797] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dominic A. Sirianni
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Lori A. Burns
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C. David Sherrill
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
13
|
Silva NJ, Machado FBC, Lischka H, Aquino AJA. π–π stacking between polyaromatic hydrocarbon sheets beyond dispersion interactions. Phys Chem Chem Phys 2016; 18:22300-10. [DOI: 10.1039/c6cp03749f] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
A slipped parallel structure of a stacked graphene flake showing a biconcave curvature.
Collapse
Affiliation(s)
- Nadeesha J. Silva
- Department of Chemistry and Biochemistry
- Texas Tech University
- Lubbock Texas 79409
- USA
| | | | - Hans Lischka
- Department of Chemistry and Biochemistry
- Texas Tech University
- Lubbock Texas 79409
- USA
- Institute for Theoretical Chemistry
| | - Adelia J. A. Aquino
- Department of Chemistry and Biochemistry
- Texas Tech University
- Lubbock Texas 79409
- USA
- School of Pharmaceutical Sciences and Technology
| |
Collapse
|
14
|
Parker TM, Sherrill CD. Assessment of Empirical Models versus High-Accuracy Ab Initio Methods for Nucleobase Stacking: Evaluating the Importance of Charge Penetration. J Chem Theory Comput 2015; 11:4197-204. [DOI: 10.1021/acs.jctc.5b00588] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Trent M. Parker
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C. David Sherrill
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
15
|
Burns LA, Marshall MS, Sherrill CD. Appointing silver and bronze standards for noncovalent interactions: a comparison of spin-component-scaled (SCS), explicitly correlated (F12), and specialized wavefunction approaches. J Chem Phys 2015; 141:234111. [PMID: 25527923 DOI: 10.1063/1.4903765] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
A systematic examination of noncovalent interactions as modeled by wavefunction theory is presented in comparison to gold-standard quality benchmarks available for 345 interaction energies of 49 bimolecular complexes. Quantum chemical techniques examined include spin-component-scaling (SCS) variations on second-order perturbation theory (MP2) [SCS, SCS(N), SCS(MI)] and coupled cluster singles and doubles (CCSD) [SCS, SCS(MI)]; also, method combinations designed to improve dispersion contacts [DW-MP2, MP2C, MP2.5, DW-CCSD(T)-F12]; where available, explicitly correlated (F12) counterparts are also considered. Dunning basis sets augmented by diffuse functions are employed for all accessible ζ-levels; truncations of the diffuse space are also considered. After examination of both accuracy and performance for 394 model chemistries, SCS(MI)-MP2/cc-pVQZ can be recommended for general use, having good accuracy at low cost and no ill-effects such as imbalance between hydrogen-bonding and dispersion-dominated systems or non-parallelity across dissociation curves. Moreover, when benchmarking accuracy is desirable but gold-standard computations are unaffordable, this work recommends silver-standard [DW-CCSD(T**)-F12/aug-cc-pVDZ] and bronze-standard [MP2C-F12/aug-cc-pVDZ] model chemistries, which support accuracies of 0.05 and 0.16 kcal/mol and efficiencies of 97.3 and 5.5 h for adenine·thymine, respectively. Choice comparisons of wavefunction results with the best symmetry-adapted perturbation theory [T. M. Parker, L. A. Burns, R. M. Parrish, A. G. Ryno, and C. D. Sherrill, J. Chem. Phys. 140, 094106 (2014)] and density functional theory [L. A. Burns, Á. Vázquez-Mayagoitia, B. G. Sumpter, and C. D. Sherrill, J. Chem. Phys. 134, 084107 (2011)] methods previously studied for these databases are provided for readers' guidance.
Collapse
Affiliation(s)
- Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Michael S Marshall
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
16
|
Waldrop JM, Song B, Patkowski K, Wang X. Accurate ab initio potential for the krypton dimer and transport properties of the low-density krypton gas. J Chem Phys 2015; 142:204307. [DOI: 10.1063/1.4921623] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan M. Waldrop
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| | - Bo Song
- Key Laboratory of Thermo-Fluid Science and Engineering, Ministry of Education, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| | - Xiaopo Wang
- Key Laboratory of Thermo-Fluid Science and Engineering, Ministry of Education, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China
| |
Collapse
|
17
|
Temelso B, Renner CR, Shields GC. Importance and Reliability of Small Basis Set CCSD(T) Corrections to MP2 Binding and Relative Energies of Water Clusters. J Chem Theory Comput 2015; 11:1439-48. [DOI: 10.1021/ct500944v] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Berhane Temelso
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Carla R. Renner
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - George C. Shields
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| |
Collapse
|
18
|
Kats D, Kreplin D, Werner HJ, Manby FR. Accurate thermochemistry from explicitly correlated distinguishable cluster approximation. J Chem Phys 2015; 142:064111. [DOI: 10.1063/1.4907591] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
- Daniel Kats
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - David Kreplin
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS,United Kingdom
| |
Collapse
|
19
|
Nascimento DR, DePrince AE. A Parametrized Coupled-Pair Functional for Molecular Interactions: PCPF-MI. J Chem Theory Comput 2014; 10:4324-31. [PMID: 26588129 DOI: 10.1021/ct500462p] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a parametrized coupled-pair functional, optimized to describe molecular interactions (PCPF-MI). The method has the same computational cost as singles and doubles configuration interaction, is size extensive, and yields energies that are stationary with respect to variations in all of its excitation coefficients. This last property facilitates the construction of density matrices and the evaluation of one- and two-electron properties. For the S22, HSG, S66, and A24 databases of van der Waals dimers, PCPF-MI computations yield interaction energies with mean unsigned errors of 0.326, 0.149, 0.214, and 0.044 kcal mol(-1), respectively, relative to benchmark computations. PCPF-MI interaction energies consistently improve upon those obtained from several other coupled-pair methods, including the averaged coupled-pair functional (ACPF), the coupled electron-pair approximation [CEPA(n), n = 0,1,3], and averaged quadratic coupled cluster (AQCC). Optimal parameters for spin-component-scaled (SCS) variants of each of these methods are also presented. SCS-CEPA(0), SCS-CEPA(1), SCS-ACPF, and SCS-PCPF-MI all perform similarly, with average errors for the four databases of roughly 0.1 kcal mol(-1). The PCPF-MI method, without additional SCS parametrization, has an average error of 0.192 kcal mol(-1) over the four databases and, when compared to the SCS-parametrized coupled-pair methods, has the added benefit that the energy is stationary with respect to variations in all its excitation coefficients.
Collapse
Affiliation(s)
- Daniel R Nascimento
- Department of Chemistry and Biochemistry, Florida State University , Tallahassee, Florida 32306-4390, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University , Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
20
|
Ahnen S, Hehn AS, Vogiatzis KD, Trachsel MA, Leutwyler S, Klopper W. Accurate computations of the structures and binding energies of the imidazole ⋯benzene and pyrrole ⋯benzene complexes. Chem Phys 2014. [DOI: 10.1016/j.chemphys.2014.05.023] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
21
|
Kalescky R, Kraka E, Cremer D. Accurate determination of the binding energy of the formic acid dimer: The importance of geometry relaxation. J Chem Phys 2014; 140:084315. [DOI: 10.1063/1.4866696] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Huang Y, Shao Y, Beran GJO. Accelerating MP2C dispersion corrections for dimers and molecular crystals. J Chem Phys 2014; 138:224112. [PMID: 23781788 DOI: 10.1063/1.4809981] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
The MP2C dispersion correction of Pitonak and Hesselmann [J. Chem. Theory Comput. 6, 168 (2010)] substantially improves the performance of second-order Møller-Plesset perturbation theory for non-covalent interactions, albeit with non-trivial computational cost. Here, the MP2C correction is computed in a monomer-centered basis instead of a dimer-centered one. When applied to a single dimer MP2 calculation, this change accelerates the MP2C dispersion correction several-fold while introducing only trivial new errors. More significantly, in the context of fragment-based molecular crystal studies, combination of the new monomer basis algorithm and the periodic symmetry of the crystal reduces the cost of computing the dispersion correction by two orders of magnitude. This speed-up reduces the MP2C dispersion correction calculation from a significant computational expense to a negligible one in crystals like aspirin or oxalyl dihydrazide, without compromising accuracy.
Collapse
Affiliation(s)
- Yuanhang Huang
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | | | | |
Collapse
|
23
|
Melicherčík M, Pitoňák M, Kellö V, Hobza P, Neogrády P. Off-Center Gaussian Functions, an Alternative Atomic Orbital Basis Set for Accurate Noncovalent Interaction Calculations of Large Systems. J Chem Theory Comput 2013; 9:5296-304. [DOI: 10.1021/ct400692b] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Miroslav Melicherčík
- Department
of Computer Science, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, 974 01 Banská Bystrica, Slovakia
| | - Michal Pitoňák
- Department
of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská Dolina, 842 15 Bratislava, Slovakia
- Computing
Center of the Slovak Academy of Sciences, Dúbravská cesta č. 9, 845 35 Bratislava, Slovakia
| | - Vladimír Kellö
- Department
of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská Dolina, 842 15 Bratislava, Slovakia
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy
of Sciences of the Czech Republic, v. v. i., Flemingovo nám. 2, 166 10 Praha 6, Czech Republic
- Department
of Physical Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| | - Pavel Neogrády
- Department
of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská Dolina, 842 15 Bratislava, Slovakia
| |
Collapse
|
24
|
Howard JC, Tschumper GS. Wavefunction methods for the accurate characterization of water clusters. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1168] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry University of Mississippi, University Mississippi USA
| |
Collapse
|
25
|
Vogiatzis KD, Klopper W. Accurate non-covalent interactions with basis-set corrections from interference-corrected perturbation theory: comparison with the S22B database. Mol Phys 2013. [DOI: 10.1080/00268976.2013.805888] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Konstantinos D. Vogiatzis
- Center for Functional Nanostructures (CFN) and Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), , P. O. Box 6980, Karlsruhe, D-76049, Germany
| | - Wim Klopper
- Center for Functional Nanostructures (CFN) and Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), , P. O. Box 6980, Karlsruhe, D-76049, Germany
| |
Collapse
|
26
|
Kalescky R, Kraka E, Cremer D. Local vibrational modes of the formic acid dimer – the strength of the double hydrogen bond. Mol Phys 2013. [DOI: 10.1080/00268976.2013.796070] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- R. Kalescky
- a Department of Chemistry , Southern Methodist University , 3215 Daniel Avenue, Dallas , Texas , 75275 , USA
| | - E. Kraka
- a Department of Chemistry , Southern Methodist University , 3215 Daniel Avenue, Dallas , Texas , 75275 , USA
| | - D. Cremer
- a Department of Chemistry , Southern Methodist University , 3215 Daniel Avenue, Dallas , Texas , 75275 , USA
| |
Collapse
|
27
|
Patkowski K. Basis set converged weak interaction energies from conventional and explicitly correlated coupled-cluster approach. J Chem Phys 2013; 138:154101. [DOI: 10.1063/1.4800981] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Cho Y, Min SK, Yun J, Kim WY, Tkatchenko A, Kim KS. Noncovalent Interactions of DNA Bases with Naphthalene and Graphene. J Chem Theory Comput 2013; 9:2090-6. [PMID: 26583555 DOI: 10.1021/ct301097u] [Citation(s) in RCA: 67] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The complexes of a DNA base bound to graphitic systems are studied. Considering naphthalene as the simplest graphitic system, DNA base-naphthalene complexes are scrutinized at high levels of ab initio theory including coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] at the complete basis set (CBS) limit. The stacked configurations are the most stable, where the CCSD(T)/CBS binding energies of guanine, adenine, thymine, and cytosine are 9.31, 8.48, 8.53, 7.30 kcal/mol, respectively. The energy components are investigated using symmetry-adapted perturbation theory based on density functional theory including the dispersion energy. We compared the CCSD(T)/CBS results with several density functional methods applicable to periodic systems. Considering accuracy and availability, the optB86b nonlocal functional and the Tkatchenko-Scheffler functional are used to study the binding energies of nucleobases on graphene. The predicted values are 18-24 kcal/mol, though many-body effects on screening and energy need to be further considered.
Collapse
Affiliation(s)
- Yeonchoo Cho
- Center for Superfunctional Materials, Department of Chemistry, Pohang University of Science and Technology , San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| | - Seung Kyu Min
- Center for Superfunctional Materials, Department of Chemistry, Pohang University of Science and Technology , San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| | - Jeonghun Yun
- Center for Superfunctional Materials, Department of Chemistry, Pohang University of Science and Technology , San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| | - Woo Youn Kim
- Department of Chemistry and KAIST Institute for NanoCentury, KAIST , Daejeon 305-701, Korea
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195, Berlin, Germany
| | - Kwang S Kim
- Center for Superfunctional Materials, Department of Chemistry, Pohang University of Science and Technology , San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| |
Collapse
|
29
|
Hollman DS, Wilke JJ, Schaefer HF. Explicitly correlated atomic orbital basis second order Møller–Plesset theory. J Chem Phys 2013; 138:064107. [DOI: 10.1063/1.4790582] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
DePrince AE, Sherrill CD. Accurate Noncovalent Interaction Energies Using Truncated Basis Sets Based on Frozen Natural Orbitals. J Chem Theory Comput 2012; 9:293-9. [PMID: 26589031 DOI: 10.1021/ct300780u] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We assess the accuracy of basis set truncations based on natural orbitals determined by second-order perturbation theory for computing noncovalent interaction energies with coupled cluster through perturbative triples [CCSD(T)]. We consider two methods for truncation: (i) the usual frozen natural orbital approach (FNO) in which the basis set truncation occurs before the iterative CCSD computation [FNO CCSD(T)] and (ii) an approach in which the truncation occurs only for the perturbative triples contribution [CCSD+FNO(T)]. The errors incurred are comparable for both methods and are small enough for the methods to be used for benchmark-quality studies of noncovalent interactions. For the FNO CCSD(T) method with a modest natural orbital occupation tolerance of 10(-5), the mean absolute error in the interaction energies for the S22 data set in an aug-cc-pVDZ basis set is only 0.012 kcal mol(-1) versus canonical CCSD(T) values. The same method exhibits a mean absolute error of 0.020 kcal mol(-1) for the S11 data set in the aug-cc-pVTZ basis set versus canonical CCSD(T) values.
Collapse
Affiliation(s)
- A Eugene DePrince
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
31
|
Platts JA, Hill JG, Riley KE, Řezáč J, Hobza P. Basis Set Dependence of Interaction Energies Computed Using Composite Post-MP2 Methods. J Chem Theory Comput 2012; 9:330-7. [DOI: 10.1021/ct300842d] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- James A. Platts
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, United
Kingdom
| | - J. Grant Hill
- School of Chemistry, University of Glasgow, Joseph Black Building, University
Ave, Glasgow, G12 8QQ, United Kingdom
| | - Kevin E. Riley
- Department of Chemistry, Xavier University of Louisiana, 1 Drexel Drive, Box
22, New Orleans, Louisiana 70125, United States
| | - Jan Řezáč
- Institute of Organic Chemistry
and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry
and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, 166 10 Prague, Czech Republic
| |
Collapse
|
32
|
Smith DGA, Patkowski K. Interactions between Methane and Polycyclic Aromatic Hydrocarbons: A High Accuracy Benchmark Study. J Chem Theory Comput 2012; 9:370-89. [DOI: 10.1021/ct3008809] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Daniel G. A. Smith
- Department of Chemistry and Biochemistry,
Auburn University,
Auburn, Alabama 36849, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry,
Auburn University,
Auburn, Alabama 36849, United States
| |
Collapse
|