1
|
Tahir MN, Shang H, Li J, Ren X. Efficient Structural Relaxation Based on the Random Phase Approximation: Applications to Water Clusters. J Phys Chem A 2024. [PMID: 39240284 DOI: 10.1021/acs.jpca.4c02411] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
We report an improved implementation for evaluating the analytical gradients of the random phase approximation (RPA) electron-correlation energy based on atomic orbitals and the localized resolution of the identity scheme. The more efficient RPA force calculations allow us to relax the structures of medium-sized water clusters. Particular attention is paid to the structures and energy orderings of the low-energy isomers of (H2O)n clusters with n = 21, 22, and 25. It is found that the RPA energy ordering of the low-energy isomers of these water clusters is rather sensitive to the basis set used. For the five low-energy isomers of (H2O)25, the RPA energy ordering still undergoes a change by increasing the basis set to the quadruple to quintuple level. The standard RPA underbinds the water clusters, and this underbinding behavior becomes more pronounced by increasing the basis size to the complete basis set (CBS) limit. The renormalized single excitation (rSE) correction remedies this underbinding, giving rise to a noticeable overbinding behavior at finite basis sets. However, as the CBS limit is approached, RPA+rSE yields an accuracy for the binding energies that is comparable to that of the best available double hybrid functionals, as demonstrated for the WATER27 test set.
Collapse
Affiliation(s)
- Muhammad N Tahir
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Honghui Shang
- University of Science and Technology of China, Hefei 230026, China
| | - Jia Li
- Shenzhen Geim Graphene Center and Institute of Materials Research, Tsinghua Shenzhen International Graduate School, Tsinghua University, Shenzhen 518055, China
| | - Xinguo Ren
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
2
|
Kokott S, Merz F, Yao Y, Carbogno C, Rossi M, Havu V, Rampp M, Scheffler M, Blum V. Efficient all-electron hybrid density functionals for atomistic simulations beyond 10 000 atoms. J Chem Phys 2024; 161:024112. [PMID: 38990115 DOI: 10.1063/5.0208103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2024] [Accepted: 06/19/2024] [Indexed: 07/12/2024] Open
Abstract
Hybrid density functional approximations (DFAs) offer compelling accuracy for ab initio electronic-structure simulations of molecules, nanosystems, and bulk materials, addressing some deficiencies of computationally cheaper, frequently used semilocal DFAs. However, the computational bottleneck of hybrid DFAs is the evaluation of the non-local exact exchange contribution, which is the limiting factor for the application of the method for large-scale simulations. In this work, we present a drastically optimized resolution-of-identity-based real-space implementation of the exact exchange evaluation for both non-periodic and periodic boundary conditions in the all-electron code FHI-aims, targeting high-performance central processing unit (CPU) compute clusters. The introduction of several new refined message passing interface (MPI) parallelization layers and shared memory arrays according to the MPI-3 standard were the key components of the optimization. We demonstrate significant improvements of memory and performance efficiency, scalability, and workload distribution, extending the reach of hybrid DFAs to simulation sizes beyond ten thousand atoms. In addition, we also compare the runtime performance of the PBE, HSE06, and PBE0 functionals. As a necessary byproduct of this work, other code parts in FHI-aims have been optimized as well, e.g., the computation of the Hartree potential and the evaluation of the force and stress components. We benchmark the performance and scaling of the hybrid DFA-based simulations for a broad range of chemical systems, including hybrid organic-inorganic perovskites, organic crystals, and ice crystals with up to 30 576 atoms (101 920 electrons described by 244 608 basis functions).
Collapse
Affiliation(s)
- Sebastian Kokott
- The NOMAD Laboratory at the Fritz Haber Institute of the Max-Planck-Gesellschaft and IRIS Adlershof of the Humboldt-Universität zu Berlin, Berlin, Germany
| | - Florian Merz
- Lenovo HPC Innovation Center, Stuttgart, Germany
| | - Yi Yao
- Thomas Lord Department of Mechanical Engineering and Material Science, Duke University, Durham, North Carolina 27708, USA
| | - Christian Carbogno
- The NOMAD Laboratory at the Fritz Haber Institute of the Max-Planck-Gesellschaft and IRIS Adlershof of the Humboldt-Universität zu Berlin, Berlin, Germany
| | - Mariana Rossi
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Ville Havu
- Department of Applied Physics, School of Science, Aalto University, Espoo, Finland
| | - Markus Rampp
- Max Planck Computing and Data Facility, 85748 Garching, Germany
| | - Matthias Scheffler
- The NOMAD Laboratory at the Fritz Haber Institute of the Max-Planck-Gesellschaft and IRIS Adlershof of the Humboldt-Universität zu Berlin, Berlin, Germany
| | - Volker Blum
- Thomas Lord Department of Mechanical Engineering and Material Science, Duke University, Durham, North Carolina 27708, USA
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
3
|
Riemelmoser S, Verdi C, Kaltak M, Kresse G. Machine Learning Density Functionals from the Random-Phase Approximation. J Chem Theory Comput 2023; 19:7287-7299. [PMID: 37800677 PMCID: PMC10601474 DOI: 10.1021/acs.jctc.3c00848] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Indexed: 10/07/2023]
Abstract
Kohn-Sham density functional theory (DFT) is the standard method for first-principles calculations in computational chemistry and materials science. More accurate theories such as the random-phase approximation (RPA) are limited in application due to their large computational cost. Here, we use machine learning to map the RPA to a pure Kohn-Sham density functional. The machine learned RPA model (ML-RPA) is a nonlocal extension of the standard gradient approximation. The density descriptors used as ingredients for the enhancement factor are nonlocal counterparts of the local density and its gradient. Rather than fitting only RPA exchange-correlation energies, we also include derivative information in the form of RPA optimized effective potentials. We train a single ML-RPA functional for diamond, its surfaces, and liquid water. The accuracy of ML-RPA for the formation energies of 28 diamond surfaces reaches that of state-of-the-art van der Waals functionals. For liquid water, however, ML-RPA cannot yet improve upon the standard gradient approximation. Overall, our work demonstrates how machine learning can extend the applicability of the RPA to larger system sizes, time scales, and chemical spaces.
Collapse
Affiliation(s)
- Stefan Riemelmoser
- Faculty
of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, A-1090 Vienna, Austria
| | - Carla Verdi
- Faculty
of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
- School
of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia
- School
of Mathematics and Physics, The University
of Queensland, Brisbane, Queensland 4072, Australia
| | - Merzuk Kaltak
- VASP
Software GmbH, Sensengasse
8/12, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty
of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
- VASP
Software GmbH, Sensengasse
8/12, A-1090 Vienna, Austria
| |
Collapse
|
4
|
Engel EA. Identification of synthesisable crystalline phases of water – a prototype for the challenges of computational materials design. CrystEngComm 2021. [DOI: 10.1039/d0ce01260b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We discuss the identification of experimentally realisable crystalline phases of water to outline and contextualise some of the diverse building blocks of a computational materials design process.
Collapse
Affiliation(s)
- Edgar A. Engel
- TCM Group
- Cavendish Laboratory
- University of Cambridge
- Cambridge CB3 0HE
- UK
| |
Collapse
|
5
|
Jana S, Patra A, Śmiga S, Constantin LA, Samal P. Insights from the density functional performance of water and water–solid interactions: SCAN in relation to other meta-GGAs. J Chem Phys 2020; 153:214116. [DOI: 10.1063/5.0028821] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Abhilash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland
| | - Lucian A. Constantin
- Istituto di Nanoscienze, Consiglio Nazionale delle Ricerche CNR-NANO, 41125 Modena, Italy
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
6
|
Greenwell C, McKinley JL, Zhang P, Zeng Q, Sun G, Li B, Wen S, Beran GJO. Overcoming the difficulties of predicting conformational polymorph energetics in molecular crystals via correlated wavefunction methods. Chem Sci 2020; 11:2200-2214. [PMID: 32190277 PMCID: PMC7059316 DOI: 10.1039/c9sc05689k] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2019] [Accepted: 01/13/2020] [Indexed: 11/21/2022] Open
Abstract
Molecular crystal structure prediction is increasingly being applied to study the solid form landscapes of larger, more flexible pharmaceutical molecules. Despite many successes in crystal structure prediction, van der Waals-inclusive density functional theory (DFT) methods exhibit serious failures predicting the polymorph stabilities for a number of systems exhibiting conformational polymorphism, where changes in intramolecular conformation lead to different intermolecular crystal packings. Here, the stabilities of the conformational polymorphs of o-acetamidobenzamide, ROY, and oxalyl dihydrazide are examined in detail. DFT functionals that have previously been very successful in crystal structure prediction perform poorly in all three systems, due primarily to the poor intramolecular conformational energies, but also due to the intermolecular description in oxalyl dihydrazide. In all three cases, a fragment-based dispersion-corrected second-order Møller-Plesset perturbation theory (MP2D) treatment of the crystals overcomes these difficulties and predicts conformational polymorph stabilities in good agreement with experiment. These results highlight the need for methods which go beyond current-generation DFT functionals to make crystal polymorph stability predictions truly reliable.
Collapse
Affiliation(s)
- Chandler Greenwell
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| | - Jessica L McKinley
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| | - Peiyu Zhang
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Qun Zeng
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Guangxu Sun
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Bochen Li
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Shuhao Wen
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Gregory J O Beran
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| |
Collapse
|
7
|
Jana S, Constantin LA, Samal P. Accurate Water Properties from an Efficient ab Initio Method. J Chem Theory Comput 2020; 16:974-987. [DOI: 10.1021/acs.jctc.9b01018] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Lucian A. Constantin
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
8
|
Sakti AW, Nishimura Y, Nakai H. Recent advances in quantum‐mechanical molecular dynamics simulations of proton transfer mechanism in various water‐based environments. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1419] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Aditya W. Sakti
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
| | - Hiromi Nakai
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering Waseda University Tokyo Japan
| |
Collapse
|
9
|
Caldeweyher E, Mewes JM, Ehlert S, Grimme S. Extension and evaluation of the D4 London-dispersion model for periodic systems. Phys Chem Chem Phys 2020; 22:8499-8512. [PMID: 32292979 DOI: 10.1039/d0cp00502a] [Citation(s) in RCA: 117] [Impact Index Per Article: 29.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present an extension of the DFT-D4 model [J. Chem. Phys., 2019, 150, 154122] for periodic systems. The main new ingredients are additional reference polarizabilities for highly-coordinated group 1-5 elements derived from pseudo-periodic electrostatically-embedded cluster calculations. To illustrate the performance of the updated method, several test cases are considered, for which we compare D4 to its predecessor D3(BJ), as well as to a comprehensive set of other dispersion-corrected methods. The largest improvements are observed for solid-state polarizabilities of 16 inorganic salts, where the D4 model achieves an unprecedented accuracy, surpassing its predecessor as well as other, computationally much more demanding approaches. For cell volumes and lattice energies of two sets of chemically diverse molecular crystals, the accuracy gain is less pronounced compared to the already excellently performing D3(BJ) method. For the challenging adsorption energies of small organic molecules on metallic as well as on ionic surfaces, DFT-D4 provides values in good agreement with experimental and/or high-level references. These results suggest the application of the proposed D4 model as a physically improved yet computationally efficient dispersion correction for standard DFT calculations as well as low-cost approaches like semi-empirical or even force-field models.
Collapse
Affiliation(s)
| | | | | | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Bonn, Germany.
| |
Collapse
|
10
|
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
|
11
|
Cheng B, Engel EA, Behler J, Dellago C, Ceriotti M. Ab initio thermodynamics of liquid and solid water. Proc Natl Acad Sci U S A 2019; 116:1110-1115. [PMID: 30610171 PMCID: PMC6347673 DOI: 10.1073/pnas.1815117116] [Citation(s) in RCA: 152] [Impact Index Per Article: 30.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Thermodynamic properties of liquid water as well as hexagonal (Ih) and cubic (Ic) ice are predicted based on density functional theory at the hybrid-functional level, rigorously taking into account quantum nuclear motion, anharmonic fluctuations, and proton disorder. This is made possible by combining advanced free-energy methods and state-of-the-art machine-learning techniques. The ab initio description leads to structural properties in excellent agreement with experiments and reliable estimates of the melting points of light and heavy water. We observe that nuclear-quantum effects contribute a crucial [Formula: see text] to the stability of ice Ih, making it more stable than ice Ic. Our computational approach is general and transferable, providing a comprehensive framework for quantitative predictions of ab initio thermodynamic properties using machine-learning potentials as an intermediate step.
Collapse
Affiliation(s)
- Bingqing Cheng
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland;
| | - Edgar A Engel
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Universität Göttingen, 37073 Göttingen, Germany
| | | | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
12
|
Goel H, Ling S, Ellis BN, Taconi A, Slater B, Rai N. Predicting vapor liquid equilibria using density functional theory: A case study of argon. J Chem Phys 2018; 148:224501. [PMID: 29907054 DOI: 10.1063/1.5025726] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase.
Collapse
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Sanliang Ling
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Breanna Nicole Ellis
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Anna Taconi
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| |
Collapse
|
13
|
Bokdam M, Lahnsteiner J, Ramberger B, Schäfer T, Kresse G. Assessing Density Functionals Using Many Body Theory for Hybrid Perovskites. PHYSICAL REVIEW LETTERS 2017; 119:145501. [PMID: 29053325 DOI: 10.1103/physrevlett.119.145501] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2017] [Indexed: 05/18/2023]
Abstract
Which density functional is the "best" for structure simulations of a particular material? A concise, first principles, approach to answer this question is presented. The random phase approximation (RPA)-an accurate many body theory-is used to evaluate various density functionals. To demonstrate and verify the method, we apply it to the hybrid perovskite MAPbI_{3}, a promising new solar cell material. The evaluation is done by first creating finite temperature ensembles for small supercells using RPA molecular dynamics, and then evaluating the variance between the RPA and various approximate density functionals for these ensembles. We find that, contrary to recent suggestions, van der Waals functionals do not improve the description of the material, whereas hybrid functionals and the strongly constrained appropriately normed (SCAN) density functional yield very good agreement with the RPA. Finally, our study shows that in the room temperature tetragonal phase of MAPbI_{3}, the molecules are preferentially parallel to the shorter lattice vectors but reorientation on ps time scales is still possible.
Collapse
Affiliation(s)
- Menno Bokdam
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Jonathan Lahnsteiner
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Tobias Schäfer
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, A-1090 Vienna, Austria
| |
Collapse
|
14
|
Pham CH, Reddy SK, Chen K, Knight C, Paesani F. Many-Body Interactions in Ice. J Chem Theory Comput 2017; 13:1778-1784. [DOI: 10.1021/acs.jctc.6b01248] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- C. Huy Pham
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Sandeep K. Reddy
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Karl Chen
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Chris Knight
- Leadership
Computing Facility, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, United States
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| |
Collapse
|
15
|
Salzmann CG, Slater B, Radaelli PG, Finney JL, Shephard JJ, Rosillo-Lopez M, Hindley J. Detailed crystallographic analysis of the ice VI to ice XV hydrogen ordering phase transition. J Chem Phys 2016; 145:204501. [DOI: 10.1063/1.4967167] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Christoph G. Salzmann
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Paolo G. Radaelli
- Department of Physics, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom
| | - John L. Finney
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Jacob J. Shephard
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Martin Rosillo-Lopez
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - James Hindley
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
16
|
Wilhelm J, Seewald P, Del Ben M, Hutter J. Large-Scale Cubic-Scaling Random Phase Approximation Correlation Energy Calculations Using a Gaussian Basis. J Chem Theory Comput 2016; 12:5851-5859. [DOI: 10.1021/acs.jctc.6b00840] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Wilhelm
- Department
of Chemistry and National Centre for Computational Design and Discovery
of Novel Materials (MARVEL), University of Zurich, 8057 Zurich, Switzerland
| | - Patrick Seewald
- Department
of Chemistry and National Centre for Computational Design and Discovery
of Novel Materials (MARVEL), University of Zurich, 8057 Zurich, Switzerland
| | - Mauro Del Ben
- Computational
Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jürg Hutter
- Department
of Chemistry and National Centre for Computational Design and Discovery
of Novel Materials (MARVEL), University of Zurich, 8057 Zurich, Switzerland
| |
Collapse
|
17
|
Liu Y, Ojamäe L. Raman and IR Spectra of Ice Ih and Ice XI with an Assessment of DFT Methods. J Phys Chem B 2016; 120:11043-11051. [PMID: 27690444 DOI: 10.1021/acs.jpcb.6b07001] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
IR and Raman spectroscopic technology can be directly used to identify the occurrence of ferroelectric ice XI in laboratory or extraterrestrial settings. The performance of 16 different DFT methods applied on the ice Ih, VIII, IX, and XI crystal phases is evaluated. Based on a selected DFT computational scheme, the IR and Raman spectra of ice Ih and XI are derived and compared. When the spectra, both IR and Raman, of ice Ih and ice XI are compared, the librational vibrations are found to be the most affected by the proton ordering. The spectroscopic fingerprint of ice XI can be used to distinguish ferroelectric ice XI from ice Ih in the universe. Furthermore, the existence of only one kind of H-bond in ice Ih is demonstrated from the overlapping subspectra for different types of H-bonded pair configurations in 16 isomers of ice Ih, which provides an illustration to the historic debate on whether one or two kinds of H-bonds existed in ice.
Collapse
Affiliation(s)
- Yuan Liu
- Department of Chemistry, IFM, Linköping University , SE-581 83 Linköping, Sweden
| | - Lars Ojamäe
- Department of Chemistry, IFM, Linköping University , SE-581 83 Linköping, Sweden
| |
Collapse
|
18
|
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
|
19
|
Komatsu K, Noritake F, Machida S, Sano-Furukawa A, Hattori T, Yamane R, Kagi H. Partially ordered state of ice XV. Sci Rep 2016; 6:28920. [PMID: 27375120 PMCID: PMC4931510 DOI: 10.1038/srep28920] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2016] [Accepted: 06/13/2016] [Indexed: 11/24/2022] Open
Abstract
Most ice polymorphs have order-disorder "pairs" in terms of hydrogen positions, which contributes to the rich variety of ice polymorphs; in fact, three recently discovered polymorphs- ices XIII, XIV, and XV-are ordered counter forms to already identified disordered phases. Despite the considerable effort to understand order-disorder transition in ice crystals, there is an inconsistency among the various experiments and calculations for ice XV, the ordered counter form of ice VI, i.e., neutron diffraction observations suggest antiferroelectrically ordered structures, which disagree with dielectric measurement and theoretical studies, implying ferroelectrically ordered structures. Here we investigate in-situ neutron diffraction measurements and density functional theory calculations to revisit the structure and stability of ice XV. We find that none of the completely ordered configurations are particular favored; instead, partially ordered states are established as a mixture of ordered domains in disordered ice VI. This scenario in which several kinds of ordered configuration coexist dispels the contradictions in previous studies. It means that the order-disorder pairs in ice polymorphs are not one-to-one correspondent pairs but rather have one-to-n correspondence, where there are n possible configurations at finite temperature.
Collapse
Affiliation(s)
- K. Komatsu
- Geochemical Research Center, Graduate School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
| | - F. Noritake
- Geochemical Research Center, Graduate School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
| | - S. Machida
- CROSS-Tokai, Research Center for Neutron Science and Technology, IQBRC Bldg, 162-1 Shirakata, Tokai, Ibaraki 319-1106, Japan
| | - A. Sano-Furukawa
- J-PARC Center, Japan Atomic Energy Agency, 2-4 Shirakata-Shirane, Tokai, Ibaraki 319-1195, Japan
| | - T. Hattori
- J-PARC Center, Japan Atomic Energy Agency, 2-4 Shirakata-Shirane, Tokai, Ibaraki 319-1195, Japan
| | - R. Yamane
- Geochemical Research Center, Graduate School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
| | - H. Kagi
- Geochemical Research Center, Graduate School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
| |
Collapse
|
20
|
Cutini M, Civalleri B, Corno M, Orlando R, Brandenburg JG, Maschio L, Ugliengo P. Assessment of Different Quantum Mechanical Methods for the Prediction of Structure and Cohesive Energy of Molecular Crystals. J Chem Theory Comput 2016; 12:3340-52. [DOI: 10.1021/acs.jctc.6b00304] [Citation(s) in RCA: 66] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michele Cutini
- Department
of Chemistry and NIS (Nanostructured Interfaces and Surfaces) Center, University of Turin, Via P. Giuria 7, 10125 Turin, Italy
| | - Bartolomeo Civalleri
- Department
of Chemistry and NIS (Nanostructured Interfaces and Surfaces) Center, University of Turin, Via P. Giuria 7, 10125 Turin, Italy
| | - Marta Corno
- Department
of Chemistry and NIS (Nanostructured Interfaces and Surfaces) Center, University of Turin, Via P. Giuria 7, 10125 Turin, Italy
| | - Roberto Orlando
- Department
of Chemistry and NIS (Nanostructured Interfaces and Surfaces) Center, University of Turin, Via P. Giuria 7, 10125 Turin, Italy
| | - Jan Gerit Brandenburg
- Mulliken
Center of Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstraße
4, 53115 Bonn, Germany
| | - Lorenzo Maschio
- Department
of Chemistry and NIS (Nanostructured Interfaces and Surfaces) Center, University of Turin, Via P. Giuria 7, 10125 Turin, Italy
| | - Piero Ugliengo
- Department
of Chemistry and NIS (Nanostructured Interfaces and Surfaces) Center, University of Turin, Via P. Giuria 7, 10125 Turin, Italy
| |
Collapse
|
21
|
Salim MA, Willow SY, Hirata S. Ice Ih anomalies: Thermal contraction, anomalous volume isotope effect, and pressure-induced amorphization. J Chem Phys 2016; 144:204503. [DOI: 10.1063/1.4951687] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael A. Salim
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - Soohaeng Yoo Willow
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| |
Collapse
|
22
|
Paier J. Hybrid Density Functionals Applied to Complex Solid Catalysts: Successes, Limitations, and Prospects. Catal Letters 2016. [DOI: 10.1007/s10562-016-1735-4] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
23
|
Abstract
Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
24
|
Brandenburg JG, Maas T, Grimme S. Benchmarking DFT and semiempirical methods on structures and lattice energies for ten ice polymorphs. J Chem Phys 2016; 142:124104. [PMID: 25833562 DOI: 10.1063/1.4916070] [Citation(s) in RCA: 71] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water in different phases under various external conditions is very important in bio-chemical systems and for material science at surfaces. Density functional theory methods and approximations thereof have to be tested system specifically to benchmark their accuracy regarding computed structures and interaction energies. In this study, we present and test a set of ten ice polymorphs in comparison to experimental data with mass densities ranging from 0.9 to 1.5 g/cm(3) and including explicit corrections for zero-point vibrational and thermal effects. London dispersion inclusive density functionals at the generalized gradient approximation (GGA), meta-GGA, and hybrid level as well as alternative low-cost molecular orbital methods are considered. The widely used functional of Perdew, Burke and Ernzerhof (PBE) systematically overbinds and overall provides inconsistent results. All other tested methods yield reasonable to very good accuracy. BLYP-D3(atm) gives excellent results with mean absolute errors for the lattice energy below 1 kcal/mol (7% relative deviation). The corresponding optimized structures are very accurate with mean absolute relative deviations (MARDs) from the reference unit cell volume below 1%. The impact of Axilrod-Teller-Muto (atm) type three-body dispersion and of non-local Fock exchange is small but on average their inclusion improves the results. While the density functional tight-binding model DFTB3-D3 performs well for low density phases, it does not yield good high density structures. As low-cost alternative for structure related problems, we recommend the recently introduced minimal basis Hartree-Fock method HF-3c with a MARD of about 3%.
Collapse
Affiliation(s)
- Jan Gerit Brandenburg
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Tilo Maas
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
25
|
Del Ben M, Hutter J, VandeVondele J. Probing the structural and dynamical properties of liquid water with models including non-local electron correlation. J Chem Phys 2015; 143:054506. [PMID: 26254660 DOI: 10.1063/1.4927325] [Citation(s) in RCA: 75] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaric-isothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance.
Collapse
Affiliation(s)
- Mauro Del Ben
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Joost VandeVondele
- Department of Materials, ETH Zurich, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland
| |
Collapse
|
26
|
Grimme S, Brandenburg JG, Bannwarth C, Hansen A. Consistent structures and interactions by density functional theory with small atomic orbital basis sets. J Chem Phys 2015; 143:054107. [PMID: 26254642 DOI: 10.1063/1.4927476] [Citation(s) in RCA: 531] [Impact Index Per Article: 59.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
A density functional theory (DFT) based composite electronic structure approach is proposed to efficiently compute structures and interaction energies in large chemical systems. It is based on the well-known and numerically robust Perdew-Burke-Ernzerhoff (PBE) generalized-gradient-approximation in a modified global hybrid functional with a relatively large amount of non-local Fock-exchange. The orbitals are expanded in Ahlrichs-type valence-double zeta atomic orbital (AO) Gaussian basis sets, which are available for many elements. In order to correct for the basis set superposition error (BSSE) and to account for the important long-range London dispersion effects, our well-established atom-pairwise potentials are used. In the design of the new method, particular attention has been paid to an accurate description of structural parameters in various covalent and non-covalent bonding situations as well as in periodic systems. Together with the recently proposed three-fold corrected (3c) Hartree-Fock method, the new composite scheme (termed PBEh-3c) represents the next member in a hierarchy of "low-cost" electronic structure approaches. They are mainly free of BSSE and account for most interactions in a physically sound and asymptotically correct manner. PBEh-3c yields good results for thermochemical properties in the huge GMTKN30 energy database. Furthermore, the method shows excellent performance for non-covalent interaction energies in small and large complexes. For evaluating its performance on equilibrium structures, a new compilation of standard test sets is suggested. These consist of small (light) molecules, partially flexible, medium-sized organic molecules, molecules comprising heavy main group elements, larger systems with long bonds, 3d-transition metal systems, non-covalently bound complexes (S22 and S66×8 sets), and peptide conformations. For these sets, overall deviations from accurate reference data are smaller than for various other tested DFT methods and reach that of triple-zeta AO basis set second-order perturbation theory (MP2/TZ) level at a tiny fraction of computational effort. Periodic calculations conducted for molecular crystals to test structures (including cell volumes) and sublimation enthalpies indicate very good accuracy competitive to computationally more involved plane-wave based calculations. PBEh-3c can be applied routinely to several hundreds of atoms on a single processor and it is suggested as a robust "high-speed" computational tool in theoretical chemistry and physics.
Collapse
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Jan Gerit Brandenburg
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
27
|
Klimeš J, Kaltak M, Maggio E, Kresse G. Singles correlation energy contributions in solids. J Chem Phys 2015; 143:102816. [PMID: 26374009 DOI: 10.1063/1.4929346] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The random phase approximation to the correlation energy often yields highly accurate results for condensed matter systems. However, ways how to improve its accuracy are being sought and here we explore the relevance of singles contributions for prototypical solid state systems. We set out with a derivation of the random phase approximation using the adiabatic connection and fluctuation dissipation theorem, but contrary to the most commonly used derivation, the density is allowed to vary along the coupling constant integral. This yields results closely paralleling standard perturbation theory. We re-derive the standard singles of Görling-Levy perturbation theory [A. Görling and M. Levy, Phys. Rev. A 50, 196 (1994)], highlight the analogy of our expression to the renormalized singles introduced by Ren and coworkers [Phys. Rev. Lett. 106, 153003 (2011)], and introduce a new approximation for the singles using the density matrix in the random phase approximation. We discuss the physical relevance and importance of singles alongside illustrative examples of simple weakly bonded systems, including rare gas solids (Ne, Ar, Xe), ice, adsorption of water on NaCl, and solid benzene. The effect of singles on covalently and metallically bonded systems is also discussed.
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
| | - Merzuk Kaltak
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Emanuele Maggio
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, 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
|
28
|
Mezei PD, Csonka GI, Ruzsinszky A, Kállay M. Construction and application of a new dual-hybrid random phase approximation. J Chem Theory Comput 2015; 11:4615-26. [PMID: 26574252 DOI: 10.1021/acs.jctc.5b00420] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The direct random phase approximation (dRPA) combined with Kohn-Sham reference orbitals is among the most promising tools in computational chemistry and applicable in many areas of chemistry and physics. The reason for this is that it scales as N(4) with the system size, which is a considerable advantage over the accurate ab initio wave function methods like standard coupled-cluster. dRPA also yields a considerably more accurate description of thermodynamic and electronic properties than standard density-functional theory methods. It is also able to describe strong static electron correlation effects even in large systems with a small or vanishing band gap missed by common single-reference methods. However, dRPA has several flaws due to its self-correlation error. In order to obtain accurate and precise reaction energies, barriers and noncovalent intra- and intermolecular interactions, we construct a new dual-hybrid dRPA (hybridization of exact and semilocal exchange in both the energy and the orbitals) and test the performance of this new functional on isogyric, isodesmic, hypohomodesmotic, homodesmotic, and hyperhomodesmotic reaction classes. We also use a test set of 14 Diels-Alder reactions, six atomization energies (AE6), 38 hydrocarbon atomization energies, and 100 reaction barrier heights (DBH24, HT-BH38, and NHT-BH38). For noncovalent complexes, we use the NCCE31 and S22 test sets. To test the intramolecular interactions, we use a set of alkane, cysteine, phenylalanine-glycine-glycine tripeptide, and monosaccharide conformers. We also discuss the delocalization and static correlation errors. We show that a universally accurate description of chemical properties can be provided by a large, 75% exact exchange mixing both in the calculation of the reference orbitals and the final energy.
Collapse
Affiliation(s)
- Pál D Mezei
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics , H-1521 Budapest, Hungary
| | - Gábor I Csonka
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics , H-1521 Budapest, Hungary
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University , Philadelphia, Pennsylvania 19122, United States
| | - Mihály Kállay
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics , P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
29
|
Li Y, Lousada CM, Soroka IL, Korzhavyi PA. Bond Network Topology and Antiferroelectric Order in Cuprice CuOH. Inorg Chem 2015; 54:8969-77. [DOI: 10.1021/acs.inorgchem.5b01030] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yunguo Li
- Department of Materials Science and Engineering and ‡Applied
Physical Chemistry, School of Chemical Science and Engineering, Royal Institute of Technology (KTH), S-100 44 Stockholm, Sweden
| | - Cláudio M. Lousada
- Department of Materials Science and Engineering and ‡Applied
Physical Chemistry, School of Chemical Science and Engineering, Royal Institute of Technology (KTH), S-100 44 Stockholm, Sweden
| | - Inna L. Soroka
- Department of Materials Science and Engineering and ‡Applied
Physical Chemistry, School of Chemical Science and Engineering, Royal Institute of Technology (KTH), S-100 44 Stockholm, Sweden
| | - Pavel A. Korzhavyi
- Department of Materials Science and Engineering and ‡Applied
Physical Chemistry, School of Chemical Science and Engineering, Royal Institute of Technology (KTH), S-100 44 Stockholm, Sweden
| |
Collapse
|
30
|
Bankura A, Santra B, DiStasio RA, Swartz CW, Klein ML, Wu X. A systematic study of chloride ion solvation in water using van der Waals inclusive hybrid density functional theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1059959] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Arindam Bankura
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | | | - Charles W. Swartz
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Michael L. Klein
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Xifan Wu
- Department of Physics, Temple University , Philadelphia, PA, USA
| |
Collapse
|
31
|
DiStasio RA, Santra B, Li Z, Wu X, Car R. The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water. J Chem Phys 2015; 141:084502. [PMID: 25173016 DOI: 10.1063/1.4893377] [Citation(s) in RCA: 186] [Impact Index Per Article: 20.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, POOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.
Collapse
Affiliation(s)
- Robert A DiStasio
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Zhaofeng Li
- Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
32
|
Fang C, Li WF, Koster RS, Klimeš J, van Blaaderen A, van Huis MA. The accurate calculation of the band gap of liquid water by means of GW corrections applied to plane-wave density functional theory molecular dynamics simulations. Phys Chem Chem Phys 2015; 17:365-75. [DOI: 10.1039/c4cp04202f] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Knowledge about the intrinsic electronic properties of water is imperative for understanding the behaviour of aqueous solutions that are used throughout biology, chemistry, physics, and industry.
Collapse
Affiliation(s)
- Changming Fang
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Wun-Fan Li
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Rik S. Koster
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Jiří Klimeš
- University of Vienna
- Faculty of Physics and Center for Computational Materials Science
- A-1090 Vienna
- Austria
| | - Alfons van Blaaderen
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Marijn A. van Huis
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| |
Collapse
|
33
|
Del Ben M, VandeVondele J, Slater B. Periodic MP2, RPA, and Boundary Condition Assessment of Hydrogen Ordering in Ice XV. J Phys Chem Lett 2014; 5:4122-8. [PMID: 26278943 DOI: 10.1021/jz501985w] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Ice XV is the hydrogen-ordered form of the ice VI phase whose structure was predicted to be Cc and ferroelectric using periodic DFT approaches. However, neutron diffraction and Raman spectroscopy data show the structure to have P1̅ symmetry and to be antiferroelectric. Recent work1 using fragment-based MP2 and CCSD(T) approaches predicts the experimental structure as the ground state. We have reconsidered this problem using fully periodic MP2 and RPA approaches and find that the ferroelectric Cc structure is the lowest energy configuration. However, ubiquitously employed tinfoil boundary conditions stabilize polar structures. We suggest that ferroelectric Cc crystals can grow within a paraelectric ice VI matrix but may become unstable once a fraction of the matrix has become hydrogen-ordered. The reduction in dielectric constant causes P1̅ and other structures with small polarization to become favored, providing a possible resolution between observation and theoretical predictions.
Collapse
Affiliation(s)
- Mauro Del Ben
- †Department of Chemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Joost VandeVondele
- ‡Department of Materials, ETH Zürich, Wolfgang-Pauli-Strasse 27, CH-8093 Zürich, Switzerland
| | - Ben Slater
- §Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
34
|
Slater B, Quigley D. Crystal nucleation: Zeroing in on ice. NATURE MATERIALS 2014; 13:670-671. [PMID: 24947779 DOI: 10.1038/nmat4017] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Affiliation(s)
- Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London W C1H 0AJ, UK
| | - David Quigley
- Department of Physics and Centre for Scientific Computing, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK
| |
Collapse
|
35
|
Geiger P, Dellago C, Macher M, Franchini C, Kresse G, Bernard J, Stern J, Loerting T. Proton Ordering of Cubic Ice Ic: Spectroscopy and Computer Simulations. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2014; 118:10989-10997. [PMID: 24883169 PMCID: PMC4032183 DOI: 10.1021/jp500324x] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/11/2014] [Revised: 04/29/2014] [Indexed: 06/03/2023]
Abstract
Several proton-disordered crystalline ice structures are known to proton order at sufficiently low temperatures, provided that the right preparation procedure is used. For cubic ice, ice Ic, however, no proton ordering has been observed so far. Here, we subject ice Ic to an experimental protocol similar to that used to proton order hexagonal ice. In situ FT-IR spectroscopy carried out during this procedure reveals that the librational band of the spectrum narrows and acquires a structure that is observed neither in proton-disordered ice Ic nor in ice XI, the proton-ordered variant of hexagonal ice. On the basis of vibrational spectra computed for ice Ic and four of its proton-ordered variants using classical molecular dynamics and ab initio simulations, we conclude that the features of our experimental spectra are due to partial proton ordering, providing the first evidence of proton ordering in cubic ice. We further find that the proton-ordered structure with the lowest energy is ferroelectric, while the structure with the second lowest energy is weakly ferroelectric. Both structures fit the experimental spectral similarly well such that no unique assignment of proton order is possible based on our results.
Collapse
Affiliation(s)
- Philipp Geiger
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Markus Macher
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Cesare Franchini
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Jürgen Bernard
- Institute
of Physical Chemistry, University of Innsbruck, Innrain 52a, 6020 Innsbruck, Austria
| | - Josef
N. Stern
- Institute
of Physical Chemistry, University of Innsbruck, Innrain 52a, 6020 Innsbruck, Austria
| | - Thomas Loerting
- Institute
of Physical Chemistry, University of Innsbruck, Innrain 52a, 6020 Innsbruck, Austria
| |
Collapse
|