Molecular dynamics of zigzag single walled carbon nanotube immersion in water

The results of enthalpy of immersion in water for finite single-walled carbon nanotubes are reported. Using molecular dynamics simulation, we discuss the relation between the value of this enthalpy and tube diameters showing that the obtained plot can be divided into three regions. The structure of water inside tubes in all three regions is discussed and it is shown that the existence of the strong maximum of enthalpy observed for tube diameter ca. 1.17 nm is due to freezing of water under confinement. The calculations of hydrogen bond statistics and water density profiles inside tubes are additionally reported to confirm the obtained results. Next, we show the results of calculation for the same tubes but containing surface carbonyl oxygen groups at pore entrances. A remarkable rise in the value of enthalpy of immersion in comparison to the initial tubes is observed. We also discuss the influence of charge distribution between oxygen and carbon atom forming surface carbonyls on the structure of confined water. It is concluded for the first time that the presence of surface oxygen atoms at the pore entrances remarkably influences the structure and stability of ice created inside nanotubes, and surface carbonyls appear to be chaotropic (i.e. structure breaking) for confined water. This effect is explained by the pore blocking leading to a decrease (compared to initial structure) in the number of confined water molecules after introduction of surface oxygen groups at pore entrances.