Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCluster Monte Carlo and numerical mean field analysis for the water liquid--liquid phase transition
145
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
By the Wolffs cluster Monte Carlo simulations and numerical minimization within a mean field approach, we study the low temperature phase diagram of water, adopting a cell model that reproduces the known properties of water in its fluid phases. Both methods allows us to study the water thermodynamic behavior at temperatures where other numerical approaches --both Monte Carlo and molecular dynamics-- are seriously hampered by the large increase of the correlation times. The cluster algorithm also allows us to emphasize that the liquid--liquid phase transition corresponds to the percolation transition of tetrahedrally ordered water molecules.
rate research
Read More
Liquid-liquid phase transition (LLPT) in supercooled water has been a long-standing controversial issue. We show simulation results of real stable first-order phase transitions between high and low density liquid (HDL and LDL)-like structures in confined supercooled water in both positive and negative pressures. These topological phase transitions originate from H-bond network ordering in molecular rotational mode after molecular exchanges are frozen. It is explained by the order parameter-dependent free energy change upon mixing liquid-like and ice-like moieties of H-bond orientations which is governed by their two- to many-body interactions. This unexplored purely H-bond orientation-driven topological phase gives mid-density and stable intermediate mixed-phase with high and low density structures. The phase diagram of supercooled water demonstrate the second and third critical points of water.
Investigations of the phase diagram of biaxial liquid crystal systems through analyses of general Hamiltonian models within the simplifications of mean-field theory (MFT), as well as by computer simulations based on microscopic models, are directed towards an appreciation of the role of the underlying molecular-level interactions to facilitate its spontaneous condensation into a nematic phase with biaxial symmetry. Continuing experimental challenges in realising such a system unambiguously, despite encouraging predictions from MFT for example, are requiring more versatile simulational methodologies capable of providing insights into possible hindering barriers within the system, typically gleaned through its free energy dependences on relevant observables as the system is driven through the transitions. The recent brief report from this group [B. Kamala Latha, et. al., Phys. Rev. E 89, 050501 (R), 2014] summarizing the outcome of detailed Monte Carlo simulations carried out employing entropic sampling technique, suggested a qualitative modification of the MFT phase diagram as the Hamiltonian is asymptotically driven towards the so-called partly-repulsive regions. It was argued that the degree of the (cross) coupling between the uniaxial and biaxial tensor components of neighbouring molecules plays a crucial role in facilitating, or otherwise, a ready condensation of the biaxial phase, suggesting that this could be a plausible f actor in explaining the experimental difficulties. In this paper, we elaborate this point further, providing additional evidences from curious variations of free-energy profiles with respect to the relevant orientational order parameters, at different temperatures bracketing the phase transitions.
Using event driven molecular dynamics simulations, we study a three dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point and a liquid-liquid critical point. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and give us insight on the mechanisms ruling the dependence of the two critical points on the potentials parameters. The soft core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.
We investigate the phase behavior of a single-component system in 3 dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature {bf 409}, 692 (2001)] that, even with no evidences of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas--low-density liquid (LDL) critical point, and the other in a gas--high-density liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the 3-parameter space of the soft-core potential and we perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.
Longitudinal and transverse sound velocities of Lennard-Jones systems are calculated at the liquid-solid coexistence using the additivity principle. The results are shown to agree well with the ``exact values obtained from their relations to excess energy and pressure. Some consequences, in particular, in the context of the Lindemanns melting rule and Stokes-Einstein relation between the self-diffusion and viscosity coefficients are discussed. Comparison with available experimental data on the sound velocities of solid argon at melting conditions is provided.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
