State-to-state photodissociation dynamics of triatomic molecules: H2O in the B band

State-to-state photodissociation dynamics of H(2)O in its B band has been investigated quantum mechanically on a new set of non-adiabatically coupled potential energy surfaces for the lowest two (1)A' states of H(2)O, which are developed at the internally contracted multi-reference configuration interaction level with the aug-cc-pVQZ basis set. Quantum dynamical calculations carried out using the Chebyshev propagator yield absorption spectra, product state distributions, branching ratios, and differential cross sections, which are in reasonably good agreement with the latest experimental results. Particular focus is placed here on the dependence of various dynamical observables on the photon energy. Detailed analyses of the dynamics have assigned the diffuse structure in absorption spectrum to short-time recurring dynamics near the HOH conical intersection. The non-adiabatic dissociation to the ground state OH product via the HOH conical intersection is facile, direct, fast, and produces rotationally hot OH(X?) products. On the other hand, the adiabatic channel on the excited state leading to the OH(A?) product is dominated by long-lived resonances, which depend sensitively on the potential energy surfaces.     

Potential surfaces and dynamics of the O(3P)+H2O(X1A1)-->OH(X2pi)+OH(X2pi) reaction

We present global potential energy surfaces for the three lowest triplet states in O(3P)+H2O(X1A1) collisions and present results of classical dynamics calculations on the O(3P)+H2O(X1A1)-->OH(X2pi)+OH(X2pi) reaction using these surfaces. The surfaces are spline-based fits of approximately 20,000 fixed geometry ab initio calculations at the complete-active-space self-consistent field+second-order perturbation theory (CASSCF+MP2) level with a O(4s3p2d1f)/H(3s2p) one electron basis set. Computed rate constants compare well to measurements in the 1000-2500 K range using these surfaces. We also compute the total, rovibrationally resolved, and differential angular cross sections at fixed collision velocities from near threshold at approximately 4 km s(-1) (16.9 kcal mol(-1) collision energy) to 11 km s(-1) (122.5 kcal mol(-1) collision energy), and we compare these computed cross sections to available space-based and laboratory data. A major finding of the present work is that above approximately 40 kcal mol(-1) collision energy rovibrationally excited OH(X2pi) products are a significant and perhaps dominant contributor to the observed 1-5 micro spectral emission from O(3P)+H2O(X1A1) collisions. Another important result is that OH(X2pi) products are formed in two distinct rovibrational distributions. The "active" OH products are formed with the reagent O atom, and their rovibrational distributions are extremely hot. The remaining "spectator" OH is relatively rovibrationally cold. For the active OH, rotational energy is dominant at all collision velocities, but the opposite holds for the spectator OH. Summed over both OH products, below approximately 50 kcal mol(-1) collision energy, vibration dominates the OH internal energy, and above approximately 50 kcal mol(-1) rotation is greater than vibrational energy. As the collision energy increases, energy is diverted from vibration to mostly translational energy. We note that the present fitted surfaces can also be used to investigate direct collisional excitation of H2O(X1A1) by O(3P) and also OH(X2pi)+OH(X2pi) collisions.     

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 and classical studies of the O(3P) + H2(v = 0-3,j = 0) --> OH + H reaction using benchmark potential surfaces

We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.     

Transition state dynamics of the OH + H2O hydrogen exchange reaction studied by dissociative photodetachment of H3O2-

Dynamics in the transition state region of the bimolecular OH + H2O-->H2O + OH hydrogen exchange reaction have been studied by photoelectron-photofragment coincidence spectroscopy of the H3O2- negative ion and its deuterated analog D3O2-. The data reveal vibrationally resolved product translational energy distributions. The total translational energy distribution shows a vibrational progression indicating excitation of the antisymmetric stretch of the water product. Electronic structure calculations at the QCISD level of theory support this analysis. Examination of the translational energy release between the neutral products reveals a dependence on the product vibrational state. These data should provide a critical test of ab initio potential energy surfaces and dynamics calculations.     

Imaging H2O photofragments in the predissociation of the HCl-H2O hydrogen-bonded dimer

The state-to-state vibrational predissociation (VP) dynamics of the hydrogen-bonded HCl-H(2)O dimer was studied following excitation of the dimer's HCl stretch by detecting the H(2)O fragment. Velocity map imaging (VMI) and resonance-enhanced multiphoton ionization (REMPI) were used to determine pair-correlated product energy distributions. Following vibrational excitation of the HCl stretch of the dimer, H(2)O fragments were detected by 2 + 1 REMPI via the C (1)B(1) (000) ← X (1)A(1) (000) transition. REMPI spectra clearly show H(2)O from dissociation produced in the ground vibrational state. The fragments' center-of-mass (c.m.) translational energy distributions were determined from images of selected rotational states of H(2)O and were converted to rotational state distributions of the HCl cofragment. The distributions were consistent with the previously measured dissociation energy of D(0) = 1334 ± 10 cm(-1) and show a clear preference for rotational levels in the HCl fragment that minimize translational energy release. The usefulness of 2 + 1 REMPI detection of water fragments is discussed.     

Differential and integral cross sections for the H + O2 --> OH + O combustion reaction

We present accurate differential and integral cross sections for the H + O2 --> OH + O reaction obtained on a newly developed ab initio potential energy surface using time-independent and time-dependent quantum mechanical methods. The product angular distributions near the reaction threshold show pronounced forward and backward peaks, reflecting the complex-forming mechanism. However, the asymmetry of these peaks suggests certain nonstatistical behaviors, presumably due to some relatively short-lived resonances. The integral cross section increases monotonically with the collision energy above a reaction threshold.     

Photofragment coincidence imaging of small I-(H2O)n clusters excited to the charge-transfer-to-solvent state

The photodissociation dynamics of small I-(H2O)n(n=2-5) clusters excited to their charge-transfer-to-solvent (CTTS) states have been studied using photofragment coincidence imaging. Upon excitation to the CTTS state, two photodissociation channels were observed. The major channel (approximately 90%) is a two-body process forming neutral I+(H2O)n photofragments, and the minor channel is a three-body process forming I+(H2O)n-1+H2O fragments. Both processes display translational energy [P(ET)] distributions peaking at ET=0 with little available energy partitioned into translation. Clusters excited to the detachment continuum rather than to the CTTS state display the same two channels with similar P(ET) distributions. The observation of similar P(ET) distributions from the two sets of experiments suggests that in the CTTS experiments, I atom loss occurs after autodetachment of the excited [I(H2O)n-]* cluster or, less probably, that the presence of the excess electron has little effect on the departing I atom.     

A reduced dimensionality quasiclassical and quantum study of the proton transfer reaction H3O(+)+H2O-->H2O+H3O(+)

We report quantum and quasiclassical calculations of proton transfer in the reaction H(3)O(+)+H(2)O in three degrees of freedom, the two OH(+) bond lengths and the OH(+)O angle. The reduced dimensional potential energy surface is obtained from the full dimensional OSS3(p) energy function of H(5)O(2) (+) [L. Ojamae, I. Shavitt, and S. J. Singer, J. Chem. Phys. 109, 5547 (1998)], with an additional long-range correction to reproduce the correct ion-molecule interaction. This surface is used to perform both quasiclassical trajectory and quantum reactive scattering calculations of the zero total angular momentum cumulative reaction probability and cross sections for initial rotational states 0, 1, and 2. Comparison of these quantities are made to assess the importance of quantum effects in this reduced dimensional reaction. Additional quasiclassical cross sections are calculated to obtain the thermal rate constant for the reaction.     

Dynamics of hyperthermal collisions of O(3P) with CO

The dynamics of O(3P) + CO collisions at a hyperthermal collision energy near 80 kcal mol-1 have been studied with a crossed molecular beams experiment and with quasi-classical trajectory calculations on computed potential energy surfaces. In the experiment, a rotatable mass spectrometer detector was used to monitor inelastically and reactively scattered products as a function of velocity and scattering angle. From these data, center-of-mass (c.m.) translational energy and angular distributions were derived for the inelastic and reactive channels. Isotopically labeled C18O was used to distinguish the reactive channel (16O + C18O 16OC + 18O) from the inelastic channel (16O + C18O 16O + C18O). The reactive 16OC molecules scattered predominantly in the forward direction, i.e., in the same direction as the velocity vector of the reagent O atoms in the c.m. frame. The c.m. translational energy distribution of the reactively scattered 16OC and 18O was very broad, indicating that 16OC is formed with a wide range of internal energies, with an average internal excitation of approximately 40% of the available energy. The c.m. translational energy distribution of the inelastically scattered C18O and 16O products indicated that an average of 15% of the collision energy went into internal excitation of C18O, although a small fraction of the collisions transferred nearly all the collision energy into internal excitation of C18O. The theoretical calculations, which extend previously published results on this system, predict c.m. translational energy and angular distributions that are in near quantitative agreement with the experimentally derived distributions. The theoretical calculations, thus validated by the experimental results, have been used to derive internal state distributions of scattered CO products and to probe in detail the interactions that lead to the observed dynamical behavior.     

Communication: quasiclassical trajectory calculations of correlated product-state distributions for the dissociation of (H2O)2 and (D2O)2

Stimulated by recent experiments [B. E. Rocher-Casterline, L. C. Ch'ng, A. K. Mollner, and H. Reisler, J. Chem. Phys. 134, 211101 (2011)], we report quasiclassical trajectory calculations of the dissociation dynamics of the water dimer, (H(2)O)(2) (and also (D(2)O)(2)) using a full-dimensional ab initio potential energy surface. The dissociation is initiated by exciting the H-bonded OH(OD)-stretch, as done experimentally for (H(2)O)(2). Normal mode analysis of the fragment pairs is done and the correlated vibrational populations are obtained by (a) standard histogram binning (HB), (b) harmonic normal-mode energy-based Gaussian binning (GB), and (c) a modified version of (b) using accurate vibrational energies obtained in the Cartesian space. We show that HB allows opening quantum mechanically closed states, whereas GB, especially via (c), gives physically correct results. Dissociation of both (H(2)O)(2) and (D(2)O)(2) mainly produces either fragment in the bending excited (010) state. The H(2)O(J) and D(2)O(J) rotational distributions are similar, peaking at J = 3-5. The computations do not show significant difference between the ro-vibrational distributions of the donor and acceptor fragments. Diffusion Monte Carlo computations are performed for (D(2)O)(2) providing an accurate zero-point energy of 7247 cm(-1), and thus, a benchmark D(0) of 1244 ± 5 cm(-1).     

Reaction dynamics simulations of the identity S(N)2 reaction H(2)O + HOOH(2)(+)--> H(2)OOH(+)+ H(2)O. Requirements for reaction and competition with proton transfer

Despite the fact that the transition structure of the gas phase S(N)2 reaction H(2)O + HOOH(2)(+)--> HOOH(2)(+)+ H(2)O is well below the reactants in potential energy, the reaction has not yet been observed by experiment. Variational transition state RRKM theory reveals a strong preference for the competing proton transfer reaction H(2)O + HOOH(2)(+)--> H(3)O(+)+ HOOH due to entropy factors. Born-Oppenheimer reaction dynamics simulations confirm these results. However, by increasing the collision energy to around 7.5 eV the probability for nucleophilic substitution increases relative to proton transfer. These observations are explained by the presence of the key common intermediate HOO(H)[dot dot dot]H-OH(2)(+) which leads to effective proton transfer, but can be avoided with increasing collision energy. However, the S(N)2 probability remains below 0.2 since successful passage through the TS requires optimum initial orientation of the reactants, excitation of the relative translational motion and good phase correlation between the O-O vibration and the motion of the incoming water.     

Water line lists close to experimental accuracy using a spectroscopically determined potential energy surface for H2(16)O, H2(17)O, and H2(18)O

Line lists of vibration-rotation transitions for the H(2) (16)O, H(2) (17)O, and H(2) (18)O isotopologues of the water molecule are calculated, which cover the frequency region of 0-20 000 cm(-1) and with rotational states up to J=20 (J=30 for H(2) (16)O). These variational calculations are based on a new semitheoretical potential energy surface obtained by morphing a high accuracy ab initio potential using experimental energy levels. This potential reproduces the energy levels with J=0, 2, and 5 used in the fit with a standard deviation of 0.025 cm(-1). Linestrengths are obtained using an ab initio dipole moment surface. That these line lists make an excellent starting point for spectroscopic modeling and analysis of rotation-vibration spectra is demonstrated by comparison with recent measurements of Lisak and Hodges [J. Mol. Spectrosc. (unpublished)]: assignments are given for the seven unassigned transitions and the intensity of the strong lines are reproduced to with 3%. It is suggested that the present procedure may be a better route to reliable line intensities than laboratory measurements.     

Femtosecond study of Cu(H(2)O) dynamics

The short-time nuclear dynamics of Cu(H(2)O) is investigated using femtosecond photodetachment-photoionization spectroscopy and time-dependent quantum wave packet calculations. The Cu(H(2)O) dynamics is initiated in the electronic ground state of the complex by electron photodetachment from the Cu(-)(H(2)O) complex, where hydrogen atoms are oriented toward Cu. Several time-resolved resonant multiphoton ionization schemes are used to probe the ensuing reorientation and dissociation. Immediately following photodetachment, the neutral complex is far from its minimum energy geometry and possesses an internal energy comparable to the Cu-H(2)O dissociation energy and undergoes both large-amplitude H(2)O motion and dissociation. Dissociation is observed to occur on three distinct time scales: 0.6, 8, and 100 ps. These results are compared to the results of time-dependent J=0 wave packet calculations, propagating the initial anion vibrational wave functions on the ground-state potential of the neutral complex. An excellent agreement is obtained between the experimental results and the ionization signals derived from the calculated probability amplitudes. Related experiments and calculations are carried out on the Cu(D(2)O) complex, with results very similar to those of Cu(H(2)O).     

Quantum/classical studies of vibrationally mediated photodissociation of Ar x H2O

Results of multiple configuration quantum/classical simulations of the dynamics of Ar x H2O photodissociation are reported. In agreement with experimental studies of Nesbitt and co-workers [J. Chem. Phys. 2000, 112, 7449], we find that the OH products emerge rotationally excited, compared to the dissociation of bare H2O. The wavelength dependence of the total cross section and the energy transfer to the argon atom are also investigated. The trends are interpreted in terms of features in the Ar x H2O A state potential surface.     

Quantum dynamics of proton migration in H2O dications: H2+ formation on ultrafast timescales

Irradiation of isolated water molecules by few-cycle pulses of intense infrared laser light can give rise to ultrafast rearrangement resulting in formation of the H(2) (+) ion. Such unimolecular reactions occur on the potential energy surface of the H(2)O(2+) dication that is accessed when peak laser intensities in the 10(15) W cm(-2) range and pulse durations as short as 9-10 fs are used; ion yields of ~1.5% relative to the H(2)O(+) ion are measured. We also study such reactions by means of time-dependent wavepacket dynamics on an ab initio potential energy surface of the dication and show that a proton, generated from O-H bond rupture, migrates towards the H-atom, and forms vibrationally excited H(2)(+) in a well-defined spatial zone.     

Crossed-beam radical-radical reaction dynamics of O(3P)+C3H3-->H(2S)+C3H2O

The radical-radical oxidation reaction, O(3P)+C3H3 (propargyl)-->H(2S)+C3H2O (propynal), was investigated using vacuum-ultraviolet laser-induced fluorescence spectroscopy in a crossed-beam configuration, together with ab initio and statistical calculations. The barrierless addition of O(3P) to C3H3 is calculated to form energy-rich addition complexes on the lowest doublet potential energy surface, which subsequently undergo direct decomposition steps leading to the major reaction products, H+C3H(2)O (propynal). According to the nascent H-atom Doppler-profile analysis, the average translational energy of the products and the fraction of the average transitional energy to the total available energy were determined to be 5.09+/-0.36 kcal/mol and 0.077, respectively. On the basis of a comparison with statistical prior calculations, the reaction mechanism and the significant internal excitation of the polyatomic propynal product can be rationalized in terms of the formation of highly activated, short-lived addition-complex intermediates and the adiabaticity of the excess available energy along the reaction coordinate.     

A theoretical and computational study of the anion, neutral, and cation Cu(H(2)O) complexes

An ab initio investigation of the potential energy surfaces and vibrational energies and wave functions of the anion, neutral, and cation Cu(H(2)O) complexes is presented. The equilibrium geometries and harmonic frequencies of the three charge states of Cu(H(2)O) are calculated at the MP2 level of theory. CCSD(T) calculations predict a vertical electron detachment energy for the anion complex of 1.65 eV and a vertical ionization potential for the neutral complex of 6.27 eV. Potential energy surfaces are calculated for the three charge states of the copper-water complexes. These potential energy surfaces are used in variational calculations of the vibrational wave functions and energies and from these, the dissociation energies D(0) of the anion, neutral, and cation charge states of Cu(H(2)O) are predicted to be 0.39, 0.16, and 1.74 eV, respectively. In addition, the vertical excitation energies, that correspond to the 4 (2)P<--4 (2)S transition of the copper atom, and ionization potentials of the neutral Cu(H(2)O) are calculated over a range of Cu(H(2)O) configurations. In hydrogen-bonded, Cu-HOH configurations, the vertical excitation and ionization energies are blueshifted with respect to the corresponding values for atomic copper, and in Cu-OH(2) configurations where the copper atom is located near the oxygen end of water, both quantities are redshifted.