A Simple Method to Determine Critical Coagulation Concentration from Electrophoretic Mobility
Marco Galli Szilard Saringer Istvan Szilagyi Gregor Trefalt
cond-mat.soft





§ INTRODUCTION

Dispersions of nano or colloidal particles attract widespread contemporary interest due their extensive use in various processes in fundamental research and in more applied disciplines, in which the phyisco-chemical properties of the particles and the media vary in a wide range <cit.>. Applications include catalysis <cit.>, energy storage <cit.>, sensing <cit.>, drug delivery <cit.> and other biomedical utilizations <cit.>, however, the stability of these dispersions is always a key issue.
A typical example is the use of particles as catalysts in liquid media, where stable dispersions (i.e., homogeneously distributed primary particles) are requested during the catalytic run, while the samples can be destabilized (i.e., particle aggregation occurs) once the reaction is terminated <cit.>. The bigger aggregates then sediment or cream according to their density and they can be removed from the system by filtration.
Moreover, aggregation of particles must be suppressed during biomedical delivery processes, since the formation of aggregates in biofluids may cause health complications such as thrombosis by blocking the veins <cit.>. Stable particle dispersions are also required in product manufacturing processes in the textile <cit.>, food <cit.> and cosmetic <cit.> industry as well as in material science, where these dispersions provide processable source of particles for building up composite materials <cit.>.
In contrast, other applications rely on destabilization of particle dispersions by induced aggregation. A typical example is secondary water treatment, in which multivalent ions and polyelectrolytes act as aggregating agents of dispersed dust particles, which sediment and hence, can be eliminated from the tanks in the aggregated form <cit.>. Nanoparticles can also be used to eliminate toxic contaminants from waters <cit.> and the removal of these particles also occurs by induced aggregation. Besides, the paper industry uses particles during production to improve certain properties of the paper and the particles undergo coaggregation with cellulose fibers in the papermaking process <cit.>.

It is evident from the above examples that colloidal stability of particles is a key factor in most of the applications. The usual quantity to estimate colloidal stability of a particle dispersion is the critical coagulation concentration (CCC) or critical coagulation ionic strength (CCIS) indicating the necessary electrolyte concentration for the destabilization of the sample <cit.>. In other words, the particle collision efficiency becomes unity at the CCC, while this efficiency decreases with decreasing the concentration of the aggregating agent and becomes unmeasurable for highly stable samples, in which the particles do not form dimers after collisions.
The accurate knowledge of the CCC is therefore a critical issue to estimate colloidal stability and to design stable or unstable particle dispersions. It can be determined by measuring aggregation rates of particles with a suitable technique <cit.>. One of the most handy methods involves light scattering either in static or dynamic mode <cit.>. In brief, time-resolved measurements are applied and the change in the scattered intensity or in the hydrodynamic radius can be used to calculate aggregation rate constants. These constants were also determined in turbidity measurements by following the changes of the transmittance <cit.> or absorbance <cit.> data in time-resolved experiments. A more advanced, but time-consuming approach is to measure inter-particle forces with the multi-particle colloidal probe technique based on atomic force microscopy <cit.>. From the primary experimental force-distance data, the rates can be calculated at different aggregating agent concentrations using appropriate theories and thus, the CCC can be determined.

The CCCs are typically strongly decreased in the presence of multivalent counterions. The correlation between valence of the counterion and CCC was observed more than 100 years ago by Schulze and Hardy <cit.>. This strong correlation is refereed to as the Schulze-Hardy rule. Later it was realized that also strongly adsorbing monovalent ions can substantially shift the CCC to lower values <cit.>. Therefore, the strong interaction of the counterion with the surface is the key parameter for determining the CCC. The CCC values often correlate with the magnitude of charge or electrokinetic potential of the particles, i.e., particle aggregation occurs at low magnitudes of electrokinetic potentials <cit.>. However, a threshold electrokinetic potential, under which the dispersion can be considered as unstable, cannot be generated for the individual systems, since aggregation processes depend on several factors including size of particles and ionic composition of the surrounding solution.

Here, we propose and implement a simple method for calculation of CCC, based on electrophoretic mobility. Electrophoretic mobility data (which are correlated with electrokinetic potentials <cit.>) of colloidal or nanoparticles measured at different electrolyte concentrations were used to calculate the CCC of various dispersions containing organic or inorganic charged particles of different shapes. The developed method relies on the Derjaguin, Landau, Verwey, and Overbeek (DLVO) theory, which takes into account electrostatic and van der Waals interactions. The calculated CCC values were compared to experimental ones determined independently by suitable techniques.



§ CALCULATING CRITICAL COAGULATION CONCENTRATION


In order to calculate the critical coagulation concentration (CCC), appropriate model has to be chosen. Here we employ DLVO theory for these calculations, since we are interested in aggregation of charged colloidal particles. In some systems non-DLVO interactions importantly shift the CCC, and in these cases the proposed approach is not applicable, as it will be seen in the results section. The DLVO theory assumes the interaction energy, V, between two charged particles to be composed of two contributions <cit.>

    V = V_ vdw + V_ dl ,

where V_ vdw and V_ dl represent the van der Waals and double-layer contributions, respectively. Van der Waals interaction between two spherical particles with radius, R, can be calculated as <cit.> 

    V_ vdW = -HR/12·1/h,

where H is the Hamaker constant and h is the surface separation distance. For the double-layer interaction Debye-Hückel superposition approximation is used

    V_ dl = 2π Rεε_0ψ_ dl^2 e^-κ h ,

where ε is the dielectric constant, ε_0 is the vacuum permittivity, ψ_ dl is the diffuse-layer potential, and κ is the inverse Debye length. The latter parameter can be calculated as

    κ^2 = 8πℓ_ B N_ A 10^3 I ,

where ℓ_ B = e_0^2/4πεε_0 k_ BT is the Bjerrum length, N_ A is the Avogadro constant, I is the ionic strength expressed in (mol/L), e_0 is the elementary charge, k_ B is the Boltzmann constant, and T is the temperature. Within the Debye-Hückel approximation the charge density and potential of the surface are connected with

    σ = εε_0 κψ_ dl .

At CCC the repulsive energy barrier vanishes, which can be mathematically written as

    V = 0    and   d V/ d h = 0 .

Combining  (<ref>)–(<ref>) with conditions in  (<ref>) permits to calculate the critical coagulation ionic strength (CCIS) <cit.>

    CCIS = 1/8πℓ_ B·( 24π/Heεε_0)^2/3σ^4/3 ,

where e = 2.7182… is the base of the natural logarithm.  (<ref>) can be use to calculate the ionic strength, which corresponds to the CCC. In order to convert the CCIS to the CCC the relation between ionic strength and concentration has to be used

    I = 1/2∑_i c_i z_i^2 ,
 
where c_i is the concentration of all ionic species in the solution and z_i is their valence. The relation between CCIS and the surface charge density given in  (<ref>) has been tested on different types of particles and different electrolytes and gives relatively accurate results <cit.>. This analysis requires to first measure the CCC and then measure the surface charge density at the CCC. Here we propose a simple numerical procedure, which is able to predict the CCC from the measurements of electrophoretic mobility as a function of salt concentration. Electrophoretic mobility, μ, is converted to electrokinetic potential, ζ, via Smoluchowski equation <cit.>

    ζ = μη/εε_0≈ψ_ dl ,

where η is the solvent viscosity. The electrokinetic potential is usually a good approximation for the diffuse layer potential especially if the potentials are relatively low <cit.>. Note that van der Waals and double-layer interaction energies written above correspond to the case of spherical colloids. As long as the dimensions of the particles are larger than the diffuse-layer thickness at the CCC, the Derjaguin approximation is valid and one can use the corresponding effective radii in the equations <cit.>. For small particles of arbitrary shape however some inaccuracies steaming from the use of Derjaguin approximation can be expected.

Based on the equations written above one can construct a simple algorithm to calculate the CCC from electrophoretic mobility vs concentration data. The pseudo-code for this algorithm is given below:




convert mobility to potential     # use  (<ref>)

convert concentration to ionic strength     # use  (<ref>)

convert potential to surface charge     # use  (<ref>)

calculate CCIS from surface charge     # use  (<ref>)

find roots for (CCIS - ionic  strength)= 0     # find where calculated ccis is equal 

     # to experimental ionic strength 



convert the resulting roots from CCIS to CCC     # use  (<ref>)






We have implemented this algorithm in python and the source code is available at <https://github.com/colloidlab/ccc-calculator> <cit.>.




§ RESULTS AND DISCUSSION


The data needed to calculate the critical coagulation concentration (CCC) is an array of electrophoretic mobility values with corresponding array of concentrations, electrolyte composition, viscosity of the solvent, and Hamaker constant for the measured system. All the analyzed data we show in this paper is taken from published work and the respective data sources are cited along the presented data sets.

Let us first explain the procedure for calculating the CCC in more detail. The available electrophoretic mobility data as a function of concentration of added salt is first interpolated by fitting an appropriate interpolating function. Typically a function describing a log-normal distribution works well in majority of the cases. An example of such an interpolation is shown in Fig. <ref>a, where electrophoretic mobility data is shown for positively and negatively charged latex particles in NaCl solutions.



With increasing salt concentration, the electrophoretic mobility for both particles tends to zero as the particle charge is screened by addition of salt. At low salt levels the particle suspensions are stable and the stability ratios are well above 1, see Fig. <ref>b. At increased concentrations one can observe the transition from stable to unstable suspensions and at this transition the stability ratio reaches unity. The concentration where this transition happens is refereed to as the CCC. One can also calculate the ionic strength at this concentration using Eq. <ref>, we refer to this value as critical coagulation ionic strength (CCIS). The interpolation functions shown in Fig. <ref> are used to calculate the CCCs by the procedure described above. The resulting calculated CCCs for amidine and sulfate latex particles are, 225 mM and 460 mM respectively. These values are marked with arrows in Fig. <ref>b and match perfectly with the experimentally measured stability ratios.

In Fig. <ref>cd results for titania nanosheets at two different pH are shown <cit.>. These particles are postively charged at pH 4.0 and negatively charged at pH 10. Again the calculated CCCs marked with arrows in Fig. <ref>d match perfectly the experimental data. The examples shown in Fig. <ref> confirm that our method for calculating CCCs gives good results for positively and negatively charged as well as organic and inorganic particles in the presence of simple monovalent electrolytes.

Let us now look at examples, where monovalent and multivalent ions strongly adsorb to particles surfaces. In the first case, shown in the top panel of Fig. <ref>, negatively charged sulfate latex particles are in contact with 1-hexyl-3-methylimidazolium (HMIM^+) or 1-octyl-3-methylimidazolium (OMIM^+) ions.



The particles are negatively charged at low concentrations and undergo charge-reversal at increased concentrations. This strong adsorption of imidazolium based ions is reflected also in stability of suspensions, see Fig. <ref>b. The particles are stable at low concentration and aggregate fast close to charge-neutralization point after which they re-stabilize due to overcharging and finally become unstable again at high salt levels. This behavior gives rise to three CCCs. Our method can predict three CCCs for each salt and they are marked with red and green arrows, respectively. Similarly, the charge reversal is observed for amidine latex particles in the presence of Fe(CN)_6^3- and Fe(CN)_6^4- multivalent anions, see Fig. <ref>c. Again multiple CCCs are present and can be predicted relatively accurately with the proposed method (see positions of the arrows in Fig. <ref>d). However, the accuracy of predictions for strongly adsorbing ions is in some cases lower. For example, in the case of HMIMCl the third CCC is predicted at ∼ 400 mM while it is observed below 100 mM. These deviations are probably related to non-DLVO interactions, which are not taken into account in our simple method, however they have been shown to be important in similar systems <cit.>. Albeit some inaccuracies, our proposed method for calculating CCCs is still able to quantitatively predict CCCs in most cases.

In the following we focus more on quantitative aspects of the proposed method. In Fig. <ref> calculated versus measured CCCs are shown for seven different types of particles at different conditions.



This analysis enables us to judge the quantitative accuracy of the proposed method. The results for different types of organic and inorganic particles as well as spherical and layered materials are shown. For each system CCCs resulting from the aggregation induced by different ions are shown. All the calculations for a specific system are done with one value for Hamaker constant. These Hamaker constants need to be adjusted because theoretical values of Hamaker constants are either not known or they can overestimate the magnitude of the van der Waals force. The Hamaker constants used for polystyrene latex particles shown in Fig. <ref>a-c are between 2 and 3· 10^-21 J. These values are in agreement with typical values for latex particles in aqueous solutions extracted from direct force measurements <cit.>. Note that however, the measured values are substantially lower as compared to the theoretical value for polystyrene latex calculated from the Liftsitz theory, which is equal to 9.0· 10^-21 J <cit.>. These deviations between the measured and calculated Hamaker constants can be explained by the surface roughness of the particles used in the experiments <cit.>. Therefore, the use of theoretical values of the Hamaker constants for calculation of CCCs leads to underestimation of CCCs, see Fig. <ref>a, b, and f. This procedure of adjusting the Hamaker constants for a given material therefore enables the determination of the effective Hamaker constant for the system at hand.

In the majority of the cases the points shown in Fig. <ref> lay very close to the diagonal line. These results confirm that our simple method is able to quantitatively predict the CCCs for a wide variety of systems. These systems include spherical, non-spherical, and platelet particles composed of different materials in the presence of simple monovalent as well as multivalent, and complex organic ions. In Fig. <ref>f also results for coated titania nanosheets are shown. In the presented case the titania nanosheets were first coated with poly(diallyldimethylammonium chloride) (PDADMAC) or poly(styrene sulfonate) (PSS) polyelectrolytes and then their stability was measured as a function of NaCl concentration at pH 4.0 and pH 10 <cit.>. Again one can quantitatively predict the CCCs of polyelectrolyte coated particles with the Hamaker constant of 1.7· 10^-20 J, which is also used for uncoated titania nanosheets. A more detailed observation of the results shown in Fig. <ref> reveals that in certain cases marked deviations between calculated and measured CCCs exist. The most prominent deviation is visible in Fig. <ref>f for titania nanosheets coated with PDADMAC at pH 10. For this system the calculation underestimates the CCC by almost an order of magnitude. In other words, the experimental system is more stable than predicted. This result is probably a consequence of some additional steric repulsion between two platelets induced by adsorbed polyelectrolytes, which is not part of our DLVO model. Other deviations include titania particles in the presence of KSCN and KNO_3 salts at pH 4.0. In these cases non-DLVO interactions are probably present due to strong specific interactions of anions with the surface. In general, for the systems below the diagonal line, non-DLVO repulsions, while for the systems above this line, non-DLVO attractions are present. Therefore, these deviations can give us a further information about the presence of the non-DLVO forces.

Finally we perform some statistical analysis on the performance of the proposed calculation method. In addition to the systems shown in Fig. <ref> additional data points were collected and are shown in Fig. <ref>a and in Table <ref>.



The data in Fig. <ref>a show that the CCCs can be predicted relatively accurately over more than four orders of magnitude. The distribution of the relative deviation of the calculated CCCs from the measured CCCs is shown Fig. <ref>b. Vast majority of the samples fall between the deviations of [-40 %, 40 %]. Few samples have deviations above 100 % and for these samples additional non-DLVO interactions substantially shift the CCC. Furthermore, the deviations in the positive side are more prominent. This observation can be explained by the fact that in majority of cases non-DLVO interactions, which considerably contribute to the shifting of the CCC, are attractive. Such attractive non-DLVO interactions are typically observed in systems, where ions strongly interact with the surface <cit.>.

By taking the absolute values of the relative deviations one can construct a distribution of deviations shown in Fig. <ref>c. The average deviation for all the sample is 38 %. However, this relatively big average deviation is due to the skewed distribution. There are few samples, for which the deviation is extremely big and these points have a big influence on the average deviation. For more than 92 % of the points the deviation is smaller that 0.75. For the points with the deviations larger than 0.75 the influence of the non-DLVO interactions is substantial and our proposed method is not applicable. If these points are omitted, the average deviation comes down to 25 %. Bearing in mind that the error of the measurements of the CCCs is typically between 10-20 % one can conclude that our proposed method performs surprisingly well.

Let us finally address the conventional wisdom that states that the suspensions lose their stability when the electrokinetic potential of particles reduces below 25 mV <cit.>. As it was already shown by some of us <cit.> this rule is not applicable for suspensions containing multivalent ions. Here we are further testing this simple rule by estimating the CCCs using 25 mV as a threshold of stability. The deviations between CCCs estimated from 25 mV rule and experimental CCCs are shown in Fig. <ref>.



While the 25 mV rule works reasonably well for the CCCs between 0.1 and 1 M, the results are significantly deviating for lower concentrations. This systematic deviation is clearly evident from the fact that the points in Fig. <ref>a follow a slope steeper than one. Furthermore, the histogram of the relative deviations presented in Fig. <ref>b shows that the majority of points have pronounced negative deviations. The average relative deviations in this case is larger than 65 %, which is much worse than our proposed method. What is even more dramatic is the systematic shift to negative deviations. The 25 mV rule is clearly unreliable for low CCCs and therefore its applicability is very limited.




§ CONCLUSIONS


In the present work, experimental electrophoretic mobility values measured in electrolyte solutions were used to calculate CCC of dispersed particles without direct determination of their aggregation rates. The developed method is based on the DLVO theory, i.e., charged particles are expected to form stable dispersions at low electrolyte concentrations due to the stabilization effect by the electrical double-layer force and the samples are destabilized above the CCC by the van der Waals attractions. For systems where the experimental Hamaker constant is not known, this parameter needs to be adjusted in order to get reliable estimations of the CCC. Therefore this method can also be used for determination of the effective Hamaker constant by measuring CCCs with different salts for the same particle system and comparing them to the calculated values.

Statistical analysis of the results and comparison to experimental CCC data determined by independent measurements revealed, that the implemented method provides CCCs, which agree very well with the experimental ones. Note that presence of non-DLVO forces may lead to certain failure of the method, as pointed out with polymer-coated particles, where steric forces provided additional stabilization of the dispersions giving rise to deviation between the calculated and experimental CCC. The method was tested for various systems including particles of different composition (polymer, metal hydroxide, and oxide), structure (sphere, sheet, and lamellar) and size (nano and colloidal particles). In addition, several types of salt solutions (ionic liquid constituents, mono and multivalent electrolytes) were used as dispersion medium. The developed method is able to determine CCC values with a deviation from the experimental values of less than 25 % for majority of samples. On the other hand a simple rule based on the assumption that the CCC occurs when electrokinetic potential reaches 25 mV is much less accurate and has large systematic deviations for small CCC values, which suggests that this assumption is not physically sound. 

 The developed algorithm is provided as free and open source software, allowing other researchers to determine CCCs solely from electrophoretic mobility data. This tool can be of special importance in systems, in which the direct measurement of the CCC is not possible due to non-ideal sample conditions (e.g., high polydispersity, size and concentration of the particles) and also in dispersions containing electrolyte mixtures, as in many industrial and environmental processes.



 










Conceptualization, I.S. and G.T.; methodology, M.G., S.S., I.S. and G.T.; software, G.T.; validation, M.G., S.S., I.S. and G.T.; resources, M.G., S.S., I.S. and G.T.; data curation, M.G., S.S., I.S. and G.T.; writing–original draft preparation, I.S. and G.T.; writing–review and editing, M.G., S.S., I.S. and G.T.; visualization, M.G., S.S., I.S. and G.T.; supervision, I.S. and G.T.; funding acquisition, I.S. and G.T. All authors have read and agreed to the published version of the manuscript.


Financial support by the Swiss National Science Foundation through grant 162420 and National Research, Development and Innovation Office (SNN131558) and the Ministry of Human Capacities of Hungary (20391-3/2018/FEKUSTRAT) is gratefully acknowledged.


The authors declare no conflict of interest. 



no


§ 






Acronyms used in the table: LDH: layered double hydroxide, TNP: titania nanoparticles, TNS: titania nanosheets, PDADMAC: poly(diallyldimethylammonium chloride), PSS: poly(styrene sulfonate), BMPL: 1-butyl-1-methylpyrrolidinium, MIM: 3-methylimidazolium, EMIM: 1-ethyl-3-methylimidazolium, BMIM: 1-butyl-3-methylimidazolium, HMIM: 1-hexyl-3-methylimidazolium, OMIM: 1-octyl-3-methylimidazolium. 



References






yes





