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.

A novel model for smectic liquid crystals: Elastic anisotropy and response to a steady-state flow.

By means of a combination of equilibrium Monte Carlo and molecular dynamics simulations and nonequilibrium molecular dynamics we investigate the ordered, uniaxial phases (i.e., nematic and smectic A) of a model liquid crystal. We characterize equilibrium behavior through their diffusive behavior and elastic properties. As one approaches the equilibrium isotropic-nematic phase transition, diffusion becomes anisotropic in that self-diffusion D⊥ in the direction orthogonal to a molecule's long axis is more hindered than self-diffusion D∥ in the direction parallel to that axis. Close to nematic-smectic A phase transition the opposite is true, D∥ < D⊥. The Frank elastic constants K1, K2, and K3 for the respective splay, twist, and bend deformations of the director field n̂ are no longer equal and exhibit a temperature dependence observed experimentally for cyanobiphenyls. Under nonequilibrium conditions, a pressure gradient applied to the smectic A phase generates Poiseuille-like or plug flow depending on whether the convective velocity is parallel or orthogonal to the plane of smectic layers. We find that in Poiseuille-like flow the viscosity of the smectic A phase is higher than in plug flow. This can be rationalized via the velocity-field component in the direction of the flow. In a sufficiently strong flow these smectic layers are not destroyed but significantly bent.

Coarse-grained treatment of the self-assembly of colloids suspended in a nematic host phase.

The complex interplay of molecular scale effects, nonlinearities in the orientational field and long-range elastic forces makes liquid-crystal physics very challenging. A consistent way to extract information from the microscopic, molecular scale up to the meso- and macroscopic scale is still missing. Here, we develop a hybrid procedure that bridges this gap by combining extensive Monte Carlo (MC) simulations, a local Landau-de Gennes theory, classical density functional theory, and finite-size scaling theory. As a test case to demonstrate the power and validity of our novel approach we study the effective interaction among colloids with Boojum defect topology immersed in a nematic liquid crystal. In particular, at sufficiently small separations colloids attract each other if the angle between their center-of-mass distance vector and the far-field nematic director is about 30°. Using the effective potential in coarse-grained two-dimensional MC simulations we show that self-assembled structures formed by the colloids are in excellent agreement with experimental data.

Defect topologies in chiral liquid crystals confined to mesoscopic channels.

We present Monte Carlo simulations in the grand canonical and canonical ensembles of a chiral liquid crystal confined to mesochannels of variable sizes and geometries. The mesochannels are taken to be quasi-infinite in one dimension but finite in the two other directions. Under thermodynamic conditions chosen and for a selected value of the chirality coupling constant, the bulk liquid crystal exhibits structural characteristics of a blue phase II. This is established through the tetrahedral symmetry of disclination lines and the characteristic simple-cubic arrangement of double-twist helices formed by the liquid-crystal molecules along all three axes of a Cartesian coordinate system. If the blue phase II is then exposed to confinement, the interplay between its helical structure, various anchoring conditions at the walls of the mesochannels, and the shape of the mesochannels gives rise to a broad variety of novel, qualitative disclination-line structures that are reported here for the first time.

Disclination lines at homogeneous and heterogeneous colloids immersed in a chiral liquid crystal.

In the present work we perform Monte Carlo simulations in the isothermal-isobaric ensemble to study defect topologies formed in a cholesteric liquid crystal due to the presence of a spherical colloidal particle. Topological defects arise because of the competition between anchoring at the colloidal surface and the local director. We consider homogeneous colloids with either local homeotropic or planar anchoring to validate our model by comparison with earlier lattice Boltzmann studies. Furthermore, we perform simulations of a colloid in a twisted nematic cell and discuss the difference between induced and intrinsic chirality on the formation of topological defects. We present a simple geometrical argument capable of describing the complex three-dimensional topology of disclination lines evolving near the surface of the colloid. The presence of a Janus colloid in a cholesteric host fluid reveals a rich variety of defect structures. Using the Frank free energy we analyze these defects quantitatively indicating a preferred orientation of the Janus colloid relative to the cholesteric helix.

Defect topologies in a nematic liquid crystal near a patchy colloid.

Using isothermal-isobaric Monte Carlo simulations we investigate defect topologies due to a spherical colloidal particle immersed in a nematic liquid crystal. Defects arise because of the competition between the preferential orientation at the colloid's surface and the far-field director ̂n(0). Considering a chemically homogeneous colloid as a special case we observe the well-known surface and saturn ring defect topologies for weak and strong perpendicular anchoring, respectively; for homogeneous, strong parallel anchoring we find a boojum defect topology that has been seen experimentally [see P. Poulin and D. A. Weitz, Phys. Rev. E 57, 626 (1998)] but not in computer simulations. We also consider a heterogeneous, patchy colloid where the liquid-crystal molecules anchor either preferentially planar or perpendicular at the surface of the colloid. For a patchy colloid we observe a boojum ring defect topology in agreement with recent experimental studies [see M. Conradi, M. Ravnik, M. Bele, M. Zorko, S. Zumer, and I. Musevic, Soft Matter 5, 3905 (2009)]. We also observe two other novel defect topologies that have not been reported thus far neither experimentally nor theoretically.

Finite-size scaling analysis of isotropic-polar phase transitions in an amphiphilic fluid.

We present Monte Carlo simulations of the isotropic-polar (IP) phase transition in an amphiphilic fluid carried out in the isothermal-isobaric ensemble. Our model consists of Lennard-Jones spheres where the attractive part of the potential is modified by an orientation-dependent function. This function gives rise to an angle dependence of the intermolecular attractions corresponding to that characteristic of point dipoles. Our data show a substantial system-size dependence of the dipolar order parameter. We analyze the system-size dependence in terms of the order-parameter distribution and a cumulant involving its first and second moments. The order parameter, its distribution, and susceptibility observe the scaling behavior characteristic of the 3D Ising universality class. Because of this scaling behavior and because all cumulants have a common intersection irrespective of system size we conclude that the IP phase transition is continuous. Considering pressures 1.3 ≤ P ≤ 3.0 we demonstrate that a line of continuous phase transitions exists which is analogous to the Curie line in systems exhibiting a ferroelectric transition. Our results are qualitatively consistent with Landau's theory of continuous phase transitions.