Full-dimensional quantum dynamics study of the H + H2O and H + HOD exchange reactions

The initial state-selected time-dependent wave packet approach is employed to study the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H, HOH' + D exchange reactions with both OH bonds in the H(2)O reactant and OH(D) bond in the HOD reactant treated as reactive bonds. The total reaction probabilities for different partial waves, as well as the integral cross sections, which are the exact CC (coupled-channel) results, are first obtained in this study for the H(2)O(HOD) reactant initially in the ground rovibrational state. Because of the shallow C(3v) minimum along the reaction path, the reaction probabilities for the three reactions present several resonance peaks, with one dominant resonance peak just above the threshold. The cross sections for the H' + HOD → HOH' + D reaction are substantially smaller than those for the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H reactions, indicating that the H'/H exchange reactions are much more favored. In the CC calculations, the resonance peaks in the reaction probabilities diminish quickly with the increase in total angular momenta J, resulting in the existence of a clear step-like feature just above the threshold in the cross sections for the title reactions, which manifests the signature of shape resonances in these reactions. In the CS calculations, the resonance peaks on reaction probabilities persist in many partial waves, and thus the resonance structures can no longer survive the partial-wave summation and are washed out completely in the CS cross sections for the title reactions.     

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.     

H2O(HOD) + Cl→HCl(DCl) + OH(OD)反应动态学的理论研究

用从头算方法, 获得了H2O + Cl→HCl + OH(R1), HOD +Cl→DCl + OH(R2), HOD + Cl→HCl + OD(R3)反应的内禀反应坐标(IRC)。根据传统过渡态、变分过渡态理论及相应的隧道效应校正, 计算了反应的速率常数。对已有实验速率常数值的R1反应, 我们计算的结果和实验一致。根据Truhlar的振动选态公式, 分别讨论了激发HOD中OH, OD振动模式对反应速率的影响,得到激发HOD中的OH振动模式将有利于产物OD + HCl生成, 和实验的结论相一致。     

Experimental and theoretical studies of charge transfer and deuterium ion transfer between D2O+ and C2H4

The charge transfer and deuterium ion transfer reactions between D(2)O(+) and C(2)H(4) have been studied using the crossed beam technique at relative collision energies below one electron volt and by density functional theory (DFT) calculations. Both direct and rearrangement charge transfer processes are observed, forming C(2)H(4) (+) and C(2)H(3)D(+), respectively. Independent of collision energy, deuterium ion transfer accounts for approximately 20% of the reactive collisions. Between 22 and 36 % of charge transfer collisions occur with rearrangement. In both charge transfer processes, comparison of the internal energy distributions of products with the photoelectron spectrum of C(2)H(4) shows that Franck-Condon factors determine energy disposal in these channels. DFT calculations provide evidence for transient intermediates that undergo H/D migration with rearrangement, but with minimal modification of the product energy distributions determined by long range electron transfer. The cross section for charge transfer with rearrangement is approximately 10(3) larger than predicted from the Rice-Ramsperger-Kassel-Marcus isomerization rate in transient complexes, suggesting a nonstatistical mechanism for H/D exchange. DFT calculations suggest that reactive trajectories for deuterium ion transfer follow a pathway in which a deuterium atom from D(2)O(+) approaches the pi-cloud of ethylene along the perpendicular bisector of the C-C bond. The product kinetic energy distributions exhibit structure consistent with vibrational motion of the D-atom in the bridged C(2)H(4)D(+) product perpendicular to the C-C bond. The reaction quantitatively transforms the reaction exothermicity into internal excitation of the products, consistent with mixed energy release in which the deuterium ion is transferred in a configuration in which both the breaking and the forming bonds are extended.     

The mechanism of proton exchange: guided ion beam studies of the reactions, H(H2O)n+ (n=1-4)+D2O and D(D2O)n+ (n=1-4)+H2O

Reactions of protonated water clusters, H(H(2)O)(n) (+) (n=1-4) with D(2)O and their "mirror" reactions, D(D(2)O)(n) (+) (n=1-4) with H(2)O, are studied using guided-ion beam mass spectrometry. Absolute reaction cross sections are determined as a function of collision energy from thermal energy to over 10 eV. At low collision energies, we observe reactions in which H(2)O and D(2)O molecules are interchanged and reactions where H-D exchange has occurred. As the collision energy is increased, the H-D exchange products decrease and the water exchange products become dominant. At high collision energies, processes in which one or more water molecules are lost from the reactant ions become important, with simple collision-induced dissociation processes, i.e., those without H-D exchange, being dominant. Threshold energies of endothermic channels are measured and used to determine binding energies of the proton bound complexes, which are consistent with those determined by thermal equilibrium measurements and previous collision-induced dissociation studies. A kinetic scheme that relies only on the ratio of isomerization and dissociation rate constants successfully accounts for the kinetic energy dependence observed in the branching ratios for H-D and water exchange products in all systems. Rice-Ramsperger-Kassel-Marcus theory and ab initio calculations confirm the feasibility and establish the details of this kinetic model.     

Exact quantum calculations of the kinetic isotope effect: cross sections and rate constants for the F+HD reaction and role of tunneling

In this paper we present integral cross sections (in the 5-220 meV collision energy range) and rate constants (in the 100-300 K range of temperature) for the F+HD reaction leading to HF+D and DF+H. The exact quantum reactive scattering calculations were carried out using the hyperquantization algorithm on an improved potential energy surface which incorporates the effects of open shell and fine structure of the fluorine atom in the entrance channel. The results reproduce satisfactorily molecular beam scattering experiments as well as chemical kinetics data for both the HF and DF channels. In particular, the agreement of the rate coefficients and the vibrational branching ratios with experimental measurements is improved with respect to previous studies. At thermal and subthermal energies, the rates are greatly influenced by tunneling through the reaction barrier. Therefore exchange of deuterium is shown to be penalized with respect to exchange of hydrogen, and the isotopic branching exhibits a strong dependence on translational energy. Also, it is found that rotational excitation of the reactant HD molecule enhances the production of HF and decreases the reactivity at the D end, obtaining insight on the reaction stereodynamics.     

The Reaction X + Cl2→XCl + Cl (X = Mu,H, D). II. Comparison of experimental data with theoretical results derived from a new potential energy surface

We consider experimental implications for the Mu + Cl₂, H + Cl₂, and D + Cl₂ reactions of the extended London—Eyring—Polanyi—Sato (LEPS) potential energy surface derived from experimental data in paper I. In the present calculations, it is necessary to make additional implicit and explicit assumptions concerning the three-dimensional (3D) nature of the potential surface, since the inversion procedure of paper I yields information only on the collinear (1D) part of the surface. We have performed accurate 1D quantum calculations of reaction probabilities, which are then transformed into 3D by an information theoretic 1D → 3D transformation incorporating a constraint to allow for angular momentum transfer effects in light+heavy—heavy atom reactions. This procedure implicitly accounts for the 3D nature of the potential surface. The calculated vibrational and vibrotational product distributions are in good agreement with those determined in thermal chemiluminescence experiments. The Sato parameters for the 1D surface also define a full 3D surface. This is used as an approximation to the true surface, and its properties are explored in 3D quasiclassical trajectory calculations. Comparison is made for the H and D reactions with available chemiluminescence, molecular beam and kinetic experimental data for differential and total reaction cross sections, energy disposal, rate coefficients and Arrhenius parameters. Some kinetic isotope effects in the Mu, H, and D reactions are discussed using vibrationally adiabatic theory. Comparison is also made with results from other calculations in the literature for the H + Cl₂ and D + Cl₂ reactions.     

Planar study of H2O+Cl<-->HO+HCl reactions

The first four dimensional (4D) quantum scattering calculations on the tetra-atomic H2O+Cl<-->HO+HCl reactions are reported. With respect to a full (6D) treatment, only the planar constraint and a fixed length for the HO spectator bond are imposed. This work explicitly accounts for the bending and local HO stretching vibrations in H2O, for the vibration of HCl and for the in-plane rotation of the H2O, HO and HCl molecules. The calculations are performed with the potential energy surface of Clary et al. and use a Born-Oppenheimer type separation between the motions of the light and the heavy nuclei. State-to-state cross sections are reported for a collision energy range 0-1.8 eV measured with respect to H2O+Cl. For the H2O+Cl reaction, present results agree with previous (3D) non planar calculations and confirm that excitation of the H2O stretching promotes more reactivity than excitation of the bending. New results are related to the rotation of the H2O molecule: the cross sections are maximal for planar rotational states corresponding to 10相似文献   

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

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

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.     

Reaction of HOD+ with NO2: effects of OD and OH stretching, bending, and collision energy on reactions on the singlet and triplet potential surfaces

Integral cross sections and product recoil velocity distributions were measured for the reaction of HOD(+) with NO(2), in which the HOD(+) reactant was prepared in its ground state and with mode-selective excitation in the 001 (OH stretch), 100 (OD stretch), and 010 (bend) modes. In addition, we measured the 300 K thermal kinetics in a selected ion flow tube reactor and report product branching ratios different from previous measurements. Reaction is found to occur on both the singlet and triplet surfaces with near-unit efficiency. At 300 K, the product branching indicates that triplet → singlet transitions occur in about 60% of triplet-coupled collisions, which we attribute to long interaction times mediated by complexes on the triplet surface. Because the collision times are much shorter in the beam experiments, the product distributions show no signs of such transitions. The dominant product on the singlet surface is charge transfer. Reactions on the triplet surface lead to NO(+), NO(2)H(+), and NO(2)D(+). There is also charge transfer, producing NO(2)(+) (a(3)B(2)); however, this triplet NO(2)(+) mostly predissociates. The NO(2)H(+)/NO(2)D(+) cross sections peak at low collision energies and are insignificant above ~1 eV due to OH/OD loss from the nascent product ions. The effects of HOD(+) vibration are mode-specific. Vibration inhibits charge transfer, with the largest effect from the bend. The NO(2)H(+)/NO(2)D(+) channels are also vibrationally inhibited, and the mode dependence reveals how energy in different reactant modes couples to the internal energy of the product ions.     

Experimental and theoretical study of the reactions between vanadium oxide cluster cations and water

Vanadium oxide cluster cations V(x)O(y)(+) (x = 2-6) are prepared by laser ablation and are reacted with D(2)O in a fast flow reactor under room temperature conditions. A time-of-flight mass spectrometer is used to detect the cluster distribution before and after the reactions. Observation of the products (V(2)O(5))(1-3)D(+) indicates the deuterium atom abstraction reaction (V(2)O(5))(1-3)(+) + D(2)O → (V(2)O(5))(1-3)D(+) + OD. In addition, significant association products (V(2)O(5))(1-3)D(2)O(+) are also observed in the experiments. Density functional theory calculations are performed to study the reaction mechanisms of V(4)O(10)(+) with H(2)O. The calculated results are in agreement with the experimental observations and indicate that H(2)O is dissociatively rather than molecularly adsorbed in V(4)O(10)H(2)O(+) complex.     

Isotope branching and tunneling in O(3P)+HD-->OH+D; OD+H reactions

The O((3)P)+HD and O((3)P)+D(2) reactions are studied using quantum scattering calculations and chemically accurate potential energy surfaces developed for the O((3)P)+H(2) system by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)]. Cross sections and rate coefficients for OH and OD products are calculated using accurate quantum methods as well as the J-shifting approximation. The J-shifting approach is found to work remarkably well for both O+HD and O+D(2) collisions. The reactions are dominated by tunneling at low temperatures and for the O+HD reaction the hydrogen atom transfer leading to the OH product dominates at low temperatures. Our result for the OH/OD branching ratio is in close agreement with previous calculations over a wide range of temperatures. The computed OH/OD branching ratios are also in close agreement with experimental results of Robie et al. [Chem. Phys. Lett. 134, 579 (1987)] at temperatures above 400 K but the theoretical results do not reproduce the rapid rise in the experimental values of the branching ratio for temperatures lower than 350 K. We believe that new measurements could resolve the long-standing discrepancy between experiment and theory for this benchmark reaction.     

Kinetics of Al + H2O reaction: theoretical study

Quantum chemical calculations were carried out to study the reaction of Al atom in the ground electronic state with H(2)O molecule. Examination of the potential energy surface revealed that the Al + H(2)O → AlO + H(2) reaction must be treated as a complex process involving two steps: Al + H(2)O → AlOH + H and AlOH + H → AlO + H(2). Activation barriers for these elementary reaction channels were calculated at B3LYP/6-311+G(3df,2p), CBS-QB3, and G3 levels of theory, and appropriate rate constants were estimated by using a canonical variational theory. Theoretical analysis exhibited that the rate constant for the Al + H(2)O → products reaction measured by McClean et al. must be associated with the Al + H(2)O → AlOH + H reaction path only. The process of direct HAlOH formation was found to be negligible at a pressure smaller than 100 atm.     

Integral and state-to-state cross sections for the reaction D + H2(νf) → HD(νf) + H: A quantum mechanical study within the infinite order sudden approximation

In this paper are presented quantum mechanical t-initial and t-average cross sections and rate constants for the reactions D + H₂(ν₁ = 0, 1) → HD(ν_f = 0, 1) + H. The calculations were done employing the infinite order sudden approximation. It was found that the t-average total cross sections overlap very nicely with the available classical cross sections. As for rate constants a reasonably good fit was found with available experimental results.     

Reaction rates for O + HD → OH + D and O + HD → OD + H

We have calculated reaction rates for the reactions O + HD → OH + D and O + DH → OD + H using improved canonical variational transition state theory and least-action ground-state transmission coefficients with an ab initio potential energy surface. The kinetic isotope effects are in good agreement with experiment. The optimized tunneling paths and properties of the variational transition states and the rate enhancement for vibrationally excited reactants are also presented and compared with those for the isotopically unsubstituted reaction O + H₂ → OH + H. The thermal reactions at low and room temperature are predicted to occur by tunneling at extended configurations, i.e., to initiate early on the reaction path and to avoid the saddle point regions. Tunneling also dominates the low and room temperature reactions for excited vibrational states, but in these cases the results are not as sensitive to the nature of the tunneling path. Overbarrier mechanisms dominate for both thermal and excited-vibrational state reactions for T > 600 K. For the excited-state reaction (with initial vibrational quantum number n > 0) a transition state switch occurs for T > 1000 K for the O + HD(n = 1) → OD + H case and for T > 1500 K for the O + DH(n = 1) → OD + H reaction, and this may be a general phenomenon for excited-state reactions at higher temperature. In the present case the switch occurs from an early variational transition state where the vibrationally adiabatic approximation is expected to be valid to a tighter variational transition state where nonadiabatic effects are probably important and should be included.     

Multidimensional separations of ubiquitin conformers in the gas phase: Relating ion cross sections to H/D exchange measurements

Investigating gas-phase structures of protein ions can lead to an improved understanding of intramolecular forces that play an important role in protein folding. Both hydrogen/deuterium (H/D) exchange and ion mobility spectrometry provide insight into the structures and stabilities of different gas-phase conformers, but how best to relate the results from these two methods has been hotly debated. Here, high-field asymmetric waveform ion mobility spectrometry (FAIMS) is combined with Fourier-transform ion cyclotron resonance mass spectrometry (FT/ICR MS) and is used to directly relate ubiquitin ion cross sections and H/D exchange extents. Multiple conformers can be identified using both methods. For the 9+ charge state of ubiquitin, two conformers (or unresolved populations of conformers) that have cross sections differing by 10% are resolved by FAIMS, but only one conformer is apparent using H/D exchange at short times. For the 12+ charge state, two conformers (or conformer populations) have cross sections differing by <1%, yet H/D exchange of these conformers differ significantly (6 versus 25 exchanges). These and other results show that ubiquitin ion collisional cross sections and H/D exchange distributions are not strongly correlated and that factors other than surface accessibility appear to play a significant role in determining rates and extents of H/D exchange. Conformers that are not resolved by one method could be resolved by the other, indicating that these two methods are highly complementary and that more conformations can be resolved with this combination of methods than by either method alone.     

H+ versus D+ transfer from HOD+ to CO2: bond-selective chemistry and the anomalous effect of bending excitation

Reactions of HOD(+) with CO(2) have been studied for HOD(+) in its ground state, and with one quantum of excitation in each of its vibrational modes: (001)--predominantly OH stretch, 0.396 eV; (010)--bend, 0.153 eV; and (100)--predominantly OD stretch, 0.293 eV. Integral cross sections and product recoil velocities were recorded for collision energies from threshold to 3 eV. The cross sections for both H(+) and D(+) transfer rise with increasing collision energy from threshold to ～1 eV, then become weakly dependent of the collision energy. All three vibrational modes enhance the total reactivity, but quite mode specifically. The H(+) transfer reaction is enhanced by OH stretch excitation, whereas OD stretch excitation has little effect. Conversely, the D(+) transfer reaction is enhanced by OD stretch excitation, while the OH stretch has little effect. Excitation of the bend strongly enhances both channels. The effects of the stretch excitations are consistent with previous studies of neutral HOD mode-selective chemistry, and can be at least qualitatively understood in terms of a late barrier to product formation. The fact that bend excitation produces the largest overall enhancement is surprising, because this is the lowest energy excitation, and is not obviously connected with the reaction coordinates for either H(+) or D(+) transfer. A rationalization in terms of the effects of water distortion on the potential surface is proposed.     

Highly efficient D2 generation by dehydrogenation of formic acid in D2O through H+/D+ exchange on an iridium catalyst: application to the synthesis of deuterated compounds by transfer deuterogenation

Deuterated compounds have received increasing attention in both academia and industrial fields. However, preparations of these compounds are limited for both economic and practical reasons. Herein, convenient generation of deuterium gas (D(2)) and the preparation of deuterated compounds on a laboratory scale are demonstrated by using a half-sandwich iridium complex with 4,4'-dihydroxy-2,2'-bipyridine. The "umpolung" (i.e., reversal of polarity) of a hydrogen atom of water was achieved in consecutive reactions, that is, a cationic H(+)/D(+) exchange reaction and anionic hydride or deuteride transfer, under mild conditions. Selective D(2) evolution (purity up to 89?%) was achieved by using HCO(2)H as an electron source and D(2)O as a deuterium source; a rhodium analogue provided HD gas (98?%) under similar conditions. Furthermore, pressurized D(2) (98?%) without CO gas was generated by using DCO(2)D in D(2)O in a glass autoclave. Transfer deuterogenation of ketones gave α-deuterated alcohols with almost quantitative yields and high deuterium content by using HCO(2)H in D(2)O. Mechanistic studies show that the H(+)/D(+) exchange reaction in the iridium hydride complex was much faster than β-elimination and hydride (deuteride) transfer.     

State-to-state reactive scattering in six dimensions using reactant-product decoupling: OH + H2 → H2O + H (J = 0)

We extend to full dimensionality a recently developed wave packet method [M. T. Cvitas? and S. C. Althorpe, J. Phys. Chem. A 113, 4557 (2009)] for computing the state-to-state quantum dynamics of AB + CD → ABC + D reactions and also increase the computational efficiency of the method. This is done by introducing a new set of product coordinates, by applying the Crank-Nicholson approximation to the angular kinetic energy part of the split-operator propagator and by using a symmetry-adapted basis-to-grid transformation to evaluate integrals over the potential energy surface. The newly extended method is tested on the benchmark OH + H(2) → H(2)O + H reaction, where it allows us to obtain accurately converged state-to-state reaction probabilities (on the Wu-Schatz-Fang-Lendvay-Harding potential energy surface) with modest computational effort. These methodological advances will make possible efficient calculations of state-to-state differential cross sections on this system in the near future.