1
|
Oudot B, Doblhoff-Dier K. Reaction barriers at metal surfaces computed using the random phase approximation: Can we beat DFT in the generalized gradient approximation? J Chem Phys 2024; 161:054708. [PMID: 39092949 DOI: 10.1063/5.0220465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2024] [Accepted: 07/14/2024] [Indexed: 08/04/2024] Open
Abstract
Reaction barriers for molecules dissociating on metal surfaces (as relevant to heterogeneous catalysis) are often difficult to predict accurately with density functional theory (DFT). Although the results obtained for several dissociative chemisorption reactions via DFT in the generalized gradient approximation (GGA), in meta-GGA, and for GGA exchange + van der Waals correlation scatter around the true reaction barrier, there is an entire class of dissociative chemisorption reactions for which GGA-type functionals collectively underestimate the reaction barrier. Little is known why GGA-DFT collectively fails in some cases and not in others, and we do not know whether other methods suffer from the same inconsistency. Here, we present barrier heights for dissociative chemisorption reactions obtained from the random phase approximation in the adiabatic-connection fluctuation-dissipation theorem (ACFDT-RPA) and from hybrid functionals with different amounts of exact exchange. By comparing the results obtained for the dissociative chemisorption reaction of H2 on Al(110) (where GGA-DFT collectively underestimates the barrier) and H2 on Cu(111) (where GGA-DFT scatters around the true barrier), we can gauge whether the inconsistent description of the systems persists for hybrid functionals and ACFDT-RPA. We find hybrid functionals to improve the relative description of the two systems, but to fall short of chemical accuracy. ACFDT-RPA improves the results further and leads to chemically accurate barriers for both systems. Together with an analysis of the density of states and the results from selected GGA, meta-GGA, and GGA exchange + van der Waals correlation functionals, these results allow us to discuss possible origins for the inconsistent behavior of GGA-based functionals for molecule-metal reaction barriers.
Collapse
Affiliation(s)
- B Oudot
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, Leiden 2300 RA, The Netherlands
| | - K Doblhoff-Dier
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, Leiden 2300 RA, The Netherlands
| |
Collapse
|
2
|
Chen D, Xu Y, Lu J, Tian Y, Li T, Jia P, Wang X, Zhang L, Hou Y, Wang L, Zhang Q, Ye Z, Lu J. Intercalation-Induced Localized Conversion Reaction in h-CuSe for Ultrafast-Rechargeable and Long-Cycling Sodium Metal Battery. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2404640. [PMID: 38775475 DOI: 10.1002/adma.202404640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2024] [Revised: 05/07/2024] [Indexed: 06/13/2024]
Abstract
Cathode materials of sodium-based batteries with high specific capacity and fast charge-discharge mode, as well as ultralong reversible cycles at wide applied temperatures, are essential for future development of advanced energy storage system. Developing transition metal selenides with intercalation features provides a new strategy for realizing the above cathode materials. Herein, this work reports a storage mechanism of sodium ion in hexagonal CuSe (h-CuSe) based on the density functional theory (DFT) guidance. This work reveals that the two-dimensional ion intercalation triggers localized redox reaction in the h-CuSe bulk phase, termed intercalation-induced localized conversion (ILC) mechanism, to stabilize the sodium storage structure by forming localized Cu7Se4 transition phase and adjusting the near-edge coordination state of the Cu sites to achieve high reversible capacity and ultra-long cycling life, while allowing rapid charge-discharge cycling over a wide temperature range.
Collapse
Affiliation(s)
- Dongliang Chen
- State Key Laboratory of Silicon and Advanced Semiconductor Materials, School of Materials Science and Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Yunkai Xu
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Jianguo Lu
- State Key Laboratory of Silicon and Advanced Semiconductor Materials, School of Materials Science and Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Yang Tian
- Zhijiang Lab, Hangzhou, 311121, China
| | - Tongtong Li
- School of Materials Science and Engineering, Zhejiang Sci-Tech University, Hangzhou, 310018, China
| | - Peng Jia
- State Key Laboratory of Metastable Materials Science and Technology, School of Materials Science and Engineering, Yanshan University, Qinhuangdao, 066004, China
| | - Xu Wang
- State Key Laboratory of Silicon and Advanced Semiconductor Materials, School of Materials Science and Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Liqiang Zhang
- State Key Laboratory of Metastable Materials Science and Technology, School of Materials Science and Engineering, Yanshan University, Qinhuangdao, 066004, China
| | - Yang Hou
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Liguang Wang
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Qinghua Zhang
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Zhizhen Ye
- State Key Laboratory of Silicon and Advanced Semiconductor Materials, School of Materials Science and Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Jun Lu
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou, 310027, China
| |
Collapse
|
3
|
Chadwick H. Characterisation of magnetic atomic and molecular beamlines for the extraction of empirical scattering-matrices. Phys Chem Chem Phys 2024; 26:19630-19645. [PMID: 38888009 PMCID: PMC11267877 DOI: 10.1039/d4cp01785d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2024] [Accepted: 06/11/2024] [Indexed: 06/20/2024]
Abstract
A recently developed magnetic molecular interferometry technique allows the experimental determination of how the amplitudes and phases of the molecular wave-function change during the collision of a gas phase molecule with a surface. This information, quantified by a scattering-matrix, provides a very stringent benchmark for developing accurate theoretical models as they can also be determined from scattering calculations and are particularly sensitive to the underlying interaction potential. However, the value of this comparison is necessarily limited by the accuracy with which an empirical scattering-matrix can be extracted from the experimental data. This paper presents the methods used to analyse the measurements and uses simulations to determine how various uncertainties in modelling the different magnetic elements which make up the beamline of the apparatus affect the accuracy with which the scattering-matrix can be extracted. It is shown that when signals have a noise level which corresponds to on the order of 1% of the oscillation amplitude, the uncertainties in the modelling do not significantly affect the ability to extract the scattering-matrix elements, with the error in the extracted values increasing to a few percent as the noise in the signals is increased to 10% of the oscillation amplitude. This therefore gives an estimate of the accuracy of the parameters that can be obtained from future measurements.
Collapse
Affiliation(s)
- Helen Chadwick
- Department of Chemistry, Faculty of Science and Engineering, Swansea University, Swansea, SA2 8PP, UK.
| |
Collapse
|
4
|
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
|
5
|
Cai Y, Michiels R, De Luca F, Neyts E, Tu X, Bogaerts A, Gerrits N. Improving Molecule-Metal Surface Reaction Networks Using the Meta-Generalized Gradient Approximation: CO 2 Hydrogenation. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:8611-8620. [PMID: 38835935 PMCID: PMC11145648 DOI: 10.1021/acs.jpcc.4c01110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/20/2024] [Revised: 05/07/2024] [Accepted: 05/09/2024] [Indexed: 06/06/2024]
Abstract
Density functional theory is widely used to gain insights into molecule-metal surface reaction networks, which is important for a better understanding of catalysis. However, it is well-known that generalized gradient approximation (GGA) density functionals (DFs), most often used for the study of reaction networks, struggle to correctly describe both gas-phase molecules and metal surfaces. Also, GGA DFs typically underestimate reaction barriers due to an underestimation of the self-interaction energy. Screened hybrid GGA DFs have been shown to reduce this problem but are currently intractable for wide usage. In this work, we use a more affordable meta-GGA (mGGA) DF in combination with a nonlocal correlation DF for the first time to study and gain new insights into a catalytically important surface reaction network, namely, CO2 hydrogenation on Cu. We show that the mGGA DF used, namely, rMS-RPBEl-rVV10, outperforms typical GGA DFs by providing similar or better predictions for metals and molecules, as well as molecule-metal surface adsorption and activation energies. Hence, it is a better choice for constructing molecule-metal surface reaction networks.
Collapse
Affiliation(s)
- Yuxiang Cai
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
- Department
of Electrical Engineering and Electronics, University of Liverpool, Liverpool L69 3GJ, U.K.
| | - Roel Michiels
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
| | - Federica De Luca
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
- Department
of ChiBioFarAM (Industrial Chemistry), ERIC aisbl and INSTM/CASPE, University of Messina, V.le F. Stagno d’Alcontres 31, Messina 98166, Italy
| | - Erik Neyts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
| | - Xin Tu
- Department
of Electrical Engineering and Electronics, University of Liverpool, Liverpool L69 3GJ, U.K.
| | - Annemie Bogaerts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
| | - Nick Gerrits
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, Antwerp, Wilrijk BE-2610, Belgium
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, Leiden 2300 RA, The Netherlands
| |
Collapse
|
6
|
van Bree RAB, Gerrits N, Kroes GJ. Dissociative chemisorption of O 2 on Al(111): dynamics on a potential energy surface computed with a non-self-consistent screened hybrid density functional approach. Faraday Discuss 2024. [PMID: 38787655 DOI: 10.1039/d3fd00165b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2024]
Abstract
Density functional theory (DFT) at the generalized gradient approximation (GGA) level is often considered the best compromise between feasibility and accuracy for reactions of molecules on metal surfaces. Recent work, however, strongly suggests that density functionals (DFs) based on GGA exchange are not able to describe molecule-metal surface reactions for which the work function of the metal surface minus the electron affinity of the molecule is less than 7 eV. Systems for which this is true exhibit an increased charge transfer from the metal to the molecule at the transition state, increasing the delocalisation of the electron density. This enlarged delocalisation can cause GGA-DFT to underestimate energy values relative to the gas-phase and thus underestimate the barrier height, similar to what has been observed for several gas-phase reactions. An example of such a molecule-metal surface system is O2 + Al(111). Following a similar strategy as for gas-phase reactions, previous work showed results of increased accuracy when using a screened hybrid DF for O2 + Al(111). However, even screened hybrid DFs are computationally expensive to use for metal surfaces. To resolve this, we test a non-self-consistent field (NSCF) screened hybrid DF approach. This approach computes screened hybrid DFT energies based on self-consistent-field (SCF) GGA electronic densities. Here, we explore the accuracy of the NSCF screened hybrid DF approach by implementing the NSCF HSE03-1/3x@RPBE DF for O2 + Al(111). We compute and analyse molecular beam sticking probabilities as well as a set of sticking probabilities for rotationally aligned O2. Our results show that the NSCF approach results in reaction probability curves that reproduce SCF results with near-chemical accuracy, suggesting that the NSCF approach can be used advantageously for exploratory purposes. An analysis of the potential energy surface and the barriers gives insight into the cause of the disagreement between the SCF and NSCF reaction probabilities and into the changes needed in theoretical modelling to further improve the description of the O2 + Al(111) system. Finally, the hole model yields fair agreement with dynamics results for the reaction probability curve, but results in an increased slope of the reaction probability curve compared to the molecular dynamics, with a shift to lower or higher energies depending on whether the vibrational energy of the molecule is included in the initial energy of the molecule or not.
Collapse
Affiliation(s)
- Robert A B van Bree
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| | - Nick Gerrits
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| | - Geert-Jan Kroes
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| |
Collapse
|
7
|
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
|
8
|
Bernard ME, Harrison I. Microcanonical treatment of HCl dissociative chemisorption on Au(111): Reactive dampening through inefficient translational energy coupling and an active surface. J Chem Phys 2024; 160:084702. [PMID: 38391017 DOI: 10.1063/5.0193675] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 02/02/2024] [Indexed: 02/24/2024] Open
Abstract
Microcanonical unimolecular rate theory is applied to Shirhatti and Wodtke's recent supersonic molecular beam experiments examining the activated dissociative chemisorption of HCl on Au(111). A precursor mediated microcanonical trapping (PMMT) model (where the surface vibrates and HCl rotations, vibration, and translation directed along the surface normal are treated as active degrees of freedom) gave dissociative sticking coefficient predictions that are several orders of magnitude higher than experimental values but in good accord with prior quantum and molecular dynamics simulations. Density functional theory (DFT) electronic structure calculations using the Perdew-Burke-Ernzerhof (PBE) functional served to fix the vibrational frequencies of the reactive transition state and the threshold energy for dissociation, E0 = 72.9 kJ/mol. To explore the possibilities of varying threshold energy, coupling to phonons, and dynamics, a three-parameter [E0, s, ɛn] dynamically biased (d-) PMMT model was fit to the experiments. A dynamical bias was introduced using an efficiency, ɛn, of normal translational energy to contribute to the active exchangeable energy capable of promoting reactivity. To achieve the low sticking probabilities observed in experiment, severe normal translational energy dampening (ɛn → 0.26) was imposed, leading to a large vibrational efficacy of ηv = εv/εn = 3.85. The optimal threshold energy for dissociation was E0 = 30.88 kJ/mol, some 40 kJ/mol below the PBE-DFT prediction, and the optimal number of Au surface oscillators was s = 1. The d-PMMT modeling indicates that HCl/Au(111) reactivity can be consistent with electronically adiabatic passage across a relatively low and late transition state that dynamically disfavors normal translational energy.
Collapse
Affiliation(s)
- Mark E Bernard
- Department of Chemistry, University of Virgina, Charlottesville, Virginia 22904, USA
| | - Ian Harrison
- Department of Chemistry, University of Virgina, Charlottesville, Virginia 22904, USA
| |
Collapse
|
9
|
Powell A, Gerrits N, Tchakoua T, Somers MF, Busnengo HF, Meyer J, Kroes GJ, Doblhoff-Dier K. Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals. J Phys Chem Lett 2024; 15:307-315. [PMID: 38169287 PMCID: PMC10788952 DOI: 10.1021/acs.jpclett.3c02972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 12/21/2023] [Accepted: 12/28/2023] [Indexed: 01/05/2024]
Abstract
Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.
Collapse
Affiliation(s)
- Andrew
D. Powell
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Theophile Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Mark F. Somers
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Heriberto F. Busnengo
- Instituto
de Física Rosario (IFIR), CONICET-UNR, 2000 Rosario, Argentina
- Facultad
de Ciencias Exatas, Ingeniería y
Agrimensura, UNR, 2000 Rosario, Argentina
| | - Jörg Meyer
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| |
Collapse
|
10
|
Xie Y, Han J, Wen L, Li Z, Xiao Y, Wang T, Yang X, Yang T. Enhanced Total Vibrational Excitation Yield in a Slow Narrow-Pulsed Hydrogen Molecular Beam. J Phys Chem Lett 2023; 14:11603-11609. [PMID: 38100090 DOI: 10.1021/acs.jpclett.3c03015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2023]
Abstract
High-efficiency excitation of a molecular beam is critical for investigating state-selected chemistry. However, achieving vibrational excitation of the entire beam for Raman-active molecules such as H2 proves extremely challenging, primarily because laser pulses are much shorter than the molecular beam. In this study, we achieve a total excitation efficiency of over 20% by employing stimulated Raman pumping (SRP) in a slow, narrow-pulsed molecular beam. Through optimizing the intensity and spot shape of the SRP lasers, we attain saturated excitation within the laser crossing region. Furthermore, by reducing the beam velocity and narrowing the beam pulse using a cold valve and a fast chopper, we significantly enhance the total excitation yield. COMSOL simulation and a newly developed model reveal that a critical velocity allows the chopper to block unexcited molecules and reserve most of the excited ones from the beam, resulting in the highest overall excitation yield. This innovative setup opens new possibilities for state-selected experiments in surface science and ion-molecule reaction dynamics, particularly involving weak transitions and pulsed lasers.
Collapse
Affiliation(s)
- Yurun Xie
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
- Institute of Advanced Science Facilities, Shenzhen, Guangdong 518107, China
| | - Jie Han
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Liping Wen
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Zhichao Li
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Yue Xiao
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Tao Wang
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Xueming Yang
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Dalian, Liaoning 116023, China
- Hefei National Laboratory, Hefei, Anhui 230088, China
| | - Tiangang Yang
- Shenzhen Key Laboratory of Energy Chemistry, Department of Chemistry, and Center for Advanced Light Source, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| |
Collapse
|
11
|
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
|
12
|
Tchakoua T, Jansen T, van Nies Y, van den Elshout RFA, van Boxmeer BAB, Poort SP, Ackermans MG, Beltrão GS, Hildebrand SA, Beekman SEJ, van der Drift T, Kaart S, Šantić A, Spuijbroek EE, Gerrits N, Somers MF, Kroes GJ. Constructing Mixed Density Functionals for Describing Dissociative Chemisorption on Metal Surfaces: Basic Principles. J Phys Chem A 2023; 127:10481-10498. [PMID: 38051300 PMCID: PMC10726370 DOI: 10.1021/acs.jpca.3c01932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
The production of a majority of chemicals involves heterogeneous catalysis at some stage, and the rates of many heterogeneously catalyzed processes are governed by transition states for dissociative chemisorption on metals. Accurate values of barrier heights for dissociative chemisorption on metals are therefore important to benchmarking electronic structure theory in general and density functionals in particular. Such accurate barriers can be obtained using the semiempirical specific reaction parameter (SRP) approach to density functional theory. However, this approach has thus far been rather ad hoc in its choice of the generic expression of the SRP functional to be used, and there is a need for better heuristic approaches to determining the mixing parameters contained in such expressions. Here we address these two issues. We investigate the ability of several mixed, parametrized density functional expressions combining exchange at the generalized gradient approximation (GGA) level with either GGA or nonlocal correlation to reproduce barrier heights for dissociative chemisorption on metal surfaces. For this, seven expressions of such mixed density functionals are tested on a database consisting of results for 16 systems taken from a recently published slightly larger database called SBH17. Three expressions are derived that exhibit high tunability and use correlation functionals that are either of the PBE GGA form or of one of two limiting nonlocal forms also describing the attractive van der Waals interaction in an approximate way. We also find that, for mixed density functionals incorporating GGA correlation, the optimum fraction of repulsive RPBE GGA exchange obtained with a specific GGA density functional is correlated with the charge-transfer parameter, which is equal to the difference in the work function of the metal surface and the electron affinity of the molecule. However, the correlation is generally not large and not large enough to obtain accurate guesses of the mixing parameter for the systems considered, suggesting that it does not give rise to a very effective search strategy.
Collapse
Affiliation(s)
- Théophile Tchakoua
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Tim Jansen
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Youri van Nies
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | | | - Bart A B van Boxmeer
- Faculty of Applied Engineering, University of Antwerp, Groenenborgerlaan 171, 2020 Antwerpen, Belgium
| | - Saskia P Poort
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Michelle G Ackermans
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Gabriel Spiller Beltrão
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Stefan A Hildebrand
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Steijn E J Beekman
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Thijs van der Drift
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Sam Kaart
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Anthonie Šantić
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Esmee E Spuijbroek
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Nick Gerrits
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
13
|
Wei Z, Martirez JMP, Carter EA. Introducing the embedded random phase approximation: H2 dissociative adsorption on Cu(111) as an exemplar. J Chem Phys 2023; 159:194108. [PMID: 37971031 DOI: 10.1063/5.0181229] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2023] [Accepted: 10/16/2023] [Indexed: 11/19/2023] Open
Abstract
The random phase approximation (RPA) as a means of treating electron correlation recently has been shown to outperform standard density functional theory (DFT) approximations in a variety of cases. However, the computational cost of the RPA is substantially more than DFT, especially when aiming to study extended surfaces. Properly accounting for sufficient surface ensemble size, Brillouin zone sampling, and vacuum separation of periodic images in standard periodic-planewave-based DFT code raises the cost to achieve converged results. Here, we show that sub-system embedding schemes enable use of the RPA for modeling heterogeneous reactions at reduced computational cost. We explore two different embedded RPA (emb-RPA) approaches, periodic emb-RPA and cluster emb-RPA. We use the (experimentally and theoretically) well-studied H2 dissociative adsorption on Cu(111) as our exemplar, and first perform full periodic RPA calculations as a benchmark. The full RPA results match well the semi-empirical barrier fit to experimental observables and others derived from high-level computations, e.g., from recent embedded n-electron valence second order perturbation theory [Zhao et al., J. Chem. Theory Comput. 16(11), 7078-7088 (2020)] and quantum Monte Carlo [Doblhoff-Dier et al., J. Chem. Theory Comput. 13(7), 3208-3219 (2017)] simulations. Among the two emb-RPA approaches tested, the cluster emb-RPA accurately reproduces the energy profile (maximum error of 50 meV along the reaction pathway) while reducing the computational cost by approximately two orders of magnitude. We therefore expect that the embedded cluster approach will enable wider RPA implementation in heterogeneous catalysis.
Collapse
Affiliation(s)
- Ziyang Wei
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, USA
| | - John Mark P Martirez
- Applied Materials and Sustainability Sciences, Princeton Plasma Physics Laboratory, Princeton, New Jersey 08540-6655, USA
| | - Emily A Carter
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, USA
- Applied Materials and Sustainability Sciences, Princeton Plasma Physics Laboratory, Princeton, New Jersey 08540-6655, USA
- Andlinger Center for Energy and the Environment and Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
14
|
Liu T, Peng T, Fu B, Zhang DH. Charge-Transfer-Controlled Quantum Dynamics of HCl Dissociation on the Ag/Au(111) Bimetallic Alloy Surface. J Phys Chem Lett 2023; 14:9713-9719. [PMID: 37877754 DOI: 10.1021/acs.jpclett.3c02556] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2023]
Abstract
Understanding polar molecule dynamics on bimetallic surfaces, especially electropositivity and electronegativity, remains a challenge. Here, we report the reactivity of HCl on a strained Ag monolayer on Au(111) using six-dimensional quantum dynamics with a new machine-learning-based potential energy surface. Surprisingly, HCl reactivity is significantly suppressed by the Ag-Au interaction despite a lower HCl+Ag/Au(111) barrier than pure Ag(111). This arises from charge transfer between Ag and Au, where electronegative Au makes the top Ag layer on Ag/Au(111) electropositive, unlike that on pure Ag(111). Electropositive Ag in HCl+Ag/Au(111) attracts Cl, yielding an unfavorable H-Cl configuration and reduced reactivity. These findings deepen our understanding of polar molecule interactions on bimetallic surfaces, highlighting the role of charge transfer in dissociative chemisorption and the implications for catalyst design in heterogeneous catalysis.
Collapse
Affiliation(s)
- Tianhui Liu
- School of Sciences, Great Bay University, Songshan Lake International Innovation Entrepreneurship Community A5, Dongguan 523000, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Tianze Peng
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Hefei National Laboratory, Hefei 230088, China
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
15
|
Li Y, Wang D. Ab Initio Molecular Dynamics Study of H 2 Dissociation Mechanisms on Cu 13 and Defective Graphene-supported Cu 13 Clusters: Active Sites, Energy Barriers and Adsorption States. Chemphyschem 2023; 24:e202300369. [PMID: 37439149 DOI: 10.1002/cphc.202300369] [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: 05/26/2023] [Revised: 07/12/2023] [Accepted: 07/13/2023] [Indexed: 07/14/2023]
Abstract
Ab initio molecular dynamics calculations were performed to study H2 dissociation mechanisms on Cu13 and defective graphene-supported Cu13 clusters. The study reveals that seven types of corresponding dissociation processes are found on the two clusters. The average dissociation energy barriers are 0.51 eV on the Cu13 cluster and 0.12 eV on the defective graphene-supported Cu13 cluster, which are lowered by ~19 % and ~81 % compared with the pristine Cu(111) surface, respectively. Furthermore, compared with the pure Cu13 cluster, the average dissociation energy barrier on the defective graphene-supported Cu13 cluster is substantially reduced by about 76 %. The preferred dissociation mechanisms on the two clusters are H2 located at a top-bridge site with the barrier heights of 0.30 eV on the Cu13 cluster and -0.31 eV on the defective graphene-supported Cu13 cluster. Analysis of the H-Cu bond interactions in the transition states shows that the antibonding-orbital center shifts upward on the defective graphene-supported Cu13 cluster compared with the one on the Cu13 cluster, which explains the reduction of the dissociation energy barrier. The average adsorption energy of dissociated H atoms is also greatly enhanced on the defective graphene-supported Cu13 cluster, about twice that on the Cu13 cluster.
Collapse
Affiliation(s)
- Yueru Li
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, Shandong, China
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, Shandong, China
| |
Collapse
|
16
|
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
|
17
|
Chen Z, Liu Z, Xu X. Accurate descriptions of molecule-surface interactions in electrocatalytic CO 2 reduction on the copper surfaces. Nat Commun 2023; 14:936. [PMID: 36807556 PMCID: PMC9941474 DOI: 10.1038/s41467-023-36695-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 02/14/2023] [Indexed: 02/22/2023] Open
Abstract
Copper-based catalysts play a pivotal role in many industrial processes and hold a great promise for electrocatalytic CO2 reduction reaction into valuable chemicals and fuels. Towards the rational design of catalysts, the growing demand on theoretical study is seriously at odds with the low accuracy of the most widely used functionals of generalized gradient approximation. Here, we present results using a hybrid scheme that combines the doubly hybrid XYG3 functional and the periodic generalized gradient approximation, whose accuracy is validated against an experimental set on copper surfaces. A near chemical accuracy is established for this set, which, in turn, leads to a substantial improvement for the calculated equilibrium and onset potentials as against the experimental values for CO2 reduction to CO on Cu(111) and Cu(100) electrodes. We anticipate that the easy use of the hybrid scheme will boost the predictive power for accurate descriptions of molecule-surface interactions in heterogeneous catalysis.
Collapse
Affiliation(s)
- Zheng Chen
- grid.8547.e0000 0001 0125 2443Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, 200433 Shanghai, People’s Republic of China
| | - Zhangyun Liu
- grid.8547.e0000 0001 0125 2443Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, 200433 Shanghai, People’s Republic of China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, 200433, Shanghai, People's Republic of China. .,Hefei National Laboratory, 230088, Hefei, P. R. China.
| |
Collapse
|
18
|
Zhu L, Hu C, Chen J, Jiang B. Investigating the Eley-Rideal recombination of hydrogen atoms on Cu (111) via a high-dimensional neural network potential energy surface. Phys Chem Chem Phys 2023; 25:5479-5488. [PMID: 36734463 DOI: 10.1039/d2cp05479e] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
As a prototypical system for studying the Eley-Rideal (ER) mechanism at the gas-surface interface, the reaction between incident H/D atoms and pre-covered D/H atoms on Cu (111) has attracted much experimental and theoretical interest. Detailed final state-resolved experimental data have been available for about thirty-years, leading to the discovery of many interesting dynamical features. However, previous theoretical models have suffered from reduced-dimensional approximations and/or omitting energy transfer to surface phonons and electrons, or the high cost of on-the-fly ab initio molecular dynamics, preventing quantitative comparisons with experimental data. Herein, we report the first high-dimensional neural network potential (NNP) for this ER reaction based on first-principles calculations including all molecular and surface degrees of freedom. Thanks to the high efficiency of this NNP, we are able to perform extensive quasi-classical molecular dynamics simulations with the inclusion of the excitation of low-lying electron-hole pairs (EHPs), which generally yield good agreement with various experimental results. More importantly, the isotopic and/or EHP effects in total reaction cross-sections and distributions of the product energy, scattering angle, and individual ro-vibrational states have been more clearly shown and discussed. This study sheds valuable light on this important ER prototype and opens a new avenue for further investigations of ER reactions using various initial conditions, surface temperatures, and coverages in the future.
Collapse
Affiliation(s)
- Lingjun Zhu
- School of Chemistry and Materials Science, 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.
| | - Ce Hu
- School of Chemistry and Materials Science, 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.
| | - Jialu Chen
- Department of Physics, City University of Hong Kong, Hong Kong, SAR, People's Republic of China
| | - Bin Jiang
- School of Chemistry and Materials Science, 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
|
19
|
Smits B, Somers MF. The quantum dynamics of H 2 on Cu(111) at a surface temperature of 925 K: Comparing state-of-the-art theory to state-of-the-art experiments 2. J Chem Phys 2023; 158:014704. [PMID: 36610948 DOI: 10.1063/5.0134817] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
State-of-the-art 6D quantum dynamics simulations for the dissociative chemisorption of H2 on a thermally distorted Cu(111) surface, using the static corrugation model, were analyzed to produce several (experimentally available) observables. The expected error, especially important for lower reaction probabilities, was quantified using wavepackets on several different grids as well as two different analysis approaches to obtain more accurate results in the region where a slow reaction channel was experimentally shown to be dominant. The lowest reaction barrier sites for different thermally distorted surface slabs are shown to not just be energetically, but also geometrically, different between surface configurations, which can be used to explain several dynamical effects found when including surface temperature effects. Direct comparison of simulated time-of-flight spectra to those obtained from state-of-the-art desorption experiments showed much improved agreement compared to the perfect lattice BOSS approach. Agreement with experimental rotational and vibrational efficacies also somewhat improved when thermally excited surfaces were included in the theoretical model. Finally, we present clear quantum effects in the rotational quadrupole alignment parameters found for the lower rotationally excited states, which underlines the importance of careful quantum dynamical analyses of this system.
Collapse
Affiliation(s)
- B Smits
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Building, 2300 RA Leiden, The Netherlands
| | - M F Somers
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Building, 2300 RA Leiden, The Netherlands
| |
Collapse
|
20
|
Liu T, Shi H, Fu B, Zhang DH. Quantum dynamics reveal different ligand effects by vibrational excitation in the dissociative chemisorption of HCl on the Au/Ag(111) surface. J Chem Phys 2022; 157:244702. [PMID: 36586991 DOI: 10.1063/5.0131503] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
The reactivity and selectivity of bimetallic surfaces are of fundamental importance in industrial applications. Here, we report the first six-dimensional (6D) quantum dynamics study for the role of surface strain and ligand effects on the reactivity of HCl on a strained pseudomorphic monolayer of Au deposited onto a Ag(111) substrate, with the aid of accurate machine learning-based potential energy surfaces. The substitute of Au into Ag changes the location of the transition state; however, the static barrier height remains roughly the same as pure Au(111). The 6D quantum dynamics calculations reveal that the surface strain due to lattice expansion slightly enhances the reactivity. The ligand effect due to electronic structure interactions between Au and Ag substantially suppresses the reactivity of HCl in the ground vibrational state but promotes the reactivity via vibrational excitation at high kinetic energies. This finding can be attributed to more close interaction with Ag atoms at the transition state close to the fcc site, as well as the tight transition-state region, making the vibrational excitation highly efficient in enhancing the reactivity. Our study quantitatively unravels the dynamical origin of reactivity control by two metals, which will ultimately provide valuable insight into the selectivity of the catalyst.
Collapse
Affiliation(s)
- Tianhui Liu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Huixia Shi
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| |
Collapse
|
21
|
Tchakoua T, Gerrits N, Smeets EWF, Kroes GJ. SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces. J Chem Theory Comput 2022; 19:245-270. [PMID: 36529979 PMCID: PMC9835835 DOI: 10.1021/acs.jctc.2c00824] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.
Collapse
Affiliation(s)
- T. Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands
| | - N. Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,PLASMANT,
Department of Chemistry, University of Antwerp, BE-2610Antwerp, Belgium
| | - E. W. F. Smeets
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,ALTEN
Nederland, Technology, Fascinatio Boulevard 582, 2909 VACapelle a/d IJssel, The Netherlands
| | - G.-J. Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,E-mail: . Phone: +31 71 527 4396
| |
Collapse
|
22
|
Koval NE, Sánchez-Portal D, Borisov AG, Díez Muiño R. Time-dependent density functional theory calculations of electronic friction in non-homogeneous media. Phys Chem Chem Phys 2022; 24:20239-20248. [PMID: 35996966 DOI: 10.1039/d2cp01972h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The excitation of low-energy electron-hole pairs is one of the most relevant processes in the gas-surface interaction. An efficient tool to account for these excitations in simulations of atomic and molecular dynamics at surfaces is the so-called local density friction approximation (LDFA). The LDFA is based on a strong approximation that simplifies the dynamics of the electronic system: a local friction coefficient is defined using the value of the electronic density for the unperturbed system at each point of the dynamics. In this work, we apply real-time time-dependent density functional theory to the problem of the electronic friction of a negative point charge colliding with spherical jellium metal clusters. Our non-adiabatic, parameter-free results provide a benchmark for the widely used LDFA approximation and allow the discussion of various processes relevant to the electronic response of the system in the presence of the projectile.
Collapse
Affiliation(s)
- Natalia E Koval
- Centro de Física de Materiales (CFM) CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain. .,Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, E-20018 San Sebastián, Spain.,CIC Nanogune BRTA, Tolosa Hiribidea 76, E-20018 San Sebastián, Spain
| | - Daniel Sánchez-Portal
- Centro de Física de Materiales (CFM) CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain. .,Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, E-20018 San Sebastián, Spain
| | - Andrei G Borisov
- Institut des Sciences Moléculaires d'Orsay (ISMO), UMR 8214, CNRS-Université Paris-Saclay, Bât. 520, F-91405 Orsay CEDEX, France
| | - Ricardo Díez Muiño
- Centro de Física de Materiales (CFM) CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain. .,Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, E-20018 San Sebastián, Spain
| |
Collapse
|
23
|
Meng Q, Chen J, Ma J, Zhang X, Chen J. Adiabatic models for the quantum dynamics of surface scattering with lattice effects. Phys Chem Chem Phys 2022; 24:16415-16436. [PMID: 35766107 DOI: 10.1039/d2cp01560a] [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/21/2022]
Abstract
In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Junbo Chen
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China. .,Xi'an Modern Chemistry Research Institute, China North Industries Group Corp., Ltd., East Zhangba Road 168, 710065 Xi'an, China
| | - Jianxing Ma
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Jun Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Yangqiao Road West 155, 350002 Fuzhou, China.,Fujian Science & Technology Innovation Laboratory for Optoelectronic Information of China, Optoelectronic Industry Base at High-tech Zone, 350108 Fuzhou, China
| |
Collapse
|
24
|
Díaz C, Gravielle MS. Grazing incidence fast atom and molecule diffraction: theoretical challenges. Phys Chem Chem Phys 2022; 24:15628-15656. [PMID: 35730987 DOI: 10.1039/d2cp01246d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This perspective article reviews the state-of-the-art of grazing incidence fast atom and molecule diffraction (GIFAD and GIFMD) simulations and addresses the main challenges that theorists, aiming to provide useful inputs in this topic, are facing. We first discuss briefly the methods used to build accurate potential energy surfaces describing the interaction between the projectile and the surface. Subsequently, we focus on the dynamics simulation methods for GIFAD, a phenomenon that has received a lot of experimental attention since 2007, when the first measurements were published. Following this experimental effort, theorists have developed and adapted a bunch of methods able to simulate, analyze and extract information from the experimental outputs. We review these methods, from the very simple ones based on classical dynamics to the full quantum ones, paying special attention to more versatile semiclassical approaches, which include quantum ingredients in the dynamics at a computational cost only slightly higher than that required in classical dynamics. Within the semiclassical framework it is possible, for example, to include in the dynamics the surface phonons and the projectile coherence, two factors that may have a relevant influence on the experimental measurements, at a reasonable computational cost. Finally, we address GIFMD, a phenomenon that has received much less attention and for which there is still a lot of room for research. We review the few examples of GIFMD available in the literature, and we discuss new phenomena associated with the molecular internal degrees of freedom, which may have some impact in other closely related fields, such as molecular reactivity on metal surfaces. Finally, we point out opened questions, raised from the comparisons between theoretical and experimental results, which claim for further experimental efforts.
Collapse
Affiliation(s)
- Cristina Díaz
- Departamento de Química Física, Facultad de CC. Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | - María Silvia Gravielle
- Instituto de Astronomía y Física del Espacio (IAFE, UBA-CONICET), Ciudad Universitaria, C1428EGA, Buenos Aires, Argentina.
| |
Collapse
|
25
|
Töpfer K, Upadhyay M, Meuwly M. Quantitative molecular simulations. Phys Chem Chem Phys 2022; 24:12767-12786. [PMID: 35593769 PMCID: PMC9158373 DOI: 10.1039/d2cp01211a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 04/30/2022] [Indexed: 11/21/2022]
Abstract
All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.
Collapse
Affiliation(s)
- Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| |
Collapse
|
26
|
Smits B, Litjens LG, Somers MF. Accurate Description of the Quantum Dynamical Surface Temperature Effects on the Dissociative Chemisorption of H 2 from Cu(111). J Chem Phys 2022; 156:214706. [DOI: 10.1063/5.0094985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurately describing surface temperature effects for the dissociative scattering of H2 on a metal surface on a quantum dynamical level is currently one of the open challenges for theoretical surface scientists. We present the first quantum dynamical (QD) simulations of hydrogen dissociating on a Cu(111) surface which accurately describe all relevant surface temperature effects, using the static corrugation model (SCM). The reaction probabilities we obtain show very good agreement with those found using quasi-classical dynamics (QCD), both for individual surface slabs and for an averaged, thus Monte-Carlo sampled, set of thermally distorted surface configurations. Rovibrationally elastic scattering probabilities show a much clearer difference between the QCD and QD results, which appears to be traceable back towards thermally distorted surface configurations with very low dissociation probabilities and underlines the importance of investigating more observables than just dissociation. By reducing the number of distorted surface atoms included in the dynamical model, we also show that only including one, or even three, surface atoms is generally not enough to accurately describe the effects of surface temperature on dissociation and elastic scattering. These results are a major step forward in accurately describing hydrogen scattering from a thermally excited Cu(111) surface, and open up a pathway to better describe reaction and scattering from other relevant crystal facets, such as stepped surfaces, at moderately elevated surface temperatures where quantum effects are expected to play a more important role in the dissociation of H2 on Cu.
Collapse
Affiliation(s)
- Bauke Smits
- Theoretical Chemistry, Leiden University Institute of Chemistry, Netherlands
| | | | - Mark F Somers
- Leiden Institute of Chemistry, Leiden University, Netherlands
| |
Collapse
|
27
|
Sundararaman R, Vigil-Fowler D, Schwarz K. Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface. Chem Rev 2022; 122:10651-10674. [PMID: 35522135 PMCID: PMC10127457 DOI: 10.1021/acs.chemrev.1c00800] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.
Collapse
Affiliation(s)
- Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th Street, Troy, New York 12180, United States
| | - Derek Vigil-Fowler
- Materials, Chemical, and Computational Science Directorate, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
28
|
Chadwick H, Somers MF, Stewart AC, Alkoby Y, Carter TJD, Butkovicova D, Alexandrowicz G. Stopping molecular rotation using coherent ultra-low-energy magnetic manipulations. Nat Commun 2022; 13:2287. [PMID: 35484103 PMCID: PMC9050693 DOI: 10.1038/s41467-022-29830-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2021] [Accepted: 03/31/2022] [Indexed: 11/21/2022] Open
Abstract
Rotational motion lies at the heart of intermolecular, molecule-surface chemistry and cold molecule science, motivating the development of methods to excite and de-excite rotations. Existing schemes involve perturbing the molecules with photons or electrons which supply or remove energy comparable to the rotational level spacing. Here, we study the possibility of de-exciting the molecular rotation of a D2 molecule, from J = 2 to the non-rotating J = 0 state, without using an energy-matched perturbation. We show that passing the beam through a 1 m long magnetic field, which splits the rotational projection states by only 10−12 eV, can change the probability that a molecule-surface collision will stop a molecule from rotating and lose rotational energy which is 9 orders larger than that of the magnetic manipulation. Calculations confirm that different rotational orientations have different de-excitation probabilities but underestimate rotational flips (∆mJ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ne$$\end{document}≠0), highlighting the importance of the results as a sensitive benchmark for further developing theoretical models of molecule-surface interactions. Manipulating the rotational motions of molecules may provide a tool for controlling chemical processes. Here the authors demonstrate that the rotation of a D2 molecule can be stopped, upon collision with a metal surface, by a magnetic field that affects the rotational levels to a much smaller extent than the energy difference upon de-excitation.
Collapse
Affiliation(s)
- Helen Chadwick
- Department of Chemistry, Faculty of Science and Engineering, Swansea University, Swansea, SA2 8PP, UK.
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, PO Box 9502, 2300 RA, Leiden, The Netherlands
| | - Aisling C Stewart
- Department of Chemistry, Faculty of Science and Engineering, Swansea University, Swansea, SA2 8PP, UK
| | - Yosef Alkoby
- Department of Chemistry, Faculty of Science and Engineering, Swansea University, Swansea, SA2 8PP, UK
| | - Thomas J D Carter
- Department of Chemistry, Faculty of Science and Engineering, Swansea University, Swansea, SA2 8PP, UK
| | - Dagmar Butkovicova
- Department of Chemistry, Faculty of Science and Engineering, Swansea University, Swansea, SA2 8PP, UK
| | - Gil Alexandrowicz
- Department of Chemistry, Faculty of Science and Engineering, Swansea University, Swansea, SA2 8PP, UK.
| |
Collapse
|
29
|
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
|
30
|
Urrego-Ortiz R, Builes S, Calle-Vallejo F. Impact of Intrinsic Density Functional Theory Errors on the Predictive Power of Nitrogen Cycle Electrocatalysis Models. ACS Catal 2022; 12:4784-4791. [PMID: 35465243 PMCID: PMC9017217 DOI: 10.1021/acscatal.1c05333] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Indexed: 01/07/2023]
Affiliation(s)
- Ricardo Urrego-Ortiz
- Departamento de Ingeniería de Procesos, Universidad EAFIT, Carrera 49 No 7 sur 50, 050022 Medellín, Colombia
| | - Santiago Builes
- Departamento de Ingeniería de Procesos, Universidad EAFIT, Carrera 49 No 7 sur 50, 050022 Medellín, Colombia
| | - Federico Calle-Vallejo
- Department of Materials Science and Chemical Physics & Institute of Theoretical and Computational Chemistry (IQTCUB), University of Barcelona, C/Martí i Franquès 1, 08028 Barcelona, Spain
| |
Collapse
|
31
|
Shi H, Liu T, Fu Y, Wu H, Fu B, Zhang DH. Fundamental invariant-neural network potential energy surface and dissociative chemisorption dynamics of N2 on rigid Ni(111). COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
32
|
Lončarić I, Alducin M, Juaristi JI. O2 on Ag(110): A puzzle for exchange-correlation functionals. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
33
|
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]
|
34
|
Xue M, Jia J, Wu H. A density functional theory study on the catalytic performance of metal (Ni, Pd) single atom, dimer and trimer for H2 dissociation. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111336] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
35
|
Shen Q, Wu J, Zhou F, Song Y, Dong W, Wang X, Wang T, Yang X. A molecular beam-surface apparatus for quantum state-resolved adsorption studies. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2022; 93:013201. [PMID: 35104941 DOI: 10.1063/5.0049178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 12/07/2021] [Indexed: 06/14/2023]
Abstract
Understanding the microscopic mechanism of molecule-surface interaction is of great importance in the study of chemical dynamics. Yet, it remains challenging to experimentally acquire quantum state resolved results, particularly the results related to different degrees of freedom of the reactants. Here, we report the design and performance of a new apparatus for molecule-surface dynamics studies, which enable the measurement of quantum state-resolved adsorption. A continuous narrow-linewidth infrared laser source and molecular beam unit were developed and employed in this new apparatus to achieve independent control on different degrees of freedom (translation, vibration, and rotation) of the molecule. Preliminary results on hydrogen and hydrogen chloride adsorption on the Cu (111) surface were also presented.
Collapse
Affiliation(s)
- Qiqi Shen
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning 116023, China
| | - Jiawei Wu
- Department of Chemical Physics, University of Science and Technology of China, Hefei 230026, China
| | - Feiyue Zhou
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning 116023, China
| | - Yunlong Song
- Department of Chemistry, College of Science, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Wenrui Dong
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning 116023, China
| | - Xingan Wang
- Department of Chemical Physics, University of Science and Technology of China, Hefei 230026, China
| | - Tao Wang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning 116023, China
| | - Xueming Yang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning 116023, China
| |
Collapse
|
36
|
Six-dimensional State-to-State Quantum Dynamics of H 2/D 2 Scattering from Cu(100): Validity of the Site-Averaging Model. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
37
|
Wei F, Smeets EWF, Voss J, Kroes GJ, Lin S, Guo H. Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Fenfei Wei
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Egidius W. F. Smeets
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Johannes Voss
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park CA 94025, USA
| | - Geert-Jan Kroes
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Sen Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
38
|
Etim UJ, Bai P, Gazit OM, Zhong Z. Low-Temperature Heterogeneous Oxidation Catalysis and Molecular Oxygen Activation. CATALYSIS REVIEWS 2021. [DOI: 10.1080/01614940.2021.1919044] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Ubong J. Etim
- Department of Chemical Engineering, Guangdong Technion-Israel Institute of Technology (GTIIT), Shantou, Guangdong, China
| | - Peng Bai
- College of Chemical Engineering, China University of Petroleum, Qingdao, China
| | - Oz M. Gazit
- Wolfson Faculty of Chemical Engineering, Technion – Israel Institute of Technology, Haifa, Israel
| | - Ziyi Zhong
- Department of Chemical Engineering, Guangdong Technion-Israel Institute of Technology (GTIIT), Shantou, Guangdong, China
- Technion Israel Institute of Technology (IIT), Haifa, Israel
| |
Collapse
|
39
|
Bracco G, Vattuone L, Smerieri M, Carraro G, Savio L, Paolini G, Benedek G, Echenique PM, Rocca M. Prominence of Terahertz Acoustic Surface Plasmon Excitation in Gas-Surface Interaction with Metals. J Phys Chem Lett 2021; 12:9894-9898. [PMID: 34609889 DOI: 10.1021/acs.jpclett.1c02669] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The current understanding of the dynamics of gas-surface interactions is that all of the energy lost in the collision is transferred to vibrations of the target. Electronic excitations were shown to play a marginal role except for cases in which the impinging particles have energies of several electronvolts. Here we show that this picture does not hold for metal surfaces supporting acoustic surface plasmons. Such loss, dressed with a vibronic structure, is shown to make up a prominent energy transfer route down to the terahertz region for Ne atoms scattering off Cu(111) and is expected to dominate for most metals. This mechanism determines, e.g., the drag force acting on telecommunication satellites, which are typically gold-plated to reduce overheating by sunshine. The electronic excitations can be unambiguously discerned from the vibrational ones under mild hyperthermal impact conditions.
Collapse
Affiliation(s)
- G Bracco
- Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, 16146 Genova, Italy
- IMEM-CNR Unità di Genova, Via Dodecaneso 33, 16146 Genova, Italy
| | - L Vattuone
- Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, 16146 Genova, Italy
- IMEM-CNR Unità di Genova, Via Dodecaneso 33, 16146 Genova, Italy
| | - M Smerieri
- IMEM-CNR Unità di Genova, Via Dodecaneso 33, 16146 Genova, Italy
| | - G Carraro
- Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, 16146 Genova, Italy
- IMEM-CNR Unità di Genova, Via Dodecaneso 33, 16146 Genova, Italy
| | - L Savio
- IMEM-CNR Unità di Genova, Via Dodecaneso 33, 16146 Genova, Italy
| | - G Paolini
- Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, 16146 Genova, Italy
| | - G Benedek
- Dipartimento di Scienza dei Materiali, Università di Milano Bicocca, Via R. Cozzi 55, 20125 Milano, Italy
- DIPC, Paseo Manuel de Lardizabal 4, 20018 Donostia, Spain
| | - P M Echenique
- DIPC, Paseo Manuel de Lardizabal 4, 20018 Donostia, Spain
| | - M Rocca
- Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, 16146 Genova, Italy
- IMEM-CNR Unità di Genova, Via Dodecaneso 33, 16146 Genova, Italy
| |
Collapse
|
40
|
Gu K, Wei F, Cai Y, Lin S, Guo H. Dynamics of Initial Hydrogen Spillover from a Single Atom Platinum Active Site to the Cu(111) Host Surface: The Impact of Substrate Electron-Hole Pairs. J Phys Chem Lett 2021; 12:8423-8429. [PMID: 34436916 DOI: 10.1021/acs.jpclett.1c02019] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The initial impulsive diffusion of hot hydrogen atoms resulted from the dissociative chemisorption of H2 at atomically dispersed Pt atoms embedded in Cu(111) is investigated using ab initio molecular dynamics. Upon dissociation, one of the two hydrogen atoms tends to roam away from the dissociation site while the other remains trapped. It is shown that the fraction of diffusion and the average diffusion length increase with the incident energy and H2 vibrational excitation, due apparently to the increased initial kinetic energy of the hot atoms. Most importantly, the strong interaction with surface electron-hole pairs, modeled using an electronic friction model, is shown to play an important role in rapid energy dissipation and significant retardation of the impulsive diffusion.
Collapse
Affiliation(s)
- Kaixuan Gu
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Fenfei Wei
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Yuhui Cai
- 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
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
41
|
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
|
42
|
Janesko BG. Replacing hybrid density functional theory: motivation and recent advances. Chem Soc Rev 2021; 50:8470-8495. [PMID: 34060549 DOI: 10.1039/d0cs01074j] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Density functional theory (DFT) is the most widely-used electronic structure approximation across chemistry, physics, and materials science. Every year, thousands of papers report hybrid DFT simulations of chemical structures, mechanisms, and spectra. Unfortunately, hybrid DFT's accuracy is ultimately limited by tradeoffs between over-delocalization and under-binding. This review summarizes these tradeoffs, and introduces six modern attempts to go beyond them while maintaining hybrid DFT's relatively low computational cost: DFT+U, self-interaction corrections, localized orbital scaling corrections, local hybrid functionals, real-space nondynamical correlation, and our rung-3.5 approach. The review concludes with practical suggestions for DFT users to identify and mitigate these tradeoffs' impact on their simulations.
Collapse
Affiliation(s)
- Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, 2800 S. University Dr, Fort Worth, TX 76129, USA.
| |
Collapse
|
43
|
Potapov A, McCoustra M. Physics and chemistry on the surface of cosmic dust grains: a laboratory view. INT REV PHYS CHEM 2021. [DOI: 10.1080/0144235x.2021.1918498] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Alexey Potapov
- Laboratory Astrophysics Group of the Max Planck Institute for Astronomy at the Friedrich Schiller University Jena, Jena, Germany
| | - Martin McCoustra
- Institute of Chemical Sciences, Heriot-Watt University, Edinburgh, UK
| |
Collapse
|
44
|
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
|
45
|
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
|
46
|
Lushchikova OV, Tahmasbi H, Reijmer S, Platte R, Meyer J, Bakker JM. IR Spectroscopic Characterization of H 2 Adsorption on Cationic Cu n+ ( n = 4-7) Clusters. J Phys Chem A 2021; 125:2836-2848. [PMID: 33787276 PMCID: PMC8054246 DOI: 10.1021/acs.jpca.0c11527] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
![]()
IR spectra of cationic
copper clusters Cun+ (n = 4–7) complexed with hydrogen molecules are recorded
via IR multiple-photon
dissociation (IRMPD) spectroscopy. To this end, the copper clusters
are generated via laser ablation and reacted with H2 and
D2 in a flow-tube-type reaction channel. The complexes
formed are irradiated using IR light provided by the free-electron
laser for intracavity experiments (FELICE). The spectra are interpreted
by making use of isotope-induced shifts of the vibrational bands and
by comparing them to density functional theory calculated spectra
for candidate structures. The structural candidates have been obtained
from global sampling with the minima hopping method, and spectra are
calculated at the semilocal (PBE) and hybrid (PBE0) functional level.
The highest-quality spectra have been recorded for [5Cu, 2H/2D]+, and we find that the semilocal functional provides better
agreement for the lowest-energy isomers. The interaction of hydrogen
with the copper clusters strongly depends on their size. Binding energies
are largest for Cu5+, which goes hand in hand
with the observed predominantly dissociative adsorption. Due to smaller
binding energies for dissociated H2 and D2 for
Cu4+, also a significant amount of molecular
adsorption is observed as to be expected according to the Evans–Polanyi
principle. This is confirmed by transition-state calculations for
Cu4+ and Cu5+, which show
that hydrogen dissociation is not hindered by an endothermic reaction
barrier for Cu5+ and by a slightly endothermic
barrier for Cu4+. For Cu6+ and Cu7+, it was difficult to draw clear conclusions
because the IR spectra could not be unambiguously assigned to structures.
Collapse
Affiliation(s)
- Olga V Lushchikova
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands
| | - Hossein Tahmasbi
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Stijn Reijmer
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands
| | - Rik Platte
- 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
| | - Joost M Bakker
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands
| |
Collapse
|
47
|
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
|
48
|
Fallaque JG, Ramos M, Busnengo HF, Martín F, Díaz C. Normal and off-normal incidence dissociative dynamics of O 2(v,J) on ultrathin Cu films grown on Ru(0001). Phys Chem Chem Phys 2021; 23:7768-7776. [PMID: 33000830 DOI: 10.1039/d0cp03979a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dissociative adsorption of molecular oxygen on metal surfaces has long been controversial, mostly due to the spin-triplet nature of its ground state, to possible non-adiabatic effects, such as an abrupt charge transfer from the metal to the molecule, or even to the role played by the surface electronic state. Here, we have studied the dissociative adsorption of O2 on CuML/Ru(0001) at normal and off-normal incidence, from thermal to super-thermal energies, using quasi-classical dynamics, in the framework of the generalized Langevin oscillator model, and density functional theory based on a multidimensional potential energy surface. Our simulations reveal a rather intriguing behavior of dissociative adsorption probabilities, which exhibit normal energy scaling at incidence energies below the reaction barriers and total energy scaling above, irrespective of the reaction channel, either direct dissociation, trapping dissociation, or molecular adsorption. We directly compare our results with existing scanning tunneling spectroscopy and microscopy measurements. From this comparison, we infer that the observed experimental behavior at thermal energies may be due to ligand and strain effects, as already found for super-thermal incidence energies.
Collapse
Affiliation(s)
- J G Fallaque
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | | | | | | | | |
Collapse
|
49
|
Rivero Santamaría A, Ramos M, Alducin M, Busnengo HF, Díez Muiño R, Juaristi JI. High-Dimensional Atomistic Neural Network Potential to Study the Alignment-Resolved O 2 Scattering from Highly Oriented Pyrolytic Graphite. J Phys Chem A 2021; 125:2588-2600. [PMID: 33734696 DOI: 10.1021/acs.jpca.1c00835] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A high dimensional and accurate atomistic neural network potential energy surface (ANN-PES) that describes the interaction between one O2 molecule and a highly oriented pyrolytic graphite (HOPG) surface has been constructed using the open-source package (aenet). The validation of the PES is performed by paying attention to static characteristics as well as by testing its performance in reproducing previous ab initio molecular dynamics simulation results. Subsequently, the ANN-PES is used to perform quasi-classical molecular dynamics calculations of the alignment-dependent scattering of O2 from HOPG. The results are obtained for 200 meV O2 molecules with different initial alignments impinging with a polar incidence angle with respect to the surface normal of 22.5° on a thermalized (110 and 300 K) graphite surface. The choice of these initial conditions in our simulations is made to perform comparisons to recent experimental results on this system. Our results show that the scattering of O2 from the HOPG surface is a rather direct process, that the angular distributions are alignment dependent, and that the final translational energy of end-on molecules is around 20% lower than that of side-on molecules. Upon collision with the surface, the molecules that are initially aligned perpendicular to the surface become highly rotationally excited, whereas a very small change in the rotational state of the scattered molecules is observed for the initial parallel alignments. The latter confirms the energy transfer dependence on the stereodynamics for the present system. The results of our simulations are in overall agreement with the experimental observations regarding the shape of the angular distributions and the alignment dependence of the in-plane reflected molecules.
Collapse
Affiliation(s)
- Alejandro Rivero Santamaría
- 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
| | - Maximiliano Ramos
- Instituto de Física Rosario, CONICET and Universidad Nacional de Rosario, Bv. 27 de Febrero 210 bis, 2000 Rosario, Argentina.,Facultad de Ciencias Exactas, Ingeniera y Agrimensura, Universidad Nacional de Rosario, Av. Pellegrini 250, S2000BTP Rosario, Argentina
| | - 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
| | - Heriberto Fabio Busnengo
- Instituto de Física Rosario, CONICET and Universidad Nacional de Rosario, Bv. 27 de Febrero 210 bis, 2000 Rosario, Argentina
| | - Ricardo Díez Muiño
- 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
| | - J Iñaki Juaristi
- 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.,Departamento de Polímeros y Materiales Avanzados: Física, Química y Tecnología, Facultad de Químicas, Universidad del País Vasco (UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain
| |
Collapse
|
50
|
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
|