Liquid polyamorphism and double criticality in a lattice gas model

Influence of disordered porous media on the anomalous properties of a simple water model

The thermodynamic, dynamic, and structural behavior of a water-like system confined in a matrix is analyzed for increasing confining geometries. The liquid is modeled by a two-dimensional associating lattice gas model that exhibits density and diffusion anomalies, similar to the anomalies present in liquid water. The matrix is a triangular lattice in which fixed obstacles impose restrictions to the occupation of the particles. We show that obstacles shorten all lines, including the phase coexistence, the critical and the anomalous lines. The inclusion of a very dense matrix not only suppresses the anomalies but also the liquid-liquid critical point.

Model of waterlike fluid under confinement for hydrophobic and hydrophilic particle-plate interaction potentials

Molecular dynamic simulations were employed to study a waterlike model confined between hydrophobic and hydrophilic plates. The phase behavior of this system is obtained for different distances between the plates and particle-plate potentials. For both hydrophobic and hydrophilic walls, there are the formation of layers. Crystallization occurs at lower temperature at the contact layer than at the middle layer. In addition, the melting temperature decreases as the plates become more hydrophobic. Similarly, the temperatures of maximum density and extremum diffusivity decrease with hydrophobicity.

Thermodynamics of supercooled water

We review the available experimental information on the thermodynamic properties of supercooled water and demonstrate the possibility of modeling these thermodynamic properties on a theoretical basis. We show that by assuming the existence of a liquid-liquid critical point in supercooled water, the theory of critical phenomena can give an accurate account of the experimental thermodynamic-property data up to a pressure of 150 MPa. In addition, we show that a phenomenological extension of the theoretical model can account for all currently available experimental data in the supercooled region, up to 400 MPa. The stability limit of the liquid state and possible coupling between crystallization and liquid-liquid separation are also discussed. It is concluded that critical-point thermodynamics describes the available thermodynamic data for supercooled water within experimental accuracy, thus establishing a benchmark for further developments in this area.

Structure and anomalous solubility for hard spheres in an associating lattice gas model

In this paper we investigate the solubility of a hard-sphere gas in a solvent modeled as an associating lattice gas. The solution phase diagram for solute at 5% is compared with the phase diagram of the original solute free model. Model properties are investigated both through Monte Carlo simulations and a cluster approximation. The model solubility is computed via simulations and is shown to exhibit a minimum as a function of temperature. The line of minimum solubility (TmS) coincides with the line of maximum density (TMD) for different solvent chemical potentials, in accordance with the literature on continuous realistic models and on the "cavity" picture.

Solution of an associating lattice-gas model with density anomaly on a Husimi lattice

We study a model of a lattice gas with orientational degrees of freedom which resemble the formation of hydrogen bonds between the molecules. In this model, which is the simplified version of the Henriques-Barbosa model, no distinction is made between donors and acceptors in the bonding arms. We solve the model in the grand-canonical ensemble on a Husimi lattice built with hexagonal plaquettes with a central site. The ground state of the model, which was originally defined on the triangular lattice, is exactly reproduced by the solution on this Husimi lattice. In the phase diagram, one gas and two liquid [high density liquid (HDL) and low density liquid (LDL)] phases are present. All phase transitions (GAS-LDL, GAS-HDL, and LDL-HDL) are discontinuous, and the three phases coexist at a triple point. A line of temperatures of maximum density in the isobars is found in the metastable GAS phase, as well as another line of temperatures of minimum density appears in the LDL phase, part of it in the stable region and another in the metastable region of this phase. These findings are at variance with simulational results for the same model on the triangular lattice, which suggested a phase diagram with two critical points. However, our results show very good quantitative agreement with the simulations, both for the coexistence loci and the densities of particles and of hydrogen bonds. We discuss the comparison of the simulations with our results.

Frustration and anomalous behavior in the Bell-Lavis model of liquid water

We have reconsidered the Bell-Lavis model of liquid water and investigated its relation to its isotropic version, the antiferromagnetic Blume-Emery-Griffiths model on the triangular lattice. Our study was carried out by means of an exact solution on the sequential Husimi cactus. We show that the ground states of both models share the same topology and that fluid phases (gas and low- and high-density liquids) can be mapped onto magnetic phases (paramagnetic, antiferromagnetic, and dense paramagnetic, respectively). Both models present liquid-liquid coexistence and several thermodynamic anomalies. This result suggests that anisotropy introduced through orientational variables play no specific role in producing the density anomaly, in agreement with a similar conclusion discussed previously following results for continuous soft core models. We propose that the presence of liquid anomalies may be related to energetic frustration, a feature common to both models.

Liquid polymorphism and density anomaly in a three-dimensional associating lattice gas

The authors investigate the phase diagram of a three-dimensional associating lattice gas (ALG) model. This model combines orientational icelike interactions and "van der Waals" that might be repulsive, representing, in this case, a penalty for distortion of hydrogen bonds. These interactions can be interpreted as two competing distances, making the connection between this model and continuous isotropic soft-core potentials. The authors present Monte Carlo studies of the ALG model showing the presence of two liquid phases, two critical points, and density anomaly.

Scaled equation of state for supercooled water near the liquid-liquid critical point

We have developed a scaled parametric equation of state to describe and predict thermodynamic properties of supercooled water. The equation of state, built on the growing evidence that the critical point of supercooled liquid-liquid water separation exists, is universal in terms of theoretical scaling fields and is shown to belong to the Ising-model class of universality. The theoretical scaling fields are postulated to be analytical combinations of the physical fields, pressure, and temperature. The equation of state enables us to accurately locate the "Widom line" (the locus of stability minima) and determine that the critical pressure is considerably lower than predicted by computer simulations.

Structural anomalies for a three dimensional isotropic core-softened potential

Using molecular dynamics simulations we investigate the structure of a system of particles interacting through a continuous core-softened interparticle potential. We found for the translational order parameter t a local maximum at a density rho(t-max) and a local minimum at rho(t-min)>rho(t-max). Between rho(t-max) and rho(t-min), the t parameter anomalously decreases upon increasing pressure. For the orientational order parameter Q(6) a maximum was observed at a density rho(t-max)<rho(Qmax)<rho(t-min). For densities between rho(Qmax) and rho(t-min), both the translational (t) and orientational (Q(6)) order parameters have anomalous behavior. We know that this system also exhibits density and diffusion anomalies. We found that the region in the pressure-temperature phase diagram of the structural anomaly englobes the region of the diffusion anomaly that is larger than the region limited by the temperature of maximum density. This cascade of anomalies (structural, dynamic, and thermodynamic) for our model has the same hierarchy as that observed for the simple point charge/extended water.

Hydrophobic Hydration in an Orientational Lattice Model

To shed light on the microscopic mechanism of hydrophobic hydration, we study a simplified lattice model for water solutions in which the orientational nature of hydrogen bonding as well as the degeneracy related to proton distribution are taken into account. Miscibility properties of the model are looked at for both polar (hydrogen bonding) and nonpolar (non-hydrogen bonding) solutes. A quasichemical solution for the pure system is reviewed and extended to include the different kinds of solute. A Monte Carlo study of our model yields a novel feature for the local structure of the hydration layer: energy correlation relaxation times for solvation water are larger than the corresponding relaxation times for bulk water. This result suggests the presence of ordering of water particles in the first hydration shell. A nonassociating model solvent, represented by a lattice gas, presents opposite behavior, indicating that this effect is a result of the directionality of the interaction. In presence of polar solutes, we find an ordered mixed pseudophase at low temperatures, indicating the possibility of closed loops of immiscibility.

Thermodynamic and dynamic anomalies for dumbbell molecules interacting with a repulsive ramplike potential

Using collision driven discrete molecular dynamics, we investigate the thermodynamics and dynamics of systems of 500 dumbbell molecules interacting by a purely repulsive ramplike discretized potential, consisting of n steps of equal size. We compare the behavior of the two systems, with n=18 and n=144 steps. Each system exhibits both thermodynamic and dynamic anomalies, a density maximum and the translational and rotational mobilities show anomalous behavior. Starting with very dense systems and decreasing the density, both mobilities first increase, reach a maximum, then decrease, reach a minimum, and finally increase; this behavior is similar to the behavior of SPC/E (Simple Point Charge-Extended) water. The regions in the pressure-temperature plane of translational and rotational mobility anomalies depend strongly on n. The product of the translational diffusion coefficient and the orientational correlation time increases with temperature, in contrast with the behavior of most liquids.

Thermodynamic and dynamic anomalies for a three-dimensional isotropic core-softened potential

Using molecular-dynamics simulations and integral equations (Rogers-Young, Percus-Yevick, and hypernetted chain closures) we investigate the thermodynamics of particles interacting with continuous core-softened intermolecular potential. Dynamic properties are also analyzed by the simulations. We show that, for a chosen shape of the potential, the density, at constant pressure, has a maximum for a certain temperature. The line of temperatures of maximum density (TMD) was determined in the pressure-temperature phase diagram. Similarly the diffusion constant at a constant temperature, D, has a maximum at a density rho(max) and a minimum at a density rho(min) < rho(max). In the pressure-temperature phase diagram the line of extrema in diffusivity is outside of the TMD line. Although this interparticle potential lacks directionality, this is the same behavior observed in simple point charge/extended water.