Effect of reagent rotation on isotopic branching in (He, HD+) collisions

A three-dimensional time-dependent quantum mechanical wave packet approach is used to calculate reaction probability (P(R)) and integral reaction cross section (sigma(R)) values for both the channels of the reaction He + HD(+) (v = 1; j = 0, 1, 2, 3) --> HeH(D)(+) + D(H), over a range of translational energy (E(trans)) on the McLaughlin-Thompson-Joseph-Sathyamurthy (MTJS) potential energy surface using centrifugal sudden approximation for nonzero total angular momentum (J) values. The reaction probability plots as a function of translational energy for different J values exhibit several oscillations, which are characteristic of the system. It is shown that HeH(+) is preferred over HeD(+) for large J values and that HeD(+) is preferred over HeH(+) for small J values for all the rotational (j) states studied. The integral reaction cross section for both the channels and therefore the isotopic branching ratio for the reaction depend strongly on j in contrast to the marginal dependence shown by earlier QCT calculations. The computed results are in overall agreement with the available experimental results.     

Isotopic branching in (He, HD+) collisions

A three-dimensional time-dependent quantum mechanical approach is used to calculate the reaction probability (P(R)) and the integral reaction cross section (sigma(R)) for both channels of the reaction He + HD+(v = 0, 1, 2, 3; j = 0) --> HeH(D)+ + D(H), over a range of translational energy (E(trans)) on two different ab initio potential energy surfaces (McLaughlin-Thompson-Joseph-Sathyamurthy and Palmieri et al.). The reaction probability plots as a function of translational energy exhibit several oscillations, which are characteristic of the system. The vibrational enhancement of the reaction probability and the integral reaction cross section values are reproduced qualitatively by our calculations, in accordance with the experimental results. The isotopic branching ratio for the reaction decreases in going from v = 0 to v = 1 and then becomes nearly v-independent in going from v = 1 to v =3 on both the surfaces.     

Influence of reagent rotation on (H-, D2) and (D-, H2) collisions: a quantum mechanical study

Time-independent quantum mechanical (TIQM) approach (helicity basis truncated at k = 2) has been used for computing differential and integral cross sections for the exchange reaction H- + D2 (v = 0, j = 0-4) --> HD + D- and D- + H2 (v = 0, j = 0-3) --> HD + H- in three dimensions on an accurate ab initio potential energy surface. It is shown that the j-weighted differential reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder at four different relative translational energies (Etrans = 0.55, 0.93, 1.16 and 1.48 eV) for (H-, D2) and at one relative translational energy (Etrans = 0.6 eV) by Haufler et al. for both (H-, D2) and (D-, H2) collisions. The j-weighted integral reaction cross section values are in good agreement with the crossed beam measurements by Zimmer and Linder in the Etrans range 0.5-1.5 eV and close to the guided ion beam results by Haufler et al. for (H-, D2) in the range 0.8-1.2 eV. Time-dependent quantum mechanical (TDQM) results obtained using centrifugal sudden approximation are reported in the form of integral reaction cross section values as a function of Etrans in the range 0.3-3.0 eV for both reactions in three dimensions on the same potential energy surface. The TDQM reaction cross section values decline more sharply than the TIQM results with increase in the initial rotational quantum number (j) for the D2 molecules in their ground vibrational state (v = 0) for (H-, D2) collisions. The computed j-weighted reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder for (H-, D2) collisions and guided ion beam results by Haufler et al. for both (H-, D2) and (D-, H2) collisions for energies below the threshold for electron detachment channel.     

Exact quantum scattering study of the H + HS reaction on a new ab initio potential energy surface H2S (3A")

We present an exact quantum dynamical study and quasi-classical trajectory (QCT) calculations for the exchange and abstraction processes for the H + HS reaction. These calculations were based on a newly constructed high-quality potential energy surface for the lowest triplet state of H(2)S ((3)A"). The ab initio single-point energies were computed using complete active space self-consistent field and multi-reference configuration interaction method with a basis set of aug-cc-pV5Z. The time-dependent wave packet (TDWP) method was used to calculate the total reaction probabilities and integral cross sections over the collision energy (E(col)) range of 0.0-2.0 eV for the reactant HS initially at the ground state and the first vibrationally excited state. It was found that the initial vibrational excitation of HS enhances both abstraction and exchange processes. In addition, a good agreement is found between QCT and TDWP reaction probabilities at the total momentum J = 0 as a function of collision energy for the H + HS (v = 0, j = 0) reaction.     

Time-dependent quantum mechanical wave packet study of the He+H(2) (+)(v,j)-->HeH(+)+H reaction

A detailed three-dimensional time-dependent quantum dynamical study of the He+H(2) (+)(v=0-3,j=0)-->HeH(+)+H reaction is reported for different vibrational v states of H(2) (+) in its ground rotational (j=0) state over a range of translational E(trans) energies on an accurate ab initio potential energy surface published by Palmieri et al. Plots of reaction probability as a function of total energy E reveal a large number of oscillations indicating the presence of a number of reactive scattering resonances. When averaged over total angular momentum J, some of the oscillations survive, indicating that they may be amenable to experimental observation. A comparison of our present results with our earlier results on the McLaughlin-Thompson-Joseph-Sathyamurthy surface and the experimental results from different research groups reveal a good deal of agreement as well as some discrepancies between theory and experiment at the level of state-selected gas phase dynamics.     

Quenching of highly vibrationally excited pyrimidine by collisions with CO2

Relaxation of highly vibrationally excited pyrimidine (C(4)N(2)H(4)) by collisions with carbon dioxide has been investigated using diode laser transient absorption spectroscopy. Vibrationally hot pyrimidine (E(')=40 635 cm(-1)) was prepared by 248-nm excimer laser excitation, followed by rapid radiationless relaxation to the ground electronic state. The nascent rotational population distribution (J=58-80) of the 00(0)0 ground state of CO(2) resulting from collisions with hot pyrimidine was probed at short times following the excimer laser pulse. Doppler spectroscopy was used to measure the CO(2) recoil velocity distribution for J=58-80 of the 00(0)0 state. Rate constants and probabilities for collisions populating these CO(2) rotational states were determined. The measured energy transfer probabilities, indexed by final bath state, were resorted as a function of DeltaE to create the energy transfer distribution function, P(E,E(')), from E(')-E approximately 1300-7000 cm(-1). P(E,E(')) is fitted to a single exponential and a biexponential function to determine the average energy transferred in a single collision between pyrimidine and CO(2) and parameters that can be compared to previously studied systems using this technique, pyrazineCO(2), C(6)F(6)CO(2), and methylpyrazineCO(2). P(E,E(')) parameters for these four systems are also compared to various molecular properties of the donor molecules. Finally, P(E,E(')) is analyzed in the context of two models, one which suggests that the shape of P(E,E(')) is primarily determined by the low-frequency out-of-plane donor vibrational modes and one which suggests that the shape of P(E,E(')) can be determined by how the donor molecule final density of states changes with DeltaE.     

A time-dependent wave packet quantum scattering study of the reaction HD+ (v = 0 - 3;j0 = 1) + He --> HeH+(HeD+) + D(H)

Time-dependent wave packet quantum scattering (TWQS) calculations are presented for HD(+) (v = 0 - 3;j(0)=1) + He collisions in the center-of-mass collision energy (E(T)) range of 0.0-2.0 eV. The present TWQS approach accounts for Coriolis coupling and uses the ab initio potential energy surface of Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. For a fixed total angular momentum J, the energy dependence of reaction probabilities exhibits quantum resonance structure. The resonances are more pronounced for low J values and for the HeH(+) + D channel than for the HeD(+) + H channel and are particularly prominent near threshold. The quantum effects are no longer discernable in the integral cross sections, which compare closely to quasiclassical trajectory calculations conducted on the same potential energy surface. The integral cross sections also compare well to recent state-selected experimental values over the same reactant and translational energy range. Classical impulsive dynamics and steric arguments can account for the significant isotope effect in favor of the deuteron transfer channel observed for HD(+)(v<3) and low translational energies. At higher reactant energies, angular momentum constraints favor the proton-transfer channel, and isotopic differences in the integral cross sections are no longer significant. The integral cross sections as well as the J dependence of partial cross sections exhibit a significant alignment effect in favor of collisions with the HD(+) rotational angular momentum vector perpendicular to the Jacobi R coordinate. This effect is most pronounced for the proton-transfer channel at low vibrational and translational energies.     

Observation of partial wave structures in the integral cross section of the S((1)D2) + H2(j = 0) reaction

The integral cross section of the S((1)D(2)) + H(2)(j = 0) → SH + H reaction has been measured for the first time at collision energies from 0.820 down to 0.078 kJ mol(-1) in a high-resolution crossed beam experiment. The excitation function obtained exhibits a non-monotonic variation with collision energy and compares well with the results of high-level quantum calculations. In particular, the structures observed in the lower energy part, where only a few partial waves contribute, can be described in terms of the sequential opening of individual channels, consistent with the theoretical calculations.     

An experimental and quasiclassical trajectory study of the rovibrationally state-selected reactions: HD+(v=0-15,j=1)+He-->HeH+(HeD+)+D(H)

The absolute integral cross sections for the formation of HeH+ and HeD+ from the collisions of HD+(v,j=1)+He have been examined over a broad range of vibrational energy levels v=0-13 at the center-of-mass collision energies (ET) of 0.6 and 1.4 eV using the vacuum ultraviolet (VUV) pulsed field ionization photoelectron secondary ion coincidence method. The ET dependencies of the integral cross sections for products HeH+ and HeD+ from HD+(v=0-4)+He collisions in the ET range of 0-3 eV have also been measured using the VUV photoionization guided ion beam mass spectrometric technique, in which vibrationally selected HD+(v) reactant ions were prepared via excitation of selected autoionization resonances of HD. At low total energies, a pronounced isotope effect is observed in absolute integral cross sections for the HeH++D and HeD++H channels with significant favoring of the deuteron transfer channel. As v is increased in the range of v=0-9, the integral cross sections of the HeH++D channel are found to approach those of HeD++H. The observed velocity distributions of products HeD+ and HeH+ are consistent with an impulsive or spectator-stripping mechanism. Detailed quasiclassical trajectory (QCT) calculations are also presented for HD+(v,j=1)+He collisions at the same energies of the experiment. The QCT calculations were performed on the most accurate ab initio potential energy surface available. If the zero-point energy of the reaction products is taken into account, the QCT cross sections for products HeH+ and HeD+ from HD+(v)+He are found to be significantly lower than the experimental results at ET values near the reaction thresholds. The agreement between the experimental and QCT cross sections improves with translational energy. Except for prethreshold reactivity, QCT calculations ignoring the zero-point energy in the products are generally in good agreement with experimental absolute cross sections. The experimental HeH+/HeD+ branching ratios for the HD+(v=0-9)+He collisions are generally consistent with QCT predictions. The observed isotope effects can be rationalized on the basis of differences in thermochemical thresholds and angular momentum conservation constraints.     

Full-dimensional time-dependent wave packet dynamics of H2 + D2 reaction

Collision induced dissociation (CID), four center reaction (4C), and single exchange reaction (SE) in H(2) (v(1) = high) + D(2) (v(2) = low) were studied by means of time-dependent wave packet approach within a full-dimensional model. Initial state-selected total reaction probabilities for the three competitive processes have been computed on two realistic global potential energy surfaces of Aguado-Suárez-Paniagua and Boothroyd-Martin-Keogh-Peterson (BMKP) with the total angular momentum J = 0. The role of both vibrationally excited and rotationally excited reagents was examined by varying the initial vibrational and rotational states. The vibrational excitation of the hot diatom gives an enhancement effect on the CID process, while the vibrational excitation of the cold diatom gives an inhibition effect. The rotational excitation of both reagents has a significant effect on the reaction process. The 4C and SE probabilities are at least one order of magnitude smaller than the CID probabilities over the energy range considered. Isotope substitution effects were also studied by substituting the collider D(2) by H(2) and HD on the BMKP potential energy surfaces. The CID process is most efficient for the H(2) + D(2) combination and least efficient for the H(2) + H(2) combination and is different for the 4C and SE processes.     

Cross sections of the O+ + H2 --> OH+ + H ion-molecule reaction and isotopic variants (D2, HD): quasiclassical trajectory study and comparison with experiments

A dynamics study [cross section and microscopic mechanism versus collision energy (E(T))] of the reaction O+ + H2 --> OH+ + H, which plays an important role in Earth's ionosphere and interstellar chemistry, was conducted using the quasiclassical trajectory method, employing an analytical potential energy surface (PES) recently derived by our group [R. Martinez et al., J. Chem. Phys. 120, 4705 (2004)]. Experimental excitation functions for the title reaction, as well as its isotopic variants with D2 and HD, were near-quantitatively reproduced in the calculations in the very broad collision energy range explored (E(T) = 0.01-6.0 eV). Intramolecular and intermolecular isotopic effects were also examined, yielding data in good agreement with experimental results. The reaction occurs via two microscopic mechanisms (direct and nondirect abstraction). The results were satisfactorily interpreted based on the reaction probability and the maximum impact parameter dependences with E(T), and considering the influence of the collinear [OHH]+ absolute minimum of the PES on the evolution from reactants to products. The agreement between theory and experiment suggests that the reaction mainly occurs through the lowest energy PES and nonadiabatic processes are not very important in the wide collision energy range analyzed. Hence, the PES used to describe this reaction is suitable for both kinetics and dynamics studies.     

Dynamics of (H(-), H(2)) collisions: a time-dependent quantum mechanical investigation on a new ab initio potential energy surface

A global analytical potential energy surface for the ground state of H(3)(-) has been constructed by fitting an analytic function to the ab initio potential energy values computed using coupled cluster singles and doubles with perturbative triples [CCSD(T)] method and Dunning's augmented correlation consistent polarized valence triple zeta basis set. Using this potential energy surface, time-dependent quantum mechanical wave packet calculations were carried out to calculate the reaction probabilities (P(R)) for the exchange reaction H(-)+H(2)(v, j)-->H(2)+H(-), for different initial vibrational (v) and rotational (j) states of H(2), for total angular momentum equal to zero. With increase in v, the number of oscillations in the P(R)(E) plot increases and the oscillations become more pronounced. While P(R) increases with increase in rotational excitation from j=0 to 1, it decreases with further increase in j to 2 over a wide range of energies. In addition, rotational excitation quenches the oscillations in P(R)(E) plots.     

An investigation of nonadiabatic interactions in Cl(2Pj) + D2 via crossed-molecular-beam scattering

We have determined limits on the cross section for both electronically nonadiabatic excitation and quenching in the Cl((2)P(j)) + D(2) system. Our experiment incorporates crossed-molecular-beam scattering with state-selective Cl((2)P(12,32)) detection and velocity-mapped ion imaging. By colliding atomic chlorine with D(2), we address the propensity for collisions that result in a change of the spin-orbit level of atomic chlorine either through electronically nonadiabatic spin-orbit excitation Cl((2)P(32)) + D(2)-->Cl(*)((2)P(12)) + D(2) or through electronically nonadiabatic spin-orbit quenching Cl(*)((2)P(12)) + D(2)-->Cl((2)P(32)) + D(2). In the first part of this report, we estimate an upper limit for the electronically nonadiabatic spin-orbit excitation cross section at a collision energy of 5.3 kcal/mol, which lies above the energy of the reaction barrier (4.9 kcal/mol). Our analysis and simulation of the experimental data determine an upper limit for the excitation cross section as sigma(NA)< or =0.012 A(2). In the second part of this paper we investigate the propensity for electronically nonadiabatic spin-orbit quenching of Cl(*) following a collision with D(2) or He. We perform these experiments at collision energies above and below the energy of the reaction barrier. By comparing the amount of scattered Cl(*) in our images to the amount of Cl(*) lost from the atomic beam we obtain the maximum cross section for electronically nonadiabatic quenching as sigma(NA)< or =15(-15) (+44) A(2) for a collision energy of 7.6 kcal/mol. Our experiments show the probability for electronically nonadiabatic quenching in Cl(*) + D(2) to be indistinguishable to that for the kinematically identical system of Cl(*) + He.     

State-resolved collisional quenching of vibrationally excited pyrazine (E(vib) = 37,900 cm(-1)) by D35Cl(v = 0)

Supercollision relaxation of highly vibrationally excited pyrazine (E(vib) = 37,900 cm(-1)) with D35Cl is investigated using high-resolution transient IR diode laser absorption spectroscopy at 4.4 microm. Highly excited pyrazine is prepared by pulsed UV excitation at 266 nm, followed by rapid radiationless decay to the ground electronic state. The rotational energy distribution of the scattered DCl (v = 0,J) molecules with J = 15-21 is characterized by T(rot) = 755+/-90 K. The relative translational energy increases as a function of rotational quantum number for DCl with T(rel) = 710+/-190 K for J = 15 and T(rel) = 1270+/-240 K for J = 21. The average change in recoil velocity correlates with the change in rotational angular momentum quantum number and highlights the role of angular momentum in energy gain partitioning. The integrated energy-transfer rate for appearance of DCl (v = 0,J = 15-21) is k(2)(int) = 7.1x10(-11) cm3 molecule(-1) s(-1), approximately one-eighth the Lennard-Jones collision rate. The results are compared to earlier energy gain measurements of CO2 and H2O.     

Guided-ion beam and theoretical study of the potential energy surface for activation of methane by W+

A guided-ion beam tandem mass spectrometer is used to study the reactions, W(+) + CH(4) (CD(4)) and [W,C,2H](+) + H(2) (D(2)), to probe the [W,C,4H](+) potential energy surface. The reaction W(+) + CH(4) produces [W,C,2H](+) in the only low-energy process. The analogous reaction in the CD(4) system exhibits a cross section with strong differences at the lowest energies caused by zero-point energy differences, demonstrating that this reaction is slightly exothermic for CH(4) and slightly endothermic for CD(4). The [W,C,2H](+) product ion reacts further at thermal energies with CH(4) to produce W(CH(2))(x)(+) (x = 2-4). At higher energies, the W(+) + CH(4) reaction forms WH(+) as the dominant ionic product with smaller amounts of WCH(3)(+), WCH(+), and WC(+) also formed. The energy dependent cross sections for endothermic formation of the various products are analyzed and allow the determination of D(0)(W(+)-CH(3)) approximately 2.31 +/- 0.10 eV, D(0)(W(+)-CH(2)) = 4.74 +/- 0.03 eV, D(0)(W(+)-CH) = 6.01 +/- 0.28 eV, and D(0)(W(+)-C) = 4.96 +/- 0.22 eV. We also examine the reverse reaction, [W,C,2H](+) + H(2) (D(2)) --> W(+) + CH(4) (CH(2)D(2)). Combining the cross sections for the forward and reverse processes yields an equilibrium constant from which D(0)(W(+)-CH(2)) = 4.72 +/- 0.04 eV is derived. Theoretical calculations performed at the B3LYP/HW+/6-311++G(3df,3p) level yield thermochemistry in reasonable agreement with experiment. These calculations help identify the structures and electronic states of the species involved and characterize the potential energy surface for the [W,C,4H](+) system.     

A crossed molecular beam imaging study of the O(1D2)+HCl-->OH+Cl(2P J=3/2, 1/2) reaction.

A crossed molecular beam study is presented for the O((1)D(2))+HCl-->OH+Cl((2)P(J)) reaction at the collision energy of 6 kcal mol(-1). State-resolved doubly differential cross sections are obtained for the Cl((2)P(J=3/2) ) and Cl*((2)P(J=1/2) ) products by velocity-map ion imaging. Both products are slightly more forward scattered, which suggests a reaction mechanism without a long-lived intermediate in the ground electronic state. A small fraction (23 %) of the energy release into the translational degree of freedom indicates strong internal excitation of the counterpart OH radical. The contribution of the electronic excited states of O--HCl to the overall reaction is also examined from the doubly differential cross sections.     

Ultracold collisions and reactions of vibrationally excited OH radicals with oxygen atoms

We report a quantum dynamics study of O + OH (v = 1, j = 0) collisions on its ground electronic state, employing two different potential energy surfaces: the DIMKP surface by Kendrick and Pack, and the XXZLG surface by Xu et al. A time-independent quantum mechanical method based on hyperspherical coordinates has been adopted for the dynamics calculations. Energy-dependent probabilities and rate coefficients are computed for the elastic, inelastic, and reactive channels over the collision energy range E(coll) = 10(-10)-0.35 eV, for J = 0 total angular momentum. Initial state-selected reaction rate coefficients are also calculated from the J = 0 reaction probabilities by applying a J-shifting approximation, for temperatures in the range T = 10(-6)-700 K. Our results show that the dynamics of the collisional process and its outcome are strongly influenced by long-range forces, and chemical reactivity is found to be sensitive to the choice of the potential energy surface. For O + OH (v = 1, j = 0) collisions at low temperatures, vibrational relaxation of OH competes with reactive scattering. Since long-range interactions can facilitate vibrational relaxation processes, we find that the DIMKP potential (which explicitly includes van der Waals dispersion terms) favours vibrational relaxation over chemical reaction at low temperatures. On the DIMKP potential in the ultracold regime, the reaction rate coefficient for O + OH (v = 1, j = 0) is found to be a factor of thirteen lower than that for O + OH (v = 0, j = 0). This significantly high reactivity of OH (v = 0, j = 0), compared to that of OH (v = 1, j = 0), is attributed to enhancement caused by the presence of a HO(2) quasibound state (scattering resonance) with energy near the O + OH (v = 0, j = 0) dissociation threshold. In contrast, the XXZLG potential does not contain explicit van der Waals terms, being just an extrapolation by a nearly constant function at large O-OH distances. Therefore, long-range potential couplings are absent in calculations using the XXZLG surface, which does not induce vibrational relaxation as efficiently as the DIMKP potential. The XXZLG potential leads to a slightly higher reactivity (a factor of 1.4 higher) for O + OH (v = 1, j = 0) compared to that for O + OH (v = 0, j = 0) at ultracold temperatures. Overall, both potential surfaces yield comparable values of reaction rate coefficients at low temperatures for the O + OH (v = 1, j = 0) reaction.     

Quantum tunneling and quantum phase interference in a [Mn(II)2Mn(III)2] single-molecule magnet

[Mn4(hmp)6(H2O)2(NO3)2](NO3)2.2.5H2O (1) has been synthesized from the reaction of 2-hydroxymethylpyridine (Hhmp) with Mn(NO3)2.4H2O in the presence of tetraethylammonium hydroxide. 1 crystallizes in the triclinic P space group with two crystallographically independent centrosymmetrical [Mn4(hmp)6(H2O)2(NO3)2]2+ complexes in the packing structure. Four Mn ions are arranged in a double-cuboidal fashion where outer Mn2+ are heptacoordinated and inner Mn3+ are hexacoordinated. dc magnetic measurements show that both Mn2+...Mn3+ and Mn3+...Mn3+ interactions are ferromagnetic with J(wb)/k(B) = +0.80(5) K, and J(bb)/k(B) = +7.1(1) K, respectively, leading to an S(T) = 9 ground state. Combined ac and dc measurements reveal the single-molecule magnet (SMM) behavior of 1 with both thermally activated and ground-state tunneling regimes, including quantum phase interference. In the thermally activated regime, the characteristic relaxation time (tau) of the system follows an Arrhenius law with tau0 = 6.7 x 10(-)(9) s and delta(eff)/k(B) = 20.9 K. Below 0.34 K, tau saturates indicating that the quantum tunneling of the magnetization becomes the dominant relaxation process as expected for SMMs. Down to 0.04 K, field dependence of the magnetization measured using the mu-SQUID technique shows the presence of very weak inter-SMM interactions (zJ'/k(B) approximately -1.5 x 10(-3) K) and allows an estimation of D/k(B) at -0.35 K. Quantum phase interference has been used to confirm the D value and to estimate the transverse anisotropic parameter to E/k(B) = +0.083 K and the ground-state tunnel splitting delta(LZ) = 3 x 10(-7) K at H(trans) = 0 Oe. These results rationalize the observed tunneling time (tau(QTM)) and the effective energy barrier (delta(eff)).