Arxiv:1808.01828V1 [Gr-Qc] 6 Aug 2018 Rvt Paradigm Gravity Hydrody- the of Is- of Exploration This Namics ﬁrst-Ever Addresses a Letter Through Our Sue Elusive

Home , Hawking radiation

arXiv:1808.01828v1 [gr-qc] 6 Aug 2018 rvt paradigm gravity hydrody- the of is- of exploration this namics first-ever addresses a Letter through Our sue elusive. hole black remained fluid has the in frame-dragging gyro- detect the freely-precessing to scope crystals, a of paramagnetic construction certain conceptual in this for phonons of spin the quantum observation induced for an [10] involving phenomenon proposal a despite How- the ever, black for frequency. precession given acoustic Lense-Thirring predictions rotating acoustic quantitative and a [10] for an- hole recently been has very phenomenon alyzed The this of hith- situations. analogue been gravity acoustic have strong [9], in unobservable holes erto black rel- rotating general for in studied ativity well gyroscopes, of precession fluid as ! [3] experiments [2], laboratory analogues are mechanical to spacetime accessible physical in now been observed gravitational have not kinematic to Thus, phenomena mildly reported [8]. water for already in recently is observed more [7], and fluids [6] viscous for [5], initially fluids super- inviscid predicted acoustic Further, (‘superresonance’), Un- radiance prediction. confirms [4] original radiation ruh’s observa- Hawking activity experimental acoustic ongoing reported of recently tion robust The a 3]. been [2, has cold and systems, matter grav- atom condensed kinematic in of experimen- phenomena, analogues itational and acoustic of theoretical studies holes, tal black such in from analogues of prediction hole the servable black and fluids, acoustic inviscid barotropic, of [1] discovery Introduction: nrilfaedagn n h Lense-Thirring the and frame-dragging Inertial atraAon nAosi lc oe rpigadFrame- and Trapping Hole: Black Acoustic an Around Bacteria idyvsossadr ud,i hspeitdt cu fo occur to predicted thus o is fluids, sign standard the analyz viscous v on theoretically mildly net superradiance, Depending the enhance acoustic viscosity. or situation, effective behaviour, active superfluid-like small yield of a either diffusivity po to dynami of the lead the inclusion up to anticipated, addition, opens In as indeed that, particles gyroscopes. show individual of further We larization’) parti of hole. possibility the black towards reveal hinting sp concentration horizon, smal hole particle acoustic vanishingly active black For the analogue of potential. acoustic perturbations velocity discern the we of fluid, turbations nematic a of aino h yrdnmc of hydrodynamics the analogue on acoustic rotating gation in recently predicted dragging akn aito fpoo excitations phonon of radiation Hawking oiae ytene ocnev reypeesn gyrosc freely-precessing conceive to need the by Motivated eai active nematic olwn po nu’ brilliant Unruh’s on up Following egnrlz nu’ approach Unruh’s generalize We . aarsn iso ieaad nvriy eu ah71 Math Belur University, Vivekananda Mission Ramakrishna eateto hsc,Peiec nvriy okt 700 Kolkata University, Presidency Physics, of Department fluids ihnteanalogue the within mtv Banerjee Amitava eai active nematic ataaah Majumdar Parthasarathi ob- ud.Wt pcfi supino barotropicity on assumption specific a With fluids. ∗ n an Koley Ratna and ocpulzn yocpsoto cieprils in fluids. particles, of such active possibility of the acoustic out study gyroscopes discern to conceptualizing and to analogues, hole perturbations, black linear on based eoiyptnili uhbakhl backgrounds, hole black perturbed such the in describing potential equation velocity the effective weak in matter a viscosity produce active directly fluids, to Remark- seen solvent is diffusivity inviscid system. for nematic even the ably, of barotropicity regard- made ing is ne- assumption in an indeed provided flows fluids, background are matic possible as analogues emerge hole to shown black system- Acoustic out carried atically. are backgrounds, parti- appropri- stationary around active ate orientation, the and density, concentration solvent cle continuity the respective for the equations and equations fluid matic behavior ‘superfluid-like’ possible shear- [14–21]. to in swim- result by leading may generated ming, forces flow the swimmers a by Such active reduction viscosity of flow. alignment shear the a with is with orientations flow, hydrody- swimmers’ fluid key the the to A of due arises interaction which the phenomena, of [11–13]. such in range scales effect broad namic time a collec- and over fascinating length phenomena of dynamical variety tive spectacular a of manifestation appear there the interactions, as particle-particle and systems, fluid-particle such In asymmetric their structure. to an due carry (‘polarization’) and such orientation diffusive Often inherently are itself. interact particles media thereby active the and pro- and media themselves birds, fluid among and in themselves fishes pel to microswimmers artificial ciemtl agents, motile Active olwn nu,lna etrain ftene- the of perturbations linear Unruh, Following dadeprmnal bevdrcnl for recently observed experimentally and ed ltapn fatv atrb h acoustic the by matter active of trapping al rfiewt nehneetcoet the to close enhancement an with profile a lc oe,w eota niin investi- incipient an report we holes, black iuiiyo h cieprils linear particles, active the of diffusivity l atce nteivsi ovn sshown is solvent inviscid the in particles soiy ftenmtcsse.I either In system. nematic the of iscosity, sblt fteruea freely-precessing as use their of ssibility eai fluids. nematic r ‡ ctmseprecdb ierper- linear by experienced acetimes psfrdtcigaosi frame- acoustic detecting for opes a aueo h retto (‘po- orientation the of nature cal h iuincecet hscan this coefficient, diffusion the f † 7,India. 073, 22 India 1202, e.g. lntn,bacteria, planktons, Dragging given in terms of the diffusion coefficient. Adjusting viscid fluids, with the total stress tensor σij (whose the latter thus paves the way for reducing shear vis- full expression is given in [18]) now containing con- cosity and confirming the possibility of simulating tributions generated by the swimming of the active an active ‘superfluid’, as also the alternative likeli- elements as well as the pressure term −Πδij hood of enhancing the net viscosity. Whereas analogue gravity studies traditionally rely on inviscid (ρ + c) (∂t + v.∇) v = ∂j σij fluids [2, 3], in view of recent work [7], [4] on rotat- Π = Π(ρ,c) = Π(ρ + c) , (1) ing acoustic black holes in mildly viscous (‘inactive’) fluids, both situations offer the opportunity for the where, the last equation is taken to characterize incipient observation of acoustic superradiance [5] in barotropicity in the nematic fluid in this Letter. In nematic fluids. this respect, since the corresponding assumption of For vanishingly small diffusivity, the perturbed ac- an equation of state relating pressure and density tive matter concentration around a vortex flow cor- of an active fluid is highly nontrivial and debated responding to a rotating acoustic black hole, exhibits [22, 23], this crucial assumption restricts the classes an intriguing profile : a sharp enhancement of con- of active systems that the analysis in this Letter can be relevant to. Under such an assumption, there is centration is manifest, close to the acoustic horizon 2 ! We shall present a heuristic derivation of this phe- a unique speed of sound given by vs ≡ dΠ/d(ρ + c). nomenon in the sequel. Such a behaviour is strongly The continuity equation for the solvent fluid is un- indicative of the possibility of externally trapping changed and is given by bacteria and other nematic active elements in a fluid ∂ ρ + ∇ · (ρv) = 0 (2) by the acoustic black hole, and may have applica- t tions in spatial control of bacterial concentrations in For the active part, we assume the swimmers in a po- biologically active fluids in general. larized state with high orientational order among the Returning to our initial motivation for investigat- polarization vectors of individual particles. In that ing active fluids, we consider with some approxi- situation, as shown in Ref. [18], we can assume |P| mations, the continuity equation corresponding to to be a constant, and work with p ≡ P/|P| instead. linear perturbations in the orientation of individual Then, the hydrodynamic equations for the variables active particles. Numerical solution of the pertur- corresponding to the active elements of such a sys- bation equation yields a dynamical orientation angle tem can be written as [18] with respect to the direction of the flow, depending ′ on the flow parameters of the rotating acoustic black ∂tc + ∇ · cv = ∂i [Dij ∂j c + λγ uklpkplpi] (3) hole. Such a strongly dynamical (time-dependent) orientation of an individual bacterium, depending and on the background flow parameters, makes this ac- [∂ + v · ∇]p + ω p = δT λu p + κ∇2p (4) tive particle a rather suitable candidate for a freely- t i ij j ij jk k j precessing gyroscope. Even though technical details ′ T are still awaited, we envisage immediate application with γ ,κ,λ constants, δij = δij − pipj being the of this phenomenon (of dynamical orientation of in- transverse projection operator , uij ≡ (∂ivj + dividual active particle) in acoustic black hole space- ∂j vi)/2 , ωij ≡ (∂ivj − ∂j vi)/2. eqn. (3) is just times), to the detection of acoustic frame-dragging the continuity equation for the concentration of ac- and measurement of the acoustic Lense-Thirring fre- tive particles, appended with an extra term to ac- quency. Since no ‘weak-field’ approximation is in- count for their diffusion in the background fluid with volved in our analysis, nematic fluids stand to en- Dij being the effective diffusion tensor governing able reliable laboratory observations of acoustic ana- orientation-dependent diffusion of active nematics logues of strong-gravity frame-dragging phenomena and is given by for the first time ! Dij = D(δij − ξpipj ) (5) Model of the Active Fluid And Realization of Ana- logue Gravity: We begin with the hydrodynamic where D is the diffusion constant and ξ is another equations appropriate to a nematic active fluid [18]. constant related to how much diffusion is influenced The total system is described by the following vari- by the alignment of active particles. The parameter ables: fluid velocity v, density of the fluid solvent |λ| controls the tendency of the swimming particles ρ, density of the active particles c and their local to align to a shear flow in the background fluid. The orientation P. At this point, we do not assume any polarization equation (4) contains the usual convec- relative velocity of the active particles with respect tive derivative for a flowing fluid and terms involving to the background fluid. The background fluid fol- the symmetric and antisymmetric derivatives of the T lows the conventional Navier-Stokes equation for in- flow velocity (ωij pj and δij λujkpk) accounting for

2 alignment of the active particles to sheared and ro- tential being ψ, so that vi = −∂iψ and uij = tating fluid flows. Finally, the contribution of elastic −∂i∂j ψ and ωij = 0. We decompose all dynami- energy cost of deforming the polarization field of the cal fields into their time-independent steady back- particles, which arise due to the aligning interactions grounds (denoted by terms with subscript 0) and between individual particles [18], is modeled in the linear fluctuations (denoted by terms with subscript last term in right hand side of Eq. (4). This model 1) : v = v0 + ǫv1, ρ = ρ0 + ǫρ1, c = c0 + ǫc1, Π= of active matter hydrodynamics can be derived phe- Π0 + ǫΠ1, ψ = ψ0 + ǫψ1, with ǫ << 1. In this anal- nomenologically or from microscopic models of the ysis, we keep the orientation field at its static back- active particle dynamics [18] and as explained above, ground, and neglect the effect of activity in Eq. (1) is quite intuitive. at order ǫ. With these assumptions and approxima- We proceed to analogue gravity through linear tions, the equation governing the velocity potential perturbations for this fluid, following ref. [2], as- ψ1 may be written as suming an irrotational flow with the velocity po-

c + ρ c + ρ −∂ 0 0 (∂ ψ + v .∇ψ ) + ∇. (c + ρ )∇ψ − 0 0 (∂ ψ + v .∇ψ )v t v2 t 1 0 1 0 0 1 v2 t 1 0 1 0 s s c = ∂ −D ∂ 0 (∂ ψ + v .∇ψ ) + λγ′p p p ∂ ∂ ψ . (6) i ij j v2 t 1 0 1 k l i k l 1 s

One can interpret Eq. (6) in terms of the stan- that respect, it is similar to the known fact [14–21] dard analogue gravity acoustic spacetime [2]. thus that flow-aligning swimmers can enhance or reduce the LHS of Eq. (6) reduces to the covariant the viscosity of the solvent by aligning to a shear 2 d’Alembertian wave operator g corresponding to flow and generating flow patterns about themselves, the acoustic metric given by reinforcing or opposing the existing shear. On one hand, these effects help impart ‘superfluidic’ prop- erty to active systems [15], which is usually an es- j µν 1 −1 −v0 sentiality [2] to realize conventional analogue gravity. g = i 2 ij i j (ρ0 + c0)vs " −v0 vs δ − v0v0 # On the other hand, these terms contribute as explic- itly Lorent-violating terms [7] in the model, thereby acting on the perturbed potential ψ1. The RHS facilitating simulations of the physics of Lorentz vi- is clearly not Lorentz invariant and thus offers new olation [24]. From studies on such systems, one con- opportunities to simulate Lorentz-violation Physics cludes, as demonstrated in [7], acoustic superradi- [24]. In particular, neglecting polarization terms in ance (‘superresonance’), i.e., amplification of sound Eq. (6) one obtains, waves at the cost of the rotational energy of the acoustic analogue black hole, is certainly possible for nematic fluids. The physical effects of such phe- ρ0 + c0 2 Dc0 2 ψ = ∇ (∂t + v · ∇)ψ (7) nomena on active systems is left as a possible future v2 g 1 v2 0 1 s s work. where we have used the relation Active Matter Concentrations in Draining Sink 2 c1 = (c0/vs )(∂t + v0 · ∇)ψ1 (8) Flow: Since the equation for the perturbed velocity potential in the active fluid is complicated in general, which follows from barotropicity. An interesting as- we shall first make several simplifying assumptions. pect of Eq. (7) is its striking similarity with the wave We neglect the contributions from the diffusion of ac- equation for ψ derived in conventional barotropic 1 tive matter by assuming D <

3 vortex profile the horizon, and a corresponding monotonic decay away from it. Furthermore, we also see that the ˆ −Arˆ + Bφ length scale of the radial decay of this oscillation v0 = (10) r amplitude can be controlled by the parameter B/A with A > 0 and B being two parameters. We also of the background flow. In general, this spatial assume a homogeneous c0 and ρ0 profile leading variation of the active matter density wave ampli- to homogeneous sound speed. This velocity profile tude comes as a direct consequence of the underly- is widely used in the analogue gravity community ing acoustic metric which controls the perturbation [2, 3, 8] and indeed vortical flows are ubiquitous in in the velocity potential ψ1. As seen from Eq.(8), active matter systems [25–30] as well. The corre- this variation can be controlled totally by manipu- sponding metric is well-known to admit an acoustic lating the background fluid flow v0. As such, this be A useful in localizing bacteria or colloidal particles or trapping horizon at r = rH = vs where the inward radial component of the fluid velocity equals sound in targeted delivery of such objects with a carefuly speed and to demonstrate superresonance [5, 6] and designed flow. This general scheme of controlling inertial frame dragging [10]. To calculate the fluctu- amount of fluctuation in active matter concentration ation in the active matter concentration c1, we shall will complement the existing methods of trapping use an explicit solution of the homogeneous wave and manipulation of active particles using, for ex- equation, which is valid near the horizon and given ample, topological defects in the orientation profiles by [5, 6] [31–33], curvature [34], deformation [35] or constric- tion [36] of the bounding walls, swirl flows [25, 37] or imB i(mφ−ωt) 2A ln(1+x) ψ1m(r, φ, t)= e e F (x) (11) shear flows [38, 39], most of which are very system- specific. where m ∈ I+ is the mode number, ω > 0 is the 2 r Active Nematic Particles as Probes of Acoustic frequency of the mode, x = 2 − 1 is the normalized rH Lense-Thirring Precession: Rotating acoustic black squared radial distance from the horizon and F (x)= hole spacetimes like that realized in DS offer the pos- 2F1(α, β; γ; −x) is the hypergeometric function with sibility of observation of acoustic inertial frame drag- 2 2 ωA Bm 2 2Bmω 2ω A Q = 2 − ,S = m + 2 − 4 , α = ging and Lense-Thirring precession [10], wihtout re- vs A vs vs Q Q sorting to any weak field approximation. In this re- − S − i ,β = S − i q,γ = 1 − iQ. This solution 2 2 2 2 spect, active particles in the vortex flow can work as is valid outside the horizon and in the low-frequency ωA ωA freely-rotating gyroscopes, probing acoustic frame- regime where the relations 2 x<mass position x of the individ-

ual particles and write their equations of motion as: 1 B=0 x˙ = v (13) B=1 B=2 which implies that the center of mass of each parti- B=3 cle drifts with the same speed as the fluid flow while 0.5 B=4 their orientations tends to align to local shear in the B=5 flow according to the Jeffrey’s equation[40] corresponding to an irrotational flow

Normalized Amplitude p˙i = βujkpk(δij − pipj ) (14) 0 0 5 10 15 20 x where β is a parameter depending on the shape and aspect ratio of individual active particles. If we write 2 p |c1(x)| = (cos θ, sin θ) then, using Jeﬀrey’s equation this FIG. 1: Active particle concentration proﬁles 2 for |c1(0)| leads to the equation for orientation angle θ as vs = 1,ω = 0.001,m = 1,A = 1. uxx − uyy 2 ∂tθ = β uxy cos2θ − sin 2θ . (15) there is an enhancement of the amplitude |c1| near 2

4 We numerically solve this equation for the DS ve- time needed to reach the horizon from where the locity profile we used above and plot the orientation particle starts. In those cases, we can use the dy- fluctuation due to the velocity perturbation in the namics of active particles to probe the correspond- fluid. The resulting Fig. 2 shows that the tempo- ing dynamics of test objects in the emergent acoustic spacetime. (a) 2.5 Concluding Remarks: In the foregoing, active nematic fluids have been shown to emerge as a novel 2 Background arena for analogue gravity phenomena. We have demonstrated several implications of this symbiotic 1.5 Total relationship between active fluids and acoustic grav-

θ 1 ity analogues, working within certain approxima- tions for simplicity. With regard to future outlook, 0.5 inclusion of the contribution of the diffusion term in 0 the right hand side of our wave equation, neglected in our analyses so far, is a high priority. While ignor- -0.5 0 4 8 ing this contribution is reasonable for the limits of r low concentration or low diffusivity which we men- (b) tioned earlier, inclusions of such terms may account 3 for simulations of more exotic phenomena. Further- more, for small D, the effect of the source term in 1.5 Background the wave equation can be evaluated perturbatively Total starting from a solution of ψ1 for the corresponding homogeneous wave equation. Analysis of such /dt 0 θ

d kinds is inspired by the fact that since the diffusion (by eqn.(5)) is dependent on the polarization -1.5 of the particles, it can be controlled experimentally if we can manipulate the alignment of the particles -3 by, e.g., the fluid flow itself as we discussed, pas- 0 4 8 r sive liquid crystals [41], intense laser beams [42] or external magnetic fields for magnetotactic bacteria FIG. 2: Dynamics of the (a) orientation θ and (b) an- [43]. This would allow us greater control over the gular velocity dθ/dt of active particles with radial coor- dispersion of the active particles over the fluid, com- dinate r; with ω = 0.001, vs = 1,m = 1,A = 1, B = plementing our proposed way of doing the same by 10, β = 1, ǫ = 0.1. The initial conditions are x = 8,y = imposing analogue spacetime structures to the flow 0,θ = 0, the ‘background’ corresponds to the dynamics perturbations. In summary, since both active mat- with v = v0 while the case with v = v0 + ǫRe[v1] is ter and analogue gravity enjoy established stature as plotted as ‘total’. robust areas of reseach in contemporary theoretical and experimental sciences, this incipient investiga- ral fluctuations in the orientation due to ψ1, rep- tion of analogue gravity in active nematic fluids of- resenting precession, are particularly prominent for fers myriad possibilities of rich and novel phenomena larger B/A ratio and are enhanced as the particles involving both biology and physics. approach the acoustic horizon. The former observa- Acknowledgements: We thank O. Ganguly and S. tion can be attributed to the fact that in our acous- Chatterjee for useful discussions. One of us (PM) tic metric, the angular momentum of the black hole 2 2 thanks C. Chakraborty, L. Crispino and J. Kunz for horizon ΩH = (Bvs )/A [10] increases with B/A and illuminating remarks and discussions. frame-dragging effects increase accordingly. This precession of the orientation of individual active particles can be experimentally observed directly. The above simulation results can thus be compared with References those observations possibly after averaging over the dynamics of a sufficiently large number of individual [1] W. G. Unruh, Phys. Rev. Lett, 46, 13511353 particles to eliminate random fluctuations in their (1981). position and orientations resulting in translational [2] C. Barcelo, S. Liberati, and M. Visser, Living and rotational diffusion. The effects would be best Rev. Relativ. 14: 3 (2011). observable for very small rotational diffusion constant (Dr) values relative to the orientation fluctua- [3] “Analogue Gravity Phenomenology Analogue 2 tions δθ due to v1, i.e., 2Drτ <<δθ where τ is the Spacetimes and Horizons, from Theory to

5 Experiment”, Editors: D. Faccio, F. Bel- [22] A. P. Solon, J. Stenhammar, R. Wittkowski, giorno, S. Cacciatori, V. Gorini, S. Liberati, M. Kardar, Y. Kafri, M. E. Cates, and J. U. Moschella, Springer (2013). Tailleur Phys. Rev. Lett. 114, 198301 (2015).

[4] J. Steinhauer, Nature Physics volume 10, 864- [23] A. P. Solon, Y. Fily, A. Baskaran, M. E. 869 (2014). Cates, Y. Kafri, M. Kardar, J. Tailleur, Na- ture Physics 11 (8), 673 (2015). [5] S. Basak and P. Majumdar, Class. Quantum Grav. 20 3907 (2003). [24] S. Liberati and L. Maccione, Ann. Rev. Nucl. Part. Sci. 59:245-267 (2009). [6] S. Basak and P. Majumdar, Class. Quantum Grav. 20 2929 (2003). [25] M. Tarama, A. M. Menzel, and H. Lowen, Phys. Rev. E 90, 032907 (2014). [7] O. Ganguly, arXiv:1705.04935v1 (2017). [26] L. Giomi, Phys. Rev. X 5, 031003 (2015). [8] T. Torres, S. Patrick, A. Coutant, M. Richartz, E. W. Tedford, and S. Weinfurtner, Nature [27] A. Doostmohammadi, M. F. Adamer, S. P. Physics volume 13, 833-836 (2017). Thampi, and J. M. Yeomans, Nature Com- munications volume 7, Article number: 10557 [9] “Einstein Gravity in a Nutshell”, A. Zee, (2016). Princeton University Press (2013). [28] M. Neef and K. Kruse, Phys. Rev. E 90, [10] C. Chakraborty, O. Ganguly, and P. Majum- 052703 (2014). dar, Ann. Phys. (Berlin) 529, 1700231 (2017). [29] H. Wioland, F. G. Woodhouse, J. Dunkel, J. [11] M. C. Marchetti, J. F. Joanny, S. Ramaswamy, O. Kessler, and R. E. Goldstein, Phys. Rev. T. B. Liverpool, J. Prost, Madan Rao, and Lett. 110, 268102 (2013). R. Aditi Simha, Rev. Mod. Phys. 85, 1143 [30] Y. Yang, F. Qiu, and G. Gompper, Phys. Rev. (2013). E 89, 012720 (2014). [12] C. Bechinger, R. Di Leonardo, H. Lowen, C. [31] M. M. Genkin, A. Sokolov, O. D. Lavrentovich, Reichhardt, G. Volpe, and G. Volpe, Rev. and I. S. Aranson, Phys. Rev. X 7, 011029 Mod. Phys. 88, 045006 (2016). (2017). [13] E. Lauga, Annu. Rev. Fluid Mech. 48:10530 [32] C. Peng, T. Turiv, Y. Guo, Q.-H. Wei, O. D. (2016). Lavrentovich, Science Vol. 354, Issue 6314, pp. [14] E. Clement, A. Lindner, C. Douarche, and H. 882-885 (2016). Auradou, Eur. Phys. J. Special Topics 225, [33] P. Guillamat, J. Ignes-Mullol and F. Sagues, 23892406 (2016). Nature Communications volume 8, Article [15] H. M. Lopez, J. Gachelin, C. Douarche, H. Au- number: 564 (2017). radou, and E. Clement, Phys. Rev. Lett. 115, [34] O. Sipos, K. Nagy, R. Di Leonardo, and P. 028301 (2015). Galajda, Phys. Rev. Lett. 114, 258104 (2015). [16] Y. Hatwalne, S. Ramaswamy, M. Rao, and R. [35] P. W. Miller, N. Stoop, and J. Dunkel, Phys. A. Simha, Phys. Rev. Lett. 92, 118101(2004). Rev. Lett. 120, 268001 (2018). [17] M. E. Cates, S. M. Fielding, D. Marenduzzo, [36] M. Potomkin, A. Kaiser, L. Berlyand1 and I. E. Orlandini, and J. Yeomans, Phys. Rev. Aranson, New J. Phys. 19 115005 (2017). Lett. 101, 068102 (2008). [37] A. Sokolov and I. S. Aranson, Nature Com- [18] L. Giomi, T. B. Liverpool,and M. C. munications volume 7, Article number: 11114 Marchetti, Phys. Rev. E 81, 051908 (2010). (2016). [19] A. Sokolov, and I. S. Aranson, Phys. Rev. [38] R. Rusconi, J. S. Guasto and R. Stocker, Na- Lett. 103, 148101 (2009). ture Physics volume 10, 212 - 217 (2014). [20] J. Gachelin, et al. Phys. Rev. Lett. 110, [39] R. N. Bearon, and A. L. Hazel, Journal of 268103 (2013). Fluid Mechanics, vol. 771, id. R3 (2015). [21] J. Slomka and J. Dunkel, Phys. Rev. Fluids [40] J. S. Guasto, R. Rusconi, and R. Stocker, 2, 043102 (2017). Annu. Rev. Fluid Mech. 44:373400 (2012).

6 [41] P. Guillamat, J. Ignes-Mullol, F. Sagues, Molecular Crystals and Liquid Crystals 646 (1), 226-234 (2017). ∗ Electronic address: [email protected] † [42] A. Bezryadina et al, Phys. Rev. Lett. 119, Electronic address: [email protected] ‡ Electronic address: [email protected] 058101 (2017). [43] P. Guillamat, J. Ignes-Mullol, F. Sagues, Pro- ceedings of the National Academy of Sciences 113 (20), 5498-5502 (2016).