Liquid Polymorphism and Density Anomaly in a Lattice Gas Model

We present a simple model for an associating liquid in which polymorphism and density anomaly are connected. Our model combines a two dimensional lattice gas with particles interacting through a soft core potential and orientational degrees of freedom represented through thermal \char`\"{}ice variables\char`\"{} . The competition between the directional attractive forces and the soft core potential leads to a phase diagram in which two liquid phases and a density anomaly are present. The coexistence line between the low density liquid and the high density liquid has a positive slope contradicting the surmise that the presence of a density anomaly implies that the high density liquid is more entropic than the low density liquid

Molecular correlations and solvation in simple fluids

We study the molecular correlations in a lattice model of a solution of a low-solubility solute, with emphasis on how the thermodynamics is reflected in the correlation functions. The model is treated in Bethe-Guggenheim approximation, which is exact on a Bethe lattice (Cayley tree). The solution properties are obtained in the limit of infinite dilution of the solute. With $h_{11}(r)$, $h_{12}(r)$, and $h_{22}(r)$ the three pair correlation functions as functions of the separation $r$ (subscripts 1 and 2 referring to solvent and solute, respectively), we find for $r \geq 2$ lattice steps that $h_{22}(r)/h_{12}(r) \equiv h_{12}(r)/h_{11}(r)$. This illustrates a general theorem that holds in the asymptotic limit of infinite $r$. The three correlation functions share a common exponential decay length (correlation length), but when the solubility of the solute is low the amplitude of the decay of $h_{22}(r)$ is much greater than that of $h_{12}(r)$, which in turn is much greater than that of $h_{11}(r)$. As a consequence the amplitude of the decay of $h_{22}(r)$ is enormously greater than that of $h_{11}(r)$. The effective solute-solute attraction then remains discernible at distances at which the solvent molecules are essentially no longer correlated, as found in similar circumstances in an earlier model. The second osmotic virial coefficient is large and negative, as expected. We find that the solvent-mediated part $W(r)$ of the potential of mean force between solutes, evaluated at contact, $r=1$, is related in this model to the Gibbs free energy of solvation at fixed pressure, $\Delta G_p^*$, by $(Z/2) W(1) + \Delta G_p^* \equiv p v_0$, where $Z$ is the coordination number of the lattice, $p$ the pressure, and $v_0$ the volume of the cell associated with each lattice site. A large, positive $\Delta G_p^*$ associated with the low solubility is thus reflected in a strong attraction (large negative $W$ at contact), which is the major contributor to the second osmotic virial coefficient. In this model, the low solubility (large positive $\Delta G_p^*$) is due partly to an unfavorable enthalpy of solvation and partly to an unfavorable solvation entropy, unlike in the hydrophobic effect, where the enthalpy of solvation itself favors high solubility, but is overweighed by the unfavorable solvation entropy.Comment: 9 pages, 2 figure

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 pressure. For the orientational order parameter, $Q_6$, was observed a maximum 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 anomaly. 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 of that one observed for the SPC/E water.Comment: 19 pages, 8 figure

Diffusion anomaly and dynamic transitions in the Bell-Lavis water model

In this paper we investigate the dynamic properties of the minimal Bell-Lavis (BL) water model and their relation to the thermodynamic anomalies. The Bell-Lavis model is defined on a triangular lattice in which water molecules are represented by particles with three symmetric bonding arms interacting through van der Waals and hydrogen bonds. We have studied the model diffusivity in different regions of the phase diagram through Monte Carlo simulations. Our results show that the model displays a region of anomalous diffusion which lies inside the region of anomalous density, englobed by the line of temperatures of maximum density (TMD). Further, we have found that the diffusivity undergoes a dynamic transition which may be classified as fragile-to-strong transition at the critical line only at low pressures. At higher densities, no dynamic transition is seen on crossing the critical line. Thus evidence from this study is that relation of dynamic transitions to criticality may be discarded

Anomalies in a waterlike model confined between plates

Using molecular dynamic simulations we study a waterlike model confined between two fixed hydrophobic plates. The system is tested for density, diffusion and structural anomalous behavior and compared with the bulk results. Within the range of confining distances we had explored we observe that in the pressure-temperature phase diagram the temperature of maximum density (TMD line), the temperature of maximum and minimum diffusion occur at lower temperatures when compared with the bulk values. For distances between the two layers below a certain threshold ,$d\le d_c$, only two layers of particles are formed, for $d\ge d_c$ three or more layers are formed. In the case of three layers the central layer stays liquid while the contact layers crystallize. This result is in agreement with simulations for atomistic models