Do you want to publish a course? Click here

Two-State Thermodynamics and the Possibility of a Liquid-Liquid Phase Transition in Supercooled TIP4P/2005 Water

99   0   0.0 ( 0 )
 Added by John Biddle
 Publication date 2016
  fields Physics
and research's language is English




Ask ChatGPT about the research

Water shows intriguing thermodynamic and dynamic anomalies in the supercooled liquid state. One possible explanation of the origin of these anomalies lies in the existence of a metastable liquid-liquid phase transition (LLPT) between two (high and low density) forms of water. While the anomalies are observed in experiments on bulk and confined water and by computer simulation studies of water-like models, the existence of a LLPT in water is still debated. Unambiguous experimental proof of the existence of a LLPT in bulk supercooled water is hampered by fast ice nucleation which is a precursor of the hypothesized LLPT. Moreover, the hypothesized LLPT, being metastable, in principle cannot exist in the thermodynamic limit (infinite size, infinite time). Therefore, computer simulations of water models are crucial for exploring the possibility of the metastable LLPT and the nature of the anomalies. In this work, we present new simulation results in the NVT ensemble for one of the most accurate classical molecular models of water, TIP4P/2005. To describe the computed properties and explore the possibility of a LLPT we have applied two-structure thermodynamics, viewing water as a non-ideal mixture of two interconvertible local structures (states). The results suggest the presence of a liquid-liquid critical point and a LLPT in this model for the simulated length and time scales. We have compared the behavior of TIP4P/2005 with other popular water-like models, namely mW and ST2, and with real water, all of which are well described by two-state thermodynamics. In view of the current debate involving different studies of TIP4P/2005, we discuss consequences of metastability and finite size in observing the liquid-liquid separation. We also address the relationship between the phenomenological order parameter of two-structure thermodynamics and the microscopic nature of the low-density structure.



rate research

Read More

We present a review on the study of metastable silicon, primarily focusing mainly on the aspects of liquid-liquid transition, critical point and phase behaviour, structural and dynamic properties of liquid phase as well as crystal nucleation. We begin with an extensive survey of the investigations of liquid silicon pursued over three decades, with salient experimental, theoretical and simulation results. Following which we present various scenarios put forward to rationalize the density and related anomalies often observed in water and other network forming liquids. After which we present the more recent investigations (both simulation and experimental works) of the phase behavior of Silicon. Since a significant part of metastable silicon work is on a classical empirical potential an important question to address is the reliability of this potential in describing the behavior of silicon. To provide a critical assessment of the applicability of classical simulation results to real silicon we present a comparison of the structural, dynamical, and thermodynamic quantities obtained from the SW potential with those from ab initio simulations and with available experimental data. We also discuss the sensitivity of the thermodynamic properties to model parameters.
One of the most promising frameworks for understanding the anomalies of cold and especially supercooled water is that of two-structure thermodynamics, in which water is viewed as a non-ideal mixture of two interconvertible local structures. The non-ideality of this mixture may give rise, at very low temperatures, to a liquid-liquid phase transition (LLPT) and a liquid-liquid critical point (LLCP), at which thermodynamic response functions diverge. Vario
We propose a unified model combining the first-order liquid-liquid and the second-order ferroelectric phase transitions models and explaining various features of the $lambda$-point of liquid water within a single theoretical framework. It becomes clear within the proposed model that not only does the long-range dipole-dipole interaction of water molecules yield a large value of dielectric constant $epsilon$ at room temperatures, our analysis shows that the large dipole moment of the water molecules also leads to a ferroelectric phase transition at a temperature close to the lambda-point. Our more refined model suggests that the phase transition occurs only in the low density component of the liquid and is the origin of the singularity of the dielectric constant recently observed in experiments with supercooled liquid water at temperature T~233K. This combined model agrees well with nearly every available set of experiments and explains most of the well-known and even recently obtained results of MD simulations.
A novel liquid-liquid phase transition has been proposed and investigated in a wide variety of pure substances recently, including water, silica and silicon. From computer simulations using the Stillinger-Weber classical empirical potential, Sastry and Angell [1] demonstrated a first order liquid-liquid transition in supercooled silicon, subsequently supported by experimental and simulation studies. Here, we report evidence for a liquid-liquid critical end point at negative pressures, from computer simulations using the SW potential. Compressibilities exhibit a growing maximum upon lowering temperature below 1500 K and isotherms exhibit density discontinuities below 1120 K, at negative pressure. Below 1120 K, isotherms obtained from constant volume-temperature simulations exhibit non-monotonic, van der Waals-like behavior signaling a first order transition. We identify Tc ~ 1120 +/- 12 K, Pc -0.60 +/- 0.15 GPa as the critical temperature and pressure for the liquid-liquid critical point. The structure of the liquid changes dramatically upon decreasing the temperature and pressure. Diffusivities vary over 4 orders of magnitude, and exhibit anomalous pressure dependence near the critical point. A strong relationship between local geometry quantified by the coordination number, and diffusivity, is seen, suggesting that atomic mobility in both low and high density liquids can usefully be analyzed in terms of defects in the tetrahedral network structure. We have constructed the phase diagram of supercooled silicon. We identify the lines of compressibility, density extrema (maxima and minima) and the spinodal which reveal the interconnection between thermodynamic anomalies and the phase behaviour of the system as suggested in previous works [2-9]
We investigate the relaxation mechanism of a supercooled tetrahedral liquid at its limit of stability using isothermal isobaric ($NPT$) Monte Carlo (MC) simulations. In similarity with systems which are far from equilibrium but near the onset of jamming [OHern et.al., Phys. Rev. Lett. {bf 93}, 165702 (2004)], we find that the relaxation is characterized by two time-scales: the decay of long-wavelength (slow) fluctuations of potential energy is controlled by the the slope $[partial (G/N)/partial phi]$ of the Gibbs free energy ($G$) at a unique value of per particle potential energy $phi = phi_{mid}$. The short-wavelength (fast) fluctuations are controlled by the bath temperature $T$. The relaxation of the supercooled liquid is initiated with a dynamical crossover after which the potential energy fluctuations are biased towards values progressively lesser than $phi_{mid}$. The dynamical crossover leads to the change of time-scale, i.e., the decay of long-wavelength potential energy fluctuations (intermediate relaxation). Because of the condition [$partial^2 (G/N)/partial phi^2 = 0$] at $phi = phi_{mid}$, the slope $[partial (G/N)/partial phi]$ has a unique value and governs the intermediate stage of relaxation, which ends just after the crossover. In the subsequent stage, there is a relatively rapid crystallization due to lack of long-wavelength fluctuations and the instability at $phi_{mid}$, i.e., the condition that $G$ decreases as configurations with potential energies lower than $phi_{mid}$ are accessed. The dynamical crossover point and the associated change in the time-scale of fluctuations is found to be consistent with the previous studies.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا