No Arabic abstract
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.
Spontaneous onset of a low temperature topologically ordered phase in a 2-dimensional (2D) lattice model of uniaxial liquid crystal (LC) was debated extensively pointing to a suspected underlying mechanism affecting the RG flow near the topological fixed point. A recent MC study clarified that a prior crossover leads to a transition to nematic phase. The crossover was interpreted as due to the onset of a perturbing relevant scaling field originating from the extra spin degree of freedom. As a counter example and in support of this hypothesis, we now consider V-shaped bent-core molecules with rigid rod-like segments connected at an assigned angle. The two segments of the molecule interact with the segments of all the nearest neighbours on a square lattice, prescribed by a biquadratic interaction. We compute equilibrium averages of different observables with Monte Carlo techniques as a function of temperature and sample size. For the chosen molecular bend angle and symmetric inter-segment interaction between neighbouirng molecules, the 2D system shows two transitions as a function of T: the higher one at T1 leads to a topological ordering of defects associated with the major molecular axis without a crossover, imparting uniaxial symmetry to the medium described by the first fundamental group of the order parameter space $pi_{1}$= $Z_{2}$ (inversion symmetry). The second at T2 leads to a medium displaying biaxial symmetry with $pi_{1}$ = Q (quaternion group). The biaxial phase shows a self-similar microscopic structure with the three axes showing power law correlations with vanishing exponents as the temperature decreases.
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.
We consider an off-lattice liquid crystal pair potential in strictly two dimensions. The potential is purely repulsive and short-ranged. Nevertheless, by means of a single parameter in the potential, the system is shown to undergo a first-order phase transition. The transition is studied using mean-field density functional theory, and shown to be of the isotropic-to-nematic kind. In addition, the theory predicts a large density gap between the two coexisting phases. The first-order nature of the transition is confirmed using computer simulation and finite-size scaling. Also presented is an analysis of the interface between the coexisting domains, including estimates of the line tension, as well as an investigation of anchoring effects.
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.
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.