CO2 and CH4 clathrate hydrates of type I were studied by means of DFT and QTAIM, in order to better understand their properties at the molecular level. Sub-cells of type I hydrates were modeled as independent rigid cages, both empty and containing guest molecules. Interaction potentials of guest molecules inside each cage, and moving from a cell to the adjacent one, were calculated using the DFT approximation B3LYP/6-311+g(d,p), considering the cases with and without long range Coulombic corrections. The selected theory level was validated by comparison of the simulated Raman spectra with the experimental ones, for the case of type I lattice at full occupation of CO2 and CH4, respectively. For this comparison, Fermi resonances of CO2 were taken into account by transforming experimental bands to the corresponding theoretical non-mixed states. On the one hand, our results confirm the validity of the theory level selected for the model. We have shown the high anisotropy of the guest-cell interaction potential of the molecules analyzed, which has implications in the formulation and use of equations of state, and in the study of transport properties as well. On the other hand, our results suggest that the concentration of guest species inside type I hydrates could be computed from the comparison of experimental and predicted Raman spectra, although there are non-trivial experimental limitations to get over for that purpose.

In this work, we use dual cage explicit atomic systems to demonstrate theoretically that direct transitions are feasible through hexagonal and pentagonal faces in type I hydrate without compromising the overall structure integrity.

Carbon dioxide (CO) molecules show a rich orientation landscape when they are enclathrated in type I hydrates. Previous studies have described experimentally their preferential orientations, and some theoretical works have explained, but only partially, these experimental results. In the present paper, we use classical molecular dynamics and electronic density functional theory to advance in the theoretical description of CO orientations within type I hydrates. Our results are fully compatible with those previously reported, both theoretical and experimental, the geometric shape of the cavities in hydrate being, and therefore, the steric constraints, responsible for some (but not all) preferential angles. In addition, our calculations also show that guest-guest interactions in neighbouring cages are a key factor to explain the remaining experimental angles. Besides the implication concerning equation of state hydrate modeling approximations, the conclusion is that these guest-guest interactions should not be neglected, contrary to the usual practice.