Lattice Boltzmann algorithm to simulate isotropic-nematic emulsions

Topology controls flow patterns in active double emulsions

Active emulsions and liquid crystalline shells are intriguing and experimentally realisable types of topological matter. Here we numerically study the morphology and spatiotemporal dynamics of a double emulsion, where one or two passive small droplets are embedded in a larger active droplet. We find activity introduces a variety of rich and nontrivial nonequilibrium states in the system. First, a double emulsion with a single active droplet becomes self-motile, and there is a transition between translational and rotational motion: both of these regimes remain defect-free, hence topologically trivial. Second, a pair of particles nucleate one or more disclination loops, with conformational dynamics resembling a rotor or chaotic oscillator, accessed by tuning activity. In the first state a single, topologically charged, disclination loop powers the rotation. In the latter state, this disclination stretches and writhes in 3D, continuously undergoing recombination to yield an example of an active living polymer. These emulsions can be self-assembled in the lab, and provide a pathway to form flow and topological patterns in active matter in a controlled way, as opposed to bulk systems that typically yield active turbulence.

Command of three-dimensional solitary waves via photopatterning

Multidimensional solitons are prevalent in numerous research fields. In orientationally ordered soft matter system, three-dimensional director solitons exemplify the localized distortion of molecular orientation. However, their precise manipulation remains challenging due to unpredictable and uncontrolled generation. Here, we utilize preimposed programmable photopatterning in nematics to control the kinetics of director solitons. This enables both unidirectional and bidirectional generation at specific locations and times, confinement within micron-scaled patterns of diverse shapes, and directed propagation along predefined trajectories. A focused dynamical model provides insight into the origins of these solitons and aligns closely with experimental observations, underscoring the pivotal role of anchoring conditions in soliton manipulation. Our findings pave the way for diverse fundamental research avenues and promising applications, including microcargo transportation and optical information processing.

Generation and Propagation of Flexoelectricity-Induced Solitons in Nematic Liquid Crystals

Solitons in nematic liquid crystals facilitate the rapid transport and sensing in microfluidic systems. Little is known about the elementary conditions needed to create solitons in nematic materials. In this study, we apply a combination of theory, computational simulations, and experiments to examine the formation and propagation of solitary waves, or "solitons", in nematic liquid crystals under the influence of an alternating current (AC) electric field. We find that these solitary waves exhibit "butterfly"-like or "bullet"-like structures that travel in the direction perpendicular to the applied electric field. Such structures propagate over long distances without losing their initial shape. The theoretical framework adopted here helps identify several key factors leading to the formation of solitons in the absence of electrostatic interactions. These factors include surface irregularities, flexoelectric polarization, unequal elastic constants, and negative anisotropic dielectric permittivity. The results of simulations are shown to be in good agreement with our own experimental observations, serving to establish the validity of the theoretical concepts and ideas advanced in this work.

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.

Hierarchical defect-induced condensation in active nematics

Topological defects play a central role in the formation and organization of various biological systems. Historically, such nonequilibrium defects have been mainly studied in the context of homogeneous active nematics. Phase-separated systems, in turn, are known to form dense and dynamic nematic bands, but typically lack topological defects. In this paper, we use agent-based simulations of weakly aligning, self-propelled polymers and demonstrate that contrary to the existing paradigm phase-separated active nematics form -1/2 defects. Moreover, these defects, emerging due to interactions among dense nematic bands, constitute a novel second-order collective state. We investigate the morphology of defects in detail and find that their cores correspond to a strong increase in density, associated with a condensation of nematic fluxes. Unlike their analogs in homogeneous systems, such condensed defects form and decay in a different way and do not involve positively charged partners. We additionally observe and characterize lateral arc-like structures that separate from a band's bulk and move in transverse direction. We show that the key control parameters defining the route from stable bands to the coexistence of dynamic lanes and defects are the total density of particles and their path persistence length. We introduce a hydrodynamic theory that qualitatively recapitulates all the main features of the agent-based model, and use it to show that the emergence of both defects and arcs can be attributed to the same anisotropic active fluxes. Finally, we present a way to artificially engineer and position defects, and speculate about experimental verification of the provided model.

Director alignment at the nematic-isotropic interface: elastic anisotropy and active anchoring

Activity in nematics drives interfacial flows that lead to preferential alignment that is tangential or planar for extensile systems (pushers) and perpendicular or homeotropic for contractile ones (pullers). This alignment is known as active anchoring and has been reported for a number of systems and described using active nematic hydrodynamic theories. The latter are based on the one-elastic constant approximation, i.e. they assume elastic isotropy of the underlying passive nematic. Real nematics, however, have different elastic constants, which lead to interfacial anchoring. In this paper, we consider elastic anisotropy in multiphase and multicomponent hydrodynamic models of active nematics and investigate the competition between the interfacial alignment driven by the elastic anisotropy of the passive nematic and the active anchoring. We start by considering systems with translational invariance to analyse the alignment at flat interfaces and, then, consider two-dimensional systems and active nematic droplets. We investigate the competition of the two types of anchoring over a wide range of the other parameters that characterize the system. The results of the simulations reveal that the active anchoring dominates except at very low activities, when the interfaces are static. In addition, we found that the elastic anisotropy does not affect the dynamics but changes the active length that becomes anisotropic. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Motility of active nematic films driven by "active anchoring"

We provide a minimal model for an active nematic film in contact with both a solid substrate and a passive isotropic fluid, and explore its dynamics in one and two dimensions using a combination of hybrid Lattice Boltzmann simulations and analytical calculations. By imposing nematic anchoring at the substrate while active flows induce a preferred alignment at the interface ("active anchoring"), we demonstrate that directed fluid flow spontaneously emerges in cases where the two anchoring types are opposing. In one dimension, our model reduces to an analogue of a loaded elastic column. Here, the transition from a stationary to a motile state is akin to the buckling bifurcation, but offers the possibility to reverse the flow direction for a given set of parameters and boundary conditions solely by changing initial conditions. The two-dimensional variant of our model allows for additional tangential instabilities, and it is found that undulations form in the interface above a threshold activity. Our results might be relevant for the design of active microfluidic geometries or curvature-guided self-assembly.

Electric-field-induced shape transition of nematic tactoids

The occurrence of new textures of liquid crystals is an important factor in tuning their optical and photonics properties. Here, we show, both experimentally and by numerical computation, that under an electric field chitin tactoids (i.e., nematic droplets) can stretch to aspect ratios of more than 15, leading to a transition from a spindlelike to a cigarlike shape. We argue that the large extensions occur because the elastic contribution to the free energy is dominated by the anchoring. We demonstrate that the elongation involves hydrodynamic flow and is reversible: the tactoids return to their original shapes upon removing the field.

Shear dynamics of an inverted nematic emulsion

Here we study theoretically the dynamics of a 2D and a 3D isotropic droplet in a nematic liquid crystal under a shear flow. We find a large repertoire of possible nonequilibrium steady states as a function of the shear rate and of the anchoring of the nematic director field at the droplet surface. We first discuss homeotropic anchoring. For weak anchoring, we recover the typical behaviour of a sheared isotropic droplet in a binary fluid, which rotates, stretches and can be broken by the applied flow. For intermediate anchoring, new possibilities arise due to elastic effects in the nematic fluid. We find that in this regime the 2D droplet can tilt and move in the flow, or tumble incessantly at the centre of the channel. For sufficiently strong anchoring, finally, one or both of the topological defects which form close to the surface of the isotropic droplet in equilibrium detach from it and get dragged deep into the nematic state by the flow. In 3D, instead, the Saturn ring associated with the normal anchoring disclination line can be deformed and shifted downstream by the flow, but remains always localized in the proximity of the droplet, at least for the parameter range we explored. Tangential anchoring in 2D leads to a different dynamic response, as the boojum defects characteristic of this situation can unbind from the droplet under a weaker shear with respect to the normal anchoring case. Our results should stimulate further experiments with inverted liquid crystal emulsions under shear, as most of the predictions can be testable in principle by monitoring the evolution of liquid crystalline orientation patterns or by tracking the position and shape of the droplet over time.

Lattice Boltzmann simulation of asymmetric flow in nematic liquid crystals with finite anchoring

Liquid crystals (LCs) display many of the flow characteristics of liquids but exhibit long range orientational order. In the nematic phase, the coupling of structure and flow leads to complex hydrodynamic effects that remain to be fully elucidated. Here, we consider the hydrodynamics of a nematic LC in a hybrid cell, where opposite walls have conflicting anchoring boundary conditions, and we employ a 3D lattice Boltzmann method to simulate the time-dependent flow patterns that can arise. Due to the symmetry breaking of the director field within the hybrid cell, we observe that at low to moderate shear rates, the volumetric flow rate under Couette and Poiseuille flows is different for opposite flow directions. At high shear rates, the director field may undergo a topological transition which leads to symmetric flows. By applying an oscillatory pressure gradient to the channel, a net volumetric flow rate is found to depend on the magnitude and frequency of the oscillation, as well as the anchoring strength. Taken together, our findings suggest several intriguing new applications for LCs in microfluidic devices.

Controlled deformation of vesicles by flexible structured media

Liquid crystalline (LC) materials, such as actin or tubulin networks, are known to be capable of deforming the shape of cells. Here, elements of that behavior are reproduced in a synthetic system, namely, a giant vesicle suspended in a LC, which we view as a first step toward the preparation of active, anisotropic hybrid systems that mimic some of the functionality encountered in biological systems. To that end, we rely on a coupled particle-continuum representation of deformable networks in a nematic LC represented at the level of a Landau-de Gennes free energy functional. Our results indicate that, depending on its elastic properties, the LC is indeed able to deform the vesicle until it reaches an equilibrium, anisotropic shape. The magnitude of the deformation is determined by a balance of elastic and surface forces. For perpendicular anchoring at the vesicle, a Saturn ring defect forms along the equatorial plane, and the vesicle adopts a pancake-like, oblate shape. For degenerate planar anchoring at the vesicle, two boojum defects are formed at the poles of the vesicle, which adopts an elongated, spheroidal shape. During the deformation, the volume of the topological defects in the LC shrinks considerably as the curvature of the vesicle increases. These predictions are confirmed by our experimental observations of spindle-like shapes in experiments with giant unilamellar vesicles with planar anchoring. We find that the tension of the vesicle suppresses vesicle deformation, whereas anchoring strength and large elastic constants promote shape anisotropy.

Stabilization of active matter by flow-vortex lattices and defect ordering

Active systems, from bacterial suspensions to cellular monolayers, are continuously driven out of equilibrium by local injection of energy from their constituent elements and exhibit turbulent-like and chaotic patterns. Here we demonstrate both theoretically and through numerical simulations, that the crossover between wet active systems, whose behaviour is dominated by hydrodynamics, and dry active matter where any flow is screened, can be achieved by using friction as a control parameter. Moreover, we discover unexpected vortex ordering at this wet-dry crossover. We show that the self organization of vortices into lattices is accompanied by the spatial ordering of topological defects leading to active crystal-like structures. The emergence of vortex lattices, which leads to the positional ordering of topological defects, suggests potential applications in the design and control of active materials.

Lattice Boltzmann simulations of a viscoelastic shear-thinning fluid

We present a hybrid lattice Boltzmann algorithm for the simulation of flow glass-forming fluids, characterized by slow structural relaxation, at the level of the Navier-Stokes equation. The fluid is described in terms of a nonlinear integral constitutive equation, relating the stress tensor locally to the history of flow. As an application, we present results for an integral nonlinear Maxwell model that combines the effects of (linear) viscoelasticity and (nonlinear) shear thinning. We discuss the transient dynamics of velocities, shear stresses, and normal stress differences in planar pressure-driven channel flow, after switching on (startup) and off (cessation) of the driving pressure. This transient dynamics depends nontrivially on the channel width due to an interplay between hydrodynamic momentum diffusion and slow structural relaxation.

Biphasic, lyotropic, active nematics

We perform dynamical simulations of a two-dimensional active nematic fluid in coexistence with an isotropic fluid. Drops of active nematic become elongated, and an effective anchoring develops at the nematic-isotropic interface. The activity also causes an undulatory instability of the interface. This results in defects of positive topological charge being ejected into the nematic, leaving the interface with a diffuse negative charge. Quenching the active lyotropic fluid results in a steady state in which phase-separating domains are elongated and then torn apart by active stirring.

Interfacial motion in flexo- and order-electric switching between nematic filled states

We consider a nematic liquid crystal, in coexistence with its isotropic phase, in contact with a substrate patterned with rectangular grooves. In such a system the nematic phase may fill the grooves without the occurrence of complete wetting. There may exist multiple (meta)stable filled states, each characterized by the type of distortion (bend or splay) in each corner of the groove and by the shape of the nematic-isotropic interface, and additionally the plateaux that separate the grooves may be either dry or wet with a thin layer of nematic. Using numerical simulations, we analyse the dynamical response of the system to an externally-applied electric field, with the aim of identifying switching transitions between these filled states. We find that order-electric coupling between the fluid and the field provides a means of switching between states where the plateaux between grooves are dry and states where they are wetted by a nematic layer, without affecting the configuration of the nematic within the groove. We find that flexoelectric coupling may change the nematic texture in the groove, provided that the flexoelectric coupling differentiates between the types of distortion at the corners of the substrate. We identify intermediate stages of the transitions, and the role played by the motion of the nematic-isotropic interface. We determine quantitatively the field magnitudes and orientations required to effect each type of transition.

Quadrupolar particles in a nematic liquid crystal: effects of particle size and shape

We investigate the effects of particle size and shape on the quadrupolar (Saturn-ring-like) defect structures formed by a nematic liquid crystal around nm-sized and mum -sized particles with spherical and spherocylindrical shapes. We also report results for the potentials of mean force in our systems, calculated using a mesoscale theory for the tensor order parameter Q of the nematic. Our results indicate that for pairs of nm-sized particles in close proximity, the nematic forms "entangled hyperbolic" defect structures regardless of the shape of the nanoparticles. In our calculations with nanoparticles we did not observe any other entangled or unentangled defect structures, in contrast to what was reported for pairs of mum -sized spherical particles. Such a finding suggests that the "entangled hyperbolic" defect structures are the most stable for pairs of nanoparticles in close proximity. For pairs of mum -sized particles, our results indicate that the nematic forms entangled "figure-of-eight" defect structures around pairs of spheres and spherocylinders. Our results suggest that the transition between "entangled hyperbolic" and figure-of-eight defect structures takes place when the diameter of the particle is between D=100 nm and 1 microm . We have also calculated the torques that develop when pairs of spherocylindrical nanoparticles in a nematic approach each other. Our calculations suggest that the nematic-mediated interactions between the nm-sized particles are fairly strong, up to 5700 k{B}T for the case of pairs of spherocylindrical nanoparticles arranged with their long axis parallel to each other. Furthermore, these interactions can make the particles to bind together at specific locations, and thus could be used to assemble the particles into ordered structures with different morphologies.

Non-isothermal model for nematic spherulite growth

A computational study of the growth of two-dimensional nematic spherulites in an isotropic phase was performed using a Landau-de Gennes-type quadrupolar tensor order parameter model for the first-order isotropic/nematic transition of 5CB (pentylcyanobiphenyl). An energy balance, taking anisotropy into account, was derived and incorporated into the time-dependent model. Growth laws were determined for two different spherulite morphologies of the form t(n), with and without the inclusion of thermal effects. Results show that incorporation of the thermal energy balance correctly predicts the transition of the growth law exponent from the volume driven regime (n=1) to the thermally limited regime (approaching n = 1/2), agreeing well with experimental observations. An interfacial nematodynamic model is used to gain insight into the interactions that result in the progression of different spherulite growth regimes.

Nanoparticles in nematic liquid crystals: Interactions with nanochannels

A mesoscale theory for the tensor order parameter Q is used to investigate the structures that arise when spherical nanoparticles are suspended in confined nematic liquid crystals (NLCs). The NLC is "sandwiched" between a wall and a small channel. The potential of mean force is determined between particles and the bottom of the channels or between several particles. Our results suggest that strong NLC-mediated interactions between the particles and the sidewalls of the channels, on the order of hundreds of k(B)T, arise when the colloids are inside the channels. The magnitude of the channel-particle interactions is dictated by a combination of two factors, namely, the type of defect structures that develop when a nanoparticle is inside a channel, and the degree of ordering of the nematic in the region between the colloid and the nanochannel. The channel-particle interactions become stronger as the nanoparticle diameter becomes commensurate with the nanochannel width. Nanochannel geometry also affects the channel-particle interactions. Among the different geometries considered, a cylindrical channel seems to provide the strongest interactions. Our calculations suggest that small variations in geometry, such as removing the sharp edges of the channels, can lead to important reductions in channel-particle interactions. Our calculations for systems of several nanoparticles indicate that linear arrays of colloids with Saturn ring defects, which for some physical conditions are not stable in a bulk system, can be stabilized inside the nanochannels. These results suggest that nanochannels and NLCs could be used to direct the assembly of nanoparticles into ordered arrays with unusual morphologies.

Steady-state hydrodynamic instabilities of active liquid crystals: hybrid lattice Boltzmann simulations

We report hybrid lattice Boltzmann (HLB) simulations of the hydrodynamics of an active nematic liquid crystal sandwiched between confining walls with various anchoring conditions. We confirm the existence of a transition between a passive phase and an active phase, in which there is spontaneous flow in the steady state. This transition is attained for sufficiently "extensile" rods, in the case of flow-aligning liquid crystals, and for sufficiently "contractile" ones for flow-tumbling materials. In a quasi-one-dimensional geometry, deep in the active phase of flow-aligning materials, our simulations give evidence of hysteresis and history-dependent steady states, as well as of spontaneous banded flow. Flow-tumbling materials, in contrast, rearrange themselves so that only the two boundary layers flow in steady state. Two-dimensional simulations, with periodic boundary conditions, show additional instabilities, with the spontaneous flow appearing as patterns made up of "convection rolls." These results demonstrate a remarkable richness (including dependence on anchoring conditions) in the steady-state phase behavior of active materials, even in the absence of external forcing; they have no counterpart for passive nematics. Our HLB methodology, which combines lattice Boltzmann for momentum transport with a finite difference scheme for the order parameter dynamics, offers a robust and efficient method for probing the complex hydrodynamic behavior of active nematics.