A quantum reaction dynamics study of the translational, vibrational, and rotational motion effects on the HD + H3+ reaction

Time-dependent, quantum reaction dynamics wavepacket approach is employed to investigate the impacts of the translational, vibrational, and rotational motion on the HD+H(3)(+) → H(2)D(+) + H(2) reaction using the Xie-Braams-Bowman potential energy surface [Z. Xie, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 122, 224307 (2005)]. We treat this five atom reaction with a seven-degree-of-freedom model by fixing one Jacobi and one torsion angle related to H(3) (+) at the lowest saddle point geometry of the potential energy surface. The initial state selected reaction probabilities show that the rotational excitations of H(+)-H(2) greatly enhance the reactivity with the reaction probabilities increased double at high rotational states compared to the ground state. However, the vibrational excitations of H(3) (+) hinder the reactivity. The ground state reaction probability shows no reaction threshold for this exoergic reaction, and as the translational energy increases, the reaction probability decreases. Furthermore, reactive resonances and zero point energy play very important roles on the reaction dynamics. The obtained integral cross section has the character of an exoergic reaction without a threshold: it decreases with the translational energy increasing. The calculated thermal rate constants using this seven-degree-of-freedom model are in agreement with a later experiment measurement.     

Collisional quenching of OD A 2Σ+ by H2: experimental and theoretical studies of the state-resolved OD X 2Π product distribution and branching fraction

We report joint experimental and theoretical studies of outcomes resulting from the nonreactive quenching of electronically excited OD?A (2)Σ(+) by H(2). The experiments utilize a pump-probe technique to detect the OD?X (2)Π product state distribution under single collision conditions. The OD?X (2)Π products are observed primarily in their lowest vibrational state (v(") = 0) with substantially less population in v(") = 1. The OD?X (2)Π products are generated with a high degree of rotational excitation, peaking at N(") = 21 with an average rotational energy of 4600 cm(-1), and a strong propensity for populating the Π(A(')) Λ-doublet component indicative of alignment of the half-filled pπ orbital in the plane of OD rotation. Branching fraction measurements show that the nonreactive channel accounts for less than 20% of quenching outcomes. Complementary classical trajectory calculations of the postquenching dynamics are initiated from representative points along seams of conical intersections between the ground and excited-state potentials of OD(A (2)Σ(+),X (2)Π) + H(2). Diabatic modeling of the initial momenta in the dynamical calculations captures the key experimental trends: OD?X (2)Π products released primarily in their ground vibrational state with extensive rotational excitation and a branching ratio that strongly favors reactive quenching. The OD?A (2)Σ(+) + H(2) results are also compared with previous studies on the quenching of OH?A (2)Σ(+) + H(2); the two experimental studies show remarkably similar rotational energy distributions for the OH and OD?X (2)Π radical products.     

Full dimensional time-dependent quantum dynamics study of the H + NH3 --> H2 + NH2 reaction

A rigorous full dimensional time-dependent wave packet method has been developed for the reactive scattering between an atom and a tetra-atomic molecule. The method has been applied to the hydrogen abstraction reaction H+NH(3)-->H(2)+NH(2). Initial state-selected total reaction probabilities are investigated for the reactions from the ground vibrational state and from four excited vibrational states of ammonia. The total reaction probabilities from two lowest "tunneling doublets" due to the inversion barrier for the umbrella bending motion of NH(3) and from two pairs of doubly degenerate vibrational states of NH(3) are also inspected. Integral cross sections and rate constants are calculated for the reaction from the ground state with the centrifugal-sudden approximation. The calculated results are compared with those from the previous seven dimensional calculations [M. Yang and J. C. Corchado, J. Chem. Phys. 126, 214312 (2007)]. This work shows that the full dimensional rate constants are a factor of 3 larger than the corresponding seven dimensional calculated values at T=200 K and are overall smaller than those obtained from the variational transition state theory in the whole temperature region. The work also reveals that nonreactive NH bonds of NH(3) cannot be treated as spectators due to the fact that three NH bonds are coupled with each other during the reaction process.     

A theoretical and experimental study on translational and internal energies of H2O and OH from the 157 nm irradiation of amorphous solid water at 90 K

The photodesorption of H(2)O in its vibrational ground state, and of OH radicals in their ground and first excited vibrational states, following 157 nm photoexcitation of amorphous solid water has been studied using molecular dynamics simulations and detected experimentally by resonance-enhanced multiphoton ionization techniques. There is good agreement between the simulated and measured energy distributions. In addition, signals of H(+) and OH(+) were detected in the experiments. These are inferred to originate from vibrationally excited H(2)O molecules that are ejected from the surface by two distinct mechanisms: a direct desorption mechanism and desorption induced by secondary recombination of photoproducts at the ice surface. This is the first reported experimental evidence of photodesorption of vibrationally excited H(2)O molecules from water ice.     

A new ab initio potential energy surface for studying vibrational relaxation in NO(v) + NO collisions

A new ab initio potential energy surface for the ground state of the NO-NO system has been calculated within a reduced dimensionality model. We find an unusually large vibrational dependence of the interaction potential which explains previous spectroscopic observations. The potential can be used to model vibrational energy transfer, and here we perform quantum scattering calculations of the vibrational relaxation of NO(v). We show that the vibrational relaxation for v = 1 is 4 orders of magnitude larger than that for the related O(2)(v) + O(2) system without having to invoke nonadiabatic mechanisms as had been suggested in the past. For highly vibrationally excited states, we predict a strong dependence of the rates on the vibrational quantum number as has been observed experimentally, although there remain important quantitative differences. The importance of a chemically bound isomer on the relaxation mechanism is analyzed, and we conclude it does not play a role for the values of v considered in the experiment. Finally, the intriguing negative temperature dependence of the vibrational relaxation rate constants observed in experiments was studied using an statistical model to include the presence of many asymptotically degenerate spin-orbit states.     

3D Generalized langevin equation approach to gas–surface reactive scattering: model H+H→H2/Si(100)-(2×1)

Classical trajectory calculation has been performed for the H+H→H₂/Si(100)-(2×1) reaction by the 3D Generalized Langevin Equation (GLE) approach. The implementation of the 3D GLE approach to the H+H→H₂/Si(100)-(2×1) reaction is presented. Reaction probabilities are calculated for given surface temperatures and given collision energies. We also calculated vibrational and rotational distributions of product H₂ molecules from the reaction. About 80% of the product hydrogen molecules are in the ground vibrational state and the remaining 20% of the products are in the excited states. The rotational state shows non-Boltzmann distribution which can be seen in the direct collision process. Vibrational and Rotational distributions are strongly related to the impact parameter. The vibrational distribution is correlated with the x-component of the impact parameter and reflects the dimer nature of the silicon (2×1) surface. Details of the dynamics involving vibrational and rotational transitions are discussed.     

Doubly excited 2 1Delta g state of Na2

The doubly excited valence (3p+3p) 2 (1)Delta(g) state of Na(2) is experimentally observed by using optical-optical double resonance spectroscopy. A single line Ar(+) laser (a total of nine lines) was used to pump the sodium dimers from thermally populated ground state X (1)Sigma(g) (+) to the intermediate B (1)Pi(u) state. Then, a single mode Ti:sapphire laser was used to probe the doubly excited 2 (1)Delta(g) state. Violet fluorescence emitted from the highly excited states (mainly 2 (3)Pi(g) or 3 (3)Pi(g) states which are transferred from 2 (1)Delta(g) state via collision) to the a (3)Sigma(u) (+) state was monitored by a filtered photomultiplier tube (PMT). A total of 582 rovibrational levels of 2 (1)Delta(g) state were observed, identified, and assigned to the vibrational and rotational quantum numbers in the range of 0< or =v< or =28 and 11< or =J< or =99, respectively. The absolute vibrational quantum number assignment was verified by comparing the totally resolved fluorescence with the calculated Franck-Condon factors between 2 (1)Delta(g) state and B (1)Pi(u) state. Dunham coefficients and Rydberg-Klein-Rees potential curve were derived from these observed quantum levels. The primary molecular constants of Na(2) 2 (1)Delta(g) state are T(e)=32 416.759(15) cm(-1), omega(e)=124.8484(36) cm(-1), B(e)=0.119 158(3) cm(-1), and R(e)=3.508 20(5) A.     

Photoionization and photodissociation dynamics of the B 1sigmau + and C 1Piu states of H2 and D2

The photoionization and photodissociation dynamics of H(2) and D(2) in selected rovibrational levels of the B (1)Sigma(u) (+) and C (1)Pi(u) states have been investigated by velocity map ion imaging. The selected rotational levels of the B (1)Sigma(u) (+) and C (1)Pi(u) states are prepared by three-photon excitation from the ground state. The absorption of fourth photon results in photoionization to produce H(2)(+) X (2)Sigma(g)(+) or photodissociation to produce a ground-state H(1s) atom and an excited H atom with n >or= 2. The H(2) (+) ion can be photodissociated by absorption of a fifth photon. The resulting H(+) or D(+) ion images provide information on the vibrational state dependence of the photodissociation angular distribution of the molecular ion. The excited H(n >or= 2) atoms produced by the neutral dissociation process can also be ionized by the absorption of a fifth photon. The resulting ion images provide insight into the excited state branching ratios and angular distributions of the neutral photodissociation process. While the experimental ion images contain information on both the ionic and neutral processes, these can be separated based on constraints imposed on the fragment translational energies. The angular distribution of the rings in the ion images indicates that the neutral dissociation of molecular hydrogen and its isotopes is quite complex, and involves coupling to both doubly excited electronic states and the dissociation continua of singly excited Rydberg states.     

Spectroscopy above the ionization threshold: Dissociative recombination of the ground vibrational level of N(2) (+)

Comprehensive theoretical calculations are reported for the dissociative recombination of the lowest vibrational level of the N(2) (+) ground state. Fourteen dissociative channels, 21 electron capture channels, and 48 Rydberg series including Rydberg states having the first excited state of the ion as core are described for electron energies up to 1.0 eV. The calculation of potential curves, electron capture and predissociation widths, cross sections and rate constants are described. The cross sections and rate constants are calculated using Multichannel Quantum Defect Theory which allows for efficient handling of the Rydberg series. The most important dissociative channel is 2(3)Π(u) followed by 4(3)Π(u). Dissociative states that do not cross the ion within the ground vibrational level turning points play a significant role in determining the cross section structure and at isolated energies can be more important than states having a favorable crossing. By accounting for autoionization, the interactions between resonances, between dissociative states, and between resonances and dissociative states it is found that the cross section can be viewed as a complex dissociative recombination spectrum in which resonances overlap and interfere. The detailed cross section exhibits a rapid variation in atomic quantum yields for small changes in the electron energy. A study of this rapid variation by future high resolution storage ring experiments is suggested. A least squares fit to the calculated rate constant from the ground vibrational level is 2.2+0.2-0.4×10(-7)×(T(e)/300)(-0.40)?cm(3)/sec for electron temperatures, T(e), between 100 and 3000 K and is in excellent agreement with experimentally derived values.     

The Al+ -H(2) cation complex: rotationally resolved infrared spectrum, potential energy surface, and rovibrational calculations

The infrared spectrum of the Al(+)-H(2) complex is recorded in the H-H stretch region (4075-4110 cm(-1)) by monitoring Al(+) photofragments. The H-H stretch band is centered at 4095.2 cm(-1), a shift of -66.0 cm(-1) from the Q(1)(0) transition of the free H(2) molecule. Altogether, 47 rovibrational transitions belonging to the parallel K(a)=0-0 and 1-1 subbands were identified and fitted using a Watson A-reduced Hamiltonian, yielding effective spectroscopic constants. The results suggest that Al(+)-H(2) has a T-shaped equilibrium configuration with the Al(+) ion attached to a slightly perturbed H(2) molecule, but that large-amplitude intermolecular vibrational motions significantly influence the rotational constants derived from an asymmetric rotor analysis. The vibrationally averaged intermolecular separation in the ground vibrational state is estimated as 3.03 A, decreasing by 0.03 A when the H(2) subunit is vibrationally excited. A three-dimensional potential energy surface for Al(+)-H(2) is calculated ab initio using the coupled cluster CCSD(T) method and employed for variational calculations of the rovibrational energy levels and wave functions. Effective dissociation energies for Al(+)-H(2)(para) and Al(+)-H(2)(ortho) are predicted, respectively, to be 469.4 and 506.4 cm(-1), in good agreement with previous measurements. The calculations reproduce the experimental H-H stretch frequency to within 3.75 cm(-1), and the calculated B and C rotational constants to within approximately 2%. Agreement between experiment and theory supports both the accuracy of the ab initio potential energy surface and the interpretation of the measured spectrum.     

Calculations on the rate of the ion-molecule reaction between NH3+ and H2

The rate coefficient for the ion-molecule reaction NH3(+) + H2 --> NH4(+) + H has been calculated as a function of temperature with the use of the statistical phase space approach. The potential surface and reaction complex and transition state parameters used in the calculation have been taken from ab initio quantum chemical calculations. The calculated rate coefficient has been found to mimic the unusual temperature dependence measured in the laboratory, in which the rate coefficient decreases with decreasing temperature until 50-100 K and then increases at still lower temperatures. Quantitative agreement between experimental and theoretical rate coefficients is satisfactory given the uncertainties in the ab initio results and in the dynamics calculations. The rate coefficient for the unusual three-body process NH3(+) + H2 + He --> NH4(+) + H + He has also been calculated as a function of temperature and the result found to agree well with a previous laboratory determination.     

Radiative association of He+ with H2 at temperatures below 100 K

The paper presents a theoretical study of the low-energy dynamics of radiative association processes in the He+ + H2 collision system. Formation of the triatomic HeH2(+) ion in its bound rotation-vibration states on the potential-energy surfaces of the ground and of the first excited electronic states is investigated. Close-coupling calculations are performed to determine detailed state-to-state characteristics (bound <-- free transition rates, radiative and dissociative widths of resonances) as well as temperature-average characteristics (rate constants, photon emission spectra) of the two-state (X <-- A) reaction He+(2S) + H2(X1sigma(g)+) --> HeH2(+)(X2A') + h nu and of the single-state (A <-- A) reaction He+(2S) + H2(X1sigma(g)+) --> HeH2(+)(A2A') + h nu. The potential-energy surfaces of the X- and A-electronic states of HeH2(+) and the dipole moment surfaces determined ab initio in an earlier work [Kraemer, Spirko, and Bludsky, Chem. Phys. 276, 225 (2002)] are used in the calculations. The rate constants k(T) as functions of temperature are calculated for the temperature interval 1 < or = T < or = 100 K. The maximum k(T) values are predicted as 3.3 x 10(-15) s(-1) cm3 for the X <-- A reaction and 2.3 x 10(-20) s(-1) cm3 for the A <-- A reaction at temperatures around 2 K. Rotationally predissociating states of the He+-H2 complex, correlating with the upsilon = 0, j = 2 state of free H2, are found to play a crucial role in the dynamics of the association reactions at low temperatures; their contribution to the k(T) function of the X <-- A reaction at T < 30 K is estimated as larger than 80%. The calculated partial rate constants and emission spectra show that in the X <-- A reaction the HeH2(+)(X) ion is formed in its highly excited vibrational states. This is in contrast with the vibrational state population of the ion when formed via the (X <-- X) reaction He(1S) + H2(+)(X2sigma(g)+) --> HeH2(+)(X2A') + h nu.     

Gas phase electronic spectrum of linear AlCCH

The electronic spectrum of the aluminium containing species AlCCH has been detected in the gas phase in the region 315-355 nm. The experiment used a mass selective resonant two-color two-photon ionization technique coupled to a laser ablation source. Structures of the AlCCH isomers have been optimized using density functional theory (DFT) and the excitation energies to the low-lying electronic excited states calculated. Based on the analysis of the observed rotational structure and the theoretical data, the spectrum is assigned to the A (1)Pi<-- X (1)Sigma(+) electronic transition of linear AlCCH. The vibronic band system is complicated by the Renner-Teller effect in the excited state. The assignment yields nu(4)' = 516.4 cm(-1) for the stretching mode in the ground X (1)Sigma(+) state and nu(4)' = 654.5 cm(-1) for A (1)Pi excited state. Molecular constants determined from the rotational analysis are B(0)' = 0.16487(14), B(0)' = 0.17845(13) and T(0) = 28 755.04 cm(-1). The experimental and theoretical data indicate a shorter Al-C bond in the A (1)Pi excited than the X (1)Sigma(+) ground state.     

A 3000 K laboratory emission spectrum of water

An emission spectrum of hot water with a temperature of about 3000 K is obtained using an oxy-acetylene torch. This spectrum contains a very large number of transitions. The spectrum, along with previous cooler laboratory emission spectra and an absorption spectrum recorded from a sunspot, is analyzed in the 500-2000 cm(-1) region. Use of a calculated variational linelist for water allows significant progress to be made on assigning transitions involving highly excited vibrational and rotational states. In particular emission from rotationally excited states up to J=42 and vibrational levels with up to eight quanta of bending motion are assigned.     

The dissociation dynamics of He...I 35Cl(B,v'=2,3) complexes with varying amounts of internal energy

The He...I (35)Cl intermolecular vibrational levels with n'=0-6 that are bound within the He+ICl(B,v'=3) potential [A. B. McCoy, J. P. Darr, D. S. Boucher, P. R. Winter, M. D. Bradke, and R. A. Loomis, J. Chem. Phys. 120, 2677 (2004)] are identified in laser-induced fluorescence experiments performed at very low temperatures within a supersonic expansion. Comparisons of the positions and intensities of these lines with the excitation spectra, calculated using potential surfaces to describe the interactions between the helium atom and ICl in its ground and excited state, assist in the assignments. Based on these comparisons the excited state potential was rescaled so that the experimental and calculated J'=0 energies agree to within the experimental uncertainties for all but the lowest, n'=0, intermolecular level. Two-laser, action, and pump-probe spectroscopy experiments indicate that the bound He...I (35)Cl(B,v'=3) intermolecular vibrational levels undergo vibrational predissociation forming rotationally excited I (35)Cl(B,v'=2,j') products with distributions that depend upon the initial intermolecular vibrational level excited. Action spectra recorded in the ICl B-X, 2-0 region while monitoring the Deltav=0, I (35)Cl(B,v'=2) channel reveal two additional dissociation mechanisms for the He...I (35)Cl(B,v') excited state complexes: rotational predissociation of discrete metastable states lying slightly above the He+I (35)Cl(B,v'=2) asymptote and direct dissociation that occurs when the linear conformer is excited to the continuum of states above the same asymptote. The rotational predissociation pathway forms I (35)Cl(B,v'=2,j') products in all of the rotational states energetically accessible. The direct dissociation mechanism yields very cold rotational product state distributions; for instance, the average rotational energy in the product state distribution measured when the linear complexes are prepared 20 cm(-1) above the dissociation limit is only 1.51 cm(-1), representing only 7.6% of the available energy.     

Spectroscopy and metastability of CO2 2+ molecular ions

The spectroscopy and metastability of the carbon dioxide doubly charged ion, the CO(2) (2+) dication, have been studied with photoionization experiments: time-of-flight photoelectron photoelectron coincidence (TOF-PEPECO), threshold photoelectrons coincidence (TPEsCO), and threshold photoelectrons and ion coincidence (TPEsCO ion coincidence) spectroscopies. Vibrational structure is observed in TOF-PEPECO and TPEsCO spectra of the ground and first two excited states. The vibrational structure is dominated by the symmetric stretch except in the TPEsCO spectrum of the ground state where an antisymmetric stretch progression is observed. All three vibrational frequencies are deduced for the ground state and symmetric stretch and bending frequencies are deduced for the first two excited states. Some vibrational structure of higher electronic states is also observed. The threshold for double ionization of carbon dioxide is reported as 37.340+/-0.010 eV. The fragmentation of energy selected CO(2) (2+) ions has been investigated with TPEsCO ion coincidence spectroscopy. A band of metastable states from approximately 38.7 to approximately 41 eV above the ground state of neutral CO(2) has been observed in the experimental time window of approximately 0.1-2.3 mus with a tendency towards shorter lifetimes at higher energies. It is proposed that the metastability is due to slow spin forbidden conversion from bound excited singlet states to unbound continuum states of the triplet ground state. Another result of this investigation is the observation of CO(+)+O(+) formation in indirect dissociative double photoionization below the threshold for formation of CO(2) (2+). The threshold for CO(+)+O(+) formation is found to be 35.56+/-0.10 eV or lower, which is more than 2 eV lower than previous measurements.     

Bimolecular reactions of vibrationally excited molecules. Roaming atom mechanism at low kinetic energies

Quasiclassical trajectory calculations have been performed for the H + H'X(v) → X + HH' abstraction and H + H'X(v) → XH + H' (X = Cl, F) exchange reactions of the vibrationally excited diatomic reactant at a wide collision energy range extending to ultracold temperatures. Vibrational excitation of the reactant increases the abstraction cross sections significantly. If the vibrational excitation is larger than the height of the potential barrier for reaction, the reactive cross sections diverge at very low collision energies, similarly to capture reactions. The divergence is quenched by rotational excitation but returns if the reactant rotates fast. The thermal rate coefficients for vibrationally excited reactants are very large, approach or exceed the gas kinetic limit because of the capture-type divergence at low collision energies. The Arrhenius activation energies assume small negative values at and below room temperature, if the vibrational quantum number is larger than 1 for HCl and larger than 3 for HF. The exchange reaction also exhibits capture-type divergence, but the rate coefficients are larger. Comparisons are presented between classical and quantum mechanical results at low collision energies. At low collision energies the importance of the exchange reaction is enhanced by a roaming atom mechanism, namely, collisions leading to H atom exchange but bypassing the exchange barrier. Such collisions probably have a large role under ultracold conditions. The calculations indicate that for roaming to occur, long-range attractive interaction and small relative kinetic energy in the chemical reaction at the first encounter are necessary, which ensures that the partners can not leave the attractive well. Large orbital angular momentum of the primary products (equivalent to large rotational excitation in a unimolecular reaction) is favorable for roaming.     

Calculation of rovibrational energy states for water vapor using the symmetric top approximation in the Pade form

Considering the water vapor molecule for ka > or = J/2 > 1 as a symmetric top, the simple analytical Pade form for vibrational rotational energy levels was proposed. The rotational spectroscopic parameters for ground and 010 vibrational states of H2O molecule in the framework of this model have been obtained with good prediction. The absorption coefficient of water vapor in the spectral range 725-925 cm(-1) has been calculated for 1000-6000 K temperatures taking into account the transitions to high-excited states up to J < or = 35 and J < or = 30 for the 000 and 010 vibrational states, respectively. It is shown that the prominent role of transitions to high-excited states is in 8-12 microm atmospheric transparency window.     

Ab initio many-electron study for the low-lying states of the alkali hydride cations in the adiabatic representation

An ab initio multireference single- and double-excitation configuration interaction (CI) study is carried out for the ground and excited electronic states of alkali-hydride cations (LiH(+), NaH(+), KH(+), RbH(+), and CsH(+)). For all alkali-metal atoms, the first inner-shell and valence electrons (nine active electrons, three for Li) are considered explicitly in the ab initio self-consistent-field and CI calculations. The adiabatic potential energy curves, radial and rotational couplings are calculated and presented. Short-range (～3 a.u.) potential wells produced by the excitation of the inner-shell electrons are found. The depths of the inner potential wells are much greater than those of the outer wells for the CsH(+) system. The computed spectroscopic constants for the long-range potential well of the 2 (2)Σ(+) state are very close to the available theoretical and experimental data. The electronic states of alkali-hydrogen cations are also compared with each other, it is found that the positions of the potential wells shift to larger internuclear distances gradually, and the depths of these potential wells become greater with increasing alkali-metal atomic number. The relationships between structures of the radial coupling matrix elements and the avoiding crossings of the potential curves are analyzed. From NaH(+) to CsH(+), radial coupling matrix elements display more and more complex structures due to the gradual decrease of energy separations for avoided crossings. Finally, the behavior of some rotational couplings is also shown.