Arxiv:1008.2681V1 [Physics.Bio-Ph]

arXiv:1008.2681v1 [physics.bio-ph] 16 Aug 2010 ioe rsrslt[2,wihdcy ihdistance with decays which [12], stresslet or dipole, radius oyfrepootoa otecl ocnrto and concentration cell the to force proportional gravitational as coarse-grained force as (Stokeslets), body equations monopoles the a Navier-Stokes force express of the instability sum to a this cells bio- of of of one Models contribution gravity, plumes buoyancy-driven of [14]. the effects convection na- consider striking in only the realized need appreciate generally To not ture. [9] density-matching cise approaches [13]. promising microorganisms most of the behavior of collective some ten- to to stress fluid central the is to stresslets sor such of suspension a from d/ηr T dis- are body their distance on drag a viscous placed the and flagella their tkse easas decays Stokeslet aissrl ly e oei hs ffcs detailed around a field effects, flow these the in of role hydrody- key knowledge As a plays 10]. surely [9, in- namics surfaces the 7] with and [6, organisms [8], of transport suspensions teraction [5], of yet behavior properties exists rheological collective consensus and of such No origins the [4]. fluid on propulsive viscous a generate mo- in appendages mechanical forces how microorganism on [3], of consensus algae general tions and and a [2] [1] now is spermatozoa bacteria there of of flagella flagella undulating rotating the the of studies tailed hc xr ontfreo h ud h thrust The fluid. the on force net swimmers no collective buoyant exert neutrally and which assume interactions typically fluid-mediated behavior away. of far the sufficiently with Theories dominating [11], mode solutions expressed singularity decaying be slowest of can superposition microorganism a as a around flow far-field far and measurements. such near-field first the the present in we both Here away. needed, is croorganisms h oc-reielzto fsimr eurspre- requires swimmers of idealization force-free The ie yavne niaigtcnqe htalwde- allow that techniques imaging in advances by Aided h iert fteSoe qain mle htthe that implies equations Stokes the of linearity The 2 R where , ietMaueeto h lwFedAon wmigMicroo Swimming Around Field Flow the of Measurement Direct ,adblnet ietefrfil o faforce a of flow far-field the give to balance and ), sue tese o nyacreto otesbedn so subleading the to correction a of only suspensions flow stresslet assumed xesof excess h einaround region the hooyo hi upnin.Wieetnieysuidt microorgani freely-swimming studied microphytes any extensively the around While detail in suspensions. measured their of rheology ASnmes 87.17.Jj,87.16.Qp,47.63.Gd describ numbers: be PACS can and stresslet, “puller” a from qualitatively wmigmcoraim raeflw htiflec hi mu their influence that flows create microorganisms Swimming ntDece,RyodE odti,NclsMce,MroPln an Polin, Marco Michel, Nicolas Goldstein, E. Raymond Drescher, Knut η stefli’ icst.Tecontribution The viscosity. fluid’s the is d F F g pr otncmaal otecell the to comparable (often apart .carteri V. g e el[4.A h o ruda around flow the As [14]. cell per /ηr ti la,i o appreciated not if clear, is it , .carteri V. ovxcarteri Volvox eateto ple ahmtc n hoeia Physics Theoretical and Mathematics Applied of Department .reinhardtii C. vrwtrlast togydmnn tkse contribut Stokeslet dominant strongly a to leads water over nvriyo abig,CmrdeC30A UK 0WA, CB3 Cambridge Cambridge, of University reyswimming freely aefaue iia ossesoso eietn partic sedimenting of suspensions to similar features have hr infiathdoyai neato slkl oocc to likely is interaction hydrodynamic significant where and Dtd coe 4 2018) 24, October (Dated: hayooa reinhardtii Chlamydomonas T r mi- as of rcrdnmc 7.Tento fna edrqie dis- requires field distances near between of to notion tinguishing and The [16] interac- [7]. cell-cell dynamics suspensions tracer to concentrated relevant in is especially field tions, near form Stokeslet the observ- the but situ- concentrations the suffices, low physical and At concentration particular interest. cell of a able the the to on organism, Λ depends given length ation vis- a the For the of than [15]. relevance smaller depth ef- scales penetration crossover unsteady on cous several ciliates, important of for become one swimmers; fects is to This relevant dis- lengths the regime. to over Stokeslet crosses tant contribution stresslet nearby the rvosy htteei itneΛ distance a is there that previously, 1] peia la( alga spherical a [18], unicel- the ( for singu- alga [17] topology force green flow competing biflagellate near-field the lular the quantify to and here larities used is stresslet. try a of that from differ can and topology well-ordered flow be the not may contributions multipole the hr tese ecito a od and hold, may description stresslet a where cinof action scmail ihasml ple”srsltol tdis- at only stresslet “puller” tances simple a doublet, with around source compatible flow is a The of stresslet. form smaller a the and in correction near-field ing ( organisms unicellular common of that than oevr h ihsmer of symmetry high the Moreover, mere component, a of Stokeslet excess its density by a dominated despite strongly is field For et h nt eaaino h aelradbd forces. body re- and topology flagellar the flow of the separation cell, finite the the flects to closer speed; swimming A eim[7 na ria hkr ohi diurnal a in both shaker, orbital an whereas on bubbling, [17] medium air TAP sterile with reinhardtii [19] C. SVM in enically ytei ftakn irsoyadfli velocime- fluid and microscopy tracking of synthesis A .carteri V. Volvox db ipetreSoeltmodel. three-Stokeslet simple a by ed R ertcly hs oshv o been not have flows these heoretically, & m erpr uhmaueet for measurements such report We sm. ∼ ∼ redpl em hsipisthat implies This term. dipole urce 7 5 R 10 u otsgicn nigi htteflow the that is finding significant most our h iue( minute The . f. hr h udvlct is velocity fluid the where , µ ulitrcin n oiythe modify and interactions tual 3 sri TX8)wsgonaeial in axenically grown was 89) UTEX (strain )adislre relative larger its and m) nagariensis Chlamydomonas R o,wt h widely- the with ion, ∼ , 200 sri V)wsgonax- grown was EVE) (strain ∼ dnTuval Idan d hayooa reinhardtii Chlamydomonas e.Teflwin flow The les. 0 r µ . lk el nissurface. its on cells -like % density 3%) Volvox )wihsisb the by swims which m) satisfying ∼ rdiffers ur rganisms 0 ∼ . % uhsmaller much 3%, eut nalead- a in results Chlamydomonas d/F T ovxcarteri Volvox . R r ∼ ∼ g ≪ %o the of 1% 5 twhich at R − r where 10%). ≪ Λ, 2

scope at 40× (NA 0.6) by exciting their chlorophyll aut- ofluorescence with a laser (635 nm, ∼ 60 mW), which also excited 1.6 µm fluorescent polystyrene microspheres (FS04F, Bangs Labs) used as tracers. Cylindrical poly- dimethylsiloxane sample chambers (5 mm radius, 0.4 mm height) were prepared, pacified, and filled following [7]. Experiments were performed at 21 ± 1◦ C, with the laser providing the only light source. We focused on a plane 150 µm inside the chamber to minimize surface effects, and recorded movies at 250 fps (Fastcam SA3, Photron). Movies were analysed with standard algorithms to track cells and tracers. For each cell swimming along the focal plane for more than 1 s (∼ 10 body lengths), we collected the instantaneous velocity of all tracers at r< 14 R, nor- malized by the swimmer’s speed. The resulting 3.3 × 106 velocity vectors were binned into a 2.5 µm square grid (shown in Fig. 4 below), and the mean of the well- resolved Gaussian in each bin was used for the flow field. In both experiments U0 indicates the swimmer veloc- ity, while u(r) and v(r)= u(r)−U0 are the velocity field in the laboratory and comoving frames respectively. A typical experimental flow field around Volvox is shown in Fig. 1(a). We fit these fields to a superpo- sition of a uniform background velocity (U0), a Stokeslet (St), a stresslet (str) and a source doublet (sd):

ASt vfit(r)= − U0 yˆ − (I + ˆrˆr) · yˆ (1) FIG. 1: (color online). Flow field of a freely swimming V. r carteri in the laboratory frame. (a,b) Magnitude and stream- Astr 2 Asd I u u − 1 − 3(y/r) ˆr − − ˆrˆr · yˆ lines of and its fitted approximation fit respectively. (c) r2 r3 3  Relative error of the fit: ǫ = ||u−ufit||/||u||. (a) and (b) have
the same colorbar, different from (c). ~g indicates gravity. where I is the unit tensor, yˆ is the upward vertical unit vector, ˆr = r/r, and r is measured from the center of the organism (xc,yc). The orientation of all multipoles is growth chamber with 16 h in artificial cool daylight (∼ fixed to be along the vertical, and we are left with six pa- ◦ ◦ 4000lux) at 28 C, and 8hin the darkat 26 C. The large rameters: (U0, ASt, Astr, Asd, xc,yc). The fits, obtained difference in organism size between Volvox and Chlamy- by minimizing the integrated squared difference between domonas required two distinct methods to measure the the model and the experimental flows, describe remark- flows they create [20]. A CCD camera (Pike, Allied Vi- ably well the experimental flow, almost down to the sur- sion Technologies) mounted on a continuously-focusable face of the organisms [see Fig. 1(b,c)]. Typical values 2 4 2 microscope (Infinivar, Infinity Optics) and connected to for the parameters are U0 ∼ 10 µm/s, ASt ∼ 10 µm /s, 6 3 a vertical motorized XY stage (Thorlabs) followed indi- Astr ∼ 10 µm /s (indicating a pusher-type stresslet), 9 4 vidual Volvox colonies as they swam upwards [21] in a Asd ∼ 10 µm /s, with the actual magnitude depend- straight line along the central axis of a 5 × 5 × 50 mm ing on the colony radius R. From the Stokeslet com- ◦ sample chamber filled with SVM at 21 ± 1 C. The stage ponent, we can calculate the average colony density as 3 −3 was controlled by a custom LabView routine. The fluid ∆ρ =6ηASt/gR , where η = 10 Pa s and g is the grav- −5 was seeded at volume fraction 10 with 1 µm nile-red itational acceleration. The dependence of both ∆ρ and polystyrene microspheres (Invitrogen) illuminated by a U0 on R (Fig. 2) compares well with previously published vertical ∼ 500 µm thick laser sheet (λ = 532 nm). Volvox data [21] obtained by different means, thereby validat- is phototactic [22] at this wavelength, and at the inten- ing the measurements and analysis procedures. Remov- sities used here it swims smoothly along the laser sheet. ing the Stokeslet contribution from the experimental flow We recorded the flow field of 19 different colonies at 30fps field [Fig. 3(a)] reveals that the near field is dominated by for ∼ 2 − 3 min each. The measured flow field v was ob- the source doublet component, with the stresslet respon- tained by particle image velocimetry (Dantec Dynamics). sible only for a slight forward-backward asymmetry [Fig. Background flows in the chamber were < 10 µm/s. 3(b,c)]. The orientation of the source doublet is opposite We observed a dilute suspension (∼ 3 × 106 cells/cm3) to that around a translating solid sphere, and is compat- of Chlamydomonas in TAP on a Nikon inverted micro- ible with a model that assigns a constant force density to 3

FIG. 2: (color online). The dependence of the excess den- sity ∆ρ and swimming speed U0 (inset) of V. carteri colonies FIG. 3: (color online). Near field around V. carteri. (a) (red stars) on their radius R are compatible with previous Magnitude, vector fields and streamlines of u after subtract- measurements (blue circles) [21], and also [25]. ing the fitted Stokeslet. Colorbar as in Fig. 1(a). (b) ||u|| along a horizontal section through the center of the organism. The average Stokeslet (green dashed line) follows the decay the colony surface [23], as well as a particular case of the of the experimental flow (blue circles) averaged over 19 differ- ent colonies (black dots). Deviations from a pure monopole “squirmer” model [24]. Average values of the parameters appear from . 5R, and can be captured adding a source dou- in Eq. 1 show that the crossover distance between the blet and a stresslet (red solid line). (c) Vertical section of the source doublet and stresslet (∼ 3 R) is beyond that at flow field u from the experiment in (a) through the center which the Stokeslet becomes the leading component of of the colony. The stresslet component is responsible for the the flow (∼ 1.5 R). This peculiar ordering of multipoles forward-backward asymmetry. Symbols as in (b). results from the high degree of anterior-posterior sym- metry of Volvox [23], and highlights the influence that a swimmer’s body plan can have on its flow field. distance [see Fig. 4(c)], are very similar to those mea- The two ∼ 12 µm long flagella of Chlamydomonas sured. Including no-slip boundary conditions on the cell beat mostly in a synchronous breast stroke at ∼ 50Hz body [26] has little effect on the velocity field as the cell- [3], pulling the cell body through the fluid at speeds drag Stokeslet nearly produces the appropriate velocity U0 ∼ 100 µm/s. Despite its ∼ 5% density excess over field on the cell surface. water [25], the gravitational Stokeslet of Chlamydomonas These flow field measurements around freely-swimming only becomes dominant at distances Λ & 35 R, as its ra- microorganisms provide the basis for a deeper under- tio of U0 to the sedimentation speed (Fg /6πηR) is much standing of a number of issues in biological fluid dynam- larger than for Volvox. Therefore this swimmer has often ics, including the interactions of microorganisms with been modelled as a puller stresslet [14]. A slight three- surfaces, with each other, and the rheology of suspen- dimensionality of the beating causes the cell to spin about sions. For example, it was recently discovered that Volvox its swimming direction at ∼ 2Hz, so ensemble-averaged colonies can form hydrodynamic bound states whose measurements of the kind presented here average out az- properties are quantitatively described by a model of in- imuthal asymmetries in the flow field. Figure 4(a) shows teracting Stokeslets near a no-slip wall [21]. The near that for r & 7R the measured flow topology begins to complete dominance of the flow field around Volvox by resemble a puller stresslet, yet flow speeds at such dis- the Stokeslet term found here provides ex post facto justi- tances are already . 1 µm/s. Closer to the organism, fication for the neglect of higher moments. Perhaps more the field becomes more complex. It includes side vortices importantly this result shows that in terms of interparti- and a flow in front of the cell body that is along the direc- cle hydrodynamic couplings a suspension of Volvox is like tion of motion, towards a stagnation point. The velocity a sedimenting suspension [27], except that the velocity field can be accurately captured, even in the near-field, of each colony is the sum of a self-propelled contribution by modeling the flow created by the pulled cell body as and mutual advection in the flow field of other spheres. a Stokeslet, distributing the thrust among two Stokeslets Elsewhere we illustrate this correspondence in detail [28]. located at the approximate positions of the two flagella, The correspondence between the measured time- and and averaging the flow over one rotation about U0. The azimuthally-averaged flow field of Chlamydomonas and streamlines of this simple extension to the force dipole the three-Stokeslet model illustrates how well such a sim- model [see Fig. 4(b)], as well as the decay of ||u(r)|| with plification captures the complex flow topology, lending 4

FIG. 4: (color online). Time- and azimuthally-averaged flow field of C. reinhardtii. (a) Streamlines (red) computed from velocity vectors (blue). The spiraling near elliptic points is an artifact of the direct integration a noisy experimental velocity field. A color scheme indicates flow speed magnitudes. (b) Streamlines of the azimuthally averaged flow of the three-Stokeslet model: flagellar thrust is distributed among two Stokeslets placed (not fitted) at the approximate flagellar position (lateral green arrows), whose sum balances drag on the cell body (central red arrow). (c) Decay of ||u(r)|| for the three directions indicated by separate colors in the inset, compared to results from the three-Stokeslet model (dashed). support to this approximation in modelling ciliary inter- Rev. Lett. 103, 148101 (2009). actions [29]. Our results indicate that the simple puller- [9] A.P. Berke, et al., Phys. Rev. Lett. 101, 038102 (2008). type description for Chlamydomonas is only valid at dis- [10] G. Li and J.X. Tang, Phys. Rev. Lett. 103, 078101 81 tances & 7 R, where the flow field is already . 1% of (2009); D.G. Crowdy and Y. Or, Phys. Rev. E , 036313 (2010). U0. We then expect interactions with other swimmers, [11] C. Pozrikidis, Boundary Integral and Singularity Meth- boundaries or tracers, to be influenced mostly by the flow ods for Linearized Viscous Flow (Cambridge University structure at shorter separations, where the full time de- Press, Cambridge, 1992). pendence of the flow may be important [30]. We are cur- [12] G.K. Batchelor, J. Fluid Mech. 41, 545 (1970). rently investigating whether similar conclusions hold for [13] R.A. Simha and S. Ramaswamy, Phys. Rev. Lett. 89, the flow field around bacteria, the prototypical “pusher” 058101 (2002); J.P. Hernandez-Ortiz, C.G. Stoltz, M.D. 95 microorganisms. Graham, Phys. Rev. Lett. , 204501 (2005); D. Sain- tillan and M.J. Shelley, Phys. Fluids 20, 123304 (2008); We thank K.C. Leptos for suggesting the use of autoflu- T.J. Pedley, J. Fluid Mech. 647, 335 (2010). orescence to track Chlamydomonas cells, S.B. Dalziel, V. [14] T.J. Pedley and J.O. Kessler, Annu. Rev. Fluid Mech. Kantsler and T.J. Pedley for discussions, D. Page-Croft 24, 313 (1992). and N. Price for technical assistance, and acknowledge [15] C. Brennen, J. Fluid Mech. 65, 799 (1974). support from the EPSRC, the BBSRC, the Marie-Curie [16] T. Ishikawa, M.P. Simmonds, T.J. Pedley, J. Fluid Mech. Program (M.P.), and the Schlumberger Chair Fund. 568, 119 (2006). C.M. Pooley, G.P. Alexander, J.M. Yeo- mans, Phys. Rev. Lett. 99, 228103 (2007). [17] E. H. Harris, The Chlamydomonas Sourcebook (Academic Press, Oxford, 2009), Vol. 1. [18] D.L. Kirk, Volvox (Cambridge University Press, Cam- 182 bridge, 1998). [1] L. Turner, W.S. Ryu, and H.C. Berg, J. Bacter. , [19] D.L. Kirk, M.M. Kirk, Dev. Biol. 96, 493 (1983). 2793 (2000). [20] See movies at http://link.aps.org [2] B.M. Friedrich, I.H. Riedel-Kruse, J. Howard, and F. 102 213 [21] K. Drescher, et al., Phys. Rev. Lett. , 168101 (2009). J¨ulicher, J. Exp. Bio. , 1226 (2010). [22] K. Drescher, et al., Proc. Natl. Acad. Sci. (USA) 107, [3] M. Polin, et al., Science 325, 487 (2009). 9 11171 (2010). [4] C. Brennen and H. Winet, Annu. Rev. Fluid Mech. , [23] M.B. Short, et al., Proc. Natl. Acad. Sci. (USA) 103, 339 (1977); E. Lauga, T.R. Powers, Rep. Prog. Phys. 8315 (2006). 72, 096601 (2009). [24] J.R. Blake, J. Fluid Mech. 46, 199 (1971). [5] C. Dombrowski et al., Phys. Rev. Lett. 93, 098103 [25] C.A. Solari, et al., Am. Nat. 167, 537 (2006). (2004). [26] C.W. Oseen, Hydrodynamik. Leipzig (1927). [6] P.T. Underhill, J.P. Hernandez-Ortiz, M. D. Graham, 28 100 [27] R.E. Caflisch, J.H.C. Luke, Phys. Fluids. , 759 (1985). Phys. Rev. Lett. , 248101 (2008). [28] V. Kantsler, et al., preprint (2010). [7] K. Leptos, et al., Phys. Rev. Lett. 103, 198103 (2009). [29] A. Vilfan and F. J¨ulicher, Phys. Rev. Lett. 96, 058102 [8] Y. Hatwalne et al., Phys. Rev. Lett. 92, 118101 (2004); 588 (2006); T. Niedermayer, B. Eckhardt and P. Lenz, Chaos T. Ishikawa and T.J. Pedley, J. Fluid Mech. , 399 18, 037128 (2008). (2007); A.W.C. Lau and T.C. Lubensky, Phys. Rev. E [30] J.S. Guasto and J.P. Gollub, private communication. 80, 011917 (2009); A. Sokolov and I.S. Aranson, Phys.