1
|
Kříž K, Schmidt L, Andersson AT, Walz MM, van der Spoel D. An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation. J Chem Inf Model 2023; 63:412-431. [PMID: 36630710 PMCID: PMC9875315 DOI: 10.1021/acs.jcim.2c01127] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2022] [Indexed: 01/12/2023]
Abstract
Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.
Collapse
Affiliation(s)
- Kristian Kříž
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - Lisa Schmidt
- Faculty
of Biosciences, University of Heidelberg, Heidelberg69117, Germany
| | - Alfred T. Andersson
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - Marie-Madeleine Walz
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - David van der Spoel
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| |
Collapse
|
2
|
Huang HH, Wang YS, Chao SD. A Minimum Quantum Chemistry CCSD(T)/CBS Data Set of Dimeric Interaction Energies for Small Organic Functional Groups: Heterodimers. ACS OMEGA 2022; 7:20059-20080. [PMID: 35722020 PMCID: PMC9201891 DOI: 10.1021/acsomega.2c01888] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 05/04/2022] [Indexed: 06/15/2023]
Abstract
We extend our previous quantum chemistry calculations of interaction energies for 31 homodimers of small organic functional groups (the SOFG-31 data set) by including 239 heterodimers with monomers selected within the SOFG-31 data set, thus resulting in the SOFG-31+239 data set. The minimum-level theoretical scheme contains (1) the basis set superposition error corrected supermolecule (BSSE-SM) approach for intermolecular interactions; (2) the second-order Møller-Plesset perturbation theory (MP2) with the Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets for the geometry optimization and correlation energy calculations; and (3) the single-point energy calculations with the coupled cluster with single, double, and perturbative triple excitations method at the complete basis set limit [CCSD(T)/CBS] using the well-tested extrapolation methods for the MP2 energy calibrations. In addition, we have performed a parallel series of energy decomposition calculations based on the symmetry adapted perturbation theory (SAPT) in order to gain chemical insights. That the above procedure cannot be further reduced has been proven to be very crucial for constructing reliable data sets of interaction energies. The calculated CCSD(T)/CBS interaction energy data can serve as a benchmark for testing or training less accurate but more efficient calculation methods, such as the electronic density functional theory. As an application, we employ a segmental SAPT model previously developed for the SOFG-31 data set to predict binding energies of large heterodimer complexes. These model energy "quanta" can be used in coarse-grained molecular dynamics simulations by avoiding large-scale calculations.
Collapse
Affiliation(s)
- Hsing-Hsiang Huang
- Institute
of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan R.O.C.
| | - Yi-Siang Wang
- School
of Chemistry & Biochemistry, Georgia
Institute of Technology, Atlanta, Georgia 30332, United States
| | - Sheng D. Chao
- Institute
of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan R.O.C.
| |
Collapse
|
3
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2913-2930. [PMID: 35412817 DOI: 10.1021/acs.jctc.2c00036] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, MALTA Consolider Team, Oviedo E-33006, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
4
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2208-2232. [PMID: 35313106 DOI: 10.1021/acs.jctc.1c01128] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- MALTA Consolider Team, Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, E-33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
5
|
Krstić M, Fink K, Sharapa DI. The Adsorption of Small Molecules on the Copper Paddle-Wheel: Influence of the Multi-Reference Ground State. Molecules 2022; 27:912. [PMID: 35164179 PMCID: PMC8840508 DOI: 10.3390/molecules27030912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Revised: 01/26/2022] [Accepted: 01/26/2022] [Indexed: 11/16/2022] Open
Abstract
We report a theoretical study of the adsorption of a set of small molecules (C2H2, CO, CO2, O2, H2O, CH3OH, C2H5OH) on the metal centers of the "copper paddle-wheel"-a key structural motif of many MOFs. A systematic comparison between DFT of different rungs, single-reference post-HF methods (MP2, SOS-MP2, MP3, DLPNO-CCSD(T)), and multi-reference approaches (CASSCF, DCD-CAS(2), NEVPT2) is performed in order to find a methodology that correctly describes the complicated electronic structure of paddle-wheel structure together with a reasonable description of non-covalent interactions. Apart from comparison with literature data (experimental values wherever possible), benchmark calculations with DLPNO-MR-CCSD were also performed. Despite tested methods show qualitative agreement in the majority of cases, we showed and discussed reasons for quantitative differences as well as more fundamental problems of specific cases.
Collapse
Affiliation(s)
- Marjan Krstić
- Institute for Theoretical Solid State Physics (TFP), Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany;
| | - Karin Fink
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany;
| | - Dmitry I. Sharapa
- Institute of Catalysis Research and Technology (IKFT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| |
Collapse
|
6
|
Loipersberger M, Bertels LW, Lee J, Head-Gordon M. Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals. J Chem Theory Comput 2021; 17:5582-5599. [PMID: 34382394 PMCID: PMC9948597 DOI: 10.1021/acs.jctc.1c00469] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N6) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.
Collapse
Affiliation(s)
| | - Luke W. Bertels
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Virginia Tech, Blacksburg, VA 24061, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Columbia University, NY
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
7
|
Towne J, Liang X, Kontogeorgis GM. Application of Quantum Chemistry Insights to the Prediction of Phase Equilibria in Associating Systems. Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.1c00072] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- John Towne
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, Søltofts Plads, Building 229, 2800 Kongens Lyngby, Denmark
| | - Xiaodong Liang
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, Søltofts Plads, Building 229, 2800 Kongens Lyngby, Denmark
| | - Georgios M. Kontogeorgis
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, Søltofts Plads, Building 229, 2800 Kongens Lyngby, Denmark
| |
Collapse
|
8
|
Chang YM, Wang YS, Chao SD. A minimum quantum chemistry CCSD(T)/CBS dataset of dimeric interaction energies for small organic functional groups. J Chem Phys 2020; 153:154301. [PMID: 33092384 DOI: 10.1063/5.0019392] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We have performed a quantum chemistry study on the bonding patterns and interaction energies for 31 dimers of small organic functional groups (dubbed the SOFG-31 dataset), including the alkane-alkene-alkyne (6 + 4 + 4 = 14, AAA) groups, alcohol-aldehyde-ketone (4 + 4 + 3 = 11, AAK) groups, and carboxylic acid-amide (3 + 3 = 6, CAA) groups. The basis set superposition error corrected super-molecule approach using the second order Møller-Plesset perturbation theory (MP2) with the Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets has been employed in the geometry optimization and energy calculations. To calibrate the MP2 calculated interaction energies for these dimeric complexes, we perform single-point calculations with the coupled cluster with single, double, and perturbative triple excitations method at the complete basis set limit [CCSD(T)/CBS] using the well-tested extrapolation methods. In order to gain more physical insights, we also perform a parallel series of energy decomposition calculations based on the symmetry adapted perturbation theory (SAPT). The collection of these CCSD(T)/CBS interaction energy values can serve as a minimum quantum chemistry dataset for testing or training less accurate but more efficient calculation methods. As an application, we further propose a segmental SAPT model based on chemically recognizable segments in a specific functional group. These model interactions can be used to construct coarse-grained force fields for larger molecular systems.
Collapse
Affiliation(s)
- Yu-Ming Chang
- Institute of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan
| | - Yi-Siang Wang
- Institute of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan
| | - Sheng D Chao
- Institute of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan
| |
Collapse
|
9
|
Gottschalk HC, Poblotzki A, Fatima M, Obenchain DA, Pérez C, Antony J, Auer AA, Baptista L, Benoit DM, Bistoni G, Bohle F, Dahmani R, Firaha D, Grimme S, Hansen A, Harding ME, Hochlaf M, Holzer C, Jansen G, Klopper W, Kopp WA, Krasowska M, Kröger LC, Leonhard K, Mogren Al-Mogren M, Mouhib H, Neese F, Pereira MN, Prakash M, Ulusoy IS, Mata RA, Suhm MA, Schnell M. The first microsolvation step for furans: New experiments and benchmarking strategies. J Chem Phys 2020; 152:164303. [DOI: 10.1063/5.0004465] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Hannes C. Gottschalk
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Anja Poblotzki
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Mariyam Fatima
- Deutsches Elektronen-Synchrotron, Notkestr. 85, 22607 Hamburg, Germany
| | | | - Cristóbal Pérez
- Deutsches Elektronen-Synchrotron, Notkestr. 85, 22607 Hamburg, Germany
| | - Jens Antony
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Alexander A. Auer
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Leonardo Baptista
- Departamento de Química e Ambiental, Universidade do Estado do Rio de Janeiro, Faculdade de Tecnologia, Resende, RJ, Brazil
| | - David M. Benoit
- Department of Physics and Mathematics, E. A. Milne Centre for Astrophysics and G. W. Gray Centre for Advanced Materials Chemistry, University of Hull, Hull HU6 7RX, United Kingdom
| | - Giovanni Bistoni
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Fabian Bohle
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Rahma Dahmani
- Université Gustave Eiffel, COSYS/LISIS, 5 Blvd. Descartes, 77454 Marne-La-Vallée, France
| | - Dzmitry Firaha
- Lehrstuhl für Technische Thermodynamik, RWTH Aachen University, 52062 Aachen, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Michael E. Harding
- Institute of Nanotechnology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Majdi Hochlaf
- Université Gustave Eiffel, COSYS/LISIS, 5 Blvd. Descartes, 77454 Marne-La-Vallée, France
| | - Christof Holzer
- Theoretical Chemistry Group, Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Georg Jansen
- Fakultät für Chemie, Universität Duisburg-Essen, Universitätsstr. 5, 45117 Essen, Germany
| | - Wim Klopper
- Theoretical Chemistry Group, Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Wassja A. Kopp
- Lehrstuhl für Technische Thermodynamik, RWTH Aachen University, 52062 Aachen, Germany
| | - Małgorzata Krasowska
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Leif C. Kröger
- Lehrstuhl für Technische Thermodynamik, RWTH Aachen University, 52062 Aachen, Germany
| | - Kai Leonhard
- Lehrstuhl für Technische Thermodynamik, RWTH Aachen University, 52062 Aachen, Germany
| | - Muneerah Mogren Al-Mogren
- Chemistry Department, Faculty of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Kingdom of Saudi Arabia
| | - Halima Mouhib
- Université Gustave Eiffel, COSYS/LISIS, 5 Blvd. Descartes, 77454 Marne-La-Vallée, France
| | - Frank Neese
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Max N. Pereira
- Departamento de Química e Ambiental, Universidade do Estado do Rio de Janeiro, Faculdade de Tecnologia, Resende, RJ, Brazil
| | - Muthuramalingam Prakash
- Department of Chemistry, SRM Institute of Science and Technology, Kattankulathur 603 203, Tamil Nadu, India
| | - Inga S. Ulusoy
- Theoretical Chemistry, Institute of Physical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Ricardo A. Mata
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Martin A. Suhm
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Melanie Schnell
- Deutsches Elektronen-Synchrotron, Notkestr. 85, 22607 Hamburg, Germany
| |
Collapse
|
10
|
Sethio D, Martins JBL, Lawson Daku LM, Hagemann H, Kraka E. Modified Density Functional Dispersion Correction for Inorganic Layered MFX Compounds (M = Ca, Sr, Ba, Pb and X = Cl, Br, I). J Phys Chem A 2020; 124:1619-1633. [PMID: 31999454 DOI: 10.1021/acs.jpca.9b10357] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
MFX (M = Ca, Ba, Sr, Pb and X = Cl, Br, I) compounds have received considerable attention due to their technological application as X-ray detectors, pressure sensors, and optical data storage materials, when doped with rare-earth ions. MFX compounds belong to the class of layered materials with a tetragonal Matlockite crystal structure, characterized by weakly stacked double-halide layers along the crystallographic c-axis. These layers predominantly determine phase transitions, elastic, and mechanical properties. However, the correct description of the lattice parameter c is a challenge for most standard DFT functionals, which tend to overestimate the lattice parameter c. Because of the weak interactions between the halide layers, dispersion-corrected functionals seem to be a better choice. We investigated 11 different inorganic layered MFX compounds for which experimental data are available, with standard and dispersion-corrected functionals to assess their performance in reproducing the lattice parameter c, structural, and vibrational properties of the MFX compounds. Our results revealed that these functionals do not describe the weak interactions between the halide layers in a balanced way. Therefore, we modified Grimme's popular DFT-D2 dispersion correction scheme in two different ways by (i) replacing the dispersion coefficients and van der Waals radii with those of noble gas atoms or (ii) increasing the van der Waals radii of the MFX atoms up to 40%. Comparison with the available experimental data revealed that the latter approach applied to the PBE (Perdew-Burke-Ernzerhof)-D2 functional with 30% increased van der Waals radii, which we coined PBE-D2* (Srvdw 1.30) is best suited to fine-tune the description of the weak interlayer interactions in MFX compounds, thus significantly improving the description of their structural, vibrational, and mechanical properties. Work is in progress applying this new, computationally inexpensive scheme to other inorganic layered compounds and periodic systems with weakly stacked layers.
Collapse
Affiliation(s)
- Daniel Sethio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry , Southern Methodist University , 3215 Daniel Avenue , Dallas , Texas 75275-0314 , United States
| | - João B L Martins
- Institute of Chemistry , University of Brasilia , Brasilia , DF 70910-900 , Brazil
| | - Latévi Max Lawson Daku
- Department of Physical Chemistry , University of Geneva , 30 Quai Ernest-Ansermet , CH-1211 Geneva 4 , Switzerland
| | - Hans Hagemann
- Department of Physical Chemistry , University of Geneva , 30 Quai Ernest-Ansermet , CH-1211 Geneva 4 , Switzerland
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry , Southern Methodist University , 3215 Daniel Avenue , Dallas , Texas 75275-0314 , United States
| |
Collapse
|
11
|
Shaw RA, Hill JG. A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach. J Chem Theory Comput 2019; 15:5352-5369. [PMID: 31465215 DOI: 10.1021/acs.jctc.9b00615] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A novel method for the accurate and efficient calculation of interaction energies in weakly bound complexes composed of a large number of molecules is presented. The new ALMO+RPAd method circumvents the prohibitive scaling of coupled cluster singles and doubles while still providing similar accuracy across a diverse range of weakly bound chemical systems. Linear-scaling procedures for the Fock build are given utilizing absolutely localized molecular orbitals (ALMOs), resulting in the a priori exclusion of basis set superposition errors. A bespoke data structure and algorithm using density fitting are described, leading to linear scaling for the storage and computation of the two-electron integrals. Electron correlation is included through a new, linear-scaling pairwise local random phase approximation approach, including exchange interactions, and decomposed into purely dispersive excitations (RPAxd). Collectively, these allow meaningful decomposition of the interaction energy into physically distinct contributions: electrostatic, polarization, charge transfer, and dispersion. Comparison with symmetry-adapted perturbation theory shows good qualitative agreement. Tests on various dimers and the S66 benchmark set demonstrate results within 0.5 kcal mol-1 of coupled cluster singles and doubles results. On a large cluster of water molecules, we achieve calculations involving over 3500 orbital and 12,000 auxiliary basis functions in under 10 min on a single CPU core.
Collapse
Affiliation(s)
- Robert A Shaw
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| | - J Grant Hill
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| |
Collapse
|
12
|
Morgante P, Peverati R. ACCDB: A collection of chemistry databases for broad computational purposes. J Comput Chem 2018; 40:839-848. [DOI: 10.1002/jcc.25761] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Revised: 10/09/2018] [Accepted: 11/11/2018] [Indexed: 01/04/2023]
Affiliation(s)
- Pierpaolo Morgante
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| | - Roberto Peverati
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| |
Collapse
|
13
|
Witte J, Neaton JB, Head-Gordon M. Push it to the limit: comparing periodic and local approaches to density functional theory for intermolecular interactions. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1542164] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Jonathon Witte
- Department of Chemistry, University of California, Berkeley, California, United States
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California, United States
| | - Jeffrey B. Neaton
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California, United States
- Department of Physics, University of California, Berkeley, California, United States
- Kavli Energy Nanosciences Institute at Berkeley, Berkeley, California, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States
| |
Collapse
|
14
|
Hänninen V, Murdachaew G, Nathanson GM, Gerber RB, Halonen L. Ab initio molecular dynamics studies of formic acid dimer colliding with liquid water. Phys Chem Chem Phys 2018; 20:23717-23725. [PMID: 30191926 DOI: 10.1039/c8cp03857k] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Ab initio molecular dynamics simulations of formic acid (FA) dimer colliding with liquid water at 300 K have been performed using density functional theory. The two energetically lowest FA dimer isomers were collided with a water slab at thermal and high kinetic energies up to 68kBT. Our simulations agree with recent experimental observations of nearly a complete uptake of gas-phase FA dimer: the calculated average kinetic energy of the dimers immediately after collision is 5 ± 4% of the incoming kinetic energy, which compares well with the experimental value of 10%. Simulations support the experimental observation of no delayed desorption of FA dimers following initial adsorption. Our analysis shows that the FA dimer forms hydrogen bonds with surface water molecules, where the hydrogen bond order depends on the dimer structure, such that the most stable isomer possesses fewer FA-water hydrogen bonds than the higher energy isomer. Nevertheless, even the most stable isomer can attach to the surface through one hydrogen bond despite its reduced hydrophilicity. Our simulations further show that the probability of FA dimer dissociation is increased by high collision energies, the dimer undergoes isomerization from the higher energy to the lowest energy isomer, and concerted double-proton transfer occurs between the FA monomers. Interestingly, proton transfer appears to be driven by the release of energy arising from such isomerization, which stimulates those internal vibrational degrees of freedom that overcome the barrier of a proton transfer.
Collapse
Affiliation(s)
- Vesa Hänninen
- Department of Chemistry, University of Helsinki, P. O. Box 55 (A.I. Virtasen aukio 1), FI-00014, Finland.
| | | | | | | | | |
Collapse
|
15
|
Mardirossian N, Head-Gordon M. Survival of the most transferable at the top of Jacob's ladder: Defining and testing the ωB97M(2) double hybrid density functional. J Chem Phys 2018; 148:241736. [PMID: 29960332 PMCID: PMC5991970 DOI: 10.1063/1.5025226] [Citation(s) in RCA: 117] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2018] [Accepted: 05/09/2018] [Indexed: 01/19/2023] Open
Abstract
A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resulting functional, ωB97M(2), is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 data points clearly surpasses that of all of the tested density functionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.
Collapse
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
16
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Atom-Centered Potentials with Dispersion-Corrected Minimal-Basis-Set Hartree–Fock: An Efficient and Accurate Computational Approach for Large Molecular Systems. J Chem Theory Comput 2018; 14:726-738. [DOI: 10.1021/acs.jctc.7b01158] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Viki Kumar Prasad
- Department
of Chemistry, University of British Columbia, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- Department
of Chemistry, University of British Columbia, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Gino A. DiLabio
- Department
of Chemistry, University of British Columbia, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
- Faculty
of Management, University of British Columbia, 1137 Alumni Avenue, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
17
|
Witte J, Neaton JB, Head-Gordon M. Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C. J Chem Phys 2017. [DOI: 10.1063/1.4986962] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
|
18
|
Mardirossian N, Head-Gordon M. Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1333644] [Citation(s) in RCA: 709] [Impact Index Per Article: 101.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
19
|
Červinka C, Fulem M. State-of-the-Art Calculations of Sublimation Enthalpies for Selected Molecular Crystals and Their Computational Uncertainty. J Chem Theory Comput 2017; 13:2840-2850. [PMID: 28437618 DOI: 10.1021/acs.jctc.7b00164] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A computational methodology for calculation of sublimation enthalpies of molecular crystals from first principles is developed and validated by comparison to critically evaluated literature experimental data. Temperature-dependent sublimation enthalpies for a set of selected 22 molecular crystals in their low-temperature phases are calculated. The computational methodology consists of several building blocks based on high-level electronic structure methods of quantum chemistry and statistical thermodynamics. Ab initio methods up to the coupled clusters with iterative treatment of single and double excitations and perturbative triples correction with an estimated complete basis set description [CCSD(T)/CBS] are used to calculate the cohesive energies of crystalline phases within a fragment-based additive scheme. Density functional theory (DFT) calculations with periodic boundary conditions (PBC) coupled with the quasi-harmonic approximation are used to evaluate the thermal contributions to the enthalpy of the solid phase. The properties of the vapor phase are calculated within the ideal-gas model using the rigid-rotor harmonic-oscillator model with correction for internal rotation using a one-dimensional hindered rotor approximation and a proper treatment of the molecular rotational degrees of freedom in the vicinity of 0 K. All individual terms contributing to the sublimation enthalpy as a function of temperature are discussed and their uncertainties estimated by comparison to critically evaluated experimental data.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology , Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Michal Fulem
- Department of Physical Chemistry, University of Chemistry and Technology , Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| |
Collapse
|
20
|
Witte J, Mardirossian N, Neaton JB, Head-Gordon M. Assessing DFT-D3 Damping Functions Across Widely Used Density Functionals: Can We Do Better? J Chem Theory Comput 2017; 13:2043-2052. [DOI: 10.1021/acs.jctc.7b00176] [Citation(s) in RCA: 56] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | | | - Jeffrey B. Neaton
- Kavli Energy
Nanosciences
Institute at Berkeley, Berkeley, California 94720, United States
| | | |
Collapse
|
21
|
Kokoschka M, Galgonek J, Vondrasek J, Hobza P. Computational methods for the description of pharmacologically relevant platinum complexes--molecular structure and bond dissociation. Phys Chem Chem Phys 2016; 18:4051-62. [PMID: 26777459 DOI: 10.1039/c5cp03391h] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cancer is after cardiovascular disease the most frequent cause of death in Europe. In 28 of 53 countries considered in this area it is already the leading cause of death and expected to gain even more importance until the year 2020. Amongst the large arsenal of different anti-cancer drugs, platinum drugs belong to the first developed anticancer drugs and still have a large impact on cancer therapy. Nevertheless therapy with platinum-anticancer drugs is accompanied by severe adverse effects caused by frequent interactions with the amino acids of different human proteins. Computational chemistry offers methods to study such interactions and even those of not yet synthesized drugs in silico. For such studies a profound knowledge of the prediction quality of various computational methods towards platinum-drug-like complexes is necessary. By this article we are aiming on delivering important accuracy information of the frequently used computational methods. Most important findings are the high performance of the double hybrid functional B2PLYP for the calculation of geometries, even in small basis sets, followed by BP86 and PBE and the still acceptable performance of the semi-empirical Method PM6-D3H4X for extremely large systems. To follow absolute energies of the dissociation process, LPNO-CEPA and B3LYP-D3 can be suggested while SCS-MP2 shows an extremely narrow standard deviation and a low maximum error, which make it an ideal candidate for relative energy calculations in the exploration of reaction mechanisms.
Collapse
Affiliation(s)
- Malte Kokoschka
- Gilead Sciences & IOCB Research Center, Department of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic. and Department of Physical Chemistry, Palacky University, 771 46 Olomouc, Czech Republic
| | - Jakub Galgonek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic
| | - Jiri Vondrasek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic
| | - Pavel Hobza
- Gilead Sciences & IOCB Research Center, Department of Computational Chemistry, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic. and Department of Physical Chemistry, Palacky University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
22
|
Mardirossian N, Head-Gordon M. How Accurate Are the Minnesota Density Functionals for Noncovalent Interactions, Isomerization Energies, Thermochemistry, and Barrier Heights Involving Molecules Composed of Main-Group Elements? J Chem Theory Comput 2016; 12:4303-25. [PMID: 27537680 DOI: 10.1021/acs.jctc.6b00637] [Citation(s) in RCA: 286] [Impact Index Per Article: 35.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
The 14 Minnesota density functionals published between the years 2005 and early 2016 are benchmarked on a comprehensive database of 4986 data points (84 data sets) involving molecules composed of main-group elements. The database includes noncovalent interactions, isomerization energies, thermochemistry, and barrier heights, as well as equilibrium bond lengths and equilibrium binding energies of noncovalent dimers. Additionally, the sensitivity of the Minnesota density functionals to the choice of basis set and integration grid is explored for both noncovalent interactions and thermochemistry. Overall, the main strength of the hybrid Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., M06-2X), barrier heights (e.g., M08-HX, M08-SO, MN15), and systems heavily characterized by self-interaction error (e.g., M06-2X, M08-HX, M08-SO, MN15), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-2X is recommended from the 10 hybrid Minnesota functionals). Similarly, the main strength of the local Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., MN15-L), barrier heights (e.g., MN12-L), and systems heavily characterized by self-interaction error (e.g., MN12-L and MN15-L), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-L is clearly the best from the four local Minnesota functionals). As an overall guide, M06-2X and MN15 are perhaps the most broadly useful hybrid Minnesota functionals, while M06-L and MN15-L are perhaps the most broadly useful local Minnesota functionals, although each has different strengths and weaknesses.
Collapse
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California , Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory , Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California , Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory , Berkeley, California 94720, United States
| |
Collapse
|
23
|
Červinka C, Fulem M, Růžička K. CCSD(T)/CBS fragment-based calculations of lattice energy of molecular crystals. J Chem Phys 2016; 144:064505. [PMID: 26874495 DOI: 10.1063/1.4941055] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
A comparative study of the lattice energy calculations for a data set of 25 molecular crystals is performed using an additive scheme based on the individual energies of up to four-body interactions calculated using the coupled clusters with iterative treatment of single and double excitations and perturbative triples correction (CCSD(T)) with an estimated complete basis set (CBS) description. The CCSD(T)/CBS values on lattice energies are used to estimate sublimation enthalpies which are compared with critically assessed and thermodynamically consistent experimental values. The average absolute percentage deviation of calculated sublimation enthalpies from experimental values amounts to 13% (corresponding to 4.8 kJ mol(-1) on absolute scale) with unbiased distribution of positive to negative deviations. As pair interaction energies present a dominant contribution to the lattice energy and CCSD(T)/CBS calculations still remain computationally costly, benchmark calculations of pair interaction energies defined by crystal parameters involving 17 levels of theory, including recently developed methods with local and explicit treatment of electronic correlation, such as LCC and LCC-F12, are also presented. Locally and explicitly correlated methods are found to be computationally effective and reliable methods enabling the application of fragment-based methods for larger systems.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Michal Fulem
- 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
| |
Collapse
|
24
|
Mardirossian N, Head-Gordon M. ωB97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation. J Chem Phys 2016; 144:214110. [DOI: 10.1063/1.4952647] [Citation(s) in RCA: 370] [Impact Index Per Article: 46.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
25
|
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
|
26
|
Řezáč J, Hobza P. Benchmark Calculations of Interaction Energies in Noncovalent Complexes and Their Applications. Chem Rev 2016; 116:5038-71. [DOI: 10.1021/acs.chemrev.5b00526] [Citation(s) in RCA: 281] [Impact Index Per Article: 35.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Řezáč
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
27
|
Møller KH, Hansen AS, Kjaergaard HG. Gas Phase Detection of the NH-P Hydrogen Bond and Importance of Secondary Interactions. J Phys Chem A 2015; 119:10988-98. [PMID: 26451467 DOI: 10.1021/acs.jpca.5b08358] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have observed the NH···P hydrogen bond in a gas phase complex. The bond is identified in the dimethylamine-trimethylphosphine complex by a red shift of the fundamental NH-stretching frequency observed using Fourier transform infrared spectroscopy (FT-IR). On the basis of the measured NH-stretching frequency red shifts, we find that P is a hydrogen bond acceptor atom similar in strength to S. Both are stronger acceptors than O and significantly weaker acceptors than N. The hydrogen bond angle, ∠NHP, is found to be very sensitive to the functional employed in density functional theory (DFT) optimizations of the complex and is a possible parameter to assess the quality of DFT functionals. Natural bonding orbital (NBO) energies and results from the topological methods atoms in molecules (AIM) and noncovalent interactions (NCI) indicate that the sensitivity is caused by the weakness of the hydrogen bond compared to secondary interactions. We find that B3LYP favors the hydrogen bond and M06-2X favors the secondary interactions leading to under- and overestimation, respectively, of the hydrogen bond angle relative to a DF-LCCSD(T)-F12a calculated angle. The remaining functionals tested, B3LYP-D3, B3LYP-D3BJ, CAM-B3LYP, and ωB97X-D, as well as MP2, show comparable contributions from the hydrogen bond and the secondary interactions and are close to DF-LCCSD(T)-F12a results.
Collapse
Affiliation(s)
- Kristian H Møller
- Department of Chemistry, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | - Anne S Hansen
- Department of Chemistry, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | - Henrik G Kjaergaard
- Department of Chemistry, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| |
Collapse
|
28
|
Theoretical description of 2D-cluster formation of nonionic surfactants at the air/water interface. Colloid Polym Sci 2015. [DOI: 10.1007/s00396-015-3630-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
29
|
Skelton AA, Khedkar VM, Fried JR. All-atom molecular dynamics simulations of an artificial sodium channel in a lipid bilayer: the effect of water solvation/desolvation of the sodium ion. J Biomol Struct Dyn 2015; 34:529-39. [PMID: 26046587 DOI: 10.1080/07391102.2015.1044473] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
All-atom molecular dynamics is used to investigate the transport of Na(+) across a 1,2-dioleoyl-sn-glycero-3-phosphocholine lipid bilayer facilitated by a diazacrown hydraphile. Specifically, the free energy of Na(+) passing through the bilayer is calculated using the adaptive biasing force method to study the free energy associated with the increase in Na(+) transport in the presence of the hydraphile molecule. The results show that water interaction greatly influences Na(+) transport through the lipid bilayer as water is pulled through the bilayer with Na(+) forming a water channel. The hydraphile causes a reduction in the free energy barrier for the transport of Na(+) through the head group part of the lipid bilayer since it complexes the Na(+) reducing the necessity for water to be complexed and, therefore, dragged through with Na(+), an energetically unfavorable process. The free energy associated with Na(+) being desolvated within the bilayer is significantly decreased in the presence of the hydraphile molecule; the hydraphile increases the number of solvation states of Na(+) that can be adopted, and this increase in the number of available configurations provides an entropic explanation for the success of the hydraphile.
Collapse
Affiliation(s)
- A A Skelton
- a School of Health Sciences , University of KwaZulu Natal , Westville Campus, Durban , South Africa
| | - V M Khedkar
- a School of Health Sciences , University of KwaZulu Natal , Westville Campus, Durban , South Africa
| | - J R Fried
- b Chemical Engineering , University of Louisville , Louisville , KY , USA
| |
Collapse
|
30
|
Wang C, Mo Y, Wagner JP, Schreiner PR, Jemmis ED, Danovich D, Shaik S. The Self-Association of Graphane Is Driven by London Dispersion and Enhanced Orbital Interactions. J Chem Theory Comput 2015; 11:1621-30. [DOI: 10.1021/acs.jctc.5b00075] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Changwei Wang
- Department
of Chemistry, School of Science, China University of Petroleum (East China), Changjianxi Road 66, 266580 Tsingtao, China
| | - Yirong Mo
- Department
of Chemistry, Western Michigan University, Kalamazoo, Michigan 49008, United States
| | - J. Philipp Wagner
- Institute
of Organic Chemistry, Justus-Liebig University, Heinrich-Buff-Ring 58, 35392 Giessen, Germany
| | - Peter R. Schreiner
- Institute
of Organic Chemistry, Justus-Liebig University, Heinrich-Buff-Ring 58, 35392 Giessen, Germany
| | - Eluvathingal D. Jemmis
- Department
of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012 India
| | - David Danovich
- Institute
of Chemistry and Lise Meitner Minerva Center for Computational Quantum
Chemistry, The Hebrew University, Jerusalem 91904, Israel
| | - Sason Shaik
- Institute
of Chemistry and Lise Meitner Minerva Center for Computational Quantum
Chemistry, The Hebrew University, Jerusalem 91904, Israel
| |
Collapse
|
31
|
Saleh G, Gatti C, Lo Presti L. Energetics of non-covalent interactions from electron and energy density distributions. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2014.10.011] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
32
|
Budzák Š, Mach P, Juhász G, Medved’ M, Kysel’ O. Theoretical study of charge transfer complexes between antithyroid thioamides and TCNE: Thermodynamics of the complex formation. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2014.10.035] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
33
|
Cunha WF, Gargano R, Garcia E, Politi JRS, Albernaz AF, Martins JBL. Rovibrational energy and spectroscopic constant calculations of CH 4 ⋯ CH 4, CH 4 ⋯ H 2 O, CH 4 ⋯ CHF 3, and H 2 O ⋯ CHF 3 dimers. J Mol Model 2014; 20:2298. [DOI: 10.1007/s00894-014-2298-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2014] [Accepted: 05/05/2014] [Indexed: 11/28/2022]
|
34
|
Theoretical study (CC2, DFT and PCM) of charge transfer complexes between antithyroid thioamides and TCNE: electronic CT transitions. J Mol Model 2014; 20:2312. [DOI: 10.1007/s00894-014-2312-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2014] [Accepted: 05/19/2014] [Indexed: 10/25/2022]
|
35
|
Huang Y, Goldey M, Head-Gordon M, Beran GJO. Achieving High-Accuracy Intermolecular Interactions by Combining Coulomb-Attenuated Second-Order Møller–Plesset Perturbation Theory with Coupled Kohn–Sham Dispersion. J Chem Theory Comput 2014; 10:2054-63. [DOI: 10.1021/ct5002329] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Yuanhang Huang
- Department
of Chemistry, University of California, Riverside, California 92521, United States
| | - Matthew Goldey
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Gregory J. O. Beran
- Department
of Chemistry, University of California, Riverside, California 92521, United States
| |
Collapse
|
36
|
Alonso M, Woller T, Martín-Martínez FJ, Contreras-García J, Geerlings P, De Proft F. Understanding the fundamental role of π/π, σ/σ, and σ/π dispersion interactions in shaping carbon-based materials. Chemistry 2014; 20:4931-41. [PMID: 24692007 DOI: 10.1002/chem.201400107] [Citation(s) in RCA: 86] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2014] [Indexed: 11/07/2022]
Abstract
Noncovalent interactions involving aromatic rings, such as π-stacking and CH/π interactions, are central to many areas of modern chemistry. However, recent studies proved that aromaticity is not required for stacking interactions, since similar interaction energies were computed for several aromatic and aliphatic dimers. Herein, the nature and origin of π/π, σ/σ, and σ/π dispersion interactions has been investigated by using dispersion-corrected density functional theory, energy decomposition analysis, and the recently developed noncovalent interaction (NCI) method. Our analysis shows that π/π and σ/σ stacking interactions are equally important for the benzene and cyclohexane dimers, explaining why both compounds have similar boiling points. Also, similar dispersion forces are found in the benzene⋅⋅⋅methane and cyclohexane⋅⋅⋅methane complexes. However, for systems larger than naphthalene, there are enhanced stacking interactions in the aromatic dimers adopting a parallel-displaced configuration compared to the analogous saturated systems. Although dispersion plays a decisive role in stabilizing all the complexes, the origin of the π/π, σ/σ, and σ/π interactions is different. The NCI method reveals that the dispersion interactions between the hydrogen atoms are responsible for the surprisingly strong aliphatic interactions. Moreover, whereas σ/σ and σ/π interactions are local, the π/π stacking are inherently delocalized, which give rise to a non-additive effect. These new types of dispersion interactions between saturated groups can be exploited in the rational design of novel carbon materials.
Collapse
Affiliation(s)
- Mercedes Alonso
- ALGC Research Group General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2 1050 Brussels (Belgium), Fax: (+32) 2-629-33-17.
| | | | | | | | | | | |
Collapse
|
37
|
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
|
38
|
Mukherjee S, Kailasam S, Bansal M, Bhattacharyya D. Energy hyperspace for stacking interaction inAU/AUdinucleotide step: Dispersion-corrected density functional theory study. Biopolymers 2013; 101:107-20. [DOI: 10.1002/bip.22289] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2013] [Revised: 05/13/2013] [Accepted: 05/13/2013] [Indexed: 12/20/2022]
Affiliation(s)
- Sanchita Mukherjee
- Biophysics Division; Saha Institute of Nuclear Physics; Kolkata 700064 India
| | - Senthilkumar Kailasam
- Molecular Biophysics Unit; Indian Institute of Science; Bangalore 560012 Karnataka India
| | - Manju Bansal
- Molecular Biophysics Unit; Indian Institute of Science; Bangalore 560012 Karnataka India
| | | |
Collapse
|
39
|
DiLabio GA, Koleini M, Torres E. Extension of the B3LYP–dispersion-correcting potential approach to the accurate treatment of both inter- and intra-molecular interactions. Theor Chem Acc 2013. [DOI: 10.1007/s00214-013-1389-x] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
40
|
Sedlak R, Janowski T, Pitoňák M, Řezáč J, Pulay P, Hobza P. The accuracy of quantum chemical methods for large noncovalent complexes. J Chem Theory Comput 2013; 9:3364-3374. [PMID: 24098094 PMCID: PMC3789125 DOI: 10.1021/ct400036b] [Citation(s) in RCA: 255] [Impact Index Per Article: 23.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
We evaluate the performance of the most widely used wavefunction, density functional theory, and semiempirical methods for the description of noncovalent interactions in a set of larger, mostly dispersion-stabilized noncovalent complexes (the L7 data set). The methods tested include MP2, MP3, SCS-MP2, SCS(MI)-MP2, MP2.5, MP2.X, MP2C, DFT-D, DFT-D3 (B3-LYP-D3, B-LYP-D3, TPSS-D3, PW6B95-D3, M06-2X-D3) and M06-2X, and semiempirical methods augmented with dispersion and hydrogen bonding corrections: SCC-DFTB-D, PM6-D, PM6-DH2 and PM6-D3H4. The test complexes are the octadecane dimer, the guanine trimer, the circumcoronene…adenine dimer, the coronene dimer, the guanine-cytosine dimer, the circumcoronene…guanine-cytosine dimer, and an amyloid fragment trimer containing phenylalanine residues. The best performing method is MP2.5 with relative root mean square deviation (rRMSD) of 4 %. It can thus be recommended as an alternative to the CCSD(T)/CBS (alternatively QCISD(T)/CBS) benchmark for molecular systems which exceed current computational capacity. The second best non-DFT method is MP2C with rRMSD of 8 %. A method with the most favorable "accuracy/cost" ratio belongs to the DFT family: BLYP-D3, with an rRMSD of 8 %. Semiempirical methods deliver less accurate results (the rRMSD exceeds 25 %). Nevertheless, their absolute errors are close to some much more expensive methods such as M06-2X, MP2 or SCS(MI)-MP2, and thus their price/performance ratio is excellent.
Collapse
Affiliation(s)
- Robert Sedlak
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University in Prague, Albertov 6, 128 43 Prague, Czech Republic
| | - Tomasz Janowski
- Department of Chemistry and Biochemistry, Fulbright College of Arts and Sciences, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Michal Pitoňák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská Dolina, 842 15 Bratislava, Slovak Republic
- Computing Center of the Slovak Academy of Sciences, Dúbravská cesta č. 9, 845 35 Bratislava, Slovak Republic
| | - Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
| | - Peter Pulay
- Department of Chemistry and Biochemistry, Fulbright College of Arts and Sciences, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Regional Center of Advanced Technologies and Materials, Department of Physical Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
41
|
Mottishaw JD, Sun H. Effects of Aromatic Trifluoromethylation, Fluorination, and Methylation on Intermolecular π–π Interactions. J Phys Chem A 2013; 117:7970-9. [DOI: 10.1021/jp403679x] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Jeffery D. Mottishaw
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United
States
| | - Haoran Sun
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United
States
| |
Collapse
|
42
|
Addicoat MA, Nishimura Y, Sato T, Tsuneda T, Irle S. Stochastic Search of Molecular Cluster Interaction Energy Surfaces with Coupled Cluster Quality Prediction. The Phenylacetylene Dimer. J Chem Theory Comput 2013; 9:3848-54. [DOI: 10.1021/ct4003515] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matthew A. Addicoat
- Department of Chemistry, Nagoya University, Furo-cho, Chikusa-ku Nagoya 464-4602,
Japan
| | - Yoshifumi Nishimura
- Department of Chemistry, Nagoya University, Furo-cho, Chikusa-ku Nagoya 464-4602,
Japan
| | - Takeshi Sato
- Photon Science Center, University of Tokyo, Tokyo 113-8656, Japan
| | - Takao Tsuneda
- Fuel Cell Nanomaterials Center, University of Yamanashi, Kofu 400-0021, Japan
| | - Stephan Irle
- Department of Chemistry, Nagoya University, Furo-cho, Chikusa-ku Nagoya 464-4602,
Japan
| |
Collapse
|
43
|
Sure R, Grimme S. Corrected small basis set Hartree-Fock method for large systems. J Comput Chem 2013; 34:1672-85. [PMID: 23670872 DOI: 10.1002/jcc.23317] [Citation(s) in RCA: 289] [Impact Index Per Article: 26.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2013] [Revised: 03/28/2013] [Accepted: 04/03/2013] [Indexed: 01/23/2023]
Abstract
A quantum chemical method based on a Hartree-Fock calculation with a small Gaussian AO basis set is presented. Its main area of application is the computation of structures, vibrational frequencies, and noncovalent interaction energies in huge molecular systems. The method is suggested as a partial replacement of semiempirical approaches or density functional theory (DFT) in particular when self-interaction errors are acute. In order to get accurate results three physically plausible atom pair-wise correction terms are applied for London dispersion interactions (D3 scheme), basis set superposition error (gCP scheme), and short-ranged basis set incompleteness effects. In total nine global empirical parameters are used. This so-called Hartee-Fock-3c (HF-3c) method is tested for geometries of small organic molecules, interaction energies and geometries of noncovalently bound complexes, for supramolecular systems, and protein structures. In the majority of realistic test cases good results approaching large basis set DFT quality are obtained at a tiny fraction of computational cost.
Collapse
Affiliation(s)
- Rebecca Sure
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115, Bonn, Germany
| | | |
Collapse
|
44
|
Remya K, Suresh CH. Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using Gaussian09. J Comput Chem 2013; 34:1341-53. [DOI: 10.1002/jcc.23263] [Citation(s) in RCA: 99] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2012] [Revised: 01/15/2013] [Accepted: 02/04/2013] [Indexed: 11/09/2022]
|
45
|
Grimme S, Steinmetz M. Effects of London dispersion correction in density functional theory on the structures of organic molecules in the gas phase. Phys Chem Chem Phys 2013; 15:16031-42. [DOI: 10.1039/c3cp52293h] [Citation(s) in RCA: 197] [Impact Index Per Article: 17.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|
46
|
Remko M, Duijnen PTV, Broer R. Structure and stability of complexes of charged structural units of heparin with arginine and lysine. RSC Adv 2013. [DOI: 10.1039/c2ra22062h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
47
|
Beran GJO, Wen S, Nanda K, Huang Y, Heit Y. Accurate and Robust Molecular Crystal Modeling Using Fragment-Based Electronic Structure Methods. Top Curr Chem (Cham) 2013; 345:59-93. [DOI: 10.1007/128_2013_502] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
48
|
Gunasinghe RN, Reuven DG, Suggs K, Wang XQ. Filled and Empty Orbital Interactions in a Planar Covalent Organic Framework on Graphene. J Phys Chem Lett 2012; 3:3048-3052. [PMID: 26292248 DOI: 10.1021/jz301304f] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The electronic characteristics of a planar covalent organic framework (COF) on graphene are investigated by means of dispersion-corrected density functional theory. The aromatic central molecule of the COF acts as an electron donor to graphene, while the linker of the COF acts as an electron acceptor. The concerted interaction between the filled orbitals of the central molecule and empty orbitals of the linker promotes the formation of planar COF networks on graphene. The calculation results are in very good agreement with experimental findings of an ordered hexagonal and square COF planar on graphene, which sheds light on the supermolecular assembly mechanism.
Collapse
Affiliation(s)
- Rosi N Gunasinghe
- Department of Physics and Center for Functional Nanoscale Materials, Clark Atlanta University, Atlanta, Georgia 30314, United States
| | - Darkeyah G Reuven
- Department of Physics and Center for Functional Nanoscale Materials, Clark Atlanta University, Atlanta, Georgia 30314, United States
| | - Kelvin Suggs
- Department of Physics and Center for Functional Nanoscale Materials, Clark Atlanta University, Atlanta, Georgia 30314, United States
| | - Xiao-Qian Wang
- Department of Physics and Center for Functional Nanoscale Materials, Clark Atlanta University, Atlanta, Georgia 30314, United States
| |
Collapse
|