# Confinement-induced accumulation and de-mixing of microscopic active-passive mixtures

Stephen Williams, Raphaël Jeanneret, Idan Tuval, Marco Polin

► To cite this version:

Stephen Williams, Raphaël Jeanneret, Idan Tuval, Marco Polin. Confinement-induced accumulation and de-mixing of microscopic active-passive mixtures. Nature Communications, 2022, 13 (1), pp.4776. 10.1038/s41467-022-32520-9. hal-03775236

HAL Id: hal-03775236

https://hal.science/hal-03775236v1

Submitted on 9 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Confinement-induced accumulation and de-mixing of microscopic active-passive mixtures

Received: 9 May 2022

Accepted: 3 August 2022

Published online: 15 August 2022

<figure><img src="image_3.png" /></figure>

Check for updates

Stephen Williams<sup>1</sup>, Raphaël Jeanneret<sup>2</sup>, Idan Tuval<sup>3,4</sup> & Marco Polin<sup>1,3,4</sup> ✉

Understanding the out-of-equilibrium properties of noisy microscale systems and the extent to which they can be modulated externally, is a crucial scientific and technological challenge. It holds the promise to unlock disruptive new technologies ranging from targeted delivery of chemicals within the body to directed assembly of new materials. Here we focus on how active matter can be harnessed to transport passive microscopic systems in a statistically predictable way. Using a minimal active-passive system of weakly Brownian particles and swimming microalgae, we show that spatial confinement leads to a complex non-monotonic steady-state distribution of colloids, with a pronounced peak at the boundary. The particles’ emergent active dynamics is well captured by a space-dependent Poisson process resulting from the space-dependent motion of the algae. Based on our findings, we then realise experimentally the de-mixing of the active-passive suspension, opening the way for manipulating colloidal objects via controlled activity fields.

Evolution has enabled living systems to achieve an exquisite control of matter at the microscopic level. From the precise positioning of chromosomes along the mitotic spindle<fnref n="1" /> to the many types of embryonic gastrulation² cells harness their internal and external motility to reach a predictable order despite the stochasticity intrinsic to the microscopic realm³. Understanding how order emerges in these active systems is a fundamental scientific challenge with the potential to bring disruptive technologies for the macroscopic control of microscopic structures. Here we address this problem within a minimal active-passive experimental model system.

From a physical perspective, living systems fall under the general category of active matter⁴, characterised by emergent phenomena including flocking⁵⁻⁷, active turbulence⁸⁻¹¹ and motility-induced phase separation¹²,¹³. As our understanding of single-species active systems progresses, attention has started to veer towards more complex cases where components with different levels of activity interact. Indeed this is often the case for biological active matter. Intracellular activity, for example, can be used for spatial organisation of passive intracellular organelles¹⁴⁻¹⁶; and clustering induced by motility differentials helps

bacterial swarms expand despite antibiotic exposure¹⁷. In order to study the emergent properties of these complex systems, an important and phenomenologically rich minimal model is one that mixes active and inert agents. Active impurities have been used to alter the dynamics of grain boundaries in colloidal crystals¹⁸,¹⁹ and favour the formation of metastable clusters in semi-dilute suspensions²⁰. Sufficiently large concentrations of both active and passive species often reveal a rich spectrum of phases depending on the interactions between constituents²¹⁻²⁹.

Active baths, where individual passive inclusions are dispersed within an active suspension, are particularly appealing. Their conceptual simplicity makes them a natural starting point to develop a statistical theory of active transport¹⁰⁻³³ with potential applications to micro-cargo delivery, micro-actuation³⁴⁻³⁶ and nutrient transport³⁷,³⁸. In general, passive particles in homogeneous and isotropic active baths display enhanced diffusion due to a continuous energy transfer from the active component via direct collisions and hydrodynamic interactions³⁹⁻⁴⁷. Designing larger objects with asymmetric shapes can then turn active diffusion into noisy active translation or rotation³⁰,⁴⁸,⁴⁹.

<footnote n="1">¹Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom. ²Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, F-75005 Paris, France. ³Departament de Física, Universitat de les Illes Balears, 07071 Palma de Mallorca, Spain. ⁴Instituto Mediterráneo de Estudios Avanzados, IMEDEA, Miquel Marques 21, 07190 Esporles, Spain. ✉ e-mail: mpolin@imedea.uib-csic.es</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1** | Experimental setup and colloidal behaviour. a Dark-field image of cells and beads within a straight microfluidic channel (here width $2W = 200$ µm). The two species are highlighted with, respectively, green ellipses and red circles. The schematic illustrates the coordinate system used throughout. b Steady-state colloidal probability distribution functions across the channel (edge to midpoint) for different $W$ values. Boundary accumulation is followed by a depleted adjacent region which gives way to a uniform distribution within the channel. c Example of a typical colloidal trajectory inside a 100 µm-width channel. Colours represent the mean frame-to-frame speed. The dynamics is composed of slow diffusive-like motion (blue sections) interspersed with fast straighter jumps (red sections).</figcaption></figure>

However the strategies for the predictable patterning and transport of passive cargo are still very limited. This is in contrast with the exquisite level of external control possible on the motility of biological and synthetic active particles themselves$^{50–55}$, such as photokinetic bacteria whose swim speed depends on the local light intensity; property that can be harnessed to pattern surfaces$^{53,54}$.

Here we show that a steady gradient of activity within active-passive suspensions can be harnessed to control the fate of generic passive particles. Quasy-2D microfluidic channels are filled with a binary suspension of polystyrene colloids and unicellular biflagellate microalgae *Chlamydomonas reinhardtii* (see Supplementary Movie 1). Confinement induces a spatially inhomogeneous and anisotropic distribution of microswimmers as a result of wall scattering$^{56–59}$, which translates into a space-dependent active noise for the colloids. We show that this produces complex non-monotonic colloidal distributions with accumulation at the boundaries and then develop an effective microscopic jump-diffusion model for the colloidal dynamics. The latter shows excellent analytical and numerical agreement with the experimental results. Finally we demonstrate how confinement with aptly designed microfluidic chips can fuel the de-mixing of active-passive suspensions, opening the way for manipulating passive colloidal objects via controlled activity fields.

## Results

**Colloid concentration is non-homogeneous under confinement**

A full description of the experimental procedures (culturing, video microscopy and data analysis) can be found in the Methods section and Supplementary Materials (data and codes also available on Zenodo$^{60}$). Figure 1a shows a detail of the main section of the first type of microfluidic setup used. These devices, between 14 µm and 20 µm thick, are composed of two reservoir chambers connected by long and straight channels of constant width $2W$ ranging from 50 to 200 µm. These are filled with a dilute mixture of *Chlamydomonas reinhardtii* (CR; strain CC-125, radius $R = 4–5$ µm) and weakly Brownian colloids (radius $a = 5 \pm 0.5$ µm) at surface fractions $\phi_{CR} \approx \phi_{col} \approx 2–3\%$ (bulk concentrations $\sim 2–3 \times 10^7$ particles/mL). The design ensures a steady concentration along the connecting channels for both species (see Supplementary Movie 1). Figure 1a shows the coordinate system

employed, which has been symmetrized with respect to the channel axis.

As typical of self-propelled particles, the algae tend to accumulate at the boundaries due to their interactions with the side walls$^{56–58}$. This appears as a significant peak in their time-averaged concentration profile at a position $y_{CR} \approx 15$ µm (Supplementary Fig. 1), roughly equivalent to the sum of the cell radius and the flagellar length. Unexpectedly, we find that also the colloids explore the available space in an inhomogeneous way. Their steady-state distribution (Fig. 1b) shows a clear peak at about one particle radius ($y_{col} = 5.9 \pm 1.5$ µm), followed by a depleted region between 10 µm and 20 µm from the wall, which roughly corresponds to the peak in algal density, before plateauing to a uniform concentration further inside the channels. This effect is also observed within circular chambers (Supplementary Fig. 2) suggesting it is a robust feature of the system. These distributions are in stark contrast with the equilibrium case (i.e., without microswimmers) for which the colloids are expected to be uniformly distributed despite spatial variations in colloidal diffusivity due to hydrodynamics$^{61}$ (Supplementary Fig. 3).

As a first step to gain insight into the colloids' experimental distributions, we characterise their dynamics within our active bath. Figure 1c shows a typical colloidal trajectory (561 s) colour-coded for the average frame-to-frame speed $v$. The dynamics can be understood as a combination of periods of slow diffusive-like displacements ($v \sim 5$ µm/s;) and fast longer jumps ($v \sim 30$ µm/s; see Supplementary Section 1). The latter is reminiscent of hydrodynamic entrainment events reported for micron-sized colloids$^{47,62}$, although for these larger particles the prominent role appears to be played by direct flagellar interactions (see Supplementary Movie 2). Irrespective of their specific origin, jump events dominate the active transport of colloids within the system. Following$^{47,63}$, we use time-correlation of displacements along individual trajectories and an estimate for the expected magnitude of frame-to-frame diffusive displacements to identify active colloidal jumps and extract their statistical properties as a function of starting position $y$ (we assume translational invariance along the channels' axis). Figure 2a shows the space-dependent distributions of waiting times between successive jumps, with colours representing the distance from the boundary. Regardless of the position, all the distributions are exponential above ~3 s but deviate from it at shorter times. This deviation from a simple Poisson process is due to the large colloidal size influencing the motility of nearby algae as reported also in the case of bacteria$^{64}$. Here, this can lead to a rapid succession of interactions with the same cell (see Supplementary Movie 3). Nevertheless, for our purposes, the waiting dynamics can still be approximated with a single effective rate $\lambda(y)$ (Fig. 2a solid lines; for further details see Supplementary Section 2). Figure 2b shows that these rates increase monotonically with increasing distance from the wall. Notice that this curve does not mirror the algal accumulation at the boundary (Supplementary Fig. 1), a reflection of the fact that proximity to the wall curtails the range of possible swimming directions that algae can have when interacting with the colloid.

Next, we look at the modulation in jump orientation and magnitude. Figure 2c shows the distribution of jump directions $P(\theta, y)$, with motion towards or away from the boundary corresponding to negative and positive values of $\theta$ respectively. Approaching the wall, the distribution develops a marked peak at $\theta = 0$ indicating a strongly anisotropic active motion preferentially parallel to the boundary. This feature reflects the anisotropy in algal dynamics that results from the interaction with the wall$^{56}$ and that can be measured up to ~100 µm from the boundary because of the persistence in cells' trajectories$^{59}$. The distribution functions of jump magnitudes are similar to those already reported for micron-sized particles$^{47,62}$, with an exponential decay above ~4 µm (Supplementary Fig. 4a). These can be used to calculate the average jump length $\langle l(y) \rangle$ (Supplementary Fig. 4b) which decreases by ~30% from the bulk level within
<figure><img src="image_5.png" /><figcaption>**Fig. 2 | Colloidal jump dynamics.** a Probability distribution function of waiting times between consecutive jumps (W = 50 µm, semilog plot) Circles: experiments. Solid lines: exponential fits. The distributions are shifted vertically for clarity. Colours code for position across the channel (red: boundary; blue: channel centre). b Characteristic encounter rates λ(y) for the four values of W, as a function of distance from the channel boundaries. For better comparison, the values have been rescaled to the cell concentration of the 100 µm-wide channel (see Supplementary Section 3). Circles: experiments. Solid lines: exponential fits (see Methods). Error bars from the fit uncertainties. c Probability distribution function P(θ, y) for the direction θ of active jumps vs. distance y from the channel boundary (W = 50 µm). d Semilog plot of active jump distributions, q(y), for different distances y from the boundary (Colours as in C). Circles: experiments. Solid lines: exponential fits. The distributions are shifted vertically for clarity.</figcaption></figure>

the first 10–15 µm from the boundary (first 5–10 µm accessible to the beads). A decrease is indeed expected due to the obvious limitations to the active movement of the colloids imposed by a nearby boundary.

As we are interested in the passive particles' dynamics across the channels, it is useful to reduce the full two-dimensional jump distribution functions to their projection $q_y(y)$ along the y axis. Here $y_j$ is the y-coordinate of the active displacement and the subscript indicates the distance from the nearest wall (see Fig. 1a for the frame of reference). Figure 2d shows that these distributions are generally well approximated by the combination of two exponentials, one each for positive and negative directions (respectively away from and towards the boundary). The exception is for the positive tails of the distributions closest to the wall, which decay slower than expected from the exponential fit. They will not be considered in the following. The distributions $q_y(y_j)$ can then be approximated as

$$q_y(y_j) = \frac{1}{L_+(y) + L_-(y)} e^{-|y_j|/L_-(y)}, 0 \le y_j, \quad (1)$$

where $L_+(y)$ are the characteristic lengths of the exponential fits to positive and negative jumps respectively (see Methods for details on how the characteristic lengths were extracted). As shown in Fig. 3-inset (red and blue symbols and solid lines) these are identical in the core of the channel ($L_+ = L_- \simeq 4$ µm) and decrease to the same small value close to the wall (-0.5 µm) as a consequence of the increasing polarisation of the active displacements along the boundary (Fig. 2c). However, the length scale for the transition between boundary and core values is shorter for $L_+(y)$ than for $L_-(y)$ (-12 µm vs. -25 µm, see Fig. 3-inset). This behaviour, common to all values of W explored (Supplementary Fig. 5), gives rise to a net drift $L_+(y)-L_-(y) > 0$ towards the bulk of the system (green symbols and solid line, Fig. 3-inset) which, as discussed below, is responsible for the depleted region observed in the experimental

<figure><img src="image_6.png" /><figcaption>**Fig. 3 | Comparison between experimental results and the jump-diffusion model.** Steady-state colloidal probability distribution functions in a $2W = 100$ µm channel: experiments (circles); full dynamics (black solid line); constant encounter rate $\lambda(y) := \text{bulk value}$ (teal dashed line); no thermal diffusion ($D_0 = 0$) (purple dashed line); homogeneous and isotropic jump size distribution ($q_y(y) := \text{bulk distribution}$) (orange solid line); analytical model (olive solid line). The analytical model curve has been shifted downward for clarity. Inset: characteristic jump size away from the boundary ($L_+(y)$, red circles) and toward the boundary ($L_-(y)$, blue circles) together with the first and second moments $m_1(y)$ and $\sqrt{m_2(y)}$ derived from them (Eq. (5)). Solid lines are fitted heuristic analytical functions used to reconstruct the colloidal distribution (see Methods). Error bars propagated from fit uncertainties.</figcaption></figure>

colloidal distributions. The difference between the $L_+$'s is likely due to a combination of local anisotropy in the distribution of swimming directions of the algae, and the stopping of colloids' active displacements by the wall.

The focus on the active jumps displayed by the passive particles can only be justified if this part of the dynamics is indeed sufficient to capture the experimental steady-state distributions. This was first tested within a 1D numerical simulation of a weakly Brownian particle (diffusivity $D_0$), moving in $y \in [0, 2W]$ and subject to a space-dependent Poisson noise of rate $\lambda(y)$ and value drawn from $q_y(y)$ (see Methods Sec. 4.5 for details on the integration scheme). Figure 3 shows that the resulting steady-state distribution (black solid curve) is in excellent agreement with the experimental one (black circles) ($2W = 100$ µm). Numerical simulations also provide a convenient way to explore which elements of the effective colloidal dynamics play the most prominent role. For example, fixing $\lambda(y)$ to the bulk value everywhere leads to a very minimal change in the spatial distribution (Fig. 3, teal dashed line), showing that the spatial dependence of the jump frequency close to the wall is not a major factor in the present case. Similarly, removing completely the background diffusion by setting $D_0 = 0$ leaves the distribution unchanged (Fig. 3, purple dashed line) consistent with the fact that, in our experimental system, colloidal transport is dominated by the Poissonian jumps. On the other hand, simulations that remove the spatial dependence of the jump distributions $q_y(y)$ and use their isotropic bulk value everywhere, show a boundary accumulation of colloids that is narrower and higher than the experimental curve and has no intermediate depletion (Fig. 3, orange solid line). This confirms that space-dependent jump anisotropy plays a key role in determining the experimental colloidal distribution.

### Minimal continuum model for the dynamics of the colloids

The numerical validation of the jump-diffusion dynamics motivates a simple analytical model for the evolution of $P_t(y)$, the probability
density of finding a colloid at position y at time t:

$$ \frac{\partial P_t(y)}{\partial t} = D_0 \frac{\partial^2 P_t(y)}{\partial y^2} - \lambda(y)P_t(y) + \int_{-\infty}^{+\infty} \lambda(y-y_j)P_t(y-y_j)q_{y-y_j}(y_j)dy_j. \quad (2) $$

In this master equation, local changes in $P_t(y)$ are due either to diffusion (first term in the r.h.s) or to the balance between active jumps from the current position or towards it from elsewhere (second and third terms in the r.h.s., respectively). This continuum equation can also be derived more formally from a stochastic description of the dynamics of single colloids, following Denisov and Bystrik⁶⁵ (see Supplementary Section 4). It is worth noting also that Eq. (2) assumes infinitely fast jumps ("teleportation"). While this might seem like a drastic approximation, it is not expected to impact the resulting spatial distributions as long as the typical jump duration is sufficiently shorter than the average waiting time between jumps. This is indeed the case in the current system (-1.2 s vs. ≥3.5 s respectively; see Fig. 2b and Supplementary Fig. 6).

The non-local term in Eq. (2) means that the model is non-tractable for generic kernels $\lambda(y)q_y(y-y_j)$. In order to make progress, we, therefore, perform a Kramers-Moyal expansion⁶⁶ and truncate the series at second-order, providing an approximation of the system at the drift-diffusion level. This type of approximation has already been used successfully to capture the essential features of run-and-tumble bacteria⁶⁷, active Brownian particles⁶⁸ and active Ornstein-Uhlenbeck particles⁶⁹. This approximation is expected to hold if the characteristic jump size is small compared to the length scale of the heterogeneities in the dynamics. In our case we have $\langle l_{bulk} \rangle \approx 5$ µm (Supplementary Figure 4b) while a length scale for the heterogeneity can be estimated for example from the $L_\pm$ curves (≥12 µm). As detailed in the Supplementary Materials (Supplementary Section 5) the second-order expansion of Eq. (2) gives a drift-diffusion equation which depends on $\lambda(y)$ and the first and second moments $m_{1,2}(y)$ of $q_y(y_j)$:

$$ \begin{aligned} \frac{\partial P_t(y)}{\partial t} &= \frac{\partial}{\partial y} \left[ D_{\text{eff}}(y) \frac{\partial P_t}{\partial y} - V_{\text{eff}}(y) P_t(y) \right] \\ D_{\text{eff}}(y) &= D_0 + \frac{\lambda(y)m_2(y)}{2} \\ V_{\text{eff}}(y) &= \lambda(y)m_1(y) - \frac{1}{2} \frac{\partial}{\partial y} [\lambda(y)m_2(y)] \\ m_n(y) &= \int_{-\infty}^{+\infty} y_j^n q_y(y_j) dy_j. \end{aligned} \quad (3) $$

Notice that $D_{\text{eff}}$ is similar to the asymptotic diffusivity obtained in⁴⁷ for homogeneous and isotropic entrainment-dominated transport, with a contribution from the jump events given by the product of jump frequency and variance of jump size. Enforcing no-flux at the boundaries one obtains a closed form for the steady-state solution:

$$ P(y) = \frac{B}{D_{\text{eff}}(y)} \exp \left( \int^{y} \frac{\lambda(y')m_1(y')}{D_{\text{eff}}(y')} dy' \right), \quad (4) $$

where B is the normalisation constant. Without a net asymmetry in the dynamics (i.e., when $m_1(y) \equiv 0$), the solution is that of a state-dependent diffusive process with Itô convention for multiplicative noise integration. In our case, however, the first moment of $q_y$ takes significant values in the first ~30 µm from the boundaries (green solid line Fig. 3-Inset) and should not be neglected.

Comparison between Eq. (4) and the experiments is facilitated by having analytical expressions for $\lambda$, $m_1$ and $m_2$. Figure 2b shows that $\lambda$ can be described heuristically by a simple exponential relaxation from the wall to the bulk values (black solid line). As for $m_1$ and $m_2$, the

<figure><img src="image_7.png" /><figcaption>Fig. 4 | De-mixing of a binary mixture. a Schematics and details of the microfluidic chamber: 200 µm-radius circular chambers decorated with side channels of width $d=7–8$ µm. Both cells ($2R=8–10$ µm) and polystyrene beads ($2a=6$ µm) are visible in the picture. b, c Time evolution of the number of particles and surface number density within the circular chamber ($N_c, \rho_c$ red solid line) and the side channels ($N_s, \rho_s$ cyan solid line) in presence of the microalgae. Within this time the total number of particles in the chamber remains constant (purple solid line). Dashed/dash-dot lines are a fit to the first order kinetics (see Supplementary Section 6). Dotted lines indicate the initial values of $N_c$ and $\rho_c$. Dashed red line corresponds to the control experiment without microalgae.</figcaption></figure>

description of $q_y(y_j)$ given by Eq. (1) implies that

$$ \begin{aligned} m_1(y) &= L_+(y) - L_-(y) \\ m_2(y) &= 2 \frac{L_+(y) + L_-(y)}{L_+(y) + L_-(y)}. \end{aligned} \quad (5) $$

Given that $L_+(y)$ themselves are well approximated by exponentials, this affords analytical approximations also for $m_{1,2}(y)$ (Fig. 2c-Inset). Armed with this description we can now compare the experimental curves with the prediction from the approximate jump-diffusion model. The olive-coloured curve in Fig. 3 shows that the theoretical distributions recapitulate extremely well the experimental ones. On one hand, this provides a justification a posteriori for the modelling approach taken above. On the other it reveals that, at a continuum level, the dynamics of colloids in an active bath can be reduced to two quantities: the first and the second moments of the active displacements. This provides an intuitive way to conceptualise the complex dynamics of colloidal particles within an active suspension.

## Induced de-mixing of the passive particles

What we learned so far on the colloids' emergent active dynamics can be harnessed to induce the de-mixing of our active-passive suspensions, as a first step towards more complex control of microscale cargo transport. The idea is to include a confining boundary which acts as a kind of selective membrane, letting the beads cross but not the algae. This is achieved by considering circular chambers decorated with 190 µm-long side channels whose width $d$ (Fig. 4a) is smaller than the swimmer diameter but larger than the bead diameter (here taken as $2a = 6$ µm $< d = 7–8$ µm $< 2R = 8–10$ µm). Due to the active dynamics of the colloids, we expect them to be pushed inside the side channels by the algae, therefore depleting the particles in the main chamber and spatially separating the active and passive species. As shown Fig. 4b (see also Supplementary Movie 4) this is exactly what we obtain, with the number of particles in the circular chamber $N_c$ (red solid line) quickly decreasing and the number in the side channels $N_s$ (cyan solid line) increasing. Notice that the total number of particles $N_t = N_c + N_s$ in the system (purple solid line) remains constant over the duration of the experiment. The process follows a first order kinetic with rates $k_{in}$.
and $k_{out}$ for transitions towards or out-of the side channels respectively (see Supplementary Section 6). The occupancy numbers follow an exponential relaxation with time-scale $1/(k_{in} + k_{out}) = 77 \pm 4$ min and stationary values $N_c/N_t = k_{out}/(k_{in} + k_{out}) = 0.64 \pm 0.02$ and $N_s/N_t = k_{in}/(k_{in} + k_{out}) = 0.36 \pm 0.02$, leading to estimates for the two rates of $k_{in} = (7.8 \pm 0.6) \times 10^{-5}$ s⁻¹ and $k_{out} = (13.9 \pm 0.8) \times 10^{-5}$ s⁻¹. The latter results from the passive dynamics of the colloids within the side channels, which depends on their thermal diffusivity, potential electrostatic interactions with the boundaries, short-range hydrodynamic interactions between the beads⁷⁰, channel design, etc⁷¹. The former, instead, is the consequence of active colloidal displacements. Assuming no colloid/wall attraction within the side channels (as confirmed by our control experiment, see below), it encapsulates the intrinsic out-of-equilibrium nature of the system, a property which is immediately clear from the significant difference in steady-state densities between the side channels ($\rho_s = 0.0295 \pm 0.007$ beads/µm²) and in the main chamber ($\rho_c = 0.0019 \pm 0.0001$ beads/µm²) (Fig. 4c). The value of $k_{in}$ can be compared with an estimate derived from the effective colloidal dynamics from Eq. (3) (see Supplementary Section 7). The mean first passage time of a colloid to the entrance of any of the side channels returns a rate of $k_{in}^{V,D} = (7.408 \pm 0.1) \times 10^{-5}$ s⁻¹ which agrees very well with the experimental value. This shows that the simplified advection-diffusion model for the passive particles is a good description not just for time-independent properties like the steady-state distribution of Fig. 3, but also for dynamic ones like the rate of active cargo transport to the side channels. It is instructive also to use the model to explore the contributions of the different parts of the dynamics. Hence we see that setting $V_{eff} = 0$ and $D_{eff} = 3.55$ µm² s⁻¹, the effective diffusivity far from the boundaries, one obtains a rate $k_{in}^D = (16.2 \pm 4.1) \times 10^{-5}$ s⁻¹. The stark difference between $k_{in}^{V,D}$ and $k_{in}^D$ suggests that, in our case, it is not sufficient to know the behaviour of the active-passive suspension in the bulk of the chamber. The behaviour close to the boundary, where the inhomogeneities are found, is essential to grasp the dynamics of the de-mixing process. Repeating the same experiment only with colloids shows indeed no appreciable filling of the side channels within the duration of the experiment (-10 h; see Supplementary Movie 5) and the particle density in the main chamber remaining constant (Fig. 4b red dashed line).

Although the microfluidic chip in Fig. 4 leads to a much higher concentration of colloids in the side channels than in the main chamber, this still involves only ~36% of the colloids. A different design, however, could maximise the fraction of de-mixed colloids instead of their concentration. For example, modifying the current design to have two identical chambers connected by a narrow gap of size $d = 7-8$ µm should produce a fraction of de-mixed colloids given by $k_{in}/(k_{in} + k_{in}^{th}) \approx 95\%$. The design could then be tailored to different needs.

## Discussion

In this study we used an active-passive system as a minimal model to investigate the active transport of passive microscopic objects. Leveraging the well-known boundary accumulation of microswimmers, we saw that the emerging interactions between active and passive components are also heterogeneous and anisotropic. In turn, they lead to a complex spatial distribution of colloids and to the possibility of de-mixing the suspension. These properties are well described by a simple advection-diffusion model for the passive species, based only on the first and second moments of their ensuing active displacements. This type of approximation, then, provides a degree of universality in the description of the dynamics of active systems including purely active and active-passive species. From the targeted delivery of chemicals, biomarkers or contrast agents at specific locations in the body to the remediation of water and soils, artificial and biological microscopic self-propelled particles hold a great potential to address critical issues in areas ranging from personalised medical care to environmental

sustainability⁷²,⁷³. Development in these areas will depend on our ability to use and control micro/nano-swimmers to perform essential basic functions, such as sensing, collecting and delivering passive cargoes in an autonomous, targeted and selective way⁷⁴. Within this perspective, we have shown here that heterogeneous activity fields can be employed to control the fate of colloidal objects in a manner similar to biological systems, where the spatial modulation of the activity can control, for example, the positioning of organelles within cells¹⁴⁻¹⁶. Although this was achieved here through simple confinement, many microswimmers possess -or can be designed to have- the ability to respond to a range of external physico-chemical cues or stimuli⁵¹⁻⁵⁴. We expect that similar phenomenologies should also emerge whenever the active particles in a mixed system autonomously generate a space-dependent concentration profile. Our work paves the way to using these to dynamically manipulate the fate of colloidal cargoes by externally altering micro-swimmers' dynamics. Advance in this area will require to understand how to predict the active displacements of passive cargoes directly from the dynamics of microswimmers⁷⁵, a theoretical development that we leave for a future study.

## Methods

### Microorganism culturing

Cultures of CR strain CC-125 were grown axenically in a Tris-Acetate-Phosphate medium⁷⁶ at 21 °C under periodic fluorescent illumination (16h light/8h dark cycles, 100 µM² s⁻¹, OSRAM Fluora). Cells were harvested at -5 × 10⁶ cells/ml in the exponentially growing phase, then centrifuged at 10 × g for 10 min and the supernatants replaced by DI-water containing the desired concentration of polystyrene colloids (Polybead® Carboxylate Microspheres, radius $a = (5 \pm 0.5)$ µm, Polysciences).

### Microfluidics and microscopy

The mixed colloid-swimmer solution was injected into polydimethylsiloxane-based microfluidic channels, previously passivated with a 0.5% w/w Pluronic-127 solution to prevent sticking of the particles to the surfaces of the chips. The microfluidic devices were then sealed with Vaseline to prevent evaporation. The channels of varying widths were observed under bright-field illumination on a Zeiss AxioVert 200M inverted microscope. A long-pass filter (cutoff wavelength 765 nm) was added to the optical path to prevent photo-tactic responses of the cells. For each channel width stacks of 80,000 images at ×10 magnification were acquired at 10 fps with an IDS UI-3370CP-NIR-GL camera.

### Data analysis: trajectories

Particle Trajectories were digitised using a standard Matlab particle tracking algorithm (The code can be downloaded at http://people.umass.edu/kilfoil/downloads.html). After the particle trajectories were reconstructed, the large jumps were isolated using a combination of move size and directional correlation along the trajectory (see detailed protocol in⁴⁷,⁶³).

### Data analysis: fitting of $L_+(y)$, $L_-(y)$ and $\lambda(y)$

The experimental curves $L_+(y)$, $L_-(y)$ and $\lambda(y)$ have been fitted by the following analytical functions:

$$L_+(y) = (L_+^w - L^b)e^{-\frac{y-y_{col}}{l^b}} + L^b, \quad (6)$$

$$L_-(y) = (L_-^w - L^b)e^{-\frac{y-y_{col}}{l^b}} + L^b, \quad (7)$$

$$\lambda(y) = (\lambda^w - \lambda^b)e^{-\frac{y-y_{col}}{l^b}} + \lambda^b, \quad (8)$$
<table><caption>**Table 1 | Fit parameters (µm) for Eqs. (6), (7) and (8) and their standard error**</caption><thead><tr><th>L<sup>b</sup></th><th>l<sup>+</sup></th><th>l<sup>-</sup></th><th>λ<sup>w</sup></th><th>λ<sup>b</sup></th><th>l<sup>a</sup></th></tr></thead><tbody><tr><td>3.70</td><td>5.05</td><td>10.8</td><td>0.237</td><td>0.315</td><td>8.20</td></tr><tr><td>0.079</td><td>0.764</td><td>0.978</td><td>0.0046</td><td>0.0040</td><td>1.278</td></tr></tbody></table>

<table><caption>**Table 2 | Limits in y<sub>j</sub> for the y-dependent distributions of jumps perpendicular to the channel boundaries**</caption><tbody><tr><td>y(µm) (bin centre)</td><td>5.0</td><td>8.3</td><td>11.7</td><td>15.0</td><td>18.4</td><td>21.7</td><td>25.0</td></tr><tr><td>lim(y) (µm)</td><td>7</td><td>12</td><td>12</td><td>15</td><td>15</td><td>18</td><td>20</td></tr><tr><td>y(µm) (bin centre)</td><td>28.4</td><td>31.7</td><td>35.0</td><td>38.4</td><td>41.7</td><td>45.0</td><td>48.4</td></tr><tr><td>lim(y) (µm)</td><td>20</td><td>20</td><td>20</td><td>20</td><td>20</td><td>20</td><td>20</td></tr></tbody></table>

with $L_+^w = L_-^w = 0.375$ µm fixed. The best fits (dotted lines in Figs. 2 and 3) are given by the parameters in Table 1.

The experimental curves $L_+(y)$ and $L_-(y)$ were obtained from the experimental jump size distributions $q_j(y_j)$ by fitting the logarithm of each side (i.e., positive and negative $y_j$) with affine functions. Because these distributions deviate from simple exponentials as we get closer to the boundary (Fig. 2d), the fits were restricted to $|y_j|<\lim(y)$ in order to extract the characteristic length scales of the initial decay of the distributions. The resulting fits are shown in Fig. 2d (solid lines). The limits for the fits of the experimental distributions were fixed as in Table 2.

## Numerical simulations

One hundred trajectories composed of 10,000 time-steps (with $\delta t = 0.1$s, equal to the acquisition period of the experiments) were simulated for each condition. At each time step an acceptance-rejection scheme is performed to decide whether the particle undergoes a jump: (i) a (uniform) random number is first drawn from the interval [0,1], (ii) if this number sits within the interval $[(1-\lambda(y)\delta t)/2; (1+\lambda(y)\delta t)/2]$, the particle undergoes a jump whose size is taken from the experimental distribution $q_j(y_j)$ (using inverse transform sampling), (iii) otherwise the particle simply undergoes Brownian motion with diffusivity $D_0 = 0.0439 \ \mu m^2 \ s^{-1}$ (the theoretical bulk diffusivity for a 5 µm-radius particle at room temperature). Upon reaching a boundary jump moves undergo a stopping boundary condition, stopping at $\pm 5 \ \mu m$ (the particle's radius) from the boundary (this is to be as faithful with our experimental observations as possible). For simplicity, diffusive moves that cross the boundary undergo a reflective boundary condition. More accurate boundary implementations for the diffusive motion do not lead to different results as the impact of thermal diffusivity in the channel experiments is negligible.

## Data availability

The main data used in this study are available in the Zenodo database with https://doi.org/10.5281/zenodo.6866462.

## Code availability

The codes used in this study are available in the Zenodo database with https://doi.org/10.5281/zenodo.6866462.

## References

1. Civelekoglu-Scholey, G. & Scholey, J. M. Mitotic force generators and chromosome segregation. *Cell. Mol. Life Sci* **67**, 2231 (2010).

2. Honerkamp-Smith, A. R. et al. The noisy basis of morphogenesis: mechanisms and mechanics of cell sheet folding inferred from developmental variability. *PLoS Biol.* **16**, 1 (2018).

3. Tsimring, L. S. Noise in biology. *Rep. Prog. Phys* **77**, 026601 (2014).

4. Ramaswamy, S. The mechanics and statistics of active matter. *Annu. Rev. Condens. Matter Phys* **1**, 323 (2010).

5. Vicsek, T., Czirk, A., Ben-Jacob, E., Cohen, I. & Shochet, O. Novel type of phase transition in a system of self-driven particles. *Phys. Rev. Lett* **75**, 1226 (1995).

6. Bricard, A., Caussin, J.-B., Desreumaux, N., Dauchot, O. & Bartolo, D. Emergence of macroscopic directed motion in populations of motile colloids. *Nature* **503**, 95 (2013).

7. Giavazzi, F., Malinverno, C., Scita, G. & Cerbino, R. Tracking-free determination of single-cell displacements and division rates in confluent monolayers. *Front. Phys* **6**, 1 (2018).

8. Dombrowski, C., Cisneros, L., Chatkaew, S., Goldstein, R. E. & Kessler, J. O. Self-concentration and large-scale coherence in bacterial dynamics. *Phys. Rev. Lett* **93**, 098103 (2004).

9. Wensink, H. et al. Meso-scale turbulence in living fluids. *Proc. Natl. Acad. Sci. USA* **109**, 14308 (2012).

10. Wu, K. T. et al. Transition from turbulent to coherent flows in confined three-dimensional active fluids. *Science* **355**, 1705.02030 (2017).

11. Doostmohammadi, A., Shendruk, T. N., Thijssen, K. & Yeomans, J. M. Onset of meso-scale turbulence in active nematics. *Nat. Commun* **8**, 1 (2017).

12. Cates, M. E. & Tailleur, J. Motility-induced phase separation. *Annu. Rev. Condens. Matter Phys* **6**, 219 (2015).

13. Grobas, I., Polin, M. & Asally, M. Swarming bacteria undergo localized dynamic phase transition to form stress-induced biofilms. *eLife* **10**, 1 (2021).

14. Almonacid, M. et al. Active diffusion positions the nucleus in mouse oocytes. *Nat. Cell Biol* **17**, 470 (2015).

15. Cadot, B., Gache, V. & Gomes, E. R. Moving and positioning the nucleus in skeletal muscle-one step at a time. *Nucleus* **6**, 373 (2015).

16. Lin, C. et al. Active diffusion and microtubule-based transport oppose myosin forces to position organelles in cells. *Nat. Commun* **7**, 11814 (2016).

17. Be'er, A. & Ariel, G. A statistical physics view of swarming bacteria. *Mov. Ecol* **7**, 1 (2019).

18. Van Der Meer, B., Filion, L. & Dijkstra, M. Fabricating large two-dimensional single colloidal crystals by doping with active particles. *Soft Matter* **12**, 3406 (2016).

19. Ramananarivo, S., Ducrot, E. & Palacci, J. Activity-controlled annealing of colloidal monolayers. *Nat. Commun* **10**, 3380 (2019).

20. Kümmel, F., Shabestari, P., Lozano, C., Volpe, G. & Bechinger, C. Formation, compression and surface melting of colloidal clusters by active particles. *Soft Matter* **11**, 6187 (2015).

21. McCandlish, S. R., Baskaran, A. & Hagan, M. F. Spontaneous segregation of self-propelled particles with different motilities. *Soft Matter* **8**, 2527 (2012).

22. Stenhammar, J., Wittkowski, R., Marenduzzo, D. & Cates, M. E. Activity-induced phase separation and self-assembly in mixtures of active and passive particles. *Phys. Rev. Lett* **114**, 018301 (2015).

23. Takatori, S. C. & Brady, J. F. A theory for the phase behavior of mixtures of active particles. *Soft Matter* **11**, 7920 (2015).

24. Wysocki, A., Winkler, R. G. & Gompper, G. Propagating interfaces in mixtures of active and passive Brownian particles. *New J. Phys* **18**, 123030 (2016).

25. Weber, S. N., Weber, C. A. & Frey, E. Binary mixtures of particles with different diffusivities demix. *Phys. Rev. Lett* **116**, 058301 (2016).

26. Ma, Z., Lei, Q. L. & Ni, R. Driving dynamic colloidal assembly using eccentric self-propelled colloids. *Soft Matter* **13**, 8940 (2017).

27. Jahanshahi, S., Lozano, C., TenHagen, B., Bechinger, C. & Löwen, H. Colloidal Brazil nut effect in microswimmer mixtures induced by motility contrast. *J. Chem. Phys* **150**, 1 (2019).

28. Ilker, E. & Joanny, J.-F. Phase separation and nucleation in mixtures of particles with different temperatures. *Phys. Rev. Res.* **2**, 23200 (2020).
29. Rogel Rodriguez, D., Alarcon, F., Martinez, R., Ramírez, J. & Valeriani, C. Phase behaviour and dynamical features of a two-dimensional binary mixture of active/passive spherical particles. *Soft Matter* **16**, 1162 (2020).

30. Kaiser, A. et al. Transport powered by bacterial turbulence. *Phys. Rev. Lett* **112**, 158101 (2014).

31. Razin, N., Voituriez, R., Elgeti, J. & Gov, N. S. Generalized Archimedean's principle in active fluids. *Phys. Rev. E* **96**, 1 (2017).

32. Pietzonka, P., Fodor, É., Lohrmann, C., Cates, M. E. & Seifert, U. Autonomous engines driven by active matter: energetics and design principles. *Phys. Rev. X* **9**, 41032 (2019).

33. Knežević, M. & Stark, H. Effective Langevin equations for a polar tracer in an active bath. *New J. Phys* **22**, 113025 (2020).

34. Palacci, J., Sacanna, S., Vatchinsky, A., Chaikin, P. M. & Pine, D. J. Photoactivated colloidal dockers for cargo transportation. *JACS* **135**, 15978 (2013).

35. Koumakis, N., Lepore, A., Maggi, C., & Leonardo, R. Di. Targeted delivery of colloids by swimming bacteria. *Nat. Commun.* **4**, 2588 (2013).

36. Wang, L. & Simmchen, J. Review: interactions of active colloids with passive tracers. *Condens. Matter* **4**, 78 (2019).

37. Mathijssen, A. J. T. M., Guzmán-Lastra, F., Kaiser, A. & Löwen, H. Nutrient transport driven by microbial active carpets. *Phys. Rev. Lett.* **121**, 248101 (2018).

38. F. Guzmán-Lastra, H. Löwen, and A. J. Mathijssen. Active carpets drive non-equilibrium diffusion and enhanced molecular fluxes. *Nat. Commun.* **12**, 1906 (2021).

39. Wu, X. L. & Libchaber, A. Particle diffusion in a quasi-two-dimensional bacterial bath. *Phys. Rev. Lett* **84**, 3017 (2000).

40. Leptos, K. C., Guasto, J. S., Gollub, J. P., Pesci, A. & Goldstein, R. E. Dynamics of enhanced tracer diffusion in suspensions of swimming eukaryotic microorganisms. *Phys. Rev. Lett* **103**, 198103 (2009).

41. Valeriani, C., Li, M., Novosel, J., Arlt, J. & Marenduzzo, D. Colloids in a bacterial bath: simulations and experiments. *Soft Matter* **7**, 5228 (2011).

42. Miño, G. et al. Enhanced diffusion due to active swimmers at a solid surface. *Phys. Rev. Lett* **106**, 048102 (2011).

43. Kurtuldu, H., Guasto, J. S., Johnson, K. A. & Gollub, J. P. Enhancement of biomixing by swimming algal cells in two-dimensional films. *Proc. Natl. Acad. Sci. USA* **108**, 10391 (2011).

44. Miño, G., Dunstan, J., Rousselet, A., Clément, E. & Soto, R. Induced diffusion of tracers in a bacterial suspension: theory and experiments. *J. Fluid Mech.* **729**, 423 (2013).

45. Jepson, A., Martinez, V. A., Schwarz-Linek, J., Morozov, A. & Poon, W. C. Enhanced diffusion of nonswimmers in a three-dimensional bath of motile bacteria. *Phys. Rev. E* **88**, 041002(R) (2013).

46. Patteson, A., Gopinath, A., Purohit, P. K. & Arratia, P. E. Particle diffusion in active fluids is non-monotonic in size. *Soft. Matter* **12**, 2365 (2016).

47. Jeanneret, R., Pushkin, D. O., Kantsler, V. & Polin, M. Entrainment dominates the interaction of microalgae with micron-sized objects. *Nat. Commun.* **7**, 12518 (2016).

48. DiLeonardo, R. et al. Bacterial ratchet motors. *Proc. Natl. Acad. Sci. USA* **107**, 9541 (2010).

49. Sokolov, A., Apodaca, M. M., Grzybowski, B. A. & Aranson, I. S. Swimming bacteria power microscopic gears. *Proc. Natl. Acad. Sci. USA* **107**, 969 (2010).

50. Palacci, J., Sacanna, S., Steinberg, A. P., Pine, D. J. & Chaikin, P. M. Living crystals of light-activated colloidal surfers. *Science* **339**, 936 (2013).

51. Kirkegaard, J. B., Bouillant, A., Marron, A. O., Leptos, K. C. & Goldstein, R. E. Aerotaxis in the closest relatives of animals. *eLife* **5**, e18109 (2016).

52. Arrieta, J., Barreira, A., Chioccioli, M., Polin, M. & Tuval, I. Phototaxis beyond turning: persistent accumulation and response acclimation

of the microalga Chlamydomonas reinhardtii. *Sci. Rep* **7**, 3447 (2017).

53. Frangipane, G. et al. Dynamic density shaping of photokinetic E. Coli. *eLife* **7**, e36608 (2018).

54. Arlt, J., Martinez, V. A., Dawson, A., Pilizota, T. & Poon, W. C. K. Painting with light-powered bacteria. *Nat. Commun.* **9**, 768 (2018).

55. Lozano, C., tenHagen, B., Löwen, H. & Bechinger, C. Phototaxis of synthetic microswimmers in optical landscapes. *Nat. Commun.* **7**, 12828 (2016).

56. Kantsler, V., Dunkel, J., Polin, M. & Goldstein, R. E. Ciliary contact interactions dominate surface scattering of swimming eukaryotes. *Proc. Natl. Acad. Sci. USA* **110**, 1187 (2013).

57. Contino, M., Lushi, E., Tuval, I., Kantsler, V. & Polin, M. Microalgae scatter off solid surfaces by hydrodynamic and contact forces. *Phys. Rev. Lett* **115**, 258102 (2015).

58. Ostapenko, T. et al. Curvature-guided motility of microalgae in geometric confinement. *Phys. Rev. Lett* **120**, 068002 (2018).

59. M. Contino, Characterisation and control of the dynamical properties of swimming microorganisms under confinement, *PhD thesis*, Warwick, 2017.

60. S. Williams, R. Jeanneret, I. Tuval, and M. Polin, Supplementary data and codes for "Confinement-induced accumulation accumulation and de-mixing of microscopic active-passive mixtures", *Zenodo https://zenodo.org/record/6866462#.YvMOx3ZBy3A* (2022).

61. Lançon, P., Batrouni, G., Lobry, L. & Ostrowsky, N. Drift without flux: Brownian walker with a space-dependent diffusion coefficient. *Europhys. Lett* **54**, 0005004 (2001).

62. Mathijssen, A. J., Jeanneret, R. & Polin, M. Universal entrainment mechanism controls contact times with motile cells. *Phys. Rev. Fluids* **3**, 033103 (2018).

63. Mosby, L. S., Polin, M. & Köster, D. V. A Python based automated tracking routine for myosin II filaments. *J. Phys. D: Appl. Phys* **53**, 304002 (2020).

64. Lagarde, A. et al. Colloidal transport in bacteria suspensions: from bacteria collision to anomalous and enhanced diffusion. *Soft Matter* **16**, 7503 (2020).

65. Denisov, S. I. & Bystrik, Y. S. Statistics of bounded processes driven by Poisson white noise. *Physica A* **515**, 38 (2019).

66. H. Risken & H. Haken, *The Fokker-Planck Equation: Methods of Solution and Applications Second Edition* (Springer, 1989).

67. Schnitzer, M. Theory of continuum random walks and application to chemotaxis. *Phys. Rev. E* **48**, 2553 (1993).

68. M. E. Cates and J. Tailleur. When are active Brownian particles and run-and-tumble particles equivalent? *Consequences for motility-induced phase separation*. *Eur. Phys. Lett.* **101**, 20010 (2013).

69. Martin, D. et al. Statistical mechanics of active Ornstein-Uhlenbeck particles. *Phys. Rev. E* **103**, 2008.12972 (2021).

70. Misiunas, K., Pagliara, S., Lauga, E., Lister, J. R. & Keyser, U. F. Nondecaying hydrodynamic interactions along narrow channels. *Phys. Rev. Lett* **115**, 038301 (2015).

71. Locatelli, E. et al. Single-file escape of colloidal particles from microfluidic channels. *Phys. Rev. Lett.* **117**, 038001 (2016).

72. Wang, J. & Gao, W. Nano/microscale motors: biomedical opportunities and challenges. *ACS Nano* **6**, 5745 (2012).

73. Gao, W. & Wang, J. The environmental impact of micro/nano-machines: a review. *ACS Nano* **8**, 3170 (2014).

74. Bechinger, C. et al. Active particles in complex and crowded environments. *Rev. Mod. Phys.* **88**, 045006 (2016).

75. Cammann, J. et al. Emergent probability fluxes in confined microbial navigation. *Proc. Natl. Acad. Sci. USA* **118**, 1 (2021).

76. E. S. Harris, *The Chlamydomonas sourcebook* (Elsevier, 2009).

## Acknowledgements

R.J. would like to thank Denis Bartolo, Jean-François Rupprecht and Jérémie Palacci for enlightening discussions. This work was
supported in part through the Ramón y Cajal Programme (RYC-2018-02534; MP), the Spanish Ministry of Science and Innovation (Grant PID2019-104232GB-100 funded by MCIN/AEI/ 10.13039/501100011033; M.P. and I.T.), the Margalida Comas Programme (PD/007/2016; R.J.), the Juan de la Cierva Incorporación Programme (IJCI-2017-32952; R.J.), the Junior Research Chair Programme (ANR-10-LABX-0010/ANR-10-IDEX-0001-02 PSL; R.J.).

**Author contributions**

R.J., I.T. and M.P. conceived and designed the experiments. R.J. and S.W. performed the experiments. R.J., S.W. and M.P. analysed the data. R.J. developed the analytical model and S.W. performed the numerical simulations. All authors discussed the results and wrote the article.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-32520-9.

**Correspondence** and requests for materials should be addressed to Marco Polin.

**Peer review information** *Nature Communications* thanks Giovanni Volpe and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**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/.

© The Author(s) 2022, corrected publication 2022