1
|
Ferretti A, Canal L, Sorodoc RA, Sinha S, Brancato G. Fine Tuning the Intermolecular Interactions of Water Clusters Using the Dispersion-Corrected Density Functional Theory. Molecules 2023; 28:molecules28093834. [PMID: 37175249 PMCID: PMC10180381 DOI: 10.3390/molecules28093834] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Revised: 04/26/2023] [Accepted: 04/28/2023] [Indexed: 05/15/2023] Open
Abstract
Dispersion-inclusive density functional theory (DFT) methods have unequivocally demonstrated improved performances with respect to standard DFT approximations for modeling large and extended molecular systems at the quantum mechanical level. Yet, in some cases, disagreements with highly accurate reference calculations, such as CCSD(T) and quantum Monte Carlo (MC) calculations, still remain. Furthermore, the application of general-purpose corrections, such as the popular Grimme's semi-classical models (DFT-D), to different Kohn-Sham exchange-correlation functionals sometimes leads to variable and inconsistent results, which recommend a careful prior evaluation. In a recent study, we proposed a simple optimization protocol for enhancing the accuracy of these DFT-D methods by following an alternative and system-specific approach. Here, adopting the same computational strategy, we show how the accurate MC intermolecular interactions of a large set of water clusters of variable sizes (i.e., 300 (H2O)n structures, n = 9, 15, 27) can be reproduced remarkably well by dispersion-corrected DFT models (i.e., B3LYP-D4, PBE-D4, revPBE(0)-D4) upon re-optimization, reaching a mean absolute error per monomer of ~0.1 kcal/mol. Hence, the obtained results support the use of this procedure for fine-tuning tailored DFT-D models for the accurate description of targeted molecular systems.
Collapse
Affiliation(s)
- Alfonso Ferretti
- Scuola Normale Superiore and CSGI, Classe di Scienze, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Istituto Nazionale di Fisica Nucleare (INFN) Sezione di Pisa, Largo Bruno Pontecorvo 3, I-56127 Pisa, Italy
| | - Laura Canal
- Dipartimento di Ingegneria Civile ed Industriale, Università di Pisa, Largo Lucio Lazzarino 2, I-56124 Pisa, Italy
| | - Robert A Sorodoc
- Dipartimento di Ingegneria Civile ed Industriale, Università di Pisa, Largo Lucio Lazzarino 2, I-56124 Pisa, Italy
| | - Sourab Sinha
- Scuola Normale Superiore and CSGI, Classe di Scienze, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Giuseppe Brancato
- Scuola Normale Superiore and CSGI, Classe di Scienze, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Istituto Nazionale di Fisica Nucleare (INFN) Sezione di Pisa, Largo Bruno Pontecorvo 3, I-56127 Pisa, Italy
| |
Collapse
|
2
|
Li C, Paesani F, Voth GA. Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature. J Chem Theory Comput 2022; 18:2124-2131. [PMID: 35263110 PMCID: PMC9059465 DOI: 10.1021/acs.jctc.1c01223] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. We find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.
Collapse
Affiliation(s)
- Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, 60637
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, Materials Science and Engineering, and San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, 60637
| |
Collapse
|
3
|
Lindsey RK, Fried LE, Goldman N. Application of the ChIMES Force Field to Nonreactive Molecular Systems: Water at Ambient Conditions. J Chem Theory Comput 2018; 15:436-447. [DOI: 10.1021/acs.jctc.8b00831] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rebecca K. Lindsey
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Laurence E. Fried
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Nir Goldman
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
- Department of Chemical Engineering, University of California, Davis, California 95616, United States
| |
Collapse
|
4
|
Koziol L, Fried LE, Goldman N. Using Force Matching To Determine Reactive Force Fields for Water under Extreme Thermodynamic Conditions. J Chem Theory Comput 2016; 13:135-146. [DOI: 10.1021/acs.jctc.6b00707] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lucas Koziol
- Physical and Life Sciences
Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Laurence E. Fried
- Physical and Life Sciences
Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Nir Goldman
- Physical and Life Sciences
Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
5
|
Azadi S, Cohen RE. Chemical accuracy from quantum Monte Carlo for the benzene dimer. J Chem Phys 2016; 143:104301. [PMID: 26374029 DOI: 10.1063/1.4930137] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
Collapse
Affiliation(s)
- Sam Azadi
- Department of Earth Science and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - R E Cohen
- London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom and Extreme Materials Initiative, Geophysical Laboratory, Carnegie Institution of Washington, Washington, D.C. 20015, USA
| |
Collapse
|
6
|
Affiliation(s)
- Matúš Dubecký
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
7
|
Gillan MJ, Alfè D, Michaelides A. Perspective: How good is DFT for water? J Chem Phys 2016; 144:130901. [DOI: 10.1063/1.4944633] [Citation(s) in RCA: 478] [Impact Index Per Article: 59.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Michael J. Gillan
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Dario Alfè
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
8
|
Gillan MJ, Alfè D, Manby FR. Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations. J Chem Phys 2015; 143:102812. [PMID: 26374005 DOI: 10.1063/1.4926444] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.
Collapse
Affiliation(s)
- M J Gillan
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - D Alfè
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - F R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
9
|
Ohto T, Usui K, Hasegawa T, Bonn M, Nagata Y. Toward ab initio molecular dynamics modeling for sum-frequency generation spectra; an efficient algorithm based on surface-specific velocity-velocity correlation function. J Chem Phys 2015; 143:124702. [DOI: 10.1063/1.4931106] [Citation(s) in RCA: 84] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tatsuhiko Ohto
- Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Kota Usui
- Max-Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Taisuke Hasegawa
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyoku, Kyoto 606-8502, Japan
| | - Mischa Bonn
- Max-Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Yuki Nagata
- Max-Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
10
|
Medders GR, Götz AW, Morales MA, Bajaj P, Paesani F. On the representation of many-body interactions in water. J Chem Phys 2015; 143:104102. [DOI: 10.1063/1.4930194] [Citation(s) in RCA: 78] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gregory R. Medders
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, USA
| | - Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
11
|
Lobanova O, Avendaño C, Lafitte T, Müller EA, Jackson G. SAFT-γ force field for the simulation of molecular fluids: 4. A single-site coarse-grained model of water applicable over a wide temperature range. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1004804] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
12
|
Gillan MJ. Many-body exchange-overlap interactions in rare gases and water. J Chem Phys 2014; 141:224106. [DOI: 10.1063/1.4903240] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Cassone G, Giaquinta PV, Saija F, Saitta AM. Effect of Electric Field Orientation on the Mechanical and Electrical Properties of Water Ices: An Ab-initio Study. J Phys Chem B 2014; 118:12717-24. [DOI: 10.1021/jp507376v] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Giuseppe Cassone
- Dipartimento
di Fisica e di Scienze della Terra, Università degli Studi di Messina, Contrada Papardo, 98166 Messina, Italy
- CNR-IPCF, Viale Ferdinando Stagno d’Alcontres
37, 98158 Messina, Italy
- UMR 7590, IMPMC,
UPMC Univ Paris 06, Sorbonne Universités, F-75005 Paris, France
- UMR 7590, IMPMC, CNRS, F-75005 Paris, France
| | - Paolo V. Giaquinta
- Dipartimento
di Fisica e di Scienze della Terra, Università degli Studi di Messina, Contrada Papardo, 98166 Messina, Italy
| | - Franz Saija
- CNR-IPCF, Viale Ferdinando Stagno d’Alcontres
37, 98158 Messina, Italy
| | - A. Marco Saitta
- UMR 7590, IMPMC,
UPMC Univ Paris 06, Sorbonne Universités, F-75005 Paris, France
- UMR 7590, IMPMC, CNRS, F-75005 Paris, France
| |
Collapse
|