No Arabic abstract
Ice-water, water-vapor interfaces and ice surface are studied by molecular dynamics simulations with the SPC/E model of water molecules having the purpose to estimate the profiles of electrostatic potential across the interfaces. We have proposed a methodology for calculating the profiles of electrostatic potential based on a trial particle, which showed good agreement for the case of electrostatic potential profile of the water-vapor interface of TIP4P model calculated in another way. The measured profile of electrostatic potential for the pure ice-water interface decreases towards the liquid bulk region, which is in agreement with simulations of preferential direction of motion of Li$^{+}$ and F$^{-}$ solute ions at the liquid side of the ice-water interface. These results are discussed in connection with the Workman-Reynolds effect.
Crystallization from a supercooled liquid initially proceeds via the formation of a small solid embryo (nucleus), which requires surmounting an activation barrier. This phenomenon is most easily studied by numerical simulation, using specialized biased-sampling techniques to overcome the limitations imposed by the rarity of nucleation events. Here, I focus on the barrier to homogeneous ice nucleation in supercooled water, as represented by the monatomic-water model, which in the bulk exhibits a complex interplay between different ice structures. I consider various protocols to identify solidlike particles on a computer, which perform well enough for the Lennard-Jones model, and compare their respective impact on the shape and height of the nucleation barrier. It turns out that the effect is stronger on the nucleus size than on the barrier height. As a by-product of the analysis, I determine the structure of the nucleation cluster, finding that the relative amount of ice phases in the cluster heavily depends on the method used for classifying solidlike particles. Moreover, the phase which is most favored during the earlier stages of crystallization may happen, depending on the nucleation coordinate adopted, to be different from the stable polymorph. Therefore, the quality of a reaction coordinate cannot be assessed simply on the basis of the barrier height obtained. I explain how this outcome is possible and why it just points out the shortcoming of collective variables appropriate to simple fluids in providing a robust method of particle classification for monatomic water.
Over the last eight years, the Visual and Infrared Mapping Spectrometer (VIMS) aboard the Cassini orbiter has returned hyperspectral images in the 0.35-5.1 micron range of the icy satellites and rings of Saturn. These very different objects show significant variations in surface composition, roughness and regolith grain size as a result of their evolutionary histories, endogenic processes and interactions with exogenic particles. The distributions of surface water ice and chromophores, i.e. organic and non-icy materials, across the saturnian system, are traced using specific spectral indicators (spectral slopes and absorption band depths) obtained from rings mosaics and disk-integrated satellites observations by VIMS.
We examine the applicability of various model profiles for the liquid/vapor interface by X-ray reflectivities on water and ethanol and their mixtures at room temperature. Analysis of the X-ray reflecivities using various density profiles shows an error-function like profile is the most adequate within experimental error. Our finding, together with recent observations from simulation studies on liquid surfaces, strongly suggest that the capillary-wave dynamics shapes the interfacial density profile in terms of the error function.
In-situ NMR spin-lattice relaxation measurements were performed on several vapor deposited ices. The measurements, which span more than 6 orders of magnitude in relaxation times, show a complex spin-lattice relaxation pattern that is strongly dependent on the growth conditions of the sample. The relaxation patterns change from multi-timescale relaxation for samples grown at temperatures below the amorphous-crystalline transition temperature to single exponential recovery for samples grown above the transition temperature. The slow-relaxation contribution seen in cold-grown samples exhibits a temperature dependence, and becomes even slower after the sample is annealed at 200K. The fast-relaxation contribution seen in these samples, does not seem to change or disappear even when heating to temperatures where the sample is evaporated. The possibility that the fast relaxation component is linked to the microporous structures in amorphous ice samples is further examined using an environmental electron scanning microscope. The images reveal complex meso-scale microporous structures which maintain their morphology up to their desorption temperatures. These findings, support the possibility that water molecules at pore surfaces might be responsible for the fast-relaxation contribution. Furthermore, the results of this study indicate that the pore-collapse dynamics observed in the past in amorphous ices using other experimental techniques, might be effectively inhibited in samples which are grown by relatively fast vapor deposition.
Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across waters solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization, and compared the MB-pol results to experimental data. We also compared rigid vs. fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system size or the rigid vs. flexible treatment of the intramolecular degrees of freedom. These results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pols application across a large swath of waters phase diagram.