Angular distributions for the F+H(2)-->HF+H reaction: the role of the F spin-orbit excited state and comparison with molecular beam experiments

We report quantum mechanical calculations of center-of-mass differential cross sections (DCS) for the F+H(2)-->HF+H reaction performed on the multistate [Alexander-Stark-Werner (ASW)] potential energy surfaces (PES) that describe the open-shell character of this reaction. For comparison, we repeat single-state calculations with the Stark-Werner (SW) and Hartke-Stark-Werner (HSW) PESs. The ASW DCSs differ from those predicted for the SW and HSW PES in the backward direction. These differences arise from nonadiabatic coupling between several electronic states. The DCSs are then used in forward simulations of the laboratory-frame angular distributions (ADs) measured by Lee, Neumark, and co-workers [J. Chem. Phys. 82, 3045 (1985)]. The simulations are scaled to match experiment over the range 12 degrees 相似文献   

Role of the F spin-orbit excited state in the F+HD reaction: contributions to the dynamical resonance

We report quantum mechanical calculations of excitation functions (relative reaction cross sections) for the F+HD reaction. We include three potential energy surfaces and an accurate treatment of all couplings (non-adiabatic, spin-orbit, and Coriolis). Comparison with experimental results [Dong, Lee, and Liu, J. Chem. Phys., 113, 3633 (2000)] show excellent agreement for the DF product channel and an improved but not perfect agreement for the HF product channel. In the former case, when weighted by the (16%) fractional population of the spin-orbit excited state (F(*)) in the beam, the overall reactivity of the F(*) is small (approximately 5%). For the HF product channel and with the same (16%) fractional weight, F(*) reactivity makes a contribution of approximately 12% in the high-energy tail of the resonance peak. As a result, averaging over the population of F spin-orbit states in the beam changes the shape of the resonance. The greater the fraction of F(*) in the beam, the less pronounced will be the resonance modulation of the reaction excitation function.     

Decay dynamics of the long-range H1Sigmag+ state of D2 and H2: experiment and theory

We present accurate experimental measurements of the lifetimes of rovibrational levels of the long-range H1Sigmag+ state for both D2 and H2, obtained directly from the observation of the time-dependent decay of the fluorescence from these excited levels. These results improve upon and extend those of Reinhold et al. [J. Chem. Phys. 112, 10754 (2000)]. Several decay pathways are open to these levels including fluorescence, predissociation, and autoionization. We present theoretical results for each of these processes, each calculated using the simplest but still appropriate level of theory. In particular, the theoretical calculations provide a quantitative explanation of the dramatic vibrational dependence of the observed lifetimes, the isotope dependence of the lifetimes for levels well localized within the H potential well and therefore not subject to significant tunneling, and an insight into the role of enhanced tunneling in autoionization. In these calculations each of the rovibrational levels of the H state is treated individually, without having to engage in a global coupled-state calculation.     

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.     

Quantum transition state theory for the full three-dimensional H+H2 reaction

A recently developed quantum transition state theory (QTST) [E. Pollak and J. L. Liao, J. Chem. Phys. 108, 2733 (1998)] for calculating thermal rate constants of chemical reactions is applied to the full three-dimensional hydrogen exchange reaction. Results are compared with other numerical results, for temperatures ranging from T=300 K to T=1500 K. The QTST rate is almost exact at high temperature and is 20% greater than the exact rate at T=300 K, where there is extensive tunneling.     

Resonance spectrum and dissociation dynamics of ozone in the 3B2 electronically excited state: experiment and theory

The rovibrational spectrum assigned to the low-lying (3)B(2) electronic state of ozone is measured with intracavity laser absorption spectroscopy. The experimental results are interpreted by means of quantum dynamical calculations on a global ab initio potential energy surface. The observed spectrum is shown to originate from the vibrational ground state in the local minimum of the (3)B(2) potential. The spectrum of short-lived resonance states in this local minimum is analyzed. Additionally, the global minimum of the surface is shown to lie in the dissociation channel in the van der Waals region. This region supports a short sequence of weakly bound vibrational states.     

The microscopic reaction dynamics and branching ratio in the F + HCOOH and F + H2CO reactions

Infrared chemiluminescence under conditions of arrested relaxation has been applied to the study of the hydrogen and deuterium abstraction reactions of HCOOH, DCOOH and H₂CO with F atoms. Two distinctly different modes of product excitation are observed, depending upon whether the reaction proceeds via the formyl or carboxyl hydrogen. Reaction at the formyl hydrogen (or deuterium) causes substantial inversion in the diatomic product internal energy distributions. The F + H₂CO and F + DCOOH reactions respectively channel 56% and 54% of the available energy into vibration in the product diatomic when they occur at the formyl site. In both cases the product energy distributions are qualitatively similar to those observed in direct reactions of triatomic systems on repulsive energy surfaces. In contrast to these, reaction at the carboxyl hydrogen of DCOOH gives an HF2 product vibrational distribution having a Boltzmann equilibrium shape at a temperature of 4300 K. The ratio of HF to DF product from the F + DCOOH study shows that reaction occurs at the carboxyl hydrogen approximately twice as often as at the formyl site. Comparison with triatomic reactions involving the same mass-combinations implies that abstraction of the formyl hydrogen occurs via single-collision, direct encounters, whereas reaction at the carboxyl site involves a long-lived complex in which extensive randomisation of the reaction exoergicity among all the product vibrational modes can occur.     

An improved calculation of the transition state for the F + H2 reaction

Using a large, balanced one-electron basis set, fully optimized reaction space (FORS) calculations to optimize the orbitals and to estimate the internal correlation energy, multi-reference configuration interaction calculations including all single and double excitations out of the FORS reference space to estimate a fraction of the external correlation energy, and the method of scaled external correlation (SEC), we calculate the interaction energy of F with H₂ in the vicinity of the saddle point for the reaction F + H₂ → HF + H. Our calculated barrier height, 1.6 kcal/mol, is considerably lower than values obtained in recent ab initio calculations, and the saddle point geometry is about 0.3 a₀ looser. This indicates that the part of the external correlation energy omitted from MR CISD calculations because of the incompleteness of the one-electron basis set and the truncation of the CI expansion, as estimated by the SEC method, has a significant effect on both the saddle point energy and its geometry.     

Comparison of theory and experiment for excited singlet states of the H2 molecule

This paper presents a survey of published and unpublished ab initio calculations of the vibrational structures of the ten lowest electronic singlet states of the hydrogen molecule up to the H(n = 1) + H(n = 2) dissociation limit. The data are based on adiabatic potential functions (clamped-nuclei electronic energies and nuclear-mass-dependent diagonal corrections). Nonadiabatic coupling has been treated ab initio within the five states. of ¹Λ symmetry (X,EF, GK, HH?) and ¹Σ I.¹Π_g. The accuracies of the theoretical energies are determined by comparisons with experimental data for H₂, HD, and D₂. The level shifts and predissociation probabilities of the excited ¹Σ states, generated by nonadiabatic coupling with the discrete and continuous vibrational structure of the ground state, and radiative properties have also been calculated.     

Mixed quantum-classical reaction path dynamics of C2H5F --> C2H4 + HF

A mixed quantum-classical method for calculating product energy partitioning based on a reaction path Hamiltonian is presented and applied to HF elimination from fluoroethane. The goal is to describe the effect of the potential energy release on the product energies using a simple model of quantized transverse vibrational modes coupled to a classical reaction path via the path curvature. Calculations of the minimum energy path were done at the B3LYP/6-311++G(2d,2p) and MP2/6-311++G** levels of theory, followed by energy-partitioning dynamics calculations. The results for the final HF vibrational state distribution were found to be in good qualitative agreement with both experimental studies and quasiclassical trajectory simulations.     

Ion-molecule chemistry of C+2: Evidence for the production of excited state dicarbon cation

Ion-molecule reactions of C⁺₂ with several neutral partners have been studied using Fourier transform mass spectrometry. The reactant ion was formed by electron impact of various neutral precursors and the internal energy content of the ion estimated using charge transfer/energy-bracketing reactions. These reactions indicate the production of a long-lived excited state ion when C⁺₂ is formed from C₂N₂ and the production of a mixture of states when formed from small hydrocarbon molecules. The observed reactions and energetics are consistent with the calculated electronic structure of this ion.     

Dynamics of the D+ + H2 and H+ + D2 reactions: a detailed comparison between theory and experiment

An extensive set of experimental measurements on the dynamics of the H(+) + D(2) and D(+) + H(2) ion-molecule reactions is compared with the results of quantum mechanical (QM), quasiclassical trajectory (QCT), and statistical quasiclassical trajectory (SQCT) calculations. The dynamical observables considered include specific rate coefficients as a function of the translational energy, E(T), thermal rate coefficients in the 100-500 K temperature range. In addition, kinetic energy spectra (KES) of the D(+) ions reactively scattered in H(+) + D(2) collisions are also presented for translational energies between 0.4 eV and 2.0 eV. For the two reactions, the best global agreement between experiment and theory over the whole energy range corresponds to the QCT calculations using a gaussian binning (GB) procedure, which gives more weight to trajectories whose product vibrational action is closer to the actual integer QM values. The QM calculations also perform well, although somewhat worse over the more limited range of translational energies where they are available (E(T) < 0.6 eV and E(T) < 0.2 eV for the H(+) + D(2) and D(+) + H(2) reactions, respectively). The worst agreement is obtained with the SQCT method, which is only adequate for low translational energies. The comparison between theory and experiment also suggests that the most reliable rate coefficient measurements are those obtained with the merged beams technique. It is worth noting that none of the theoretical approaches can account satisfactorily for the experimental specific rate coefficients of H(+) + D(2) for E(T)≤ 0.2 eV although there is a considerable scatter in the existing measurements. On the whole, the best agreement with the experimental laboratory KES is obtained with the simulations carried out using the state resolved differential cross sections (DCSs) calculated with the QCT-GB method, which seems to account for most of the observed features. In contrast, the simulations with the SQCT data predict kinetic energy spectra (KES) considerably cooler than those experimentally determined.     

Cumulative reaction probabilities and transition state properties: a study of the F + H2 reaction and its deuterated isotopic variants

A comparative quantum mechanical (QM) and quasiclassical trajectory (QCT) study of the cumulative reaction probabilities (CRPs) is presented in this work for the F + H(2) reaction and its isotopic variants for low values of the total angular momentum J. The agreement between the two sets of calculations is very good with the exception of some features whose origin is genuinely QM. The agreement also extends to the CRP resolved in the helicity quantum number k. The most remarkable feature is the steplike structure, which becomes clearly distinct when the CRPs are resolved in odd and even rotational states j. The analysis of these steps shows that each successive increment is due to the opening of the consecutive rovibrational states of the H(2) or D(2) molecule, which, in this case, nearly coincide with those of the transition state. Moreover, the height of each step reflects the number of helicity states compatible with a given J and j values, thus indicating that the various helicity states for a specific j have basically the same contribution to the CRPs at a given total energy. As a consequence, the dependence with k of the reactivity is practically negligible, suggesting very small steric restrictions for any possible orientation of the reactants. This behavior is in marked contrast to that found in the D + H(2) reaction, wherein a strong k dependence was found in the threshold and magnitude of the CRP. The advantages of a combined QCT and QM approaches to the study of CRPs are emphasized in this work.     

The formation of highly excited H in the reaction H2+(v) + H2 → H + H

A quasiclassical trajectory surface hopping method has been used to study H(v) + H₂ → H + H for v = 0, 3, 7, 10, 13, and 17 with an emphasis on determining the H internal energy and angular momentum distributions for high v. For v = 13 and 17, significant cross sections are found for producing H at energies above its dissociation energy. An average metastable H lifetime of 11.5 ps for v = 13 and 4.7 ps for v = 17 is found, but there is also a much longer lived component to the lifetime distributions that is more important for v = 13 than for v = 17. Some of the longer lived metastables correspond to high angular momentum orbiting states of H, but other sources of metastability are also present.     

Collision theory treatment of the H + H2 exchange reaction

The rate constant of the three-dimensional H + H₂ reaction is calculated accurately in the temperature range 100–600 K. By adjusting the value of an effective collision diameter one obtains good agreement of the theory with experimental data. The precise corrections to the simplest approximate kinetic theories are also computed.     

Ground state structures and excited state dynamics of pyrrole-water complexes: ab initio excited state molecular dynamics simulations

Structures of the ground state pyrrole-(H2O)n clusters are investigated using ab initio calculations. The charge-transfer driven femtosecond scale dynamics are studied with excited state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field method for pyrrole-(H2O)n clusters. Upon the excitation of these clusters, the charge density is located over the farthest water molecule which is repelled by the depleted pi-electron cloud of pyrrole ring, resulting in a highly polarized complex. For pyrrole-(H2O), the charge transfer is maximized (up to 0.34 a.u.) around approximately 100 fs and then oscillates. For pyrrole-(H2O)2, the initial charge transfer occurs through the space between the pyrrole and the pi H-bonded water molecule and then the charge transfer takes place from this water molecule to the sigma H-bonded water molecule. The total charge transfer from the pyrrole to the water molecules is maximized (up to 0.53 a.u.) around approximately 100 fs.     

On the dynamics of the H+ +D2(v=0,j=0)-->HD+D + reaction: a comparison between theory and experiment

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.     

An improved potential energy surface for the F+H2 reaction

A new ground state potential energy surface has been developed for the F+H₂ reaction. Using the UCCSD(T) method, ab initio calculations were performed for 786 geometries located mainly in the exit channel of the reaction. The new data was used to correct exit channel errors that have become apparent in the potential energy surface of Stark and Werner [J. Chem. Phys. 104 (1996) 6515]. While the entrance channel and saddlepoint properties of the Stark–Werner surface are unchanged on the new potential, the exit channel behavior is more satisfactory. The exothermicity on the new surface is much closer to the experimental value. The new surface also greatly diminishes the exit channel van der Waals well that was too pronounced on the Stark–Werner surface. Several preliminary dynamical scattering calculations were carried out using the new surface for total angular momentum equal to zero for F+H₂ and F+HD. It is found that gross features of the reaction dynamics are quite similar to those predicted by the Stark–Werner surface, in particular the reactive resonance for F+HD and F+H₂ survive. However, the most of the exit channel van der Waals resonances disappear on the new surface. It is predicted that the differential cross-sections at low collision energy for the F+H₂ reaction may be drastically modified from the predictions based on the Stark–Werner surface.     

Quantum dynamics study of H+NH3-->H2+NH2 reaction

We report in this paper a quantum dynamics study for the reaction H+NH3-->NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].