Erratum: “Comparative dynamics of Cl(2P) and O(3P) interactions with a hydrocarbon surface” [J. Chem. Phys. 112, 5975 (2000)]

Theoretical study of the O(3P) + C2H6 reaction based on a new ab initio-based global potential energy surface

A new analytical potential energy surface was developed for the first time for the nine-body O(3P) + C2H6 hydrogen abstraction reaction, named PES-2020, which was fitted to explicitly-correlated high-level electronic structure calculations. This surface simulates the topography of the reactive system, from reactants to products, OH(v,j) + C2H5. The adiabatic energy of reaction, ΔHr(0 K) = -2.33 kcal mol-1, reproduces the experimental evidence, and the barrier height, 10.70 kcal mol-1, agrees with the ab initio calculations used as input. In addition, an intermediate complex in the exit channel is observed, which is stabilized with respect to the products of the reaction. Based on PES-2020 a dynamics study was carried out, where quasi-classical trajectory calculations were performed for collision energies in the range of 7.0-60.0 kcal mol-1, which covers high collision energy regions. The reaction cross section increases with collision energy; the largest fraction of available energy is deposited as translational energy (44-66%), and the scattering distribution evolves from backward to forward with collision energy. These findings reproduce previous theoretical calculations using electronic structure calculations of lower levels. However, where these previous studies failed, viz. in rotational and vibrational OH(v,j) distributions, PES-2020 reproduces practically quantitatively the experimental evidence, i.e., cold vibration and rotation, the rotational distribution peaking at j = 1-3 depending on the collision energy. In sum, this behaviour is typical of gas-phase hydrogen abstraction reactions with direct mechanism and high reaction barrier.

Scattering of hyperthermal nitrogen atoms from the Ag(111) surface

Measurements on scattering of hyperthermal N atoms from the Ag(111) surface at temperatures of 500, 600, and 730 K are presented. The scattered atoms have a two-component angular distribution. One of the N components is very broad. In contrast, scattered Ar atoms exhibit only a sharp, single-component angular distribution. There are noteworthy features in the angle-resolved energy of the scattered N when compared with Ar. Taking into account the relative masses involved, N atoms lose significantly more energy at the surface than Ar. However, there is a preferential loss mechanism that predominantly affects low-energy N atoms with small total scattering angle trajectories. The results are interpreted in terms of probing of different interaction potentials: strongly attractive and almost purely repulsive, and spin-state changes during the interaction of N with the surface appear probable.

Theoretical Investigation of Hyperthermal Reactions at the Gas−Liquid Interface: O (3P) and Squalane

Hyperthermal collisions (5 eV) of ground-state atomic oxygen [O ((3)P)] with a liquid-saturated hydrocarbon, squalane (C(30)H(62)), have been studied using QM/MM hybrid "on-the-fly" direct dynamics. The surface structure of the liquid squalane is obtained from a classical molecular dynamics simulation using the OPLS-AA force field. The MSINDO semiempirical Hamiltonian is combined with OPLS-AA for the QM/MM calculations. In order to achieve a more consistent and efficient simulation of the collisions, we implemented a dynamic partitioning of the QM and MM atoms in which atoms are assigned to QM or MM regions based on their proximity to "seed" (open-shell) atoms that determine where bond making/breaking can occur. In addition, the number of seed atoms is allowed to increase or decrease as time evolves so that multiple reactive events can be described. The results show that H abstraction is the most important process for all incident angles, with H elimination, double H abstraction, and C-C bond cleavage also being important. A number of properties of these reactive channels, as well as inelastic nonreactive scattering, are investigated, including angular and translational energy distributions, the effect of incident collision angle, variation with depth of the reactive event within the liquid, with the reaction site on the hydrocarbon, and the effect of dynamics before and after reaction (direct reaction versus trapping reaction-desorption).

Quantum-state resolved reaction dynamics at the gas-liquid interface: Direct absorption detection of HF(v,J) product from F(P2)+Squalane

Exothermic reactive scattering of F atoms at the gas-liquid interface of a liquid hydrocarbon (squalane) surface has been studied under single collision conditions by shot noise limited high-resolution infrared absorption on the nascent HF(v,J) product. The nascent HF(v,J) vibrational distributions are inverted, indicating insufficient time for complete vibrational energy transfer into the surface liquid. The HF(v=2,J) rotational distributions are well fit with a two temperature Boltzmann analysis, with a near room temperature component (T(TD) approximately equal to 290 K) and a second much hotter scattering component (T(HDS) approximately equal to 1040 K). These data provide quantum state level support for microscopic branching in the atom abstraction dynamics corresponding to escape of nascent HF from the liquid surface on time scales both slow and fast with respect to rotational relaxation.

Molecular Dynamics Study to Identify the Reactive Sites of a Liquid Squalane Surface

Molecular dynamics simulations of liquid squalane, C30H62, were performed, focusing in particular on the liquid-vacuum interface. These theoretical studies were aimed at identifying potentially reactive sites on the surface, knowledge of which is important for a number of inelastic and reactive scattering experiments. A united atom force field (Martin, M. G.; Siepmann, J. I. J. Phys. Chem. B 1999, 103, 4508-4517) was used, and the simulations were analyzed with respect to their interfacial properties. A modest but clearly identifiable preference for methyl groups to protrude into the vacuum has been found at lower temperatures. This effect decreases when going to higher temperatures. Additional simulations tracking the flight paths of projectiles directed at a number of randomly chosen surfaces extracted from the molecular dynamics simulations were performed. The geometrical parameters for these calculations were chosen to imitate a typical abstraction reaction, such as the reaction between ground-state oxygen atoms and hydrocarbons. Despite the preference for methyl groups to protrude further into the vacuum, Monte Carlo tracking simulations suggest, on geometric grounds, that primary and secondary hydrogen atoms are roughly equally likely to react with incoming gas-phase atoms. These geometric simulations also indicate that a substantial fraction of the scattered products is likely to undergo at least one secondary collision with hydrocarbon side chains. These results help to interpret the outcome of previous measurements of the internal and external energy distribution of the gas-phase OH products of the interfacial reaction between oxygen atoms and liquid squalane.

Direct Gas−Liquid Interfacial Dynamics: The Reaction between O(3P) and a Liquid Hydrocarbon

We report the first measurements of internal energy distributions of the OH produced via a direct mechanism, isolated from other components on the basis of time-of-flight, in the interfacial reaction between gas-phase O((3)P) atoms and the liquid hydrocarbon squalane, C(30)H(62). O((3)P) atoms were generated by laser photolysis of NO(2) above the liquid. Resulting hydroxyl radicals that escape from the surface were detected by laser-induced fluorescence. Time-of-flight profiles demonstrate that the kinetic energy of the fastest OH (nu' = 1) is lower than that of (nu' = 0). Rotational distributions were measured at the rising edge of their appearance for both OH (nu' = 0) and (nu' = 1). They were found to differ substantially more than at the peak of their profiles. They were also less dependent on the bulk liquid temperature. We conclude that the new data confirm strongly that at least two mechanisms contribute to the production of OH. The higher-velocity component has translational and rotational energy distributions, observed cleanly for the first time, consistent with a direct mechanism. The close correspondence of these rotational distributions to those from the corresponding homogeneous gas-phase reaction of O((3)P) with smaller hydrocarbons suggests a very similar, near collinear direct abstraction. This is accompanied by a slower component with kinetic energy and rotational (but not vibrational) distributions reflecting the temperature of the liquid, consistent with a distinct trapping-desorption mechanism.

The effects of surface temperature on the gas-liquid interfacial reaction dynamics of O(3P)+squalane

OH/OD product state distributions arising from the reaction of gas-phase O(3P) atoms at the surface of the liquid hydrocarbon squalane C30H62/C30D62 have been measured. The O(3P) atoms were generated by 355 nm laser photolysis of NO2 at a low pressure above the continually refreshed liquid. It has been shown unambiguously that the hydroxyl radicals detected by laser-induced fluorescence originate from the squalane surface. The gas-phase OH/OD rotational populations are found to be partially sensitive to the liquid temperature, but do not adapt to it completely. In addition, rotational temperatures for OH/OD(v'=1) are consistently colder (by 34+/-5 K) than those for OH/OD(v'=0). This is reminiscent of, but less pronounced than, a similar effect in the well-studied homogeneous gas-phase reaction of O(3P) with smaller hydrocarbons. We conclude that the rotational distributions are composed of two different components. One originates from a direct abstraction mechanism with product characteristics similar to those in the gas phase. The other is a trapping-desorption process yielding a thermal, Boltzmann-like distribution close to the surface temperature. This conclusion is consistent with that reached previously from independent measurements of OH product velocity distributions in complementary molecular-beam scattering experiments. It is further supported by the temporal profiles of OH/OD laser-induced fluorescence signals as a function of distance from the surface observed in the current experiments. The vibrational branching ratios for (v'=1)/(v'=0) for OH and OD have been found to be (0.07+/-0.02) and (0.30+/-0.10), respectively. The detection of vibrationally excited hydroxyl radicals suggests that secondary and/or tertiary hydrogen atoms may be accessible to the attacking oxygen atoms.

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.