Quasiclassical trajectory study of the collision-induced dissociation dynamics of Ar + CH3SH+ using an ab initio interpolated potential energy surface

Classical trajectory calculations have been performed to investigate the collision-induced dissociation (CID) of the CH(3)SH(+) cation with Ar atoms. A new intramolecular potential energy surface for the CH(3)SH(+) cation is evaluated by interpolation of 3000 ab initio data points calculated at the MP2/6-311G(d,p) level of theory. The new potential energy surface includes seven accessible dissociation channels of the cation. The present QCT calculations show that migration of hydrogen atoms, leading to the rearrangement CH(3)SH(+) <--> CH(2)SH(2)(+), is significant at the collision energies considered (6.5-34.7 eV) and that the formation of CH(3)(+), CH(3)S(+), and CH(2)(+) cations takes place primarily by a "shattering" mechanism in which the products are formed just after the collision. The theoretical product abundances are found to be in qualitative agreement with the experimental data. However, at high collision energies, the calculated total cross sections for the formation of CH(3)(+) and CH(2)SH(+) cations are noticeably larger than the experimental determinations. Several features of the dynamics of the CID processes are discussed.     

Trajectory dynamics study of collision-induced dissociation of the Ar + CH4 reaction at hyperthermal conditions: vibrational excitation and isotope substitution

We investigate the role of vibrational energy excitation of methane and two deuterated species (CD(4) and CH(2)D(2)) in the collision-induced dissociation (CID) process with argon at hyperthermal energies. The quasi-classical trajectory method has been applied, and the reactive Ar + CH(4) system has been modeled by using a modified version of the CH(4) potential energy surface of Duchovic et al. (J. Phys. Chem. 1984, 88, 1339) and the Ar-CH(4) intermolecular potential function obtained by Troya (J. Phys. Chem. A 2005, 109, 5814). This study clearly shows that CID is markedly enhanced with vibrational excitation and, to a lesser degree, with collision energy. In general, CID increases by exciting stretch vibrational modes of the reactant molecule. For the direct dissociation of CH(4), however, the CID cross sections appear to be essentially independent of which vibrational mode is initially excited. In all situations studied, the CID cross sections are always greater for the Ar + CD(4) reaction than for the Ar + CH(4) one, the Ar + CH(2)D(2) being an intermediate situation. A detailed analysis of the energy transfer processes, including their relation with CID, is also presented.     

An ab initio based global potential energy surface describing CH5+ --> CH3+ + H2

A full-dimensional, ab initio based potential energy surface (PES) for CH(5)(+), which can describe dissociation is reported. The PES is a precise fit to 36173 coupled-cluster [CCSD(T)] calculations of electronic energies done using an aug-cc-pVTZ basis. The fit uses a polynomial basis that is invariant with respect to permutation of the five H atoms, and thus describes all 120 equivalent minima. The rms fitting error is 78.1 cm(-1) for the entire data set of energies up to 30,000 cm(-1) and a normal-mode analysis of CH(5)(+) also verifies the accuracy of the fit. Two saddle points have been located on the surface as well and compared with previous theoretical work. The PES dissociates correctly to the fragments CH(3)(+) + H(2) and the equilibrium geometry and normal-mode analyses of these fragments are also presented. Diffusion Monte Carlo calculations are done for the zero-point energies of CH(5)(+) (and some isotopologs) as well as for the separated fragments of CH(5)(+), CH(3)(+) + H(2) and those of CH(4)D(+), CH(3)(+) + HD and CH(2)D(+) + H(2). Values of D(0) are reported for these dissociations. A molecular dynamics calculation of CH(4)D(+) dissociation at one total energy is also performed to both validate the applicability of the PES for dynamics studies as well as to test a simple classical statistical prediction of the branching ratio of the dissociation products.     

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.     

Dissociation of energy selected Sn(CH3)4(+), Sn(CH3)3Cl+, and Sn(CH3)3Br+ ions: evidence for isolated excited state dynamics

The dissociation dynamics of Sn(CH(3))(4)(+), Sn(CH(3))(3)Cl(+), and Sn(CH(3))(3)Br(+) were investigated by threshold photoelectron photoion spectrometry using an electron imaging apparatus (iPEPICO) at the Swiss Light Source. The tetramethyltin ion was found to dissociate via Sn(CH(3))(4)(+) → Sn(CH(3))(3)(+) + CH(3) → Sn(CH(3))(2)(+) + 2CH(3), while the trimethyltin halide ions dissociated via methyl loss at low energies, and a competitive halogen loss at somewhat higher energies. The 0 K methyl loss onset for the three ions was found to be 9.410 ± 0.020 eV, 10.058 ± 0.020 eV, and 9.961 ± 0.020 eV, respectively. Statistical theory could not reproduce the observed onsets for the halogen loss steps in the halotrimethyltin ions. The halide loss signal as a function energy mimicked the excited state threshold photoelectron spectrum, from which we conclude that the halide loss from these ions takes place on an isolated excited state potential energy surface, which we describe by time dependent density functional calculations. The sequential loss of a second methyl group in the Sn(CH(3))(4)(+) ion, observed at about 3 eV higher energies than the first one, is also partially non-statistical. The derived product energy distribution resulting from the loss of the first methyl group is two-component with about 50% being statistical and the remainder associated with high translational energy products that peak at 2 eV. Time dependent DFT calculations show that a dissociative ?B state lies in the vicinity of the experimental measurements. We thus propose that 50% of the Sn(CH(3))(4)(+) ions produced in this energy range internally convert to the ?X state, on which they dissociate statistically, while the remainder dissociate directly from the repulsive ?B state leading to high kinetic energy products.     

On the parallel mechanism of the dissociation of energy-selected P(CH3)3+ ions

Energy selected trimethyl phosphine ions were prepared by threshold photoelectron photoion coincidence (TPEPICO) spectroscopy. This ion dissociates via H, CH(3), and CH(4) loss, the latter two involving hydrogen transfer steps. The ion time-of-flight distribution and the breakdown diagram are analyzed in terms of the statistical RRKM theory, which includes tunneling. Ab initio and DFT calculations provide the vibrational frequencies required for the RRKM modeling. CH(3) loss could produce both the P(CH(3))(2)(+) by a simple bond dissociation step, and the more stable HP(CH(2))CH(3)(+) ion by a hydrogen transfer step. Quantum chemical calculations are extensively used to uncover the reaction scheme, and they strongly suggest that the latter product is exclusively formed via an isomerization step in the energy range of the experiment. The data analysis, which includes modeling with the trimethyl phosphine thermal energy distribution, provides accurate onset energies for both H (E(0K) = 1024.1 +/- 3.5 kJ/mol) and CH(3) (E(0K) = 1024.8 +/- 3.5 kJ/mol) loss reactions. From this analysis, we conclude that the Delta(f)H(298K) degrees [HP(CH(2))(CH(3))(+)] = 783 +/- 8 kJ/mol and Delta(f)H(298K) degrees [P(CH(2))(CH(3))(2)(+)] = 711 +/- 8 kJ/mol.     

Potential energy surface for activation of methane by Pt(+): a combined guided ion beam and DFT study

A guided-ion beam tandem mass spectrometer is used to study the reactions of Pt(+) with methane, PtCH(2)(+) with H(2) and D(2), and collision-induced dissociation of PtCH(4)(+) and PtCH(2)(+) with Xe. These studies experimentally probe the potential energy surface for the activation of methane by Pt(+). For the reaction of Pt(+) with methane, dehydrogenation to form PtCH(2)(+) + H(2) is exothermic, efficient, and the only process observed at low energies. PtH(+), formed in a simple C-H bond cleavage, dominates the product spectrum at high energies. The observation of a PtH(2)(+) product provides evidence that methane activation proceeds via a (H)(2)PtCH(2)(+) intermediate. Modeling of the endothermic reaction cross sections yields the 0 K bond dissociation energies in eV (kJ/mol) of D(0)(Pt(+)-H) = 2.81 +/- 0.05 (271 +/- 5), D(0)(Pt(+)-2H) = 6.00 +/- 0.12 (579 +/- 12), D(0)(Pt(+)-C) = 5.43 +/- 0.05 (524 +/- 5), D(0)(Pt(+)-CH) = 5.56 +/- 0.10 (536 +/- 10), and D(0)(Pt(+)-CH(3)) = 2.67 +/- 0.08 (258 +/- 8). D(0)(Pt(+)-CH(2)) = 4.80 +/- 0.03 eV (463 +/- 3 kJ/mol) is determined by measuring the forward and reverse reaction rates for Pt(+) + CH(4) right harpoon over left harpoon PtCH(2)(+) + H(2) at thermal energy. We find extensive hydrogen scrambling in the reaction of PtCH(2)(+) with D(2). Collision-induced dissociation (CID) of PtCH(4)(+), identified as the H-Pt(+)-CH(3) intermediate, with Xe reveals a bond energy of 1.77 +/- 0.08 eV (171 +/- 8 kJ/mol) relative to Pt(+) + CH(4). The experimental thermochemistry is favorably compared with density functional theory calculations (B3LYP using several basis sets), which also establish the electronic structures of these species and provide insight into the reaction mechanism. Results for the reaction of Pt(+) with methane are compared with those for the analogous palladium system and the differences in reactivity and mechanism are discussed.     

Structure of the oxygen-rich cluster cation Al2O7+ and its reactivity toward methane and water

The oxygen-rich cluster Al(2)O(7)(+) is generated in the gas phase and investigated with respect to both its structure and its reactivity toward small, inert molecules using Fourier-transform ion cyclotron resonance (FT-ICR) mass spectrometry and DFT-based calculations. Al(2)O(7)(+) reacts with CH(4) under ambient conditions via hydrogen atom transfer (HAT), and with H(2)O a ligand exchange occurs which gives rise to the evaporation of two O(2) molecules. The resulting product ion Al(2)O(4)H(2)(+) is also capable of abstracting a hydrogen atom from both H(2)O and CH(4). As indicated in the H(2)O/2O(2) ligand exchange and supported by collision-induced dissociation (CID) experiments, two O(2) units constitute structural elements of Al(2)O(7)(+). Further insight is provided by DFT calculations, performed at the unrestricted B3LYP/TZVP level, and reaction mechanisms are suggested on the basis of both the experimental and theoretical results.     

Electronic and vibrational spectroscopy of intermediates in methane-to-methanol conversion by CoO+

At room temperature, cobalt oxide cations directly convert methane to methanol with high selectivity but very low efficiency. Two potential intermediates of this reaction, the [HO-Co-CH(3)](+) insertion intermediate and [H(2)O-Co=CH(2)](+) aquo-carbene complex are produced in a laser ablation source and characterized by electronic and vibrational spectroscopy. Reaction of laser-ablated cobalt cations with different organic precursors seeded in a carrier gas produces the intermediates, which subsequently expand into vacuum and cool. Ions are extracted into a time-of-flight mass spectrometer and spectra are measured via photofragment spectroscopy. Photodissociation of [HO-Co-CH(3)](+) in the visible and via infrared multiple photon dissociation (IRMPD) makes only Co(+) + CH(3)OH, while photodissociation of [H(2)O-Co=CH(2)](+) produces CoCH(2)(+) + H(2)O. The electronic spectrum of [HO-Co-CH(3)](+) shows progressions in the excited state Co-C stretch (335 cm(-1)) and O-Co-C bend (90 cm(-1)); the IRMPD spectrum gives ν(OH) = 3630 cm(-1). The [HO-Co-CH(3)](+)(Ar) complex has been synthesized and its vibrational spectrum measured in the O-H stretching region. The resulting spectrum is sharper than that obtained via IRMPD and gives ν(OH) = 3642 cm(-1). Also, an improved potential energy surface for the reaction of CoO(+) with methane has been developed using single point energies calculated by the CBS-QB3 method for reactants, intermediates, transition states and products.     

Direct dynamics classical trajectory simulations of the O+ + CH4 reaction at hyperthermal energies

A Born-Oppenheimer direct dynamics simulation of the O(+) + CH(4) reaction dynamics at hyperthermal energies has been carried out with the PM3 (ground quartet state) Hamiltonian. Calculations were performed at various collision energies ranging from 0.5 to 10 eV with emphasis on high energy collisions where this reaction is relevant to materials erosion studies in low Earth orbit and geosynchronous Earth orbit. Charge transfer to give CH(4)(+) is the dominant channel arising from O(+) + CH(4) collisions in this energy range, but most of the emphasis in our study is on collisions that lead to reaction. All energetically accessible reaction channels were found, including products containing carbon-oxygen bonds, which is in agreement with the results of recent experiments. After correcting for compensating errors in competing reaction channels, our excitation functions show quantitative agreement with experiment (for which absolute magnitudes of cross sections are available) at high collision energies (several eV). More detailed properties, such as translational and angular distributions, show qualitative agreement. The opacity function reveals a high selectivity for producing OH(+) at high impact parameters, CH(3)(+)/CH(2)(+)/H(2)O(+) at intermediate impact parameters, and H(2)CO(+)/HCO(+)/CO(+) at small impact parameters. Angular distributions for CH(3)(+)/CH(2)(+)/OH(+) are forward scattered at high collision energies which implies the importance of direct reaction mechanisms, while reaction complexes play an important role at lower energies, especially for the H(2)O(+) product. Finally, we find that the nominally spin-forbidden product CH(3)(+) + OH can be produced by a spin-allowed pathway that involves the formation of the triplet excited product CH(3)(+)(?(3)E). This explains why CH(3)(+) can have a high cross section, even at very low collision energies. The results of this work suggest that the PM3 method may be applied directly to the study of O(+) reactions with small alkane molecules and polymer surfaces.     

The effects of collision energy, vibrational mode, and vibrational angular momentum on energy transfer and dissociation in NO2+-rare gas collisions: an experimental and trajectory study

A combined experimental and trajectory study of vibrationally state-selected NO2+ collisions with Ne, Ar, Kr, and Xe is presented. Ne, Ar, and Kr are similar in that only dissociation to the excited singlet oxygen channel is observed; however, the appearance energies vary by approximately 4 eV between the three rare gases, and the variation is nonmonotonic in rare gas mass. Xe behaves quite differently, allowing efficient access to the ground triplet state dissociation channel. For all four rare gases there are strong effects of NO2+ vibrational excitation that extend over the entire collision energy range, implying that vibration influences the efficiency of collision to internal energy conversion. Bending excitation is more efficient than stretching; however, bending angular momentum partially counters the enhancement. Direct dynamics trajectories for NO2+ + Kr reproduce both the collision energy and vibrational state effects observed experimentally and reveal that intracomplex charge transfer is critical for the efficient energy transfer needed to drive dissociation. The strong vibrational effects can be rationalized in terms of bending, and to a lesser extent, stretching distortion enhancing transition to the Kr+ -NO2 charge state.     

Understanding the complex dissociation dynamics of energy selected dichloroethylene ions: neutral isomerization energies and heats of formation by imaging photoelectron-photoion coincidence

The dissociative photoionization of 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2) has been investigated at high energy and mass resolution using the imaging photoelectron photoion coincidence instrument at the Swiss Light Source. The asymmetric Cl-atom loss ion time-of-flight distributions were fitted to obtain the dissociation rates in the 10(3) s(-1) < k < 10(7) s(-1) range as a function of the ion internal energy. The results, supported by ab initio calculations, show that all three ions dissociate to the same C(2v) symmetry ClC═CH(2)(+) product ion. The 0 K onset energies thus establish the relative heats of formation of the neutral isomers, that is, the isomerization energies. The experimental rate constants, k(E), as well as ab initio calculations indicate an early isomerization transition state and no overall reverse barrier to dissociation. The major high energy channels are the parallel HCl loss and the sequential ClC═CH(2)(+) → HCCH(+) + Cl process, the latter in competition with a ClC═CH(2)(+) → ClCCH(+) + H reaction. A parallel C(2)H(2)Cl(2)(+) → C(2)HCl(2)(+) + H channel also weakly asserts itself. The 0 K onset energy for the sequential Cl loss reaction suggests no barrier to the production of the most stable acetylene ion product; thus the sequential Cl-atom loss is preceded by a ClC═CH(2)(+) → HC(Cl)CH(+) reorganization step with a barrier lower than that of the second Cl-atom loss. The breakdown diagram corresponding to this sequential dissociation reveals the internal energy distribution of the first C(2)H(2)Cl(+) daughter ion, which is determined by the kinetic energy release in the first, Cl loss reaction at high excess energies. At low kinetic energy release, this distribution corresponds to the predicted two translational degrees of freedom, whereas at higher energies, the excess energy partitioning is characteristic of only one translational degree of freedom. New Δ(f)H(o)(298K) of 3.7, 2.5, and 0.2 ± 1.75 kJ mol(-1) are proposed for 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2), respectively, and the proton affinity of ClCCH is found to be 708.6 ± 2.5 kJ mol(-1).     

Theoretical investigation of rare gas hydride cations: HRgN2+ (Rg=He, Ar, Kr, and Xe)

Rare gas containing protonated nitrogen cations, HRgN(2)(+) (Rg=He, Ar, Kr, and Xe), have been predicted using quantum computational methods. HRgN(2)(+) ions exhibit linear structure (C(∞v) symmetry) at the minima and show planar structure (C(s) symmetry) at the transition state. The stability is determined by computing the energy differences between the predicted ions and its various unimolecular dissociation products. Analysis of energy diagram indicates that HXeN(2)(+) is thermodynamically stable with respect to dissociated products while HHeN(2)(+), HArN(2)(+), and HKrN(2)(+) ions are metastable with small barrier heights. Moreover, the computed intrinsic reaction coordinate analysis also confirms that the minima and the 2-body global dissociation products are connected through transition states for the metastable ions. The coupled-cluster theory computed dissociation energies corresponding to the 2-body dissociation (HN(2)(+) + Rg) is -288.4, -98.3, -21.5, and 41.4 kJ mol(-1) for HHeN(2)(+), HArN(2)(+), HKrN(2)(+), and HXeN(2)(+) ions, respectively. The dissociation energies are positive for all the other channels implying that the predicted ions are stable with respect to other 2- and 3-body dissociation channels. Atoms-in-molecules analysis indicates that predicted ions may be best described as HRg(+)N(2). It should be noted that the energetic of HXeN(2)(+) ion is comparable to that of the experimentally observed stable mixed cations, viz. (RgHRg')(+). Therefore, it may be possible to prepare and characterize HXeN(2)(+) ions in an electron bombardment matrix isolation technique.     

Overtone-induced dissociation and isomerization dynamics of the hydroxymethyl radical (CH2OH and CD2OH). I. A theoretical study

The dissociation of the hydroxymethyl radical, CH(2)OH, and its isotopolog, CD(2)OH, following the excitation of high OH stretch overtones is studied by quasi-classical molecular dynamics calculations using a global potential energy surface (PES) fitted to ab initio calculations. The PES includes CH(2)OH and CH(3)O minima, dissociation products, and all relevant barriers. Its analysis shows that the transition states for OH bond fission and isomerization are both very close in energy to the excited vibrational levels reached in recent experiments and involve significant geometry changes relative to the CH(2)OH equilibrium structure. The energies of key stationary points are refined using high-level electronic structure calculations. Vibrational energies and wavefunctions are computed by coupled anharmonic vibrational calculations. They show that high OH-stretch overtones are mixed with other modes. Consequently, trajectory calculations carried out at energies about ~3000 cm(-1) above the barriers reveal that despite initial excitation of the OH stretch, the direct OH bond fission is relatively slow (10 ps) and a considerable fraction of the radicals undergoes isomerization to the methoxy radical. The computed dissociation energies are: D(0)(CH(2)OH → CH(2)O + H) = 10,188 cm(-1), D(0)(CD(2)OH → CD(2)O + H) = 10,167 cm(-1), D(0)(CD(2)OH → CHDO + D) = 10,787 cm(-1). All are in excellent agreement with the experimental results. For CH(2)OH, the barriers for the direct OH bond fission and isomerization are: 14,205 and 13,839 cm(-1), respectively.     

Activation of the C-I and C-OH bonds of 2-iodoethanol by gas phase silver cluster cations yields subvalent silver-iodide and -hydroxide cluster cations

The gas phase ion-molecule reactions of silver cluster cations (Ag(n)(+)) and silver hydride cluster cations (Ag(m)H(+)) with 2-iodoethanol have been examined using multistage mass spectrometry experiments in a quadrupole ion trap mass spectrometer. These clusters exhibit size selective reactivity: Ag(2)H(+), Ag(3)(+), and Ag(4)H(+) undergo sequential ligand addition only, while Ag(5)(+) and Ag(6)H(+) also promote both C-I and C-OH bond activation of 2-iodoethanol. Collision induced dissociation (CID) of Ag(5)HIO(+), the product of C-I and C-OH bond activation by Ag(5)(+), yielded Ag(4)OH(+), Ag(4)I(+) and Ag(3)(+), consistent with a structure containing AgI and AgOH moieties. Ag(6)H(+) promotes both C-I and C-OH bond activation of 2-iodoethanol to yield the metathesis product Ag(6)I(+) as well as Ag(6)H(2)IO(+). The metathesis product Ag(6)I(+) also promotes C-I and C-OH bond activation.DFT calculations were carried out to gain insights into the reaction of Ag(5)(+) with ICH(2)CH(2)OH by calculating possible structures and their energies for the following species: (i) initial adducts of Ag(5)(+) and ICH(2)CH(2)OH, (ii) the subsequent Ag(5)HIO(+) product, (iii) CID products of Ag(5)HIO(+). Potential adducts were probed by allowing ICH(2)CH(2)OH to bind in different ways (monodentate through I, monodentate through OH, bidentate) at different sites for two isomers of Ag(5)(+): the global minimum "bowtie" structure, 1, and the higher energy trigonal bipyramidal isomer, 2. The following structural trends emerged: (i) ICH(2)CH(2)OH binds in a monodentate fashion to the silver core with little distortion, (ii) ICH(2)CH(2)OH binds to 1 in a bidentate fashion with some distortion to the silver core, and (iii) ICH(2)CH(2)OH binds to 2 and results in a significant distortion or rearrangement of the silver core. The DFT calculated minimum energy structure of Ag(5)HIO(+) consists of an OH ligated to the face of a distorted trigonal bipyramid with I located at a vertex, while those for both Ag(4)X(+) (X = OH, I) involve AgX bound to a Ag(3)(+) core. The calculations also predict the following: (i) the ion-molecule reaction of Ag(5)(+) and ICH(2)CH(2)OH to yield Ag(5)HIO(+) is exothermic by 34.3 kcal mol(-1), consistent with the fact that this reaction readily occurs under the near thermal experimental conditions, (ii) the lowest energy products for fragmentation of Ag(5)HIO(+) arise from loss of AgI, consistent with this being the major pathway in the CID experiments.     

Effect of the Ar-Ni(s) potential on the cross section for Ar+CH4/Ni{111} collision-induced desorption and the need for a more accurate CH4/Ni{111} potential

In a previous paper [L. Sun, P. de Sainte Claire, O. Meroueh, and W. L Hase, J. Chem. Phys. 114, 535 (2001)], a classical trajectory simulation was reported of CH(4) desorption from Ni{111} by Ar-atom collisions. At an incident angle theta(i) of 60 degrees (with respect to the surface normal), the calculated collision-induced desorption (CID) cross sections are in excellent agreement with experiment. However, for smaller incident angles the calculated cross sections are larger than the experimental values and for normal collisions, theta(i)=0 degrees , the calculated cross sections are approximately a factor of 2 larger. This trajectory study used an analytic function for the Ar+Ni(s) intermolecular potential which gives an Ar-Ni{111} potential energy minimum which is an order of magnitude too deep. In the work reported here, the previous trajectory study is repeated with an Ar+Ni(s) analytic intermolecular potential which gives an accurate Ar-Ni{111} potential energy minimum and also has a different surface corrugation than the previous potential. Though there are significant differences between the two Ar+Ni(s) analytic potentials, they have no important effects on the CID dynamics and the cross sections reported here are nearly identical to the previous values. Zero-point energy motions of the surface and the CH(4)-Ni(s) intermolecular modes are considered in the simulation and they are found to have a negligible effect on the CID cross sections. Calculations of the intermolecular potential between CH(4) and a Ni atom, at various levels of theory, suggest that there are substantial approximations in the ab initio calculation used to develop the CH(4)+Ni{111} potential. The implication is that the differences between the trajectory and experimental CID cross sections may arise from an inaccurate CH(4)+Ni{111} potential used in the trajectory simulation.     

Dissociative photoionization of 1,3-butadiene: experimental and theoretical insights

The vacuum-ultraviolet photoionization and dissociative photoionization of 1,3-butadiene in a region ～8.5-17 eV have been investigated with time-of-flight photoionization mass spectrometry using tunable synchrotron radiation. The adiabatic ionization energy of 1,3-butadiene and appearance energies for its fragment ions, C(4)H(5)(+), C(4)H(4)(+), C(4)H(3)(+), C(3)H(3)(+), C(2)H(4)(+), C(2)H(3)(+), and C(2)H(2)(+), are determined to be 9.09, 11.72, 13.11, 15.20, 11.50, 12.44, 15.15, and 15.14 eV, respectively, by measurements of photoionization efficiency spectra. Ab initio molecular orbital calculations have been performed to investigate the reaction mechanism of dissociative photoionization of 1,3-butadiene. On the basis of experimental and theoretical results, seven dissociative photoionization channels are proposed: C(4)H(5)(+) + H, C(4)H(4)(+) + H(2), C(4)H(3)(+) + H(2) + H, C(3)H(3)(+) + CH(3), C(2)H(4)(+) + C(2)H(2), C(2)H(3)(+) + C(2)H(2) + H, and C(2)H(2)(+) + C(2)H(2) + H(2). Channel C(3)H(3)(+) + CH(3) is found to be the dominant one, followed by C(4)H(5)(+) + H and C(2)H(4)(+) + C(2)H(2). The majority of these channels occur via isomerization prior to dissociation. Transition structures and intermediates for those isomerization processes were also determined.     

The OH* + CH3SH reaction: support for an addition-elimination mechanism from ab initio calculations

Several intermediates for the CH(3)SH + OH(*) --> CH(3)S(*) + H(2)O reaction were identified using MP2(full) 6-311+g(2df,p) ab initio calculations. An adduct, CH(3)S(H)OH(*), I, with electronic energy 13.63 kJ mol(-1) lower than the reactants, and a transition state, II(double dagger), located 5.14 kJ mol(-1) above I, are identified as the entrance channel for an addition-elimination reaction mechanism. After adding zero-point and thermal energies, DeltaH(r,298) ( degrees )(reactants --> I) = -4.85 kJ mol(-1) and DeltaH(298) (double dagger)(I --> II(double dagger)) = +0.10 kJ mol(-1), which indicates that the potential energy surface is broad and flat near the transition state. The calculated imaginary vibrational frequency of the transition state, 62i cm(-1), is also consistent with an addition-elimination mechanism. These calculations are consistent with experimental observations of the OH(*) + CH(3)SH reaction that favored an addition-elimination mechanism rather than direct hydrogen atom abstraction. An alternative reaction, CH(3)SH + OH(*) --> CH(3)SOH + H(*), with DeltaH(r,298) ( degrees ) = +56.94 kJ mol(-1) was also studied, leading to a determination of DeltaH(f,298) ( degrees )(CH(3)SOH) = -149.8 kJ mol(-1).     

State-selected dynamics of the complex-forming bimolecular reaction Cl- +CH3 Cl'-->ClCH3+Cl'-: a four-dimensional quantum scattering study

Time-independent quantum scattering calculations have been carried out on the Walden inversion S(N)2 reaction Cl(-)+CH(3)Cl(')(v(1),v(2),v(3))-->ClCH(3)(v(1) ('),v(2) ('),v(3) ('))+Cl('-). The two C-Cl stretching modes (quantum numbers v(3) and v(3) (')) and the totally symmetric internal modes of the methyl group (C-H stretching vibration, v(1) and v(1) ('), and inversion bending vibration, v(2) and v(2) (')) are treated explicitly. A four-dimensional coupled cluster potential energy surface is employed. The scattering problem is formulated in hyperspherical coordinates using the exact Hamiltonian and exploiting the full symmetry of the problem. Converged state-selected reaction probabilities and product distributions have been calculated up to 6100 cm(-1) above the vibrational ground state of CH(3)Cl, i.e., up to initial vibrational excitation (2,0,0). In order to extract all scattering resonances, the energetic grid was chosen to be very fine, partly down to a resolution of 10(-12) cm(-1). Up to 2500 cm(-1) translational energy, initial excitation of the umbrella bending vibration, (0,1,0), is more efficient for reaction than exciting the C-Cl stretching mode, (0,0,1). The combined excitation of both vibrations results in a synergic effect, i.e., a considerably higher reaction probability than expected from the sum of both independent excitations, even higher than (0,0,2) up to 1500 cm(-1) translational energy. Product distributions show that the umbrella mode is strongly coupled to the C-Cl stretching mode and cannot be treated as a spectator mode. The reaction probability rises almost linearly with increasing initial excitation of the umbrella bending mode. The effect with respect to the C-Cl stretch is five times larger for more than two quanta in this mode, and in agreement with previous work saturation is found. Exciting the high-frequency C-H stretching mode, (1,0,0), yields a large increase for small energies [more than two orders of magnitude larger than (0,0,0)], while for translational energies higher than 2000 cm(-1), it becomes a pure spectator mode. For combined initial excitations including the symmetric C-H stretch, the spectator character of the latter is even more pronounced. However, up to more than 1500 cm(-1) translational energy, the C-H vibration does not behave adiabatically during the course of reaction, because only 20% of the initial energy is found in the same mode of the product molecule. The distribution of resonance widths and peak heights is discussed, and it is found that individual resonances pertinent to intermediate complexes Cl(-)...CH(3)Cl show product distributions independent of the initial vibrational state of the reactant molecule. The relatively high reactivity, of resonance states with respect to excitation of any mode, found in previous work is confirmed in the present calculations. However, reactivity of intermediate states and reactivity with respect to initial vibrational excitation have to be distinguished. There is a strong mixing between the vibrational states reflected in numerous avoided crossings of the hyperspherical adiabatic curves.