Molecular dynamic simulations were performed for ice I h with a free surface by using four water models, SPC/E, TIP4P, TIP4P/Ice and TIP4P/2005. The behavior of the basal plane, the primary prismatic plane and of the secondary prismatic plane when exposed to vacuum was analyzed. We observe the formation of a thin liquid layer at the ice surface at temperatures below the melting point for all models and the three planes considered. For a given plane it was found that the thickness of a liquid layer was similar for different water models, when the comparison is made at the same undercooling with respect to the melting point of the model. The liquid layer thickness is found to increase with temperature. For a fixed temperature it was found that the thickness of the liquid layer decreases in the following order: the basal plane, the primary prismatic plane, and the secondary prismatic plane. For the TIP4P/Ice model, a model reproducing the experimental value of the melting temperature of ice, the first clear indication of the formation of a liquid layer appears at about -100 Celsius for the basal plane, at about -80 Celsius for the primary prismatic plane and at about -70 Celsius for the secondary prismatic plane.

No abstract

Cite as : C. Vega, M. Martin-Conde and A.Patrykiejew, Mol.Phys., 104, 3583, (2006). Molecular dynamic simulations were performed for ice I h with a free surface. The simulations were carried out at several temperatures and each run lasted more than 7ns. At high temperatures the ice melts. It is demonstrated that the melting process starts at the surface and propagates to the bulk of the ice block. Already at the temperatures below the melting point, we observe a thin liquid layer at the ice surface, but the block of ice remains stable along the run. As soon as the temperature reaches the melting point the entire ice block melts. Our results demonstrate that, unlike in the case of conventional simulations in the NpT ensemble, overheating of the ice I h with a free surface does not occur. That allows to estimate the melting point of ice at zero pressure.We applied the method to the following models of water: SPC/E, TIP4P, TIP4P/Ew, TIP4P/Ice and TIP4P/2005, and found good agreement between the melting temperatures obtained by this procedure and the values obtained either from free energy calculations or from direct simulations of the ice/water interface.

No abstract

We present a density functional theory of nonuniform ionic fluids. This theory is based on the application of the electrostatic contribution to the free energy functional arising from mean spherical approximation for a bulk restricted primitive model and from the energy route bulk equation of state. In order to employ this functional we define a reference fluid and additional averaged densities, according to the approach introduced by Gillespie, Nonner and Eisenberg [J. Phys.: Condens. Matter 14, 12129 (2002)]. In the case of bulk systems the proposed theory reduces to the mean spherical approximation equation of state, arising from the energy route and thus it predicts the first-order phase transition. We use this theory to investigate the effects of confinement on the liquid-vapor equilibria. Two cases are considered, namely an electrolyte confined to the pore with uncharged walls and with charged walls. The dependence of the capillary evaporation diagrams on the pore width and on the electrostatic potential is determined.