Time-dependent density functional study of the electronic potential energy curves and excitation spectrum of the oxygen molecule

Orbital energies, ionization potentials, molecular constants, potential energy curves, and the excitation spectrum of O(2) are calculated using time-dependent density functional theory (TDDFT) with Tamm-Dancoff approximation (TDA). The calculated negative highest occupied molecular orbital energy (-epsilon(HOMO)) is compared with the energy difference ionization potential for five exchange correlation functionals consisting of the local density approximation (LDAxc), gradient corrected Becke exchange plus Perdew correlation (B(88X)+P(86C)), gradient regulated asymptotic correction (GRAC), statistical average of orbital potentials (SAOP), and van Leeuwen and Baerends asymptotically correct potential (LB94). The potential energy curves calculated using TDDFT with the TDA at internuclear distances from 1.0 to 1.8 A are divided into three groups according to the electron configurations. The 1pi(u) (4)1pi(g) (2) electron configuration gives rise to the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) states; the 1pi(u) (3)1pi(g) (3) electron configuration gives rise to the c (1)Sigma(u) (-), C (3)Delta(u), and A (3)Sigma(u) (+) states; and the B (3)Sigma(u) (-), A (1)Delta(u), and f (1)Sigma(u) (+) states are determined by the mixing of two or more electron configurations. The excitation spectrum of the oxygen molecule, calculated with the aforementioned exchange correlation functionals, shows that the results are quite sensitive to the choice of functional. The LDAxc and the B(88X)+P(86C) functionals produce similar spectroscopic patterns with a single strongly absorbing band positioned at 19.82 and 19.72 eV, respectively, while the asymptotically corrected exchange correlation functionals of the SAOP and the LB94 varieties yield similar excitation spectra where the computed strongly absorbing band is located at 16.09 and 16.42 eV, respectively. However, all of the exchange correlation functionals yield only one strongly absorbing band (oscillator strength greater than 0.1) in the energy interval of 0-20 eV, which is assigned to a X (3)Sigma(g) (-) to (3)Sigma(u) (-) transition. Furthermore, the oxygen molecule has a rich spectrum in the energy range of 14-20 eV and no spin allowed absorption bands are predicted to be observed in the range of 0-6 eV.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.     

Photodissociation channels for N2O near 130 nm studied by product imaging

The photodissociation of N(2)O at wavelengths near 130 nm has been investigated by velocity-mapped product imaging. In all, five dissociation channels have been detected, leading to the following products: O((1)S)+N(2)(X (1)Sigma), N((2)D)+NO(X (2)Pi), N((2)P)+NO(X (2)Pi), O((3)P) + N(2)(A (3)Sigma(+) (u)), and O((3)P) + N(2)(B (3)Pi(g)). The most significant channel is to the products O((1)S) + N(2)(X(1)Sigma), with strong vibrational excitation in the N(2). The O((3)P) + N(2)(A,B):N((2)D,(2)P) + NO branching ratio is measured to be 1.4 +/- 0.5, while the N(2)(A) + O((3)P(J)):N(2)(B) + O((3)P(J)) branching ratio is determined to be 0.84+/-0.09. The spin-orbit distributions for the O((3)P(J)), N((2)P(J)), and N((2)D(J)) products were also determined. The angular distributions of the products are in qualitative agreement with excitation to the N(2)O(D (1)Sigma(+)) state, with participation as well by the (3)Pi(v) state.     

Cluster-enhanced X-O2 photochemistry (X=CH3I, C3H6, C6H12, and Xe)

The effect of a local environment on the photodissociation of molecular oxygen is investigated in the van der Waals complex X-O(2) (X=CH(3)I, C(3)H(6), C(6)H(12), and Xe). A single laser operating at wavelengths around 226 nm is used for both photodissociation of the van der Waals complex and simultaneous detection of the O((3)P(J),J=2,1,0) atom photoproduct via (2+1) resonance enhanced multiphoton ionization. The kinetic energy distribution (KED) and angular anisotropy of the product O atom recoil in this dissociation are measured using the velocity map imaging technique configured for either full ("crush") or partial ("slice") detection of the three-dimensional O((3)P(J)) atom product Newton sphere. The measured KED and angular anisotropy reveal a distinct difference in the mechanism of O atom generation from an X-O(2) complex compared to a free O(2) molecule. The authors identify two one-photon excitation pathways, the relative importance of which depends on IPx, the ionization potential of the X partner. One pathway, observed for all complexes independent of IPx, involves a direct transition to the perturbed covalent state X-O(2)(A'(3)Delta(u)) with excitation localized on the O(2) subunit. The predominantly perpendicular character of this channel relative to the laser polarization detection, together with data on the structure of the complex, allows us to confirm that X partner induced admixing of an X(+)-O(2) (-) charge transfer (CT) state is the perturbing factor resulting in the well-known enhancement of photoabsorption within the Herzberg continuum of molecular oxygen. The second excitation pathway, observed for X-O(2) complexes with X=CH(3)I and C(3)H(6), involves direct excitation into the (3)(X(+)-O(2) (-)) CT state of the complex. The subsequent photodissociation of this CT state by the same laser pulse gives rise to the superoxide anion O(2) (-), which then photodissociates, providing fast (0.69 eV) O atoms with a parallel image pattern. Products from the photodissociation of singlet oxygen O(2)(b (1)Sigma(g) (+)) are also observed when the CH(3)I-O(2) complex was irradiated. Potential energy surfaces (PES) for the ground and relevant excited states of the X-O(2) complex have been constructed for CH(3)I-O(2) using the results of CASSCF calculations for the ground and CT states of the complex as well as literature data on PES of the subunits. These model potential energy surfaces allowed us to interpret all of the observed O((3)P(J)) atom production channels.     

Characterization of singlet ground and low-lying electronic excited states of phosphaethyne and isophosphaethyne

The singlet ground ((approximate)X(1)Sigma1+) and excited (1Sigma-,1Delta) states of HCP and HPC have been systematically investigated using ab initio molecular electronic structure theory. For the ground state, geometries of the two linear stationary points have been optimized and physical properties have been predicted utilizing restricted self-consistent field theory, coupled cluster theory with single and double excitations (CCSD), CCSD with perturbative triple corrections [CCSD(T)], and CCSD with partial iterative triple excitations (CCSDT-3 and CC3). Physical properties computed for the global minimum ((approximate)X(1)Sigma+HCP) include harmonic vibrational frequencies with the cc-pV5Z CCSD(T) method of omega1=3344 cm(-1), omega2=689 cm(-1), and omega3=1298 cm(-1). Linear HPC, a stationary point of Hessian index 2, is predicted to lie 75.2 kcal mol(-1) above the global minimum HCP. The dissociation energy D0[HCP((approximate)X(1)Sigma+)-->H(2S)+CP(X2Sigma+)] of HCP is predicted to be 119.0 kcal mol(-1), which is very close to the experimental lower limit of 119.1 kcal mol(-1). Eight singlet excited states were examined and their physical properties were determined employing three equation-of-motion coupled cluster methods (EOM-CCSD, EOM-CCSDT-3, and EOM-CC3). Four stationary points were located on the lowest-lying excited state potential energy surface, 1Sigma- -->1A", with excitation energies Te of 101.4 kcal mol(-1) (1A"HCP), 104.6 kcal mol(-1)(1Sigma-HCP), 122.3 kcal mol(-1)(1A" HPC), and 171.6 kcal mol(-1)(1Sigma-HPC) at the cc-pVQZ EOM-CCSDT-3 level of theory. The physical properties of the 1A" state with a predicted bond angle of 129.5 degrees compare well with the experimentally reported first singlet state ((approximate)A1A"). The excitation energy predicted for this excitation is T0=99.4 kcal mol(-1) (34 800 cm(-1),4.31 eV), in essentially perfect agreement with the experimental value of T0=99.3 kcal mol(-1)(34 746 cm(-1),4.308 eV). For the second lowest-lying excited singlet surface, 1Delta-->1A', four stationary points were found with Te values of 111.2 kcal mol(-1) (2(1)A' HCP), 112.4 kcal mol(-1) (1Delta HPC), 125.6 kcal mol(-1)(2(1)A' HCP), and 177.8 kcal mol(-1)(1Delta HPC). The predicted CP bond length and frequencies of the 2(1)A' state with a bond angle of 89.8 degrees (1.707 A, 666 and 979 cm(-1)) compare reasonably well with those for the experimentally reported (approximate)C(1)A' state (1.69 A, 615 and 969 cm(-1)). However, the excitation energy and bond angle do not agree well: theoretical values of 108.7 kcal mol(-1) and 89.8 degrees versus experimental values of 115.1 kcal mol(-1) and 113 degrees. of 115.1 kcal mol(-1) and 113 degrees.     

Ab initio study on the decomposition of first excited state HOOO radicals

The HOOO radical plays a crucial role in atmospheric processes involving the OH radical and O(2) molecule. We present an ab initio molecular orbital theory study on the decomposition reaction of the first excited state HOOO((2)A') with respect to OH and O(2). The geometries and harmonic vibrational frequencies of all stationary points are calculated at the CASSCF and MRCI levels of theory in conjunction with the 6-31+G(d,p) basis set. The potential energy profile of the decomposition reaction is studied at the CASSCF/6-31+G(d,p) level of theory, in which the complete valence orbitals and electrons are included in the active space. The energies of the potential energy profile are further refined at the CASPT2 and MRCI levels of the theory. Additionally, we have determined the interesting reaction process: the HOOO((2)A') radical with C(s) symmetry does not dissociate to OH((2)Pi) and O(2)((3)Sigma(-)(g)) directly as this is forbidden by orbital symmetry, but dissociates to OH((2)Pi) and O(2)((3)Sigma(-)(g)) via the change in symmetry from C(s) to C(infinity v) symmetry with a low barrier.     

Photodissociation dynamics of N2O at 130 nm: the N2(A 3Sigmau +,B 3Pig)+O(3PJ = 2,1,0) channels

Oxygen Rydberg time-of-flight spectroscopy was used to study the vacuum ultraviolet photodissociation dynamics of N(2)O near 130 nm. The O((3)P(J)) products were tagged by excitation to high-n Rydberg levels and subsequently field ionized at a detector. In agreement with previous work, we find that O((3)P(J)) formation following excitation to the repulsive N(2)O D((1)Sigma(+)) state produces the first two electronically excited states of the N(2) counterfragment, N(2)(A (3)Sigma(u) (+)) and N(2)(B (3)Pi(g)). The O((3)P(J)) translational energy distribution reveals that the overall branching ratio between N(2)(A (3)Sigma(u) (+)) and N(2)(B (3)Pi(g)) formation is approximately 1.0:1.0 for J = 1 and 2, with slightly less N(2)(B (3)Pi(g)) produced in coincidence with O((3)P(0)). The angular distributions were found to be independent of J and highly anisotropic, with beta = 1.5+/-0.2.     

Coupled-states statistical investigation of vibrational and rotational relaxation of OH(2pi) by collisions with atomic hydrogen

We report state-to-state cross sections and thermal rate constants for vibrational and rotational relaxation of OH(2pi) by collision with H atoms. The cross sections are calculated by the coupled-states (CS) statistical method including the full open-shell character of the OH + H system. Four potential energy surfaces (PESs) ((1,3)A' and (1,3)A') describe the interaction of OH(X2pi) with H atoms. Of these, three are repulsive, and one (1A') correlates with the deep H2O well. Consequently, rotationally and ro-vibrationally inelastic scattering of OH in collisions with H can occur by scattering on the repulsive PESs, in a manner similar to the inelastic scattering of OH by noble gas atoms, or by collisions which enter the H2O well and then reemerge. At 300 K, we predict large (approximately 1 x 10(-10) cm3 molecule(-1) s(-1)) vibrational relaxation rates out of both v = 2 and v = 1, comparable to earlier experimental observations. This anomalously fast relaxation results from capture into the H2O complex. There exists a significant propensity toward formation of OH in the pi(A') lambda-doublet level. We also report state-resolved cross sections and rate constants for rotational excitation within the OH v = 0 manifold. Collisional excitation from the F1 to the F2 spin-orbit manifold leads to an inverted lambda-doublet population.     

Mechanisms of the reactions of W AND W+ with H2O: computational studies

The mechanism of the reactions of W and W(+) with the water molecule have been studied for several lower-lying electronic states of tungsten centers at the CCSD(T)/6-311G(d,p)+SDD and B3LYP/6-31G(d,p)+SDD levels of theory. It is shown that these reactions are essentially multistate processes, during which lower-lying electronic states of the systems cross several times. They start with the formation of initial prereaction M(H(2)O) complexes with M-H(2)O bonding energies of 9.6 and 48.2 kcal/mol for M = W and W(+), followed by insertion of the metal center into an O-H bond with 20.0 and 53.3 kcal/mol barriers for neutral and cationic systems, respectively. The overall process of M + H(2)O --> t-HM(OH) is calculated to be highly exothermic, 48.4 and 48.8 kcal/mol for M = W and W(+). From the HM(OH) intermediate the reaction may proceed via several different channels, among which the stepwise HM(OH) --> HMO + H --> (H)(2)MO and concerted HM(OH) --> (H)(2)MO pathways are more favorable and can compete (energetically) with each other. For the neutral system (M = W), the concerted process is the most favorable, whereas for the charged system (M = W(+)), the stepwise pathway is slightly more favorable. From the energetically most favorable intermediate (H)(2)MO the reactions proceed via H(2)-molecule formation with a 53.1 kcal/mol activation barrier for the neutral system. For the cationic system, H-H formation and dissociation is an almost barrierless process. The overall reaction of W and W(+) with the water molecule leading to H(2) + MO formation is found to be exothermic by 48.2 and 39.8 kcal/mol, respectively. In the gas phase with the collision-less conditions the reactions W((7)S) + H(2)O --> H(2) + WO((3)Sigma(+)), and W(+)((6)D) + H(2)O --> H(2) + WO(+)((4)Sigma(+)) are expected to proceed via a 10.4 and 5.1 kcal/mol overall energy barrier corresponding to the first O-H dissociation at the TS1. On the basis of these PESs, we predict kinetic rate constants for the reactions of W and W(+) with H(2)O.     

Ab initio potential energy surfaces and nonadiabatic collision dynamics in H(+)+O(2) system

The adiabatic potential energy surfaces for the lowest five electronic states of (3)A" symmetry for the H(+)+O(2) collision system have been obtained at the multireference configuration interaction level of accuracy using Dunning's correlation consistent polarized valence triple zeta basis set. The radial nonadiabatic coupling terms and the mixing angle between the lowest two electronic states (1 (3)A" and 2 (3)A"), which adiabatically correlate in the asymptotic limit to H((2)S)+O(2) (+)(X (2)Pi(g)) and H(+)+O(2)(X (3)Sigma(g)(-)), respectively, have been computed using ab initio procedures at the same level of accuracy to yield the corresponding quasidiabatic potential energy matrix. The computed strengths of the vibrational coupling matrix elements reflect the trend observed for inelastic vibrational excitations of O(2) in the experiments at collision energy of 9.5 eV. The quantum dynamics has been preformed on the newly obtained coupled quasidiabatic potential energy surfaces under the vibrational close-coupling rotational infinite-order sudden framework at the experimental collision energy of 9.5 eV. The present theoretical results for vibrational elastic/inelastic excitations of O(2) are in overall good agreement with the available experimental data obtained from the proton energy-loss spectra in molecular beam experiments [F. A. Gianturco et al., J. Phys. B 14, 667 (1981)]. The results for the complementary charge transfer processes are also presented at this collision energy.     

An ab initio study of the low-lying electronic states of S3

Accurate calculations of the low-lying singlet and triplet electronic states of thiozone, S(3), have been carried out using large multireference configuration interaction wave functions. Cuts of the full potential energy surfaces along the stretching and bending coordinates have been presented, together with the vertical excitation spectra. The strong experimentally observed absorption around 395 nm is assigned to the 1 (1)B(2) state, which correlates to ground state products. Absorption at wavelengths shorter than 260 nm is predicted to lead to singlet excited state products, S(2) (a (1)Delta(g))+S((1)D). The spectroscopic properties of the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) electronic states of the S(2) radical have also been accurately characterized in this work. The investigations of the low-lying electronic states were accompanied by accurate ground state coupled cluster calculations of the thermochemistry of both S(2) and S(3) using large correlation consistent basis sets with corrections for core-valence correlation, scalar relativity, and atomic spin-orbit effects. Resulting values for D(0)(S(2)+S) and SigmaD(0) for S(3) are predicted to be 61.3 and 162.7 kcal/mol, respectively, with conservative uncertainties of +/-1 kcal/mol. Analogous calculations predict the C(2v)-D(3h) (open-cyclic) isomerization energy of S(3) to be 4.4+/-0.5 kcal/mol.     

Observation and calculation of the Cs2 2 3Delta(1g) and b 3Pi(0u) states

The Cs(2) 2 (3)Delta(1g) and b (3)Pi(0u) states have been observed by infrared-infrared double resonance spectroscopy for the first time. 221 2 (3)Delta(1g)<--A (1)Sigma(u) (+)<--X (1)Sigma(g) (+) double resonance lines have been assigned to transitions into the 2 (3)Delta(1g) v=6-13 vibrational levels. Resolved fluorescence into the b (3)Pi(0u) v(')=0-48 levels has been recorded. Molecular constants and potential energy curves are determined by the global fit of the entire set of the experimental data. Theoretical potential energy curves of the 2 (3)Delta(g) and b (3)Pi(u) states have been determined in the framework of the pseudopotential method and are compared with the experimental results.     

Ab initio potential-energy surfaces of O2(X3Sigmag -, a1Deltag, b1Sigmag +) +O2 (X3Sigmag -, a1Deltag, b1Sigmag +): mechanism of quenching of O2 (a 1Deltag)

Ab initio computational studies were carried out in order to explore the possible mechanisms of quenching of O(2)(a (1)Delta(g)) by O(2)(X (3)Sigma(g) (-)): the self-quenching of O(2)(a (1)Delta(g)) and other energy-transfer processes involving two O(2) molecules. All eighteen states arising from two O(2) molecules in the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) states are considered. After scans at the state-averaged complete active space self-consistent field method to identify possible regions of crossing between states belonging to different asymptotes, complete active state second-order perturbation theory high-symmetry optimization and low-symmetry scans established that four different minima on the seams of crossing (MSXs), arising between the a+a manifold and the X+b manifold and responsible for self-quenching: O(2)(a (1)Delta(g))+O(2)(a (1)Delta(g))-->O(2)(X (3)Sigma(g) (-))+O(2)(b (1)Sigma(g) (+)), have coplanar C(2h) or C(2v) symmetries and are only 0.45 eV barrier relative to the a+a asymptote and energetically easily accessible. The rate constant for this process was estimated based on the Landau-Zener formalism. The MSXs for quenching of O(2)(a (1)Delta(g)) by the ground state O(2)(X (3)Sigma(g) (-)):O(2)(a (1)Delta(g))+O(2)(X (3)Sigma(g) (-))-->O(2)(X (3)Sigma(g) (-))+O(2)(X (3)Sigma(g) (-)) require higher energies and the process is not likely to be important.     

Spectroscopic investigation of Al2N and its anion via negative ion photoelectron spectroscopy

Negative ion photoelectron spectroscopy was used to elucidate the electronic and geometric structure of the gaseous Al2N/Al2N- molecules, using photodetachment wavelengths of 416 nm (2.977 eV), 355 nm (3.493 eV), and 266 nm (4.661 eV). Three electronic bands are observed and assigned to the X2Sigma(u)+ <-- X1Sigma(g)+, A2Pi(u) <-- X1Sigma(g)+, and B2Sigma(g)+ <-- X1Sigma(g)+ electronic transitions, with the caveat that one or both excited states may be slightly bent. With the aid of density functional theory calculations and Franck-Condon spectral simulations, we determine the adiabatic electron affinity of Al2N, 2.571 +/- 0.008 eV, along with geometry changes upon photodetachment, vibrational frequencies, and excited-state term energies. Observation of excitation of the odd vibrational levels of the antisymmetric stretch (nu3) suggests a breakdown of the Franck-Condon approximation, caused by the vibronic coupling between the X2Sigma(u)+ and B2Sigma(g)+ electronic states through the nu3 mode.     

Quantum calculations of the O(3P)+H2-->OH+H reaction

Quantum scattering calculations are reported for the O(3P)+H2(v=0,1) reaction using chemically accurate potential energy surfaces of 3A' and 3A" symmetry. We present state-to-state reaction cross sections and rate coefficients as well as thermal rate coefficients for the title reaction using accurate quantum calculations. Our calculations yield reaction cross sections that are in quantitative accord with results of recent crossed molecular beam experiments. Comparisons with results obtained using the J-shifting calculations show that the J-shifting approximation is quite reliable for this system. Thermal rate coefficients from the exact calculations and the J-shifting approximation agree remarkably well with experimental results. Our calculations also reproduce the markedly different OH(v'=0)/OH(v'=1) branching in O(3P)+H2(v=1) reaction, observed in experiments that use different O(3P) atom sources. In particular, we show that the branching ratio is a strong function of the kinetic energy of the O(3P) atom.     

Negative-ion photoelectron spectroscopy, gas-phase acidity, and thermochemistry of the peroxyl radicals CH(3)OO and CH(3)CH(2)OO

Methyl, methyl-d(3), and ethyl hydroperoxide anions (CH(3)OO(-), CD(3)OO(-), and CH(3)CH(2)OO(-)) have been prepared by deprotonation of their respective hydroperoxides in a stream of helium buffer gas. Photodetachment with 364 nm (3.408 eV) radiation was used to measure the adiabatic electron affinities: EA[CH(3)OO, X(2)A' '] = 1.161 +/- 0.005 eV, EA[CD(3)OO, X(2)A' '] = 1.154 +/- 0.004 eV, and EA[CH(3)CH(2)OO, X(2)A' '] = 1.186 +/- 0.004 eV. The photoelectron spectra yield values for the term energies: Delta E(X(2)A' '-A (2)A')[CH(3)OO] = 0.914 +/- 0.005 eV, Delta E(X(2)A' '-A (2)A')[CD(3)OO] = 0.913 +/- 0.004 eV, and Delta E(X(2)A' '-A (2)A')[CH(3)CH(2)OO] = 0.938 +/- 0.004 eV. A localized RO-O stretching mode was observed near 1100 cm(-1) for the ground state of all three radicals, and low-frequency R-O-O bending modes are also reported. Proton-transfer kinetics of the hydroperoxides have been measured in a tandem flowing afterglow-selected ion flow tube (FA-SIFT) to determine the gas-phase acidity of the parent hydroperoxides: Delta(acid)G(298)(CH(3)OOH) = 367.6 +/- 0.7 kcal mol(-1), Delta(acid)G(298)(CD(3)OOH) = 367.9 +/- 0.9 kcal mol(-1), and Delta(acid)G(298)(CH(3)CH(2)OOH) = 363.9 +/- 2.0 kcal mol(-1). From these acidities we have derived the enthalpies of deprotonation: Delta(acid)H(298)(CH(3)OOH) = 374.6 +/- 1.0 kcal mol(-1), Delta(acid)H(298)(CD(3)OOH) = 374.9 +/- 1.1 kcal mol(-1), and Delta(acid)H(298)(CH(3)CH(2)OOH) = 371.0 +/- 2.2 kcal mol(-1). Use of the negative-ion acidity/EA cycle provides the ROO-H bond enthalpies: DH(298)(CH(3)OO-H) = 87.8 +/- 1.0 kcal mol(-1), DH(298)(CD(3)OO-H) = 87.9 +/- 1.1 kcal mol(-1), and DH(298)(CH(3)CH(2)OO-H) = 84.8 +/- 2.2 kcal mol(-1). We review the thermochemistry of the peroxyl radicals, CH(3)OO and CH(3)CH(2)OO. Using experimental bond enthalpies, DH(298)(ROO-H), and CBS/APNO ab initio electronic structure calculations for the energies of the corresponding hydroperoxides, we derive the heats of formation of the peroxyl radicals. The "electron affinity/acidity/CBS" cycle yields Delta(f)H(298)[CH(3)OO] = 4.8 +/- 1.2 kcal mol(-1) and Delta(f)H(298)[CH(3)CH(2)OO] = -6.8 +/- 2.3 kcal mol(-1).     

Quasi-classical trajectory study of Si+O2-->SiO+O reaction

Quasi-classical trajectory calculations for the Si(3P)+O2(X 3Sigmag-)-->SiO(X 1Sigma+)+O(1D) reaction have been carried out using the analytical ground 1A' potential energy surface (PES) recently reported by Dayou and Spielfiedel [J. Chem. Phys. 119, 4237 (2003)]. The reaction has been studied for a wide range of collision energies (0.005-0.6 eV) with O2 in its ground rovibrational state. The barrierless PES leads to a decrease of the total reaction cross section with increasing collision energy. It has been brought to evidence that the reaction proceeds through different reaction mechanisms whose contributions to reactivity are highly dependent on the collision energy range. At low collision energy an abstraction mechanism occurs involving the collinear SiOO potential well. The associated short-lived intermediate complex leads to an inverted vibrational distribution peaked at v'=3 and low rotational excitation of SiO(v',j') with a preferentially backward scattering. At higher energies the reaction proceeds mainly through an insertion mechanism involving the bent and linear OSiO deep potential wells and associated long-lived intermediate complexes, giving rise to nearly statistical energy disposals into the product modes and a forward-backward symmetry of the differential cross section.     

Direct dynamics simulations of O(3P) + HCl at hyperthermal collision energies

The dynamics of the O(3P) + HCl reaction at hyperthermal collision energies were investigated using the quasiclassical trajectory method. Stationary points on the OClH 3A" and 3A' potential energy surfaces (PESs) were also examined. The lowest transition state leading to OCl + H on the 3A" surface is 2.26 eV above the reagents at the CCSD(T)/cc-pVTZ level of theory. This saddle point is bent and product-like. Direct dynamics calculations at the MP2/cc-pVTZ level of theory were used to investigate the excitation functions for OH + Cl, OCl + H, and O + H + Cl formation. OCl is formed mainly from small-impact-parameter collisions, and the OCl + H excitation function peaks around 5 eV, where it is similar in magnitude to the OH + Cl excitation function. The shape of the OCl + H excitation function is discussed, and features are identified that should be general to hyperthermal collision dynamics.     

Application of mean-field and surface-hopping approaches for interrogation of the Xe3+ molecular ion photoexcitation dynamics

The dissociation dynamics of the excited Xe(3) (+) molecular ion through the Pi(12)(u) and Pi(12)(g) conical intersection was interrogated by computational simulation in which no adjustable parameters were used. The electronic ground and excited state potential energy surfaces were generated by the diatomics-in-molecules method, and the Ehrenfest mean-field and Tully surface-hopping approaches treated the nonadiabatic interactions. Reproduction of the experimental spectrum of the symmetric photofragmentation as a function of excitation energy was obtained within the region of interest (2.5-3.75 eV), with the exception of a 0.25 eV width on the red side of the spectral apex. Good agreement was obtained with the experimental dissociated photofragment kinetic energy spectra. It was determined that the greatest contribution to the nonadiabatic coupling between the two states originated from the bending vibrational mode of the molecule in the Sigma(12)(u), ground electronic state before excitation.