Stochastic rotation dynamics for nematic liquid crystals

Entangled nematic disclinations using multi-particle collision dynamics

Colloids dispersed in nematic liquid crystals form topological composites in which colloid-associated defects mediate interactions while adhering to fundamental topological constraints. Better realising the promise of such materials requires numerical methods that model nematic inclusions in dynamic and complex scenarios. We employ a mesoscale approach for simulating colloids as mobile surfaces embedded in a fluctuating nematohydrodynamic medium to study the kinetics of colloidal entanglement. In addition to reproducing far-field interactions, topological properties of disclination loops are resolved to reveal their metastable states and topological transitions during relaxation towards ground state. The intrinsic hydrodynamic fluctuations distinguish formerly unexplored far-from-equilibrium disclination states, including configurations with localised positive winding profiles. The adaptability and precision of this numerical approach offers promising avenues for studying the dynamics of colloids and topological defects in designed and out-of-equilibrium situations.

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.

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.

Multiparticle Collision Dynamics for Ferrofluids

Detailed studies of the intriguing field-dependent dynamics and transport properties of confined flowing ferrofluids require efficient mesoscopic simulation methods that account for fluctuating ferrohydrodynamics. Here, we propose such a new mesoscopic model for the dynamics and flow of ferrofluids, where we couple the multi-particle collision dynamics method as a solver for the fluctuating hydrodynamics equations to the stochastic magnetization dynamics of suspended magnetic nanoparticles. This hybrid model is validated by reproducing the magnetoviscous effect in Poiseuille flow, obtaining the rotational viscosity in quantitative agreement with theoretical predictions. We also illustrate the new method for the benchmark problem of flow around a square cylinder. Interestingly, we observe that the length of the recirculation region is increased whereas the drag coefficient is decreased in ferrofluids when an external magnetic field is applied, compared with the field-free case at the same effective Reynolds number. The presence of thermal fluctuations and the flexibility of this particle-based mesoscopic method provides a promising tool to investigate a broad range of flow phenomena of magnetic fluids and could also serve as an efficient way to simulate solvent effects when colloidal particles are immersed in ferrofluids.

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.

Multiparticle collision dynamics simulations of a squirmer in a nematic fluid

We study the dynamics of a squirmer in a nematic liquid crystal using the multiparticle collision dynamics (MPCD) method. A recently developed nematic MPCD method [Phys. Rev. E 99, 063319 (2019)] which employs a tensor order parameter to describe the spatial and temporal variations of the nematic order is used to simulate the suspending anisotropic fluid. Considering both nematodynamic effects (anisotropic viscosity and elasticity) and thermal fluctuations, in the present study, we couple the nematic MPCD algorithm with a molecular dynamics (MD) scheme for the squirmer. A unique feature of the proposed method is that the nematic order, the fluid, and the squirmer are all represented in a particle-based framework. To test the applicability of this nematic MPCD-MD method, we simulate the dynamics of a spherical squirmer with homeotropic surface anchoring conditions in a bulk domain. The importance of anisotropic viscosity and elasticity on the squirmer's speed and orientation is studied for different values of self-propulsion strength and squirmer type (pusher, puller or neutral). In sharp contrast to Newtonian fluids, the speed of the squirmer in a nematic fluid depends on the squirmer type. Interestingly, the speed of a strong pusher is smaller in the nematic fluid than for the Newtonian case. The orientational dynamics of the squirmer in the nematic fluid also shows a non-trivial dependence on the squirmer type. Our results compare well with existing experimental and numerical data. The full particle-based framework could be easily extended to model the dynamics of multiple squirmers in anisotropic fluids.

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

Nematic multiparticle collision dynamics is used to simulate disclination ring defects around spherocylinders suspended in a liquid crystal. A solvent-solute interaction potential is integrated over a short-time scale by an auxiliary molecular dynamics procedure that updates the translational and angular coordinates of the spherocylinders. For suspended particles with length in the range ∼(60,160)nm and a fixed aspect ratio, this method is able to simulate static defects reported previously in the literature. It also simulates orientation fluctuations of the elongated colloids that exhibit a broad distribution and a slow relaxation rate. Finally, a nematic host driven from equilibrium by shear flow is simulated, and the consequent dynamic behavior of the colloid-defect pair is studied. Defects under shear present significant structural transformations from chairlike disclination rings to extended defects that cover most of the cylindrical surface of the colloid. This effect results from the hydrodynamic torque on the nematic field caused by the distorted flow around the spherocylinder, and it is present for small Reynolds and Ericksen numbers of order unity.

Multi-particle collision dynamics with a non-ideal equation of state. I

The method of multi-particle collision dynamics (MPCD) and its different implementations are commonly used in the field of soft matter physics to simulate fluid flow at the micron scale. Typically, the coarse-grained fluid particles are described by the equation of state of an ideal gas, and the fluid is rather compressible. This is in contrast to conventional fluids, which are incompressible for velocities much below the speed of sound, and can cause inhomogeneities in density. We propose an algorithm for MPCD with a modified collision rule that results in a non-ideal equation of state and a significantly decreased compressibility. It allows simulations at less computational costs compared to conventional MPCD algorithms. We derive analytic expressions for the equation of state and the corresponding compressibility as well as shear viscosity. They show overall very good agreement with simulations, where we determine the pressure by simulating a quiet bulk fluid and the shear viscosity by simulating a linear shear flow and a Poiseuille flow.

Dynamics control of an in-plane-switching liquid crystal cell using heterogeneous substrates

The control of surface anchoring strength can be achieved by using heterogeneous substrates. In contrast to conventional substrates that control the anchoring strength by using temperature or chemical processes, heterogeneous substrates provide surface anchoring to liquid crystal molecules by using a mixed composition of (1) a zero anchoring surface and (2) planar-anchoring patches. To study the dynamics of in-plane-switching liquid crystal displays (IPS-LCDs) under external fields, a new particle-based numerical algorithm is developed to simulate both nematic liquid crystals and heterogeneous surfaces. This new method allows us to create different heterogeneous surfaces easily by adopting predefined distributions of numerical particles. The generated effective anchoring strength from the heterogeneous surface is thus calculated, and the dynamical response is found to be similar to that of conventional homogeneously processed substrates. The results suggest that the use of a heterogeneous LCD cell is a suitable alternative for creating desirable LCD substrates, for which chemical/temperature dependence can be transferred to a more controllable configurational dependence. Interestingly, we found master curves in the peak transmittance/recovery time phase space, and they appeared to be dependent solely on the cell thickness. This discovery clarifies the fundamental optical dynamics of IPS-LCD cells.

Multiparticle collision dynamics for tensorial nematodynamics

Liquid crystals establish a nearly unique combination of thermodynamic, hydrodynamic, and topological behavior. This poses a challenge to their theoretical understanding and modeling. The arena where these effects come together is the mesoscopic (micron) scale. It is then important to develop models aimed at capturing this variety of dynamics. We have generalized the particle-based multiparticle collision dynamics (MPCD) method to model the dynamics of nematic liquid crystals. Following the Qian-Sheng theory [Phys. Rev. E 58, 7475 (1998)1063-651X10.1103/PhysRevE.58.7475] of nematics, the spatial and temporal variations of the nematic director field and order parameter are described by a tensor order parameter. The key idea is to assign tensorial degrees of freedom to each MPCD particle, whose mesoscopic average is the tensor order parameter. This nematic MPCD method includes backflow effect, velocity-orientation coupling, and thermal fluctuations. We validate the applicability of this method by testing (i) the nematic-isotropic phase transition, (ii) the flow alignment of the director in shear and Poiseuille flows, and (iii) the annihilation dynamics of a pair of line defects. We find excellent agreement with existing literature. We also investigate the flow field around a force dipole in a nematic liquid crystal, which represents the leading-order flow field around a force-free microswimmer. The anisotropy of the medium not only affects the magnitude of velocity field around the force dipole, but can also induce hydrodynamic torques depending on the orientation of dipole axis relative to director field. A force dipole experiences a hydrodynamic torque when the dipole axis is tilted with respect to the far-field director. The direction of hydrodynamic torque is such that the pusher- (or puller-) type force dipole tends to orient along (or perpendicular to) the director field. Our nematic MPCD method can have far-reaching implications not only in modeling of nematic flows, but also to study the motion of colloids and microswimmers immersed in an anisotropic medium.

Electroconvection of pure nematic liquid crystals without free charge carriers

We consider electroconvection as a response of nematic liquid crystals to an external electric AC field, in the absence of free charge carriers. Previous experimental and theoretical results emphasized charge carriers as a necessary precondition of electroconvection because free-charges in the fluid can respond to an external electric field. Therefore, ionized molecules are considered as responsible for the driving of electroconvective flows. In experiments, finite conductivity is achieved by adding charge-carrying dye molecules or in non-dyed liquid crystals by impurities of the samples. The phenomenon of electroconvection is explained by the Carr-Helfrich theory, supported by numerical simulations. In the present paper, we show that electroconvection may occur also in pure nematic liquid crystals. By means of particle-based numerical simulations we found that bound charges emerge by alignment of polarized liquid crystal molecules in response to the external electric field. In our simulations we could reproduce the characteristic features of electroconvection, such as director-flow patterns, the phase-transition in the voltage-frequency diagram, and dislocation climb/glide motion, which are well known from experiments and hydrodynamic simulations under the assumption of free charge carriers.

Particle-based modeling of heterogeneous chemical kinetics including mass transfer

Connecting the macroscopic world of continuous fields to the microscopic world of discrete molecular events is important for understanding several phenomena occurring at physical boundaries of systems. An important example is heterogeneous catalysis, where reactions take place at active surfaces, but the effective reaction rates are determined by transport limitations in the bulk fluid and reaction limitations on the catalyst surface. In this work we study the macro-micro connection in a model heterogeneous catalytic reactor by means of stochastic rotation dynamics. The model is able to resolve the convective and diffusive interplay between participating species, while including adsorption, desorption, and reaction processes on the catalytic surface. Here we apply the simulation methodology to a simple straight microchannel with a catalytic strip. Dimensionless Damkohler numbers are used to comment on the spatial concentration profiles of reactants and products near the catalyst strip and in the bulk. We end the discussion with an outlook on more complicated geometries and increasingly complex reactions.

Multi-particle collision dynamics algorithm for nematic fluids

Research on transport, self-assembly and defect dynamics within confined, flowing liquid crystals requires versatile and computationally efficient mesoscopic algorithms to account for fluctuating nematohydrodynamic interactions. We present a multi-particle collision dynamics (MPCD) based algorithm to simulate liquid-crystal hydrodynamic and director fields in two and three dimensions. The nematic-MPCD method is shown to successfully reproduce the features of a nematic liquid crystal, including a nematic-isotropic phase transition with hysteresis in 3D, defect dynamics, isotropic Frank elastic coefficients, tumbling and shear alignment regimes and boundary condition-dependent order parameter fields.