Density-functional theory of inhomogeneous systems of hard spherocylinders

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.

Computational aspects of the smectization process in liquid crystals: An example study of a perfectly aligned two-dimensional hard-boomerang system

A replica method for calculation of smectic liquid crystal properties within the Onsager theory has been presented and applied to an exemplary case of two-dimensional perfectly aligned needlelike boomerangs. The method allows one to consider the complete influence of the interaction terms in contrast to the Fourier expansion method which uses mostly first or second order terms of expansion. The program based on the replica algorithm is able to calculate a single representative layer as an equivalent set of layers, depending on the size of the considered width of the sample integration interval. It predicts successfully smectic density distributions, energies, and layer thicknesses for different types of layer arrangement-of the antiferroelectric or of the smectic A order type. Specific features of the algorithm performance and influence of the numerical accuracy on the physical properties are presented. Future applications of the replica method to freely rotating molecules are discussed.

Fundamental measure theory for non-spherical hard particles: predicting liquid crystal properties from the particle shape

Density functional theory (DFT) for hard bodies provides a theoretical description of the effect of particle shape on inhomogeneous fluids. We present improvements of the DFT framework fundamental measure theory (FMT) for hard bodies and validate these improvements for hard spherocylinders. To keep the paper self-contained, we first discuss the recent advances in FMT for hard bodies that lead to the introduction of fundamental mixed measure theory (FMMT) in our previous paper (2015 Europhys. Lett. 109 26003). Subsequently, we provide an efficient semi-empirical alternative to FMMT and show that the phase diagram for spherocylinders is described with similar accuracy in both versions of the theory. Finally, we present a semi-empirical modification of FMMT whose predictions for the phase diagram for spherocylinders are in excellent quantitative agreement with computer simulation results.

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.

Neural-network approach to modeling liquid crystals in complex confinement

Finding the structure of a confined liquid crystal is a difficult task since both the density and order parameter profiles are nonuniform. Starting from a microscopic model and density-functional theory, one has to either (i) solve a nonlinear, integral Euler-Lagrange equation, or (ii) perform a direct multidimensional free energy minimization. The traditional implementations of both approaches are computationally expensive and plagued with convergence problems. Here, as an alternative, we introduce an unsupervised variant of the multilayer perceptron (MLP) artificial neural network for minimizing the free energy of a fluid of hard nonspherical particles confined between planar substrates of variable penetrability. We then test our algorithm by comparing its results for the structure (density-orientation profiles) and equilibrium free energy with those obtained by standard iterative solution of the Euler-Lagrange equations and with Monte Carlo simulation results. Very good agreement is found and the MLP method proves competitively fast, flexible, and refinable. Furthermore, it can be readily generalized to the richer experimental patterned-substrate geometries that are now experimentally realizable but very problematic to conventional theoretical treatments.

Simulation and theory of hybrid aligned liquid crystal films

We present a study of the effects of nanoconfinement on a system of hard Gaussian overlap particles interacting with planar substrates through the hard-needle-wall potential, extending earlier work by two of us [D. J. Cleaver and P. I. C. Teixeira, Chem. Phys. Lett. 338, 1 (2001)]. Here, we consider the case of hybrid films, where one of the substrates induces strongly homeotropic anchoring, while the other favors either weakly homeotropic or planar anchoring. These systems are investigated using both Monte Carlo simulation and density-functional theory, the latter implemented at the level of Onsager's second-virial approximation with Parsons-Lee rescaling. The orientational structure is found to change either continuously or discontinuously depending on substrate separation, in agreement with earlier predictions by others. The theory is seen to perform well in spite of its simplicity, predicting the positional and orientational structure seen in simulations even for small particle elongations.

Smectic, nematic, and isotropic phases in binary mixtures of thin and thick hard spherocylinders

A second-virial Onsager theory, based on Parsons-Lee rescaling and suitably extended to deal with multicomponent systems and smectic phases, has been used to calculate the phase diagram of a collection of binary mixtures of thin and thick hard spherocylinders. In particular, two types of phase diagrams are investigated. First, a number of binary mixtures where the two components have the same total length have been considered; in addition, the phase diagram of a binary mixture where the two components have the same volume has been calculated. For the particles of one of the two components, the length of the cylindrical part and the diameter have always been set equal to 5 and 1, respectively. Spherocylinders of the same total length and different diameter tend to demix considerably as soon as the diameter ratio deviates from unity. This happens especially at high pressures, when at least the phase richer in the thicker component is smectic. In the case where the two components have equal volumes, demixing is further increased due to the disparity not only in particle diameter but also in particle lengths. The incorporation of inhomogeneous layered phases is seen to alter significantly the phase diagrams calculated if only homogeneous phases are allowed, since transitions to a smectic phase often preempt those to a nematic or an isotropic phase. The apparent versatility of the recent experimental techniques suggests that the phase diagram features predicted by the theory might be also observed in real systems.

Dynamical density functional theory for anisotropic colloidal particles

We generalize the formalism of dynamical density functional theory for translational Brownian dynamics toward that of anisotropic colloidal particles which perform both translational and rotational Brownian motion. Using a mean-field approximation for the density functional and a Gaussian-segment model for the rod interaction, the dynamical density functional theory is then applied to a concentrated rod suspension in a confined slab geometry made by two parallel soft walls. The walls are either expanded or compressed and the relaxation behavior is investigated for an equilibrated starting configuration. We find distinctly different orientational ordering during expansion and compression. During expansion we observe preferential parallel ordering of the rods relative to the wall while during compression there is homeotropic ordering perpendicular to the wall. We find a nonexponential relaxation behavior in time. Furthermore, an external field which aligns the rods perpendicular to the walls is turned on or switched off and similar differences in the relaxational dynamics are found. Comparing the theoretical predictions to Brownian dynamics computer simulation data, we find good agreement.

Lennard-Jones sticks: A new model for linear molecules

We consider the anisotropic interaction between two line segments consisting of a homogeneous distribution of Lennard-Jones centers. The potential energy of such a pair cannot be expressed in closed form. However, we show that it may be approximated in a way that renders this intuitively appealing model competitive both for simulations and theory.

Microstructure and Self-Assembly of Inhomogeneous Rigid Rodlike Chains between Two Neutral Surfaces: A Hybrid Density Functional Approach

We use a hybrid density functional approach to investigate the microstructure and self-assembly of inhomogeneous rigid rodlike chains between two neutral surfaces, i.e., two hard walls. In the calculation, the rodlike molecule is modeled as a rigid rod linearly connected by the tangent sphere beads. The hybrid method combines a single-chain Monte Carlo (MC) simulation for the ideal-gas part of Helmholtz energy and a DFT approach for the excess Helmholtz energy. The DFT approach includes a modified fundamental measure theory for the excluded-volume effect, the first order thermodynamics perturbation theory for chain connectivity, and the mean field approximation for the van der Waals attraction. We investigate the effect of the chain length (i.e., aspect ratio) of the rodlike molecule and the separation between two surfaces on the microstructure and self-assembly of inhomogeneous rigid rodlike chains. For the athermal systems, the rodlike chain fluids present a smaller partitioning coefficient compared to the flexible chain fluids. For the thermal systems, lamellar thin films formed by the rigid rodlike molecules perpendicular to the neutral surface are observed. The effects of the head-head interaction and the separation on the self-assembly of the rodlike chain fluids in the slit are investigated.

Elasticity of a system with noncentral potentials

We derive expressions for determination of the stress and the elastic constants in systems composed of particles interacting via noncentral two-body potentials as thermal averages of products of first and second partial derivatives of the interparticle potentials and components of the interparticle separation vectors. These results are adapted to hard potentials, where the stress and the elastic constants are expressed as thermal averages of the components of normals to contact surfaces between the particles and components of vectors separating their centers. The averages require knowledge of the simultaneous contact probabilities of two pairs of particles. We apply the expressions to particles for which a contact function can be defined, and demonstrate the feasibility of the method by computing the stress and the elastic constants of a two-dimensional system of hard ellipses using Monte Carlo simulations.

Self-assembly of rod-coil block copolymers

We present a self-consistent field theory model for the self-assembly behavior of rod-coil block copolymers. The orientational interactions between the rods were modeled through a Maier-Saupe interaction, while the enthalpic interactions between rods and coils were modeled through a standard Flory-Huggins approach. We outline a "real-space" numerical approach to solve the self-consistent field equations for such rod-coil block copolymers. A major focus of our work is upon the nonlamellar phases observed in the experiments on such polymers. To develop a physical understanding of these phases and their regimes of occurrence, we compute the two-dimensional phase diagram for our model. The latter shows significant departures from the one-dimensional phase diagram, but matches qualitatively with the existing experimental results. We also present scaling arguments that rationalize the numerical results for the self-assembly behavior.

Capillary effects in a confined smectic phase of hard spherocylinders: influence of particle elongation

A system of hard rods confined into a pore with slit geometry (two parallel planar substrates) is studied theoretically in the regime of high packing fraction. In this regime the bulk system exhibits a nematic phase as well as a smectic-A (spatially layered) phase. When the system is confined, strong commensuration effects between the layer spacing and the pore width bring about a rich phenomenology, with a phase diagram showing layering and capillary transitions. The latter include capillary smectization transitions whereby a confined smectic phase occurs at conditions of saturation different from those of the corresponding bulk fluid. These transitions are seen to be intimately connected with layering transitions involving discontinuous changes in the number of layers inside the pore. This rich phenomenology is obtained by use of a sophisticated density-functional, Onsager-theory-based approach, especially suited to deal with strongly inhomogeneous fluids. The theory allows for a unified description of ordering and phase behavior of the fluid in confined geometry, and permits us to correlate the above behavior with the wetting properties of the fluid on a single substrate.

Parsons−Lee and Monte Carlo Study of Soft Repulsive Nematogens

A general approach based on the Parsons-Lee theory for soft repulsive molecular fluids is employed to investigate the nematogenic behavior of prolate thermotropic liquid crystals over a broad temperature range. The theory is solved for the particular case of the Kihara soft repulsive spherocylinder model, which is mapped into an effective hard core interaction with a temperature-dependent molecular diameter, expected to resemble the average size and shape of the soft molecules at a given temperature. The reduction of the effective molecular diameter with temperature in the Kihara soft repulsive fluid implies implicitly an increase of the elongation of the molecule and induces the stabilization of the nematic phase at smaller effective packing fractions, contrary to what is found for other fluid models. The rationalization of this effect in terms of excluded volume steric arguments is corroborated by the good general agreement between the Parsons-Lee approach and Monte Carlo simulations for the equation of state of the fluid in the vicinity of the isotropic-nematic transition.

Diffusion and viscosity of a calamitic liquid crystal model studied by computer simulation

We report a molecular dynamics simulation study on an ensemble of rod-like particles, each composed of nine soft spheres held rigidly along a line. We have calculated translational mean square displacements and velocity autocorrelation functions in the fluid phases exhibited by the model, i.e., smectic A, nematic and isotropic. These quantities have then been used to compute diffusion coefficients. In addition, we have calculated viscosities in the nematic and isotropic phases. Despite its crude nature, the model is capable of providing a faithful reproduction of many features of the transport behavior observed in real liquid-crystalline materials. The simulation results have been compared with the predictions of the modified affine transformation theory, finding only a fair agreement.

Smectic phase in a system of hard ellipsoids with isotropic attractive interactions

The smectic phase is studied for a thermotropic fluid model consisting of aligned hard ellipsoids with superimposed square-well attractive interactions of variable range. The system is analyzed using a density functional theory in which the hard-core contributions to the free-energy functional are treated within a nonlocal weighted density approximation and the attractive contributions are considered at a mean-field level. In the absence of attractions the model reduces, under appropriate scaling, to a fluid of hard spheres and therefore does not exhibit smectic ordering. It is shown that above a certain value of the square-well range, smectic ordering is stable relative to the nematic state at densities well inside the fluid region. The nematic-smectic-A transition is found to be continuous at high temperatures and first order at low temperatures, these two regimes being separated by a tricritical point at an intermediate temperature. These predictions have been confirmed by computer simulation of the model fluid. The results highlight that smectic ordering can be stabilized by coupling anisotropic short-range repulsions with the isotropic contribution of the soft attractive interactions. By increasing the pressure, the range of stability of the smectic phase is seen to decrease. At sufficiently high pressure, the smectic phase is suppressed, and the solid phase dominates. Our calculations show that smectic ordering is no longer stable if the range of the attractions is made too long ranged.

Smectic ordering of homogeneous semiflexible polymers

A self-consistent-field theory for fluids of homogeneous wormlike polymers exhibiting a one-dimensional spatial variation is presented. We have extended the treatment of excluded-volume effects by adding an effective interaction term which describes the excluded volume between wormlike cylindrical segments and terminal (or end) segments of the polymer molecules. This enables us to find a smectic-A phase in the case of homogeneous semiflexible polymers. Using this framework, we have investigated the occurrence of smectic-A, nematic, and isotropic phases in the second-virial (Onsager) approximation. Phase diagrams are calculated for systems characterized by different rigidities (i.e., persistence lengths). For the case of infinitely rigid molecules, the nematic-smectic transition appears to be mostly second order. Systems of semiflexible molecules exhibit mainly a first-order smectic-nematic transition, and their isotropic-nematic-smectic triple points are accessed for different rigidity values. The nematic-smectic transition line is in good agreement with previous analytical calculations, which were also performed assuming the second-virial approximation. However, the values of the volume fraction at the nematic-smectic transition are large compared with computer simulation results, indicating limitations of the second-virial approximation.

Numerical study of a calamitic liquid-crystal model: phase behavior and structure

We have studied an idealized calamitic liquid-crystal model, consisting of a linear rigid array of nine soft repulsive spheres, employing both theory and molecular dynamics simulation. The phase behavior (which includes crystalline, smectic, nematic, and isotropic phases) and structure of a collection of these rodlike particles have been determined by molecular dynamics simulation in an isothermal-isobaric ensemble. The liquid crystalline part of the phase diagram has been compared to that emerging from an Onsager-type density-functional theory. We have found a fair agreement between theory and computer simulation results, with a similar accuracy for the smectic to nematic and nematic to isotropic phase transitions.

Effect of particle geometry on phase transitions in two-dimensional liquid crystals

Using a version of density-functional theory which combines Onsager approximation and fundamental-measure theory for spatially nonuniform phases, we have studied the phase diagram of freely rotating hard rectangles and hard discorectangles. We find profound differences in the phase behavior of these models, which can be attributed to their different packing properties. Interestingly, bimodal orientational distribution functions are found in the nematic phase of hard rectangles, which cause a certain degree of biaxial order, albeit metastable with respect to spatially ordered phases. This feature is absent in discorectangles, which always show unimodal behavior. This result may be relevant in the light of recent experimental results which have confirmed the existence of biaxial phases. We expect that some perturbation of the particle shapes (either a certain degree of polydispersity or even bimodal dispersity in the aspect ratios) may actually destabilize spatially ordered phases thereby stabilizing the biaxial phase.

Capillary smectization and layering in a confined liquid crystal

Using density-functional theory, we have analyzed the phase behavior of a model liquid crystal confined between two parallel, planar surfaces (i.e., the so-called slit pore). As a result of confinement, a rich phase behavior arises. The complete liquid-crystal phase diagram of the confined fluid is mapped out as a function of wall separation and chemical potential. Strong commensuration effects in the film with respect to wall separation lead to enhanced smectic ordering, which gives capillary smectization (i.e., formation of a smectic phase in the pore), or frustrated smectic ordering, which suppresses capillary smectization. These effects also produce layering transitions. Our nonlocal density-functional-based analysis provides a unified picture of all the above phenomena.

Liquid–crystal phase diagrams of binary mixtures of hard spherocylinders

We have built the liquid crystal phase diagram of several binary mixtures of freely rotating hard spherocylinders employing a second-order virial density functional theory with Parsons scaling, suitably generalized to deal with mixtures and smectic phases. The components have the same diameter and aspect ratio of moderate value, typical of many mesogens. Attention has been paid to smectic-smectic demixing and the types of arrangement that rods can adopt in layered phases. Results are shown to depend on the aspect ratio of the individual components and on the ratio of their lengths. Smectic phases are seen not to easily mix together at sufficiently high pressures. Layered phases where the longer rods are the majority component have a smectic-A structure. In the opposite case, a smectic-A(2) phase is obtained where the shorter particles populate the layers and the longer ones prefer to stay parallel to the latter in the interlayer region.

Bulk inhomogeneous phases of anisotropic particles: a fundamental measure functional study of the restricted orientations model

The phase diagram of prolate and oblate particles in the restricted orientations approximation (Zwanzig model) is calculated. Transitions to different inhomogeneous phases (smectic, columnar, oriented, or plastic solid) are studied through minimization of the fundamental measure functional (FMF) of hard parallelepipeds. The study of parallel hard cubes (PHC's) as a particular case is also included motivated by recent simulations of this system. As a result a rich phase behavior is obtained which include, apart from the usual liquid crystal phases, a very peculiar phase (called here discotic smectic) which was already found in the only existing simulation of the model, and which turns out to be stable because of the restrictions imposed on the orientations. The phase diagram is compared at a qualitative level with simulation results of other anisotropic particle systems.

Effects of wetting and anchoring on capillary phenomena in a confined liquid crystal

A fluid of hard spherocylinders of length-to-breadth ratio L/D=5 confined between two identical planar, parallel walls--forming a pore of slit geometry--has been studied using a version of the Onsager density-functional theory. The walls impose an exclusion boundary condition over the particle's centers of mass, while at the same time favoring a particular anchoring at the walls, either parallel or perpendicular to the substrate. We observe the occurrence of a capillary transition, i.e., a phase transition associated with the formation of a nematic film inside the pore at a chemical potential different from micro(b)-the chemical potential at the bulk isotropic-nematic transition. This transition terminates at an Ising-type surface critical point. In line with previous studies based on the macroscopic Kelvin equation and the mesoscopic Landau-de Gennes approach, our microscopic model indicates that the capillary transition is greatly affected by the wetting and anchoring properties of the semi-infinite system, i.e., when the fluid is in contact with a single wall or, equivalently, the walls are at a very large distance. Specifically, in a situation where the walls are preferentially wetted by the nematic phase in the semi-infinite system, one has the standard scenario with the capillary transition taking place at chemical potentials less than micro(b) (capillary nematization transition or capillary ordering transition). By contrast, if the walls tend to orientationally disorder the fluid, the capillary transition may occur at chemical potentials larger than micro(b), in what may be called a capillary isotropization transition or capillary disordering transition. Moreover, the anchoring transition that occurs in the semi-infinite system may affect very decisively the confinement properties of the liquid crystal and the capillary transitions may become considerably more complicated.

Wetting properties of a hard-spherocylinder fluid on a substrate

A density-functional theory is used to analyze the wetting properties of a fluid made up of hard spherocylinders of a length-to-breadth ratio L/D=5 on a model substrate. The substrate imposes an exclusion boundary condition over the molecular centres of mass, while at the same time favoring a definite molecular orientation, either parallel or perpendicular to the substrate, in a region next to the substrate. The wetting properties of this system are seen to depend on the strength with which the substrate orients the molecules: as the latter is increased, wetting by nematic phase is followed by a region of partial wetting which then leads to reentrant wetting by nematic. The two wetting transitions correspond to wetting films with nematic director perpendicular and parallel to the substrate, respectively. Also, in the region of partial wetting, an anchoring transition occurs in the substrate-nematic interface between two different director configurations (parallel and perpendicular to the substrate). Finally, a metastable wetting transition by isotropic is also obtained. This model considerably enriches the wetting phenomenology of the hard-spherocylinder fluid on substrates, of which only the pure hard wall, with no surface control parameter available, has been considered so far.

Density-functional study of the nematic-isotropic interface of hard spherocylinders

The Somoza-Tarazona density-functional theory is applied to the isotropic-nematic interface of hard spherocylinders with length (L)-to-diameter (D) ratios in the range L/D=5-20. Properties such as the density and orientational order-parameter profiles and the variation of interfacial tension with bulk nematic tilt angle agree qualitatively with results of previous studies at larger values of L/D using both computer simulation and the Onsager second-virial approximation. The minimum interfacial tension is obtained at a tilt angle of 90 degrees. For values of L/D approximately 5, it is found that the Onsager approximation predicts a spurious minimum in the interfacial tension at small tilt angles.