J Phys Chem C Nanomater Interfaces
November 2024
Using ab initio molecular dynamics (AIMD) simulations, based on density functional theory that also accounts for van der Waals interactions, we study the oxidation of gas-phase CO on MoSe with a Se vacancy and oxygen coverage of 0.125 ML. In the equilibrium configuration, one of the O atoms is adsorbed on the vacancy and the other one atop one Se atom.
View Article and Find Full Text PDFJ Phys Chem Lett
February 2023
Hydrogen molecules dissociate on clean W(110) surfaces. This reaction is progressively inhibited as the tungsten surface is precovered with oxygen. We use density functional theory and ab initio molecular dynamics to rationalize, at the atomic scale, the influence of the adsorbed O atoms on the H dissociation process.
View Article and Find Full Text PDFThe 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.
View Article and Find Full Text PDFThe classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction. This approach was recently applied to the heterogeneous gas-surface reaction between H2 in its internal ground state and Pd(111) surface e.g.
View Article and Find Full Text PDFWhen elementary reactive processes occur at such low energies that only a few states of reactants and/or products are available, quantum effects strongly manifest and the standard description of the dynamics within the classical framework fails. We show here, for H scattering on Pd(111), that by pseudoquantizing in the spirit of Bohr the relevant final actions of the system, along with adequately treating the diffraction-mediated trapping of the incoming wave, classical simulations achieve an unprecedented agreement with state-of-the-art quantum dynamics calculations.
View Article and Find Full Text PDFPhys Chem Chem Phys
September 2017
We study the dissociative adsorption dynamics of N on clean bcc Fe(111) surfaces. We base our theoretical analysis on a multidimensional potential energy surface built from density functional theory. The dissociative sticking probability is computed by means of quasi-classical trajectory calculations.
View Article and Find Full Text PDFWe present an extensive theoretical study of N adsorption mechanisms on an Fe(111) surface. We combine the static analysis of a six-dimensional potential energy surface (6D-PES), based on ab initio density functional theory (DFT) calculations for the system, with quasi-classical trajectory (QCT) calculations to simulate the adsorption dynamics. There are four molecular adsorption states, usually called γ, δ, α, and ε, arising from our DFT calculations.
View Article and Find Full Text PDFAccurately modeling surface temperature and surface motion effects is necessary to study molecule-surface reactions in which the energy dissipation to surface phonons can largely affect the observables of interest. We present here a critical comparison of two methods that allow to model such effects, namely, the ab initio molecular dynamics (AIMD) method and the generalized Langevin oscillator (GLO) model, using the dissociation of N2 on W(110) as a benchmark. AIMD is highly accurate as the surface atoms are explicitly part of the dynamics, but this advantage comes with a large computational cost.
View Article and Find Full Text PDFPhys Chem Chem Phys
July 2015
We study the adsorption dynamics of N2 on an expanded monolayer of Fe grown pseudomorphically on W(110). To this aim we have performed molecular dynamics simulations in a six-dimensional potential energy surface calculated within density functional theory. Our results show that N2 dissociation on this surface is a highly activated process with an energy barrier of around 1.
View Article and Find Full Text PDFUsing ab initio molecular dynamics (AIMD) calculations, we investigate the role of the van der Waals (vdW) interaction in the dissociative adsorption of N2 on W(110). Hitherto, existing classical dynamics calculations performed on six-dimensional potential energy surfaces based on density functional theory (DFT), and the semi-local PW91 and RPBE [Hammer et al. Phys.
View Article and Find Full Text PDFWe 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.
View Article and Find Full Text PDFWe study the dynamics of transient hot H atoms on Pd(100) that originated from dissociative adsorption of H2. The methodology developed here, denoted AIMDEF, consists of ab initio molecular dynamics simulations that include a friction force to account for the energy transfer to the electronic system. We find that the excitation of electron-hole pairs is the main channel for energy dissipation, which happens at a rate that is five times faster than energy transfer into Pd lattice motion.
View Article and Find Full Text PDFWe have studied survival and rotational excitation probabilities of H(2)(v(i) = 1, J(i) = 1) and D(2)(v(i) = 1, J(i) = 2) upon scattering from Cu(111) using six-dimensional (6D) adiabatic (quantum and quasi-classical) and non-adiabatic (quasi-classical) dynamics. Non-adiabatic dynamics, based on a friction model, has been used to analyze the role of electron-hole pair excitations. Comparison between adiabatic and non-adiabatic calculations reveals a smaller influence of non-adiabatic effects on the energy dependence of the vibrational deexcitation mechanism than previously suggested by low-dimensional dynamics calculations.
View Article and Find Full Text PDFPhys Chem Chem Phys
May 2012
We study the adsorption dynamics of N(2) on the Fe(110) surface. Classical molecular dynamics calculations are performed on top of a six-dimensional potential energy surface calculated within density functional theory. Our results show that N(2) dissociation on this surface is a highly activated process that takes place along a very narrow reaction path with an energy barrier of around 1.
View Article and Find Full Text PDFWe investigate the role played by electron-hole pair and phonon excitations in the interaction of reactive gas molecules and atoms with metal surfaces. We present a theoretical framework that allows us to evaluate within a full-dimensional dynamics the combined contribution of both excitation mechanisms while the gas particle-surface interaction is described by an ab initio potential energy surface. The model is applied to study energy dissipation in the scattering of N(2) on W(110) and N on Ag(111).
View Article and Find Full Text PDFFirst-principles calculations are performed to obtain the dielectric function and loss spectra of bulk PdH(x). Hydrogen concentrations between x = 0 and 1 are considered. The calculated spectra are dominated by a broad peak that redshifts in energy with x.
View Article and Find Full Text PDFThe non-reactive scattering of N(2) from the W(110) surface is studied with six dimensional (6D) classical dynamics and two distinct potential energy surfaces (PES). Here, we use the PESs calculated with density functional theory and two different exchange-correlation functionals, the PW91 [J. E.
View Article and Find Full Text PDFThe role of electron-hole pair excitations in the dynamics of N(2) on W(100) and W(110) is evaluated using a theoretical model that accounts for the six-dimensionality of the problem in the whole calculation. The six-dimensional potential energy surface is determined in each case from an extensive grid of energies calculated with density functional theory. Dissipative effects due to electron-hole pair excitations are introduced in the classical dynamics equations through a friction force.
View Article and Find Full Text PDFWe study the dissociative dynamics of O(2) molecules on the Ag(100) surface. Initially, the impinging molecules are either in the spin-triplet ground state or in the spin-singlet excited state. The molecule-surface interaction is obtained in each case by constructing the six-dimensional potential energy surface (PES) from the interpolation of the energies calculated with spin-polarized and non-spin-polarized density functional theories, respectively.
View Article and Find Full Text PDFWe study the dissociative adsorption of N(2) on W(100) and W(110) by means of density functional theory and classical dynamics. Working with a full six-dimensional adiabatic potential energy surface (PES), we find that the theoretical results of the dynamical problem strongly depend on the choice of approximate exchange-correlation functional for the determination of the PES. We consider the Perdew-Wang-91 [Perdew et al.
View Article and Find Full Text PDFThe wave nature of particles is rarely observed, in part because of their very short de Broglie wavelengths in most situations. However, even with wavelengths close to the size of their surroundings, the particles couple to their environment (for example, by gravity, Coulomb interaction, or thermal radiation). These couplings shift the wave phases, often in an uncontrolled way, and the resulting decoherence, or loss of phase integrity, is thought to be a main cause of the transition from quantum to classical behavior.
View Article and Find Full Text PDFThe six-dimensional potential energy surface for the dissociation of N2 molecules on the W(110) surface has been determined by density functional calculations and interpolated using the corrugation reducing procedure. Examination of the resulting six-dimensional potential energy surface shows that nonactivated paths are available for dissociation. In spite of this, the dissociation probability goes to a very small value when the impact energy goes to zero and increases with increasing energy, a behavior usually associated with activated systems.
View Article and Find Full Text PDFLow-energy N2 molecules easily dissociate on W(100) but not on W(110). In this Letter, the six-dimensional potential energy surface for the dissociation of N2 molecules on W(110) has been determined by density functional calculations. Results are compared to those of N2 dissociation on W(100).
View Article and Find Full Text PDF