Monte Carlo simulations of chiral and achiral nematic droplets: thermal quenches and role of the elastic constants

We used Monte Carlo simulations with an elastically anisotropic pairwise potential to simulate the equilibrium structure of chiral and achiral nematic liquid crystals in spherical confinement under perpendicular and tangential strong anchoring. The stable structures are obtained by a quenching process, where the temperature is a parameter of the system. We calculate and map the distribution of elastic energy corresponding to splay, twist, and bend. Furthermore, we study the effect of changing the ratio of the elastic constants for achiral and low chirality droplets and found a diverse scenario of possible structures. The method described here allows the introduction of temperature and elastic anisotropy naturally within its framework. It is simpler to implement and faster to execute than the Landau-de Gennes (LdG) method used for simulating liquid crystal droplets.