Hard-body models of bulk liquid crystals

Ordering properties of anisotropic hard bodies in one-dimensional channels

The phase behavior and structural properties of hard anisotropic particles (prisms and dumbbells) are examined in one-dimensional channels using the Parsons-Lee (PL) theory, and the transfer-matrix and neighbor-distribution methods. The particles are allowed to move freely along the channel, while their orientations are constrained such that one particle can occupy only two or three different lengths along the channel. In this confinement setting, hard prisms behave as an additive mixture, while hard dumbbells behave as a non-additive one. We prove that all methods provide exact results for the phase properties of hard prisms, while only the neighbor-distribution and transfer-matrix methods are exact for hard dumbbells. This shows that non-additive effects are incorrectly included into the PL theory, which is a successful theory of the isotropic-nematic phase transition of rod-like particles in higher dimensions. In the one-dimensional channel, the orientational ordering develops continuously with increasing density, i.e., the system is isotropic only at zero density, while it becomes perfectly ordered at the close-packing density. We show that there is no orientational correlation in the hard prism system, while the hard dumbbells are orientationally correlated with diverging correlation length at close packing. On the other hand, positional correlations are present for all the systems, the associated correlation length diverging at close packing.

Elastic response of colloidal smectic liquid crystals: Insights from microscopic theory

Elongated colloidal rods at sufficient packing conditions are known to form stable lamellar or smectic phases. Using a simplified volume-exclusion model, we propose a generic equation of state for hard-rod smectics that is robust against simulation results and is independent of the rod aspect ratio. We then extend our theory by exploring the elastic properties of a hard-rod smectic, including the layer compressibility (B) and bending modulus (K_{1}). By introducing weak backbone flexibility we are able to compare our predictions with experimental results on smectics of filamentous virus rods (fd) and find quantitative agreement between the smectic layer spacing, the out-of-plane fluctuation strength, as well as the smectic penetration length λ=sqrt[K_{1}/B]. We demonstrate that the layer bending modulus is dominated by director splay and depends sensitively on lamellar out-of-plane fluctuations that we account for on the single-rod level. We find that the ratio between the smectic penetration length and the lamellar spacing is about two orders of magnitude smaller than typical values reported for thermotropic smectics. We attribute this to the fact that colloidal smectics are considerably softer in terms of layer compression than their thermotropic counterparts while the cost of layer bending is of comparable magnitude.

Liquid-crystal composites of carbon nanotubes in a magnetic field: Bridging continuum theory and a molecular-statistical approach

We propose an approach combining the continuum theory and molecular-statistical approach for a suspension of carbon nanotubes based on a negative diamagnetic anisotropy liquid crystal. Using the continuum theory, we show that in the case of an infinite sample in suspension it is possible to observe peculiar magnetic Fréedericksz-like transitions between three nematic phases: planar, angular, and homeotropic with different mutual orientations of liquid-crystal and nanotube directors. The transition fields between these phases are found analytically as functions of material parameters of the continuum theory. To account for the effects associated with temperature changes, we propose a molecular-statistical approach that allows obtaining the equations of orientational state for the orientation angles of the main axes of the nematic order, i.e., the liquid-crystal and carbon-nanotube directors in a similar form as was obtained within the continuum theory. Thus, it is possible to relate the parameters of the continuum theory, such as the surface-energy density of a coupling between molecules and nanotubes, to the parameters of the molecular-statistical model and the order parameters of the liquid crystal and carbon nanotubes. This approach allows determining the temperature dependencies of the threshold fields of transitions between different nematic phases, which is impossible in the framework of the continuum theory. In the framework of the molecular-statistical approach we predict the existence of an additional direct transition between the planar and homeotropic nematic phases of the suspension, which cannot be described based on the continuum theory. As the main results, the magneto-orientational response of the liquid-crystal composite is studied and a possible biaxial orientational ordering of the nanotubes in the magnetic field is shown.

Evidence of T-type structures of hard square boards in capillary confinement

We employ Onsager's second virial density functional theory combined with the Parsons-Lee theory within the restricted orientation (Zwanzig) approximation to examine the phase structure of hard square boards of dimensions (L×D×D) uniaxially confined in narrow slabs. Depending on the wall-to-wall separation (H), we predict a number of distinctly different capillary nematic phases, including a monolayer uniaxial or biaxial planar nematic, homeotropic with a variable number of layers, and a T-type structure. We determine that the favored phase is homotropic, and we observe first-order transitions from the homeotropic structure with n layers to n+1 layers as well as from homeotropic surface anchoring to a monolayer planar or T-type structure involving both planar and homeotropic anchoring at the pore surface. By increasing the packing fraction, we further demonstrate a reentrant homeotropic-planar-homeotropic phase sequence in a particular range (i.e., H/D=1.1 and 0.25≤L/D<0.26). We find that the T-type structure is more stable when the pore is wide enough with respect to the planar phase. The enhanced stability of the mixed-anchoring T-structure is unique for square boards and becomes manifest at pore width exceeding L+D. More specifically, the biaxial T-type structure emerges directly from the homeotropic state without intervention of a planar layer structure as observed for other convex particle shapes.

Reduced-variance orientational distribution functions from torque sampling

We introduce a method to sample the orientational distribution function in computer simulations. The method is based on the exact torque balance equation for classical many-body systems of interacting anisotropic particles in equilibrium. Instead of the traditional counting of events, we reconstruct the orientational distribution function via an orientational integral of the torque acting on the particles. We test the torque sampling method in two- and three-dimensions, using both Langevin dynamics and overdamped Brownian dynamics, and with two interparticle interaction potentials. In all cases the torque sampling method produces profiles of the orientational distribution function with better accuracy than those obtained with the traditional counting method. The accuracy of the torque sampling method is independent of the bin size, and hence it is possible to resolve the orientational distribution function with arbitrarily small angular resolutions.

Spontaneous Twisting of Achiral Hard Rod Nematics

Since Onsager's seminal work, hard rods have been taken as a prototype of nematic liquid crystals, characterized by uniaxial order and a uniform director field as a ground state. Here, using Onsager theory to calculate the free energy in the presence of arbitrary deformations, we find that hard rod nematics have an intrinsic tendency to twist around their ordering axis (double twist), driven by a mechanism in which the orientational fluctuations of particles play a key role. The anisotropic hard core potential used here is arguably the simplest form of interaction able to originate spontaneous breaking of mirror symmetry in a 3D fluid. Our results are discussed in relation to the recent discovery of a double twisted ground state in cylindrically confined lyotropic chromonic liquid crystals.

Multiphase Coexistence in Binary Hard Colloidal Mixtures: Predictions from a Simple Algebraic Theory

A general theoretical framework is proposed to quantify the thermodynamic properties of multicomponent hard colloidal mixtures. This framework is used to predict the phase behavior of mixtures of rods with spheres and rods with plates taking into account (liquid) crystal phases of both components. We demonstrate a rich and complex range of phase behaviors featuring a large variety of different multiphase coexistence regions, including two five-phase coexistence regions for hard rod/sphere mixtures, and even a six-phase equilibrium for hard rod/plate dispersions. The various multiphase coexistences featured in a particular mixture are in line with a recently proposed generalized phase rule and can be tuned through subtle variations of the particle shape and size ratio. Our approach qualitatively accounts for certain multiphase equilibria observed in rod/plate mixtures of clay colloids and will be a useful guide in tuning the phase behavior of shape-disperse mixtures in general.

Anomalous phase behavior of quasi-one-dimensional attractive hard rods

We study a two-state model of attractive hard rods using the transfer matrix method, where the centers of the particles are confined to a straight line, but the orientations of the rods can be parallel or perpendicular to the confining line. The rods are modeled as hard rectangles with length L and width D and decorated with attractive sites at both ends of the rectangles. We find that the particles align parallel to the line and form long chains at low densities, while they turn out of the line and form a Tonks gas at high densities. With increasing the stickiness between the rods, the structural change between parallel and perpendicular states becomes stronger and the pressure vs density curve becomes almost a horizontal line at the transition pressure. We show that such a behavior is reminiscent of the first-order phase transition. This manifests in the validity of the lever rule of the phase transitions for very sticky cases.

Active microrheology of colloidal suspensions of hard cuboids

By performing dynamic Monte Carlo simulations, we investigate the microrheology of isotropic suspensions of hard-core colloidal cuboids. In particular, we infer the local viscoelastic behavior of these fluids by studying the dynamics of a probe spherical particle that is incorporated in the host phase and is dragged by an external force. This technique, known as active microrheology, allows one to characterize the microscopic response of soft materials upon application of a constant force, whose intensity spans here three orders of magnitude. By tuning the geometry of cuboids from oblate to prolate as well as the system density, we observe different responses that are quantified by measuring the effective friction perceived by the probe particle. The resulting friction coefficient exhibits a linear regime at forces that are much weaker and larger than the thermal forces, whereas a nonlinear, force-thinning regime is observed at intermediate force intensities.

Interplay of orientational order and roughness in simulated thin film growth of anisotropically interacting particles

Roughness and orientational order in thin films of anisotropic particles are investigated using kinetic Monte Carlo simulations on a cubic lattice. Anisotropic next-neighbor interactions between the lattice particles were chosen to mimic the effects of shape anisotropy in the interactions of disk- or rodlike molecules with van der Waals attractions. Increasing anisotropy leads first to a preferred orientation in the film (which is close to the corresponding equilibrium transition) while the qualitative mode of roughness evolution (known from isotropic systems) does not change. At strong anisotropies, an effective step-edge (Ehrlich-Schwoebel) barrier appears and a nonequilibrium roughening effect is found, accompanied by reordering in the film which can be interpreted as the nucleation and growth of domains of lying-down disks or rods. The information on order and roughness is combined into a diagram of dynamic growth modes.

Phase behavior of hard circular arcs

By using Monte Carlo numerical simulation, this work investigates the phase behavior of systems of hard infinitesimally thin circular arcs, from an aperture angle θ→0 to an aperture angle θ→2π, in the two-dimensional Euclidean space. Except in the isotropic phase at lower density and in the (quasi)nematic phase, in the other phases that form, including the isotropic phase at higher density, hard infinitesimally thin circular arcs autoassemble to form clusters. These clusters are either filamentous, for smaller values of θ, or roundish, for larger values of θ. Provided the density is sufficiently high, the filaments lengthen, merge, and straighten to finally produce a filamentary phase while the roundels compact and dispose themselves with their centers of mass at the sites of a triangular lattice to finally produce a cluster hexagonal phase.

Second-virial theory for shape-persistent living polymers templated by disks

Living polymers composed of noncovalently bonded building blocks with weak backbone flexibility may self-assemble into thermoresponsive lyotropic liquid crystals. We demonstrate that the reversible polymer assembly and phase behavior can be controlled by the addition of (nonadsorbing) rigid colloidal disks which act as an entropic reorienting "template" onto the supramolecular polymers. Using a particle-based second-virial theory that correlates the various entropies associated with the polymers and disks, we demonstrate that small fractions of discotic additives promote the formation of a polymer nematic phase. At larger disk concentrations, however, the phase is disrupted by collective disk alignment in favor of a discotic nematic fluid in which the polymers are dispersed antinematically. We show that the antinematic arrangement of the polymers generates a nonexponential molecular-weight distribution and stimulates the formation of oligomeric species. At sufficient concentrations the disks facilitate a liquid-liquid phase separation which can be brought into simultaneously coexistence with the two fractionated nematic phases, providing evidence for a four-fluid coexistence in reversible shape-dissimilar hard-core mixtures without cohesive interparticle forces. We stipulate the conditions under which such a phenomenon could be found in experiment.

Topology of Orientational Defects in Confined Smectic Liquid Crystals

We propose a general formalism to characterize orientational frustration of smectic liquid crystals in confinement by interpreting the emerging networks of grain boundaries as objects with a topological charge. In a formal idealization, this charge is distributed in pointlike units of quarter-integer magnitude, which we identify with tetratic disclinations located at the end points and nodes. This coexisting nematic and tetratic order is analyzed with the help of extensive Monte Carlo simulations for a broad range of two-dimensional confining geometries as well as colloidal experiments, showing how the observed defect networks can be universally reconstructed from simple building blocks. We further find that the curvature of the confining wall determines the anchoring behavior of grain boundaries, such that the number of nodes in the emerging networks and the location of their end points can be tuned by changing the number and smoothness of corners, respectively.

Multiphase Coexistences in Rod-Polymer Mixtures

Using recently derived analytical equations of state for hard rod dispersions, we predict the phase behavior of athermal rod-polymer mixtures with free volume theory. The rods are modeled as hard spherocylinders, while the nonadsorbing polymer chains are described as penetrable hard spheres. It is demonstrated that all of the different types of phase states that are stable for pure colloidal rod dispersions can coexist with any combination of these phases if polymers are added, depending on the concentrations, rod aspect ratio, and polymer-rod size ratio. This includes novel two-, three-, and four-phase coexistences and isostructural coexistences between dilute and concentrated phases of the same kind, even for the more ordered (liquid) crystal phases. This work provides insight into the conditions at which particular multiphase coexistences are expected for well-defined model colloidal rod-polymer mixtures. We provide a quantitative map detailing the various types of isostructural coexistences, which confirms an early qualitative hypothesis by Bolhuis et al. ( J. Chem. Phys. 107, 1997 1551).

From n-layer planar ordering to the monolayer homeotropic structure of confined hard rods: The effect of shape anisotropy and wall-to-wall separation

Using the Parsons-Lee theory, we examined the effect of shape anisotropy and the wall-to-wall separation (H) on the phase behavior of the hard parallelepiped rods with dimensions L, D, and D (L>D) in such narrow slitlike pores which only one homeotropic layer can form. The phase structures, including biaxiality, planar nematic layering transition as well as planar to homeotropic, were studied for some separations in the range 2.5D≤H≤10.0D for H-D≤L<H.

Orientational ordering in a fluid of hard kites: A density-functional-theory study

Using density-functional theory we theoretically study the orientational properties of uniform phases of hard kites-two isosceles triangles joined by their common base. Two approximations are used: scaled particle theory and a new approach that better approximates third virial coefficients of two-dimensional hard particles. By varying some of their geometrical parameters, kites can be transformed into squares, rhombuses, triangles, and also very elongated particles, even reaching the hard-needle limit. Thus, a fluid of hard kites, depending on the particle shape, can stabilize isotropic, nematic, tetratic, and triatic phases. Different phase diagrams are calculated, including those of rhombuses, and kites with two of their equal interior angles fixed to 90^{∘}, 60^{∘}, and 75^{∘}. Kites with one of their unequal angles fixed to 72^{∘}, which have been recently studied via Monte Carlo simulations, are also considered. We find that rhombuses and kites with two equal right angles and not too large anisometry stabilize the tetratic phase but the latter stabilize it to a much higher degree. By contrast, kites with two equal interior angles fixed to 60^{∘} stabilize the triatic phase to some extent, although it is very sensitive to changes in particle geometry. Kites with the two equal interior angles fixed to 75^{∘} have a phase diagram with both tetratic and triatic phases, but we show the nonexistence of a particle shape for which both phases are stable at different densities. Finally, the success of the new theory in the description of orientational order in kites is shown by comparing with Monte Carlo simulations for the case where one of the unequal angles is fixed to 72^{∘}. These particles also present a phase diagram with stable tetratic and triatic phases.

Dense packings of hard circular arcs

This work investigates dense packings of congruent hard infinitesimally thin circular arcs in the two-dimensional Euclidean space. It focuses on those denotable as major whose subtended angle θ∈(π,2π]. Differently than those denotable as minor whose subtended angle θ∈[0,π], it is impossible for two hard infinitesimally thin circular arcs with θ∈(π,2π] to arbitrarily closely approach once they are arranged in a configuration, e.g., on top of one another, replicable ad infinitum without introducing any overlap. This makes these hard concave particles, in spite of being infinitesimally thin, most densely pack with a finite number density. This raises the question as to what are these densest packings and what is the number density that they achieve. Supported by Monte Carlo numerical simulations, this work shows that one can analytically construct compact closed circular groups of hard major circular arcs in which a specific, θ-dependent, number of them (counter) clockwise intertwine. These compact closed circular groups then arrange on a triangular lattice. These analytically constructed densest-known packings are compared to corresponding results of Monte Carlo numerical simulations to assess whether they can spontaneously turn up.

Statistical Physics Approach to Liquid Crystals: Dynamics of Mobile Potts Model Leading to Smectic Phase, Phase Transition by Wang-Landau Method

We study the nature of the smectic-isotropic phase transition using a mobile 6-state Potts model. Each Potts state represents a molecular orientation. We show that with the choice of an appropriate microscopic Hamiltonian describing the interaction between individual molecules modeled by a mobile 6-state Potts spins, we observe the smectic phase dynamically formed when we cool the molecules from the isotropic phase to low temperatures (T). In order to elucidate the order of the transition and the low-T properties, we use the high-performance Wang-Landau flat energy-histogram technique. We show that the smectic phase goes to the liquid (isotropic) phase by melting/evaporating layer by layer starting from the film surface with increasing T. At a higher T, the whole remaining layers become orientationally disordered. The melting of each layer is characterized by a peak of the specific heat. Such a succession of partial transitions cannot be seen by the Metropolis algorithm. The successive layer meltings/evaporations at low T are found to have a first-order character by examining the energy histogram. These results are in agreement with experiments performed on some smectic liquid crystals.

Demixing and tetratic ordering in some binary mixtures of hard superellipses

We examine the fluid phase behavior of binary mixtures of hard superellipses using the scaled particle theory. The superellipse is a general two-dimensional convex object that can be tuned between the elliptical and rectangular shapes continuously at a given aspect ratio. We find that the shape of the particle affects strongly the stability of isotropic, nematic, and tetratic phases in the mixture even if the side lengths of both species are fixed. While the isotropic-isotropic demixing transition can be ruled out using the scaled particle theory, the first order isotropic-nematic and the nematic-nematic demixing transition can be stabilized with strong fractionation between the components. It is observed that the demixing tendency is strongest in small rectangle-large ellipse mixtures. Interestingly, it is possible to stabilize the tetratic order at lower densities in the mixture of hard squares and rectangles where the long rectangles form a nematic phase, while the squares stay in the tetratic order.

Algebraic equations of state for the liquid crystalline phase behavior of hard rods

Based on simplifications of previous numerical calculations [H. Graf and H. Löwen, Phys. Rev. E 59, 1932 (1999)1063-651X10.1103/PhysRevE.59.1932], we propose algebraic free energy expressions for the smectic-A liquid crystal phase and the crystal phases of hard spherocylinders. Quantitative agreement with simulations is found for the resulting equations of state. The free energy expressions can be used to straightforwardly compute the full phase behavior for all aspect ratios and to provide a suitable benchmark for exploring how attractive interrod interactions mediate the phase stability through perturbation approaches such as free-volume or van der Waals theory.

Orientational ordering and layering of hard plates in narrow slitlike pores

We examine the ordering behavior of hard platelike particles in a very narrow, slitlike pore using the Parsons-Lee density functional theory and the restricted orientation approximation. We observe that the plates are orientationally ordered and align perpendicularly (face-on) to the walls at low densities, a first-order layering transition occurs between uniaxial nematic structures having n and n+1 layers at intermediate densities, and even a phase transition between a monolayer with parallel (edge-on) orientational order and n layers with a perpendicular one can be detected at high densities. In addition to this, the edge-on monolayer is usually biaxial nematic, and a uniaxial-biaxial nematic phase transition can be also seen at very high densities.

Biaxial, Twist-bend, and Splay-bend Nematic Phases of Banana-shaped Particles Revealed by Lifting the "Smectic Blanket"

We perform an extensive computational study on the phase behavior of hard banana-shaped particles, and show that biaxial, twist-bend, and splay-bend nematic phases are metastable with respect to a smectic phase for a system of hard bent spherocylinders. However, if the smectic phase is destabilized-either by polydispersity in the particle length or by curvature in the particle shape-stable biaxial, twist-bend, and splay-bend nematic phases are obtained. This provides a unified and consistent picture on the subtle role of particle shape on the phase behavior of bent rods.

Principal component analysis of the excluded area of two-dimensional hard particles

The excluded area between a pair of two-dimensional hard particles with given relative orientation is the region in which one particle cannot be located due to the presence of the other particle. The magnitude of the excluded area as a function of the relative particle orientation plays a major role in the determination of the bulk phase behavior of hard particles. We use principal component analysis (PCA) to identify the different types of excluded areas corresponding to randomly generated two-dimensional hard particles modeled as non-self-intersecting polygons and star lines (line segments radiating from a common origin). Only three principal components are required to have an excellent representation of the value of the excluded area as a function of the relative particle orientation for sufficiently anisotropic particles. Independent of the particle shape, the minimum value of the excluded area is always achieved when the particles are antiparallel to each other. The property that affects the value of the excluded area most strongly is the elongation of the particle shape. PCA identifies four limiting cases of excluded areas with one to four global minima at equispaced relative orientations. We study selected particle shapes using Monte Carlo simulations.

Kinetic theory of pattern formation in mixtures of microtubules and molecular motors

In this study we formulate a theoretical approach, based on a Boltzmann-like kinetic equation, to describe pattern formation in two-dimensional mixtures of microtubular filaments and molecular motors. Following the previous work by Aranson and Tsimring [Phys. Rev. E 74, 031915 (2006)PLEEE81539-375510.1103/PhysRevE.74.031915] we model the motor-induced reorientation of microtubules as collision rules, and devise a semianalytical method to calculate the corresponding interaction integrals. This procedure yields an infinite hierarchy of kinetic equations that we terminate by employing a well-established closure strategy, developed in the pattern-formation community and based on a power-counting argument. We thus arrive at a closed set of coupled equations for slowly varying local density and orientation of the microtubules, and study its behavior by performing a linear stability analysis and direct numerical simulations. By comparing our method with the work of Aranson and Tsimring, we assess the validity of the assumptions required to derive their and our theories. We demonstrate that our approximation-free evaluation of the interaction integrals and our choice of a systematic closure strategy result in a rather different dynamical behavior than was previously reported. Based on our theory, we discuss the ensuing phase diagram and the patterns observed.

Hard Competition: Stabilizing the Elusive Biaxial Nematic Phase in Suspensions of Colloidal Particles with Extreme Lengths

We use computer simulations to study the existence and stability of a biaxial nematic N_{b} phase in systems of hard polyhedral cuboids, triangular prisms, and rhombic platelets, characterized by a long (L), medium (M), and short (S) particle axis. For all three shape families, we find stable N_{b} states provided the shape is not only close to the so-called dual shape with M=sqrt[LS] but also sufficiently anisotropic with L/S>9,11,14,23 for rhombi, (two types of) triangular prisms, and cuboids, respectively, corresponding to anisotropies not considered before. Surprisingly, a direct isotropic-N_{b} transition does not occur in these systems due to a destabilization of N_{b} by a smectic (for cuboids and prisms) or a columnar (for platelets) phase at small L/S or by an intervening uniaxial nematic phase at large L/S. Our results are confirmed by a density functional theory provided the third virial coefficient is included and a continuous rather than a discrete (Zwanzig) set of particle orientations is taken into account.

Structure and Interfacial Tension of a Hard-Rod Fluid in Planar Confinement

The structural properties and interfacial tension of a fluid of rodlike hard-spherocylinder particles in contact with hard structureless flat walls are studied by means of Monte Carlo simulation. The calculated surface tension between the rod fluid and the substrate is characterized by a nonmonotonic trend as a function of the bulk concentration (density) over the range of isotropic bulk concentrations. As suggested by earlier theoretical studies, a surface-ordering scenario is confirmed by our simulations: the local orientational order close to the wall changes from uniaxial to biaxial nematic when the bulk concentration reaches about 85% of the value at the onset of the isotropic-nematic phase transition. The surface ordering coincides with a wetting transition whereby the hard wall is wetted by a nematic film. Accurate values of the fluid-solid surface tension, the adsorption, and the average particle-wall contact distance are reported (over a broad range of densities into the dense nematic region for the first time), which can serve as a useful benchmark for future theoretical and experimental studies on confined rod fluids. The simulation data are supplemented with predictions from second-virial density functional theory, which are in good qualitative agreement with the simulation results.

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.

Isotropic-nematic transition for hard rods on a three-dimensional cubic lattice

Using grand-canonical Monte Carlo (GCMC) simulations, we investigate the isotropic-nematic phase transition for hard rods of size L×1×1 on a three-dimensional cubic lattice. We observe such a transition for L≥6. For L=6, the nematic state has a negative order parameter, reflecting the co-occurrence of two dominating orientations. For L≥7, the nematic state has a positive order parameter, corresponding to the dominance of one orientation. We investigate rod lengths up to L=25 and find evidence for a very weakly first-order isotropic-nematic transition, while we cannot completely rule out a second-order transition. It was not possible to detect a density jump at the transition, despite using large systems containing several 10^{5} particles. The probability density distributions P(Q) from the GCMC simulations near the transition are very broad, pointing to strong fluctuations. Our results complement earlier results on the demixing (pseudonematic) transition for an equivalent system in two dimensions, which is presumably of Ising type and occurs for L≥7. We compare our results to lattice fundamental measure theory (FMT) and find that FMT strongly overestimates nematic order and consequently predicts a strong first-order transition. The rod packing fraction of the nematic coexisting states, however, agree reasonably well between FMT and GCMC.

Isotropic-nematic phase transition of uniaxial variable softness prolate and oblate ellipsoids

Onsager's theory of the isotropic-nematic phase separation of rod shaped particles is generalized to include particle softness and attractions in the anisotropic interparticle force field. The procedure separates a scaled radial component from the angular integral part, the latter being treated in essentially the same way as in the original Onsager formulation. Building on previous treatments of more idealised hard-core particle models, this is a step toward representing more realistic rod-like systems and also allowing temperature (and in principle specific chemical factors) to be included at a coarse grained level in the theory. The focus of the study is on the coexisting concentrations and associated coexistence properties. Prolate and oblate ellipsoids are considered in both the small and very large aspect ratio limits. Approximations to the terms in the angular integrals derived assuming the very large (prolate) and very small (oblate) aspect ratios limits are compared with the formally exact treatment. The approximation for the second virial coefficient matches the exact solution for aspect ratios above about 20 for the prolate ellipsoids and less than ca. 0.05 for the oblate ellipsoids from the numerical evaluation of the angular integrals. The temperature dependence of the coexistence density could be used to help determine the interaction potential of two molecules. The method works at temperatures above a certain threshold temperature where the second virial coefficient is positive.

Colloidal liquid crystals in square confinement: isotropic, nematic and smectic phases

We report on the confinement of colloidal liquid crystals in three dimensional chambers with a square footprint. To this end we use colloidal silica rods and exploit their relatively large density difference with respect to the dispersing solvent to study isotropic, nematic and smectic phases confined into a single chamber. Combining laser scanning confocal microscopy and soft-lithography techniques enables us to characterize the configurations down to the single particle level. We will focus on the smectic phase and compare to recent theories and simulations.

Phase behaviour and correlations of parallel hard squares: from highly confined to bulk systems

We study a fluid of two-dimensional parallel hard squares in bulk and under confinement in channels, with the aim of evaluating the performance of fundamental-measure theory (FMT). To this purpose, we first analyse the phase behaviour of the bulk system using FMT and Percus-Yevick (PY) theory, and compare the results with molecular dynamics and Monte Carlo simulations. In a second step, we study the confined system and check the results against those obtained from the transfer matrix method and from our own Monte Carlo simulations. Squares are confined to channels with parallel walls at angles of 0° or 45° relative to the diagonals of the parallel hard squares, respectively, which allows for an assessment of the effect of the external-potential symmetry on the fluid structural properties. In general FMT overestimates bulk correlations, predicting the existence of a columnar phase (absent in simulations) prior to crystallization. The equation of state predicted by FMT compares well with simulations, although the PY approach with the virial route is better in some range of packing fractions. The FMT is highly accurate for the structure and correlations of the confined fluid due to the dimensional crossover property fulfilled by the theory. Both density profiles and equations of state of the confined system are accurately predicted by the theory. The highly non-uniform pair correlations inside the channel are also very well described by FMT.

Entropy-driven formation of chiral nematic phases by computer simulations

Predicting the macroscopic chiral behaviour of liquid crystals from the microscopic chirality of the particles is highly non-trivial, even when the chiral interactions are purely entropic in nature. Here we introduce a novel chiral hard-particle model, namely particles with a twisted polyhedral shape and obtain a stable fully entropy-driven cholesteric phase by computer simulations. By slightly modifying the triangular base of the particle, we are able to switch from a left-handed prolate (calamitic) to a right-handed oblate (discotic) cholesteric phase using the same right-handed twisted particle model. Furthermore, we show that not only prolate and oblate chiral nematic phases, but also other novel entropy-driven phases, namely chiral blue phases, chiral nematic phases featuring both twist and splay deformations, chiral biaxial nematic phases with one of the axes twisted, can be obtained by varying particle biaxiality and chirality. Our results allow to identify general guidelines for the stabilization of these phases.

Entropy-Driven Chiral Order in a System of Achiral Bent Particles

Why should achiral particles organize into a helical structure? Here, using theory and molecular dynamics simulations we show that at high concentration crescent-shaped particles interacting through a purely repulsive potential form the twist-bend nematic phase, which features helical order of the twofold symmetry axes of particles, with doubly degenerate handedness. Spontaneous breaking of the chiral symmetry is driven by the entropic gain that derives from the decrease in excluded volume in the helical arrangement. Crucial to this purpose is the concave shape of particles. This study is based on a general formulation of the Onsager theory, which includes biaxiality and polarity of phase and particles, in addition to the space modulation of order. Molecular dynamics simulations corroborate the theoretical predictions and provide further insights into the structure of the helical phase.

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.

Elasticity of nematic phases with fundamental measure theory

In a previous publication [R. Wittmann, M. Marechal, and K. Mecke, Europhys. Lett. 109, 26003 (2015)], we introduced fundamental mixed measure theory (FMMT) for mixtures of anisotropic hard bodies, which shows that earlier results with an empirical parameter are inaccurate. Now we provide a deeper insight into the background of this theory in integral geometry. We study the Frank elastic coefficients in the nematic phase of the hard spherocylinder fluid. The framework of FMMT provides us with the required direct correlation function without additional input of an equation of state. A series representation of the mixed measure gives rise to closed analytical formulas for the elastic constants that only depend on the density, order parameters, and the particle geometry, pointing out a significant advantage of our geometry-based approach compared to other density functionals. Our elastic coefficients are in good agreement with computer simulations and increase with the density and the nematic order parameter. We confirm earlier mean-field predictions in the limits of low orientational order and infinitely long rods.