1
|
Oudot B, Doblhoff-Dier K. Reaction barriers at metal surfaces computed using the random phase approximation: Can we beat DFT in the generalized gradient approximation? J Chem Phys 2024; 161:054708. [PMID: 39092949 DOI: 10.1063/5.0220465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2024] [Accepted: 07/14/2024] [Indexed: 08/04/2024] Open
Abstract
Reaction barriers for molecules dissociating on metal surfaces (as relevant to heterogeneous catalysis) are often difficult to predict accurately with density functional theory (DFT). Although the results obtained for several dissociative chemisorption reactions via DFT in the generalized gradient approximation (GGA), in meta-GGA, and for GGA exchange + van der Waals correlation scatter around the true reaction barrier, there is an entire class of dissociative chemisorption reactions for which GGA-type functionals collectively underestimate the reaction barrier. Little is known why GGA-DFT collectively fails in some cases and not in others, and we do not know whether other methods suffer from the same inconsistency. Here, we present barrier heights for dissociative chemisorption reactions obtained from the random phase approximation in the adiabatic-connection fluctuation-dissipation theorem (ACFDT-RPA) and from hybrid functionals with different amounts of exact exchange. By comparing the results obtained for the dissociative chemisorption reaction of H2 on Al(110) (where GGA-DFT collectively underestimates the barrier) and H2 on Cu(111) (where GGA-DFT scatters around the true barrier), we can gauge whether the inconsistent description of the systems persists for hybrid functionals and ACFDT-RPA. We find hybrid functionals to improve the relative description of the two systems, but to fall short of chemical accuracy. ACFDT-RPA improves the results further and leads to chemically accurate barriers for both systems. Together with an analysis of the density of states and the results from selected GGA, meta-GGA, and GGA exchange + van der Waals correlation functionals, these results allow us to discuss possible origins for the inconsistent behavior of GGA-based functionals for molecule-metal reaction barriers.
Collapse
Affiliation(s)
- B Oudot
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, Leiden 2300 RA, The Netherlands
| | - K Doblhoff-Dier
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, Leiden 2300 RA, The Netherlands
| |
Collapse
|
2
|
Sah MK, Naskar K, Adhikari S, Smits B, Meyer J, Somers MF. On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K. J Chem Phys 2024; 161:014306. [PMID: 38953445 DOI: 10.1063/5.0217639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2024] [Accepted: 06/13/2024] [Indexed: 07/04/2024] Open
Abstract
We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.
Collapse
Affiliation(s)
- Mantu Kumar Sah
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Bauke Smits
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
3
|
Sahoo SJ, Xu Q, Lei X, Staros D, Iyer GR, Rubenstein B, Suryanarayana P, Medford AJ. Self-Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces. Chemphyschem 2024; 25:e202300688. [PMID: 38421371 DOI: 10.1002/cphc.202300688] [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: 09/22/2023] [Revised: 02/23/2024] [Accepted: 02/25/2024] [Indexed: 03/02/2024]
Abstract
The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as "Jacob's ladder." This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBE α ${\alpha }$ framework with α ${\alpha }$ being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.
Collapse
Affiliation(s)
| | - Qimen Xu
- Georgia Institute of Technology, Atlanta, GA
- National Supercomputing Center, Shenzhen, People's Republic of China
| | | | - Daniel Staros
- Department of Chemistry, Brown University, Providence, RI
| | - Gopal R Iyer
- Department of Chemistry, Brown University, Providence, RI
| | | | | | | |
Collapse
|
4
|
Powell A, Gerrits N, Tchakoua T, Somers MF, Busnengo HF, Meyer J, Kroes GJ, Doblhoff-Dier K. Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals. J Phys Chem Lett 2024; 15:307-315. [PMID: 38169287 PMCID: PMC10788952 DOI: 10.1021/acs.jpclett.3c02972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 12/21/2023] [Accepted: 12/28/2023] [Indexed: 01/05/2024]
Abstract
Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.
Collapse
Affiliation(s)
- Andrew
D. Powell
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Theophile Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Mark F. Somers
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Heriberto F. Busnengo
- Instituto
de Física Rosario (IFIR), CONICET-UNR, 2000 Rosario, Argentina
- Facultad
de Ciencias Exatas, Ingeniería y
Agrimensura, UNR, 2000 Rosario, Argentina
| | - Jörg Meyer
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| |
Collapse
|
5
|
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
|
6
|
Tchakoua T, Powell AD, Gerrits N, Somers MF, Doblhoff-Dier K, Busnengo HF, Kroes GJ. Simulating Highly Activated Sticking of H 2 on Al(110): Quantum versus Quasi-Classical Dynamics. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:5395-5407. [PMID: 36998253 PMCID: PMC10041643 DOI: 10.1021/acs.jpcc.3c00426] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 02/28/2023] [Indexed: 06/19/2023]
Abstract
We evaluate the importance of quantum effects on the sticking of H2 on Al(110) for conditions that are close to those of molecular beam experiments that have been done on this system. Calculations with the quasi-classical trajectory (QCT) method and with quantum dynamics (QD) are performed using a model in which only motion in the six molecular degrees of freedom is allowed. The potential energy surface used has a minimum barrier height close to the value recently obtained with the quantum Monte Carlo method. Monte Carlo averaging over the initial rovibrational states allowed the QD calculations to be done with an order of magnitude smaller computational expense. The sticking probability curve computed with QD is shifted to lower energies relative to the QCT curve by 0.21 to 0.05 kcal/mol, with the highest shift obtained for the lowest incidence energy. Quantum effects are therefore expected to play a small role in calculations that would evaluate the accuracy of electronic structure methods for determining the minimum barrier height to dissociative chemisorption for H2 + Al(110) on the basis of the standard procedure for comparing results of theory with molecular beam experiments.
Collapse
Affiliation(s)
- Theophile Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Andrew D. Powell
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, 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
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Heriberto F. Busnengo
- Instituto
de Física Rosario (IFIR), CONICET-UNR, Bv. 27 de Febrero 210 bis, 2000 Rosario, Argentina
- Facultad
de Ciencias Exactas, Ingeniería y
Agrimensura, UNR, Av.
Pellegrini 250, 2000 Rosario, Argentina
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
7
|
Hoggan PE. Quantum Monte Carlo method describing supported metal catalysis: Ni(111)/alumina decomposing methane as a route to hydrogen. ADVANCES IN QUANTUM CHEMISTRY 2023. [DOI: 10.1016/bs.aiq.2023.01.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/09/2023]
|
8
|
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
|
9
|
Gao Z, Cheng X, Zhao M. Ab initio diffusion quantum Monte Carlo study of the structural and electronic properties of small Lithium-Chloride LinCl (0,1+) (n = 1–7) clusters. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113781] [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]
|
10
|
Iyer GR, Rubenstein BM. Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry. J Phys Chem A 2022; 126:4636-4646. [PMID: 35820033 DOI: 10.1021/acs.jpca.2c01957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H2O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.
Collapse
Affiliation(s)
- Gopal R Iyer
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
11
|
Sundararaman R, Vigil-Fowler D, Schwarz K. Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface. Chem Rev 2022; 122:10651-10674. [PMID: 35522135 PMCID: PMC10127457 DOI: 10.1021/acs.chemrev.1c00800] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.
Collapse
Affiliation(s)
- Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th Street, Troy, New York 12180, United States
| | - Derek Vigil-Fowler
- Materials, Chemical, and Computational Science Directorate, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
12
|
Auerbach DJ, Tully JC, Wodtke AM. Chemical dynamics from the gas‐phase to surfaces. ACTA ACUST UNITED AC 2021. [DOI: 10.1002/ntls.10005] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Daniel J. Auerbach
- Institut für physikalische Chemie Georg‐August Universität Göttingen Göttingen Germany
- Abteilung für Dynamik an Oberflächen Max‐Planck‐Institut für biophysikalische Chemie Göttingen Germany
| | - John C. Tully
- Department of Chemistry Yale University New Haven Connecticut USA
| | - Alec M. Wodtke
- Institut für physikalische Chemie Georg‐August Universität Göttingen Göttingen Germany
- Abteilung für Dynamik an Oberflächen Max‐Planck‐Institut für biophysikalische Chemie Göttingen Germany
| |
Collapse
|
13
|
Morales-Silva MA, Jordan KD, Shulenburger L, Wagner LK. Frontiers of stochastic electronic structure calculations. J Chem Phys 2021; 154:170401. [PMID: 34241059 DOI: 10.1063/5.0053674] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In recent years there has been a rapid growth in the development and application of new stochastic methods in electronic structure. These methods are quite diverse, from many-body wave function techniques in real space or determinant space to being used to sum perturbative expansions. This growth has been spurred by the more favorable scaling with the number of electrons and often better parallelization over large numbers of central processing unit (CPU) cores or graphical processing units (GPUs) than for high-end non-stochastic wave function based methods. This special issue of the Journal of Chemical Physics includes 33 papers that describe recent developments and applications in this area. As seen from the articles in the issue, stochastic electronic structure methods are applicable to both molecules and solids and can accurately describe systems with strong electron correlation. This issue was motivated, in part, by the 2019 Telluride Science Research Center workshop on Stochastic Electronic Structure Methods that we organized. Below we briefly describe each of the papers in the special issue, dividing the papers into six subtopics.
Collapse
Affiliation(s)
- Miguel A Morales-Silva
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Kenneth D Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Luke Shulenburger
- HEDP Theory Department, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Lucas K Wagner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
14
|
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: 33] [Impact Index Per Article: 11.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
|
15
|
Gerrits N, Smeets EWF, Vuckovic S, Powell AD, Doblhoff-Dier K, Kroes GJ. Density Functional Theory for Molecule-Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not. J Phys Chem Lett 2020; 11:10552-10560. [PMID: 33295770 PMCID: PMC7751010 DOI: 10.1021/acs.jpclett.0c02452] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
While density functional theory (DFT) is perhaps the most used electronic structure theory in chemistry, many of its practical aspects remain poorly understood. For instance, DFT at the generalized gradient approximation (GGA) tends to fail miserably at describing gas-phase reaction barriers, while it performs surprisingly well for many molecule-metal surface reactions. GGA-DFT also fails for many systems in the latter category, and up to now it has not been clear when one may expect it to work. We show that GGA-DFT tends to work if the difference between the work function of the metal and the molecule's electron affinity is greater than ∼7 eV and to fail if this difference is smaller, with sticking of O2 on Al(111) being a spectacular example. Using dynamics calculations we show that, for this system, the DFT problem may be solved as done for gas-phase reactions, i.e., by resorting to hybrid functionals, but using screening at long-range to obtain a correct description of the metal. Our results suggest the GGA error in the O2 + Al(111) barrier height to be functional driven. Our results also suggest the possibility to compute potential energy surfaces for the difficult-to-treat systems with computationally cheap nonself-consistent calculations in which a hybrid functional is applied to a GGA density.
Collapse
Affiliation(s)
- Nick Gerrits
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Egidius W. F. Smeets
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Stefan Vuckovic
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Andrew D. Powell
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|