Quantum dynamics of H2, D2, and HD in the small dodecahedral cage of clathrate hydrate: evaluating H2-water nanocage interaction potentials by comparison of theory with inelastic neutron scattering experiments

We have performed rigorous quantum five-dimensional (5D) calculations and analysis of the translation-rotation (T-R) energy levels of one H(2), D(2), and HD molecule inside the small dodecahedral (H(2)O)(20) cage of the structure II clathrate hydrate, which was treated as rigid. The H(2)- cage intermolecular potential energy surface (PES) used previously in the molecular dynamics simulations of the hydrogen hydrates [Alavi et al., J. Chem. Phys. 123, 024507 (2005)] was employed. This PES, denoted here as SPC/E, combines an effective, empirical water-water pair potential [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and electrostatic interactions between the partial charges placed on H(2)O and H(2). The 5D T-R eigenstates of HD were calculated also on another 5D H(2)-cage PES denoted PA-D, used by us earlier to investigate the quantum T-R dynamics of H(2) and D(2) in the small cage [Xu et al., J. Phys. Chem. B 110, 24806 (2006)]. In the PA-D PES, the hydrogen-water pair potential is described by the ab initio 5D PES of the isolated H(2)-H(2)O dimer. The quality of the SPC/E and the PA-D H(2)-cage PESs was tested by direct comparison of the T-R excitation energies calculated on them to the results of two recent inelastic neutron scattering (INS) studies of H(2) and HD inside the small clathrate cage. The translational fundamental and overtone excitations, as well as the triplet splittings of the j=0-->j=1 rotational transitions, of H(2) and HD in the small cage calculated on the SPC/E PES agree very well with the INS results and represent a significant improvement over the results computed on the PA-D PES. Our calculations on the SPC/E PES also make predictions about several spectroscopic observables for the encapsulated H(2), D(2), and HD, which have not been measured yet.     

Hydrogen adsorbed in a metal organic framework-5: coupled translation-rotation eigenstates from quantum five-dimensional calculations

We report rigorous quantum five-dimensional (5D) calculations of the coupled translation-rotation (T-R) eigenstates of a H(2) molecule adsorbed in metal organic framework-5 (MOF-5), a prototypical nanoporous material, which was treated as rigid. The anisotropic interactions between H(2) and MOF-5 were represented by the analytical 5D intermolecular potential energy surface (PES) used previously in the simulations of the thermodynamics of hydrogen sorption in this system [Belof et al., J. Phys. Chem. C 113, 9316 (2009)]. The global and local minima on this 5D PES correspond to all of the known binding sites of H(2) in MOF-5, three of which, α-, β-, and γ-sites are located on the inorganic cluster node of the framework, while two of them, the δ- and ε-sites, are on the phenylene link. In addition, 2D rotational PESs were calculated ab initio for each of these binding sites, keeping the center of mass of H(2) fixed at the respective equilibrium geometries; purely rotational energy levels of H(2) on these 2D PESs were computed by means of quantum 2D calculations. On the 5D PES, the three adjacent γ-sites lie just 1.1 meV above the minimum-energy α-site, and are separated from it by a very low barrier. These features allow extensive wave function delocalization of even the lowest translationally excited T-R eigenstates over the α- and γ-sites, presenting significant challenges for both the quantum bound-state calculations and the analysis of the results. Detailed comparison is made with the available experimental data.     

One and two hydrogen molecules in the large cage of the structure II clathrate hydrate: quantum translation-rotation dynamics close to the cage wall

We have performed a rigorous theoretical study of the quantum translation-rotation (T-R) dynamics of one and two H2 and D2 molecules confined inside the large hexakaidecahedral (5(12)6(4)) cage of the sII clathrate hydrate. For a single encapsulated H2 and D2 molecule, accurate quantum five-dimensional calculations of the T-R energy levels and wave functions are performed that include explicitly, as fully coupled, all three translational and the two rotational degrees of freedom of the hydrogen molecule, while the cage is taken to be rigid. In addition, the ground-state properties, energetics, and spatial distribution of one and two p-H2 and o-D2 molecules in the large cage are calculated rigorously using the diffusion Monte Carlo method. These calculations reveal that the low-energy T-R dynamics of hydrogen molecules in the large cage are qualitatively different from that inside the small cage, studied by us recently. This is caused by the following: (i) The large cage has a cavity whose diameter is about twice that of the small cage for the hydrogen molecule. (ii) In the small cage, the potential energy surface (PES) for H2 is essentially flat in the central region, while in the large cage the PES has a prominent maximum at the cage center, whose height exceeds the T-R zero-point energy of H2/D2. As a result, the guest molecule is excluded from the central part of the large cage, its wave function localized around the off-center global minimum. Peculiar quantum dynamics of the hydrogen molecule squeezed between the central maximum and the cage wall manifests in the excited T-R states whose energies and wave functions differ greatly from those for the small cage. Moreover, they are sensitive to the variations in the hydrogen-bonding topology, which modulate the corrugation of the cage wall.     

Hydrogen molecule in the small dodecahedral cage of a clathrate hydrate: quantum five-dimensional calculations of the coupled translation-rotation eigenstates

We report quantum five-dimensional (5D) calculations of the energy levels and wave functions of the hydrogen molecule, para-H2 and ortho-H2, confined inside the small dodecahedral (H2O)20 cage of the sII clathrate hydrate. All three translational and the two rotational degrees of freedom of H2 are included explicitly, as fully coupled, while the cage is treated as rigid. The 5D potential energy surface (PES) of the H2-cage system is pairwise additive, based on the high-quality ab initio 5D (rigid monomer) PES for the H2-H2O complex. The bound state calculations involve no dynamical approximations and provide an accurate picture of the quantum 5D translation-rotation dynamics of H2 inside the cage. The energy levels are assigned with translational (Cartesian) and rotational quantum numbers, based on calculated root-mean-square displacements and probability density plots. The translational modes exhibit negative anharmonicity. It is found that j is a good rotational quantum number, while the threefold degeneracy of the j = 1 level is lifted completely. There is considerable translation-rotation coupling, particularly for excited translational states.     

Hydrogen molecule in the small dodecahedral cage of a clathrate hydrate: quantum translation-rotation dynamics at higher excitation energies

Higher-lying five-dimensional translation-rotation (T-R) eigenstates of a single p-H2 and o-D2 molecule confined inside the small dodecahedral (512) cage of the structure II clathrate hydrate are calculated rigorously, as fully coupled, with the cage assumed to be rigid. The calculations cover the excitation energies up to and beyond the j=2 rotational level of the free molecule, 356 cm(-1) for H2 and 179 cm(-1) for D2. It is found that j is a good quantum number for all the T-R states of p-H2, j=0 and j=2, considered. The same is not true for o-D2, where a number of T-R states in the neighborhood of the j=2 level show significant mixing of j=0 and j=2 rotational basis functions. The 5-fold degeneracy of the j=2 level of p-H2 is lifted completely due to the anisotropy of the cage environment, as is the 3-fold degeneracy of the j=1 level of o-H2 studied by us previously. Pure translational mode excitations with up to four quanta display negative anharmonicity, which was observed earlier for the translational fundamentals and their first overtones. The issues of assigning the combination states of p-H2 with excitations of two or all three translational modes, and of the strength of the mode coupling as a function of the excitation energy, are studied carefully for a range of quantum numbers. The average T-R energy of the encapsulated p-H2 is calculated as a function of temperature from 0 to 150 K.     

Quantum dynamics of coupled translational and rotational motions of H2 inside C60

We report rigorous quantum calculations of the translation-rotation (T-R) eigenstates of the H2 molecule in C60. The resulting level structure can be explained in terms of a few dominant features. These include the coupling between the orbital and the rotational angular momenta of H2 to give the total angular momentum lambda, and the splitting of the sevenfold degeneracy of T-R levels with lambda=3 by the nonsphericity of C60, according to the rules of the icosahedral I h group.     

Nonadiabatic quantum reactive scattering of the OH(A 2Σ+) + D2

The seams of conical intersection exist between the ground (1 (2)A(')) and the first-excited (2 (2)A(')) electronic potential energy surfaces (PESs) of OH(A (2)Σ(+),X (2)Π) + H(2) system. This intersection induces the nonadiabatic quenching of OH(A (2)Σ(+)) by D(2). We present nonadiabatic quantum dynamics study for OH(A (2)Σ(+)) + D(2) on new five-dimensional coplanar PESs. The ab initio calculations of PESs are based on multireference configuration interaction (MRCI)/aug-cc-pVQZ level. A back-propagation neural network is utilized to fit the PESs and nonadiabatic coupling. High degrees of rotational excitation of quenched OH(X (2)Π) products are found in nonreactive quenching channel, and the quenched D(2) products are vibrationally excited up to quantum number v(2) (')=8. The theoretical results of nonadiabatic time-dependent wave-packet calculation are in good agreement with the existing experimental data.     

Time dependent quantum dynamics study of the O++H2(v=0,j=0)-->OH++H ion-molecule reaction and isotopic variants (D2,HD)

The time dependent real wave packet method using the helicity decoupling approximation was used to calculate the cross section evolution with collision energy (excitation function) of the O++H2(v=0,j=0)-->OH++H reaction and its isotopic variants with D2 and HD, using the best available ab initio analytical potential energy surface. The comparison of the calculated excitation functions with exact quantum results and experimental data showed that the present quantum dynamics approach is a very useful tool for the study of the selected and related systems, in a quite wide collision energy interval (approximately 0.0-1.1 eV), involving a much lower computational cost than the quantum exact methods and without a significant loss of accuracy in the cross sections.     

Mean potential statistical theory of the H+ + D2 → HD + D+ reaction

The reactive collision process H(+) + D(2)(ν = 0, j = 0) → HD + D(+) is theoretically analyzed for collision energies ranging from threshold up to 1.3 eV. It is assumed that the reaction takes place via formation of a collision complex. In calculations, a statistical theory is used, based on a mean isotropic potential deduced from a full potential energy surface. Calculated integral cross sections, opacity functions, and rotational distributions of the HD products are compared with recent statistical and quantum mechanical calculations performed using a full potential energy surface. Satisfactory agreement between the results obtained using the two statistical methods is found, both of which however overestimate the existing quantum mechanical predictions. The effects due to the presence of identical particles are also discussed.     

Unraveling the dynamics of the C(3P,1D) + C2H2 reactions by the crossed molecular beam scattering technique

A detailed investigation of the dynamics of the reactions of ground- and excited-state carbon atoms, C(3P) and C(1D), with acetylene is reported over a wide collision energy range (3.6-49.1 kJ mol-1) using the crossed molecular beam (CMB) scattering technique with electron ionization mass spectrometric detection and time-of-flight (TOF) analysis. We have exploited the capability of (a) generating continuous intense supersonic beams of C(3P, 1D), (b) crossing the two reactant beams at different intersection angles (45, 90, and 135 degrees ) to attain a wide range of collision energies, and (c) tuning the energy of the ionizing electrons to low values (soft ionization) to suppress interferences from dissociative ionization processes. From angular and TOF distribution measurements of products at m/z=37 and 36, the primary reaction products of the C(3P) and C(1D) reactions with C2H2 have been identified to be cyclic (c)-C3H + H, linear (l)-C3H + H, and C3 + H2. From the data analysis, product angular and translational energy distributions in the center-of-mass (CM) system for both the linear and cyclic C3H isomers as well as the C3 product from C(3P) and for l/c-C3H and C3 from C(1D) have been derived as a function of collision energy from 3.6 to 49.1 kJ mol-1. The cyclic/linear C3H ratio and the C3/(C3 + c/l-C3H) branching ratios for the C(3P) reaction have been determined as a function of collision energy. The present findings have been compared with those from previous CMB studies using pulsed beams; here, a marked contrast is noted in the CM angular distributions for both C3H- and C3-forming channels from C(3P) and their trend with collision energy. Consequently, the interpretation of the reaction dynamics derived in the present work contradicts that previously proposed from the pulsed CMB studies. The results have been discussed in the light of the available theoretical information on the relevant triplet and singlet C3H2 ab initio potential energy surfaces (PESs). In particular, the branching ratios for the C(3P) + C2H2 reaction have been compared with the available theoretical predictions (approximate quantum scattering calculations and quasiclassical trajectory calculations on ab initio triplet PESs and, very recent, statistical calculations on ab initio triplet PESs as well as on ab initio triplet/singlet PESs including nonadiabatic effects, that is, intersystem crossing). While the experimental branching ratios have been corroborated by the statistical predictions, strong disagreement has been found with the results of the dynamical calculations. The astrophysical implications of the present results have been noted.     

Translation-rotation energy levels of one H2 molecule inside the small, medium and large cages of the structure H clathrate hydrate

We report quantum dynamics calculations of the translation-rotation energy levels of one hydrogen molecule inside the small, medium and large cages of the structure H clathrate hydrate. The calculations are performed using the multiconfiguration time-dependent Hartree (MCTDH) method. Some low-lying states are computed for para-H(2), ortho-H(2), para-D(2), ortho-D(2) and HD, by block improved relaxation.     

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.     

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.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.     

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.     

Quantum dynamics of Br + HD reaction

In this article we report the results of three-dimensional time-dependent quantum wavepacket calculations carried out for the Br + HD( v = 0, j = 0) reaction in the collision energy range 0.0-1.2 eV. An accurate potential energy surface computed by Kurosaki was used for the dynamical calculations. Both reactive channels, BrH + D and BrD + H, show vibrational enhancement of the reaction cross sections. For the three initial vibrational states considered, the production of BrD channel dominates over that of BrH for the considered collision energy range. The two arrangement channels exhibit different initial rotational state dependence. The cross section for the formation of BrD is almost independent of j whereas the same for the formation of BrH increases with increase in j. A comparison with the results on an e-LEPS surface shows that the two surfaces behave very differently with respect to the cross section for the initial rotational states.     

Direct MP2 molecular dynamics studies of H atom reaction with CD4 and CH4

The mechanism and dynamics of the H + CD₄ → CD₃ + HD (I) and H + CH₄ → CH₃ + H₂ (II) reactions have been investigated by electronic structure methods. The minimum‐energy path and vibrational frequencies along the intrinsic reaction coordinate are calculated at MP2/cc‐pVDZ level. Energy distributions of the products are also obtained by the direct classical trajectory calculations at the MP2/ cc‐pVDZ level. It is found that most of the available energy appears as product translational energy, and very little of the available energy is partitioned into internal excitation of the HD (H₂) product for reaction I (II), which is in agreement with the experimental evidence. The results indicate that the experimental results could be reproduced by the direct MP2 molecular dynamics calculations. The rotational state distributions of the products show the HD (H₂) products are formed with lower rotational quantum numbers than the CD₃ (CH₃) products. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

D_2在Ni(100)表面离解吸附的量子动力学研究

用量子含时波包法研究了Ｄ２在镍表面顶－桥位上离解吸附量子动力学．计算了不同入射动能及初始振转态的离解几率．讨论了分子的同核对称性、转动取向和振动激发对离解几率的影响，并与其他理论计算结果做了比较．     

Three-state trajectory surface hopping studies of the photodissociation dynamics of formaldehyde on ab initio potential energy surfaces

Full-dimensional, three-state, surface hopping calculations of the photodissociation dynamics of formaldehyde are reported on ab initio potential energy surfaces (PESs) for electronic states S(1), T(1), and S(0). This is the first such study initiated on S(1) with ab initio-calculated spin-orbit couplings among the three states. We employ previous PESs for S(0) and T(1), and a new PES for S(1), which we describe here, as well as new spin-orbit couplings. The time-dependent electronic state populations and the branching ratio of radical products produced from S(0) and T(1) states and that of total radical products and molecular products at three total energies are calculated. Details of the surface hopping dynamics are described, and a novel pathway for isomerization on T(1) via S(0) is reported. Final translational energy distributions of H + HCO products from S(0) and T(1) are also reported as well as the translational energy distribution and final rovibrational distributions of H(2) products from the molecular channel. The present results are compared to previous trajectory calculations initiated from the global minimum of S(0). The roaming pathway leading to low rotational distribution of CO and high vibrational population of H(2) is observed in the present calculations.