Dynamics-dependent density distribution in active suspensions
Jochen Arlt Vincent A Martinez Angela Dawson Teuta Pilizota Wilson C K Poon







j.arlt@ed.ac.uk
Corresponding author
School of Physics & Astronomy, The University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom


School of Physics & Astronomy, The University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom


School of Physics & Astronomy, The University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom


School of Biological Sciences and Centre for Synthetic and Systems Biology, The University of Edinburgh, Alexander Crum Brown Road, Edinburgh EH9 3FF, United Kingdom


School of Physics & Astronomy, The University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom



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 behavior 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. We verify this prediction quantitatively by constructing bacteria that swim with an intensity-dependent speed when illuminated. 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.  






Dynamics-dependent density distribution in active suspensions
    Wilson C K Poon
    today
=============================================================



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 

    ρ (z) = ρ(0) e^-z/z_0,

where z_0 encodes the equality of diffusive and sedimentation fluxes:

    z_0 = D_0/v_ s,

with v_ s a particle's sedimentation speed and D_0 its thermal diffusivity. For spheres of radius a in a liquid of viscosity η_0 at temperature T with density Δ lower that of the particles v_ s = 2ga^2Δ/9η_0 and
D_0 = k_ BT/6πη_0 a where k_ B is Boltzmann's constant, so that z_0 = 3k_ BT/4π ga^3 Δ. The verification of Eq. <ref> demonstrated the granularity of matter <cit.>. 

Now suppose η_0 = η_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 D_0(z) cannot appear in the equilibrium density distribution, and Eq. <ref> still holds because the η-dependence cancels out in Eq. <ref>. 
 
Active colloids <cit.>, particles that dissipate energy to propel themselves, form an important class of active matter <cit.>. 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 𝐯 for time τ_ run, then changes direction (`tumbles') instantaneously to run at 𝐯^' such that |𝐯^'| = |𝐯| but with randomised direction, so that at times ≫τ_ run it behaves as a random walker. Suppose the run speed of such particles is spatially dependent, v(𝐫). Solving a Fokker-Planck equation for the coarse-grained kinetics of RTPs, Tailleur and Cates <cit.> predicted that the resulting density distribution should be ρ(𝐫) = ρ(0)v(0)/v(𝐫), with 0 an arbitrary origin. (A purely mechanical derivation is also possible <cit.>.) The appearance of the particles' dynamics, v(𝐫), in this formula contrasts starkly with the sedimentation equilibrium of passive colloids with spatially-dependent diffusivity D_0(x), for which the D-independent Eq. <ref> still holds. 

When restricted to 1D, this result becomes

    ρ(x) = ρ(0)v(0)/v(x),

which was first derived for non-interacting random walkers by Schnitzer <cit.>. Tailleur and Cates show that it is valid for interacting RTPs whose run speed can be expressed as v(x) = v[ρ(x)] <cit.>. Moreover, under quite general conditions, ρ v = constant also holds (provided translational diffusion is negligible) for active Brownian particles (ABPs) <cit.>, which reorient gradually due to their rotational diffusivity D_ rot, losing directional memory after a persistence time of ∼ D_ rot^-1. At  t ≫ D_ rot^-1, an ABP is again a random walker. 

Qualitatively, Eq. <ref> was the basis of recent demonstrations of templated self assembly using light-activated motile bacteria <cit.>. In a spatially-varying illumination pattern, cells accumulate in the darker regions, generating contrast. Quantitatively, however, Eq. <ref> has remained unverified by experiments to date. 
Indeed, as its theoretical derivations do not explicitly take account of hydrodynamic interactions, it is unclear to what extent it is applicable to such systems.
Here, we quantitatively investigate eq. <ref> with the same light-activated Escherichia coli bacteria used previously to demonstrate templated active self assembly <cit.>.

Each E. coli is an ≈2×1 spherocylinder with ≈ 7–10 helical flagella powered by rotary motors <cit.>. When all flagella rotate counterclockwise (seen from behind), they bundle and propel the cell. Every τ_ run∼1 or so, one or more flagella reverse and unbundle, causing a change in direction: wild-type (WT) cells are RTPs <cit.>. At a typical average speed v̅≳10, they random walk with a persistence length l_p ∼v̅τ_ run∼10. Deleting the cheY gene prevents tumbling; cells become ABPs with D_ rot^-1∼10 and l_p ∼v̅ D_ rot^-1∼100 <cit.>, so that at times ≫ D_ rot^-1 cells random walk with D_ eff∼v̅^2 D_ rot^-1∼ 10^3 □. 

Our E. coli mutants carried a plasmid expressing proteorhodopsin (PR), which pumps protons in green light <cit.>. Cells suspended in nutrient-free motility buffer were sealed into 20 high compartments and imaged using 10× phase contrast microscopy.  After some minutes, v̅ dropped abruptly to zero upon oxygen exhaustion <cit.>. Thereafter, cells only swam when illuminated in green  <cit.>, with an average speed v̅ that increased with the light intensity, ℐ, saturating at v_ max beyond some ℐ_0 <cit.>. These are living analogs of synthetic light-activated active colloids <cit.>.

We projected a quasi-1D stepped illumination pattern 

    ℐ(x,y) = {[  ℐ_+ < ℐ_0,   x<0; ℐ_- < ℐ_+,   x >0 ].

on a field of these cells. This generated a spatial pattern of swimming speed, v̅(x,y), and cell density, ρ(x,y), which we measured using spatially-resolved differential dynamic microscopy (sDDM) <cit.>. Averaging over y gives v̅(x) and ρ(x), which allows a direct test of Eq. <ref>, provided that this light pattern, Eq. <ref>, generates a corresponding sharp pattern of cell run speeds:

    v̅(x) = {[        v̅_+  ,  x<0; v̅_- < v̅_+,   x >0 ]. .

This requires cells that respond rapidly to changes in ℐ, which was found not to be the case <cit.> for previous PR-expressing mutants with otherwise intact metabolism <cit.>. Indeed, a recent attempt to verify Eq. <ref> using PR-bearing E. coli found instead (in our notation) ρ = (a/v̅) + b with positive constants a and b. The latter was ascribed to a long τ_ stop, which led to memory effects <cit.>. 






We achieved rapid response by deleting the unc gene cluster encoding the F_1F_o-ATPase membrane protein complex from a parent K12-derived ΔcheY mutant, giving a fast-responding smooth-swimmer, AD10 <cit.>. 
In fully-oxygenated phosphate motility buffer (MB) at optical density OD=1, v̅≈30 and a fraction β≲ 20% of cells were non-motile. When illuminated anaerobically, v̅_ max = 28 and τ_ stop≪1, compared to a τ_ stop of many minutes in the parent strain without unc deletion <cit.>. (Details of other strains we used are given in the methods section.) 





§ RESULTS


We first used a digital mirror device <cit.> to project a binary (bright-dark) spatial intensity pattern ℐ(x,y) spelling out `UoE', Fig. <ref>(a, inset), onto a field of cells that had been uniformly illuminating for some time, so v̅ was initially constant in space. We used sDDM to measure v̅(x,y), β(x,y) and ρ(x,y) in 64× 64 (pixel)^2 tiles (see SI for details). The projected ℐ(x,y) was rapidly replicated in a pattern of v̅(x,y), Fig. <ref>(a). A similar ρ(x,y) pattern soon forms, Fig. <ref>(b), so that this effect can be used for templated self assembly <cit.>. Given that higher cell densities occur in darker regions with lower swimming speed, Eq. <ref> is clearly qualitatively correct <cit.>. We now proceed to test it quantitatively.

Strictly speaking, a non-interacting limit does not exist for bacterial suspensions <cit.>. Cells interact hydrodynamically at any concentration, although simulations show that swimmers behave effectively as non-interacting when ρ/ρ_c ≲ 0.1, where ρ_c is the density for the onset of collective behaviour. We observed collective motion in our E. coli suspension at OD ≳ 10, corresponding to a cell body volume fraction of ϕ≳ 1.4%, consistent with a previous estimate of 2% <cit.>, so that a quasi-non-interacting limit is reached at OD ≲ 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 chambers, see <cit.>) 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× 10^9 cells/ml; ϕ≈ 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 to achieve saturation speed  <cit.>. The light was then attenuated to ℐ_-, the level of the darker half of the target pattern, Eq. <ref>, for 5 to determine v̅_- = 6.5(2). Returning the intensity to its initial level, we waited another 5 for the swimming speed to return to v̅_+ = 13.2 ± 0.2. 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. <ref>, by reducing the intensity in the x > 0 half plane. 






Figure <ref>(a) shows the mean swimming speed averaged over y tiles, ⟨v̅(x,y) ⟩_y = v̅(x), normalised to the whole-sample-averaged speed, ⟨ v ⟩, plotted against x at 30 after switching on the stepped pattern. A stepped speed pattern developed, Fig. <ref>a (∘).

If Eq. <ref> is valid, we expect the swimmer density to obey

    ρ^ s_+(x) v̅_+(x) = ρ^ s_-(x) v̅_-(x),

where `±' subscripts having their obvious meanings. If the density of non-motile cells is constant throughout the experiment (see SI for justification), i.e.

    ρ_±^ nm = β_0ρ_0,

we can write the total cell density on the two sides of x= 0 as

    ρ_± = ρ_±^ s  + ρ_±^ nm = ρ_±^ s  + β_0 ρ_0.

Finally, the average cell density is  

    ρ_0 = 1/2(ρ_+ + ρ_-).

Equations <ref>-<ref> together predict the density of motile cells in the two half planes: 

    ρ_±^ s = 2ρ_0(1-β_0)/1 + v_±/v_∓.











We calculated the swimmer density in our experiments from the measured total cell density ρ(x) and non-motile fraction β(x) using ρ^ s(x) = ρ(x)[1-  β(x)], and normalised it by the whole-sample-averaged swimmer density. This data 30 after the imposition of the stepped intensity pattern is also stepped, Fig. <ref>a ( □), with the theoretical predictions from Eq. <ref> using the measured average v_± as inputs, Fig. <ref>a (dotted line), giving a reasonable account of the step amplitude. A more sophisticated version of this model which takes the measured shape of v̅(x) into account (see SI for details) is able to capture the amplitude of the jump in ρ^ s(x) at x = 0 even more precisely, Fig. <ref>a (dashed line).


The product ρ^ sv̅(x) normalised to the whole sample average, Fig. <ref>a ( ×), is indeed constant for |x| ≲200, verifying Eq. <ref>, which is the application of Eq. <ref> to our conditions. 
However, there are systematic deviations from constancy at |x| ≳200.  
One possible explanation is the emergence of collective motion with associated local nematic ordering <cit.>, which would invalidate the derivation of Eq. <ref>. However, we only observed collective motion at OD ≳ 10. 
Instead, the deviations of ρ^ sv̅ from constancy at |x| ≳200 is a kinetic effect. Figure <ref>b shows the time evolution of the normalised ρ^ sv̅(x) at different x. Steady state was reached rapidly for |x| ≲200, but was not reached by 30 at the extremes of our observation window, |x| ≳600. Given their effective diffusivity D_ eff∼ 10^3 □, cells at the extremities of our compartment take ≫30 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, Figure <ref>b suggests that ρ^ sv̅ = 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. <ref>(a, b). The data are noisier, but show the same trends. In the vicinity of x = 0, ρ^ sv̅≈ constant. To highlight the behaviour in the two 90-wide stripes of tiles bordering x = 0, we plot  ρ^ sv̅(x) at t = 30 for these two stripes against each other for a number of independent experiments, Fig. <ref> (∙). In all cases, ρ^ sv̅(x) = constant for these central stripes to within experimental uncertainties. The ratio of ρ^ sv̅(x) in these two stripes plotted against the ratio of the swimming speed on the two half planes, Fig. <ref> inset (∙), 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 ρ^ sv̅ is constant across the two central stripe of tiles on either side of x = 0, Fig. <ref>, where we are certain that steady state has been reached, verifying Eq. <ref> up ρ≈ 8 × 10^9 cells/ml. 






 §.§ 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 ∼ 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. <ref> 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 repeat generated by a digital mirror device <cit.>, the swimming speed of bacteria changed from a uniform distribution to a square-wave distribution almost instantaneously (in ≲1), Fig. <ref>(a). 
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.

Fig. <ref>b shows the normalised swimmer density ρ^ s(x)/⟨ρ^ s⟩ after 6 of patterned illumination, together with v̅(x)/⟨ v ⟩ and their product. While the data are  again somewhat noisy because of the low average cell density (OD = 1), it is clear that ρ^ sv̅= constant to within one standard deviation, which directly verifies Eq. <ref>. 









 §.§ 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 (Fig. <ref>).
This led to a spatial variation in the non-motile density (Fig. <ref>), which considerably complicates the interpretation and analysis of such experiments (see SI 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 Figs. <ref> & <ref>). From the q dependence of our DDM data <cit.> 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 behavior. This may be due to `energy taxis' <cit.>. The validity of Eq. <ref> depends on the assumption that the tumbling rate is independent of swimming direction <cit.>, so that motility wild types cannot be used to test this result.






§ SUMMARY AND CONCLUSIONS


Equation <ref> 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 RTPs <cit.>, while its `strong' form was later derived both for RTPs <cit.> and ABPs <cit.>. Taken together, our experiments using stepped and stripped light patterns give strong evidence that Eq. <ref> holds at 1 ≤ OD ≤ 8 (0.15%≲ϕ≲ 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 range <cit.>, even though collective motion is not observed until OD = 10. Thus, our results verify the `strong' form of Eq. <ref> for ABPs. 
 


The qualitative validity of Eq. <ref>, 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' <cit.>. Indeed, in a recent demonstration of how to perform bacterial painting with multiple shades of graded intensity levels <cit.>, there was attempt at checking the correctness of Eq. <ref> 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. <ref> 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 which is too complex to fulfill the assumptions leading to eq. <ref>, as evidenced by our findings for the run-and-tumble strain.
However, our experimental method of spatially resolved DDM can reliably quantify swimming speed and relative density (along with many other parameters) even in such cases. As such it can provide new insights into a wide variety of sytems displaying spatially varying dynamics, from biological taxis studies <cit.> 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. <ref> 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 themselves <cit.> has yet to be explored experimentally. On the other hand, it has recently been suggested theoretically <cit.> and demonstrated experimentally <cit.> 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. 





§ MATERIALS AND METHODS







 §.§ Cells

We constructed 3 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 F_1F_0-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 before <cit.>. 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 SI table <ref>).

The 2 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 35ml tryptone broth (TB) and grown for 4 to an optical density of OD_600≈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 K_2HPO_4, 3.8 mM KH_2PO_4, 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 sample chambers (SC-20-01-08-B, Leja, NL) of dimensions ≈6×10×20, 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 exhausted <cit.>. This happened in ≈ 20 min at OD≈ 1 (≈ 10^9 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 intensities <cit.>. 




 §.§ 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 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) overlapping with the absorption peak of our PR <cit.> and illuminated the sample in a trans-illumination geometry. By illuminating an area ≫ 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 diameter circle, covering almost all of the sample chamber. Under these conditions the cell density is conserved, thus simplifying theoretical modeling. We used a thin sheet polariser imaged onto the sample plane to attenuate the intensity on half of the sample. A digital mirror device <cit.> projected the periodic pattern onto a ≈2.9 diameter area of the sample. 




 §.§ Differential Dynamic Microscopy


DDM measures (v̅,  β) averaged over 10^4 - 10^5 cells under our conditions <cit.>. From ≈30 of wide-field, low-magnification movies, one extracts the power spectrum of the difference between pairs of images delayed by time τ, g(q⃗,τ), where q⃗ is the spatial frequency vector. Under suitable conditions and for isotropic motion, the intermediate scattering function f(q, τ), the q^ th mode of the density autocorrelation function, is given by:


    g(q,τ) = A(q)[1 - f(q,τ)] + B(q)  .
 

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 4 key motility parameters: the mean 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 μm^-1 to give, e.g. v̅=<v(q)>_q and β=<β(q)>_q.

In a dilute system whose structure factor S(q) ≈ 1, A(q) is proportional to the sample density <cit.> and can therefore be used to determine relative densities by ρ_1/ρ_0 = < A_1(q)/A_0(q) >_q  <cit.>. Note that ratioing the A(q)s removes their strong q-dependence. 




Spatially-resolved DDM 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 (v̅, β, ρ, ρ^ 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 ≥ 64. We therefore chose p = 64, corresponding to 90 in the sample (see SI <ref> 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), β_x,y(q), and ρ_x,y(q)/ρ_0=A(q)_x,y/A_0(q), where A_0(q) is measured from the same sample under uniform illumination (i.e. just before switching to a structured light pattern) averaged over 0.5≤ q ≤ 1.5 μm^-1. 
The upper q limit is somewhat lower than what typical for whole-movie analysis <cit.> due to non-systematic failure of fitting at higher q values, presumably due to noise or windowing artefacts <cit.>.





 §.§ Data availability
 The data presented here is available on the Edinburgh DataShare repository <cit.>.









 §.§ Acknowledgements

We 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 Contributions

J.A. and V.A.M. contributed equally to this work. WCKP initiated the work. TP and AD designed mutants constructed by AD. JA and VAM performed experiments, analysed and interpreted data with TP and WCKP, and wrote manuscript with WCKP.



 §.§ Competing Interests
 The authors declare that they have no competing financial interests.



 §.§ Correspondence
 Correspondence and requests for materials should be addressed to J. Arlt (email: j.arlt@ed.ac.uk).

