Clathrate hydrates hold considerable promise as safe and economical materials for hydrogen storage. Here we present a quantum mechanical study of H 2 and D 2 diffusion through a hexagonal face shared by two large cages of clathrate hydrates over a wide range of temperatures. Path integral molecular dynamics simulations are used to compute the free-energy profiles for the diffusion of H 2 and D 2 as a function of temperature. Ring polymer molecular dynamics rate theory, incorporating both exact quantum statistics and approximate quantum dynamical effects, is utilized in the calculations of the H 2 and D 2 diffusion rates in a broad temperature interval. We find that the shape of the quantum free-energy profiles and their height relative to the classical free energy barriers at a given temperature, as well as the rate of diffusion, are profoundly affected by competing quantum effects: above 25 K, zero-point energy (ZPE) perpendicular to the reaction path for diffusion between cavities decreases the quantum rate compared to the classical rate, whereas at lower temperatures tunneling outcompetes the ZPE and as result the quantum rate is greater than the classical rate.

We report a theoretical study of the frequency shift (redshift) of the stretching fundamental transition of an H molecule confined inside the small dodecahedral cage of the structure II clathrate hydrate and its dependence on the condensed-phase environment. In order to determine how much the hydrate water molecules beyond the confining small cage contribute to the vibrational frequency shift, quantum five-dimensional (5D) calculations of the coupled translation-rotation eigenstates are performed for H in the v=0 and v=1 vibrational states inside spherical clathrate hydrate domains of increasing radius and a growing number of water molecules, ranging from 20 for the isolated small cage to over 1900. In these calculations, both H and the water domains are treated as rigid. The 5D intermolecular potential energy surface (PES) of H inside a hydrate domain is assumed to be pairwise additive. The H-HO pair interaction, represented by the 5D (rigid monomer) PES that depends on the vibrational state of H, v=0 or v=1, is derived from the high-quality ab initio full-dimensional (9D) PES of the H-HO complex [P. Valiron et al., J. Chem. Phys. 129, 134306 (2008)]. The H vibrational frequency shift calculated for the largest clathrate domain considered, which mimics the condensed-phase environment, is about 10% larger in magnitude than that obtained by taking into account only the small cage. The calculated splittings of the translational fundamental of H change very little with the domain size, unlike the H j = 1 rotational splittings that decrease significantly as the domain size increases. The changes in both the vibrational frequency shift and the j = 1 rotational splitting due to the condensed-phase effects arise predominantly from the HO molecules in the first three complete hydration shells around H.

In the supramolecular complex H2@C60, the lightest of molecules, H2, is encapsulated inside the most highly symmetric molecule C60. The elegance and apparent simplicity of H2@C60 conceal highly intricate quantum dynamics of the coupled translational and rotational motions of the guest molecule in a nearly spherical nanoscale cavity, which embodies some of the most fundamental concepts of quantum mechanics. Here we present the first rigorous and highly accurate quantum calculations of the inelastic neutron scattering (INS) spectra of this prototypical endohedral fullerene complex and their temperature dependence. The calculations enable complete assignment of the recently reported experimental INS spectra of H2@C60 measured at several temperatures. We also derive a new and unexpected selection rule for the INS spectroscopy of H2 in a near-spherical confinement, which explains why the INS transitions between certain translation-rotation eigenstates of H2 in C60 have zero intensity and do not appear in the spectra.

We have performed high-resolution inelastic neutron scattering (INS) measurements on binary hydrogen clathrate hydrates exhibiting the hexagonal structure (sH). Two samples, differing only in the ortho/para fraction of hydrogen, were prepared using heavy water and methyl tert-butyl ether as the promoter in its perdeuterated form. The INS spectrum of the translation-rotation (TR) excitations of the guest H2 molecule was obtained by subtracting the very weak signal due to the D2O lattice modes. By means of a subtraction procedure, it has been possible to obtain separately the spectra of caged p-H2 and o-H2. sH clathrates are comprised of three distinct types of cages, two of which, differing in shape and size, are each occupied by one H2 molecule only. Both contribute to the measured INS spectrum which is, therefore, rather complex and challenging to assign unambiguously. To assist with the interpretation, the INS spectra are calculated accurately utilizing the quantum methodology which incorporates the coupled five-dimensional TR energy levels and wave functions of the H2 molecule confined in each type of nanocage. The computed INS spectra are highly realistic and reflect the complexity of the coupled TR dynamics of the guest H2 in the anisotropic confining environment. The simulated INS spectra of p-H2 and o-H2 in the small and medium cages are compared with the experimental data, and are indispensable for their interpretation.