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 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.     

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.     

H2, HD, and D2 inside C60: coupled translation-rotation eigenstates of the endohedral molecules from quantum five-dimensional calculations

We have performed rigorous quantum five-dimensional (5D) calculations of the translation-rotation (T-R) energy levels and wave functions of H(2), HD, and D(2) inside C(60). This work is an extension of our earlier investigation of the quantum T-R dynamics of H(2)@C(60) [M. Xu et al., J. Chem. Phys. 128, 011101 (2008)] and uses the same computational methodology. Two 5D intermolecular potential energy surfaces (PESs) were employed, differing considerably in their well depths and the degree of confinement of the hydrogen molecule. Our calculations revealed pronounced sensitivity of the endohedral T-R dynamics to the differences in the interaction potentials, and to the large variations in the masses and the rotational constants of H(2), HD, and D(2). The T-R levels vary significantly in their energies and ordering on the two PESs, as well as from one isotopomer to another. Nevertheless, they all display the same distinctive patterns of degeneracies, which can be qualitatively understood and assigned in terms the model which combines the isotropic three-dimensional harmonic oscillator, the rigid rotor, and the coupling between the orbital and the rotational angular momenta of H(2)/HD/D(2). The quantum number j associated with the rotation of H(2), HD, and D(2) was found to be a good quantum number for H(2) and D(2) on both PESs, while most of the T-R levels of HD exhibit strong mixing of two or more rotational basis functions with different j values.     

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.     

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.     

Interaction potential and infrared absorption of endohedral H2 in C60

We have measured the temperature dependence of the infrared spectra of a hydrogen molecule trapped inside a C(60) cage, H(2)@C(60), in the temperature range from 6 to 300 K and analyzed the excitation spectrum by using a five-dimensional model of a vibrating rotor in a spherical potential. The electric dipole moment is induced by the translational motion of endohedral H(2) and gives rise to an infrared absorption process where one translational quantum is created or annihilated, ΔN = ±1. Some fundamental transitions, ΔN = 0, are observed as well. The rotation of endohedral H(2) is unhindered but coupled to the translational motion. The isotropic and translation-rotation coupling part of the potential are anharmonic and different in the ground and excited vibrational states of H(2). The vibrational frequency and the rotational constant of endohedral H(2) are smaller than those of H(2) in the gas phase. The assignment of lines to ortho- and para-H(2) is confirmed by measuring spectra of a para enriched sample of H(2)@C(60) and is consistent with the earlier interpretation of the low temperature infrared spectra [Mamone et al., J. Chem. Phys. 130, 081103 (2009)].     

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.     

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.     

New ab initio potential energy surface and the vibration-rotation-tunneling levels of (H2O)2 and (D2O)2

We report a new full-dimensional potential energy surface (PES) for the water dimer, based on fitting energies at roughly 30,000 configurations obtained with the coupled-cluster single and double, and perturbative treatment of triple excitations method using an augmented, correlation consistent, polarized triple zeta basis set. A global dipole moment surface based on Moller-Plesset perturbation theory results at these configurations is also reported. The PES is used in rigorous quantum calculations of intermolecular vibrational frequencies, tunneling splittings, and rotational constants for (H2O)2 and (D2O)2, using the rigid monomer approximation. Agreement with experiment is excellent and is at the highest level reported to date. The validity of this approximation is examined by comparing tunneling barriers within that model with those from fully relaxed calculations.     

Rotational state specific dissociation dynamics of HOD → H + OD via two-photon excitation to the C̃ electronic state

The dissociation dynamics of HOD via two-photon excitation to the C? state have been investigated using the H-atom Rydberg tagging time-of-flight (TOF) technique. The H-atom action spectrum for the C? ← X? transition shows resolved rotational structure. Product translational energy distributions and angular distributions have also been recorded for the H + OD channel for three excited levels each with k(a)′ = 2. From these distributions, quantum state distributions and angular anisotropy parameters (β2 and β4) for the OD product were determined. These results are consistent with the nonadiabatic predissociation picture illustrated in the one-photon dissociation process for H2O. The heterogeneous dissociation pathway via Coriolis coupling is the dominant dissociation process in the present study. A high proportion of the total available energy is deposited into the rotational energy of the OD product. The anisotropic recoil distributions reveal the distinctive contributions from the alignment of the excited states and the dissociation process. Comparisons are also made between the results for HOD and H2O via the equivalent rotational transitions. The OH bond energy, D(o)(H?OD), of the HOD molecule is also determined to be 41283.0 ± 5 cm(-1).     

Imaging H2O photofragments in the predissociation of the HCl-H2O hydrogen-bonded dimer

The state-to-state vibrational predissociation (VP) dynamics of the hydrogen-bonded HCl-H(2)O dimer was studied following excitation of the dimer's HCl stretch by detecting the H(2)O fragment. Velocity map imaging (VMI) and resonance-enhanced multiphoton ionization (REMPI) were used to determine pair-correlated product energy distributions. Following vibrational excitation of the HCl stretch of the dimer, H(2)O fragments were detected by 2 + 1 REMPI via the C (1)B(1) (000) ← X (1)A(1) (000) transition. REMPI spectra clearly show H(2)O from dissociation produced in the ground vibrational state. The fragments' center-of-mass (c.m.) translational energy distributions were determined from images of selected rotational states of H(2)O and were converted to rotational state distributions of the HCl cofragment. The distributions were consistent with the previously measured dissociation energy of D(0) = 1334 ± 10 cm(-1) and show a clear preference for rotational levels in the HCl fragment that minimize translational energy release. The usefulness of 2 + 1 REMPI detection of water fragments is discussed.     

Discrete variable approaches to tetratomic molecules: part II: application to H2O2 and H2CO

We have carried out large-scale calculations for accurate vibrational energy levels of formaldehyde and hydrogen peroxide. The discrete variable representations of the radial and angular coordinates are employed together with the contraction scheme resulting from several diagonalization/truncation steps. The global potential energy surface due to Carter et al. [J. Mol. Spectrosc. 90 (1997) 729] is used for H2CO and due to Koput et al. [J. Phys. Chem. A 102 (1998) 6325] for H2O2. For both molecules, the calculated vibrational energy levels are characterized by combining vibrationally averaged geometries and expectation values of rotational constants with several adiabatic projection schemes for automatic quantum number assignments. The energy levels of H2CO involving the excited v2 and v3 vibrations appear as resonances beyond the zero-order picture consisting of uncoupled 3D stretching and 2D bending modes. The torsional energy levels of H2O2 are studied in great detail and different energy patterns occurring below and above the cis barrier are discussed. Our full dimensional calculations for H2O2 have shown that the OH triad levels, 2vOH, are symmetry adapted local mode states.     

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相似文献   

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.     

Molecular dynamics with quantum transitions study of the vibrational relaxation of the HOD bend fundamental in liquid D2O

The molecular dynamics with quantum transitions method is used to study the vibrational relaxation of the HOD bend fundamental in liquid D(2)O. All of the vibrational bending degrees of freedom of the HOD and D(2)O molecules are described by quantum mechanics, while the remaining translational and rotational degrees of freedom are described classically. The effect of the coupling between the rotational and vibrational degrees of freedom of the deuterated water molecules is analyzed. A kinetic mechanism based on three steps is proposed in order to interpret the dynamics of the system. It is shown that intermolecular vibrational energy transfer plays an important role in the relaxation process and also that the transfer of energy into the rotational degrees of freedom is favored over the transfer of energy into the translational motions. The thermalization of the system after the relaxation is reached in a shorter time scale than that of the recovery of the hydrogen bond network. The relaxation and equilibration times obtained compare well with experimental and previous theoretical results.     

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

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

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.     

Product angular distributions in the ultraviolet photodissociation of N2O

The angular distribution of products from the ultraviolet photodissociation of nitrous oxide yielding O((1)D) and N(2)(X Σ(g)(+)(1)) was investigated using classical trajectory calculations. The calculations modeled absorption only to the 2(1)A(') electronic state but used surface-hopping techniques to model nonadiabatic transitions to the ground electronic state late in the dissociation. Observed values of the anisotropy parameter β, which decrease as the product N(2) rotational quantum number j increases, could be well reproduced. The relatively low observed β values arise principally from nonaxial recoil due to the very strong bending forces present in the excited state. In the main part of the product rotational distribution near 203 nm, an unusual dynamical effect produces the decrease in β with increasing j; nonaxial recoil effects remain approximately constant while higher j product molecules arise from parent molecules that had their transition dipole moments aligned more closely along the molecular axis. In both low and high j tails of the rotational distribution, the variations in β with j are caused by changes in the extent of nonaxial recoil. In the high-j tail, additional torque present on the ground state potential energy surface following nonadiabatic transitions causes both the additional rotational excitation and the lower β values.     

Rovibrational product state distribution for inelastic H+D2 collisions

Experimental measurements of rovibrational product state distributions for the inelastic scattering process H + D2(nu=0,j)-->H + D2(nu' = 1,2,j') are presented and compared with the results of quasiclassical and quantum mechanical calculations. Agreement between theory and experiment is almost quantitative. Two subtle trends are found: the relative amount of energy in product rotational excitation decreases slightly with increasing collision energy and increases slightly with increasing product vibrational excitation. These trends are the reverse of what has been found for reactive scattering in which the opposite trends are much more pronounced.