Core-softened fluids, water-like anomalies, and the liquid-liquid critical points

Molecular dynamics simulations are used to examine the relationship between water-like anomalies and the liquid-liquid critical point in a family of model fluids with multi-Gaussian, core-softened pair interactions. The core-softened pair interactions have two length scales, such that the longer length scale associated with a shallow, attractive well is kept constant while the shorter length scale associated with the repulsive shoulder is varied from an inflection point to a minimum of progressively increasing depth. The maximum depth of the shoulder well is chosen so that the resulting potential reproduces the oxygen-oxygen radial distribution function of the ST4 model of water. As the shoulder well depth increases, the pressure required to form the high density liquid decreases and the temperature up to which the high-density liquid is stable increases, resulting in the shift of the liquid-liquid critical point to much lower pressures and higher temperatures. To understand the entropic effects associated with the changes in the interaction potential, the pair correlation entropy is computed to show that the excess entropy anomaly diminishes when the shoulder well depth increases. Excess entropy scaling of diffusivity in this class of fluids is demonstrated, showing that decreasing strength of the excess entropy anomaly with increasing shoulder depth results in the progressive loss of water-like thermodynamic, structural and transport anomalies. Instantaneous normal mode analysis was used to index the overall curvature distribution of the fluid and the fraction of imaginary frequency modes was shown to correlate well with the anomalous behavior of the diffusivity and the pair correlation entropy. The results suggest in the case of core-softened potentials, in addition to the presence of two length scales, energetic, and entropic effects associated with local minima and curvatures of the pair interaction play an important role in determining the presence of water-like anomalies and the liquid-liquid phase transition.     

Effects of the attractive interactions in the thermodynamic, dynamic, and structural anomalies of a two length scale potential

Using molecular dynamic simulations, we study a system of particles interacting through a continuous core-softened potentials consisting of a hard core, a shoulder at closest distances, and an attractive well at further distance. We obtain the pressure-temperature phase diagram of this system for various depths of the tunable attractive well. Since this is a two length scale potential, density, diffusion, and structural anomalies are expected. We show that the effect of increasing the attractive interaction between the molecules is to shrink the region in pressure in which the density and the diffusion anomalies are present. If the attractive forces are too strong, particle will be predominantly in one of the two length scales and no density of diffusion anomaly is observed. The structural anomalous region is present for all the cases.     

Effect of pressure on liquid crystal dimers

We report high pressure investigations on a homologous series of liquid crystalline dimeric molecules in which the terminal chain length is kept constant but the length of the methylene spacer connecting the two mesogenic units is varied. We find that unlike the nematic–isotropic transition temperature and associated entropy change, there is no alternation in the slope of the nematic–isotropic phase boundary in the pressure–temperature plane as a function of the length of the spacer group in the molecule. By applying the Clausius–Clapeyron equation we conclude that the volume change at the transition should exhibit a strong odd–even effect. Measurements on the shortest homologue of the series, which is non-mesomorphic, show the expected result that the application of pressure induces mesomorphism in non-mesomorphic compounds.     

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.     

Pressure effects on the transitions between disordered phases in supercooled liquid silicon

We investigate the pressure effects on the transitions between the disordered phases in supercooled liquid silicon through Monte Carlo simulations and efficient methods to compute free energies. Our calculations, using an environment dependent interatomic potential for Si, indicate that at zero pressure the liquid-liquid phase transition, between the high density liquid and the low density liquid, occurs at a temperature 325K below melting. We found that the liquid-liquid transition temperature decreases with increasing pressure, following the liquid-solid coexistence curve. As pressure increases, the liquid-liquid coexistence curve approaches the region where the glass transition between the low density liquid and the low density amorphous takes place. Above 5 GPa, our calculations show that the liquid-liquid transition is suppressed by the glassy dynamics of the system. We also found that above 5 GPa, the glass transition temperature is lower than that at lower pressures, suggesting that under these conditions the glass transition occurs between the high density liquid and the high density amorphous.     

Effect of pressure on liquid crystal dimers

We report high pressure investigations on a homologous series of liquid crystalline dimeric molecules in which the terminal chain length is kept constant but the length of the methylene spacer connecting the two mesogenic units is varied. We find that unlike the nematic-isotropic transition temperature and associated entropy change, there is no alternation in the slope of the nematic-isotropic phase boundary in the pressure-temperature plane as a function of the length of the spacer group in the molecule. By applying the Clausius-Clapeyron equation we conclude that the volume change at the transition should exhibit a strong odd-even effect. Measurements on the shortest homologue of the series, which is non-mesomorphic, show the expected result that the application of pressure induces mesomorphism in non-mesomorphic compounds.     

Phase diagram of a two-dimensional system with anomalous liquid properties

Using Monte Carlo simulation techniques, we calculate the phase diagram for a square-shoulder square-well potential in two dimensions that has been previously shown to exhibit liquid anomalies consistent with a metastable liquid-liquid critical point. We consider the liquid, gas, and five crystal phases, and find that all the melting lines are first order, despite a small range of metastability. One melting line exhibits a temperature maximum, as well as a pressure maximum that implies inverse melting over a small range in pressure.     

On the mechanism of hydrophobic association of nanoscopic solutes

The hydration behavior of two planar nanoscopic hydrophobic solutes in liquid water at normal temperature and pressure is investigated by calculating the potential of mean force between them at constant pressure as a function of the solute-solvent interaction potential. The importance of the effect of weak attractive interactions between the solute atoms and the solvent on the hydration behavior is clearly demonstrated. We focus on the underlying mechanism behind the contrasting results obtained in various recent experimental and computational studies on water near hydrophobic solutes. The length scale where crossover from a solvent separated state to the contact pair state occurs is shown to depend on the solute sizes as well as on details of the solute-solvent interaction. We find the mechanism for attractive mean forces between the plates is very different depending on the nature of the solute-solvent interaction which has implications for the mechanism of the hydrophobic effect for biomolecules.     

Entropy, diffusivity, and structural order in liquids with waterlike anomalies

The excess entropy, defined as the difference between the entropies of the liquid and the ideal gas under identical density and temperature conditions, is studied as a function of density and temperature for liquid silica and a two-scale ramp potential, both of which are known to possess waterlike liquid state anomalies. The excess entropy for both systems is evaluated using a fairly accurate pair correlation approximation. The connection between the excess entropy and the density and diffusional anomalies is demonstrated. Using the pair correlation approximation to the excess entropy, it can be shown that if the energetically favorable local geometries in the low and high density limits have different symmetries, then a structurally anomalous regime can be defined in terms of orientational and translational order parameters, as in the case of silica and the two-scale ramp system but not for the one-scale ramp liquid. Within the category of liquids with waterlike anomalies, we show that the relationship between the macroscopic entropy and internal energy is sufficient to distinguish between those with local anisotropy and consequent open packings at low densities and those with isotropic interactions but multiple length scales. Since it is straightforward to evaluate the pair correlation entropy and internal energy from simulations or experimental data, such plots should provide a convenient means to diagnose the existence as well as type of anomalous behavior in a range of liquids, including ionic and intermetallic melts and complex fluids with ultrasoft repulsions.     

Polymorphism in simple liquids: a Gibbs ensemble Monte Carlo study

We perform Gibbs ensemble Monte Carlo (GEMC) simulations of a one-component system of hard spheres with a repulsive shoulder and an attractive well. We show the existence of two distinct liquid-gas and liquid-liquid phase equilibria. The GEMC estimate of the critical parameters, as following from an interpolation of the binodal points, is only slightly influenced by finite size effects. The liquid-gas critical temperature and pressure are lower than those of the liquid-liquid phase separation. A discussion of our findings in comparison with those of previous numerical studies is also presented.     

Critical and phase-equilibrium properties of an ab initio based potential model of methanol and 1-propanol using two-phase molecular dynamics simulations

Two-phase molecular dynamics simulations employing a Monte Carlo volume sampling method were performed using an ab initio based force field model parameterized to reproduce quantum-mechanical dimer energies for methanol and 1-propanol at temperatures approaching the critical temperature. The intermolecular potential models were used to obtain the binodal vapor-liquid phase dome at temperatures to within about 10 K of the critical temperature. The efficacy of two all-atom, site-site pair potential models, developed solely from the energy landscape obtained from high-level ab initio pair interactions, was tested for the first time. The first model was regressed from the ab initio landscape without point charges using a modified Morse potential to model the complete interactions; the second model included point charges to separate Coulombic and dispersion interactions. Both models produced equivalent phase domes and critical loci. The model results for the critical temperature, density, and pressure, in addition to the sub-critical equilibrium vapor and liquid densities and vapor pressures, are compared to experimental data. The model's critical temperature for methanol is 77 K too high while that for 1-propanol is 80 K too low, but the critical densities are in good agreement. These differences are likely attributable to the lack of multi-body interactions in the true pair potential models used here.     

Evaporation of Lennard-Jones fluids

Evaporation and condensation at a liquid/vapor interface are ubiquitous interphase mass and energy transfer phenomena that are still not well understood. We have carried out large scale molecular dynamics simulations of Lennard-Jones (LJ) fluids composed of monomers, dimers, or trimers to investigate these processes with molecular detail. For LJ monomers in contact with a vacuum, the evaporation rate is found to be very high with significant evaporative cooling and an accompanying density gradient in the liquid domain near the liquid/vapor interface. Increasing the chain length to just dimers significantly reduces the evaporation rate. We confirm that mechanical equilibrium plays a key role in determining the evaporation rate and the density and temperature profiles across the liquid/vapor interface. The velocity distributions of evaporated molecules and the evaporation and condensation coefficients are measured and compared to the predictions of an existing model based on kinetic theory of gases. Our results indicate that for both monatomic and polyatomic molecules, the evaporation and condensation coefficients are equal when systems are not far from equilibrium and smaller than one, and decrease with increasing temperature. For the same reduced temperature T/T(c), where T(c) is the critical temperature, these two coefficients are higher for LJ dimers and trimers than for monomers, in contrast to the traditional viewpoint that they are close to unity for monatomic molecules and decrease for polyatomic molecules. Furthermore, data for the two coefficients collapse onto a master curve when plotted against a translational length ratio between the liquid and vapor phase.     

Bulk and interfacial properties of binary polymer mixtures

A microscopic density functional theory is used to investigate a binary mixture of polymers, built of freely jointed tangent hard spheres. The difference in the chain length and in the segment diameter of polymers gives rise to a demixing transition. We evaluate the bulk fluid phase equilibria (binodal) and the limit of stability of a mixed state (spinodal) for selected systems, and analyze the decay of the critical packing fraction, critical mole fraction, and critical pressure with an increase of the chain length. The bulk results are subsequently used in the calculations of the density profiles across the fluid-fluid interface. The obtained profiles are smooth and do not exhibit any oscillations on the length scale of the segment diameter. Upon approaching the critical point the interfacial tension vanishes as (Deltarho)3, where Deltarho is the difference between bulk densities of one component in bulk phases rich and poor in that species. This indicates that the microscopic density functional theory applied here is of a mean-field type.     

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.     

Surface phase behavior in Gibbs monolayers of bis(ethylene glycol) mono-n-tetradecyl ether at the air-water interface

We present the adsorption kinetics and surface morphology of the adsorbed monolayers of bis(ethylene glycol) mono-n-tetradecyl ether (C14E2) by film balance and Brewster angle microscopy. A cusp point followed by a plateau region in the pressure (pi)-time (t) adsorption isotherm indicates a first-order phase transition in the coexistence region between a lower density liquid expanded (LE) phase and a higher density liquid condensed (LC) phase. A variety of condensed phase domains surrounded by the homogeneous LE phase are observed just after the appearance of the phase transition. The domains are of a spiral or striplike structure at lower temperatures. This characteristic shape of the domains is because of strong dipole-dipole repulsion between the molecules. At 18 degrees C, the domains are found to be quadrant structures. A slight increase in subphase temperature (around 1 degrees C) brings about a quadrant-to-circular shape transition in the domains. The circular domains return to quadrant structures as the subphase temperature is lowered. The domains completely disappear when the temperature is increased beyond 19 degrees C, suggesting that the critical temperature for the condensed domain formation is 19 degrees C. Above this temperature, the hypothetical surface pressure necessary for the phase transition exceeds the actual surface pressure attainable from a solution of concentration greater than or equal to the critical micelle concentration. An increase in molecular motion with increasing temperature results in a higher degree of chain flexibility. As a result, the molecules cannot accumulate in the condensed phase form when the subphase temperature is above 19 degrees C.     

Structure and phase behavior of a confined nanodroplet composed of the flexible chain molecules

A polymer density functional theory has been employed for investigating the structure and phase behaviors of the chain polymer, which is modelled as the tangentially connected sphere chain with an attractive interaction, inside the nanosized pores. The excess free energy of the chain polymer has been approximated as the modified fundamental measure-theory for the hard spheres, the Wertheim's first-order perturbation for the chain connectivity, and the mean-field approximation for the van der Waals contribution. For the value of the chemical potential corresponding to a stable liquid phase in the bulk system and a metastable vapor phase, the flexible chain molecules undergo the liquid-vapor transition as the pore size is reduced; the vapor is the stable phase at small volume, whereas the liquid is the stable phase at large volume. The wide liquid-vapor coexistence curve, which explains the wide range of metastable liquid-vapor states, is observed at low temperature. The increase of temperature and decrease of pore size result in a narrowing of liquid-vapor coexistence curves. The increase of chain length leads to a shift of the liquid-vapor coexistence curve towards lower values of chemical potential. The coexistence curves for the confined phase diagram are contained within the corresponding bulk liquid-vapor coexistence curve. The equilibrium capillary phase transition occurs at a higher chemical potential than in the bulk phase.     

Global thermodynamics of hydrophobic cavitation, dewetting, and hydration

Pure water experimental and simulation results are combined to predict the thermodynamics of cavity formation, spanning atomic to macroscopic length scales, over the entire ambient liquid temperature range. The resulting cavity equation of state is used to quantify dewetting excess contributions to cavity formation thermodynamics and construct a thermodynamic perturbation theory of hydrophobic hydration. Predictions are compared with large cavity simulations and experimental rare-gas hydration thermodynamics data (for He, Ne, Ar, Kr, Xe, and Rn). Key findings include the strong temperature dependence of the critical length scale for hydrophobic dewetting and the evaluation of fundamental solute-solvent interaction contributions to rare-gas hydration chemical potentials.     

Phase behavior of polyampholytes from charged hard-sphere chain model

A molecular thermodynamic theory is developed for polyampholytes from the coarse-grained charged hard-sphere chain model. The phase behavior of polyampholytes with variations in sequence and chain length is satisfactorily predicted by the theory, consistent with simulation results and experimental observations. At a fixed chain length, the phase envelope expands as the sequence of charge distribution becomes less random. With increasing chain length, the phase envelope expands for diblock and random polyampholytes, but shrinks for zwitterionic polyampholytes. The predicted critical temperature, density, and pressure exhibit scaling relations with chain length for all the three (diblock, random, and zwitterionic) polyampholytes.     

Thermal response of a microgel system

Recent experimental studies [Z. Wu, B. Zhou, Z.B. Hu, Phys. Rev. Lett. 90 (2003) 048304] on an uncharged aqueous poly-N-isopropylacrylamide (PNIPAM) dispersion have shown that this microgel system is sensitive to temperature. This system was also experimentally found to be modeled quite well by microgel particles interacting via a hard-sphere repulsive plus an inverse power (temperature-dependent) attractive potential. To understand theoretically this thermally responsive PNIPAM dispersion, we apply a novel approach [G.F. Wang, S.K. Lai, Phys. Rev. E 70 (2004) 051402] to calculate its thermodynamic phase diagram. Differing from the conventional method in which the boundaries of the coexisting phases are the ultimate target, the present work places emphasis on crosshatching colloidal domains which include the homogeneous phase (gas, liquid or solid), two coexisting phases and perhaps also multi-phases in coexistence. Strategically, this was done by treating the coexisting phases as one composite system whose Helmholtz free energy density is written as the sum of constituent free energy densities each of which is weighed by its respective volume proportion. We show here that by minimizing the composite system's free energy density the phase-diagram domains can all be determined in addition to the phase boundaries customarily obtained by imposing the conditions of equal pressure and equal chemical potential. Also, we present the theoretically predicted phase diagram of PNIPAM dispersion and compare it with the one observed experimentally.