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.     

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.     

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.     

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.     

Considering multiple occupancy of cavities in clathrate hydrate phase equilibrium calculations

One of the major assumptions of the original van der Waals–Platteeuw (vdWP) model is the single occupancy of hydrate cavities. In this work, the vdWP model is modified to also account for multiple occupancies of hydrate cavities by small molecules. The developed model is evaluated by calculating the hydrate equilibrium conditions with either oxygen or nitrogen as guest molecules in pure form, as well as mixtures of nitrogen and propane (molecules of these pure gases and those in (nitrogen + propane) have double occupancy in large cavities of structure II up to a certain concentration of propane). The results of this modified model show good agreement with the experimental data reported in the literature.     

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.     

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.     

Effect of small cage guests on hydrogen bonding of tetrahydrofuran in binary structure II clathrate hydrates

Molecular dynamics simulations of the pure structure II tetrahydrofuran clathrate hydrate and binary structure II tetrahydrofuran clathrate hydrate with CO(2), CH(4), H(2)S, and Xe small cage guests are performed to study the effect of the shape, size, and intermolecular forces of the small cages guests on the structure and dynamics of the hydrate. The simulations show that the number and nature of the guest in the small cage affects the probability of hydrogen bonding of the tetrahydrofuran guest with the large cage water molecules. The effect on hydrogen bonding of tetrahydrofuran occurs despite the fact that the guests in the small cage do not themselves form hydrogen bonds with water. These results indicate that nearest neighbour guest-guest interactions (mediated through the water lattice framework) can affect the clathrate structure and stability. The implications of these subtle small guest effects on clathrate hydrate stability are discussed.     

Regularity of low-lying quantum eigenstates in a classically mixing system

Low-lying quantum energy eigenstates of the classically mixing stadium system are shown to be highly regular and to be well described by an adiabatic separable Hamiltonian. The results clearly demonstrate quantum regularity in an energy regime of extreme classical chaos.     

Accuracy of gates in a quantum computer based on vibrational eigenstates

A model is developed to study the properties of a quantum computer that uses vibrational eigenstates of molecules to implement the quantum information bits and shaped laser pulses to apply the quantum logic gates. Particular emphasis of this study is on understanding how the different factors, such as properties of the molecule and of the pulse, can be used to affect the accuracy of quantum gates in such a system. Optimal control theory and numerical time-propagation of vibrational wave packets are employed to obtain the shaped pulses for the gates NOT and Hadamard transform. The effects of the anharmonicity parameter of the molecule, the target time of the pulse and of the penalty function are investigated. Influence of all these parameters on the accuracy of qubit transformations is observed and explained. It is shown that when all these parameters are carefully chosen the accuracy of quantum gates reaches 99.9%.     

On the pressure-induced phase transformation in the structure II clathrate hydrate of tetrahydrofuran

A high-pressure phase of the clathrate hydrate of tetrahydrofuran was prepared by freezing a liquid phase of overall composition THF · 7 H₂O under a pressure of 3.0 kbar, or by pressurizing the solid structure II THF hydrate of 255K to 3.4 kbar. Unfortunately, the products recovered at 77K were always mixed phase materials as shown by X-ray powder diffraction. A number of diffraction lines could be indexed in terms of the cubic structure I hydrate with a slightly expanded lattice parameter, 12.08 Å, giving some support to Dyadin's idea that the high pressure phase transition involves a conversion of Structure II to Structure I. Other phases observed in the recovered product include Ice IX and amorphous materials. The reversion of the high pressure sample to the structure II hydrate was followed by differential scanning calorimetry. At ambient pressure, the high pressure sample converts slowly back to Structure II hydrate event at 77K.NRCC No. 35786.     

Ordering of guest-molecule dipoles in the structure I clathrate hydrate of trimethylene oxides

The rotational mobility of encaged trimethylene oxide (TMO) molecules was studied down to 1.8°K by sub-MHz dielectric measurements of the structure I H₂O clathrate and by proton magnetic resonance measurements of the corresponding D₂O clathrate. The results indicate that below a transitional temperature range about 105°K most TMO dipoles assume parallel alignment along the $\text{4}$ axes of the cages. Below the transition the proton second moment suggests the presence of hindered rotation of TMO about its polar axis until the rigid-lattice condition is reached below 5°K. Some residual very broad dielectric absorption (activation energy 2.1 kcal/mole) persists to very low temperatures. Guest-guest and guest-host interaction energies are calculated for simple models     

Dynamic (hyper)polarizabilities of the sulphur dioxide molecule: coupled cluster calculations including vibrational corrections

In this work we report results for dynamical (hyper)polarizabilities of the sulphur dioxide molecule with inclusion of vibrational corrections. The electronic contributions were computed analytically at the single and double coupled cluster level through response theories for the frequencies 0, 0.0239, 0.0428, 0.0656, 0.0720, and 0.0886 hartree. Contributions of the connected triple excitations to the dynamic electronic properties were also estimated through the multiplicative correction scheme. Vibrational corrections were calculated by means of the perturbation theoretical method. The results obtained show that the zero point vibrational correction is very small for all properties studied while the pure vibrational correction is relevant for the dc-Pockels effect, intensity dependent refractive index, and dc-Kerr effect. For these nonlinear optical processes, the pure vibrational corrections represent approximately 75%, 13%, and 6% of the corresponding electronic contributions for the higher frequencies quoted. The results presented for the polarizability are in good agreement with experimental values available in the literature. For the hyperpolarizabilities we have not obtained experimental results with precision sufficient for comparison.     

129Xe nuclear magnetic resonance in the clathrate hydrate of xenon

The sensitivity of the nuclear magnetic shielding of the Xe atom to the physical environment makes it possible to “see” well-separated. ¹²⁹Xe resonances from Xe atoms in the small and large clathrate deuteriohydrate cages as well as in the gas with which the hydrate is in equilibrium. The relative occupancy of the small cages is suostantially less than predicted oy existing models of guest-host interactions in Xe gas hydrate.     

The chemical shifts of Xe in the cages of clathrate hydrate Structures I and II

We report, for the first time, a calculation of the isotropic NMR chemical shift of 129Xe in the cages of clathrate hydrates Structures I and II. We generate a shielding surface for Xe in the clathrate cages by quantum mechanical calculations. Subsequently this shielding surface is employed in canonical Monte Carlo simulations to find the average isotropic Xe shielding values in the various cages. For the two types of cages in clathrate hydrate Structure I, we find the intermolecular shielding values [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-214.0 ppm, and [sigma(Xe@5(12)6(2) cage)-sigma(Xe atom)]=-146.9 ppm, in reasonable agreement with the values -242 and -152 ppm, respectively, observed experimentally by Ripmeester and co-workers between 263 and 293 K. For the 5(12) and 5(12)6(4) cages of Structure II we find [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-206.7 ppm, and [sigma(Xe@5(12)6(4) cage)-sigma(Xe atom)]=-104.7 ppm, also in reasonable agreement with the values -225 and -80 ppm, respectively, measured in a Xe-propane type II mixed clathrate hydrate at 77 and 220-240 K by Ripmeester et al.     

Laser Raman spectroscopy of THF clathrate hydrate in the temperature range 90-300 K

In situ Raman and Fourier transform infrared (FTIR-NIR) spectroscopic studies on tetrahydrofuran (THF-C(4)H(8)O) clathrate hydrate (CH) were reported. The Raman results in lattice (64 cm(-1)), ring breathing and C-H stretching mode regions are in conformity with earlier reports, while the FTIR (NIR) studies in second order mode region were reported for the first time. Comparison of the results indicate that the band assigned to ring breathing mode around 922 cm(-1) (in Raman) and corresponding second order mode in NIR around 4295 cm(-1) broadens and shifts in enclathrated THF. The ring breathing mode at lower temperatures (T<120 K) is highly asymmetric and splits into two and are due to different host-guest interactions at lower temperatures.     

Growth mechanism of a gas clathrate hydrate from a dilute aqueous gas solution: a molecular dynamics simulation of a three-phase system

A molecular dynamics simulation of a three-phase system including a gas clathrate, liquid water, and a gas was carried out at 298 K and high pressure in order to investigate the growth mechanism of the clathrate from a dilute aqueous gas solution. The simulation indicated that the clathrate grew on interfaces between the clathrate and the liquid water, after transfer of the gas molecules from the gas phase to the interfaces. The results suggest a two-step process for growth: first, gas molecules are arranged at cage sites, and second, H(2)O molecules are ordered near the gas molecules. The results also suggest that only the H(2)O molecules, which are surrounded or sandwiched by the gas molecules, form the stable polygons that constitute the cages of the clathrate. In addition, the growth of the clathrate from a concentrated aqueous gas solution was also simulated, and the results suggested a growth mechanism in which many H(2)O and gas molecules correctively form the structure of the clathrate. The clathrate grown from the concentrated solution contained some empty cages, whereas the formation of empty cages was not observed during the growth from the dilute solution. The results obtained by both simulations are compared with the results of an experimental study, and the growth mechanism of the clathrate in a real system is discussed.