![Arxiv:2002.12893V3 [Cond-Mat.Soft] 7 Sep 2020](http://data.docslib.org/img/7783fbcd6b3fa33e2874464a6f65f6a8-1.webp)
Emergence of dynamic vortex glasses in disordered polar active fluids
Am´elieChardac,1 Suraj Shankar,2 M. Cristina Marchetti,3 and Denis Bartolo1 1Univ. Lyon, ENS de Lyon, Univ. Claude Bernard, CNRS, Laboratoire de Physique, F-69342, Lyon, France 2Department of Physics, Harvard University, Cambridge, MA 02318, USA 3Department of Physics, University of California Santa Barbara, Santa Barbara, CA 93106, USA (Dated: September 8, 2020) In equilibrium, disorder conspires with topological defects to redefine the ordered states of matter in systems as diverse as crystals, superconductors and liquid crystals. Far from equilibrium, however, the consequences of quenched disorder on active condensed matter remain virtually uncharted. Here, we reveal a state of strongly disordered active matter with no counterparts in equilibrium: a dynamical vortex glass. Combining high-content microfluidic experiments and theory, we show how colloidal flocks collectively cruise through disorder without relaxing the topological singularities of their flows. The resulting state is highly dynamical but the flow patterns, shaped by a finite density of frozen vortices, are stationary and exponentially degenerated. Quenched isotropic disorder acts as a random gauge field turning active liquids into dynamical vortex glasses. We argue that this robust mechanism should shape the collective dynamics of a broad class of disordered active matter, from synthetic active nematics to collections of living cells exploring heterogeneous media.
From a physicist perspective, flocks are ensembles of of self-propelled particles can steadily cruise though dis- living or synthetic motile units collectively moving along order maintaining local orientational order and without a common emerging direction [1–4]. They realize one of relaxing the topological singularities of their flows. The the most robust ordered states of matter observed over associated flow patterns are exponentially degenerated five orders of magnitude in scale and in systems as di- and shaped by amorphous ensembles of frozen topologi- verse as motility assays, self-propelled colloids, shaken cal defects, yielding a dynamical state akin to the static grains and actual flocks of birds [3,5–10]. The quiet vortex-glass phase of dirty superconductors and random- flows of flocks are in stark contrast with the spatio- gauge magnets [27–29]. Building a theory of flock hy- temporal chaos consistently reported and predicted in drodynamics beyond the spin-wave approximation, we active nematic liquid crystals, another abundant form elucidate the emergence and stabilization of topological of ordered active matter realized in biological tissues, vortices by quenched disorder. Finally, we discuss the swimming cells, cellular extracts and shaken rods [2, 11]. universality of the dynamical vortex glass phase beyond Active nematics do not support any form of long range the specifics of polar active matter and colloidal flocks. order [4, 12]. Their structure is continuously bent and destroyed by the proliferation and annihilation of sin- gularities in their local orientation: topological defects DISORDERED FLOCKS ON A CHIP [11, 13–15]. Unlike in active nematics, topological de- fects in flocking matter are merely transient excitations Our experiments are based on the model system we in- which annihilate rapidly and allow uniaxial order to ex- troduced in [6]. We use Quincke rotation to power inani- tend over system-spanning scales [4]. mate polystyrene beads of radius a = 2.4 µm [30, 31], and This idyllic view of the ordered phases of active liq- turn them into colloidal rollers self-propelling at constant uids is limited, however, to pure systems. Disorder is speed 1 mm/s in circular microfluidic chambers of ra- known to profoundly alter the stability of topological de- dius R∼= 1.5 mm, see Fig.1A and Movie S1. Before fects and the corresponding ordered states in equilibrium addressing the impact of disorder on their collective dy- condensed matter [16–18], but its role in active fluids re- namics, it is worth recalling the phenomenology of the mains virtually uncharted territory. All previous stud- pure ordered phase [32, 33]. arXiv:2002.12893v3 [cond-mat.soft] 7 Sep 2020 ies [19–26], including our own early experiments [22], As the rollers are set in motion, they interact and self- have been limited to weak disorder and smooth per- assemble into a spontaneously flowing liquid equivalently turbations around topologically trivial states. Unlike referred to as a flock, a Toner–Tu fluid, or a polar liq- in equilibrium, disorder-induced topological excitations uid [2, 34] shown in Movie S1. The resulting flow pat- have been out of reach of any active-matter experiment, terns are initially isotropic and marred by a number of theory and simulations. 1 topological defects clearly visible in the Schlieren pat- In this article we show how isotropic disorder generi- terns± of Fig.1A and Movie S2. However, the velocity- cally challenges the extreme robustness of flocking mat- alignment interactions responsible for flocking motion pe- ter to topological defects. We map the full phase be- nalize flow distortions thereby causing the attraction and havior of colloidal flocks in disordered environments and annihilation of defects of opposite charge. The resulting reveal an unanticipated state of active matter: a dynam- lively coarsening dynamics lasts few tens of seconds and ical vortex glass. In dynamical vortex glasses, millions ultimately yields pristine azimuthal flows with long-range 2
a b c
Solvent inlet Outlet 0 1000 1500 2000 2500 3000 500 3000 0 Colloid inlet 0
3000 500 d 2500 2500 1000
2000 2000 1500
1500 1500 2000
1000
1000 2500
500
500 3000
0 0 500 1000 1500 2000 2500 3000 4 (%) 14 34
10 400400 600 400 200 200200 400 200 Time (s)
order in two dimensions, Fig.1A. ing the obstacle fraction Φo from 0 % to 38 % keeping the system below the Lorentz localization transition of To comprehensively investigate how disorder alters individual rollers [35, 36]. All experiments are system- the phases of polar active matter, we perform high- atically performed for several disorder realizations and content experiments on the microfluidic device sketched initial conditions, see also Methods. in Fig.1B. On a single chip, we quantify the flows of disordered roller fluids in 32 different disordered geome- Let us first discuss the impact of disorder on a polar tries. Disorder is implemented in 3 mm wide circular liquid deep in the flocking phase, for a colloid number 2 chambers decorated by random arrays of isotropic ob- density ρ0 = 7 0.5 10− . For small disorder, Φo < 9 %, stacles of diameter 10 µm that repel the rollers at a dis- the obstacles± hardly× disrupt the collective flow, as 1 tance, while leaving their speed unaltered, see Fig.1C defects formed at the onset of collective motion quickly± and Movie S3. We replicate the experiments increas- annihilate. The system coarsens to an ordered flocking 3 state which, in the disk geometry, takes the form of a nition of is natural in a circular geometry where the macroscopic vortex centered around a single +1 defect, spatial averageP of p vanishes, even in pure systems. Be- see Fig.1 and Movie S4. In stark contrast, increasing low Φc, Supplementary Fig. S1 clearly indicates that Φo above 9 %, a very distinct type of organization (`) converges to a finite value over a finite length scale emerges. Disorder essentially arrests the coarsening by ξP , signaling long range polar order. Conversely, deep in pinning 1 defects. While the global flow vanishes, theP meandering regime (`) vanishes over a finite scale. we find± locally correlated flows through meandering Polar liquids and meandersP are therefore two genuinely patterns bent by disorder (Fig.1D and Movies S5 and distinct phases of active matter. S6). Particles stream coherently through the meanders, The transition between these two dynamical phases is resulting in a finite correlation length for orientational captured by the global polarization (` = 2R), order, but the structure of the flow network itself is where R is the chamber radius. We canP ≡ see P in Fig.2B static and fixed in space, as shown in Fig.1D and Movie that decays weakly and linearly with Φ for small disor- P o S6. Further increasing the obstacle fraction (above der, but drops sharply to 0 at the critical value Φc defined 26 %), we recover the liquid-gas transition reported from the proliferation of quenched topological defects in in [22]. In the gas phase, particles are motile, but there Fig1F. The bifurcation of in Fig.2B suggests a crit- is no orientational order as defect pairs continuously ical scenario. This hypothesisP is further supported by unbind and annihilate as seen in Figs.1D,1E. The Fig.2C which reveals a divergence, or at the very least two transitions between the three dynamical states are a drastic increase of the correlation length ξ at the on- unambiguously determined by the variations of (Φ ), set of meandering motion. The polarizationP alone, N o the mean number of topological defects in the velocity however, does not distinguish the meander fromP the gas field, plotted in Fig.1F (the detection of the topological phase. We thus introduce the Edwards–Anderson pa- defects is detailed in Supplementary Note 3). The rameter EA = vˆ(r, t) vˆ(r, t + T ) r,t,T that quanti- emergence of meandering flows is signaled by the linear fies the temporalQ h persistence· of thei emergent→∞ flows [40]. increase of as Φ exceeds the critical value Φ = 9 %, The finite non-zero value of for Φ < Φ < 26 % N o c EA c o whereas the loss of local orientational order saturates (Fig.2B) confirms the persistenceQ of polar order along (Φ ) above Φ = 26 %. the meanders. We also find that the continuous tran- N o o sition, or crossover, between the meander and the gas phase, where EA 0, coincides with that identified from the topologicalQ ∼ defect statistics in Fig.1F. PHASE BEHAVIOUR OF DISORDERED POLAR Taking advantage of our high-content microfluidic ACTIVE MATTER experiments, we measure and EA in more than 400 experiments performed varyingP theQ roller fraction from 3 1 10− to 1.5 10− , and Φ from 0 to 38 % for multiple In order to quantitatively distinguish the dynamics in × o the three regimes, we inspect the statistics of the time- disorder realizations and initial conditions. The resulting phase diagram shown in Fig.2D firmly establishes that averaged polarization defined as p(r) = vˆ(r, t) t, where the unit vector vˆ(r, t) is the instantaneoush orientationi disorder and topological defects do not merely offset of the velocity field, see Methods. The distributions of the flocking transition as reported in [22], but conspire p(r) in Fig.2A show that, for small disorder, the flow is to bend the dynamics of polar liquids into amorphous uniformly polarized along the azimuthal direction, which and singular flows while preserving local flocking motion. reflects nearly perfect polar order. Conversely, above Φc global polar order is suppressed, and the polarization is isotropically distributed on the scale of the system. The strong localization of the distribution on the unit circle MEANDER FLOWS ARE DYNAMICAL VORTEX demonstrates, however, that meanders persistently dis- GLASSES tort the streamlines without arresting the local flows. Above Φo = 26 %, in the active-gas phase, the typical We now establish that the meanders realize a dynam- polarization vanishes, polar order melts and flows are ical vortex glass. To do so, we compare the flow ori- suppressed at all scales. entations of fifty replicas of the very same experiment We now demonstrate that the two flowing patterns cor- (Φo = 15 %). We solely vary the initial conditions, keep- respond to two distinct dynamical phases. We first stress ing the colloid fraction and disorder realization iden- that polar liquids with genuine long range orientational tical. We then quantify the resemblance between the order survive at finite disorder (see also Supplementary flow patterns by introducing a measure of the local over- Note 6). Unlike active nematics, deconfining the system lap qαβ(r) = pα(r) pβ(r) between replicas α and β. does not suppress the stationary vortex patterns shown Examples of overlap· maps are shown in Fig.3A for in Fig.1D[37–39]. This crucial result follows from the four different pairs of experiments, see also Supplemen- finite-size analysis of the polarization order parameter tary Fig. S5. The replicated flows are identical when p (`) = p (r) 2 + p (r) 2, where ` is the size of the q (r) = +1, opposite when q (r) = 1 and orthogo- P h r ir h θ ir αβ αβ − region where spatial averaging is performed. This defi- nal when qαβ(r) = 0. A simple inspection of the maps
3 4
0 1 -1 30
4
10 -3 30 aa -1 10 30 1 6 1 6511 6
1 6 0.1 9 d Polar liquid Meanders Gas 4
4 4 48 30
0.2 0 3 4
2 30 20 20.3 2 2
6 2 20 1 1 30
-1 -1 -1 0.4-1 -1 -1 5 30 0 00 0 0 -1 0 1 -1 0 1 -1-1 -1 0 0 11 1 20
FIG. 2. Suppression of long-range orientational order. a, PDF of the polarization fields P(pr,p✓) in the three dynamical phases state. Note thatFIG. the 2. distributionSuppression is isotropic of long-range but peaked along orientational the unit circle in order the meander. A, PDF phase: of spontaneousthe polarization flows fields P(pr, pθ) in the three dynamical are locally preserved butphases. global Note polar order that is the suppressed distribution by disorder. is isotropicb, Finite-size but analysis peaked of along the polarization the unit order circle parameter in the meander phase: spontaneous flows are (`) measured in square boxes of width `. The polarization decays exponentially towards a finite value in the polar-liquid phase P locally preserved but global polar order is suppressed by disorder. B, Variation of the global polarization and Edwards– ( o < c) demonstrating the existence of genuine long-range orientational order. No polar order is observed above c.The P chamber diameter 2RAndersonis 3 mm. c, VariationEA order of theparameters. global polarization In polarand liquids Edwards–Anderson> 0, in meandersEA order parameters.= 0 but TheEA > 0, and in a gas both and EA Q P P Q P Q P Q polar liquid phase isvanish. defined byC,> Correlation0, the meander length phase of by the=0, polarizationEA > 0andthegasphaseby field defined fromEA the=0. exponentiald, Correlation decay of (`) plotted in Supplementary P P Q Q length of the polarization field defined from the exponential decay of (`)assketchedin(b). The cusp of the ⇠ ( o) curve at P Figure S1. The cusp of the ξ (Φo) curveP at Φc hints towards a critical dynamicalP transition. The accuracy in the measurement c hints towards a critical dynamical transition. P 2 of ξ/2R is of the order of 10− . D, Phase diagram of polar active matter in isotropic disorder. The phase boundaries are equivalently determined from the number of topological defects (Fig.1F), or from the variations of and EA. In practice, we define the boundary separating the polar-liquid from the meander phase as > 0.45 and the boundaryP betweenQ the meander P distinct phases of activeand gasmatter. phases as < 0.45 and QourEA high-content> 0.2. The microfludic geometry experiments, of the phase we measure boundaries does not crucially depend on these criteria. The extentP of the meander phaseand becomesin 400 experimentsvanishingly performed small only at at XX the di↵ onseterentP of the flocking transition. The dashed The transition between these two dynamical phases is QEA captured by rectangle(` =2R), indicates where R theis the region chamber of the phaseroller densities, diagram and explored YY obstacle in [22] fractions where the for meander multiple phase was missed. Orientational order radius. We showP ⌘ inPis Fig. generically2c that lostdecays due weakly to the and emergencedisorder of realizationsmeander patterns. and initial Theconditions. gray rectangle The resulting corresponds to the series of experiments P linearly with o fordiscussed small disorder, in the butmain drops text. sharply phase diagram shown in Fig. 4a firmly establishes that to 0 at the critical value c defined from the prolifer- disorder and topological defects do not merely o↵set the ation of quenched topological defects in Fig 1e. The flocking transition as reported in [14] but generically con- bifurcation of in Fig. 2c suggests a critical scenario. spire to bend the polar liquid into an amorphous meander P phase where flocking motion is locally preserved. We This hypothesis isreadily further supported indicates by that, Fig. 2dwhich although all experiments corre- do not merely differ from each other by continuous dis- reveals a divergence, or at the very least a drastic in- now establish that the meanders realize a dynamical vor- crease of the correlationspond length to overall⇠ at differentthe onset of patterns, me- tex a few glass. macroscopicTo do so, we re-comparetortions. the flow orientations The replicas of are actually topologically inequiva- P andering motion.gions The polarization are virtuallyalone, identical however, from onefiftyreplica replicas ofto the another. very same experimentlent, as demonstrated ( o = 15 %). by the distribution of the difference does not distinguishMore the meander quantitatively, fromP the we gas plot phase. in Fig.We3C solely the vary distribution the initial conditions,in the keepingnumber the ofcol- topological-defect between each pair of We thus introduce the Edwards–Anderson parameter loid fraction and disorder realization identical. We then of the global overlap αβ = qαβ(r) r, and compare it to replicas ( α β). The regions of the flow patterns EA = vˆ(r,t) vˆ(r,t+ T ) r,t,T that quantifiesQ theh quantifyi the resemblance between the flowN patterns− N by Q h · that ofi the!1 macroscopic vortex obtainedintroducing in the the local flocking overlap q that(r)= arep ( robustr) p (r to)be- variations in the initial conditions are temporal persistence of the emergent flows [24]. The fi- ↵ ↵ · nite non-zero valuestate. of for The two< distributions< 26 % (Fig. 2c) are markedlytween the replicas different.↵ and In . Examplestherefore of overlapshaped maps by pinned defects of identical charge and QEA c o confirms the persistencethe flocking of polar order state, along P( theαβ mean-) is composedare shown of twoin Fig. symmet-4a for four diorientation↵erent pairs which of experi- persist in all replicas. In other words, Q ments, see also Supplementary Fig. S5. The replicated ders. We also findric that peaks the continuous at αβ = transition,1, as or global circulation of either defects are strongly correlated not only in their spatial crossover, between the meander andQ the gas phase,± where flows are identical when q↵ (r) = +1, opposite when handedness are equally probable.q Conversely,(r)= 1 and in orthogonalthe me- whenlocation,q (r) = but 0. A also simple in their flow direction. This becomes EA 0, coincides with that identified from the topo- ↵ ↵ logicalQ ⇠ defect statisticsander in Fig. phase,1e. Taking P( αβ advantage) is not of symmetric.inspection It of the is a maps broad readily indicatesevident uponthat all inspecting exper- the statistics of local overlap. Gaussian distributionQ centered around a positive mean Unlike the global overlap distribution, P(qαβ(r)) value. In other words, even though disorder is isotropic is bimodal and sharply peaks at qαβ(r) = 1, see and homogeneous, the obstacles bias the orientation of Fig.3E. Simply put, the meanders in each pair± of the flow field over macroscopic regions of space, as a re- replicas differ by the reversal of the flow orientation sult reversing the orientation of the flow does not yield over a finite fraction of space. This area fraction is an equally probable pattern. The bias in the flow struc- given by 1 [P(qαβ = +1) P(qαβ = 1)] 0.5, see ture is clearly visible in the replica-averaged velocity field also Supplementary− Note III.− We can then− estimate∼ the shown in Fig3B. typical extent of the compact regions where the flow Crucially, the large width of the overlap distribution can flip sign from one replica to another by measuring demonstrates a broad range of meandering patterns that the correlation length of the local overlap, or of the 5
p
q (r)
0.60.6 100 100 0 100 0 100100 0 00 100100 0.40.4
150 150 -0.5 150 -0.5 150150 -0.5 -0.5-0.5150150 0.20.2
200 200 200 200200 200200 0 0 50 100 150 200 50 100 150 200 50 100 150 200 0.15 5050 100100 150150 200200 5050 100100 150150 200200 c 0 d e 0.3 f 1 10 Overlap 0.2 0.20.2 0.1 0.2 Polarization 0.20.2 0.8 0.20.20.20.20.20.20.10.1 0.1 0.1
0.10.1 )) 0.08 0.1 0 0.2 0.10.10.10.10.10.10 r
100100150150200200 ( 0.6 -2 000 0.06 100100100150150200200 10 000000 ↵ 0100100100100100100150150150150150150200200200200200200 q 100 150 2000.05 0.4
0.02 0.2
10-4 0 0 0 0 -1 0 1 -100 0 100 -1-1 00 1 0 ⇠ 500 1000 1500 q (r)
FIG. 3. Meanders as dynamical vortex glasses. A, Maps of the overlap field qαβ (r) for 4 out of 1225 different pairs of replicas (α, β). Scalebar: 1 mm. B, Magnitude of the replica-averaged polarization field. Disorder robustly sets the flow orientation in macroscopic regions (light colors). Scalebar: 1 mm. C, Overlap distributions in the polar liquid and meander phases. In the polar liquid phase (light orange) the distribution is symmetric and bimodal. Clockwise and counterclockwise vortices are equiprobable. In the meander phase, the distribution is a wide Gaussian centered on 0.46. D, Distribution of the difference between the number of topological defects in each pair of replicas. The distribution is a Gaussian with a standard deviation σ = 18 3 much larger than the uncertainty defined by the standard deviation of the topological charge (which should be constrained± to equal 1 in a circular geometry). Inset: distribution of the number of defects in each replica. The replicated flows are topologically distinct. E, Probability density of the local overlap qαβ (r). The distribution is asymmetric but sharply peaked at 1. Locally the active fluid either flows along identical or opposite directions in each replica of the same ± experiment. F, Decay of the spatial correlation of qαβ (r) averaged over all replicas and of the replica-averaged polarization field plotted in (B). Both decorrelations occur over the same finite distance ξ = 100 10 µm. All panels correspond to experiments Q ± performed deep in the meander phase Φo = 15 %.
replica-averaged flow shown in Fig.3F. We find that the performed in periodic lattices of obstacles, see Supple- overlap decorrelates exponentially over a finite distance mentary Fig. S7 and Movie S7. A systematic investi- ξ = 100 10 µm. This finite correlation length implies gation reveals that no meander forms in the absence of thatQ the meanders± explore a conformational landscape disorder, see Supplementary Note 5. The spatial hetero- including a number of steady states that scales as geneities in the local obstacle fraction φo(r) are essential 2 2(R/ξQ) . The exponential degeneracy of the flows to generate the meandering flow patterns. At lowest or- demonstrates∼ that the meander phase is the dynamical der in gradients, disorder thus minimally alters Toner–Tu analogue of a topological defect glass. It is strongly hydrodynamics for the conserved roller density ρ and the reminiscent of the vortex glass phases found in disor- flow velocity v [1, 34], dered flux lines in superconductors and in random-gauge ∂ ρ = ∇ (ρv), (1) Heisenberg magnets [28, 29, 41]. t − · 2 2 ∂tv + λv ∇v = (a2 a4 v )v + K v β∇ρ βo∇φo , · − | | ∇ − − (2) All the transport coefficients are taken to be constant DYNAMICAL VORTEX GLASS AS A GAUGE positive parameters. When ∇φo = 0, Eq.2 reduces GLASS to the standard Toner–Tu equations, which quantita- tively predicts the vortical flows of Fig.1A[32]. The We now make a precise connection between the dy- terms a , a set the mean value of the flow velocity p2 4 namical vortex glass we observe here in active fluids and (v0 = a2/a4) deep in the homogeneous ordered phase, its static counterpart in superconductors. while λ controls self-advection and the elastic constant K To do so, we construct a hydrodynamic theory of disor- penalizes flow distortions. The βo term reflects the cou- dered polar flows, building on a final set of experiments pling to disorder, and provides a direct qualitative expla- 6 nation for the bias of the overlap distribution in Fig.3C. disorder potential φo into an effective random gauge Although the obstacles are isotropic and, on average, ho- vector field . The local circulation of then de- A A mogeneously distributed, φo(r) acts as a random pres- fines a quenched topological charge in the background sure field that locally drives the flows towards the same Q(r) = zˆ [∇ (r)]. obstacle-depleted regions in all replicas. For weak disor- For a given· × disorder A configuration, the flow stream- der (Φo < Φc), the ordered flows are only marginally per- lines on large scales must satisfy the nonlinear constraint turbed by the obstacles. As detailed in Supplementary in Eq.4, which can be used to obtain the steady state Note 6, we can then perform a simple linearized analysis distribution of θ to be P[θ(r)] limT 0 exp( E/T ) with about the polarized state, deep in the order phase where an auxiliary temperature T and∝ an effective→ − energy (see v = v0(cos θ, sin θ). We readily obtain the steady state Supplementary Note 7 for details) orientational correlator in Fourier space K Z E[θ(r)] = d2r ∇θ (r) 2 ,. (5) β2q2q2 2 o x y 2 | − A | θq = 2 2 2 2 2 2 2 4 Φo , (3) [(c v0qx c qy) + v0K qxq ] k − ⊥ The zero temperature ground state controls the flow con- figurations of the disordered flock (T 0). This com- where the overbar denotes a disorder average. c = λv0 → k pletes our mapping of the nonlinear meandering flows in and c = √βρ are the longitudinal and transverse sound 0 a polar fluid on a disordered substrate to the well known speeds⊥ respectively, and ρ is the mean colloidal roller 0 problem of finding the ground states of a random-gauge density. Velocity fluctuations hence remain finite on large XY model [27, 29, 42]. scales signaling true long-range order. The mean po- The effective energy in Eq.5 captures the competition larization merely decays linearly with Φ , in agreement o between disorder generating meanders dictated by and with our experimental measurements in the small disor- elastic restoring forces penalizing these distortions.A This der limit (Φ < Φ ), see Fig.2C and Supplementary Note o c allows a Kosterlitz-Thouless style argument for prolifer- 6. ating free vortices detailed in Supplementary Note 8. We For strong disorder, we have to go beyond this spin- stress, however, on the crucial role of activity in our anal- wave theory and account for both nonlinearities and the ysis. Although we obtain an equilibrium-like mapping in presence of topological defects. Unlike in equilibrium fer- Eq.5, activity is crucial in generating the gauge field romagnets and active-nematic films, 2D colloidal flocks via self-advection. Noting a the vortex-core size, the support genuine long range order, thereby ruling out the elasticA penalty of an isolated defect K ln(R/a) can be classical Kosterlitz-Thouless picture of defect unbinding. offset by the energy gain KΦ (β∼/λv2) ln(R/a) from Topological defects in clean polar flocks only participate o o 0 optimally screening out the∼ background charge. Their in a rapid coarsening process by annihilating quickly to balance hence predicts a critical threshold Φ λv2/β yield pristine ordered flow, as evidenced in Fig.1A. Dis- c 0 o beyond which pinned vortices proliferate, in agreement∝ order changes this picture in important ways. Quenched with our measurements of the number of frozen defects fluctuations due to the random obstacles provides an ad- at the onset of of the vortex glass formation, Fig.1F. ditional contribution that competes with orientational A more sophisticated analysis including renormaliza- elasticity to uncover a new active defect-unbinding sce- tion corrections to the elastic constant is presented in nario. For simplicity, we fix the colloidal roller density to Supplementary Note 8. This rationalizes the existence be a constant (ρ = ρ ); the sole hydrodynamic mode, is 0 of the dynamical vortex glass phase along with its phase therefore the orientation θ. As detailed in Supplementary transition as being driven by disorder driven unbinding Notes 7 and 8, retaining the full nonlinear dependence on of 1 vortices. θ in Toner-Tu equations, we find that at long-wavelengths ± the convective nonlinearities solely compete with the ran- dom pressure term in Eq.2. To gain some insight on CONCLUDING REMARKS this competition, let us first consider a local minimum of φo(r) at the origin, the inward disorder pressure must be balanced by a centrifugal kinematic pressure λv2/r In conclusion, combining high-content experiments associated with an azimuthal flow. In a more∼ general and theory we have shown how dynamical vortex glasses setting, the competition between convection and random generically emerge when the polar order of flocking mat- pressure constrain the angular fluctuations to obey, ter competes with isotropic disorder. We emphasize the generality of our predictions. Both the topological-defect β stabilisation scenario and the extensive complexity of ∇ o ∇ θ = 2 zˆ P φo (r) , (4) the flow patterns solely rely on stable uniaxial order in a −λv0 × · ≡ A fluid assembled from motile units. Beyond the specifics, where P = 1 vˆvˆ. Remarkably, this quenched distri- of collodial active matter, we therefore expect dynamical bution of phase− distortions is reminiscent of random- vortex glasses to emerge in a host of active materials gauge XY models, and flux-line glasses in type II su- ranging from confined microtubule nematics [11], to con- perconductors, see e.g. [27–29]. Eq. (4) reveals that ac- centrated bacteria suspensions [37], and cell tissues [43] tive self-advection nonlinearly transforms the quenched cruising through disorder. Establishing a quantitative 7 theory of their amorphous patterns remains a formidable inject sequentially a colloidal suspension or a colloid-free challenge. solvent. If not specified otherwise, the stationary flows are mea- Acknowledgements. We acknowledge support sured in the first chamber 120 s after the application of from ANR program WTF and Idex program ToRe. the DC field. This waiting time is more than one or- M. C. M. was supported by the US National Science der of magnitude larger that the flows’ relaxation time, Foundation through award DMR-1938187. S. S. is sup- see Supplementary Note I. For the replica experiments ported by the Harvard Society of Fellows. We thank we proceed as follows: the chambers are filled with the D. Carpentier, O. Dauchot, W. T. M. Irvine, and colloids at the desired area fraction. We motorize the col- A. Morin for insightful comments and suggestions. loids, wait for 120 s and film their motion for 5 s. We then Data availability The data that support the plots switch the voltage off, and switch it on again repeating within this paper and other findings of this study are the same procedure fifty times in a row. available from the corresponding author upon request. Author Contributions D. B. conceived the project. A. C. and D. B. designed the experiments. A. C. B. From Lagrangian trajectories to Eulerian flow performed the experiments and analyzed data. A. C. fields and D. B. interpreted and discussed the experiments. M. C. M. and S. S. performed the theory. All authors In order to track of the trajectories of the rollers, we discussed the results and wrote the manuscript. A. C. image them with a Nikon AZ100 microscope with a 4.8X and S. S. contributed equally. magnification and record videos with a CMOS camera Author Information Correspondence and requests (Basler Ace) at 190 fps. All measurements are systemat- for materials should be addressed to D. B. (email: ically repeated three times for different initial conditions [email protected]). and same disorder configuration. If not specified other- Competing interests The authors declare no compet- wise, we measure all quantities reported in the main text ing interests. after the ensemble of rollers has reached its stationary state.
METHODS 1. Lagrangian trajectories A. Quincke rollers experiments We detect the position of all the rollers with a sub- The experimental setup is similar to that described pixel accuracy using the algorithm introduced by Lu et in [44]. We disperse polystyrene colloids of radius a = al in [46]. We then reconstruct their trajectories over the 2.4 µm (Thermo Scientific G0500) in a solution of hex- whole 3 mm wide circular chambers using the Crocker 2 and Grier algorithm [47] with the MATLAB routine adecane including 5.5 10− wt % of dioctyl sulfosucci- nate sodium salt (AOT).× We then inject the solution in available at [48]. We define the individual roller veloci- microfluidic chambers made of two electrodes spaced by a ties from their displacements over two subsequent frames (time interval: δt = 5.3 ms): v (t) = r (t + δt) r (t), 25 µm-thick scotch tape. The electrodes are glass slides, i i − i coated with indium tin oxide (Solems, ITOSOL30, thick- where ri(t) and vi(t) are respectively the position and ness: 80 nm). We let the colloids sediment on the bottom velocity of particule i at time t. The accuracy of the po- electrode and apply a DC voltage of 120 V. The result- sition measurements is of the order of 0.1 µm, inducing ing electric field triggers the so-called Quincke electro- an accuracy of the order of 40 µm/s for individual speed rotation and causes the colloids to roll at a constant measurements. When powered with an electric field E of magnitude 120 V, all colloids roll at a constant speed: speed v0 = 0.8 mm/s [6, 45]. The microfluidic device is sketched in .1b. We confine the colloidal rollers inside v0 = 0.80 0.04 mm/s. (6) circular chambers of radius R = 1.5 mm including ran- ± dom lattices of circular posts of radius 10 µm. Both the In addition, when isolated, their direction of motion obstacles and the confining disks are made of a 2 µm-thick freely diffuses on the unit circle with a rotational dif- layer of insulating photoresist resin (Microposit S1818) fusivity DR defined as the exponential decorrelation rate patterned by means of conventional UV-Lithography as of the velocity orientation in an isotropic phase: explained in [22]. The patterns are lithographed on the 1 bottom electrode. Note that the distribution of the ob- DR = 2.2 0.1 s− . (7) ± stacle centers corresponds to a planar Poisson process, the circular posts can therefore overlap as see in Fig.1b and Movie S3. The experiments reported in the main 2. Eulerian fields text correspond to thirty different microfluidic chambers including obstacle fractions ranging from 0 % to 38 %. We Building on these Lagrangian measurements, we recon- fine tune the mean roller fraction ρ0 using a T-junction to struct the instantaneous Eulerian velocity fields v(r, t) as 8 follows. We average the instantaneous roller velocities in To compute the polarization order parameter from the 76.4 µm 76.4 µm binning windows arranged on a square instantaneous velocity field, we first average the radial lattice with× a lattice spacing of 15.3 µm. Given the roller and azimuthal components of the polarization field density, each PIV window typically averages the veloc- p(r) ˆv(r, t) over square boxes of size `. We then ≡ h it ity of 25 rollers. We systematically checked that none of compute the spatial average: pB ( pr(r) r, pθ(r) r), our results crucially depends on the specific choice of the in each box B. Finally, (`) corresponds≡ h toi theh normi of P size of the binning windows. The polarization and over- pB averaged over all boxes B. lap fields are computed with the same spatial resolution from v(r, t).
[1] John Toner, Yuhai Tu, and Sriram Ramaswamy, “Hydro- [15] Suraj Shankar, Sriram Ramaswamy, M Cristina dynamics and phases of flocks,” Annals of Physics 318, Marchetti, and Mark J Bowick, “Defect unbinding in 170 – 244 (2005), special Issue. active nematics,” Physical review letters 121, 108002 [2] M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. (2018). Liverpool, J. Prost, Madan Rao, and R. Aditi Simha, [16] John A Blackman and J Tag¨ue˜na, Disorder in condensed “Hydrodynamics of soft active matter,” Rev. Mod. Phys. matter physics: a volume in honour of Roger Elliott (Ox- 85, 1143–1189 (2013). ford University Press, USA, 1991). [3] Andrea Cavagna and Irene Giardina, “Bird flocks as con- [17] Pierre-Gilles De Gennes, Superconductivity of metals and densed matter,” Annu. Rev. Condens. Matter Phys. 5, alloys (CRC Press, 2018). 183–207 (2014). [18] George W Crabtree and David R Nelson, “Vortex physics [4] Hugues Chat´e,“Dry aligning dilute active matter,” An- in high-temperature superconductors,” Physics Today nual Review of Condensed Matter Physics 11, 189–212 50, 38–45 (1997). (2020). [19] Oleksandr Chepizhko, Eduardo G. Altmann, and Fer- [5] Volker Schaller, Christoph Weber, Christine Semmrich, nando Peruani, “Optimal noise maximizes collective mo- Erwin Frey, and Andreas R Bausch, “Polar patterns of tion in heterogeneous media,” Phys. Rev. Lett. 110, driven filaments,” Nature 467, 73 (2010). 238101 (2013). [6] Antoine Bricard, Jean-Baptiste Caussin, Nicolas [20] David A Quint and Ajay Gopinathan, “Topologically in- Desreumaux, Olivier Dauchot, and Denis Bartolo, duced swarming phase transition on a 2d percolated lat- “Emergence of macroscopic directed motion in popula- tice,” Physical biology 12, 046008 (2015). tions of motile colloids,” Nature 503 (2013). [21] Clemens Bechinger, Roberto Di Leonardo, Hartmut [7] Jing Yan, Ming Han, Jie Zhang, Cong Xu, Erik Lui- L¨owen, Charles Reichhardt, Giorgio Volpe, and Gio- jten, and Steve Granick, “Reconfiguring active particles vanni Volpe, “Active particles in complex and crowded by electrostatic imbalance,” Nature materials 15, 1095 environments,” Reviews of Modern Physics 88, 045006 (2016). (2016). [8] Jie Zhang, Erik Luijten, Bartosz A. Grzybowski, and [22] Alexandre Morin, Nicolas Desreumaux, Jean-Baptiste Steve Granick, “Active colloids with collective mobility Caussin, and Denis Bartolo, “Distortion and destruc- status and research opportunities,” Chem. Soc. Rev. 46, tion of colloidal flocks in disordered environments,” Na- 5551–5569 (2017). ture Physics 13, 63 (2017). [9] Julien Deseigne, S´ebastien L´eonard, Olivier Dauchot, [23] Rakesh Das, Manoranjan Kumar, and Shradha Mishra, and Hugues Chat´e,“Vibrated polar disks: spontaneous “Polar flock in the presence of random quenched rota- motion, binary collisions, and collective dynamics,” Soft tors,” Phys. Rev. E 98, 060602 (2018). Matter 8, 5629–5639 (2012). [24] CJ Olson Reichhardt and Charles Reichhardt, [10] Nitin Kumar, Harsh Soni, Sriram Ramaswamy, and “Avalanche dynamics for active matter in heteroge- AK Sood, “Flocking at a distance in active granular mat- neous media,” New Journal of Physics 20, 025002 ter,” Nature Communications 5 (2014). (2018). [11] Amin Doostmohammadi, Jordi Ign´es-Mullol, Julia M [25] John Toner, Nicholas Guttenberg, and Yuhai Tu, Yeomans, and Francesc Sagu´es,“Active nematics,” Na- “Swarming in the dirt: Ordered flocks with quenched ture communications 9, 1–13 (2018). disorder,” Phys. Rev. Lett. 121, 248002 (2018). [12] Suraj Shankar, Sriram Ramaswamy, and M Cristina [26] John Toner, Nicholas Guttenberg, and Yuhai Tu, Marchetti, “Low-noise phase of a two-dimensional active “Hydrodynamic theory of flocking in the presence of nematic system,” Physical Review E 97, 012707 (2018). quenched disorder,” Phys. Rev. E 98, 062604 (2018). [13] Tim Sanchez, Daniel T. N. Chen, Stephen J. DeCamp, [27] Michael Rubinstein, Boris Shraiman, and David R Michael Heymann, and Zvonimir Dogic, “Spontaneous Nelson, “Two-dimensional xy magnets with random motion in hierarchically assembled active matter,” Na- dzyaloshinskii-moriya interactions,” Physical Review B ture 491, 431–434 (2012). 27, 1800 (1983). [14] Luca Giomi, Mark J Bowick, Xu Ma, and M Cristina [28] Thomas Nattermann and Stefan Scheidl, “Vortex-glass Marchetti, “Defect annihilation and proliferation in ac- phases in type-ii superconductors,” Advances in Physics tive nematics,” Physical review letters 110, 228101 49, 607–704 (2000). (2013). [29] David Carpentier and Pierre Le Doussal, “Topological 9
transitions and freezing in xy models and coulomb gases “Transition from turbulent to coherent flows in confined with quenched disorder: renormalization via traveling three-dimensional active fluids,” Science 355, eaal1979 waves,” Nuclear Physics B 588, 565–629 (2000). (2017). [30] G. Quincke, “Uber rotationen im constanten electrischen [39] Achini Opathalage, Michael M Norton, Michael PN Ju- felde,” Annalen der Physik , 417–486 (1896). niper, Blake Langeslay, S Ali Aghvami, Seth Fraden, and [31] Oleg D Lavrentovich, “Active colloids in liquid crystals,” Zvonimir Dogic, “Self-organized dynamics and the transi- Current opinion in colloid & interface science 21, 97–109 tion to turbulence of confined active nematics,” Proceed- (2016). ings of the National Academy of Sciences 116, 4788–4797 [32] Antoine Bricard, Jean-Baptiste Caussin, Debasish Das, (2019). Charles Savoie, Vijayakumar Chikkadi, Kyohei Shitara, [40] Andrea Cavagna, “Supercooled liquids for pedestrians,” Oleksandr Chepizhko, Fernando Peruani, David Saintil- Physics Reports 476, 51–124 (2009). lan, and Denis Bartolo, “Emergent vortices in popula- [41] Matthew P. A. Fisher, “Vortex-glass superconductivity: tions of colloidal rollers,” Nature communications 6, 7470 A possible new phase in bulk high-tc oxides,” Phys. Rev. (2015). Lett. 62, 1415–1418 (1989). [33] Delphine Geyer, David Martin, Julien Tailleur, and De- [42] David Carpentier and Pierre Le Doussal, “Disordered xy nis Bartolo, “Freezing a flock: Motility-induced phase models and coulomb gases: renormalization via traveling separation in polar active liquids,” Physical Review X 9, waves,” Physical review letters 81, 2558 (1998). 031043 (2019). [43] G Duclos, C Blanch-Mercader, V Yashunsky, G Salbreux, [34] John Toner and Yuhai Tu, “Long-range order in a two- J-F Joanny, J Prost, and Pascal Silberzan, “Spontaneous dimensional dynamical XY model: How birds fly to- shear flow in confined cellular nematics,” Nature physics gether,” Phys. Rev. Lett. 75, 4326–4329 (1995). 14, 728–732 (2018). [35] Maria Zeitz, Katrin Wolff, and Holger Stark, “Active [44] Delphine Geyer, Alexandre Morin, and Denis Bartolo, brownian particles moving in a random lorentz gas,” The “Sounds and hydrodynamics of polar active fluids,” Na- European Physical Journal E 40, 23 (2017). ture materials 17, 789–793 (2018). [36] Alexandre Morin, David Lopes Cardozo, Vijayakumar [45] J. R. Melcher and G. I. Taylor, “Electrohydrodynamics: Chikkadi, and Denis Bartolo, “Diffusion, subdiffusion, A review of the role of interfacial shear stresses,” Annual and localization of active colloids in random post lat- Review of Fluid Mechanics 1, 111–146 (1969). tices,” Phys. Rev. E 96, 042611 (2017). [46] Peter J Lu, Peter A Sims, Hidekazu Oki, James B [37] Hugo Wioland, Enkeleida Lushi, and Raymond E Gold- Macarthur, and David A Weitz, “Target-locking acqui- stein, “Directed collective motion of bacteria under chan- sition with real-time confocal (tarc) microscopy,” Optics nel confinement,” New Journal of Physics 18, 075002 express 15, 8702–8712 (2007). (2016). [47] John C Crocker and David G Grier, “Methods of digital [38] Kun-Ta Wu, Jean Bernard Hishamunda, Daniel TN video microscopy for colloidal studies,” Journal of colloid Chen, Stephen J DeCamp, Ya-Wen Chang, Alberto and interface science 179, 298–310 (1996). Fern´andez-Nieves, Seth Fraden, and Zvonimir Dogic, [48] D. Blair and E. Dufresne, “The matlab particle tracking code repository,”.