A QM/MM model for O( 3 P) reaction with an alkyl thiolate self-assembled monolayer

Computational approaches to dissociative chemisorption on metals: towards chemical accuracy

We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (E_b, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate E_b can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol^-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for E_b for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH₄. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of E_b for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate E_b for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O₂, H₂O, NH₃, CO₂, and CH₃OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.

Reactive Scattering as a Chemically Specific Analytical Probe of Liquid Surfaces

In this Perspective, we highlight some recent progress in the reactive scattering of "chemical probe" species such as atoms or small radicals from liquid surfaces. We emphasize in particular the evolution of this area from purely dynamical studies of the scattering mechanism. The mechanistic understanding that has now been gained is sufficiently mature to allow the same methods to be used as an effective analytical tool. The use of this approach to measure liquid-surface composition and structure is illustrated through the scattering of O((3)P) atoms from a common, imidazolium-based family of ionic liquids.

Collisions of noble gases with supercooled sulfuric acid-water solutions

The collisions of hyperthermal noble gases (He, Ne, Ar, Kr, Xe) with supercooled binary sulfuric acid-water mixtures (57-77 wt%) were explored in the temperature range between 210 and 240 K. The experiments were performed by directing a molecular beam of the respective gases onto a continuously renewed liquid surface and monitoring the velocity of the scattered molecules by mass spectrometry. Depending on the initial translational energies and molecular masses, we observe both inelastic scattering from the surface as well as thermalization followed by subsequent desorption. The experiments indicate that the repulsive momentum transfer in the inelastic scattering channel increases with increasing mass of the impinging gas, while it is only weakly affected by the initial velocities. The final energy of the thermally desorbing atoms can always be approximated by a Maxwell-Boltzmann distribution equal to the liquid bulk phase temperature. The influence of the binary composition of the liquid phase is only noticeable in the case of Ne, whilst this dependence diminishes for gases with molecular masses >or=40 amu. The probability of thermalisation relative to inelastic scattering increases with the bulk phase temperature, independent of the molecular masses of the colliding gas. In contrast, the fractional energy transfer during collision does not increase with temperature, except for Neon. These results can be interpreted in the model framework of hard-sphere collisions of noble gases with the surface, during which water and sulfuric acid molecules interact independently with the impinging gas.

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).

Theoretical Study of the Dynamics of Ar Collisions with C2H6 and C2F6 at Hyperthermal Energy

We present a classical-trajectory study of the dynamics of high-energy (5-12 eV) collisions between Ar atoms and the C2H6 and C2F6 molecules. We have constructed the potential-energy surfaces for these systems considering separately the Ar-molecule interactions (intermolecular potential) and the interactions within the molecule (intramolecular potential). The intermolecular surfaces consist of pairwise empirical potentials derived from high-accuracy ab initio calculations. The intramolecular potentials for C2H6 and C2F6 are described using specific-reaction-parameters semiempirical Hamiltonians and are calculated "on the fly", i.e., while the trajectories are evolving. Trajectory analysis shows that C2F6 absorbs more energy than C2H6 and is more susceptible to collision-induced dissociation (CID). C-C bond-breakage processes are more important than C-H or C-F bond breakage at the energies explored in this work. Analysis of the reaction mechanism for CID processes indicates that, although C-C breakage is mostly produced by side-on collisions, head-on collisions are more efficient in producing C-F or C-H dissociation. Our results suggest that high-energy collisions between closed-shell species of the natural low-Earth-orbit environment and spacecraft can contribute to the observed degradation of polymers that coat spacecraft surfaces.

Diffusion on a Self-Assembled Monolayer: Molecular Modeling of a Bound + Mobile Lubricant

The diffusion of tricresyl phosphate molecules on an octadecyltrichlorosilane self-assembled monolayer (SAM) was characterized using molecular dynamics simulations. The simulations predict that when placed on the top of a close-packed SAM, the molecules remain mobile on the surface with an isotropic diffusion activation energy of approximately 9 kJ/mol. In contrast, an anisotropic barrier that results from chain tilt within the SAM is predicted for diffusion into a defect created by reducing the alkane chain length within a cylinderical region of the surface. Once in the defect, the molecules become trapped by embedding part of the molecule into the side of the SAM.

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.

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.

Ab initio and analytic intermolecular potentials for Ar-CF4

Ab initio calculations at the CCSD(T) level of theory were performed to characterize the Ar + CF4 intermolecular potential. Potential energy curves were calculated with the aug-cc-pVTZ basis set, and with and without a correction for basis set superposition error (BSSE). Additional calculations were performed with other correlation consistent basis sets to extrapolate the Ar-CF4 potential energy minimum to the complete basis set (CBS) limit. Both the size of the basis set and BSSE have substantial effects on the Ar + CF4 potential. Calculations with the aug-cc-pVTZ basis set, and with a BSSE correction, appear to give a good representation of the BSSE corrected potential at the CBS limit. In addition, MP2 theory is found to give potential energies in very good agreement with those determined by the much higher level CCSD(T) theory. Two model analytic potential energy functions were determined for Ar + CF4. One is a fit to the aug-cc-pVTZ calculations with a BSSE correction. The second was derived by fitting an average BSSE corrected potential, which is an average of the CCSD(T)/aug-cc-pVTZ potentials with and without a BSSE correction. These analytic functions are written as a sum of two-body potentials and excellent fits to the ab initio potentials are obtained by representing each two-body interaction as a Buckingham potential.

Reaction of 5−40 eV Ions with Self-Assembled Monolayers

The reaction of 5-40 eV O(+) and Ne(+) ions with alkanethiolate and semifluorinated alkanethiolate self-assembled monolayers (SAMs) is studied under ultrahigh vacuum (UHV) conditions. Whereas Ne(+) simply sputters fragments from the surface, O(+) can also abstract surface atoms and break C-C bonds in both the hydrocarbon and fluorocarbon SAM chains. Isotopic labeling experiments reveal that O(+) initially abstracts hydrogen atoms from the outermost two carbon atoms on an alkanethiolate SAM chain. However, the position of the isotopic label quickly becomes scrambled along the chain as the SAM is damaged through continuous ion bombardment. Scanning tunneling microscopy (STM) monitors changes in the SAM conformational structure at various stages during 5 eV ion bombardment. STM images indicate that O(+) reacts less efficiently with dodecanethiolate molecules packed internally within a structural domain than it does with molecules adsorbed at domain boundaries or near defect sites. STM images recorded after Ne(+) bombardment suggest that Ne(+) attacks the SAM exclusively near the domain boundaries. Taken collectively, these experiments advance our understanding of the degradation pathways suffered by polymeric satellite materials in the low-earth orbit (LEO) space environment.

Ab initioand direct quasiclassical-trajectory study of the F+CH4→HF+CH3 reaction

We present an electronic structure and dynamics study of the F+CH4-->HF+CH3 reaction. CCSD(T)/aug-cc-pVDZ geometry optimizations, harmonic-frequency, and energy calculations indicate that the potential-energy surface is remarkably isotropic near the transition state. In addition, while the saddle-point F-H-C angle is 180 degrees using MP2 methods, CCSD(T) geometry optimizations predict a bent transition state, with a 153 degrees F-H-C angle. We use these high-quality ab initio data to reparametrize the parameter-model 3 (PM3) semiempirical Hamiltonian so that calculations with the improved Hamiltonian and employing restricted open-shell wave functions agree with the higher accuracy data. Using this specific-reaction-parameter PM3 semiempirical Hamiltonian (SRP-PM3), we investigate the reaction dynamics by propagating quasiclassical trajectories. The results of our calculations using the SRP-PM3 Hamiltonian are compared with experiments and with the estimates of two recently reported potential-energy surfaces. The trajectory calculations using the SRP-PM3 Hamiltonian reproduce quantitatively the measured HF vibrational distributions. The calculations also agree with the experimental HF rotational distributions and capture the essential features of the excitation function. The results of the SRP semiempirical Hamiltonian developed here clearly improve over those using the two prior potential-energy surfaces and suggest that reparametrization of semiempirical Hamiltonians is a promising strategy to develop accurate potential-energy surfaces for reaction dynamics studies of polyatomic systems.

Surface Reactions of Molecular and Atomic Oxygen with Carbon Phosphide Films

The surface reactions of atomic and molecular oxygen with carbon phosphide films have been studied using X-ray photoelectron spectroscopy (XPS) and atomic force microscopy (AFM). Carbon phosphide films were produced by ion implantation of trimethylphosphine into polyethylene. Atmospheric oxidation of carbon phosphide films was dominated by phosphorus oxidation and generated a carbon-containing phosphate surface film. This oxidized surface layer acted as an effective diffusion barrier, limiting the depth of phosphorus oxidation within the carbon phosphide film to < 3 nm. The effect of atomic oxygen (AO) exposure on this oxidized carbon phosphide layer was subsequently probed in situ using XPS. Initially AO exposure resulted in a loss of carbon atoms from the surface, but increased the surface concentration of phosphorus atoms as well as the degree of phosphorus oxidation. For more prolonged AO exposures, a highly oxidized phosphate surface layer formed that appeared to be inert toward further AO-mediated erosion. By utilizing phosphorus-containing hydrocarbon thin films, the phosphorus oxides produced during exposure to AO were found to desorb at temperatures >500 K under vacuum conditions. Results from this study suggest that carbon phosphide films can be used as AO-resistant surface coatings on polymers.

Quasiclassical Trajectory Study of the O(3P) + CH4→ OH + CH3Reaction with a Specific Reaction Parameters Semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.

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.

Transition state dynamics and a QM/MM model for the Cl– + C2H5Cl SN2 reaction

A MP2/6-31G* direct dynamics simulation is used to study the dynamics of the central barrier [Cl-C2H5-Cl]– for the Cl– + C2H5 SN2 reaction. The majority of the trajectories move off the central barrier to form the Cl––C2H5Cl complex and appear to undergo efficient IVR as assumed by RRKM theory. However, some of the trajectories move directly to products without forming the complex, a non-RRKM result. A hydrogen atom link-atom QM/MM model is described for studying the dynamics of [X-CH2R-Y]– central barriers with the -R substituent. The model is used to calculate vibrational frequencies for the [Cl-C2H5-Cl]– central barrier.Key words: SN2 reaction dynamics, RRKM theory, QM/MM model, central barrier dynamics, direct dynamics classical trajectories.

Energetics, transition states, and intrinsic reaction coordinates for reactions associated with O(3P) processing of hydrocarbon materials

Electronic structure calculations based on multiconfiguration wave functions are used to investigate a set of archetypal reactions relevant to O(3P) processing of hydrocarbon molecules and surfaces. These include O(3P) reactions with methane and ethane to give OH plus methyl or ethyl radicals, O(3P) + ethane to give CH3O + CH3, and secondary reactions of the OH product radical with ethane and the ethyl radical. Geometry optimization is carried out with CASSCF/cc-pVTZ for all reactions, and with CASPT2/cc-pVTZ for O(3P) + methane/ethane. Single-point energy corrections are applied with CASPT2, CASPT3, and MRCI + Q with the cc-pVTZ and cc-pVQZ basis sets, and the energies extrapolated to the complete basis set limit (CBL). Where comparison of computed barriers and energies of reaction with experiment is possible, the agreement is good to excellent. The best agreement (within experimental error) is found for MRCI + Q/CBL applied to O(3P) + methane. For the other reactions, CASPT2/CBL and MRCI + Q/CBL predictions differ from experiment by 1-5 kcal/mol for 0 K enthalpies of reaction, and are within 1 kcal/mol of the best-estimate experimental range of 0 K barriers for O(3P) + ethane and OH + ethane. The accuracy of MRCI + Q/CBL is limited mainly by the quality of the active space. CASPT2/CBL barriers are consistently lower than MRCI + Q/CBL barriers with identical reference spaces.

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.