Nonstandard cages in the formation process of methane clathrate: Stability, structure, and spectroscopic implications from first-principles

• Published in: The Journal of Chemical Physics
• Main Authors: Tang, Lingli, Su, Yan, Liu, Yuan, Zhao, Jijun, Qiu, Ruifeng
• Other Authors: National Natural Science Foundation of China
• Format: Article in Journal/Newspaper
• Language: English
• Published: AIP Publishing 2012
• Subjects: Physical and Theoretical Chemistry; General Physics and Astronomy; Methane hydrate
• Online Access: http://dx.doi.org/10.1063/1.4728157; https://pubs.aip.org/aip/jcp/article-pdf/doi/10.1063/1.4728157/13968358/224508_1_online.pdf

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∼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 13C 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.