# Transient electrophoretic current in a nonpolar solvent

## Abstract

The transient electric current of surfactants dissolved in a nonpolar solvent is investigated both experimentally and theoretically in the parallel-plate geometry. Due to a low concentration of free charges the cell can be completely polarized by an external voltage of several volts. In this state, all the charged micelles are compacted against the electrodes. After the voltage is set to zero the reverse current features a sharp discharge spike and a broad peak. This shape and its variation with the compacting voltage are reproduced in a one-dimensional drift-diffusion model. The model reveals the broad peak is formed by a competition between an increasing number of charges drifting back to the middle of the cell and a decreasing electric field that drives the motion. After complete polarization is achieved, the shape of the peak stops evolving with further increase of the compacting voltage. The spike-peak separation time grows logarithmically with the charge content in the bulk. The time peak is a useful measure of the micelle mobility. Time integration of the peak yields the total charge in the system. By measuring its variation with temperature, the activation energy of bulk charge generation has been found to be 0.126 eV.

## I Introduction

Electrical conduction in nonpolar fluids is of interest in relation to electrophoretic displays that use nonpolar-based solutions as the working medium. Muerau1978 (); Novotny1979b (); Jacobson1998 () When an appropriate surfactant is added to a nonpolar solvent its molecules form inverse micelles: globular aggregates with hydrophobic surfaces and hydrophilic interiors. Eicke1980 () The micelles are believed to be responsible for the solution’s nonzero electrical conductivity by stabilizing positive and negative charges and keeping them from recombining. Morrison1993 () In addition, the micelles have been suggested to charge colloidal particles by ionizing their surface groups and preferentially adsorbing on the surfaces. Roberts2008 ()

A useful way to gain insight into the microscopic properties of inverse micelles is through a detailed analysis of the transient current curves. Novotny1979 (); Denat1982b (); Novotny1986 (); Dikarev1997 (); Bazant2004a (); Kim2005a (); Bert2006a (); Beunis2007a (); Beunis2007b (); Beunis2008 (); Prieve2008 () Typically, in such an experiment a single-surfactant solution is placed between two parallel electrodes, although more complex geometries have been used. Kim2005a () The voltage is stepped from zero to a finite value and then either back to zero or to some other value. The current recorded in the external circuit is an integral measure of the complex movement of charges inside the cell. By assuming the solution to be a symmetric electrolyte, the concentration of charged micelles, their mobility, diffusion coefficient, and the Stokes diameter can be inferred.

Most of the previous work has focused on the sometimes very detailed quantitative analysis of the forward transient current. Bazant2004a (); Beunis2008 () However, in the early studies it was noticed that the reverse transients often exhibit a nonmonotonic behavior. Novotny1979b (); Novotny1986 () Although the reverse peak has been observed in subsequent studies, Bert2005 (); Bert2006a (); Beunis2007b () it has not been analyzed or explained in sufficient detail. The present paper aims to contribute to this area.

## Ii Experiment

Poly-isobuthylene succinimide (OLOA 11000, Shevron Oronite Company LLC, 72% active component)
was dissolved in the hydrocarbon fluid Isopar^{®} M (Exxon Mobil
Chemical, dielectric constant ) by sonication at different concentrations.
No attempt was made to dry the components and the water content was not controlled.
Parallel-plate cells were made from pieces of polyethylene terephthalate (PET) coated with conducting
indium tin oxide (Sheldahl, Inc.). One side of the cell was photo-patterned with a square grid of
SU8 spacer walls whereas the other side remained flat. The width of the walls was 10 m and
the grid period was 230 m. The height of the walls (10 or 30 m) defined the cell
thickness. The active device area was about 10 cm. The voltage profile was provided by a DC
power supply (Agilent Technologies). The current was recorded by a low-noise current preamplifier
(Stanford Research Systems) at a sampling rate of 1000 Hz. Both instruments were controlled by a
customized LabView application (National Instruments Corp.). To suppress the Faradaic processes,
the measurements were done at below-room temperatures of C, C, and
C. The temperature was controlled by a convective environmental chamber
(Autronic-Melchers GmbH) with accuracy of C.

## Iii Observations

Transient currents of a 10 -thick 0.5 wt.% OLOA-11000 solution in Isopar M are shown in Fig. 1(a). The voltage between two parallel electrodes was stepped to a nonzero value between 0.25 and 10.0 V for 15 s, and then brought back to 0 V for another 14 seconds. As observed by other investigators Novotny1979 (); Novotny1986 (); Bert2005 (); Bert2006a (); Beunis2007b () the currents possess the following characteristic features: a fast decay in the forward direction, then a sharp spike followed by a broad peak in the reverse direction. Initially, the forward decay is fast but then it slows down. Phenomenologically, the forward transient curve is well described by an circuit with a decreasing number of charge carriers (increasing ). This issue is not considered in the present work. The forward current continues to evolve with increasing voltage: it grows larger in magnitude and shorter in duration. In contrast, the reverse current quickly saturates with voltage. The and 10.0 V curves are nearly identical. Figure 1(b) shows a magnified view of the reverse peak area. The peak develops between 0.5 and 0.75 V and saturates by V. This suggests 3 V cause complete polarization of the cell: all available charges are compacted near the electrodes and the charge cloud no longer changes with a further increase in voltage.

Variation of the transient current with the cell thickness has been investigated. The reverse current for two thicknesses 10 m and 30 m are displayed in Fig. 1(c). As expected intuitively, in the thicker cell, the reverse peak develops more slowly. The ratio of the peak times, , is close to the ratio of the cell thicknesses squared, . This scaling is discussed in the following sections. Variation of the current with temperature has also been investigated. Figure 1(d) compares the reverse transients for the three temperatures of C, C, and C. As the temperature decreases the peak shifts to longer times and reduces in height.

## Iv Model

To develop an understanding of the observed behavior, a one-dimensional drift-diffusion (Nernst-Planck-Poisson) model has been applied. The model is formulated here in a dimensionless form. Bazant2004a () Given the distance between two parallel electrodes , the diffusion coefficient of charged micelles (both positive and negative), and the physical distance and time , the dimensionless distance and time are introduced as and . The electrostatic potential is measured in units of the thermal voltage: . The micelles are assumed to be single-charged. The concentration of micelles is measured in units of the equilibrium initial concentration of one sign, . The micelle mobility is related to the diffusion coefficient by the Einstein relation . In such units, the micelle flux is given by a sum of the drift and diffusion contributions

(1) |

In the simplest version of the model, the dissociation-recombination processes in the bulk Silver1965 (); Novotny1986 (); Strubbe2006 () and the Faradaic processes at the electrodes Bockris1970 (); Novotny1986 () are neglected. (Both are small for the experimental times and temperatures studied in this work.) As a result, the fluxes satisfy the local conservation laws

(2) |

and the blocking boundary conditions

(3) |

These conditions also imply global conservation of positive and negative charges. The electrostatic potential is governed by the Poisson equation

(4) |

and boundary conditions

(5) |

The step-function boundary conditions are defined as

(6) |

In Eq. (4), is a dimensionless parameter which is an equivalent measure of the total charge concentration. It can also be expressed as , where is the Debye length of the undisturbed solution. Another parameter is the magnitude of the voltage step . The model is fully characterized by just two dimensionless parameters, and . As such, it allows for complete numerical investigation. The initial conditions are uniform concentrations and a zero electrostatic potential:

(7) |

Once is known, the external current (in amperes) is computed as the time derivative of the electric field at the electrode (as follows from the Gauss law):

(8) |

where is the electrode area.

The Cauchy problem (1)-(7) has been solved numerically using MATLAB (The MathWorks, Inc). The solution for several values of and is shown in Fig. 2. The model transient possesses the same qualitative features as the experimental one: a fast monotonic decay in the forward direction, then a sharp spike followed by a broad peak in the reverse direction. The reverse current saturates with . The shape of the reverse peak slowly evolves with increasing .

For the following discussion it is instructive to rewrite formula (8) to make the relationship between the external current and the distribution of charges inside the cell more explicit. Integrating by parts and applying the Poisson equation (4) one obtains for the potential difference between the electrodes

(9) |

Note that due to global electroneutrality the electric fields near the two electrodes are always equal, or . Substituting (9) into (8) one obtains

(10) |

The dimensional form of the last formula reads

(11) |

where the dot means a partial derivative with respect to the real time . The external current comprises two different contributions. The first term is charging and discharging of the geometric capacitor. (The coefficient by is recognized as the geometric capacitance of the cell .) The second term is an additional current associated with the movement of ions inside the cell. The first process is typically much faster than the second, and results in sharp charging/discharging initial spikes.

## V Origin of the reverse peak

In addition to the external current, the numerical solution provides time-dependent distributions of the charges and the fields inside the cell, some of which are shown in Fig. 3. The origin of the reverse peak can be analyzed in sufficient detail, which is now described.

During the forward phase the external electric field polarizes the cell by compacting the positive and negative ions against their respective electrodes. The equilibrium concentration profile is determined by the balance between the electrophoretic drift and diffusion when it reaches an equilibrium, cf. Fig. 3(a), . As the field is turned off, the electrodes are brought to the same potential on a time scale much shorter than the characteristic diffusion time of the ions. A new potential distribution is established almost instantaneously after the switching at , cf. Fig. 3(c), . This fast geometric capacitance discharge produces the sharp initial spike.

Just after the switching most ions are still being pushed against the electrodes by the local electric field. However, the electrostatic potential of both electrodes is now zero. By symmetry, the potential in the cell center is also zero. (In a more general case of the asymmetric electrolyte or cell, the zero-potential point will not be in the center of the cell. This fact does not change the following arguments.) Therefore the potential assumes a local maximum somewhere near the negative electrode (and a local minimum near the positive electrode). This zero-field point must be within the charge cloud since by Poisson equation only uncompensated charge can change the curvature of the potential. As a result, to the right of the maximum point there is some positive charge, which is of the order of the geometric capacitance times the minimal voltage of total polarization. (The estimate follows from the Gauss law.) Since the electric field is also positive there, the charge-field product creates a flux of positive ions toward the center of the cell. A similar drift of the negative ions away from the positive electrode is taking place at the opposite side of the cell. Closer to the electrodes, the drift and diffusion fluxes continue to effectively compensate each other.

Initially, the integral ionic flux is small and, according to Eq. (12), the external current is also small. However, the external current is a flow of the electrode charge, and as a result the latter decreases. The reduction of the electrode charge moves the zero-field point deeper into the charge cloud thereby releasing more ions for drift-diffusion toward the bulk. The electric field increases inside the cell for the same reason. The integral ionic flux grows and so does the external current. Thus the current (i.e. time decrement of the electrode charge) is linearly related to the electrode charge itself. As a result, the current grows approximately exponentially, at least for short times when the positive and negative charge clouds do not overlap. As the clouds expand, positive and negative ions begin to mix in the middle of the cell reducing the space-charge and the electric field. The charge distribution becomes more uniform which dampens the internal fluxes and the external current. The entire process can be viewed as a competition between the increasing number of participating ions and decreasing driving fields. In the early stages of the reversal the increase in the ion number dominates, and the overall current increases with time. Later, the ion number levels off but the driving electric field subsides as does the current.

The variation of the reverse current with the charge parameter is now discussed. As can be seen in Fig. 2, larger s require larger voltages to reach current saturation. In addition, the peak maximum shifts to longer times. The voltage increase can be understood from the following considerations. The reverse current saturates when the cell is fully polarized which occurs when . Here is the Gouy-Chapman double-layer capacitance and . At low charge densities, the polarization voltage scales as . At moderate and high densities the exponential enhancement of the capacitance is most important, and the scaling crosses over to . It is essential that the polarization voltage scales sublinearly with the density. The above estimates do not take into account the significant effects of the limited charge content in a finite cell on the double layer capacitance and polarization voltage. This issue was thoroughly analyzed by Verschueren et al. Verschueren2008 ()

As mentioned above, the reverse drift starts with a small amount of charge and then accelerates exponentially by pulling ions from the cloud near the electrodes. The process continues until the charge “stocks” are depleted. The depletion time approximately corresponds to the peak time in the reverse current. Therefore the peak time should scale logarithmically with the total amount of charge in the system, . This scaling is confirmed by numerical calculations, as shown in Fig. 4. Notice the envelope of the curves’ leading edges is roughly an exponential.

## Vi Micelle properties

The broad peak is the most conspicuous feature of the reverse transient current in the high-voltage regime. In this section it will be shown how useful physical information can be extracted from the time and shape of the peak.

Diffusivity. As can be seen in Fig. 4, the weak logarithmic dependence places the peak time within the interval for a large spread of the density parameter (a factor of 30). This suggests a simple method of estimating the ion’s diffusivity: measure the reverse peak time in seconds, then with the accuracy of 50%. For example, for the data of Fig. 1(d), the peak times are 2.00, 1.65, and 1.40 sec for C, C, and C. Using , one obtains the respective diffusivities of OLOA11000 in Isopar M at these temperatures as , 6.1, and 7.1 . By the same argument, the peak time should scale approximately quadratically with the cell thickness , cf. Fig. 1(c).

A more accurate method involves taking into account the changes of the peak shape with . A convenient measure of the current asymmetry is its skewness. Figure 5 shows the peak time at saturation versus skewness, as calculated from the model data of Fig. 4. To avoid ambiguity in determining the skewness the data set is limited to a symmetric interval around the peak time. Once the skewness is calculated for an experimental reverse current the dimensionless peak time can be found in Fig. 5. Then the micelle diffusivity is . The skewness of the experimental data of Fig. 1 are given in the fourth column of Table 1. The next column contains the respective dimensionless peak times determined from Fig. 5. They are 11-13 % larger than the “back-of-the-envelope” value 0.1, which corrects the diffusivities to their final values given in the sixth column of the table.

, | , C | , s | Skewness | , m/s | , s/kg | , mPas | , nm | , nC | , nmol/l | |

10 | 5.0 | 2.00 | -0.28 | 0.115 | 5.91 | 12.0 | 108.2 | 11.23 | ||

10 | 10.0 | 1.65 | -0.32 | 0.116 | 4.81 | 12.3 | 117.0 | 12.15 | ||

10 | 15.0 | 1.40 | -0.38 | 0.118 | 4.07 | 12.3 | 129.9 | 13.49 |

Micelles Stokes mobility and Stokes diameter. Once the diffusivity is known the Stokes mobility and Stokes diameter

(13) |

can be easily determined. Here is the absolute temperature in K, and is the viscosity of the host fluid. The viscosity of Isopar M has been measured for . Note that the Stokes mobility differs from the more familiar electrophoretic mobility by the elementary charge factor, . The values of , and for three temperatures are given in Table 1. Both and show significant temperature dependence while the Stokes diameter does not. This suggests that the bulk of the temperature variation in the micelle mobility is due to the corresponding change of the viscosity of the host fluid. A similar conclusion was reached by Novotny. Novotny1986 () The overall estimate of the micelle diameter is nm.

A similar analysis of the sample at C resulted in a slightly larger mobility of s/kg and a slightly smaller micelle diameter of 11.8 nm. Note that the reverse transient current, cf. Fig. 1(c), is more skewed than the curves. The corresponding skewness parameter, -0.76, is about twice as large.

Equilibrium charge density. The saturated reverse transient curve allows an independent determination of the total amount of charge in the system. The Faradaic processes are suppressed at low temperatures and the bulk charge generation is too slow to affect the overall charge dynamics on the time scale of seconds. Strubbe2006 () Under these assumptions, the total amount of the charge of each polarity is conserved. This can be determined from the experimental transient curves. The argument is based on Eq. (11). By integrating the second term from the moment just after the voltage turn off (i.e., just after the discharge spike) to infinity, one obtains

(14) |

The contribution vanishes since , the micelles are fully mixed, and the integrand is identically zero. The volumetric densities are those just after the voltage turn off. In the ideal case of complete compaction, the densities are delta functions: , . Performing the integration, the right hand side of Eq. (14) becomes . Since is the cell volume, the last expression is the total negative charge of the cell. In practice, the micelles are never compacted into a delta function; one can only speak about different degrees of approaching the limit. (In the experiment presented here, the compaction is believed to be almost complete, within a few percent of the limit.) The current integral provides the lower limit for the total charge:

(15) |

Once are known, the equilibrium concentrations of charged micelles are found by dividing by the cell volume. The values of and are given in the last two columns of Table 1. The order of magnitude is mol/l mol/m, which is similar to the values reported in the literature. Kim2005a ()

Activation energy. The temperature dependence of the equilibrium charge density allows us to estimate the activation energy of the charge generation process, which is of fundamental importance in dielectric fluids. Fitting the values from Table 1, the activation energy is determined as = 0.126 eV = 2.91 kcal/mol. This value is close to the one reported by Novotny Novotny1986 () for Aerosol OT surfactant (0.140 eV), which suggest that the ionization mechanisms in the two systems might be similar.

## Vii Summary

In summary, transient currents in a weakly conducting solution of poly-isobuthylene succinimide in Isopar M have been investigated. The limited charge content allows the complete polarization of the cell at low and medium voltages. Under these conditions, the bulk of the cell is devoid of mobile charges since bulk generation is a relatively slow process. The ability to control the bulk charge density makes this system a unique physical laboratory.

The spike-peak structure is the robust feature of the experimental and model transient currents. The development of the peak signals the onset of cell polarization. The reverse current is nonmonotonic because even after the voltage turn off, the electrodes carry a significant amount of charge. Initially, the latter keeps a significant portion of the mobile charges inside the cell close to the electrodes. As a result, the internal charges are released into the bulk not all at once, but gradually. At high external voltages the peak shape saturates. All of these features can be reproduced in a simple one-dimensional drift-diffusion model. Such a model has been formulated in dimensionless units and solved numerically for a range of the two model parameters. The dimensionless peak time increases logarithmically with the charge density and lies in the interval for a wide range of parameters. Measured in real seconds, the peak time scales quadratically with the cell gap and inversely with the micelle mobility. The scaling provides a simple way to estimate the micelle diffusivity as . Taking into account the skewness of the peak allows a more accurate determination of the micelle parameters, which have been summarized in Table 1. The Stokes diameter of the micelles is estimated to be nm. From the temperature dependence of the total charge density, the activation energy of 0.126 eV has been inferred.

###### Acknowledgements.

The authors wish to thank Gregg Combs for providing the current measurement hardware and software, Casey Lohrman for viscosity measurements, Tim Koch, Richard Henze and Kenneth Abbott for support and encouragement, and all the colleagues at Hewlett-Packard for numerous discussions on the subject of this paper.### References

- P. Mürau and B. Singer, J. Appl. Phys. 49, 4820 (1978).
- V. Novotny and M. A. Hopper, J. Electrochem. Soc. 126, 2211 (1979).
- B. Comiskey, J. D. Albert, H. Yoshizawa, and J. Jacobson, Nature 394, 253 (1998).
- H.-F. Eicke, “Micelles,” (Springer, Berlin, Heidelberg, 1980) Chap. 2, Surfactants in Nonpolar Solvents. Aggregation and Micellization, pp. 85–145.
- I. D. Morrison, Colloids and Surfaces A 71, 1 (1993).
- G. S. Roberts, R. Sanchez, R. Kemp, T. Wood, and P. Bartlett, Langmuir 24, 6530 (2008).
- V. Novotny and M. A. Hopper, J. Electrochem. Soc. 126, 925 (1979).
- A. Denat, B. Gosse, and J. P. Gosse, J. Electrostatics 12, 197 (1982).
- V. Novotny, J. Electrochem. Soc. 133, 1629 (1986).
- B. N. Dikarev, G. G. Karasev, V. I. Bolshakov, R. G. Romanets, and I. V. Potapov, J. Electrostatics 40-41, 147 (1997).
- M. Z. Bazant, K. Thornton, and A. Ajdari, Phys. Rev. E 70, 021506 (2004).
- J. Kim, J. L. Anderson, S. Garoff, and L. J. M. Schlangen, Langmuir 21, 8620 (2005).
- T. Bert, H. D. Smet, F. Beunis, and K. Neyts, Displays 27, 50 (2006).
- F. Beunis, F. Strubbe, K. Neyts, and A. R. M. Verschueren, Applied Physics Letters 90, 182103 (2007).
- F. Beunis, F. Strubbe, M. Marescaux, K. Neyts, and A. R. M. Verschueren, Applied Physics Letters 91, 182911 (2007).
- F. Beunis, F. Strubbe, M. Marescaux, J. Beeckman, K. Neyts, and A. R. M. Verschueren, Phys. Rev. E 78, 011502 (2008).
- D. C. Prieve, J. D. Hoggard, R. Fu, P. J. Sides, and R. Bethea, Langmuir 24, 1120 (2008).
- T. Bert, H. D. Smet, F. Beunis, and F. Strubbe, Opto-Electronics Review 13, 281 (2005).
- M. Silver, J. Chem. Phys. 42, 1011 (1965).
- F. Strubbe, A. R. M. Verschueren, L. J. M. Schlangen, F. Beunis, and K. Neyts, Journal of Colloid and Interface Science 300, 396 (2006).
- J. O. Bockris and A. K. N. Reddy, Modern Electrochemistry, Vol. 2 (Plenum, New York, 1970).
- A. R. M. Verschueren, P. H. L. Notten, L. J. M. Schlangen, F. Strubbe, F. Beunis, and K. Neyts, Journal of Physical Chemistry B 112, 13038 (2008).