1
|
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
|
2
|
Stark W, Westermayr J, Douglas-Gallardo OA, Gardner J, Habershon S, Maurer RJ. Machine Learning Interatomic Potentials for Reactive Hydrogen Dynamics at Metal Surfaces Based on Iterative Refinement of Reaction Probabilities. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:24168-24182. [PMID: 38148847 PMCID: PMC10749455 DOI: 10.1021/acs.jpcc.3c06648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 11/12/2023] [Accepted: 11/15/2023] [Indexed: 12/28/2023]
Abstract
The reactive chemistry of molecular hydrogen at surfaces, notably dissociative sticking and hydrogen evolution, plays a crucial role in energy storage and fuel cells. Theoretical studies can help to decipher underlying mechanisms and reaction design, but studying dynamics at surfaces is computationally challenging due to the complex electronic structure at interfaces and the high sensitivity of dynamics to reaction barriers. In addition, ab initio molecular dynamics, based on density functional theory, is too computationally demanding to accurately predict reactive sticking or desorption probabilities, as it requires averaging over tens of thousands of initial conditions. High-dimensional machine learning-based interatomic potentials are starting to be more commonly used in gas-surface dynamics, yet robust approaches to generate reliable training data and assess how model uncertainty affects the prediction of dynamic observables are not well established. Here, we employ ensemble learning to adaptively generate training data while assessing model performance with full uncertainty quantification (UQ) for reaction probabilities of hydrogen scattering on different copper facets. We use this approach to investigate the performance of two message-passing neural networks, SchNet and PaiNN. Ensemble-based UQ and iterative refinement allow us to expose the shortcomings of the invariant pairwise-distance-based feature representation in the SchNet model for gas-surface dynamics.
Collapse
Affiliation(s)
- Wojciech
G. Stark
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Julia Westermayr
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | | | - James Gardner
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Scott Habershon
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Reinhard J. Maurer
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
- Department
of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| |
Collapse
|
3
|
Wagner J, Grabnic T, Sibener SJ. STM Visualization of N 2 Dissociative Chemisorption on Ru(0001) at High Impinging Kinetic Energies. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:18333-18342. [PMID: 36366757 PMCID: PMC9639351 DOI: 10.1021/acs.jpcc.2c05770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Revised: 10/06/2022] [Indexed: 06/16/2023]
Abstract
This paper examines the reactive surface dynamics of energy- and angle-selected N2 dissociation on a clean Ru(0001) surface. Presented herein are the first STM images of highly energetic N2 dissociation on terrace sites utilizing a novel UHV instrument that combines a supersonic molecular beam with an in situ STM that is in-line with the molecular beam. Atomically resolved visualization of individual N2 dissociation events elucidates the fundamental reactive dynamics of the N2/Ru(0001) system by providing a detailed understanding of the on-surface dissociation dynamics: the distance and angle between nitrogen atoms from the same dissociated N2 molecule, site specificity and coordination of binding on terrace sites, and the local evolution of surrounding nanoscopic areas. These properties are precisely measured over a range of impinging N2 kinetic energies and angles, revealing previously unattainable information about the energy dissipation channels that govern the reactivity of the system. The experimental results presented in this paper provide insight into the fundamental N2 dissociation mechanism that, in conjunction with ongoing theoretical modeling, will help determine the role of dynamical processes such as energy transfer to surface phonons and nonadiabatic excitation of electron-hole pairs (ehps). These results will not only help uncover the underlying chemistry and physics that give rise to the unique behavior of this activated dissociative chemisorption system but also represent an exciting approach to studying reaction dynamics by pairing the angstrom-level spatiotemporal resolution of an in situ STM with nonequilibrium fluxes of reactive gases generated in a supersonic molecular beam to access highly activated chemical dynamics and observe the results of individual reaction events.
Collapse
|
4
|
Litman Y, Pós ES, Box CL, Martinazzo R, Maurer RJ, Rossi M. Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. II. Benchmarks and applications. J Chem Phys 2022; 156:194107. [PMID: 35597654 DOI: 10.1063/5.0088400] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], we presented the ring-polymer instanton with explicit friction (RPI-EF) method and showed how it can be connected to the ab initio electronic friction formalism. This framework allows for the calculation of tunneling reaction rates that incorporate the quantum nature of the nuclei and certain types of non-adiabatic effects (NAEs) present in metals. In this paper, we analyze the performance of RPI-EF on model potentials and apply it to realistic systems. For a 1D double-well model, we benchmark the method against numerically exact results obtained from multi-layer multi-configuration time-dependent Hartree calculations. We demonstrate that RPI-EF is accurate for medium and high friction strengths and less accurate for extremely low friction values. We also show quantitatively how the inclusion of NAEs lowers the crossover temperature into the deep tunneling regime, reduces the tunneling rates, and, in certain regimes, steers the quantum dynamics by modifying the tunneling pathways. As a showcase of the efficiency of this method, we present a study of hydrogen and deuterium hopping between neighboring interstitial sites in selected bulk metals. The results show that multidimensional vibrational coupling and nuclear quantum effects have a larger impact than NAEs on the tunneling rates of diffusion in metals. Together with Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], these results advance the calculations of dissipative tunneling rates from first principles.
Collapse
Affiliation(s)
- Y Litman
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - E S Pós
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - C L Box
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - R Martinazzo
- Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - R J Maurer
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - M Rossi
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
5
|
Martinazzo R, Burghardt I. Quantum Dynamics with Electronic Friction. PHYSICAL REVIEW LETTERS 2022; 128:206002. [PMID: 35657868 DOI: 10.1103/physrevlett.128.206002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Revised: 12/19/2021] [Accepted: 04/08/2022] [Indexed: 06/15/2023]
Abstract
A theory of electronic friction is developed using the exact factorization of the electronic-nuclear wave function. No assumption is made regarding the electronic bath, which can be made of independent or interacting electrons, and the nuclei are treated quantally. The ensuing equation of motion for the nuclear wave function is a nonlinear Schrödinger equation including a friction term. The resulting friction kernel agrees with a previously derived mixed quantum-classical result by Dou et al., [Phys. Rev. Lett. 119, 046001 (2017)]PRLTAO0031-900710.1103/PhysRevLett.119.046001, except for a pseudomagnetic contribution in the latter that is here removed. More specifically, it is shown that the electron dynamics generally washes out the gauge fields appearing in the adiabatic dynamics. However, these are fully re-established in the typical situation where the electrons respond rapidly on the slow time scale of the nuclear dynamics (Markov limit). Hence, we predict Berry's phase effects to be observable also in the presence of electronic friction. Application to a model vibrational relaxation problem proves that the proposed approach represents a viable way to account for electronic friction in a fully quantum setting for the nuclear dynamics.
Collapse
Affiliation(s)
- Rocco Martinazzo
- Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
- Istituto di Scienze e Tecnologie Chimiche "Giulio Natta", CNR, via Golgi 19, 20133 Milano, Italy
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt/Main, Germany
| |
Collapse
|
6
|
Gardner J, Douglas-Gallardo OA, Stark WG, Westermayr J, Janke SM, Habershon S, Maurer RJ. NQCDynamics.jl: A Julia package for nonadiabatic quantum classical molecular dynamics in the condensed phase. J Chem Phys 2022; 156:174801. [PMID: 35525649 DOI: 10.1063/5.0089436] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in the condensed phase. To facilitate effective development, code sharing, and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package, which provides a framework for established and emerging methods for performing semiclassical and mixed quantum-classical dynamics in the condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods, which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical and mixed quantum-classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the molecular dynamics with electronic friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems.
Collapse
Affiliation(s)
- James Gardner
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Oscar A Douglas-Gallardo
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Wojciech G Stark
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Julia Westermayr
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Svenja M Janke
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Reinhard J Maurer
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
7
|
Zhou X, Meng G, Guo H, Jiang B. First-Principles Insights into Adiabatic and Nonadiabatic Vibrational Energy-Transfer Dynamics during Molecular Scattering from Metal Surfaces: The Importance of Surface Reactivity. J Phys Chem Lett 2022; 13:3450-3461. [PMID: 35412832 DOI: 10.1021/acs.jpclett.2c00593] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Energy transfer is ubiquitous during molecular collisions and reactions at gas-surface interfaces. Of particular importance is vibrational energy transfer because of its relevance to bond forming and breaking. In this Perspective, we review recent first-principles studies on vibrational energy-transfer dynamics during molecular scattering from metal surfaces at the state-to-state level. Taking several representative systems as examples, we highlight the intrinsic correlation between vibrational energy transfer in nonreactive scattering and surface reactivity and how it operates in both electronically adiabatic and nonadiabatic pathways. Adiabatically, the presence of a dissociation barrier softens a bond in the impinging molecule and increases its couplings with other molecular modes and surface phonons. In the meantime, the stronger interaction between the molecule and the surface also changes the electronic structure at the barrier, resulting in an increase of nonadiabatic effects. We further discuss future prospects toward a more quantitative understanding of this important surface dynamical process.
Collapse
Affiliation(s)
- Xueyao Zhou
- Department of Chemical Physics, School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Gang Meng
- Department of Chemical Physics, School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Bin Jiang
- Department of Chemical Physics, School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
8
|
Formulation of temperature dependent effective Hartree potential incorporating quadratic over linear molecular DOFs-surface modes couplings and its effect on quantum dynamics of D2 (v = 0, j = 0)/D2 (v = 0, j = 2) on Cu(111) metal surface. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
Dutta J, Mandal S, Adhikari S, Spiering P, Meyer J, Somers MF. Effect of surface temperature on quantum dynamics of H 2 on Cu(111) using a chemically accurate potential energy surface. J Chem Phys 2021; 154:104103. [PMID: 33722025 DOI: 10.1063/5.0035830] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The effect of surface atom vibrations on H2 scattering from a Cu(111) surface at different temperatures is being investigated for hydrogen molecules in their rovibrational ground state (v = 0, j = 0). We assume weakly correlated interactions between molecular degrees of freedom and surface modes through a Hartree product type wavefunction. While constructing the six-dimensional effective Hamiltonian, we employ (a) a chemically accurate potential energy surface according to the static corrugation model [M. Wijzenbroek and M. F. Somers, J. Chem. Phys. 137, 054703 (2012)]; (b) normal mode frequencies and displacement vectors calculated with different surface atom interaction potentials within a cluster approximation; and (c) initial state distributions for the vibrational modes according to Bose-Einstein probability factors. We carry out 6D quantum dynamics with the so-constructed effective Hamiltonian and analyze sticking and state-to-state scattering probabilities. The surface atom vibrations affect the chemisorption dynamics. The results show physically meaningful trends for both reaction and scattering probabilities compared to experimental and other theoretical results.
Collapse
Affiliation(s)
- Joy Dutta
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Souvik Mandal
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Paul Spiering
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- 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
| |
Collapse
|
13
|
Spiering P, Shakouri K, Behler J, Kroes GJ, Meyer J. Orbital-Dependent Electronic Friction Significantly Affects the Description of Reactive Scattering of N 2 from Ru(0001). J Phys Chem Lett 2019; 10:2957-2962. [PMID: 31088059 PMCID: PMC6558642 DOI: 10.1021/acs.jpclett.9b00523] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2019] [Accepted: 05/15/2019] [Indexed: 05/18/2023]
Abstract
Electron-hole pair (ehp) excitation is thought to substantially affect the dynamics of molecules on metal surfaces, but it is not clear whether this can be better addressed by orbital-dependent friction (ODF) or the local density friction approximation (LDFA). We investigate the effect of ehp excitation on the dissociative chemisorption of N2 on and its inelastic scattering from Ru(0001), which is the benchmark system of highly activated dissociation, with these two different models. ODF is in better agreement with the best experimental estimates for the reaction probabilities than LDFA, yields results for vibrational excitation in better agreement with experiment, but slightly overestimates the translational energy loss during scattering. N2 on Ru(0001) is thus the first system for which the ODF and LDFA approaches are shown to yield substantially different results for easily accessible experimental observables, including reaction probabilities.
Collapse
Affiliation(s)
- Paul Spiering
- Gorlaeus Laberatories, Leiden Institute of Chemistry , Leiden University , P.O. Box 9502, 2300 RA Leiden , The Netherlands
| | - Khosrow Shakouri
- Gorlaeus Laberatories, Leiden Institute of Chemistry , Leiden University , P.O. Box 9502, 2300 RA Leiden , The Netherlands
| | - Jörg Behler
- Universität Göttingen , Institut für Physikalische Chemie, Theoretische Chemie, Tammannstr. 6 , 37077 Göttingen , Germany
| | - Geert-Jan Kroes
- Gorlaeus Laberatories, Leiden Institute of Chemistry , Leiden University , P.O. Box 9502, 2300 RA Leiden , The Netherlands
| | - Jörg Meyer
- Gorlaeus Laberatories, Leiden Institute of Chemistry , Leiden University , P.O. Box 9502, 2300 RA Leiden , The Netherlands
| |
Collapse
|
14
|
Jiang B, Guo H. Dynamics in reactions on metal surfaces: A theoretical perspective. J Chem Phys 2019; 150:180901. [PMID: 31091904 DOI: 10.1063/1.5096869] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.
Collapse
Affiliation(s)
- Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
15
|
Füchsel G, Zhou X, Jiang B, Juaristi JI, Alducin M, Guo H, Kroes GJ. Reactive and Nonreactive Scattering of HCl from Au(111): An Ab Initio Molecular Dynamics Study. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:2287-2299. [PMID: 30740194 PMCID: PMC6366682 DOI: 10.1021/acs.jpcc.8b10686] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 12/19/2018] [Indexed: 05/20/2023]
Abstract
The HCl + Au(111) system has recently become a benchmark for highly activated dissociative chemisorption, which presumably is strongly affected by electron-hole pair excitation. Previous dynamics calculations, which were based on density functional theory at the generalized gradient approximation level (GGA-DFT) for the molecule-surface interaction, have all overestimated measured reaction probabilities by at least an order of magnitude. Here, we perform ab initio molecular dynamics (AIMD) and AIMD with electronic friction (AIMDEF) calculations employing a density functional that includes the attractive van der Waals interaction. Our calculations model the simultaneous and possibly synergistic effects of surface temperature, surface atom motion, electron-hole pair excitation, the molecular beam conditions of the experiments, and the van der Waals interaction on the reactivity. We find that reaction probabilities computed with AIMDEF and the SRP32-vdW functional still overestimate the measured reaction probabilities, by a factor 18 for the highest incidence energy at which measurements were performed (≈2.5 eV). Even granting that the experiment could have underestimated the sticking probability by about a factor three, this still translates into a considerable overestimation of the reactivity by the current theory. Likewise, scaled transition probabilities for vibrational excitation from ν = 1, j = 1 to ν = 2 are overestimated by the AIMDEF theory, by factors 3-8 depending on the initial conditions modeled. Energy losses to the surface and translational energy losses are, however, in good agreement with experimental values.
Collapse
Affiliation(s)
- Gernot Füchsel
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- Institut
für Chemie und Biochemie—Physikalische und Theoretische
Chemie, Freie Universität Berlin, Takustraße3, 14195 Berlin, Germany
- E-mail: (G.F.)
| | - Xueyao Zhou
- Hefei
National Laboratory for Physical Science at the Microscale, Department
of Chemical Physics, School of Chemistry and Materials, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Bin Jiang
- Hefei
National Laboratory for Physical Science at the Microscale, Department
of Chemical Physics, School of Chemistry and Materials, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - J. Iñaki Juaristi
- Departamento
de Física de Materiales, Facultad
de Químicas (UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain
- Centro
de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
- Donostia
International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Maite Alducin
- Centro
de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
- Donostia
International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University
of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- E-mail: . Phone: +31 (0)71 527
4396 (G.-J.K.)
| |
Collapse
|
16
|
Zhang Y, Maurer RJ, Guo H, Jiang B. Hot-electron effects during reactive scattering of H 2 from Ag(111): the interplay between mode-specific electronic friction and the potential energy landscape. Chem Sci 2019; 10:1089-1097. [PMID: 30774906 PMCID: PMC6346630 DOI: 10.1039/c8sc03955k] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Accepted: 11/07/2018] [Indexed: 01/29/2023] Open
Abstract
The breakdown of the Born-Oppenheimer approximation gives rise to nonadiabatic effects in gas-surface reactions at metal surfaces. However, for a given reaction, it remains unclear which factors quantitatively determine whether these effects measurably contribute to surface reactivity in catalysis and photo/electrochemistry. Here, we systematically investigate hot electron effects during H2 scattering from Ag(111) using electronic friction theory. We combine first-principles calculations of tensorial friction by time-dependent perturbation theory based on density functional theory and an analytical neural network representation, to overcome the limitations of existing approximations and explicitly simulate mode-specific nonadiabatic energy loss during molecular dynamics. Despite sizable hot-electron-induced energy loss, no measurable nonadiabatic effects can be found for H2 scattering on Ag(111). This is in stark contrast to previous reports for vibrationally excited H2 scattering on Cu(111). By detailed analysis of the two systems, we attribute this discrepancy to a subtle interplay between the magnitude of electronic friction along intramolecular vibration and the shape of the potential energy landscape that controls the molecular velocity at impact. On the basis of this characterization, we offer guidance for the search of highly nonadiabatic surface reactions.
Collapse
Affiliation(s)
- Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale , Department of Chemical Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China .
| | - Reinhard J Maurer
- Department of Chemistry and Centre for Scientific Computing , University of Warwick , Gibbet Hill Road , Coventry , CV4 7AL , UK .
- Department of Chemistry , Yale University , New Haven , Connecticut 06520 , USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , USA
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale , Department of Chemical Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China .
| |
Collapse
|
17
|
Maurer RJ, Zhang Y, Guo H, Jiang B. Hot electron effects during reactive scattering of H2 from Ag(111): assessing the sensitivity to initial conditions, coupling magnitude, and electronic temperature. Faraday Discuss 2019; 214:105-121. [DOI: 10.1039/c8fd00140e] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We use an analytical representation of electronic friction for H2 on Ag(111) to assess the validity and robustness of the MDEF method based on TDPT.
Collapse
Affiliation(s)
- Reinhard J. Maurer
- Department of Chemistry
- Centre for Scientific Computing
- University of Warwick
- Coventry
- UK
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Hua Guo
- Department of Chemistry and Chemical Biology
- University of New Mexico
- Albuquerque
- USA
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| |
Collapse
|
18
|
Spiering P, Wijzenbroek M, Somers MF. An improved static corrugation model. J Chem Phys 2018; 149:234702. [DOI: 10.1063/1.5058271] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Affiliation(s)
- P. Spiering
- Leiden University, Leiden,
Zuid-holland 2300 RA, The Netherlands
| | | | - M. F. Somers
- Leiden University, Leiden,
Zuid-holland 2300 RA, The Netherlands
| |
Collapse
|
19
|
Shakouri K, Behler J, Meyer J, Kroes GJ. Analysis of Energy Dissipation Channels in a Benchmark System of Activated Dissociation: N 2 on Ru(0001). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2018; 122:23470-23480. [PMID: 30364480 PMCID: PMC6196344 DOI: 10.1021/acs.jpcc.8b06729] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2018] [Revised: 09/18/2018] [Indexed: 05/20/2023]
Abstract
The excitation of electron-hole pairs in reactive scattering of molecules at metal surfaces often affects the physical and dynamical observables of interest, including the reaction probability. Here, we study the influence of electron-hole pair excitation on the dissociative chemisorption of N2 on Ru(0001) using the local density friction approximation method. The effect of surface atom motion has also been taken into account by a high-dimensional neural network potential. Our nonadiabatic molecular dynamics simulations with electronic friction show that the reaction of N2 is more strongly affected by the energy transfer to surface phonons than by the energy loss to electron-hole pairs. The discrepancy between the computed reaction probabilities and experimental results is within the experimental error both with and without friction; however, the incorporation of electron-hole pairs yields somewhat better agreement with experiments, especially at high collision energies. We also calculate the vibrational efficacy for the N2 + Ru(0001) reaction and demonstrate that the N2 reaction is more enhanced by exciting the molecular vibrations than by adding an equivalent amount of energy into translation.
Collapse
Affiliation(s)
- Khosrow Shakouri
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Behler
- Institut
für Physikalische Chemie, Theoretische Chemie, Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Jörg Meyer
- 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
|
20
|
Nour Ghassemi E, Somers M, Kroes GJ. Test of the Transferability of the Specific Reaction Parameter Functional for H 2 + Cu(111) to D 2 + Ag(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2018; 122:22939-22952. [PMID: 30344838 PMCID: PMC6189907 DOI: 10.1021/acs.jpcc.8b05658] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Revised: 09/04/2018] [Indexed: 06/08/2023]
Abstract
The accurate description of the dissociative chemisorption of a molecule on a metal surface requires a chemically accurate description of the molecule-surface interaction. Previously, it was shown that the specific reaction parameter approach to density functional theory (SRP-DFT) enables accurate descriptions of the reaction of dihydrogen with metal surfaces in, for instance, H2 + Pt(111), H2 + Cu(111), and H2 + Cu(100). SRP-DFT likewise allowed a chemically accurate description of dissociation of methane on Ni(111) and Pt(111), and the SRP functional for CH4 + Ni(111) was transferable to CH4 + Pt(111), where Ni and Pt belong to the same group. Here, we investigate whether the SRP density functional derived for H2 + Cu(111) also gives chemically accurate results for H2 + Ag(111), where Ag belongs to the same group as Cu. To do this, we have performed quasi-classical trajectory calculations using the six-dimensional potential energy surface of H2 + Ag(111) within the Born-Oppenheimer static surface approximation. The computed reaction probabilities are compared with both state-resolved associative desorption and molecular beam sticking experiments. Our results do not yet show transferability, as the computed sticking probabilities and initial-state selected reaction probabilities are shifted relative to experiment to higher energies by about 2-3 kcal/mol. The lack of transferability may be due to the different character of the SRP functionals for H2 + Cu and CH4 + group 10 metals, the latter containing a van der Waals correlation functional and the former not.
Collapse
|
21
|
Montemore MM, Hoyt R, Kolesov G, Kaxiras E. Reaction-Induced Excitations and Their Effect on Surface Chemistry. ACS Catal 2018. [DOI: 10.1021/acscatal.8b03266] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Matthew M. Montemore
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Robert Hoyt
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Grigory Kolesov
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Efthimios Kaxiras
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
22
|
Spiering P, Meyer J. Testing Electronic Friction Models: Vibrational De-excitation in Scattering of H 2 and D 2 from Cu(111). J Phys Chem Lett 2018; 9:1803-1808. [PMID: 29528648 PMCID: PMC5890313 DOI: 10.1021/acs.jpclett.7b03182] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 03/12/2018] [Indexed: 05/29/2023]
Abstract
At present, molecular dynamics with electronic friction (MDEF) is the workhorse model to go beyond the Born-Oppenheimer approximation in modeling dynamics of molecules at metal surfaces. Concomitant friction coefficients can be calculated with either the local density friction approximation (LDFA) or orbital-dependent friction (ODF), which, unlike LDFA, accounts for anisotropy while relying on other approximations. Due to the computational cost of ODF, extensive high-dimensional MDEF trajectory calculations of experimentally measurable observables have hitherto only been performed based on LDFA. We overcome this limitation with a continuous neural-network-based representation. In our first application to the scattering of vibrationally excited H2 and D2 from Cu(111), we predict up to three times higher vibrational de-excitation probabilities with ODF than with LDFA. These results indicate that anisotropic electronic friction can be important for specific molecular observables. Future experiments can test for this "fingerprint" of different approximations underlying state-of-the-art MDEF.
Collapse
Affiliation(s)
- Paul Spiering
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
23
|
Dou W, Miao G, Subotnik JE. Born-Oppenheimer Dynamics, Electronic Friction, and the Inclusion of Electron-Electron Interactions. PHYSICAL REVIEW LETTERS 2017; 119:046001. [PMID: 29341745 DOI: 10.1103/physrevlett.119.046001] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2017] [Indexed: 06/07/2023]
Abstract
We present a universal expression for the electronic friction as felt by a set of classical nuclear degrees of freedom (DOFs) coupled to a manifold of quantum electronic DOFs; no assumptions are made regarding the nature of the electronic Hamiltonian and electron-electron repulsions are allowed. Our derivation is based on a quantum-classical Liouville equation for the coupled electronic-nuclear motion, followed by an adiabatic approximation whereby electronic transitions are assumed to equilibrate faster than nuclear movement. The resulting form of friction is completely general, but does reduce to previously published expressions for the quadratic Hamiltonian (i.e., Hamiltonians without electronic correlation). At equilibrium, the second fluctuation-dissipation theorem is satisfied and the frictional matrix is symmetric. To demonstrate the importance of electron-electron correlation, we study electronic friction within the Anderson-Holstein model, where a proper treatment of electron-electron interactions shows signatures of a Kondo resonance and a mean-field treatment is completely inadequate.
Collapse
Affiliation(s)
- Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Gaohan Miao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
24
|
Maurer RJ, Jiang B, Guo H, Tully JC. Mode Specific Electronic Friction in Dissociative Chemisorption on Metal Surfaces: H_{2} on Ag(111). PHYSICAL REVIEW LETTERS 2017; 118:256001. [PMID: 28696728 DOI: 10.1103/physrevlett.118.256001] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Indexed: 05/11/2023]
Abstract
Electronic friction and the ensuing nonadiabatic energy loss play an important role in chemical reaction dynamics at metal surfaces. Using molecular dynamics with electronic friction evaluated on the fly from density functional theory, we find strong mode dependence and a dominance of nonadiabatic energy loss along the bond stretch coordinate for scattering and dissociative chemisorption of H_{2} on the Ag(111) surface. Exemplary trajectories with varying initial conditions indicate that this mode specificity translates into modulated energy loss during a dissociative chemisorption event. Despite minor nonadiabatic energy loss of about 5%, the directionality of friction forces induces dynamical steering that affects individual reaction outcomes, specifically for low-incidence energies and vibrationally excited molecules. Mode-specific friction induces enhanced loss of rovibrational rather than translational energy and will be most visible in its effect on final energy distributions in molecular scattering experiments.
Collapse
Affiliation(s)
- Reinhard J Maurer
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Bin Jiang
- Department of Chemical Physics, University of Science & Technology of China, Hefei, Anhui 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - John C Tully
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
25
|
Rittmeyer SP, Ward DJ, Gütlein P, Ellis J, Allison W, Reuter K. Energy Dissipation during Diffusion at Metal Surfaces: Disentangling the Role of Phonons versus Electron-Hole Pairs. PHYSICAL REVIEW LETTERS 2016; 117:196001. [PMID: 27858423 DOI: 10.1103/physrevlett.117.196001] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Indexed: 06/06/2023]
Abstract
Helium spin echo experiments combined with ab initio based Langevin molecular dynamics simulations are used to quantify the adsorbate-substrate coupling during the thermal diffusion of Na atoms on Cu(111). An analysis of trajectories within the local density friction approximation allows the contribution from electron-hole pair excitations to be separated from the total energy dissipation. Despite the minimal electronic friction coefficient of Na and the relatively small mass mismatch to Cu promoting efficient phononic dissipation, about (20±5)% of the total energy loss is attributable to electronic friction. The results suggest a significant role of electronic nonadiabaticity in the rapid thermalization generally relied upon in adiabatic diffusion theories.
Collapse
Affiliation(s)
- Simon P Rittmeyer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| | - David J Ward
- Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge CB3 0HE, United Kingdom
| | - Patrick Gütlein
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| | - John Ellis
- Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge CB3 0HE, United Kingdom
| | - William Allison
- Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge CB3 0HE, United Kingdom
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, 85747 Garching, Germany
| |
Collapse
|
26
|
Wijzenbroek M, Helstone D, Meyer J, Kroes GJ. Dynamics of H2 dissociation on the close-packed (111) surface of the noblest metal: H2 + Au(111). J Chem Phys 2016; 145:144701. [DOI: 10.1063/1.4964486] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mark Wijzenbroek
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Darcey Helstone
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- 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
|
27
|
Luo X, Jiang B, Juaristi JI, Alducin M, Guo H. Electron-hole pair effects in methane dissociative chemisorption on Ni(111). J Chem Phys 2016; 145:044704. [DOI: 10.1063/1.4959288] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xuan Luo
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Bin Jiang
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - J. Iñaki Juaristi
- Centro de Física de Materiales CFM/MPC(CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center DIPC, P. Manuel de Lardizabal 4, 20018 San Sebastián, Spain
- Departamento de Física de Materiales, Facultad de Químicas, Universidad del País Vasco (UPV/EHU), Apartado 1072, 20080 San Sebastián, Spain
| | - Maite Alducin
- Centro de Física de Materiales CFM/MPC(CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center DIPC, P. Manuel de Lardizabal 4, 20018 San Sebastián, Spain
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
28
|
Jiang B, Alducin M, Guo H. Electron-Hole Pair Effects in Polyatomic Dissociative Chemisorption: Water on Ni(111). J Phys Chem Lett 2016; 7:327-31. [PMID: 26732612 DOI: 10.1021/acs.jpclett.5b02737] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
The influence of electron-hole pairs in dissociative chemisorption of a polyatomic molecule (water) on metal surfaces is assessed for the first time using a friction approach. The atomic local density dependent friction coefficients computed based on a free electron gas embedding model are employed in classical molecular dynamics simulations of the water dissociation dynamics on rigid Ni(111) using a recently developed nine dimensional interaction potential energy surface for the system. The results indicate that nonadiabatic effects are relatively small and they do not qualitatively alter the mode specificity in the dissociation.
Collapse
Affiliation(s)
- Bin Jiang
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
- Department of Chemical Physics, University of Science and Technology of China , Hefei, Anhui 230026, China
| | - Maite Alducin
- Centro de Física de Materiales Centro Mixto, CFM/MPC (CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center DIPC, P. Manuel de Lardizabal 4, 20018 San Sebastián, Spain
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
| |
Collapse
|
29
|
Kroes GJ, Díaz C. Quantum and classical dynamics of reactive scattering of H2 from metal surfaces. Chem Soc Rev 2016; 45:3658-700. [DOI: 10.1039/c5cs00336a] [Citation(s) in RCA: 120] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
State-of-the-art theoretical models allow nowadays an accurate description of H2/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry
- Gorlaeus Laboratories
- Leiden University
- 2300 RA Leiden
- The Netherlands
| | - Cristina Díaz
- Departamento de Química
- Módulo 13
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
| |
Collapse
|
30
|
Diesing D, Hasselbrink E. Chemical energy dissipation at surfaces under UHV and high pressure conditions studied using metal–insulator–metal and similar devices. Chem Soc Rev 2016; 45:3747-55. [DOI: 10.1039/c5cs00932d] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thin film metal heterostructures have allowed new light to be shed on the dissipation of chemical energy into electric excitations on metal surfaces.
Collapse
Affiliation(s)
- Detlef Diesing
- Fakultät f. Chemie
- Universität Duisburg-Essen
- 45141 Essen
- Germany
| | | |
Collapse
|
31
|
Kroes GJ. Toward a Database of Chemically Accurate Barrier Heights for Reactions of Molecules with Metal Surfaces. J Phys Chem Lett 2015; 6:4106-14. [PMID: 26722785 DOI: 10.1021/acs.jpclett.5b01344] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Being able to calculate reaction barrier heights to within chemical accuracy (errors < 1 kcal/mol) is crucial to the accurate modeling of chemical reactions. Although accurate databases exist that can help theorists with benchmarking new electronic structure theories on gas-phase chemical reactions, no such databases exist for reactions of molecules with metal surfaces. Nonetheless, most chemicals are made in heterogeneously catalyzed processes, of which many take place over metal particles. Presently, barrier heights for molecule-metal surface reactions have been determined with chemical accuracy for only two systems, that is, H2 + Cu(111) and H2 + Cu(100). This has been done with semiempirically determined density functionals, which were fitted through comparisons of dynamics results with molecular beam sticking probabilities. The prospects of extending the database with chemically accurate data for other molecule-metal reactions, either with the use of semiempirical density functional theory or with first-principles theory, are discussed.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
32
|
Wijzenbroek M, Klein DM, Smits B, Somers MF, Kroes GJ. Performance of a Non-Local van der Waals Density Functional on the Dissociation of H2 on Metal Surfaces. J Phys Chem A 2015; 119:12146-58. [DOI: 10.1021/acs.jpca.5b06008] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mark Wijzenbroek
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA, Leiden, The Netherlands
| | - David M. Klein
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA, Leiden, The Netherlands
| | - Bauke Smits
- 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
|
33
|
Öberg H, Gladh J, Marks K, Ogasawara H, Nilsson A, Pettersson LGM, Öström H. Indication of non-thermal contribution to visible femtosecond laser-induced CO oxidation on Ru(0001). J Chem Phys 2015; 143:074701. [PMID: 26298142 DOI: 10.1063/1.4928646] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
We studied CO oxidation on Ru(0001) induced by 400 nm and 800 nm femtosecond laser pulses where we find a branching ratio between CO oxidation and desorption of 1:9 and 1:31, respectively, showing higher selectivity towards CO oxidation for the shorter wavelength excitation. Activation energies computed with density functional theory show discrepancies with values extracted from the experiments, indicating both a mixture between different adsorbed phases and importance of non-adiabatic effects on the effective barrier for oxidation. We simulated the reactions using kinetic modeling based on the two-temperature model of laser-induced energy transfer in the substrate combined with a friction model for the coupling to adsorbate vibrations. This model gives an overall good agreement with experiment except for the substantial difference in yield ratio between CO oxidation and desorption at 400 nm and 800 nm. However, including also the initial, non-thermal effect of electrons transiently excited into antibonding states of the O-Ru bond yielded good agreement with all experimental results.
Collapse
Affiliation(s)
- H Öberg
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - J Gladh
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - K Marks
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - H Ogasawara
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - A Nilsson
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - L G M Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - H Öström
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
34
|
Rittmeyer SP, Meyer J, Juaristi JI, Reuter K. Electronic Friction-Based Vibrational Lifetimes of Molecular Adsorbates: Beyond the Independent-Atom Approximation. PHYSICAL REVIEW LETTERS 2015; 115:046102. [PMID: 26252696 DOI: 10.1103/physrevlett.115.046102] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2014] [Indexed: 05/11/2023]
Abstract
We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) type charge partitioning scheme accounts for intramolecular contributions and overcomes the systematic underestimation of the nonadiabatic losses obtained within the prevalent independent-atom approximation. The quantitative agreement obtained with theoretical and experimental benchmark data suggests the LDFA-AIM scheme as an efficient and reliable approach to account for electronic dissipation in ab initio molecular dynamics simulations of surface chemical reactions.
Collapse
Affiliation(s)
- Simon P Rittmeyer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - J Iñaki Juaristi
- Departamento de Física de Materiales, Facultad de Químicas, UPV/EHU, Apartado 1072, 20080 San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 5, 20018 San Sebastián, Spain
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| |
Collapse
|
35
|
Mandal S, Sahoo T, Ghosh S, Adhikari S. The effect of phonon modes and electron–hole pair couplings on molecule–surface scattering processes. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2015. [DOI: 10.1142/s0219633615500285] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The effect of phonon modes and electron–hole pair (elhp) couplings at different surface temperature on D 2(v = 0, 1; j = 0)– Cu (111) collision has been explored by assuming weakly correlated interactions between molecular Degrees of Freedoms (DOFs) with surface modes and elhp excitations through a Hartree product type wavefunction, where the initial state distributions for the phonon modes and the elhp couplings are incorporated by using Bose–Einstein and Fermi–Dirac probability factors, respectively. We carry out four (4D⊗2D)- and six (6D)- dimensional quantum dynamics on such an effective Hamiltonian, and depict the calculated sticking/transition probabilities and energy transfer from molecule to the surface. The phonon modes slightly affect the sticking probability by broadening the profile, but the transition probability are substantially changed with respect to the rigid surface. On the contrary, the inclusion of elhp coupling along with phonon modes does not change the results much compared to the only phonon case.
Collapse
Affiliation(s)
- Souvik Mandal
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Tapas Sahoo
- Department of Chemical Physics, Wiseman Institute of Science, Rehovot 76100, Israel
| | - Sandip Ghosh
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Satrajit Adhikari
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| |
Collapse
|
36
|
Saalfrank P, Juaristi JI, Alducin M, Blanco-Rey M, Muiño RD. Vibrational lifetimes of hydrogen on lead films: An ab initio molecular dynamics with electronic friction (AIMDEF) study. J Chem Phys 2014; 141:234702. [DOI: 10.1063/1.4903309] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
37
|
Nattino F, Genova A, Guijt M, Muzas AS, Díaz C, Auerbach DJ, Kroes GJ. Dissociation and recombination of D2 on Cu(111): Ab initio molecular dynamics calculations and improved analysis of desorption experiments. J Chem Phys 2014; 141:124705. [DOI: 10.1063/1.4896058] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Francesco Nattino
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Alessandro Genova
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Marieke Guijt
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Alberto S. Muzas
- Departamento de Química Módulo 13, Universitad Autónoma de Madrid, 28049 Madrid, Spain
| | - Cristina Díaz
- Departamento de Química Módulo 13, Universitad Autónoma de Madrid, 28049 Madrid, Spain
| | - Daniel J. Auerbach
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
38
|
Wijzenbroek M, Kroes GJ. The effect of the exchange-correlation functional on H2 dissociation on Ru(0001). J Chem Phys 2014; 140:084702. [DOI: 10.1063/1.4865946] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Boereboom JM, Wijzenbroek M, Somers MF, Kroes GJ. Towards a specific reaction parameter density functional for reactive scattering of H2 from Pd(111). J Chem Phys 2013; 139:244707. [DOI: 10.1063/1.4851355] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Mondal A, Wijzenbroek M, Bonfanti M, Díaz C, Kroes GJ. Thermal Lattice Expansion Effect on Reactive Scattering of H2 from Cu(111) at Ts = 925 K. J Phys Chem A 2013; 117:8770-81. [DOI: 10.1021/jp4042183] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Arobendo Mondal
- Indian Institute of Science Education and Research Kolkata (IISER-K),
741252, Nadia, West Bengal, India
| | - Mark Wijzenbroek
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, Post Office Box 9502, 2300 RA Leiden, The Netherlands
| | - Matteo Bonfanti
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, Post Office Box 9502, 2300 RA Leiden, The Netherlands
| | - Cristina Díaz
- Departamento de Química, Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, Post Office Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
41
|
Tremblay JC. A unifying model for non-adiabatic coupling at metallic surfaces beyond the local harmonic approximation: From vibrational relaxation to scanning tunneling microscopy. J Chem Phys 2013; 138:244106. [DOI: 10.1063/1.4811150] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
|
42
|
Füchsel G, Schimka S, Saalfrank P. On the Role of Electronic Friction for Dissociative Adsorption and Scattering of Hydrogen Molecules at a Ru(0001) Surface. J Phys Chem A 2013; 117:8761-9. [DOI: 10.1021/jp403860p] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Gernot Füchsel
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Straße 24-25, D-14476
Potsdam-Golm, Germany
| | - Selina Schimka
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Straße 24-25, D-14476
Potsdam-Golm, Germany
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Straße 24-25, D-14476
Potsdam-Golm, Germany
| |
Collapse
|
43
|
Schindler B, Diesing D, Hasselbrink E. Electronic Excitations in the Course of the Reaction of H with Coinage and Noble Metal Surfaces: A Comparison. Z PHYS CHEM 2013. [DOI: 10.1524/zpch.2013.0408] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
44
|
|
45
|
Sementa L, Wijzenbroek M, van Kolck BJ, Somers MF, Al-Halabi A, Busnengo HF, Olsen RA, Kroes GJ, Rutkowski M, Thewes C, Kleimeier NF, Zacharias H. Reactive scattering of H2 from Cu(100): Comparison of dynamics calculations based on the specific reaction parameter approach to density functional theory with experiment. J Chem Phys 2013; 138:044708. [DOI: 10.1063/1.4776224] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|
46
|
|
47
|
Non-adiabatic Effects at Surfaces Simulated with TDDFT Molecular Dynamics. DYNAMICS OF GAS-SURFACE INTERACTIONS 2013. [DOI: 10.1007/978-3-642-32955-5_12] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
|
48
|
Kroes GJ. Towards chemically accurate simulation of molecule-surface reactions. Phys Chem Chem Phys 2012; 14:14966-81. [PMID: 23037951 DOI: 10.1039/c2cp42471a] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This perspective addresses four challenges facing theorists whose aim is to make quantitatively accurate predictions for reactions of molecules on metal surfaces, and suggests ways of meeting these challenges, focusing on dissociative chemisorption reactions of H(2), N(2), and CH(4). Addressing these challenges is ultimately of practical importance to a more accurate description of overall heterogeneously catalysed reactions, which play a role in the production of more than 90% of man-made chemicals. One challenge is to describe the interaction of a molecule with a metal surface with chemical accuracy, i.e., with errors in reaction barrier heights less than 1 kcal mol(-1). In this framework, the potential of a new implementation of specific reaction parameter density functional theory (SRP-DFT) will be discussed, with emphasis on applications to reaction of H(2) with metal surfaces. Two additional challenges are to come up with improved descriptions of the effects of phonons and electron-hole pairs on reaction of molecules like N(2) on metal surfaces. Phonons can be tackled using sudden approximations in quantum dynamics, and through Ab Initio Molecular Dynamics (AIMD) calculations using classical dynamics. To additionally achieve an accurate description of the effect of electron-hole pair excitation on dissociative chemisorption within a classical dynamics framework, it may be possible to combine AIMD with electronic friction. The fourth challenge we will consider is how to achieve an accurate quantum mechanical description of the dissociative chemisorption of a polyatomic molecule, like methane, on a metal surface. A method of potential interest is the Multi-Configuration Time-Dependent Hartree (MCTDH) method.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
49
|
Muzas AS, Juaristi JI, Alducin M, Muiño RD, Kroes GJ, Díaz C. Vibrational deexcitation and rotational excitation of H2 and D2 scattered from Cu(111): Adiabatic versus non-adiabatic dynamics. J Chem Phys 2012; 137:064707. [DOI: 10.1063/1.4742907] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Martin-Gondre L, Bocan G, Alducin M, Juaristi J, Díez Muiño R. Energy dissipation channels in the adsorption of N on Ag(111). COMPUT THEOR CHEM 2012. [DOI: 10.1016/j.comptc.2012.03.009] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|