Hydrodynamics Defines the Stable Swimming Direction of Spherical Squirmers in a Nematic Liquid Crystal

Microswimmers Knead Nematics into Cholesterics

The hydrodynamic stresses created by active particles can destabilize orientational order present in the system. This is manifested, for example, by the appearance of a bend instability in active nematics or in quasi-two-dimensional living liquid crystals consisting of swimming bacteria in thin nematic films. Using large-scale hydrodynamics simulations, we study a system consisting of spherical microswimmers within a three-dimensional nematic liquid crystal. We observe a spontaneous chiral symmetry breaking, where the uniform nematic state is kneaded into a continuously twisting state, corresponding to a helical director configuration akin to a cholesteric liquid crystal. The transition arises from the hydrodynamic coupling between the liquid crystalline elasticity and the swimmer flow fields, leading to a twist-bend instability of the nematic order. It is observed for both pusher (extensile) and puller (contractile) swimmers. Further, we show that the liquid crystal director and particle trajectories are connected: in the cholesteric state the particle trajectories become helicoidal.

Simulating dynamics of ellipsoidal particles using lattice Boltzmann method

Anisotropic particles are often encountered in different fields of soft matter and complex fluids. In this work, we present an implementation of the coupled hydrodynamics of solid ellipsoidal particles and the surrounding fluid using the lattice Boltzmann method. A standard link-based mechanism is used to implement the solid-fluid boundary conditions. We develop an implicit method to update the position and orientation of the ellipsoid. This exploits the relations between the quaternion which describes the ellipsoid's orientation and the ellipsoid's angular velocity to obtain a stable and robust dynamic update. The proposed algorithm is validated by looking at four scenarios: (i) the steady translational velocity of a spheroid subject to an external force in different orientations, (ii) the drift of an inclined spheroid subject to an imposed force, (iii) three-dimensional rotational motions in a simple shear flow (Jeffrey's orbits), and (iv) developed fluid flows and self-propulsion exhibited by a spheroidal microswimmer. In all cases the comparison of numerical results shows good agreement with known analytical solutions, irrespective of the choice of the fluid properties, geometrical parameters, and lattice Boltzmann model, thus demonstrating the robustness of the proposed algorithm.

Direct numerical simulations of a microswimmer in a viscoelastic fluid

This study presents the application of the smoothed profile (SP) method to perform direct numerical simulations for the motion of both passive and active "squirming" particles in Newtonian and viscoelastic fluids. We found that fluid elasticity has a significant impact on both the transient behavior and the steady-state velocity of the particles. Specifically, we observe that the swirling flow generated by the squirmer's surface velocity significantly enhances their swimming speed as the Weissenberg number increases, regardless of the swimming type. Furthermore, we find that pushers outperform pullers in Oldroyd-B fluids, suggesting that the speed of a squirmer depends on the swimmer type. To understand the physical origin of the phenomenon of swirling flow enhancing the swimming speed, we investigate the velocity field and polymer conformation around non-swirling and swirling neutral squirmers in viscoelastic fluids. Our investigation reveals that the velocity field around the neutral swirling squirmers exhibits pusher-like extensional flow characteristics, as well as an asymmetric polymer conformation distribution, which gives rise to this increased propulsion. This is confirmed by the investigation of the force on a fixed squirmer, which revealed that the polymer stress, particularly its diagonal components, plays a critical role in enhancing the swimming speed of swirling squirmers in viscoelastic fluids. Additionally, our results demonstrate that the maximum swimming speeds of swirling squirmers occur at an intermediate value of the fluid viscosity ratio for all swimmer types. These findings have important implications for understanding the behavior of particles and micro-organisms in complex fluids.

Symbiotic dynamics in living liquid crystals

An amalgam of nematic liquid crystals and active matter, referred to as living liquid crystals, is a promising self-healing material with futuristic applications for targeted delivery of information and microcargo. We provide a phenomenological model to study the symbiotic pattern dynamics in this contemporary system using the Toner-Tu model for active matter (AM), the Landau-de Gennes free energy for liquid crystals (LCs), and an experimentally motivated coupling term that favours coalignment of the active and nematic components. Our extensive theoretical studies unfold two novel steady states, chimeras and solitons, with sharp regions of distinct orientational order that sweep through the coupled system in synchrony. The induced dynamics in the passive nematic is unprecedented. We show that the symbiotic dynamics of the AM and LC components can be exploited to induce and manipulate order in an otherwise disordered system.

Spontaneous vortex formation by microswimmers with retarded attractions

Collective states of inanimate particles self-assemble through physical interactions and thermal motion. Despite some phenomenological resemblance, including signatures of criticality, the autonomous dynamics that binds motile agents into flocks, herds, or swarms allows for much richer behavior. Low-dimensional models have hinted at the crucial role played in this respect by perceived information, decision-making, and feedback, implying that the corresponding interactions are inevitably retarded. Here we present experiments on spherical Brownian microswimmers with delayed self-propulsion toward a spatially fixed target. We observe a spontaneous symmetry breaking to a transiently chiral dynamical state and concomitant critical behavior that do not rely on many-particle cooperativity. By comparison with the stochastic delay differential equation of motion of a single swimmer, we pinpoint the delay-induced effective synchronization of the swimmers with their own past as the key mechanism. Increasing numbers of swimmers self-organize into layers with pro- and retrograde orbital motion, synchronized and stabilized by steric, phoretic, and hydrodynamic interactions. Our results demonstrate how even most simple retarded interactions can foster emergent complex adaptive behavior in small active-particle ensembles.

Anchoring-induced nonmonotonic velocity versus temperature dependence of motile bacteria in a lyotropic nematic liquid crystal

The elastic and viscous properties of lyotropic chromonic liquid crystals have a very sharp, often exponential temperature dependence. Self-propelled bacteria swimming in this viscoelastic medium induce director deformations which can strongly influence their velocity, and we study the temperature behavior of their motility in the whole range of the nematic phase. We observe experimentally that, with increasing temperature, while the viscosity drops exponentially and the frequency of the flagellum rotation grows linearly, the swimmers' speed first conventionally increases but then, above some crossover temperature, slows down and at the same time bacteria-induced director distortions become visible. It is shown that the physics behind this temperature-driven effect is in a sharp rise in the ability of the bacterium's flagellum to induce director deformations. As temperature increases, the splay and bend elastic constants sharply decrease and the anchoring extrapolation length of the flagellum surface gets shorter and shorter. At the crossover temperature the resulting effective anchoring effect dominates the fast dropping viscosity and the distortion strengthens. As a result, a fraction of the torque the flagellum applies for the propulsion is spent for the elastic degrees of freedom, which results in a bacterium slowdown. To find the director distortions, the flagellum is presented as a collection of anchoring-induced elastic monopoles, and the bacterium velocity is found from the balance of the energy spent for the propulsion and the viscous drag and nematodynamic dissipation.

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. The patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

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.

Towards an analytical description of active microswimmers in clean and in surfactant-covered drops

Geometric confinements are frequently encountered in the biological world and strongly affect the stability, topology, and transport properties of active suspensions in viscous flow. Based on a far-field analytical model, the low-Reynolds-number locomotion of a self-propelled microswimmer moving inside a clean viscous drop or a drop covered with a homogeneously distributed surfactant, is theoretically examined. The interfacial viscous stresses induced by the surfactant are described by the well-established Boussinesq-Scriven constitutive rheological model. Moreover, the active agent is represented by a force dipole and the resulting fluid-mediated hydrodynamic couplings between the swimmer and the confining drop are investigated. We find that the presence of the surfactant significantly alters the dynamics of the encapsulated swimmer by enhancing its reorientation. Exact solutions for the velocity images for the Stokeslet and dipolar flow singularities inside the drop are introduced and expressed in terms of infinite series of harmonic components. Our results offer useful insights into guiding principles for the control of confined active matter systems and support the objective of utilizing synthetic microswimmers to drive drops for targeted drug delivery applications.

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.

Anomalous Diffusion of Motile Colloids Dispersed in Liquid Crystals

We study the superdiffusion of driven colloidal particles dispersed in a nematic liquid crystal. While motion is ballistic in the driving direction, our experiments show that transversal fluctuations become superdiffusive depending on the topological defect pattern around the inclusions. The phenomenon can be reproduced with different driving methods and propulsion speeds, while it is strongly dependent on particle size and temperature. We propose a mechanism based on the geometry of the liquid crystal backflow around the inclusions to justify the persistence of thermal fluctuations and to explain the observed temperature and particle size dependence of the superdiffusive behavior based on material and geometrical parameters.

Long-time anomalous swimmer diffusion in smectic liquid crystals

The dynamics of self-locomotion of active particles in aligned or liquid crystalline fluids strongly deviates from that in simple isotropic media. We explore the long-time dynamics of a swimmer moving in a three-dimensional smectic liquid crystal and find that the mean-square displacement transverse to the director exhibits a distinct logarithmic tail at long times. The scaling is distinctly different from that in an isotropic or nematic fluid and hints at the subtle but important role of the director fluctuation spectrum in governing the long-time motility of active particles. Our findings are based on a generic hydrodynamic theory and Brownian dynamics computer simulation of a three-dimensional soft mesogen model.

Simulating squirmers with multiparticle collision dynamics

Multiparticle collision dynamics is a modern coarse-grained simulation technique to treat the hydrodynamics of Newtonian fluids by solving the Navier-Stokes equations. Naturally, it also includes thermal noise. Initially it has been applied extensively to spherical colloids or bead-spring polymers immersed in a fluid. Here, we review and discuss the use of multiparticle collision dynamics for studying the motion of spherical model microswimmers called squirmers moving in viscous fluids.

Hydrodynamic interaction of a self-propelling particle with a wall : Comparison between an active Janus particle and a squirmer model

Using lattice Boltzmann simulations we study the hydrodynamics of an active spherical particle near a no-slip wall. We develop a computational model for an active Janus particle, by considering different and independent mobilities on the two hemispheres and compare the behaviour to a standard squirmer model. We show that the topology of the far-field hydrodynamic nature of the active Janus particle is similar to the standard squirmer model, but in the near-field the hydrodynamics differ. In order to study how the near-field effects affect the interaction between the particle and a flat wall, we compare the behaviour of a Janus swimmer and a squirmer near a no-slip surface via extensive numerical simulations. Our results show generally a good agreement between these two models, but they reveal some key differences especially with low magnitudes of the squirming parameter [Formula: see text]. Notably the affinity of the particles to be trapped at a surface is increased for the active Janus particles when compared to standard squirmers. Finally, we find that when the particle is trapped on the surface, the velocity parallel to the surface exceeds the bulk swimming speed and scales linearly with [Formula: see text].