Dynamics of defects around anisotropic particles in nematic liquid crystals under shear

Lock-key microfluidics: simulating nematic colloid advection along wavy-walled channels

Liquid crystalline media mediate interactions between suspended particles and confining geometries, which not only has potential to guide patterning and bottom-up colloidal assembly, but can also control colloidal migration in microfluidic devices. However, simulating such dynamics is challenging because nemato-elasticity, diffusivity and hydrodynamic interactions must all be accounted for within complex boundaries. We model the advection of colloids dispersed in flowing and fluctuating nematic fluids confined within 2D wavy channels. A lock-key mechanism between homeotropic colloids and troughs is found to be stronger for planar anchoring on the wavy walls compared to homeotropic anchoring on the wavy walls due to the relative location of the colloid-associated defects. Sufficiently large amplitudes result in stick-slip trajectories and even permanent locking of colloids in place. These results demonstrate that wavy walls not only have potential to direct colloids to specific docking sites but also to control site-specific resting duration and intermittent elution.

Particle-based and continuum models for confined nematics in two dimensions

We use the particle-based stochastic multi-particle collision dynamics (N-MPCD) algorithm to simulate confined nematic liquid crystals in regular two-dimensional polygons such as squares, pentagons and hexagons. We consider a range of values of the nematicities, U, and simulation domain sizes, R, that canvass nano-sized polygons to micron-sized polygons. We use closure arguments to define mappings between the N-MPCD parameters and the parameters in the continuum deterministic Landau-de Gennes framework. The averaged N-MPCD configurations agree with those predicted by Landau-de Gennes theory, at least for large polygons. We study relaxation dynamics or the non-equilibrium dynamics of confined nematics in polygons, in the N-MPCD framework, and the kinetic traps bear strong resemblance to the unstable saddle points in the Landau-de Gennes framework. Finally, we study nematic defect dynamics inside the polygons in the N-MPCD framework and the finite-size effects slow down the defects and attract them to polygon vertices. Our work is a comprehensive comparison between particle-based stochastic N-MPCD methods and deterministic/continuum Landau-de Gennes methods, and this comparison is essential for new-age multiscale theories.

Defect-influenced particle advection in highly confined liquid crystal flows

We study the morphology of the Saturn ring defect and director structure around a colloidal particle with normal anchoring conditions and within the flow of the nematic host phase through a rectangular duct of comparable size to the particle. The changes in the defect structures and director profile influence the advection behaviour of the particle, which we compare to that in a simple Newtonian host phase. These effects lead to a non-monotonous dependence of the differential velocity of particle and fluid, also known as retardation ratio, on the Ericksen number.

Active nematic liquid crystals simulated by particle-based mesoscopic methods

Two Multi-particle collision dynamics algorithms that simulate nematic liquid crystals are generalised to reproduce active behaviour. One of the algorithms is due to Shendruk and Yeomans and is based on particles that carry an orientation vector ordered by a mean-field energy [T. N. Shendruk and J. M. Yeomans, Soft Matter, 2015, 11, 5101]. In the other algorithm, due to Mandal and Mazza, particles possess an order parameter tensor which evolves according to the Qian-Sheng model of nematohydrodynamics [S. Mandal and M. G. Mazza, Phys. Rev. E, 2019, 99, 063319]. For both methods activity is incorporated through a force proportional to the divergence of the local average order parameter tensor. Both implementations produce disclination curves in the nematic fluid that undergo nucleation and self-annihilation dynamics. Topological defects are found to be consistent with those observed in recent experiments of three-dimensional active nematics. Results permit to compare the length-scales over which the different nematic Multi-particle collision dynamics methods operate. The structure and dynamics of the orientation and flow fields agree with those obtained recently in numerical studies of continuum three-dimensional active nematics. Overall, our results open the opportunity to use mesoscopic particle-based approaches to study active liquid crystals in situations such as nonequilibrium states driven by flow or colloidal particles in active anisotropic solvents.

Mesoscopic simulations of active nematics

Coarse-grained, mesoscale simulations are invaluable for studying soft condensed matter because of their ability to model systems in which a background solvent plays a substantial role but is not the primary interest. Such methods generally model passive solvents; however, far-from-equilibrium systems may also be composed of complex solutes suspended in an active fluid. Yet, few coarse-grained simulation methods exist to model an active medium. We introduce an algorithm to simulate active nematics, which builds on multiparticle collision dynamics (MPCD) for passive fluctuating nematohydrodynamics by introducing dipolar activity in the local collision operator. Active nematic MPCD (AN-MPCD) simulations not only exhibit the key characteristics of active nematic turbulence but, as a particle-based algorithm, also reproduce crucial attributes of active particle models. Thus, mesoscopic AN-MPCD is an approach that bridges microscopic and continuum descriptions, allowing simulations of composite active-passive systems.

Interactions of charged microrods in chiral nematic liquid crystals

We study the pair interaction of charged silica microrods in chiral nematic liquid crystals and show that the microrods with homeotropic surface anchoring form a bound state due to the competing effect of electrostatic (Coulomb) and elastic interactions. The robustness of the bound state is demonstrated by applying external electrical and mechanical forces that perturbs their equilibrium position as well as orientation. In the bound state we have measured the correlated thermal fluctuations of the position, using two-particle cross-correlation spectroscopy that uncovers their hydrodynamic interaction. These findings reveal unexplored aspects of liquid-crystal dispersions which are important for understanding the assembly and dynamics of nano- and microparticles in chiral nematic liquid crystals.

Nonequilibrium Dynamics of a Magnetic Nanocapsule in a Nematic Liquid Crystal

Colloidal particles in nematic liquid crystals show a beautiful variety of complex phenomena with promising applications. Their dynamical behaviour is determined by topology and interactions with the liquid crystal and external fields. Here, a nematic magnetic nanocapsule reoriented periodically by time-varying magnetic fields is studied using numerical simulations. The approach combines Molecular Dynamics to resolve solute–solvent interactions and Nematic Multiparticle Collision Dynamics to incorporate nematohydrodynamic fields and fluctuations. A Saturn ring defect resulting from homeotropic anchoring conditions surrounds the capsule and rotates together with it. Magnetically induced rotations of the capsule can produce transformations of this topological defect, which changes from a disclination curve to a defect structure extending over the surface of the capsule. Transformations occur for large magnetic fields. At moderate fields, elastic torques prevent changes of the topological defect by tilting the capsule out from the rotation plane of the magnetic field.