ABSTRACT
The two-dimensional Mercedes-Benz model of water has been studied by molecular simulations over a wide range of thermodynamic conditions as an attempt to locate the supercooled region where a liquid-liquid separation and, potentially, also other structures may occur. Both the correlation functions and a number of local structure factors have been used to identify different structural arrangements. These include, in addition to the hexatic phase, also the hexagon, pentagon, and quadruplet arrangements. All these structures result from the competition between the hydrogen bonding and Lennard-Jones interactions and their effect at different temperatures and pressures. Based on the obtained results, an attempt is made to sketch a (rather complex) phase diagram of the model.
ABSTRACT
Knowledge of the vapor-liquid equilibrium (VLE) properties of molten salts is important in the design of thermal energy storage systems for solar power and nuclear energy production applications. The high temperatures involved make their experimental determination problematic, and the development of both macroscopic thermodynamic correlations and predictive molecular-based methodologies are complicated by the requirement to appropriately incorporate the chemically reacting vapor-phase species. We derive a general thermodynamic-based atomistic simulation framework for molten salt VLE prediction and show its application to NaCl. Its input quantities are temperature-dependent ideal-gas free energy data for the vapor phase reactions and density and residual chemical potential data for the liquid. If these are not available experimentally, the former may be predicted using standard electronic structure software, and the latter may be predicted by means of classical atomistic simulation methodology. The framework predicts the temperature dependence of vapor pressure, coexisting phase densities, vapor phase composition, and vaporization enthalpy. It also predicts the concentrations of vapor phase species present in minor amounts (such as the free ions), quantities that are extremely difficult to measure experimentally. We furthermore use the results to obtain an approximation to the complete VLE binodal dome and the critical properties. We verify the framework for molten NaCl, for which experimentally based density and chemical potential data are available in the literature. We then apply it to the analysis of NaCl simulation data for two commonly used atomistic force fields. The framework can be readily extended to molten salt mixtures and to ionic liquids.
ABSTRACT
A general and direct computational scheme to locate the surface separating arbitrarily shaped domains made up of molecules (or any other particles) has been developed and is described and illustrated for several, both artificial and physical examples. The proposed scheme consists of two modules: (i) triangulation and (ii) assignment of simplices to domains. Three different triangulation methods are employed, viz., the Delaunay triangulation, regular triangulation, and quasi-triangulation. In the triangulated system, the assignment step is carried out in two different ways, one based on the characteristic metric of a particular triangulation procedure and the other on the concept of a touching sphere. Some of the combinations of the triangulation and assignment steps lead to methods already used by others to find interfacial or surface molecules, namely the alpha-shape-based method of Usabiaga nad Duque [Phys. Rev. E 79 (2009) 046709] and GITIM of Sega et al. [J. Chem. Phys. 138 (2013) 044110]. The resulting surface is defined not only as a discrete set of particles, but it is build up of facets of the triangulation forming a broken line in two dimensions or a polyhedral surface in three dimensions. Individual molecular layers are identified in a very straightforward manner, starting with the interfacial layer itself and proceeding into the interior of the phase. The proposed scheme is illustrated first by identifying border molecules of pre-sampled domains of several shapes in a plane and then applied to five physically meaningful examples: thin films, near critical water, liquid water slab in an electric field, liquid water at a solid wall, and water at condition of electric-field-induced jetting. Performance of the considered methods is critically assessed. Treatment of domains forming percolating clusters through periodic boundary conditions is also described along with the determination of their periodicity and dimensionality.
