Phase behavior of shape-changing spheroids

We introduce a simple model for a biaxial nematic liquid crystal. This consists of hard spheroids that can switch shape between prolate (rodlike) and oblate (platelike) subject to an energy penalty Δε. The spheroids are approximated as hard Gaussian overlap particles and are treated at the level of Onsager's second-virial description. We use both bifurcation analysis and a numerical minimization of the free energy to show that, for additive particle shapes, (i) there is no stable biaxial phase even for Δε=0 (although there is a metastable biaxial phase in the same density range as the stable uniaxial phase) and (ii) the isotropic-to-nematic transition is into either one of two degenerate uniaxial phases, rod rich or plate rich. We confirm that even a small amount of shape nonadditivity may stabilize the biaxial nematic phase.