Density-functional, Landau, and Onsager theories of the isotropic-nematic transition of hard ellipsoids

On Using the BMCSL Equation of State to Renormalize the Onsager Theory Approach to Modeling Hard Prolate Spheroidal Liquid Crystal Mixtures

Modifications to the traditional Onsager theory for modeling isotropic-nematic phase transitions in hard prolate spheroidal systems are presented. Pure component systems are used to identify the need to update the Lee-Parsons resummation term. The Lee-Parsons resummation term uses the Carnahan-Starling equation of state to approximate higher-order virial coefficients beyond the second virial coefficient employed in Onsager's original theoretical approach. As more exact ways of calculating the excluded volume of two hard prolate spheroids of a given orientation are used, the division of the excluded volume by eight, which is an empirical correction used in the original Lee-Parsons resummation term, must be replaced by six to yield a better match between the theoretical and simulation results. These modifications are also extended to binary mixtures of hard prolate spheroids using the Boublík-Mansoori-Carnahan-Starling-Leland (BMCSL) equation of state.

Emergence of biaxial nematic phases in solutions of semiflexible dimers

We investigate the isotropic, uniaxial nematic and biaxial nematic phases, and the transitions between them, for a model lyotropic mixture of flexible molecules consisting of two rigid rods connected by a spacer with variable bending stiffness. We apply density-functional theory within the Onsager approximation to describe strictly excluded-volume interactions in this athermal model and to self-consistently find the orientational order parameters dictated by its complex symmetry, as functions of the density. Earlier work on lyotropic ordering of rigid bent-rod molecules is reproduced and extended to show explicitly the continuous phase transition at the Landau point, at a critical bend angle of 36^{∘}. For flexible dimers with no intrinsic biaxiality, we find that a biaxial nematic phase can nevertheless form at a sufficiently high density and low bending stiffness. For bending stiffness κ>0.86k_{B}T, this biaxial phase manifests as dimer bending fluctuations occurring preferentially in one plane. When the dimers are more flexible, κ<0.86k_{B}T, the modal shape of the fluctuating dimer is a V with an acute opening angle, and one of the biaxial order parameters changes sign, indicating a rotation of the directors. These two regions are separated by a narrow strip of uniaxial nematic in the phase diagram, which we generate in terms of the spacer stiffness and particle density.

Density functional theory for dense nematic liquid crystals with steric interactions

The celebrated work of Onsager on hard particle systems, based on the truncated second order virial expansion, is valid at relatively low volume fractions for large aspect ratio particles. While it predicts the isotropic-nematic phase transition, it does not provide a realistic equation of state in that the pressure remains finite for arbitrarily high densities. In this work, we derive a mean field density functional form of the Helmholtz free energy for nematics with hard core repulsion. In addition to predicting the isotropic-nematic transition, the model provides a more realistic equation of state. The energy landscape is much richer, and the orientational probability distribution function in the nematic phase possesses a unique feature-it vanishes on a nonzero measure set in orientation space.

Explicit excluded volume of cylindrically symmetric convex bodies

We represent explicitly the excluded volume V(e){B(1),B(2)} of two generic cylindrically symmetric, convex rigid bodies, B(1) and B(2), in terms of a family of shape functionals evaluated separately on B(1) and B(2). We show that V(e){B(1),B(2)} fails systematically to feature a dipolar component, thus making illusory the assignment of any shape dipole to a tapered body in this class. The method proposed here is applied to cones and validated by a shape-reconstruction algorithm. It is further applied to spheroids (ellipsoids of revolution), for which it shows how some analytic estimates already regarded as classics should indeed be emended.

Hard-body models of bulk liquid crystals

Hard models for particle interactions have played a crucial role in the understanding of the structure of condensed matter. In particular, they help to explain the formation of oriented phases in liquids made of anisotropic molecules or colloidal particles and continue to be of great interest in the formulation of theories for liquids in bulk, near interfaces and in biophysical environments. Hard models of anisotropic particles give rise to complex phase diagrams, including uniaxial and biaxial nematic phases, discotic phases and spatially ordered phases such as smectic, columnar or crystal. Also, their mixtures exhibit additional interesting behaviours where demixing competes with orientational order. Here we review the different models of hard particles used in the theory of bulk anisotropic liquids, leaving aside interfacial properties and discuss the associated theoretical approaches and computer simulations, focusing on applications in equilibrium situations. The latter include one-component bulk fluids, mixtures and polydisperse fluids, both in two and three dimensions, and emphasis is put on liquid-crystal phase transitions and complex phase behaviour in general.

Spinodal decomposition of a binary magnetic fluid confined to a surface

In our previous work [J. Chem. Phys. 136, 024502 (2012)], we reported a demixing phase transition of a quasi-two-dimensional, binary Heisenberg fluid mixture driven by the ferromagnetic interactions of the magnetic species. Here, we present a theoretical study for the time-dependent coarsening occurring within the two-phase region in the density-concentration plane, also known as spinodal decomposition. Our investigations are based on dynamical density functional theory (DDFT). The particles in the mixture are modeled as Gaussian soft spheres on a two-dimensional surface, where one component carries a classical spin of Heisenberg type. To investigate the two-phase region, we first present a linear stability analysis with respect to small, harmonic density perturbations. Second, to capture nonlinear effects, we calculate time-dependent structure factors by combining DDFT with Percus' test particle method. For the growth of the average domain size l during spinodal decomposition with time t, we observe a power-law behavior l∝t^{δ_{α}} with δ_{m}≃0.333 for the magnetic species and δ_{n}≃0.323 for the nonmagnetic species.

Isotropic and nematic phases of the hard spheroid fluid: a virial approach

We use a high level virial expansion to investigate the properties of the isotropic and nematic phases of the hard spheroid fluid. We use the Monte Carlo techniques described previously to calculate the virial coefficients up to seventh order and we represent the dependence of these coefficients on particle orientations via a spherical harmonic expansion. The expansion coefficients are determined using Lebedev quadrature which carries out the angular integration required exactly. For fairly spherical spheroids (1/3<a/b<3, where a and b are the semiaxes of the spheroid) the virial expansion appears to be convergent in both the isotropic and nematic phases up to high density. For more aspherical particles, the predictions varied in a more erratic fashion upon increasing the number of virials. The use of a finite angular basis set to represent the orientationally dependent virial coefficients could be a contributory cause to this.

Molecular-field-theory approach to the Landau theory of liquid crystals: uniaxial and biaxial nematics

Nematic liquid crystal phase diagrams in temperature-biaxiality space are usually complex. We construct a Landau theory based on the analogous molecular-field theory for orthorhombic biaxial nematic fluids. A formal procedure yields coefficients (some of which, unusually, can be tensorial) in this Landau expansion, correctly predicts the complete set of invariants formed from the ordering tensors, and avoids ad hoc parametrization of the molecular biaxiality. By regularizing the Landau expansion to avoid unwanted order parameter divergences at low temperatures, we predict phase behavior over the whole range of biaxiality. The resulting phase diagrams have the same topology as those of molecular-field theory.

Fluids of hard natural and Gaussian ellipsoids: A comparative study by integral equation theories

The hard Gaussian overlap (HGO) model for ellipsoids is compared to the hard ellipsoid of revolution (HER) model, in the isotropic fluid phase and within the framework of the Percus-Yevick (PY) and hypernetted chain (HNC) integral equation theories. The former model is often used in place of the latter in many approximate theories. Since the HGO model slightly overestimates the contact distance when the two ellipsoids are perpendicular to each other, it leads to small differences in the Mayer function of the two models, but nearly none in the integrals of these functions and particularly for the second virial coefficients. However, it leads to notable differences in the pair correlation functions, as obtained by the Percus-Yevick and the hypernetted chain theories, especially at high densities. The prediction of the stability of the isotropic phase with respect to orientational order, at high densities, is notably influenced by these small differences. Both theories predict that, for same aspect ratios, the HGO model overestimates the ordering, when compared to the HER model. This explains why the PY approximation predicts ordering for the HGO model with aspect ratio of 1:3, while it does not for the HER model, in accordance with the very first integral equation results obtained for this system, and at variance with many opposite claims from subsequent publications that used the HGO model in place of the HER model.

Direct correlation functions of binary mixtures of hard Gaussian overlap molecules

We study the direct correlation function (DCF) of a classical fluid mixture of nonspherical molecules. The components of the mixture are two types of hard ellipsoidal molecules with different elongations, interacting through the hard Gaussian overlap (HGO) model. Two different approaches are used to calculate the DCFs of this fluid, and the results are compared. Here, the Pynn approximation [J. Chem. Phys. 60, 4579 (1974)] is extended to calculate the DCF of the binary mixtures of HGO molecules, then we use a formalism based on the weighted density functional theory introduced by Chamoux and Perera [J. Chem. Phys. 104, 1493 (1996)]. These results are fairly in agreement with each other. The pressure of this system is also calculated using the Fourier zero components of the DCF. The results are in agreement with the Monte Carlo molecular simulation.

Anomalous viscoelasticity near the isotropic-nematic phase transition in liquid crystals

Recent optical Kerr effect experiments have shown that orientational relaxation of nematogens shows a pronounced slow down of the response function at intermediate times and also a power law decay near the isotropic-nematic (I-N) transition. In many aspects, this behavior appears to be rather similar to the ones observed in the supercooled liquid near-glass transition. We have performed molecular dynamics simulations of model nematogens (Gay-Berne with aspect ratio 3) to explore the viscoelasticity near the I-N transition and also investigated the correlation of viscoelasticity (if any) with orientational relaxation. It is found that although the viscosity indeed undergoes a somewhat sharper than normal change near the I-N transition, it is not characterized by any divergence-like behavior (like the ones observed in the supercooled liquid). The rotational friction, on the other hand, shows a much sharper rise as the I-N transition is approached. Interestingly, the probability distribution of the amplitude of the three components of the stress tensor shows anisotropy near the I-N transition-similar anisotropy has also been seen in the deeply supercooled liquid. Frequency dependence of viscosity shows several unusual behaviors: (a) There is a weak, power law dependence on frequency [eta(')(omega) approximately omega(-alpha)] at low frequencies and (b) there is a rapid increase in the sharp peak observed in eta(')(omega) in the intermediate frequency on approach to the I-N transition density. These features can be explained from the stress-stress time correlation function. The angular velocity correlation function also exhibits a power law decay in time. The reason for this is discussed.

In search of temporal power laws in the orientational relaxation near isotropic–nematic phase transition in model nematogens

Recent Kerr relaxation experiments by Gottke et al. have revealed the existence of a pronounced temporal power law decay in the orientational relaxation near the isotropic-nematic phase transition (INPT) of nematogens of rather small aspect ratio, kappa (kappa approximately 3-4). We have carried out very long (50 ns) molecular dynamics simulations of model (Gay-Berne) prolate ellipsoids with aspect ratio 3 in order to investigate the origin of this power law. The model chosen is known to undergo an isotropic to nematic phase transition for a range of density and temperature. The distance dependence of the calculated angular pair correlation function correctly shows the emergence of a long range correlation as the INPT is approached along the density axis. In the vicinity of INPT, the single particle second rank orientational time correlation function exhibits power law decay, (t(-alpha)) with exponent alpha approximately 2/3. More importantly, we find the sudden appearance of a pronounced power-law decay in the collective part of the second rank orientational time correlation function at short times when the density is very close to the transition density. The power law has an exponent close to unity, that is, the correlation function decays almost linearly with time. At long times, the decay is exponential-like, as predicted by Landau-de Gennes mean field theory. Since Kerr relaxation experiments measure the time derivative of the collective second rank orientational pair correlation function, the simulations recover the near independence of the signal on time observed in experiments. In order to capture the microscopic essence of the dynamics of pseudonematic domains inside the isotropic phase, we introduce and calculate a dynamic orientational pair correlation function (DOPCF) obtained from the coefficients in the expansion of the distinct part of orientational van Hove time correlation function in terms of spherical harmonics. The DOPCF exhibits power law relaxation when the pair separation length is below certain critical length. The orientational relaxation of a local director, defined in terms of the sum of unit vectors of all the ellipsoidal molecules, is also found to show slow power law relaxation over a long time scale. These results have been interpreted in terms of a newly developed mode coupling theory of orientational dynamics near the INPT. In the present case, the difference between the single particle and the collective orientational relaxation is huge which can be explained by the frequency dependence of the memory kernel, calculated from the mode coupling theory. The relationship of this power law with the one observed in a supercooled liquid near its glass transition temperature is explored.