No Arabic abstract
Using grand-canonical Monte Carlo simulations, we investigate the phase diagram of hard rods of length $L$ with additional contact (sticky) attractions on square and cubic lattices. The phase diagram shows a competition between gas-liquid and ordering transitions (which are of demixing type on the square lattice for $L ge 7$ and of nematic type on the cubic lattice for $L ge 5$). On the square lattice, increasing attractions initially lead to a stabilization of the isotropic phase. On the cubic lattice, the nematic transition remains of weak first order upon increasing the attractions. In the vicinity of the gas-liquid transition, the coexistence gap of the nematic transition quickly widens. These features are different from nematic transitions in the continuum.
Using overdamped Brownian dynamics simulations we investigate the isotropic-nematic (IN) transition of self-propelled rods in three spatial dimensions. For two well-known model systems (Gay-Berne potential and hard spherocylinders) we find that turning on activity moves to higher densities the phase boundary separating an isotropic phase from a (nonpolar) nematic phase. This active IN phase boundary is distinct from the boundary between isotropic and polar-cluster states previously reported in two-dimensional simulation studies and, unlike the latter, is not sensitive to the system size. We thus identify a generic feature of anisotropic active particles in three dimensions.
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 report phase separation and liquid-crystal ordering induced by scalar activity in a system of Soft Repulsive Spherocylinders (SRS) of aspect ratio $L/D = 5 $. Activity was introduced by increasing the temperature of half of the SRS (labeled textit{`hot}) while maintaining the temperature of the other half constant at a lower value (labeled textit{`cold}). The difference between the two temperatures scaled by the lower temperature provides a measure of the activity. Starting from different equilibrium initial phases, we find that activity leads to segregation of the hot and cold particles. Activity also drives the cold particles through a phase transition to a more ordered state and the hot particles to a state of less order compared to the initial equilibrium state. The cold components of a homogeneous isotropic (I) structure acquire nematic (N) and, at higher activity, crystalline (K) order. Similarly, the cold zone of a nematic initial state undergoes smectic (Sm) and crystal ordering above a critical value of activity while the hot component turns isotropic. We find that the hot particles occupy a larger volume and exert an extra kinetic pressure, confining, compressing and provoking an ordering transition of the cold-particle domains.
Isotropic-Nematic and Nematic-Nematic transitions from a homogeneous base state of a suspension of high aspect ratio, rod-like magnetic particles are studied for both Maier-Saupe and the Onsager excluded volume potentials. A combination of classical linear stability and asymptotic analyses provides insight into possible nematic states emanating from both the isotropic and nematic non-polarized equilibrium states. Local analytical results close to critical points in conjunction with global numerical results (Bhandar, 2002) yields a unified picture of the bifurcation diagram and provides a convenient base state to study effects of external orienting fields.
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.