Active transport in a channel: stabilisation by flow or thermodynamics

Channel instability in binary mixtures with differential diffusivity

Rupture dynamics and pinch-off phenomena are fundamental for understanding instabilities in fluid dynamics and biological systems. In this study, we investigate the rupture of two-dimensional, channel-like configurations in a binary mixture of particles with differential diffusivities. Through computational simulations, we analyze the evolution of this instability under various conditions, identifying key parameters such as aspect ratio, particle density, and drift strength that influence the system's stability. While its behavior resembles the Plateau-Rayleigh instability (PRI), the underlying mechanism differs fundamentally, as PRI is restricted to three-dimensional systems. Interestingly, similar instabilities have been observed in chiral fluids, further supporting the existence of a novel instability mechanism unique to two-dimensional non-equilibrium systems. Our results suggest that this phenomenon is not a finite-size effect, but rather an intrinsic property of systems with differential diffusivities, offering new insights into pattern formation and instability dynamics in active matter.

Solute dispersion in pre-turbulent confined active nematics

We investigate the dispersion of solutes in active nematic fluids confined to narrow channels based on simulations of nematohydrodynamics. The study focuses on two pre-turbulent regimes: oscillatory flow, with net mass flux, and dancing flow, without net flux. Non-diffusive tracers exhibit markedly different behaviors in oscillatory and dancing flows. By contrast, the hydrodynamic dispersion of solutes driven by active flows, both in the oscillatory and dancing flows, are similar and can be described by an extension of the Taylor-Aris law. This study contributes to our understanding of micromixing in active flows both in nature and in applications.

Vertex model with internal dissipation enables sustained flows

Complex tissue flows in epithelia are driven by intra- and inter-cellular processes that generate, maintain, and coordinate mechanical forces. There has been growing evidence that cell shape anisotropy, manifested as nematic order, plays an important role in this process. Here we extend an active nematic vertex model by replacing substrate friction with internal viscous dissipation, dominant in epithelia not supported by a substrate or the extracellular matrix, which are found in many early-stage embryos. When coupled to cell shape anisotropy, the internal viscous dissipation allows for long-range velocity correlations and thus enables the spontaneous emergence of flows with a large degree of spatiotemporal organisation. We demonstrate sustained flow in epithelial sheets confined to a channel, providing a link between the cell-level vertex model of tissue dynamics and continuum active nematics, whose behaviour in a channel is theoretically understood and experimentally realisable. Our findings also show a simple mechanism that could account for collective cell migration correlated over distances large compared to the cell size, as observed during morphogenesis.

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.

Periodic orbits, pair nucleation, and unbinding of active nematic defects on cones

Geometric confinement and topological constraints present promising means of controlling active materials. By combining analytical arguments derived from the Born-Oppenheimer approximation with numerical simulations, we investigate the simultaneous impact of confinement together with curvature singularity by characterizing the dynamics of an active nematic on a cone. Here, the Born-Oppenheimer approximation means that textures can follow defect positions rapidly on the timescales of interest. Upon imposing strong anchoring boundary conditions at the base of a cone, we find a rich phase diagram of multidefect dynamics, including exotic periodic orbits of one or two +1/2 flank defects, depending on activity and nonquantized geometric charge at the cone apex. By characterizing the transitions between these ordered dynamical states, we present detailed understanding of (i) defect unbinding, (ii) defect absorption, and (iii) defect pair nucleation at the apex. Numerical simulations confirm theoretical predictions of not only the nature of the circular orbits but also defect unbinding from the apex.

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.

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.

Viscoelastic confinement induces periodic flow reversals in active nematics

We use linear stability analysis and hybrid lattice Boltzmann simulations to study the dynamical behavior of an active nematic confined in a channel made of viscoelastic material. We find that the quiescent, ordered active nematic is unstable above a critical activity. The transition is to a steady flow state for high elasticity of the channel surroundings. However, below a threshold elastic modulus, the system produces spontaneous oscillations with periodic flow reversals. We provide a phase diagram that highlights the region where time-periodic oscillations are observed and explain how they are produced by the interplay of activity and viscoelasticity. Our results suggest experiments to study the role of viscoelastic confinement in the spatiotemporal organization and control of active matter.

Spontaneous flows and dynamics of full-integer topological defects in polar active matter

Polar active matter of self-propelled particles sustain spontaneous flows through the full-integer topological defects. We study theoretically the incompressible flow profiles around ±1 defects induced by polar and dipolar active forces. We show that dipolar forces induce vortical flows around the +1 defect, while the flow around the -1 defect has an 8-fold rotational symmetry. The vortical flow changes its chirality near the +1 defect core in the absence of the friction with a substrate. We show analytically that the flow induced by polar active forces is vortical near the +1 defect and is 4-fold symmetric near the -1 defect, while it becomes uniform in the far-field. For a pair of oppositely charged defects, this polar flow contributes to a mutual interaction force that depends only on the orientation of the defect pair relative to the background polarization, and that enhances defect pair annihilation. This is in contradiction with the effect of dipolar active forces which decay inversely proportional with the defect separation distance. As such, our analyses reveals a long-ranged mechanism for the pairwise interaction between topological defects in polar active matter.

Data-Driven Discovery of Active Nematic Hydrodynamics

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatiotemporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

Flow transitions and length scales of a channel-confined active nematic

We perform lattice Boltzmann simulations of an active nematic fluid confined in a two-dimensional channel to study the range of flow states that are stabilised by the confinement: unidirectional flow, oscillatory flow, the dancing state, localised active turbulence and fully-developed active turbulence. We analyse the flows in Fourier space, and measure a range of different length scales which describe the flows. We argue that the different states occur as a result of flow instabilities inherent to the system. As a consequence the characteristic length scale for oscillatory flow, the dancing state and localised active turbulence is set by the channel width. Fully-developed active turbulence occurs only when the channel width is larger than the intrinsic, active length scale of the bulk fluid. The results clarify why the activity number is a control parameter for the flow transitions.

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.

Activity pulses induce spontaneous flow reversals in viscoelastic environments

Complex interactions between cellular systems and their surrounding extracellular matrices are emerging as important mechanical regulators of cell functions, such as proliferation, motility and cell death, and such cellular systems are often characterized by pulsating actomyosin activities. Here, using an active gel model, we numerically explore spontaneous flow generation by activity pulses in the presence of a viscoelastic medium. The results show that cross-talk between the activity-induced deformations of the viscoelastic surroundings and the time-dependent response of the active medium to these deformations can lead to the reversal of spontaneously generated active flows. We explain the mechanism behind this phenomenon based on the interaction between the active flow and the viscoelastic medium. We show the importance of relaxation time scales of both the polymers and the active particles and provide a phase space over which such spontaneous flow reversals can be observed. Our results suggest new experiments investigating the role of controlled pulses of activity in living systems ensnared in complex mircoenvironments.

Confinement-Induced Self-Pumping in 3D Active Fluids

Two dimensional active fluids display a transition from turbulent to coherent flow upon decreasing the size of the confining geometry. A recent experiment suggests that the behavior in three dimensions is remarkably different; emergent flows transition from turbulence to coherence upon increasing the confinement height to match the width. Using a simple hydrodynamic model of a suspension of extensile rodlike units, we provide the theoretical explanation for this puzzling behavior. Furthermore, using extensive numerical simulations supported by theoretical arguments, we map out the conditions that lead to coherent flows and elucidate the critical role played by the aspect ratio of the confining channel. The mechanism that we identify applies to a large class of symmetries and propulsion mechanisms, leading to a unified set of design principles for self-pumping 3D active fluids.

Flow States and Transitions of an Active Nematic in a Three-Dimensional Channel

We use active nematohydrodynamics to study the flow of an active fluid in a 3D microchannel, finding a transition between active turbulence and regimes where there is a net flow along the channel. We show that the net flow is only possible if the active nematic is flow aligning and that, in agreement with experiments, the appearance of the net flow depends on the aspect ratio of the channel cross section. We explain our results in terms of when the hydrodynamic screening due to the channel walls allows the emergence of vortex rolls across the channel.

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.

Active phase field crystal systems with inertial delay and underdamped dynamics

Active matter systems often are well approximated as overdamped, meaning that any inertial momentum is immediately dissipated by the environment. On the other hand, especially for macroscopic systems but also for many mesoscopic ones particle mass can become relevant for the dynamics. For such systems we recently proposed an underdamped continuum model which captures translationally inertial dynamics via two contributions. First, convection and second a damping time scale of inertial motion. In this paper, we ask how both of these features influence the collective behavior compared to overdamped dynamics by studying the example of the active phase field crystal model. We first focus on the case of suppressed convection to study the role of the damping time. We quantify that the relaxation process to the steady collective motion state is considerably prolonged with damping time due to the increasing occurrence of transient groups of circularly moving density peaks. Finally, we illustrate the fully underdamped case with convection. Instead of collective motion of density peaks we then find a coexistence of constant high and low density phases reminiscent of motility-induced phase separation.

Mean field approach of dynamical pattern formation in underdamped active matter with short-ranged alignment and distant anti-alignment interactions

Many active matter systems, especially on the microscopic scale, are well approximated as overdamped, meaning that any inertial momentum is immediately dissipated by the environment. On the other hand, especially for macroscopic active systems but also for many mesoscopic systems the time scale of translational inertial motion can become large enough to be relevant for the dynamics. This raises the question how collective dynamics and the resulting states in active matter are influenced by inertia. Therefore, we propose a coarse-grained continuum model for underdamped active matter based on a mean field description for passive systems. Furthermore, as an example, we apply the model to a system with interactions that support an alignment on short distances and an anti-alignment on longer length scales as known in the context of pattern formation due to orientational interactions. Our numerical calculations of the under- and overdamped dynamics both predict a structured laning state. However, activity induced convective flows that are only present in the underdamped model destabilize this state when the anti-alignment is weakened, leading to a collective motion state which does not occur in the overdamped limit. A turbulent transition regime between the two states can be characterized by strong density fluctuations and the absence of global ordering.

Propagation of active nematic-isotropic interfaces on substrates

Motivated by results for the propagation of active-passive interfaces of bacterial Serratia marcescens swarms [Nat. Commun., 2018, 9, 5373], we used a hydrodynamic multiphase model to investigate the propagation of interfaces of active nematics on substrates. We characterized the active nematic phase of the model through the calculation of the spatial and temporal auto correlation functions and the energy spectrum and discussed its description of the statistical dynamics of the swarms reported in the experiment. We then studied the propagation of circular and flat active-passive interfaces. We found that the closing time of the circular passive domain decays quadratically with the activity and that the structure factor of the flat interface is similar to that reported for the swarms, with an activity dependent exponent. Finally, the effect of the substrate friction was investigated. We found an activity dependent threshold, above which the turbulent active nematic forms isolated islands that shrink until the system becomes isotropic and below which the active nematic expands, with a well defined propagating interface. We also found that the interface becomes static in the presence of a friction gradient.

Active nematic-isotropic interfaces in channels

We use numerical simulations to investigate the hydrodynamic behavior of the interface between nematic (N) and isotropic (I) phases of a confined active liquid crystal. At low activities, a stable interface with constant shape and velocity is observed separating the two phases. For nematics in homeotropic channels, the velocity of the interface at the NI transition increases from zero (i) linearly with the activity for contractile systems and (ii) quadratically for extensile ones. Interestingly, the nematic phase expands for contractile systems while it contracts for extensile ones, as a result of the active forces at the interface. Since both activity and temperature affect the stability of the nematic, for active nematics in the stable regime the temperature can be tuned to observe static interfaces, providing an operational definition for the coexistence of active nematic and isotropic phases. At higher activities, beyond the stable regime, an interfacial instability is observed for extensile nematics. In this regime defects are nucleated at the interface and move away from it. The dynamics of these defects is regular and persists asymptotically for a finite range of activities. We used an improved hybrid model of finite differences and the lattice Boltzmann method with a multi-relaxation-time collision operator, the accuracy of which allowed us to characterize the dynamics of the distinct interfacial regimes.