Year: 2024
Type: article
Abstract: ABSTRACT The dynamic stability of various methane‐water clathrate‐like cages (CH4@(H2O)n, n = 16, 18, 20, 22), has been analyzed explicitly considering thermal effects by means of ab initio M06‐2X/6–31+G*/PCM calculations, which make use of Gaussian basis functions. Starting from the equilibrium filled cage structures, classical, dynamic reaction coordinate (DRC) on the Born–Oppenheimer surface, and semiclassical, Born–Oppenheimer plus harmonic zero‐point energy surface (BOMD), molecular dynamics have been carried out. Water molecules have a high tendency to orient covalent O–H bonds tangentially to the hydrophobic surface, thus clathrate‐like arrangements are an acceptable model to fully hydrate methane. If the cage size is such as to minimize core repulsion, due to electron cloud overlap, and to maximize host–guest van der Waals attractions, the clathrate‐like structures have a life‐time of two picoseconds in classical DRC simulations. The inclusion of quantum kinetic energy in BOMD simulations results in less structured cages with a reduced amount of hydrogen bond network. The preferential tangential orientation of the O‐H bonds is largely maintained, although few of them point toward the methane for a very short time in BOMD simulations. The reduced configurational space of water molecules hydrating hydrophobic moiety is highlighted, thus any satisfactory molecular modeling has to account for it.
DOI: 10.1002/qua.27487
Author: Giuseppe Lanza
Institution: University of Catania
Language: English
Subfield: Atomic and Molecular Physics, and Optics
Field: Physics and Astronomy
Domain: Physical Sciences
Open Access status: hybrid
Estimated APC paid: $4,330
Funder: Università di Catania
OpenAlex ID: W4402732793