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.     

Full-dimensional quantum dynamics of vibrationally highly excited NHD2

We report on full-dimensional vibrational quantum dynamics of the highly excited ammonia isotopologue NHD(2) using a newly developed potential energy surface and the MCTDH program package. The calculations allow to realistically simulate an infrared laser induced stereomutation reaction at the pyramidal nitrogen atom in the femtosecond time domain. Our results allow for a thorough qualitative and quantitative understanding of infrared photoinduced stereomutation kinetics, the underlying quantum dynamics, and the reaction mechanisms. Comparison is made with a previous, reduced dimensionality study of the same reaction [R. Marquardt, M. Quack, I. Thanopulos, and D. Luckhaus, J. Chem. Phys. 118, 643 (2003)], and it is shown that slight variances of reduced spaces lead to significantly different kinetics. Because the quantum dynamics depends subtly on variances of reduced spaces, reduced dimensionality treatments are not reliable even for qualitative predictions of the stereomutation kinetics. The first direct comparison between the Multiconfigurational Time Dependent Hartree [M. H. Beck, A. Ja?ckle, G. A. Worth et al., Phys. Rep. 324, 1 (2000)] and Unimolecular Reactions Induced by Monochromatic Infrared Radiation [M. Quack and E. Sutcliffe, QCPE Bulletin 6, 98 (1986)] program packages on a specific, four dimensional quantum dynamical problem allows for their full validation in the present work.     

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.     

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.     

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H(5)O(2) (+) using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh et al. [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.     

Quantum dynamics of rovibrational transitions in H2-H2 collisions: internal energy and rotational angular momentum conservation effects

We present a full dimensional quantum mechanical treatment of collisions between two H(2) molecules over a wide range of energies. Elastic and state-to-state inelastic cross sections for ortho-H(2)?+ para-H(2) and ortho-H(2)?+ ortho-H(2) collisions have been computed for different initial rovibrational levels of the molecules. For rovibrationally excited molecules, it has been found that state-to-state transitions are highly specific. Inelastic collisions that conserve the total rotational angular momentum of the diatoms and that involve small changes in the internal energy are found to be highly efficient. The effectiveness of these quasiresonant processes increases with decreasing collision energy and they become highly state-selective at ultracold temperatures. They are found to be more dominant for rotational energy exchange than for vibrational transitions. For non-reactive collisions between ortho- and para-H(2) molecules for which rotational energy exchange is forbidden, the quasiresonant mechanism involves a purely vibrational energy transfer albeit with less efficiency. When inelastic collisions are dominated by a quasiresonant transition calculations using a reduced basis set involving only the quasiresonant channels yield nearly identical results as the full basis set calculation leading to dramatic savings in computational cost.     

Anatomy of the S(1D) + H2 reaction: the dynamics on two new potential energy surfaces from quantum dynamics calculations

We present exact quantum integral and differential cross sections for the title reaction from a time-dependent wavepacket method which takes account of all Coriolis couplings. We employ two new potential energy surfaces fitted using the double many-body expansion (DMBE) method. The difference between the two surfaces is that for the first the data was extrapolated to the complete basis set limit (CBS) and for the second the data was corrected semi-empirically (SEC). While the DMBE/CBS surface is, on first impressions, regarded as the most accurate, our results show that this surface gives consistent smaller cross section when compared to previous results employing an earlier surface, named Ho after its first author. We also find that the DMBE/CBS surface features an unphysical barrier for contracted H(2) distances which explains the smaller results. The DMBE/SEC surface, which is based on the same data, does not show the same barrier and the results compare much better to previous theoretical results as well as those from experiment. While we find that overall the differential cross sections from the DMBE/SEC surface are forward scattered, which is in line with experiment, the cross sections do not rise steeply enough with decreasing energy showing that this surface is not sufficiently attractive at low energies. We find this is due to a shallow van der Waals well present for the Ho surface but not on the DMBE surfaces.     

Coherence parameters for He+(2p) capture and H(2p) excitation in He2+-H(1s) collisions

Semiclassical coupled channel calculations have been carried out for the collision system He²⁺-H(1s) in the velocity range 0.15–3.0 a.u. (impact energies 0.5–225 keV/amu) in order to study capture probabilities and alignment and orientation parameters for the dominant He⁺(n=2) channels. A 14-state AO basis set calculation has been combined with an analytical treatment of the asymptotic collision region. For impact velocities about and abovev=0.6 a.u. a strong propensity for resonance capture into an oriented He⁺(2p) state with the same sense of rotation as the collisional rotation of the internuclear axis is predicted together with a very smooth behaviour of the alignment angle as function of impact parameter. Eikonal method calculations of differential capture cross sections predict that the left/right orientation asymmetry will prevail in differential scattering experiments. The resulting total cross sections for capture into specificnl-substates (n=2, 3) and the total light polarisation parameter for He⁺(2p) capture compare well with previous work. Finally we report H(2s,2p) excitation cross sections, probabilties and H(2p) alignment and orientation parameters, following the established propensity rule for orientation in H(2p) excitation.     

Reduced dimensionality quantum calculations of resonances in H + H2

Reduced dimensionality quantum reaction probabilities are reported for the H+H₂ reaction using the ab initio potential surface of Liu, Siegbahn, Truhlar and Horowitz. Resonances are found for the ground and first two excited adiabatic bending states of H₃. Comparison of the resonance energies with the new coupled states calculations of Colton and Schatz shows good agreement. Additional resonances are reported for energies greater than those considered in the coupled states studies.     

Six-dimensional quantum dynamics of dissociative chemisorption of H2 on Ru(0001)

Six-dimensional quantum dynamics calculations on dissociative chemisorption of H(2) on Ru(0001) are performed. The six-dimensional potential energy surface is generated using density functional theory. Two different generalized gradient approximations are used, i.e., RPBE and PW91, to allow the results to be compared. The dissociation probability for normally incident H(2) on a clean Ru(0001) surface is calculated. Large differences between the reaction probabilities calculated using the RPBE and PW91 are seen, with the PW91 results showing a much narrower reaction probability curve and a much higher reactivity. Using the reaction probabilities and assuming normal energy scaling reaction rates are generated for temperatures between 300 and 800 K. The rate generated using the PW91 results is higher by about a factor 5 than the rate based on the RPBE results in the range of temperatures relevant to ammonia production.     

Full-dimensional quantum calculations of vibrational spectra of six-atom molecules. I. Theory and numerical results

Two quantum mechanical Hamiltonians have been derived in orthogonal polyspherical coordinates, which can be formed by Jacobi and/or Radau vectors etc., for the study of the vibrational spectra of six-atom molecules. The Hamiltonians are expressed in an explicit Hermitian form in the spatial representation. Their matrix representations are described in both full discrete variable representation (DVR) and mixed DVR/nondirect product finite basis representation (FBR) bases. The two-layer Lanczos iteration algorithm [H.-G. Yu, J. Chem. Phys. 117, 8190 (2002)] is employed to solve the eigenvalue problem of the system. A strategy regarding how to carry out the Hamiltonian-vector products for a high-dimensional problem is discussed. By exploiting the inversion symmetry of molecules, a unitary sequential 1D matrix-vector multiplication algorithm is proposed to perform the action of the Hamiltonian on the wavefunction in a symmetrically adapted DVR or FBR basis in the azimuthal angular variables. An application to the vibrational energy levels of the molecular hydrogen trimer (H2)3 in full dimension (12D) is presented. Results show that the rigid-H2 approximation can underestimate the binding energy of the trimer by 27%. Finally, it is demonstrated that the two-layer Lanczos algorithm is also capable of computing the eigenvectors of the system with minor effort.     

Renner-Teller quantum dynamics of the N(2D) + H2-->NH + H reaction

We present the Born-Oppenheimer (BO) and Renner-Teller (RT) quantum dynamics of the reaction (14)N((2)D)+(1)H(2)(X (1)Sigma(g) (+))-->NH(X (3)Sigma(-))+H((2)S), considering the NH(2) electronic states X (2)B(1) and A (2)A(1). These states correlate to the same (2)Pi(u) linear species, are coupled by RT nonadiabatic effects, and give NH(X (3)Sigma(-))+H and NH(a (1)Delta)+H, respectively. We develop the Hamiltonian matrix elements in the R embedding of the Jacobi coordinates and in the adiabatic electronic representation, using the permutation-inversion symmetry, and taking into account the nuclear-spin statistics. Collision observables are calculated via the real wave-packet (WP) and flux methods, using the potential-energy surfaces of Santoro et al. [J. Phys. Chem. A 106, 8276 (2002)]. WP snapshots show that the reaction proceeds via an insertion mechanism, and that the RT-WP avoids the A (2)A(1) potential barrier, jumping from the excited to the ground surface and giving mainly the NH(X (3)Sigma(-)) products. X (2)B(1) BO probabilities and cross sections show large tunnel effects and are approximately four to ten times larger than the A (2)A(1) ones. This implies a BO rate-constant ratio k(X (2)B(1))k(A (2)A(1)) approximately 10(5) at 300 K, i.e., a negligible BO formation of NH(a (1)Delta). When H(2) is rotationally excited, RT couplings reduce slightly the X (2)B(1) reaction observables, but enhance strongly the A (2)A(1) reactivity. These couplings are important at all collision energies, reduce the collision threshold, and increase remarkably reaction probabilities and cross sections. The RT k(A (2)A(1)) is thus approximately 3.3 order of magnitude larger than the BO value, and degeneracy-averaged, initial-state-resolved rate constants increase by approximately 13% and by approximately 47% at 300 and 500 K, respectively. Owing to an overestimation of the X (2)B(1) potential barrier, the calculated thermal rate is too low with respect to that observed, but we obtain a good agreement by shifting down the calculated cross section.     

Vibrationally inelastic collisions in H+ +CO system: comparing quantum calculations with experiments

State-resolved cross beam experiments [H. Udseth et al., J. Chem. Phys. 60, 3051 (1974); J. Krutein and F. Linder, J. Chem. Phys. 71, 599 (1979); G. Niedner-Schatteburg and J. P. Toennies, Adv. Chem. Phys. LXXXII, 553 (1992)], coupled with proton energy loss spectroscopy for the inelastic scattering of H(+) from CO in the collision range of 10-30 eV show very low vibrational excitation of the target molecule. Stimulated by the experimentally observed low vibrational inelasticity in the system the ground and the first two low-lying excited electronic potential-energy surfaces have been computed using the ab initio multireference configuration interaction method. Quantum dynamics has been performed on the ground potential energy surface in the framework of vibrational close-coupling rotational infinite-order sudden approximation. The various computed dynamical attributes such as differential and integral cross sections, and average vibrational energy transfer are analyzed in detail, and compared successfully with the available experimental results.     

Three dimensional quantum dynamics of (H-, H2) and its isotopic variants

We present the results of a time-dependent quantum mechanical investigation using centrifugal sudden approximation in the form of reaction probability as a function of collision energy (E(trans)) in the range 0.3-3.0 eV for a range of total angular momentum (J) values and the excitation function sigma(E(trans)) for the exchange reaction H(-) + H(2) (v = 0, j = 0) --> H(2) + H(-) and its isotopic variants in three dimensions on an accurate ab initio potential energy surface published recently (J. Chem. Phys. 2004, 121, 9343). The excitation function results are shown to be in excellent agreement with those obtained from crossed beam measurements by Zimmer and Linder for H(-) + D(2) collisions for energies below the threshold for electron detachment channel and somewhat larger than the most recent results of Haufler et al. for (H(-), D(2)) and (D(-), H(2)) collisions.     

Off-normal incidence dissociative sticking of H2 on Cu(100) studied using six-dimensional quantum calculations

Six-dimensional quantum calculations of the sticking probability for H2 hitting a Cu(100) surface with off-normal incidence are presented. The multiconfiguration time-dependent Hartree approach is employed for an efficient wave-packet propagation. The sticking probability is calculated for different initial momenta parallel to the surface. In contrast with the picture described in the literature, the sticking probability was found to depend on the parallel momentum. The results are explained by the topology of the potential-energy surface, which shows significant corrugation with a moderate variation of the barrier height with the surface site.     

A quantum wave packet dynamics study of the N(2D) + H2 reaction

We report a dynamics study of the reaction N((2)D) + H(2) (v=0, j=0-5) --> NH + H using the time-dependent quantum wave packet method and a recently reported single-sheeted double many-body expansion potential energy surface for NH(2)(1(2)A' ') which has been modeled from accurate ab initio multireference configuration-interaction calculations. The calculated probabilities for (v=0, j=0-5) are shown to display resonance structures, a feature also visible to some extent in the calculated total cross sections for (v=0, j=0). A comparison between the calculated centrifugal-sudden and coupled-channel reaction probabilities validate the former approximation for the title system. Rate constants calculated using a uniform J-shifting scheme and averaged over a Boltzmann distribution of rotational states are shown to be in good agreement with the available experimental values. Comparisons with other theoretical results are also made.     

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.     

Full-dimensional quantum dynamics study of vinylidene-acetylene isomerization: a scheme using the normal mode Hamiltonian

Full-dimensional quantum dynamics calculations of vinylidene-acetylene isomerization are performed and the state-specific resonance decay lifetimes of vinylidene(-d(2)) are computed. The theoretical scheme is a combination of several methods: normal coordinates are chosen to describe the nuclear motion of vinylidene, with both the parity and permutation symmetry exploited; phase space optimization in combination with physical considerations is used to generate an efficient discrete variable representation; the reaction coordinate is defined by us according to the three most relevant normal coordinates, along which a kind of optimal complex absorbing potential is imposed; the preconditioned inexact spectral transform method combined with an efficient preconditioner is employed to extract the energies and lifetimes of vinylidene. The overall computation is efficient. The computed energy levels generally agree with experiment well, and several state-specific lifetimes are reported for the first time.     

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcalmol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively.     

Exploring the dynamics of reaction C(3P) + C2H4 with crossed beam/photoionization experiments and quantum chemical calculations

We investigated the title reaction at collision energy 3.5 kcal mol(-1) in a crossed molecular beam apparatus using undulator radiation as an ionization source. Time-of-flight (TOF) spectra of product C(3)H(3) were measured in laboratory angles from 20° to 100° using two photoionization energies 9.5 and 11.6 eV. These two sets of experimental data exhibit almost the same TOF distributions and laboratory angular distributions. From the best simulation, seven angle-specific kinetic-energy distributions and a nearly isotropic angular distribution are derived for product channel C(3)H(3) + H that has an average kinetic-energy release of 15.5 kcal mol(-1), corresponding to an average internal energy of 33.3 kcal mol(-1) in C(3)H(3). Furthermore, TOF spectra of product C(3)H(3) were measured at laboratory angle 52° with ionizing photon energies from 7 to 12 eV. The appearance of TOF spectra remains almost the same, indicating that a species exclusively contributes to product C(3)H(3); the species is identified as H(2)CCCH (propargyl) based on the ionization energy of 8.6 ± 0.2 eV and the maximal kinetic-energy release of 49 kcal mol(-1). Theoretical calculations indicate that the rapid inversion mechanism and rotation in intermediate H(2)CCCH(2) can result in a forward-backward symmetric angular distribution for product C(3)H(3) + H. The present work avoids the interference of reactions of C((1)D) and C(2) radicals with C(2)H(4) and rules out the probability of production of other isomers like c-C(3)H(3) and H(3)CCC proposed in the previous work at least at the investigated collision energy.