1
|
Sitkiewicz S, Zaleśny R, Ramos-Cordoba E, Luis JM, Matito E. How Reliable Are Modern Density Functional Approximations to Simulate Vibrational Spectroscopies? J Phys Chem Lett 2022; 13:5963-5968. [PMID: 35735354 PMCID: PMC9251762 DOI: 10.1021/acs.jpclett.2c01278] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
We show that properties of molecules with low-frequency modes calculated with density functional approximations (DFAs) suffer from spurious oscillations along the nuclear displacement coordinate due to numerical integration errors. Occasionally, the problem can be alleviated using extensive integration grids that compromise the favorable cost-accuracy ratio of DFAs. Since spurious oscillations are difficult to predict or identify, DFAs are exposed to severe performance errors in IR and Raman intensities and frequencies or vibrational contributions to any molecular property. Using Fourier spectral analysis and digital signal processing techniques, we identify and quantify the error due to these oscillations for 45 widely used DFAs. LC-BLYP and BH&H are revealed as the only functionals showing robustness against the spurious oscillations of various energy, dipole moment, and polarizability derivatives with respect to a nuclear displacement coordinate. Given the ubiquitous nature of molecules with low-frequency modes, we warrant caution in using modern DFAs to simulate vibrational spectroscopies.
Collapse
Affiliation(s)
- Sebastian
P. Sitkiewicz
- Donostia
International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain
- Polimero
eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Unibertsitatea UPV/EHU, P.K. 1072, 20080 Donostia, Euskadi, Spain
| | - Robert Zaleśny
- Faculty
of Chemistry, Wrocław University of
Science and Technology, Wyb. Wyspiańskiego 27, PL−50370 Wrocław, Poland
| | - Eloy Ramos-Cordoba
- Donostia
International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain
- Polimero
eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Unibertsitatea UPV/EHU, P.K. 1072, 20080 Donostia, Euskadi, Spain
| | - Josep M. Luis
- Institut
de Química Computacional i Catàlisi (IQCC) and Departament
de Química, Universitat de Girona, 17003 Girona, Catalonia, Spain
| | - Eduard Matito
- Donostia
International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain
- Ikerbasque
Foundation for Science, Plaza Euskadi 5, 48009 Bilbao, Euskadi, Spain
| |
Collapse
|
2
|
Hinsch JJ, Liu J, Wang Y. Reinvestigating oxygen adsorption on Ag(111) by using strongly constrained and appropriately normed semi-local density functional with the revised Vydrov van Voorhis van der Waals force correction. J Chem Phys 2021; 155:234704. [PMID: 34937376 DOI: 10.1063/5.0073407] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
While density functional theory (DFT) at the generalized gradient approximation (GGA) level has made great success in catalysis, it fails in some important systems such as the adsorption of the oxygen molecule on the Ag(111) surface. Previous DFT studies at the GGA level revealed theoretical inconsistencies on the adsorption energies and dissociation barriers of O2 on Ag(111) in comparison with the experimental conclusion. In this study, the strongly constrained and appropriately normed-revised Vydrov van Voorhis van der Waals correction functional (SCAN-rVV10) method at the meta-GGA level with the nonlocal van der Waals (vdW) force correction was used to reinvestigate the adsorption properties of O2 on the Ag(111) surface. The SCAN-rVV10 results successfully confirm the experimental observation that both molecular and dissociative adsorptions can exist for oxygen on Ag(111). The calculated adsorption energy for the physisorption state and the relevant dissociation energy barrier are close to the experimental data. It demonstrates that SCAN-rVV10 can outperform functionals at the GGA level for O2/Ag(111). Therefore, our findings suggest that SCAN-rVV10 can be the desired method for systems where the correct description of intermediate-ranged vdW forces is essential, such as the physisorption of small molecules on the solid surface.
Collapse
Affiliation(s)
- Jack J Hinsch
- School of Environment and Science, Centre for Catalysis and Clean Energy, Griffith University, Gold Coast QLD 4222, Australia
| | - Junxian Liu
- School of Environment and Science, Centre for Catalysis and Clean Energy, Griffith University, Gold Coast QLD 4222, Australia
| | - Yun Wang
- School of Environment and Science, Centre for Catalysis and Clean Energy, Griffith University, Gold Coast QLD 4222, Australia
| |
Collapse
|
3
|
Sparrow ZM, Ernst BG, Joo PT, Lao KU, DiStasio RA. NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts. J Chem Phys 2021; 155:184303. [PMID: 34773949 DOI: 10.1063/5.0068862] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we present NENCI-2021, a benchmark database of ∼8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation-π and anion-π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion-π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (Eint) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The Eint values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) Eint value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent-these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation-π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the Eint values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive Eint values, and maximum errors of ±0.2-0.3 kcal/mol (i.e., ∼±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.
Collapse
Affiliation(s)
- Zachary M Sparrow
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Brian G Ernst
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Paul T Joo
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Ka Un Lao
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
4
|
Price AJA, Bryenton KR, Johnson ER. Requirements for an accurate dispersion-corrected density functional. J Chem Phys 2021; 154:230902. [PMID: 34241263 DOI: 10.1063/5.0050993] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Post-self-consistent dispersion corrections are now the norm when applying density-functional theory to systems where non-covalent interactions play an important role. However, there is a wide range of base functionals and dispersion corrections available from which to choose. In this work, we opine on the most desirable requirements to ensure that both the base functional and dispersion correction, individually, are as accurate as possible for non-bonded repulsion and dispersion attraction. The base functional should be dispersionless, numerically stable, and involve minimal delocalization error. Simultaneously, the dispersion correction should include finite damping, higher-order pairwise dispersion terms, and electronic many-body effects. These criteria are essential for avoiding reliance on error cancellation and obtaining correct results from correct physics.
Collapse
Affiliation(s)
- Alastair J A Price
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd., Halifax, Nova Scotia B3H 4R2, Canada
| | - Kyle R Bryenton
- Department of Physics and Atmospheric Science, Dalhousie University, 6310 Coburg Rd., Halifax, Nova Scotia B3H 4R2, Canada
| | - Erin R Johnson
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd., Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
5
|
Chakraborty D, Berland K, Thonhauser T. Next-Generation Nonlocal van der Waals Density Functional. J Chem Theory Comput 2020; 16:5893-5911. [PMID: 32786912 DOI: 10.1021/acs.jctc.0c00471] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The fundamental ideas for a nonlocal density functional theory-capable of reliably capturing van der Waals interactions-were already conceived in the 1990s. In 2004, a seminal paper introduced the first practical nonlocal exchange-correlation functional called vdW-DF, which has become widely successful and laid the foundation for much further research. However, since then, the functional form of vdW-DF has remained unchanged. Several successful modifications paired the original functional with different (local) exchange functionals to improve performance, and the successor vdW-DF2 also updated one internal parameter. Bringing together different insights from almost 2 decades of development and testing, we present the next-generation nonlocal correlation functional called vdW-DF3, in which we change the functional form while staying true to the original design philosophy. Although many popular functionals show good performance around the binding separation of van der Waals complexes, they often result in significant errors at larger separations. With vdW-DF3, we address this problem by taking advantage of a recently uncovered and largely unconstrained degree of freedom within the vdW-DF framework that can be constrained through empirical input, making our functional semiempirical. For two different parameterizations, we benchmark vdW-DF3 against a large set of well-studied test cases and compare our results with the most popular functionals, finding good performance in general for a wide array of systems and a significant improvement in accuracy at larger separations. Finally, we discuss the achievable performance within the current vdW-DF framework, the flexibility in functional design offered by vdW-DF3, as well as possible future directions for nonlocal van der Waals density functional theory.
Collapse
Affiliation(s)
- D Chakraborty
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina 27109, United States.,Center for Functional Materials, Wake Forest University, Winston-Salem, North Carolina 27109, United States
| | - K Berland
- Faculty of Science and Technology, Norwegian University of Life Sciences, 1430 Ås, Norway
| | - T Thonhauser
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina 27109, United States.,Center for Functional Materials, Wake Forest University, Winston-Salem, North Carolina 27109, United States
| |
Collapse
|
6
|
Brémond É, Ciofini I, Sancho-García JC, Adamo C. Double-Hybrid Functionals and Tailored Basis Set: Fullerene (C60) Dimer and Isomers as Test Cases. J Phys Chem A 2019; 123:10040-10046. [DOI: 10.1021/acs.jpca.9b06536] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Éric Brémond
- ITODYS, UMR CNRS 7086, Université de Paris, 15 rue J.-A. de Baïf, F-75013 Paris, France
| | - Ilaria Ciofini
- Institut de Recherche de Chimie Paris, PSL Research University, CNRS, Chimie ParisTech, 11, rue Pierre et Marie Curie, F-75005 Paris, France
| | | | - Carlo Adamo
- Institut de Recherche de Chimie Paris, PSL Research University, CNRS, Chimie ParisTech, 11, rue Pierre et Marie Curie, F-75005 Paris, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
7
|
Garcı́a JS, Brémond É, Campetella M, Ciofini I, Adamo C. Small Basis Set Allowing the Recovery of Dispersion Interactions with Double-Hybrid Functionals. J Chem Theory Comput 2019; 15:2944-2953. [DOI: 10.1021/acs.jctc.8b01203] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Juan Sanz Garcı́a
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, 11, rue Pierre et Marie Curie, F-75005 Paris, France
| | - Éric Brémond
- Univ Paris Diderot, Sorbonne Paris Cité, ITODYS, UMR CNRS 7086, 15 rue J.-A. de Baïf, F-75013 Paris, France
| | - Marco Campetella
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, 11, rue Pierre et Marie Curie, F-75005 Paris, France
| | - Ilaria Ciofini
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, 11, rue Pierre et Marie Curie, F-75005 Paris, France
| | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, 11, rue Pierre et Marie Curie, F-75005 Paris, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
8
|
Li M, Reimers JR, Dobson JF, Gould T. Faraday cage screening reveals intrinsic aspects of the van der Waals attraction. Proc Natl Acad Sci U S A 2018; 115:E10295-E10302. [PMID: 30327347 PMCID: PMC6217410 DOI: 10.1073/pnas.1811569115] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
General properties of the recently observed screening of the van der Waals (vdW) attraction between a silica substrate and silica tip by insertion of graphene are predicted using basic theory and first-principles calculations. Results are then focused on possible practical applications, as well as an understanding of the nature of vdW attraction, considering recent discoveries showing it competing against covalent and ionic bonding. The traditional view of the vdW attraction as arising from pairwise-additive London dispersion forces is considered using Grimme's "D3" method, comparing results to those from Tkatchenko's more general many-body dispersion (MBD) approach, all interpreted in terms of Dobson's general dispersion framework. Encompassing the experimental results, MBD screening of the vdW force between two silica bilayers is shown to scale up to medium separations as 1.25 de/d, where d is the bilayer separation and de is its equilibrium value, depicting antiscreening approaching and inside de Means of unifying this correlation effect with those included in modern density functionals are urgently required.
Collapse
Affiliation(s)
- Musen Li
- International Centre for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
- Department of Physics, Shanghai University, Shanghai 200444, China
| | - Jeffrey R Reimers
- International Centre for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China;
- Department of Physics, Shanghai University, Shanghai 200444, China
- School of Mathematical and Physical Sciences, University of Technology Sydney, Ultimo, NSW 2007, Australia
| | - John F Dobson
- School of Natural Sciences, Griffith University, Nathan, QLD 4111, Australia
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan, QLD 4111, Australia
| | - Tim Gould
- School of Natural Sciences, Griffith University, Nathan, QLD 4111, Australia;
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan, QLD 4111, Australia
| |
Collapse
|
9
|
Walters ET, Mohebifar M, Johnson ER, Rowley CN. Evaluating the London Dispersion Coefficients of Protein Force Fields Using the Exchange-Hole Dipole Moment Model. J Phys Chem B 2018; 122:6690-6701. [DOI: 10.1021/acs.jpcb.8b02814] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Evan T. Walters
- Department of Chemistry, Memorial University of Newfoundland, St. John’s A1C 5S7, Newfoundland and Labrador, Canada
| | - Mohamad Mohebifar
- Department of Chemistry, Memorial University of Newfoundland, St. John’s A1C 5S7, Newfoundland and Labrador, Canada
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie University, Halifax B3H 4R2, Nova Scotia, Canada
| | - Christopher N. Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John’s A1C 5S7, Newfoundland and Labrador, Canada
| |
Collapse
|