Cholesteric and screw-like nematic phases in systems of helical particles

Stability of the chiral crystal phase and breakdown of the cholesteric phase in mixtures of active-passive chiral rods

In this study, we aim to explore the effect of chirality on the phase behavior of active helical particles driven by two-temperature scalar activity. We first calculate the equation of state of soft helical particles of various intrinsic chiralities using molecular dynamics (MD) simulation. In equilibrium, the emergence of various liquid crystal (LC) phases such as nematic (N), cholesteric , smectic (Sm) and crystal (K) crucially depends on the presence of walls that induce planar alignment. Next, we introduce activity through the two-temperature model: keep increasing the temperature of half of the helical particles (labeled as 'hot' particles) while maintaining the temperature of the other half at a lower value (labeled as 'cold' particles). Starting from a homogeneous isotropic (I) phase, we find the emergence of 2-TIPS: two temperature-induced phase separations between the hot and cold particles. We also observe that the cold particles undergo an ordering transition to various LC phases even in the absence of a wall. This observation reveals that the hot-cold interface in the active system plays the role of a wall in the equilibrium system by inducing an alignment direction for the cold particles. However, in the case of a cholesteric phase, we observe that activity destabilizes the phase by inducing smectic ordering in the cold zone while an isotropic structure in the hot zone. The smectic ordering in the cold zone eventually transforms to a chiral crystal phase with high enough activity.

Janus helices: From fully attractive to hard helices

The phase diagram of hard helices differs from its hard rods counterpart by the presence of chiral "screw" phases stemming from the characteristic helical shape, in addition to the conventional liquid crystal phases also found for rod-like particles. Using extensive Monte Carlo and Molecular Dynamics simulations, we study the effect of the addition of a short-range attractive tail representing solvent-induced interactions to a fraction of the sites forming the hard helices, ranging from a single-site attraction to fully attractive helices for a specific helical shape. Different temperature regimes exist for different fractions of the attractive sites, as assessed in terms of the relative Boyle temperatures, that are found to be rather insensitive to the specific shape of the helical particle. The temperature range probed by the present study is well above the corresponding Boyle temperatures, with the phase behaviour still mainly entropically dominated and with the existence and location of the various liquid crystal phases only marginally affected. The pressure in the equation of state is found to decrease upon increasing the fraction of attractive beads and/or on lowering the temperature at fixed volume fraction, as expected on physical grounds. All screw phases are found to be stable within the considered range of temperatures with the smectic phase becoming more stable on lowering the temperature. By contrast, the location of the transition lines do not display a simple dependence on the fraction of attractive beads in the considered range of temperatures.

Conformational degrees of freedom and stability of splay-bend ordering in the limit of a very strong planar anchoring

We study the self-organization in a monolayer (a two-dimensional system) of flexible planar trimer particles. The molecules are made up of two mesogenic units linked by a spacer, all of which are modeled as hard needles of the same length. Each molecule can dynamically adopt two conformational states: an achiral bent-shaped (cis-) and a chiral zigzag (trans-) one. Using constant pressure Monte Carlo simulations and Onsager-type density functional theory (DFT), we show that the system consisting of these molecules exhibits a rich spectrum of liquid crystalline phases. The most interesting observation is the identification of stable smectic splay-bend (S_{SB}) and chiral smectic-A (S_{A}^{*}) phases. The S_{SB} phase is also stable in the limit, where only cis- conformers are allowed. The second phase that occupies a considerable portion of the phase diagram is S_{A}^{*} with chiral layers, where the chirality of the neighboring layers is of opposite sign. The study of the average fractions of the trans- and cis- conformers in various phases shows that while in the isotropic phase all fractions are equally populated, the S_{A}^{*} phase is dominated by chiral conformers (zigzag), but the achiral conformers win in the smectic splay-bend phase. To clarify the possibility of stabilization of the nematic splay-bend (N_{SB}) phase for trimers, the free energy of the N_{SB} and S_{SB} phases is calculated within DFT for the cis- conformers, for densities where simulations show stable S_{SB}. It turns out that the N_{SB} phase is unstable away from the phase transition to the nematic phase, and its free energy is always higher than that of S_{SB}, down to the transition to the nematic phase, although the difference in free energies becomes extremely small when approaching the transition.

The thermodynamic origins of chiral twist in monolayer assemblies of rod-like colloids

The propagation of chirality across scales is a common but poorly understood phenomenon in soft matter. Here, using computer simulations, we study twisted monolayer assemblies formed by both chiral and achiral rod-like particles in the presence of non-adsorbing polymer and characterise the thermodynamic driving forces responsible for the twisting. We observe assemblies with both like and inverted chirality relative to the rods and show that the preferred twist is already determined during the initial stage of the self-assembly. Depending on the geometry of the constituent rods, the chiral twist is regulated by either the entropy gain of the polymer, or of the rods, or both. This can include important contributions from changes in both the surface area and volume of the monolayer and from rod fluctuations perpendicular to the monolayer. These findings can deepen our understanding of why chirality propagates and of how to control it.

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.

A versatile simulation method for studying phase behavior and dynamics in colloidal rod and rod-polymer suspensions

Here, we present an implicit-solvent model for dynamic simulations of hard-rod and rod-polymer suspensions. Individual rods are represented by a rigid linear chain consisting of overlapping spheres which interact through a pseudohard-core potential based on the cut-and-shifted Mie (generalized Lennard-Jones) potential with exponents (50, 49). In the rod-polymer suspensions, the polymers are modeled as freely interpenetrable spheres with respect to each other, while there is the pseudohard-core repulsion between the polymer and rod spheres. Dynamic simulations with this model are carried out with a dissipative particle dynamics (DPD) thermostat-each sphere is put in a larger DPD sphere and thus interacts with others via additional pairwise frictional and random forces-which captures the effects of Brownian forces due to the solvent while conserving local momentum. The phase behavior of these models, obtained from continuous compression and expansion simulations, reproduces previous predictions based on theoretical calculations and Monte Carlo simulations. Our method is suited to study dynamic processes in these suspensions, including nucleation and self-assembly, and can be readily extended to colloidal particles of different shapes and chemistry.

Self-Assembled Chiral Photonic Crystals from a Colloidal Helix Racemate

Chiral crystals consisting of microhelices have many optical properties, while presently available fabrication processes limit their large-scale applications in photonic devices. Here, by using a simplified simulation method, we investigate a bottom-up self-assembly route to build up helical crystals from the smectic monolayer of a colloidal helix racemate. With increasing the density, the system undergoes an entropy-driven cocrystallization by forming crystals of various symmetries with different helical shapes. In particular, we identify two crystals of helices arranged in binary honeycomb and square lattices, which are essentially composed of two sets of opposite-handed chiral crystals. Photonic calculations show that these chiral structures can have large complete photonic band gaps. In addition, in the self-assembled chiral square crystal, we also find dual polarization band gaps that selectively forbid the propagation of circularly polarized light of a specific handedness along the helical axis direction. The self-assembly process in our proposed system is robust, suggesting possibilities of using chiral colloids to assemble photonic metamaterials.