Perturbative density functional methods for cholesteric liquid crystals

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.

Chiral shape fluctuations and the origin of chirality in cholesteric phases of DNA origamis

Lyotropic cholesteric liquid crystal phases are ubiquitously observed in biological and synthetic polymer solutions, characterized by a complex interplay between thermal fluctuations and entropic and enthalpic forces. The elucidation of the link between microscopic features and macroscopic chiral structure, and of the relative roles of these competing contributions on phase organization, remains a topical issue. Here, we provide theoretical evidence of a previously unidentified mechanism of chirality amplification in lyotropic liquid crystals, whereby phase chirality is governed by fluctuation-stabilized helical deformations in the conformations of their constituent molecules. Our results compare favorably to recent experimental studies of DNA origami assemblies and demonstrate the influence of intramolecular mechanics on chiral supramolecular order, with potential implications for a broad class of experimentally relevant colloidal systems.

Morphological analysis of chiral rod clusters from a coarse-grained single-site chiral potential

We present a coarse-grained single-site potential for simulating chiral interactions, with adjustable strength, handedness, and preferred twist angle. As an application, we perform basin-hopping global optimisation to predict the favoured geometries for clusters of chiral rods. The morphology phase diagram based upon these predictions has four distinct families, including previously reported structures for potentials that introduce chirality based on shape, such as membranes and helices. The transition between these two configurations reproduces some key features of experimental results for fd bacteriophage. The potential is computationally inexpensive, intuitive, and versatile; we expect it will be useful for large scale simulations of chiral molecules. For chiral particles confined in a cylindrical container we reproduce the behaviour observed for fusilli pasta in a jar. Hence this chiropole potential has the capability to provide insight into structures on both macroscopic and molecular length scales.

Manipulation of cholesteric liquid crystal phase behavior and molecular assembly by molecular chirality

Molecular simulation is used to study the effect of molecular chirality on liquid crystalline phase transition and molecular assembly behavior. Based on a flexible chain (FCh) model with helical arrangement of side beads, the phase behavior of FCh models with various molecular chiralities are studied as functions of pressure (or density). By modifying the molecular chirality of FCh, we can manipulate the relative stability of the nematic and cholesteric phases continuously; and we found that increasing molecular chirality may destabilize cholesteric order due to the effective reduction of chiral interactions. A semismectic phase is identified in the high-density region, in which the two-dimensional fluid layers overlap due to shift alignment formed by FCh particles. The global phase diagram of the FCh model is constructed and the potential energy surface is calculated to elucidate the formation of cholesteric phase in terms of two-body interactions.

Polymorphism of stable collagen fibrils

Collagen fibrils are versatile self-assembled structures that provide mechanical integrity within mammalian tissues. The radius of collagen fibrils vary widely depending on experimental conditions in vitro or anatomical location in vivo. Here we explore the variety of thermodynamically stable fibril configurations that are available. We use a liquid crystal model of radial collagen fibril structure with a double-twist director field. Using a numerical relaxation method we show that two dimensionless parameters, the ratio of saddle-splay to twist elastic constants k24/K22 and the ratio of surface tension to chiral strength [small gamma, Greek, tilde] ≡ γ/(K22q), largely specify both the scaled fibril radius and the associated surface twist of equilibrium fibrils. We find that collagen fibrils are the stable phase with respect to the cholesteric phase only when the reduced surface tension is small, [small gamma, Greek, tilde] ⪅ 0.2. Within this stable regime, collagen fibrils can access a wide range of radii and associated surface twists. Remarkably, we find a maximal equilibrium surface twist of 0.33 rad (19°). Our results are compatible with corneal collagen fibrils, and we show how the large surface twist can explain the narrow distribution of corneal fibril radii. Conversely, we show how small surface twist is required for the thermodynamic stability of tendon fibrils in the face of considerable polydispersity of radius.

Hierarchical bounding structures for efficient virial computations: Towards a realistic molecular description of cholesterics

We detail the application of bounding volume hierarchies to accelerate second-virial evaluations for arbitrary complex particles interacting through hard and soft finite-range potentials. This procedure, based on the construction of neighbour lists through the combined use of recursive atom-decomposition techniques and binary overlap search schemes, is shown to scale sub-logarithmically with particle resolution in the case of molecular systems with high aspect ratios. Its implementation within an efficient numerical and theoretical framework based on classical density functional theory enables us to investigate the cholesteric self-assembly of a wide range of experimentally relevant particle models. We illustrate the method through the determination of the cholesteric behavior of hard, structurally resolved twisted cuboids, and report quantitative evidence of the long-predicted phase handedness inversion with increasing particle thread angles near the phenomenological threshold value of 45°. Our results further highlight the complex relationship between microscopic structure and helical twisting power in such model systems, which may be attributed to subtle geometric variations of their chiral excluded-volume manifold.