Structure of molecular liquids: hard rod-disk mixtures

Bulk phase behavior of binary hard platelet mixtures from density functional theory

We investigate isotropic-isotropic, isotropic-nematic, and nematic-nematic phase coexistence in binary mixtures of circular platelets with vanishing thickness, continuous rotational degrees of freedom, and radial size ratios lambda up to 5. A fundamental measure density functional theory, previously used for the one-component model, is presented and results are compared against those from Onsager theory as a benchmark. For lambda<or=1.7 the system displays isotropic-nematic phase coexistence with a widening of the biphasic region for increasing values of lambda . For size ratios lambda>or=2, we find that demixing into two nematic states becomes stable and an isotropic-nematic-nematic triple point can occur. Fundamental measure theory gives a smaller isotropic-nematic biphasic region than Onsager theory and locates the transition at lower densities. Furthermore, nematic-nematic demixing occurs over a larger range of compositions at a given value of lambda than found in Onsager theory. Both theories predict the same topologies of the phase diagrams. The partial nematic order parameters vary strongly with composition and indicate that the larger particles are more strongly ordered than the smaller particles.

Structure and stability of isotropic states of hard platelet fluids

We study the thermodynamics and the pair structure of hard, infinitely thin, circular platelets in the isotropic phase. Monte Carlo simulation results indicate a rich spatial structure of the spherical expansion components of the direct correlation function, including nonmonotonical variation of some of the components with density. Integral equation theory is shown to reproduce the main features observed in simulations. The hypernetted chain closure, as well as its extended versions that include the bridge function up to second and third order in density, perform better than both the Percus-Yevick closure and Verlet bridge function approximation. Using a recent fundamental measure density functional theory, an analytic expression for the direct correlation function is obtained as the sum of the Mayer bond and a term proportional to the density and the intersection length of two platelets. This is shown to give a reasonable estimate of the structure found in simulations, but to fail to capture the nonmonotonic variation with density. We also carry out a density functional stability analysis of the isotropic phase with respect to nematic ordering and show that the limiting density is consistent with that where the Kerr coefficient vanishes. As a reference system, we compare to simulation results for hard oblate spheroids with small, but nonzero elongations, demonstrating that the case of vanishingly thin platelets is approached smoothly.