Spontaneous mirror-symmetry breaking induces inverse energy cascade in 3D active fluids

Vortex Pattern Stabilization in Thin Films Resulting from Shear Thickening of Active Suspensions

The need for structuring on micrometer scales is abundant, for example, in view of phononic applications. We here outline a novel approach based on the phenomenon of active turbulence on the mesoscale. As we demonstrate, a shear-thickening carrier fluid of active microswimmers intrinsically stabilizes regular vortex patterns of otherwise turbulent active suspensions. The fluid self-organizes into a periodically structured nonequilibrium state. Introducing additional passive particles of intermediate size leads to regular spatial organization of these objects. Our approach opens a new path toward functionalization through patterning of thin films and membranes.

Pattern formation by turbulent cascades

Fully developed turbulence is a universal and scale-invariant chaotic state characterized by an energy cascade from large to small scales at which the cascade is eventually arrested by dissipation1-6. Here we show how to harness these seemingly structureless turbulent cascades to generate patterns. Pattern formation entails a process of wavelength selection, which can usually be traced to the linear instability of a homogeneous state7. By contrast, the mechanism we propose here is fully nonlinear. It is triggered by the non-dissipative arrest of turbulent cascades: energy piles up at an intermediate scale, which is neither the system size nor the smallest scales at which energy is usually dissipated. Using a combination of theory and large-scale simulations, we show that the tunable wavelength of these cascade-induced patterns can be set by a non-dissipative transport coefficient called odd viscosity, ubiquitous in chiral fluids ranging from bioactive to quantum systems8-12. Odd viscosity, which acts as a scale-dependent Coriolis-like force, leads to a two-dimensionalization of the flow at small scales, in contrast with rotating fluids in which a two-dimensionalization occurs at large scales4. Apart from odd viscosity fluids, we discuss how cascade-induced patterns can arise in natural systems, including atmospheric flows13-19, stellar plasma such as the solar wind20-22, or the pulverization and coagulation of objects or droplets in which mass rather than energy cascades23-25.

Topology, vorticity, and limit cycle in a stabilized Kuramoto-Sivashinsky equation

A noisy stabilized Kuramoto-Sivashinsky equation is analyzed by stochastic decomposition. For values of the control parameter for which periodic stationary patterns exist, the dynamics can be decomposed into diffusive and transverse parts which act on a stochastic potential. The relative positions of stationary states in the stochastic global potential landscape can be obtained from the topology spanned by the low-lying eigenmodes which interconnect them. Numerical simulations confirm the predicted landscape. The transverse component also predicts a universal class of vortex-like circulations around fixed points. These drive nonlinear drifting and limit cycle motion of the underlying periodic structure in certain regions of parameter space. Our findings might be relevant in studies of other nonlinear systems such as deep learning neural networks.

Incompressible Polar Active Fluids with Quenched Random Field Disorder in Dimensions d>2

We present a hydrodynamic theory of incompressible polar active fluids with quenched random field disorder. This theory shows that such fluids can overcome the disruption caused by the quenched disorder and move coherently, in the sense of having a nonzero mean velocity in the hydrodynamic limit. However, the scaling behavior of this class of active systems cannot be described by linearized hydrodynamics in spatial dimensions between 2 and 5. Nonetheless, we obtain the exact dimension-dependent scaling exponents in these dimensions.

Helical flow states in active nematics

We show that confining extensile nematics in three-dimensional (3D) channels leads to the emergence of two self-organized flow states with nonzero helicity. The first is a pair of braided antiparallel streams-this double helix occurs when the activity is moderate, anchoring negligible, and reduced temperature high. The second consists of axially aligned counter-rotating vortices-this grinder train arises between spontaneous axial streaming and the vortex lattice. These two unanticipated helical flow states illustrate the potential of active fluids to break symmetries and form complex but organized spatiotemporal structures in 3D fluidic devices.

Spontaneous chiral symmetry breaking in a random driven chemical system

Living systems have evolved to efficiently consume available energy sources using an elaborate circuitry of chemical reactions which, puzzlingly, bear a strict restriction to asymmetric chiral configurations. While autocatalysis is known to promote such chiral symmetry breaking, whether a similar phenomenon may also be induced in a more general class of configurable chemical systems—via energy exploitation—is a sensible yet underappreciated possibility. This work examines this question within a model of randomly generated complex chemical networks. We show that chiral symmetry breaking may occur spontaneously and generically by harnessing energy sources from external environmental drives. Key to this transition are intrinsic fluctuations of achiral-to-chiral reactions and tight matching of system configurations to the environmental drives, which together amplify and sustain diverged enantiomer distributions. These asymmetric states emerge through steep energetic transitions from the corresponding symmetric states and sharply cluster as highly-dissipating states. The results thus demonstrate a generic mechanism in which energetic drives may give rise to homochirality in an otherwise totally symmetrical environment, and from an early-life perspective, might emerge as a competitive, energy-harvesting advantage.

“A hallmark of living systems is their homochirality, the selection of specific mirror symmetry in their molecules. Here, the authors show that chiral symmetry can be spontaneously broken in complex, random chemical systems via exploitation of environmental energy sources – a possible mechanism for the emergence of homochirality in life.”

Active Extensile Stress Promotes 3D Director Orientations and Flows

We use numerical simulations and linear stability analysis to study an active nematic layer where the director is allowed to point out of the plane. Our results highlight the difference between extensile and contractile systems. Contractile stress suppresses the flows perpendicular to the layer and favors in-plane orientations of the director. By contrast extensile stress promotes instabilities that can turn the director out of the plane, leaving behind a population of distinct, in-plane regions that continually elongate and divide. This supports extensile forces as a mechanism for the initial stages of layer formation in living systems, and we show that a planar drop with extensile (contractile) activity grows into three dimensions (remains in two dimensions). The results also explain the propensity of disclination lines in three dimensional active nematics to be of twist type in extensile or wedge type in contractile materials.

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.

Global potential, topology, and pattern selection in a noisy stabilized Kuramoto-Sivashinsky equation

We formulate a general method to extend the decomposition of stochastic dynamics developed by Ao et al. [J. Phys. Math. Gen. 37, L25-L30 (2004)] to nonlinear partial differential equations which are nonvariational in nature and construct the global potential or Lyapunov functional for a noisy stabilized Kuramoto-Sivashinsky equation. For values of the control parameter where singly periodic stationary solutions exist, we find a topological network of a web of saddle points of stationary states interconnected by unstable eigenmodes flowing between them. With this topology, a global landscape of the steady states is found. We show how to predict the noise-selected pattern which agrees with those from stochastic simulations. Our formalism and the topology might offer an approach to explore similar systems, such as the Navier Stokes equation.

Condensate formation and multiscale dynamics in two-dimensional active suspensions

The collective effects of microswimmers in active suspensions result in active turbulence, a spatiotemporally chaotic dynamics at mesoscale, which is characterized by the presence of vortices and jets at scales much larger than the characteristic size of the individual active constituents. To describe this dynamics, Navier-Stokes-based one-fluid models driven by small-scale forces have been proposed. Here, we provide a justification of such models for the case of dense suspensions in two dimensions (2D). We subsequently carry out an in-depth numerical study of the properties of one-fluid models as a function of the active driving in view of possible transition scenarios from active turbulence to large-scale pattern, referred to as condensate, formation induced by the classical inverse energy cascade in Newtonian 2D turbulence. Using a one-fluid model it was recently shown [M. Linkmann et al., Phys. Rev. Lett 122, 214503 (2019)10.1103/PhysRevLett.122.214503] that two-dimensional active suspensions support two nonequilibrium steady states, one with a condensate and one without, which are separated by a subcritical transition. Here, we report further details on this transition such as hysteresis and discuss a low-dimensional model that describes the main features of the transition through nonlocal-in-scale coupling between the small-scale driving and the condensate.

Learning the Physics of Pattern Formation from Images

Using a framework of partial differential equation-constrained optimization, we demonstrate that multiple constitutive relations can be extracted simultaneously from a small set of images of pattern formation. Examples include state-dependent properties in phase-field models, such as the diffusivity, kinetic prefactor, free energy, and direct correlation function, given only the general form of the Cahn-Hilliard equation, Allen-Cahn equation, or dynamical density functional theory (phase-field crystal model). Constraints can be added based on physical arguments to accelerate convergence and avoid spurious results. Reconstruction of the free energy functional, which contains nonlinear dependence on the state variable and differential or convolutional operators, opens the possibility of learning nonequilibrium thermodynamics from only a few snapshots of the dynamics.

Phase Transition to Large Scale Coherent Structures in Two-Dimensional Active Matter Turbulence

The collective motion of microswimmers in suspensions induce patterns of vortices on scales that are much larger than the characteristic size of a microswimmer, attaining a state called bacterial turbulence. Hydrodynamic turbulence acts on even larger scales and is dominated by inertial transport of energy. Using an established modification of the Navier-Stokes equation that accounts for the small-scale forcing of hydrodynamic flow by microswimmers, we study the properties of a dense suspension of microswimmers in two dimensions, where the conservation of enstrophy can drive an inverse cascade through which energy is accumulated on the largest scales. We find that the dynamical and statistical properties of the flow show a sharp transition to the formation of vortices at the largest length scale. The results show that 2D bacterial and hydrodynamic turbulence are separated by a subcritical phase transition.

Twist-induced crossover from two-dimensional to three-dimensional turbulence in active nematics

While studies of active nematics in two dimensions have shed light on various aspects of the flow regimes and topology of active matter, three-dimensional properties of topological defects and chaotic flows remain unexplored. By confining a film of active nematics between two parallel plates, we use continuum simulations and analytical arguments to demonstrate that the crossover from quasi-two-dimensional (quasi-2D) to three-dimensional (3D) chaotic flows is controlled by the morphology of the disclination lines. For small plate separations, the active nematic behaves as a quasi-2D material, with straight topological disclination lines spanning the height of the channel and exhibiting effectively 2D active turbulence. Upon increasing channel height, we find a crossover to 3D chaotic flows due to the contortion of disclinations above a critical activity. Above this critical activity highly contorted disclination lines and disclination loops are formed. We further show that these contortions are engendered by twist perturbations producing a sharp change in the curvature of disclinations.

Anomalous Chained Turbulence in Actively Driven Flows on Spheres

Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence.

Discontinuous Transition from Direct to Inverse Cascade in Three-Dimensional Turbulence

Inviscid invariants of flow equations are crucial in determining the direction of the turbulent energy cascade. In this work we investigate a variant of the three-dimensional Navier-Stokes equations that shares exactly the same ideal invariants (energy and helicity) and the same symmetries (under rotations, reflections, and scale transforms) as the original equations. It is demonstrated that the examined system displays a change in the direction of the energy cascade when varying the value of a free parameter which controls the relative weights of the triadic interactions between different helical Fourier modes. The transition from a forward to inverse cascade is shown to occur at a critical point in a discontinuous manner with diverging fluctuations close to criticality. Our work thus supports the observation that purely isotropic and three-dimensional flow configurations can support inverse energy transfer when interactions are altered and that inside all turbulent flows there is a competition among forward and backward transfer mechanisms which might lead to multiple energy-containing turbulent states.