Measurements and simulations of high energy O(3P) + Ar(1S) angular scattering: single and multi-collision regimes

We present differential angular cross sections for O(3P) + Ar(1S) scattering at collision energies near 90 kcal mol(-1) (approximately 8 km s(-1) relative velocity) from molecular beam measurements and high-level theoretical calculations. Beams of hyperthermal O(3P) are now being used to investigate novel gas-phase and gas-surface chemistries, and the comparison of theory and measurements on this simple system will be a stringent test of the experimental methodology. Potential energy curves were generated for O(3P) + Ar(1S) using a large cc-pVQZ basis within a valence multi-configuration plus perturbation theory treatment. These curves were then used in quantum scattering calculations to generate differential cross sections. Agreement between experiment and theory is excellent. In addition to these comparisons, the cross sections were used in direct simulation Monte Carlo calculations to investigate effects of increasing the Ar flux above the "single-collision" regime. As the Ar flux increases, the observed differential angular cross sections change in two ways. In addition to the main "single-scatter" peak along the incident O-atom beam direction, a secondary O-atom peak appears in the direction of the incident Ar beam, and the multiple-scattered O-atom translational energy starts to reflect the energy of the relatively slow moving Ar beam.     

Scattering of Xe from graphite

Recently a series of experimental measurements for the scattering of Xe atoms from graphite has been reported for both energy-resolved spectra and angular distributions. This system is of fundamental interest because the projectile Xe atoms are considerably more massive than the carbon atoms making up the graphite surface. These measurements were initially analyzed using the hard cubes model and molecular dynamics simulations, and both treatments indicated that the scattering process was a single collision in which the incoming Xe atom interacted strongly with a large number of carbon atoms in the outermost graphite layer. In this work we analyze the data using a single scattering theory that has been shown to explain a number of other experiments on molecular beam scattering from surfaces. These calculations confirm that the scattering process is a single collision with an effective surface mass that is substantially larger than that of the basic graphite ring.     

Experiments and simulations of hyperthermal Xe interacting with an ordered 1-decanethiol/Au(111) monolayer: penetration followed by high-energy, directed ejection

A study of the interaction of hyperthermal Xe with a well-ordered standing-up phase of 1-decanethiol adsorbed on Au(111) is presented. Experimentally, double-differential measurements were made of the postcollision Xe kinetic energy as a function of incident and final angles. These experiments are compared to classical trajectory calculations. The results show the two expected channels: direct-inelastic scattering from the surface and accommodated Xe due to trapping-desorption. There is also evidence of a further interaction mechanism. This involves the penetration of the atom deep into the channels between the aligned chains of the monolayer. When the collision energy has been dissipated, the implanted Xe is expelled as the chains return to their equilibrium positions. The expelled Xe leaves the surface with an energy much higher than expected for trapping-desorption, and with an angular-intensity distribution peaked close to the direction of the 1-decanethiol chain orientation. For this reason, we call this new scattering mechanism directed ejection.     

Crossed beams and theoretical studies of the dynamics of hyperthermal collisions between Ar and ethane

Crossed molecular beams experiments and classical trajectory calculations have been used to study the dynamics of Ar+ethane collisions at hyperthermal collision energies. Experimental time-of-flight and angular distributions of ethane molecules that scatter into the backward hemisphere (with respect to their original direction in the center-of-mass frame) have been collected. Translational energy distributions, derived from the time-of-flight distributions, reveal that a substantial fraction of the collisions transfer abnormally large amounts of energy to internal excitation of ethane. The flux of the scattered ethane molecules increased only slightly from directly backward scattering to sideways scattering. Theoretical calculations show angular and translational energy distributions which are in reasonable agreement with the experimental results. These calculations have been used to examine the microscopic mechanism for large energy transfer collisions ("supercollisions"). Collinear ("head-on") or perpendicular ("side-on") approaches of Ar to the C-C axis of ethane do not promote energy transfer as much as bent approaches, and collisions in which the H atom is "sandwiched" in a bent Ar...H-C configuration lead to the largest energy transfer. The sensitivity of collisional energy transfer to the intramolecular potential energy of ethane has also been examined.     

Angular intensity distribution of a molecular oxygen beam scattered from a graphite surface

The scattering of the oxygen molecule from a graphite surface has been studied using a molecular beam scattering technique. The angular intensity distributions of scattered oxygen molecules were measured at incident energies from 291 to 614 meV with surface temperatures from 150 to 500 K. Every observed distribution has a single peak at a larger final angle than the specular angle of 45° which indicates that the normal component of the translation energy of the oxygen molecule is lost by the collision with the graphite surface. The amount of the energy loss by the collision has been roughly estimated as about 30-41% based on the assumption of the tangential momentum conservation during the collision. The distributions have also been analyzed with two theoretical models, the hard cubes model and the smooth surface model. These results indicate that the scattering is dominated by a single collision event of the particle with a flat surface having a large effective mass. The derived effective mass of the graphite surface for the incoming oxygen is 9-12 times heavier than that of a single carbon atom, suggesting a large cooperative motion of the carbon atoms in the topmost graphene layer.     

Hyperthermal reactions of O and O2 with a hydrocarbon surface: direct C-C bond breakage by O and H-atom abstraction by O2

A C-C bond-breaking reaction has been observed when a beam containing hyperthermal oxygen was directed at a continuously refreshed saturated hydrocarbon liquid (squalane) surface. The dynamics of this C-C bond-breaking reaction have been investigated by monitoring time-of-flight and angular distributions of the volatile product, OCH3 or H2CO. The primary product is believed to be the methoxy radical, OCH3, but if this radical is highly internally excited, then it may undergo secondary dissociation to form formaldehyde, H2CO. Either the primary or the secondary product may scatter directly into the gas phase before thermal equilibrium with the surface is reached, or they may become trapped on the surface and desorb in thermal equilibrium with the surface. Direct, single-collision scattering events that produce a C-C bond-breaking product are described with a kinematic picture that allows the determination of the effective surface mass encountered by an incident O atom, the atom-surface collision energy in the center-of-mass frame, and the fraction of the center-of-mass collision energy that goes into translation of the scattered gaseous product and the recoiling surface fragment. The dynamical behavior of the C-C bond-breaking reaction is compared with that of the H-atom abstraction reaction, which was the subject of an earlier study. Another reaction, H-atom abstraction by O2 (which is present in the hyperthermal beam), has also been observed, and the dynamics of this reaction are compared with the inelastic scattering dynamics of O2 and the dynamics of H-atom abstraction by O. The dynamics involving direct inelastic and reactive scattering of O2 are also described in terms of a kinematic picture where the incident O2 molecule is viewed as interacting with a local region of the surface.     

Theoretical studies of hyperthermal O(3P) collisions with hydrocarbon self-assembled monolayers

We present a dynamics study of inelastic and reactive scattering processes in collisions of hyperthermal (5 eV) O(3P) atoms with a hydrocarbon self-assembled monolayer (SAM). Molecular-dynamics simulations are carried out using a quantum mechanics/molecular mechanics (QM/MM) interaction potential that uses a high quality semiempirical Hamiltonian for the QM part and the MM3 force field for the MM part. A variety of products coming from reaction are identified, including H abstraction to generate OH, O atom addition to the SAM with subsequent elimination of H atoms, and direct C-C breakage. The C-C breakage mechanism provides a pathway for significant surface mass loss in single reactive events whereas the O addition-H elimination channel leads to surface oxidation. Reaction probabilities, product energy, and angular distributions are examined to gain insight on polymer erosion in low Earth orbit conditions and on fundamentals of inelastic and reactive hyperthermal gas-surface interactions.     

Dynamics of energy transfer in collisions of O(3P) atoms with a 1-decanethiol self-assembled monolayer surface

Chemical dynamics simulations are reported of energy transfer in collisions of O(3P) atoms with a 300 K 1-decanethiol self-assembled monolayer (H-SAM) surface. The simulations are performed with a nonreactive potential energy surface, developed from PMP2/aug-cc-pVTZ calculations of the O(3P) + H-SAM intermolecular potential, and the simulation results represent the energy transfer dynamics in the absence of O(3P) reaction. Collisions energies E(i) of 0.12, 2.30, 11.2, 75.0, and 120.5 kcal/mol and incident angles theta(i) of 15, 30, 45, 60, and 75 degrees were considered in the study (theta(i) = 0 degrees is the surface normal). The translational energy distribution of the scattered O(3P) atoms, P(E(f)), may be deconvoluted into Boltzmann and non-Boltzmann components, with the former fraction identified as f(B). The trajectories are also analyzed in terms of three types; that is, direct scattering from and physisorption on the top of the H-SAM and penetration of the H-SAM. There are three energy regimes in the scattering dynamics. For the low E(i) values of 0.12 and 2.30 kcal/mol, physisorption is important and both f(B) and the average final translational energy of the scattered O(3P) atom, E(f), are nearly independent of the incident angle. The dynamics is much different for hyperthermal energies of 75.0 and 120.5 kcal/mol, where penetration of the surface is important. For hyperthermal collisions, the penetration probability decreases as theta(i) is increased, with a significant transition between theta(i) of 60 and 75 degrees . Hyperthermal penetration occurs upon initial surface impact and is more probable if the impinging O(3P) atom may move down a channel between the chains. For E(i) = 120.5 kcal/mol, 90% of the trajectories penetrate at theta(i) = 15 degrees , while only 3% penetrate at theta(i) = 75 degrees. For the former theta(i), the energy transfer to the surface is efficient with E(f) = 4.04 kcal/mol, but for the latter theta(i), E(f) = 85.3 kcal/mol! Particularly interesting penetrating trajectories are those in which O(3P) is trapped in the H-SAM for times exceeding 60 ps, linger near the Au substrate, and strike the Au substrate and scatter directly. For E(i) = 11.2 kcal/mol, there is a transition between the scattering dynamics for the low and hyperthermal collision energies. Additional detail in the energy transfer dynamics is obtained from the final polar and azimuthal angles, the residence time on/in the H-SAM, the minimum height with respect to the Au substrate, and the number of inner turning points in the O-atom's velocity. Calculated values of E(f) vs the final polar angle, theta(f), are in qualitative agreement with experiment. The O(3P) + H-SAM nonreactive energy transfer dynamics, for E(i) of 11.2 kcal/mol and lower, are very similar to previously reported Ne + H-SAM simulations.     

Scattering of hyperthermal argon atoms from clean and D-covered Ru(0001) surfaces

Hyperthermal Ar atoms were scattered from a Ru(0001) surface held at temperatures of 180, 400 and 600 K, and from a Ru(0001)-(1×1)D surface held at 114 and 180 K. The resultant angular intensity and energy distributions are complex. The in-plane angular distributions have narrow (FWHM ≤ 10°) near-specular peaks and additional off-specular features. The energy distributions show an oscillatory behavior as a function of outgoing angle. In comparison, scattered Ar atoms from a Ag(111) surface exhibit a broad angular intensity distribution and an energy distribution that qualitatively tracks the binary collision model. The features observed for Ru, which are most evident when scattering from the clean surface at 180 K and from the Ru(0001)-(1×1)D surface, are consistent with rainbow scattering. The measured TOF profiles cannot be adequately described with a single shifted Maxwell-Boltzmann distribution. They can be fitted by two components that exhibit complex variations as a function of outgoing angle. This suggests at least two significantly different site and∕or trajectory dependent energy loss processes at the surface. The results are interpreted in terms of the stiffness of the surface and highlight the anomalous nature of the apparently simple hcp(0001) ruthenium surface.     

Calculation of inelastic helium atom scattering from H2/NaCl(001)

The one-phonon inelastic low energy helium atom scattering theory is adapted to cases where the target monolayer is a p(1 × 1) commensurate square lattice. Experimental data for para-H(2)/NaCl(001) are re-analyzed and the relative intensities of energy loss peaks in the range 6 to 9 meV are determined. The case of the H(2)/NaCl(001) monolayer for 26 meV scattering energy is computationally challenging and difficult because it has a much more corrugated surface than those in the previous applications for triangular lattices. This requires a large number of coupled channels for convergence in the wave-packet-scattering calculation and a long series of Fourier amplitudes to represent the helium-target potential energy surface. A modified series is constructed in which a truncated Fourier expansion of the potential is constrained to give the exact value of the potential at some key points and which mimics the potential with fewer Fourier amplitudes. The shear horizontal phonon mode is again accessed by the helium scattering for small misalignment of the scattering plane relative to symmetry axes of the monolayer. For 1° misalignment, the calculated intensity of the longitudinal acoustic phonon mode frequently is higher than that of the shear horizontal phonon mode in contrast to what was found at scattering energies near 10 meV for triangular lattices of Ar, Kr, and Xe on Pt(111).     

Predicting stability limits for pure and doped dicationic noble gas clusters undergoing coulomb explosion: A parallel tempering based study

We have used a replica exchange Monte‐Carlo procedure, popularly known as Parallel Tempering, to study the problem of Coulomb explosion in homogeneous Ar and Xe dicationic clusters as well as mixed Ar–Xe dicationic clusters of varying sizes with different degrees of relative composition. All the clusters studied have two units of positive charges. The simulations reveal that in all the cases there is a cutoff size below which the clusters fragment. It is seen that for the case of pure Ar, the value is around 95 while that for Xe it is 55. For the mixed clusters with increasing Xe content, the cutoff limit for suppression of Coulomb explosion gradually decreases from 95 for a pure Ar to 55 for a pure Xe cluster. The hallmark of this study is this smooth progression. All the clusters are simulated using the reliable potential energy surface developed by Gay and Berne (Gay and Berne, Phys. Rev. Lett. 1982, 49, 194). For the hetero clusters, we have also discussed two different ways of charge distribution, that is one in which both positive charges are on two Xe atoms and the other where the two charges are at a Xe atom and at an Ar atom. The fragmentation patterns observed by us are such that single ionic ejections are the favored dissociating pattern. © 2017 Wiley Periodicals, Inc.     

Vibrational excitation of SF6 scattering from graphite

We have observed vibrationally excited sulfur hexafluoride molecules in direct inelastic scattering from hot graphite surfaces. The vibrational temperature for the scattered flux has been determined by probing the effect of internal temperature on electron-induced fragmentation observed in mass spectra. The vibrational excitation depends on incident translational energy, E_tr, and a maximum temperature increase of 50 K is reached in direct scattering at E_tr = 2.5 eV. No effect of surface temperature has been observed at 950–1400 K. Inelastic angular distributions are reproduced by a collision complex model, and the experimental results are related to existing models for vibrational excitation.     

Scattering of O2 from AL111

Experimental results are presented for the scattering of well-defined beams of molecular oxygen incident on clean Al(111). The data consist of scattered angular distributions measured as a function of incident angle, and for fixed incident angle, the dependence on surface temperature of the angular distributions. The measurements are interpreted in terms of a scattering theory that treats the exchange of energy between the translational and rotational motions of the molecule and the phonons of the surface using classical dynamics. The dependence of the measured angular distributions on incident beam angle and temperature is well explained by the theory. Rotational excitation and quantum excitation of the O(2) internal stretching mode are briefly discussed.     

Ti2: accurate determination of the dissociation energy from matrix resonance Raman spectra and chemical interaction with noble gases

UV-visible and resonance Raman spectra of Ti(2) isolated in Ar, Kr, and Xe matrices at temperatures of 10 K were measured by using the 514 nm line of an Ar ion laser. The data show that the Ti(2) molecule interacts strongly with Xe, leading to a significant weakening of the Ti[bond]Ti bond strength. The f(Ti[bond]Ti) force constant decreases in the series Ar>Kr>Xe, from 232.8 Nm(-1) in Ar and 225.5 Nm(-1) in Kr to 199.7 Nm(-1) in Xe. Additional experiments in an Ar matrix containing 2 % of Xe indicate the formation of a molecule of the formula Ti(2)Xe. Our spectra for Ti(2) in an Ar matrix give evidence for several previously not observed members of the Stokes progression. The sum of experimental data allows for an improved estimation of the dissociation energy on the basis of a LeRoy-Bernstein-Lam analysis. A dissociation energy of 1.18 eV was derived from this analysis. The UV-visible data give evidence of the vibrational levels of an excited state of Ti(2).     

Collisions of slow polyatomic ions with surfaces: The scattering method and results

Surface-induced dissociation (SID) and reactions following impact of well-defined ion beams of polyatomic cations C₂H₅OH⁺, CH₄⁺, and CH₅⁺ (and its deuterated variants) at several incident angles and energies with self-assembled monolayers (SAM), carbon surfaces, and hydrocarbon covered stainless steel were investigated by the scattering method. Energy transfer and partitioning of the incident projectile energy into internal excitation of the projectile, translational energy of products, and energy transferred into the surface were deduced from the mass spectra and the translational energy and angular distributions of the product ions. Conversion of ion impact energy into internal energy of the recoiling ions peaked at about 17% of the incident energy for the perfluoro-hydrocarbon SAM, and at about 6% for the other surfaces investigated. Ion survival probability is about 30–50 times higher for closed-shell ions than for open-shell radical cations (e.g., 12% for CD₅⁺ versus 0.3% for CD₄⁺, at the incident angle of 60° with respect to the surface normal). Contour velocity plots for inelastic scattering of CD₅⁺ from hydrocarbon-coated and hydrocarbon-free highly oriented pyrolytic graphite (HOPG) surfaces gave effective masses of the surface involved in the scattering event, approximately matching that of an ethyl group (or two methyl groups) and four to five carbon atoms, respectively. Internal energy effects in impacting ions on SID were investigated by comparing collision energy resolved mass spectra (CERMS) of methane ions generated in a low pressure Nier-type electron impact source versus those generated in a Colutron source in which ions undergo many collisions prior to extraction and are essentially vibrationally relaxed. This comparison supports the hypothesis that internal energy of incident projectile ions is fully available to drive their dissociation following surface impact.     

The interaction of hyperthermal argon atoms with CO-covered Ru(0001): scattering and collision-induced desorption

Hyperthermal Ar atoms were scattered under grazing incidence (θ(i) = 60°) from a CO-saturated Ru(0001) surface held at 180 K. Collision-induced desorption involving the ejection of fast CO (～1 eV) occurs. The angularly resolved in-plane CO desorption distribution has a peak along the surface normal. However, the angular distribution varies with the fractional coverage of the surface. As the total CO coverage decreases, the instantaneous desorption maximum shifts to larger outgoing angles. The results are consistent with a CO desorption process that involves lateral interaction with neighboring molecules. Furthermore, the data indicate that the incident Ar cannot readily penetrate the saturated CO overlayer. Time-of-flight measurements of scattered Ar exhibit two components-fast and slow. The slow component is most evident when scattering from the fully covered surface. The ratio and origin of these components vary with the CO coverage.     

Semiclassical perturbation theory for diffraction in heavy atom surface scattering

The semiclassical perturbation theory formalism of Hubbard and Miller [J. Chem. Phys. 78, 1801 (1983)] for atom surface scattering is used to explore the possibility of observation of heavy atom diffractive scattering. In the limit of vanishing ? the semiclassical theory is shown to reduce to the classical perturbation theory. The quantum diffraction pattern is sensitive to the characteristics of the beam of incoming particles. Necessary conditions for observation of quantum diffraction are derived for the angular width of the incoming beam. An analytic expression for the angular distribution as a function of the angular and momentum variance of the incoming beam is obtained. We show both analytically and through some numerical results that increasing the angular width of the incident beam leads to decoherence of the quantum diffraction peaks and one approaches the classical limit. However, the incoherence of the beam in the parallel direction does not destroy the diffraction pattern. We consider the specific example of Ar atoms scattered from a rigid LiF(100) surface.     

Potential energy surfaces for the benzene–rare gas systems

A harmonic expansion representation of the intermolecular interaction has been exploited to obtain the potential energy surface (PES) for the C₆H₆–He, –Ne, –Ar, –Kr and –Xe systems in an analytical form. Basic data employed are binding energy, equilibrium distance and long-range attraction predicted by a semi-empirical method for selected configurations of the complexes. For those favorable cases where additional information are available the proposed PESs exhibit features in good agreement with those derived from spectroscopy and scattering experiments and/or ab initio calculations. The availability of realistic PESs expressed in an analytical form opens new perspectives of calculations in molecular dynamics and spectroscopic simulations where the benzene molecule and rare gas atoms are involved.     

The N2Ar potential energy surface

The potential energy surface for the N₂Ar system has been obtained assuming a spherical average interaction previously reported from this laboratory. The angular dependence has been assessed by a combined analysis of the integral and differential scattering cross sections and sonic spectroscope data. The potential energy surface is given via a parametric model. A similar potential energy surface for O₂Ar has been obtained with the same procedure. This surface is an improvement of an earlier one, because it reproduces the differential total cross sections recently measured.