Phase diagram of the uniaxial and biaxial soft–core Gay–Berne model

Classical molecular dynamics simulations have been used to explore the phase diagrams for a family of attractive–repulsive soft–core Gay–Berne models [R. Berardi, C. Zannoni, J. S. Lintuvuori, and M. R. Wilson, J. Chem. Phys. 131, 174107 (2009)] and determine the effect of particle softness, i.e., of a moderately repulsive short–range interaction, on the order parameters and phase behaviour of model systems of uniaxial and biaxial ellipsoidal particles. We have found that isotropic, uniaxial, and biaxial nematic and smectic phases are obtained for the model. Extensive calculations of the nematic region of the phase diagram show that endowing mesogenic particles with such soft repulsive interactions affect the stability range of the nematic phases, and in the case of phase biaxiality it also shifts it to lower temperatures. For colloidal particles, stabilised by surface functionalisation, (e.g., with polymer chains), we suggest that it should be possible to tune liquid crystal behaviour to increase the range of stability of uniaxial and biaxial phases (by varying solvent quality). We calculate second virial coefficients and show that they are a useful means of characterising the change in effective softness for such systems. For thermotropic liquid crystals, the introduction of softness in the interactions between mesogens with overall biaxial shape (e.g., through appropriate conformational flexibility) could provide a pathway for the actual chemical synthesis of stable room – temperature biaxial nematics.