Activity induced isotropic-polar transition in active liquid crystals

Maria Grazia Giordano,1 Francesco Bonelli,2 Livio Nicola Carenza,1 Giuseppe Gonnella,1 and Giuseppe Negro1 1Dipartimento di Fisica, Università degli Studi di Bari and INFN, via Amendola 173, Bari, I-70126, Italy 2Dipartimento di Meccanica, Matematica e Management, DMMM, Politecnico di Bari, 70125, Bari, Italy Active fluids are intrinsically out-of-equilibrium systems due to the internal energy injection of the active constituents. We show here that a transition from a motion-less isotropic state towards a flowing polar one can be possibly driven by the sole active injection through the action of polar- hydrodynamic interactions in absence of an ad hoc free-energy which favors the development of an ordered phase. In particular, we propose an analytical argument and we perform lattice Boltzmann simulations where the appearance of large temporal fluctuations in the polar fraction of the system is observed at the transition point. Moreover, we make use of a scale-to-scale analysis to unveil the energy transfer mechanism, proving that elastic absorption plays a relevant role in the overall dynamics of the system, contrary to what reported in previous works on the usual active gel theory where this term could be factually neglected.

INTRODUCTION The physical origin of such motility-induced transition from the isotropic to an ordered state is to be related Active fluids intrinsically evolve out of equilibrium due to the intricate nemato-hydrodynamic behavior and the to internal energy injection [1, 2]. Systems of biological connected rheological properties of liquid-crystal suspen- origin – such as bacterial [3, 4] and cytoskeletal [5–7] sions. For instance, Markovich et al. [55] have recently suspensions – or synthetic realizations, as Janus parti- shown, analytically and numerically, that an external cles [8] and polyacrylic acid hydrogels [9], convert chemi- shear can induce a first-order phase transition from the cal energy stored in some reservoir into non-equilibrium isotropic to the ordered state in a (passive) polar liq- stresses which result in self-sustained flows [6, 10, 11]. uid crystal. The mechanism at base of this effect lies Recent research has proved that the mutual interactions in the renormalization of the bulk constants describing between hydrodynamics and active energy input may the ordering properties of the polarization, due to the give rise to non-equilibrium dynamical states connected strain effect of the mechanically imposed shear flow, thus to a plethora of unexpected behaviors as, for instance, allowing for the development of a polar phase in a re- super-fluidic and negative viscosity states [12–18], spon- gion of the parameter space which would otherwise be taneous flows [19–27] and active turbulence regime [3, 28– isotropic. This is of intereset for our work, since the 32]. swimming mechanism of active liquid crystals is able to Even with intrinsic structural differences, a common generate local shear flows which may possibly alter the feature to many active systems is the breaking of isotropic threshold of the orientational transition. Indeed, exten- symmetry. This may either occur at the level of the in- sile swimmers are able to expel fluid at their ends and dividual constituents [33, 34] – since their bodies often draw it across their body acting on the surrounding en- exhibit an elongated or filamentous structure – or due to vironment as out-warding force dipoles, while contractile the direction of motion [3, 35–37], resulting into the emer- swimmers act as in-warding force dipoles, leading to a gence of local patterns with polar/nematic or even hex- reversed pattern [12]. In both cases, such active flow atic symmetry [38] characterized by the proliferation of circulation is able to substantially affect the overall dy- topological defects whose behavior strongly differs from namics of the system [52, 57] and its rheological proper- their passive counter-part [6, 39–45]. ties [15, 16, 18, 56]. However, experiments on concentrated suspensions of Our main objective in this Letter is to investigate the bacteria and microtubules show that polar and nematic possibility that the flows internally generated by polar order only emerges in presence of activity, while isotropic active constituents may induce orientational order even symmetry is restored in the passive limit [33, 46, 47], thus in absence of an ad hoc free-energy, analogously to what suggesting that the thermodynamic ground state of a pas- happens in the case of an external shear. A similar study sive suspension is basically disordered [48] and activity- has been carried out by Thampi et al. [58] for the case of induced visco-elastic couplings may alone drive a transi- a system with nematic symmetry, proving that activity arXiv:2010.14124v1 [cond-mat.soft] 27 Oct 2020 tion to an ordered state by effective aligning interactions alone is actually able to affect the degree of nematic or- among the individual constituents. This is somehow in dering. They found flow-driven macroscopic fluctuations contrast with the typical approach of many models for ac- of the concentration and phase-separation, resulting into tive fluids – as the well-known active gel theory [49, 50] – a dynamical mixture of extensile and contractile regions. where the existence of an equilibrium free-energy favour- We will show that an isotropic-ordered phase transition ing an ordered equilibrium configuration in absence of also develops in polar systems, exhibiting yet a different activity is postulated a priori [29, 51–54]. nature due to different symmetry features. By means of 2 numerical simulations, we will show that activity is able The evolution equation for the polarization field is to give rise to a first-order isotropic-polar (I-P) transi- given by the Ericksen-Leslie equation [2, 60] adapted for tion for extensile arrow-like systems and for contractile the treatment of a vector order parameter and reads as disk-like ones, in presence of nemato-hydrodynamic in- follows teractions. In addition, we perform a scale-to-scale anal- 1 ysis in Fourier space to investigate the intricate energy- dtP = ξD P + h, (4) · Γ transfer mechanism in our model, finding that the most relevant contribution to energy absorption is the polar where we denoted with dtP the sum of the material one, contrary to what reported in analogous investiga- derivative ∂tP + (v ) P and the co-rotational deriva- · ∇ tions [29–31, 59] on usual active gel theory. tive Ω P . Here, Dαβ = (∂βvα + ∂αvβ)/2 and Ωαβ = · (∂βvα ∂αvβ)/2 are respectively the symmetric and anti- symmetric− parts of the strain tensor v. The alignment MODEL parameter ξ controls the aspect ratio∇ of the suspended particles (ξ > 0 for rod-like particles and ξ < 0 for disk- We consider a polar fluid in a 2d geometry with mass line ones) as well as the response to a shear flow. In density ρ and velocity . The ordering properties of the particular, when ξ > 1, the polarization field uniformly | | suspended particles are encoded in the polarization field orients at the Leslie angle with respect to the imposed P , while the concentration of the nutrient is denoted with flow [55, 62]. In this case, the is said to be φ. The equations which rules the hydrodynamic of the flow-aligning. Conversely, when ξ < 1, the liquid crys- | | system are tal rotates under the effect of shear and is customarily addressed as flow tumbling. The rotational viscosity Γ v = 0 (1) measures the relevance of the hydrodynamic interactions ∇ · ρ(∂tv + v v) = p + σ (2) with respect to the driving effect of the molecular field · ∇ −∇ ∇ · h = δ /δP , where is a suitable free energy which where the first one is the solenoidal condition for the flow we introduce− F to guaranteeF the thermodynamic stability field, enforcing the incompressibility of the system, and of the system. To do so, we consider an expansion in the second one is the Navier-Stokes equation. Here, p is power of P 2 up to sixth order: the hydrodynamic pressure, while σ is the stress tensor. | | Z  k  This is in turn the sum of three contributions. The first = dr α P 2 + β P 4 + γ P 6 + P ( P )2 (5) is the usual viscous stress F | | | | | | 2 ∇

visc where kP is the elastic constant of the polar phase and α, σαβ = η(∂αvβ + ∂βvα), β and γ are three bulk phenomenological parameters with where η is the nominal viscosity of the fluid and it is γ > 0 to ensure stability. When both α and β are positive responsible for energy dissipation. The polar stress [60] or null (see continuous orange line in Fig. 1, showing the possible behaviors of the bulk potential relevant for pol 1 ξ σ = (Pαhβ Pβhα) (Pαhβ + Pβhα) the following discussion), the only stable solution is the αβ 2 2 − − trivial one at P = 0, corresponding to the isotropic phase. kP ∂αPγ ∂βPγ , − This condition remains unaltered even for positive values provides a non-linear coupling between the dynamics of of α, as long as β > β∗∗ = √3αγ (see continuous purple − the fluid and the evolution of the polarization field. The line in Fig. 1), where β∗∗ is a critical value corresponding explicit structure of the molecular field h and the role to the appearance of a meta-stable ordered region at P = 6 of the flow-aligning parameter ξ will be discussed in the 0 (yellow line in Fig. 1). Finally for β < βcr = √4αγ − following. Finally, the active stress tensor the minimum corresponding to the polar phase becomes absolute.  2  act P In the usual active gel theory, the parameters are usu- σαβ = ζφ PαPβ | | δαβ − − 2 ally chosen to ensure the free-energy to have a global minimum corresponding to non-null polarization. For is a phenomenological term [61] responsible for the lo- the case of a polar gel, this can be obtained by mak- cal energy input due to the active constituents whose ing use of a double-well potential by setting α < 0, rate is proportional to the activity parameter ζ –positive β > 0 and γ = 0 [51]. Conversely, here we will only for extensile swimmers and negative for contractile ones– consider cases with α > 0 and β > β∗∗, whose equi- and to the nutrient concentration φ which in turn evolves librium configuration is the isotropic one. However, it through an advection-diffusion equation turns out that visco-elastic interactions arising due to 2 activity-induced flows, may lead to a significant renor- ∂tφ + v φ = D φ, (3) · ∇ ∇ malization of the bulk coefficients, leading to a situation where D is the diffusion constant. where polar order is also observed. To show this, we start 3

a 1.0 0.004 b 1.0 , α = 0 β = 0 0.8 α = 0, β = 0, ζ > 0 0.8

F α > 0, β = β∗∗ ζ =0.014 0.6 0.6 0.003 , ζ =0.016 α > 0 β > β∗∗ ζ =0.018 α > 0, β > β∗∗, ζ > 0 ζ =0.024 0.4 0.4 ζ =0.030 Polar Fraction

0.002 0.2 0.2

0.0 0.0 0.0 0.2 0.4 0.6 0.8 1.0 1 2 3 4 t 106 0.001 × 101 c 1.0 d × 10 10 0.8 1.5 0 8 8 Bulk Free Energy Density 0.6 6 Fluctuations 6 Fluctuations Mean Mean 1.0 0.001 0.4 4 4 Polar Fraction − 0 0.2 0.4 0.6 0.8 1 1.2 Kinetic Energy Polarization Magnitude ~ 0.2 2 2 0.5 Polar Fraction Fluctions (%)

P Kinetic Energy Fluctions (%) | | 0 0 0.0 0 1 2 0.0 0.5 1.0 1.5 2.0 (ζ ζ )/ζ (ζ ζ )/ζ Figure 1. Effective free energy in active polar gels. − cr cr − cr cr Continuous lines show the functional dependence of the free energy from the polarization field. For positive α and γ and Figure 2. Uniform nutrient concentration. (a) Velocity β > β∗∗, the only stable solution is P = 0 corresponding field superimposed to the contour plot of the magnitude of 2 to the isotropic phase, while at β = β∗∗ the free-energy ex- P on a 64 fraction of the system for ξ = 1.1 and ζ = 0.015 hibit an inflection point at P =6 0. Dashed lines represent ((ζ−ζcr)/ζcr = 0.5). (b) Polar fraction temporal evolution for the corresponding free energies after the renormalization of some values of the activity ζ. Panels (c) and (d) respectively the coefficient of the quartic term β −→ β˜, in accordance to show the behavior of the polar fraction and kinetic energy Eq. (8). and their fluctuations, at varying ζ across the isotropic-polar transition (ζcr = 0.01). Temportal fluctuations are computed here and elsewhere as percent deviation from the mean. from the assumption –which we will verify a posteriori in numerical simulations– of local mechanical equilibrium σvisc + σpol + σact 0. Retaining only the lowest bulk lines in Fig. 1. This may either occur for extensile rod- order contributions≈ to the molecular field ( P ) allows like particles or contractile disk-like ones. In particular, us to find a functional relation between the flow∼ structure we observe that the predicted nature of the transition is ∗∗ and the polarization first order, hence in the range βcr < β˜(ζ, ξ) < β we expect phase coexistence of polar (active) and isotropic ζφ  P 2  αξ (passive) regions. The transition becomes second order Dαβ PαPβ | | δαβ PαPβ. (6) ≈ 2η − 2 − η in the limit α, β = 0.

By substituting this relation into Eq. (4) we obtain:   UNIFORM NUTRIENT CONCENTRATION ξ ζφ 2 1 dtP = αξ P P + h. (7) η 4 − | | Γ We numerically integrate the hydrodynamic equations This new term has cubic dependence and induce a renor- Eqs. (1)-(2) on a squared computational grid of size malization of the coefficient of the quartic term in Eq. (5) L = 256 by means of a predictor-corrector LB model as follows with periodic boundary conditions. The dynamics of the order parameter φ and P is solved through a predictor- ξΓ ζφ  β β˜ = β αξ , (8) corrector finite-difference algorithm, implementing first- −→ − 4η 4 − order upwind scheme and fourth order accurate stencils for space derivatives [50, 63]. Except otherwise stated, so that for strong enough active forcing we set η = 5/3,D = 4 10−4 and ξ = 1.1, while the × 2 polarization field is initialized with random orientation 4αξ 16η  p  ζξ > + β + Θαγ , (9) and . Moreover we choose φ φΓ P = 1 α = γ = 0.005 β = | 0|.0065 > β∗∗, deep in the isotropic region of the the renormalized free energy may either develop a local parameter− space (see Fig. 1, continuous purple line). minimum at P = 0, (for Θ = 3 corresponding to the We first consider the case of an active polar gel with critical condition6 β˜ < β∗∗) or a global one (for Θ = 4 uniform nutrient concentration (i.e. φ = 1 in the whole ˜ corresponding to β < βcr), as suggested by the dashed volume and constant in time). For 0 < ζ . ζcr = 0.01, 4

2 10− a b × 1.5 Active injection energy as shown in Fig. 2(d), and invade the whole com- Viscous dissipation Polar absorption Hydrodynamic transfer putational domain, pushing the system into a polar state 1.0 with ψP & 0.9 and negligible fluctuations. This state, 0.5 characterized by vortical regions both in the velocity and the polarization field, can be considered the analogue of 0.0 the active turbulent regime for the usual active gel theory. Power spectrum 0.5 Hence, our results suggest that activity alone is able to − drive a transition from a quiescent isotropic state at low 1.0 − activity towards a dynamical polar state where rotational 100 101 102 k symmetry is restored at a statistical level due to the chaotic vortical patterning of the polarization field. Such Figure 3. Dynamical nutrient concentration. Panel (a) transition develops in an extended range of the activity shows a typical non-equilibrium steady-state configuration of parameters (0.01 . ζ . 0.02) and exhibits the typical the polarization field (colored vectors) superimposed to the features of a first-order transition, including coexistence contour plot of the concentration fluctuations δφ = φ − φ0, between ordered and disordered regions and large fluctu- for ξ = 1.1 and ζ = 0.015. Notice that polar/isotropic regions correspond to negative/positive values of δφ. Panel (b) shows ations. These values are of the same order of magnitude that the energy balance in Fourier space (see Eq. (10)) satisfies and roughly in agreement with those obtained through the hypothesis of local mechanical equilibrium. Eq. (9) which predicts the polar metastability region in the range 0.020 . ζ . 0.022 (1 < (ζ ζcr)/ζcr < 1.2). A further increment of the activity eventually− leads the sys- the rate of active injection is not strong enough to affect tem in a region of metastability for the isotropic phase the relaxation towards the ground state and the polariza- where the effective minimum for P is favored and the tion rapidly disappears, leaving the system in an isotropic polar fraction saturates. motion-less state. As activity is increased over the critical threshold ζcr, a dynamical non equilibrium state sets up, character- DYNAMICAL NUTRIENT CONCENTRATION ized by the coexistence of disordered/passive regions and polar/active domains with elongated shape, where is P We shall now treat the concentration of the nutrient φ consistently different from 0 (see for instance panels (a) as a dynamical field and let it evolve from its uniform ini- and (b) of Fig. 2). These evolve in time under the fu- tial state ( φ =1) according to Eq. (3). TFor small active eling action of the activity and continuously shrink and doping theh systemi ends up in an isotropic quiescent state, enlarge, giving rise to large fluctuations of the mean polar while for ζ > ζcr = 0.012 fluctuating polar domains begin fraction ψP of the system (Fig. 2(c)) computed as to populate the system. Their dynamics is accompanied Z by fluctuations in the concentration field caused by the 1 ¯ ψP (t) = ϑ( P P )dr advective effect of active flows. Interestingly, tends to L2 h | | − i φ lower in those regions populated by the polarization (see where ϑ( ) is the Heaviside step function and P¯ = 0.2 Fig. 3). This can be easily understood in terms of the is a suitable· threshold to distinguish polar regions from relation between P and the flow field. Indeed Eq. (6) isotropic ones and ... denotes the time average. More suggests that ∂αvβ (ζφ 2αξ)PαPβ/η hence, when the specifically, regionsh of welli aligned liquid crystal undergo condition in Eq. (9)∼ is fulfilled,− we get from Eq. (3) that 2 a banding instability. In the usual active gel theory this ∂tφ φP < 0. Therefore, the concentration field is mechanism is at the the base of the production of topo- advected∼ − towards isotropic regions, leading to an effec- logical defects [42, 44]. However, in our model, strong tive increment in the rate of active injection which favors distortions of the polarization pattern are not sustainable in turn the formation of polar regions. This is at the since they are associated with a reduction of the polar or- base of the observed φ and P patterns in ours simula- der, which results in the thinning of the polar domains tions, where thin stripes of polarization are interspaced and eventually in their disappearance, due to the diffu- with isotropic regions (see Fig. 3a). The presence of a 2 sive elastic term kP P in the molecular field h. dynamical concentration field has the important effect Nevertheless, these∼ − regions∇ of deformation also act as to reduce the temporal fluctuations in the polar fraction act a source of momentum –since the active force f ψP at the isotropic/polar transition (see comparison in PP – and produce jets of fluid which advect and shear∝ Fig. 4c) which now develops in a more narrow range of ac- ∇· the liquid crystal, counteracting elastic effects and re- tivities, signalled by a sharp increase of ψP in correspon- inforcing the polar order. This mechanism rapidly be- dence of the critical activity ζcr (see red lines in Fig. 2c comes more and more dominant as extensile activity and Fig.4a). is increased over 2ζcr = 0.02. Activity-induced flows A feature that can enlighten the dynamical role of po- strengthen, causing a sharp increase of the total kinetic larization comes from the analysis of different contribu- 5

1 a 10− tions to energy transfer. We shall now consider a scale- × b 2.5 to-scale balance equation for the kinetic energy in Fourier 8 space: 2.0 6 1.5 visc pol act extensile extensile ∂tEk + k = + + , (10) contractile contractile k k k 4 T S S S 1.0 Polar Fraction where the energy spectrum v 2 and Kinetic Energy Ek = k Ω/2 ... Ω 2 0.5 stands for the average in Fourier spaceh| over| i shells ofh equali ∗ 0 0.0 momentum ( k = k). k = vk Jk Ω represents the 0 1 2 0 1 2 | | T h · i (ζ ζ )/ζ (ζ ζ )/ζ advective transfer rate, being Jk the hydrodynamic flux − cr cr − cr cr 101 c d 1.5 1.0 p+v v in the k-space. The terms k are source/sink × −∇ ·∇ S 1.0 terms for the kinetic energy. In particular, 1.0 0.8

0.8 Fraction Polar 2 0.5 visc ∗ visc 8π η 2 0.6 0.6 Dynamical nutrient concentration ξ k = 2πi vk k : σk Ω/L = 2 k Ek 0.0 S h ⊗ i − L Uniform nutrient concentration 0.4 0.4 0.5 gives the rate at which energy is dissipated by viscosity − 0.2 0.2 Polar Fraction Fluctions 1.0 and − 0.0 1.5 0.0 pol ∗ pol 0 1 2 − 5 0 5 − 2 k = 2πi vk k : σk Ω/L (ζ ζcr)/ζcr ζ 10− S h ⊗ i − × gives the rate at which the polarization field absorbs en- Figure 4. Polar phase diagram. Panels (a) and (b) re- ergy to be maintained in its non-equilibrium state. Fi- spectively show the behavior of the polar fraction and the nally, mean kinetic energy at varying activity for an extensile sys- tem with and a contractile one with − . The act ∗ act ξ = 1.1 ξ = 1.1 k = 2πi vk k : σk Ω/L critical extensile (contractile) threshold is . S h ⊗ i ζcr = 0.012 (0.015) Panel (c) shows the comparison of the time fluctuations of the gives the rate at which energy is injected in the system polar fraction between the case with uniform and dynamical by the activity. nutrient concentration. Panel (d) shows the contour plot of Panel (b) of Fig. 3 shows the energy balance for the the polar fraction in the ξ − ζ plane. The dashed yellow line case at ζ = 0.03. As expected, the active power spec- represent the theoretically predicted transition line of Eq. (9). act trum injects energy over all scales ( k > 0) and is counter-balanced by viscous dissipationS and polarization absorption, while the non-linear hydrodynamic advective state for any value of activity, proving that the align- term , in accordance to the value of the Reynolds ment term D P is crucial to the development of the k 0 ∼ · numberT Re≈ = ρl∗v∗/η (where l∗ and v∗ are respectively nemato-hydrodynamic instability discussed in this paper. a typical scale and velocity of the flow [59]) which never By considering progressively higher values of ξ we find | | exceeds 10−1 in our simulations. We found that the polar that the transition towards the ordered state develops at term gives here a greater contribution to energy absorp- smaller and smaller values of ζ . This mechanism dras- | | tion than viscosity. This is in contrast with respect to tically depends on the sign of both ξ and ζ, as suggested previous studies which accounted for the energy transfer by Eq. (9). Indeed, the phase-diagram in Fig. 4d shows in active systems, where the polar term could be factu- that polar regions are only observed for either extensile ally neglected [29–31, 59]. The explanation lies in the rod-like particles (ζ > 0, ξ > 0) or contractile disk-like fact that while in the usual active gel approach, the po- ones (ζ < 0, ξ < 0). The observed position of the tran- larization field can be assumed close to equilibrium with sition line in the ξ ζ plane is in good agreement with − the only relevant contribution arising from deformations the theoretical prediction of Eq. (9) (see dashed line in 2 (h kP P ) in the present case activity must provide Fig. 4d). Moreover, the response of the system is roughly enough∼ − energy∇ to overcome the energy barrier due to bulk symmetric under a sign change in both the aligning pa- terms, which would drive the system into an isotropic rameter and the activity, as suggested by the comparisons state in absence of active forcing. between the polar fraction and the kinetic energy for the Power spectra behavior exhibits a localized balance in extensile and contractile case at ξ = 1.1 respectively | | Fourier space, in agreement with our assumption of lo- shown in Fig. 4(a) and (b). cal mechanical equilibrium. Moreover, the most relevant scales λ, at which the rate of active injection reaches its maximum for the case in Fig. 3, range in 64 < λ < 13 CONCLUSION (with λ = L/k) corresponding to the typical size of the structures developed by the polarization field. In this Letter we have investigated the dynamics of Finally, we comment on the effect of the alignment an active polar gel in absence of a free-energy favoring parameter ξ. At ξ = 0 we do not observe a polar a polar state, contrary to most previous work based on 6 active gel theory where an ordered phase is obtained by [7] P. Guillamat, Ž. Kos, J. Hardoüin, J. Ignés-Mullol, opportunely setting the parameters of some free-energy M. Ravnik, and F. Sagués. Taming active turbulence functional. We found that nemato-hydrodynamic cou- with patterned soft interfaces. Nat. Comm., 8, 12 2017. plings may induce a transition from the motion-less and [8] S. Ebbens, D.A. Gregory, G. Dunderdale, J.R. Howse, Y. Ibrahim, T.B. Liverpool, and R. Golestanian. Elec- energetically favored disordered phase towards a self- trokinetic effects in catalytic platinum-insulator janus sustained flowing polar state evolving far from equilib- swimmers. Europhys. Lett., 106:5, 2014. rium. In our model the rate of energy absorption due [9] P.A. Korevaar, C.N. Kaplan, A. Grinthal, R.M. Rust, to the polarization cannot be neglected –as we showed and J. Aizenberg. Non-equilibrium signal integration in through a scale-to-scale analysis in Fourier space– con- hydrogels. Nat. Comm., 11:386, 2020. trary to what reported by previous works on energy [10] H. Wioland, F.G. Woodhouse, J. Dunkel, J.O. Kessler, transfer in active gels. On these assumptions, we pro- and R.E. Goldstein. Confinement stabilizes a bacte- rial suspension into a spiral vortex. Phys. Rev. Lett., posed an analytical argument which suggests such transi- 110:268102, Jun 2013. tion to be first-order. This picture was confirmed through [11] J. Dunkel, S. Heidenreich, K. Drescher, H.H. Wensink, numerical experiments which showed that coexistence of M. Bär, and R.E. Goldstein. of bacterial polar and isotropic regions occurs over a large range of turbulence. Phys. Rev. Lett., 110:228102, 2013. activities at the transition, accompanied by important [12] Y. Hatwalne, S. Ramaswamy, M. Rao, and R.A. Simha. temporal fluctuations in the polar fraction of the sys- of active-particle suspensions. Phys. Rev. Lett., tem. Moreover, in presence of a diffusive field describing 92:118101, 2004. [13] T.B. Liverpool and M.C. Marchetti. Rheology of active the concentration of the nutrient, the transition becomes filament solutions. Phys. Rev. Lett., 97:268101, 2006. sharper and fluctuations are considerably reduced. [14] Salima Rafaï, Levan Jibuti, and Philippe Peyla. Effective We believe that results of our paper provide a solid viscosity of microswimmer suspensions. Phys. Rev. Lett., ground for a deeper understanding of the motility of ac- 104:098102, 2010. tive gels and furnish a starting point to study the effects [15] S. Guo, D. Samanta, Y. Peng, X. Xu, and X. Cheng. of confinement. In this case, the presence of boundaries Symmetric shear banding and swarming vortices in bac- either soft or rigid may result into dynamical states which terial superfluids. Proc. Natl. Acad. Sci., 2018. [16] G. Negro, L.N. Carenza, A. Lamura, A. Tiribocchi, and are commonly observed in experiments [33]. G. Gonnella. Rheology of active polar emulsions: from We warmly thank Julia Yeomans and Davide Maren- linear to unidirectional and unviscid flow, and intermit- duzzo for discussions from which this work originated. tent viscosity. Soft Matter, 15:8251–8265, 2019. The work has been performed under the Project HPC- [17] Jérémie Gachelin, Gastón Miño, Hélène Berthet, Anke EUROPA3 (INFRAIA-2016-1-730897), with the support Lindner, Annie Rousselet, and Éric Clément. Non- of the EC Research Innovation Action under the H2020 newtonian viscosity of escherichia coli suspensions. Phys. Programme; in particular, the authors gratefully ac- Rev. Lett., 110:268103, 2013. [18] H.M. López, J. Gachelin, C. Douarche, H. Auradou, and knowledges the support and the computer resources pro- E. Clément. Turning Bacteria Suspensions into Superflu- vided by EPCC. ids. Phys. Rev. Lett., 115:028301, 2015. [19] S. Sankararaman and S. Ramaswamy. Instabilities and waves in thin films of living fluids. Phys. Rev. Lett., 102:118107, 2009. [20] R. Voituriez, J.F. Joanny, and J. Prost. Spontaneous flow [1] S. Ramaswamy. The mechanics and statistics of active transition in active polar gels. E.P.L., 70:404, 2005. matter. Annu. Rev. Condens. Matter Phys., 1:323, 2010. [21] F.C. Keber, E. Loiseau, T. Sanchez, S.J. DeCamp, [2] M.C. Marchetti, J.F. Joanny, S. Ramaswamy, T.B. Liv- L. Giomi, M.J. Bowick, M.C. Marchetti, Z. Dogic, and erpool, J.J. Prost, M. Rao, and R.A. Simha. Hydrody- A.R. Bausch. Topology and dynamics of active nematic namics of soft active matter. Rev. Mod. Phys., 85:1143, vesicles. Science, 345(6201):1135–1139, 2014. 2013. [22] L.N. Carenza, G. Gonnella, D. Marenduzzo, and G. Ne- [3] H.H. Wensink, J. Dunkel, S. Heidenreich, K. Drescher, gro. Rotation and propulsion in 3d active chiral droplets. H. Lowen R.E. Goldstein, and J.M. Yeomans. Meso-scale Proc. Natl. Acad. Sci., 116(44):22065–22070, 2019. turbulence in living fluids. Proc. Natl. Acad. Sci., 109, [23] L.N. Carenza, G. Gonnella, D. Marenduzzo, and G. Ne- 2012. gro. Chaotic and periodical dynamics of active chiral [4] V. Bratanov, F. Jenko, and E. Frey. New class of droplets. Physica A, 559:125025, 2020. turbulence in active fluids. Proc. Natl. Acad. Sci., [24] L. Giomi, M.C. Marchetti, and T.B Liverpool. Complex 112(49):15048–15053, 2015. spontaneous flows and concentration banding in active [5] P.M. Bendix, G.H. Koenderink, D. Cuvelier, Z. Dogic, polar films. Phys. Rev. Lett., 101:198101, 2008. B.N. Koeleman, W.M. Brieher, C.M. Field, L. Mahade- [25] E. Tjhung, D. Marenduzzo, and M.E. Cates. Spon- van, and D.A. Weitz. A quantitative analysis of contrac- taneous symmetry breaking in active droplets provides tility in active cytoskeletal protein networks. Biophys. J., a generic route to motility. Proc. Natl. Acad. Sci., 94:117960, 2008. 109(31):12381–12386, 2012. [6] S.J. DeCamp, G.S. Redner, A. Baskaran, M.F. Hagan, [26] L. Giomi and A. DeSimone. Spontaneous division and and Z. Dogic. Orientational order of motile defects in motility in active nematic droplets. Phys. Rev. Lett., active nematics. Nat. Mater., 14:11110, 2015. 112:147802, 2014. 7

[27] F. Bonelli, L.N. Carenza, G. Gonnella, D. Marenduzzo, Self-organized dynamics and the transition to turbulence E. Orlandini, and A. Tiribocchi. Lamellar ordering, of confined active nematics. Proc. Nat. Acad. Sci., droplet formation and phase inversion in exotic active 116(11):4788–4797, 2019. emulsions. Sci. Rep., 9, 2019. [47] E. Lushi, H. Wioland, and R.E. Goldstein. Fluid flows [28] C.W. Wolgemuth. Collective swimming and the dynam- created by swimming bacteria drive self-organization ics of bacterial turbulence. Biophys. J., 95(4):1564 – in confined suspensions. Proc. Nat. Acad. Sci., 1574, 2008. 111(27):9733–9738, 2014. [29] L. Giomi. Geometry and topology of turbulence in active [48] É. Fodor, C. Nardini, M.E. Cates, J. Tailleur, P. Visco, nematics. Physical Review X, 5, 2015. and F. van Wijland. How far from equilibrium is active [30] J. Urzay, A. Doostmohammadi, and J.M. Yeomans. matter? Phys. Rev. Lett., 117:038103, 2016. Multi-scale statistics of turbulence motorized by active [49] W. Kung, M.C. Marchetti, and K. Saunders. Hydrody- matter. Journal of Fluid Mechanics, 822, 07 2017. namics of polar liquid crystals. Phys. Rev. E, 73:031708, [31] R. Alert, J.-F. Joanny, and J. Casademunt. Universal 2006. scaling of active nematic turbulence. Nat. Phys., 3 2020. [50] L.N. Carenza, G. Gonnella, A. Lamura, G. Negro, and [32] L.N. Carenza, L. Biferale, and G. Gonnella. Multiscale A. Tiribocchi. Lattice boltzmann methods and active control of active emulsion dynamics. Phys. Rev. Fluids, fluids. Eur. Phys. J. E, 42(6):81, 2019. 5:011302, Jan 2020. [51] K. Kruse, J.F. Joanny, F. Jülicher, J. Prost, and K. Seki- [33] T. Sanchez, D.T.N. Chen, S.J. Decamp, M. Heymann, moto. Asters, vortices, and rotating spirals in active gels and Z. Dogic. Spontaneous motion in hierarchically as- of polar filaments. Phys. Rev. Lett., 92:078101, 2004. sembled active matter. Nature, 491:431–434, 2012. [52] M.E. Cates, S.M. Fielding an D. Marenduzzo, E. Or- [34] F. Nédélec, T. Surrey, A.C. Maggs, and S. Leibler. landini, and J.M. Yeomans. Shearing active gels close Self-organization of microtubules and motors. Nature, to the isotropic-nematic transition. Phys. Rev. Lett., 389:305, 1997. 101:068102, 2008. [35] D. Saintillan and M. J. Shelley. Orientational order and [53] M.L Blow, S.P. Thampi, and J.M Yeomans. Biphasic, instabilities in suspensions of self-locomoting rods. Phys. lyotropic, active nematics. Phys. Rev. Lett., 113:248303, Rev. Lett, 99:058102, 2007. 2014. [36] T. Vicsek, A. Czirók, E. Ben-Jacob, I. Cohen, and [54] T.N. Shendruk, A. Doostmohammadi, K. Thijssen, and O. Shochet. Novel type of phase transition in a system J.M. Yeomans. Dancing disclinations in confined active of self-driven particles. Phys. Rev. Lett., 75:1226, 1995. nematics. Soft Matter, 2017. [37] A. Baskaran and M.C. Marchetti. Enhanced diffusion [55] T. Markovich, E. Tjhung, and M.E. Cates. Shear-induced and ordering of self-propelled rods. Phys. Rev. Lett., first-order transition in polar liquid crystals. Phys. Rev. 101:268101, 2008. Lett., 122:088004, Feb 2019. [38] L.F. Cugliandolo, P. Digregorio, G. Gonnella, and [56] L.N. Carenza, G. Gonnella, A. Lamura, D. Marenduzzo, A. Suma. Phase coexistence in two-dimensional pas- G. Negro, and A. Tiribocchi. Soft channel formation and sive and active dumbbell systems. Phys. Rev. Lett., symmetry breaking in exotic active emulsions. Sci. Rep., 119:268002, 2017. 10, 2020. [39] G. Tóth, C. Denniston, and J.M Yeomans. Hydrodynam- [57] L. Giomi, T.B. Liverpool, and M.C. Marchetti. Sheared ics of topological defects in nematic liquid crystals. Phys. active fluids: Thickening, thinning, and vanishing viscos- Rev. Lett., 88:105504, 2002. ity. Phys. Rev. E, 81:051908, 2010. [40] S. Zhou, A. Sokolov, O. Lavrentovich, and I.S. Aran- [58] P.T. Sumesh, A. Doostmohammadi, R. Golestanian, and son. Living liquid crystals. Proc. Natl. Acad. Sci., J.M. Yeomans. Intrinsic free energy in active nematics. 111(4):1265–1270, 2014. EPL, 112(2):28004, 2015. [41] E. Putzig, G.S. Redner, A. Baskaran, and A. Baskaran. [59] L.N. Carenza, L. Biferale, and G. Gonnella. Cascade Instabilities, defects, and defect ordering in an over- or not cascade? Energy Transfer and Elastic Effects in damped active nematic. Soft Matter, 12:3854–3859, 2016. Active Nematics. Under Review, 2020. [42] Luca Giomi, Mark J. Bowick, Xu Ma, and M. Cristina [60] P.G. de Gennes and J. Prost. The physics of liquid crys- Marchetti. Defect annihilation and proliferation in active tals. The International series of monographs on physics nematics. Phys. Rev. Lett., 110:228101, 2013. 83 Oxford science publications. , [43] J. Elgeti, M.E. Cates, and D. Marenduzzo. Defect hydro- 2nd ed edition, 1993. dynamics in 2d polar active fluids. Soft Matter, 7:3177– [61] R.A. Simha and S. Ramaswamy. Hydrodynamic fluc- 3185, 2011. tuations and instabilities in ordered suspensions of self- [44] A. Doostmohammadi, T.N. Shendruk, K. Thijssen, and propelled particles. Phys. Rev. Lett., 89:058101, 2002. J.M. Yeomans. Onset of meso-scale turbulence in active [62] R.G. Larson. The structure and rheology of complex flu- nematics. Nat. Comm., 8, 2017. ids. Topics in chemical engineering. Oxford University [45] G. Negro, L.N. Carenza, P. Digregorio, G. Gonnella, and Press, 1999. A. Lamura. Morphology and flow patterns in highly [63] C. Denniston, E. Orlandini, and J.M. Yeomans. Lattice asymmetric active emulsions. Physica A, 503:464 – 475, boltzmann simulations of liquid crystal hydrodynamics. 2018. Phys. Rev. E, 63:056702, 2001. [46] A. Opathalage, M.M. Norton, M.P.N. Juniper, B. Langeslay, S.A. Aghvami, S. Fraden, and Z. Dogic.