Regulating Aggregation of Colloidal Particles in an Electro-Osmotic Micropump

Unrestricted particle transport through microfluidic channels is of paramount importance to a wide range of applications, including lab-on-a-chip devices. In this article, we study using video microscopy the electro-osmotic aggregation of colloidal particles at the opening of a micrometer-sized silica channel in presence of a salt gradient. Particle aggregation eventually leads to clogging of the channel, which may be undone by a time-adjusted reversal of the applied electric potential. We numerically model our system via the Stokes-Poisson-Nernst-Planck equations in a geometry that approximates the real sample. This allows us to identify the transport processes induced by the electric field and salt gradient and to provide evidence that a balance thereof leads to aggregation. We further demonstrate experimentally that a net flow of colloids through the channel may be achieved by applying a square-waveform electric potential with an appropriately tuned duty cycle. Our results serve to guide the design of microfluidic and nanofluidic pumps that allow for controlled particle transport and provide new insights for anti-fouling in ultra-filtration.


Introduction
Micropores and microporous structures are omnipresent building blocks of applied fluid mechanics and systems for industrial fluid processing 1 . Introducing colloidal suspensions into such structures introduces an extra level of complexity, but also is of high technological relevance in filtration, separation, or assisted concentration of colloids 2 . Examples of this can be found in fields ranging from oil recovery 3 to food processing 4 . The interaction between channel walls and the (often charged) colloidal particles (CPs) initiates a range of multi-scale processes 5 , complicating the rational design of fluid processing and lab-on-chip devices.
Electro-osmotic pumping (EOPing) is a common technique for circulating fluids in microfluidic channels [6][7][8][9] . It is the method of choice 5,10 for nanofluidic channels with a cross section smaller than 1 µm 2 . The reason for this is that in EOPs the average flow velocity is independent of the channel diameter 11,12 , unlike pressure-driven flow, where the average flow velocity decreases with the square of the channel diameter. EOP-based nanoporous membranes have been applied in various methods of analyti-plied potential for a short period of time; as much as an order of magnitude shorter than the period associated with clogging. This observation enabled us to formulate a procedure to avoid clogging through the application of an asymmetric square waveform potential with appropriately tuned duty cycle and subsequently demonstrate its efficacy. We showed this time-dependent potential can indeed result in a net CP flow through the channel. This insight, combined with our understanding of the underlying interplay between electric-field-and fluid-driven transport processes will prove valuable for the future design of non-clogging EOPs for use in, e.g., lab-on-a-chip devices, microfluidic nanoparticle manipulation, and industrial processes such as ultrafiltration.

Electrically Driven Colloid Transport
In this section, we place our work into context by briefly reviewing recent work on the electrically driven transport of CPs through micro-and nanofluidic channels. The two dominant effects are: (i) fluid flow through the channel, be that pressure driven or osmotic, and (ii) phoretic transport of CPs. A typical feature of transport is the dominance of one of these over the other, while a balance -opposing with comparable magnitude -of these two can lead to CP-aggregation and eventual clogging.
A common way of inducing flow through a charged nano-or micropore containing an electrolyte is to impose an electric field (E-field) over the pore 9 . The E-field acts on the charge excess in the double layer screening the charge on the surface of the pore. As a result, an electro-osmotic flow of the fluid is induced in the channel, i.e., an EOP is formed. Such a pump is made more effective by the dielectric contrast between the fluid and the (insulating) material of the pore, which causes the E-field to decay predominantly through the pore 17 . That is, the density of field lines that go through the material of the pore is negligible compared to that going through the pore opening. This in turn leads to a strong electro-osmotic pumping effect.
A change in dielectric properties and the associated gradients in the E-field can also be exploited to trap/separate particles. For example, microfabricated post, ridges, and hurdles have been used in this manner [18][19][20][21] . Zhu et al. similarly observed focusing and separation processes of microparticles and cells in spiral and serpentine microchannels [22][23][24][25] . Dielectric differences can further lead to induced-charge electro-osmosis (ICEO) -a nonlinear electro-osmotic slip velocity that occurs when an applied Efield acts on the ionic charge it induces around a polarizable (potentially uncharged) surface 26 . This can play a critical role near the sharp corners found at the opening of a microfluidic channel, where the change in electric field is significant. Squires and Bazant 26 argued that this results in fluid vortices and localized particle aggregation. Such vortices were first studied by Thamida and Chang 27,28 and later applied by Zehavi et al. to induce rapid particle accumulation and trapping 29,30 . Similar non-linear electrokinetic transport phenomena have also been reported by using isolated symmetrical sharp corners 31,32 or bipolar electrodes inside a straight microchannel [33][34][35] .
Other approaches to transporting particles through or trapping them in/near pores use linear electro-osmotic effects in combination with a secondary driving mechanism. For example, Cevheri and Yoda 36,37 have reported that CPs were concentrated at microchannel wall and form bands when the pressure driven flow opposes electro-osmotic flow. Lochab et al. 38 demonstrated that the cross-stream migration of colloids gives rise to varying spatial distributions within microchannels. They observed CPaggregation both in fluid bulk and near the wall of the microchannel by using the inertial effect 39 in microfluidics. Huang et al. 40 demonstrated that continuous pressure-driven flow and electrophoresis could lead to the separation and focusing of proteins and other charged molecules in a channel. Similarly, Stein et al. 41 revealed that DNA-aggregation may be achieved through the interplay of electro-phoresis, electro-osmosis, and the unique statistical properties of confined polymers. Lastly, Rempfer et al. reported the use of a microcapillary pipette to concentrate λ -phage DNA at the tip of, and its subsequent delivery into, the capillary using a combination of electro-osmotic flow, pressure-driven flow, and electro-phoresis 42,43 .
These linear electro-osmotic and electro-phoretic effects may be enhanced by imposing a salt gradient over the pore. For example, imposing a salt gradient led to a substantial improvement of the capture rate of DNA into a nanofabricated SiN pore [44][45][46] . A salt gradient can also induce ionic diffusio-osmosis and -phoresis, i.e., a tendency of charged objects to move in a salt gradient. Rabinowitz et al. 47 showed that similar to ICEO, this combination can lead to nanoscale fluid vortexes and non-linear electro-osmotic flow around a pipette tip. Recently, Shin et al. 48 used a combination of fluid flow and diffusio-phoresis to accumulate CPs at the outlet of a microchannel; the particles may be orders of magnitudes smaller than the pore size. These authors further showed that flow-induced particle clustering can cover the pore and lead to irreversible clogging when left unperturbed for a substantial time (∼ 20 minutes).
In the above examples, direct-current (DC) E-fields were predominantly used. However, alternating-current (AC) fields may be applied to create greater sensitivity or selectivity. For instance, Xuan and coworkers 49-55 demonstrated selective concentration of particles and control of their transport into a reservoirmicrochannel junction using biased alternating E-fields. In addition, reversal of an outwardly directed E-field can be used to concentrate DNA near a pipette tip, when combined with electro-osmotic flow, electro-phoresis, and dielectro-phoretic effects 56,57 . Lastly, AC fields are also suited to remove or prevent clogging. Harrison et al. 58 reported a high AC voltage inducedcharge electro-osmosis that can concentrate particles inside channel reservoirs and avoid channel clogging. Such clogging may also be overcome by using well-designed geometries, like pillar arrays 59 and diamond-shaped gradient nanopillar arrays 60,61 .
To summarize, there is a wealth of knowledge on the transport, concentration, and trapping of CPs in nano-and microfluidic channels and the use of geometry and AC fields to prevent clogging. Our work explores the combination of AC electric fields and imposed salt gradients, to contribute to the current level of understanding and we will come back to the various features, introduced here, in our discussion section. Therein, we will also argue that our numerical approach provides the essential insights to understand and qualitatively explain our experimental findings and we identify the innovative aspects of our research.

Experimental Results
We used two measurement configurations, vertical and horizontal, which are schematically depicted in Fig. 1. Both feature two liquid reservoirs containing DI (de-ionized) water on one side and a sodium-chloride (NaCl) solution on the opposite side. The reservoirs are connected via a silica glass capillary (10 mm in length) containing the microchannel of 5 µm in diameter (length to diameter ratio of 2×10 3 ). An electrical potential was applied at the two reservoirs using copper (Cu) electrodes, which resulted in an electrokinetic flow. We prepared for a typical measurement by filling the top (T) and left (L) reservoirs, respectively, with a solution of high salt concentration ranging from 0.25 mM to 10 mM NaCl. The respective bottom (B) and right (R) reservoirs were filled with a solution of lower salt concentration or DI water. The fluorescent CPs (200 nm diameter) were added to the liquid in both reservoirs before placement. Horizontal configuration to observe the motion of the CPs inside the channel, 'L' and 'R' labels denote left and right, respectively. The two reservoirs (green) at either end of silica capillary (blue) contained aqueous CP (yellow) solutions and were connected to a waveform generator via immersed copper electrodes. The illuminated area is shown in red. from this region fluorescent light coming from the CPs (yellow arrows) is collected through the microscope objective.
The vertical setup is shown in Fig. 1(a). Here, the capillary was held vertically connecting two volumes of liquids. The T reservoir consisted of a polymer tube containing 80 µL of aqueous solution. A droplet on a glass cover-slip acted as the B reservoir. We recorded the particles motion at the B opening of the microchannel using an inverted microscope in fluorescence mode. Fig. 1(b) depicts the horizontal measurement configuration. Here, the glass capillary was placed horizontally on top of the microscope objective inside the immersion oil, which is necessary to obtain an aberration free image of the particles inside the channel. The two reservoirs at either end of the capillary (L and R, respectively) contained 40 µL-droplets of the aqueous CP solution. The equal size of the droplets also minimizes the hydrostatic pressure difference between the two reservoirs. Because of the length of the channel, the flow caused by any residual pressure difference is orders of magnitudes slower than that caused by osmotic pressure and electro-osmotic forces. CP traces were recorded through the microscope objective focused on the channel inside the capillary. In Fig. 2, we present a typical measurement of the total fluorescence signal at the B opening of the vertical channel. Upon applying the potential, we observed the intensity of the fluorescence emission from CPs at the B opening to increase rapidly, before saturating. This dynamics corresponded to CP aggregation, as can be appreciated from the insets of Fig. 2, which contain fluorescence still images made during this process. The ESI † contains a movie that shows the associated dynamics of aggregation. Next, we examined if the CP aggregate could be removed by inverting the applied potential. We thus applied an asymmetric square-waveform potential to our system, with a period of 2.5 s, a duty cycle of 80%, and an amplitude of 30 V. Fig. 3(a) shows the fluorescence intensity, which is an indicator of the CP accumulation, around opening B over 4 periods of this applied potential. Fig. 3(b,c) depict representative images of the CP intensity at the B opening of the channel after clogging and clearing, re- spectively. The ESI † contains a movie that shows the clearance process resulting from the potential reversal.
From the fluorescence signal it is clear that the CPs started to aggregate when applying positive potential (the potential decreases from opening T to opening B) across the two reservoirs. The aggregation process took a few seconds, after which the fluorescence intensity remained roughly constant. When the potential was reversed, the CPs were quickly dispersed away from the bottom opening of the channel. Interestingly, the time of unclogging is, in all our observations, including a wide range of applied potentials and salt concentrations in the T reservoir, a factor of 10 smaller than the time of clogging. The corresponding measured time-constants are plotted in Fig. S6 of the ESI † .
This last observation proved key to keeping the channel from permanent clogging, while simultaneously obtaining a net transport of CPs, as we will describe later. In relation to this, it should also be noted that when the positive potential is applied for too long, we were not always able to open the pathway by reversing the EOP direction. Therefore, to avoid irreversible clogging, the applied potential must be reversed at an appropriate moment.
We studied the impact of CP aggregation on the ionic flow through the channel by measuring the electric current passing through it. In Fig. 4, we plot the ionic current across the channel while applying square-waveform potentials of varying amplitude. Note that in all cases, there appears to be a dip in the current when the positive voltage is just applied, which saturates to a constant value. Upon switching to negative voltages, the current reversed and its magnitude decreased over the time interval. From this observation, we conclude that the clogging is predominantly affecting the CP transport and has a limited effect on the current.

Measuring Net Particle Transport
In the previous section, we have shown that an appropriately tuned alternating potential could be used to aggregate and disaggregate CPs at the pore opening. However, this is in itself not sufficient to show that net particle transport can be obtained. In this section, we therefore consider the trajectories of CPs inside the channel. Obtaining these trajectories necessitated a switch to our (effectively equivalent) horizontal geometry, which is shown in Fig. 1(b). Here, the L, R openings are equivalent to the T, B openings of our vertical setup. The microscopy data is, however, taken at the mid-point of the channel, instead of the opening. The tracks were recorded after applying an asymmetric square-waveform electric potential with period of 1 s, duty cycle of 90%, and amplitude of 30 V for approximately 10 minutes. The salt concentration in the L reservoir was C L = 5 mM; the R reservoir was filled with DI water. The sketch to the right of the graph illustrates how a kymograph trace can be related to a particle path.
Following the no-clogging recipe proposed in the previous section, we applied a square-waveform periodic potential with period of 1 s, duty cycle of 90%, and amplitude of 30 V. The salt concentration in the L reservoir was C L = 5 mM, the R reservoir was filled with DI water. In this case, even after 10 minutes from the start of the measurement, no clogging was observed and the net motion of the CPs through the channel was preserved. This net flow is clearly present in the kymograph in Fig. 5 that depicts particle tracks in our field of view for one period of the AC field; the ESI † contains a movie showing the corresponding tracks in real time. Clearly, the particles can be transported through the channel. At the same time, we also measured a non-zero ionic current which similar to the measurements in the vertical configuration, see Fig. S5 of the ESI † .
We note in passing that by applying a DC potential for about 10 minutes, to verify clogging, we could also observe a different form of particle stagnation, akin to what has been recently reported by Shin et al. 48 . During this time a cluster of fluorescent CPs slowly grew into the imaging area in the middle of the channel, eventually filling the whole field of view. This clogging phenomenon starts from inside the channel due to a separate mechanism that is the subject of a follow-up study 62 . Once the channel had clogged in this manner, we were unable to open the pathway via EOP in the opposite direction, even at a greater applied voltage. Therefore, to avoid clogging in presence of CPs, applying a constant potential across the channel for times exceeding a few seconds has to be avoided at any moment during the experiment.

Numerical Modeling
We modeled our system numerically to understand the processes that cause CP clogging near the B opening. Here, we will briefly outline the model to aid our discussion, full details are provided in our follow-up study 62 . We solved the fully coupled Stokes-Poisson-Nernst-Planck (SPNP) equations for a stationary response using the COMSOL Multiphysics 5.4 suite that uses the finiteelement method. This means that for the ions: advection by the fluid, thermal diffusion, and migration in an electric field are accounted for. From the obtained fluid velocity, salt concentration profiles, and electrostatic potential, the velocity of the CPs was computed, which incorporated advection, electro-phoresis, and ionic diffusio-phoresis, respectively. Here, we made a pointparticle approximation, i.e., the presence of the particle does not significantly perturb its environment. This is a necessary reduction to keep our calculation tractable due to the wide separation of length scales involved. Fig. 6 Sketch of the axisymmetric (about the central red dashed line) microchannel geometry used in our numerical modeling. The inset shows a zoom-in of a part of the geometry that is of interest. The setup consists of a top and bottom hemisphere, radius R, with respective salt concentrations C T and C B and electrostatic potential V app (red) and ground (blue). These are separated by a divide and connected via cylindrical pore (grey), length L, inner diameter D, and a small rounding of the edges with radius r. The setup is filled with suspending fluid with relative permittivity ε f and the divide and pore are solids with relative permittivity ε p . The boundary color black indicates a no-slip, no-penetration, nosurface-charge boundary. On the boundaries marked yellow there is an osmotic surface slip velocity for the fluid.
We use the geometry depicted in Fig. 6 to perform these calculations with a set of parameters that best approximate the experimental setup. The sketch in Fig. 6 is not to scale. Specifically, the pore has an inner diameter of D = 2.5 µm, length of L = 10 mm, thickness of t = 50 µm, the corners are rounded with a radius of r = 500 nm, and the divide between the reservoirs has a width of 500 µm. The latter proved sufficient to ensure that most of the field decayed through the pore, as in the experiment. We ensured that we obtained mesh-independent results using the meshing ap-proaches laid out in refs. 42,43 , also see Fig. S7 of the ESI † .
An electrostatic potential V app and ground applied in the respective reservoirs, placed sufficiently far away from the pore to limit finite-size effects. On the walls, except for the inner channel walls, we use no-surface-charge boundary conditions and the fluid is treated as a dielectric medium with relative permittivity ε f . The divide and pore are solids with relative permittivity ε p . The fluid experiences a no-normal-stress condition at the reservoir edges and no-slip conditions on all the black surfaces of Fig. 6. On the yellow surfaces an effective slip velocity is imposed that captures ionic diffusio-osmosis and electro-osmosis. As a result, it is not necessary in our model to explicitly resolve the Debye screening layer. On the outer-edge of the pore the effective slip velocity smoothly transitions to the no-slip boundary condition on the divide, as indicated by the gradient from yellow to black. This improves the numerical stability of our SPNP solver.
Our system parameter choices are temperature T = 300 K, fluid dynamic viscosity η = 10 −3 Pa s and relative permittivity ε r = 80. The divide and pore have ε p = 4. For the ions we use Na + and Cl − with diffusion coefficients D + = 1.3 · 10 −9 m 2 s −1 and D − = 2.0 · 10 −9 m 2 s −1 , respectively. We varied V app from −100 to 100 V and take the zeta potential of the wall for the computation of the effective slip velocity to be −75 mV. The effective velocity of the CPs is computed using a zeta potential of −25 mV. The salt concentration in the B reservoir was fixed to C B = 0.5 mM, whilst the salt concentration in the T reservoir was set to C T = 7.5 mM. Fig. 7(a,b) shows the CP velocity around the bottom opening of the channel. Note that the speed is predominantly directed inward close to the pore axis, when the salt concentration in the T reservoir is C T = 7.5 mM. Right at the periphery of the pore opening, however, the particle flow is directed outward, see Fig. 7(c). The particle trajectories around the pore mouth are part of a larger vortex, not shown here, a vortex that is also present in the fluid velocity field, see Fig. S8 of the ESI † . Fig. 7b and the zoom-in in Fig. 7d shows that when there is no salt gradient and the E-field is reversed, the vortex disappears. This results in a reversal of the direction of CP motion near the wall compared to the center line, which is the origin of the observed vortexes around the pore. These plots further show that the salt gradient has almost no direct influence on the motion of the CPs. However, as can be appreciated by contrasting this result with Fig. 7(f), the presence of a salt gradient strongly modulates the way the E-field is shaped and decays throughout the pore, thereby having a dominant, albeit indirect effect. N.B. If only the E-field had been flipped in Fig. 7(f), graph (e) would be speed-wise inverted.  Fig. 6. The solid black curves shows the total speed, the red dashed curves show the advective contribution, the blue dashed curves show the electro-phoretic component, and green curve shows the effect of ionic-diffusio-phoresis. The inset to (e) shows that the latter has a small contribution; there is no contribution in (f). All scale bars are 2 µm.

Discussion
The combined experimental and numerical results lead us to conclude that CP aggregation and clogging of the pore can occur via a combination of advection, driven by electro-osmotic flow, and electro-phoresis. Both of these effects are strongly modulated by the presence of a salt gradient, which was found to be a necessary ingredient to the clogging of the pore opening. The second condition is that the applied potential is upstream with the salt concentration gradient. In this case, clogging of the channel occurs at the low-salt-concentration side.
The clogging of the opening of the channel was further found to only strongly affect the transport of CPs, as no appreciable change of ionic current was observed after CP aggregation. To explain this observation, we estimate the Dukhin length based on the definition l Du = σ /(eC), with σ the surface charge density, e the elementary charge, and C the bulk (monovalent) ion concentration 63 . Assuming a salt concentration of 0.1 mM for the DI water (an upper estimate) and a surface charge density of σ = 0.3e/nm 2 (typical number for Silica glass), the Dukhin length at the B reservoir is roughly 5 µm and hence comparable with the channel diameter. This estimation indicates that the surface conductance plays a major role in the measured ionic current, while the CPs mainly influence the bulk conductivity.
The slight change in electrical conductivity of the channel under negative bias (see Fig. 4) is followed by an equally long drop in the current signal after reversing the bias. This non-stationary behavior can be explained by considering the plug flow of liquid from the low concentration reservoir into the opening B of the channel, which is subsequently pushed out upon reversing the bias. The lower salt concentration in the B reservoir temporarily decreases the electric conductance of the channel, hence the decrease in the measured current. The duration of this negative bias is, however, not long enough for reaching an equilibrium with the higher salt concentration inside the channel. Therefore, under positive bias, the same portion of low-salt concentration liquid is pushed out of the channel, causing a short dip in conductance at beginning of each cycle. This is the reason for providing results with different 'effective' salt gradients in our Numerical Modeling section, as will become clear shortly.
Our finite-element based calculations revealed that under the application of a salt gradient, vortex like motion of the particles occurs close to the opening. This effect is reminiscent of the results by Rabinowitz et al. 47 , as well as the type of motion that emerges due to ICEO 26 . It suggests that the initial CP aggregation at the pore opening may be attributed to this vortex-like motion, as depicted in Fig. 7(e). A notion that is supported by a lack of aggregation when there is no salt gradient in the experiment, combined with the absence of a balance in the two driving components in the numerical calculations. Furthermore, reversing the bias temporarily reduces the salt gradient at the opening of the channel by pulling in the low-salinity solution. As a result, the vortex-like driving force on the CPs is replaced by another flow pattern that contain no nodes, similar to the pattern expected from stationary results for no salt gradient, which is depicted in Fig. 7(f). This outward flow can result in clearing all the CPs aggregated at the opening.
The main outcome of our investigation is that the flow pattern at the opening of the channel strongly depends on the saltgradient, as identified by our experiments and verified by numerical modelling. Since reaching the stationary state can take several minutes, the exact calculation of influx and outflux of CPs will require time-dependent modeling of the flow pattern. However, the separation of length scales relevant to our problem proved prohibitive, in spite of the significant computational power at our dis-posal. Such modeling will therefore be left for future work. Nevertheless, from these stationary calculations we can qualitatively understand the competition between different transport processes on CPs, which can direct us in choosing proper experimental parameters.

Summary and Outlook
Summarizing, we have experimentally and numerically studied the electro-osmotic aggregation of nanoparticles at the opening of a micrometer-sized silica channel in the presence of a salt gradient. Under the application of a DC electric field, aggregation of CPs at the pore opening occurred, which led to clogging of the channel. Based on numerical modelling, we explain the formation of the aggregate in terms of a balance between advective and phoretic forces. The clogging may be undone for sufficiently small aggregates through reversal of the potential. This insight allowed us to formulate a recipe for continuous nanoparticle transport by applying a square-waveform electric potential with an appropriately tuned duty cycle. Such a practical recipe is relevant to applications requiring continuous EOP, e.g., sensing using membrane nanopores.
A full understanding of the time-dependence of this process is still lacking and will be the focus of future efforts. However, considering the wide use of EOP in micro/nanofluidic devices, we are confident that our present results will already provide a new handle on the design and operation of such devices. The description we have provided for the underlying physical processes can further be useful for engineering situations where controlled aggregation can be advantageous, for example, for concentrating proteins in a well-defined location or for inducing crystallization by the long-period aggregation of (nano)particles.

Preparation of the Channel and Fluorescent Colloids
Silica capillary (microchannel): The Silica glass capillary (Fused Silica Capillary, CM Scientific, Silsden, UK) used in our experiments had an inner-channel diameter of 5 µm and an outer diameter of 150 µm. It was coated with a 12 µm thick protective layer on its outer surface. In order to receive the fluorescent light from our CPs inside the channel, we removed the opaque protective layer and wiped the capillary surface using a tissue wet with ethanol and isopropanol before constructing the sample cell. Finally, we cut the capillary to 10 mm using a diamond cutter.

Preparation of the sample:
The vertical configuration We put the prepared capillary into the top reservoir which had a 0.5 mm diameter hole in it. We subsequently added optical epoxy (NOA61, Thorlabs) in the gap between the capillary and the hole, and used a UV-light source (CS2010, Thorlabs) cure the glue, thereby fixing the capillary to the top reservoir. The horizontal configuration We put the prepared capillary in the middle of the plastic microscope slide which contains a rectangular opening of 1 × 2 mm 2 area. We used the optical epoxy (Griffon combi fast, Bolton Adhesives) between capillary and slide at the edge of the hole to fixed the capillary and to prevent the solution in the reservoirs from mixing with the immersion oil on top of the objective.

Fluorescent Colloids:
The fluorescent CPs (sicastar-greenF, micromod Partikeltechnologie GmbH) were produced by hydrolysis of orthosilicates and related compounds. They have a hydrophilic surface with terminal Si-OH-groups. The CPs were monodisperse with a mean diameter of 200 nm and nonporous with with a density of 50 mg/ml. They emit florescent light with wavelength of 510 nm when they are excited by a laser with a wavelength of 485 nm. The concentration of the CPs was controlled by diluting the stock solution with DI water. In our experiments, the CPs were suspended in different salt solutions, prepared by diluting a stock solution of 100 mM NaCl (≥ 99.5%, Sigma-Aldrich) to their final salt concentrations.

Microscopy Measurement Setup
Light from a solid-state single-mode laser (Coherent OBIS, 488 nm) was used to illuminate the observation area. In the vertical setup, we observed the CPs around the bottom opening, at a distance of roughly 5 µm from the end of the channel. This enabled us to study the channel cross section. In the horizontal setup, we instead covered in the middle of the channel in a direction along the axis, allowing us to see CP traces. In both cases, we used a home-built microscope with a 100 × 1.30 NA Olympus oil-immersion objective, 100 × overall magnification, and an effective field of view of around 200 ×200 µm. The fluorescent light from the diffusing particles was recorded using a sCMOS camera (Hamamatsu Orca Flash4).

Measurement preparation procedures
A small drop of solution without CPs was placed in the top reservoir for the vertical configuration setup or on the one end of the capillary for the horizontal setup. We then waited for 2 to 3 minutes until the liquid was pulled into the whole channel by capillary forces. After this, another drop of solution was placed on the other end of the capillary. An electrical potential with a maximum amplitude of 100 V was applied to our setup using Cu electrodes connected in series to a commercial current amplifier (DLPCA-200, FEMTO), an oscilloscope, and a data acquisition card (NI 6216-USB, National Instruments). The combination of these allowed us to simultaneously view and record ionic currents. By monitoring the ionic currents, we could confirm that the two reservoirs at either side of the channel were connected through the capillary. After establishing the electric connection, the solution in the top and bottom reservoir was replaced with drops of suspension containing the CPs. Square-waveform potentials in the range −100 V to +100 V with a range of duty cycles were applied across the channel, whilst the particle motion was recorded through the microscope.

Conflicts of interest
There are no conflicts to declare.