摘要

Endohedral CH4@(H2O)(n) (n = 16, 18, 20, 22, 24) clusters with standard and nonstandard cage configurations containing four-, five-, six-, seven-membered rings were generated by spiral algorithm and were systematically explored using DFT-D methods. The geometries of all isomers were optimized in vacuum and aqueous solution. In vacuum, encapsulation of methane molecules can stabilize the hollow (H2O)(n) cage by 2.31 similar to 5.44 kcal/mol; but the endohedral CH4@(H2O)(n) cages are still less stable than the pure (H2O)(n) clusters. Aqueous environment could promote the stabilities of the hollow (H2O) n cages as well as the CH4@(H2O)(n) clusters, and the CH4@(H2O)(n) clusters possess larger stabilization energies with regard to the pure (H2O) (n) clusters except for n = 24. The lowest energy structures of the CH4@(H2O)(20) and CH4@(H2O)(24) cages are identical to the building units in the crystalline sI clathrate hydrate. All of the low-energy cages (including both regular and irregular ones) have large structural similarity and can be connected by "dimer-insertion" operation and Stone-Wales transformation. Our calculation also showed that in the range of cluster size n = 16-24, the relative energies of cage isomers tend to decrease with increasing number of the adjacent pentagons in the oxygen skeleton structures. In addition to the regular endohedral CH4@(H2O)(20) and CH4@(H2O)(24) cage structures, some nonstandard CH4@(H2O)(n) (n = 18, 20, 22, 24) cages have lower energies and might appear during nucleation process of methane hydrate. For the methane molecules in these low-energy cage isomers, we found that the C-H symmetric stretching frequencies show a red-shift trend and the C-13 NMR chemical shifts generally move toward negative values as the cavity size increases. These theoretical results are comparable to the available experimental data and might help experimental identification of the endohedral water cages during nucleation.