Molecular Dynamics Simulation of the Influence of Nanoscale Structure on Water Wetting and Condensation.

Recent advances in the microfabrication technology have made it possible to control surface properties at micro- and nanoscale levels. Functional surfaces drastically change wettability and condensation processes that are essential for controlling of heat transfer. However, the direct observation of condensation on micro- and nanostructure surfaces is difficult, and further understanding of the effects of the microstructure on the phase change is required. In this research, the contact angle of droplets with a wall surface and the initial condensation process were analyzed using a molecular dynamics simulation to investigate the impact of nanoscale structures and their adhesion force on condensation. The results demonstrated the dependence of the contact angle of the droplets and condensation dynamics on the wall structure and attractive force of the wall surface. Condensed water droplets were adsorbed into the nanostructures and formed a water film in case of a hydrophilic surface.

A molecular dynamics study of guest-host hydrogen bonding in alcohol clathrate hydrates.

Clathrate hydrates are typically stabilized by suitably sized hydrophobic guest molecules. However, it has been experimentally reported that isomers of amyl-alcohol C5H11OH can be enclosed into the 5(12)6(4) cages in structure II (sII) clathrate hydrates, even though the effective radii of the molecules are larger than the van der Waals radii of the cages. To reveal the mechanism of the anomalous enclathration of hydrophilic molecules, we performed ab initio and classical molecular dynamics simulations (MD) and analyzed the structure and dynamics of a guest-host hydrogen bond for sII 3-methyl-1-butanol and structure H (sH) 2-methyl-2-butanol clathrate hydrates. The simulations clearly showed the formation of guest-host hydrogen bonds and the incorporation of the O-H group of 3-methyl-1-butanol guest molecules into the framework of the sII 5(12)6(4) cages, with the remaining hydrophobic part of the amyl-alcohol molecule well accommodated into the cages. The calculated vibrational spectra of alcohol O-H bonds showed large frequency shifts due to the strong guest-host hydrogen bonding. The 2-methyl-2-butanol guests form strong hydrogen bonds with the cage water molecules in the sH clathrate, but are not incorporated into the water framework. By comparing the structures of the alcohols in the hydrate phases, the effect of the location of O-H groups in the butyl chain of the guest molecules on the crystalline structure of the clathrate hydrates is indicated.

Observation of interstitial molecular hydrogen in clathrate hydrates.

The current knowledge and description of guest molecules within clathrate hydrates only accounts for occupancy within regular polyhedral water cages. Experimental measurements and simulations, examining the tert-butylamine + H2 + H2O hydrate system, now suggest that H2 can also be incorporated within hydrate crystal structures by occupying interstitial sites, that is, locations other than the interior of regular polyhedral water cages. Specifically, H2 is found within the shared heptagonal faces of the large (4(3)5(9)6(2)7(3)) cage and in cavities formed from the disruption of smaller (4(4)5(4)) water cages. The ability of H2 to occupy these interstitial sites and fluctuate position in the crystal lattice demonstrates the dynamic behavior of H2 in solids and reveals new insight into guest-guest and guest-host interactions in clathrate hydrates, with potential implications in increasing overall energy storage properties.

Water proton configurations in structures I, II, and H clathrate hydrate unit cells.

Position and orientation of water protons need to be specified when the molecular simulation studies are performed for clathrate hydrates. Positions of oxygen atoms in water are experimentally determined by X-ray diffraction analysis of clathrate hydrate structures, but positions of water hydrogen atoms in the lattice are disordered. This study reports a determination of the water proton coordinates in unit cell of structure I (sI), II (sII), and H (sH) clathrate hydrates that satisfy the ice rules, have the lowest potential energy configuration for the protons, and give a net zero dipole moment. Possible proton coordinates in the unit cell were chosen by analyzing the symmetry of protons on the hexagonal or pentagonal faces in the hydrate cages and generating all possible proton distributions which satisfy the ice rules. We found that in the sI and sII unit cells, proton distributions with small net dipole moments have fairly narrow potential energy spreads of about 1 kJ∕mol. The total Coulomb potential on a test unit charge placed in the cage center for the minimum energy∕minimum dipole unit cell configurations was calculated. In the sI small cages, the Coulomb potential energy spread in each class of cage is less than 0.1 kJ∕mol, while the potential energy spread increases to values up to 6 kJ∕mol in sH and 15 kJ∕mol in the sII cages. The guest environments inside the cages can therefore be substantially different in the sII case. Cartesian coordinates for oxygen and hydrogen atoms in the sI, sII, and sH unit cells are reported for reference.

Vibrational modes of methane in the structure H clathrate hydrate from ab initio molecular dynamics simulation.

Vibrational spectra of guest molecules in clathrate hydrates are frequently measured to determine the characteristic signatures of the molecular environment and dynamical properties of guest-host interactions. Here, we present results of our study on the vibrational frequencies of methane molecules in structure H clathrate hydrates, namely, in the 5(12) and 4(3)5(6)6(3) cages, as the frequencies of stretching vibrational modes in these environments are still unclear. The vibrational spectra of methane molecules in structure H clathrate hydrate were obtained from ab initio molecular dynamics simulation and computed from Fourier transform of autocorrelation functions for each distinct vibrational mode. The calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules were found to be lower in the 4(3)5(6)6(3) cages than in the 5(12) cages (3.8 cm(-1) for symmetric stretching and 6.0 cm(-1) for asymmetric stretching). The C-H bond length and average distance between methane molecules and host-water molecules in 4(3)5(6)6(3) cages were slightly longer than those in the 5(12) cages.

Molecular vibrations of methane molecules in the structure I clathrate hydrate from ab initio molecular dynamics simulation.

Vibrational frequencies of guest molecules in clathrate hydrates reflect the molecular environment and dynamical behavior of molecules. A detailed understanding of the mechanism for the vibrational frequency changes of the guest molecules in the clathrate hydrate cages is still incomplete. In this study, molecular vibrations of methane molecules in a structure I clathrate hydrate are calculated from ab initio molecular dynamics simulation. The vibrational spectra of methane are computed by Fourier transform of autocorrelation functions, which reveal distinct separation of each vibrational mode. Calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules are lower in the large cages than in the small cages (8 and 16 cm(-1) for symmetric and asymmetric stretching, respectively). These changes are closely linked with the C-H bond length. The vibrational frequencies for the bending and rocking vibrational modes nearly overlap in each of the cages.