No Arabic abstract
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.
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.
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.
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 consider a monomer-dimer system with a strong attractive dimer-dimer interaction that favors alignment. In 1979, Heilmann and Lieb conjectured that this model should exhibit a nematic liquid crystal phase, in which the dimers are mostly aligned, but do not manifest any translational order. We prove this conjecture for large dimer activity and strong interactions. The proof follows a Pirogov-Sinai scheme, in which we map the dimer model to a system of hard-core polymers whose partition function is computed using a convergent cluster expansion.
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.