How Adsorbed Oxygen Atoms Inhibit Hydrogen Dissociation on Tungsten Surfaces.

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 H2 dissociation process. The reaction probability is calculated for kinetic energies below 300 meV and different O nominal coverages. We show that the adsorbed O atoms act as repulsive centers that modulate the dynamics of the impinging H2 molecules by closing dissociation pathways. In agreement with existing experimental information, H2 dissociation is absent for an O coverage of half a monolayer. The results show that the influence of O adsorbates on the dissociation dynamics on W(110) goes much beyond the blocking of possible H adsorption sites. Adsorbed O atoms create a sort of chemical shield at the surface that prevents further approach and dissociation of the H2 molecules.

Statistical properties of quantum probability fluctuations in complex-forming chemical reactions.

The energy dependence of quantum complex-forming reaction probabilities is well known to involve sharp fluctuations, but little seems to be known about their amplitudes. We develop here, for triatomic reactions, an analytical approach of their statistical distribution. This approach shows that the fluctuation amplitudes depend essentially on the number of available quantum states in the reagent and product channels. Moreover, the more numerous the product states, the more efficiently the fluctuations of their populations compensate each other when they add up to give the reaction probability. The predictions of our approach appear to be in good quantitative agreement with quantum scattering calculations for the prototypical reaction H+ + H2.

When Classical Trajectories Get to Quantum Accuracy: The Scattering of H2 on Pd(111).

When 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 H2 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.

Theoretical Study of Barrierless Chemical Reactions Involving Nearly Elastic Rebound: The Case of S(1D) + X2, X = H, D.

For some values of the total angular momentum consistent with reaction, the title processes involve nonreactive trajectories proceeding through a single rebound mechanism during which the internal motion of the reagent diatom is nearly unperturbed. When such paths are in a significant amount, the classical reaction probability is found to be markedly lower than the quantum mechanical one. This finding was recently attributed to an unusual quantum effect called diffraction-mediated trapping, and a semiclassical correction was proposed in order to take into account this effect in the classical trajectory method. In the present work, we apply the resulting approach to the calculation of opacity functions as well as total and state-resolved integral cross sections (ICSs) and compare the values obtained with exact quantum ones, most of which are new. As the title reactions proceed through a deep insertion well, mean potential statistical calculations are also presented. Seven values of the collision energy, ranging from 30 to 1127 K, are considered. Two remarkable facts stand out: (i) The corrected classical treatment strongly improves the accuracy of the opacity function as compared to the usual classical treatment. When the entrance transition state is tight, however, those trajectories crossing it with a bending vibrational energy below the zero point energy must be discarded. (ii) The quantum opacity function, particularly its cutoff, is finely reproduced by the statistical approach. Consequently, the total ICS is also very well described by the two previous approximate methods. These, however, do not predict state-resolved ICSs with the same accuracy, proving thereby that (i) one or several genuine quantum effects involved in the dynamics are missed by the corrected classical treatment and (ii) the dynamics are not fully statistical.

Dynamics of H2 Eley-Rideal abstraction from W(110): sensitivity to the representation of the molecule-surface potential.

Dynamics of the Eley-Rideal (ER) abstraction of H2 from W(110) is analyzed by means of quasi-classical trajectory calculations. Simulations are based on two different molecule-surface potential energy surfaces (PES) constructed from Density Functional Theory results. One PES is obtained by fitting, using a Flexible Periodic London-Eyring-Polanyi-Sato (FPLEPS) functional form, and the other by interpolation through the corrugation reducing procedure (CRP). Then, the present study allows us to elucidate the ER dynamics sensitivity on the PES representation. Despite some sizable discrepancies between both H+H/W(110) PESs, the obtained projectile-energy dependence of the total ER cross sections are qualitatively very similar ensuring that the main physical ingredients are captured in both PES models. The obtained distributions of the final energy among the different molecular degrees of freedom barely depend on the PES model, being most likely determined by the reaction exothermicity. Therefore, a reasonably good agreement with the measured final vibrational state distribution is observed in spite of the pressure and material gaps between theoretical and experimental conditions.

Surface temperature effects on the dynamics of N2 Eley-Rideal recombination on W(100).

Quasiclassical trajectories simulations are performed to study the influence of surface temperature on the dynamics of a N atom colliding a N-preadsorbed W(100) surface under normal incidence. A generalized Langevin surface oscillator scheme is used to allow energy transfer between the nitrogen atoms and the surface. The influence of the surface temperature on the N(2) formed molecules via Eley-Rideal recombination is analyzed at T = 300, 800, and 1500 K. Ro-vibrational distributions of the N(2) molecules are only slightly affected by the presence of the thermal bath whereas kinetic energy is rather strongly decreased when going from a static surface model to a moving surface one. In terms of reactivity, the moving surface model leads to an increase of atomic trapping cross section yielding to an increase of the so-called hot atoms population and a decrease of the direct Eley-Rideal cross section. The energy exchange between the surface and the nitrogen atoms is semi-quantitatively interpreted by a simple binary collision model.

High-harmonic transient grating spectroscopy of NO2 electronic relaxation.

We study theoretically and experimentally the electronic relaxation of NO(2) molecules excited by absorption of one ∼400 nm pump photon. Semiclassical simulations based on trajectory surface hopping calculations are performed. They predict fast oscillations of the electronic character around the intersection of the ground and first excited diabatic states. An experiment based on high-order harmonic transient grating spectroscopy reveals dynamics occurring on the same time scale. A systematic study of the detected transient is conducted to investigate the possible influence of the pump intensity, pump wavelength, and rotational temperature of the molecules. The quantitative agreement between measured and predicted dynamics shows that, in NO(2), high harmonic transient grating spectroscopy encodes vibrational dynamics underlying the electronic relaxation.

Dynamical reaction pathways in Eley-Rideal recombination of nitrogen from W(100).

The scattering of atomic nitrogen over a N-pre-adsorbed W(100) surface is theoretically described in the case of normal incidence off a single adsorbate. Dynamical reaction mechanisms, in particular Eley-Rideal (ER) abstraction, are scrutinized in the 0.1-3.0 eV collision energy range and the influence of temperature on reactivity is considered between 300 and 1500 K. Dynamics simulations suggest that, though non-activated reaction pathways exist, the abstraction process exhibits a significant collision energy threshold (0.5 eV). Such a feature, which has not been reported so far in the literature, is the consequence of a repulsive interaction between the impinging and the pre-adsorbed nitrogens along with a strong attraction towards the tungsten atoms. Above threshold, the cross section for ER reaction is found one order of magnitude lower than the one for hot-atoms formation. The abstraction process involves the collision of the impinging atom with the surface prior to reaction but temperature effects, when modeled via a generalized Langevin oscillator model, do not affect significantly reactivity.

The O((1)D) + H2 (X (1)Σ+, v, j) → OH(X (2)Π, v', j') + H((2)S) reaction at low collision energy: when a simple statistical description of the dynamics works.

In this communication, we highlight that statistical approaches for chemical reactions describe reasonably well the low energy dynamics of the title process. Consequently, such methods prove to be valuable to compute rate constants from low to room temperatures. Results are compared with experiment and recent precise quantum wave packet calculations [J. Phys. Chem. A, 2009, 113, 5285].

Dynamics simulation of N(2) scattering onto W(100,110) surfaces: A stringent test for the recently developed flexible periodic London-Eyring-Polanyi-Sato potential energy surface.

An efficient method to construct the six dimensional global potential energy surface (PES) for two atoms interacting with a periodic rigid surface, the flexible periodic London-Eyring-Polanyi-Sato model, has been proposed recently. The main advantages of this model, compared to state-of-the-art interpolated ab initio PESs developed in the past, reside in its global nature along with the small number of electronic structure calculations required for its construction. In this work, we investigate to which extent this global representation is able to reproduce the fine details of the scattering dynamics of N(2) onto W(100,110) surfaces reported in previous dynamics simulations based on locally interpolated PESs. The N(2)/W(100) and N(2)/W(110) systems are chosen as benchmarks as they exhibit very unusual and distinct dissociative adsorption dynamics although chemically similar. The reaction pathways as well as the role of dynamic trapping are scrutinized. Besides, elastic/inelastic scattering dynamics including internal state and angular distributions of reflected molecules are also investigated. The results are shown to be in fair agreement with previous theoretical predictions.

Transformation from angle-action variables to Cartesian coordinates for polyatomic reactions.

The transformation from angle-action variables to Cartesian coordinates is an important step of the semiclassical description of bimolecular collisions and photofragmentations. The basic reason is that dynamical conditions corresponding to molecular beam experiments are ideally generated in angle-action variables, whereas the classical equations of motion are ideally solved in Cartesian coordinates by standard numerical approaches. To our knowledge, this transformation is available in the literature only for atom-diatom arrangements. The goal of the present work is to derive it for diatom-polyatom ones. The analogous transformation for any type of arrangement may then be straightforwardly deduced from that presented here.

Absorption and emission lineshapes and ultrafast solvation dynamics of NO in parahydrogen.

We perform a theoretical study on the electronic spectroscopy of dilute NO impurity embedded in parahydrogen (p-H(2)). Absorption and emission lineshapes for the A (2)Sigma(+)<--X (2)Pi Rydberg transition of NO in parahydrogen have been previously measured and simulated, which yielded results for the NO/p-H(2) ground and excited state pair potentials [L. Bonacina et al., J. Chem. Phys. 125, 054507 (2006)]. Using these potentials, we performed molecular dynamics simulation, theoretical statistical mechanical calculations of absorption and emission lineshapes, and both equilibrium and nonequilibrium solvation correlation functions for NO chromophore in parahydrogen. Theory was shown to be in good agreement with simulation. Linear response treatment of solvation dynamics was shown to break down due to a dramatic change in the solute-solvent microstructure upon solute excitation to the Rydberg state and the concomitant increase of the solute size.

Mean potential phase space theory of chemical reactions.

A nonconventional application of phase space theory to the insertion reactions A+H(2), with A=C((1)D) and S((1)D), is presented. Instead of approximating the potential energies of interaction between separated fragments by their isotropic long-range contributions, as in the original theory, the latter are replaced by the accurate potential energies averaged with respect to Jacobi angles. The integral and differential cross sections obtained from this mean potential phase space theory (MPPST) turn out to be in very satisfying agreement with the benchmark predictions of the time-independent and time-dependent statistical quantum methods. The formal and numerical simplicity of MPPST with respect to any approach combining statistical assumptions and dynamical calculations makes it a promising tool for studying indirect polyatomic reactions.

On the theory of complex-forming chemical reactions: effect of parity conservation on the polarization of differential cross sections.

For complex-forming triatomic reactions such as the prototypical insertion reactions intensively studied in the last few years, quantum mechanical differential cross sections (DCS) present sharp forward/backward polarization peaks when the reagent rotational angular momentum quantum number j is zero. Moreover, the size of the peaks decreases rapidly with increasing j values so that for j = 3, they are no longer visible. In contrast, the polarization peaks are always missing in the classical mechanical DCSs. Apart from the peaks, however, the quantum and classical DCSs are usually in good agreement. In a recent rapid communication, we showed that the fundamental reason for the previous differences in the quantum and classical scenarios is that parity conservation leads in quantum mechanics to an angular momentum constraint without equivalent in classical mechanics. We also proposed a parity-restoring approximation leading to an accurate semi-classical description of the peaks. While only the main lines of the demonstration were given in the communication, we report here the whole developments. We also analyse why the peaks disappear when the reagent diatom is rotationally excited. As a by-product of the previous developments, we finally discuss the possibility of a general statistico-dynamical semiclassical approach.

Classical treatment of molecular collisions: striking improvement of the description of recoil energy distributions using Gaussian weighted trajectories.

The Gaussian weighting (GW) procedure, recently used in the classical treatment of molecular collisions, is a practical way of taking into account quantization of product vibrational actions. The goal of this brief communication is to show that the GW procedure may drastically improve the predictions of the recoil energy distribution between final fragments, an observable frequently measured in molecular beam experiments.

Parity conservation and polarization of differential cross sections in complex-forming chemical reactions.

For complex-forming chemical reactions, such as atom-diatom insertion reactions, quantum scattering and quantum statistical calculations usually predict sharp forward/backward peaks in the Differential Cross Sections (DCS). Conversely, the corresponding classical calculations are unable to reproduce these peaks. We show here that the basic reason for such an intriguing failure is that parity conservation is ignored in classical mechanics. A by-product of the analysis is a simple parity-restoring approximation that might significantly increase the ability of classical mechanics to describe DCSs over the whole angular range for the title processes.

The ultrafast structural response of solid parahydrogen: a complementary experimental/simulation investigation.

We present a complete characterization, based on femtosecond pump-probe spectroscopy and molecular dynamics simulations, of the ultrafast dynamics of electronic bubble formation in solid parahydrogen upon impulsive excitation of impurity-doped sites, which correlate with the lowest Rydberg state of the NO impurity. The high temporal resolution of the experiment allows us to identify three time scales in the structural dynamics. A first ultrafast expansion (<150 fs), associated with the release of approximately 80% of the excess energy available to the system after excitation, is accompanied by a transient narrowing of the spatial distribution of the first shell of H2 molecules around the impurity. In a subsequent stage (up to approximately 800 fs), the cavity expansion slows down, and energy starts to flow irreversibly into the crystal. Finally, the lattice undergoes a slow structural reorganization at the impurity site (5-10 ps). A weak low-frequency recurrence, probably associated with an elastic response of the crystal, is observed at approximately 10 ps. The absence of polarization dependence indicates that the dynamics is largely dominated by translational (radial) motions of the molecules surrounding NO and not by the rotational motion of the impurity. Molecular dynamics simulations with temperature corrections, to mimic zero-point fluctuations, fully support the experimental results and show that the bubble model is suited to describe the dynamics of the system. It appears that the response of the medium around the impurity at short times is typical of a liquid solvent rather than that of a solid.

Validity of phase space theory for atom-diatom insertion reactions.

Phase space theory (PST) is applied to the calculation of state-resolved integral and differential cross sections for the complex-forming atom-diatom insertion reactions A + H(2) --> AH(2) --> AH + H with A = C((1)D), S((1)D), O((1)D), and N((2)D). In the asymptotic channels, vibration motion is quantized while rotation and translation motions are treated classically. The approach is compared to exact quantum scattering calculations and quantum statistical models. Given the simplicity of PST, the agreement with the previous much more refined treatments is very satisfying. Although PST is a well-established theory, this work is, to our knowledge, the first such systematic comparison of its predictions with accurate quantum scattering and quantum statistical calculations.

Lattice response of quantum solids to an impulsive local perturbation.

The lattice response of solid para-H2 to an impulsive electronic excitation was studied using femtosecond pump-probe spectroscopy. The evolution of an electronic bubble in the crystal, created upon excitation of the A(3ssigma) Rydberg state of an NO impurity, was followed in real time, with a resolution of 100 fs. The experimental results, interpreted in connection with molecular dynamics simulations with quantum corrections, indicate the presence of three stages in the dynamics: a sub-100 fs "adiabatic" phase, a 0.5-1 ps phase, corresponding to the interaction of the first with the next shells driven by the bubble expansion, and a 5 ps phase, corresponding to a slow rearrangement of the environment surrounding the impurity. These findings indicate that the lattice response in solid para-H2 resembles that of a liquid.