A molecular dynamics study of a bi-disperse liquid crystal mixture using a generalized Gay–Berne potential

Fractionation in Gay-Berne liquid crystal mixtures

We present a constant-pressure molecular dynamics simulation study of the phase behavior of binary (50:50) Gay-Berne liquid crystal mixtures consisting of elongated particles with different lengths (LA>LB) and equal diameters. We focus on systems at dense liquid-state conditions. Considering three mixtures characterized by different values of LA(B) and different length ratios q=LB/LA<1, we find complex fluid-fluid phase behavior resulting from the interplay between nematic, smectic-A-type, or smectic-B-type orientational ordering, on the one hand, and demixing into two phases of different composition (fractionation), on the other hand. The driving "forces" of demixing transitions are the temperature and the length ratio. Indeed, in the system characterized by the largest value of q (q=0.86) orientational order occurs already in mixed states, whereas full fractionation is found at q=0.71. The two resulting states are either of type smectic-B-nematic (intermediate temperatures) or smectic-B-smectic-B (low temperatures). In the intermediate case q=0.80 we observe a stepwise ordering and demixing behavior on cooling the system from high temperatures. Moreover, our results show that the stability range of (partially) nematic structures in mixtures of sufficiently small q can be significantly larger than in the pure counterparts, in qualitative agreement with experimental observations.

Lennard-Jones sticks: A new model for linear molecules

We consider the anisotropic interaction between two line segments consisting of a homogeneous distribution of Lennard-Jones centers. The potential energy of such a pair cannot be expressed in closed form. However, we show that it may be approximated in a way that renders this intuitively appealing model competitive both for simulations and theory.

Two-dimensional chiral model for liquid crystals, bent hard needles: A Monte Carlo simulation

The liquid crystalline behavior of a two dimensional (2D) model of hard needles bent into a "zigzag shape" is studied. This model, originally designed to study two dimensional chiral segregation, also shows liquid crystalline behavior and has some anomalous features which are contrasted in relation to the following: (i) Most of the microscopical models used to study liquid crystals have a symmetry axis that coincides with a molecular axis; (ii) in three-dimensions, chiral molecules can form cholesteric instead of nematic phases; (iii) the smectic phase is usually found when attractions are present or at least when the molecules have finite volume. Despite the fact that the present 2D model does not have any of these characteristics, numerical evidence is found for the occurrence of nematic and smectic phases. Since these molecules are athermal, infinitely repulsive, and infinitesimally thin, the liquid crystalline characteristics are attributed to excluded volume effects. To determine the mesophases of the model, both nematic and smectic order parameters as well as distribution functions are computed.