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 confinement on the liquid-liquid phase transition of supercooled water

We report on an observation of the phase transition between two liquid phases of supercooled confined water in simulations. The temperature of the liquid-liquid transition of water at zero pressure slightly decreases due to confinement in the hydrophobic pore. The hydrophilic confinement affects this temperature in the opposite direction and shifts the critical point of the liquid-liquid transition to a higher pressure. As a result, in a strongly hydrophilic pore the liquid-liquid phase transition becomes continuous at zero pressure, indicating the shift of its critical point from negative to a positive pressure. These findings indicate that experimental studies of water confined in the pores of various hydrophobicity/hydrophilicity may clarify the location of the liquid-liquid critical point of bulk water.     

The putative liquid-liquid transition is a liquid-solid transition in atomistic models of water

We use numerical simulation to examine the possibility of a reversible liquid-liquid transition in supercooled water and related systems. In particular, for two atomistic models of water, we have computed free energies as functions of multiple order parameters, where one is density and another distinguishes crystal from liquid. For a range of temperatures and pressures, separate free energy basins for liquid and crystal are found, conditions of phase coexistence between these phases are demonstrated, and time scales for equilibration are determined. We find that at no range of temperatures and pressures is there more than a single liquid basin, even at conditions where amorphous behavior is unstable with respect to the crystal. We find a similar result for a related model of silicon. This result excludes the possibility of the proposed liquid-liquid critical point for the models we have studied. Further, we argue that behaviors others have attributed to a liquid-liquid transition in water and related systems are in fact reflections of transitions between liquid and crystal.     

乙腈-水的液液、固液相平衡和过量Gibbs自由能

测制了xCH~3CN+(1-x)H~2O的液液、固液平衡相图, 此系液液分层的简单低共熔混合物类相图。低共熔点温度为227.44K, 组成x=0.955。最高临界溶解温度为271.0K, 临界组成x=0.35, 临界指数n=2.64。两液相与冰的平衡温度为263.07K。计算出体系在263.07K的过量Gibbs自由能G~m^E, 液液分层时G~m^E的最大值为1174J·mol^-^1。     

Solid-liquid transitions of sodium chloride at high pressures

We investigate solid-liquid transitions in NaCl at high pressures using molecular dynamics simulations, including the melting curve and superheating/supercooling as well as solid-solid and liquid-liquid transitions. The first-order B1-B2 (NaCl-CsCl type) transition in solid is observed at high pressures besides continuous liquid structure transitions, which are largely analogous to the B1-B2 transition in solid. The equilibrium melting temperatures (T(m)) up to megabar pressure are obtained from the solid-liquid coexistence technique and the superheating-supercooling hysteresis method. Lindemann's vibrational and Born's mechanical instabilities are found upon melting. The Lindemann frequency is calculated from the vibrational density of states. The Lindemann parameter (fractional root-mean-squared displacement) increases with pressure and approaches a constant asymptotically, similar to the Lennard-Jones system. However, the Lindemann melting relation holds for both B1 and B2 phases to high accuracy as for the Lennard-Jonesium. The B1 and B2 NaCl solids can be superheated by 0.18T(m) and 0.24T(m), and the NaCl liquid, supercooled by 0.22T(m) and 0.32T(m), respectively, at heating or cooling rates of 1 K/s and 1 K/ps. The amount of maximum superheating or supercooling and its weak pressure dependence observed for NaCl are in accord with experiments on alkali halides and with simulations on the Lennard-Jones system and Al.     

Liquid-liquid transition in supercooled water investigated by interaction with LiCl and Xe

The hypothesis that supercooled water consists of two distinct liquid phases has been explored on the basis of their ability to hydrate nonpolar (Xe) and electrolytic (LiCl) species. Xe incorporated in the bulk of amorphous solid water survives in the deeply supercooled regime above the glass-transition temperature of 136 K and is finally dehydrated at 165 K, whereas LiCl dissolves only in the liquid phase appearing above 165 K. The second liquid phase connects with normal water as inferred from high (poor) solubility of LiCl (Xe). This result also suggests that decoupling of translational diffusion and viscosity in the deeply supercooled regime is caused by domain structures of the two liquid phases formed during a possible liquid-liquid transition.     

Two liquid phases of water in the deeply supercooled region and their roles in crystallization and formation of LiCl solution

The properties of supercooled liquid water and the mechanism of crystallization in it were investigated using time-of-flight secondary ion mass spectrometry and reflection absorption infrared spectroscopy. The self-diffusion of the water molecules commences at 136 K, and then the liquid-liquid phase transition occurs at 160-165 K. The latter is evidenced not only by the occurrence of fluidity but also by the formation of a LiCl solution. The infrared absorption band also changes drastically above 160 K due to crystallization of water (on the Au film) and the formation of LiCl solution (on the LiCl film). The immediate crystallization and dissolution of LiCl are thought to be characteristic of normal water that is created in a deeply supercooled region, indicating that viscous liquid water (T > 136 K) is transformed into supercooled liquid water at around 160 K. The crystallization kinetics is different between these two phases because the former (latter) involves nuclear growth (spontaneous nucleation). Without nuclei, crystallization is quenched below 160 K in the present experiment. It is suggested that the viscous liquid phase coexists at the surface or grain boundaries of metastable ice Ic.     

Phosphorus: first principle simulation of a liquid-liquid phase transition

We report a Car-Parrinello molecular dynamics study of the liquid-liquid phase transition in phosphorus. We employed a gradient corrected density functional (B-LYP) to describe the electronic structure and performed simulations at constant pressure. Upon increasing pressure we observed, along the 1500 K isotherm, a structural transition converting the molecular P4 liquid into an atomic liquid with a network structure. Our calculations suggest this transition to be first order with a discontinuous density increase accompanied by an insulator into metal transition. The transition pressure is significantly higher than obtained by employing the less accurate local density functional (LDA) [Morishita, Phys. Rev. Lett. 87, 105701 (2001)], which matches the experimental value for the pressure. We argue why the LDA result should be considered fortuitous. The change of the calculated structure factor upon the transition shows the same trend as experimentally observed. Analysis of the structural changes during the phase transition revealed that a chain of linked and opened up ("butterfly") P4 molecules may serve as a seed triggering the transition from the molecular to the network phase.     

Effect of temperature on the surface phase behavior of n-hexadecyl dihydrogen phosphate in adsorption layers at the air-water interface

We present the adsorption kinetics and the surface phase behavior of n-hexadecyl dihydrogen phosphate (n-HDP) at the air-water interface by film balance and Brewster angle microscopy (BAM). A phase diagram, which shows a triple point at about 25.8 degrees C, is constructed by measuring the surface pressure (pi)-time (t) adsorption isotherms. Below 25.8 degrees C, each of the pi-t curves shows a plateau at about zero surface pressure indicating the existence of a first-order phase transition. The BAM observation confirms the order of this phase transition by presenting two-surface phases during this plateau. However, the BAM observation also shows clearly another second-order phase transition from an isotropic phase to a mosaic-textured liquid condensed (LC) phase. The initial phase is a gas (G) phase. Considering the peculiarity of the middle phase, we suggest this phase as an intermediate (I) phase. Above the triple point, the pi-t curves predict the existence of two-step first-order phase transitions. Similar to the results at lower temperatures, the BAM images show two-surface phases during these first-order phase transitions together with a second-order phase transition from an isotropic phase to an LC phase. These transitions are classified as a first-order G-LE (liquid expanded) phase transition, which is followed by another first-order LE-I phase transition. The second-order phase transition is an I-LC phase transition. Contrary to these results, at 36 degrees C both the pi-t measurements and the BAM observation present only two first-order phase transitions, which are G-LE at zero surface pressure and LE-LC transition at higher surface pressure. The shape of the domains during the main transitions shows a peculiar change from a circular at 20 degrees C to an elongated at 24 degrees C and finally to a circular shape at 36 degrees C. Such a change in the domain shapes has been explained considering the dehydration effect at higher temperatures as well as the nature of phases.     

Free energy of the solid C60 fullerene orientational order-disorder transition

The free energies of the orientationally ordered crystal phase of C60 at low temperatures and the disordered crystal phase at high temperatures are calculated to an accuracy of +/-0.05 kJ/mol using the expanded ensemble Monte Carlo method with the potential model of Sprik et al. [J. Phys. Chem. 96, 2027 (1992)]. The order-disorder transition temperature at zero pressure is determined directly from these free energies, and is found to be consistent with the abrupt changes in configurational energy and unit cell size also found in simulation. A modification of the potential results in predictions of the transition temperature of 257 K and the entropy change of 18.1 J/mol K at this transition, which are in good agreement with the experimental values of 260 K and 19 J/mol K, respectively. The orientational distinguishability in the ordered phase and the indistinguishability in the disordered phase lead to a contribution to the entropy difference of k ln 60, with 60 being the symmetry number of C60. This quantum mechanical correction is important for the accurate prediction of the phase transition properties of the C60 crystals.     

Thermodynamics of the liquid states of Langmuir monolayers

The liquid states and the liquid-liquid equilibrium of surfactant molecules forming an interphase between air and water have been considered using Monte Carlo computer simulations. Specifically, the expanded and compressed liquid phases observed for surfactant molecules were characterized as a function of pressure and temperature. Simple modified beadlike potentials were implemented in order to describe the interparticle forces between the hydrophobic and hydrophilic portions of surfactant molecules at the air/water interface. A simulation box was defined such that the monolayer was exposed to an externally applied lateral pressure in a modified isothermal-isobaric ensemble, whereas the water bath was modeled in a canonical ensemble. The simulation resembles the experimental setup used to measure lateral pressure (Pi) versus area isotherms obtained with Langmuir troughs. The applied lateral pressure-surface area phase diagram clearly showed the coexistence of the expanded and compressed liquid phases within certain temperature and pressure ranges. Distribution functions of distances and enthalpies for the monolayer were computed to clearly identify each liquid phase and the coexistence region.     

Phase-ordering kinetics of the liquid-liquid transition in single-component molecular liquids

Recently it has been revealed that even a single-component liquid can have more than two liquid states. The transition between these liquid states is called the "liquid-liquid transition." Most known liquid-liquid transitions occur at temperatures and pressures which are difficult to access experimentally, so the physical nature of the transition, particularly the kinetics, has remained elusive. However, the recent discovery of liquid-liquid transitions in molecular liquids opens up a possibility to study the kinetics in detail. Here, we report the first phase field simulation on the kinetics of a liquid-liquid transition and its direct comparison with experimental results of the molecular liquids. Both nucleation-growth-type and spinodal-decomposition-type liquid-liquid transformation observed experimentally are well reproduced by numerical simulation based on a two-order-parameter model of liquid that regards the liquid-liquid transition as the cooperative formation of locally favored structures. Thus, phase field calculations may allow us to predict the kinetics of liquid-liquid transitions and the resulting spatiotemporal change of various physical properties of the liquid, such as density and refractive index.     

On the use of Bennett's acceptance ratio method in multi-canonical-type simulations

A common strategy for mapping coexistence curves is to employ multi-canonical (MUCA) sampling to simulate along a macrostate path connecting two phases. Central to this approach is the task of accurately calculating the importance weights used in the MUCA procedure, which are needed for both effective sampling and accurate determination of phase boundaries. The purpose of this study is to develop a strategy for determining the importance weights that is built upon Bennett's optimized acceptance ratio method. This approach is shown to be closely related to transition matrix schemes, and is used to compute the vapor-liquid equilibrium of a Lennard-Jones fluid and the liquid-liquid equilibrium of a n-hexane/n-perfluorohexane mixture. For the Lennard-Jones system, the importance weights as a function of the number of particles "N" (at fixed temperature and volume) are obtained by using Bennett's method to estimate free energy differences between N and N+1 particle systems over the desired range of N values. In this application, the method is found to perform slightly better than a related transition matrix scheme. For the n-hexane/n-perfluorohexane liquid mixture, the method is designed to obtain weights as a function of composition (for fixed temperature, pressure, and total number of particles); in this case, the method is found to outperform the Gibbs ensemble approach.     

Dissociation behavior of C2H6 hydrate at temperatures below the ice point: melting to liquid water followed by ice nucleation

The dissociation of C(2)H(6) hydrate particles by slow depressurization at temperatures slightly below the ice melting point was studied using optical microscopy and Raman spectroscopy. Visual observations and Raman measurements revealed that ethane hydrates can be present as a metastable state at pressures lower than the dissociation pressures of the three components: ice, hydrate, and free gas. However, they decompose into liquid water and gas phases once the system pressure drops to the equilibrium boundary for supercooled water, hydrate, and free gas. Structural analyses of obtained Raman spectra indicate that structures of the metastable hydrates and liquid water from the hydrate decay are fundamentally identical to those of the stable hydrates and supercooled water without experience of the hydration. These results imply a considerably high energy barrier for the direct hydrate-to-ice transition. Water solidification, probably induced by dynamic nucleation, was also observed during melting.     

On the role of frustration on the glass transition and polyamorphism of mesoscopically heterophase liquids

The model of heterophase fluctuations is developed accounting frustration of the mesoscopic solidlike fluctuons. Within the framework of this model, the glass transition and polyamorphous transformations are considered. It is shown that the frustration increases the temperature range in which the heterophase liquid state exists. the upper and lower boundaries of this temperature range are determined. These boundaries separate different phase states-amorphous solid, heterophase liquid, and fluid phases. Polyamorphous liquid-liquid transitions in the liquid are investigated. Frustration can call forth continuous fluid-solid phase transformation avoiding the first- or second-order phase transition. Conditions under which the first-order phase transition fraction takes place are formulated. Two scenarios of the first-order liquid-liquid polyamorphous transformation are described. As an example the glacial phase formation and the first-order liquid-liquid phase transition in triphenyl phosphate are considered and discussed. Impact of frustration on the liquid crystallization and crystallinity of the glassy state is studied.     

Liquid-liquid phase transitions in supercooled water studied by computer simulations of various water models

Liquid-liquid and liquid-vapor coexistence regions of various water models were determined by Monte Carlo (MC) simulations of isotherms of density fluctuation-restricted systems and by Gibbs ensemble MC simulations. All studied water models show multiple liquid-liquid phase transitions in the supercooled region: we observe two transitions of the TIP4P, TIP5P, and SPCE models and three transitions of the ST2 model. The location of these phase transitions with respect to the liquid-vapor coexistence curve and the glass temperature is highly sensitive to the water model and its implementation. We suggest that the apparent thermodynamic singularity of real liquid water in the supercooled region at about 228 K is caused by an approach to the spinodal of the first (lowest density) liquid-liquid phase transition. The well-known density maximum of liquid water at 277 K is related to the second liquid-liquid phase transition, which is located at positive pressures with a critical point close to the maximum. A possible order parameter and the universality class of liquid-liquid phase transitions in one-component fluids are discussed.     

Direct Raman evidence for a weak continuous phase transition in liquid water

This paper presents the Raman depolarization ratio of degassed ultrapure water as a function of temperature, in the range 303.4-314.4 K (30.2-41.2 degrees C). The pressure of the sample was the vapor pressure of water at the measurement temperature. The data provide a direct indication of the existence of a phase transition in the liquid at 307.7 K, 5.8 kPa (34.6 degrees C, 0.057 atm). The minimum in the heat capacity, C(p)(), of water occurs at 34.5 degrees C, 1.0 atm (J. Res. Natl. Bur. Stand. 1939, 23, 197(1)). The minimum in C(p)() is shallow, and the transition is a weak-continuous phase transition. The pressure coefficient of the viscosity of water changes sign as pressure increases for temperatures below 35 degrees C (Nature 1965, 207, 620(2)). The viscosity minimum tracks the liquid phase transition in the P, T plane where it connects with the minimum in the freezing point of pure water in the same plane (Proc. Am. Acad. Arts Sci. 1911-12, 47, 441(3)). Previously we argued (J. Chem. Phys. 1998, 109, 7379(4)) that the minimum in the pressure coefficient of viscosity signaled the elimination of three-dimensional connectivity in liquid water. These observations coupled with recent measurements of the coordination shell of water near 300 K (Science 2004, 304, 995(5)) suggest that the structural component that changes during the phase transition is tetrahedrally coordinated water. At temperatures above the transition, there is no tetrahedrally coordinated water in the liquid and locally planar water structures dominate the liquid structure. Water is a structured liquid with distinct local structures that vary with temperature. Furthermore, liquid water has a liquid-liquid phase transition near the middle of the normal liquid range.     

Liquid-solid transition in fully ionized hydrogen at ultra-high pressures

We study the phase diagram of an effective ion model of fully ionized hydrogen at ultra-high pressure. We assume that the protons interact with a screened Coulomb potential derived from a static linear response theory. This model accurately reproduces the physical properties of hydrogen for densities greater than g/ρ(m)=10 cm(3) corresponding to the range of the coupling parameter r(s) ? 0.6. The pressure range, P ? 20 TPa, is well beyond present experimental limitations. Assuming classical protons, we find that the zero temperature enthalpy of the perfect bcc crystal is slightly lower than for other structures at g/ρ(m)=12.47 cm(3) while the fcc structure gains stability at higher density. Using Monte Carlo calculations, we compute the free energy of various phases and locate the melting transition versus density. We find that on melting, bcc is energetically favored with respect to fcc over the entire range investigated. In the solid phase the system undergoes a structural transition from bcc at higher temperature to fcc at lower temperature. The free energy difference between these two structures is very small so that obtaining a quantitative estimate of this second transition line requires accuracy beyond that provided by our method. We estimate the effect of proton zero point motion on the bcc melting line for hydrogen, deuterium, and tritium by a path integral Monte Carlo method. Although zero point effects on hydrogen are large, since the two competing phases (bcc and liquid) have locally similar environments, the effect on the melting line is small; the melting temperature for hydrogen is lowered by about 10% with respect to the classical value.     

Temperature Effect on Liquid-Liquid Phase Equilibria of Aqueous Solutions of Fluorobenzene in Ethanol at the Pressure 101.2 kPa

Solubilities and liquid-liquid equilibrium (LLE) data for (water + ethanol + fluorobenzene) solutions are determined at T = (288.15, 298.15, 308.15) K and atmospheric pressure. The titration method is used to construct the binodal curves. All measured solubilities and liquid-liquid equilibrium data are reproduced quantitatively by empirical equations and the NRTL activity coefficient model. The correlated data are consistent to better than 0.3% in phase composition. The studied temperatures have a small effect on the phase equilibrium properties of the investigated systems. The influence of water concentration on the distribution coefficient of ethanol in the coexistent liquid phases and on the fluorobenzene selectivity is important. Practically, fluorobenzene may well be used to produce absolute alcohol.     

Roles of deeply supercooled ethanol in crystallization and solvation of LiI

The mechanisms of glass-liquid transition and crystallization of amorphous solid ethanol were investigated through detailed analyses of the interaction with LiI using time-of-flight secondary ion mass spectrometry and reflection absorption infrared spectroscopy. The LiI species adsorbed on the surface are incorporated into the bulk of ethanol at temperatures higher than 100 K, concomitantly with the reorganization of the ethanol molecules at the surface. This behavior is explicable by self-diffusion of the ethanol molecules as a consequence of the glass-liquid transition. The resulting liquid is a distinct phase, as revealed from the similarity of the IR absorption band to that of amorphous solid ethanol rather than liquid ethanol. The liquid-liquid phase transition occurs at 130 K, and a supercooled liquid ethanol is formed, as evidenced by formation of a metastable LiI solution when ethanol is deposited on the LiI film. The supercooled liquid ethanol is unstable, so that it crystallizes immediately at 130 K on the Ni(111) substrate. The film morphology changes continuously, even after crystallization, and the film abruptly becomes smoother before film evaporation. This behavior implies that crystallization is not completed and that a liquidlike phase coexists.