Abstract
Self-propelled colloids constitute an important class of intrinsically non-equilibrium matter. Typically, such a particle moves ballistically at short times, but eventually changes its orientation, and displays random-walk behaviour in the long-time limit. Theory predicts that if the velocity of non-interacting swimmers varies spatially in 1D, v(x), then their density Ï(x) satisfies Ï(x)â=âÏ(0)v(0)/v(x), where xâ=â0 is an arbitrary reference point. Such a dependence of steady-state Ï(x) on the particle dynamics, which was the qualitative basis of recent work demonstrating how to âpaintâ with bacteria, is forbidden in thermal equilibrium. Here we verify this prediction quantitatively by constructing bacteria that swim with an intensity-dependent speed when illuminated and implementing spatially-resolved differential dynamic microscopy (sDDM) for quantitative analysis over millimeter length scales. Applying a spatial light pattern therefore creates a speed profile, along which we find that, indeed, Ï(x)v(x)â=âconstant, provided that steady state is reached.
Similar content being viewed by others
Introduction
Einstein predicted and Perrin verified that, in a gravitational field, the equilibrium (number) density of a dilute dispersion of colloidal particles varied with height z according to
where z0 encodes the equality of diffusive and sedimentation fluxes:
with vs a particleâs sedimentation speed and D0 its thermal diffusivity. For spheres of radius a in a liquid of viscosity \(\eta _0\) at temperature T with density Î lower that of the particles \(v_{\mathrm{s}} = 2ga^2{\mathrm{\Delta }}/9\eta _0\) and \(D_0 = k_{\mathrm{B}}T/6\pi \eta _0a\) where kB is Boltzmannâs constant, so that \(z_0 = 3k_{\mathrm{B}}T/4\pi ga^3{\mathrm{\Delta }}\). The verification of Eq. (1) demonstrated the granularity of matter1.
Now suppose \(\eta _0 = \eta _0(z)\), so that the particle dynamics also varies in space: \(D_0 = D_0(z)\). Nevertheless, statistical mechanics stipulates that the spatially dependent dynamical coefficient D0(z) cannot appear in the equilibrium density distribution, and Eq. (1) still holds because the \(\eta\)-dependence cancels out in Eq. (2).
Active colloids2, particles that dissipate energy to propel themselves, form an important class of active matter3,4. Such dissipative states of matter, which include all living organisms, are intrinsically non-equilibrium, and give rise to new physics. Consider a system of run-and-tumble particles (RTPs). An RTP self propels (ârunsâ) at velocity v for time \(\tau _{{\mathrm{run}}}\), then changes direction (âtumblesâ) instantaneously to run at vâ² such that \(|{\mathbf{v}}\prime | = |{\mathbf{v}}|\) but with randomised direction, so that at times \(\gg\)Ïrun it behaves as a random walker. Suppose the run speed of such particles is spatially dependent, v(r). Solving a Fokker-Planck equation for the coarse-grained kinetics of RTPs, Tailleur and Cates5,6 predicted that the resulting density distribution should be \(\rho ({\mathbf{r}}) = \rho (0)v(0)/v({\mathbf{r}})\), with 0 an arbitrary origin. A purely mechanical derivation is also possible7. The appearance of the particlesâ dynamics, v(r), in this formula contrasts starkly with the sedimentation equilibrium of passive colloids with spatially dependent diffusivity D0(x), for which the D-independent Eq. (1) still holds.
When restricted to 1D, this result becomes
which was first derived for non-interacting random walkers by Schnitzer8. Later, Tailleur and Cates showed that it is valid for interacting RTPs whose run speed can be expressed as \(v(x) = v[\rho (x)]\)5. Moreover, under quite general conditions, Ïvâ=âconstant also holds (provided translational diffusion is negligible) for active Brownian particles (ABPs)9, which reorient gradually due to their rotational diffusivity Drot, losing directional memory after a persistence time of ~\(D_{{\mathrm{rot}}}^{ - 1}\). At \(t \gg D_{{\mathrm{rot}}}^{ - 1}\), an ABP is again a random walker.
Qualitatively, Eq. (3) was the basis of recent demonstrations of templated self assembly using light-activated motile bacteria10,11. In a spatially varying illumination pattern, cells accumulate in the darker regions, generating contrast. Quantitatively, however, Eq. (3) has remained unverified by experiments to date. Moreover, its theoretical derivations do not include hydrodynamic interactions, so that its applicability to real systems is open to doubt.
Here, we investigate Eq. (3) quantitatively with the same light-activated Escherichia coli bacteria used previously to demonstrate templated active self assembly10. Our strainâs rapid response to changes in light intensity allows us to impose well defined changes in local swimming speed. Differential dynamic microscopy (DDM)12 has previously been shown to reliably measure swimming speed and density of motile bacteria on uniform samples10,13. Here we develop and demonstrate a spatially resolved implementation (sDDM) for quantitative analysis of \(\bar v(x,y)\) and Ï(x, y) over larger length scales (~mm) with a resolution of â100âμm. We use sDDM to confirm that indeed Ï(x)v(x)â=âconstant, provided that steady state is reached.
Results
Light controlled E. coli and spatially resolved analysis
Each E. coli is an â2âÃâ1âμm spherocylinder with â7â10âμm helical flagella powered by rotary motors14. When all flagella rotate counterclockwise (seen from behind), they bundle and propel the cell. Every Ïrunâ~â1âs or so, one or more flagella reverse and unbundle, causing a change in direction: wild-type (WT) cells are RTPs15. At a typical average speed \(\bar v \gtrsim 10\)âμmâsâ1, they random walk with a persistence length \(l_p \sim \bar v\tau _{{\mathrm{run}}} \sim 10\)âμm. Deleting the cheY gene prevents tumbling; cells become ABPs with \(D_{{\mathrm{rot}}}^{ - 1} \sim 10\)âs and \(l_p \sim \bar vD_{{\mathrm{rot}}}^{ - 1} \sim 100\)âμm16, so that at times \(\gg D_{{\mathrm{rot}}}^{ - 1}\) cells random walk with \(D_{{\mathrm{eff}}} \sim \bar v^2D_{{\mathrm{rot}}}^{ - 1} \sim 10^3\)âμm2âsâ1.
Our E. coli mutants carried a plasmid expressing proteorhodopsin (PR), which pumps protons when exposed to green light17. Cells suspended in nutrient-free motility buffer were sealed into 20âμm high compartments and imaged using 10à phase contrast microscopy. After some minutes, \(\bar v\) dropped abruptly to zero upon oxygen exhaustion10. Thereafter, cells only swam when illuminated in green10,18,19, with an average speed \(\bar v\) that increased with the light intensity, \({\cal{I}}\), saturating at vmax beyond some \({\cal{I}}_0\)10. These are living analogues of synthetic light-activated active colloids20,21.
We first used a digital mirror device10 to project a binary (bright-dark) spatial intensity pattern \({\cal{I}}(x,y)\) spelling out âUoEâ (inset of Fig. 1a) onto a field of cells that had been uniformly illuminated for some time, so \(\bar v\) was initially constant in space. We used sDDM to measure \(\bar v(x,y)\), β(x, y) and Ï(x, y) in 64âÃâ64â(pixel)2 tiles (see Methods section and Supplementary Note 1 for details). The projected \({\cal{I}}(x,y)\) was rapidly replicated in a pattern of \(\bar v(x,y)\) (Fig. 1a). A similar Ï(x, y) pattern soon forms (Fig. 1b) so that this effect can be used for templated self assembly10,11. Given that higher cell densities occur in darker regions with lower swimming speed, Eq. (3) is clearly qualitatively correct10,11. We now proceed to test it quantitatively.
We projected a quasi-1D stepped illumination pattern
on a field of these cells. This generated a spatial pattern of swimming speed, \(\bar v(x,y)\), and cell density, Ï(x, y), which we quantified using sDDM. Averaging over y gives \(\bar v(x)\) and Ï(x), which allows a direct test of Eq. (3), provided that this light pattern (Eq. (4)) generates a corresponding sharp pattern of cell run speeds:
This requires cells that respond rapidly to changes in \({\cal{I}}\), which was found not to be the case10 for previous PR-expressing mutants with otherwise intact metabolism18. Indeed, a recent attempt to verify Eq. (3) using PR-bearing E. coli found instead (in our notation) \(\rho = (a/\bar v) + b\) with positive constants a and b. The latter was ascribed to a long Ïstop, which led to memory effects11.
We achieved rapid response by deleting the unc gene cluster encoding the F1Fo-ATPase membrane protein complex from a parent K12-derived ÎcheY mutant, giving a fast-responding smooth-swimmer, AD1010. In fully-oxygenated phosphate motility buffer (MB) at optical density ODâ=â1, \(\bar v \approx 30\)âμmâsâ1 and a fraction \(\beta \lesssim 20\%\) of cells were non-motile. When illuminated anaerobically, \(\bar v_{{\mathrm{max}}} = 28\)âμmâsâ1 and \(\tau _{{\mathrm{stop}}} \ll 1\)âs, compared with a Ïstop of many minutes in the parent strain without unc deletion10. (Details of other strains we used are given in the methods section.)
Strictly speaking, a non-interacting limit does not exist for bacterial suspensions22. Cells interact hydrodynamically at any concentration, although simulations show that swimmers behave effectively as non-interacting when Ï/\(\rho _c \lesssim 0.1\), where Ïc is the density for the onset of collective behaviour. We observed collective motion in our E. coli suspension at \({\mathrm{OD}} \gtrsim 10\), corresponding to a cell body volume fraction of \(\phi \gtrsim 1.4\%\), consistent with a previous estimate of 2%22,23, so that a quasi-non-interacting limit is reached at \({\mathrm{OD}} \lesssim 1\). It was not possible to work below this limit because of an increasing fraction of cells trapped in circular trajectories (due to hydrodynamic interactions with walls of the sample chamber24) that did not explore the whole sample compartment, hindering relaxation towards a steady state. We therefore worked at ODââ¥â1. We report first data for ODâ=â5 (Ïâââ5âÃâ109âcells/ml; \(\phi \approx 0.7\%\)) before discussing ODâ=â1, where the data are noisier due to lower cell numbers.
Stepped light pattern at ODâ=â5
A field of AD10 cells rendered stationary by oxygen exhaustion was uniformly illuminated for â20âmin to achieve saturation speed10. The light was then attenuated to \({\cal{I}}_ -\), the level of the darker half of the target pattern (Eq. (4)) for 5âmin to determine \(\bar v_ - = 6.5(2)\)âμmâsâ1. Returning the intensity to its initial level, we waited another 5âmin for the swimming speed to return to \(\bar v_ + = 13.2(2)\)âμmâsâ1. We measured the cell density Ï0 and non-motile fraction β0 of this high-speed uniform sample, and then switched on a stepped pattern (Eq. (4)) by reducing the intensity in the xâ>â0 half plane.
Figure 2a shows the mean swimming speed averaged over y tiles, \(\langle \bar v(x,y)\rangle _y = \bar v(x)\), normalised to the whole-sample-averaged speed, \(\langle \bar v\rangle\), plotted against x at 30âmin after switching on the stepped pattern. A stepped speed pattern was developed (Fig. 2a, black circles).
If Eq. (3) is valid, we expect the swimmer density to obey
where â±â subscripts having their obvious meanings. If the density of non-motile cells is constant throughout the experiment (see Supplementary Note 2 for justification), i.e.
we can write the total cell density on the two sides of xâ=â0 as
Finally, the average cell density is
Equations (6)â(9) together predict the density of motile cells in the two half-planes:
We calculated the swimmer density in our experiments from the measured total cell density Ï(x) and non-motile fraction β(x) using \(\rho ^{\mathrm{s}}(x) = \rho (x)[1 - \beta (x)]\), and normalised it by the whole-sample-averaged swimmer density. This function 30âmin after the imposition of the stepped intensity pattern is also stepped (Fig. 2a, blue squares), with the theoretical predictions from Eq. (10), using the measured average v± as inputs (Fig. 2a, dotted line), giving a reasonable account of the step amplitude. A more sophisticated version of this model which takes the measured shape of \(\bar v(x)\) into account (see Supplementary Note 2 for details) is able to capture the amplitude of the jump in Ïs(x) at xâ=â0 even more precisely (Fig. 2a, dashed line).
The product \(\rho ^{\mathrm{s}}\bar v(x)\) normalised to the whole sample average (Fig. 2a, red xâs), is indeed constant for \(|x| \lesssim 200\)âμm, verifying Eq. (6), which is the application of Eq. (3) to our conditions. However, there are systematic deviations from constancy at \(|x| \gtrsim 200\)âμm. One possible explanation is the emergence of collective motion with associated local nematic ordering25, which would invalidate the derivation of Eq. (3). However, we only observed collective motion at \({\mathrm{OD}} \gg 10\). Instead, the deviation of \(\rho ^{\mathrm{s}}\bar v\) from constancy at \(|x| \gtrsim 200\)âμm is a kinetic effect. Figure 2b shows the time evolution of the normalised \(\rho ^{\mathrm{s}}\bar v(x)\) at different x. Steady state was reached rapidly for \(|x| \lesssim 200\)âμm, but was not reached by 30âmin at the extremes of our observation window, \(|x| \gtrsim 600\)âμm. Given their effective diffusivity Deffâ~â103âμm2âsâ1, cells at the extremities of our compartment take \(\gg\)30âmin to sufficiently sample both speed regions, preventing the attainment of steady state within our observational time window. This leads to the deviations between observed and predicted Ïs(x) away from xâ=â0. Nevertheless, Fig. 2b suggests that \(\rho ^{\mathrm{s}}\bar v\)â=âconstant should be attained at all x in the long-time limit.
Stepped pattern at other cell densities
Measurements and model predictions for the lower ODâ=â1 are shown in Fig. 3. The data are noisier, but show the same trends. In the vicinity of xâ=â0, \(\rho ^{\mathrm{s}}\bar v\)âââconstant. To highlight the behaviour in the two 90âμm-wide stripes of tiles bordering xâ=â0, we plot \(\rho ^{\mathrm{s}}\bar v(x)\) at tâ=â30âmin for these two stripes against each other for a number of independent experiments (Fig. 4, black solid circles). In all cases, \(\rho ^{\mathrm{s}}\bar v(x)\)â=âconstant for these central stripes to within experimental uncertainties. The ratio of \(\rho ^{\mathrm{s}}\bar v(x)\) in these two stripes plotted against the ratio of the swimming speed on the two half-planes (Fig. 4 inset, black solid circles) is consistent with this claim.
We performed experiments using the stepped light pattern at other cell densities and also using an additional smooth swimming strain (DM1). In all cases up to ODâ=â8, we find that \(\rho ^{\mathrm{s}}\bar v (x)\) is constant across the two central stripe of tiles on either side of xâ=â0 (Fig. 4), where we are certain that a steady state has been reached, verifying Eq. (3) up to Ïâââ8âÃâ109âcells/ml. Spatial maps and time evolution of spatial profiles for our highest density sample (ODâ=â8) are shown in Supplementary Figs. 6 and 7, respectively.
Measurements using a periodic light pattern
The complicating factor so far is the slow global convergence towards Ïs(x)v(x)â=âconst, so that steady state will only be reached in approximately hours. Experiments on such time scales are impractical due to mechanical and biological stability issues. Thus, we only have direct evidence for the validity of Eq. (3) in the vicinity of the intensity step at xâ=â0. This suggests that the use of a series of thin stripes would give more clear-cut results unencumbered by kinetic issues. We found that this was indeed the case.
In response to the imposition of a one-dimensional square-wave illumination pattern of brighter-darker stripes with 540âμm repeat generated by a digital mirror device10, the swimming speed of bacteria changed from a uniform distribution to a square-wave distribution almost instantaneously (in \(\lesssim\)1âs) (Fig. 5a). This in turn modified the cell density (initially uniform at ODâ=â1), which approached a steady state much more quickly. This is possible not only because of the length scale reduction, but also because swimmers can enter (say) a high-intensity region from low-intensity regions on both the left and right.
Figure 5b shows the normalised swimmer density \(\rho ^{\mathrm{s}}(x)\)/\(\langle \rho ^{\mathrm{s}}\rangle\) after 15âmin of patterned illumination, together with \(\bar v(x)\)/\(\langle v\rangle\) and their product. While the data are again somewhat noisy because of the low average cell density (ODâ=â1), it is clear that \(\rho ^{\mathrm{s}}\bar v (x)\)â=âconstant to within one standard deviation, which directly verifies Eq. (3).
Experiments with β(I) dependency
Interestingly, experiments using low light intensities (which gave low swimming speeds) proved less successful, because at very low intensities we found a noticeably higher percentage of non-motile cells in the sample than at high light intensities (Supplementary Fig. 8). This led to a spatial variation in the non-motile density (Supplementary Fig. 9), which considerably complicates the interpretation and analysis of such experiments (see Supplementary Note 3 for details).
Experiments using light-powered run-and-tumble strain
We end by explaining why we did not use motility wild-type (run and tumble) strains for our experiments. Their motion randomises much more rapidly than smooth swimming mutants, which would have significantly alleviated the non-steady-state issue for the stepped intensity pattern. However, we found that AD4, a PR-bearing motility wild type, gathered near xâ=â0, on the darker side of the intensity step (see Supplementary Figs. 10 and 11 and Supplementary Note 4 for details). From the q-dependence of our DDM data13 we can deduce that the tumbling rate increases noticeably as cells swim from light to dark, whereas cells swimming from dark to light do not show any obvious change in their tumbling behaviour. This may be due to âenergy taxisâ26. The validity of Eq. (3) depends on the assumption that the tumbling rate is independent of swimming direction5,8, so that motility wild types cannot be used to test this result.
Discussion
Equation (3) is one of only a handful of exact predictions to date on the statistical mechanics of active particle systems. Its âweakâ form, for non-interacting systems was derived for RTPs8, while its âstrongâ form was later derived both for RTPs5 and ABPs9. Taken together, our experiments using stepped and stripped light patterns give strong evidence that Eq. (3) holds at 1ââ¤âODââ¤â8 (\(0.15\% \lesssim \phi \lesssim 1.2\%\)) for smooth swimming E. coli whenever we can be sure that steady state has been reached, either in the vicinity of xâ=â0 in the stepped pattern or throughout the stripped pattern. Our swimmers are interacting throughout our concentration range22, even though collective motion is not observed until ODâ>>â10. Thus, our results verify the âstrongâ form of Eq. (3) for ABPs.
The qualitative validity of Eq. (3), viz., that cells gather where they swim slower, or, equivalently for our cells, where the light intensity is lower, has already been assumed and utilised in recent work deploying such cells in smart (or reconfigurable) templated self assembly, or âpainting with bacteriaâ10,11. Indeed, in a recent demonstration of how to perform bacterial painting with multiple shades of graded intensity levels11, there was attempt at checking the correctness of Eq. (3) en passant, which, however, was unsuccessful because of a high number of non-motile cells and the long stopping time of their strain, the latter producing memory effects. Our success in verifying Eq. (3) shows that carefully quantifying and subtracting the non-motile fraction and the use of a strain of bacteria with very short stopping time are essential ingredients in such an experiment. Indeed, without careful design most ârealâ active systems are likely to display dynamic behaviour that is too complex to fulfil the assumptions leading to Eq. (3), as evidenced by our findings for the wild type strain with illumination-dependent tumbling rate. Note, however, that âsmart templated active self assemblyâ10 using photo-activated swimmers is possible for any relationship of Ïâv in which the increase in one variable necessitates the decrease in the other.
Our experiments would not have been possible without spatially resolved differential dynamic microscopy (sDDM), which can reliably quantify swimming speed and relative density (along with many other parameters) over a wide range of length scales and cell concentrations. Its adoption can therefore provide new insights into a variety of systems displaying spatially varying dynamics, from biological taxis27 to collective motion.
Throughout, we have focussed on steady-state effects, although the consideration of time dependence proved crucial in interpreting apparent systematic deviations from the prediction of Eq. (3) for imposed stepped intensity patterns. Time-dependent effects are, of course, interesting in their own right. Thus, the response of active particles to a time-dependent topographic landscape that is self-assembled by the cells themselves9 has yet to be explored experimentally. On the other hand, it has recently been suggested theoretically28 and demonstrated experimentally29 that travelling-wave light fields can be exploited for transporting and rectifying light-activated swimmers. Exploitation of these and other opportunities should open up new fields of fundamental studies and applications.
Methods
Strains and sample preparation
We constructed three different strains of E. coli using plasmids expressing SAR86 γ-proteorhodopsin (a gift from Jan Liphardt, UC Berkley). These strains were designed to exhibit a fast response to changes in light intensity. This was achieved by deleting the unc gene cluster, so that the F1F0-ATPase membrane protein complex cannot work in reverse in the dark to generate proton motive force to power swimming. The detailed molecular biology and strain characterisation have been reported before10. AD4 is a WT (run-and-tumble) swimmer derived from AB1157, whereas DM1 and AD10 are smooth swimming strains derived from RP437 and AB1157, respectively (see Supplementary Table 1). The two smooth swimming strains behaved similarly, although AD10 achieved a much higher swimming speed than DM1 and was also more efficiently powered by light. Therefore we mostly used AD10, with some additional data acquired using DM1.
Overnight cultures were grown aerobically in 10âmL Luria-Bertani Broth (LB) using an orbital shaker at 30â°C and 200ârpm. A fresh culture was inoculated as 1:100 dilution of overnight grown cells in 35âml tryptone broth (TB) and grown for 4âh to an optical density of OD600âââ0.2. The production of proteorhodopsin (PR) was induced by adding arabinose to a concentration of 1âmM as well as the necessary cofactor all-trans-retinal to 10âμM to the growth medium. Cells were incubated under the same conditions for a further hour to allow protein expression to take place and then transferred to motility buffer (MB, pHâ=â7.0, 6.2âmM K2HPO4, 3.8âmM KH2PO4, 67âmM NaCl and 0.1âmM EDTA). Single filtration (0.45âμm HATF filter, Millipore) was used to prepare high density stock solutions (ODâââ8) which were diluted with MB to the desired cell concentration.
The samples were loaded into commercial 2âμL sample chambers (SC-20-01-08-B, Leja, NL) of dimensions ââ6âÃâ10âmmâÃâ20âμm, where cells predominantly swim in the (x, y) (imaging) plane, but have enough room to âovertakeâ each other in all three spatial dimensions. The chamber was then sealed using vaseline to stop air flow, so that swimming stopped once dissolved oxygen was exhausted19. This happened within 10âmin at ODâââ1 (â109âcells/ml or 0.2% volume fraction of cell bodies). Thereafter, we controlled the activity of the cells by illuminating with green light of various intensities10.
Experimental setup
The samples were observed using a Nikon TE2000 inverted microscope with a PF 10Ã, N.A. 0.3 phase contrast objective. Time series of movies (~40âs duration at 100âframes per second) were recorded using a CMOS camera (MC 1362, Mikrotron). A long-pass filter (RG630, Schott Glass) in the bright-field light path ensured that the imaging light did not activate PR. The light controlling bacterial swimming was provided by an LED (Sola SE II, Lumencor) whose intensity was set via a computer interface. The LED light was filtered to a green wavelength range (510â560ânm) overlapping with the absorption peak of our PR18 and illuminated the sample in a trans-illumination geometry. By illuminating an area much larger than the field of view of the objective, we minimised the loss of swimmers over time. If only a small region of the sample is illuminated, the density of swimmers continuously drops, because they reach the illumination boundaries and accumulate there (no lightâ=âno swimming). For the stepped pattern experiment we uniformly illuminated a â7âmm diameter circle, covering almost all of the sample chamber. Under these conditions the cell density is conserved, thus simplifying theoretical modelling. We used a thin sheet polariser imaged onto the sample plane to attenuate the intensity on half of the sample. A digital mirror device10 projected the periodic pattern onto a â2.9âmm diameter area of the sample.
Differential dynamic microscopy
DDM measures \((\bar v,\beta )\) averaged over 104â105âcells under our conditions13,19,30. From â40âs of wide-field, low-magnification movies, one extracts the power spectrum of the difference between pairs of images delayed by time Ï, \(g(\vec q,\tau )\), where \(\vec q\) is the spatial frequency vector. Under suitable conditions and for isotropic motion, the intermediate scattering function f(q, Ï), the qth mode of the density autocorrelation function, is given by
Here, B(q) relates to the background noise and A(q) is the signal amplitude. Fitting f(q, Ï) to a suitable swimming model of E. coli yields four key motility parameters: the mean \(\bar v(q)\), and width Ï(q) of the speed distribution P(v), the non-motile fraction β(q), and the diffusion coefficient of non-motile cells D(q), as a function of q. All of these should, ideally, be q-independent. In practice, there is some q-variation. We typically averaged the fitting parameters over \(0.5 < q < 2.2\,\upmu {\mathrm{m}}^{ - 1}\) to give, e.g. \(\bar v = \langle v(q)\rangle _q\) and \(\beta = \langle \beta (q)\rangle _q\).
In a dilute system whose structure factor \(S(q) \approx 1\), A(q) is proportional to the sample density31,32 and can therefore be used to determine relative densities by \(\rho _1/\rho _0 = \langle A_1(q)/A_0(q)\rangle _q\)10. Note that ratioing the A(q)s removes their strong q-dependence.
Spatially resolved DDM (sDDM) is in principle straightforward: the above algorithm simply needs to be implemented on pâÃâpâ(pixel)2 sub-movies. In practice, care is required in choosing the minimum p for which meaningful results can be obtained. We do this by illuminating a field of cells uniformly, measuring \((\bar v,\beta ,\rho ,\rho ^{\mathrm{s}}v)\) from individual pâÃâp tiles in the steady state, and obtaining the probability distribution of these parameters. Under our imaging conditions, we found that these distributions became p-independent when \(p \ge 64\). We therefore chose \(p = 64\), corresponding to 90âμm in the sample (see Supplementary Note 1 for details).
A full 512âÃâ512 movie yields g(q, Ï)s at 512/2â=â256 distinct q values. We divide it into 64 sub-movies of size 64âÃâ64â(pixel)2. This yields 8âÃâ8âÃâ(64/2)â=â2048âg(q, Ï)s to be fitted to give for each sub-movie \(v_{x,y}(q)\), \(\beta _{x,y}(q)\) and \(\rho _{x,y}(q)/\rho _0 = A(q)_{x,y}/A_0(q)\), where A0(q) is measured from the same sample under uniform illumination (i.e. just before switching to a structured light pattern). These were averaged over \(0.5 \le q \le 1.5\,\upmu {\mathrm{m}}^{ - 1}\). The upper q limit is somewhat lower than what is typical for whole-movie analysis13 due to non-systematic failure of fitting at higher q values, presumably due to noise or windowing artefacts33.
Data availability
The data presented here is available on the Edinburgh DataShare repository34.
References
Perrin, J. Atoms (Ox Bow Press, Woodbridge, 1990). (Original: Les Atomes, Librairie Félix Alcan, Paris, 1913).
Poon, W. C. K. in Physics of Complex Colloids (eds Bechinger, C., Sciortino, F. & Ziherl, P.) 317â386 (Società Italiana di Fisica, Bologna, 2013).
Ramaswamy, S. The mechanics and statistics of active matter. Annu. Rev. Condens. Matter Phys. 1, 323â345 (2010).
Marchetti, M. C. et al. Hydrodynamics of soft active matter. Rev. Mod. Phys. 85, 1143 (2014).
Tailleur, J. & Cates, M. E. Statistical mechanics of interacting run-and-tumble bacteria. Phys. Rev. Lett. 100, 218103 (2008).
Tailleur, J. & Cates, M. E. Sedimentation, trapping and rectification of dilute bacteria. Europhys. Lett. 86, 60002 (2009).
Takatori, S. C. & Brady, J. F. Towards a thermodynamics of active matter. Phys. Rev. E 91, 032117 (2015).
Schnitzer, M. J. Theory of continuum random walks and application to chemotaxis. Phys. Rev. E 48, 2553â2568 (1993).
Stenhammar, J., Wittkowski, R., Marenduzzo, D. & Cates, M. E. Light-induced self-assembly of active rectification devices. Sci. Adv. 2, e1501850 (2016).
Arlt, J., Martinez, V. A., Dawson, A., Pilizota, T. & Poon, W. C. K. Painting with lightâpowered bacteria. Nat. Commun. 9, 768 (2018).
Frangipane, G. et al. Dynamic density shaping of photokinetic E. coli. eLife 7, 1â6 (2018).
Cerbino, R. & Trappe, V. Differential dynamic microscopy: probing wave vector dependent dynamics with a microscope. Phys. Rev. Lett. 100, 188102 (2008).
Martinez, V. A. et al. Differential dynamic microscopy: a high-throughput method for characterizing the motility of microorganisms. Biophys. J. 103, 1637â1647 (2012).
Berg, H. C. E. coli in Motion. (Springer, Berlin, 2003).
Berg, H. C. & Brown, D. A. Chemotaxis in Escherichia coli analysed by three-dimensional tracking. Nature 239, 500â504 (1972).
Wu, M., Roberts, J. W., Kim, S., Koch, D. L. & DeLisa, M. P. Collective bacterial dynamics revealed using a three-dimensional population-scale defocused particle tracking technique. Appl. Environ. Microbiol. 72, 4987â4994 (2006).
Béjà , O. et al. Bacterial rhodopsin: evidence for a new type of phototrophy in the sea. Science 289, 1902â1906 (2000).
Walter, J. M., Greenfield, D., Bustamante, C. & Liphardt, J. Light-powering Escherichia coli with proteorhodopsin. Proc. Natl Acad. Sci. USA 104, 2408â2412 (2007).
Schwarz-Linek, J. et al. Escherichia coli as a model active colloid: a practical introduction. Colloids Surf. B 137, 2â16 (2016).
Buttinoni, I., Volpe, G., Kümmel, F., Volpe, G. & Bechinger, C. Active Brownian motion tunable by light. J. Phys. Condens. Matter 24, 284129 (2012).
Palacci, J., Sacanna, S., Steinberg, A. P., Pine, D. J. & Chaikin, P. M. Living crystals of light-activated colloidal surfers. Science 339, 936â940 (2013).
Stenhammar, J., Nardini, C., Nash, R. W., Marenduzzo, D. & Morozov, A. Role of correlations in the collective behavior of microswimmer suspensions. Phys. Rev. Lett. 119, 028005 (2017).
Gachelin, J., Rousselet, A., Lindner, A. & Clement, E. Collective motion in an active suspension of Escherichia coli bacteria. New J. Phys. 16, 025003 (2014).
Lauga, E., Di Luzio, W. R., Whitesides, G. M. & Stone, H. A. Swimming in circles: motion of bacteria near solid boundaries. Biophys. J. 90, 400â412 (2006).
Cisneros, L. H., Kessler, J. O., Ganguly, S. & Goldstein, R. E. Dynamics of swimming bacteria: transition to directional order at high concentration. Phys. Rev. E 83, 061907 (2011).
Schweinitzer, T. & Josenhans, C. Bacterial energy taxis: a global strategy? Arch. Microbiol. 192, 507â520 (2010).
Rosko, J., Martinez, V. A., Poon, W. C. K. & Pilizota, T. Osmotaxis in Escherichia coli through changes in motor speed. Proc. Natl Acad. Sci. USA 114, E7969âE7976 (2017).
Maggi, C., Angelani, L., Frangipane, G. & Di Leonardo, R. Currents and flux-inversion in photokinetic active particles. Soft Matter. 12, 1â5 (2018).
Koumakis, N. et al. Self-assembly and transport of motile bacteria with time-dependent light fields. https://arxiv.org/abs/1811.09133 (2018).
Wilson, L. G. et al. Differential dynamic microscopy of bacterial motility. Phys. Rev. Lett. 106, 018101 (2011).
Reufer, M., Martinez, V. A., Schurtenberger, P. & Poon, W. C. K. Differential dynamic microscopy for anisotropic colloidal dynamics. Langmuir 28, 4618â4624 (2012).
Lu, P. J. et al. Characterizing concentrated, multiply scattering, and actively driven fluorescent systems with confocal differential dynamic microscopy. Phys. Rev. Lett. 108, 218103 (2012).
Giavazzi, F., Edera, P., Lu, P. J. & Cerbino, R. Image windowing mitigates edge effects in Differential Dynamic Microscopy. Eur. Phys. J. E 40, 1â8 (2017).
Arlt, J., Martinez, V. A., Dawson, A., Pilizota, T. & Poon, W. C. K. Dynamics-dependent density distribution in active suspensions. https://doi.org/10.7488/ds/2550 (University of Edinburgh, 2019).
Acknowledgements
The authors were funded by the EPSRC (EP/J007404/1) and the ERC (AdG 340877 PHYSAPS). We thank Mike Cates, Julien Tailleur, Alexander Morozov and Nick Koumakis for helpful discussions, Jan Liphardt for a gift of PR plasmids and Dario Miroli for E. coli strain DM1.
Author information
Authors and Affiliations
Contributions
J.A. and V.A.M. contributed equally to this work. W.C.K.P. initiated the work. T.P. and A.D. designed mutants constructed by A.D. J.A. and V.A.M. performed experiments, analysed and interpreted data with T.P. and W.C.K.P., and wrote manuscript with W.C.K.P.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Journal peer review information: Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Arlt, J., Martinez, V.A., Dawson, A. et al. Dynamics-dependent density distribution in active suspensions. Nat Commun 10, 2321 (2019). https://doi.org/10.1038/s41467-019-10283-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-019-10283-0
This article is cited by
-
Colloidal transport by light induced gradients of active pressure
Nature Communications (2023)
-
Rectification and confinement of photokinetic bacteria in an optical feedback loop
Nature Communications (2022)
-
Time irreversibility in active matter, from micro to macro
Nature Reviews Physics (2022)
-
Feedback-controlled active brownian colloids with space-dependent rotational dynamics
Nature Communications (2020)