# Emergent probability fluxes in confined microbial navigation

Edited by David A. Weitz, Harvard University, Cambridge, MA, and approved August 8, 2021 (received for review December 1, 2020)

## Significance

Motile microorganisms commonly live in porous media comprising microhabitats filled with interfaces of complex shape. On such small scales, the interactions with these interfaces, rather than external gradients, dominate their motion in the search for favorable living conditions. We demonstrate with experiments and theory that the geometry of confining interfaces shapes the topology of the most likely, average trajectory, leading to directed fluxes of probability that are not exclusively localized at the near-wall region. Employing this principle allows us to actively shape a microbe’s average direction of movement, which could be of use in the design of topological transport mechanisms for microfluidic environments.

## Abstract

When the motion of a motile cell is observed closely, it appears erratic, and yet the combination of nonequilibrium forces and surfaces can produce striking examples of organization in microbial systems. While most of our current understanding is based on bulk systems or idealized geometries, it remains elusive how and at which length scale self-organization emerges in complex geometries. Here, using experiments and analytical and numerical calculations, we study the motion of motile cells under controlled microfluidic conditions and demonstrate that probability flux loops organize active motion, even at the level of a single cell exploring an isolated compartment of nontrivial geometry. By accounting for the interplay of activity and interfacial forces, we find that the boundary’s curvature determines the nonequilibrium probability fluxes of the motion. We theoretically predict a universal relation between fluxes and global geometric properties that is directly confirmed by experiments. Our findings open the possibility to decipher the most probable trajectories of motile cells and may enable the design of geometries guiding their time-averaged motion.

### Sign up for PNAS alerts.

Get alerts for new articles, or get an alert when an article is cited.

The presence of hidden order and regularities in living systems—seemingly intractable from the point of view of mathematics and physics—was first intuited by Erwin Schrödinger (1). More recently, there has been a renaissance of discoveries of physical principles governing living matter (2, 3). In systems of motile microorganisms, activity and geometry often conspire to create organized collective states (4–15). The robustness of these results invites the question: At what level does such order start emerging? Specifically, is there a lower bound in either the number of participating microbes or the available space for regularities to affect the activity of the cells? While the dynamics of active microswimmers, which propel themselves—e.g., by the periodic beating of one or multiple flagella (16, 17)—are often studied in idealized bulk situations, these microorganisms live in the proximity of interfaces, prosper in wet soil, inhabit porous rocks, and generally encounter complex boundaries regularly (18–20). When colliding with a solid boundary, a microswimmer may interact with it through hydrodynamic interactions mediated by the fluid (21, 22), through steric interactions (23), or through a combination of the two (24, 25). Studies of puller-type microswimmers like

*Chlamydomonas reinhardtii*suggest that, for such microswimmers, steric interactions dominate the dynamics in the presence of such boundaries (23, 26–28).*C. reinhardtii*naturally lives in wet soil (29), an environment dominated by fluid–solid interfaces. A common observation is that such microswimmers are more likely to be found at or in close proximity to the boundary, rather than within the interior (30–33), showing a preference for regions of high wall curvature (28, 30, 34).Here, we report experiments on a single

*C. reinhardtii*cell confined in a compartment with varying boundary curvature, analytical calculations, and simulations that model the microswimmer as an asymmetric dumbbell undergoing active Brownian motion and interacting sterically with the compartment wall. We find that when the compartment’s boundary exhibits nonconstant curvature, the accessible space is partitioned by loops of probability flux that direct and organize the cell’s motion. This finding becomes evident when the nonequilibrium fluxes are extracted from the analysis of experimental and simulation data and also estimated through the Fokker–Planck description of our system. The cell’s front–back asymmetry results in a torque, which, in turn, causes an active reorientation when interacting with the wall. Such shape asymmetry proves crucial for the quantitative agreement of experiments and simulations, making steric effects dominant over hydrodynamics in situations of strong confinement for such flagellated microorganisms. We show that the loops of probability flux are directly linked to the boundary’s gradient of curvature and can be quantified by a dimensionless number $\chi $ accounting for the relevant ratio of length scales.We conduct experiments on single biflagellated

*C. reinhardtii*cells, wild-type strain SAG 11-32b, confined within a quasi-two-dimensional (quasi-2D) microfluidic compartment in the absence of any inlets or outlets. We employ optical bright-field microscopy and particle-tracking techniques to extract cell trajectories over extended time periods (28). Experiments are performed in elliptical compartments with eccentricities varying in the range 0.5 to 0.9 and accessible areas of $\left(7-200\right)\times 1{0}^{3}$ $\mathrm{\mu}{\mathrm{m}}^{2}$. The height of all compartments is ∼22 $\mathrm{\mu}$m, i.e., slightly larger than one cell diameter, comprising cell body and flagella, such that the cell’s motion is confined in two dimensions (*Materials and Methods*).We statistically investigate the average motion by tracking the cell’s position over long trajectories. Fig. 1

*A*shows the overlaid positions of a single*C. reinhardtii*cell confined within a sealed elliptical compartment. While a certain preference to travel alongside the boundary is noticeable, frequent excursions across the microfluidic microhabitat are evident. Periods of time when the cell travels along the boundary are interrupted by rapid reorientations, crossing of the compartments, or curved arcs. The short-time motion of the cell is influenced by the stochasticity associated to the biological motors powering the flagella and their coordination (16, 21, 35) and the small-scale hydrodynamic fluctuations of the fluid stirred by the swimming cell.Fig. 1.

To investigate the nonequilibrium dynamics of the cell’s motion, we compute the trajectory transition rates as follows. Starting from the cell’s trajectories obtained via particle tracking, we divide the space in the compartment with a square grid and determine the probabilities to traverse a certain box in either direction by averaging over all passages of the trajectory on that box. By counting the directed transitions between adjacent boxes $\alpha $ and $\beta $, we can extract the net transition rates ${w}_{\alpha ,\beta}$ from which the $x$ and $y$ components of the flux $\mathit{j}$ are computed (

*SI Appendix*, section I) (36, 37). This analysis reveals loops in the experimental probability fluxes, underlying the run-and-tumble-like motion (21) of the cell within the compartment (Fig. 1*B*). Because of the symmetry of the system and the absence of any preferential direction in the cell’s motion, the loops are symmetrically placed within the ellipse. The fluxes form a strong equatorial current pointing directly toward the ellipse’s apices (the two points on the boundary with the highest curvature) and are then redirected along the boundary away from the apices. These loops are a manifestation of the inherent nonequilibrium nature of the cell’s active motion, which explicitly breaks the invariance of the dynamics between a transition from one state to another and the inverse transition.Numerical simulations of our system lend insight into the properties of these flux loops. We employ a minimal mechanistic model of

*C. reinhardtii*cells as asymmetric dumbbells (*SI Appendix*, Fig. S1), representing the characteristic fore–aft asymmetry of the cell’s body and of a larger space spanned by the stroke-averaged flagella (30, 38, 39). Note that other microswimmer geometries produce fluxes as well, though without matching experiments quantitatively (*SI Appendix*, section V). The translational dynamics of the dumbbell are governed by an overdamped Langevin equation: $\text{d}\mathit{r}/\text{d}t={v}_{0}\mathit{e}+{\mu}_{w}{\mathit{F}}_{w}+\mathit{\eta}$, where $\mathit{r}$ is the position of the geometric center of the dumbbell, ${v}_{0}$ is the self-propulsion velocity, $\mathit{e}$ is the orientation of the dumbbell, ${\mu}_{w}$ is the mobility, and ${\mathit{F}}_{w}$ is the force stemming from steric wall interactions. Furthermore, $\mathit{\eta}$ is a Gaussian white noise to account for translational diffusion with coefficient ${D}_{T}={k}_{B}T{\mu}_{w}$. The orientational dynamics are governed by the equation: $\text{d}\mathit{e}/\text{d}t=\left({\mathit{T}}_{w}/{\tau}_{w}+\mathit{\xi}\right)\times \mathit{e}$, where ${\mathit{T}}_{w}$ is the torque acting at the wall, ${\tau}_{w}$ is the rotational drag coefficient, and $\mathit{\xi}$ is a Gaussian white noise accounting for rotational diffusion with coefficient ${D}_{R}$. All parameters entering the equations of motion were either directly measured in our experiments or extracted from the literature (*Materials and Methods*).Fig. 1

*B*,*Lower*also shows the results of our simulations for the probability flux for a single cell moving according to the above Langevin equations within an elliptical compartment with the same area and eccentricity as in the experiments. We find the analogous structure of flux loops emerging from the numerical results. Strong directional fluxes point from the central region of the compartment toward the apices and then move away along the boundary, closing the flux loops. For any point on the compartment boundary, the fluxes $\mathit{j}$ can be decomposed in a component tangential and one normal to the boundary. We choose the positive normal direction to point into the compartment and the positive tangential direction to always point down the gradient of curvature. The flux analysis reveals that, in the bulk, the cell is most likely to be on a path leading toward a region of high curvature $\kappa $, as indicated by the negative normal fluxes in Fig. 1*C*; after colliding with the wall, a smaller reorientation is needed to escape again at regions of low curvature, whereas escaping at high curvature requires larger turns, which are statistically less likely. This leads to a net statistical flux along the wall toward regions of lower curvature, resulting in a rise of tangential flux with decreasing $\kappa $. The strong, outward equatorial fluxes (Fig. 1*B*) meet the boundary and turn up or down along regions of lower curvature. At those turning points (close to the ellipse apices), the normal fluxes change sign, while the tangential flux become strongest. The quantitative agreement between experiments and simulations in Fig. 1*C*confirms the fundamental role of the force and torque applied by the boundary on the cell.For a more general understanding of the origin of the flux loops, we now turn to a continuum-mechanics approach of the probability flow by computing the Fokker–Planck equation for the system, which readswhere the gorque ${\mathit{G}}_{w}$ is a more convenient way to treat the torque ${\mathit{T}}_{w}=\mathit{e}\times {\mathit{G}}_{w}$. A moment expansion of the probability distribution function $p=p\left(\mathit{r},\mathit{e},t\right)$ in terms of density $\rho \left(\mathit{r}\right)$ and polarization $\mathit{P}\left(\mathit{r}\right)$ allows us to identify the nonequilibrium flux in position space ${\mathit{j}}_{r}={v}_{0}\mathit{P}+{\mu}_{w}{\mathit{F}}_{w}\rho -{D}_{T}\nabla \rho $. The probability flux ${\mathit{j}}_{r}$ obeys a solenoidal condition $\nabla \cdot {\mathit{j}}_{r}=0$ (where the vorticity $\omega =\nabla \times {\mathit{j}}_{r}$ is a function of the forces ${\mathit{F}}_{w}$ and gorques ${\mathit{G}}_{w}$ exerted by the boundary on the swimming cell. Additionally, because of the anisotropic shape of the microswimmer, the vorticity crucially couples with the gradient of curvature $\nabla \kappa $.

$$\begin{array}{ll}\hfill \frac{\partial p}{\partial t}=& -\nabla \cdot \left({v}_{0}\mathit{e}+{\mu}_{w}{\mathit{F}}_{w}-{D}_{T}\nabla \right)p\hfill \\ \hfill & -\mathit{e}\times \frac{\partial}{\partial \mathit{e}}\cdot \left(\frac{1}{{\tau}_{w}}\mathit{e}\times {\mathit{G}}_{w}-{D}_{R}\mathit{e}\times \frac{\partial}{\partial \mathit{e}}\right)p,\hfill \end{array}$$

[1]

*SI Appendix*, section II), which is obeyed in both the experimental and simulated fluxes (Fig. 1*B*). We can solve the dynamics of the probability loops by introducing a stream function ${\mathit{j}}_{r}=\nabla \times \mathit{\psi}$, which satisfies the governing Poisson equation$$\mathrm{\Delta}\psi =-\omega \left({\mathit{F}}_{w},{\mathit{G}}_{w},\nabla \kappa \right),$$

[2]

The flux equation [

**2]**with elliptical boundary conditions can be solved analytically (*SI Appendix*, section III). Fig. 2 shows the resulting fluxes. The qualitative features of the nonequilibrium fluxes found in experiments and simulations are reproduced by the solution of Eq.**2**. Four symmetrically placed flux loops emerge in the region close to the apices. This fact points to a general connection between fluxes and the boundary’s curvature. We elucidate this relationship in the following.Fig. 2.

We can deduct a quantitative relation of the fluxes from the above arguments. From our general arguments (

*SI Appendix*, section II), the steady-state expression for the polarization in close proximity of the wall ${\mathit{P}}_{\mathrm{w}\mathrm{a}\mathrm{l}\mathrm{l}}\approx -{\left(2\nabla \cdot {\mathit{F}}_{w}\right)}^{-1}\left({v}_{0}\nabla +{{\tau}_{w}}^{-1}{\mathit{G}}_{w}\right)\rho $. Recalling the definition of ${\mathit{j}}_{r}$, and the fact that close to the boundary the probability density is proportional to the curvature $\rho \propto \kappa $ (28)—i.e., the cell is more likely to spend more time in regions of high local curvature—we find that, generally, the flux depends on the curvature $\kappa $ and its gradient $\nabla \kappa $. In a circular domain, where the curvature is constant, symmetry prevents the emergence of fluxes (*SI Appendix*, section IV), even though the system is out of equilibrium. A local dependence of curvature breaks the symmetry and allows for nonzero fluxes. Thus, we expect the fluxes to depend on $\nabla \kappa $, but not directly on $\kappa $, as observed in our experiments and simulations.We now generalize our arguments to more complex shapes by inferring the following effective rules: 1) Because flux loops are generated by curvature gradients, the number of flux loops equals the number of zero crossings of curvature gradient, $\nabla \kappa =0$; 2) the magnitude of a flux loop is proportional to the integrated change in curvature $\int \left|\frac{\partial \kappa}{\partial l}\right|\mathrm{d}l$ along a portion of the boundary with arc length $l$; and 3) the number of stagnation points in the flux ($\left|{\mathit{j}}_{r}\right|\approx 0$) is at most one for every two flux loops. These predictions are in line with the topological structure of the fluxes for simulations of our model in compartments of growing complexity with multiple lobes and points of both positive and negative curvature, as shown in Fig. 3.

Fig. 3.

Although Eq. which gives an effective measure for the strength of the nonequilibrium fluxes within the region $S$ (we choose a strip of $20\hspace{0.17em}\mathrm{\mu}$m along the boundary where the fluxes are clearly distinguishable from statistical noise) and where self-propulsion velocity ${v}_{0}$ is used to nondimensionalize $\mathcal{F}$. Through $\mathit{j}$, $\mathcal{F}$ inherits the dependence on $\nabla \kappa \sim \partial \kappa /\partial l$. To capture the global characteristics of the geometry, we are naturally led to define the dimensionless numbercomparing the global change in curvature over the total perimeter $\int \mathrm{d}l$ of the boundary with the typical area of the cell ${A}_{c}$ calculated by squaring the swimmer’s length (

**2**reveals the importance of forces acting at the boundaries, the nonequilibrium nature of the probability distribution imposes a nonlocal spatial distribution of the fluxes. It is then natural to consider the integral of the flux over an area of the compartment$$\mathcal{F}\equiv \frac{1}{{v}_{0}}{\int}_{S}\left|\mathit{j}\right|\text{d}S,$$

[3]

$$\chi \equiv {A}_{c}\frac{\int \left|\frac{\partial \kappa}{\partial l}\right|\mathrm{d}l}{\int \mathrm{d}l},$$

[4]

*Materials and Methods*). We find that the integral flux strength $\mathcal{F}$ depends uniquely on $\chi $ for both experiments and simulations of elliptical compartments with various areas $A$ and eccentricities $\u03f5$ (Fig. 4). The cross-over in the integrated fluxes at $\chi \approx 1{0}^{-3}$ corresponds to the point at which the fluxes can be effectively distinguished from noise, given the available statistics. In fact, to quantify the effect of statistical noise, we consider simulations in circular compartments—where fluxes are strictly absent—and use the same amount of statistics as elliptical compartments; we find an average value ${\mathcal{F}}_{\text{noise}}=5\times 1{0}^{-3}$, marked with a horizontal gray line; the shaded area marks the standard deviation (SD) derived from simulations of circular compartments with the same areas as the ellipses.Fig. 4.

For values of scaled curvature gradients $\chi <1{0}^{-3}$, the fluxes are rather weak, whereas for $\chi \ge 1{0}^{-3}$, strong fluxes emerge and exhibit closed loops, such that their impact on the global dynamics within the compartment is much stronger. The near collapses of both experiments and simulations for all shapes confirm that the boundary’s geometry determines the strength and shape of flux loops within the interior of the compartment.

Taken together, our results show that the boundary of the confining domain imposes a robust topology of loops of probability fluxes at the level of a single active cell. Our experimental and theoretical results demonstrate the intimate connection between the geometric properties of the boundary and the interior of a compartment confining a motile cell. The ensuing probability fluxes impose an organizing structure to the whole compartment’s interior, which statistically guides the cell’s motion. Our study shows that

*C. reinhardtii*cells are very efficient at exploring the available space and that simple geometric features can leave an imprint on the cells’ overall motion in their microhabitat.Harnessing the motion of microorganisms is a promising direction of technological development in active matter (25, 40). Improved efficiency of micromachines will require a better understanding of how microbial cells navigate complex environments and interact with their boundaries. Inducing a statistical bias to directional motion, as a consequence of the nonequilibrium nature of the motion and confining boundaries, might, in fact, help in producing efficient microdevices, even at the scale of a single cell.

## Materials and Methods

### Cell Cultivation.

All experiments were performed by using the wild-type

*Chlamydomonas reinhardtii*strain SAG 11-32b (provided by the Göttingen Algae Culture Collection; SAG). The cells were cultivated axenically in Tris–acetate–phosphate medium in a Memmert IPP110Plus incubator on a 12-h-day–12-h-night cycle. The temperature of the cell cultures in the incubator was kept at 24 $\xb0$C during the day (white light; light intensity of about $3\times 1{0}^{19}$ photons/(m^{2}s) in the center of the incubator) and reduced to a temperature of 22 $\xb0$C during the night (no illumination), respectively. Prior to every experiment, about 50 mL of cell culture was centrifuged (Eppendorf 5804R) for about 10 min at an acceleration of 100 $g$ at ambient temperature. About 40 to 45 mL of the suspension was then removed, and the remaining 5 to 10 mL was allowed to relax for 90 to 120 min in the incubator to ensure that the cells had a sufficient amount of time to eventually regrow their flagella. Finally, this suspension was diluted with cell-culture medium to very low cell concentrations in order to enhance the probability of capturing precisely one single cell per compartment.### Microfluidics.

The experimental chambers were composed of 2D arrays of stand-alone elliptical microfluidic compartments exhibiting a height of about 22 $\mathrm{\mu}$m. These compartment arrays were manufactured by employing soft-lithography techniques using a curable elastomer (polydimethylsiloxane [PDMS]; Sylgard 184 Elastomer Kit, Dow Inc.) and master structures, which we produced by means of photolithography techniques in a clean room. Before filling the chambers with the cell culture, the PDMS-based microfluidic device and a glass microscope slide were treated with air plasma (Pico; Diener Electronic) for 30 s to render their surfaces hydrophilic. A droplet of the diluted cell suspension was placed onto the feature side of the microfluidic device, such that the compartment array was completely filled, and, subsequently, the glass slide was placed atop and gently pressed to tightly seal the experimental compartment. The confining surfaces at the top and bottom prevented the organism from turning out of the plane, enabling a 2D description.

### Microscopy.

For time-resolved cell imaging, we employed bright-field microscopy (Olympus IX-81 inverted microscope) in controlled light conditions (closed box). During all experiments, the microfluidic compartment was illuminated by using a narrow red-light bandpass filter (671 nm, full-width at half-maximum = 10 nm) in order to exclude any photoactive response of the cell, including phototaxis (41, 42) and light-induced adhesion to surfaces (43, 44). A Canon 600D camera (at 25 frames per second; resolution: 1,920 $\times $ 1,080 pixels) was used to record videos of single cells swimming in isolated elliptical compartments for about 5 to 30 min each. In order to increase the statistics substantially, these single-cell recordings were independently repeated three to eight times for each chamber geometry.

### Image Processing and Cell Tracking

All videos were converted into sequences of 8-bit grayscale images with improved contrast using a custom-made MATLAB algorithm. The compartment boundaries were manually identified to restrict the cell tracking to the region available to the motion of the cell. Finally, 2D cell tracking was performed by using MATLAB based on the protocol developed by Crocker and Grier (45).

### Numerical Model and Simulation Parameters

The

*C. reinhardtii*cells are modeled as asymmetric dumbbells (*SI Appendix*, Fig. S1) with a large sphere in front and a smaller sphere in the back, representing the fore–aft asymmetry of body and appendages (38, 39). The equation of motion for the position $\mathit{r}$ of the active dumbbell is given by$$\frac{\text{d}\mathit{r}}{\text{d}t}={v}_{0}\mathit{e}+{\mu}_{w}{\mathit{F}}_{w}+\mathit{\eta}.$$

[5]

Here, ${v}_{0}=100\hspace{0.17em}\mathrm{\mu}\text{m/s}$ is the self-propulsion speed of the cell, and $\mathit{\eta}$ is a Gaussian white noise with correlator $\u27e8\mathit{\eta}\left(t\right)\mathit{\eta}\left({t}^{\prime}\right)\u27e9=2{k}_{\mathrm{B}}T{\mu}_{w}\mathbf{1}\delta \left(t-{t}^{\prime}\right)$ and translational diffusion coefficient ${k}_{\mathrm{B}}T{\mu}_{w}=250\hspace{0.17em}\mathrm{\mu}$m

^{2}/s (both numerical values are taken from ref. 28). The term ${\mathit{F}}_{w}$ accounts for steric wall interactions of the dumbbell and was computed by ${\mathit{F}}_{w}={\mathit{F}}_{1}+{\mathit{F}}_{2}$ with ${\mathit{F}}_{\alpha}=-\nabla {U}_{\alpha}\left(r\right)$, $\alpha =\mathrm{1,2}$, where 1 and 2 refer to the large and small sphere of the dumbbell, respectively. To compute the respective steric forces, we used the Weeks–Chandler–Anderson potential (46)$${U}_{\alpha}\left(d\right)/\left({k}_{\mathrm{B}}T\right)=4{U}_{0}\left[{\left(\frac{{a}_{\alpha}}{d}\right)}^{12}-{\left(\frac{{a}_{\alpha}}{d}\right)}^{6}\right]+{U}_{0}\hspace{0.17em},$$

[6]

if $d<{2}^{1/6}{a}_{\alpha}$, and zero otherwise, where $d$ is the distance of the sphere $\alpha \in \left\{\mathrm{1,2}\right\}$ to the wall of the compartment. The radii of the large (1) and small (2) circles of the dumbbell are ${a}_{1}=5\hspace{0.17em}\mathrm{\mu}$m and ${a}_{2}=2.5\hspace{0.17em}\mathrm{\mu}$m, and we use $\epsilon =10$ to obtain a sufficiently strong screening (the values are taken from ref. 28).

The orientation $\mathit{e}$ is defined as the unit vector pointing from the small to the large circle of the dumbbell. The equation of motion for the orientation is given by

$$\frac{\text{d}\mathit{e}}{\text{d}t}=\left({\mathit{T}}_{w}/{\tau}_{w}+\mathit{\xi}\right)\times \mathit{e}.$$

[7]

Here, $\mathit{\xi}$ is a Gaussian white noise with correlator $\u27e8\mathit{\xi}\left(t\right)\mathit{\xi}\left({t}^{\prime}\right)\u27e9=\frac{2{k}_{\mathrm{B}}T}{{\tau}_{p}}\mathbf{1}\delta \left(t-{t}^{\prime}\right)$ and rotational diffusion coefficient ${\tau}_{p}/{k}_{\mathrm{B}}T=2\hspace{0.17em}\text{s}$ (23). The torque acting at the wall is computed by ${\mathit{T}}_{w}={\mathit{T}}_{1}+{\mathit{T}}_{2}$, where we use ${\mathit{T}}_{1}=\left({\mathit{r}}_{1}-\mathit{r}\right)\times {\mathit{F}}_{1}=l\left(\mathit{e}\times {\mathit{F}}_{1}\right)/2$, ${\mathit{T}}_{2}=-l\left(\mathit{e}\times {\mathit{F}}_{2}\right)/2$, and $l=5\hspace{0.17em}\mathrm{\mu}$m. For the shear time at the wall, we use ${\tau}_{w}/{k}_{\mathrm{B}}T=0.15\hspace{0.17em}\text{s}$ (23). In addition to Eq.

**7**, we include a run-and-tumble motion of the cell. Each tumbling event is instantaneous, and the time between each tumbling event is sampled for an exponential distribution with mean ${\tau}_{p}/{k}_{\mathrm{B}}T$. The relative tumbling angle ${\varphi}_{\mathrm{t}\mathrm{u}\mathrm{m}\mathrm{b}\mathrm{l}\mathrm{e}}$ is drawn from a Gaussian distribution with an SD of 0.1 and a mean of $\pi /2$.## Data Availability

Motile cell trajectory data, analysis and simulation code have been deposited at zenodo.org (DOI: 10.5281/zenodo.5497768) (47).

## Acknowledgments

We thank the Göttingen Algae Culture Collection for providing the

*C. reinhardtii*wild-type strain SAG 11-32b. We thank Stephan Herminghaus for helpful discussions.## Supporting Information

Appendix (PDF)

- Download
- 3.61 MB

Movie S1.

Video of a single

*Chlamydomonas reinhardtii*cell (strain SAG 11-32b) swimming within a stand-alone elliptical microfluidic compartment with a major semi-axis of 157 μm and a minor semi-axis of 63 μm.- Download
- 14.90 MB

## References

1

E. Schrödinger,

*What Is Life? The Physical Aspect of the Living Cell*(Cambridge University Press, Cambridge, UK, 1944).2

T. Vicsek, A. Zafeiris, Collective motion.

*Phys. Rep.***517**, 71–140 (2012).3

M. C. Marchetti et al., Hydrodynamics of soft active matter.

*Rev. Mod. Phys.***85**, 1143 (2013).4

I. H. Riedel, K. Kruse, J. Howard, A self-organized vortex array of hydrodynamically entrained sperm cells.

*Science***309**, 300–303 (2005).5

A. Sokolov, I. S. Aranson, Physical properties of collective motion in suspensions of bacteria.

*Phys. Rev. Lett.***109**, 248109 (2012).6

F. G. Woodhouse, R. E. Goldstein, Spontaneous circulation of confined active suspensions.

*Phys. Rev. Lett.***109**, 168105 (2012).7

H. Wioland, F. G. Woodhouse, J. Dunkel, J. O. Kessler, R. E. Goldstein, Confinement stabilizes a bacterial suspension into a spiral vortex.

*Phys. Rev. Lett.***110**, 268102 (2013).8

R. Großmann, P. Romanczuk, M. Bär, L. Schimansky-Geier, Pattern formation in active particle systems due to competing alignment interactions.

*Eur. Phys. J. Spec. Top.***224**, 1325 (2015).9

E. Lushi, H. Wioland, R. E. Goldstein, Fluid flows created by swimming bacteria drive self-organization in confined suspensions.

*Proc. Natl. Acad. Sci. U.S.A.***111**, 9733–9738 (2014).10

H. Wioland, F. G. Woodhouse, J. Dunkel, R. E. Goldstein, Ferromagnetic and antiferromagnetic order in bacterial vortex lattices.

*Nat. Phys.***12**, 341–345 (2016).11

H. Wioland, E. Lushi, R. E. Goldstein, Directed collective motion of bacteria under channel confinement.

*New J. Phys.***18**, 075002 (2016).12

K. Beppu et al., Geometry-driven collective ordering of bacterial vortices.

*Soft Matter***13**, 5038–5043 (2017).13

M. Theillard, R. Alonso-Matilla, D. Saintillan, Geometric control of active collective motion.

*Soft Matter***13**, 363–375 (2017).14

G. Frangipane et al., Invariance properties of bacterial random walks in complex structures.

*Nat. Commun.***10**, 2442 (2019).15

H. H. Wensink, V. Kantsler, R. E. Goldstein, J. Dunkel, Controlling active self-assembly through broken particle-shape symmetry.

*Phys. Rev. E Stat. Nonlin. Soft Matter Phys.***89**, 010302 (2014).16

K. Y. Wan, R. E. Goldstein, Coordinated beating of algal flagella is mediated by basal coupling.

*Proc. Natl. Acad. Sci. U.S.A.***113**, E2784–E2793 (2016).17

T. J. Böddeker, S. Karpitschka, C. T. Kreis, Q. Magdelaine, O. Bäumchen, Dynamic force measurements on swimming

*Chlamydomonas*cells using micropipette force sensors.*J. R. Soc. Interface***17**, 20190580 (2020).18

W. C. Ghiorse, J. T. Wilson, “Microbial ecology of the terrestrial subsurface,” in

*Advances in Applied Microbiology*, A. I. Laskin, ed. (Elsevier, Amsterdam, Netherlands, 1988), vol. 33, pp. 107–172.19

M. G. van der Heijden, R. D. Bardgett, N. M. van Straalen, The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems.

*Ecol. Lett.***11**, 296–310 (2008).20

G. M. Gadd, Metals, minerals and microbes: Geomicrobiology and bioremediation.

*Microbiology (Reading)***156**, 609–643 (2010).21

M. Polin, I. Tuval, K. Drescher, J. P. Gollub, R. E. Goldstein,

*Chlamydomonas*swims with two “gears” in a eukaryotic version of run-and-tumble locomotion.*Science***325**, 487–490 (2009).22

M. Contino, E. Lushi, I. Tuval, V. Kantsler, M. Polin, Microalgae scatter off solid surfaces by hydrodynamic and contact forces.

*Phys. Rev. Lett.***115**, 258102 (2015).23

V. Kantsler, J. Dunkel, M. Polin, R. E. Goldstein, Ciliary contact interactions dominate surface scattering of swimming eukaryotes.

*Proc. Natl. Acad. Sci. U.S.A.***110**, 1187–1192 (2013).24

J. Elgeti, G. Gompper, Microswimmers near surfaces.

*Eur. Phys. J. Spec. Top.***225**, 2333 (2016).25

S. E. Hulme et al., Using ratchets and sorters to fractionate motile cells of

*Escherichia coli*by length.*Lab Chip***8**, 1888–1895 (2008).26

R. Ledesma-Aguilar, J. M. Yeomans, Enhanced motility of a microswimmer in rigid and elastic confinement.

*Phys. Rev. Lett.***111**, 138101 (2013).27

L. Caprini, U. M. B. Marconi, Active particles under confinement and effective force generation among surfaces.

*Soft Matter***14**, 9044–9054 (2018).28

T. Ostapenko et al., Curvature-guided motility of microalgae in geometric confinement.

*Phys. Rev. Lett.***120**, 068002 (2018).29

E. H. Harris,

*The Chlamydomonas Sourcebook: Introduction to Chlamydomonas and its Laboratory Use*(Academic Press, New York, 2009),**vol. 1**.30

A. Wysocki, J. Elgeti, G. Gompper, Giant adsorption of microswimmers: Duality of shape asymmetry and wall curvature.

*Phys. Rev. E Stat. Nonlin. Soft Matter Phys.***91**, 050302 (2015).31

Y. Ibrahim, T. B. Liverpool, How walls affect the dynamics of self-phoretic microswimmers.

*Eur. Phys. J. Spec. Top.***225**, 1843 (2016).32

J. Elgeti, G. Gompper, Wall accumulation of self-propelled spheres.

*EPL***101**, 48003 (2013).33

J. Elgeti, G. Gompper, Run-and-tumble dynamics of self-propelled particles in confinement.

*EPL***109**, 58003 (2015).34

Y. Fily, A. Baskaran, M. F. Hagan, Dynamics of self-propelled particles under strong confinement.

*Soft Matter***10**, 5609–5617 (2014).35

V. F. Geyer, F. Jülicher, J. Howard, B. M. Friedrich, Cell-body rocking is a dominant mechanism for flagellar synchronization in a swimming alga.

*Proc. Natl. Acad. Sci. U.S.A.***110**, 18058–18063 (2013).36

R. K. P. Zia, B. Schmittmann, Probability currents as principal characteristics in the statistical mechanics of non-equilibrium steady states.

*J. Stat. Mech.***2007**, P07012 (2007).37

C. Battle et al., Broken detailed balance at mesoscopic scales in active biological systems.

*Science***352**, 604–607 (2016).38

A. M. Roberts, F. M. Deacon, Gravitaxis in motile micro-organisms: The role of fore–aft body asymmetry.

*J. Fluid Mech.***452**, 405–423 (2002).39

A. M. Roberts, Mechanisms of gravitaxis in

*Chlamydomonas*.*Biol. Bull.***210**, 78–80 (2006).40

R. Di Leonardo et al., Bacterial ratchet motors.

*Proc. Natl. Acad. Sci. U.S.A.***107**, 9541–9545 (2010).41

P. Berthold et al., Channelrhodopsin-1 initiates phototaxis and photophobic responses in Chlamydomonas by immediate light-induced depolarization.

*Plant Cell***20**, 1665–1677 (2008).42

K. W. Foster et al., A rhodopsin is the functional photoreceptor for phototaxis in the unicellular eukaryote

*Chlamydomonas*.*Nature***311**, 756–759 (1984).43

C. T. Kreis, M. Le Blay, C. Linne, M. M. Makowski, O. Bäumchen, Adhesion of

*Chlamydomonas*microalgae to surfaces is switchable by light.*Nat. Phys.***14**, 45–49 (2018).44

C. T. Kreis, A. Grangier, O. Bäumchen, In vivo adhesion force measurements of

*Chlamydomonas*on model substrates.*Soft Matter***15**, 3027–3035 (2019).45

J. C. Crocker, D. G. Grier, Methods of digital video microscopy for colloidal studies.

*J. Coll. Interf. Sci.***179**, 298–310 (1996).46

J. D. Weeks, D. Chandler, H. C. Andersen, Role of repulsive forces in determining the equilibrium structure of simple liquids.

*J. Chem. Phys.***54**, 5237 (2003).47

J. Cammann, et al., Emergent probability fluxes in confined microbial navigation: Data and code, Zenodo. https://zenodo.org/record/5497768. Deposited 9 September 2021.

## Information & Authors

### Information

#### Published in

#### Classifications

#### Copyright

Copyright © 2021 the Author(s). Published by PNAS. This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).

#### Data Availability

Motile cell trajectory data, analysis and simulation code have been deposited at zenodo.org (DOI: 10.5281/zenodo.5497768) (47).

#### Submission history

**Accepted**: August 8, 2021

**Published online**: September 23, 2021

**Published in issue**: September 28, 2021

#### Keywords

#### Acknowledgments

We thank the Göttingen Algae Culture Collection for providing the

*C. reinhardtii*wild-type strain SAG 11-32b. We thank Stephan Herminghaus for helpful discussions.#### Notes

This article is a PNAS Direct Submission.

### Authors

#### Competing Interests

The authors declare no competing interest.

## Metrics & Citations

### Metrics

#### Citation statements

#### Altmetrics

### Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

#### Cited by

Loading...

## View Options

### View options

#### PDF format

Download this article as a PDF file

DOWNLOAD PDF### Get Access

#### Login options

Check if you have access through your login credentials or your institution to get full access on this article.

Personal login Institutional Login#### Recommend to a librarian

Recommend PNAS to a Librarian#### Purchase options

Purchase this article to get full access to it.