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
|
Gu K, Lin S. Advances in the Dynamics of Adsorbate Diffusion on Metal Surfaces: Focus on Hydrogen and Oxygen. Chemphyschem 2024; 25:e202400083. [PMID: 38511509 DOI: 10.1002/cphc.202400083] [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: 01/28/2024] [Revised: 02/26/2024] [Accepted: 03/21/2024] [Indexed: 03/22/2024]
Abstract
Adsorbates on metal surfaces are typically formed from the dissociative chemisorption of molecules occurring at gas-solid interfaces. These adsorbed species exhibit unique diffusion behaviors on metal surfaces, which are influenced by their translational energy. They play crucial roles in various fields, including heterogeneous catalysis and corrosion. This review examines recent theoretical advancements in understanding the diffusion dynamics of adsorbates on metal surfaces, with a specific emphasis on hydrogen and oxygen atoms. The diffusion processes of adsorbates on metal surfaces involve two energy transfer mechanisms: surface phonons and electron-hole pair excitations. This review also surveys new theoretical methods, including the characterization of the electron-hole pair excitation within electronic friction models, the acceleration of quantum chemistry calculations through machine learning, and the treatment of atomic nuclear motion from both quantum mechanical and classical perspectives. Furthermore, this review offers valuable insights into how energy transfer, nuclear quantum effects, supercell sizes, and the topography of potential energy surfaces impact the diffusion behavior of hydrogen and oxygen species on metal surfaces. Lastly, some preliminary research proposals are presented.
Collapse
Affiliation(s)
- Kaixuan Gu
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou, 350002, China
| | - Sen Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou, 350002, China
| |
Collapse
|
3
|
Liu XY, Wang SR, Fang WH, Cui G. Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation. J Chem Theory Comput 2024; 20:3426-3439. [PMID: 38656202 DOI: 10.1021/acs.jctc.4c00068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Sheng-Rui Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
4
|
Gerrits N, Jackson B, Bogaerts A. Accurate Reaction Probabilities for Translational Energies on Both Sides of the Barrier of Dissociative Chemisorption on Metal Surfaces. J Phys Chem Lett 2024; 15:2566-2572. [PMID: 38416779 PMCID: PMC10926167 DOI: 10.1021/acs.jpclett.3c03408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Revised: 02/16/2024] [Accepted: 02/26/2024] [Indexed: 03/01/2024]
Abstract
Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule-metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.
Collapse
Affiliation(s)
- Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, Post Office
Box 9502, 2300 RA Leiden, Netherlands
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, BE-2610, Wilrijk, Antwerp, Belgium
| | - Bret Jackson
- Department
of Chemistry, University of Massachusetts
Amherst, Amherst, Massachusetts 01003, United States
| | - Annemie Bogaerts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, BE-2610, Wilrijk, Antwerp, Belgium
| |
Collapse
|
5
|
Zhang L, Zuo J, Suleimanov YV, Guo H. Ring Polymer Molecular Dynamics Approach to Quantum Dissociative Chemisorption Rates. J Phys Chem Lett 2023; 14:7118-7125. [PMID: 37531595 DOI: 10.1021/acs.jpclett.3c01848] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2023]
Abstract
A ring polymer molecular dynamics (RPMD) method is proposed for the calculation of the dissociative chemisorption rate coefficient on surfaces. The RPMD rate theory is capable of handling quantum effects such as the zero-point energy and tunneling in dissociative chemisorption, while it relies on classical trajectories for the simulation. Applications to H2 dissociative chemisorption are demonstrated. For the highly activated process on Ag(111), strong deviations from Arrhenius behavior are found at low temperatures and attributed to tunneling. On Pt(111), where the dissociation has a barrierless pathway, the RPMD rate coefficient is found to agree with the experimentally derived thermal sticking coefficient within a factor of 2 over a large temperature range. Significant quantum effects are also found.
Collapse
Affiliation(s)
- Liang Zhang
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Junxiang Zuo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yury V Suleimanov
- American Association for the Advancement of Science, 1200 New York Ave NW, Washington, D.C. 20005, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
6
|
Li C, Li Y, Jiang B. First-principles surface reaction rates by ring polymer molecular dynamics and neural network potential: role of anharmonicity and lattice motion. Chem Sci 2023; 14:5087-5098. [PMID: 37206404 PMCID: PMC10189860 DOI: 10.1039/d2sc06559b] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 04/05/2023] [Indexed: 08/04/2023] Open
Abstract
Elementary gas-surface processes are essential steps in heterogeneous catalysis. A predictive understanding of catalytic mechanisms remains challenging due largely to difficulties in accurately characterizing the kinetics of such steps. Experimentally, thermal rates for elementary surface reactions can now be measured using a novel velocity imaging technique, providing a stringent testing ground for ab initio rate theories. Here, we propose to combine ring polymer molecular dynamics (RPMD) rate theory with state-of-the-art first-principles-determined neural network potential to calculate surface reaction rates. Taking NO desorption from Pd(111) as an example, we show that the harmonic approximation and the neglect of lattice motion in the commonly-used transition state theory overestimates and underestimates the entropy change during the desorption process, respectively, leading to opposite errors in rate coefficient predictions and artificial error cancellations. Including anharmonicity and lattice motion, our results reveal a generally neglected surface entropy change due to significant local structural change during desorption and obtain the right answer for the right reasons. Although quantum effects are found to be less important in this system, the proposed approach establishes a more reliable theoretical benchmark for accurately predicting the kinetics of elementary gas-surface processes.
Collapse
Affiliation(s)
- Chen Li
- Key Laboratory of Precision and Intelligent Chemistry, 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 Anhui 230026 China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures, Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University Shanghai 200444 China
| | - Bin Jiang
- Key Laboratory of Precision and Intelligent Chemistry, 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 Anhui 230026 China
| |
Collapse
|
7
|
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
|
8
|
Young TA, Johnston-Wood T, Zhang H, Duarte F. Reaction dynamics of Diels-Alder reactions from machine learned potentials. Phys Chem Chem Phys 2022; 24:20820-20827. [PMID: 36004770 DOI: 10.1039/d2cp02978b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent advances in the development of reactive machine-learned potentials (MLPs) promise to transform reaction modelling. However, such methods have remained computationally expensive and limited to experts. Here, we employ different MLP methods (ACE, NequIP, GAP), combined with automated fitting and active learning, to study the reaction dynamics of representative Diels-Alder reactions. We demonstrate that the ACE and NequIP MLPs can consistently achieve chemical accuracy (±1 kcal mol-1) to the ground-truth surface with only a few hundred reference calculations. These strategies are shown to enable routine ab initio-quality classical and quantum dynamics, and obtain dynamical quantities such as product ratios and free energies from non-static methods. For ambimodal reactions, product distributions were found to be strongly dependent on the QM method and less so on the type of dynamics propagated.
Collapse
Affiliation(s)
- Tom A Young
- Chemistry Research Laboratory, 12 Mansfield Road, Oxford, OX1 3TA, UK.
| | | | - Hanwen Zhang
- Chemistry Research Laboratory, 12 Mansfield Road, Oxford, OX1 3TA, UK.
| | - Fernanda Duarte
- Chemistry Research Laboratory, 12 Mansfield Road, Oxford, OX1 3TA, UK.
| |
Collapse
|
9
|
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
|
10
|
Marjollet A, Inhester L, Welsch R. Initial state-selected scattering for the reactions H + CH4/CHD3 and F + CHD3 employing ring polymer molecular dynamics. J Chem Phys 2022; 156:044101. [DOI: 10.1063/5.0076216] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- A. Marjollet
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
- Department of Physics, Universität Hamburg, Notkestr. 9-11, 22607 Hamburg, Germany
| | - L. Inhester
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - R. Welsch
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
11
|
Theoretical Description of Water from Single-Molecule to Condensed Phase: a Review of Recent Progress on Potential Energy Surfaces and Molecular Dynamics. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2201005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
12
|
Bhowmick S, Hernández MI, Campos-Martínez J, Suleimanov YV. Isotopic separation of helium through graphyne membranes: a ring polymer molecular dynamics study. Phys Chem Chem Phys 2021; 23:18547-18557. [PMID: 34612392 DOI: 10.1039/d1cp02121d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Microscopic-level understanding of the separation mechanism for two-dimensional (2D) membranes is an active area of research due to potential implications of this class of membranes for various technological processes. Helium (He) purification from the natural resources is of particular interest due to the shortfall in its production. In this work, we applied the ring polymer molecular dynamics (RPMD) method to graphdiyne (Gr2) and graphtriyne (Gr3) 2D membranes having variable pore sizes for the separation of He isotopes, and compare for the first time with rigorous quantum calculations. We found that the transmission rate through Gr3 is many orders of magnitude greater than Gr2. The selectivity of either isotope at low temperatures is a consequence of a delicate balance between the zero-point energy effect and tunneling of 4He and 3He. In particular, a remarkable tunneling effect is reported on the Gr2 membrane at 10 K, leading to a much larger permeation of the lighter species as compared to the heavier isotope. RPMD provides an efficient approach for studying the separation of He isotopes, taking into account quantum effects of light nuclei motions at low temperatures, which classical methods fail to capture.
Collapse
Affiliation(s)
- Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | | | | | | |
Collapse
|
13
|
Zhou X, Zhang Y, Yin R, Hu C, Jiang B. Neural Network Representations for Studying
Gas‐Surface
Reaction Dynamics: Beyond the
Born‐Oppenheimer
Static Surface Approximation
†. CHINESE J CHEM 2021. [DOI: 10.1002/cjoc.202100303] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Rongrong Yin
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Ce Hu
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| |
Collapse
|
14
|
Li C, Liu Q, Zhang L, Li Y, Jiang B. Ring polymer molecular dynamics in gas-surface reactions: tests on initial sampling and potential energy landscape. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1941367] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Chen Li
- 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, People’s Republic of China
| | - Qinghua Liu
- 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, People’s Republic of China
| | - Liang Zhang
- 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, People’s Republic of China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai, People’s Republic of China
| | - 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, People’s Republic of China
| |
Collapse
|
15
|
Bünermann O, Kandratsenka A, Wodtke AM. Inelastic Scattering of H Atoms from Surfaces. J Phys Chem A 2021; 125:3059-3076. [PMID: 33779163 PMCID: PMC8154602 DOI: 10.1021/acs.jpca.1c00361] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Revised: 03/09/2021] [Indexed: 11/29/2022]
Abstract
We have developed an instrument that uses photolysis of hydrogen halides to produce nearly monoenergetic hydrogen atom beams and Rydberg atom tagging to obtain accurate angle-resolved time-of-flight distributions of atoms scattered from surfaces. The surfaces are prepared under strict ultrahigh vacuum conditions. Data from these experiments can provide excellent benchmarks for theory, from which it is possible to obtain an atomic scale understanding of the underlying dynamical processes governing H atom adsorption. In this way, the mechanism of adsorption on metals is revealed, showing a penetration-resurfacing mechanism that relies on electronic excitation of the metal by the H atom to succeed. Contrasting this, when H atoms collide at graphene surfaces, the dynamics of bond formation involving at least four carbon atoms govern adsorption. Future perspectives of H atom scattering from surfaces are also outlined.
Collapse
Affiliation(s)
- Oliver Bünermann
- Institute
for Physical Chemistry, Georg-August-University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
- Department
of Dynamics at Surfaces, Max-Planck Institute
for Biophysical Chemistry, Am Faßberg 11, 37077 Göttingen, Germany
- International
Center for Advanced Studies of Energy Conversion, Georg-August University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Alexander Kandratsenka
- Department
of Dynamics at Surfaces, Max-Planck Institute
for Biophysical Chemistry, Am Faßberg 11, 37077 Göttingen, Germany
| | - Alec M. Wodtke
- Institute
for Physical Chemistry, Georg-August-University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
- Department
of Dynamics at Surfaces, Max-Planck Institute
for Biophysical Chemistry, Am Faßberg 11, 37077 Göttingen, Germany
- International
Center for Advanced Studies of Energy Conversion, Georg-August University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
16
|
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
|
17
|
Jiang H, Tao X, Kammler M, Ding F, Wodtke AM, Kandratsenka A, Miller TF, Bünermann O. Small Nuclear Quantum Effects in Scattering of H and D from Graphene. J Phys Chem Lett 2021; 12:1991-1996. [PMID: 33596383 DOI: 10.1021/acs.jpclett.0c02933] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We study nuclear quantum effects in H/D sticking to graphene, comparing scattering experiments at near-zero coverage with classical, quantized, and transition-state calculations. The experiment shows H/D sticking probabilities that are indistinguishable from one another and markedly smaller than those expected from a consideration of zero-point energy shifts of the chemisorption transition state. Inclusion of dynamical effects and vibrational anharmonicity via ring-polymer molecular dynamics (RPMD) yields results that are in good agreement with the experimental results. RPMD also reveals that nuclear quantum effects, while modest, arise primarily from carbon and not from H/D motion, confirming the importance of a C atom rehybridization mechanism associated with H/D sticking on graphene.
Collapse
Affiliation(s)
- Hongyan Jiang
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
| | - Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Marvin Kammler
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
| | - Feizhi Ding
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Alec M Wodtke
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
- Institut für Physikalische Chemie, Georg-Georg-August-Universität Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Georg-August University of Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
| | - Alexander Kandratsenka
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Oliver Bünermann
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
- Institut für Physikalische Chemie, Georg-Georg-August-Universität Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Georg-August University of Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
| |
Collapse
|
18
|
Chen Q, Hu X, Guo H, Xie D. Theoretical H + O 3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties. Phys Chem Chem Phys 2021; 23:3300-3310. [PMID: 33506830 DOI: 10.1039/d0cp05771a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. The theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O3 → OD + O2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.
Collapse
Affiliation(s)
- Qixin Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Xixi Hu
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| |
Collapse
|
19
|
Zhang L, Jiang B. State-to-state quantum dynamics of H2O/HOD scattering from Cu(111): Mode- and bond-selective vibrational energy transfer. J Chem Phys 2020; 153:214702. [DOI: 10.1063/5.0030490] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Liang Zhang
- 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, Anhui 230026, China
| | - 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, Anhui 230026, China
| |
Collapse
|
20
|
Li J, Zhao B, Xie D, Guo H. Advances and New Challenges to Bimolecular Reaction Dynamics Theory. J Phys Chem Lett 2020; 11:8844-8860. [PMID: 32970441 DOI: 10.1021/acs.jpclett.0c02501] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.
Collapse
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
21
|
Marjollet A, Welsch R. Nuclear quantum effects in state-selective scattering from ring polymer molecular dynamics. J Chem Phys 2020; 152:194113. [DOI: 10.1063/5.0004179] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrien Marjollet
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- Department of Physics, Universität Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
22
|
Jiang B, Guo H. Erratum: Dynamics of Water Dissociative Chemisorption on Ni(111): Effects of Impact Sites and Incident Angles [Phys. Rev. Lett. 114, 166101 (2015)]. PHYSICAL REVIEW LETTERS 2020; 124:149901. [PMID: 32338964 DOI: 10.1103/physrevlett.124.149901] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Accepted: 03/25/2020] [Indexed: 06/11/2023]
Abstract
This corrects the article DOI: 10.1103/PhysRevLett.114.166101.
Collapse
|
23
|
Tao X, Shushkov P, Miller TF. Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches. J Chem Phys 2020; 152:124117. [DOI: 10.1063/1.5144307] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Philip Shushkov
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
24
|
Zhu L, Zhang Y, Zhang L, Zhou X, Jiang B. Unified and transferable description of dynamics of H2 dissociative adsorption on multiple copper surfaces via machine learning. Phys Chem Chem Phys 2020; 22:13958-13964. [DOI: 10.1039/d0cp02291h] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Schematic of the developed neural network potential energy surface enabling a unified and transferable description of dynamics of H2 dissociative adsorption on multiple copper surfaces.
Collapse
Affiliation(s)
- Lingjun Zhu
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Liang Zhang
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| |
Collapse
|