1
|
Farahvash A, Willard AP. A theory of phonon-induced friction on molecular adsorbates. Proc Natl Acad Sci U S A 2024; 121:e2400589121. [PMID: 39052839 PMCID: PMC11295025 DOI: 10.1073/pnas.2400589121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 06/07/2024] [Indexed: 07/27/2024] Open
Abstract
In this manuscript, we provide a general theory for how surface phonons couple to molecular adsorbates. Our theory maps the extended dynamics of a surface's atomic vibrational motions to a generalized Langevin equation, and by doing so captures these dynamics in a single quantity: the non-Markovian friction. The different frequency components of this friction are the phonon modes of the surface slab weighted by their coupling to the adsorbate degrees of freedom. Using this formalism, we demonstrate that physisorbed species couple primarily to acoustic phonons while chemisorbed species couple to dispersionless local vibrations. We subsequently derive equations for phonon-adjusted reaction rates using transition state theory and demonstrate that these corrections improve agreement with experimental results for CO desorption rates from Pt(111).
Collapse
Affiliation(s)
- Ardavan Farahvash
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA02139
| | - Adam P. Willard
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA02139
| |
Collapse
|
2
|
Sah MK, Naskar K, Adhikari S, Smits B, Meyer J, Somers MF. On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K. J Chem Phys 2024; 161:014306. [PMID: 38953445 DOI: 10.1063/5.0217639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2024] [Accepted: 06/13/2024] [Indexed: 07/04/2024] Open
Abstract
We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.
Collapse
Affiliation(s)
- Mantu Kumar Sah
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Bauke Smits
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
3
|
Farahvash A, Agrawal M, Peterson AA, Willard AP. Modeling Surface Vibrations and Their Role in Molecular Adsorption: A Generalized Langevin Approach. J Chem Theory Comput 2023; 19:6452-6460. [PMID: 37682532 DOI: 10.1021/acs.jctc.3c00473] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
The atomic vibrations of a solid surface can play a significant role in the reactions of surface-bound molecules, as well as their adsorption and desorption. Relevant phonon modes can involve the collective motion of atoms over a wide array of length scales. In this paper, we demonstrate how the generalized Langevin equation can be utilized to describe these collective motions weighted by their coupling to individual sites. Our approach builds upon the generalized Langevin oscillator (GLO) model originally developed by Tully. We extend the GLO by deriving parameters from atomistic simulation data. We apply this approach to study the memory kernel of a model platinum surface and demonstrate that the memory kernel has a bimodal form due to coupling to both low-energy acoustic modes and high-energy modes near the Debye frequency. The same bimodal form was observed across a wide variety of solids of different elemental compositions, surface structures, and solvation states. By studying how these dominant modes depend on the simulation size, we argue that the acoustic modes are frozen in the limit of macroscopic lattices. By simulating periodically replicated slabs of various sizes, we quantify the influence of phonon confinement effects in the memory kernel and their concomitant effect on simulated sticking coefficients.
Collapse
Affiliation(s)
- Ardavan Farahvash
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Mayank Agrawal
- School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| | - Andrew A Peterson
- School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| | - Adam P Willard
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
4
|
Martin-Barrios R, Hertl N, Galparsoro O, Kandratsenka A, Wodtke AM, Larrégaray P. H atom scattering from W(110): A benchmark for molecular dynamics with electronic friction. Phys Chem Chem Phys 2022; 24:20813-20819. [PMID: 36004823 PMCID: PMC9472596 DOI: 10.1039/d2cp01850k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Molecular dynamics with electronic friction (MDEF) at the level of the local density friction approximation (LDFA) has been applied to describe electronically non-adiabatic energy transfer accompanying H atom collisions with many solid metal surfaces. When implemented with full dimensional potential energy and electron density functions, excellent agreement with experiment is found. Here, we compare the performance of a reduced dimensional MDEF approach involving a simplified description of H atom coupling to phonons to that of full dimensional MDEF calculations known to yield accurate results. Both approaches give remarkably similar results for H atom energy loss distributions with a 300 K W(110) surface. At low surface temperature differences are seen; but, quantities like average energy loss are still accurately reproduced. Both models predict similar conditions under which H atoms that have penetrated into the subsurface regions could be observed in scattering experiments. Molecular dynamics with electronic friction (MDEF) at the level of the local density friction approximation (LDFA) has been applied to describe electronically non-adiabatic energy transfer accompanying H atom collisions with many solid metal surfaces.![]()
Collapse
Affiliation(s)
- Raidel Martin-Barrios
- Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR5255, F-33400, France. .,Dynamical processes in Atomic and Molecular Systems (DynAMoS), Facultad de Física, Universidad de la Habana, La Habana, 10400, Cuba
| | - Nils Hertl
- Max-Planck Institut für multidisziplinäre Naturwissenschaften, Am Faßberg 11, Göttingen, Germany. .,Institut für physikalische Chemie, Georg-August-Universität, Tammannstraße 6, Göttingen, Germany
| | - Oihana Galparsoro
- Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia Saila, Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU) Lardizabal Pasealekua 3, 20018, Donostia-San Sebastián, Spain
| | - Alexander Kandratsenka
- Max-Planck Institut für multidisziplinäre Naturwissenschaften, Am Faßberg 11, Göttingen, Germany. .,Institut für physikalische Chemie, Georg-August-Universität, Tammannstraße 6, Göttingen, Germany
| | - Alec M Wodtke
- Max-Planck Institut für multidisziplinäre Naturwissenschaften, Am Faßberg 11, Göttingen, Germany. .,Institut für physikalische Chemie, Georg-August-Universität, Tammannstraße 6, Göttingen, Germany
| | | |
Collapse
|
5
|
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
|
6
|
Kang K, Shakouri K, Kroes GJ, Kleyn AW, Meyer J. Three-dimensional Langevin dynamics of N atom scattering from N-covered Ag(111). Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111557] [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]
|
7
|
Hu C, Lin Q, Guo H, Jiang B. Influence of supercell size on Gas-Surface Scattering: A case study of CO scattering from Au(1 1 1). Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111423] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
8
|
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
|
9
|
Hertl N, Martin-Barrios R, Galparsoro O, Larrégaray P, Auerbach DJ, Schwarzer D, Wodtke AM, Kandratsenka A. Random Force in Molecular Dynamics with Electronic Friction. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2021; 125:14468-14473. [PMID: 34267855 PMCID: PMC8273891 DOI: 10.1021/acs.jpcc.1c03436] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 06/10/2021] [Indexed: 06/13/2023]
Abstract
Originally conceived to describe thermal diffusion, the Langevin equation includes both a frictional drag and a random force, the latter representing thermal fluctuations first seen as Brownian motion. The random force is crucial for the diffusion problem as it explains why friction does not simply bring the system to a standstill. When using the Langevin equation to describe ballistic motion, the importance of the random force is less obvious and it is often omitted, for example, in theoretical treatments of hot ions and atoms interacting with metals. Here, friction results from electronic nonadiabaticity (electronic friction), and the random force arises from thermal electron-hole pairs. We show the consequences of omitting the random force in the dynamics of H-atom scattering from metals. We compare molecular dynamics simulations based on the Langevin equation to experimentally derived energy loss distributions. Despite the fact that the incidence energy is much larger than the thermal energy and the scattering time is only about 25 fs, the energy loss distribution fails to reproduce the experiment if the random force is neglected. Neglecting the random force is an even more severe approximation than freezing the positions of the metal atoms or modelling the lattice vibrations as a generalized Langevin oscillator. This behavior can be understood by considering analytic solutions to the Ornstein-Uhlenbeck process, where a ballistic particle experiencing friction decelerates under the influence of thermal fluctuations.
Collapse
Affiliation(s)
- Nils Hertl
- Max-Planck-Institut
für Biophysikalische Chemie, Am Faßberg 11, 37077 Göttingen, Germany
- Institut
für Physikalische Chemie, Georg-August-Universität Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
| | - Raidel Martin-Barrios
- Université
de Bordeaux, 351 Cours
de la Libération, 33405 Talence, France
- CNRS, 351 Cours de la Libération, 33405 Talence, France
- Universidad
de La Habana, San Lázaro
y L, CP 10400 La
Habana, Cuba
| | - Oihana Galparsoro
- Max-Planck-Institut
für Biophysikalische Chemie, Am Faßberg 11, 37077 Göttingen, Germany
- Institut
für Physikalische Chemie, Georg-August-Universität Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
| | - Pascal Larrégaray
- Université
de Bordeaux, 351 Cours
de la Libération, 33405 Talence, France
- CNRS, 351 Cours de la Libération, 33405 Talence, France
| | - Daniel J. Auerbach
- Max-Planck-Institut
für Biophysikalische Chemie, Am Faßberg 11, 37077 Göttingen, Germany
| | - Dirk Schwarzer
- Max-Planck-Institut
für Biophysikalische Chemie, Am Faßberg 11, 37077 Göttingen, Germany
| | - Alec M. Wodtke
- Max-Planck-Institut
für Biophysikalische Chemie, Am Faßberg 11, 37077 Göttingen, Germany
- Institut
für Physikalische Chemie, Georg-August-Universität Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
| | - Alexander Kandratsenka
- Max-Planck-Institut
für Biophysikalische Chemie, Am Faßberg 11, 37077 Göttingen, Germany
| |
Collapse
|
10
|
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
|
11
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (Eb, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate Eb can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for Eb for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH4. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of Eb for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate Eb for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O2, H2O, NH3, CO2, and CH3OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
12
|
Fujiwara K, Shibahara M. Atomic-scale thermal manipulation with adsorbed atoms on a solid surface at a liquid-solid interface. Sci Rep 2019; 9:13202. [PMID: 31519938 PMCID: PMC6744397 DOI: 10.1038/s41598-019-49677-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2019] [Accepted: 08/29/2019] [Indexed: 11/10/2022] Open
Abstract
Modulating thermal transport through interfaces is one of the central issues in nanoscience and nanotechnology. This study examined thermal transport between atoms adsorbed on a solid surface and a liquid phase based on non-equilibrium molecular dynamics. The heat flux was detected at sub-atomic spatial resolution, yielding a two-dimensional map of local heat flux in the vicinity of the adsorbed atoms on the surface. Based on the detected heat flux, the possibility of atomic-scale thermal manipulation with the adsorbed atoms was examined by varying the interaction strengths between the liquid molecules and atoms adsorbed on the surface. The results of the local heat flux at the single-atom scale clearly showed effects of the adsorbed atoms on the thermal transport through the liquid-solid interface; they can significantly enhance the heat flux at the single-atom scale using degrees of freedom normal to the macroscopic temperature gradient. The effect was especially evident for a low wettability surface, which provides key information on local enhancement at the single-atom scale of the thermal transport through a liquid-solid interface.
Collapse
Affiliation(s)
- Kunio Fujiwara
- Center for Atomic and Molecular Technologies, Osaka University, 2-1 Yamadaoka, Suita, Osaka, 565-0871, Japan.
| | - Masahiko Shibahara
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka, 565-0871, Japan
| |
Collapse
|
13
|
Jiang B, Guo H. Dynamics in reactions on metal surfaces: A theoretical perspective. J Chem Phys 2019; 150:180901. [PMID: 31091904 DOI: 10.1063/1.5096869] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.
Collapse
Affiliation(s)
- Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
14
|
Serwatka T, Tremblay JC. Stochastic wave packet approach to nonadiabatic scattering of diatomic molecules from metals. J Chem Phys 2019; 150:184105. [PMID: 31091890 DOI: 10.1063/1.5092698] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this contribution, we present a quantum dynamical approach to study inelastic scattering of diatomic molecules from metal surfaces at normal incidence. The dissipative dynamics obeys a stochastic Schrödinger equation describing the time-evolution of the system as a piecewise deterministic process. Energy exchange between the molecular vibrational degrees of freedom and the metal electrons is represented using operators in tensor product form, which are coupled via anharmonic transition rates calculated from first-order perturbation theory. Full dimensional observables are obtained by averaging over simulations in 4D-including the internal stretch, the distance to the surface, and the orientation angles-at different surface sites. The method is applied to the state-resolved scattering of vibrationally excited NO from Au(111), revealing important channels for quantized energy relaxation.
Collapse
Affiliation(s)
- T Serwatka
- Institut für Chemie und Biochemie, Freie Universität Berlin, D-14195 Berlin, Germany
| | - J C Tremblay
- Institut für Chemie und Biochemie, Freie Universität Berlin, D-14195 Berlin, Germany
| |
Collapse
|
15
|
Zhou X, Jiang B. A modified generalized Langevin oscillator model for activated gas-surface reactions. J Chem Phys 2019; 150:024704. [DOI: 10.1063/1.5078541] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale, 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, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
16
|
Luo X, Zhou X, Jiang B. Effects of surface motion and electron-hole pair excitations in CO2 dissociation and scattering on Ni(100). J Chem Phys 2018; 148:174702. [DOI: 10.1063/1.5025029] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Xuan Luo
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale, 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, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
17
|
Peña-Torres A, Busnengo HF, Juaristi JI, Larregaray P, Crespos C. Dynamics of N2 sticking on W(100): the decisive role of van der Waals interactions. Phys Chem Chem Phys 2018; 20:19326-19331. [DOI: 10.1039/c8cp03515f] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The reactive dynamics of N2 on W(100) has been investigated by means of quasi-classical trajectory calculations using an interpolated six-dimensional potential energy surface (PES) based on density functional theory energies obtained employing the vdW-DF2 functional.
Collapse
Affiliation(s)
| | - H. Fabio Busnengo
- Instituto de Física de Rosario (CONICET-UNR) and Facultad de Ciencias Exactas
- Ingeniería y Agrimensura
- Universidad Nacional de Rosario
- 2000 Rosario
- Argentina
| | - J. Iñaki Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
- Departamento de Física de Materiales
- Facultad de Químicas (UPV/EHU)
| | | | | |
Collapse
|
18
|
Lončarić I, Füchsel G, Juaristi JI, Saalfrank P. Strong Anisotropic Interaction Controls Unusual Sticking and Scattering of CO at Ru(0001). PHYSICAL REVIEW LETTERS 2017; 119:146101. [PMID: 29053313 DOI: 10.1103/physrevlett.119.146101] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/05/2017] [Indexed: 06/07/2023]
Abstract
Complete sticking at low incidence energies and broad angular scattering distributions at higher energies are often observed in molecular beam experiments on gas-surface systems which feature a deep chemisorption well and lack early reaction barriers. Although CO binds strongly on Ru(0001), scattering is characterized by rather narrow angular distributions and sticking is incomplete even at low incidence energies. We perform molecular dynamics simulations, accounting for phononic (and electronic) energy loss channels, on a potential energy surface based on first-principles electronic structure calculations that reproduce the molecular beam experiments. We demonstrate that the mentioned unusual behavior is a consequence of a very strong rotational anisotropy in the molecule-surface interaction potential. Beyond the interpretation of scattering phenomena, we also discuss implications of our results for the recently proposed role of a precursor state for the desorption and scattering of CO from ruthenium.
Collapse
Affiliation(s)
- Ivor Lončarić
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Gernot Füchsel
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - J I Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
- Departamento de Física de Materiales, Facultad de Químicas, Universidad del País Vasco (UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain
- Donostia International Physics Center DIPC, P. Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Strasse 24-25, D-14476 Potsdam, Germany
| |
Collapse
|
19
|
Kroes GJ, Juaristi JI, Alducin M. Vibrational Excitation of H 2 Scattering from Cu(111): Effects of Surface Temperature and of Allowing Energy Exchange with the Surface. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2017; 121:13617-13633. [PMID: 28729891 PMCID: PMC5510092 DOI: 10.1021/acs.jpcc.7b01096] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2017] [Revised: 05/23/2017] [Indexed: 05/11/2023]
Abstract
In scattering of H2 from Cu(111), vibrational excitation has so far defied an accurate theoretical description. To expose the causes of the large discrepancies with experiment, we investigate how the feature due to vibrational excitation (the "gain peak") in the simulated time-of-flight spectrum of (v = 1, j = 3) H2 scattering from Cu(111) depends on the surface temperature (Ts) and the possibility of energy exchange with surface phonons and electron-hole pairs (ehp's). Quasi-classical dynamics calculations are performed on the basis of accurate semiempirical density functionals for the interaction with H2 + Cu(111). The methods used include the quasi-classical trajectory method within the Born-Oppenheimer static surface model, the generalized Langevin oscillator (GLO) method incorporating energy transfer to surface phonons, the GLO + friction (GLO+F) method also incorporating energy exchange with ehp's, and ab initio molecular dynamics with electronic friction (AIMDEF). Of the quasi-classical methods tested, comparison with AIMDEF suggests that the GLO+F method is accurate enough to describe vibrational excitation as measured in the experiments. The GLO+F calculations also suggest that the promoting effect of raising Ts on the measured vibrational excitation is due to an electronically nonadiabatic mechanism. However, by itself, enabling energy exchange with the surface by modeling surface phonons and ehp's leads to reduced vibrational excitation, further decreasing the agreement with experiment. The simulated gain peak is quite sensitive to energy shifts in calculated vibrational excitation probabilities and to shifts in a specific experimental parameter (the chopper opening time). While the GLO+F calculations allow important qualitative conclusions, comparison to quantum dynamics results suggests that, with the quasi-classical way of describing nuclear motion and the present box quantization method for assigning the final vibrational state, the gain peak is not yet described with quantitative accuracy. Ways in which this problem might be resolved in the future are discussed.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- E-mail:
| | - J. I. Juaristi
- Departamento
de Física de Materiales, Facultad de Químicas, Universidad del País Vasco/Euskal Herriko Unibertsitatea
(UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain
- Centro
de Física de Materiales (CFM/MPC), Consejo Superior de Investigaciones Científicas (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
| | - M. Alducin
- Centro
de Física de Materiales (CFM/MPC), Consejo Superior de Investigaciones Científicas (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
| |
Collapse
|
20
|
Strak P, Sakowski K, Kempisty P, Krukowski S. Dissipation of the excess energy of the adsorbate-thermalization via electron transfer. Phys Chem Chem Phys 2017; 19:9149-9155. [PMID: 28318002 DOI: 10.1039/c7cp00235a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new scenario for the thermalization process of adsorbates at solid surfaces is proposed. The scenario is based on the existence of an electric dipole layer in which the electron wavefunctions extend over the positive ions, creating a strong local electric field which drags the electrons into the solid interior and repels the positive ions. During adsorption the electrons tunnel into the solid interior, conveying the excess energy. The positive ions are retarded by the field, losing the excess kinetic energy, and are located smoothly into the adsorption sites. In such a scheme, the excess energy is not dissipated locally, avoiding melting or the creation of defects which is in accordance with experiments. The scenario is supported by ab initio calculation results, including density function theory of the slabs representing the AlN surface and the Schrodinger equation for the time evolution of hydrogen-like atoms at the solid surface.
Collapse
Affiliation(s)
- Pawel Strak
- Institute of High Pressure Physics, Polish Academy of Sciences, Sokołowska 29/37, 01-142 Warsaw, Poland.
| | - Konrad Sakowski
- Institute of High Pressure Physics, Polish Academy of Sciences, Sokołowska 29/37, 01-142 Warsaw, Poland.
| | - Pawel Kempisty
- Institute of High Pressure Physics, Polish Academy of Sciences, Sokołowska 29/37, 01-142 Warsaw, Poland. and Center for Integrated Research of Future Electronics (CIRFE), Institute of Materials and Systems for Sustainability (IMaSS), Nagoya University, Nagoya 464-8603, Japan
| | - Stanislaw Krukowski
- Institute of High Pressure Physics, Polish Academy of Sciences, Sokołowska 29/37, 01-142 Warsaw, Poland.
| |
Collapse
|
21
|
Dou W, Subotnik JE. Electronic friction near metal surfaces: A case where molecule-metal couplings depend on nuclear coordinates. J Chem Phys 2017. [DOI: 10.1063/1.4965823] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
22
|
Nosir MA, Martin-Gondre L, Bocan GA, Díez Muiño R. Dissociative adsorption dynamics of nitrogen on a Fe(111) surface. Phys Chem Chem Phys 2017; 19:24626-24635. [DOI: 10.1039/c7cp03701e] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dissociative adsorption dynamics of N2 on clean Fe(111) surfaces is theoretically investigated by means of quasi-classical trajectory calculations based on a multidimensional potential energy surface built from density functional theory.
Collapse
Affiliation(s)
- M. A. Nosir
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- Donostia-San Sebastián
- Spain
- Donostia International Physics Center (DIPC)
- Donostia-San Sebastián
| | | | - G. A. Bocan
- CONICET and Centro Atómico Bariloche (CNEA)
- Argentina
| | - R. Díez Muiño
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- Donostia-San Sebastián
- Spain
- Donostia International Physics Center (DIPC)
- Donostia-San Sebastián
| |
Collapse
|
23
|
Nosir MA, Martin-Gondre L, Bocan GA, Díez Muiño R. Adsorption dynamics of molecular nitrogen at an Fe(111) surface. Phys Chem Chem Phys 2017; 19:7370-7379. [DOI: 10.1039/c6cp07174k] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present an extensive theoretical study of N2 adsorption mechanisms on an Fe(111) surface.
Collapse
Affiliation(s)
- M. A. Nosir
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- Donostia-San Sebastián
- Spain
- Donostia International Physics Center (DIPC)
- Donostia-San Sebastián
| | - L. Martin-Gondre
- Institut UTINAM – Université Bourgogne Franche-Comté
- Besançon
- France
| | - G. A. Bocan
- CONICET and Centro Atómico Bariloche (CNEA)
- Argentina
| | - R. Díez Muiño
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- Donostia-San Sebastián
- Spain
- Donostia International Physics Center (DIPC)
- Donostia-San Sebastián
| |
Collapse
|
24
|
Ramos M, Díaz C, Martínez AE, Busnengo HF, Martín F. Dissociative and non-dissociative adsorption of O2 on Cu(111) and CuML/Ru(0001) surfaces: adiabaticity takes over. Phys Chem Chem Phys 2017; 19:10217-10221. [DOI: 10.1039/c7cp00753a] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Adiabatic molecular spin-quenching during the approach of O2 to Cu(111) and CuML/Ru(0001) surfaces.
Collapse
Affiliation(s)
- M. Ramos
- Instituto de Fsica Rosario
- CONICET
- and Universidad Nacional de Rosario
- Rosario
- Argentina
| | - C. Díaz
- Departamento de Qumica Módulo 13
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
- Condensed Matter Physics Center (IFIMAC)
| | - A. E. Martínez
- Instituto de Fsica Rosario
- CONICET
- and Universidad Nacional de Rosario
- Rosario
- Argentina
| | - H. F. Busnengo
- Instituto de Fsica Rosario
- CONICET
- and Universidad Nacional de Rosario
- Rosario
- Argentina
| | - F. Martín
- Departamento de Qumica Módulo 13
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
- Condensed Matter Physics Center (IFIMAC)
| |
Collapse
|
25
|
Ma L, Li X, Liu C. The derivation and approximation of coarse-grained dynamics from Langevin dynamics. J Chem Phys 2016; 145:204117. [DOI: 10.1063/1.4967936] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Affiliation(s)
- Lina Ma
- Department of Mathematics, The Pennsylvania State University, University Park, Pennsylvania 16802-6400, USA
| | - Xiantao Li
- Department of Mathematics, The Pennsylvania State University, University Park, Pennsylvania 16802-6400, USA
| | - Chun Liu
- Department of Mathematics, The Pennsylvania State University, University Park, Pennsylvania 16802-6400, USA
| |
Collapse
|
26
|
Ma L, Li X, Liu C. From generalized Langevin equations to Brownian dynamics and embedded Brownian dynamics. J Chem Phys 2016. [DOI: 10.1063/1.4962419] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lina Ma
- Department of Mathematics, The Pennsylvania State University, University Park, Pennsylvania 16802-6400, USA
| | - Xiantao Li
- Department of Mathematics, The Pennsylvania State University, University Park, Pennsylvania 16802-6400, USA
| | - Chun Liu
- Department of Mathematics, The Pennsylvania State University, University Park, Pennsylvania 16802-6400, USA
| |
Collapse
|
27
|
Kajita S. Green's function nonequilibrium molecular dynamics method for solid surfaces and interfaces. Phys Rev E 2016; 94:033301. [PMID: 27739703 DOI: 10.1103/physreve.94.033301] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2016] [Indexed: 06/06/2023]
Abstract
This study presents a comprehensive procedure to calculate the exact dynamic Green's function of a harmonic semi-infinite solid and the time trajectories of the atoms, in the framework of the Green's function molecular dynamics. This Green's function properly describes the energy dissipation caused by excitations of the surface phonons, and the simulated atoms serve as well-defined thermo- and barostats for the nonequilibrium surface and interface systems. Moreover, the use of the exact dynamic Green's function coupled with a fast convolution algorithm significantly improves both the accuracy and the computing speed. The presented method is applied to a diamond (001) surface, and the results demonstrate that the properties of the nonreflecting boundary, the thermal fluctuations, and the energy dissipations involving long-wavelength phonons are correctly reproduced. These distinctive performances potentially allow us to reveal the nonequilibrium phenomena in a wide spectrum of applications such as catalysis, thermal transport, fracture mechanics, mechanochemistry, and tribology.
Collapse
Affiliation(s)
- Seiji Kajita
- Istituto Nanoscienze, CNR-Consiglio Nazionale delle Ricerche, Via Campi 213A, I-41125 Modena, Italy and Toyota Central R&D Labs., Inc., 41-1, Yokomichi, Nagakute, Aichi 480-1192, Japan
| |
Collapse
|
28
|
Ouyang W, Dou W, Jain A, Subotnik JE. Dynamics of Barrier Crossings for the Generalized Anderson–Holstein Model: Beyond Electronic Friction and Conventional Surface Hopping. J Chem Theory Comput 2016; 12:4178-83. [DOI: 10.1021/acs.jctc.6b00533] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Wenjun Ouyang
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amber Jain
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
29
|
Dou W, Subotnik JE. A many-body states picture of electronic friction: The case of multiple orbitals and multiple electronic states. J Chem Phys 2016; 145:054102. [DOI: 10.1063/1.4959604] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
30
|
Luo X, Jiang B, Juaristi JI, Alducin M, Guo H. Electron-hole pair effects in methane dissociative chemisorption on Ni(111). J Chem Phys 2016; 145:044704. [DOI: 10.1063/1.4959288] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xuan Luo
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Bin Jiang
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - J. Iñaki Juaristi
- Centro de Física de Materiales CFM/MPC(CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center DIPC, P. Manuel de Lardizabal 4, 20018 San Sebastián, Spain
- Departamento de Física de Materiales, Facultad de Químicas, Universidad del País Vasco (UPV/EHU), Apartado 1072, 20080 San Sebastián, Spain
| | - Maite Alducin
- Centro de Física de Materiales CFM/MPC(CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center DIPC, P. Manuel de Lardizabal 4, 20018 San Sebastián, Spain
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
31
|
Nattino F, Galparsoro O, Costanzo F, Díez Muiño R, Alducin M, Kroes GJ. Modeling surface motion effects in N2 dissociation on W(110): Ab initio molecular dynamics calculations and generalized Langevin oscillator model. J Chem Phys 2016; 144:244708. [DOI: 10.1063/1.4954773] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Francesco Nattino
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Oihana Galparsoro
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Université de Bordeaux, F-33400 Talence, France
- CNRS, ISM, UMR 5255, F-33400 Talence, France
| | - Francesca Costanzo
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Ricardo Díez Muiño
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Maite Alducin
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
32
|
Jiang B, Alducin M, Guo H. Electron-Hole Pair Effects in Polyatomic Dissociative Chemisorption: Water on Ni(111). J Phys Chem Lett 2016; 7:327-31. [PMID: 26732612 DOI: 10.1021/acs.jpclett.5b02737] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
The influence of electron-hole pairs in dissociative chemisorption of a polyatomic molecule (water) on metal surfaces is assessed for the first time using a friction approach. The atomic local density dependent friction coefficients computed based on a free electron gas embedding model are employed in classical molecular dynamics simulations of the water dissociation dynamics on rigid Ni(111) using a recently developed nine dimensional interaction potential energy surface for the system. The results indicate that nonadiabatic effects are relatively small and they do not qualitatively alter the mode specificity in the dissociation.
Collapse
Affiliation(s)
- Bin Jiang
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
- Department of Chemical Physics, University of Science and Technology of China , Hefei, Anhui 230026, China
| | - Maite Alducin
- Centro de Física de Materiales Centro Mixto, CFM/MPC (CSIC-UPV/EHU), P. Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center DIPC, P. Manuel de Lardizabal 4, 20018 San Sebastián, Spain
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
| |
Collapse
|
33
|
Kroes GJ, Díaz C. Quantum and classical dynamics of reactive scattering of H2 from metal surfaces. Chem Soc Rev 2016; 45:3658-700. [DOI: 10.1039/c5cs00336a] [Citation(s) in RCA: 120] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
State-of-the-art theoretical models allow nowadays an accurate description of H2/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry
- Gorlaeus Laboratories
- Leiden University
- 2300 RA Leiden
- The Netherlands
| | - Cristina Díaz
- Departamento de Química
- Módulo 13
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
| |
Collapse
|
34
|
Lončarić I, Alducin M, Juaristi JI. Molecular dynamics simulation of O2 adsorption on Ag(110) from first principles electronic structure calculations. Phys Chem Chem Phys 2016; 18:27366-27376. [DOI: 10.1039/c6cp05199e] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
State of the art simulations show that the physisorption state could be important for O2/Ag(110) adsorption.
Collapse
Affiliation(s)
- Ivor Lončarić
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
| | - M. Alducin
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
- Donostia International Physics Center DIPC
- 20018 Donostia-San Sebastián
| | - J. I. Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
- Donostia International Physics Center DIPC
- 20018 Donostia-San Sebastián
| |
Collapse
|
35
|
Abstract
When a classical master equation (CME) is used to describe the nonadiabatic dynamics of a molecule at metal surfaces, we show that in the regime of reasonably strong molecule-metal couplings, the CME can be reduced to a Fokker-Planck equation with an explicit form of electronic friction. For a single metal substrate at thermal equilibrium, the electronic friction and random force satisfy the fluctuation-dissipation theorem. When we investigate the time scale for an electron transfer (ET) event between the molecule and metal surface, we find that the ET rates show a turnover effect (just as in Kramer's theory) as a function of frictional damping.
Collapse
Affiliation(s)
- Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Abraham Nitzan
- School of Chemistry, The Sackler Faculty of Science, Tel Aviv University, Tel Aviv 69978, Israel
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
36
|
Angular distributions and rovibrational excitation of N2 molecules recombined on N-covered Ag(111) by the Eley–Rideal mechanism. Catal Today 2015. [DOI: 10.1016/j.cattod.2014.06.028] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
37
|
Dombrowski E, Peterson E, Del Sesto D, Utz A. Precursor-mediated reactivity of vibrationally hot molecules: Methane activation on Ir(111). Catal Today 2015. [DOI: 10.1016/j.cattod.2014.10.025] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
38
|
Fujiwara K, Shibahara M. Local pressure components and interfacial tensions of a liquid film in the vicinity of a solid surface with a nanometer-scale slit pore obtained by the perturbative method. J Chem Phys 2015; 142:094702. [PMID: 25747094 DOI: 10.1063/1.4913495] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
A classical molecular dynamics simulation was conducted for a liquid-solid interfacial system with a nanometer-scale slit pore in order to reveal local thermodynamic states: local pressure components and interfacial tensions of a liquid film in the vicinity of the slit. The simulation also examined the transition mechanism between the two states of the liquid film: (a) liquid film on the slit and (b) liquid film in the slit, based on the local thermodynamic quantities from a molecular point of view. An instantaneous expression of the local pressure components and interfacial tensions, which is based on a volume perturbation, was presented to investigate time-dependent phenomena in molecular dynamics simulations. The interactions between the particles were described by the 12-6 Lennard-Jones potential, and effects of the fluid-solid interaction intensity on the local pressure components and interfacial tensions of the fluid in the vicinity of the slit were examined in detail by the presented perturbative method. The results revealed that the local pressure components tangential to the solid surface in the vicinity of the 1st fluid layer from the solid surface are different in a two dimensional plane, and the difference became pronounced in the vicinity of the corner of the slit, for cases where the fluid-solid interaction intensities are relatively strong. The results for the local interfacial tensions of the fluid inside the slit suggested that the local interfacial tensions in the vicinity of the 2nd and 3rd layers of the solid atoms from the entrance of the slit act as a trigger for the transition between the two states under the influence of a varying fluid-solid interaction.
Collapse
Affiliation(s)
- K Fujiwara
- R & D Department, SCREEN Holdings Co., Ltd., 322 Furukawa-cho, Hazukashi, Fushimi-ku, Kyoto, Kyoto 612-8486, Japan
| | - M Shibahara
- Department of Mechanical Engineering, Graduate School of Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
39
|
Martin-Gondre L, Juaristi JI, Blanco-Rey M, Díez Muiño R, Alducin M. Influence of the van der Waals interaction in the dissociation dynamics of N2 on W(110) from first principles. J Chem Phys 2015; 142:074704. [DOI: 10.1063/1.4908060] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- L. Martin-Gondre
- Institut UTINAM-CNRS UMR 6213, Université de Franche-Comté, 16 Route de Gray, 25030 Besançon Cedex, France
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - J. I. Juaristi
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Departamento de Física de Materiales, Facultad de Químicas (UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - M. Blanco-Rey
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Departamento de Física de Materiales, Facultad de Químicas (UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain
| | - R. Díez Muiño
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - M. Alducin
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| |
Collapse
|
40
|
Goikoetxea I, Alducin M, Díez Muiño R, Juaristi JI. The dynamics of adsorption and dissociation of N2 in a monolayer of iron on W(110). Phys Chem Chem Phys 2015; 17:19432-45. [DOI: 10.1039/c5cp02051d] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Does N2 adsorption increase on strained Fe monolayers?
Collapse
Affiliation(s)
- I. Goikoetxea
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
| | - M. Alducin
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
- Donostia International Physics Center (DIPC)
- 20018 Donostia-San Sebastián
| | - R. Díez Muiño
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
- Donostia International Physics Center (DIPC)
- 20018 Donostia-San Sebastián
| | - J. I. Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
- Donostia International Physics Center (DIPC)
- 20018 Donostia-San Sebastián
| |
Collapse
|
41
|
Sharia O, Holzgrafe J, Park N, Henkelman G. Rare event molecular dynamics simulations of plasma induced surface ablation. J Chem Phys 2014; 141:074706. [DOI: 10.1063/1.4892841] [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] Open
Affiliation(s)
- Onise Sharia
- Department of Chemistry and the Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712-0165, USA
| | - Jeffrey Holzgrafe
- Department of Chemistry and the Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712-0165, USA
| | - Nayoung Park
- Department of Chemistry and the Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712-0165, USA
| | - Graeme Henkelman
- Department of Chemistry and the Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712-0165, USA
| |
Collapse
|
42
|
Goikoetxea I, Juaristi JI, Díez Muiño R, Alducin M. Surface strain improves molecular adsorption but hampers dissociation for N2 on the Fe/W(110) surface. PHYSICAL REVIEW LETTERS 2014; 113:066103. [PMID: 25148338 DOI: 10.1103/physrevlett.113.066103] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2014] [Indexed: 06/03/2023]
Abstract
We compare the adsorption dynamics of N(2) on the unstrained Fe(110) and on a 10% expanded Fe monolayer grown on W(110) by performing classical molecular dynamics simulations that use potential energy surfaces calculated with density functional theory. Our results allow us to understand why, experimentally, the molecular adsorption of N(2) is observed on the strained layer but not on Fe(110). Surprisingly, we also find that while surface strain favors the molecular adsorption of N(2) it seems, on the contrary, to impede the dissociative adsorption. This result contrasts with previous examples for which strain is found to modify equally the energetics of chemisorption and dissociation.
Collapse
Affiliation(s)
- I Goikoetxea
- Humboldt Universität zu Berlin, Institut für Chemie, Unter den Linden 6, D-10009 Berlin, Germany and Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - J I Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain and Departamento de Física de Materiales, Facultad de Químicas (UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain and Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - R 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 and Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - M Alducin
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain and Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| |
Collapse
|
43
|
Fujiwara K, Shibahara M. Local pressure components and interfacial tension at a liquid-solid interface obtained by the perturbative method in the Lennard-Jones system. J Chem Phys 2014; 141:034707. [DOI: 10.1063/1.4890036] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Pétuya R, Larrégaray P, Crespos C, Busnengo HF, Martínez AE. Dynamics of H2 Eley-Rideal abstraction from W(110): Sensitivity to the representation of the molecule-surface potential. J Chem Phys 2014; 141:024701. [DOI: 10.1063/1.4885139] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- R. Pétuya
- Université de Bordeaux, ISM, CNRS UMR 5255, 33405 Talence Cedex, France
- CNRS, ISM, UMR5255, F-33400 Talence, France
| | - P. Larrégaray
- Université de Bordeaux, ISM, CNRS UMR 5255, 33405 Talence Cedex, France
- CNRS, ISM, UMR5255, F-33400 Talence, France
| | - C. Crespos
- Université de Bordeaux, ISM, CNRS UMR 5255, 33405 Talence Cedex, France
- CNRS, ISM, UMR5255, F-33400 Talence, France
| | - H. F. Busnengo
- Instituto de Física Rosario (IFIR) CONICET-UNR. Ocampo y Esmeralda (2000) Rosario, Argentina
| | - A. E. Martínez
- Instituto de Física Rosario (IFIR) CONICET-UNR. Ocampo y Esmeralda (2000) Rosario, Argentina
| |
Collapse
|
45
|
Goikoetxea I, Meyer J, Juaristi JI, Alducin M, Reuter K. Role of physisorption states in molecular scattering: a semilocal density-functional theory study on O2/Ag(111). PHYSICAL REVIEW LETTERS 2014; 112:156101. [PMID: 24785056 DOI: 10.1103/physrevlett.112.156101] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/23/2014] [Indexed: 06/03/2023]
Abstract
We simulate the scattering of O2 from Ag(111) with classical dynamics simulations performed on a six-dimensional potential energy surface calculated within semilocal density-functional theory. The enigmatic experimental trends that originally required the conjecture of two types of repulsive walls, arising from a physisorption and chemisorption part of the interaction potential, are fully reproduced. Given the inadequate description of the physisorption properties in semilocal density-functional theory, our work casts severe doubts on the prevalent notion to use molecular scattering data as indirect evidence for the existence of such states.
Collapse
Affiliation(s)
- I Goikoetxea
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain
| | - J Meyer
- Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstrasse 4, D-85747 Garching, Germany
| | - J I Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain and Departamento de Física de Materiales, Facultad de Químicas, UPV/EHU, Apartado 1072, E-20080 San Sebastián, Spain and Donostia International Physics Center DIPC, Paseo Manuel de Lardizabal 4, E-20018 San Sebastián, Spain
| | - M Alducin
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain and Donostia International Physics Center DIPC, Paseo Manuel de Lardizabal 4, E-20018 San Sebastián, Spain
| | - K Reuter
- Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstrasse 4, D-85747 Garching, Germany
| |
Collapse
|
46
|
Gamallo P, Prats H, Sayós R. ReaxFF molecular dynamics simulations of CO collisions on an O-preadsorbed silica surface. J Mol Model 2014; 20:2160. [PMID: 24633769 DOI: 10.1007/s00894-014-2160-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2013] [Accepted: 01/27/2014] [Indexed: 11/28/2022]
Abstract
A quasiclassical trajectory dynamics study was performed for carbon monoxide collisions over an oxygen preadsorbed β-cristobalite (001) surface. A reactive molecular force field (ReaxFF) was used to model the potential energy surface. The collisions were performed fixing several initial conditions: CO rovibrational states (v = 0-5 and j = 0, 20, 35), collision energies (0.05 ≤ E(col) ≤ 2.5 eV), incident angles (θ(v) = 0°, 45°) and surface temperatures (T(surf) = 300 K, 900 K). The principal elementary processes were the molecular reflection and the non-dissociative molecular adsorption. CO₂ molecules were also formed in minor extension via an Eley-Rideal reaction although some of them were finally retained on the surface. The scattered CO molecules tend to be translationally colder and internally hotter (rotationally and vibrationally) than the initial ones. The present study supports that CO + O(ad) reaction should be less important than O + O(ad) reaction over silica for similar initial conditions of reactants, in agreement with experimental data.
Collapse
Affiliation(s)
- Pablo Gamallo
- Departament de Química Física and Institut de Química Teòrica i Computacional (IQTC-UB), University of Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | | | | |
Collapse
|
47
|
De Luca S, Todd BD, Hansen JS, Daivis PJ. A new and effective method for thermostatting confined fluids. J Chem Phys 2014; 140:054502. [DOI: 10.1063/1.4862544] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
|
48
|
|
49
|
Nakamura Y, Carlson A, Amberg G, Shiomi J. Dynamic wetting at the nanoscale. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:033010. [PMID: 24125347 DOI: 10.1103/physreve.88.033010] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/16/2013] [Indexed: 06/02/2023]
Abstract
Although the capillary spreading of a drop on a dry substrate is well studied, understanding and describing the physical mechanisms that govern the dynamics remain challenging. Here we study the dynamics of spreading of partially wetting nanodroplets by combining molecular dynamics simulations and continuum phase field simulations. The phase field simulations account for all the relevant hydrodynamics, i.e., capillarity, inertia, and viscous stresses. By coordinated continuum and molecular dynamics simulations, the macroscopic model parameters are extracted. For a Lennard-Jones fluid spreading on a planar surface, the liquid slip at the solid substrate is found to be significant, in fact crucial for the motion of the contact line. Evaluation of the different contributions to the energy transfer shows that the liquid slip generates dissipation of the same order as the bulk viscous dissipation or the energy transfer to kinetic energy. We also study the dynamics of spreading on a substrate with a periodic nanostructure. Here it is found that a nanostructure with a length scale commensurate with molecular size completely inhibits the liquid slip. The dynamic spreading is thus about 30% slower on a nanostructured surface compared to one that is atomically smooth.
Collapse
Affiliation(s)
- Yoshinori Nakamura
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | | | | | | |
Collapse
|
50
|
Si C, Wang XD, Chen M, Duan YY. Effect of wall adsorption on the nano-droplet evaporation in a nano-channel: A molecular dynamics investigation. Colloids Surf A Physicochem Eng Asp 2013. [DOI: 10.1016/j.colsurfa.2013.07.008] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|