Chemokinesis-driven accumulation of active colloids in low-mobility regions of fuel gradients

Many motile cells exhibit migratory behaviors, such as chemotaxis (motion up or down a chemical gradient) or chemokinesis (dependence of speed on chemical concentration), which enable them to carry out vital functions including immune response, egg fertilization, and predator evasion. These have inspired researchers to develop self-propelled colloidal analogues to biological microswimmers, known as active colloids, that perform similar feats. Here, we study the behavior of half-platinum half-gold (Pt/Au) self-propelled rods in antiparallel gradients of hydrogen peroxide fuel and salt, which tend to increase and decrease the rods’ speed, respectively. Brownian Dynamics simulations, a Fokker–Planck theoretical model, and experiments demonstrate that, at steady state, the rods accumulate in low-speed (salt-rich, peroxide-poor) regions not because of chemotaxis, but because of chemokinesis. Chemokinesis is distinct from chemotaxis in that no directional sensing or reorientation capabilities are required. The agreement between simulations, model, and experiments bolsters the role of chemokinesis in this system. This work suggests a novel strategy of exploiting chemokinesis to effect accumulation of motile colloids in desired areas.

Chemotaxis vs. chemokinesis. Self-propulsion at the microscale is abundant and important in biology. Many eukaryotic and prokaryotic cells propel themselves in biological fluids to perform critical functions. Examples include sperm cells, which navigate through cervical mucus to fuse with and fertilize an egg 1 , and leukocytes, whose migration is essential for coordinated immune responses and often occurs over long distances 2,3 . Two common migratory behaviors exhibited by motile cells are taxis and kinesis. Taxis refers to the phenomenon in which the direction of an organism's motion is determined by the non-uniform distribution of a physical quantity (i.e., the cell moves up or down a gradient). In contrast, kinesis occurs when an organism's speed, either translational or rotational, depends on the spatial distribution of a quantity 2 . Unlike taxis, kinesis does not imply a direction and is generally random.
In the natural world, taxis manifests in many forms, including chemotaxis (migration up or down a chemical concentration gradient), aerotaxis (oxygen gradient), phototaxis (light intensity gradient), gravitaxis (gravitational potential gradient), and others [4][5][6] . Taxis drives many cell migration behaviors known to be essential for various physiological and pathological processes. For example, motile bacteria execute chemotaxis to find nutrients and evade predators 7 . Aerobic bacteria perform aerotaxis to find the oxygen they need to survive 8 . Leukocytes reach sites of infection by orienting toward higher concentrations of chemicals secreted at these sites 2,9 . To perform taxis, cells employ a coordinated series of sensing and signaling processes. For example, Escherichia coli bacteria possess transmembrane receptors that bind attractants (such as glucose) and repellents (such as phenol) and thereby detect spatial or temporal differences in chemical concentration. These differences are relayed via a signaling protein to the flagella, which execute runs (straight-line motions) and tumbles (i.e., orientation changes), adjusting the cell's direction according to its needs. In addition to spatial gradients, bacteria can also respond to temporal gradients. For example, Salmonella typhimurium bacteria change their tumbling frequency, effectively modifying direction, in response to a change from one uniform concentration of a chemoattractant to another 10 . This suggests that the detection of a spatial gradient by bacteria may involve detection of a temporal concentration variation, which the cell experiences by moving through a spatial concentration Self-propelled microparticles: vehicles for artificial chemotaxis? Since the early 2000s, extensive research efforts have been devoted to developing synthetic analogues to biological microswimmers. These selfpropelled particles, often called active colloids, typically range in size from 0.1 to 10 μm and have the unique ability to harvest energy from their surroundings and convert it to motion 21 . They accomplish this via several mechanisms including bubble propulsion 22 , conversion of external electric or magnetic fields 23,24 , ultrasound 25,26 , or incident light 27 into motion, and phoretic mechanisms 28 in which the particle moves as a result of self-generated gradients in temperature, chemical concentration, or electric potential. The first active colloids to be widely studied were half-platinum half-gold (Pt/Au) bimetallic rods (2 μm long, between 200 and 400 nm diameter) that self-propel with the Pt end forward in hydrogen peroxide (H 2 O 2 ) solutions 29 . Several mechanistic analyses 30,31 , including some from our group [32][33][34] , established that electrochemical reactions on the Pt and Au surfaces generate an electric field in the rod's vicinity, which exerts a propulsive force on the charged rods (commonly called self-electrophoresis 30,32,35,36 ). Figure 1 shows a schematic of a Pt/Au rod with the electrochemical reactions powering the motion.
The speed of Pt/Au rods increases with local H 2 O 2 concentration, and the rods thus exhibit a positive chemokinetic response to H 2 O 2 . Some studies report a linear dependence of speed on concentration while others predict nonlinear relationships of varying types 29,32,33,37,38 . The general consensus is that the dependence is linear at low to moderate concentrations, and slightly levels off at high concentrations because of saturation of available reaction sites on the Pt and Au surfaces 37,38 . Other active colloid designs also exhibit a chemokinetic response to their fuel (typically H 2 O 2 ), including bubble-propelled particles 39 and Pt/polystyrene or Pt/silica "Janus" particles 40,41 . At higher fuel concentrations, speeds are usually faster because the reactions responsible for motion proceed more rapidly, leading to stronger propulsive forces.
In contrast to H 2 O 2 , the addition of electrolytes tends to decrease the speed of Pt/Au rods. This was first reported in 2006 for sodium nitrate and lithium nitrate by Paxton et al. 42 , who found that at constant H 2 O 2 concentration, speed is roughly inversely proportional to solution conductivity. (H 2 O 2 is required for electrolyteinduced speed reductions to be present, since without H2O2 there would be no self-propelled motion to impede.) The only exception to this trend is electrolytes that contain silver, which tend to increase the speed of the rods 43 . solutions. Pt and Au catalyze the indicated electrochemical reactions, resulting in an electron current through the rod (from Pt to Au) and a corresponding ion current in the solution, and create dipolar charge density (ρ e ) and electric field (E) distributions around the rod. The electric field exerts a propulsive force on the negativelycharged rod, propelling it with the Pt end forward (left to right) 32,33 . The speed of Pt/Au rods increases as the H 2 O 2 concentration is increased, while the addition of a nonreacting electrolyte tends to decrease the selfpropelled speed 34,42 .
Electrolytes not containing silver cause significant reductions in rod speed even at modest concentration, such that millimolar concentrations are often enough to effectively eliminate self-propulsion 42 . To elucidate the physical mechanism for this phenomenon, we simulated the motion of a single rod in varying salt solutions 34 and found that the self-generated electric field magnitude decreases markedly upon addition of salt (and the concomitant increase in conductivity). This trend is robust and has been characterized for a wide range of electrolytes 43 .
Several groups have attempted to realize autonomous taxis in artificial active colloids 44 . Some studies [45][46][47][48] have reported chemotaxis in chemically-powered particles, effectively making two simultaneous claims: (1) the particles move by consuming a chemical fuel that influences their speed (chemokinesis); (2) the particles move up gradients of that same fuel (chemotaxis). As previously stated, chemotaxis and chemokinesis can coexist in some organisms 19 ; however, given the signaling and sensing required for chemotaxis in living cells 49 , the question naturally arises as to how synthetic active colloids could achieve "intentional" motion along gradients. The argument is often made that when an active colloid happens to move up a fuel gradient (i.e., in a direction of increasing concentration), it accelerates and therefore advances farther toward increasing concentration than if it were moving toward lower concentration. This scenario is plausible but does not address the scenario in which a particle encounters a local maximum in chemical concentration. When a chemokinetic active colloid reaches a concentration maximum, it will quickly disperse because of the maximal motility there. In contrast, chemotactic cells do not disperse upon reaching a chemoattractant source; rather, their directional motion ceases. For example, chemotactic bacteria seek out and accumulate at nutrient maxima 50 . Sperm cells accumulate in follicular fluid containing a chemoattractant 51 , and there is a strong correlation between this in vitro accumulation and egg fertilizability 52 . Neutrophils use chemotaxis to find and accumulate at wound sites 53 . A recent numerical study considered a self-propelled particle in a linear fuel gradient and assumed a linear dependence of speed on chemical concentration; the main conclusion was that, although initially the particle may move in either direction, it would ultimately end up in the region of lowest motility 55 . As noted by Popescu et al. in a systematic study of Janus swimmers in concentration gradients 54 , "chemokinesis alone cannot lead to chemotaxis. " Some theoretical analyses have been undertaken of phoretic active colloids to identify the conditions in which they could exhibit chemotaxis, positing a rich variety of single-particle and collective behaviors 56,57 . It is wellknown from the work of Golestanian and others that the response of a phoretic self-propelled particle depends on its phoretic mobility and surface activity 35 . Popescu et al. theoretically analyzed the response of a phoretic self-propelled particle to a concentration gradient and concluded that, while a particle's chemokinetic response to a fuel depends primarily on its average phoretic mobility (in the case of Pt/Au rods, this is the electrophoretic mobility, which depends on zeta potential 58 ), chemotaxis primarily depends on the mismatch in the phoretic mobility between the front and back halves of the particle 54 . The zeta potential of metallic nanorod particles in aqueous solutions has been measured to be consistently negative (attributed to preferential adsorption of anions) 59 , and is around -40 mV for Pt/Au rods 42 . Since the particle is conducting, our previous theoretical analyses have predicted that the zeta potential of metallic nanorods differs negligibly between the Pt and Au surfaces 33 . Thus, the theoretical basis for Pt/Au active colloids to exhibit chemotaxis of this type has not been definitively established.
Other demonstrations of artificial taxis are based on different physical mechanisms. Silica/carbon two-faced "Janus" microspheres, when immersed in a binary mixture of water and 2,6-lutidine and irradiated with 532-nm light, asymmetrically de-mix the solution surrounding the particle, establishing a lutidine concentration gradient that produces a diffusiophoretic propulsive force on the particle 60 . In a light intensity gradient, a torque on the particle arises that selectively orients particles toward low-intensity regions, exemplifying artificial negative phototaxis. In another example, droplets of a liquid crystal move by generating surfactant gradients, driving Marangoni flows that propel the droplets 61 . The droplets evinced chemotaxis by successfully navigating a microfluidic maze with the chemoattractant source at the exit. Other examples include asymmetric artificial liposomes made from copolymer mixtures, which undergo runs and tumbles akin to bacteria and are propelled by enzymatic reactions 62 ; these particles were shown to be promising for blood-brain barrier crossing, a crucial step in many biomedical applications. Finally, gravitaxis was demonstrated in asymmetrically-coated Janus microspheres 63 ; here, the mismatch between the densities of the polystyrene colloids and metallic coatings led to a torque orienting the particles along the gravitational direction with the metallic side down.
In contrast to taxis, a lesser-explored strategy is to exploit the chemokinetic response that most active colloids naturally exhibit to bring about their accumulation in desired areas. One recent example study 64 demonstrated polymeric self-propelled particles that swell in response to decreases in pH. In a pH gradient, particles accumulate in the most acidic regions because the PVP swells in response to the local pH decrease, leading to an increased drag profile and reduced swimming speed there. Given that many active colloids naturally exhibit chemokinesis, a systematic characterization of the behavior of active colloids in static concentration gradients is warranted.
In this paper, we show that Pt/Au rods undergo chemokinesis-driven accumulation in steady-state concentration gradients of hydrogen peroxide (H 2 O 2 ) and potassium chloride (KCl) salt. We use a microfluidic device to generate steady-state linear gradients of H 2 O 2 and KCl. The rods are initially distributed randomly in these gradients and exhibit a positive and negative orthokinetic response to H 2 O 2 and KCl, respectively. The distribution of rods reaches a pseudo-equilibrium and concentrates in regions with low H 2 O 2 and high KCl concentrations, where their effective diffusivity is minimized. The experiments show good agreement with Brownian dynamics simulations as well as a theoretical model based on a 1-D Fokker-Planck equation, which models the particles as a continuous substance undergoing enhanced diffusion with spatially-dependent diffusivity. The simulations, PDE model, and experiments point to a simple explanation that chemokinesis can lead to accumulation of active colloids in chemical gradients. This accumulation is distinct from chemotaxis in that it does not require sensing of, or a direct response to, temporal or spatial gradients in concentration.

Methods
The experiments were conducted in a microfluidic device fabricated from polycarbonate, glass, double-adhesive Mylar, and nitrocellulose membranes (see Fig. 2a). All layers are patterned using a CO 2 laser ablator (Universal Laser Systems, Scottsdale, AZ). The lower and upper surfaces of the microchannels are formed from a glass microscope slide and poly(methyl methacrylate) (PMMA) sheets, respectively. Nitrocellulose membranes (Whatman, Maidstone, UK), with average pore diameter of 0.4 μm, form the boundaries between the main channel and the side channels and are sealed to the upper and lower channel surfaces using 200-μm-thick double sided adhesive Mylar sheets cut in the same pattern as the membranes. The device is assembled using metal bolts through the acrylic upper surface and an additional PMMA superstructure below the glass lower surface. The center microfluidic channel is 15 mm long, 400 μm wide, and 450 μm deep. The outer channels are 600 μm wide and 450 μm deep. The assembled device is depicted schematically in Fig. 2a.
The rods are fabricated using a templated electrodeposition procedure 65 that is widely used for synthesis of self-propelled rod-shaped particles 29,31,32 . Briefly, platinum and gold are sequentially electrodeposited into the cylindrical 200-nm-diameter pores of an anodic aluminum oxide (AAO) membrane. After metal electrodeposition is complete, the AAO membrane is chemically etched, and the rods are centrifuged and resuspended in pure water. The length of the Pt and Au segments is each approximately 1 μm, yielding an overall rod length of 2 μm. The rods are imaged in the microfluidic device using optical microscopy with a 20 × objective (NA = 0.45, Nikon TE2000, Japan) and a CCD camera (Coolsnap HQ, Photometric, Tucson, AZ) at 2 frames per second. The motion of the rods is tracked and analyzed using a custom MATLAB-based particle tracking algorithm. In each frame, the positions of the rods' centers were calculated from the intensity-weighted centers of the rod images. The rods' centers at each time were paired using an optical flow algorithm. In a typical experiment, there are between  ; and (f) rods centered at x/w = 0.9 (D eff /D B = 6.2). (g) The measured effective diffusivity (blue circles) varies roughly linearly with x. The error bars represent one standard deviation of the measured diffusivities and the red solid line is a least-squares linear fit that is used in the models. Panels (c)-(g) originally appeared in a previous review by our group 28  www.nature.com/scientificreports/ 350 and 600 rods in the field of view. The width of the channel is divided into 20 segments. The mean squared displacement (MSD) for each rod is tracked for 100 frames. An effective diffusivity for each of the 20 segments is determined by averaging the MSDs for all rods starting in the respective segments at the start of the 100 frames. We drive flow of 30% H 2 O 2 in deionized water at 5 μL/min through the left channel, which serves as an H 2 O 2 source and KCl sink. We drive 100 μM KCl + 1 μM fluorescein dye (the dye is used to visualize the concentration gradient using optical microscopy) through the right channel, which serves as a KCl source and H 2 O 2 sink. The flows in the side channels are driven with syringe pumps. As described in the supplementary information, the maximum and minimum concentrations of KCl and fluorescein in the main channel are calculated as C min,KCl = 45.5 μM, and C max,KCl = 54.5 μM. C min,Fl = 455 nM, and C max,Fl = 545 nM. As shown in Fig. 2b,

Theory.
Complementing the experiments, we model chemokinesis using Brownian dynamics (BD) simulations as well as a mass conservation partial differential equation (PDE). The BD simulations are based on the standard Langevin model that neglects inertial effects, collisions between rods, and interparticle hydrodynamic and electrostatic interactions. We assume the motion of the rods is confined to two dimensions, since the negatively buoyant rods sediment to a quasi-2-D fluid layer a few micrometers above the bottom surface, where they remain for the duration of the experiment. Many of the materials from which active colloids are made are denser than water, and thus 2-D (or even 1-D) models of ensembles of active colloids are common 40,55,66-68 .
The differential equations governing the position x and y and orientation θ of each rod are the Langevin equations 55,66,67,69 , where ẋ and ẏ are (respectively) the translational speeds in the x and y directions, u is the translational speed of the rod (which depends on x), and ξ x , ξ y , and ξ θ are stochastic noise terms in the indicated directions. This basic system of Langevin equations has been used to model quasi-2-D swarms of chemokinetic active colloids in several previous studies 55, [66][67][68]70,71 . The right-hand side of the third Langevin Eq. (3) includes only the stochastic term, ξ θ , reflecting our assumption that the rods are not "self-polarizing 66 . " In other words, we assume that the orientation of the rods varies in time because of rotational diffusion alone, and that the rotational diffusivity is approximately independent of H 2 O 2 or KCl concentration. The assumption of no self-polarization has been made in several previous models of chemokinetic microswimmers 40,55,70,72 . For many phoretic self-propelled particles that show a strong chemokinetic response, the rotational diffusivity is primarily governed by rotational Brownian motion and does not depend strongly on stimulus 67,72,73 . For Janus active colloids whose front and back mobilities are nearly equal, the systematic analysis by Popescu et al. found that the reorienting torque, and thus the chemotactic response, are negligible 54 .
The BD simulations do not specify an explicit dependence of speed on either H 2 O 2 or KCl concentration. As discussed above, the dependence of speed on both concentrations has been measured by several groups, including ours 29,32,37,42,43 . Although the qualitative trends are well-established (speed increases approximately linearly with H2O2 concentration and with the reciprocal of KCl concentration), different studies report different values of the rod speed for the same H 2 O 2 concentration 29,32,37 . In this work, we found that even the same batch of rods showed different average speeds from one day to another under nominally identical conditions. In experiments, ensemble averages of rod particle tracking showed a roughly linear dependence of rod speed on horizontal position x 74 . As a result, in the BD simulations presented here, we assumed a linear variation of u with x, bounded by maximum and minimum speeds, V max and V min , which were specified inputs for each simulation run. The minimum value of V min we considered was 1 μm/s and the maximum value of V max was 25 μm/s. This is a similar range to the speeds measured in the experiments 74 .
The stochastic displacements ξ x , ξ y , and ξ θ are independent random variables with zero ensemble average, �ξ x � = �ξ y � = �ξ θ � = 0, and their standard deviations are tuned such that the mean-squared displacements (without self-propulsion) are consistent with Brownian motion with translational diffusivity D and rotational diffusivity D rot , where brackets indicate an ensemble average over the population of rods, D is the Brownian diffusion coefficient of the rods, and D rot is the Brownian rotational diffusivity of the rods. For the purposes of tuning the noise terms, the translational and rotational Brownian diffusivities are estimated from the Stokes-Einstein equations for translational and rotational diffusion of a spheroidal particle 75 . For the dimensions considered here, D ≈ 1 (2) y = usinθ + ξ y , www.nature.com/scientificreports/ μm 2 /s and D rot ≈ 1.38 rad 2 /s. The boundary conditions employed for the BD simulations enforce no flux of rods (i.e., the channel is closed and rods cannot escape) and elastic boundaries. That is, rods are assumed to undergo perfectly elastic collisions with the walls and the overall kinetic energy of the rod is the same before and after the collision. The simulations track 10,000 rods over 6,000 s with time step Δt = 0.1 s. The BD simulations and experiments are complemented with a theoretical model. This model assumes the swarm of rods to be a continuously distributed substance that obeys the mass conservation equation where n is the number density of rods and the flux of rods J is defined by the modified Fokker-Planck equation 76-78 , where D eff is the effective diffusivity of the rods (which depends, in turn, on the spatially-dependent rod speed), n is the number density of the rods, and α represents the Itô-Stratonovich convention coefficient. The Fokker-Planck equation has been used to model chemokinesis of active colloids in activation gradients in a variety of situations 55,66,67 . The supporting information (SI) provides a full derivation of Eq. (2). The appropriate value of α for a random-walk process is 0.5 as discussed by in detail by Schnitzer 76 . Equation (6) is solved in one dimension with no-flux boundary conditions, J = 0, i.e. D eff ∇n = −αn∇D eff . The only input into the models is the variation of the effective diffusivity of the rods with space, which is obtained from a least-squares linear fit (r 2 = 0.84) to the experimental data in Fig. 2g. We solve Eqs. (6) and (7) using a second-order centered-difference method in space and a 1 st -order backward Euler method in time. We assume a uniform distribution of rods for the initial condition and impose no-flux boundary conditions on the channel walls. Further details on the numerical method are provided in the SI.

Results and discussion
Since the rods' velocity increases with H 2 O 2 concentration 10,13,14,27 and decreases with electrolyte concentration 34,42 , the linear static gradients in H 2 O 2 and KCl concentrations result in a spatial variation in the translational speeds of the rods. Figure 2c shows the pathlines of individual rods in the center channel as a function of position, x/w, over 100 s. The rods on the left (x/w < 0.5, high H 2 O 2 and low KCl concentration) move faster, resulting in longer pathlines (green) relative to the ones on the right (x/w > 0.5, locally low H 2 O 2 and high KCl concentration) over the same time interval (red). We use the rod pathlines to calculate the MSD, L 2 , of the rods for the lifetime of the experiment.  Fig. 2d. In the channel center, x/w = 0.6, the effective diffusivity decreases to D eff /D B = 30.7 (Fig. 2e). On the right-hand side where the salt is highly concentrated and there is little peroxide, x/w = 0.9, the effective diffusivity reduces to roughly six times the rods' thermal diffusivity. Figure 2g shows that the effective diffusivity of the rods decreases roughly linearly with space, x/w. This gradient arises from the gradient in velocity of the individual rods. Figure 3a shows an instantaneous micrograph of the rods in the center channel at t = 0, showing a relatively uniform initial distribution of rods. We plot the number of rods as a function of position in Fig. 3c for the experiments, simulations, and model. Corresponding to Fig. 3a, at t = 0 the rods are uniformly dispersed throughout the channel. After 12 min, the rods accumulate on the right due to chemokinesis, as shown in Fig. 3b. At this time, the rod number distribution reaches a pseudo-equilibrium state in which the rods are concentrated in the high-salt, low peroxide region on the right as shown in Fig. 3b. The temporal evolution of chemokinesis in both simulations and experiments is included as a video in the supplementary information (SI). The SI also includes a description of control experiments that showed no migration of rods in a KCl gradient alone (no peroxide), effectively eliminating diffusiophoresis due to the KCl gradient as an alternative mechanism for the rods' accumulation. We plot both the BD and PDE model predictions for the number density as a function of space in Fig. 3c, which show good agreement with experimental results. We also note the qualitative resemblance between the quasi-steady-state distributions (solid lines) in Fig. 3c and the steady-state probability density function in Fig. 2a of Ghosh et al. 55 , who also found a roughly inverse relationship between speed and steady-state accumulation.
To quantify the asymmetry in the accumulation of rods at steady state, we define a chemotactic index, CI. Our use of the term "chemotactic index" should not be construed as a claim that the rods exhibit chemotaxis here; rather, we use this term because it has been used in previous literature 79,80 to denote accumulation of microswimmers in a given region. Indeed, in this system accumulation results not from chemotaxis, but from chemokinesis. The CI is defined as the number of rods in the high-salt, low-peroxide concentration region (right hand side) divided by the number of rods in the low-salt, high-peroxide concentration region (left hand side) 81 . Specifically, where N R and N L are the number of rods in the rightmost and leftmost regions of the main channel, respectively. The width of the rightmost and leftmost regions is defined to be 1/10 th of the channel width. Figure 3d shows the CI as a function of time. The CI is initially unity and increases as CI(t)/CI final = 1 − e −t/τ 1 + CI final −1 over (6) ∂n ∂t + ∇ · J = 0, www.nature.com/scientificreports/ τ ~ 4 min until it reaches a pseudo-equilibrium CI ~ 5. Cases with large diffusivity gradient, dD/dx, and small minimum diffusivity, D min , tend to attain higher equilibrium CI values. Although Fig. 3c shows good agreement among the experiments, simulations, and model after 12 min, Fig. 3d shows some disparity between the experimental estimates of CI and the BD and PDE predictions. Agreement between experiments and predictions may be improved by introducing a weak dependence of rotational diffusivity on position. Specifically, if D rot were to increase with increasing H 2 O 2 concentration and/or with decreasing KCl concentration, this would effectively increase the predicted steady-state CI. Rods in high-salt, low-H 2 O 2 regions may reorient less often (and thus be less likely to escape the low-motility region, where they accumulate) and rods in high-H 2 O 2 , low-salt regions would reorient more often (and thus escape more frequently to lowmotility regions). In the PDE model, this would correspond to a case in which 0.5 < α < 1 74 .
Some of the disparity between the experiments and model predictions may also be attributed to hydrodynamic and electrostatic interactions among rods, which were ignored in both the Fokker-Planck and BD models. In the low-speed (salt-rich, peroxide-poor) regions, rods accumulate over time and therefore come close to each other more often. More detailed computational models (beyond the scope of this work) could reveal that the non-spherical nature of the rods, combined with their finite size, could produce hydrodynamic effects in the form of additional self-propulsion velocity terms and additional torques. Additionally, each rod generates its own electric field to propel itself, and so the interactions among the self-generated fields of neighboring rods could also play a role. For simplicity we neglect these complications here and show that while self-polarization may be occurring to a limited extent, it is not necessary to explain the observed accumulation.
To demonstrate the time scale and repeatability of the chemokinesis-based accumulation, we conducted a separate experiment in which the directions of the H 2 O 2 and KCl gradients were reversed several times, and in response the rods accumulated on alternating sides of the channel. The results are depicted in Fig. 4 (note the horizontal x-axis and vertical time axis). The channel was subdivided into horizontal bins and the number of rods in each bin was measured at each time. Initially (t = 0), the KCl reservoir was at x = 0 and the H 2 O 2 reservoir at x = w, so that H 2 O 2 concentration increased from left to right and KCl concentration increased from right to left. This arrangement resulted in accumulation of the rods near x = 0, where the rods' motility was minimized, which is visible after 800 s. At t = 800 s, the gradients were switched such that the H 2 O 2 concentration increased from right to left and KCl concentration increased from left to right. Over a 20-min period, the particles underwent a net chemokinetic accumulation on the opposite side of the channel. At t = 3100 s and t = 4800 s, the gradients were switched again, resulting in the same chemokinetic accumulation response in regions of minimal motility. Red dashes indicate the times when the gradient directions were switched. This experiment demonstrates the repeatability of chemokinesis-induced accumulation and shows that the time scale for chemokinesis to occur is on the order of 1500-2000s.
One can use standard diffusion theory to estimate the timescale over which the chemokinesis-based accumulation occurs. For Fickian diffusion in one dimension with diffusivity D, the root-mean-square displacement (the square root of the MSD), L, scales with the square root of time t as Let L be the width of the channel, in our case 400 μm. As discussed, the minimum and maximum values of effective diffusivity in the experiments range from roughly D min = 6D B to D max = 120D B (see Fig. 2g), where D B = 1 μm 2 /s. Inserting numerical values, we estimate that accumulation should begin to become significant after approximately t = L 2 /2D max = 666 seconds (about 11 min). As shown in Fig. 3b,c, accumulation is quantifiably present after 12 min, in rough agreement with this prediction.
Using the PDE model, we compiled a phase map of the pseudo-equilibrium CI as a function of D min and dD/ dx. The theoretical phase map along with the experimental results is shown in Fig. 5a. The phase map exhibits three distinct regions: (1) at low diffusivity gradient magnitude, there is no net accumulation and the rods exhibit nearly uniform velocity/diffusivities over the entire space; (2) when the effective diffusivity of the rods dips below the Brownian diffusivity D min < D B , this represents the case where the particles aggregate (and stop moving) or reach a sticky boundary; (3) the physical regime where chemokinesis occurs, is that for which D min > D B and |dD/ dx|> 0. The experiments are performed for a variety of effective diffusivity gradients dD/dx, but all have a minimum effective diffusivity near D B . If we normalize the effective diffusivity gradient by the minimum diffusivity, we can collapse the entire phase map onto a single line as shown in Fig. 5b, which compares the predictions of the PDE model (solid line) to experimental measurements (squares) and BD simulations (circles). The experiments, model, and BD simulations show good agreement and indicate that the CI increases with increasing diffusivity gradient. This is expected because stronger gradients result in greater net drift of high-diffusivity particles into regions with low diffusivities.
Since the addition of electrolytes increases the viscosity of an aqueous solution, we sought to rule out viscosity as a potential contributor to the accumulation of nanorods in KCl-rich regions of the channel. The viscosity of aqueous KCl solutions was measured and tabulated by Grimes et al. 82 By interpolating this data (presented in the SI), we found that the viscosity of a KCl solution at the maximum concentration considered in our experiments (approximately 54.5 µM) exceeds that of pure water by 0.00014% (less than 1 part per thousand). Thus, we conclude that KCl-induced viscosity increases play a negligible role in the chemokinetic accumulation of rods described in this work.

Summary
We have shown that Pt/Au self-propelled rods accumulate in low-motility regions when immersed in static fuel gradients because of their chemokinetic response to H 2 O 2 and salt. This chemokinesis-induced accumulation is similar to that observed in various prokaryotic and eukaryotic cells 20 . We demonstrated good agreement between experiments and both Brownian Dynamics simulations and a conservation equation based on the Fokker-Planck Equation.
Chemokinesis-based accumulation has the potential to be leveraged for a variety of applications including environmental remediation, targeted drug delivery, self-healing materials, and more.
Although this study has focused on Pt/Au rods, it is important to note that our findings are expected to apply generally to a variety of self-propelled particles that exhibit an orthokinetic response to fuel concentration. The key components of this system necessary for chemokinesis-driven accumulation, namely orthokinetic response In each case, the rods preferentially accumulate in high-KCl, low-H2O2 regions (i.e., where the rods' speed is minimized) within 1000 s. www.nature.com/scientificreports/ to fuel and random motion trajectories, are observed in a wide variety of active colloid systems. In any active colloid system in which the particle trajectory varies in a stochastic manner, the particles exhibit an orthokinetic response, and the gradient is reasonably strong, one can expect chemokinesis-driven accumulation will be observed and the phase diagram in Fig. 5a will be applicable.
In the present work, we have assumed klinokinesis to be negligible. That is, we have assumed the rods' rotational velocity is independent of chemical concentration and is governed exclusively by rotational diffusion. Likewise, the simulations assumed pure orthokinesis, and the agreement between the simulations and experiments appears to validate this assumption. However, there are active colloid systems, such as Quincke rollers 83 or Pt/Au-based microgears 84 , in which the rotational speed clearly depends on stimulus intensity (e.g. electric field magnitude or local fuel concentration), and the collective dynamics of these systems are the subject of ongoing study 85 . A logical extension of our work would be to study an active colloid system that exhibits both orthokinesis and klinokinesis, and to elucidate the role of each in determining accumulation of particles in certain locations. Experiments are shown as colored squares, where the color corresponds to the chemotactic index. The chemotactic index increases with decreasing minimum diffusivity and increasing diffusivity gradient. Larger gradients result in a higher asymmetry in diffusion (higher on left) and smaller minimum diffusivities result in rods that tend to remain at the right end where salt concentrations are large. The white dashed line delineates a region where there is negligible accumulation (CI ≈ 1) because the gradients in diffusivity are too small to result in asymmetric diffusion. A dashed black line shows the rod's Brownian diffusivity. Large CI values are attained below this diffusivity because the rods effectively become spatially fixed due to their lack of motility. This corresponds to rods aggregating and losing motility or becoming stuck to a physical boundary. (b) For linear effective diffusivity gradients, the contour plot can be collapsed onto to a single line of CI as function of dD * /(dx * D * min ). Open squares show experimental data, while black dots show results from Brownian Dynamics simulations. For the BD simulation data, the abscissa is varied by changing the minimum and maximum speeds, which in general simultaneously varies the effective diffusivity gradient and the minimum effective diffusivity in the system.