1
|
Voss J. Machine learning for accuracy in density functional approximations. J Comput Chem 2024; 45:1829-1845. [PMID: 38668453 DOI: 10.1002/jcc.27366] [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: 10/30/2023] [Revised: 02/16/2024] [Accepted: 03/25/2024] [Indexed: 07/21/2024]
Abstract
Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. Here, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.
Collapse
Affiliation(s)
- Johannes Voss
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, Menlo Park, California, USA
| |
Collapse
|
2
|
Cai Y, Michiels R, De Luca F, Neyts E, Tu X, Bogaerts A, Gerrits N. Improving Molecule-Metal Surface Reaction Networks Using the Meta-Generalized Gradient Approximation: CO 2 Hydrogenation. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:8611-8620. [PMID: 38835935 PMCID: PMC11145648 DOI: 10.1021/acs.jpcc.4c01110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/20/2024] [Revised: 05/07/2024] [Accepted: 05/09/2024] [Indexed: 06/06/2024]
Abstract
Density functional theory is widely used to gain insights into molecule-metal surface reaction networks, which is important for a better understanding of catalysis. However, it is well-known that generalized gradient approximation (GGA) density functionals (DFs), most often used for the study of reaction networks, struggle to correctly describe both gas-phase molecules and metal surfaces. Also, GGA DFs typically underestimate reaction barriers due to an underestimation of the self-interaction energy. Screened hybrid GGA DFs have been shown to reduce this problem but are currently intractable for wide usage. In this work, we use a more affordable meta-GGA (mGGA) DF in combination with a nonlocal correlation DF for the first time to study and gain new insights into a catalytically important surface reaction network, namely, CO2 hydrogenation on Cu. We show that the mGGA DF used, namely, rMS-RPBEl-rVV10, outperforms typical GGA DFs by providing similar or better predictions for metals and molecules, as well as molecule-metal surface adsorption and activation energies. Hence, it is a better choice for constructing molecule-metal surface reaction networks.
Collapse
Affiliation(s)
- Yuxiang Cai
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
- Department
of Electrical Engineering and Electronics, University of Liverpool, Liverpool L69 3GJ, U.K.
| | - Roel Michiels
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
| | - Federica De Luca
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
- Department
of ChiBioFarAM (Industrial Chemistry), ERIC aisbl and INSTM/CASPE, University of Messina, V.le F. Stagno d’Alcontres 31, Messina 98166, Italy
| | - Erik Neyts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
| | - Xin Tu
- Department
of Electrical Engineering and Electronics, University of Liverpool, Liverpool L69 3GJ, U.K.
| | - Annemie Bogaerts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
| | - Nick Gerrits
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, Leiden 2300 RA, The Netherlands
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Trepte K, Voss J. Data-driven and constrained optimization of semi-local exchange and nonlocal correlation functionals for materials and surface chemistry. J Comput Chem 2022; 43:1104-1112. [PMID: 35474584 DOI: 10.1002/jcc.26872] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Revised: 04/07/2022] [Accepted: 04/11/2022] [Indexed: 11/08/2022]
Abstract
Reliable predictions of surface chemical reaction energetics require an accurate description of both chemisorption and physisorption. Here, we present an empirical approach to simultaneously optimize semi-local exchange and nonlocal correlation of a density functional approximation to improve these energetics. A combination of reference data for solid bulk, surface, and gas-phase chemistry and physical exchange-correlation model constraints leads to the VCML-rVV10 exchange-correlation functional. Owing to the variety of training data, the applicability of VCML-rVV10 extends beyond surface chemistry simulations. It provides optimized gas phase reaction energetics and an accurate description of bulk lattice constants and elastic properties.
Collapse
Affiliation(s)
- Kai Trepte
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, Menlo Park, California, USA
| | - Johannes Voss
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, Menlo Park, California, USA
| |
Collapse
|
6
|
Li WL, Lininger CN, Chen K, Vaissier Welborn V, Rossomme E, Bell AT, Head-Gordon M, Head-Gordon T. Critical Role of Thermal Fluctuations for CO Binding on Electrocatalytic Metal Surfaces. JACS AU 2021; 1:1708-1718. [PMID: 34723274 PMCID: PMC8549055 DOI: 10.1021/jacsau.1c00300] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2021] [Indexed: 06/01/2023]
Abstract
This work considers the evaluation of density functional theory (DFT) when comparing against experimental observations of CO binding trends on the strong binding Pt(111) and intermediate binding Cu(111) and for weak binding Ag(111) and Au(111) surfaces important in electrocatalysis. By introducing thermal fluctuations using appropriate statistical mechanical NVT and NPT ensembles, we find that the RPBE and B97M-rV DFT functionals yield qualitatively better metal surface strain trends and CO enthalpies of binding for Cu(111) and Pt(111) than found at 0 K, thereby correcting the overbinding by 0.2 to 0.3 eV to yield better agreement with the enthalpies determined from experiment. The importance of dispersion effects are manifest for the weak CO binding Ag(111) and Au(111) surfaces at finite temperatures in which the RPBE functional does not bind CO at all, while the B97M-rV functional shows that the CO-metal interactions are a mixture of chemisorbed and physisorbed species with binding enthalpies that are within ∼0.05 eV of experiment. Across all M(111) surfaces, we show that the B97M-rV functional consistently predicts the correct atop site preference for all metals due to thermally induced surface distortions that preferentially favor the undercoordinated site. This study demonstrates the need to fully account for finite temperature fluctuations to make contact with the binding enthalpies from surface science experiments and electrocatalysis applications.
Collapse
Affiliation(s)
- Wan-Lu Li
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
| | - Christianna N. Lininger
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
| | - Kaixuan Chen
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
| | - Valerie Vaissier Welborn
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 26067, United States
| | - Elliot Rossomme
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
| | - Alexis T. Bell
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical
and Biomolecular Engineering, and Department of
Bioengineering, University of California, Berkeley, California 94720, United States
| |
Collapse
|
7
|
Smeets EF, Kroes GJ. Performance of Made Simple Meta-GGA Functionals with rVV10 Nonlocal Correlation for H 2 + Cu(111), D 2 + Ag(111), H 2 + Au(111), and D 2 + Pt(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2021; 125:8993-9010. [PMID: 34084265 PMCID: PMC8162760 DOI: 10.1021/acs.jpcc.0c11034] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Revised: 04/08/2021] [Indexed: 06/12/2023]
Abstract
Accurately modeling heterogeneous catalysis requires accurate descriptions of rate-controlling elementary reactions of molecules on metal surfaces, but standard density functionals (DFs) are not accurate enough for this. The problem can be solved with the specific reaction parameter approach to density functional theory (SRP-DFT), but the transferability of SRP DFs among chemically related systems is limited. We combine the MS-PBEl, MS-B86bl, and MS-RPBEl semilocal made simple (MS) meta-generalized gradient approximation (GGA) (mGGA) DFs with rVV10 nonlocal correlation, and we evaluate their performance for the hydrogen (H2) + Cu(111), deuterium (D2) + Ag(111), H2 + Au(111), and D2 + Pt(111) gas-surface systems. The three MS mGGA DFs that have been combined with rVV10 nonlocal correlation were not fitted to reproduce particular experiments, nor has the b parameter present in rVV10 been reoptimized. Of the three DFs obtained the MS-PBEl-rVV10 DF yields an excellent description of van der Waals well geometries. The three original MS mGGA DFs gave a highly accurate description of the metals, which was comparable in quality to that obtained with the PBEsol DF. Here, we find that combining the three original MS mGGA DFs with rVV10 nonlocal correlation comes at the cost of a slightly less accurate description of the metal. However, the description of the metal obtained in this way is still better than the descriptions obtained with SRP DFs specifically optimized for individual systems. Using the Born-Oppenheimer static surface (BOSS) model, simulations of molecular beam dissociative chemisorption experiments yield chemical accuracy for the D2 + Ag(111) and D2 + Pt(111) systems. A comparison between calculated and measured E 1/2(ν, J) parameters describing associative desorption suggests chemical accuracy for the associative desorption of H2 from Au(111) as well. Our results suggest that ascending Jacob's ladder to the mGGA rung yields increasingly more accurate results for gas-surface reactions of H2 (D2) interacting with late transition metals.
Collapse
Affiliation(s)
- Egidius
W. F. Smeets
- Gorlaeus Laboratories, Leiden
Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Gorlaeus Laboratories, Leiden
Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
8
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (Eb, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate Eb can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for Eb for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH4. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of Eb for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate Eb for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O2, H2O, NH3, CO2, and CH3OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
9
|
Lininger CN, Gauthier JA, Li WL, Rossomme E, Welborn VV, Lin Z, Head-Gordon T, Head-Gordon M, Bell AT. Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals. Phys Chem Chem Phys 2021; 23:9394-9406. [DOI: 10.1039/d0cp03821k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
We assess four DFT functionals, RTPSS, RPBE, SCAN and B97M-rV, for surface interactions. We find that B97M-rV predicts the correct site preference for CO binding on Ag and Au while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.
Collapse
Affiliation(s)
- Christianna N. Lininger
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Department of Chemical and Biomolecular Engineering
| | - Joseph A. Gauthier
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Department of Chemical and Biomolecular Engineering
| | - Wan-Lu Li
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Kenneth S. Pitzer Center for Theoretical Chemistry
| | - Elliot Rossomme
- Kenneth S. Pitzer Center for Theoretical Chemistry
- Department of Chemistry
- University of California
- Berkeley
- USA
| | - Valerie Vaissier Welborn
- Kenneth S. Pitzer Center for Theoretical Chemistry
- Department of Chemistry
- University of California
- Berkeley
- USA
| | - Zhou Lin
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Kenneth S. Pitzer Center for Theoretical Chemistry
| | - Teresa Head-Gordon
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Department of Chemical and Biomolecular Engineering
| | - Martin Head-Gordon
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Kenneth S. Pitzer Center for Theoretical Chemistry
| | - Alexis T. Bell
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Department of Chemical and Biomolecular Engineering
| |
Collapse
|
10
|
Kabalan L, Kowalec I, Catlow CRA, Logsdail AJ. A computational study of the properties of low- and high-index Pd, Cu and Zn surfaces. Phys Chem Chem Phys 2021; 23:14649-14661. [PMID: 34212951 DOI: 10.1039/d1cp01602d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report a detailed Density Functional Theory (DFT) based investigation of the structure and stability of bulk and surface structures for the Group 10-12 elements Pd, Cu and Zn, considering the effect of the choice of exchange-correlation density functional and computation parameters. For the initial bulk structures, the lattice parameter and cohesive energy are calculated, which are then augmented by calculation of surface energies and work functions for the lower-index surfaces. Of the 22 density functionals considered, we highlight the mBEEF density functional as providing the best overall agreement with experimental data. The optimal density functional choice is applied to the study of higher index surfaces for the three metals, and Wulff constructions performed for nanoparticles with a radius of 11 nm, commensurate with nanoparticle sizes commonly employed in catalytic chemistry. For Pd and Cu, the low-index (111) facet is dominant in the constructed nanoparticles, covering ∼50% of the surface, with (100) facets covering a further 10 to 25%; however, non-negligible coverage from higher index (332), (332) and (210) facets is also observed for Pd, and (322), (221) and (210) surfaces are observed for Cu. In contrast, only the (0001) and (10-10) facets are observed for Zn. Overall, our results highlight the need for careful validation of computational settings before performing extensive density functional theory investigations of surface properties and nanoparticle structures of metals.
Collapse
Affiliation(s)
- Lara Kabalan
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, Wales, UK.
| | - Igor Kowalec
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, Wales, UK.
| | - C Richard A Catlow
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, Wales, UK. and Department of Chemistry, University College London, 20 Gordon Street, London, WC1E 6BT, UK and UK Catalysis Hub, Research Complex at Harwell, Rutherford Appleton Laboratory, Didcot, OX11 OFA, UK
| | - Andrew J Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, Wales, UK.
| |
Collapse
|
11
|
Powell AD, Kroes GJ, Doblhoff-Dier K. Quantum Monte Carlo calculations on dissociative chemisorption of H2 + Al(110): Minimum barrier heights and their comparison to DFT values. J Chem Phys 2020; 153:224701. [DOI: 10.1063/5.0022919] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrew D. Powell
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Katharina Doblhoff-Dier
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| |
Collapse
|
12
|
Castro-Latorre P, Miranda-Rojas S, Mendizabal F. Theoretical exploration of the forces governing the interaction between gold-phthalocyanine and gold surface clusters. RSC Adv 2020; 10:3895-3901. [PMID: 35492636 PMCID: PMC9049280 DOI: 10.1039/c9ra07959a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Accepted: 01/17/2020] [Indexed: 12/03/2022] Open
Abstract
Here we aim to explore the nature of the forces governing the adsorption of gold-phthalocyanine on gold substrates. For this, we designed computational models of metal-free phthalocyanine and gold-phthalocyanine deposited over a gold metallic surface represented by cluster models of different sizes and geometries. Thereby, we were able to determine the role of the metal center and of the size of the substrate in the interaction process. For this purpose, we worked within the framework provided by density functional theory, were the inclusion of the semi-empirical correction of the dispersion forces of Grimme's group was indispensable. It has been shown that the interaction between molecules and surfaces is ruled by van der Waals attractive forces, which determine the stabilization of the studied systems and their geometric properties. Their contribution was characterized by energy decomposition analysis and through the visualization of the dispersion interactions by means of the NCI methodology. Moreover, calculations of Density of States (DOS) showed that the molecule-surface system displays a metal-organic interface evidenced by changes in their electronic structure, in agreement with a charge transfer process found to take place between the interacting parts.
Collapse
Affiliation(s)
- Pablo Castro-Latorre
- Departamento de Química, Facultad de Ciencias, Universidad de Chile Casilla 653 Santiago Chile
| | - Sebastián Miranda-Rojas
- Departamento de Ciencias Químicas, Facultad de Ciencias Exactas, Universidad Andres Bello Avenida República 275 Santiago Chile
| | - Fernando Mendizabal
- Departamento de Química, Facultad de Ciencias, Universidad de Chile Casilla 653 Santiago Chile
| |
Collapse
|
13
|
Smeets EF, Voss J, Kroes GJ. Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H 2 + Cu(111) and H 2 + Ag(111). J Phys Chem A 2019; 123:5395-5406. [PMID: 31149824 PMCID: PMC6600505 DOI: 10.1021/acs.jpca.9b02914] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2019] [Revised: 05/09/2019] [Indexed: 01/08/2023]
Abstract
Specific reaction parameter density functionals (SRP-DFs), which can describe dissociative chemisorption reactions on metals to within chemical accuracy, have so far been based on exchange functionals within the generalized gradient approximation (GGA) and on GGA correlation functionals or van der Waals correlation functionals. These functionals are capable of describing the molecule-metal surface interaction accurately, but they suffer from the general GGA problem that this can be done only at the cost of a rather poor description of the metal. Here, we show that it is possible also to construct SRP-DFs for H2 dissociation on Cu(111) based on meta-GGA functionals, introducing three new functionals based on the "made-simple" (MS) concept. The exchange parts of the three functionals (MS-PBEl, MS-B86bl, and MS-RPBEl) are based on the expressions for the PBE, B86b, and RPBE exchange functionals. Quasi-classical trajectory (QCT) calculations performed with potential energy surfaces (PESs) obtained with the three MS functionals reproduce molecular beam experiments on H2, D2 + Cu(111) with chemical accuracy. Therefore, these three non-empirical functionals themselves are also capable of describing H2 dissociation on Cu(111) with chemical accuracy. Similarly, QCT calculations performed on the MS-PBEl and MS-B86bl PESs reproduced molecular beam and associative desorption experiments on D2, H2 + Ag(111) more accurately than was possible with the SRP48 density functional for H2 + Cu(111). Also, the three new MS functionals describe the Cu, Ag, Au, and Pt metals more accurately than the all-purpose Perdew-Burke-Ernzerhof (PBE) functional. The only disadvantage we noted of the new MS functionals is that, as found for the example of H2 + Cu(111), the reaction barrier height obtained by taking weighted averages of the MS-PBEl and MS-RPBEl functionals is tunable over a smaller range (9 kJ/mol) than possible with the standard GGA PBE and RPBE functionals (33 kJ/mol). As a result of this restricted tunability, it is not possible to construct an SRP-DF for H2 + Ag(111) on the basis of the three examined MS meta-GGA functionals.
Collapse
Affiliation(s)
- Egidius
W. F. Smeets
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Johannes Voss
- SLAC
National Accelerator Laboratory, SUNCAT Center Interface Science &
Catalysis, 2575 Sand
Hill Rd, Menlo Park, California 94025, United States
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
14
|
Mahlberg D, Sakong S, Forster-Tonigold K, Groß A. Improved DFT Adsorption Energies with Semiempirical Dispersion Corrections. J Chem Theory Comput 2019; 15:3250-3259. [DOI: 10.1021/acs.jctc.9b00035] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- David Mahlberg
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
| | - Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
| | - Katrin Forster-Tonigold
- Helmholtz Institute Ulm (HIU) for Electrochemical Energy Storage, 89069 Ulm, Germany
- Karlsruhe Institute of Technology (KIT), P.O. Box
3640, 76021 Karlsruhe, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
- Helmholtz Institute Ulm (HIU) for Electrochemical Energy Storage, 89069 Ulm, Germany
| |
Collapse
|
15
|
Jana S, Patra A, Samal P. Assessing the performance of the Tao-Mo semilocal density functional in the projector-augmented-wave method. J Chem Phys 2018; 149:044120. [DOI: 10.1063/1.5040786] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Abhilash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|