Arxiv:1907.04278V1 [Cond-Mat.Soft] 9 Jul 2019 As Schematically Shown in Fig

Nonlinear Electrophoresis of Highly Charged Nonpolarizable Particles Soichiro Tottori, Karolis Misiunas, and Ulrich F. Keyser∗ Cavendish Laboratory, Department of Physics, University of Cambridge, JJ Thomson Avenue, Cambridge, CB3 0HE, United Kingdom Douwe Jan Bonthuis Institute of Theoretical and Computational Physics, Graz University of Technology, 8010 Graz, Austria (Dated: July 10, 2019) Nonlinear field dependence of electrophoresis in high fields has been investigated theoretically, yet experimental studies have failed to reach consensus on the effect. In this work, we present a systematic study on the nonlinear electrophoresis of highly charged submicron particles in applied electric fields of up to several kV/cm. First, the particles are characterized in the low-field regime at different salt concentrations and the surface charge density is estimated. Subsequently, we use microfluidic channels and video tracking to systematically characterize the nonlinear response over a range of field strengths. Using velocity measurements on the single particle level, we prove that nonlinear effects are present at electric fields and surface charge densities that are accessible in practical conditions. Finally, we show that nonlinear behavior leads to unexpected particle trapping in channels. Electrophoresis has been widely used for sensing, filtra- trokinetics found for polarizable particles [20], is not well tion, manipulation, and the separation of molecules and understood for intermediate to large κa in contrast to the particles, particularly recently using micro- and nanoflu- case of small κa in nonpolar electrolytes [21–26]. The- idic devices [1–4]. On these length scales, the applied oretical analyses have predicted that the electrophoretic electric field may reach kV/cm because the confinement velocity becomes nonlinear to the applied field at a mod- ∼ focuses the electric field, even with a moderate applied erate to high field β = aE/φth & 1 [14, 27–29]. However, voltage to the system. The electrophoretic velocity vep a well-controlled experiment in this regime has been lack- of a rigid nonpolarizable sphere in the limit of low zeta ing. Indeed, there are experimental reports showing that potential ζ and low applied electric field E is particle electrophoretic velocity is linear [30] and nonlinear [27] under similar conditions. Even for the nonlinear 2 v = f(κa)ζE, (1) case, the properties of particles were not clearly char- ep 3η acterized. Moreover, the zeta potential was reported to where and η are the permittivity and viscosity of the be severalfold lower than the value used for fitting one fluid, respectively [5]. f(κa) is Henry’s function, where of the models, and the difference was then attributed to 1 κ− and a are the Debye length and particle radius, the divergence of the slip plane from the actual particle respectively [6, 7]. Within a canonical model (homo- surface [31, 32]. Here, we resolve the conflicting results geneous and η, no specific ion-surface interactions, found in the literature by systematic experiments in com- no hydrodynamic slip) and using a low zeta potential bination with detailed simulations. ( ζ < φth = kBT/e 26 mV), ζ calculated using Eq. (1) We used high-speed video tracking to study the non- corresponds| | to the≈ homogeneous electrostatic potential linear electrophoresis of highly charged particles in mi- at the particle surface. In real systems, however, the crofluidic channels. We characterized particle mobility zeta potentials of particles may exceed φth, in which case with different salt concentration, and subsequently mea- the direct action of the electric field on the double layer sured the nonlinear electrophoretic velocity under high as well as the advection of ions by the flow result in an field up to β 3. The experimental data were compared ≈ asymmetric shape of both coion and counterion clouds, with the values obtained with coupled Stokes-Poisson- arXiv:1907.04278v1 [cond-mat.soft] 9 Jul 2019 as schematically shown in Fig. 1(a) [8, 9]. For interme- Nernst-Planck (SPNP) simulations, showing that high diate values of κa, the double-layer distortion, known as field particle electrophoresis experiments can be quan- relaxation effect, leads to a nonlinear dependence of vep titatively explained with our SPNP model. We also on ζ [7, 8, 10]. Therefore, in case ζ & φth, the model in- demonstrate that particles can be trapped in the channel cluding relaxation effect [11, 12] needs| | to be used instead at high fields due to this nonlinear effect in combination of Eq. (1) to extract the surface potential that reproduces with the electro-osmotic flow resulting from the channel vep within the canonical model. wall. In addition, the relaxation effect is field dependent, Microfluidic channels were made out of polydimethyl- complicating the interpretation of experimental results siloxane (PDMS) cast onto microfabricated molds [33]. even further (see Supplemental Material [13]). This non- The mold for the “wide” channel (w = 12.7 µm, h = linear effect, which is different from the nonlinear elec- 5.2 µm, l = 100 µm) was fabricated with SU-8 photore- 2 (a) (b) Laser Doppler Microelectrophoresis (Zetasizer Nano ZSP, E Malvern). All particles were diluted in KCl solution Ag/AgCl electrode to the final concentration (volume fraction) of approx- ζ < ϕth + + + 5 5 4 + imately 2.5 10− , 1 10− , and 2.5–10 10− for + tube KCl solution × × × + + narrow and wide microchannel experiments and Laser + + PDMS Doppler Microelectrophoresis measurement, respectively. ζ ϕ E ≳ th + + In Fig. 1(c), the electrophoretic mobility of PS particles + PS particle - + - in 10 mM KCl is plotted as a function of pH. The pH was + + access channel PDMS-coated + glass + - objective adjusted by adding small quantities of KOH or HCl. The lens + + 8 2 1 1 observed negative mobility of 7 10− m V− s− is ≈ − × (c) -4 due the negative surface charge associated with the sul- ) -1 PS 620 nm, 10 mM KCl fate end groups on PS, originating from the decomposi- s -1 tion of the initiator used during particle synthesis (ini- V 2 + tiator: K2S2O8, end group: -OSO3−K ). The reported m -6 -8 values vary but typically are in the range of several tens 2 (10 of mC/m . As the dissociation constant pKa of sulfate ep μ end groups is below 2 [35–37], no significant change in -8 electrophoretic mobility was observed for pH 4. We 2 4 6 8 10 decided that no buffer was necessary in the following≥ ex- pH periments, and pH values were all kept in the range of 6–9. FIG. 1. (a) Schematic illustrations of the counterion cloud around particles with low and high zeta potentials under an Coupled SPNP simulations were performed using a fi- electric field E. (b) Schematic illustration of our experimental nite element method (COMSOL Multiphysics ver. 4.4). system. Microfluidic chip is connected to two open reservoirs. The particle was fixed in a large enough space (= 100a) Particle motion in channels is imaged using a camera via the with open boundary conditions. The surface charge den- objective. (c) Bulk electrophoretic mobility of PS particle sity was set at a constant value on the particle surface. (620 nm in diameter) as a function of pH. Dotted line is the The electro-osmotic flow generated far away from the averaged value of µep for pH 4. ≥ particle was used to calculate the electrophoretic mobility of the particle [5]. More details on the simulations can be found in Supplemental Material [13]. sist. The “narrow” channels (w = 750 nm, h = 750 nm, In Fig. 2, the electrophoretic mobilities µep of var- l = 10 µm) were made by depositing platinum onto silicon ious particles measured with Laser Doppler Micro- substrates with focused ion beam. The access channels electrophoresis are plotted as a function of KCl concen- for wide and narrow channels were made out of SU-8 tration c. The field used in the measurement was be- (h = 100 µm) and AZ 9260 (h = 12 µm), respectively. low 0.025 kV/cm (β 0.05); thus nonlinear effects were The bottom glass slides were coated with a thin layer of negligible. The electrophoretic≈ mobility curves of PS PDMS, then plasma bonded to the PDMS channels [34]. particles display negative peaks at about 10 mM, in a The overview of our experimental system is shown in good agreement to previous literature values for other Fig. 1(b). The channels were connected to the two reser- highly charged particles [35, 38–40]. The negative peak voirs through access channels and tubes. The reservoir of PMMA particles, which had lower charge density, was height was adjusted to eliminate the pressure difference found at around 1 mM. between both sides before each experiment. Voltage ( ≤ In order to estimate the surface charge density from 35 V) was applied through Ag/AgCl electrodes connected these data, we fitted the approximate model including the to a source measure unit (Keithley 2450). The particle relaxation effect by Chen and Keh [12] (valid for κa & 20) motion was recorded at maximum 2000 fps and tracked to the experimental data with a surface charge density as by custom-written software after experiments. Conver- the only fitting parameter (see Supplemental Material for sion coefficients k from voltage to field inside a channel, more details [13]). The extracted surface charge densities E = kV , were calibrated for each microfluidic chip de- of PS 620 nm, PS 370 nm, PMMA 520 nm, and PS-COOH sign by measuring the current I through the system, as 380 nm are 51.2, 41.7, 13.1, and 11.2 mC/m2, re- k = I/(ρAV ), where ρ and A are the conductivity of the spectively.− Note that− these− surface charge− densities are electrolyte and total channel cross sectional area, respec- not necessarily the bare surface charge densities because tively.

Load more