1
|
Francisco H, Cancio AC, Trickey SB. Reworking the Tao-Mo Exchange-Correlation Functional. III. Improved Deorbitalization Strategy and Faithful Deorbitalization. J Phys Chem A 2024; 128:6010-6018. [PMID: 38994657 DOI: 10.1021/acs.jpca.4c02635] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
We present a deorbitalization of the recent simplified, regularized Tao-Mo exchange functional ( J. Chem. Phys. 2023, 159, 214102) that is faithful to the parent functional. That is a major gain relative to our earlier deorbitalization which did poorly on molecular heats of formation ( J. Chem. Phys. 2023, 159, 214103). The improvement arises from augmentation of the Mejía-Rodríguez and Trickey deorbitalization strategy ( Phys. Rev. A 2017, 96, 052512) to use a smoothed replacement for the reduced density Laplacian (conventionally denoted q) obtained from that Laplacian itself. The augmentation also rationalizes the improvement obtained from the cutoff of q < 0 that was poorly understood at the time of the previous paper. The new scheme yields deorbitalized chemical region indicators that are much closer to those from the parent, orbital-dependent functional than were obtainable from the previous deorbitalization. It also replicates the good 3d elemental magnetization of the parent functional reasonably well.
Collapse
Affiliation(s)
- H Francisco
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, United States
| | - A C Cancio
- Department of Physics and Astronomy, Ball State University, Muncie, Indiana 47306, United States
| | - S B Trickey
- Quantum Theory Project, Department of Physics and Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
2
|
Fischer M, Brauer J. Studying the adsorption of emerging organic contaminants in zeolites with dispersion-corrected density functional theory calculations: From numbers to recommendations. ChemistryOpen 2024; 13:e202300273. [PMID: 38385822 PMCID: PMC11230941 DOI: 10.1002/open.202300273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Revised: 01/10/2024] [Indexed: 02/23/2024] Open
Abstract
Adsorption energies obtained from dispersion-corrected density functional theory (DFT) calculations show a considerable dependence on the choice of exchange-correlation functional and dispersion correction. A number of investigations have employed different approaches to compute adsorption energies of small molecules in zeolites, using reference values from high-level calculations and/or experiments. Such comparative studies are lacking for larger functional organic molecules such as pharmaceuticals or personal care products, despite their potential relevance for applications, e. g., in contaminant removal or drug delivery. The present study aims to fill this gap by comparing adsorption energies and, for selected cases, equilibrium structures of emerging organic contaminants adsorbed in MOR- and FAU-type all-silica zeolites. A total of 13 dispersion-corrected DFT approaches are compared, including methods using a pairwise dispersion correction as well as non-local van der Waals density functionals. While absolute values of adsorption energies vary widely, qualitative trends across the set of zeolite-guest combinations are not strongly dependent on the choice of functional. For selected cluster models, DFT adsorption energies are compared to reference values from coupled cluster (DLPNO-CCSD(T)) calculations. Although all DFT approaches deliver systematically more negative adsorption energies than the coupled cluster reference, this tendency is least pronounced for the rev-vdW-DF2 functional.
Collapse
Affiliation(s)
- Michael Fischer
- Crystallography and GeomaterialsFaculty of GeosciencesUniversity of BremenKlagenfurter Straße 2–428359BremenGermany
- Bremen Center for Computational Materials Science and MAPEX Center for Materials and ProcessesUniversity of Bremen28359BremenGermany
| | - Jakob Brauer
- Crystallography and GeomaterialsFaculty of GeosciencesUniversity of BremenKlagenfurter Straße 2–428359BremenGermany
- Bremen Center for Computational Materials Science and MAPEX Center for Materials and ProcessesUniversity of Bremen28359BremenGermany
| |
Collapse
|
3
|
Desmarais JK, Maul J, Civalleri B, Erba A, Vignale G, Pittalis S. Spin Currents via the Gauge Principle for Meta-Generalized Gradient Exchange-Correlation Functionals. PHYSICAL REVIEW LETTERS 2024; 132:256401. [PMID: 38996240 DOI: 10.1103/physrevlett.132.256401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Accepted: 05/15/2024] [Indexed: 07/14/2024]
Abstract
The prominence of density functional theory in the field of electronic structure computation stems from its ability to usefully balance accuracy and computational effort. At the base of this ability is a functional of the electron density: the exchange-correlation energy. This functional satisfies known exact conditions that guide the derivation of approximations. The strongly constrained and appropriately normed (SCAN) approximation stands out as a successful, modern, example. In this Letter, we demonstrate how the SU(2) gauge invariance of the exchange-correlation functional in spin current density functional theory allows us to add an explicit dependence on spin currents in the SCAN functional (here called JSCAN)-and similar meta-generalized-gradient functional approximations-solely invoking first principles. In passing, a spin-current dependent generalization of the electron localization function (here called JELF) is also derived. The extended forms are implemented in a developer's version of the crystal23 program. Applications on molecules and materials confirm the practical relevance of the extensions.
Collapse
Affiliation(s)
| | | | | | | | | | - Stefano Pittalis
- Istituto Nanoscienze, Consiglio Nazionale delle Ricerche, Via Campi 213A, I-41125 Modena, Italy
| |
Collapse
|
4
|
Riedel ZW, Shoemaker DP. Design Rules, Accurate Enthalpy Prediction, and Synthesis of Stoichiometric Eu 3+ Quantum Memory Candidates. J Am Chem Soc 2024; 146:2113-2121. [PMID: 38214913 DOI: 10.1021/jacs.3c11615] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
Stoichiometric Eu3+ compounds have recently shown promise for building dense, optically addressable quantum memory as the cations' long nuclear spin coherence times and shielded 4f electron optical transitions provide reliable memory platforms. Implementing such a system, though, requires ultranarrow, inhomogeneous linewidth compounds. Finding this rare linewidth behavior within a wide range of potential chemical spaces remains difficult, and while exploratory synthesis is often guided by density functional theory (DFT) calculations, lanthanides' 4f electrons pose unique challenges for stability predictions. Here, we report DFT procedures that reliably reproduce known phase diagrams and correctly predict two experimentally realized quantum memory candidates. We are the first to synthesize the double perovskite halide Cs2NaEuF6. It is an air-stable compound with a calculated band gap of 5.0 eV that surrounds Eu3+ with mononuclidic elements, which are desirable for avoiding inhomogeneous linewidth broadening. We also analyze computational database entries to identify phosphates and iodates as the next generation of chemical spaces for stoichiometric quantum memory system studies. This work identifies new candidate platforms for exploring chemical effects on quantum memory candidates' inhomogeneous linewidth while also providing a framework for screening Eu3+ compound stability with DFT.
Collapse
Affiliation(s)
- Zachary W Riedel
- Department of Materials Science and Engineering, University of Illinois Urbana─Champaign, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois Urbana─Champaign, Urbana, Illinois 61801, United States
| | - Daniel P Shoemaker
- Department of Materials Science and Engineering, University of Illinois Urbana─Champaign, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois Urbana─Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
5
|
Liu H, Bai X, Ning J, Hou Y, Song Z, Ramasamy A, Zhang R, Li Y, Sun J, Xiao B. Assessing r2SCAN meta-GGA functional for structural parameters, cohesive energy, mechanical modulus, and thermophysical properties of 3d, 4d, and 5d transition metals. J Chem Phys 2024; 160:024102. [PMID: 38189614 DOI: 10.1063/5.0176415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Accepted: 12/17/2023] [Indexed: 01/09/2024] Open
Abstract
The recent development of accurate and efficient semilocal density functionals on the third rung of Jacob's ladder of density functional theory, such as the revised regularized strongly constrained and appropriately normed (r2SCAN) density functional, could enable rapid and highly reliable prediction of the elasticity and temperature dependence of thermophysical parameters of refractory elements and their intermetallic compounds using the quasi-harmonic approximation (QHA). Here, we present a comparative evaluation of equilibrium cell volumes, cohesive energy, mechanical moduli, and thermophysical properties (Debye temperature and thermal expansion coefficient) for 22 transition metals using semilocal density functionals, including the local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE) and PBEsol generalized gradient approximations (GGAs), and the r2SCAN meta-GGA. PBEsol and r2SCAN deliver the same level of accuracies for structural, mechanical, and thermophysical properties. PBE and r2SCAN perform better than LDA and PBEsol for calculating cohesive energies of transition metals. Among the tested density functionals, r2SCAN provides an overall well-balanced performance for reliably computing cell volumes, cohesive energies, mechanical properties, and thermophysical properties of various 3d, 4d, and 5d transition metals using QHA. Therefore, we recommend that r2SCAN could be employed as a workhorse method to evaluate thermophysical properties of transition metal compounds and alloys in high throughput workflows.
Collapse
Affiliation(s)
- Haoliang Liu
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Xue Bai
- State Key Laboratory for Mechanical Behavior of Materials, Xi'an Jiaotong University, Xi'an 710049, China
| | - Jinliang Ning
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - Yuxuan Hou
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Zifeng Song
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Akilan Ramasamy
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - Ruiqi Zhang
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - Yefei Li
- State Key Laboratory for Mechanical Behavior of Materials, Xi'an Jiaotong University, Xi'an 710049, China
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - Bing Xiao
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| |
Collapse
|
6
|
Kovács P, Blaha P, Madsen GKH. Origin of the success of mGGAs for bandgaps. J Chem Phys 2023; 159:244118. [PMID: 38156634 DOI: 10.1063/5.0179260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Accepted: 12/05/2023] [Indexed: 01/03/2024] Open
Abstract
One of the well-known limitations of Kohn-Sham density functional theory is the tendency to strongly underestimate bandgaps. Meta-generalized gradient approximations (mGGAs), which include the kinetic energy density in the functional form, have been shown to significantly alleviate this deficiency. In this study, we explore the mechanisms responsible for this improvement from the angle of the underlying local densities. We find that the highest occupied and lowest unoccupied states are distinct in the space of the underlying descriptors. The gap opening is compared to a simple scaling of the local density approximation, and two mechanisms responsible for opening the mGGA gaps are identified. First of all, the relatively large negative derivative of the functional form with respect to reduced kinetic energy tends to elevate the lowest unoccupied state. Second, the curvature of functional, which ensures that it is bounded, tends to lower the highest occupied state. Remarkably, these two mechanisms are found to be transferable over a large and diverse database of compounds.
Collapse
Affiliation(s)
- Péter Kovács
- Institute of Materials Chemistry, Technical University of Vienna, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Peter Blaha
- Institute of Materials Chemistry, Technical University of Vienna, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Georg K H Madsen
- Institute of Materials Chemistry, Technical University of Vienna, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| |
Collapse
|
7
|
Tchakoua T, Jansen T, van Nies Y, van den Elshout RFA, van Boxmeer BAB, Poort SP, Ackermans MG, Beltrão GS, Hildebrand SA, Beekman SEJ, van der Drift T, Kaart S, Šantić A, Spuijbroek EE, Gerrits N, Somers MF, Kroes GJ. Constructing Mixed Density Functionals for Describing Dissociative Chemisorption on Metal Surfaces: Basic Principles. J Phys Chem A 2023; 127:10481-10498. [PMID: 38051300 PMCID: PMC10726370 DOI: 10.1021/acs.jpca.3c01932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
The production of a majority of chemicals involves heterogeneous catalysis at some stage, and the rates of many heterogeneously catalyzed processes are governed by transition states for dissociative chemisorption on metals. Accurate values of barrier heights for dissociative chemisorption on metals are therefore important to benchmarking electronic structure theory in general and density functionals in particular. Such accurate barriers can be obtained using the semiempirical specific reaction parameter (SRP) approach to density functional theory. However, this approach has thus far been rather ad hoc in its choice of the generic expression of the SRP functional to be used, and there is a need for better heuristic approaches to determining the mixing parameters contained in such expressions. Here we address these two issues. We investigate the ability of several mixed, parametrized density functional expressions combining exchange at the generalized gradient approximation (GGA) level with either GGA or nonlocal correlation to reproduce barrier heights for dissociative chemisorption on metal surfaces. For this, seven expressions of such mixed density functionals are tested on a database consisting of results for 16 systems taken from a recently published slightly larger database called SBH17. Three expressions are derived that exhibit high tunability and use correlation functionals that are either of the PBE GGA form or of one of two limiting nonlocal forms also describing the attractive van der Waals interaction in an approximate way. We also find that, for mixed density functionals incorporating GGA correlation, the optimum fraction of repulsive RPBE GGA exchange obtained with a specific GGA density functional is correlated with the charge-transfer parameter, which is equal to the difference in the work function of the metal surface and the electron affinity of the molecule. However, the correlation is generally not large and not large enough to obtain accurate guesses of the mixing parameter for the systems considered, suggesting that it does not give rise to a very effective search strategy.
Collapse
Affiliation(s)
- Théophile Tchakoua
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Tim Jansen
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Youri van Nies
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | | | - Bart A B van Boxmeer
- Faculty of Applied Engineering, University of Antwerp, Groenenborgerlaan 171, 2020 Antwerpen, Belgium
| | - Saskia P Poort
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Michelle G Ackermans
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Gabriel Spiller Beltrão
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Stefan A Hildebrand
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Steijn E J Beekman
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Thijs van der Drift
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Sam Kaart
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Anthonie Šantić
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Esmee E Spuijbroek
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Nick Gerrits
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
8
|
Francisco H, Cancio AC, Trickey SB. Reworking the Tao-Mo exchange-correlation functional. I. Reconsideration and simplification. J Chem Phys 2023; 159:214102. [PMID: 38038713 DOI: 10.1063/5.0167868] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2023] [Accepted: 10/16/2023] [Indexed: 12/02/2023] Open
Abstract
The revised, regularized Tao-Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn-Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao-Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew-Burke-Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization.
Collapse
Affiliation(s)
- H Francisco
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - A C Cancio
- Center for Computational Nanoscience, Department of Physics and Astronomy, Ball State University, Muncie, Indiana 47306, USA
| | - S B Trickey
- Quantum Theory Project, Department of Physics and Department of Chemistry, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
9
|
Francisco H, Cancio AC, Trickey SB. Reworking the Tao-Mo exchange-correlation functional. II. De-orbitalization. J Chem Phys 2023; 159:214103. [PMID: 38038712 DOI: 10.1063/5.0167873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2023] [Accepted: 11/02/2023] [Indexed: 12/02/2023] Open
Abstract
In Paper I [H. Francisco, A. C. Cancio, and S. B. Trickey, J. Chem. Phys. 159, 214102 (2023)], we gave a regularization of the Tao-Mo exchange functional that removes the order-of-limits problem in the original Tao-Mo form and also eliminates the unphysical behavior introduced by an earlier regularization while essentially preserving compliance with the second-order gradient expansion. The resulting simplified, regularized (sregTM) functional delivers performance on standard molecular and solid state test sets equal to that of the earlier revised, regularized Tao-Mo functional. Here, we address de-orbitalization of that new sregTM into a pure density functional. We summarize the failures of the Mejía-Rodríguez and Trickey de-orbitalization strategy [Phys. Rev. A 96, 052512 (2017)] when used with both versions. We discuss how those failures apparently arise in the so-called z' indicator function and in substitutes for the reduced density Laplacian in the parent functionals. Then, we show that the sregTM functional can be de-orbitalized somewhat well with a rather peculiarly parameterized version of the previously used deorbitalizer. We discuss, briefly, a de-orbitalization that works in the sense of reproducing error patterns but that apparently succeeds by cancelation of major qualitative errors associated with the de-orbitalized indicator functions α and z, hence, is not recommended. We suggest that the same issue underlies the earlier finding of comparatively mediocre performance of the de-orbitalized Tao-Perdew-Staroverov-Scuseri functional. Our work demonstrates that the intricacy of such two-indicator functionals magnifies the errors introduced by the Mejía-Rodríguez and Trickey de-orbitalization approach in ways that are extremely difficult to analyze and correct.
Collapse
Affiliation(s)
- H Francisco
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - A C Cancio
- Department of Physics and Astronomy, Ball State University, Muncie, Indiana 47306, USA
| | - S B Trickey
- Quantum Theory Project, Department of Physics and Department of Chemistry, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
10
|
Ricard TC, Zhu X, Iyengar SS. Capturing Weak Interactions in Surface Adsorbate Systems at Coupled Cluster Accuracy: A Graph-Theoretic Molecular Fragmentation Approach Improved through Machine Learning. J Chem Theory Comput 2023. [PMID: 38019639 DOI: 10.1021/acs.jctc.3c00955] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2023]
Abstract
The accurate and efficient study of the interactions of organic matter with the surface of water is critical to a wide range of applications. For example, environmental studies have found that acidic polyfluorinated alkyl substances, especially perfluorooctanoic acid (PFOA), have spread throughout the environment and bioaccumulate into human populations residing near contaminated watersheds, leading to many systemic maladies. Thus, the study of the interactions of PFOA with water surfaces became important for the mitigation of their activity as pollutants and threats to public health. However, theoretical study of the interactions of such organic adsorbates on the surface of water, and their bulk concerted properties, often necessitates the use of ab initio methods to properly incorporate the long-range electronic properties that govern these extended systems. Notable theoretical treatments of "on-water" reactions thus far have employed hybrid DFT and semilocal DFT, but the interactions involved are weak interactions that may be best described using post-Hartree-Fock theory. Here, we aim to demonstrate the utility of a graph-theoretic approach to molecular fragmentation that accurately captures the critical "weak" interactions while maintaining an efficient ab initio treatment of the long-range periodic interactions that underpin the physics of extended systems. We apply this graph-theoretical treatment to study PFOA on the surface of water as a model system for the study of weak interactions seen in the wide range of surface interactions and reactions. The approach divides a system into a set of vertices, that are then connected through edges, faces, and higher order graph theoretic objects known as simplexes, to represent a collection of locally interacting subsystems. These subsystems are then used to construct ab initio molecular dynamics simulations and for computing multidimensional potential energy surfaces. To further improve the computational efficiency of our graph theoretic fragmentation method, we use a recently developed transfer learning protocol to construct the full system potential energy from a family of neural networks each designed to accurately model the behavior of individual simplexes. We use a unique multidimensional clustering algorithm, based on the k-means clustering methodology, to define our training space for each separate simplex. These models are used to extrapolate the energies for molecular dynamics trajectories at PFOA water interfaces, at less than one-tenth the cost as compared to a regular molecular fragmentation-based dynamics calculation with excellent agreement with couple cluster level of full system potential energies.
Collapse
Affiliation(s)
- Timothy C Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Xiao Zhu
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
11
|
Erba A, Desmarais JK, Casassa S, Civalleri B, Donà L, Bush IJ, Searle B, Maschio L, Edith-Daga L, Cossard A, Ribaldone C, Ascrizzi E, Marana NL, Flament JP, Kirtman B. CRYSTAL23: A Program for Computational Solid State Physics and Chemistry. J Chem Theory Comput 2023; 19:6891-6932. [PMID: 36502394 PMCID: PMC10601489 DOI: 10.1021/acs.jctc.2c00958] [Citation(s) in RCA: 33] [Impact Index Per Article: 33.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Indexed: 12/14/2022]
Abstract
The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.
Collapse
Affiliation(s)
- Alessandro Erba
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jacques K. Desmarais
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Silvia Casassa
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Bartolomeo Civalleri
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Lorenzo Donà
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Ian J. Bush
- STFC
Rutherford Appleton Laboratory, Chilton Didcot, Oxfordshire OX11 0QX, United Kingdom
| | - Barry Searle
- SFTC
Daresbury Laboratory, Daresbury, Cheshire WA4 4AD, United Kingdom
| | - Lorenzo Maschio
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Loredana Edith-Daga
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Alessandro Cossard
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Chiara Ribaldone
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Eleonora Ascrizzi
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Naiara L. Marana
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université
de Lille, CNRS, UMR 8523 — PhLAM — Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department
of Chemistry and Biochemistry, University
of California, Santa
Barbara, California 93106, United States
| |
Collapse
|
12
|
Maciaszek M, Žalandauskas V, Silkinis R, Alkauskas A, Razinkovas L. The application of the SCAN density functional to color centers in diamond. J Chem Phys 2023; 159:084708. [PMID: 37642256 DOI: 10.1063/5.0154319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Accepted: 07/31/2023] [Indexed: 08/31/2023] Open
Abstract
Detailed characterization of deep-level color centers requires understanding their electronic and atomic structure, which is most commonly investigated utilizing the Kohn-Sham density functional theory. Standard semilocal functionals based on the generalized gradient approximation (GGA) are inclined toward an imprecise quantitative description of defects' electronic structure. Hybrid functionals provide an improved prediction of electronic properties, albeit at a much higher computational cost. In this work, we test the newly developed Strongly Constrained and Appropriately Normed (SCAN) family of meta-GGA density functionals for selected color centers in diamond. In particular, we study nitrogen-, silicon-, germanium-, and tin-vacancy centers that have been recently investigated for their use in quantum technological applications. We show that SCAN and its derivatives, the rSCAN and r2SCAN functionals, significantly improve the calculated energies of optical transitions within the delta-self-consistent-field approach, almost reaching the accuracy of the hybrid Heyd-Scuseria-Ernzerhof (HSE) functional. In the case of the NV- center, we also show that the SCAN family of functionals improves the description of the adiabatic potential energy surfaces compared to both GGA and hybrid functionals, improving calculated luminescence lineshapes. As a result of these findings, we recommend using the SCAN family of functionals as a promising alternative for studying color centers in solids.
Collapse
Affiliation(s)
- Marek Maciaszek
- Center for Physical Sciences and Technology (FTMC), Vilnius LT-10257, Lithuania
- Faculty of Physics, Warsaw University of Technology, Koszykowa 75, 00-662 Warsaw, Poland
| | | | - Rokas Silkinis
- Center for Physical Sciences and Technology (FTMC), Vilnius LT-10257, Lithuania
| | - Audrius Alkauskas
- Center for Physical Sciences and Technology (FTMC), Vilnius LT-10257, Lithuania
| | - Lukas Razinkovas
- Center for Physical Sciences and Technology (FTMC), Vilnius LT-10257, Lithuania
| |
Collapse
|
13
|
Swathilakshmi S, Devi R, Sai Gautam G. Performance of the r 2SCAN Functional in Transition Metal Oxides. J Chem Theory Comput 2023. [PMID: 37329316 DOI: 10.1021/acs.jctc.3c00030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
We assess the accuracy and computational efficiency of the recently developed meta-generalized gradient approximation (metaGGA) functional, restored regularized strongly constrained and appropriately normed (r2SCAN), in transition metal oxide (TMO) systems and compare its performance against SCAN. Specifically, we benchmark the r2SCAN-calculated oxidation enthalpies, lattice parameters, on-site magnetic moments, and band gaps of binary 3d TMOs against the SCAN-calculated and experimental values. Additionally, we evaluate the optimal Hubbard U correction required for each transition metal (TM) to improve the accuracy of the r2SCAN functional, based on experimental oxidation enthalpies, and verify the transferability of the U values by comparing against experimental properties on other TM-containing oxides. Notably, including the U-correction with r2SCAN increases the lattice parameters, on-site magnetic moments, and band gaps of TMOs, apart from an improved description of the ground state electronic state in narrow band gap TMOs. The r2SCAN and r2SCAN+U calculated oxidation enthalpies follow the qualitative trends of SCAN and SCAN+U, with r2SCAN and r2SCAN+U predicting marginally larger lattice parameters, smaller magnetic moments, and lower band gaps compared to SCAN and SCAN+U, respectively. We observe the overall computational time (i.e., for all ionic+electronic steps) required for r2SCAN(+U) to be lower than SCAN(+U). Thus, the r2SCAN(+U) framework can offer a reasonably accurate description of the ground state properties of TMOs with better computational efficiency than SCAN(+U).
Collapse
Affiliation(s)
- S Swathilakshmi
- Department of Materials Engineering, Indian Institute of Science, Bengaluru 560012, India
| | - Reshma Devi
- Department of Materials Engineering, Indian Institute of Science, Bengaluru 560012, India
| | | |
Collapse
|
14
|
Garrity KF, Choudhary K. Fast and accurate prediction of material properties with three-body tight-binding model for the periodic table. PHYSICAL REVIEW MATERIALS 2023; 7:10.1103/physrevmaterials.7.044603. [PMID: 37274125 PMCID: PMC10236923 DOI: 10.1103/physrevmaterials.7.044603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Parametrized tight-binding models fit to first-principles calculations can provide an efficient and accurate quantum mechanical method for predicting properties of molecules and solids. However, well-tested parameter sets are generally only available for a limited number of atom combinations, making routine use of this method difficult. Furthermore, many previous models consider only simple two-body interactions, which limits accuracy. To tackle these challenges, we develop a density functional theory database of nearly 1 000 000 materials, which we use to fit a universal set of tight-binding parameters for 65 elements and their binary combinations. We include both two-body and three-body effective interaction terms in our model, plus self-consistent charge transfer, enabling our model to work for metallic, covalent, and ionic bonds with the same parameter set. To ensure predictive power, we adopt a learning framework where we repeatedly test the model on new low-energy crystal structures and then add them to the fitting data set, iterating until predictions improve. We distribute the materials database and tools developed in this paper publicly.
Collapse
|
15
|
Cytter Y, Nandy A, Duan C, Kulik HJ. Insights into the deviation from piecewise linearity in transition metal complexes from supervised machine learning models. Phys Chem Chem Phys 2023; 25:8103-8116. [PMID: 36876903 DOI: 10.1039/d3cp00258f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023]
Abstract
Virtual high-throughput screening (VHTS) and machine learning (ML) with density functional theory (DFT) suffer from inaccuracies from the underlying density functional approximation (DFA). Many of these inaccuracies can be traced to the lack of derivative discontinuity that leads to a curvature in the energy with electron addition or removal. Over a dataset of nearly one thousand transition metal complexes typical of VHTS applications, we computed and analyzed the average curvature (i.e., deviation from piecewise linearity) for 23 density functional approximations spanning multiple rungs of "Jacob's ladder". While we observe the expected dependence of the curvatures on Hartree-Fock exchange, we note limited correlation of curvature values between different rungs of "Jacob's ladder". We train ML models (i.e., artificial neural networks or ANNs) to predict the curvature and the associated frontier orbital energies for each of these 23 functionals and then interpret differences in curvature among the different DFAs through analysis of the ML models. Notably, we observe spin to play a much more important role in determining the curvature of range-separated and double hybrids in comparison to semi-local functionals, explaining why curvature values are weakly correlated between these and other families of functionals. Over a space of 187.2k hypothetical compounds, we use our ANNs to pinpoint DFAs for which representative transition metal complexes have near-zero curvature with low uncertainty, demonstrating an approach to accelerate screening of complexes with targeted optical gaps.
Collapse
Affiliation(s)
- Yael Cytter
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| |
Collapse
|
16
|
Song S, Vuckovic S, Kim Y, Yu H, Sim E, Burke K. Extending density functional theory with near chemical accuracy beyond pure water. Nat Commun 2023; 14:799. [PMID: 36781855 PMCID: PMC9925738 DOI: 10.1038/s41467-023-36094-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Accepted: 01/13/2023] [Indexed: 02/15/2023] Open
Abstract
Density functional simulations of condensed phase water are typically inaccurate, due to the inaccuracies of approximate functionals. A recent breakthrough showed that the SCAN approximation can yield chemical accuracy for pure water in all its phases, but only when its density is corrected. This is a crucial step toward first-principles biosimulations. However, weak dispersion forces are ubiquitous and play a key role in noncovalent interactions among biomolecules, but are not included in the new approach. Moreover, naïve inclusion of dispersion in HF-SCAN ruins its high accuracy for pure water. Here we show that systematic application of the principles of density-corrected DFT yields a functional (HF-r2SCAN-DC4) which recovers and not only improves over HF-SCAN for pure water, but also captures vital noncovalent interactions in biomolecules, making it suitable for simulations of solutions.
Collapse
Affiliation(s)
- Suhwan Song
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea ,grid.266093.80000 0001 0668 7243Department of Chemistry, University of California, Irvine, CA 92697 USA
| | - Stefan Vuckovic
- grid.472716.10000 0004 1758 7362Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy ,grid.12380.380000 0004 1754 9227Departments of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Youngsam Kim
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea
| | - Hayoung Yu
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722, Korea.
| | - Kieron Burke
- grid.266093.80000 0001 0668 7243Department of Chemistry, University of California, Irvine, CA 92697 USA ,grid.266093.80000 0001 0668 7243Departments of Physics & Astronomy, University of California, Irvine, CA 92697 USA
| |
Collapse
|
17
|
Gygi F. All-Electron Plane-Wave Electronic Structure Calculations. J Chem Theory Comput 2023; 19:1300-1309. [PMID: 36757291 PMCID: PMC9979607 DOI: 10.1021/acs.jctc.2c01191] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
We demonstrate the use of the plane wave basis for all-electron electronic structure calculations. The approach relies on the definition of an analytic, norm-conserving, regularized Coulomb potential, and a scalable implementation of the plane wave method capable of handling large energy cutoffs (up to 80 kRy in the examples shown). The method is applied to the computation of electronic properties of isolated atoms as well as the diamond and silicon crystals, MgO, solid argon, and a configuration of 64 water molecules extracted from a first-principles molecular dynamics simulation. The computed energies, band gaps, ionic forces, and stress tensors provide reference results for the validation of pseudopotentials and/or localized basis sets. A calculation of the all-electron band structure of diamond and silicon using the SCAN meta-GGA density functional allows for a validation of calculations based on pseudopotentials derived using the PBE exchange-correlation functional. In the case of (H2O)64, the computed ionic forces provide a reference from which the errors incurred in pseudopotential calculations and in localized Gaussian basis sets calculations can be estimated.
Collapse
|
18
|
Perez I. Ab initio methods for the computation of physical properties and performance parameters of electrochemical energy storage devices. Phys Chem Chem Phys 2023; 25:1476-1503. [PMID: 36602004 DOI: 10.1039/d2cp03611h] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
With the rapid development of electric vehicles and mobile technologies, there is a high demand for electrochemical energy storage devices and electrochemical energy conversion devices. Devices meeting these needs include metal-ion batteries (MIBs), supercapacitors (SCs), electrochromic devices (ECDs), and multifunctional devices such as electrochromic batteries and supercapatteries. Currently, the goal has been the enhancement of operational parameters and physical properties that results in a higher performance of these devices. In the case of batteries, SCs, and supercapatteries, scientists seek to improve the equilibrium voltage, energy density, power, capacitance, and charge rate. In the case of ECDs, the focus is on improvement of the optical modulation and coloration efficiency. However, synthesis and characterization of new materials, or of materials with optimized properties, is time consuming and highly expensive. Computational simulation of materials can expedite the experimental endeavor by modelling novel atomic structures and predicting device performance. This is possible using ab initio theories and applying physical principles that allow us to understand the underlying mechanisms governing the behavior of materials in these devices. Taking as a point of departure density functional theory (DFT), in this review, we discuss the first principles methods used for the computation of physical properties and performance parameters of electrochemical energy storage devices. A wide coverage of DFT is given, dealing with the strengths and weaknesses of the most popular functionals used in the field of electrochemical energy storage. With these tools, ab initio methods for the computation of basic properties such as effective mass, mobility, optical band gap, transmissivity, conductivity (ionic and electronic), and criteria for structure stability (cohesive energy, formation energy, adsorption energy, and phonon frequency) are addressed. We also highlight the first principles techniques for the calculation of performance parameters in MIBs, SCs, and ECDs.
Collapse
Affiliation(s)
- Israel Perez
- National Council of Science and Technology (CONACYT)-Department of Physics and Mathematics, Institute of Engineering and Technology, Universidad Autonoma de Ciudad Juarez, Av. del Charro 450 Col. Romero Partido, C.P. 32310, Juarez, Chihuahua, Mexico.
| |
Collapse
|
19
|
Zhang C, Verma P, Wang J, Liu Y, He X, Wang Y, Truhlar DG, Liu Z. Performance of Screened-Exchange Functionals for Band Gaps and Lattice Constants of Crystals. J Chem Theory Comput 2023; 19:311-323. [PMID: 36520598 DOI: 10.1021/acs.jctc.2c00822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Kohn-Sham density functional theory is the most widely used method for electronic structure calculations of solid-state systems. The screened-exchange functionals developed following the influential work of Scuseria and co-workers in 2003-2006 have significantly improved the accuracy of the predictions of solid-state properties. This work assesses six screened-exchange density functionals for the prediction of 60 band gaps (database BG60) and 68 lattice constants (database LC68). The band gaps are calculated with both consistently calculated lattice constants and experimental lattice constants. Results for the nonlocal screened-exchange functionals are compared with those for six widely used or recently developed local functionals. The results show that all the screened-exchange functionals have smaller mean absolute errors (MAEs) than any of the local functionals. All the functionals except HLE17 overestimate (on average) the lattice constants, and M06-SX gives the best performance among the compared functionals, with a MAE of 0.051 Å. All the functionals underestimate (on average) the band gaps, and M06-SX outperforms all other functionals, with a MAE of 0.47 eV. M06-SX also has the lowest root-mean-squared error for both LC68 and BG60. For the subdatabases of BG60, M06-SX shows better performance for ionic crystals and systems with large band gaps, while HSE12s gives better results for semiconductors and systems with small band gaps. Overall, M06-SX shows the best performance for solid-state systems, followed by N12-SX and HSE12s. The best-performing local functionals are M06-L, revM06-L, and HLE17 for band gaps and M06-L and revM06-L for lattice constants. We found that M06-SX, revM06-L, and N12-SX not only are well optimized for a broad array of chemical properties but also have very good performance for the databases in this paper, making them well-suited for applications involving heterogeneous chemistry.
Collapse
Affiliation(s)
- Cheng Zhang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Jiaxu Wang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| | - Yiwei Liu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| | - Ying Wang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zhonghua Liu
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| |
Collapse
|
20
|
Manian A, Hudson RJ, Ramkissoon P, Smith TA, Russo SP. Interexcited State Photophysics I: Benchmarking Density Functionals for Computing Nonadiabatic Couplings and Internal Conversion Rate Constants. J Chem Theory Comput 2023; 19:271-292. [PMID: 36490305 DOI: 10.1021/acs.jctc.2c00888] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We present the first benchmarking study of nonadiabatic matrix coupling elements (NACMEs) calculated using different density functionals. Using the S1 → S0 transition in perylene solvated in toluene as a case study, we calculate the photophysical properties and corresponding rate constants for a variety of density functionals from each rung of Jacob's ladder. The singlet photoluminescence quantum yield (sPLQY) is taken as a measure of accuracy, measured experimentally here as 0.955. Important quantum chemical parameters such as geometries, absorption, emission, and adiabatic energies, NACMEs, Hessians, and transition dipole moments were calculated for each density functional basis set combination (data set) using density functional theory based multireference configuration interaction (DFT/MRCI) and compared to experiment where possible. We were able to derive simple relations between the TDDFT and DFT/MRCI photophysical properties; with semiempirical damping factors of ∼0.843 ± 0.017 and ∼0.954 ± 0.064 for TDDFT transition dipole moments and energies to DFT/MRCI level approximations, respectively. NACMEs were dominated by out-of-plane derivative components belonging to the center-most ring atoms with weaker contributions from perturbations along the transverse and longitudinal axes. Calculated theoretical spectra compared well to both experiment and literature, with fluorescence lifetimes between 7.1 and 12.5 ns, agreeing within a factor of 2 with experiment. Internal conversion (IC) rates were then calculated and were found to vary wildly between 106-1016 s-1 compared with an experimental rate of the order 107 s-1. Following further testing by mixing data sets, we found a strong dependence on the method used to obtain the Hessian. The 5 characterized data sets ranked in order of most promising are PBE0/def2-TZVP, ωB97XD/def2-TZVP, HCTH407/TZVP, PBE/TZVP, and PBE/def2-TZVP.
Collapse
Affiliation(s)
- Anjay Manian
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne3000, Australia
| | - Rohan J Hudson
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Melbourne, Parkville3010, Australia
| | - Pria Ramkissoon
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Melbourne, Parkville3010, Australia
| | - Trevor A Smith
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Melbourne, Parkville3010, Australia
| | - Salvy P Russo
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne3000, Australia
| |
Collapse
|
21
|
Tchakoua T, Gerrits N, Smeets EWF, Kroes GJ. SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces. J Chem Theory Comput 2022; 19:245-270. [PMID: 36529979 PMCID: PMC9835835 DOI: 10.1021/acs.jctc.2c00824] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.
Collapse
Affiliation(s)
- T. Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands
| | - N. Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,PLASMANT,
Department of Chemistry, University of Antwerp, BE-2610Antwerp, Belgium
| | - E. W. F. Smeets
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,ALTEN
Nederland, Technology, Fascinatio Boulevard 582, 2909 VACapelle a/d IJssel, The Netherlands
| | - G.-J. Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,E-mail: . Phone: +31 71 527 4396
| |
Collapse
|
22
|
Goldzak T, Wang X, Ye HZ, Berkelbach TC. Accurate thermochemistry of covalent and ionic solids from spin-component-scaled MP2. J Chem Phys 2022; 157:174112. [PMID: 36347707 PMCID: PMC9637026 DOI: 10.1063/5.0119633] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Accepted: 10/11/2022] [Indexed: 11/14/2022] Open
Abstract
We study the performance of spin-component-scaled second-order Møller-Plesset perturbation theory (SCS-MP2) for the prediction of the lattice constant, bulk modulus, and cohesive energy of 12 simple, three-dimensional covalent and ionic semiconductors and insulators. We find that SCS-MP2 and the simpler scaled opposite-spin MP2 (SOS-MP2) yield predictions that are significantly improved over the already good performance of MP2. Specifically, when compared to experimental values with zero-point vibrational corrections, SCS-MP2 (SOS-MP2) yields mean absolute errors of 0.015 (0.017) Å for the lattice constant, 3.8 (3.7) GPa for the bulk modulus, and 0.06 (0.08) eV for the cohesive energy, which are smaller than those of leading density functionals by about a factor of two or more. We consider a reparameterization of the spin-scaling parameters and find that the optimal parameters for these solids are very similar to those already in common use in molecular quantum chemistry, suggesting good transferability and reliable future applications to surface chemistry on insulators.
Collapse
Affiliation(s)
- Tamar Goldzak
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Xiao Wang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
23
|
Kovacs P, Tran F, Blaha P, Madsen G. What is the optimal mGGA exchange functional for solids? J Chem Phys 2022; 157:094110. [DOI: 10.1063/5.0098787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The space of GGA and mGGA exchange approximations is systematically explored bytraining 25 new functionals to produce accurate lattice parameter, cohesive energy andband gap predictions. The trained functionals are used to reproduce previous knowledge ina data-driven way and to understand the accuracy tradeoff between the mentioned proper-ties. The functionals are compared to notable mGGA functionals to analyze how changesin the enhancement factor maps influence the accuracy of the predictions. Some of thetrained functionals are found to perform on par with specialized functionals for band gaps,while outperforming them on the other two properties. The error surface of our trainedfunctionals can serve as a soft-limit of what mGGA functionals can achieve.
Collapse
Affiliation(s)
- Peter Kovacs
- Vienna University of Technology Faculty of Technical Chemistry, Austria
| | - Fabien Tran
- Institute of Materials Chemistry, Vienna University of Technology, Austria
| | - Peter Blaha
- Materials Chemistry, Vienna University of Technology Faculty of Technical Chemistry, Austria
| | | |
Collapse
|
24
|
Jana S, Constantin LA, Smiga S, Samal P. Solid-state performance of a meta-GGA screened hybrid density functional constructed from Pauli kinetic enhancement factor dependent semilocal exchange hole. J Chem Phys 2022; 157:024102. [DOI: 10.1063/5.0096674] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The semilocal form of the exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, that depends only on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is used subsequently to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional.
Collapse
Affiliation(s)
- Subrata Jana
- Department of Chemistry and Biochemistry, The Ohio State University, United States of America
| | - Lucian A. Constantin
- Istituto di Nanoscienze, Consiglio Nazionale delle Ricerche CNR-NANO, 41125 Modena, Italy, Italy
| | - Szymon Smiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University Institute of Physics, Poland
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, India
| |
Collapse
|
25
|
Khan I, Ullah I, Haq IU, Ali A, Dahshan A, Ali Z, Ahmad I. Structural and optoelectronic properties of CsLnZnTe3 (Ln = La, Pr, Nd and Sm). J RARE EARTH 2022. [DOI: 10.1016/j.jre.2022.03.021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
26
|
Schmidt J, Wang HC, Cerqueira TFT, Botti S, Marques MAL. A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals. Sci Data 2022; 9:64. [PMID: 35236866 PMCID: PMC8891291 DOI: 10.1038/s41597-022-01177-w] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Accepted: 01/25/2022] [Indexed: 11/09/2022] Open
Abstract
In the past decade we have witnessed the appearance of large databases of calculated material properties. These are most often obtained with the Perdew-Burke-Ernzerhof (PBE) functional of density-functional theory, a well established and reliable technique that is by now the standard in materials science. However, there have been recent theoretical developments that allow for increased accuracy in the calculations. Here, we present a dataset of calculations for 175k crystalline materials obtained with two functionals: geometry optimizations are performed with PBE for solids (PBEsol) that yields consistently better geometries than the PBE functional, and energies are obtained from PBEsol and from SCAN single-point calculations at the PBEsol geometry. Our results provide an accurate overview of the landscape of stable (and nearly stable) materials, and as such can be used for reliable predictions of novel compounds. They can also be used for training machine learning models, or even for the comparison and benchmark of PBE, PBEsol, and SCAN.
Collapse
Affiliation(s)
- Jonathan Schmidt
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120, Halle (Saale), Germany
| | - Hai-Chen Wang
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120, Halle (Saale), Germany
| | - Tiago F T Cerqueira
- CFisUC, Department of Physics, University of Coimbra, Rua Larga, 3004-516, Coimbra, Portugal
| | - Silvana Botti
- Institut für Festkörpertheorie und -optik and European Theoretical Spectroscopy Facility, Friedrich-Schiller-Universität Jena, D-07743, Jena, Germany
| | - Miguel A L Marques
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120, Halle (Saale), Germany.
| |
Collapse
|
27
|
Kuklin MS, Eklund K, Linnera J, Ropponen A, Tolvanen N, Karttunen AJ. Structural Properties and Magnetic Ground States of 100 Binary d-Metal Oxides Studied by Hybrid Density Functional Methods. Molecules 2022; 27:molecules27030874. [PMID: 35164135 PMCID: PMC8838575 DOI: 10.3390/molecules27030874] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Revised: 01/23/2022] [Accepted: 01/24/2022] [Indexed: 02/04/2023] Open
Abstract
d-metal oxides play a crucial role in numerous technological applications and show a great variety of magnetic properties. We have systematically investigated the structural properties, magnetic ground states, and fundamental electronic properties of 100 binary d-metal oxides using hybrid density functional methods and localized basis sets composed of Gaussian-type functions. The calculated properties are compared with experimental information in all cases where experimental data are available. The used PBE0 hybrid density functional method describes the structural properties of the studied d-metal oxides well, except in the case of molecular oxides with weak intermolecular forces between the molecular units. Empirical D3 dispersion correction does not improve the structural description of the molecular oxides. We provide a database of optimized geometries and magnetic ground states to facilitate future studies on the more complex properties of the binary d-metal oxides.
Collapse
|
28
|
Furness JW, Kaplan AD, Ning J, Perdew JP, Sun J. Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals. J Chem Phys 2022; 156:034109. [DOI: 10.1063/5.0073623] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- James W. Furness
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - Aaron D. Kaplan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Jinliang Ning
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - John P. Perdew
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| |
Collapse
|
29
|
Kovács P, Tran F, Hanbury A, Madsen GKH. Similarity Clustering for Representative Sets of Inorganic Solids for Density Functional Testing. J Chem Theory Comput 2022; 18:441-447. [PMID: 34919396 PMCID: PMC8757462 DOI: 10.1021/acs.jctc.1c00536] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Indexed: 11/30/2022]
Abstract
Benchmarking DFT functionals is complicated since the results highly depend on which properties and materials were used in the process. Unwanted biases can be introduced if a data set contains too many examples of very similar materials. We show that a clustering based on the distribution of density gradient and kinetic energy density is able to identify groups of chemically distinct solids. We then propose a method to create smaller data sets or rebalance existing data sets in a way that no region of the meta-GGA descriptor space is overrepresented, yet the new data set reproduces average errors of the original set as closely as possible. We apply the method to an existing set of 44 inorganic solids and suggest a representative set of seven solids. The representative sets generated with this method can be used to make more general benchmarks or to train new functionals.
Collapse
Affiliation(s)
- Péter Kovács
- Institute
of Materials Chemistry, Technical University
of Vienna, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Fabien Tran
- Institute
of Materials Chemistry, Technical University
of Vienna, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Allan Hanbury
- Institute
for Information Systems Engineering, Technical
University of Vienna, Favoritenstrasse 9-11/194, A-1040 Vienna, Austria
| | - Georg K. H. Madsen
- Institute
of Materials Chemistry, Technical University
of Vienna, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| |
Collapse
|
30
|
Understanding the effect of the exchange-correlation functionals on methane and ethane formation over ruthenium catalysts. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2203049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
31
|
Wu L, Tang SY, Zhou S. Computational Study on Homolytic Bond Energies of the Ag-X (X = C, O, and H) Complexes and Hammett-Type Analysis of Reactivity. ACS OMEGA 2021; 6:34904-34911. [PMID: 34963973 PMCID: PMC8697619 DOI: 10.1021/acsomega.1c05563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Accepted: 11/18/2021] [Indexed: 06/14/2023]
Abstract
Thirty-seven calculation methods were benchmarked against the available experimental bond lengths and energies data regarding the Ag-X bonds. The theoretical protocol PBE0/VDZ//ωB97x-D/mVTZ was found to be capable of accurately predicting the homolytic bond dissociation energies (BDEs) of Ag-X complexes with a precision of 1.9 kcal/mol. With the available method in hand, a wide range of different Ag-X BDEs were estimated. BDE(Ag-CH2X), BDE(Ag-PhX), BDE(Ag-OPhX), and BDE(Ag-OCOPhX) (X = NH2, OMe, Me, H, Cl, and NO2) were found to be in the ranges of 27-47, 51-54, 19-39, and 64-70 kcal/mol, respectively. Subsequently, Hammett-type analysis was carried out with reactivity parameters. Good positive linear relationships were found for BDE of Ag-O bands and decarboxylation barriers of Ag-OCOPhX with the Hammett constant σ. It suggested that electron-donating substituents could promote either the homolytic cleavage of the Ag-OPhX bond to undergo a radical process or Ag-OCOPhX decarboxylation. Moreover, ligand effects on Ag-H bonds were investigated using BDE(Ag-H) and related NPA charges on Ag. In the case of P-ligands, carbene ligands, and other small molecule ligands (i.e., CO, CO2, and H2O), a good negative linear relationship was found. In contrast, N-ligands could have a reverse effect. Understanding the intrinsic relationships of BDE(Ag-X) with related reactivity parameters might help gain insights into the structure-reactivity relationships in Ag-X-assisted C-H activation/decarboxylation.
Collapse
Affiliation(s)
- Lei Wu
- Zhejiang
Provincial Key Laboratory of Advanced Chemical Engineering Manufacture
Technology, College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310027, P. R. China
- SINOPEC
Research Institute of Safety Engineering, Qingdao 266000, P. R. China
- State
Key Laboratory of Safety and Control for Chemicals, Qingdao 266000, P. R. China
| | - Shi-Ya Tang
- SINOPEC
Research Institute of Safety Engineering, Qingdao 266000, P. R. China
- State
Key Laboratory of Safety and Control for Chemicals, Qingdao 266000, P. R. China
| | - Shaodong Zhou
- Zhejiang
Provincial Key Laboratory of Advanced Chemical Engineering Manufacture
Technology, College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310027, P. R. China
| |
Collapse
|
32
|
Wei F, Smeets EWF, Voss J, Kroes GJ, Lin S, Guo H. Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Fenfei Wei
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Egidius W. F. Smeets
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Johannes Voss
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park CA 94025, USA
| | - Geert-Jan Kroes
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Sen Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
33
|
Duan C, Chen S, Taylor MG, Liu F, Kulik HJ. Machine learning to tame divergent density functional approximations: a new path to consensus materials design principles. Chem Sci 2021; 12:13021-13036. [PMID: 34745533 PMCID: PMC8513898 DOI: 10.1039/d1sc03701c] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Accepted: 09/01/2021] [Indexed: 01/17/2023] Open
Abstract
Virtual high-throughput screening (VHTS) with density functional theory (DFT) and machine-learning (ML)-acceleration is essential in rapid materials discovery. By necessity, efficient DFT-based workflows are carried out with a single density functional approximation (DFA). Nevertheless, properties evaluated with different DFAs can be expected to disagree for cases with challenging electronic structure (e.g., open-shell transition-metal complexes, TMCs) for which rapid screening is most needed and accurate benchmarks are often unavailable. To quantify the effect of DFA bias, we introduce an approach to rapidly obtain property predictions from 23 representative DFAs spanning multiple families, “rungs” (e.g., semi-local to double hybrid) and basis sets on over 2000 TMCs. Although computed property values (e.g., spin state splitting and frontier orbital gap) differ by DFA, high linear correlations persist across all DFAs. We train independent ML models for each DFA and observe convergent trends in feature importance, providing DFA-invariant, universal design rules. We devise a strategy to train artificial neural network (ANN) models informed by all 23 DFAs and use them to predict properties (e.g., spin-splitting energy) of over 187k TMCs. By requiring consensus of the ANN-predicted DFA properties, we improve correspondence of computational lead compounds with literature-mined, experimental compounds over the typically employed single-DFA approach. Machine learning (ML)-based feature analysis reveals universal design rules regardless of density functional choices. Using the consensus among multiple functionals, we identify robust lead complexes in ML-accelerated chemical discovery.![]()
Collapse
Affiliation(s)
- Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584.,Department of Chemistry, Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Shuxin Chen
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584.,Department of Chemistry, Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge MA 02139 USA +1-617-253-4584
| |
Collapse
|
34
|
Keith JA, Vassilev-Galindo V, Cheng B, Chmiela S, Gastegger M, Müller KR, Tkatchenko A. Combining Machine Learning and Computational Chemistry for Predictive Insights Into Chemical Systems. Chem Rev 2021; 121:9816-9872. [PMID: 34232033 PMCID: PMC8391798 DOI: 10.1021/acs.chemrev.1c00107] [Citation(s) in RCA: 190] [Impact Index Per Article: 63.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Indexed: 12/23/2022]
Abstract
Machine learning models are poised to make a transformative impact on chemical sciences by dramatically accelerating computational algorithms and amplifying insights available from computational chemistry methods. However, achieving this requires a confluence and coaction of expertise in computer science and physical sciences. This Review is written for new and experienced researchers working at the intersection of both fields. We first provide concise tutorials of computational chemistry and machine learning methods, showing how insights involving both can be achieved. We follow with a critical review of noteworthy applications that demonstrate how computational chemistry and machine learning can be used together to provide insightful (and useful) predictions in molecular and materials modeling, retrosyntheses, catalysis, and drug design.
Collapse
Affiliation(s)
- John A. Keith
- Department
of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Valentin Vassilev-Galindo
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Bingqing Cheng
- Accelerate
Programme for Scientific Discovery, Department
of Computer Science and Technology, 15 J. J. Thomson Avenue, Cambridge CB3 0FD, United Kingdom
| | - Stefan Chmiela
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Michael Gastegger
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587, Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea
- Max-Planck-Institut für Informatik, 66123 Saarbrücken, Germany
- Google Research, Brain Team, 10117 Berlin, Germany
| | - Alexandre Tkatchenko
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
35
|
Jana S, Behera SK, Śmiga S, Constantin LA, Samal P. Accurate density functional made more versatile. J Chem Phys 2021; 155:024103. [PMID: 34266258 DOI: 10.1063/5.0051331] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We propose a one-electron self-interaction-free correlation energy functional compatible with the order-of-limit problem-free Tao-Mo (TM) semilocal functional (regTM) [J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016) and Patra et al., J. Chem. Phys. 153, 184112 (2020)] to be used for general purpose condensed matter physics and quantum chemistry. The assessment of the proposed functional for large classes of condensed matter and chemical systems shows its improvement in most cases compared to the TM functional, e.g., when applied to the relative energy difference of MnO2 polymorphs. In this respect, the present exchange-correction functional, which incorporates the TM technique of the exchange hole model combined with the slowly varying density correction, can achieve broad applicability, being able to solve difficult solid-state problems.
Collapse
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Sushant Kumar Behera
- 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
|
36
|
Huran AW, Wang HC, San-Miguel A, Marques MAL. Atomically Thin Pythagorean Tilings in Two Dimensions. J Phys Chem Lett 2021; 12:4972-4979. [PMID: 34014100 DOI: 10.1021/acs.jpclett.1c00903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We perform a theoretical study of an atomically thin, two-dimensional layer obtained by positioning atoms at the vertices of the classical Pythagorean tiling. This leads to an unusual geometrical pattern that is only stable for the three halogens Cl, Br, and I. In this Pythagorean structure, halogen atoms are arranged in strongly bound diatomic units that bind together by weaker electrostatic bonds. The energy of these phases is competitive with those of the low-temperature phase of the halogens and the two-dimensional layer obtained by exfoliating it. The Pythagorean layers are semiconducting, with an unusual band structure composed of very mobile holes and extremely heavy electrons. They are also soft, exhibiting small values of the elastic constants and a very low energy flexural mode. Analysis of the allowed Raman transitions reveals breathing-like modes that might be used to fingerprint, experimentally, the Pythagorean structure. Finally, we present a series of substrates that, due to lattice matching and compatible symmetry, can be used to stabilize these peculiar two-dimensional layers.
Collapse
Affiliation(s)
- Ahmad W Huran
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120 Halle (Saale), Germany
| | - Hai-Chen Wang
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120 Halle (Saale), Germany
| | - Alfonso San-Miguel
- Institut Lumière Matière, Univ. Lyon, Université Claude Bernard Lyon 1, CNRS, F-69622 Lyon, France
| | - Miguel A L Marques
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120 Halle (Saale), Germany
| |
Collapse
|
37
|
Peng D, Nakashima PNH. Measuring Density Functional Parameters from Electron Diffraction Patterns. PHYSICAL REVIEW LETTERS 2021; 126:176402. [PMID: 33988394 DOI: 10.1103/physrevlett.126.176402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/17/2020] [Revised: 01/13/2021] [Accepted: 03/15/2021] [Indexed: 06/12/2023]
Abstract
We integrate density functional theory (DFT) into quantitative convergent-beam electron diffraction (QCBED) to create a synergy between experiment and theory called QCBED-DFT. This synergy resides entirely in the electron density which, in real materials, gives rise to the experimental CBED patterns used by QCBED-DFT to refine DFT model parameters. We use it to measure the Hubbard energy U for two strongly correlated electron systems, NiO and CeB_{6} (U=7.4±0.6 eV for d orbitals in NiO and U=3.0±0.6 eV for f orbitals in CeB_{6}), and the boron position parameter x for CeB_{6} (x=0.1992±0.0003). In verifying our measurements, we demonstrate an accuracy test for any modeled electron density.
Collapse
Affiliation(s)
- Ding Peng
- Department of Materials Science and Engineering, Monash University, Victoria 3800, Australia
| | - Philip N H Nakashima
- Department of Materials Science and Engineering, Monash University, Victoria 3800, Australia
| |
Collapse
|
38
|
Kraisler E, Hodgson MJP, Gross EKU. From Kohn-Sham to Many-Electron Energies via Step Structures in the Exchange-Correlation Potential. J Chem Theory Comput 2021; 17:1390-1407. [PMID: 33595312 PMCID: PMC8363072 DOI: 10.1021/acs.jctc.0c01093] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Accurately
describing excited states within Kohn–Sham (KS)
density functional theory (DFT), particularly those which induce ionization
and charge transfer, remains a great challenge. Common exchange-correlation
(xc) approximations are unreliable for excited states owing, in part,
to the absence of a derivative discontinuity in the xc energy (Δ),
which relates a many-electron energy difference to the corresponding
KS energy difference. We demonstrate, analytically and numerically,
how the relationship between KS and many-electron energies leads to
the step structures observed in the exact xc potential in four scenarios:
electron addition, molecular dissociation, excitation of a finite
system, and charge transfer. We further show that steps in the potential
can be obtained also with common xc approximations, as simple as the
LDA, when addressed from the ensemble perspective. The article therefore
highlights how capturing the relationship between KS and many-electron
energies with advanced xc approximations is crucial for accurately
calculating excitations, as well as the ground-state density and energy
of systems which consist of distinct subsystems.
Collapse
Affiliation(s)
- Eli Kraisler
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| | - M J P Hodgson
- Department of Physics, Durham University, South Road, Durham DH1 3LE, United Kingdom.,Max-Planck-Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - E K U Gross
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| |
Collapse
|
39
|
Li HX, Zhou XY, Wang YC, Jiang H. Theoretical study of Na+ transport in the solid-state electrolyte Na3OBr based on deep potential molecular dynamics. Inorg Chem Front 2021. [DOI: 10.1039/d0qi00921k] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Deep potential molecular dynamics is used to study Na+ transport in Na3OBr.
Collapse
Affiliation(s)
- Han-Xiao Li
- Beijing National Laboratory of Molecular Science
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
- China
| | - Xu-Yuan Zhou
- Beijing National Laboratory of Molecular Science
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
- China
| | - Yue-Chao Wang
- Beijing National Laboratory of Molecular Science
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
- China
| | - Hong Jiang
- Beijing National Laboratory of Molecular Science
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
- China
| |
Collapse
|
40
|
GABRIEL JOSHUAJ, PAULSON NOAHH, DUONG THIENC, TAVAZZA FRANCESCA, BECKER CHANDLERA, CHAUDHURI SANTANU, STAN MARIUS. Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review. JOM (WARRENDALE, PA. : 1989) 2021; 73:10.1007/s11837-020-04436-6. [PMID: 34511862 PMCID: PMC8431950 DOI: 10.1007/s11837-020-04436-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2020] [Accepted: 10/05/2020] [Indexed: 06/13/2023]
Abstract
The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.
Collapse
Affiliation(s)
- JOSHUA J. GABRIEL
- Applied Materials Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - NOAH H. PAULSON
- Applied Materials Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - THIEN C. DUONG
- Energy and Global Security, Argonne National Laboratory, Lemont, IL 60439, USA
| | - FRANCESCA TAVAZZA
- Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| | - CHANDLER A. BECKER
- Office of Data and Informatics, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| | - SANTANU CHAUDHURI
- Manufacturing Science and Engineering, Energy and Global Security, Argonne National Laboratory, Lemont, IL 60439, USA
- Civil, Materials, and Environmental Engineering, University of Illinois at Chicago, Chicago, IL 60607, USA
| | - MARIUS STAN
- Applied Materials Division, Argonne National Laboratory, Lemont, IL 60439, USA
| |
Collapse
|
41
|
Jana S, Śmiga S, Constantin LA, Samal P. Generalizing Double-Hybrid Density Functionals: Impact of Higher-Order Perturbation Terms. J Chem Theory Comput 2020; 16:7413-7430. [PMID: 33205659 PMCID: PMC7735712 DOI: 10.1021/acs.jctc.0c00823] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Indexed: 11/28/2022]
Abstract
Connections between the Görling-Levy (GL) perturbation theory and the parameters of double-hybrid (DH) density functional are established via adiabatic connection formalism. Moreover, we present a more general DH density functional theory, where the higher-order perturbation terms beyond the second-order GL2 one, such as GL3 and GL4, also contribute. It is shown that a class of DH functionals including previously proposed ones can be formed using the present construction. Based on the proposed formalism, we assess the performance of higher-order DH and long-range corrected DH formed on the Perdew-Burke-Ernzerhof (PBE) semilocal functional and second-order GL2 correlation energy. The underlying construction of DH functionals based on the generalized many-body perturbation approaches is physically appealing in terms of the development of the non-local forms using more accurate and sophisticated semilocal functionals.
Collapse
Affiliation(s)
- Subrata Jana
- 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
- Consiglio
Nazionale delle Ricerche CNR-NANO, Istituto
di Nanoscienze, 41125 Modena, Italy
| | - Prasanjit Samal
- School
of Physical Sciences, National Institute
of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
42
|
Patra A, Jana S, Samal P. A way of resolving the order-of-limit problem of Tao–Mo semilocal functional. J Chem Phys 2020; 153:184112. [DOI: 10.1063/5.0025173] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
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
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
43
|
Mejía-Rodríguez D, Trickey SB. Spin-Crossover from a Well-Behaved, Low-Cost meta-GGA Density Functional. J Phys Chem A 2020; 124:9889-9894. [DOI: 10.1021/acs.jpca.0c08883] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Daniel Mejía-Rodríguez
- Center for Molecular Magnetic Quantum Materials, Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, United States
| | - S. B. Trickey
- Center for Molecular Magnetic Quantum Materials, Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
44
|
|
45
|
Furness JW, Kaplan AD, Ning J, Perdew JP, Sun J. Accurate and Numerically Efficient r 2SCAN Meta-Generalized Gradient Approximation. J Phys Chem Lett 2020; 11:8208-8215. [PMID: 32876454 DOI: 10.1021/acs.jpclett.0c02405] [Citation(s) in RCA: 251] [Impact Index Per Article: 62.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The recently proposed rSCAN functional [ J. Chem. Phys. 2019 150, 161101] is a regularized form of the SCAN functional [ Phys. Rev. Lett. 2015 115, 036402] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN.
Collapse
Affiliation(s)
- James W Furness
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| | - Aaron D Kaplan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Jinliang Ning
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| | - John P Perdew
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| |
Collapse
|
46
|
Yamamoto Y, Salcedo A, Diaz CM, Alam MS, Baruah T, Zope RR. Assessing the effect of regularization on the molecular properties predicted by SCAN and self-interaction corrected SCAN meta-GGA. Phys Chem Chem Phys 2020; 22:18060-18070. [PMID: 32760934 DOI: 10.1039/d0cp02717k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent regularization of the SCAN meta-GGA functional (rSCAN) has simplified the numerical complexities of the SCAN functional, alleviating SCAN's stringent demand on the numerical integration grids to some extent. The regularization of rSCAN, however, results in the breaking of some constraints such as the uniform electron gas limit, the slowly varying density limit, and coordinate scaling of the iso-orbital indicator. Here, we assess the effects of regularization on the electronic, structural, vibrational, and magnetic properties of molecules by comparing the SCAN and rSCAN predictions. The properties studied include atomic energies, atomization energies, ionization potentials, electron affinities, barrier heights, infrared intensities, dissociation and reaction energies, spin moments of molecular magnets, and isomer ordering of water clusters. Our results show that rSCAN requires less dense numerical grids and gives very similar results to those of SCAN for all properties examined with the exception of atomization energies, which are worsened in rSCAN. We also examine the performance of self-interaction-corrected (SIC) rSCAN with respect to SIC-SCAN using the Perdew-Zunger (PZ) SIC method. The PZSIC method uses orbital densities to compute one-electron self-interaction errors and places an even more stringent demand on numerical grids. Our results show that SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work with numerical grids that are on average half or less as dense as that needed for SIC-SCAN.
Collapse
Affiliation(s)
- Yoh Yamamoto
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA.
| | - Alan Salcedo
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA.
| | - Carlos M Diaz
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA. and Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Md Shamsul Alam
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA. and Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA. and Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R Zope
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA. and Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|
47
|
Tran A, Tranchida J, Wildey T, Thompson AP. Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys. J Chem Phys 2020; 153:074705. [PMID: 32828092 DOI: 10.1063/5.0015672] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials' design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum-niobium-titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.
Collapse
Affiliation(s)
- Anh Tran
- Optimization and Uncertainty Quantification, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Julien Tranchida
- Computational Multiscale, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Tim Wildey
- Optimization and Uncertainty Quantification, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Aidan P Thompson
- Computational Multiscale, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| |
Collapse
|
48
|
Otero-de-la-Roza A, Johnson ER. Application of XDM to ionic solids: The importance of dispersion for bulk moduli and crystal geometries. J Chem Phys 2020; 153:054121. [PMID: 32770899 DOI: 10.1063/5.0015133] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Dispersion corrections are essential in the description of intermolecular interactions; however, dispersion-corrected functionals must also be transferrable to hard solids. The exchange-hole dipole moment (XDM) model has demonstrated excellent performance for non-covalent interactions. In this article, we examine its ability to describe the relative stability, geometry, and compressibility of simple ionic solids. For the specific cases of the cesium halides, XDM-corrected functionals correctly predict the energy ranking of the B1 and B2 forms, and a dispersion contribution is required to obtain this result. Furthermore, for the lattice constants of the 20 alkali halides, the performance of XDM-corrected functionals is excellent, provided that the base functional's exchange enhancement factor properly captures non-bonded repulsion. The mean absolute errors in lattice constants obtained with B86bPBE-XDM and B86bPBE-25X-XDM are 0.060 Å and 0.039 Å, respectively, suggesting that delocalization error also plays a minor role in these systems. Finally, we considered the calculation of bulk moduli for alkali halides and alkaline-earth oxides. Previous claims in the literature that simple generalized gradient approximations, such as PBE, can reliably predict experimental bulk moduli have benefited from large error cancellations between neglecting both dispersion and vibrational effects. If vibrational effects are taken into account, dispersion-corrected functionals are quite accurate (4 GPa-5 GPa average error), again, if non-bonded repulsion is correctly represented. Careful comparisons of the calculated bulk moduli with experimental data are needed to avoid systematic biases and misleading conclusions.
Collapse
Affiliation(s)
- A Otero-de-la-Roza
- Departamento de Química Física y Analítica and MALTA Consolider Team, Facultad de Química, Universidad de Oviedo, 33006 Oviedo, Spain
| | - Erin R Johnson
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd., Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
49
|
Xu X, Jiang H. Anion order in perovskite oxynitrides AMO 2N (A = Ba, Sr, Ca; M = Ta, Nb): a first-principles based investigation. RSC Adv 2020; 10:24410-24418. [PMID: 35516175 PMCID: PMC9055086 DOI: 10.1039/d0ra03681a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Accepted: 06/13/2020] [Indexed: 11/21/2022] Open
Abstract
Perovskite-type oxynitrides have attracted a lot of research interest as emerging functional materials with promising wide applications. The ordering of O/N anions in perovskite oxynitrides plays an important role in determining their physical properties, while it is still challenging to characterize the actual anion order in a particular material and understand the underlying physics. In this work, we have investigated anion order in a series of perovskite oxynitrides AMO2N (A = Ba, Sr, Ca; M = Ta, Nb) through first-principles calculations and the cluster-expansion-model-based Monte Carlo simulations. In terms of cluster correlation functions, it can be explicitly demonstrated that short-range anion order is present in all these perovskite oxynitrides. In addition, the anion order varies with the temperature of thermal equilibrium and depends on the cation type. Special quasi-ordered structures are then constructed as representative structures by taking the calculated anion order at finite temperature into consideration and their band gaps and dielectric tensors are predicted by first-principles calculations and compared to experimental values.
Collapse
Affiliation(s)
- Xi Xu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| | - Hong Jiang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| |
Collapse
|
50
|
Hyldgaard P, Jiao Y, Shukla V. Screening nature of the van der Waals density functional method: a review and analysis of the many-body physics foundation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:393001. [PMID: 32213670 DOI: 10.1088/1361-648x/ab8250] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2020] [Accepted: 03/23/2020] [Indexed: 05/03/2023]
Abstract
We review the screening nature and many-body physics foundation of the van der Waals density functional (vdW-DF) method [Berland Ket al2015Rep. Prog. Phys.78066501], a systematic approach to construct truly nonlocal exchange-correlation energy density functionals. To that end we define and focus on a class of consistent vdW-DF versions that adhere to the Lindhard screening logic of the full method formulation. The consistent-exchange vdW-DF-cx version [Berland K and Hyldgaard P 2014Phys. Rev. B89035412] and its spin extension [Thonhauser Tet al2015Phys. Rev. Lett.115136402] represent the first examples of this class; in general, consistent vdW-DFs reflect a concerted expansion of a formal recast of the adiabatic-connection formula [Hyldgaard Pet al2014Phys. Rev. B90075148], an exponential summation of contributions to the local-field response, and the Dyson equation. We argue that the screening emphasis is essential because the exchange-correlation energy reflects an effective electrodynamics set by a long-range interaction. Two consequences are that (1) there are, in principle, no wiggle room in how one balances exchange and correlation, for example, in vdW-DF-cx, and that (2) consistent vdW-DFs have a formal structure that allows them to incorporate vertex-correction effects, at least in the case of levels that experience recoil-less interactions (for example, near the Fermi surface). We explore the extent to which the strictly nonempirical vdW-DF-cx formulation can serve as a systematic extension of the constraint-based semilocal functionals. For validation, we provide a complete survey of vdW-DF-cx performance for broad molecular processes, for the full set of 55 benchmarks in GMTKN55 [Goerigk Let al2017Phys. Chem. Chem. Phys.1932184] and comparing to the quantum-chemistry calculations that are summarized in that paper. We also provide new vdW-DF-cx results for metal surface energies and work functions that we compare to experiment. Finally, we use the screening insight to separate the vdW-DF nonlocal-correlation term into pure-vdW-interaction and local-field-susceptibility effects and present tools to compute and map the binding signatures.
Collapse
Affiliation(s)
- Per Hyldgaard
- Department of Microtechnology and Nanoscience (MC2), Chalmers University of Technology, SE-41296 Gothenburg, Sweden
| | - Yang Jiao
- Department of Microtechnology and Nanoscience (MC2), Chalmers University of Technology, SE-41296 Gothenburg, Sweden
| | - Vivekanand Shukla
- Department of Microtechnology and Nanoscience (MC2), Chalmers University of Technology, SE-41296 Gothenburg, Sweden
| |
Collapse
|