Frequency Response of Induced-Charge Electrophoretic Metallic Janus Particles

Abstract

The ability to manipulate and control active microparticles is essential for designing microrobots for applications. This paper describes the use of electric and magnetic fields to control the direction and speed of induced-charge electrophoresis (ICEP) driven metallic Janus microrobots. A direct current (DC) magnetic field applied in the direction perpendicular to the electric field maintains the linear movement of particles in a 2D plane. Phoretic force spectroscopy (PFS), a phase-sensitive detection method to detect the motions of phoretic particles, is used to characterize the frequency-dependent phoretic mobility and drag coefficient of the phoretic force. When the electric field is scanned over a frequency range of 1 kHz–1 MHz, the Janus particles exhibit an ICEP direction reversal at a crossover frequency at ~30 kH., Below this crossover frequency, the particle moves in a direction towards the dielectric side of the particle, and above this frequency, the particle moves towards the metallic side. The ICEP phoretic drag coefficient measured by PFS is found to be similar to that of the Stokes drag. Further investigation is required to study microscopic interpretations of the frequency at which ICEP mobility switched signs and the reason why the magnitudes of the forward and reversed modes of ICEP are so different.

1. Introduction

Microrobots are considered as a potential future workforce. Some examples of such applications include the use of active Janus particles to enhance optical resolution for measurements of molecular interactions in biological samples [1], enabling the optimization of transportation or navigation [2], self-assembly and formation of microscopic smart materials [3,4], serving as cargo movers for medicine delivery [5], and functioning as micromanipulators or micromixers [6]. Today, robotic devices are made at increasing smaller scales, reaching that of colloidal particles [5]. Among such endeavors are creative efforts to make active colloidal particles that convert energy provided by an external source to kinetic energy in order to move persistently [7,8]. Such active colloids can be powered by a variety of different mechanisms, such as self-thermophoresis [9], chemical decomposition [10], electric fields [11], and magnetic fields [12]. In this work, we apply an AC electric field and use induced-charge-electrophoresis (ICEP) to drive particle motion and address the possibilities of controlling such colloid-based microrobots by their electric frequency response.

Control of active colloids by external fields has been demonstrated in terms of precise alignment [13], directed control of micro motors [6], and cargo transport [3] for potential microrobot applications. Gangwal et al. [14] used a uniform distributed alternating current (AC) electric field in the kHz frequency range to produce a propulsive metal-dielectric Janus particle via ICEP. A theory for ICEP was first introduced by Squires and Bazant [15]. Based on the theory, an AC electric field could induce different electroosmotic flows around the dielectric surface and around the metallic surface, causing the particle to move. According to the model, the flow around the metallic surface was stronger than the flow around the dielectric surface; the pushing force on the metal surface was higher than the force on the dielectric surface. The net force pushes the particle in the direction from the metal side toward the dielectric side (normal direction), shown in Figure 1.

According to the model by Squires and Bazant [15], the ICEP mobility should decay to zero as the frequency of the AC electric field exceeds the characteristic times for the ions in the electric double layers to respond to the applied electric field. In Mano et al.’s experiments [16], however, a phoretic motion reversal was observed during a frequency scan, i.e., the Janus particle was found to be moving in a direction led by the metallic side. One possible explanation for such a mobility reversal was the reversal of induced-charge electroosmosis (ICEO) flow [15], i.e., the external electric field-induced flows of liquid around the two sides of the particles changed direction [17]. In an alternating current electrophoresis (ACEO) setting [18], the additional nonlinear mechanisms including Faradic reactions [19,20] and steric effects [17] due to ion crowding within the electrical double layers (EDLs), were considered to explain the phenomena. In traveling wave electroosmosis (TWEOF), the motility reversal could be expressed as the electrical body force [21,22]. All these explanations neither can be applied to the Janus particle system [21,22] nor need unreasonable ionic density [17,19,20].

A better understanding of the ICEP motility reversal requires better experiments that can determine the crossover frequency and the phoretic drag coefficient accurately. Using traditional DC measurements to determine the crossover frequency is difficult because the Brownian motion often swamps measurements of the null phoretic motion near the crossover frequency [23]. It is impossible to determine the drag coefficient of any phoretic motion by DC measurements, because the phoretic motion, absent of acceleration, is subject to zero net force, as the drag force perfectly opposes the phoretic force [24].

It is, however, possible to determine the drag coefficient by an AC measurement from the phase delay of the particle’s phoretic motion relative to that of the harmonically varying electrophoretic force, as demonstrated in earlier work [24]. We have used such a detection method to measure the crossover frequency and phoretic force of dielectrophoresis with high precision [23,25]. Briefly, our detection method, phoretic force spectroscopy (PFS), measures the harmonic response of an optically trapped phoretic particle in an amplitude-modulated AC electric field. A lock-in amplifier analyzes the motion of the particle to yield the particle motion’s amplitudes, and phase lags relative to the phase of the harmonic field. The crossover frequency is the frequency at which the phase lag changes by 180 degrees. In this paper, we describe how we use PFS to study the frequency dependence and phoretic drag coefficient of ICEP.

2. Materials and Methods 2.1. Fabrication of Metal-Dielectric Janus Particles

We made metallic Janus particles by depositing a thin film of metal on a monolayer of dry silica particles (SS05N, Bangs Laboratories, Inc., Fishers, IN, USA) [26] with a process shown in Figure 2a. A droplet of 3 micron diameter silica particle suspension was deposited on a glass substrate to create a monolayer by a vibration-assisted convective deposition method. The slide covered with a monolayer of particles was loaded into an e-beam evaporator (Indel E-beam Evaporator, International Delta Systems, LLC, Tucson, AZ, USA) to receive first a coating of 50 nm Ni and then a coating of 10 nm Au on the top half of the particle surface. The slide was then placed in a water-filled centrifuge tube and sonicated in a sonication bath (Branson 1510, Branson Ultrasonics Co., Danbury, CT, USA) for 6 h to release the particles into the water for further treatment. Schematics of the Janus particles and a scanning electron microscope (SEM) (JEOL JXA-8900 EPMA Microprobe, JEOL USA, Inc., Pleasanton, CA, USA) micrograph of the Janus particles are shown in Figure 2b.

2.2. Application of AC Electric Field to Drive Janus Particles Based on ICEP

Figure 2c shows the schematic of the sample chamber in which the particle’s motions were examined. A pair of indium-tin-oxide (ITO) coated glass slides, separated by a 50 um spacer of polycarbonate tape (3M 980, 3M, St. Paul, MN, USA), were fixed together by wrapping with the same tape. Copper tapes were attached to the ITO coating on each glass slide to form an electrical contact. An AC electric signal was applied across the pair of ITO slides to create a uniform electric field in the vertical direction. In the presence of the AC electric field, the Janus particles were driven into drift motions in the horizontal direction by ICEP. Since Janus particles settled very close to the bottom of the sample chamber, the particles moved in a two dimensional (2D) plane.

2.3. Application of a Magnetic Field to Fix the Direction of the ICEP Driven Phoretic Motion

While the vertically applied electric field exerted a torque on the Janus particles and aligned the metal cap edge parallel to the E field, the polar axis of the particles still had the freedom to orient freely in the plane perpendicular to the E field. Since the metal cap of the Janus particles exhibited magnetic dipole moment in the polar direction, an external magnetic field was applied to align the Janus particle orientation in the 2D plane perpendicular to the E field. The ICEP force was aligned with the polar direction of the Janus particles, and a magnetic field was used to control the direction of particle movement, as shown by Lin et al. [27] and Han et al. [28] to control and manipulate the metallic Janus particles, respectively. In this work, the magnetic field was generated by four 15 mm × 6.5 mm × 3 mm rectangular commercial Neodymium magnets (rectangular magnets, Theodora LLC, Seattle, WA, USA). The four magnets were divided into two groups, with two magnets in each group. Two groups of magnets were placed 10 cm apart and with the S and N poles facing each other in the horizontal plane. The sample chamber was placed at the center of the two groups of magnets. The intensity of the magnetic field at the location of the sample chamber was measured by Gauss meter (SJ200, Guilin Senjie Technology Co., Ltd., Guilin, China) to be1.0±0.1mT.

2.4. Position Detection of Phoretic Particles in 2D Using Image Analysis

The particle’s motion was measured by video imaging using a 20× objective lens on an inverted microscope (Olympus IX-81, Olympus Corporation, Tokyo, Japan). The videos of the particle movement were recorded by a charge-coupled device (CCD) camera (Sony Sscm256, Sony Corp., Tokyo, Japan) at a frame rate of 1 frame per second for a typical duration of a few minutes. The video recordings were analyzed using a particle tracking program (MOSAIC Suit on ImageJ, Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG), Dresden, Germany.) [29] to extract the trajectories of particles in the experimental chamber. Mean-square-displacement (MSD), defined as<|x^(t)2−x^(t+_t0)2|>, was calculated by taking the average of the MSD for all particles in the field of view. The experimental MSD vs. time for 3 µm metallic Janus particles at varying applied voltage over a fixed gap distance (50µm ) between the electrodes is shown in Figure 3.

2.5. Trapping and Manipulation of a Janus Particle

A 1064 nm wavelength laser (Nd:YVO₄ 1064nm diode-pumped solid-state CW laser, Spectra-Physics, Santa Clara, CA, USA) was used to create an optical trap, as shown in Figure 4a. The trap confines the Janus particle in a quadratic potential well [30]. The displacement of the particle in a trap, too small for video analysis to pick up, was tracked by a tracking beam (980 nm CL-2000 diode-pumped crystal laser, CrystaLaser, Reno, NV, USA) aligned collinearly with the 1064 nm laser. The tracking beam was received by a quadrant photodiode (QPD) that measured the position of the particle to give 5 nm spatial resolution at 1 kHz sampling rate. We used this setup to measure coefficients of the Stokes’ drag [31] An AC voltage was applied using a function generator (Stanford Research DS345, Stanford Research Systems, Sunnyvale, CS, USA) across the electrodes in the experimental chamber. The electric field produced by the electric potential across the electrodes drives ICEP and makes the particles an active Brownian particle. Active driving, trapping, and tracking of a metal-dielectric Janus particle was achieved in the same setup. A magnetic field was added to control the ICEP motion, so the particle moved linearly.

2.6. Phoretic Force Spectroscopy

Phoretic force spectroscopy (PFS) was used to measure the phoretic force of the Janus particles in the electric field [23]. Over a range of frequencies of the ICEP driving field, the motion of a particle in an optical trap was analyzed with an amplitude-modulated AC electric field. By using a lock-in amplifier, we measured the amplitudes and phase lags of an oscillated particle with the modulation signal as a reference. Not only the amplitude of the motion but also the phase relative to the driving force were obtained. The directional change of ICEP was detected when a 180-degree phase change was observed.

The lock-in amplifier (Stanford Research SR830 DSP Dual Phase Lock-In Amplifier, Stanford ResearchSystems, Sunnyvale, CA, USA) required a signal input and a reference input, as shown in Figure 4b. The reference was from the modulation frequency output of the function generator (Stanford Research, DS345, Stanford Research Systems, Sunnyvale, CA, USA). The signal was from the quadrant photodiode (homemade QPD) which measured the position of the particle. Here, the signal was amplified by a sensing amplifier (On-Trak OT301, On-Trak Photonics, Irvine, CA, USA) before it went into the lock-in amplifier. The amplitude (A) and the phase lag (δ) of the particle motion were measured by the lock-in amplifier.

The motion of the particle was analyzed as follows. Janus particles are assumed to have a constant speed,v, in an AC electric field with amplitudeE. The speed is proportional to the square of the AC electric field strength,v=β^E2, whereβis ICEP mobility. In the water, the drag force equals the ICEP forces,_Fe=_ζICEPv=_ζICEPβ^E2, where_ζICEPis the drag coefficient. At any instance, tweezer force, drag force, and resistance are in balance, so the Langevin equation of the position, x, of the particle is:

dxdt+_kotx_ζICEP=β^E2

where_kotis the spring constant of the optical trap. In our experiment, the electric field with 30% depth modulation wasE=_E0[0.85+0.15cos(_ΩMt)]cos(_ΩBt), where the_ΩMis amplitude modulation frequency. Equation (1) can be rewritten to be

dxdt+_kotx_ζICEP=β_E02 ^{[0.85+0.15cos(_ΩMt)]2}cos^(_ΩBt)2=β_E02{0.832+0.255cos(_ΩMt)+0.01125cos(2_ΩMt)}(12+12cos(2_ΩBt))=β_E02{0.416+0.1775cos(_ΩMt)+0.005625cos(2_ΩMt)+0.416cos(2_ΩBt)+12[0.1775cos((2_ΩB−_ΩM)t)+0.1775cos((2_ΩB+_ΩM)t)]+12[0.005625cos(2(_ΩB−_ΩM)t)+0.005625cos(2(_ΩB+_ΩM)t)]}

The amplitudes of the high-frequency terms (higher than_ΩB ) are small due to water damping [32]. Thus, we can ignore the higher frequency term and only consider the first three terms (in boldface in Equation (2)). The steady-state solution for Equation (2),x(t)has a term of a DC offset, a first harmonic term and a second harmonic term, each with amplitudeA(_ΩM), andA′(2_ΩM)and phase delayδ(_ΩM)andδ′(2_ΩM), respectively, i.e.,

x(t)=_{xDC offset}+A(_ΩM)^{ei(_ΩMt−δ)}+^A′(2_ΩM)^{ei(2_ΩMt−δ′)}

The frequency-dependent amplitudes A, A′ and the phase delaysδ,δ′are

A(Ω)=0.255β_{E02 ζICEPkot2}+^{(_{ζICEP ΩM})2}

^A′(2Ω)=0.01125β_{E02 ζICEPkot2}+^{(2_{ζICEP ΩM})2}

δ(Ω)=^tan−1_{ζICEP ΩMkot}

^δ′(2Ω)=^tan−12_{ζICEP ΩMkot}

Since the second harmonic (2_ΩM) terms have a much smaller amplitude than that of the fundamental frequency terms we use the_ΩMterm for further calculation. The two parameters of interest of this study, the mobilityβand the phoretic coefficient_ζICEP, can be determined in terms of the experimental measurables,A(_ΩM),δ(_ΩM)and_kot. According to Equations (4) and (6) above, we have the following relationships:

β=A(_ΩM)_ΩM0.255_E02sinδ(_ΩM)

_ζICEP=_kotΩMtanδ(_ΩM)

Although theA(_ΩM)andδ(_ΩM)are a function of modulation frequency, according to our model, theβ and_ζICEPshould be independent of the modulation frequency. Theβ and_ζICEPwere tested to be the same with modulation frequency from 1 Hz to 3 Hz. However, theA(_ΩM)andδ(_ΩM)decreased to about zero when_ΩMwas larger than 3 Hz. Thus, we set the modulation frequency_ΩMto be 1 Hz in our experiment for the best signal and use a lock-in amplifier to measureA(_ΩM)andδ(_ΩM)from which we calculated the mobilityβand the phoretic coefficient_ζICEPof ICEP-driven metallic Janus particles

The force diagram of a metallic Janus particle under the modulated AC electric field was shown in Figure 4b. The optical force,_Fot, was generated from the optical trap. The ICEP force was generated from the AC electric field. Since the electric field and magnetic field were applied to fix the orientation of the Janus particle, the direction of the ICEP force was fixed. The average ICEP force was non-zero due to the 30% amplitude modulation. The particle oscillated around an offset that was predicted by Equation (3). The 30% amplitude modulation generated an offset which became the balance point between the force from the trap and the force from ICEP.

We conducted an independent experiment (passive microrheology [33]) to measure the trap stiffness_kotusing the same optically trapped particle without the presence of the AC electric field. According to the equal partition relationship, we have

12_kot ^x2=12_kBT

where⟨^x2⟩is the mean-squared displacement of the particle in the trap,_kBthe Boltzmann constant, andTthe absolute temperature. We determined_kotby measuring the mean-squared displacement⟨^x2⟩.

2.7. Stokes’ Drag Coefficient of the Particle near the Bottom of the ITO Glass Chamber

We used the active microrheology to determined the Stokes’ drag coefficient,_ζStokes [31]. We applied an oscillatory optical force to a trapped bead and analyzed the oscillatory motion of the particle to determine the Stokes’ drag coefficient of the particle in 2D near the bottom of the glass chamber. The trapped bead was forced to oscillate along the x-direction by the oscillating tweezers driven by the piezoelectric lead zirconate titanate (PZT)-controlled mirror. The equation of motion of the trapped bead can be written as

_ζStokesdxdt=_kot[_A0cos(ωt)−x]

where_A0is the maximum displacement of the trap andωis the angular velocity of the oscillation. The steady-state solution of the motion is:

x=D(ω)cos(ωt+δ″)

whereD(ω)is the amplitude of displacement andδ″is the phase lag. Both values were from the lock-in amplifier. The frequency-dependent amplitudes D and the phase delaysδ″ are

D(ω)=_{A0 kotkot2}+^{(_ζStokesΩ)2}

δ″(ω)=^tan−1_ζStokesΩ_kot

Thus,

_ζStokesω_kot=tanδ′

Since we already know_kotfrom above, we can determine the Stokes drag coefficient_ζStokesof the Janus particle near the bottom of the glass plate.

We estimate the distance between the particle bottom and the ITO surface with the Faxen equation. Faxen law describes the drag coefficient change near a flat surface due to the non-slipping boundary condition, which can be written as [34]

_ζStokes6πηr=11−916rh+18^(rh)3−45256^(rh)4−116^(rh)5

whereηis the viscosity of water, r is the radius of Janus particle and h is the distance between the center of particle and the surface. Since we already know_ζStokes, we can determine the gap (h − r) between the particle bottom and the ITO surface.

3. Results and Discussions 3.1. ICEP Movement of an Unconfined Janus Particle in 2D

Due to the high-mass density and micron sizes, these Janus particles tended to sediment to the bottom of the sample chamber. These particles did not stick to the substrate presumably due to Coulomb repulsion between the same negatively charged surface of the indium–tin–oxide (ITO) coated glass slide and that of the Janus particle surfaces. The electric repulsion was sufficiently strong to elevate the particles at a small distance above the bottom surface. We estimated the electrical repulsive force using Derjaguin–Landau–Verwey–Overbeek (DLVO) theory [35] and Derjaguin and superposition approximations [36],

F=4πεκr_{ϕJanuseff ϕITOeff}exp(κz)

whereεis the permittivity of water,κis one over the Debye screening length, r is the radius of Janus particle,_ϕJanuseff and _ϕITOeffare the effective potential of the Janus particle and ITO surface, respectively;zis the distance between the bottom surface of the Janus particle and the ITO surface, and

_{ϕJanus, ITOeff}=(4_kBTe)tanh(e_ϕJanus,ITO4_kBT)

We measured the zeta potential of the Janus particle (_ϕJanus)to be−17.01±1.13mV (ZetaPlus Brookhaven instrument). We used a published value for the zeta potential of ITO surface (_ϕITO)in DI water to be−32.7±0.2 m W [37] to calculate effective potential [36]. We used the conductance of DI water 167 ms (ZetaPlus, Brookhaven Instruments Corp., Holtsville, NY, USA) to determine the Debye screening length to be 5.55 nm. The permittivity of water was80.2_ε0at 25 °C, where_ε0 is the permittivity of vacuum [38]. The gravitational force of 3μmJanus particle was 0.23 pN, calculated with∆ρVg, where∆ρis the mass density difference of silica (2.65 g/^cm3 ) [39] and water (1.00 g/^cm3),Vwas the volume of particle, and g was gravity acceleration. Thus, the gap distance between the particle’s bottom surface and the ITO surface was47±8 nm, which was the balance point of the particle at which the electric repulsive force and gravitational force were equal to each other.

In the absence of an electric field, Janus particles moved as Brownian particles in a 2-dimensional (2D) plane defined by the glass substrate. Video analysis of the particles’ motion in 2D yielded a 2D mean-squared-displacement (MSD) shown in Figure 3. The MSD curve in the absence of the applied electric field was characteristic of a Brownian particle. Using_Dt=d(MSD)4dtand the measured MSD curve, we determined the diffusivity to be0.053±0.008 μ^m2/s in the absence of the applied electric field. With this diffusivity, we calculated the drag coefficient using the Stokes-Einstein relation [40]. The drag coefficient was determined to be7.86±1.37×¹⁰⁻⁸ Ns/m . Using Faxen’s law [34], we determined the distance between the particle’s bottom surface and the ITO surface to be12.8±0.2 nm; the gap size was about 4 times smaller compared to the distance estimated by the balance between the electrostatic repulsion and the weight of the particle. The discrepancy was probably caused by the overestimation of the zeta potentials of the particle surface, and the ITO surface since during the expriment, the pH value of the water could be ~5.6 instead of 7.0, which would reduce the thickness of the EDL of both surfaces.

When an AC electric field was applied, the Janus particles underwent induced-charge dielectrophoresis (ICEP). Modeling the ICEP driven particle as an active swimmer, its MSD vs. time curve was predicted by Marchetti [41] as Equation (19).

MSD=4_Dtt+2_{v02 τr}[t−_τr(1−^e−t_τr)]

According to the equation, a log-log plot of the MSD vs. time curve would yield the active ICEP speed from the intercept of the portion where the slope is 2. Equation (19) could be deducted toMSD=_v02 ^t2, during the time interval between_τrand_vo2/4_Dt. This deduction works when the_vo2/4_Dtwas larger than_τr. The vertical axis intercept of the ballistic motion section of the curve was2ln(_v0), and_v0was calculated from the experimental intercept value.

We examined the size-dependent speed of the Janus particles under the same experimental conditions. The relationship between ICEP active velocity and applied voltage at 5 kHz was shown in Figure 5a. The speed of 3μm, 10μm, and 20μ m of metallic Janus particles were nearly proportional to the square of applied voltage as predicted theoretically by Squires and Bazant [15] for normal ICEP at low frequency. The metallic Janus particle driven by ICEP with movement towards the dielectric side had a speed of:

v(ICEP)=964ϵaη(1+δ)^E2

whereϵis the electric permittivity,ηis the viscosity of the bulk solvent,a is the radius of the particle, E is the electric field, and δ is the ratio between the capacitances of the compact and diffuse layer in the electric double layer. With a uniform thickness in the experimental cell, the ICEP velocity should be proportional to the square of the voltage according to this theory, which agreed with Figure 5a, as well as Gangwal et al.’s experiment data [14]. The size dependence, which was predicted to be linear with particle diameter, however, was not in agreement with our experimental results. In this work, measured v(ICEP) ratios of 3:3.9:6.6, shown in the Figure 5a insert, were lower than the theoretical prediction, which was 3:10:20. It is possible that the discrepancy in size dependence for the largest particles of the group might be due to increased hydrodynamic drag between the bottom of the particle and the glass substrate since gap sizes decrease as the particle sizes increase in our experiment. It is also possible that the discrepancies might be caused by the fact that the parameters δ in these three cases were different.

The frequency-dependent ICEP velocity was measured and shown in Figure 5b. There was a characteristic frequency_ffat which the ICEP velocity was at the maximum in the normal direction (from the metal side toward the dielectric side). As the frequency increases, we observed that the velocity first decreases gradually to zero, then the direction of motion reverses at a frequency_fc, the crossover frequency. The negative velocity (the reverse direction) increased as the frequency continuously increased until a peak was reached at a frequency_fr, after which the ICEP velocity reaches zero at a very high frequency.

According to Squires et al. [15], the ICEP movement is generated by slip flow around the particle which appears in a certain range of the driving frequency,_ωe≤ω≤_ωp. The lower bound_ωeis the minimium charging frequency of the electrode. Here,_ωe=DλLand D is the ion diffusion coefficient, λ is the Debye length, and 2L is the electrode separation. The upper bound,_ωpis the maximum formation frequency of induced ionic cloud screening on the particle, with_ωp=DλR, and R is the particle radius. For our experiment the_ωe=1 kHz, and for the 3μm particle,_ωp=30 kHz . In Figure 5b, our experiment shows there was a peak for3 μmparticles in velocity between 1 kHz–5 kHz. The larger 10μmand 20μmparticles archived maximum speed at even lower frequencies. The frequency range predicted by Squires does include these frequencies where maximum speeds occurred in our experiments. When testing the frequency response of 10μm and 20μm diameter particles at frequencies below 1 kHz, Joule heating was found to not be negligible; water vapor bubbles were observed to form at these low frequencies. For frequencies in the range of 100 kHz to 1 MHz, reversed ICEP was observed for the 3μm particle, where the particle moved towards the metal side. The reversed ICEP had a much lower speed at the same applied voltage (2 V). The 3μm Janus particle seemed to have a higher reversed ICEP speed in comparison to 10μm and 20μm particles. The crossover frequencies for these particles were difficult to determine by video imaging methods because video analysis was not sensitive enough for particles moving at low speeds and when displacements were small.

3.2. ICEP Movements in 1D of a Phoretic Particle Confined in a Quadratic Potential The use of PFS to detect the phoretic motion requires the orientation of the Janus particle to be fixed by a magnetic field, and the amplitude of the AC electric field be modulated. Two permanent magnets applied the magnetic field in a direction perpendicular to the vertical AC electric field. Janus particles with a ferromagnetic Ni layer in the metal cap were used to align the direction of the particle’s largest magnetic moment to the B field. The magnetic field fixed the direction of one axis, and the AC electric field fixed the direction of the other axis of the magnetic Janus particle, thus maintaining the orientation of particles in-plane. With the application of the AC electric field and the DC magnetic fields, the orientation of the particle was fixed in a direction perpendicular to both fields. We used PFS to determine the crossover frequency and phoretic drag coefficient of the phoretic Janus particle which moved in a linear motion. According to Equation (3), particle motion had two harmonic components and an offset. The lock-in amplifier calculated the amplitude and the phase lag of the particle motion with a frequency component in the first harmonic of the amplitude modulation (AM) frequency. Using PFS, the crossover frequency was measured at the point when the phase lag signal shifted 180 degrees. Using Equations (8) and (9), we calculated the drag coefficient and mobility of the ICEP-driven motion under a range of base wave frequencies (1 kHz–1 MHz).

Figure 6a shows a plot of the phase shift vs. frequency. A 180-degree phase shift between 20 kHz to 30 kHz indicated the normal direction of motion changed into the reverse direction; the frequency at which the 180-degree phase shift was identified as the crossover frequency. Figure 6b shows the speed of the Janus particle as calculated byv=_ΩMA(_ΩM), i.e., the maximum speed during the oscillating movement under amplitude-modulated ICEP. At frequencies lower than the crossover frequency (~27 kHz), the particle moved towards the dielectric side, defined as the positive velocity (normal direction). At frequencies higher than the crossover frequency of 27 kHz, the particle moved with the metal side in the front. It is informative to compare the crossover frequency of the 3μm Janus particle with other experiments. Our crossover frequency for the 3μm Janus particle was 30 kHz. Mano et al.’s experiment on a 3μ m rotational Janus particle showed a cross-over frequency at 22 kHz [16]. Suzuki et al.’s work, also using a 3μ m Janus particle, showed a cross-over frequency at around 30 kHz [42]. The reversal of ICEO often showed smaller crossover frequencies of around 5 kHz [17]. The frequency scales, in the range 5 kHz to 30 kHz, revealed a length scale of 0.2~0.6μmusing the length scale defined byL=_Dion/_fc. The ionic diffusivity,_Dion, of K⁺ or Cl⁻ is about2×^{10−9 m2}/s [39]. The 30-microsecond time scale revealed by the crossover frequency was on the same order of magnitude of the diffusion time for ions to diffuse the perimeter of the Janus particle, suggesting that the relaxation of the electrokinetic flows in the electric double layer must play a significant role.

Frequency-dependent ICEP mobilityβ calculated according to Equations (8) and (9), is shown in Figure 6c. The frequency when the speed approached zero in Figure 6c is in agreement with the frequency at which the phase delay changed by 180 degrees. The frequency at which the reversed ICEP speed occurred was the same as that observed in the frequency-dependentβor the frequency-dependent speed. Theβvs. frequency curve looked different from that of the speed vs. frequency because of thesinδterm in the denominator of the right side of Equation (8). Using amplitude-modulated ICEP, the oscillation of the driving force had a phase lag (δ)relative to the particle motion. The in-phase term of the driving force was0.1275_E02sinδ, which was used to calculate theβ. The out-off-phase term of the driving force was0.1275_E02cosδ, which was used to balance the optical trap force. Mobilityβ obtained by the AM modulated, and the phoretic force spectroscopy agreed with the value of the 2D measurements for particles not confined by the optical trapping. The maximum speed in the normal direction was at a frequency of 1 kHz–5 kHz, and the peak frequency for the reverse direction was around 200 kHz. We don’t know exactly why the reversed ICEP had much slower mobility than that of the forward motion. A few papers [14,15,16,17,18,19,20] that mentioned the frequency-dependent reversal of electroosmotic flows did not propose mechanisms that would explain the reversal.

We determined the phoretic drag coefficientξ according to Equation (9), shown in Figure 6d. The drag coefficients were found to be independent of frequency except a few outliers obtained when the particle had very low moving speeds (near crossover frequency) with large error bars. Taking the average over all the frequency, we determined the ICEP drag to be7.5±2.7×¹⁰⁻⁸ N·s/m . The Stokes drag coefficient, measured by oscillating tweezers [29,41] for the same Janus particle in the absence of ICEP, was7.9±1.4×¹⁰⁻⁸ N·s/m, similar to the ICEP drag coefficient.

4. Conclusions This paper reports a study that used phoretic force spectroscopy (PFS) to determine the crossover frequency, the maximum phoretic speeds in the normal and the reverse direction, ICEP mobility, and the ICEP drag coefficient. The direction of the ICEP-driven metallic-dielectric Janus particles was in the normal direction at low frequencies and the reversed direction at high frequencies. Phoretic mobility was proportional to the square of an applied electric field. Our data are in agreement with prior experimental studies. However, our PFS detection method provided better accuracy in determining the crossover frequencies. The 30-microsecond time scale revealed by the crossover frequency is on the same order of magnitude as the diffusion time for ions to diffuse the perimeter of the Janus particle, suggesting the relaxation of the electrokinetic flows in the electric double layer must play a role. This study found the ICEP drag coefficients of the forward and the reversed motion were similar, and both are comparable to the Stokes drag. Unfortunately, large error bars in our data did not permit more detailed analysis to distinguish the ICEP drag from the Stokes drag. Microscopic interpretations of both the frequency at which ICEP mobility switched signs could benefit from the accurately determined crossover frequencies. The reason for the large difference between the magnitudes of ICEP mobility in the forward and reversed modes of ICEP requires further study. With the demonstrated ability to determine the high-resolution frequency-dependent response function of the metallic Janus particle by a phase-sensitive detection method, this study can help broaden the application areas and enhance the performance of current applications of colloidal microrobots. For example, a more accurately determined crossover frequency can provide more precise sorting or separation of particles by exposing unsorted particles to judiciously controlled frequencies in a microfluidic channel setting. Metallic Janus particles coated with layers of different magnetic properties exposed in cleverly arranged magnetic fields could be used to transport different kinds of medicines to different locations. Combining the use of electric fields with different frequencies and/or an added magnetic field might also be useful to control the self-assembly of active colloids to achieve complex and time-varying microstructures.

Figure 1. Reproduction of a schematic by Gangwal et al. of a particle in one-half cycle of an AC electric field in the stable configuration. The electric double layer on the gold side (black hemisphere) is more strongly polarized and thus drives a stronger induced-charge-electroosmosis slip (arrows) than the polystyrene side, resulting in induced-charge electrophoresis (ICEP) motion in the direction of the dielectric side. Reprinted with permission [14].

Figure 2.(a) Fabrication process of Janus particles. Silica particles were deposited from an aqueous suspension to form a sub-monolayer on a glass substrate, followed by E-beam metal evaporation deposition of Ni and Au. (b) On the left is a cartoon showing the shape and geometry of resultant Janus particles, and on the right is a scanning electron microscope (SEM) micrograph of actual particles. (c) A depiction of a trapped Janus particle in the sample chamber. The directions of the particle's motion, the E and B fields, are shown relative to the sample chamber.

Figure 3. Mean-squared-displacement (MSD) of ICEP-driven metallic Janus particle under 5 kHz AC electric field. The extrapolated dash lines are from ballistic motion (the slope of 2) to the active diffusive motion (the slope of 1). The red vertical dotted line represents the particle's rotational relaxation time.

Figure 4. Phoretic force spectroscopy (PFS) set-up and force diagram. (a) Lightpath diagram of PFS. A 1064 nm laser was used to trap a Janus particle. A 980 nm tracking laser was co-aligned with a 1064 nm trapping laser. The 980 nm tracking laser was detected by a quadrant photodiode (QPD). A magnet provided an in-plane magnetic field to further align the Janus Particles (insert). (b) Signal chain diagram. A lock-in amplifier analyzed the signal of particle position with the modulation signal as a reference. (c) Force diagram of an optically trapped Janus particle under a modulated AC electric field. Uot is the potential of a particle in the trap. FICEP (active force) and Fot (trap force) are in balance at a non-center spot. Electric field (E) is in the z-direction, and the magnetic field is in the y-direction. The orientation of the Janus particle is in the x-direction. The particle oscillates, and the oscillation frequency is the same as the frequency of the modulation signal.

Figure 5. ICEP of a Janus particle without confinement. (a) ICEP speed of 3μm, 10μm, and 20μm Janus particles under different applied voltage at 5 kHz. Insert:βvs. the diameter of the particle. The dashed line is the prediction using Equation (20). (b) The velocity of 3μm, 10μm, and 20μm diameter particles under the potential difference 2 V over a 50μm gap between the electrodes and with different frequencies from 100 Hz to 1 MHz. The positive velocity (normal direction) is defined as moving toward the dielectric side.

Figure 6. Experimental results of PFS at 0.1 V/μm of the carrier frequency of the electric field from 1 kHz to 1 MHz with 1 Hz amplitude modulation frequency at 30% modulation amplitude. The vertical dash rad lines(fc)represent the crossover frequency at 27 kHz. The vertical dash green line(fr)represents the maximum reversal velocity at 220 kHz. The vertical dash purple line(ff)represents the maximum forward velocity at 1-5 kHz. (a) Phase delays of a 3μm Janus particle. (b) Speeds of the same 3μm Janus particle. (c) Comparison of the Beta values of the Janus particles measured with trapping (Equation (3)) and without trapping. (d) Drag coefficients determined from the phase delay data.

Author Contributions

Conceptualization, C.S. and H.D.O.; methodology, C.S., J.F.G. and H.D.O.; software, C.S.; data curation, C.S., Z.J. and L.L.; writing-original draft preparation, C.S., Z.J. and H.D.O.-Y.; writing-review and editing, C.S., L.L., Z.J., J.F.G. and H.D.O.-Y.; supervision, H.D.O.-Y.; funding acquisition, H.D.O.-Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was supported in part by the Emulsion Polymers Institute and the Department of Physics of Lehigh University. The authors thank Thitiporn Kaewpetch of the Department of Chemical and Biomolecular Engineering for assistance in preparing the monolayer particle arrays and Theadosia Kurniawan for critical reading of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

1. Choi, J.; Zhao, Y.; Zhang, D.; Chien, S.; Lo, Y.H. Patterned Fluorescent Particles as Nanoprobes for the Investigation of Molecular Interactions. Nano Lett. 2003, 3, 995-1000.

2. Yang, Y.; Bevan, M.A. Optimal Navigation of Self-Propelled Colloids. ACS Nano 2018, 12, 10712-10724.

3. Solovev, A.A.; Sanchez, S.; Pumera, M.; Mei, Y.F.; Schmidt, O.C. Magnetic Control of Tubular Catalytic Microbots for the Transport, Assembly, and Delivery of Micro-Objects. Adv. Funct. Mater. 2010, 20, 2430-2435.

4. Walther, A.; Müller, A.H.E. Janus Particles: Synthesis, Self-Assembly, Physical Properties, and Applications. Chem. Rev. 2013, 113, 5194-5261.

5. Bogue, R. The Development of Medical Microrobots: A Review of Progress. Ind. Rob. 2008, 35, 294-299.

6. Shields, C.W.; Velev, O.D. The Evolution of Active Particles: Toward Externally Powered Self-Propelling and Self-Reconfiguring Particle Systems. Chem 2017, 3, 539-559.

7. Boymelgreen, A.M.; Balli, T.; Miloh, T.; Yossifon, G. Active Colloids as Mobile Microelectrodes for Unified Label-Free Selective Cargo Transport. Nat. Commun. 2018, 9, 1-8.

8. Ebbens, S.J. Active Colloids: Progress and Challenges towards Realising Autonomous Applications. Curr. Opin. Colloid Interface Sci. 2016, 21, 14-23.

9. Jiang, H.R.; Yoshinaga, N.; Sano, M. Active Motion of a Janus Particle by Self-Thermophoresis in a Defocused Laser Beam. Phys. Rev. Lett. 2010, 105, 268302.

10. Uspal, W.E.; Popescu, M.N.; Dietrich, S.; Tasinkevych, M. Self-Propulsion of a Catalytically Active Particle near a Planar Wall: From Reflection to Sliding and Hovering. Soft Matter 2015, 11, 434-438.

11. Kilic, M.S.; Bazant, M.Z. Induced-Charge Electrophoresis near a Wall. Electrophoresis 2011, 32, 614-628.

12. Sanchez, S.; Solovev, A.A.; Harazim, S.M.; Schmidt, O.G. Microbots Swimming in the Flowing Streams of Microfluidic Channels. J. Am. Chem. Soc. 2011, 133, 701-703.

13. Zhang, L.; Zhu, Y. Directed Assembly of Janus Particles under High Frequency Ac-Electric Fields: Effects of Medium Conductivity and Colloidal Surface Chemistry. Langmuir 2012, 28, 13201-13207.

14. Gangwal, S.; Cayre, O.J.; Bazant, M.Z.; Velev, O.D. Induced-Charge Electrophoresis of Metallodielectric Particles. Phys. Rev. Lett. 2008, 100, 058302.

15. Squires, T.M.; Bazant, M.Z. Breaking Symmetries in Induced-Charge Electro-Osmosis and Electrophoresis. J. Fluid Mech. 2006, 560, 65-101.

16. Mano, T.; Delfau, J.B.; Iwasawa, J.; Sano, M. Optimal Run-And-Tumble-Based Transportation of a Janus Particle with Active Steering. Proc. Natl. Acad. Sci. USA 2017, 114, E2580-E2589.

17. Storey, B.D.; Edwards, L.R.; Kilic, M.S.; Bazant, M.Z. Steric Effects on Ac Electro-Osmosis in Dilute Electrolytes. Phys. Rev. E 2008, 77, 036317.

18. Bazant, M.Z.; Kilic, M.S.; Storey, B.D.; Ajdari, A. Towards an Understanding of Induced-Charge Electrokinetics at Large Applied Voltages in Concentrated Solutions. Adv. Colloid Interface Sci. 2009, 152, 48-88.

19. Ng, W.Y.; Lam, Y.C.; Rodríguez, I. Experimental Verification of Faradaic Charging in Ac Electrokinetics. Biomicrofluidics 2009, 3, 022405.

20. González, A.; Ramos, A.; García-Sánchez, P.; Castellanos, A. Effect of the Combined Action of Faradaic Currents and Mobility Differences in Ac Electro-Osmosis. Phys. Rev. E. 2010, 81, 016320.

21. Ramos, A.; Morgan, H.; Green, N.G.; González, A.; Castellanos, A. Pumping of Liquids with Traveling-Wave Electroosmosis. J. Appl. Phys. 2005, 97, 084906.

22. Yang, H.; Jiang, H.; Shang, D.; Ramos, A.; Garcia-Sanchez, P. Experiments on Traveling-Wave Electroosmosis: Efect of Electrolyte Conductivity. IEEE Trans. Dielectr. Electr. Insul. 2009, 16, 417-423.

23. Wei, M.T.; Junio, J.; Ou-Yang, D.H. Direct Measurements of the Frequency-Dependent Dielectrophoresis Force. Biomicrofluidics 2009, 3, 012003.

24. Liao, M.-J.; Wei, M.-T.; Xu, S.-X.; Daniel Ou-Yang, H.; Sheng, P. Chinese Physics B Non-Stokes Drag Coefficient in Single-Particle Electrophoresis: New Insights on a Classical Problem. Chin. Phys. B 2019, 28, 084701.

25. Park, H.; Wei, M.T.; Ou-Yang, H.D. Dielectrophoresis Force Spectroscopy for Colloidal Clusters. Electrophoresis 2012, 33, 2491-2497.

26. Muangnapoh, T.; Weldon, A.L.; Gilchrist, J.F. Enhanced Colloidal Monolayer Assembly via Vibration-Assisted Convective Deposition. Appl. Phys. Lett. 2013, 103, 181603.

27. Lin, C.H.; Chen, Y.L.; Jiang, H.R. Orientation-Dependent Induced-Charge Electrophoresis of Magnetic Metal-Coated Janus Particles with Different Coating Thicknesses. RSC Adv. 2017, 7, 46118-46123.

28. Han, M.; Yan, J.; Granick, S.; Luijten, E. Effective Temperature Concept Evaluated in an Active Colloid Mixture. Proc. Natl. Acad. Sci. USA 2017, 114, 7513-7518.

29. Sbalzarini, I.F.; Koumoutsakos, P. Feature Point Tracking and Trajectory Analysis for Video Imaging in Cell Biology. J. Struct. Biol. 2005, 151, 182-195.

30. Shen, C.; Ou-Yang, H.-C.D. The Far-from- Equilibrium Fluctuation of an Active Brownian Particle in an Optical Trap. In Optical Trapping and Optical Micromanipulation XVI, Proceedings of the SPIE Nanoscience + Engineering, San Diego, CA, USA, 11-15 August 2019; Dholakia, K., Spalding, G.C., Eds.; SPIE: Bellingham, WA, USA, 2019; Volume 11083, p. 63.

31. Valentine, M.T.; Dewalt, L.E.; Ou-Yang, H.D. Forces on a Colloidal Particle in a Polymer Solution: A Study Using Optical Tweezers. J. Phys. Condens. Matter 1996, 8, 9477-9482.

32. Wang, J.; Wei, M.T.; Cohen, J.A.; Ou-Yang, H.D. Mapping Alternating Current Electroosmotic Flow at the Dielectrophoresis Crossover Frequency of a Colloidal Probe. Electrophoresis 2013, 34, 1915-1921.

33. Zia, R.N. Active and Passive Microrheology: Theory and Simulation. Annu. Rev. Fluid Mech. 2018, 50, 371-405.

34. Ha, C.; Ou-Yang, H.D.; Pak, H.K. Direct Measurements of Colloidal Hydrodynamics near Flat Boundaries Using Oscillating Optical Tweezers. Phys. A Stat. Mech. Appl. 2013, 392, 3497-3504.

35. BUTLER, J.A.V. Theory of the Stability of Lyophobic Colloids. Nature 1948, 162, 315-316.

36. Behrens, S.H.; Plewa, J.; Grier, D.G. Measuring a Colloidal Particle's Interaction with a Flat Surface under Nonequilibrium Conditions Total Internal Reflection Microscopy with Absolute Position Information. Eur. Phys. J. E 2003, 10, 115-121.

37. Kline, T.R.; Chen, G.; Walker, S.L. Colloidal Deposition on Remotely Controlled Charged Micropatterned Surfaces in a Parallel-Plate Flow Chamber. Langmuir 2008, 24, 9381-9385.

38. Kaatze, U. Complex Permittivity of Water as a Function of Frequency and Temperature. J. Chem. Eng. Data 1989, 34, 371-374.

39. Lide, D.R. CRC Handbook of Chemistry and Physics; CRC Press: Boca Raton, FL, USA, 2004; ISBN 0-8493-0487-3.

40. Kholodenko, A.L.; Douglas, J.F. Generalized Stokes-Einstein Equation for Spherical Particle Suspensions. Phys. Rev. E 1995, 51, 1081-1090.

41. Marchetti, M.C.; Fily, Y.; Henkes, S.; Patch, A.; Yllanes, D. Minimal Model of Active Colloids Highlights the Role of Mechanical Interactions in Controlling the Emergent Behavior of Active Matter. Curr. Opin. Colloid Interface Sci. 2016, 21, 34-43.

42. Suzuki, R.; Jiang, H.-R.; Sano, M. Validity of Fluctuation Theorem on Self-Propelling Particles. arXiv 2011, arXiv:1104.5607.