Active nematics with anisotropic friction: the decisive role of the flow aligning parameter

Active nematic ratchet in asymmetric obstacle arrays

We numerically investigate the effect of an asymmetric periodic obstacle array in a two-dimensional active nematic. We find that activity in conjunction with the asymmetry leads to a ratchet effect or unidirectional flow of the fluid along the asymmetry direction. The directional flow is still present even in the active turbulent phase when the gap between obstacles is sufficiently small. We demonstrate that the dynamics of the topological defects transition from flow mirroring to smectic-like as the gap between obstacles is made smaller, and explain this transition in terms of the pinning of negative winding number defects between obstacles. This also leads to a nonmonotonic ratchet effect magnitude as a function of obstacle size, so that there is an optimal obstacle size for ratcheting at fixed activity.

Controlling active turbulence by activity patterns

By patterning activity in space, one can control active turbulence. To show this, we use Doi's hydrodynamic equations of a semidilute solution of active rods. A linear stability analysis reveals the resting isotropic fluid to be unstable above an absolute pusher activity. The emergent activity-induced paranematic state displays active turbulence, which we characterize by different quantities including the energy spectrum, which shows the typical power-law decay with exponent -4. Then, we control the active turbulence by a square lattice of circular spots where activity is switched off. In the parameter space lattice constant versus surface-to-surface distance of the spots, we identify different flow states. Most interestingly, for lattice constants below the vorticity correlation length and for spot distances smaller than the nematic coherence length, we observe a multi-lane flow state, where flow lanes with alternating flow directions are separated by a street of vortices. The flow pattern displays pronounced multistability and also appears transiently at the transition to the isotropic active-turbulence state. At larger lattice constants a trapped vortex state is identified with a non-Gaussian vorticity distribution due to the low flow vorticity at the spots. It transitions to conventional active turbulence for increasing spot distance.

Vortex Lattices in Active Nematics with Periodic Obstacle Arrays

We numerically model a two-dimensional active nematic confined by a periodic array of fixed obstacles. Even in the passive nematic, the appearance of topological defects is unavoidable due to planar anchoring by the obstacle surfaces. We show that a vortex lattice state emerges as activity is increased, and that this lattice may be tuned from "ferromagnetic" to "antiferromagnetic" by varying the gap size between obstacles. We map the rich variety of states exhibited by the system as a function of distance between obstacles and activity, including a pinned defect state, motile defects, the vortex lattice, and active turbulence. We demonstrate that the flows in the active turbulent phase can be tuned by the presence of obstacles, and explore the effects of a frustrated lattice geometry on the vortex lattice phase.

Elasticity tunes mechanical stress localization around active topological defects

Mechanical stresses are increasingly found to be associated with various biological functionalities. At the same time, topological defects are being identified across a diverse range of biological systems and are points of localized mechanical stress. It is therefore important to ask how mechanical stress localization around topological defects is controlled. Here, we use continuum simulations of nonequilibrium, fluctuating and active nematics to explore the patterns of stress localization, as well as their extent and intensity around topological defects. We find that by increasing the orientational elasticity of the material, the isotropic stress pattern around topological defects is changed substantially, from a stress dipole characterized by symmetric compression-tension regions around the core of the defect, to a localized stress monopole at the defect position. Moreover, we show that elastic anisotropy alters the extent and intensity of the stresses, and can result in the dominance of tension or compression around defects. Finally, including both nonequilibrium fluctuations and active stress generation, we find that the elastic constant tunes the relative effect of each, leading to the flipping of tension and compression regions around topological defects. This flipping of the tension-compression regions only by changing the elastic constant presents an interesting, simple, way of switching the dynamic behavior in active matter by changing a passive material property. We expect these findings to motivate further exploration tuning stresses in active biological materials by varying material properties of the constituent units.

Vorticity phase separation and defect lattices in the isotropic phase of active liquid crystals

We use numerical simulations and linear stability analysis to study the dynamics of an active liquid crystal film on a substrate in the regime where the passive system would be isotropic. Extensile activity builds up local orientational order and destabilizes the quiescent isotropic state above a critical activity, eventually resulting in spatiotemporal chaotic dynamics akin to the one observed ubiquitously in the nematic state. Here we show that tuning substrate friction yields a variety of emergent structures at intermediate activity, including lattices of flow vortices with associated regular arrangements of topological defects and a new state where flow vortices trap pairs of +1/2 defect that chase each other's tail. These chiral units spontaneously pick the sense of rotation and organize in a hexagonal lattice, surrounded by a diffuse flow of opposite rotation to maintain zero net vorticity. The length scale of these emergent structures is set by the screening length of the flow, controlled by the shear viscosity η and the substrate friction Γ, and can be captured by simple mode selection of the vortical flows. We demonstrate that the emergence of coherent structures can be interpreted as a phase separation of vorticity, where friction plays a role akin to that of birth/death processes in breaking conservation of the phase separating species and selecting a characteristic scale for the patterns. Our work shows that friction provides an experimentally accessible tuning parameter for designing controlled active flows.

Measuring response functions of active materials from data

From flocks of birds to biomolecular assemblies, systems in which many individual components independently consume energy to perform mechanical work exhibit a wide array of striking behaviors. Methods to quantify the dynamics of these so-called active systems generally aim to extract important length or time scales from experimental fields. Because such methods focus on extracting scalar values, they do not wring maximal information from experimental data. We introduce a method to overcome these limitations. We extend the framework of correlation functions by taking into account the internal headings of displacement fields. The functions we construct represent the material response to specific types of active perturbation within the system. Utilizing these response functions we query the material response of disparate active systems composed of actin filaments and myosin motors, from model fluids to living cells. We show we can extract critical length scales from the turbulent flows of an active nematic, anticipate contractility in an active gel, distinguish viscous from viscoelastic dissipation, and even differentiate modes of contractility in living cells. These examples underscore the vast utility of this method which measures response functions from experimental observations of complex active systems.

Friction-mediated phase transition in confined active nematics

Using a minimal continuum model, we investigate the interplay between circular confinement and substrate friction in active nematics. Upon increasing the friction from low to high, we observe a dynamical phase transition from a circulating flow phase to an anisotropic flow phase in which the flow tends to align perpendicular to the nematic director at the boundary. We demonstrate that both the flow structure and dynamic correlations in the latter phase differ from those of an unconfined, active turbulent system and may be controlled by the prescribed nematic boundary conditions. Our results show that substrate friction and geometric confinement act as valuable control parameters in active nematics.

Flow patterns and defect dynamics of active nematic liquid crystals under an electric field

The effects of an electric field on the flow patterns and defect dynamics of two-dimensional active nematic liquid crystals are numerically investigated. We found that field-induced director reorientation causes anisotropic active turbulence characterized by enhanced flow perpendicular to the electric field. The average flow speed and its anisotropy are maximized at an intermediate field strength. Topological defects in the anisotropic active turbulence are localized and show characteristic dynamics with simultaneous creation of two pairs of defects. A laning state characterized by stripe domains with alternating flow directions is found at a larger field strength near the transition to the uniformly aligned state. We obtained periodic oscillations between the laning state and active turbulence, which resembles an experimental observation of active nematics subject to anisotropic friction.

Coupling the topological defect phase to the extrinsic curvature in nematic shells

In two dimensional nematics, topological defects are point like singularities with both a charge and a phase. We study the topological defects within curved nematic textures on the surface of a cylinder. This allows us to isolate the effect of extrinsic curvature on the structure of the topological defect. By minimizing the energy associated with distortions in the nematic director around the core of a defect, we show that the phase of the topological defect is coupled to the orientation of the cylinder. This coupling depends on the relative energetic cost associated with splay, bend and twist distortions of the nematic director. We identify a bistability in the phase of the defects when twist deformations dominate. Finally, we show a similar effect for integer charge topological defects.

Topography-induced large-scale antiparallel collective migration in vascular endothelium

Collective migration of vascular endothelial cells is central for embryonic development, angiogenesis, and wound closure. Although physical confinement of cell assemblies has been shown to elicit specific patterns of collective movement in various cell types, endothelial migration in vivo often occurs without confinement. Here we show that unconfined endothelial cell monolayers on microgroove substrates that mimic the anisotropic organization of the extracellular matrix exhibit a specific type of collective movement that takes the form of a periodic pattern of antiparallel cell streams. We further establish that the development of these streams requires intact cell-cell junctions and that stream sizes are particularly sensitive to groove depth. Finally, we show that modeling the endothelial cell sheet as an active fluid with the microgrooves acting as constraints on cell orientation predicts the occurrence of the periodic antiparallel cell streams as well as their lengths and widths. We posit that in unconfined cell assemblies, physical factors that constrain or bias cellular orientation such as anisotropic extracellular matrix cues or directed flow-derived shear forces dictate the pattern of collective cell movement.

The physical environment dictates the emergence of specific patterns of collective cell migration. Here, authors show that unconfined endothelial monolayers on microgroove substrates exhibit an original pattern of antiparallel cell streams.

Role of Advective Inertia in Active Nematic Turbulence

Suspensions of active agents with nematic interactions exhibit complex spatiotemporal dynamics such as mesoscale turbulence. Since the Reynolds number of microscopic flows is very small on the scale of individual agents, inertial effects are typically excluded in continuum theories of active nematic turbulence. Whether active stresses can collectively excite inertial flows is currently unclear. To address this question, we investigate a two-dimensional continuum theory for active nematic turbulence. In particular, we compare mesoscale turbulence with and without the effects of advective inertia. We find that inertial effects can influence the flow already close to the onset of the turbulent state and, moreover, give rise to large-scale fluid motion for strong active driving. A detailed analysis of the kinetic energy budget reveals an energy transfer to large scales mediated by inertial advection. While this transfer is small in comparison to energy injection and dissipation, its effects accumulate over time. The inclusion of friction, which is typically present in experiments, can compensate for this effect. The findings suggest that the inclusion of inertia and friction may be necessary for dynamically consistent theories of active nematic turbulence.

Submersed micropatterned structures control active nematic flow, topology, and concentration

Coupling between flows and material properties imbues rheological matter with its wide-ranging applicability, hence the excitement for harnessing the rheology of active fluids for which internal structure and continuous energy injection lead to spontaneous flows and complex, out-of-equilibrium dynamics. We propose and demonstrate a convenient, highly tunable method for controlling flow, topology, and composition within active films. Our approach establishes rheological coupling via the indirect presence of fully submersed micropatterned structures within a thin, underlying oil layer. Simulations reveal that micropatterned structures produce effective virtual boundaries within the superjacent active nematic film due to differences in viscous dissipation as a function of depth. This accessible method of applying position-dependent, effective dissipation to the active films presents a nonintrusive pathway for engineering active microfluidic systems.

Properties of twisted topological defects in 2D nematic liquid crystals

Topological defects are one of the most conspicuous features of liquid crystals. In two dimensional nematics, they have been shown to behave effectively as particles with both charge and orientation, which dictate their interactions. Here, we study "twisted" defects that have a radially dependent orientation. We find that twist can be partially relaxed through the creation and annihilation of defect pairs. By solving the equations for defect motion and calculating the forces on defects, we identify four distinct elements that govern the relative relaxational motion of interacting topological defects, namely attraction, repulsion, co-rotation and co-translation. The interaction of these effects can lead to intricate defect trajectories, which can be controlled by setting relevant timescales.

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.

Memory effects, arches and polar defect ordering at the cross-over from wet to dry active nematics

We use analytic arguments and numerical solutions of the continuum, active nematohydrodynamic equations to study how friction alters the behaviour of active nematics. Concentrating on the case where there is nematic ordering in the passive limit, we show that, as the friction is increased, memory effects become more prominent and +1/2 topological defects leave increasingly persistent trails in the director field as they pass. The trails are preferential sites for defect formation and they tend to impose polar order on any new +1/2 defects. In the absence of noise and for high friction, it becomes very difficult to create defects, but trails formed by any defects present at the beginning of the simulations persist and organise into parallel arch-like patterns in the director field. We show aligned arches of equal width are approximate steady state solutions of the equations of motion which co-exist with the nematic state. We compare our results to other models in the literature, in particular dry systems with no hydrodynamics, where trails, arches and polar defect ordering have also been observed.

Spectral energy analysis of bulk three-dimensional active nematic turbulence

We perform energy spectrum analysis of the active turbulence in a 3D bulk active nematic using continuum numerical modelling. Specifically, we calculate the spectra of the two main energy contributions - kinetic energy and nematic elastic energy - and combine this with the geometrical analysis of the nematic order and flow fields, based on direct defect tracking and calculation of autocorrelations. We show that the active nematic elastic energy is concentrated at scales corresponding to the effective defect-to-defect separation, scaling with activity as ∼ζ0.5, whereas the kinetic energy is largest at somewhat larger scales of typically several 100 nematic correlation lengths. Nematic biaxiality is shown to have no role in active turbulence at most length scales, but can affect the nematic elastic energy by an order of magnitude at scales of the active defect core size. The effect of an external aligning field on the 3D active turbulence is explored, showing a transition from an effective active turbulent to an aligned regime. The work is aimed at providing a contribution towards understanding active turbulence in general three-dimensions, from the perspective of main energy-relevant mechanisms at different length scales of the system.

Spontaneous helical flows in active nematics lying on a cylindrical surface

Within the framework of the two-dimensional Ericksen-Leslie model, we explore the effect of geometric confinement on the spontaneous flow of active nematic gels. The nematic particles are assumed to flow on a cylindrical surface, while a degenerate tangential anchoring is enforced. Using the linear approximation of the motion equations, we show that there is a close interplay among extrinsic curvature, flow, director alignment, and activity. We find that the extrinsic curvature promotes the director alignment parallel to the cylindrical axis and is responsible for raising the critical threshold with respect to the flat case. Our analysis reveals a very rich scenario where the key quantities are the activity coefficient, the tumbling parameter, and the anisotropic viscosity ratio. Thus, solutions can exhibit a double periodicity in both the azimuthal and axial variables. As a consequence, the velocity field can make a finite angle with the cylinder axis and the active flow winds on the surface with a helical pattern, while the director oscillates around the cylinder generators. Our results can be validated on thin layers of nematic gels placed between two concentric cylinders and suggest which material properties are most suited for the design of active microfluidic devices.