Calculating the helical twisting power of dopants in a liquid crystal by computer simulation

Pre-Smectic Ordering and the Unwinding Helix in Monte Carlo Simulations of Cholesteric Liquid-Crystals

In this paper, molecular chirality is studied for liquid-crystal fluids represented by hard rods with the addition of an attractive chiral dispersion term. Chiral forces between molecular pairs are assumed to be long-ranged and are described in terms of the pseudotensor of Goossens [W. J. A. Goossens, Mol. Cryst. Liq. Cryst. 1971, 12, 237-244]. Following Varga and Jackson [S. Varga and G. Jackson, Chem. Phys. Lett. 2003, 377, 6-12], this is combined with a hard-spherocylinder core. We investigate the relationship between molecular chirality and the helical pitch of the system, which occurs in the absence of full three-dimensional periodic boundary conditions. The dependence of the wavenumber of this pitch on the thermodynamic variables, temperature, and density is measured. We also explore the use of a novel surface boundary interaction model. As a result of this approach, we are able to lower the temperature of the system without the occurrence of nematic droplets, which would interfere with the formation of a uniaxial pitch. Regarding the theoretical predictions of Wensink and Jackson [H. H. Wensink and G. Jackson, J. Chem. Phys. 2009, 130, 234911], on the one hand, we have qualitative agreement with the observed non-monotonic density dependence of the wavenumber. Initially increasing with density, the wavenumber reaches a maximum, before falling as the density moves toward the point of phase transition from cholesteric to smectic. However, further analysis for shorter rods, in the presence of novel boundary conditions, reveals some disagreement with the theory, at least in this case; the unwinding of the cholesteric helix in the cholesteric phase occurs simultaneously with subtle increases in smectic ordering. These pre-smectic fluctuations have not been accounted for so far in theories on cholesterics but turn out to play a key role in controlling the pitch of cholesteric phases of rod-shaped mesogens with a small to moderate aspect ratio.

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.

Simulating the pitch sensitivity of twisted nematics of patchy rods

Stiff, elongated biomolecules such as filamentous viruses, DNA or cellulose nanocrystals are known to form liquid crystals often exhibiting a helical supramolecular organization. Little is known about the microscopic origin, size and handedness of the helical pitch in these, so-called cholesteric phases. Experimental observations in chiral lyotropics suggest that long-ranged chiral forces of electrostatic origin acting between the mesogens are responsible for such organization. Using large-scale computer simulation we study the sensitivity of the pitch imparted by soft microscopic helices and confirm that the helical sense is sensitive to a change of packing fraction, magnitude of the molecular pitch and amplitude of the chiral interactions. In particular, we find evidence that the cholesteric helix sense may change spontaneously upon variation of particle density, at fixed molecular chirality. These pitch inversions have been reported in recent theoretical studies but simulation evidence remains elusive. We rationalize these sudden changes in the supramolecular helical symmetry on the basis of detailed measurements of the mean-torque generated by the twisting of the helices. The simulation methodology employed does not require confining the twisted nematic in a slab geometry and allows for a simultaneous measurement of the pitch and the twist elastic constant. We find that the twist elastic constant increases almost linearly with density suggesting that twisted nematic shows no signs of anomalous stiffening due to pre-smectic fluctuations at higher packing fraction.

Chiral amplification in a cyanobiphenyl nematic liquid crystal doped with helicene-like derivatives

The addition of a chiral non-racemic dopant to a nematic liquid crystal (LC) has the effect of transferring the molecular chirality to the phase organization and a chiral nematic phase is formed. This molecular chirality amplification in the LC provides a unique possibility for investigating the relationship between molecular structure, intermolecular interactions, and mesoscale organization. It is known that axially chiral or helical-shaped molecules with reduced conformational disorder are good candidates for high helical twisting power derivatives. In particular, biaryl derivatives are known to be efficient chiral inducers in biaryl nematic mesophases. In this paper, we focus on a new series of helicene-like molecules of known absolute configuration. We have integrated cholesteric pitch measurements with geometry optimization by DFT calculations and analysis of the twisting ability by the Surface Chirality model to shed light on the structural features responsible for the analogies and differences exhibited by these derivatives. The investigation of these dopants with well-defined geometry, by virtue of the low conformational freedom, and the substituents variously distributed around the core, allows us to extend our knowledge of the molecular origin of the chirality amplification in liquid crystals and to confirm the simple relationship “molecular P-helicity” → “cholesteric P-handedness” for helical-shaped helicene-like derivatives.

Calculations of helical twisting powers from intermolecular torques

We present a Monte Carlo molecular simulation method that calculates the helical twisting power of a chiral molecule by sampling intermolecular torques. The approach is applied to an achiral nematic liquid crystalline system, composed of Gay-Berne particles, that is doped with chiral molecules. Calculations are presented for six chiral dopant molecules and the results show a good correlation with the sign and magnitude of experimentally determined helical twisting powers.

Simulations using hard particles

A short review is given of recent progress in the computer simulation of liquid crystal phases using hard particles. Emphasis is placed on the richness of phase behaviour that may result from the effects of molecular size and shape alone, and on the role of simulations in testing modern theories of liquid crystal phase transitions, structure and dynamics. Two specific examples are treated in detail: the simulation of twisted nematic liquid crystals, allowing a direct calculation of the twist elastic constant and the helical twisting power of chiral dopant molecules; and the recent quantitative explanation of diffusive behaviour in isotropic and nematic liquids using kinetic theory.