Quasiclassical trajectory calculations comparing the reactivity and dynamics of symmetric and asymmetric stretch and the role of the bending mode excitations of methane in the Cl + CH4 reaction

To analyze the effects of the symmetric (nu(1)) and asymmetric (nu(3)) stretch mode excitations and the role played by the "umbrella" bending (nu(4)) mode excitation in the reactivity and the dynamics of the gas-phase Cl+CH(4) reaction, an exhaustive dynamics study was performed. Quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were used in this work on an analytical potential energy surface previously developed by Espinosa-Garcia et al. [J. Chem. Phys. (to be published)]. First, with respect to the reactivity, we found that the nu(1) mode excitation is more reactive than the nu(3) mode by a factor of 1.20, in agreement with the experimental tendency between these modes. The inclusion of the nu(4) bending mode practically does not affect this relative reactivity, (nu(1+)nu(4))(nu(3+)nu(4)) = 1.16. Second, with respect to the dynamics (rotovibrational and angular distributions of the products), the two stretch modes, nu(1) and nu(3), give very similar pictures, reproducing the experimental behavior, and the nu(4) "umbrella" mode does not affect the dynamics. The satisfactory reproduction (always qualitatively acceptable and sometimes even quantitatively) of a great variety of experimental data by the QCT study presented here lends confidence to the potential energy surface constructed by Espinosa-Garcia et al. [J. Chem. Phys. (to be published)].     

Quasi-classical trajectory calculations analyzing the reactivity and dynamics of asymmetric stretch mode excitations of methane in the H + CH4 reaction

An exhaustive dynamics study was performed at two collision energies, 1.52 and 2.20 eV, analyzing the effects of the asymmetric (nu3) stretch mode excitation in the reactivity and dynamics of the gas-phase H + CH4 reaction. Quasi-classical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were performed on an analytical potential energy surface previously developed by our group. First, strong coupling between different vibrational modes in the entry channel was observed, indicating that energy can flow between these modes, and therefore that they do not preserve their adiabatic character along the reaction path; i.e., the reaction is nonadiabatic. Second, we found that the reactant vibrational excitation has a significant influence on the vibrational and rotational product distributions. With respect to the vibrational distribution, our results confirm the purely qualitative experimental evidence, although the theoretical results presented here are also quantitative. The rotational distributions are predictive, because no experimental data have been reported. Third, with respect to the reactivity, we found that the nu3 mode excitation by one quantum is more reactive than the ground state by a factor of about 2, independently of the collision energy, and in agreement with the experimental measurement of 3.0 +/- 1.5. Fourth, the state-to-state angular distributions of the products reproduce the experimental behavior at 1.52 eV, where the CH3 products scatter sideways and backward. At 2.20 eV this experimental information is not available, and therefore the results reported here are again predictive. The satisfactory reproduction of a great variety of experimental data by the present QCT study lends confidence to the potential energy surface constructed by our group and to those results whose accuracy cannot be checked by comparison with experiment.     

Vibrationally controlled chemistry: mode- and bond-selected reaction of CH3D with Cl

Selective vibrational excitation controls the competition between C-H and C-D bond cleavage in the reaction of CH(3)D with Cl, which forms either HCl + CH(2)D or DCl + CH(3). The reaction of CH(3)D molecules with the first overtone of the C-D stretch (2nu(2)) excited selectively breaks the C-D bond, producing CH(3) exclusively. In contrast, excitation of either the symmetric C-H stretch (nu(1)), the antisymmetric C-H stretch (nu(4)), or a combination of antisymmetric stretch and CH(3) umbrella bend (nu(4) + nu(3)) causes the reaction to cleave only a C-H bond to produce solely CH(2)D. Initial preparation of C-H stretching vibrations with different couplings to the reaction coordinate changes the rate of the H-atom abstraction reaction. Excitation of the symmetric C-H stretch (nu(1)) of CH(3)D accelerates the H-atom abstraction reaction 7 times more than excitation of the antisymmetric C-H stretch (nu(4)) even though the two lie within 80 cm(-1) of the same energy. Ab initio calculations and a simple theoretical model help identify the dynamics behind the observed mode selectivity.     

Effect of bending and torsional mode excitation on the reaction Cl+CH4-->HCl+CH3

A beam containing CH(4), Cl(2), and He is expanded into a vacuum chamber where CH(4) is prepared via infrared excitation in a combination band consisting of one quantum of excitation each in the bending and torsional modes (nu(2)+nu(4)). The reaction is initiated by fast Cl atoms generated by photolysis of Cl(2) at 355 nm, and the resulting CH(3) and HCl products are detected in a state-specific manner using resonance-enhanced multiphoton ionization (REMPI). By comparing the relative amplitudes of the action spectra of Cl+CH(4)(nu(2)+nu(4)) and Cl+CH(4)(nu(3)) reactions, we determine that the nu(2)+nu(4) mode-driven reaction is at least 15% as reactive as the nu(3) (antisymmetric stretch) mode-driven reaction. The REMPI spectrum of the CH(3) products shows no propensity toward the formation of umbrella bend mode excited methyl radical, CH(3)(nu(2)=1), which is in sharp distinction to the theoretical expectation based on adiabatic correlations between CH(4) and CH(3). The rotational distribution of HCl(v=1) products from the Cl+CH(4)(nu(2)+nu(4)) reaction is hotter than the corresponding distribution from the Cl+CH(4)(nu(3)) reaction, even though the total energies of the two reactions are the same within 4%. An explanation for this enhanced rotational excitation of the HCl product from the Cl+CH(4)(nu(2)+nu(4)) reaction is offered in terms of the projection of the bending motion of the CH(4) reagent onto the rotational motion of the HCl product. The angular distributions of the HCl(nu=0) products from the Cl+CH(4)(nu(2)+nu(4)) reaction are backward scattered, which is in qualitative agreement with theoretical calculation. Overall, nonadiabatic product vibrational correlation and mode specificity of the reaction indicate that either the bending mode or the torsional mode or both modes are strongly coupled to the reaction coordinate.     

Role of the C-H stretch mode excitation in the dynamics of the Cl + CHD3 reaction: a quasi-classical trajectory calculation

To analyze the effect of the C-H stretch mode excitation on the dynamics of the Cl + CHD3 gas-phase abstraction reaction, an exhaustive state-to-state dynamics study was performed. This reaction can evolve along two channels: H-abstraction, CD3 + ClH, and D-abstraction, CHD2 + ClD. On an analytical potential energy surface constructed previously by our group, named PES-2005, quasi-classical trajectory calculations were performed at a collision energy of 0.18 eV, including corrections to avoid zero-point energy leakage along the trajectories. First, strong coupling between different vibrational modes in the entry valley was observed; i.e., the reaction is vibrationally nonadiabatic. Second, for the ground-state CHD3(nu=0) reaction, the diatomic fragments appeared in their ground states, and the H- and D-abstraction reactions showed similar reactivities. However, when the reactivity per atom is considered, the H is three times more reactive than the D atom. Third, when the C-H stretch mode is excited by one quantum, CHD3(nu1=1), the H-abstraction is strongly favored, and the C-H stretch excitation is maintained in the product CHD2(nu1=1) + ClD channel; i.e., the reaction shows mode selectivity, reproducing the experimental evidence, and also the reactivity of the vibrational ground state is increased, in agreement with experiment. Fourth, the state-to-state angular distributions of the CD3 and CHD2 products showed the products to be practically sideways for the reactant ground state, while the C-H excitation yielded a more forward scattering, reproducing the experimental data. The role of the zero-point energy correction was also analyzed, and we find that the dynamics results are very sensitive on how the ZPE issue is treated. Finally, a comparison is made with the similar H + CHD3(nu1=0,1) and Cl + CH4(nu1=0,1) reactions.     

State-to-state dynamics analysis of the F + CHD3 reaction: a quasiclassical trajectory study

An exhaustive state-to-state dynamics study was performed to analyze the F + CHD3 --> FD(nu', j') + CHD2(nu) gas-phase abstraction reaction. Quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were performed at different collision energies on an analytical potential energy surface (PES-2006) recently developed by our group. Whereas the CHD2 coproduct appears vibrationally and rotationally cold, most of the available energy appears as FD(nu') product vibrational energy, peaking at nu' = 2 and nu' = 3, with the population in the latter level growing as the energy increases. The excitation function rises from the threshold of the reaction and then levels off at higher energies, with the maximum contribution from the FD(nu' = 3) level. The state-specific FD(nu') scattering distributions correlated with the coproduct CHD2 in the nu4 = 2 and nu3 = 1 states, at different collision energies, show a steady change from backward to forward scattering as the energy increases. This similar behavior for the two coproduct vibrational states, nu4 = 2 and nu3 = 1, agrees qualitatively with the experimental measurements. Comparison with theoretical and experimental results for the isotopic analogues, F + CH4 and F + CD4, shows that the title reaction presents a direct mechanism, similar to the perdeuterated reaction, but contrasts with that of the F + CH4 reaction. These results for the dynamics of different isotopic variants, always in qualitative and sometimes in quantitative agreement with experiment, show the capacity of the PES-2006 surface to correctly describe the title reaction, even though there are differences that could be due to deficiencies of the PES but also to the known limitations of the classical treatment in the QCT method.     

Resonance Raman and density functional theory investigation of the photodissociation dynamics of the A-band absorption of (E)-beta-nitrostyrene in cyclohexane solution

Resonance Raman spectra were obtained for (E)-beta-nitrostyrene in cyclohexane solution with excitation wavelengths in resonance with the charge transfer (CT)-band absorption spectrum. These spectra indicate that the Franck-Condon region photodissociation dynamics have multidimensional character with motion predominantly along the nominal NO(2) symmetric stretch mode (nu(14)), the nominal C=C stretch mode (nu(8)), the nominal benzene ring stretch mode (nu(9)), accompanied by a smaller amount of motion along the nominal ONO symmetric bend/benzene ring stretch mode (nu(34)), the nominal CCH in-plane bending mode (nu(20)), the nominal HC=CH in-plane bending mode (nu(18)), the nominal NO(2) asymmetric stretch mode (nu(11)), the nominal C-N stretch/benzene ring breathing mode (nu(27)), and the nominal CCC trigonal bending mode (nu(25)). A preliminary resonance Raman intensity analysis was done and these results for (E)-beta-nitrostyrene were compared to results previously reported for several nitrobenzene and trans-stilbene compounds. The differences and similarities between the CT-band resonance Raman spectra and vibrational reorganizational energies for (E)-beta-nitrostyrene relative to those for nitrobenzene and trans-stilbene were briefly discussed.     

Infrared absorption spectra of vinyl radicals isolated in solid Ne

Irradiation of samples of solid Ne near 3.0 K containing ethene (C(2)H(4)) with vacuum ultraviolet radiation at 120 nm from synchrotron yielded new spectral lines at 3141.0, 2953.6, 2911.5, 1357.4, 677.1, 895.3, and 857.0 cm(-1). These features are assigned to alpha-CH stretching (nu(1)), CH(2) antisymmetric stretching (nu(2)), CH(2) symmetric stretching (nu(3)), CH(2)-bending (nu(5)), HCCH cis bending (nu(7)), CH(2) out-of-plane bending (nu(8)), and alpha-CH out-of-plane bending (nu(9)) modes of C(2)H(3), respectively, based on results of (13)C- and D-isotopic experiments and quantum-chemical calculations. These calculations using density-functional theory (B3LYP and PW91PW91/aug-cc-pVTZ) predict vibrational wavenumbers, IR intensities, and isotopic ratios of vinyl radical that agree satisfactorily with our experimental results.     

Potential energy surface, kinetics, and dynamics study of the Cl+CH4-->HCl+CH3 reaction

A modified and recalibrated potential energy surface for the gas-phase Cl+CH4-->HCl+CH3 reaction is reported and tested. It is completely symmetric with respect to the permutation of the four methane hydrogen atoms and is calibrated with respect to updated experimental and theoretical stationary point properties and experimental forward thermal rate constants. From the kinetics point of view, the forward and reverse thermal rate constants and the activation energies were calculated using the variational transition-state theory with semiclassical transmission coefficients over a wide temperature range of 150-2500 K. The theoretical results reproduce the available experimental data, with a small curvature of the Arrhenius plot which indicates the role of tunneling in this hydrogen abstraction reaction. A dynamics study was also performed on this PES using quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found a noticeable internal energy in the coproduct methyl radical, both in the ground-state [CH4 (v=0)] and vibrationally excited [CH4 (v=1)] reactions. This CH3 internal energy was directly precluded in some experiments or oversimplified in previous theoretical studies using pseudotriatomic models. Second, our QCT calculations give HCl rotational distributions slightly hotter than those in experiment, but correctly describing the experimental trend of decreasing the HCl product rotation excitation in going from HCl (v'=0) to HCl (v'=1) for the CH4 (v=1) reaction. Third, the state specific scattering distributions present a reasonable agreement with experiment, although they tend to make the reaction more forward and backward scattered than found experimentally probably because of the hotter rotational distribution and the deficiencies of the QCT methods.     

Communication: Experimental and theoretical investigations of the effects of the reactant bending excitations in the F + CHD3 reaction

The effects of the reactant bending excitations in the F+CHD(3) reaction are investigated by crossed molecular beam experiments and quasiclassical trajectory (QCT) calculations using a high-quality ab initio potential energy surface. The collision energy (E(c)) dependence of the cross sections of the F+CHD(3)(v(b)=0,1) reactions for the correlated product pairs HF(v('))+CD(3)(v(2)=0,1) and DF(v('))+CHD(2)(v(4)=0,1) is obtained. Both experiment and theory show that the bending excitation activates the reaction at low E(c) and begins to inactivate at higher E(c). The experimental F+CHD(3)(v(b)=1) excitation functions display surprising peak features, especially for the HF(v(')=3)+CD(3)(v(2)=0,1) channels, indicating reactive resonances (quantum effects), which cannot be captured by quasiclassical calculations. The reactant state-specific QCT calculations predict that the v(5)(e) bending mode excitation is the most efficient to drive the reaction and the v(6)(e) and v(5)(e) modes enhance the DF and HF channels, respectively.     

Imaging the state-specific vibrational predissociation of the C2H2-NH3 hydrogen-bonded dimer

The state-to-state vibrational predissociation (VP) dynamics of the hydrogen-bonded ammonia-acetylene dimer were studied following excitation in the asymmetric CH stretch. Velocity map imaging (VMI) and resonance-enhanced multiphoton ionization (REMPI) were used to determine pair-correlated product energy distributions. Following vibrational excitation of the asymmetric CH stretch fundamental, ammonia fragments were detected by 2 + 1 REMPI via the B1E' <-- X1A1' and C'1A1' <-- X1A1' transitions. The fragments' center-of-mass (c.m.) translational energy distributions were determined from images of selected rotational levels of ammonia with one or two quanta in the symmetric bend (nu2 umbrella mode) and were converted to rotational-state distributions of the acetylene co-fragment. The latter is always generated with one or two quanta of bending excitation. All the distributions could be fit well when using a dimer dissociation energy of D0 = 900 +/- 10 cm(-1). Only channels with maximum translational energy <150 cm(-1) are observed. The rotational excitation in the ammonia fragments is modest and can be fit by temperatures of 150 +/- 50 and 50 +/- 20 K for 1nu2 and 2nu2, respectively. The rotational distributions in the acetylene co-fragment pair-correlated with specific rovibrational states of ammonia appear statistical as well. The vibrational-state distributions, however, show distinct state specificity among channels with low translational energy release. The predominant channel is NH3(1nu2) + C2H2(2nu4 or 1nu4 + 1nu5), where nu4 and nu5 are the trans- and cis-bend vibrations of acetylene, respectively. A second observed channel, with much lower population, is NH3(2nu2) + C2H2(1nu4). No products are generated in which the ammonia is in the vibrational ground state or the asymmetric bend (1nu4) state, nor is acetylene ever generated in the ground vibrational state or with CC stretch excitation. The angular momentum (AM) model of McCaffery and Marsh is used to estimate impact parameters in the internal collisions that give rise to the observed rotational distributions. These calculations show that dissociation takes place from bent geometries, which can also explain the propensity to excite fragment bending levels. The low recoil velocities associated with the observed channels facilitate energy exchange in the exit channel, which results in statistical-like fragment rotational distributions.     

Reaction dynamics and vibrational spectroscopy of CH3D molecules with both C-H and C-D stretches excited

State-resolved reactions of CH3D molecules containing both C-H and C-D stretching excitation with Cl atoms provide new vibrational spectroscopy and probe the consumption and disposal of vibrational energy in the reactions. The vibrational action spectra have three different components, the combination of the C-H symmetric stretch and the C-D stretch (nu1 + nu2), the combination of the C-D stretch and the C-H antisymmetric stretch (nu2 + nu4), and the combination of the C-D stretch and the first overtone of the CH3 bend (nu2 + 2nu5). The simulation for the previously unanalyzed (nu2 + nu4) state yields a band center of nu0 = 5215.3 cm(-1), rotational constants of A = 5.223 cm(-1) and B = 3.803 cm(-1), and a Coriolis coupling constant of zeta = 0.084. The reaction dynamics largely follow a spectator picture in which the surviving bond retains its initial vibrational excitation. In at least 80% of the reactive encounters of vibrationally excited CH3D with Cl, cleavage of the C-H bond produces CH2D radicals with an excited C-D stretch, and cleavage of the C-D bond produces CH3 radicals with an excited C-H stretch. Deviations from the spectator picture seem to reflect mixing in the initially prepared eigenstates and, possibly, collisional coupling during the reaction.     

Intramolecular vibrational energy redistribution in the highly excited fluoroform molecule: a quantum mechanical study using the multiconfiguration time-dependent Hartree algorithm

The present paper is devoted to a detailed study of the intramolecular vibrational energy redistribution in fluoroform initiated by a local mode excitation of the CH stretch [nnu(CH) (n=1,...,4)]. All nine internal degrees of freedom are explicitly taken into account and the full quantum mechanical simulation is performed by means of the multiconfiguration time-dependent Hartree algorithm. The existence of different time scales considerably complicates the dynamics. The mode-to-mode energy transfer is analyzed by calculating the evolution of the partial energies of all vibrational modes. This study emphasizes the crucial role played by the two-dimensional FCH bending modes which act as an energy reservoir. The fast energy flow into these bending modes significantly hinders an energy flow from the CH chromophore. Finally, our results are compared with those obtained previously with the wave operator sorting algorithm approach.     

Searching for resonances in the reaction Cl+CH4-->HCl+CH3: quantum versus quasiclassical dynamics and comparison with experiments

A quantum-mechanical (QM) and quasiclassical trajectory (QCT) study was performed on the title reaction, using a pseudotriatomic ab initio based surface. Probabilities and integral cross sections present some clear peaks versus the collision energy E(col), which we assign to Feshbach resonances of the transition state, where the light H atom oscillates between the heavy Cl and CH(3) groups. For ground-state reactants, reactivity is essentially of quantum origin (QCT observables and oscillations are smaller, or much smaller, than QM ones), and the calculated integral cross section and product distributions are in reasonable agreement with the experiment. The reaction occurs through an abstraction mechanism, following both a direct and an indirect mechanism. The quasiclassical trajectory calculations show the participation of a short-lived collision complex in the microscopic reaction mechanism. Finally, QCT differential cross sections of Cl+CH(4)-->HCl (nu(')=0 and 1)+CH(3) oscillate versus E(col), whereas experimentally this only occurs for HCl (nu(')=1). This theoretical result and other oscillating properties found here could, however, be related to the existence of a Feshbach resonance for the production of HCl (nu(')=1), as suggested by experimentalists.     

Resonance Raman intensity analysis of the excited state proton transfer dynamics of 2-nitrophenol in the charge-transfer band absorption

Resonance Raman spectra were obtained for 2-nitrophenol in cyclohexane solution with excitation wavelengths in resonance with the charge-transfer (CT) proton transfer band absorption. These spectra indicate that the Franck-Condon region photodissociation dynamics have multidimensional character with motion along more than 15 normal modes: the nominal CCH bend+CC stretch nu(12) (1326 cm(-1)), the nominal CCC bend nu(23) (564 cm(-1)), the nominal CO stretch+NO stretch+CC stretch nu(14) (1250 cm(-1)), the nominal CCH bend+CC stretch+COH bend nu(15) (1190 cm(-1)); the nominal CCH bend+CC stretch nu(17) (1134 cm(-1)), the nominal CCC bend+CC stretch nu(22) (669 cm(-1)), the nominal CCN bend nu(27) (290 cm(-1)), the nominal NO(2) bend+CC stretch nu(21) (820 cm(-1)), the nominal CCO bend+CNO bend nu(25) (428 cm(-1)), the nominal CC stretch nu(7) (1590 cm(-1)), the nominal NO stretch nu(8) (1538 cm(-1)), the nominal CCC bend+NO(2) bend nu(20) (870 cm(-1)), the nominal CC stretch nu(6) (1617 cm(-1)), the nominal COH bend+CC stretch nu(11) (1382 cm(-1)), nominal CCH bend+CC stretch nu(9) (1472 cm(-1)). A preliminary resonance Raman intensity analysis was done and the results for 2-nitrophenol were compared to previously reported results for nitrobenzene, p-nitroaniline, and 2-hydroxyacetophenone. The authors briefly discuss the differences and similarities in the CT-band absorption excitation of 2-nitrophenol relative to those of nitrobenzene, p-nitroaniline, and 2-hydroxyacetophenone.     

Vibrational predissociation dynamics of methane-Ar: an ab initio approach

We calculated the cross sections for vibrational predissociation of methane-Ar induced by excitation of the methane nu 3 mode with the aid of an ab initio CH4-Ar potential depending explicitly on the nu 3 and nu 1 normal coordinates of the CH4 monomer. We found that dissociation into CH4 fragments excited in the nu 1 mode, a V-->V' process with very low kinetic energy release, strongly dominates over direct dissociation into Ar and ground state CH4, and is responsible for the line broadening observed experimentally. The (observed and calculated) strong variation of the line widths for the Van der Waals levels excited in combination with the nu 3 mode (giving states of A, F and E symmetry) is related to the opening up of appropriate nu 1 dissociation channels and the occurrence of rotational resonances in the nu 1 continuum in the energy range of the quasi-bound nu 3 levels.     

Six-dimensional potential energy surface and rovibrational energies of the HCCN radical in the ground electronic state

We report large-scale quantum mechanical calculations for the HCCN radical in its ground electronic state. A six-dimensional potential energy surface based on MR-ACPF/cc-pVQZ ab initio energy points is developed and adjusted to reproduce experimental findings for and nu1 of HCCN. Rovibrational energy levels of HCCN and DCCN are computed for total rotational angular momentum J = 0-4 by making use of combined (functional + point wise) coordinate representations together with contraction schemes resulting from several diagonalization/truncation steps. The classical barrier to linearity is determined to be 287 cm(-1). Spectroscopic parameters are calculated for low lying states and compared with available experimental data. Energy patterns attributed to the nu4 bending mode and to the quasilinear nu5 bending mode are identified. It has been also found that nu2 and nu3 + (nu4(1),nu5(1))(0,0) are coupled in HCCN, while the mixing between nu3 and (2nu4(0), 2nu5(0))(0,0) is seen in DCCN.     

Collaborative control of branching ratio in the O + HO2 → OH + O2 reaction via vibrational and rotational excitation

To understand the effect of different vibrational and rotational modes of reactant to enhance the reactivity of the O + HO₂ → OH + O₂ reaction, we revisited this important atmospheric reaction. We report here a quasi-classical trajectory (QCT) study of the reaction dynamics on a recently developed full-dimensional potential energy surface (PES). Our previous work has indicated that this reaction has two pathways, the H abstraction (HA) channel and the O abstraction (OA) channel, which lead to totally different product energy distribution. In this work, we identified that the vibrational excitation of the OH stretching (v₁) mode of HO₂ is the switch of the HA channel at low collision energy; meanwhile, the rotational excitation can also greatly change the branching ratio of the two pathways. With the excitation of v₁ mode, the original negligible HA channel controlled by the tight transition state becomes quite important. This work presents an approach to control the branching ratio via collaboration between vibrational and rotational excitation and will enrich the knowledge of the O + HO₂ reaction in atmospheric chemistry and physics.