Measuring liquid crystal elastic constants with free energy perturbations

On the elusive saddle-splay and splay-bend elastic constants of nematic liquid crystals

The elastic behavior of nematics is commonly described in terms of the three so-called bulk deformation modes, i.e., splay, twist, and bend. However, the elastic free energy contains also other terms, often denoted as saddle-splay and splay-bend, which contribute, for instance, in confined systems. The role of such terms is controversial, partly because of the difficulty of their experimental determination. The saddle-splay (K24) and splay-bend (K13) elastic constants remain elusive also for theories; indeed, even the possibility of obtaining unambiguous microscopic expressions for these quantities has been questioned. Here, within the framework of Onsager theory with Parsons-Lee correction, we obtain microscopic estimates of the deformation free energy density of hard rod nematics in the presence of different director deformations. In the limit of a slowly changing director, these are directly compared with the macroscopic elastic free energy density. Within the same framework, we derive also closed microscopic expressions for all elastic coefficients of rodlike nematics. We find that the saddle-splay constant K24 is larger than both K11 and K22 over a wide range of particle lengths and densities. Moreover, the K13 contribution comes out to be crucial for the consistency of the results obtained from the analysis of the microscopic deformation free energy density calculated for variants of the splay deformation.

Structures, thermodynamics and dynamics of topological defects in Gay-Berne nematic liquid crystals

Topological defects are a ubiquitous phenomenon across different physical systems. A better understanding of defects can be helpful in elucidating the physical behaviors of many real materials systems. In nematic liquid crystals, defects exhibit unique optical signatures and can segregate impurities, showing their promise as molecular carriers and nano-reactors. Continuum theory and simulations have been successfully applied to link static and dynamical behaviors of topological defects to the material constants of the underlying nematic. However, further evidence and molecular details are still lacking. Here we perform molecular dynamics simulations of Gay-Berne particles, a model nematic, to examine the molecular structures and dynamics of +1/2 defects in a thin-film nematic. Specifically, we measure the bend-to-splay ratio K₃/K₁ using two independent, indirect measurements, showing good agreement. Next, we study the annihilation event of a pair of ±1/2 defects, of which the trajectories are consistent with experiments and hydrodynamic simulations. We further examine the thermodynamics of defect annihilation in an NVE ensemble, leading us to correctly estimate the elastic modulus by using the energy conservation law. Finally, we explore effects of defect annihilation in regions of nonuniform temperature within these coarse-grained molecular models which cannot be analysed by existing continuum level simulations. We find that +1/2 defects tend to move toward hotter areas and their change of speed in a temperature gradient can be quantitatively understood through a term derived from the temperature dependence of the elastic modulus. As such, our work has provided molecular insights into structures and dynamics of topological defects, presented unique and accessible methods to measure elastic constants by inspecting defects, and proposed an alternative control parameter of defects using temperature gradient.

Characterization of optically thin cells and experimental liquid crystals

The current development of new liquid crystal devices often requires the use of thin cells and new experimental materials. Characterizing these devices and materials with optical methods can be challenging if (1) the total phase lag is small ("thin cells") or (2) the liquid crystal optical and dielectric properties are only partially known. We explore the limitations of these two challenges for efficient characterization and assessment of new, to the best of our knowledge, liquid crystal devices. We show that it is possible to extract a wealth of liquid crystal parameters even for cells with a phase lag of ΔΦ≈π, such as E7 liquid crystal in a 1.5 µm cell, using cross-polarized intensity measurements. The reliability of the optical method is also demonstrated for liquid crystals without precise values of dielectric or refractive index coefficients.

Elastic constants of biological filamentous colloids: estimation and implications on nematic and cholesteric tactoid morphologies

Biological liquid crystals, originating from the self-assembly of biological filamentous colloids, such as cellulose and amyloid fibrils, show a complex lyotropic behaviour that is extremely difficult to predict and characterize. Here we analyse the liquid crystalline phases of amyloid fibrils, and sulfated and carboxylated cellulose nanocrystals and measure their Frank-Oseen elastic constants K1, K2 and K3 by four different approaches. The first two approaches are based on the benchmark of the predictions of: (i) a scaling form and (ii) a variational form of the Frank-Oseen energy functional with the experimental critical volumes at order-order liquid crystalline transitions of the tactoids. The third and the fourth methods imply: (iii) the direct scaling equations of elastic constants and (iv) a molecular theory predicting the elastic constants from the experimentally accessible contour length distributions of the filamentous colloids. These three biological systems exhibit diverse liquid crystalline behaviour, governed by the distinct elastic constants characterizing each colloid. Differences and similarities among the three systems are highlighted and interpreted based on the present analysis, providing a general framework to study dispersed liquid crystalline systems.

Interplay of Particle Morphology and Director Distortions in Nematic Fluids

The existing microscopic theories for elasticity of nematics are challenged by recent findings on systems, whether bent molecules or semiflexible polymers, which do not comply with the model of rigid rodlike particles. Here, we propose an extension of Onsager-Straley second-virial theory, based on a model for the orientational distribution function that, through explicit account of the director profile along a particle, changes in the presence of deformations. The elastic constants reveal specific effects of particle morphology, which are not captured by the existing theories. This paves the way to microscopic modeling of the elastic properties of semiflexible liquid crystal polymers, which is a longstanding issue.

Electric field-induced crossover from 3D to 2D topological defects in a nematic liquid crystal: experimental verification

A substrate was patterned with two pairs of half-integer strength topological defects, (+1/2, +1/2) and (+1/2, -1/2). In a sufficiently thick cell, a disclination line runs in an arch above the substrate connecting the two half integer defects within each pair. The director around the disclination line for the like-sign pair must rotate in 3D, whereas for the opposite-sign defect pair the director lies in the xy-plane parallel to the substrate. For a negative dielectric anisotropy nematic, an electric field applied normal to the substrate drives the director into the xy-plane, forcing the arch of the disclination line of the like-sign pair to become extended along the z-axis. For sufficiently large field the arch splits, resulting in two nearly parallel disclination lines traversing the cell from one substrate to the other. The opposite-sign defect pair is largely unaffected by the electric field as the director already lies in the xy-plane. Experimental results are presented, which are consistent with numerical simulations.

Novel elastic response in twist-bend nematic models

Bent-shaped liquid crystals have attracted significant attention recently due to their novel mesostructure and the intriguing behavior of their elastic constants, which are strongly anisotropic and have an unusual temperature dependence. Though theories explain the onset of the twist-bend nematic phase (NTB) through spontaneous symmetry breaking concomitant with transition to a negative bend (K3) elastic constant, this has not been observed as yet in experiments. There, the small bend elastic constant has a strongly non-monotonic temperature dependence, which first increases after crossing the isotropic (I)-nematic (N) transition, then dips near the nematic (N)-twist-bend (NTB) transition before it increases again as the transition is crossed. The molecular mechanisms responsible for this exotic behavior are unclear. Here, we utilize density of states algorithms in Monte Carlo simulation applied to a variant of the Lebwohl-Lasher model which includes bent-shaped-like interactions to analyze the mechanism behind elastic response in this novel mesostructure and understand the temperature dependence of its Frank-Oseen elastic constants.

Elasticity and yielding of mesophases of block copolymers in water-oil mixtures

Amphiphilic block copolymers self-assemble at the water/oil interface to form different mesomorphic structures, such as lamellar, micellar cubic, normal hexagonal, and reverse hexagonal structures. Usually, these structures are polycrystalline and the value of their elastic modulus depends on the average orientation of their constituent's single crystals. We provide a model to predict the elastic modulus and yielding of mesophases from their characteristic length and intermicellar interactions. Shear modulus of each structure is calculated as a function of deformation (strain). Zero-shear modulus, G0, depends on the inverse of the intermicellar distance with a power law model. The power law index for each structure is approximately n + 2 where n is the degree of confinement in the mesophase: 1 for lamellar, 2 for both normal and reverse hexagonal, and 3 for micellar cubic structures. Rheological properties of different mesophases of Pluronic P84 in the presence of water and p-xylene are used as a case study. The model is found to be in good agreement with experimental data in the linear viscoelastic region. When compared to experimental data, the yield strain value obtained from the model is one order of magnitude higher than the limit of the linear viscoelastic regime and close to the strain at the cross-over point of storage and loss moduli. Frequency sweep measurements are done to characterize the relaxation and cooperative model behaviors of each mesophase structure.

Effects of chiral dopants on double-twist configurations of lyotropic chromonic liquid crystals in a cylindrical cavity

We investigate how chiral dopants affect the chiral symmetry breaking of lyotropic chromonic liquid crystals (LCLCs) focusing on the double-twist (DT) director configurations in a cylindrical capillary. LCLCs of unusual elastic properties tend to exhibit chiral director configurations under confinement despite the absence of intrinsic chirality. The DT director configuration in a cylindrical cavity with a degenerate planar anchoring, resulting from the large saddle-splay-to-twist elastic modulus ratio, is a representative example. Here we start by reexamining the DT configuration of nematic disodium cromoglycate in a cylindrical capillary and estimate the ratio of saddle splay to bend modulus K_{24}/K_{3}=0.5±0.1. Additionally, we study the DT configurations of the chiral nematic LCLCs with chiral dopants. The DT configuration becomes homochiral when the dopant concentration surpasses the critical concentration. We characterize these chiral DT configurations and provide a theoretical model on their energetics. Finally, we observe how the enantiomeric excess of chiral dopants determines the director configuration when dopants of two different handednesses are mixed.

Fluctuations and phase transitions of uniaxial and biaxial liquid crystals using a theoretically informed Monte Carlo and a Landau free energy density

In this work, we explore fluctuations during phase transitions of uniaxial and biaxial liquid crystals using a phenomenological free energy functional. We rely on a continuum-level description of the liquid crystal ordering with a tensorial parameter and a temperature dependent Landau polynomial expansion of the tensor's invariants. The free energy functional, over a three-dimensional periodic domain, is integrated with a Gaussian quadrature and minimized with a theoretically informed Monte Carlo method. We reconstruct analytical phase diagrams, following Landau and Doi's notations, to verify that the free energy relaxation reaches the global minimum. Importantly, our relaxation method is able to follow the thermodynamic behavior provided by other non-phenomenological approaches; we predict the first order character of the isotropic-nematic transition, and we identify the uniaxial-biaxial transition as second order. Finally, we use a finite-size scaling method, using the nematic susceptibility, to calculate the transition temperatures for 4-Cyano-4'-pentylbiphenyl (5CB) and N-(4-methoxybenzylidene)-4-butylaniline (MBBA). Our results show good agreement with experimental values, thereby validating our minimization method. Our approach is an alternative towards the relaxation of temperature dependent continuum-level free energy functionals, in any geometry, and can incorporate complicated elastic and surface energy densities.

The interplay between activity and filament flexibility determines the emergent properties of active nematics

Active nematics are microscopically driven liquid crystals that exhibit dynamical steady states characterized by the creation and annihilation of topological defects. Motivated by differences between previous simulations of active nematics based on rigid rods and experimental realizations based on semiflexible biopolymer filaments, we describe a large-scale simulation study of a particle-based computational model that explicitly incorporates filament semiflexibility. We find that energy injected into the system at the particle scale preferentially excites bend deformations, reducing the apparent filament bend modulus. The emergent characteristics of the active nematic depend on activity and flexibility only through this activity-renormalized bend 'modulus', demonstrating that apparent values of material parameters, such as the Frank 'constants', depend on activity. Thus, phenomenological parameters within continuum hydrodynamic descriptions of active nematics must account for this dependence. Further, we present a systematic way to estimate these parameters from observations of deformation fields and defect shapes in experimental or simulation data.

Self-assembly of semiflexible polymers confined to thin spherical shells

Confinement effects are critical for stiff macromolecules in biological cells, vesicles, and other systems in soft matter. For these molecules, the competition between the packing entropy and the enthalpic cost of bending is further shaped by strong confinement effects. Through coarse-grained molecular dynamics simulations, we explore the self-assembly of semiflexible polymers confined in thin spherical shells for various chain lengths, chain stiffnesses, and shell thicknesses. Here, we focus on the case where the contour and persistence length of the polymers are comparable to the radius of the confining cavity. The range of ordered structures is analyzed using several order parameters to elucidate the nature of orientational ordering in different parameter regimes. Previous simulations have revealed the emergence of bipolar and quadrupolar topological defects on the surface when the entire cavity was filled with a concentrated polymer solution [Phys. Rev. Lett., 2017, 118, 217803]. In contrast, spherical shell confinement restricts the appearance of a bipolar order. Instead, only the extent of the quadrupolar order changes with chain stiffness, as evidenced by the relative motion of topological defects. In the case of monolayers, we observe a nematic to smectic transition accompanied by a change in the nematic grain-size distribution as the contour length was decreased.

In Silico Measurement of Elastic Moduli of Nematic Liquid Crystals

Experiments on confined droplets of the nematic liquid crystal 5CB have questioned long-established bounds imposed on the elastic free energy of nematic systems. This elasticity, which derives from molecular alignment within nematic systems, is quantified through a set of moduli which can be difficult to measure experimentally and, in some cases, can only be probed indirectly. This is particularly true of the surfacelike saddle-splay elastic term, for which the available experimental data indicate values on the cusp of stability, often with large uncertainties. Here, we demonstrate that all nematic elastic moduli, including the saddle-splay elastic constant k_{24}, may be calculated directly from atomistic molecular simulations. Importantly, results obtained through in silico measurements of the 5CB elastic properties demonstrate unambiguously that saddle-splay elasticity alone is unable to describe the observed confined morphologies.

The Impact of Colloidal Surface-Anchoring on the Smectic A Phase

Liquid-crystalline phases are known for their unique properties, i.e., the combination of fluidity and long-range orientational and/or positional order. The presence of a colloidal particle gives rise to perturbations of this order locally. These perturbations are the origin of intercolloidal forces driving the colloidal self-assembly in a directed manner. Hence, the understanding of these perturbations is the first step in understanding and controlling the self-assembly process. Here, we perform Monte Carlo simulations to investigate the perturbations of orientational and positional order in a smectic A phase caused by a spherical colloid. We model the host phase via an interaction potential that reproduces characteristic features of phase behavior, structure, dynamics, and elasticity [S. Püschel-Schlotthauer et al. J. Chem. Phys. 2016, 145, 164903]. For strong homeotropic anchoring conditions, we find a Saturn ring defect and an onion structure in the smectic A phase; the latter has never been reported for colloids so far. For strong planar anchoring conditions, we find Boojum defects that become elongated at low temperature, similar to what is observed in experiments. However, for weak planar anchoring conditions, a double surface ring defect is exhibited in the smectic A phase.

Retrieving the saddle-splay elastic constant K24 of nematic liquid crystals from an algebraic approach

The physics of light interference experiments is well established for nematic liquid crystals. Using well-known techniques, it is possible to obtain important quantities, such as the differential scattering cross section and the saddl-splay elastic constant K24. However, the usual methods to retrieve the latter involve adjusting of computational parameters through visual comparisons between the experimental light interference pattern or a (2) H-NMR spectral pattern produced by an escaped-radial disclination, and their computational simulation counterparts. To avoid such comparisons, we develop an algebraic method for obtaining of saddle-splay elastic constant K24. Considering an escaped-radial disclination inside a capillary tube with radius R0 of tens of micrometers, we use a metric approach to study the propagation of the light (in the scalar wave approximation), near the surface of the tube and to determine the light interference pattern due to the defect. The latter is responsible for the existence of a well-defined interference peak associated to a unique angle [Formula: see text] . Since this angle depends on factors such as refractive indexes, curvature elastic constants, anchoring regime, surface anchoring strength and radius R0, the measurement of [Formula: see text] from the interference experiments involving two different radii allows us to algebraically retrieve K24. Our method allowed us to give the first reported estimation of K24 for the lyotropic chromonic liquid crystal Sunset Yellow FCF: K 24 = 2.1 pN.

Elastic response and phase behavior in binary liquid crystal mixtures

Utilizing density-of-states simulations, we perform a full mapping of the phase behavior and elastic responses of binary liquid crystalline mixtures represented by the multicomponent Lebwohl-Lasher model. Our techniques are able to characterize the complete phase diagram, including nematic-nematic phase separation predicted by mean-field theories, but previously not observed in simulations. Mapping this phase diagram permits detailed study of elastic properties across the miscible nematic region. Importantly, we observe for the first time local phase separation and disordering driven by the application of small linear perturbations near the transition temperature and more significantly through nonlinear stresses. These findings are of key importance in systems of blended nematics which contain particulate inclusions, or are otherwise confined.

Relevance of saddle-splay elasticity in complex nematic geometries

We demonstrate the relevance of saddle-splay elasticity in nematic liquid crystalline fluids in the context of complex surface anchoring conditions and the complex geometrical confinement. Specifically, nematic cells with patterns of surface anchoring and colloidal knots are shown as examples where saddle-splay free energy contribution can have a notable role which originates from nonhomogeneous surface anchoring and the varying surface curvature. Patterned nematic cells are shown to exhibit various (meta)stable configurations of nematic field, with relative (meta)stability depending on the saddle-splay. We show that for high enough values of saddle-splay elastic constant K24 a previously unstable conformation can be stabilised, more generally indicating that the saddle-splay can reverse or change the (meta)stability of various nematic structures affecting their phase diagrams. Furthermore, we investigate saddle-splay elasticity in the geometry of highly curved boundaries - the colloidal particle knots in nematic - where the local curvature of the particles induces complex spatial variations of the saddle-splay contributions. Finally, a nematic order parameter tensor based saddle-splay invariant is shown, which allows for the direct calculation of saddle-splay free energy from the Q-tensor, a possibility very relevant for multiple mesoscopic modelling approaches, such as Landau-de Gennes free energy modelling.

Chiral structures and defects of lyotropic chromonic liquid crystals induced by saddle-splay elasticity

An experimental and theoretical study of lyotropic chromonic liquid crystals (LCLCs) confined in cylinders with degenerate planar boundary conditions elucidates LCLC director configurations. When the Frank saddle-splay modulus is more than twice the twist modulus, the ground state adopts an inhomogeneous escaped-twisted configuration. Analysis of the configuration yields a large saddle-splay modulus, which violates Ericksen inequalities but not thermodynamic stability. Lastly, we observe point defects between opposite-handed domains, and we explain a preference for point defects over domain walls.

Fluctuation spectra in polymer nematics and Frank elastic constants: a coarse-grained modelling study

Monte Carlo simulations of uniaxial nematic polymer melts are performed, based on a discrete worm-like chain model combined with soft, anisotropic non-bonded potentials. Different chain lengths are considered, the contour length of the longest being an order of magnitude larger than the persistence length. From equilibrated melt configurations, density and director fluctuation spectra are calculated and compared with analytical predictions available in literature. The latter typically correspond to hydrodynamic treatments of correlations and assume that there is no chain backfolding along the nematic director. Nevertheless, it is demonstrated that the analytical theories capture several features of the spectra obtained in the current simulations, where moderate backfolding of polymer chains is observed. Based on the available analytical expressions for density and director fluctuation spectra, material properties, such as Frank elastic constants, are extracted. Their dependence on polymerisation degree is studied and found to reproduce theoretically expected trends. For example, evidence is provided that the splay constant increases linearly with chain length, when effects of hairpins are negligible.

Basis function sampling: a new paradigm for material property computation

Wang-Landau sampling, and the associated class of flat histogram simulation methods have been remarkably helpful for calculations of the free energy in a wide variety of physical systems. Practically, convergence of these calculations to a target free energy surface is hampered by reliance on parameters which are unknown a priori. Here, we derive and implement a method built upon orthogonal functions which is fast, parameter-free, and (importantly) geometrically robust. The method is shown to be highly effective in achieving convergence. An important feature of this method is its ability to attain arbitrary levels of description for the free energy. It is thus ideally suited to in silico measurement of elastic moduli and other material properties related to free energy perturbations. We demonstrate the utility of such applications by applying our method to calculate the Frank elastic constants of the Lebwohl-Lasher model of liquid crystals.

Reversible Switching of Liquid Crystalline Order Permits Synthesis of Homogeneous Populations of Dipolar Patchy Microparticles

The spontaneous positioning of colloids on the surfaces of micrometer-sized liquid crystalline droplets and their subsequent polymerization offers the basis of a general and facile method for the synthesis of patchy microparticles. The existence of multiple local energetic minima, however, can generate kinetic traps for colloids on the surfaces of the liquid crystal (LC) droplets and result in heterogeneous populations of patchy microparticles. To address this issue, here we demonstrate that adsorbate-driven switching of the internal configurations of LC droplets can be used to sweep colloids to a single location on the LC droplet surfaces, thus resulting in the synthesis of homogeneous populations of patchy microparticles. The surface-driven switching of the LC can be triggered by addition of surfactant or salts, and permits the synthesis of dipolar microparticles as well as "Janus-like" microparticles. By using magnetic colloids, we illustrate the utility of the approach by synthesizing magnetically-responsive patchy microdroplets of LC with either dipolar or quadrupolar symmetry that exhibit distinct optical responses upon application of an external magnetic field.