The Scientific World Journal

Volume 2014 (2014), Article ID 567560, 15 pages

http://dx.doi.org/10.1155/2014/567560

## Translational Dielectric Friction on a Chain of Charged Spheres

Faculty of Sciences of Monastir, UR Electrolytes, University of Monastir, 5000 Monastir, Tunisia

Received 30 August 2013; Accepted 1 October 2013; Published 29 January 2014

Academic Editors: T. Bancewicz and T. Takayanagi

Copyright © 2014 Sondès Boughammoura and Jalel M'halla. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We have proved in details that the dielectric friction remains the principal frictional effect for a stretched polyion modeled as a chain of charged spheres, whereas, in the case of Manning’s model (infinite thread with a continuous distribution of charge), this friction effect is nonexistent. According to this chain model, it is therefore possible to detect by conductivity measurements any transition from a coiled configuration (ellipsoidal model) to a stretched configuration during dilution process. We have also underlined the important interdependence between the dielectric friction and the ionic condensation of the counterions, in order to distinguish between the Ostwald regime and the Manning regime for which the degree of condensation is practically constant in a large range of concentrations.

#### 1. Introduction

It was generally assumed that, when the concentration of a polyelectrolyte is sufficiently low, screen effect due to free counterions of ionic atmosphere is relatively weak so that the intrarepulsion between charged monomers inside each flexible chain of structural length prevents the collapse of the polyion and consequently the stretched configuration becomes the most probable at high dilution. This prevention is enhanced in the case of polyelectrolytes which present a concentration regime satisfying the Ostwald dilution principle, since in this case the degree of ionic condensation () of counterions on charged monomers (partial neutralization of the structural charge of the polyion) decreases with dilution [1–4]. It follows that continuous transition from stretched (or rod-like) conformation to coiled shape could be observed by increasing the concentration of counterions: , or by decreasing the permittivity of the solvent. Indeed, since the concentration effect changes the apparent charge, , the shape, and the size of polyions and since, in their turn, these parameters govern the different frictional processes, it therefore results in more or less sharp variation of the mobility of the polyions with the concentration. In previous studies [3, 5], the various friction effects on polyions have been classified to five types: (a) the hydrodynamic friction which depends on the viscosity of the solvent, , the size and the shape of the polyion, and it is quantitatively evaluated by the value of the hydrodynamic equivalent conductivity of the polyion, ; (b) the electrophoretic friction which expresses the hydrodynamic friction on the ionic atmosphere of the polyion, and it is quantified by the electrophoretic conductibility increment, ; (c) the ionic friction (or ionic relaxation effect) due to the local field caused by the polarization of the ionic atmosphere by the external field during its relaxation so that the frictional force acting on the polyion is equal to , the intensity of this effect is evaluated by the ionic friction coefficient ; (d) the translational dielectric friction effect [3, 6] due to the perturbation of the polarization of solvent molecules around the moving polyion. This effect is evaluated by the coefficient , where is the local dielectric field. Note that is proportional to the square power of the degree of dissociation: .

The relative importance of each friction contribution depends on the concentration regime and on the conformation of the polyion. However, it is possible to express formally, in the general case, the equivalent conductivity of the polyion in terms of the different friction contributions as follows [3, 6]: For ellipsoidal polyions (see Figure 1) of focuses and , of interfocuses distance: , of minor axis , and major axis , obeying Ostwald regime, so that , when , we have generalized the Debye-Onsager theory concerning simple electrolytes for the calculation of the electrophoretic increment and the ionic friction coefficient [6]. On the other hand, we have generalized the Boyd-Zwanzig’s approach [7] concerning simple spherical ions, to the case of ellipsoidal polyions [6], in order to evaluate the dielectric friction coefficient at infinite dilution . The different contributions have been expressed in terms of the degree of dissociation of counterions , the geometric parameters of the ellipsoidal polyion and , the effective radius of the counterions , the charge numbers of the polyion and the counterions respectively, , , and also in terms of the minor axis of the ellipsoidal ionic atmosphere surrounding the polyion, , (see Figure 1) which depends on the concentration of free counterions and therefore on the Debye-MSA [8–11] screen parameters and (in ): is the Avogadro number and is the Bjerrum length, so that the hydrodynamic and the electrophoretic contributions are respectively given by is the Faraday, and are, respectively, the mean radius of the polyion and of its ionic atmosphere defined as follows: where is the so-called “generating function,” characterizing the conformation of the polyion (and its ionic atmosphere) which can vary from the spherical shape () to the cylindrical configuration () [3, 6]: Note that measures also the electrostatic capacitance of the ellipsoidal polyion in c.g.s.u.e units [6] and measures also the inverse of the electrostatic capacitance of the ellipsoidal capacitor constituted by the polyion and its ionic atmosphere “Gouy capacitance” (see Figure 1).

The ionic friction due to the perturbation of the ionic atmosphere during its relaxation is expressed in terms of the coefficient as follows: Note that in the limiting case of Manning’s model (i.e., , ), the ionic friction coefficient .

It follows according to (2), (3), (5), and (8) that the contributions and , relative to ellipsoidal polyions, vanish at infinite dilution , , , so that
The general expression of the dielectric friction coefficient is [6]
and are, respectively, the static and the high-frequency dielectric constants of the solvent. The apparent radius is a function of the eccentricity of the polyion so that for and for [6]. For spherical polyion (), , and (10) becomes in this case identical to Boyd-Zwanzig’s equation. Note that the dielectric friction coefficient at infinite dilution is in general of the order of few percent and *decreases with elongation * of the polyion so that the conductibility of the polyion is essentially governed by the ionic condensation process and the hydrodynamic friction .

In the case of stretched polyions obeying Manning’s regime , is negligible, whereas and remain relatively important and especially quasi-independent on the concentration so that .

However, in a recent work we have observed a *sharp increase* of the dielectric friction due to the* elongation* of some polyions during the *dilution* process [12–16]. We have concluded that (10) predicting the *decrease* of the dielectric friction coefficient with the elongation is only valid for stretched polyions characterized by a *continuous* linear distribution of their apparent charge .

In fact, the real structure of a stretched polyion is rather similar to a chain of charged spheres than to a rod uniformly charged. In other terms, the stretched polyion could be modeled as a succession of charged rigid spherical monomers “pearls” or “groups” of radius and ± charge. The charge distribution of the polyion is therefore discontinuous “necklace.” The distance separating two successive groups is equal to . Note that this model converges toward the rod-like Manning’s model when the ratio . On the other hand, for coiled conformation, the polyion could be treated as a rigid ellipsoid uniformly charged of and focuses, length, minor axis, and () charge.

The essential difference between the chain configuration and the ellipsoidal conformation concerns the importance of the dielectric friction. In the first case, the total dielectric friction results from the superposition of the various dielectric frictions acting on the successive charged groups so that the corresponding dielectric friction coefficient is proportional to the ratio , whereas, in the second case, the friction coefficient is proportional to ()^{2}. Consequently, an important variation of the dielectric friction is expected during this conformation transition.

The main objective of the present work is therefore to explain in details why the transition from coiled configuration to completely stretched chain is accompanied by a sharp increase of the dielectric friction on the moving polyion, whereas such variation is undetectable according to the model of Manning.

In order to achieve this objective progressively and completely we have organized the rest of the paper as follows. In Section 2 we present an new simplified (heuristic) derivation of the expression of the dielectric frictional force acting on a moving spherical charge, different from those of Zwanzig [7], Hubbard and Douglas [17], and Hubbard and Onsager [18] which has the advantage of underlining the physical significance of the process of dielectric friction. In Section 3 we give a brief recall of the general expression of according to the time-correlation function formalism [6] in terms of the displacement field created by the polyion, the memory function () relating the polarization around the polyion to , and the key integral allowing the direct calculation of via and . In Section 4 we apply the previous approach to calculate in details the explicit expression of and therefore the frictional force in the case of a chain of charged spheres. In Section 5 we discuss the coupling between the dielectric friction effect and the ionic condensation processes on the basis of a generalization of the Fuoss’s approach. Section 6 gives the explicit expression of the variation of the conductibility of a polyion with dilution. Section 7 concludes with a brief discussion of the results using some experiments.

#### 2. A Simplified Derivation of the Expression of the Dielectric Frictional Force on a Moving Charged Rigid Sphere

Figure 2 represents a moving rigid sphere along the -axis, of radius, charge, and velocity . represents its center at time “” and its center at a previous instant so that the distance is equal to . The moving charge creates at each point of the dielectric medium, a time dependent electric displacement which in its turn polarizes the solvent molecules. We design by the induced polarization at the surface of the sphere and by the induced polarization at the surface of the sphere at . Now, according to the theory of dielectric mediums, the superficial charge density of the dielectric molecules at the interface between the solvent and the charged sphere at time is related to the orthogonal component of the induced polarization at the surface of the sphere as follows: On the same way the charge density is equal to On the other hand, according to time-correlation function formalism, can be expressed via a temporal convolution integral, as a sum of linear responses to the successive orthogonal components of at the different anterior instants : According to Figure 2, is the associated after-effect function representing both electronic relaxation and rotational diffusion of dipolar molecules [7]: is the relaxation time of the solvent molecules and and are the static and high-frequency dielectric constants of the solvent. For water at 25°C, and . Notice that as the ratio , the exponential function exp() decreases sharply with time and therefore the temporal correlation between and vanishes rapidly when .

The expression of is more ordinary because it is induced by the displacement field which is constant for all so that Slowly moving particle is characterized by the condition , so that during the period of correlation . Now, according to (15), cos() is always 0; it follows therefore from (13) that for , for all cos(). In the case of speedy particle so that , (13) and (14) imply that for all 2 and for all 2. Note also that for all .

In both cases, the dielectric superficial charges and present an axial symmetry around the -axis, and consequently they are created at the center of the moving sphere, and, at time , a reacting dielectric relaxation field is directed along the -axis. According to the principle of superposition, its component can be expressed as a sum of two local fields, respectively, and : The relation between the component and is obviously The integration is done over the surface of the moving sphere at time , that is, for .

On the other hand, the expression of the component due to is more subtle and depends on the movement of the charge. Indeed, if , is inside the sphere of center and of radius; therefore, the image of the charge (which is a punctual charge concentrated in ) is situated inside the cavity of center and of radius , and consequently: In contrast, for , is outside the cavity, and therefore Now, the component of the resulting dielectric frictional force acting on the moving charged sphere due to the superficial charges and is therefore equal to Calculation of according to the different above equations leads to the following expressions which depend on the condition imposed on the velocity of the polyion.

(a) For , Therefore, , and according to (19) and (22), The denominator must be always , so that when cos, , because . If we substitute the function cos() by the new variable , we get Now the integral can be decomposed as a sum of two integrals: The integration over involves elementary functions so that the above integrals exist in the literature. Since for and , we have , the integration leads to the following exact results: so that Consequently, (25) can be simplified as follows:

(b) For .

In this case, when we have , and integration leads to the following exact results: with so that Consequently, (31) can be simplified as follows: On the other hand, the use of (21), (22), and (34) leads to And therefore Finally, the general expression of the resulting dielectric force: which is valid for and for is This result is identical to original Zwanzig’s result [7].

In the limit of low velocity that is , that is, for a slow particle
it is interesting to note that we can physically interpret this last result according to the linear response theory combined to a dimensional analysis as follows. Indeed, first, frictional effect is absent for immobile particle (); therefore, we expect that must be proportional to ; second, there is no dielectric friction if the depolarization of the solvent molecules around the particle is instantaneous (); consequently, is also proportional to the relaxation time ; on the other hand, there is no relaxation effect if the solvent is dielectrically saturated (); consequently, must be proportional to (); finally, is an electric force acting on a sphere of charge and radius , and it is in principle proportional to ()^{2}()^{−2}; now as is also proportional to the length , it follows that must be proportional to ()^{−3} so that it sharply decreases with the size of the particle.

More rigorous derivations of taking into account hydrodynamic motion of the solvent were performed successively by Zwanzig [7], Hubbard and Douglas [17], and Hubbard and Onsager [18]. More recent approach taking into account molecular correlation between solvent molecules and moving particle was presented by Wolynes [19]. However, if the charged sphere is assumed to be a conductor of large radius , then hydrodynamic effects become small and all theories converge to Zwanzig’s original result.

Note finally that (38) has been generalized to the case of ellipsoidal polyion of charge characterized by an apparent radius which depends on the eccentricity so that when [6]; the result is

#### 3. General Expression of the Dielectric Frictional Force

The previous expression giving the dielectric frictional force acting on an ellipsoidal (or spherical) polyion of center , and focuses, minor axis, major axis, length (), and charge, can be found by applying the following general method [6, 14].

As explained above, the *time dependent* polarization at a position around the polyion ( point; see Figure 3) is induced by the displacement field due to the () charge of the moving polyion. In its turn, this induced polarization creates a reacting dielectric relaxation field in and therefore exerts a dielectric frictional force () back on the polyion directed along the axis. Indeed, according to the dielectric theory, the local charge density of the dielectric molecules into the element of volume at is equal to so that
with
The explicit development of the previous relations according to adequate boundary conditions enables us to reduce the expression of the dielectric frictional force to the following general form which is in fact valid for a charged macroion of any shape:
is the so-called key integral [6] related to the components , , and of as follows:
According to Figure 3, , and are the Cartesian coordinates of the vector radius of module ; is the vector radius of module . is the distance covered by the center of the moving polyion during the time with the velocity , and .

The above triple integration is effectuated over the whole volume of the dielectric medium; consequently, the main mathematical difficulty comes from the limiting conditions imposed on the interdependent lower boundaries of integration () defining the *finite* surface region surrounding the volume from which the dielectric medium must be excluded.

Demonstration of (42) and (43) was first achieved by Zwanzig in its original work [7] in the case of a moving spherical charge and then by the authors of [6] in the case of a moving ellipsoidal polyion.

As indication, we give below the explicit formulas corresponding to the components , , and , created by an ellipsoidal polyion: The parameters , and are defined in Figure 3; and is the linear charge density of the polyion. The introduction of the explicit expressions of , , and , into (43) leads to the following expression of the key integral: Therefore, integration of according to (42) leads to the same results obtained in paragraph 2 ((38) and (39)).

In the next paragraph, we will generalize this approach based on the notion of the key integral , in order to calculate in details the dielectric frictional force acting on a moving chain of charged spheres.

#### 4. Expression of the Dielectric Frictional Force on a Moving Chain of Charged Spheres

Figure 4 represents a moving polyion as a succession of identical charged rigid spheres (or “groups”) of centers of radius and charge. The distance of separation between two successive charged spheres is defined by .

indicates the fixed Cartesian reference frame so that the chain moves along with a velocity . The relative position of a charged sphere “” of center is defined compared to an *arbitrary* origin which coincides with the center of an unspecified sphere “” so that
and are, respectively, the coordinates of and . Now, if , we can write
If , then
with the obvious condition
is an arbitrary point inside the dielectric medium and and are the distances of separation between and, respectively, and at time . On the other hand, is the distance of separation between and at time (i.e., the center of the red sphere represented in Figure 4), so that is equal to the distance of separation between the center of the sphere “” at time and its center at .

According to Figure 4, in which , the distance is equal to For simplification we have used the new parameters Now, the sphere “” of charge creates in the point at an electrical displacement field Therefore, the expressions of the components , and are It is also important to note that the denominator which appears in (53)–(55) must be strictly positive.

According to the principle of superposition, the total displacement field due to total charge of the moving polyion in at is equal to the sum of the local fields created by the successive spherical charges :
It is clear that the components of created by the reference charge in at are calculated from (52)–(55) by imposing . We can also develop the equation of superposition given by (52) as follows:
Now, in order to simplify calculations of the dielectric frictional force according to the general method exposed in Section 3 on the basis of (42) and (43), we will define for each number (or ) a corresponding *cross* key integral as follows:
with . The triple integration in this equation is done over the whole volume except the volume of the spherical charge . The explicit expression of this integral in terms of the spherical coordinates and is obtained after replacing the volume element of the dielectric by . The result is
The insertion of the expressions of , , and given by (53)–(55) into the above relations leads to
We can decompose the second integral over the variable as a sum of three elementary integrals:
At this stage, we must distinguish two conditions.

(a) , so that and ; thus, Therefore

(b) , so that and ; thus, Therefore This last equation means that polarized solvent molecules inside the sphere of radius do not participate in the dielectric friction. (Note that we have obtained a similar result in the case of an ellipsoidal polyion of center and length [6]. All occur as if during its translational motion, the polyion turns around its center so that the solvent molecules inside the sphere of center and radius do not take part in the process of dielectric friction).

Finally, by introducing (63) and (65) into (60) and after simple integration we obtain the following final expression of the cross key integral :
It is important to underline that represents a *cross * effect. Indeed, it is proportional to the dielectric frictional force acting on the reference charge , due to the polarizationof the solvent by the field created by the spherical charge “”:
In the case of a spherical charge “” characterized by the condition , we define in the same way a key integral related to its corresponding force :
with the condition .

In contrast, the *self*-key integral is proportional to the *self*-dielectric frictional force acting on the reference charge which is induced by its *own* charge . is therefore Zwanzig’s integral relative to a moving sphere of charge and radius; its expression is therefore analogous to (45):
Consequently, the total dielectric force acting on the reference spherical charge which is induced by the different spherical charges of the chain is given by
The resulting dielectric force acting on the polyion is obtained via the total key integral by a summation over the reference spherical charges :
Recall that and are related by the condition . According to (66), (68), and (69), the key integrals , , and are, respectively, proportional to , , and ; one can therefore regroup the different terms of the above sums over and in order to transform them into a sum of couples (, ) of the form ; each couple has a “degeneracy” equal to (). Consider
Table 1 gives as an example the matrix of the various terms of the key integral of a chain of eight () charged spheres.

In the approximation of the linear response, must be proportional to () with the condition ; its expression can therefore be simplified as follows:
As the sums and converge rapidly, respectively, to
the final result for large is therefore
We now define the “*interference factor*” as follows:
The final expression of the dielectric force acting on the charged chain of the polyion is obtained according to (76), (71), and (16) and by using the relation
Note that when the distance of separation between two successive charges is large compared to the radius of the charged groups (i.e., and ), the mutual influence between charged spheres vanishes, so that the solvent molecules surrounding each sphere “” are polarized essentially by the displacement field caused by its own charge (*no interference*). Consequently, the dielectric force: undergone by each group “” is therefore reduced to Zwanzig’s force , and the total force acting on the polyion is obtained by superposition of these Zwanzig’s forces:
It is also important to realize that in all cases , which means that the total force is always opposed by the movement of the polyion; however, this braking effect is attenuated by the effect of interference, as it is shown by the ratio
It is now possible to derive the expression of the dielectric friction coefficient defined previously by . is the local field caused by the polarization of the solvent molecules around the moving polyion and it is related to the frictional force as follows:
() is the effective (apparent) electric charge of the polyion equal to
is the effective charge of any charged sphere of the chain. On the other hand, the velocity of the polyion is related to its electrical mobility and to the external field (directed along the -axis), according to
Introduction of (83) into (78)–(82) leads to the general expression of the dielectric friction coefficient :

#### 5. Coupling between the Dielectric Friction Effect and the Ionic Condensation Processes

Because of the ionic condensation effect, the charged groups of the polyion are partially neutralized by the counterions “” of charge and effective radius. The degree of condensation () depends on the configuration of the polyion. Recall that according to Manning’s rod-like limiting model () [5], the degree of ionic condensation () is independent of the counterions concentration so that ; where is the Bjerrum length, is the Boltzmann constant. In fact, experimental conductivity results are not in general in conformity with Manning’s theory [1–5]. In particular, we have proved that for more realistic polyelectrolyte models the degree of dissociation obeys in general Ostwald’s principle of dilution and consequently when . The degree was calculated on the basis of the two-state model [6] (double layer) proposed by Dobrynin and Rubinstein [20] (see Figure 1) and by the generalization of the theory of ionic association of Fuoss [8, 9], so that for spherical polyion (), and for an infinite chain (). The result of an ellipsoidal polyion () of concentration and of volume is [6, 12] The degree () of ionic condensation on a charged chain is obtained from this relation by replacing by the structural length and by . On the other hand, as , the charge is thus equal to As indicated in Section 1, the electrical mobility of the polyion is a complex function of , and its expression can be written in the following general form: The insertion of these last equations into the expression of the coefficient of dielectric friction leads after linearization to Now, if we replace the relaxation time by its Debye’s expression: , with , the general explicit expression of becomes This equation shows that the coefficient of dielectric friction corresponding to a chain of charged sphere is as expected and in a first approximation proportional to number of charged groups and to the square power of the degree of dissociation .

Comparison between the given by (89) and given by (10) leads to This ratio explains the important increase of the dielectric friction during the transition from an ellipsoid configuration toward the stretched chain configuration. Indeed, as and , , and therefore Note finally that according to (89), increases with dilution () in the case of polyions obeying Oswald’s regime () and it remains sensibly constant for polyions obeying Manning’s regime for which is constant.

#### 6. Variation of the Conductibility of a Polyion with Dilution

As indicated in the previous sections, the dilution increases the degree of ionic dissociation and the mean radius of the ionic atmosphere and modifies the configuration of the polyion. However, for polyions obeying Manning’s regime, the degree “” remains sensibly constant, so that both frictional coefficients and are also constants. In this case, it is therefore easier to test the importance of the dielectric friction effect on stretched polyions by comparing the conductibility of such polyelectrolytes to the one calculated in absence of dielectric friction. Recall that according to limiting Manning’s model the dielectric friction is nonexistent, ; indeed, the moving polyion is assumed to be an infinite thread with a continuous distribution of charge, so that it seems immobile and therefore the polarization of its surrounding solvent molecules is not perturbed by the movement.

As an illustration, Figure 5 compares the variation (with the total concentration ) of the experimental conductibility and (blue points) of sodium chondroitin sulfate in water at 25°C to the following three theoretical conductibility types: , and . The first one is calculated in absence of dielectric friction. The second is calculated by taking into account the dielectric friction but in absence of interference. The last one is calculated with dielectric friction in presence of interference.

The structural characteristics of this polyion are [6, 7] as follows. is the structural charge number. Å is the structural length. Å is the charge-to-charge distance. Å is the cylindrical radius of the polyion chain.

Sodium chondroitin sulfate is one of peculiar polyelectrolytes for which the behavior of ionic condensation in aqueous solution is compatible to the model of Manning [13, 14], so that remains practically constant equal to .

The equivalent conductivity of this polyelectrolyte is therefore equal to and depend on the concentration of the free counterions because of the braking effects due to the ionic atmosphere, giving rise to the electrophoretic effect () and to the ionic relaxation effect (). The explicit expressions of and in terms of , , , and are the following: is the Faraday. The screen parameter is given by (3) Calculations are effectuated with , cp, Å, ,