Computer simulations of biaxial nematics

Simulating the Lyotropic Phase Behavior of a Partially Self-Complementary DNA Tetramer

DNA oligomers in solution have been found to develop liquid crystal phases via a hierarchical process that involves Watson-Crick base pairing, supramolecular assembly into columns of duplexes, and long-range ordering. The multiscale nature of this phenomenon makes it difficult to quantitatively describe and assess the importance of the various contributions, particularly for very short strands. We performed molecular dynamics simulations based on the coarse-grained oxDNA model, aiming to depict all of the assembly processes involved and the phase behavior of solutions of the DNA GCCG tetramers. We find good quantitative matching to experimental data at both levels of molecular association (thermal melting) and collective ordering (phase diagram). We characterize the isotropic state and the low-density nematic and high-density columnar liquid crystal phases in terms of molecular order, size of aggregates, and structure, together with their effects on diffusivity processes. We observe a cooperative aggregation mechanism in which the formation of dimers is less thermodynamically favored than the formation of longer aggregates.

Mean-field model for a mixture of biaxial nematogens and dipolar nanoparticles

We analyze a mean-field model for mixtures involving biaxial nematogens and dipolar nanoparticles, taking into account not only orientational and isotropic pair interactions between nematogens but also orientational nematogen-nanoparticle interactions. We determine bulk equilibrium phase diagrams for a wide range of interaction strengths, identifying in each case the effect of the nanoparticles on the stability of nematic phases and on the appearance of multicritical points. Special attention is given to the limit of a low concentration of nanoparticles, in which their effect on the temperatures of both the first-order uniaxial-isotropic and the continuous biaxial-uniaxial transitions is investigated in detail.

Phases of the hard-plate lattice gas on a three-dimensional cubic lattice

We study the phase diagram of a lattice gas of 2×2×1 hard plates on the three-dimensional cubic lattice. Each plate covers an elementary plaquette of the cubic lattice, with the constraint that a site can belong to utmost one plate. We focus on the isotropic system, with equal fugacities for the three orientations of plates. We show, using grand canonical Monte Carlo simulations, that the system undergoes two phase transitions when the density of plates is increased: the first from a disordered fluid phase to a layered phase, and the second from the layered phase to a sublattice-ordered phase. In the layered phase, the system breaks up into disjoint slabs of thickness two along one spontaneously chosen Cartesian direction, corresponding to a twofold (Z_{2}) symmetry breaking of translation symmetry along the layering direction. Plates with normals perpendicular to this layering direction are preferentially contained entirely within these slabs, while plates straddling two adjacent slabs have a lower density, thus breaking the symmetry between the three types of plates. We show that the slabs exhibit two-dimensional power-law columnar order even in the presence of a nonzero density of vacancies. In contrast, interslab correlations of the two-dimensional columnar order parameter decay exponentially with the separation between the slabs. In the sublattice-ordered phase, there is twofold symmetry breaking of lattice translation symmetry along all three Cartesian directions. We present numerical evidence that the disordered to layered transition is continuous and consistent with universality class of the three-dimensional O(3) model with cubic anisotropy, while the layered to sublattice transition is first-order in nature.

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.

External field induced defect transformation in circular confined Gay-Berne liquid crystals

Normally, defects in two-dimensional, circular, confined liquid crystals can be classified into four types based on the position of singularities formed by liquid crystal molecules, i.e., the singularities located inside the circle, at the boundary, outside the circle, and outside the circle at infinity. However, it is considered difficult for small aspect ratio liquid crystals to generate all these four types of defects. In this study, we use molecular dynamics simulation to investigate the defect formed in Gay-Berne, ellipsoidal liquid crystals, with small aspect ratios confined in a circular cavity. As expected, we only find two types of defects (inside the circle and at the boundary) in circular, confined, Gay-Berne ellipsoids under static conditions at various densities, aspect ratios, and interactions between the wall and liquid crystals. However, when introducing an external field to the system, four types of defects can be observed. With increasing the strength of the external field, the singularities in the circular, confined system change from the inside to the boundary and the outside, and the farthest position that the singularities can reach depends on the strength of the external field. We further introduce an alternating, triangular wave, external field to the system to check if we can observe the transformation of different defects within an oscillating period. We find that the position of the singularities greatly depends on the oscillating intensity and oscillating period. By changing the oscillating intensity and oscillating period of the external field, the defect types can be adjusted, and the transformation between different defects can be easily observed. This provides a feasible way to modulate liquid crystal defects and investigate the transformation between different defects.

Molecular Simulation Approaches to the Study of Thermotropic and Lyotropic Liquid Crystals

Over the last decade, the availability of computer time, together with new algorithms capable of exploiting parallel computer architectures, has opened up many possibilities in molecularly modelling liquid crystalline systems. This perspective article points to recent progress in modelling both thermotropic and lyotropic systems. For thermotropic nematics, the advent of improved molecular force fields can provide predictions for nematic clearing temperatures within a 10 K range. Such studies also provide valuable insights into the structure of more complex phases, where molecular organisation may be challenging to probe experimentally. Developments in coarse-grained models for thermotropics are discussed in the context of understanding the complex interplay of molecular packing, microphase separation and local interactions, and in developing methods for the calculation of material properties for thermotropics. We discuss progress towards the calculation of elastic constants, rotational viscosity coefficients, flexoelectric coefficients and helical twisting powers. The article also covers developments in modelling micelles, conventional lyotropic phases, lyotropic phase diagrams, and chromonic liquid crystals. For the latter, atomistic simulations have been particularly productive in clarifying the nature of the self-assembled aggregates in dilute solution. The development of effective coarse-grained models for chromonics is discussed in detail, including models that have demonstrated the formation of the chromonic N and M phases.

Systematic bottom-up molecular coarse-graining via force and torque matching using anisotropic particles

We derive a systematic and general method for parametrizing coarse-grained molecular models consisting of anisotropic particles from fine-grained (e.g. all-atom) models for condensed-phase molecular dynamics simulations. The method, which we call anisotropic force-matching coarse-graining (AFM-CG), is based on rigorous statistical mechanical principles, enforcing consistency between the coarse-grained and fine-grained phase-space distributions to derive equations for the coarse-grained forces, torques, masses, and moments of inertia in terms of properties of a condensed-phase fine-grained system. We verify the accuracy and efficiency of the method by coarse-graining liquid-state systems of two different anisotropic organic molecules, benzene and perylene, and show that the parametrized coarse-grained models more accurately describe properties of these systems than previous anisotropic coarse-grained models parametrized using other methods that do not account for finite-temperature and many-body effects on the condensed-phase coarse-grained interactions. The AFM-CG method will be useful for developing accurate and efficient dynamical simulation models of condensed-phase systems of molecules consisting of large, rigid, anisotropic fragments, such as liquid crystals, organic semiconductors, and nucleic acids.

Simulating the nematic-isotropic phase transition of liquid crystal model via generalized replica-exchange method

The nematic-isotropic (NI) phase transition of 4-cyano-4'-pentylbiphenyl was simulated using the generalized replica-exchange method (gREM) based on molecular dynamics simulations. The effective temperature is introduced in the gREM, allowing for the enhanced sampling of configurations in the unstable region, which is intrinsic to the first-order phase transition. The sampling performance was analyzed with different system sizes and compared with that of the temperature replica-exchange method (tREM). It was observed that gREM is capable of sampling configurations at sufficient replica-exchange acceptance ratios even around the NI transition temperature. A bimodal distribution of the order parameter at the transition region was found, which is in agreement with the mean-field theory. In contrast, tREM is ineffective around the transition temperature owing to the potential energy gap between the nematic and isotropic phases.

Topological phase transitions in two-dimensional bent-core liquid crystal models

Two-dimensional liquid crystal (LC) models of interacting V-shaped bent-core molecules with two rigid rodlike identical segments connected at a fixed angle (θ=112^{∘}) are investigated. The model assigns equal biquadratic tensor coupling among constituents of the interacting neighboring molecules on a square lattice, allowing for reorientations in three dimensions (d=2, n=3). We find evidence of two temperature-driven topological transitions mediated by point disclinations associated with the three ordering directors, condensing the LC medium successively into uniaxial and biaxial phases. With Monte Carlo simulations, temperature dependencies of the system energy, specific heat, orientational order parameters, topological order parameters, and densities of unbound topological defects of the different ordering directors are computed. The high-temperature transition results in topological ordering of disclinations of the primary director, imparting uniaxial symmetry to the phase. The low-temperature transition precipitates simultaneous topological ordering of defects of the remaining directors, resulting in biaxial symmetry. The correlation functions, quantifying spatial variations of the orientational correlations of the molecular axes show exponential decays in the high-temperature (disordered) phase, and power-law decays in the low-temperature (biaxial) phase. Differing temperature dependencies of the topological parameters point to a significant degree of cross coupling among the uniaxial and biaxial tensors of interacting molecules. This simplified Hamiltonian leaves θ as the only free model parameter, and the system traces a θ-dependent trajectory in a plane of the phenomenological parameter space.

Lattice Monte Carlo study of orientational order in a confined system of biaxial particles: Effect of an external electric field

In this work we have used lattice Monte Carlo to determine the orientational order of a system of biaxial particles confined between two walls inducing perfect order and subjected to an electric field perpendicular to the walls. The particles are set to interact with their nearest neighbors through a biaxial version of the Lebwohl-Lasher potential. A particular set of values for the molecular reduced polarizabilities defining the potential used was considered; the Metropolis sampling algorithm was used in the Monte Carlo simulations. The relevant order parameters were determined in the middle plane of the sample and for some cases across the whole thickness of the sample. We have determined the temperature-electric field phase diagram for this system and found, as expected, five different system configurations corresponding to three different mesophases. At low temperatures and low fields the system finds itself in an undistorted biaxial phase. On increasing the field at low temperatures, a Freedericksz transition takes place and the secondary directors reorient towards the field while the primary director stays undistorted and parallel to the walls. On increasing the field further, a second Freedericksz transition occurs and the primary director orients also towards the field direction. The orientational order measured at the field strengths tested is not affected by the field. On increasing the temperature, a transition to a uniaxial phase occurs and within the range of this phase a field increase leads also to a Freedericksz transition where the main director reorients towards the field. At higher temperature a transition to a disordered phase is found. We have performed finite size scaling analysis for the Freedericksz critical fields and found that they scale with the distance L between the walls as L^{-1} as expected from continuum theory. From these fields we have also determined the temperature dependence of two elastic constant ratios. Critical exponents and critical temperatures for the order parameter and the correlation length for the biaxial-uniaxial phase transition and the uniaxial to disordered phase transition were also determined by finite size scaling and are discussed.

Quantification of nematic cell polarity in three-dimensional tissues

How epithelial cells coordinate their polarity to form functional tissues is an open question in cell biology. Here, we characterize a unique type of polarity found in liver tissue, nematic cell polarity, which is different from vectorial cell polarity in simple, sheet-like epithelia. We propose a conceptual and algorithmic framework to characterize complex patterns of polarity proteins on the surface of a cell in terms of a multipole expansion. To rigorously quantify previously observed tissue-level patterns of nematic cell polarity (Morales-Navarrete et al., eLife 2019), we introduce the concept of co-orientational order parameters, which generalize the known biaxial order parameters of the theory of liquid crystals. Applying these concepts to three-dimensional reconstructions of single cells from high-resolution imaging data of mouse liver tissue, we show that the axes of nematic cell polarity of hepatocytes exhibit local coordination and are aligned with the biaxially anisotropic sinusoidal network for blood transport. Our study characterizes liver tissue as a biological example of a biaxial liquid crystal. The general methodology developed here could be applied to other tissues and in-vitro organoids.

Topological transition in a two-dimensional three-vector model: A non-Boltzmann Monte Carlo study

Two-dimensional three-vector (d=2,n=3) lattice model of a liquid crystal (LC) system with order parameter space (R) described by the fundamental group Π_{1}(R)=Z_{2} was recently investigated based on non-Boltzmann Monte Carlo simulations. Its results indicated that the system did not undergo a topological transition condensing to a low temperature critical state as was reported earlier. Instead, a crossover to a nematic phase was observed, induced by the onset of a competing relevant length scale. This mechanism is further probed here by assigning a more restrictive R symmetry with Π_{1}(R)=Q (the discrete and non-Abelian group of quaternions), thus engaging the three spin degrees in the formation of point topological defects (disclinations). The results reported here indicate that such a choice of symmetry of the Hamiltonian with suitable model parameters leads to a defect-mediated transition to a low-temperature phase with topological order. It is characterized by a line of critical points with quasi-long-range order of its three spin degrees. The associated temperature-dependent power-law exponent decreases progressively and vanishes linearly as temperature tends to zero. The high-temperature disordered phase shows exponential spin correlations and their temperature-dependent lengths exhibit an essential singular divergence as the system approaches the topological transition point. Biaxial LC models have the required R symmetry owing to their tensor orientational orders and are suggested to serve as prototype examples to exhibit topological transition in (d=2,n=3) lattice models.

Ultrafast isomerization-induced cooperative motions to higher molecular orientation in smectic liquid-crystalline azobenzene molecules

The photoisomerization of molecules is widely used to control the structure of soft matter in both natural and synthetic systems. However, the structural dynamics of the molecules during isomerization and their subsequent response are difficult to elucidate due to their complex and ultrafast nature. Herein, we describe the ultrafast formation of higher-orientation of liquid-crystalline (LC) azobenzene molecules via linearly polarized ultraviolet light (UV) using ultrafast time-resolved electron diffraction. The ultrafast orientation is caused by the trans-to-cis isomerization of the azobenzene molecules. Our observations are consistent with simplified molecular dynamics calculations that revealed that the molecules are aligned with the laser polarization axis by their cooperative motion after photoisomerization. This insight advances the fundamental chemistry of photoresponsive molecules in soft matter as well as their ultrafast photomechanical applications.

Phase behavior of hard C_{2h}-symmetric particle systems

Using Monte Carlo numerical simulation, this work sketches the phase diagram of systems of certain hard C_{2h}-symmetric particles, formed by gluing two aligned and displaced hard spherocylinders with a cylindrical-length-to-diameter ratio realistically, if viewed not only from the lyotropic colloidal liquid-crystal side but also from the thermotropic low-molecular-mass liquid-crystal side, equal to 5, as a function of the displacement. Several distinctive phases are observed, such as a nonperiodic smectic-B-like phase, a nonperiodic smectic-H-like phase, a smectic-C phase, and a short-layer-spacing uniaxial smectic-A phase but no biaxial nematic phase.

Experimental Conditions for the Stabilization of the Lyotropic Biaxial Nematic Mesophase

Nematic phases are some of the most common phases among the lyotropic liquid crystalline structures. They have been widely investigated during last decades. In early studies, two uniaxial nematic phases (discotic, ND, and calamitic, NC) were identified. After the discovery of the third one, named biaxial nematic phase (NB) in 1980, however, some controversies in the stability of biaxial nematic phases began and still continue in the literature. From the theoretical point of view, the existence of a biaxial nematic phase is well established. This review aims to bring information about the historical development of those phases considering the early studies and then summarize the recent studies on how to stabilize different nematic phases from the experimental conditions, especially, choosing the suitable constituents of lyotropic mixtures.

Generic Model for Lamellar Self-Assembly in Conjugated Polymers: Linking Mesoscopic Morphology and Charge Transport in P3HT

We develop a generic coarse-grained model of soluble conjugated polymers, capable of describing their self-assembly into a lamellar mesophase. Polymer chains are described by a hindered-rotation model, where interaction centers represent entire repeat units, including side chains. We introduce soft anisotropic nonbonded interactions to mimic the potential of mean force between atomistic repeat units. The functional form of this potential reflects the symmetry of the molecular order in a lamellar mesophase. The model can generate both nematic and lamellar (sanidic smectic) molecular arrangements. We parametrize this model for a soluble conjugated polymer poly(3-hexylthiophene) (P3HT) and demonstrate that the simulated lamellar mesophase matches morphologies of low molecular weight P3HT, experimentally observed at elevated temperatures. A qualitative charge-transport model allows us to link local chain conformations and mesoscale order to charge transport. In particular, it shows how coarsening of lamellar domains and chain extension increase the charge carrier mobility. By modeling large systems and long chains, we can capture transport between lamellar layers, which is due to rare, but thermodynamically allowed, backbone bridges between neighboring layers.

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.

Complex free-energy landscapes in biaxial nematic liquid crystals and the role of repulsive interactions: A Wang-Landau study

General quadratic Hamiltonian models, describing the interaction between liquid-crystal molecules (typically with D_{2h} symmetry), take into account couplings between their uniaxial and biaxial tensors. While the attractive contributions arising from interactions between similar tensors of the participating molecules provide for eventual condensation of the respective orders at suitably low temperatures, the role of cross coupling between unlike tensors is not fully appreciated. Our recent study with an advanced Monte Carlo technique (entropic sampling) showed clearly the increasing relevance of this cross term in determining the phase diagram (contravening in some regions of model parameter space), the predictions of mean-field theory, and standard Monte Carlo simulation results. In this context, we investigated the phase diagrams and the nature of the phases therein on two trajectories in the parameter space: one is a line in the interior region of biaxial stability believed to be representative of the real systems, and the second is the extensively investigated parabolic path resulting from the London dispersion approximation. In both cases, we find the destabilizing effect of increased cross-coupling interactions, which invariably result in the formation of local biaxial organizations inhomogeneously distributed. This manifests as a small, but unmistakable, contribution of biaxial order in the uniaxial phase. The free-energy profiles computed in the present study as a function of the two dominant order parameters indicate complex landscapes. On the one hand, these profiles account for the unusual thermal behavior of the biaxial order parameter under significant destabilizing influence from the cross terms. On the other, they also allude to the possibility that in real systems, these complexities might indeed be inhibiting the formation of a low-temperature biaxial order itself-perhaps reflecting the difficulties in their ready realization in the laboratory.

Dense crystalline packings of ellipsoids

An ellipsoid, the simplest nonspherical shape, has been extensively used as a model for elongated building blocks for a wide spectrum of molecular, colloidal, and granular systems. Yet the densest packing of congruent hard ellipsoids, which is intimately related to the high-density phase of many condensed matter systems, is still an open problem. We discover an unusual family of dense crystalline packings of self-dual ellipsoids (ratios of the semiaxes α:sqrt[α]:1), containing 24 particles with a quasi-square-triangular (SQ-TR) tiling arrangement in the fundamental cell. The associated packing density ϕ exceeds that of the densest known SM2 crystal [ A. Donev et al., Phys. Rev. Lett. 92, 255506 (2004)10.1103/PhysRevLett.92.255506] for aspect ratios α in (1.365, 1.5625), attaining a maximal ϕ≈0.75806... at α=93/64. We show that the SQ-TR phase derived from these dense packings is thermodynamically stable at high densities over the aforementioned α range and report a phase diagram for self-dual ellipsoids. The discovery of the SQ-TR crystal suggests organizing principles for nonspherical particles and self-assembly of colloidal systems.

Multi-step modeling of liquid crystals using ab initio molecular packing and hybrid quantum mechanics/molecular mechanics simulations

A density functional theory (DFT) based multi-step simulation method is used to characterize the detailed molecular structure and inter/intra- molecular interactions of two benchmark liquid crystals (LC) 5CB, 8CB and a novel tri-biphenyl ring bent core LC material. The method uses hybrid DFT at the B3LYP/6-31G* level to obtain molecular structure and Raman data. These results are fed to a crystal packing simulation to find possible crystal structures. A pico-second quantum mechanics/molecular mechanics (QM/MM) simulation model is built for the selected structures with lower overall energy as well as optimal density. The stabilized crystal structures are then extended into a super cell, heated and simulated using a mixed force field and nano-second molecular dynamics (MD). The described simulation process sequence provides predictions of molecular Raman spectrum, LC density, isotropic depolarization ratio, ratio of differential polarizability, order parameters, molecular structures, and rotating Raman spectrum of the different mesophases. The Raman spectra, order parameters and depolarization ratios all agree well with existing experimental and previous simulation results. The study of the novel tri-biphenyl ring bent core LC system shows that the ratio of differential polarizability depends on intra-molecular interactions. The findings presented in this manuscript contribute to the on-going efforts to establish links between LC molecular structures and their properties, including optical behavior.

Maier-Saupe model for a mixture of uniaxial and biaxial molecules

We introduce shape variations in a liquid-crystalline system by considering an elementary Maier-Saupe lattice model for a mixture of uniaxial and biaxial molecules. Shape variables are treated in the annealed (thermalized) limit. We analyze the thermodynamic properties of this system in terms of temperature T, concentration c of intrinsically biaxial molecules, and a parameter Δ associated with the degree of biaxiality of the molecules. At the mean-field level, we use standard techniques of statistical mechanics to draw global phase diagrams, which are shown to display a rich structure, including uniaxial and biaxial nematic phases, a reentrant ordered region, and many distinct multicritical points. Also, we use the formalism to write an expansion of the free energy in order to make contact with the Landau-de Gennes theory of nematic phase transitions.

Reexamination of the mean-field phase diagram of biaxial nematic liquid crystals: Insights from Monte Carlo studies

Investigations of the phase diagram of biaxial liquid-crystal systems through analyses of general Hamiltonian models within the simplifications of mean-field theory (MFT), as well as by computer simulations based on microscopic models, are directed toward an appreciation of the role of the underlying molecular-level interactions to facilitate its spontaneous condensation into a nematic phase with biaxial symmetry. Continuing experimental challenges in realizing such a system unambiguously, despite encouraging predictions from MFT, for example, are requiring more versatile simulational methodologies capable of providing insights into possible hindering barriers within the system, typically gleaned through its free-energy dependences on relevant observables as the system is driven through the transitions. The recent paper from this group [Kamala Latha et al., Phys. Rev. E 89, 050501(R) (2014)], summarizing the outcome of detailed Monte Carlo simulations carried out employing an entropic sampling technique, suggested a qualitative modification of the MFT phase diagram as the Hamiltonian is asymptotically driven toward the so-called partly repulsive regions. It was argued that the degree of (cross) coupling between the uniaxial and biaxial tensor components of neighboring molecules plays a crucial role in facilitating a ready condensation of the biaxial phase, suggesting that this could be a plausible factor in explaining the experimental difficulties. In this paper, we elaborate this point further, providing additional evidence from curious variations of free-energy profiles with respect to the relevant orientational order parameters, at different temperatures bracketing the phase transitions.

Nanosecond electro-optics of a nematic liquid crystal with negative dielectric anisotropy

We study a nanosecond electro-optic response of a nematic liquid crystal in a geometry where an applied electric field E modifies the tensor order parameter but does not change the orientation of the optic axis (director N ̂). We use a nematic with negative dielectric anisotropy with the electric field applied perpendicularly to N ̂. The field changes the dielectric tensor at optical frequencies (optic tensor) due to the following mechanisms: (a) nanosecond creation of the biaxial orientational order, (b) uniaxial modification of the orientational order that occurs over time scales of tens of nanoseconds, and (c) the quenching of director fluctuations with a wide range of characteristic times up to milliseconds. We develop a model to describe the dynamics of all three mechanisms. We design the experimental conditions to selectively suppress the contributions from the quenching of director fluctuations (c) and from the biaxial order effect (a) and thus, separate the contributions of the three mechanisms in the electro-optic response. As a result, the experimental data can be well fitted with the model parameters. The analysis provides a rather detailed physical picture of how the liquid crystal responds to a strong electric field on a time scale of nanoseconds. The paper provides a useful guidance in the current search for the biaxial nematic phase. Namely, the temperature dependence of the biaxial susceptibility allows one to estimate the temperature of the potential uniaxial-to-biaxial phase transition. An analysis of the quenching of director fluctuations indicates that on a time scale of nanoseconds, the classic model with constant viscoelastic material parameters might reach its limit of validity. The effect of nanosecond electric modification of the order parameter can be used in applications in which one needs to achieve ultrafast (nanosecond) changes in optical characteristics, such as birefringence.

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.

Model of hard spheroplatelets near a hard wall

A system of hard spheroplatelets near an impenetrable wall is studied in the low-density Onsager approximation. Spheroplatelets have optimal shape between rods and plates, and the direct transition from the isotropic to biaxial nematic phase is present. A simple local approximation for the one-particle distribution function is used. Analytical results for the surface tension and the entropy contributions are derived. The density and the order-parameter profiles near the wall are calculated. The preferred orientation of the short molecule axes is perpendicular to the wall. Biaxiality close to the wall can appear only if the phase is biaxial in the bulk.

Detection of an intermediate biaxial phase in the phase diagram of biaxial liquid crystals: entropic sampling study

We investigate the phase sequence of biaxial liquid crystals, based on a general quadratic model Hamiltonian over the relevant parameter space, with a Monte Carlo simulation which constructs equilibrium ensembles of microstates, overcoming possible (free) energy barriers (combining entropic and frontier sampling techniques). The resulting phase diagram qualitatively differs from the universal phase diagram predicted earlier from mean-field theory (MFT), as well as the Monte Carlo simulations with the Metropolis algorithm. The direct isotropic-to-biaxial transition predicted by the MFT is replaced in certain regions of the space by the onset of an additional intermediate biaxial phase of very low order, leading to the sequence N(B)-N(B1)-I. This is due to inherent barriers to fluctuations of the components comprising the total energy, and may explain the difficulties in the experimental realization of these phases.

Molecular field theory for biaxial smectic A liquid crystals

Thermotropic biaxial nematic phases seem to be rare, but biaxial smectic A phases less so. Here we use molecular field theory to study a simple two-parameter model, with one parameter promoting a biaxial phase and the second promoting smecticity. The theory combines the biaxial Maier-Saupe and McMillan models. We use alternatively the Sonnet-Virga-Durand (SVD) and geometric mean approximations (GMA) to characterize molecular biaxiality by a single parameter. For non-zero smecticity and biaxiality, the model always predicts a ground state biaxial smectic A phase. For a low degree of smectic order, the phase diagram is very rich, predicting uniaxial and biaxial nematic and smectic phases, with the addition of a variety of tricritical and tetracritical points. For higher degrees of smecticity, the region of stability of the biaxial nematic phase is restricted and eventually disappears, yielding to the biaxial smectic phase. Phase diagrams from the two alternative approximations for molecular biaxiality are similar, except inasmuch that SVD allows for a first-order isotropic-biaxial nematic transition, whereas GMA predicts a Landau point separating isotropic and biaxial nematic phases. We speculate that the rarity of thermotropic biaxial nematic phases is partly a consequence of the presence of stabler analogous smectic phases.

Nanosecond electro-optic switching of a liquid crystal

Electrically induced reorientation of nematic liquid crystal (NLC) molecules caused by dielectric anisotropy of the material is a fundamental phenomenon widely used in modern technologies. Its Achilles heel is a slow (millisecond) relaxation from the field-on to the field-off state. We present an electro-optic effect in an NLC with a response time of about 30 ns to both the field-on and field-off switching. This effect is caused by the electric field induced modification of the order parameters and does not require reorientation of the optic axis (director).

Calamitic and antinematic orientational order produced by the generalized Straley lattice model

We consider here a classical model, consisting of D_{2h}-symmetric particles in a three-dimensional simple-cubic lattice; the pair potential is isotropic in orientation space, and restricted to nearest neighbors. The simplest potential model is written in terms of the squares of the scalar products between unit vectors describing the three interacting arms of the molecules, as proposed in previous literature. Two predominant antinematic couplings of equal strength (+1) are perturbed by a comparatively weaker calamitic one, parameterized by a coupling constant -z ranging in [-1,0]. This choice rules out thermodynamically stable phases endowed with macroscopic biaxiality. The antinematic terms favor states with the corresponding molecular axes mutually orthogonal. Although the low-temperature phase of the special case with null calamitic term (PP0) is uniaxial and antinematically ordered, in the general case presented here both Monte Carlo and molecular-field approaches show that, for z close to zero, the models exhibit a low-temperature uniaxial nematic phase, followed by an antinematic one, and finally by the orientationally disordered one. On the other hand, for sufficiently large values of z, we only find evidence of uniaxial calamitic behavior, as expected by following the limiting cases.

Biaxial order parameter in the homologous series of orthogonal bent-core smectic liquid crystals

The fundamental parameter of the uniaxial liquid crystalline state that governs nearly all of its physical properties is the primary orientational order parameter (S) for the long axes of molecules with respect to the director. The biaxial liquid crystals (LCs) possess biaxial order parameters depending on the phase symmetry of the system. In this paper we show that in the first approximation a biaxial orthogonal smectic phase can be described by two primary order parameters: S for the long axes and C for the ordering of the short axes of molecules. The temperature dependencies of S and C are obtained by the Haller's extrapolation technique through measurements of the optical birefringence and biaxiality on a nontilted polar antiferroelectric (Sm-AP(A)) phase of a homologous series of LCs built from the bent-core achiral molecules. For such a biaxial smectic phase both S and C, particularly the temperature dependency of the latter, are being experimentally determined. Results show that S in the orthogonal smectic phase composed of bent cores is higher than in Sm-A calamatic LCs and C is also significantly large.

Free-energy power expansion for orientationally ordered phases: energy and entropy

We propose an approach to the description of orientational phase transitions that utilizes the specific features of orientational energy and entropy. The approach is applied to building a model for nematic phases in materials with nonpolar parallelepiped-type molecules with symmetry D(2h). The model operates with complex order parameters, generalizes the Landau-de Gennes theory, and predicts the existence of a biaxial nematic phase for the fourth-order expansion of free energy.

Antinematic orientational order produced by an extreme case of the generalized Straley lattice model

We address here a special, extreme case of the quadratic pair interaction potential between classical, D(2h)-symmetric particles (the generalized Straley model) on a three-dimensional simple cubic lattice. The model involves predominant antinematic couplings and it has been studied by Monte Carlo simulation and a molecular field treatment. The obtained results show a second-order transition between the isotropic phase and the low-temperature one, exhibiting uniaxial antinematic order.

Depletion-induced biaxial nematic states of boardlike particles

With the aim of investigating the stability conditions of biaxial nematic liquid crystals, we study the effect of adding a non-adsorbing ideal depletant on the phase behavior of colloidal hard boardlike particles. We take into account the presence of the depletant by introducing an effective depletion attraction between a pair of boardlike particles. At fixed depletant fugacity, the stable liquid-crystal phase is determined through a mean-field theory with restricted orientations. Interestingly, we predict that for slightly elongated boardlike particles a critical depletant density exists, where the system undergoes a direct transition from an isotropic liquid to a biaxial nematic phase. As a consequence, by tuning the depletant density, an easy experimental control parameter, one can stabilize states of high biaxial nematic order even when these states are unstable for pure systems of boardlike particles.

Tetrahedratic mesophases, chiral order, and helical domains induced by quadrupolar and octupolar interactions

We present an exhaustive account of phases and phase transitions that can be stabilized in the recently introduced generalized Lebwohl-Lasher model with quadrupolar and octupolar microscopic interactions [L. Longa, G. Pająk, and T. Wydro, Phys. Rev. E 79, 040701(R) (2009)]. A complete mean-field analysis of the model, along with Monte Carlo simulations allows us to identify four distinct classes of the phase diagrams with a number of multicritical points where, in addition to the standard uniaxial and biaxial nematic phases, the other nematic like phases are stabilized. These involve, among the others, tetrahedratic (T), nematic tetrahedratic (N(T)), and chiral nematic tetrahedratic (N(T)(*)) phases of global T(d), D(2d), and D(2) symmetry, respectively. Molecular order parameters and correlation functions in these phases are determined. We conclude with generalizations of the model that give a simple molecular interpretation of macroscopic regions with opposite optical activity (ambidextrous chirality), observed, e.g., in bent-core systems. An estimate of the helical pitch in the N(T)(*) phase is also given.

Influence of dipole-dipole correlations on the stability of the biaxial nematic phase in the model bent-core liquid crystal

A molecular theory of biaxial nematic ordering in the system of bent-core molecules has been developed in the two-particle cluster approximation which enables one to take into account short-range polar correlations determined by both electrostatic dipole-dipole interaction and polar molecular shape. All orientational order parameters and short-range correlation functions are calculated numerically as functions of temperature in the uniaxial and in the biaxial nematic phases, and the results are compared with the ones obtained in the mean-field approximation and in the cluster approximation but without taking into consideration the dipole-dipole interaction. It is shown that short-range polar correlations and, in particular, the dipole-dipole correlations dramatically increase the temperature of the transition into the biaxial nematic phase and enhancing its stability range. The results are also very sensitive to the value of the opening angle of a model bent-core molecule.

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.

Biaxial nematic phase in the Maier-Saupe model for a mixture of discs and cylinders

We analyze the global phase diagram of a Maier-Saupe lattice model with the inclusion of shape-disordered degrees of freedom to mimic a mixture of oblate and prolate molecules (discs and cylinders). In the neighborhood of a Landau multicritical point, solutions of the statistical problem can be written as a Landau-de Gennes expansion for the free energy. If the shape-disordered degrees of freedom are quenched, we confirm the existence of a biaxial nematic structure. If orientational and disorder degrees of freedom are allowed to thermalize, this biaxial solution becomes thermodynamically unstable. Also, we use a two-temperature formalism to mimic the presence of two distinct relaxation times, and show that a slight departure from complete thermalization is enough to stabilize a biaxial nematic phase.