Reexamining the phase diagram of the Gay-Berne fluid

This work reexamines and updates earlier investigations on the phase behaviour of the Gay-Berne liquid crystal model, concentrating on the effect of varying temperature. Constant volume and constant pressure Monte Carlo simulations are combined for systems consisting of N = 500 molecules along different isotherms over the reduced temperature range 0.60 less than or equal toT less than or equal to1.25. As in previous simulation studies of the model, the study identifies nematic and smectic B phases on compressing the isotropic fluid, the particular phase sequence depending on temperature. The nematic phase is found to be stable with respect to the isotropic phase for reduced temperatures T greater than or equal to0.75. In the temperature range 0.75 less than or equal toT less than or equal to1.25, the phase boundaries of the isotropic-nematic transition are obtained by computing the Helmholtz free energy of both phases from thermodynamic integration. From the simulation data, the relative volume change at the isotropic-nematic transition is about 2%, and this value appears to be rather insensitive to changes in temperature. On compressing the nematic phase, the Gay-Berne fluid undergoes a strong first-order transition to the smectic B phase. This transition is studied by using constant pressure simulation, and the coexistence properties are estimated from the limits of mechanical stability of the nematic phase. Larger relative volume changes are found at the transition than those suggested by previous studies, with typical values increasing up to 10.5% as the temperature is decreased. The results are consistent with the existence of strong coupling between nematic and smectic order parameters. For temperatures T less than or equal to0.70 the nematic phase is no longer stable, and the phase sequence isotropic-smectic B is observed. Therefore, the Gay-Berne model exhibits an isotropic-nematic-smectic B triple point. Extrapolating the present simulation data, this triple point is located approximately at reduced temperature T-INB approximate to 0.70 and reduced pressure P-INB approximate to 1.825.