## Tunable active rotational diffusion in swimming droplets

Adrien Izzet, 1 Pepijn Moerman, 2 Jan Groenewold, 2,3 Jérôme Bibette, 4 and Jasna Bruji´ c 1, ∗

1

2 Debye Institute for Nanomaterials, Department of Chemistry, Utrecht University, Netherlands 3 Guangdong Provincial Key Laboratory of Optical Information Materials and Technology and Institute of Electronic Paper Displays, South China Academy of Advanced Optoelectronics, South China Normal University, Guangzhou 510006, P. R. China. 4 Laboratoire Colloïdes et Matériaux Divisés, CNRS UMR 8231, Chemistry Biology &amp; Innovation, ESPCI Paris, PSL Research University, 10 rue Vauquelin, 75005 Paris, France

Center for Soft Matter Research, Physics Department, New York University, New York, NY 10003, USA

(Dated: August 5, 2019)

Here we characterize the motility of athermal swimming droplets within the framework of active rotational diffusion. Just like active colloids, their trajectories can be modeled with a constant velocity V and a slow angular diffusion, but the random changes in direction are not thermally driven. Instead, V is determined by the interfacial tension gradient along the droplet surface, while local micellar fluctuations lead to droplet reorientation with a persistence time τ . We show that the origin of locomotion is the difference in the critical micellar concentration ∆ CMC in the front and the back of the droplet. Tuning this parameter by salt controls V from 3 -15 diameters d / s . Surfactant concentration has little effect on speed, but leads to a dramatic decrease in τ over four orders of magnitude. The corresponding range of the effective diffusion constant De extends beyond the realm of synthetic or living swimmers, in which V is limited by fuel consumption and τ is set by temperature or biological activity, respectively. Our tunable swimmers are ideal candidates for the study of the departure from equilibrium to high levels of activity, on both the single particle level and their collective behavior, including the motility-induced phase separation (MIPS).

Most active particles, such as Janus colloids [1], are anisotropic and their orientation determines the swimming direction. They consume fuel preferentially on one side of the particle, which results in chemical gradients that drive and direct self-propulsion. On short time scales, these swimmers move in a straight line, but they are small enough that thermal diffusion randomizes the swimming direction on long time scales. The transition between these two regimes is characterized by the persistence time, and they are modeled as Active Brownian Particles (ABPs) [2, 3]. Unlike Janus particles, active droplets are chemically isotropic, yet they spontaneously break symmetry as they dissolve to create swollen micelles [4]. At a sufficiently large Peclet number, the isotropic state becomes unstable to small perturbations, resulting in self-sustained motion in a given direction [5]. Autochemotaxis of the droplet is explained by an interfacial gradient along the surface, which induces a Marangoni stress and drives the flow of surfactants. This general mechanism is true for all droplet swimmers: direct emulsions made of oil [6] or nematic liquid crystals [7-10], as well as inverse emulsions [11, 12].

In the case of inverse emulsions, adding salt reduces the dissolution rate of the water into the oil and slows down the droplets, while increasing surfactant concentration speeds them up [11, 12]. It is thought that droplet motion is sufficient to cause the asymmetry needed to sustain a surfactant gradient, as the droplet encounters empty micelles in the front and leaves filled ones at the back [5, 13].

In this letter, we present an alternative origin for the driving force of Marangoni flows in direct emulsions. The difference between the CMC of the surfactant in pure water and in the presence of oil, ∆ CMC, creates a concentration gradient of monomers from the front to the back of a droplet, as it swims and dissolves, as shown in the schematic in Fig. 1. The re- sulting Marangoni flow gives rise to and sustains speeds up to 15 d / s , an order of magnitude faster than any known swimming droplets. We demonstrate that adding salt decreases ∆ CMC proportionally to the swimming speed, validating the mechanism. Conversely, we show that adding surfactants above the CMC has little effect on speed, since ∆ CMC is independent of surfactant concentration.

FIG. 1. Schematic of the mechanism at the origin of self-sustained droplet motion. Spontaneous symmetry breaking arises due to a higher concentration of monomers in equilibrium with the aqueous solution (on the right) compared to those facing DEP-saturated water (on the left). This effect arises because of the difference in the CMC in the two environments, and it causes a Marangoni flow, shown by the arrows, which is enhanced by droplet motion.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-c834bb66-ee57-4efd-852a-faf280b24f2f.jpg)

The direction of droplet motion is fixed by the orientation of the surfactant gradient, which varies in time due to local micellar fluctuations, as seen in the trajectories in Fig. 2. Because swimming droplets are athermal, ranging in size from 10 -300 µ m, their reorientation cannot be a thermally driven process and their motility is more complex than that of active colloids.

Depending on their size and fuel concentration, previous studies have empirically characterized trajectories of ordinary emulsions as ballistic (straight) or diffusive (random) [10, 12], and those of liquid crystal droplets as helical [8] and even os-

FIG. 2. Individual trajectories of 30 µm diameter droplets, lasting four seconds, swimming in aqueous solutions with increasing bulk SDS concentrations. The motion transitions from ballistic to diffusive, with more frequent turns leading to a progressively lower effective diffusion constant, which is still orders of magnitude higher than that of thermal rotational diffusion.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-9f903907-31d3-487a-a497-4d73c0c5c69a.jpg)

cillatory [10, 14]. Despite all these advances, a quantitative analysis of the motility profiles as a function of experimental control parameters, as well as their theoretical interpretation, is still lacking in the literature.

Here we show that single-particle dynamics in direct emulsions can be described by active rotational diffusion, in which velocity V and persistence time τ are independently tuned by the salt and surfactant concentration, respectively. Since the microscopic mechanism that underlies turning is active rather than thermal, we are able to tune τ over four orders of magnitude. As a result, the effective diffusion constant is millions of times higher than that of passive rotational diffusion. Swimming droplets therefore belong to the class of ABPs, but with a much wider effective temperature range than active colloids or even bacteria [15, 16]. This system opens the possibility to study collective motion, MIPS, and giant fluctuations, which are predicted to occur at very high levels of activity [17]. While the ABP model holds for small droplets &lt; 50 µm, we note that more exotic effects, including self-avoidance [13] and oscillations [10], are observed for larger emulsions, and are beyond the scope of this study.

Droplets of diethyl phthalate (DEP) are produced in a custom-built glass microfluidic setup with a coaxial co-flow geometry [18-20]. The continuous phase is initially a 10mM sodium dodecyl sulphate (SDS) solution saturated with DEP to prevent autophoretic motion and stabilize the droplets. To initiate swimming, 0 5 µL of the emulsion is injected using a . micropipette into a circular glass sample cell (1cm wide and 1mm deep), filled with SDS solution at the desired concentration. Droplet trajectories are observed using a Nikon Eclipse T microscope equipped with a 4 i × Plan Apochromat λ Nikon objective ( NA = 0 20). The movies are recorded using a Zyla . sCMOS Andor camera operating at 40 frames per second. To compare droplets of the same size, we crop the trajecto-

FIG. 3. Instantaneous speed V at a constant [SDS] = 10 mM sharply decreases with [NaCl]. This trend is theoretically captured by Eq. 1 with ∆ CMCvalues independently measured by pendant drop experiments. Above [NaCl] = 70 mM, corresponding to ∆ CMC = 0 5mM, . droplets cease to swim. Dashed line is computed from the difference between empirical fits of CMC vs [SDS] [23]).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-29460ff7-5271-43c5-aa2c-bed88d3cf148.jpg)

ries in a time window corresponding to a droplet diameter of 30 ± 5 µm.

More quantitatively, the droplet velocity is given by /vector V = M ∇ || C , where M is the phoretic mobility and ∇ || C is the average concentration gradient along the surface of the droplet. It has been shown that Marangoni effects are dominant compared to diffusiophoresis in swimming droplets [24]. The velocity is thus given by:

Our system therefore consists of droplets of DEP suspended in an aqueous solution of SDS. This oil is slightly soluble in water (1mM), which lowers the CMC of SDS from 8mM in pure water to 5mM in the vicinity of the droplets. This decrease corresponds to a decrease in ∆ G mic = 2 3 . kBT log ( ∆ CMC ) = 1 1 . kBT for the micellization energy of one monomer of SDS [21, 22]. The resulting DEP swollen micelles are thermodynamically more stable than the empty ones. Therefore, an infinitesimal perturbation in the velocity of the droplet leads to an anisotropy in surfactant concentration, on the order of ∆ CMC, on the droplet interface: the part facing DEP-saturated water sheds more surface active monomers into solution than the part facing the empty micelles in the bulk, as shown in Fig. 1. This asymmetry induces a Marangoni stress that forces the surface flow of surfactants. At high enough dissolution rates compared to molecular diffusion, the interfacial gradient undergoes a swimming instability and the droplets self-propel [5].

$$s \, & \text{un} \, \text{given} \, \text{by} \colon \\ \vec { V } = \frac { r } { 2 \eta _ { 0 } + 3 \eta _ { i } } \left ( - \, \frac { \partial \gamma } { \partial c } \right ) \nabla _ { | | } C, \quad \text{(1)} \\ r - 1 & \text{$1$} \, \text{im} \, \text{the radine} \, of the \, \text{don} \, \text{let} \, \text{ $n$} \, - 1 $} \text{$3$} \, \text{$2$} \, \text{is}$$

where r = 15µm is the radius of the droplet, η i = 13mPa.s is the viscosity of the DEP oil, η o = 0 89mPa.s is the viscosity of . the solute, ∂γ / ∂ c is the rate of change of surface tension with SDS concentration at the interface, and ∇ || C = ∆ CMC 2 / r is controlled by ∆ CMC.

Using the pendant drop method, we measure the surface tension as a function of surfactant concentration, in the presence and absence of DEP, and find the difference in the CMC to be ∆ CMC = 3mM [25]. For the experimentally determined speed V ∼ 450 µm/s, shown in Fig. 3, we then calculate ∂γ / ∂ c = 0 01N/m/M. Comparing this value to the tension . gradient in the range of ∆ CMC at the interface with pure water [25], we find that it is two orders of magnitude smaller than the expected value. This discrepancy is consistent with the numerical prediction of Michelin [5] that one can only achieve 1 -10% of the theoretical swimming speed in the Janus limit, depending on the Peclet number that relates the advective and diffusive surfactant transport. Alternatively, if the process is modeled as an adsorption at the interface, such that -∇ γ = kBT λ ∇ C [24], the measured speed predicts an interaction length λ = 5nm, which is in good agreement with the diameter of a SDS micelle of 3 7nm [26, 27], as proposed . by Izri et al. .

To test the hypothesis that the driving force for propulsion is ∆ CMCat the interface, we add sodium chloride (NaCl) to decrease the drive [28] and measure the resulting instantaneous velocity at a fixed surfactant concentration of 10mM. First, we estimate ∆ CMCby fitting the surface tension as a function of SDS and salt concentration using the equation proposed in [29] to reveal a gradual decrease in the CMC in both pure water and DEP-saturated water [25]. Nevertheless, the difference in the CMC decreases because the CMC decays faster in pure water.

Assuming all other terms in Eq. 1 remain the same, the change in ∆ CMCaccurately predicts the sharp decrease in the average instantaneous speed shown in Fig. 3. This assumption is corroborated by the fact that the gradients ∂γ / ∂ c below the CMC are largely independent of salt [25]. If we consider the surface adsorption picture, results are also consistent since the micelle size, which sets the characteristic length scale λ , hardly changes as a function of salt [27]. An important difference between experiment and theory is that droplets do not swim below ∆ CMC ≈ 0 5 mM, while the theoretical predic-. tion approaches a constant speed. This result shows that the surfactant gradient achieved by motility alone is not sufficient for self-sustained droplet motion, contrary to previous belief. Direct control of the speed via ∆ CMC not only presents a useful toolbox, but also lends support to the idea that micellar equilibration at the interface is at the origin of Marangoni flows.

Next, we investigate the effect of adding surfactants on the motility of individual droplets. Representative trajectories over four seconds are shown for various SDS concentrations in Fig. 2. Slightly above the CMC, droplets swim with a constant speed in a given direction. This speed is only approximately constant since droplets decelerate at a rate of 10 µ m/s 2 , but this effect is negligible over the course of the trajectory. As the concentration increases, droplets change their direction more frequently, and therefore explore less space. To quantify the motility, we compute the mean square displacement (MSD) ∆ L 2 of the trajectories, shown on a log-log

FIG. 4. Mean square displacement curves, averaged over tens of trajectories, show a decrease in the crossover time τ as a function of SDS concentration. The model of rotational diffusion in Eq. 2 gives an excellent fit to the data and the grey area corresponds to the 95% confidence interval.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-737e4c3a-7553-4965-a289-f4837a8d38d4.jpg)

plot in Fig. 4. We find that swimming droplets are ballistic at short times, where the slope is two, and they become diffusive at long times, where the slope is one. The crossover time τ between the two regimes decreases with increasing SDS concentration.

These results can be fit by the persistent random walk model [30],

$$\Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Philon _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi }
 \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Psi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Psi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi. } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi }\Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Psigma \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \ Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \psi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi_ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Psi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Psi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Psi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \Phi } \Phi \Phi \Phi \Phi \Phi _ { \Phi } \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi } \Phi _ { \Phi } \Phi } \Phi _ { \Phi } \$$

where the fitting parameters V and τ are shown in Figs. 5a,b. We find that V is largely independent of surfactant concentration, and in good agreement with direct measurements of the average instantaneous speed along the trajectories, as well as previous measurements on small droplets [10]. The speed is therefore also independent of the dissolution rate, which increases with SDS concentration, as shown in Fig. 5c. This result suggests that the equilibration timescale of the monomers at the interface with the bulk is faster than micellar solubilization. Droplet locomotion is governed by ∆ CMC and not the bulk micellar concentration, in support of the proposed mechanism for propulsion.

The surfactant concentration causes a decrease in τ that is dramatic: over one decade of SDS concentration, the fitting parameter τ in Eq. 2 decreases over four orders of magnitude, giving rise to a power law exponent of -4 in the fit in Fig. 5b. We compare the persistence time with the characteristic decorrelation time of the velocity-velocity correlation function C ( ∆ t ) in Fig. 6, averaged over the whole trajectory [31],

$$\Phi _ { \Phi } & = \left \langle \frac { \vec { v } ( t ) \cdot \vec { v } ( t + \Delta t ) } { | \vec { v } ( t ) | | \vec { v } ( t + \Delta t ) | } \right \rangle _ { t }. & & ( 3 )$$

FIG. 5. Effect of surfactant concentration on motility. a) Instantaneous speed V , estimated directly from trajectories (circles) or the MSD fits of Eq. 2 (squares), is largely independent of [SDS]. b) Persistence time τ , from the fits and the values extracted from the velocity correlation function C ( ∆ t ) , sharply decays with [SDS] with a heuristic power law fit τ = 6 × 10 6 [SDS] -4 . c) The rate of droplet dissolution increases nonlinearly as a function of SDS concentration over the experimental range.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-5087058e-99ec-4905-b44c-5282aeae5704.jpg)

The data is fit reasonably well with an exponential decay, C ( ∆ t ) = exp ( -∆ t / τ ) , to obtain the characteristic times over which the swimming direction randomizes, also shown in Fig. 4b. The agreement between the two estimates of τ confirm the trend that droplets turn more frequently as more surfactant is added. This result differs from other swimmers, in which the persistence time is tuned by temperature, i.e. particle size [32], or by biological processes in the case of living organisms [16]. This limits the range of reported τ to 0 1 . -1s for colloids with thermally activated particle sizes. The measured active rotational diffusion constants De = V 2 τ are on timescales that are at least one million times faster than thermal diffusion with D ≈ 10 -2 µm 2 /s.

the droplet, eventually leading to chaos in the hydrodynamics of the droplet. Further studies on the spatial heterogeneities in the advective flows are needed to rationalize the turning frequency characterized by the rapid power law decay in the persistence time in Fig. 4b.

The statistics of this active turning process are captured by the distributions of the instantaneous speed and reorientation angle as a function of SDS concentration in Fig. 7. For a given δ t = 25ms, the displayed position at time t + δ t along the trajectory informs on both the speed and change in direction with respect to the previous timestep. Since the position at time t lies at the origin, positions on the positive x-axis indicate straight motion, while particles that turn result in an off-axis position. As the SDS concentration increases, the characteristic time τ decreases to approach δ t . When τ /greatermuch δ t , we observe ballistic trajectories focused on a single point in panel (a), where the noise is on the order of a pixel size. Intermediate concentrations give rise to noise in both speed and reorientation angle, but the droplets rarely turn back on themselves (panel (b)). Eventually, at τ ≈ δ t droplets sample the space with equal probability,similar to Brownian diffusion, as indicated by the color map in panel (c).

To explain this active rotational diffusion, we estimate that the increase in local micellar fluctuations surrounding the droplets creates gradients that are insufficient to reorient the droplets according to Eq.1. Instead of a microscopic mechanism for turning, it has been shown theoretically and numerically that a shift in motility from a Marangoni flow to diffusiophoresis [33-35] gives rise to the spontaneous breaking of symmetry of the inner and outer flow fields surrounding

In this paper, we have presented evidence that the difference in the CMC at the front and back of the droplet is at the origin of droplet motion, which is why the average instantaneous velocity is independent of surfactant concentration above the CMC, but sensitive to changes in salt. Conversely, the frequency with which droplets change direction drastically increases with surfactant concentration, allowing us to independently tune V and τ in the experiments. This allows one to achieve effective diffusion constants that are millions of times higher than those induced by passive rotational diffusion of active particles. These swimmers broaden the range of active particle sizes to the athermal realm. Such a system promises to be useful in tuning the departure from equilibrium [36, 37] and understanding the resulting dynamics [2, 38, 39]. The self-assembly of droplets into equilibrium colloidomers [40] can now be extended to the assembly of non-equilibrium architectures, akin to living systems [41].

The authors gratefully acknowledge insights by Eric Vanden Eijnden, Katherine Newhall, Alexander Grosberg, Sascha Hilgenfeldt, Willem Kegel, and Alfons Van Blaaderen, as well as useful conversations with Micha Kornreich and Rodrigo Guerra. This work was supported by the Materials Research Science and Engineering Center (MRSEC) program of the National Science Foundation under Grants No. NSF DMR1420073, No. NSF PHY17-48958, and No. NSF DMR1710163.

∗ email: jb2929@nyu.edu ; https://wp.nyu.edu/brujiclab/

[1] J. R. Howse, R. a. L. Jones, A. J. Ryan, T. Gough, R. Vafabakhsh, and R. Golestanian, Phys. Rev. Lett. 99 , 8 (2007).

FIG. 6. Normalized auto-correlation of the velocity for each concentration. Each dataset is computed from averaging the C ( ∆ t ) of experimental trajectories at a given concentration. Solid lines are fit with an exponential decay following equation 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-67d0bf62-1a63-47ef-8aa8-862439c13dd6.jpg)

- [2] A. P. Solon, M. E. Cates, and J. Tailleur, Eur. Phys. J. ST 224 , 1231 (2015).
- [4] C. A. Weber, D. Zwicker, F. Jülicher, and C. F. Lee, Rep. Prog. Phys. 82 , 064601 (2019).
- [3] M. E. Cates and J. Tailleur, Europhys. Lett. 101 , 20010 (2013).
- [5] S. Michelin, E. Lauga, and D. Bartolo, Physics of Fluids 25 , 061701 (2013).
- [7] K. Peddireddy, P. Kumar, S. Thutupalli, S. Herminghaus, and C. Bahr, Langmuir 28 , 12426 (2012).
- [6] P. G. Moerman, H. W. Moyses, E. B. Van Der Wee, D. G. Grier, A. Van Blaaderen, W. K. Kegel, J. Groenewold, and J. Brujic, Phys. Rev. E 96 , 032607 (2017).
- [8] C. Krüger, G. Klös, C. Bahr, and C. C. Maass, Phys. Rev. Lett. 117 , 048003 (2016).
- [10] M. Suga, S. Suda, M. Ichikawa, and Y. Kimura, Phys. Rev. E 97 , 1 (2018).
- [9] N. Ueno, T. Banno, A. Asami, Y. Kazayama, Y. Morimoto, T. Osaki, S. Takeuchi, H. Kitahata, and T. Toyota, Langmuir 33 , 5393 (2017).
- [11] S. Thutupalli, R. Seemann, and S. Herminghaus, New Journal of Physics 13 , 073021 (2011).
- [13] C. Jin, C. Krüger, and C. C. Maass, Proc. Natl. Acad. Sci. U.S.A. 114 , 5089 (2017).
- [12] Z. Izri, M. N. Van Der Linden, S. Michelin, and O. Dauchot, Phys. Rev. Lett. 113 , 1 (2014).
- [14] S. Herminghaus, C. C. Maass, C. Krüger, S. Thutupalli,
- L. Goehring, and C. Bahr, Soft Matter 10 , 7008 (2014).
- [16] A. E. Patteson, A. Gopinath, M. Goulian, and P. E. Arratia, Sci. Rep. 5 , 15761 EP (2015).
- [15] J. Saragosti, P. Silberzan, and A. Buguin, PLOS ONE 7 , 1 (2012).
- [17] J. U. Klamser, S. C. Kapfer, and W. Krauth, Nat. Commun. 9 , 1 (2018).
- [19] C. H. Choi, D. a. Weitz, and C. S. Lee, Advanced Materials 25 , 2536 (2013).
- [18] C. Cramer, P. Fischer, and E. J. Windhab, Chemical Engineering Science 59 , 3045 (2004).
- [20] M. F. Haase and J. Brujic, Angewandte Chemie - International Edition 10.1002/anie.201406040 (2014).
- [22] M. J. Rosen and J. T. Kunjappu, Surfactants and Interfacial Phenomena: Fourth Edition (John Wiley &amp; Sons Inc., 2012).
- [21] J. N. Israelachvili, Intermolecular and Surface Forces: Third Edition (Academic Press by Elsevier, 2011).
- [23] S. A. Morton, D. J. Keffer, A. N. Davis, and R. M. Counce, Separation Science and Technology 43 , 310 (2008).
- [25] See Supplementary Material.
- [24] J. Anderson, Ann. Rev. Fluid Mech. 21 , 61 (1989).
- [26] S. Hayashi and S. Ikeda, J. Phys. Chem. 84 , 744 (1980).
- [28] B. Principles, Emulsion Science , Springer Tracts in Modern Physics, Vol. 181 (Springer Berlin Heidelberg, Berlin, Heidelberg, 2002).
- [27] G. Duplâtre, M. F. Ferreira Marques, and M. da Graça Miguel, J. Phys. Chem. 100 , 16608 (1996).
- [29] B. von Szyszkowski, Zeitschrift für Physikalische Chemie 64 , 385 (1908).
- [31] S. Burov, S. M. A. Tabei, T. Huynh, M. P. Murrell, L. H. Philipson, S. a. Rice, M. L. Gardel, N. F. Scherer, and A. R. Dinner, Proc. Natl. Acad. Sci. U.S.A. 110 , 19689 (2013).
- [30] Y. Fily and M. C. Marchetti, Phys. Rev. Lett. 108 , 1 (2012).
- [32] J. Palacci, S. Sacanna, A. P. Steinberg, D. J. Pine, and P. M. Chaikin, Science 579 (2012).
- [34] M. Morozov and S. Michelin, J. Chem. Phys. 150 , 044110 (2019).
- [33] M. Schmitt and H. Stark, Physics of Fluids 28 , 012106 (2016).
- [35] Chamolly, Alexander and Lauga, Eric, Eur. Phys. J. E 42 , 88 (2019).
- [37] E. Fodor, C. Nardini, M. E. Cates, J. Tailleur, P. Visco, and F. van Wijland, Phys. Rev. Lett. 117 , 038103 (2016).
- [36] G. Gonnella, D. Marenduzzo, A. Suma, and A. Tiribocchi, C. R. Physique 16 , 316 (2015).
- [38] H. Stark, Acc. Chem. Research 51 , 2681 (2018).
- [40] A. McMullen, M. Holmes-Cerfon, F. Sciortino, A. Y. Grosberg, and J. Brujic, Phys. Rev. Lett. 121 , 138002 (2018).
- [39] S. Thutupalli, D. Geyer, R. Singh, R. Adhikari, and H. A. Stone, Proc. Natl. Acad. Sci. U.S.A. , 201718807 (2018).
- [41] Z. Zeravcic, V. N. Manoharan, and M. P. Brenner, Rev. Mod. Phys. 89 , 1 (2017).

FIG. 7. Displacements of the droplet between consecutive time-steps ( i.e. ∆ t = 25ms), for three surfactant concentrations. Colors represent the density at which this location in the displacement field is hit along the trajectory.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-2d48483c-5e88-42a3-8ed3-05346df70740.jpg)

## Supplementary material for Tunable active rotational diffusion in swimming droplets

Adrien Izzet, 1 Pepijn Moerman, 2 Jan Groenewold, 2,3 Jérôme Bibette, 4 and Jasna Bruji´ c 1

1 Center for Soft Matter Research, Physics Department, New York University, New York, NY 10003, USA 2 Debye Institute for Nanomaterials, Department of Chemistry, Utrecht University, Netherlands 3 Guangdong Provincial Key Laboratory of Optical Information Materials and Technology and Institute of Electronic Paper Displays, South China Academy of Advanced Optoelectronics, South China Normal University, Guangzhou 510006, P. R. China. 4 Laboratoire Colloïdes et Matériaux Divisés, CNRS UMR 8231, Chemistry Biology &amp; Innovation, ESPCI Paris, PSL Research University, 10 rue Vauquelin, 75005 Paris, France

## SURFACE TENSION MEASUREMENTS

FIG. 1. Surface tension versus sodium dodecyl sulfate (SDS) concentration at the air/aqueous solution interface (closed symbols), diethylphthalate (DEP) and DEP-saturated aqueous solution interface (closed symbols), and silicone oil (350 cst . )/ aqueous solution interface (blue squares). Although the surface tension is lower for the silicone oil than the air interface, the interfacial tension gradient does not change significantly. We therefore assume that the transient DEP/pure aqueous solution interface at the front of the droplet exhibits a similar gradient of 0 01N/m/M stated in the manuscript, which is the dominant driving force for the droplets. As a function of [NaCl], fits to the von Szyszkowski . equation (dashed lines) reveal that the surface tension gradient is independent of salt concentration for a given interface. On the other hand, the CMC (intersection with the horizontal line) decreases with salt for both the aqueous and the DEP-saturated interface, but the difference ∆ CMC is progressively smaller. We thus use these estimates of ∆ CMC to fit the velocity decrease with [NaCl] in Fig. 3 in the manuscript.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00581v1-figure-10b1263d-55e1-4a61-901a-54e69636f00f.jpg)
Adrien Izzet Pepijn Moerman Jan Groenewold Jérôme Bibette Jasna Brujić
Tunable active rotational diffusion in swimming droplets
Here we characterize the motility of athermal swimming droplets within the
framework of active rotational diffusion. Just like active colloids, their
trajectories can be modeled with a constant velocity $V$ and a slow angular
diffusion, but the random changes in direction are not thermally driven.
Instead, $V$ is determined by the interfacial tension gradient along the
droplet surface, while local micellar fluctuations lead to droplet
reorientation with a persistence time $\tau$. We show that the origin of
locomotion is the difference in the critical micellar concentration $\Delta$CMC
in the front and the back of the droplet. Tuning this parameter by salt
controls $V$ from $3-15$ diameters $d/s$. Surfactant concentration has little
effect on speed, but leads to a dramatic decrease in $\tau$ over four orders of
magnitude. The corresponding range of the effective diffusion constant $D_e$
extends beyond the realm of synthetic or living swimmers, in which $V$ is
limited by fuel consumption and $\tau$ is set by temperature or biological
activity, respectively. Our tunable swimmers are ideal candidates for the study
of the departure from equilibrium to high levels of activity, on both the
single particle level and their collective behavior, including the
motility-induced phase separation (MIPS).