Quasiclassical trajectory study of energy transfer and collision-induced dissociation in hyperthermal Ar + CH4 and Ar + CF4 collisions

We present a study of energy transfer in collisions of Ar with methane and perfluoromethane at hyperthermal energies (E(coll) = 4-10 eV). Quasiclassical trajectory calculations of Ar + CX(4) (X = H, F) collisions indicate that energy transfer from reagents' translation to internal modes of the alkane molecule is greatly enhanced by fluorination. The reasons for the enhancement of energy transfer upon fluorination are shown to emerge from a decrease in the hydrocarbon vibrational frequencies of the CX(4) molecule with increasing the mass of the X atom, and to an increase of the steepness of the Ar-CX(4) intermolecular potential. At high collision energies, we find that the cross section of Ar + CF(4) collisions in which the amount of energy transfer is larger than needed to break a C-F bond is at least 1 order of magnitude larger than the cross sections of Ar + CH(4) collisions producing CH(4) with energy above the dissociation limit. In addition, collision-induced dissociation is detected in short time scales in the case of the fluorinated species at E(coll) = 10 eV. These results suggest that the cross section for degradation of fluorinated hydrocarbon polymers under the action of nonreactive hyperthermal gas-phase species might be significantly larger than that of hydrogenated hydrocarbon polymers. We also illustrate a practical way to derive intramolecular potential energy surfaces for bond-breaking collisions by improving semiempirical Hamiltonians based on grids of high-quality ab initio calculations.     

Quasiclassical trajectory study of the vibrational quenching of hydroxyl radicals through collision with O atoms

The collisional removal of vibrationally excited OH radicals by O atoms is studied by the quasiclassical trajectory method. To evaluate the effect of different topological features on the scattering processes two different global potential energy surfaces, DMBE IV and TU, are used. Results for reactive, exchange, and inelastic scattering probabilities are reported for central collisions (with zero total angular momentum) with a fixed relative translational energy for vibrational levels of OH ranging from nu=1 to v=8. Vibrational state distributions of product molecules are also compared on the two potential energy surfaces. Both surfaces predict higher probabilities for reaction than for exchange or inelastic scattering. The vibrational state distributions of the product diatomic molecules are different on the two surfaces. In particular, the two surfaces give substantially different probabilities for multiquantum OH vibrational relaxation transitions OH(v)+O-->OH(v')+O.     

Quasiclassical trajectory study of the Cl + Hl → HCl + I reaction

Monte Carlo quasiclassical trajectory calculations have been carried out for the reaction Cl + Hl → HCl + I for 300, 1000, and 2000 K. A semi-empirical potential-energy surface (London equation) was obtained by “transfering” parameters from surfaces computed for other reaction systems. The computed results are in general accord with experimental measurements. Thermal rate coefficients, differential scattering cross sections, and product vibrational and rotational distributions were computed for the three temperatures. Angular scattering distributions are in agreement with experiment only at elevated temperatures.     

Quasiclassical trajectory study of the collision-induced dissociation of CH3SH+ + Ar

Quasiclassical trajectory calculations were carried out to study the dynamics of energy transfer and collision-induced dissociation (CID) of CH(3)SH(+) + Ar at collision energies ranging from 4.34 to 34.7 eV. The relative abundances calculated for the most relevant product ions are found to be in good agreement with experiment, except for the lowest energies investigated. In general, the dissociation to form CH(3)(+) + SH is the dominant channel, even though it is not among the energetically favored reaction pathways. The results corroborate that this selective dissociation observed upon collisional activation arises from a more efficient translational to vibrational energy transfer for the low-frequency C-S stretching mode than for the high-frequency C-H stretching modes, together with weak couplings between the low- and high-frequency modes of vibration. The calculations suggest that CID takes place preferentially by a direct CH(3)(+) + SH detachment, and more efficiently when the Ar atom collides with the methyl group-side of CH(3)SH(+).     

Classical trajectory study of rotational excitation of HD by collisions with He

The results of a classical trajectory study of rotational excitation of HD by collisions with He are presented. The calculations spanned the relative translational energy range 0.26 – 1.80 kcal/mole for transitions from j = 0 HD and 0.52 – 1.80 kcal/mole for transitions from j = 1 HD. Three methods of defining initial states were employed and the results of each method are compared with the accurate close coupling calculations of Green.     

Quasiclassical trajectory study of reactive and dissociative processes in H2+H2: comparison with quantum-mechanical calculations

Quasiclassical trajectory calculations have been carried out for H(2)(v(1)=high)+H(2)(v(2)=low) collisions within a three degrees of freedom model where five different geometries of the colliding complex were considered. Within this approach, probabilities for different competitive processes are studied: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. The purpose is to compare in detail with equivalent quantum-mechanical wave packet calculations [Bartolomei et al., J. Chem. Phys 122, 064305 (2005)], especially the behavior of the probabilities near reaction thresholds. Quasiclassical calculations compare quite well with the quantum-mechanical ones for collision induced dissociation as well as for the four center reaction, although quantum effects become very important near thresholds, particularly for lower v(1)'s and for the four center process. Less quantitative agreement is found for reactive dissociation and three-body complex formation. It is found that most quantum effects are due to differences between quantum and classical vibrational distributions of H(2)(v(1)=high). Zero point energy violation has been found in the classical reactive-dissociative probabilities. Extension of these findings to full-dimensional treatments is examined.     

Quasiclassical trajectory study of the CH3++HD-->CH2D++H2 Reaction

A full dimensional ab initio potential energy surface for the CH5+ system based on coupled cluster electronic structure calculations and capable of describing the dissociation of methonium ion into methyl cation and molecular hydrogen (J. Phys. Chem. A 2006, 110, 1569) is used in quasiclassical trajectory calculations of the reaction CH3++HD-->CH2D++H2 for low collision energies of relevance to astrochemistry. Cross sections for the exchange are obtained at several relative translational energies and a fit to the energy dependence of the cross sections is used to obtain the rate constant at temperatures between 10 and 50 K. The calculated rate constant at 10 K agrees well with the previously reported experimental value. Internal energy distributions of the products are presented and discussed in the context of zero-point energy "noncompliance".     

Quasiclassical trajectory study of the F + CH4 reaction dynamics on a dual-level interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the F + CH4 reactive system has been constructed using the interpolation method of Collins and co-workers. The ab initio calculations have been performed using second-order M?ller-Plesset (MP2) perturbation theory to build the initial PES. Scaling all correlation (SAC) methodology has been employed to improve the ab initio calculations and to construct a dual-level PES. Using this PES, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations and internal energy distributions has been carried out for the F + CH4 and F + CD4 reactions and the theoretical results have been compared with the available experimental data.     

Classical trajectory study of collisions of Ar with alkanethiolate self-assembled monolayers: potential-energy surface effects on dynamics

We have investigated collisions between Ar and alkanethiolate self-assembled monolayers (SAMs) using classical trajectory calculations with several potential-energy surfaces. The legitimacy of the potential-energy surfaces is established through comparison with molecular-beam data and ab initio calculations. Potential-energy surfaces used in previous work overestimate the binding of Ar to the SAM, leading to larger energy transfer than found in the experiments. New calculations, based on empirical force fields that better reproduce ab initio calculations, exhibit improved agreement with the experiments. In particular, polar-angle-dependent average energies calculated with explicit-atom potential-energy surfaces are in excellent agreement with the experiments. Polar- and azimuthal-angle-dependent product translational energies are examined to gain deeper insight into the dynamics of Ar+SAM collisions.     

Quasiclassical trajectory study of colinear Cl− + HBr → HCl + Br−

Energy distributions in the products of the ion-molecule reaction Cl^? + HBr → HC1 + Br^? have been studied using quasiclassical trajectories on a semi-empirical collinear potential energy surface. Vibrational energy is favored in the products. While some trajectories are long-lived, the kinematic factor of the light central atom prevents effective energy redistribution.     

Quasiclassical trajectory study of the rotational distribution for the O+NO(v = 0) fundamental vibrational excitation

Quasiclassical trajectories have been run to study the fundamental one‐quantum vibrational transition formed from collisions of ground‐state nitric oxide with atomic oxygen at temperatures of 500, 750, and 1000 K. Two adiabatic potential energy surfaces of different symmetry (²A′ and ²A″ of NO₂) have been utilized. The rate constant for the title process is given along with the rotational distributions, and the results shown to corroborate previous atmospheric models that describe the nascent state by a Maxwell–Boltzmann distribution at the local temperature. © 2011 Wiley Peiodicals, Inc. Int J Chem Kinet 43: 345–352, 2011     

Ion-pair formation in fast collisions of metastable argon with iodine

Differential ion-pair production cross sections for Ar^*—I₂ collisions at center of mass energies of 25–133 eV are quite similar to analogous alkali—I₂ results despite the addition of competing channels (Penning ionization and excitation transfer) and possible multiple electronic surface crossings.     

Quasiclassical trajectory state-to-state cross sections for energy transfer in Ar+ OH(υ = 9,J = 0,4, and 8) collisions

State-to-state energy transfer cross sections have been computed for Ar + OH(υ₁ = 9; J_i = 0,4, and 8) at initial relative translational energy 0.2 eV using quasiclassical trajectories. The results show that a relatively small amount of initial rotational excitation has a significant effect on the energy transfer. The energy transfer for J_i = 8 is dominated by transitions for which the vibrational and rotational energy changes are such that the translational energy changes are minimized.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

A quasiclassical trajectory study of vibrational energy transfer in collisions involving intermolecular attraction of moderate strength

Monte Carlo selected, quasiclassical trajectories have been computed on six potential energy hypersurfaces possessing potential minima or “wells” up to 50 kJ mol^?1 deep. The aim of the investigation has been to examine how vibrational energy transfer in A + BC(υ = 1) collisions is promoted by intermolecular attraction of moderate strength. Here results are reported for the mass combination m_A = 20 u, m_B = 1 u, m_C = u. The results show that even quite slight intermolecular attraction can enhance energy transfer, as long as the attraction does not just depend on the separation of A from the center-of-mass of BC. The mean loss of vibrational energy does not depend only the well depth but also on its “location” (in particular, the difference in r_BC at the minimum and in isolated BC) and on the angular anisotropy of the potential. Large transfers of energy do not occur only in complex-forming collisions; indeed, a high fraction of trajectories on all surfaces are direct but show similar transfer of energy as in the more complex trajectories on the same surface. The results of the calculations are discussed in relation to the mechanisms and rates of vibrational relaxation in collisions between radicals and between species. such as HF + HF, capable of forming hydrogen bonds.     

Quasiclassical trajectory studies of H(D)+HBr(DBr) abstraction and exchange reactions

The abstraction and exchange reaction dynamics for H(D)+HBr(DBr) systems have been investigated on three LEPS potential-energy surfaces whose features are in accord with the surface topography suggested by recent molecular-beam and thermal experiments (abstraction barrier less than 1.0 kcal/mole, exchange reaction barriers of =5.0 kcal/mole, and no attractive wells with a depth greater than 0.209 kcal/mole). The surfaces differ primarily in the magnitude of the abstraction barrier which varies from 0.19 to 1.01 kcal/mole. Reaction cross sections have been computed on each surface as a function of relative collision energy from the results of 139000 quasiclassical trajectoris. Comparison of these results with measured relative abstraction cross sections suggests that the true abstraction barrier is very small, perhaps between 0.0 and 0.25 kcal/mole. However, thermal rate coefficients computed on the - best- surface at 300 K are about a factor of 2 larger than the most recently measured values. The calculated (H,D)/(D,H) isotope ratio at 300 K lies between the two reported experimental results. The computed thermal activation energy for abstraction is 835 cal/mole, which is in good agreement with a very early measurements but a factor of 2.5 less than the most recently reported experimental result. These results suggest that the molecular-beam and thermal rate measurements are inconsistent. The average fraction of the available energy which is partitioned into internal product modes <f_E> is found to be nearly independent of relative collision energy and the small topographical differences present in the potential surfaces used in these calculations. We find <f_E> = 0.40. In all reactions, the differential scattering cross sections are peaked in the backward direction for the molecular products, indicating a rebound mechanism.     

Quasiclassical trajectory study of formaldehyde unimolecular dissociation: H(2)CO-->H2 + CO, H + HCO

We report quasiclassical trajectory calculations of the dynamics of the two reaction channels of formaldehyde dissociation on a global ab initio potential energy surface: the molecular channel H(2)CO-->H(2) + CO and the radical H(2)CO-->H + HCO. For the molecular channel, it is confirmed that above the threshold of the radical channel a second, intramolecular hydrogen abstraction pathway is opened to produce CO with low rotation and vibrationally hot H(2). The low-j(CO) and high-nu(H(2) ) products from the second pathway increase with the total energy. The competition between the molecular and radical pathways is also studied. It shows that the branching ratio of the molecular products decreases with increasing energy, while the branching ratio of the radical products increases. The results agree well with very recent velocity-map imaging experiments of Suits and co-workers and solves a mystery first posed by Moore and co-workers. For the radical channel, we present the translational energy distributions and HCO rotation distributions at various energies. There is mixed agreement with the experiments of Wittig and co-workers, and this provides an indirect confirmation of their speculation that the triplet surface plays a role in the formation of the radical products.     

Second-order theories for electron loss in fast collisions with He atoms

The momentum distribution of projectile electrons ejected in collisions with light targets is calculated within the second-order Born approximation for direct ionisation and within the electron impact approximation and the impulse approximation for electron capture to the target continuum. From comparison with available experimental data it is found that for forward emission angles the electron is well described by a projectile eigenstate, while at backward angles a target final state is more appropriate. At all angles the inclusion of simultaneous target excitation is very important.