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

Hierarchy of anomalies in the two-dimensional Mercedes-Benz model of water

We investigate by Monte Carlo simulations density, diffusion, and structural anomalies of the simple two-dimensional Mercedes-Benz (MB) model of water, which is a very simple toy model for explaining the origin of water properties. MB water molecules are modeled as two-dimensional Lennard-Jones disks, with three orientation-dependent hydrogen-bonding arms, arranged as in the MB logo. The model is in a way also a variance of silica-like models. Beside the known thermodynamic anomaly for the model we also found diffusion and structural anomalies and map out the cascade of density, structural, pair entropy, and diffusivity anomalies for MB model. The orientational order parameters with three and six-fold symmetry were determined and maximum for each one observed. The anomalies occur in hierarchy order, which is a slight variation of the hierarchy order in real water. The diffusion anomaly region is the innermost in the hierarchy while for water it is the density anomaly region.

Waterlike glass polyamorphism in a monoatomic isotropic Jagla model

We perform discrete-event molecular dynamics simulations of a system of particles interacting with a spherically-symmetric (isotropic) two-scale Jagla pair potential characterized by a hard inner core, a linear repulsion at intermediate separations, and a weak attractive interaction at larger separations. This model system has been extensively studied due to its ability to reproduce many thermodynamic, dynamic, and structural anomalies of liquid water. The model is also interesting because: (i) it is very simple, being composed of isotropically interacting particles, (ii) it exhibits polyamorphism in the liquid phase, and (iii) its slow crystallization kinetics facilitate the study of glassy states. There is interest in the degree to which the known polyamorphism in glassy water may have parallels in liquid water. Motivated by parallels between the properties of the Jagla potential and those of water in the liquid state, we study the metastable phase diagram in the glass state. Specifically, we perform the computational analog of the protocols followed in the experimental studies of glassy water. We find that the Jagla potential calculations reproduce three key experimental features of glassy water: (i) the crystal-to-high-density amorphous solid (HDA) transformation upon isothermal compression, (ii) the low-density amorphous solid (LDA)-to-HDA transformation upon isothermal compression, and (iii) the HDA-to-very-high-density amorphous solid (VHDA) transformation upon isobaric annealing at high pressure. In addition, the HDA-to-LDA transformation upon isobaric heating, observed in water experiments, can only be reproduced in the Jagla model if a free surface is introduced in the simulation box. The HDA configurations obtained in cases (i) and (ii) are structurally indistinguishable, suggesting that both processes result in the same glass. With the present parametrization, the evolution of density with pressure or temperature is remarkably similar to the corresponding experimental measurements on water. Our simulations also suggest that the Jagla potential may reproduce features of the HDA-VHDA transformations observed in glassy water upon compression and decompression. Snapshots of the system during the HDA-VHDA and HDA-LDA transformations reveal a clear segregation between LDA and HDA but not between HDA and VHDA, consistent with the possibility that LDA and HDA are separated by a first order transformation as found experimentally, whereas HDA and VHDA are not. Our results demonstrate that a system of particles with simple isotropic pair interactions, a Jagla potential with two characteristic length scales, can present polyamorphism in the glass state as well as reproducing many of the distinguishing properties of liquid water. While most isotropic pair potential models crystallize readily on simulation time scales at the low temperatures investigated here, the Jagla potential is an exception, and is therefore a promising model system for the study of glass phenomenology.

Breakdown of excess entropy scaling for systems with thermodynamic anomalies

This paper presents a simulation study of the applicability of the Rosenfeld entropy scaling to the systems which cannot be approximated by the effective hard spheres. Three systems are studied: the Herzian spheres, the Gauss core model, and a soft repulsive shoulder potential. These systems demonstrate diffusion anomalies at low temperatures: the diffusion coefficient increases with increasing density or pressure. It is shown that for the first two systems belonging to a class of bounded potentials, the Rosenfeld scaling formula is valid only in the infinite-temperature limit where there are no anomalies. For the soft repulsive shoulder potential, the scaling formula is valid already at sufficiently low temperatures, however, out of the anomaly range.

Cluster formation, waterlike anomalies, and re-entrant melting for a family of bounded repulsive interaction potentials

We introduce a family of bounded repulsive potentials, which we call the cut ramp potential, obtained by cutting a linear ramp potential at different heights. We find that for the uncut ramp potential the system shows a region of anomalous re-entrant melting (a negative slope of the melting line in the temperature-pressure phase diagram), with waterlike anomalies in the same pressure range. At high pressure the melting line recovers a positive slope, a feature that we associate with the formation of clusters of particles separated by a more or less density-independent distance, the cluster separation, which is approximately equal to the ramp width sigma1. As the ramp is cut at lower and lower heights, the region of anomalous behavior shrinks and eventually disappears while at the same time the formation of clusters becomes more favored, as it is energetically less unfavorable for particles to "climb up" the ramp. We relate the occurrence of anomalous behavior to the reduced efficacy of the soft repulsive length scale with increasing pressure. The clustering phenomenon partially restores this efficacy, giving rise to an approximately constant distance sigma1 between the clusters. Our results may be useful to better understand the phase behavior of macromolecules as well as that of substances with nondirectional interactions that are capable of displaying liquid polymorphism.

Heterogeneities in confined water and protein hydration water

We report recent efforts to understand a broad range of experiments on confined water and protein hydration water, many initiated by a collaboration between workers at the University of Messina and MIT-the editors of this special issue. Preliminary calculations are not inconsistent with one tentative interpretation of these experiments as resulting from the system passing from the high-temperature high-pressure 'HDL' side of the Widom line (where the liquid might display non-Arrhenius behavior) to the low-temperature low-pressure 'LDL' side of the Widom line (where the liquid might display Arrhenius behavior). The Widom line-defined to be the line in the pressure-temperature plane where the correlation length has its maximum-arises if there is a critical point. Hence, interpreting the Messina-MIT experiments in terms of a Widom line is of potential relevance to testing, experimentally, the hypothesis that water displays a liquid-liquid critical point.

Unusual phase behavior of one-component systems with two-scale isotropic interactions

We study the phase behavior of systems of particles interacting through pair potentials with a hard core plus a soft repulsive component. We consider several different forms of soft repulsion, including a square shoulder, a linear ramp and a quasi-exponential tail. The common feature of these potentials is the presence of two repulsive length scales, which may be the origin of unusual phase behaviors such as polyamorphism both in the equilibrium liquid phase and in the glassy state, water-like anomalies in the liquid state and anomalous melting at very high pressures.

Generalized Rosenfeld scalings for tracer diffusivities in not-so-simple fluids: mixtures and soft particles

Rosenfeld [Phys. Rev. A 15, 2545 (1977)] originally noticed that casting the transport coefficients of simple monatomic equilibrium fluids in a specific dimensionless form makes them approximately single-valued functions of excess entropy. This observation has predictive value because, while the transport coefficients of dense fluids can be difficult to estimate from first principles, the excess entropy can often be accurately predicted from liquid-state theory. In this work, we use molecular simulations to investigate whether Rosenfeld's observation is a special case of a more general scaling law relating the tracer diffusivities of particles in mixtures to the excess entropy. Specifically, we study the tracer diffusivities, static structure, and thermodynamic properties of a variety of one- and two-component model fluid systems with either additive or nonadditive interactions of the hard-sphere or Gaussian-core form. The results of the simulations demonstrate that the effects of mixture concentration and composition, particle-size asymmetry and additivity, and strength of the interparticle interactions in these fluids are consistent with an empirical scaling law relating the excess entropy to a dimensionless (generalized Rosenfeld) form of tracer diffusivity, which we introduce here. The dimensionless form of the tracer diffusivity follows from knowledge of the intermolecular potential and the transport/thermodynamic behavior of fluids in the dilute limit. The generalized Rosenfeld scaling requires less information and provides more accurate predictions than either Enskog theory or scalings based on the pair-correlation contribution to the excess entropy. As we show, however, it also suffers from some limitations especially for systems that exhibit significant decoupling of individual component tracer diffusivities.

A monatomic system with a liquid-liquid critical point and two distinct glassy states

We study the glass transition (GT) in a model system that exhibits the presence of more than one liquid or glassy state ("polyamorphism") using molecular dynamics simulations. We study the Jagla model [E. A. Jagla, J. Chem. Phys. 111, 8980 (1999)], a two-scale spherically symmetric ramp potential with both attractive and repulsive interactions. The Jagla model is particularly interesting since, depending on its parametrization, it predicts two phases ("polyamorphism") not only in the glassy state but also in equilibrium as a liquid-liquid phase transition (LLPT). The Jagla model may also be useful in understanding a recent observation of polyamorphism in metallic glasses containing cerium. We use a parametrization for which crystallization can be avoided and the GT and LLPT are clearly separated, providing a unique opportunity to study the effects of the LLPT on the GT. We follow the experimental protocol employed in the classical differential scanning calorimetry experiments used to characterize the GT, cooling and heating the system through the GT and calculating the constant-pressure specific heat C(P) and the thermal expansion coefficient alpha(P). At pressures below and well above the LLPT, the same basic GT phenomenology of metallic glasses is observed, i.e., a single peak in C(P) (typical of ergodicity restoration) occurs upon heating across the GT. At pressures above the LLPT, a second peak in C(P) develops at higher temperature above the GT. This second peak in C(P) arises from the presence of a Widom line T(W) defined as the locus of maximum correlation length in the one-phase region above the liquid-liquid critical point (LLCP). The behavior of alpha(P) is different across the GT and Widom line. Near the GT temperature T(g), alpha(P) displays a small peak upon heating, which makes a negligible contribution to the C(P) peak. On the other hand, near T(W), alpha(P) displays a much larger peak, which makes a substantial contribution to the C(P) peak at higher temperature. We find that T(g) is almost independent of pressure for each of the two coexisting liquids, but shows an apparent discontinuity upon crossing the LLPT line, to a lower value for the higher-entropy phase. We compare the entropies of both phases, and the corresponding temperature dependencies, with those of the crystal phase. We also study the dependence of the GT on heating rate and find that for pressures below the LLCP, slow heating results in crystallization, as occurs in laboratory experiments. Regarding the thermal expansion properties of the Jagla model, we study the interplay of the density minimum recently observed in confined water and the GT.

Anomalous structure and dynamics of the Gaussian-core fluid

It is known that there are thermodynamic states for which the Gaussian-core fluid displays anomalous properties such as expansion upon isobaric cooling (density anomaly) and increased single-particle mobility upon isothermal compression (self-diffusivity anomaly). Here, we investigate how temperature and density affect its short-range translational structural order, as characterized by the two-body excess entropy. We find that there is a wide range of conditions for which the short-range translational order of the Gaussian-core fluid decreases upon isothermal compression (structural order anomaly). As we show, the origin of the structural anomaly is qualitatively similar to that of other anomalous fluids (e.g., water or colloids with short-range attractions) and is connected to how compression affects static correlations at different length scales. Interestingly, we find that the self-diffusivity of the Gaussian-core fluid obeys a scaling relationship with the two-body excess entropy that is very similar to the one observed for a variety of simple liquids. One consequence of this relationship is that the state points for which structural, self-diffusivity, and density anomalies of the Gaussian-core fluid occur appear as cascading regions on the temperature-density plane; a phenomenon observed earlier for models of waterlike fluids. There are, however, key differences between the anomalies of Gaussian-core and waterlike fluids, and we discuss how those can be qualitatively understood by considering the respective interparticle potentials of these models. Finally, we note that the self-diffusivity of the Gaussian-core fluid obeys different scaling laws depending on whether the two-body or total excess entropy is considered. This finding, which deserves more comprehensive future study, appears to underscore the significance of higher-body correlations for the behavior of fluids with bounded interactions.

Waterlike hierarchy of anomalies in a continuous spherical shouldered potential

We investigate by molecular dynamics simulations a continuous isotropic core-softened potential with attractive well in three dimensions, introduced by Franzese [J. Mol. Liq. 136, 267 (2007)], that displays liquid-liquid coexistence with a critical point and waterlike density anomaly. Besides the thermodynamic anomalies, here we find diffusion and structural anomalies. The anomalies, not observed in the discrete version of this model, occur with the same hierarchy that characterizes water. We discuss the differences in the anomalous behavior of the continuous and the discrete model in the framework of the excess entropy, calculated within the pair correlation approximation.

Structural anomalies of fluids: origins in second and higher coordination shells

Compressing or cooling a fluid typically enhances its static interparticle correlations. However, there are notable exceptions. Isothermal compression can reduce the translational order of fluids that exhibit anomalous waterlike trends in their thermodynamic and transport properties, while isochoric cooling (or strengthening of attractive interactions) can have a similar effect on fluids of particles with short-range attractions. Recent simulation studies by Yan [Phys. Rev. E 76, 051201 (2007)] on the former type of system and Krekelberg [J. Chem. Phys. 127, 044502 (2007)] on the latter provide examples where such structural anomalies can be related to specific changes in second and more distant coordination shells of the radial distribution function. Here, we confirm the generality of this microscopic picture through analysis, via molecular simulation and integral equation theory, of coordination shell contributions to the two-body excess entropy for several related model fluids which incorporate different levels of molecular resolution. The results suggest that integral equation theory can be an effective and computationally inexpensive tool for assessing, based on the pair potential alone, whether new model systems are good candidates for exhibiting structural (and hence thermodynamic and transport) anomalies.

Connection of translational and rotational dynamical heterogeneities with the breakdown of the Stokes-Einstein and Stokes-Einstein-Debye relations in water

We study the Stokes-Einstein (SE) and the Stokes-Einstein-Debye (SED) relations, Dt=kBT/6pietaR and Dr=kBT/8pietaR3, where Dt and Dr are the translational and rotational diffusivity, respectively, T is the temperature, eta the viscosity, kB the Boltzmann constant, and R the "molecular" radius. Our results are based on molecular dynamics simulations of the extended simple point charge model of water. We find that both the SE and SED relations break down at low temperature. To explore the relationship between these breakdowns and dynamical heterogeneities (DHs), we also calculate the SE and SED relations for subsets of the 7% "fastest" and 7% "slowest" molecules. We find that the SE and SED relations break down in both subsets, and that the breakdowns occur on all scales of mobility. Thus these breakdowns appear to be generalized phenomena, in contrast with a view where only the most mobile molecules are the origin of the breakdown of the SE and SED relations, embedded in an inactive background where these relations hold. At low temperature, the SE and SED relations in both subsets of molecules are replaced with "fractional" SE and SED relations, Dt approximately (tau/T)-xit and Dr approximately (tau/T)-xir, where xit approximately 0.84(<1) and xir approximately 0.75(<1). We also find that there is a decoupling between rotational and translational motion, and that this decoupling occurs in both the fastest and slowest subsets of molecules. Further, we find that, the decoupling increases upon cooling, but that the probability of a molecule being classified as both translationally and rotationally fastest also increases. To study the effect of time scale for SE and SED breakdown and decoupling, we introduce a time-dependent version of the SE and SED relations, and a time-dependent function that measures the extent of decoupling. Our results suggest that both the decoupling and SE and SED breakdowns originate at a time scale corresponding to the end of the cage regime, when diffusion starts. This is also the time scale when the DHs are more relevant. Our work also demonstrates that selecting DHs on the basis of translational or rotational motion more strongly biases the calculation of diffusion constants than other dynamical properties such as relaxation times.

Phase behavior of attractive and repulsive ramp fluids: Integral equation and computer simulation studies

Using computer simulations and a thermodynamically self-consistent integral equation we investigate the phase behavior and thermodynamic anomalies of a fluid composed of spherical particles interacting via a two-scale ramp potential (a hard core plus a repulsive and an attractive ramp) and the corresponding purely repulsive model. Both simulation and integral equation results predict a liquid-liquid demixing when attractive forces are present, in addition to a gas-liquid transition. Furthermore, a fluid-solid transition emerges in the neighborhood of the liquid-liquid transition region, leading to a phase diagram with a somewhat complicated topology. This solidification at moderate densities is also present in the repulsive ramp fluid, but in this case inhibits the fluid-fluid separation.

Theoretical prediction of multiple fluid-fluid transitions in monocomponent fluids

The authors use the analytical equation of state obtained by the discrete perturbation theory [A. L. Benavides and A. Gil-Villegas, Mol. Phys. 97, 1225 (1999)] to study the phase diagram of fluids with discrete spherical potentials formed by a repulsive square-shoulder plus an attractive square-well interaction (SS+SW). This interaction is characterized by the usual energy and size parameters plus three dimensionless parameters: two of them measuring the widths of the SS and the SW and the third the relative height of the SS. The matter of interest is that, for certain values of the interaction parameters, the SS+SW systems exhibit more than one first-order fluid-fluid transition. The evidence that several real substances (such as water, phosphorus, carbon, and silica, among others) exhibit an extra liquid-liquid transition has drawn interest into the study of interactions responsible for this behavior. The simple SS+SW fluid is one of the systems that, in spite of being spherically symmetric, shows multiple fluid-fluid transitions. In this work the authors investigate systematically the effect on the phase diagram of varying the interaction parameters. The use of an analytical free-energy equation gives a clear thermodynamic picture of the emergence of different types of critical points, throwing new light on the phase behavior of these fluids and thus clarifying previous results obtained by other techniques. The interplay of attractive and repulsive forces with several scale lengths produces very rich phase diagrams, including cases with three critical points. The region of the interaction-parameter space where multiple critical points appear is mapped for various families of interactions.

Local Structure and Thermodynamics of a Core-Softened Potential Fluid: Theory and Simulation

Phase behavior and structural properties of homogeneous and inhomogeneous core-softened (CS) fluid consisting of particles interacting via the potential, which combines the hard-core repulsion and double attractive well interaction, are investigated. The vapour-liquid coexistence curves and critical points for various interaction ranges of the potential are determined by discrete molecular dynamics simulations to provide guidance for the choice of the bulk density and potential parameters for the study of homogeneous and inhomogeneous structures. Spatial correlations in the homogeneous CS system are studied by the Ornstein-Zernike integral equation in combination with the modified hypernetted chain (MHNC) approximation. The local structure of CS fluid subjected to diverse external fields maintaining the equilibrium with the bulk CS fluid are studied on the basis of a recently proposed third order+second order perturbation density functional approximation (DFA). The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. Reasonable agreement between the results of both methods proves that the DFA theory applied in this work is a convenient theoretical tool for the investigation of the CS fluid, which is practically applicable for modeling numerous real systems.

Thermodynamics and dynamics of the two-scale spherically symmetric Jagla ramp model of anomalous liquids

Using molecular dynamics simulations, we study the Jagla model of a liquid which consists of particles interacting via a spherically symmetric two-scale potential with both repulsive and attractive ramps. This potential displays anomalies similar to those found in liquid water, namely expansion upon cooling and an increase of diffusivity upon compression, as well as a liquid-liquid (LL) phase transition in the region of the phase diagram accessible to simulations. The LL coexistence line, unlike in tetrahedrally coordinated liquids, has a positive slope, because of the Clapeyron relation, corresponding to the fact that the high density phase (HDL) is more ordered than low density phase (LDL). When we cool the system at constant pressure above the critical pressure, the thermodynamic properties rapidly change from those of LDL-like to those of HDL-like upon crossing the Widom line. The temperature dependence of the diffusivity also changes rapidly in the vicinity of the Widom line, namely the slope of the Arrhenius plot sharply increases upon entering the HDL domain. The properties of the glass transition are different in the two phases, suggesting that the less ordered phase is fragile, while the more ordered phase is strong, which is consistent with the behavior of tetrahedrally coordinated liquids such as water silica, silicon, and BeF2.