1
|
Pham KN, Modrzejewski M, Klimeš J. Contributions beyond direct random-phase approximation in the binding energy of solid ethane, ethylene, and acetylene. J Chem Phys 2024; 160:224101. [PMID: 38856055 DOI: 10.1063/5.0207090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2024] [Accepted: 05/22/2024] [Indexed: 06/11/2024] Open
Abstract
The random-phase approximation (RPA) includes a subset of higher than second-order correlation-energy contributions, but stays in the same complexity class as the second-order Møller-Plesset perturbation theory (MP2) in both Gaussian-orbital and plane-wave codes. This makes RPA a promising ab initio electronic structure approach for the binding energies of molecular crystals. Still, some issues stand out in practical applications of RPA. Notably, compact clusters of nonpolar molecules are poorly described, and the interaction energies strongly depend on the reference single-determinant state. Using the many-body expansion of the binding energy of a crystal, we investigate those issues and the effect of beyond-RPA corrections. We find the beneficial effect of quartic-scaling exchange and non-ring coupled-cluster doubles corrections. The nonadditive interactions in compact trimers of molecules are improved by using the self-consistent Hartree-Fock orbitals instead of the usual Kohn-Sham states, but this kind of orbital input also leads to underestimated dimer energies. Overall, a substantial improvement over the RPA with a renormalized singles approach is possible at a modest quartic-scaling cost, which encourages further research into additional RPA corrections.
Collapse
Affiliation(s)
- Khanh Ngoc Pham
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
2
|
Joshi P, Voora VK. Generalized perturbative singles corrections to the random phase approximation method: Impact on noncovalent interaction energies of closed- and open-shell dimers. J Chem Phys 2024; 160:044104. [PMID: 38258929 DOI: 10.1063/5.0180526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Accepted: 12/29/2023] [Indexed: 01/24/2024] Open
Abstract
The post-Kohn-Sham (KS) random phase approximation (RPA) method may provide a poor description of interaction energies of weakly bonded molecules due to inherent density errors in approximate KS functionals. To overcome these errors, we develop a generalized formalism to incorporate perturbative singles (pS) corrections to the RPA method using orbital rotations as a perturbation parameter. The pS schemes differ in the choice of orbital-rotation gradient and Hessian. We propose a pS scheme termed RPA singles (RPAS)[Hartree-Fock (HF)] that uses the RPA orbital-rotation gradient and time-dependent HF Hessian. This correction reduces the errors in noncovalent interaction energies of closed- and open-shell dimers. For the open-shell dimers, the RPAS(HF) method leads to a consistent error reduction by 50% or more compared to the RPA method for the cases of hydrogen-bonding, metal-solvent, carbene-solvent, and dispersion interactions. We also find that the pS corrections are more important in error reduction compared to higher-order exchange corrections to the RPA method. Overall, for open shells, the RPAS(HF)-corrected RPA method provides chemical accuracy for noncovalent interactions and is more reliable than other perturbative schemes and dispersion-corrected density functional approximations, highlighting its importance as a reliable beyond-RPA correction.
Collapse
Affiliation(s)
- Pulkit Joshi
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| |
Collapse
|
3
|
Yang S, Zhang IY, Ren X. Developing correlation-consistent numeric atom-centered orbital basis sets for krypton: Applications in RPA-based correlated calculations. J Chem Phys 2024; 160:024112. [PMID: 38193553 DOI: 10.1063/5.0174952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 12/04/2023] [Indexed: 01/10/2024] Open
Abstract
Localized atomic orbitals are the preferred basis set choice for large-scale explicit correlated calculations, and high-quality hierarchical correlation-consistent basis sets are a prerequisite for correlated methods to deliver numerically reliable results. At present, numeric atom-centered orbital (NAO) basis sets with valence correlation consistency (VCC), designated as NAO-VCC-nZ, are only available for light elements from hydrogen (H) to argon (Ar) [Zhang et al., New J. Phys. 15, 123033 (2013)]. In this work, we extend this series by developing NAO-VCC-nZ basis sets for krypton (Kr), a prototypical element in the fourth row of the periodic table. We demonstrate that NAO-VCC-nZ basis sets facilitate the convergence of electronic total-energy calculations using the Random Phase Approximation (RPA), which can be used together with a two-point extrapolation scheme to approach the complete basis set limit. Notably, the Basis Set Superposition Error (BSSE) associated with the newly generated NAO basis sets is minimal, making them suitable for applications where BSSE correction is either cumbersome or impractical to do. After confirming the reliability of NAO basis sets for Kr, we proceed to calculate the Helmholtz free energy for Kr crystal at the theoretical level of RPA plus renormalized single excitation correction. From this, we derive the pressure-volume (P-V) diagram, which shows excellent agreement with the latest experimental data. Our work demonstrates the capability of correlation-consistent NAO basis sets for heavy elements, paving the way toward numerically reliable correlated calculations for bulk materials.
Collapse
Affiliation(s)
- Sixian Yang
- Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, China
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Igor Ying Zhang
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai, Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xinguo Ren
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
4
|
Cieśliński D, Tucholska AM, Modrzejewski M. Post-Kohn-Sham Random-Phase Approximation and Correction Terms in the Expectation-Value Coupled-Cluster Formulation. J Chem Theory Comput 2023; 19:6619-6631. [PMID: 37774375 PMCID: PMC10569055 DOI: 10.1021/acs.jctc.3c00496] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Indexed: 10/01/2023]
Abstract
Using expectation-value coupled-cluster theory and many-body perturbation theory (MBPT), we formulate a series of corrections to the post-Kohn-Sham (post-KS) random-phase approximation (RPA) energy. The beyond-RPA terms are of two types: those accounting for the non-Hartree-Fock reference and those introducing the coupled-cluster doubles non-ring contractions. The contributions of the former type, introduced via the semicanonical orbital basis, drastically reduce the binding strength in noncovalent systems. The good accuracy is recovered by the attractive third-order doubles correction referred to as Ec2g. The existing RPA approaches based on KS orbitals neglect most of the proposed corrections but can perform well thanks to error cancellation. The proposed method accounts for every contribution in the state-of-the-art renormalized second-order perturbation theory (rPT2) approach but adds additional terms which initially contribute in the third order of MBPT. The cost of energy evaluation scales as noniterative O ( N 4 ) in the implementation with low-rank tensor decomposition. The numerical tests of the proposed approach demonstrate accurate results for noncovalent dimers of polar molecules and for the challenging many-body noncovalent cluster of CH4···(H2O)20.
Collapse
Affiliation(s)
- Dominik Cieśliński
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | | | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
5
|
Ruan S, Jackson KA, Ruzsinszky A. Spin-crossover complexes: Self-interaction correction vs density correction. J Chem Phys 2023; 158:064303. [PMID: 36792493 DOI: 10.1063/5.0128950] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Complexes containing a transition metal atom with a 3d4-3d7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi-Löwdin self-interaction correction. Considering that Hartree-Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree-Fock method and the Perdew-Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. We find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.
Collapse
Affiliation(s)
- Shiqi Ruan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Koblar A Jackson
- Physics Department and Science of Advanced Materials Ph. D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, USA
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
6
|
Schäfer T, Daelman N, López N. Cerium Oxides without U: The Role of Many-Electron Correlation. J Phys Chem Lett 2021; 12:6277-6283. [PMID: 34212726 PMCID: PMC8397342 DOI: 10.1021/acs.jpclett.1c01589] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 06/29/2021] [Indexed: 05/14/2023]
Abstract
Electron transfer with changing occupation in the 4f subshell poses a considerable challenge for quantitative predictions in quantum chemistry. Using the example of cerium oxide, we identify the main deficiencies of common parameter-dependent one-electron approaches, such as density functional theory (DFT) with a Hubbard correction, or hybrid functionals. As a response, we present the first benchmark of ab initio many-electron theory for electron transfer energies and lattice parameters under periodic boundary conditions. We show that the direct random phase approximation clearly outperforms all DFT variations. From this foundation, we, then, systematically improve even further. Periodic second-order Møller-Plesset perturbation theory meanwhile manages to recover standard hybrid functional values. Using these approaches to eliminate parameter bias allows for highly accurate benchmarks of strongly correlated materials, the reliable assessment of various density functionals, and functional fitting via machine-learning.
Collapse
Affiliation(s)
- Tobias Schäfer
- Institute
for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Nathan Daelman
- Institute
of Chemical Research of Catalonia, The Barcelona
Institute of Science and Technology, 43007 Tarragona, Spain
| | - Núria López
- Institute
of Chemical Research of Catalonia, The Barcelona
Institute of Science and Technology, 43007 Tarragona, Spain
| |
Collapse
|
7
|
Modrzejewski M, Yourdkhani S, Śmiga S, Klimeš J. Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage. J Chem Theory Comput 2021; 17:804-817. [PMID: 33445879 DOI: 10.1021/acs.jctc.0c00966] [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/22/2022]
Abstract
The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH4···(H2O)20, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.
Collapse
Affiliation(s)
- Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland.,Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Sirous Yourdkhani
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziądzka 5, 87-100 Toruń, Poland
| | - Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
8
|
Riemelmoser S, Kaltak M, Kresse G. Plane wave basis set correction methods for RPA correlation energies. J Chem Phys 2020; 152:134103. [DOI: 10.1063/5.0002246] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Stefan Riemelmoser
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Merzuk Kaltak
- VASP Software GmbH, Sensengasse 8/17, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria
| |
Collapse
|
9
|
Klimeš J, Tew DP. Efficient and accurate description of adsorption in zeolites. J Chem Phys 2019; 151:234108. [PMID: 31864262 DOI: 10.1063/1.5123425] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.
Collapse
Affiliation(s)
- Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - David P Tew
- Max-Planck-Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
10
|
Modrzejewski M, Yourdkhani S, Klimeš J. Random Phase Approximation Applied to Many-Body Noncovalent Systems. J Chem Theory Comput 2019; 16:427-442. [DOI: 10.1021/acs.jctc.9b00979] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcin Modrzejewski
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Sirous Yourdkhani
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
11
|
Patra A, Jana S, Myneni H, Samal P. Laplacian free and asymptotic corrected semilocal exchange potential applied to the band gap of solids. Phys Chem Chem Phys 2019; 21:19639-19650. [PMID: 31468039 DOI: 10.1039/c9cp03356d] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
It is well known that the modified semilocal exchange potentials explicitly designed for the study of solid-state band gaps are very successful in describing these properties. These exchange potentials are in principle designed either from a spherically averaged exchange hole or by satisfying the exact asymptotic conditions. In this present attempt, we use the recently developed novel technique of density matrix expansion to construct the model exchange hole potential. The proposed exchange hole potential is free from the Laplacian of density and generalized through the coordinate transformation. An improvement in the exchange energies of atoms using this potential is shown. The salient feature of the proposed semilocal potential is that it can be used within generalized Kohn-Sham formalism because of its Laplacian free representation. This modified potential is used in the framework of TBMBJ [Phys. Rev. Lett., 2009, 102, 226401] to calculate the band gaps of materials. The comparison and assessment of the newly constructed Laplacian free, asymptotically corrected semilocal potential to address the band gap problem show good agreement with the experimental band gaps and diversify the studies done in the same direction.
Collapse
Affiliation(s)
- Abhilash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India.
| | - Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India.
| | - Hemanadhan Myneni
- Department of Physics, Temple University, Philadelphia, PA 19122-1801, USA
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India.
| |
Collapse
|
12
|
Hummel F, Grüneis A, Kresse G, Ziesche P. Screened Exchange Corrections to the Random Phase Approximation from Many-Body Perturbation Theory. J Chem Theory Comput 2019; 15:3223-3236. [PMID: 30901204 DOI: 10.1021/acs.jctc.8b01247] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The random phase approximation (RPA) systematically overestimates the magnitude of the correlation energy and generally underestimates cohesive energies. This originates in part from the complete lack of exchange terms that would otherwise cancel Pauli exclusion principle violating (EPV) contributions. The uncanceled EPV contributions also manifest themselves in form of an unphysical negative pair density of spin-parallel electrons close to electron-electron coalescence. We follow considerations of many-body perturbation theory to propose an exchange correction that corrects the largest set of EPV contributions, while having the lowest possible computational complexity. The proposed method exchanges adjacent particle/hole pairs in the RPA diagrams, considerably improving the pair density of spin-parallel electrons close to coalescence in the uniform electron gas (UEG). The accuracy of the correlation energy is comparable to other variants of second-order screened exchange (SOSEX) corrections although it is slightly more accurate for the spin-polarized UEG. Its computational complexity scales as [Formula: see text] or [Formula: see text] in orbital space or real space, respectively. Its memory requirement scales as [Formula: see text].
Collapse
Affiliation(s)
- Felix Hummel
- Institute for Theoretical Physics , Technische Universität Wien , Wiedner Hauptstraße 8-10/136 , 1040 Vienna , Austria
| | - Andreas Grüneis
- Institute for Theoretical Physics , Technische Universität Wien , Wiedner Hauptstraße 8-10/136 , 1040 Vienna , Austria
| | - Georg Kresse
- University of Vienna , Faculty of Physics and Center for Computational Materials Sciences , Sensengasse 8/12 , 1090 Vienna , Austria
| | - Paul Ziesche
- Max Planck Institute for the Physics of Complex Systems Nöthnitzer Straße 38 , 01187 Dresden , Germany
| |
Collapse
|
13
|
Ajala AO, Voora V, Mardirossian N, Furche F, Paesani F. Assessment of Density Functional Theory in Predicting Interaction Energies between Water and Polycyclic Aromatic Hydrocarbons: from Water on Benzene to Water on Graphene. J Chem Theory Comput 2019; 15:2359-2374. [PMID: 30860827 DOI: 10.1021/acs.jctc.9b00110] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The interactions of water with polycyclic aromatic hydrocarbons, from benzene to graphene, are investigated using various exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. The accuracy of the different functionals is assessed through comparisons with random phase approximation (RPA) and coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] calculations. Diffusion Monte Carlo (DMC) data reported in the literature are also used for comparison. Relatively large variations are found in interaction energies predicted by different DFT models, with GGA functionals underestimating the interaction strength for configurations with the water oxygen pointing toward the aromatic molecules. The meta-GGA B97M-rV and range-separated hybrid, meta-GGA ωB97M-V functionals provide nearly quantitative agreement with CCSD(T) values for the water-benzene, water-coronene, and water-circumcoronene dimers, while RPA and DMC predict interaction energies that differ by up to ∼1 kcal/mol and ∼0.4 kcal/mol from the corresponding CCSD(T) values, respectively. Similar trends among GGA, meta-GGA, and hybrid functionals are observed for larger polycyclic aromatic hydrocarbons. By performing absolutely localized molecular orbital energy decomposition analyses (ALMO-EDA), it is found that, independently of the number of carbon atoms and exchange-correlation functional, the dominant contributions to the interaction energies between water and polycyclic aromatic hydrocarbon molecules are the electrostatic and dispersion terms while polarization and charge transfer effects are negligibly small. Calculations carried out with GGA and meta-GGA functionals indicate that, as the number of carbon atoms increases, the interaction energies slowly converge to the corresponding values obtained for an infinite graphene sheet.
Collapse
Affiliation(s)
- Adeayo O Ajala
- Department of Chemistry and Biochemistry , University of California San Diego , La Jolla , California 92093 , United States
| | - Vamsee Voora
- Department of Chemistry , University of California Irvine , Irvine , California 92697 , United States
| | - Narbe Mardirossian
- Division of Chemistry and Chemical Engineering , California Institute of Technology , 1200 E. California Boulevard , Pasadena , California 91125 , United States
| | - Filipp Furche
- Department of Chemistry , University of California Irvine , Irvine , California 92697 , United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry , University of California San Diego , La Jolla , California 92093 , United States.,Materials Science and Engineering , University of California San Diego , La Jolla , California 92093 , United States.,San Diego Supercomputer Center , University of California San Diego , La Jolla , California 92093 , United States
| |
Collapse
|
14
|
Brandenburg JG, Zen A, Fitzner M, Ramberger B, Kresse G, Tsatsoulis T, Grüneis A, Michaelides A, Alfè D. Physisorption of Water on Graphene: Subchemical Accuracy from Many-Body Electronic Structure Methods. J Phys Chem Lett 2019; 10:358-368. [PMID: 30615460 DOI: 10.1021/acs.jpclett.8b03679] [Citation(s) in RCA: 64] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Wet carbon interfaces are ubiquitous in the natural world and exhibit anomalous properties, which could be exploited by emerging technologies. However, progress is limited by lack of understanding at the molecular level. Remarkably, even for the most fundamental system (a single water molecule interacting with graphene), there is no consensus on the nature of the interaction. We tackle this by performing an extensive set of complementary state-of-the-art computer simulations on some of the world's largest supercomputers. From this effort a consensus on the water-graphene interaction strength has been obtained. Our results have significant impact for the physical understanding, as they indicate that the interaction is weaker than predicted previously. They also pave the way for more accurate and reliable studies of liquid water at carbon interfaces.
Collapse
Affiliation(s)
- Jan Gerit Brandenburg
- Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom
| | - Andrea Zen
- Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom
| | - Martin Fitzner
- Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom
| | - Benjamin Ramberger
- University of Vienna , Faculty of Physics and Center for Computational Materials Sciences , Sensengasse 8/12 , 1090 Wien , Austria
| | - Georg Kresse
- University of Vienna , Faculty of Physics and Center for Computational Materials Sciences , Sensengasse 8/12 , 1090 Wien , Austria
| | - Theodoros Tsatsoulis
- Institute for Theoretical Physics , Vienna University of Technology , Wiedner Hauptstrasse 8-10 , 1040 Vienna , Austria
- Max Planck Institute for Solid State Research , Heisenbergstrasse 1 , 70569 Stuttgart , Germany
| | - Andreas Grüneis
- Institute for Theoretical Physics , Vienna University of Technology , Wiedner Hauptstrasse 8-10 , 1040 Vienna , Austria
- Max Planck Institute for Solid State Research , Heisenbergstrasse 1 , 70569 Stuttgart , Germany
| | - Angelos Michaelides
- Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom
- Department of Earth Sciences , University College London , Gower Street , London WC1E 6BT , United Kingdom
- Dipartimento di Fisica Ettore Pancini , Università di Napoli Federico II , Monte S. Angelo, I-80126 Napoli , Italy
| |
Collapse
|
15
|
Jin Y, Su NQ, Yang W. Renormalized Singles Green's Function for Quasi-Particle Calculations beyond the G 0 W 0 Approximation. J Phys Chem Lett 2019; 10:447-452. [PMID: 30609900 DOI: 10.1021/acs.jpclett.8b03337] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Quasi-particle energies and band gaps in particular are critical for investigating novel materials. Commonly used density functional approximations (DFAs) systematically underestimate band gaps, and GW approximation is the established method of choice for good accuracy and reliability. However, G0 W0 has some undesired dependence on the DFA, while self-consistent GW (sc GW) is expensive and not consistent in accuracy improvement. Here a simple and efficient GRS W0 approach has been developed: a subspace diagonalization of the Hartree-Fock (HF) Hamiltonian with the DFA density matrix provides the new reference Green's function GRS that incorporates the effect of all single excitation contributions to the self-energy, thereby essentially eliminating the starting-point dependence. Calculations for molecules and large band gap solids demonstrate the significant improvement over G0 W0 and greatly reduced dependence on the initial DFA. GRS W0 approach also improve results for other bulks over G0 W0, but to a lesser extent, which could be due to the limitations in current implementation for bulks. The results demonstrate that to achieve good accuracy, it is not necessary to use hybrid DFA, which is expensive for bulks. This work should be greatly significant in making GW a more robust approach.
Collapse
Affiliation(s)
- Ye Jin
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States
| | - Neil Qiang Su
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States
| | - Weitao Yang
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment , South China Normal University , Guangzhou 510006 , China
| |
Collapse
|
16
|
Stöhr M, Van Voorhis T, Tkatchenko A. Theory and practice of modeling van der Waals interactions in electronic-structure calculations. Chem Soc Rev 2019; 48:4118-4154. [PMID: 31190037 DOI: 10.1039/c9cs00060g] [Citation(s) in RCA: 68] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.
Collapse
Affiliation(s)
- Martin Stöhr
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg.
| | | | | |
Collapse
|
17
|
Tuning the balance between dispersion and entropy to design temperature-responsive flexible metal-organic frameworks. Nat Commun 2018; 9:4899. [PMID: 30464249 PMCID: PMC6249296 DOI: 10.1038/s41467-018-07298-4] [Citation(s) in RCA: 65] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2018] [Accepted: 10/25/2018] [Indexed: 11/09/2022] Open
Abstract
Temperature-responsive flexibility in metal-organic frameworks (MOFs) appeals to the imagination. The ability to transform upon thermal stimuli while retaining a given crystalline topology is desired for specialized sensors and actuators. However, rational design of such shape-memory nanopores is hampered by a lack of knowledge on the nanoscopic interactions governing the observed behavior. Using the prototypical MIL-53(Al) as a starting point, we show that the phase transformation between a narrow-pore and large-pore phase is determined by a delicate balance between dispersion stabilization at low temperatures and entropic effects at higher ones. We present an accurate theoretical framework that allows designing breathing thermo-responsive MOFs, based on many-electron data for the dispersion interactions and density-functional theory entropy contributions. Within an isoreticular series of materials, MIL-53(Al), MIL-53(Al)-FA, DUT-4, DUT-5 and MIL-53(Ga), only MIL-53(Al) and MIL-53(Ga) are proven to switch phases within a realistic temperature range. Rational design of metal organic frameworks (MOFs) with shape-memory nanopores is a formidable challenge. Here the authors use an accurate theoretical approach to design thermo-responsive MOFs based on a balance of van der Waals and entropy contributions.
Collapse
|
18
|
Jin Y, Zhang D, Chen Z, Su NQ, Yang W. Generalized Optimized Effective Potential for Orbital Functionals and Self-Consistent Calculation of Random Phase Approximations. J Phys Chem Lett 2017; 8:4746-4751. [PMID: 28895734 PMCID: PMC6209318 DOI: 10.1021/acs.jpclett.7b02165] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
A new self-consistent procedure for calculating the total energy with an orbital-dependent density functional approximation (DFA), the generalized optimized effective potential (GOEP), is developed in the present work. The GOEP is a nonlocal Hermitian potential that delivers the sets of occupied and virtual orbitals and minimizes the total energy. The GOEP optimization leads to the same minimum as does the orbital optimization. The GOEP method is promising as an effective optimization approach for orbital-dependent functionals, as demonstrated for the self-consistent calculations of the random phase approximation (RPA) to the correlation functionals in the particle-hole (ph) and particle-particle (pp) channels. The results show that the accuracy in describing the weakly interacting van der Waals systems is significantly improved in the self-consistent calculations. In particular, the important single excitations contribution in non-self-consistent RPA calculations can be captured self-consistently through the GOEP optimization, leading to orbital renormalization, without using the single excitations in the energy functional.
Collapse
Affiliation(s)
- Ye Jin
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Du Zhang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Zehua Chen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Neil Qiang Su
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
19
|
Al-Hamdani YS, Rossi M, Alfè D, Tsatsoulis T, Ramberger B, Brandenburg JG, Zen A, Kresse G, Grüneis A, Tkatchenko A, Michaelides A. Properties of the water to boron nitride interaction: From zero to two dimensions with benchmark accuracy. J Chem Phys 2017; 147:044710. [DOI: 10.1063/1.4985878] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yasmine S. Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Mariana Rossi
- Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, 14195 Berlin, Germany
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Theodoros Tsatsoulis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Wien, Austria
| | - Jan Gerit Brandenburg
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AH, United Kingdom
| | - Andrea Zen
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Wien, Austria
| | - Andreas Grüneis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
20
|
Zhang D, Su NQ, Yang W. Accurate Quasiparticle Spectra from the T-Matrix Self-Energy and the Particle-Particle Random Phase Approximation. J Phys Chem Lett 2017; 8:3223-3227. [PMID: 28654275 DOI: 10.1021/acs.jpclett.7b01275] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The GW self-energy, especially G0W0 based on the particle-hole random phase approximation (phRPA), is widely used to study quasiparticle (QP) energies. Motivated by the desirable features of the particle-particle (pp) RPA compared to the conventional phRPA, we explore the pp counterpart of GW, that is, the T-matrix self-energy, formulated with the eigenvectors and eigenvalues of the ppRPA matrix. We demonstrate the accuracy of the T-matrix method for molecular QP energies, highlighting the importance of the pp channel for calculating QP spectra.
Collapse
Affiliation(s)
- Du Zhang
- Department of Chemistry, Duke University , Durham, North Carolina 27708, United States
| | - Neil Qiang Su
- Department of Chemistry, Duke University , Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University , Durham, North Carolina 27708, United States
- Key Laboratory of Theoretical Chemistry of Environment, School of Chemistry and Environment, South China Normal University , Guangzhou 510006, China
| |
Collapse
|
21
|
Ramberger B, Schäfer T, Kresse G. Analytic Interatomic Forces in the Random Phase Approximation. PHYSICAL REVIEW LETTERS 2017; 118:106403. [PMID: 28339256 DOI: 10.1103/physrevlett.118.106403] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/04/2016] [Indexed: 06/06/2023]
Abstract
We discuss that in the random phase approximation (RPA) the first derivative of the energy with respect to the Green's function is the self-energy in the GW approximation. This relationship allows us to derive compact equations for the RPA interatomic forces. We also show that position dependent overlap operators are elegantly incorporated in the present framework. The RPA force equations have been implemented in the projector augmented wave formalism, and we present illustrative applications, including ab initio molecular dynamics simulations, the calculation of phonon dispersion relations for diamond and graphite, as well as structural relaxations for water on boron nitride. The present derivation establishes a concise framework for forces within perturbative approaches and is also applicable to more involved approximations for the correlation energy.
Collapse
Affiliation(s)
- Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Vienna, Austria
| | - Tobias Schäfer
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Vienna, Austria
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Vienna, Austria
| |
Collapse
|
22
|
Hermann J, DiStasio RA, Tkatchenko A. First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications. Chem Rev 2017; 117:4714-4758. [PMID: 28272886 DOI: 10.1021/acs.chemrev.6b00446] [Citation(s) in RCA: 255] [Impact Index Per Article: 36.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.
Collapse
Affiliation(s)
- Jan Hermann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University , Ithaca, New York 14853, United States
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany.,Physics and Materials Science Research Unit, University of Luxembourg , L-1511 Luxembourg, Luxembourg
| |
Collapse
|
23
|
Masur O, Schütz M, Maschio L, Usvyat D. Fragment-Based Direct-Local-Ring-Coupled-Cluster Doubles Treatment Embedded in the Periodic Hartree–Fock Solution. J Chem Theory Comput 2016; 12:5145-5156. [DOI: 10.1021/acs.jctc.6b00651] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Oliver Masur
- Institut
für Physikalische und Theoretische Chemie, Universität Regensburg, 93049 Regensburg, Germany
| | - Martin Schütz
- Institut
für Physikalische und Theoretische Chemie, Universität Regensburg, 93049 Regensburg, Germany
| | - Lorenzo Maschio
- Dipartimento
di Chimica, and NIS (Nanostructured Interfaces and Surfaces) Centre, Università di Torino, via Giuria 5, I-10125 Torino, Italy
| | - Denis Usvyat
- Institut
für Physikalische und Theoretische Chemie, Universität Regensburg, 93049 Regensburg, Germany
| |
Collapse
|
24
|
Klimeš J. Lattice energies of molecular solids from the random phase approximation with singles corrections. J Chem Phys 2016; 145:094506. [DOI: 10.1063/1.4962188] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jiří Klimeš
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-18223 Prague 8, Czech Republic and Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|