Dynamics of charge transfer states on metal surfaces: the competition between reactivity and quenching

The dynamics of excited states of adsorbates on surfaces caused by charge transfer is studied. Both negative and positive charge transfer processes are possible. In particular we are interested in positive charge transfer from a metal surface to molecular or atomic oxygen adsorbed on the surface. Once the negatively charged oxygen on the surface loses an electron it becomes chemically activated. The ability of this species to react depends on the quenching time or back transfer. The analysis of these processes is based on a set of diabatic potential energy surfaces each representing a different charged oxygen species. The dynamics is followed by solving the multichannel time-dependent Schr?dinger equation or Liouville von Neumann equation. Due to the nonadiabatic character of these reactions large isotope effects are predicted.     

Oxygen dissociation on palladium and gold core/shell nanoparticles

Oxygen dissociation reaction on gold, palladium, and gold‐palladium core/shell nanoparticles was investigated with plane wave basis set, density functional theory. Bader population analysis of charge and electron distribution was employed to understand the change of catalytic activity as a function of the nanopaticle composition. The nanoparticles’ electronic properties were investigated and the degree of core/shell charge polarization was estimated for each composition. It was found that surface polarization plays an important role in the catalysis of the initial step of electrophile reactions such as oxygen dissociation. We have investigated the O₂ adsorption energy on each nanoparticle and the activation barrier for the oxygen dissociation reaction as a function of the nanoparticle structure. Furthermore, we have investigated the influence of surface geometry, that is., surface bond lengths on the catalytic activity. We have compared the electronic and the geometry effects on the oxygen activation and dissociation. Our design rules for core/shell nanoparticles offer an effective method for control of the surface catalytic activity. Palladium and gold are often used as catalysts in synthetic chemistry. First‐principles calculations elucidate the mechanisms that control the surface reactivity of gold, palladium, and gold‐palladium core shell nanoparticles in oxygen dissociation reactions. Oxygen dissociation is promoted on the gold surface of gold/palladium core‐shell nanoparticles by favorable electron transfer from the core to the shell. Such core‐shell electronic effects can be used for fine‐tuning the nanoparticles catalytic activity.     

Quantum wave packet study of the H+ + D2 reaction on diabatic potential energy surfaces

The exact three-dimensional nonadiabatic quantum dynamics calculations were carried out for the title reaction by a time-dependent wave packet approach based on a newly constructed diabatic potential energy surface (Kamisaka et al. J. Chem. Phys. 2002, 116, 654). Three processes including those of reactive charge transfer, nonreactive charge transfer, and reactive noncharge transfer were investigated to determine the initial state-resolved probabilities and reactive cross sections. The results show that a large number of resonances can be observed in the calculated probabilities due to the deep well on adiabatic ground surface and the dominant process is the reactive noncharge-transfer process. Some interesting dynamical features such as v-dependent and j-dependent behaviors of the probabilities are also revealed. In addition, a good agreement has been achieved in the comparison between the calculated quantum cross sections from the ground rovibrational initial state and the experimental measurement data.     

Attachment cross sections of protonated water clusters

The attachment of water molecules onto size selected protonated water clusters has been experimentally investigated. Absolute attachment cross sections are measured as a function of cluster size, collision energy, and initial cluster temperature. Although thermal evaporation is ruled out in our experiment, attachment cross sections become significantly smaller than hard sphere cross sections as the collision energy increases. This feature is attributed to a transition from adiabatic to nonadiabatic regime. It is shown to be due to a dynamical effect: as the collision duration becomes shorter than the typical time required for collision energy redistribution into clusters internal energy, the attachment probability is reduced. We relate this typical time to the period of the main surface vibrational mode excited by the collisions. This hypothesis is further supported by results obtained with deuterated water clusters.     

Influence of nuclear exchange on nonadiabatic electron processes in H(+)+H2 collisions

H(+)+H(2) collisions are studied by means of a semiclassical approach that explicitly accounts for nuclear rearrangement channels in nonadiabatic electron processes. A set of classical trajectories is used to describe the nuclear motion, while the electronic degrees of freedom are treated quantum mechanically in terms of a three-state expansion of the collision wavefunction. We describe electron capture and vibrational excitation, which can also involve nuclear exchange and dissociation, in the E = 2-1000 eV impact energy range. We compare dynamical results obtained with two parametrizations of the potential energy surface of H(3)(+) ground electronic state. Total cross sections for E > 10 eV agree with previous results using a vibronic close-coupling expansion, and with experimental data for E < 10 eV. Additionally, some prototypical features of both nuclear and electron dynamics at low E are discussed.     

Molecular dynamics of hydrogen dissociation on an oxygen covered Pt(111) surface

The dissociation of hydrogen on a Pt(111) surface covered with a p(2 x 2) oxygen phase was investigated using quasiclassical, six-dimensional molecular dynamics. The potential energy surface (PES) used in these simulations was obtained by an iterative novelty sampling algorithm. Compared to molecular beam experiments performed under similar conditions, the simulations give an accurate prediction of the reaction probability via a direct dissociation pathway. When compared to previously reported reaction probability curves for the clean Pt(111) surface, we find that the presence of an oxygen overlayer inhibits the direct pathway to hydrogen dissociation. This inhibition is a function of incident energy and cannot be described by a simple site blocking model. An indirect pathway to dissociation, which was observed in experiments, is not properly captured by the PES. Spatially resolved "reaction maps" indicate that the preferred site for hydrogen dissociation on an oxygen covered Pt surface is the top site of the Pt atom farthest from the adsorbed oxygen atom.     

Stability of atomic oxygen chemisorption on Pt‐alloy surfaces

A density functional theory calculation is used to investigate the atomic oxygen (O) stability over platinum (Pt) and Pt‐based alloy surfaces. Here, the stability is connected with the preferential adsorption sites for O chemisorptions and the adsorption energy. Thus, the interaction mechanism between atomic O and metal surfaces is studied by using charge transfer analysis. In this present paper, atomic structure and binding energy of oxygen adsorption on the Pt(111) are in a very good agreement with experiment and previous density functional theory calculations. Furthermore, we obtained that the addition of ruthenium (Ru) and molybdenum (Mo) on the pure Pt surface enhances the adsorption energy. Our charge transfer analysis shows that the largest charge transfer contributing to the metal‐O bonding formation is observed in the case of O/PtRuMo surface followed by O/PtRu surface. This is in consistency with metal d‐orbital characteristic, where Mo has much more empty d‐orbital than Ru in correspondence to accept electrons from atomic oxygen. Copyright © 2016 John Wiley & Sons, Ltd.     

The surface temperature dependence of the inelastic scattering and dissociation of hydrogen molecules from metal surfaces

High-dimensional, wave packet calculations have been carried out to model the surface temperature dependence of rovibrationally inelastic scattering and dissociation of hydrogen molecules from the Cu(111) surface. Both the molecule and the vibrating surface are treated fully quantum-mechanically. It is found, in agreement with experimental data, that the surface temperature dependence of a variety of dynamical processes has an Arrhenius form with an activation energy dependent on molecular translational energy and on the initial and final molecular states. The activation energy increases linearly with decreasing translational energy below the threshold energy. Above threshold the behavior is more complex. A quasianalytical model is proposed that faithfully reproduces the Arrhenius law and the translational energy dependence of the activation energy. In this model, it is essential to include quantized energy transfer between the surface and the molecule. It further predicts that for any process characterized by a large energy barrier and multiphonon excitation, the linear change in activation energy up to threshold has slope-1. This explains successfully the universal nature of the unit slope found experimentally for H2 and D2 dissociation on Cu.     

A novel conical intersection topography and its consequences: the 1, 2 2A conical intersection seam of the vinoxy radical

A region of the 1, 2 2A seam of accidental conical intersections in the vinoxy radical exhibits a novel topography which has important consequences for both upper-state to lower-state and lower-state to upper-state nonadiabatic transitions. The impact of this topography on these nonadiabatic transitions is described. We also considered the possibility that this conical intersection seam provides a dynamical bottleneck to the photodissociation of vinoxy to H+ketene by facilitating nonadiabatic recrossing. Our analysis of the conical topographies and the proximity of the conical intersections to the transition state for dissociation to H+ketene does not support nonadiabatic recrossing as an effective dynamical bottleneck blocking the H+ketene channel.     

Product angular distributions in the ultraviolet photodissociation of N2O

The angular distribution of products from the ultraviolet photodissociation of nitrous oxide yielding O((1)D) and N(2)(X Σ(g)(+)(1)) was investigated using classical trajectory calculations. The calculations modeled absorption only to the 2(1)A(') electronic state but used surface-hopping techniques to model nonadiabatic transitions to the ground electronic state late in the dissociation. Observed values of the anisotropy parameter β, which decrease as the product N(2) rotational quantum number j increases, could be well reproduced. The relatively low observed β values arise principally from nonaxial recoil due to the very strong bending forces present in the excited state. In the main part of the product rotational distribution near 203 nm, an unusual dynamical effect produces the decrease in β with increasing j; nonaxial recoil effects remain approximately constant while higher j product molecules arise from parent molecules that had their transition dipole moments aligned more closely along the molecular axis. In both low and high j tails of the rotational distribution, the variations in β with j are caused by changes in the extent of nonaxial recoil. In the high-j tail, additional torque present on the ground state potential energy surface following nonadiabatic transitions causes both the additional rotational excitation and the lower β values.     

Electronic nonadiabatic effects in the adsorption of hydrogen atoms on metals

The time-dependent, mean-field Newns-Anderson model for a spin-polarized adsorbate approaching a metallic surface is solved in the wide-band limit. Equations for the time evolution of the occupation of the spin dependent adsorbate states and for the nonadiabatic and nearly adiabatic adsorbate-surface energy transfer rates are derived. Numerical solutions are obtained using characteristic parameters derived from density functional theory calculations for the H/Cu(111) system. The time evolution of the model system is shown to be strongly nonadiabatic in the vicinity of the transition point between spin-polarized and nonpolarized ground states. Away from the spin transition the nonadiabatic energy transfer is in close agreement with the nearly adiabatic limit. Near the transition, nonadiabatic effects are large and the nearly adiabatic approximation fails.     

甲烷在O/Ni(100)表面上的反应动力学研究

甲烷在Ｏ／Ｎｉ（１００）表面上的反应动力学研究俞华根，程极源（中国科学院成都有机化学研究所，成都６１００１５）关键词甲烷，活化解离，预吸附氧，Ｎｉ（１００）表面，分子动力学，势能面甲烷在金属催化剂表面活化解离是重要的催化反应，受到了高度重视．近年来，...     

Application of mean-field and surface-hopping approaches for interrogation of the Xe3+ molecular ion photoexcitation dynamics

The dissociation dynamics of the excited Xe(3) (+) molecular ion through the Pi(12)(u) and Pi(12)(g) conical intersection was interrogated by computational simulation in which no adjustable parameters were used. The electronic ground and excited state potential energy surfaces were generated by the diatomics-in-molecules method, and the Ehrenfest mean-field and Tully surface-hopping approaches treated the nonadiabatic interactions. Reproduction of the experimental spectrum of the symmetric photofragmentation as a function of excitation energy was obtained within the region of interest (2.5-3.75 eV), with the exception of a 0.25 eV width on the red side of the spectral apex. Good agreement was obtained with the experimental dissociated photofragment kinetic energy spectra. It was determined that the greatest contribution to the nonadiabatic coupling between the two states originated from the bending vibrational mode of the molecule in the Sigma(12)(u), ground electronic state before excitation.     

Collision induced dissociation of doubly-charged ions: Coulomb explosion vs. neutral loss in [Ca(urea)](2+) gas phase unimolecular reactivity via chemical dynamics simulations

In this paper we report different theoretical approaches to study the gas-phase unimolecular dissociation of the doubly-charged cation [Ca(urea)](2+), in order to rationalize recent experimental findings. Quantum mechanical plus molecular mechanical (QM/MM) direct chemical dynamics simulations were used to investigate collision induced dissociation (CID) and rotational-vibrational energy transfer for Ar + [Ca(urea)](2+) collisions. For the picosecond time-domain of the simulations, both neutral loss and Coulomb explosion reactions were found and the differences in their mechanisms elucidated. The loss of neutral urea subsequent to collision with Ar occurs via a shattering mechanism, while the formation of two singly-charged cations follows statistical (or almost statistical) dynamics. Vibrational-rotational energy transfer efficiencies obtained for trajectories that do not dissociate during the trajectory integration were used in conjunction with RRKM rate constants to approximate dissociation pathways assuming complete intramolecular vibrational energy redistribution (IVR) and statistical dynamics. This statistical limit predicts, as expected, that at long time the most stable species on the potential energy surface (PES) dominate. These results, coupled with experimental CID from which both neutral loss and Coulomb explosion products were obtained, show that the gas phase dissociation of this ion occurs by multiple mechanisms leading to different products and that reactivity on the complicated PES is dynamically complex.     

Nonadiabatic time-dependent wave packet study of the D+ + H2 reaction system

A theoretical investigation on the nonadiabatic processes of the D(+) + H(2) reaction system has been carried out by means of exact three-dimensional nonadiabatic time-dependent wave packet calculations with an extended split operator scheme (XSOS). The diabatic potential energy surface newly constructed by Kamisaka et al. (J. Chem. Phys. 2002, 116, 654) was employed in the calculations. This study provided quantum cross sections for three competing channels of the reactive charge transfer, the nonreactive charge transfer, and the reactive noncharge transfer, which contrasted markedly to many previous quantum theoretical reports on the (DH(2))(+) system restricted to the total angular momentum J = 0. These quantum theoretical cross sections derived from the ground rovibrational state of H(2) show wiggling structures and an increasing trend for both the reactive charge transfer and the nonreactive charge transfer but a decreasing trend for the reactive noncharge transfer throughout the investigated collision energy range 1.7-2.5 eV. The results also show that the channel of the reactive noncharge transfer with the largest cross section is the dominant one. A further investigation of the v-dependent behavior of the probabilities for the three channels revealed an interesting dominant trend for the reactive charge transfer and the nonreactive charge transfer at vibrational excitation v = 4 of H(2). In addition, the comparison between the centrifugal sudden (CS) and exact calculations showed the importance of the Coriolis coupling for the reactive system. The computed quantum cross sections are also compared with the experimental measurement results.     

Statistical mechanics of energy transfer in gas–surface scattering: A dynamical Lie algebraic approach

The dynamical Lie algebraic (DLA) method is used to describe statistical mechanics of energy transfer in rotationally inelastic molecule–surface scattering. Statistical average values of an observable for the scattering system are calculated in terms of density operator formalism in statistical mechanics. Employing a cubic expansion procedure of molecule–surface interaction potential leads to generation of a dynamical Lie algebra. Thus these statistical average values as a function of the group parameters can be obtained analytically in this formulation. The group parameters can be found from solving a set of coupled nonlinear differential equations. The DLA method, which has no need for determination of transition probabilities in advance as made routinely in the calculation, offers an efficient alternative to the method for computing the statistical average values. This method is much less computationally intensive because most of calculations can be analytically carried out. The average final rotational energies and their dependence on the main dynamic variables and the average interaction potential are presented for the rotationally inelastic scattering of NO molecules from a flat, static Ag(111) surface. Direct comparison is made between the predictions of this model calculation and experiment. The model reproduces well the degree of rotational excitation and correlation between the average final translational and the average rotational energies. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Accurate ab initio determination of the adiabatic potential energy function and the Born-Oppenheimer breakdown corrections for the electronic ground state of LiH isotopologues

High level ab initio potential energy functions have been constructed for LiH in order to predict vibrational levels up to dissociation. After careful tests of the parameters of the calculation, the final adiabatic potential energy function has been composed from: (a) an ab initio nonrelativistic potential obtained at the multireference configuration interaction with singles and doubles level including a size-extensivity correction and quintuple-sextuple ζ extrapolations of the basis, (b) a mass-velocity-Darwin relativistic correction, and (c) a diagonal Born-Oppenheimer (BO) correction. Finally, nonadiabatic effects have also been considered by including a nonadiabatic correction to the kinetic energy operator of the nuclei. This correction is calculated from nonadiabatic matrix elements between the ground and excited electronic states. The calculated vibrational levels have been compared with those obtained from the experimental data [J. A. Coxon and C. S. Dickinson, J. Chem. Phys. 134, 9378 (2004)]. It was found that the calculated BO potential results in vibrational levels which have root mean square (rms) deviations of about 6-7 cm(-1) for LiH and ～3 cm(-1) for LiD. With all the above mentioned corrections accounted for, the rms deviation falls down to ～1 cm(-1). These results represent a drastic improvement over previous theoretical predictions of vibrational levels for all isotopologues of LiH.     

Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation

We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schro?dinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used in the diabatic basis to measure the importance of the diabatic couplings. The method is tested on three model problems introduced by Tully and on a two-surface model of dissociation of NaI.     

Electron kinetic energies from vibrationally promoted surface exoemission: evidence for a vibrational autodetachment mechanism

We report kinetic energy distributions of exoelectrons produced by collisions of highly vibrationally excited NO molecules with a low work function Cs dosed Au(111) surface. These measurements show that energy dissipation pathways involving nonadiabatic conversion of vibrational energy to electronic energy can result in electronic excitation of more than 3 eV, consistent with the available vibrational energy. We measured the dependence of the electron energy distributions on the translational and vibrational energy of the incident NO and find a clear positive correlation between final electron kinetic energy and initial vibrational excitation and a weak but observable inverse dependence of electron kinetic energy on initial translational energy. These observations are consistent with a vibrational autodetachment mechanism, where an electron is transferred to NO near its outer vibrational turning point and ejected near its inner vibrational turning point. Within the context of this model, we estimate the NO-to-surface distance for electron transfer.