Helical flow states in active nematics

Majorana quasiparticles and topological phases in 3D active nematics

Quasiparticles are low-energy excitations with important roles in condensed matter physics. An intriguing example is provided by Majorana quasiparticles, which are equivalent to their antiparticles. Despite being implicated in neutrino oscillations and topological superconductivity, their experimental realizations remain very rare. Here, we propose a purely classical realization of Majorana fermions in terms of three-dimensional disclination lines in active nematics. The underlying reason is the well-known equivalence, in 3D, between a [Formula: see text] local defect profile and a [Formula: see text] profile, which acts as its antiparticle. The mapping also requires proving that defect profiles transform as spinors, and activity is needed to overcome the elastic cost associated with these excitations, so they spontaneously appear in steady state. We combine topological considerations and numerics to show that active nematics under confinement spontaneously create in their interior topologically charged disclination lines and loops, akin to Majorana quasiparticles with finite momentum. Within a long channel, the phenomenology we observe resembles that of the Kitaev chain, as Majorana-like states appear near the boundaries, while a delocalized topological excitation arises in the form of a chiral disclination line. The analogy between 3D nematic defects and topological quasiparticles further suggests that active turbulence can be viewed as a topological phase, where defects percolate to form delocalized topological quasiparticles similar to those observed in the channel. We propose that three-dimensional active disclinations can be used to probe the physics of Majorana spinors at much larger scale than that for which they were originally introduced, potentially facilitating their experimental study.

Spontaneous flows and quantum analogies in heterogeneous active nematic films

Incorporating the inherent heterogeneity of living systems into models of active nematics is essential to provide a more realistic description of biological processes such as bacterial growth, cell dynamics and tissue development. Spontaneous flow of a confined active nematic is a fundamental feature of these systems, in which the role of heterogeneity has not yet been considered. We therefore determine the form of spontaneous flow transition for an active nematic film with heterogeneous activity, identifying a correspondence between the unstable director modes and solutions to Schrödinger's equation. We consider both activity gradients and steps between regions of distinct activity, finding that such variations can change the signature properties of the flow. The threshold activity required for the transition can be raised or lowered, the fluid flux can be reduced or reversed and interfaces in activity induce shear flows. In a biological context fluid flux influences the spread of nutrients while shear flows affect the behaviour of rheotactic microswimmers and can cause the deformation of biofilms. All the effects we identify are found to be strongly dependent on not simply the types of activity present in the film but also on how they are distributed.

Active nematics in corrugated channels

Active nematic fluids exhibit complex dynamics in both bulk and in simple confining geometries. However, complex confining geometries could have substantial impact on active spontaneous flows. Using multiparticle collision dynamics simulations adapted for active nematic particles, we study the dynamic behaviour of an active nematic fluid confined in a corrugated channel. The transition from a quiescent state to a spontaneous flow state occurs from a weak swirling flow to a strong coherent flow due to the presence of curved-wall induced active flows. We show that the active nematic fluid flows in corrugated channels can be understood in two different ways: (i) as the result of an early or delayed flow transition when compared with that in a flat-walled channel of appropriate width and (ii) boundary-induced active flows in the corrugations providing an effective slip velocity to the coherent flows in the bulk. Thus, our work illustrates the crucial role of corrugations of the confining boundary in dictating the flow transition and flow states of active fluids.

Active Darcy's Law

While bacterial swarms can exhibit active turbulence in vacant spaces, they naturally inhabit crowded environments. We numerically show that driving disorderly active fluids through porous media enhances Darcy's law. While purely active flows average to zero flux, hybrid active/driven flows display greater drift than purely pressure-driven flows. This enhancement is nonmonotonic with activity, leading to an optimal activity to maximize flow rate. We incorporate the active contribution into an active Darcy's law, which may serve to help understand anomalous transport of swarming in porous media.

Flow states of two dimensional active gels driven by external shear

Using a minimal hydrodynamic model, we theoretically and computationally study the Couette flow of active gels in straight and annular two-dimensional channels subject to an externally imposed shear. The gels are isotropic in the absence of externally- or activity-driven shear, but have nematic order that increases with shear rate. Using the finite element method, we determine the possible flow states for a range of activities and shear rates. Linear stability analysis of an unconfined gel in a straight channel shows that an externally imposed shear flow can stabilize an extensile fluid that would be unstable to spontaneous flow in the absence of the shear flow, and destabilize a contractile fluid that would be stable against spontaneous flow in the absence of shear flow. These results are in rough agreement with the stability boundaries between the base shear flow state and the nonlinear flow states that we find numerically for a confined active gel. For extensile fluids, we find three kinds of nonlinear flow states in the range of parameters we study: unidirectional flows, oscillatory flows, and dancing flows. To highlight the activity-driven spontaneous component of the nonlinear flows, we characterize these states by the average volumetric flow rate and the wall stress. For contractile fluids, we only find the linear shear flow and a nonlinear unidirectional flow in the range of parameters that we studied. For large magnitudes of the activity, the unidirectional contractile flow develops a boundary layer. Our analysis of annular channels shows how curvature of the streamlines in the base flow affects the transitions among flow states.

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.

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.