1
|
Trinter F, Miteva T, Weller M, Hartung A, Richter M, Williams JB, Gatton A, Gaire B, Sartor J, Landers AL, Berry B, Ben-Itzhak I, Sisourat N, Stumpf V, Gokhberg K, Dörner R, Jahnke T, Weber T. Ultrafast temporal evolution of interatomic Coulombic decay in NeKr dimers. Chem Sci 2022; 13:1789-1800. [PMID: 35282626 PMCID: PMC8827086 DOI: 10.1039/d1sc04630f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2021] [Accepted: 12/28/2021] [Indexed: 11/21/2022] Open
Abstract
We investigate interatomic Coulombic decay in NeKr dimers after neon inner-valence photoionization [Ne+(2s-1)] using a synchrotron light source. We measure with high energy resolution the two singly charged ions of the Coulomb-exploding dimer dication and the photoelectron in coincidence. By carefully tracing the post-collision interaction between the photoelectron and the emitted ICD electron we are able to probe the temporal evolution of the state as it decays. Although the ionizing light pulses are 80 picoseconds long, we determine the lifetime of the intermediate dimer cation state and visualize the contraction of the nuclear structure on the femtosecond time scale.
Collapse
Affiliation(s)
- F Trinter
- Institut für Kernphysik, Goethe-Universität 60438 Frankfurt am Main Germany .,Molecular Physics, Fritz-Haber-Institut der Max-Planck-Gesellschaft 14195 Berlin Germany
| | - T Miteva
- Laboratoire de Chimie Physique Matière et Rayonnement, UMR 7614, Sorbonne Université, CNRS 75005 Paris France
| | - M Weller
- Institut für Kernphysik, Goethe-Universität 60438 Frankfurt am Main Germany .,Lawrence Berkeley National Laboratory, Chemical Sciences Division Berkeley California 94720 USA
| | - A Hartung
- Institut für Kernphysik, Goethe-Universität 60438 Frankfurt am Main Germany
| | - M Richter
- Institut für Kernphysik, Goethe-Universität 60438 Frankfurt am Main Germany
| | - J B Williams
- Department of Physics, University of Nevada Reno Nevada 89557 USA
| | - A Gatton
- Lawrence Berkeley National Laboratory, Chemical Sciences Division Berkeley California 94720 USA .,Department of Physics, Auburn University Auburn Alabama 36849 USA
| | - B Gaire
- Lawrence Berkeley National Laboratory, Chemical Sciences Division Berkeley California 94720 USA
| | - J Sartor
- Department of Physics, Auburn University Auburn Alabama 36849 USA
| | - A L Landers
- Department of Physics, Auburn University Auburn Alabama 36849 USA
| | - B Berry
- J. R. Macdonald Laboratory, Department of Physics, Kansas State University Manhattan Kansas 66506 USA
| | - I Ben-Itzhak
- J. R. Macdonald Laboratory, Department of Physics, Kansas State University Manhattan Kansas 66506 USA
| | - N Sisourat
- Laboratoire de Chimie Physique Matière et Rayonnement, UMR 7614, Sorbonne Université, CNRS 75005 Paris France
| | - V Stumpf
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg 69120 Heidelberg Germany
| | - K Gokhberg
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg 69120 Heidelberg Germany
| | - R Dörner
- Institut für Kernphysik, Goethe-Universität 60438 Frankfurt am Main Germany
| | - T Jahnke
- European XFEL GmbH 22869 Schenefeld Germany
| | - T Weber
- Lawrence Berkeley National Laboratory, Chemical Sciences Division Berkeley California 94720 USA
| |
Collapse
|
2
|
Řezáč J. Non-Covalent Interactions Atlas Benchmark Data Sets 5: London Dispersion in an Extended Chemical Space. Phys Chem Chem Phys 2022; 24:14780-14793. [DOI: 10.1039/d2cp01602h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Non-Covalent Interactions Atlas (www.nciatlas.org) has been extended with two data sets of benchmark interaction energies in complexes dominated by London dispersion. The D1200 data set of equilibrium geometries provides...
Collapse
|
3
|
Jerabek P, Smits OR, Mewes JM, Peterson KA, Schwerdtfeger P. Solid Oganesson via a Many-Body Interaction Expansion Based on Relativistic Coupled-Cluster Theory and from Plane-Wave Relativistic Density Functional Theory. J Phys Chem A 2019; 123:4201-4211. [DOI: 10.1021/acs.jpca.9b01947] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Paul Jerabek
- Department for Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung (KOFO), Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Odile R. Smits
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, 0632 Auckland, New Zealand
| | - Jan-Michael Mewes
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, 0632 Auckland, New Zealand
| | - Kirk A. Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164, United States
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, 0632 Auckland, New Zealand
| |
Collapse
|
4
|
Abstract
The potential energy surfaces for the interior rotation of a series of pairs of noble gas atoms encapsulated in the C70 cavity have been explored.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore
- Classe di Scienze
- 56126 Pisa
- Italy
| | - Albeiro Restrepo
- Instituto de Química
- Universidad de Antioquia UdeA
- Medellín
- Colombia
| |
Collapse
|
5
|
Knochenmuss R, Sinha RK, Poblotzki A, Den T, Leutwyler S. Intermolecular dissociation energies of hydrogen-bonded 1-naphthol complexes. J Chem Phys 2018; 149:204311. [DOI: 10.1063/1.5055720] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Richard Knochenmuss
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Rajeev K. Sinha
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Anja Poblotzki
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Takuya Den
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Samuel Leutwyler
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| |
Collapse
|
6
|
Knochenmuss R, Maity S, Balmer F, Müller C, Leutwyler S. Intermolecular dissociation energies of 1-naphthol·n-alkane complexes. J Chem Phys 2018; 149:034306. [DOI: 10.1063/1.5034110] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Richard Knochenmuss
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Surajit Maity
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Franziska Balmer
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Charlotte Müller
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Samuel Leutwyler
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| |
Collapse
|
7
|
Vuckovic S, Gori-Giorgi P, Della Sala F, Fabiano E. Restoring Size Consistency of Approximate Functionals Constructed from the Adiabatic Connection. J Phys Chem Lett 2018; 9:3137-3142. [PMID: 29787273 PMCID: PMC5994725 DOI: 10.1021/acs.jpclett.8b01054] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2018] [Accepted: 05/22/2018] [Indexed: 05/24/2023]
Abstract
Approximate exchange-correlation functionals built by modeling in a nonlinear way the adiabatic connection (AC) integrand of density functional theory have many attractive features, being virtually parameter-free and satisfying different exact properties, but they also have a fundamental flaw: they violate the size-consistency condition, crucial to evaluate interaction energies of molecular systems. We show that size consistency in the AC-based functionals can be restored in a very simple way at no extra computational cost. Results on a large set of benchmark molecular interaction energies show that functionals based on the interaction strength interpolation approximations are significantly more accurate than second-order perturbation theory.
Collapse
Affiliation(s)
- Stefan Vuckovic
- Department
of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling,
FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Paola Gori-Giorgi
- Department
of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling,
FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Fabio Della Sala
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center
for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| | - Eduardo Fabiano
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center
for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| |
Collapse
|
8
|
Knochenmuss R, Sinha RK, Leutwyler S. Intermolecular dissociation energies of dispersively bound complexes of aromatics with noble gases and nitrogen. J Chem Phys 2018; 148:134302. [DOI: 10.1063/1.5019432] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Richard Knochenmuss
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Rajeev K. Sinha
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Samuel Leutwyler
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| |
Collapse
|
9
|
Affiliation(s)
- Robert A. Shaw
- Department of Chemistry, University of Sheffield, Sheffield, UK
| | - J. Grant Hill
- Department of Chemistry, University of Sheffield, Sheffield, UK
| |
Collapse
|
10
|
Roch LM, Baldridge KK. General optimization procedure towards the design of a new family of minimal parameter spin-component-scaled double-hybrid density functional theory. Phys Chem Chem Phys 2017; 19:26191-26200. [PMID: 28930316 DOI: 10.1039/c7cp04125j] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A general optimization procedure towards the development and implementation of a new family of minimal parameter spin-component-scaled double-hybrid (mSD) density functional theory (DFT) is presented. The nature of the proposed exchange-correlation functional establishes a methodology with minimal empiricism. This new family of double-hybrid (DH) density functionals is demonstrated using the PBEPBE functional, illustrating the optimization procedure to the mSD-PBEPBE method, and the performance characteristics shown for a set of non-covalent complexes covering a broad regime of weak interactions. With only two parameters, mSD-PBEPBE and its cost-effective counterpart, RI-mSD-PBEPBE, show a mean absolute error of ca. 0.4 kcal mol-1 averaged over 66 weak interacting systems. Following a successive 2D-grid refinement for a CBS extrapolation of the coefficients, the optimization procedure can be recommended for the design and implementation of a variety of additional DH methods using any of the plethora of currently available functionals.
Collapse
Affiliation(s)
- Loïc M Roch
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland.
| | | |
Collapse
|
11
|
Roch LM, Baldridge KK. Dispersion-Corrected Spin-Component-Scaled Double-Hybrid Density Functional Theory: Implementation and Performance for Non-covalent Interactions. J Chem Theory Comput 2017; 13:2650-2666. [PMID: 28537392 DOI: 10.1021/acs.jctc.7b00220] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The implementation of 300 combinations of generalized gradient approximation/local density approximation exchange-correlation dispersion-corrected spin-component-scaled double-hybrid (DSD) density functional theory (DFT) methods has been carried out and the performance assessed against several DFT and post-Hartree-Fock methods, enabling further advancements toward the long-standing challenge of accurate prediction of interaction energies and associated properties. The resulting framework is flexible and has been further extended to include the resolution of identity (RI) approximation for solving the critical four-center two-electron repulsion integrals in the basis of the Kohn-Sham orbitals for cost effectiveness. To evaluate the performance of this set of new cost-effective methods, denoted as RI-DSD-DFTs, seven validation data sets were designed to cover a broad range of non-covalent interactions with characteristic stabilizing contributions. Inclusion of the perturbative treatment of correlation effects is shown to significantly improve the description of weak interactions. The set of DSD-DFTs provide interaction energies with root-mean-square deviations and mean absolute errors within 0.5 kcal/mol. The cost-effective RI-DSD-DFT counterparts deviate by less than 0.18 kcal/mol on average with only 2% of the computational cost.
Collapse
Affiliation(s)
- Loïc M Roch
- Department of Chemistry, University of Zurich , Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Kim K Baldridge
- Department of Chemistry, University of Zurich , Winterthurerstrasse 190, 8057 Zurich, Switzerland.,School of Pharmaceutical Science and Technology, University of Tianjin , 92 Weijin Road, Nankai District, Tianjin 3000072, P. R. China
| |
Collapse
|
12
|
Zhao Y, Truhlar DG. Benchmark Databases for Nonbonded Interactions and Their Use To Test Density Functional Theory. J Chem Theory Comput 2015; 1:415-32. [PMID: 26641508 DOI: 10.1021/ct049851d] [Citation(s) in RCA: 800] [Impact Index Per Article: 88.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
We present four benchmark databases of binding energies for nonbonded complexes. Four types of nonbonded interactions are considered: hydrogen bonding, charge transfer, dipole interactions, and weak interactions. We tested 44 DFT methods and 1 WFT method against the new databases; one of the DFT methods (PBE1KCIS) is new, and all of the other methods are from the literature. Among the tested methods, the PBE, PBE1PBE, B3P86, MPW1K, B97-1, and BHandHLYP functionals give the best performance for hydrogen bonding. MPWB1K, MP2, MPW1B95, MPW1K, and BHandHLYP give the best performances for charge-transfer interactions, and MPW3LYP, B97-1, PBE1KCIS, B98, and PBE1PBE give the best performance for dipole interactions. Finally, MP2, B97-1, MPWB1K, PBE1KCIS, and MPW1B95 give the best performance for weak interactions. Overall, MPWB1K is the best of all the tested DFT methods, with a relative error (highly averaged) of only 11%, and MPW1K, PBE1PBE, and B98 are the best of the tested DFT methods that do not contain kinetic energy density. Moving up the rungs of Jacob's ladder for nonempirical DFT, PBE improves significantly over the LSDA, and TPSS improve slightly over PBE (on average) for nonbonded interactions.
Collapse
Affiliation(s)
- Yan Zhao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
13
|
Shirkov L, Makarewicz J. Does DFT-SAPT method provide spectroscopic accuracy? J Chem Phys 2015; 142:064102. [DOI: 10.1063/1.4907204] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
14
|
Patra A, Xiao B, Perdew JP. Short-Range Cut-Off of the Summed-Up van der Waals Series: Rare-Gas Dimers. Top Curr Chem (Cham) 2015. [DOI: 10.1007/128_2015_625] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
|
15
|
Bousquet D, Brémond E, Sancho-García JC, Ciofini I, Adamo C. Non-parametrized functionals with empirical dispersion corrections: A happy match? Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1602-6] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
16
|
Zobel JP, Kryzhevoi NV, Pernpointner M. Communication: electron transfer mediated decay enabled by spin-orbit interaction in small krypton/xenon clusters. J Chem Phys 2014; 140:161103. [PMID: 24784242 DOI: 10.1063/1.4873134] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
In this work we study the influence of relativistic effects, in particular spin-orbit coupling, on electronic decay processes in KrXe2 clusters of various geometries. For the first time it is shown that inclusion of spin-orbit coupling has decisive influence on the accessibility of a specific decay pathway in these clusters. The radiationless relaxation process is initiated by a Kr 4s ionization followed by an electron transfer from xenon to krypton and a final second ionization of the system. We demonstrate the existence of competing electronic decay pathways depending in a subtle way on the geometry and level of theory. For our calculations a fully relativistic framework was employed where omission of spin-orbit coupling leads to closing of two decay pathways. These findings stress the relevance of an adequate relativistic description for clusters with heavy elements and their fragmentation dynamics.
Collapse
Affiliation(s)
- J Patrick Zobel
- Physikalisch-Chemisches Institut, Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| | - Nikolai V Kryzhevoi
- Physikalisch-Chemisches Institut, Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| | - Markus Pernpointner
- Physikalisch-Chemisches Institut, Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
17
|
Niimi K, Nakayama A, Ono Y, Taketsugu T. Vibrational Shifts of HXeCl in Matrix Environments. J Phys Chem A 2014; 118:380-7. [DOI: 10.1021/jp411298p] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Keisuke Niimi
- Graduate
School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-0810, Japan
| | - Akira Nakayama
- Catalysis
Research Center, Hokkaido University, Sapporo 001-0021, Japan
| | - Yuriko Ono
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
| |
Collapse
|
18
|
Bousquet D, Brémond E, Sancho-García JC, Ciofini I, Adamo C. Is There Still Room for Parameter Free Double Hybrids? Performances of PBE0-DH and B2PLYP over Extended Benchmark Sets. J Chem Theory Comput 2013; 9:3444-52. [DOI: 10.1021/ct400358f] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Diane Bousquet
- Laboratoire d’Electrochimie,
Chimie des Interfaces et Modélisation pour l’Energie,
CNRS UMR-7575, Chimie ParisTech, 11 rue
P. et M. Curie, F-75231 Paris Cedex 05, France
| | - Eric Brémond
- Laboratoire d’Electrochimie,
Chimie des Interfaces et Modélisation pour l’Energie,
CNRS UMR-7575, Chimie ParisTech, 11 rue
P. et M. Curie, F-75231 Paris Cedex 05, France
| | | | - Ilaria Ciofini
- Laboratoire d’Electrochimie,
Chimie des Interfaces et Modélisation pour l’Energie,
CNRS UMR-7575, Chimie ParisTech, 11 rue
P. et M. Curie, F-75231 Paris Cedex 05, France
| | - Carlo Adamo
- Laboratoire d’Electrochimie,
Chimie des Interfaces et Modélisation pour l’Energie,
CNRS UMR-7575, Chimie ParisTech, 11 rue
P. et M. Curie, F-75231 Paris Cedex 05, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris,
France
| |
Collapse
|
19
|
Cornaton Y, Franck O, Teale AM, Fromager E. Analysis of double-hybrid density functionals along the adiabatic connection. Mol Phys 2013. [DOI: 10.1080/00268976.2013.783640] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Yann Cornaton
- a Laboratoire de Chimie Quantique, Institut de Chimie , CNRS/Université de Strasbourg , Strasbourg , France
| | - Odile Franck
- a Laboratoire de Chimie Quantique, Institut de Chimie , CNRS/Université de Strasbourg , Strasbourg , France
| | - Andrew M. Teale
- b School of Chemistry , University of Nottingham , University Park, Nottingham , United Kingdom
- c Department of Chemistry, Centre for Theoretical and Computational Chemistry , University of Oslo , Blindern, Oslo , Norway
| | - Emmanuel Fromager
- a Laboratoire de Chimie Quantique, Institut de Chimie , CNRS/Université de Strasbourg , Strasbourg , France
| |
Collapse
|
20
|
Austin A, Petersson GA, Frisch MJ, Dobek FJ, Scalmani G, Throssell K. A Density Functional with Spherical Atom Dispersion Terms. J Chem Theory Comput 2012; 8:4989-5007. [PMID: 26593191 DOI: 10.1021/ct300778e] [Citation(s) in RCA: 362] [Impact Index Per Article: 30.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A new hybrid density functional, APF, is introduced, which avoids the spurious long-range attractive or repulsive interactions that are found in most density functional theory (DFT) models. It therefore provides a sound baseline for the addition of an empirical dispersion correction term, which is developed from a spherical atom model (SAM). The APF-D empirical dispersion model contains nine adjustable parameters that were selected based on a very small training set (15 noble gas dimers and 4 small hydrocarbon dimers), along with two computed atomic properties (ionization potential and effective atomic polarizability) for each element. APF-D accurately describes a large portion of the potential energy surfaces of complexes of noble gas atoms with various diatomic molecules involving a wide range of elements and of dimers of small hydrocarbons, and it reproduces the relative conformational energies of organic molecules. The accuracy for these weak interactions is comparable to that of CCSD(T)/aug-cc-pVTZ calculations. The accuracy in predicting the geometry of hydrogen bond complexes is competitive with other models involving DFT and empirical dispersion.
Collapse
Affiliation(s)
- Amy Austin
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| | - George A Petersson
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| | - Michael J Frisch
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States.,Gaussian, Inc. , 340 Quinnipiac Street Building 40, Wallingford, Connecticut 06492, United States
| | - Frank J Dobek
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| | - Giovanni Scalmani
- Gaussian, Inc. , 340 Quinnipiac Street Building 40, Wallingford, Connecticut 06492, United States
| | - Kyle Throssell
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| |
Collapse
|
21
|
Kullie O, Saue T. Range-separated density functional theory: A 4-component relativistic study of the rare gas dimers He2, Ne2, Ar2, Kr2, Xe2, Rn2 and Uuo2. Chem Phys 2012. [DOI: 10.1016/j.chemphys.2011.06.024] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
22
|
Devdariani AZ, Belyaev AK, Alekseev AB, Liebermann HP, Buenker RJ. Lifetime of weakly bound He(21 S)-Ne molecules. Temperature dependence. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY B 2012. [DOI: 10.1134/s1990793111060066] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
23
|
Goerigk L, Kruse H, Grimme S. Benchmarking Density Functional Methods against the S66 and S66x8 Datasets for Non-Covalent Interactions. Chemphyschem 2011; 12:3421-33. [DOI: 10.1002/cphc.201100826] [Citation(s) in RCA: 259] [Impact Index Per Article: 19.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2011] [Indexed: 11/11/2022]
|
24
|
Martinez-Casado R, Mallia G, Harrison NM. Ab initio calculation of the MgO(100) interaction with He and Ne: a HF + MP2 and HF + MP2(B3LYP) comparison. Chem Commun (Camb) 2011; 47:11630-2. [PMID: 21952116 DOI: 10.1039/c1cc14623h] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Second order Rayleigh Schrödinger perturbation theory is applied to calculate the correlation energy contribution to the London dispersion interaction to approximate the interaction of the He and Ne with the MgO(100) surface; single particle orbitals using either Hartree-Fock theory or hybrid-exchange density functional theory are used as the reference state.
Collapse
Affiliation(s)
- Ruth Martinez-Casado
- Thomas Young Centre, Department of Chemistry, Imperial College London, South Kensington, London, SW7 2AZ, UK.
| | | | | |
Collapse
|
25
|
Tonner R, Frenking G, Lein M, Schwerdtfeger P. Packed to the Rafters: Filling up C60 with Rare Gas Atoms. Chemphyschem 2011; 12:2081-4. [DOI: 10.1002/cphc.201100360] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2011] [Indexed: 11/09/2022]
|
26
|
Arion T, Mucke M, Förstel M, Bradshaw AM, Hergenhahn U. Interatomic Coulombic decay in mixed NeKr clusters. J Chem Phys 2011; 134:074306. [PMID: 21341843 DOI: 10.1063/1.3552082] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tiberiu Arion
- Max-Planck-Institut für Plasmaphysik, EURATOM Association, Garching bei München, Germany
| | | | | | | | | |
Collapse
|
27
|
Goerigk L, Grimme S. A thorough benchmark of density functional methods for general main group thermochemistry, kinetics, and noncovalent interactions. Phys Chem Chem Phys 2011; 13:6670-88. [PMID: 21384027 DOI: 10.1039/c0cp02984j] [Citation(s) in RCA: 1310] [Impact Index Per Article: 100.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Affiliation(s)
- Lars Goerigk
- Theoretische Organische Chemie, Organisch-Chemisches Institut der Universität Münster, Corrensstraße 40, D-48149 Münster, Germany
| | | |
Collapse
|
28
|
Goerigk L, Grimme S. Efficient and Accurate Double-Hybrid-Meta-GGA Density Functionals-Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions. J Chem Theory Comput 2010; 7:291-309. [PMID: 26596152 DOI: 10.1021/ct100466k] [Citation(s) in RCA: 847] [Impact Index Per Article: 60.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We present an extended and improved version of our recently published database for general main group thermochemistry, kinetics, and noncovalent interactions [J. Chem. Theory Comput. 2010, 6, 107], which is dubbed GMTKN30. Furthermore, we suggest and investigate two new double-hybrid-meta-GGA density functionals called PTPSS-D3 and PWPB95-D3. PTPSS-D3 is based on reparameterized TPSS exchange and correlation contributions; PWPB95-D3 contains reparameterized PW exchange and B95 parts. Both functionals contain fixed amounts of 50% Fock-exchange. Furthermore, they include a spin-opposite scaled perturbative contribution and are combined with our latest atom-pairwise London-dispersion correction [J. Chem. Phys. 2010, 132, 154104]. When evaluated with the help of the Laplace transformation algorithm, both methods scale as N(4) with system size. The functionals are compared with the double hybrids B2PLYP-D3, B2GPPLYP-D3, DSD-BLYP-D3, and XYG3 for GMTKN30 with a quadruple-ζ basis set. PWPB95-D3 and DSD-BLYP-D3 are the best functionals in our study and turned out to be more robust than B2PLYP-D3 and XYG3. Furthermore, PWPB95-D3 is the least basis set dependent and the best functional at the triple-ζ level. For the example of transition metal carbonyls, it is shown that, mainly due to the lower amount of Fock-exchange, PWPB95-D3 and PTPSS-D3 are better applicable than the other double hybrids. Finally, we discuss in some detail the XYG3 functional [Proc. Nat. Acad. Sci. U.S.A. 2009, 106, 4963], which makes use of B3LYP orbitals and electron densities. We show that it is basically a highly nonlocal variant of B2PLYP and that its partially good performance is mainly due to a larger effective amount of perturbative correlation compared to other double hybrids. We finally recommend the PWPB95-D3 functional in general chemistry applications.
Collapse
Affiliation(s)
- Lars Goerigk
- Theoretische Organische Chemie, Organisch-Chemisches Institut der Universität Münster, Corrensstraβe 40, and NRW Graduate School of Chemistry, Wilhelm-Klemm-Straβe 10, D-48149 Münster, Germany
| | - Stefan Grimme
- Theoretische Organische Chemie, Organisch-Chemisches Institut der Universität Münster, Corrensstraβe 40, and NRW Graduate School of Chemistry, Wilhelm-Klemm-Straβe 10, D-48149 Münster, Germany
| |
Collapse
|
29
|
Faßhauer E, Kryzhevoi NV, Pernpointner M. Possible electronic decay channels in the ionization spectra of small clusters composed of Ar and Xe: A four-component relativistic treatment. J Chem Phys 2010; 133:014303. [DOI: 10.1063/1.3462246] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Wander MC, Clark AE. Gradient fit functions for two-body potential energy surfaces based upon a harmonic series. MOLECULAR SIMULATION 2010. [DOI: 10.1080/08927020903433739] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
31
|
Devdariani A, Belyaev A, Alekseyev A, Liebermann HP, Buenker R. Ab initio study of the lifetime of weakly bound He(21S)–Ne molecules. Mol Phys 2010. [DOI: 10.1080/00268971003610226] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
32
|
Sato T, Nakai H. Density functional method including weak interactions: Dispersion coefficients based on the local response approximation. J Chem Phys 2010; 131:224104. [PMID: 20001021 DOI: 10.1063/1.3269802] [Citation(s) in RCA: 154] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
A new method to calculate the atom-atom dispersion coefficients in a molecule is proposed for the use in density functional theory with dispersion (DFT-D) correction. The method is based on the local response approximation due to Dobson and Dinte [Phys. Rev. Lett. 76, 1780 (1996)], with modified dielectric model recently proposed by Vydrov and van Voorhis [J. Chem. Phys. 130, 104105 (2009)]. The local response model is used to calculate the distributed multipole polarizabilities of atoms in a molecule, from which the dispersion coefficients are obtained by an explicit frequency integral of the Casimir-Polder type. Thus obtained atomic polarizabilities are also used in the damping function for the short-range singularity. Unlike empirical DFT-D methods, the local response dispersion (LRD) method is able to calculate the dispersion energy from the ground-state electron density only. It is applicable to any geometry, free from physical constants such as van der Waals radii or atomic polarizabilities, and computationally very efficient. The LRD method combined with the long-range corrected DFT functional (LC-BOP) is applied to calculations of S22 weakly bound complex set [Phys. Chem. Chem. Phys. 8, 1985 (2006)]. Binding energies obtained by the LC-BOP+LRD agree remarkably well with ab initio references.
Collapse
Affiliation(s)
- Takeshi Sato
- Research Institute for Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| | | |
Collapse
|
33
|
Yousaf KE, Brothers EN. Applications of Screened Hybrid Density Functionals with Empirical Dispersion Corrections to Rare Gas Dimers and Solids. J Chem Theory Comput 2010; 6:864-72. [DOI: 10.1021/ct900536n] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Kazim E. Yousaf
- Science Program, Texas A&M at Qatar, Texas A&M Engineering Building, Education City, Doha, Qatar
| | - Edward N. Brothers
- Science Program, Texas A&M at Qatar, Texas A&M Engineering Building, Education City, Doha, Qatar
| |
Collapse
|
34
|
Baranowska A, Capelo SB, Fernández B. New basis sets for the evaluation of interaction energies: an ab initio study of the He–He, Ne–Ne, Ar–Ar, He–Ne, He–Ar and Ne–Ar van der Waals complex internuclear potentials and ro-vibrational spectra. Phys Chem Chem Phys 2010; 12:13586-96. [DOI: 10.1039/c0cp00535e] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
35
|
Goll E, Werner HJ, Stoll H. A short-range gradient-corrected density functional in long-range coupled-cluster calculations for rare gas dimers. Phys Chem Chem Phys 2009; 7:3917-23. [PMID: 19810319 DOI: 10.1039/b509242f] [Citation(s) in RCA: 181] [Impact Index Per Article: 12.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A previously proposed scheme for coupling short-range (sr) density functionals with wavefunction-based long-range (lr) ab initio methods has been extended by (a) developing a new gradient-corrected sr functional of the Perdew-Burke-Ernzerhof (PBE) type and (b) introducing coupled-cluster (CC) approaches (CC with single and double excitations (CCSD), and with additional perturbative triples (CCSD(T))) at the ab initio side. The results show that mixing-in of lr-ab initio correlation helps to remove deficiencies of currently used density functionals for the treatment of van-der-Waals interactions. Compared to full ab initio calculations, the basis set dependence is weaker so that the accuracy of the mixed results surpasses that of the ab initio ones for basis sets of triple-zeta quality.
Collapse
Affiliation(s)
- Erich Goll
- Institut für Theoretische Chemie, Universität Stuttgart, D-70550 Stuttgart, Germany.
| | | | | |
Collapse
|
36
|
Pirani F, Brizi S, Roncaratti LF, Casavecchia P, Cappelletti D, Vecchiocattivi F. Beyond the Lennard-Jones model: a simple and accurate potential function probed by high resolution scattering data useful for molecular dynamics simulations. Phys Chem Chem Phys 2008; 10:5489-503. [PMID: 18956082 DOI: 10.1039/b808524b] [Citation(s) in RCA: 192] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Scattering data, measured for rare gas-rare gas systems under high angular and energy resolution conditions, have been used to probe the reliability of a recently proposed interaction potential function, which involves only one additional parameter with respect to the venerable Lennard-Jones (LJ) model and is hence called Improved Lennard-Jones (ILJ). The ILJ potential eliminates most of the inadequacies at short- and long-range of the LJ model. Further reliability tests have been performed by comparing calculated vibrational spacings with experimental values and calculated interaction energies at short-range with those obtained from the inversion of gaseous transport properties. The analysis, extended also to systems involving ions, suggests that the ILJ potential model can be used to estimate the behavior of unknown systems and can help to assess the different role of the leading interaction components. Moreover, due to its simple formulation, the physically reliable ILJ model appears to be particularly useful for molecular dynamics simulations of both neutral and ionic systems.
Collapse
Affiliation(s)
- Fernando Pirani
- Dipartimento di Chimica, Università di Perugia, 06123, Perugia, Italy.
| | | | | | | | | | | |
Collapse
|
37
|
Civalleri B, Middlemiss D, Orlando R, Wilson C, Ugliengo P. Testing the combination of Hartree–Fock exchange and Wilson–Levy correlation for weakly bonded extended systems. Chem Phys Lett 2008. [DOI: 10.1016/j.cplett.2007.11.095] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
38
|
Housden MP, Pyper NC. The noble gas dimers as a probe of the energetic contributions of dispersion and short-range electron correlation in weakly bound systems. Mol Phys 2007. [DOI: 10.1080/00268970701618432] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
39
|
Zhao Y, Truhlar DG. Comparative DFT study of van der Waals complexes: rare-gas dimers, alkaline-earth dimers, zinc dimer, and zinc-rare-gas dimers. J Phys Chem A 2007; 110:5121-9. [PMID: 16610834 DOI: 10.1021/jp060231d] [Citation(s) in RCA: 565] [Impact Index Per Article: 33.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Recent interest in the application of density functional theory prompted us to test various functionals for the van der Waals interactions in the rare-gas dimers, the alkaline-earth metal dimers, zinc dimer, and zinc-rare-gas dimers. In the present study, we report such tests for 18 DFT functionals, including both some very recent functionals and some well-established older ones. We draw the following conclusions based on the mean errors in binding energies and complex geometries: (1) B97-1 gives the best performance for predicting the geometry of rare-gas dimers, whereas M05-2X and B97-1 give the best energetics for rare-gas dimers. (2) PWB6K gives the best performance for the prediction of the geometry of the alkaline-earth metal dimers, zinc dimers, and zinc-rare-gas dimers. M05-2X gives the best energetics for the metal dimers, whereas B97-1 gives the best energetics for the zinc-rare-gas dimers. (3) The M05 functional is unique in providing good accuracy for both covalent transition-metal dimers and van der Waals metal dimers. (4) The combined mean percentage unsigned error in geometries and energetics shows that M05-2X and MPWB1K are the overall best methods for the prediction of van der Waals interactions in metal and rare-gas van der Waals dimers.
Collapse
Affiliation(s)
- Yan Zhao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | | |
Collapse
|
40
|
Sato T, Tsuneda T, Hirao K. Long-range corrected density functional study on weakly bound systems: Balanced descriptions of various types of molecular interactions. J Chem Phys 2007; 126:234114. [PMID: 17600411 DOI: 10.1063/1.2747243] [Citation(s) in RCA: 126] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
The long-range correction scheme for the density functional theory, combined with a van der Waals functional, is examined for its applicability to a wide variety of weakly bonded complexes including dispersion, dipole-induced dipole, dipole-dipole, and hydrogen-bonded complexes. The present method is able to describe geometries and binding energies of all these complexes accurately. Explicit inclusion of the long-range exchange and dispersion interactions is found to be important for the balanced description of various kinds of weak interactions. The present method is a promising alternative for high-level ab initio methods in calculations of large and complex systems, because it gives equally correct descriptions for various types of molecular interactions with much less computational cost.
Collapse
Affiliation(s)
- Takeshi Sato
- Department of Applied Chemistry, School of Engineering, University of Tokyo, Tokyo 113-8656, Japan.
| | | | | |
Collapse
|
41
|
Piris M, Lopez X, Ugalde JM. Dispersion interactions within the Piris natural orbital functional theory: The helium dimer. J Chem Phys 2007; 126:214103. [PMID: 17567186 DOI: 10.1063/1.2743019] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The authors have investigated the description of the dispersion interaction within the Piris natural orbital functional (PNOF) theory. The PNOF arises from an explicit antisymmetric approach for the two-particle cumulant in terms of two symmetric matrices, Delta and Lambda. The functional forms of these matrices are obtained from the generalization of the two-particle system expressions, except for the off-diagonal elements of Delta. The mean value theorem and the partial sum rule obtained for the off-diagonal elements of Delta provide a prescription for deriving practical functionals. In particular, the previous employed approximation {Jpp/2} for the mean values {Jp*} affords several molecular properties but it is incapable to account for dispersion effects. In this work, the authors analyze a new approach for Jp* obtained by factorization of the matrix Delta within the bounds on its off-diagonal elements imposed by the positivity conditions of the two-particle reduced density matrix. Additional terms for the matrix elements of Lambda proportional to the square root of the holes are again introduced to describe properly the occupation numbers of the lowest occupied levels. The authors have found that the cross products between weakly occupied orbitals must be removed from the functional form of Lambda to obtain a correct long-range asymptotic behavior. The PNOF is used to predict the binding energy as well as the equilibrium distance of the helium dimer. The results are compared with the full configuration-interaction calculations and the corresponding experimental data.
Collapse
Affiliation(s)
- M Piris
- Kimika Fakultatea, Euskal Herriko Unibertsitatea, Donostia International Physics Center (DIPC), P.K. 1072, 20080 Donostia, Euskadi, Spain
| | | | | |
Collapse
|
42
|
Gerber IC, Angyán JG. London dispersion forces by range-separated hybrid density functional with second order perturbational corrections: The case of rare gas complexes. J Chem Phys 2007; 126:044103. [PMID: 17286458 DOI: 10.1063/1.2431644] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
A satisfactory account of the van der Waals (vdW) (London dispersion) forces is, in general not possible by the Kohn-Sham method using standard local, semilocal generalized gradient approximation (GGA), or meta-GGA density functionals. The recently proposed range-separated hybrid (RSH) approach, supplemented by second order perturbational corrections (MP2) to include long-range dynamic correlation effects, offers a physically consistent, seamless description of dispersion forces. It is based on a rigorous generalization of the Kohn-Sham method, where long-range exchange and correlation effects are treated by wave function methods, while short-range electron exchange and correlation are handled by local or semilocal functionals. The method is tested on a series of rare gas dimers in comparison with standard wave function theory and density functional theory approaches. In contrast to the most successful exchange correlation functionals, which describe at best the vdW minimum, the RSH+MP2 approach is valid also in the asymptotic region and the potential curve displays the correct 1/R(6) behavior at large internuclear separations. In contrast to usual MP2 calculations, the basis set superposition error is considerably reduced, making RSH+MP2 an ideal tool for exploring the potential energy surface of weakly bound molecular complexes.
Collapse
Affiliation(s)
- I C Gerber
- Laboratoire de Cristallographie et de Modélisation des Matériaux Minéraux et Biologiques, UMR 7036, CNRS, Nancy Université, B.P. 239, F-54506 Vandoeuvre-lès-Nancy, France
| | | |
Collapse
|
43
|
Bouazza MT, Bouledroua M. Diffusion coefficient of krypton atoms in helium gas at low and moderate temperatures. Mol Phys 2007. [DOI: 10.1080/00268970601129050] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- M. T. Bouazza
- Physics Department, Badji Mokhtar University, B.P. 12, Annaba 23000, Algeria
| | - M. Bouledroua
- Physics Department, Badji Mokhtar University, B.P. 12, Annaba 23000, Algeria
| |
Collapse
|
44
|
Grimme S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J Comput Chem 2006; 27:1787-99. [PMID: 16955487 DOI: 10.1002/jcc.20495] [Citation(s) in RCA: 13637] [Impact Index Per Article: 757.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.
Collapse
Affiliation(s)
- Stefan Grimme
- Theoretische Organische Chemie, Organisch-Chemisches Institut der Universität Münster, Corrensstrasse 40, D-48149 Münster, Germany.
| |
Collapse
|
45
|
Mori-Sánchez P, Wu Q, Yang W. Orbital-dependent correlation energy in density-functional theory based on a second-order perturbation approach: Success and failure. J Chem Phys 2005; 123:62204. [PMID: 16122290 DOI: 10.1063/1.1904584] [Citation(s) in RCA: 96] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have developed a second-order perturbation theory (PT) energy functional within density-functional theory (DFT). Based on PT with the Kohn-Sham (KS) determinant as a reference, this new ab initio exchange-correlation functional includes an exact exchange (EXX) energy in the first order and a correlation energy including all single and double excitations from the KS reference in the second order. The explicit dependence of the exchange and correlation energy on the KS orbitals in the functional fits well into our direct minimization approach for the optimized effective potential, which is a very efficient method to perform fully self-consistent calculations for any orbital-dependent functionals. To investigate the quality of the correlation functional, we have applied the method to selected atoms and molecules. For two-electron atoms and small molecules described with small basis sets, this new method provides excellent results, improving both second-order Moller-Plesset expression and any conventional DFT results significantly. For larger systems, however, it performs poorly, converging to very low unphysical total energies. The failure of PT based energy functionals is analyzed, and its origin is traced back to near degeneracy problems due to the orbital- and eigenvalue-dependent algebraic structure of the correlation functional. The failure emerges in the self-consistent approach but not in perturbative post-EXX calculations, emphasizing the crucial importance of self-consistency in testing new orbital-dependent energy functionals.
Collapse
Affiliation(s)
- Paula Mori-Sánchez
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | | | | |
Collapse
|
46
|
Grimme S. Accurate description of van der Waals complexes by density functional theory including empirical corrections. J Comput Chem 2005; 25:1463-73. [PMID: 15224390 DOI: 10.1002/jcc.20078] [Citation(s) in RCA: 2994] [Impact Index Per Article: 157.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
An empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT-D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R(-6). The use of pure, gradient-corrected density functionals (BLYP and PBE), together with the resolution-of-the-identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene. H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies.
Collapse
Affiliation(s)
- Stefan Grimme
- Organisch-Chemisches Institut der Universität Münster, Corrensstrasse 40, D-48149 Münster, Germany.
| |
Collapse
|
47
|
Furche F, Van Voorhis T. Fluctuation-dissipation theorem density-functional theory. J Chem Phys 2005; 122:164106. [PMID: 15945671 DOI: 10.1063/1.1884112] [Citation(s) in RCA: 141] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Using the fluctuation-dissipation theorem (FDT) in the context of density-functional theory (DFT), one can derive an exact expression for the ground-state correlation energy in terms of the frequency-dependent density response function. When combined with time-dependent density-functional theory, a new class of density functionals results that use approximations to the exchange-correlation kernel fxc as input. This FDT-DFT scheme holds promise to solve two of the most distressing problems of conventional Kohn-Sham DFT: (i) It leads to correlation energy functionals compatible with exact exchange, and (ii) it naturally includes dispersion. The price is a moderately expensive O(N6) scaling of computational cost and a slower basis set convergence. These general features of FDT-DFT have all been recognized previously. In this paper, we present the first benchmark results for a set of molecules using FDT-DFT beyond the random-phase approximation (RPA)-that is, the first such results with fxc not equal to 0. We show that kernels derived from the adiabatic local-density approximation and other semilocal functionals suffer from an "ultraviolet catastrophe," producing a pair density that diverges at small interparticle distance. Nevertheless, dispersion interactions can be treated accurately if hybrid functionals are employed, as is demonstrated for He2 and HeNe. We outline constraints that future approximations to fxc should satisfy and discuss the prospects of FDT-DFT.
Collapse
Affiliation(s)
- Filipp Furche
- Institut für Physikalische Chemie, Universität Karlsruhe, Kaiserstrasse 12, Karlsruhe 76128, Germany
| | | |
Collapse
|
48
|
|
49
|
Tao J, Perdew JP. Test of a nonempirical density functional: Short-range part of the van der Waals interaction in rare-gas dimers. J Chem Phys 2005; 122:114102. [PMID: 15836196 DOI: 10.1063/1.1862242] [Citation(s) in RCA: 100] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It is known that the nonempirical generalized gradient approximation (GGA) of Perdew, Burke, and Ernzerhof (PBE) provides a much more realistic description of the short-range part of the van der Waals (vdW) interaction than does the local spin density (LSD) approximation. In the present work, the ability of the higher-level nonempirical meta-GGA of Tao, Perdew, Staroverov, and Scuseria (TPSS) to describe vdW interaction is tested self-consistently in ten rare-gas dimers with Z< or =36. The one-parameter hybrid version (TPSSh) of the TPSS exchange-correlation functional is also included in this test. Calculations show that both TPSS and TPSSh functionals correctly yield vdW bonds in these dimers and significantly improve the prediction of bond lengths, binding energies, and harmonic vibrational frequencies over LSD. The rather close agreement of TPSS with PBE for these dimers confirms a principle of the TPSS construction: preservation of the PBE large-gradient behavior. More importantly, it suggests that TPSS can serve as a platform on which to construct a still-higher level of nonempirical functionals. Compared with the PBE GGA, TPSS, and TPSSh yield a slightly weaker binding. As for normally bonded molecules, TPSSh yields the most accurate vibrational frequencies. The typically too-long bond lengths and too-small binding energies of TPSS meta-GGA suggest the need for some long-range vdW interaction correction even in this class of systems. The effect of basis-set superposition error on the calculated properties of these vdW systems is investigated. We also show that the relatively strong anharmonic effects in the rare-gas dimers are described remarkably well by the Morse potential.
Collapse
Affiliation(s)
- Jianmin Tao
- Department of Physics and Quantum Theory Group, Tulane University, New Orleans, Louisiana 70118, USA
| | | |
Collapse
|
50
|
Seletskaia T, Osetsky Y, Stoller RE, Stocks GM. Magnetic interactions influence the properties of helium defects in iron. PHYSICAL REVIEW LETTERS 2005; 94:046403. [PMID: 15783579 DOI: 10.1103/physrevlett.94.046403] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/12/2004] [Indexed: 05/24/2023]
Abstract
Density functional theory calculations of He defect properties in iron have shown an unexpected influence of magnetism arising from the defect's electronic structure. In contrast with previous work that neglected such effects, the results indicate that the tetrahedral position is energetically more favorable for the He interstitial than the octahedral site. This may have significant implications for He clustering and bubble nucleation, which will impact material performance in future fusion reactors. These results provide the basis for development of improved atomistic models.
Collapse
Affiliation(s)
- Tatiana Seletskaia
- Metals and Ceramics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6138, USA
| | | | | | | |
Collapse
|