The ionic conductivity and viscous flow data of , , have been collected in a large temperature range, below and above their glass transition temperatures (). A microscopic model is proposed, assuming that the ionic displacement would result from the migration of interstitial positively charged cationic pairs whose concentration is an activated function of temperature. Below , their migration is also an activated mechanism, but a “free volume” would prevail above this temperature. This discontinuity in the migration mechanism justifies a Dienes-Macedo-Litovitz (DML) relationship to be representative of conductivity data above and an Arrhenius law below. According to this model, the enthalpy deduced by the fit of high temperature data using a DML equation would correspond to the charge carrier formation, whose migration enthalpy, below , could be deduced by the difference between the activation energy measured in the Arrhenius domain and the charge carrier formation enthalpy. To reduce the number of adjustable parameters numerical values were physically justified. We also applied a complete test for conductivity below , using the so-called weak electrolyte model, splitting activation enthalpy into formation and migration enthalpies and also explaining the variation of pre-exponential term of conductivity with composition.

1. Introduction

For glass-forming mixtures in the solid or supercooled liquid state, ionic transport due to alkali cations strongly depends on temperature . The variations of the conductivity-temperature product in an Arrhenius representation show two distinct behaviors. At the lowest temperatures, that product follows an activated relationship: where and are constants; is the gas constant.

At higher temperatures, experimental data obey another temperature behavior: where , and are constants. This expression will be deduced in more details below.

Such equation suggests an asymptotic decrease of conductivity towards and is similar to the empirical Vogel-Fulcher-Tammann-Hesse (VFTH) relationship originally established to describe the viscosity-temperature dependence of molten silicates [13]: where (the glass transition temperature) is an empirical constant at which the viscosity diverges, and has the same meaning of the preexponential constant of the Arrhenius expression: where and are constants. The regime of viscosity below follows (4) and is called isostructural (i.e., the viscosity of the glassy state where the structure is frozen). Unfortunately (4) was not used because there is no data in such temperature range.

Some authors have correlated with the Kauzmann temperature. The and parameters are erroneously denominated “activation energies.” In fact, the activation enthalpy for viscous flow in VFTH depends on and in the form , considering a fixed temperature . As described below, both parameters are related to “free volumes.”

In the high temperature domain a hybrid viscosity form associating (3) and (4) has been found: Macedo and Litovitz [4] determined the average activation energy for viscous flow from the probability of forming a hole. Their ideas resulted in an advantageous equation, that combines both the Arrhenius and VFTH equations with the free volume model [5, 6]. Some years before, and following analogous considerations, Dienes [7] obtained a similar expression, but using explicitly the equation below, that includes the VFTH term: where is another constant.

In this work we attribute (5) to Dienes, Macedo and Litovitz (DML). The temperature at which the transport properties change from the Arrhenius law to one governed by (2) or (3) is near . At this temperature, transport and mechanical properties change from solid-like to liquid-like.

Both viscosity and conductivity data in large temperature ranges are not commonly available since different experimental setups are used below and above . Nevertheless some results have been published, considering the same batches [8], something rare in the literature. The ionic conductivity and viscous flow behavior for sodium germanate glasses and supercooled liquids at wide temperature range (from below to over ) are illustrated in Figures 1 and 2, respectively. Lines correspond to adjustable models explained below.

The purpose of this paper is to verify the validity of two microscopic models, trying to find a connection between conduction and viscous flow processes. These models could ascribe a physical meaning to the different ionic conductivity parameters , , and in (2) with viscosity parameters , , and from (5)—considering the same in both cases. We also verified the validity of such obtained parameters from conductivity models above and below . Careful attention was done to predictable and values. From such conductivity values is possible to assess the activation enthalpy formation, that will be considered on application of the weak electrolyte theory for conductivity data below , as presented below. We also will take care of the conductivity dependence with composition, obtaining the resulting charge carrier concentration as well as mobilities for each glass system.

2. Experimental Procedure

Samples were melted in platinum crucible. In fact, it is easier to measure conductivity data below due to different experimental setups: the design of conductivity cells is different below and above the glass transition temperature. Below , glass samples are generally discs of 1 to 3 cm in diameter and 1 to 5 mm thickness. Graphite electrodes were deposited on the two faces in this case study. Electrical measurements were made by d.c. measurements by means of a.c. bridge method using alumina cell with wire electrodes. For conductivity values lower than  S/cm and an applied voltage below 200 mV, the electrode polarization is low, and no significant differences can be noticed between the two techniques. At higher temperatures, over , they are contained in platinum crucibles and two platinum electrodes immersed in the melt. In general the cell is previously calibrated at room temperature with a KCl solution. To avoid electrode polarization the melt resistance is determined in a frequency range in which the resistance is not frequency dependant, generally around 10 kHz.

For obvious technological reasons, the viscosity of a glass or glass-forming melt is an important characteristic justifying the large number of available experimental data. The common viscosity unit is the Poise but the SI unit of viscosity is the with the correspondence Poise. From an experimental point of view, viscosities are currently determined between and . Over a value of the melt has a solid-like behavior and is called a glass, and, below this value, designed as a liquid or supercooled liquid. Few methods are most frequently used to measure viscosity at low or high viscosity. At low viscosity, up to , the rotating cylinder, counterbalanced method, or crucible techniques appear to be the most reliable. Higher viscosities (>) are measured by the rate at which a glass rod elongates under a constant stress or even by penetration/beam-bending methods. The viscosity was measured in this sodium germanate system by counterbalanced method with Pt body and crucible for and by penetration method for higher values.

was obtained from dilatometric measurements with error . More details can be found elsewhere [8]. All adjustments were performed using a Levenberg-Marquardt nonlinear fitting using the Origin software. Table 1 gives the corresponding experimental results, according to theoretical explanation given in this work.

3. Ionic Transport below the Glass Transition Temperature

Below , the usual approach developed for ionic crystals prevails: charge carrier formation may be described by partial dissociation of ionic pair equivalent to the formation of a Frenkel defect, and their migration in a glass is caused by an indirect interstitial mechanism [9], as briefly explained below. The two processes are activated and justify an Arrhenius law for conductivity dependence with (1).

In this temperature range, and for all ionic conductive glasses, the representation of experimental data of the conductivity-temperature product, , in the Arrhenius coordinates results in straight lines whose extrapolation towards infinite temperature converges to corresponding values of the preexponential term . Temperature conductivities of sodium germanate glasses under study spread over 12 orders of magnitude only due to the narrow variation of the conductivity activation energy (), that also depends upon the nature and the concentration of the added alkali cation. Figure 1 shows such conduction results in wide temperature and composition range.

For cationic glasses and regardless of the conduction mechanism, electrical transport can be expressed as the product of three terms, the charge carriers concentration, , their electrical mobility, , and their charge, (by means of Faraday constant):

Since the relative dielectric constant of inorganic silicate glasses is low ionic species are strongly associated. For instance most monovalent cations will be associated with nonbridging oxygen atoms. Such an associated cation should be regarded as in a normal position, hence defining a regular cationic site. Nevertheless, thermal vibrations allow some partial dissociation which leads to the formation of point defects as it is the case in crystalline structures. When a cation leaves its “normal position,” it will be in what could be considered an “interstitial position” and its previously occupied site will become a “vacancy.” This defect formation in the glass structure is equivalent to the formation of a Frenkel type defect in an ionic crystal. In oxide based glasses for which alkali concentration, , is over 10% atom (as in this case), “normal” sites are close enough (3 Å to 5 Å) and a leaving cation is now expected to share a neighboring nonbridging oxygen with another cation. This combination of two cations surrounding a nonbridging oxygen can be described as an “interstitial pair defect.”

The concentration of these interstitial pairs, , may be very small compared with the total concentration of alkali cations, . In other words, glasses are weak electrolytes. In that case, the chemical equilibrium between alkali cations in regular sites and interstitial positions leads to the following relationship: where is the free energy associated with the simultaneous formation of an interstitial pair and a cationic “vacancy,” is the formation enthalpy, and is the formation entropy. The second step is to describe the defect migration caused by the electric field. Thus, the mobility, , of the charge carrier (the interstitial cation) in the electric field is expressed by the general relationship:

In this expression, is the mean distance between two cationic sites, the attempt frequency and is characterized by a migration free energy with its corresponding enthalpy and entropy .

Using (6), (7), and (8), the cationic conductivity in the low temperature range is given by

This relationship reduces to the experimentally observed the Arrhenius law (1) by identification of the experimental value of to

Equation (9) is the usual expression proposed for the intrinsic ionic conductivity-temperature dependence in ionic crystals. The successive steps of a cationic displacement according to this description are illustrated in [10].

Equation (9) interestingly gives the opportunity to calculate a theoretical value for the preexponential term, . Assuming a homogeneous distribution of the alkali cations through the glass, the jump distance, , is the mean distance between two nonbridging oxygen atoms and is related to the total alkali concentration []. For classical oxide glasses, this distance typically varies from 3 Å to 5 Å. The attempt frequency, , may be obtained from the far-infrared absorption which occurs between 500 cm−1 and 100 cm−1 in oxide glasses [11]. Corresponding values of are around 1013 Hz. Assuming a low value for the entropy terms, and , reasonable values of the preexponential term, , lie in the 104-105 (K·S/cm) range in remarkable agreement with most experimental data [12].

From the model above, (10), based on defect crystal theory for ionic transport in solids, the activation enthalpy for conduction is the sum of two enthalpy terms. Contrary to what can be done for ionic conductive materials presenting temperature-dependent intrinsic and extrinsic domains, the migration and formation enthalpies cannot be separated without further experiments and assumptions. In the following paragraphs, we will assume that remains the same at high temperatures and obtain it from viscous flow results considering the same batches.

Fortunately, as far as the author knows, the only theory that can predict the and composition dependences below is the Ravaine-Souquet model [13, 14], presented below.

3.1. Ravaine and Souquet Model (Weak Electrolyte Theory)

The weak electrolyte theory has been extended to glassy electrolytes by Ravaine and Souquet [13, 14] to interpret the important variations in conductivity with composition.

In the classical approaches of ionic transport in glasses, all cations, in the present case all Na+, are supposed to simultaneously migrate. Since the conductivity variations with composition are of several orders of magnitude larger than the total concentration in cations, these variations in conductivity with composition () should be the result of large variations of their mobility with the composition, which cannot be justified from a structural point of view.

Ravaine and Souquet [13, 14] have then supposed that the number of instantaneous charge carriers do not correspond to all the cations present in the glass but only to a small fraction of them who temporally escape from their stable cationic site to participate in conduction process. On a long time scale, that is much over the charge carriers life time, all cations participate to ionic transport, but the number of instantaneous charge carriers is far lower than the total number of cations. This hypothesis is similar to ionic transport by defects as previously described for ionic crystals. In the present case, the concentration of instantaneous charge carriers can be identified to interstitial pair concentration . These interstitial pairs can also be seen as sites of higher energy than the bound lowest energy sites as suggested by Martin and Angell [15] such intermediate sites would temporarily “store” a dissociated cation and make it “available” for ionic conduction.

Thus, Ravaine and Souquet have done the simple hypothesis that the charge carriers mobility could be considered, in a first approximation, as constant with composition in a given glass system and that the large variations in ionic conductivity with composition is the result of the large variations in instantaneous charge carriers concentration with composition. The previous considerations developed remain valid, and the hypothesis of a constant mobility with composition means that all terms in the mobility expression in (9) can be considered as almost constant at constant temperature as a function of .

To estimate, in the system, the number, or at least the variation of instantaneous charge carriers with composition, thermodynamic solutions are applied to this glass-forming system, considered as a “solution” in which the modifier is the solute and the “solvent.” The charge carriers are obtained by partial dissociation of in this low dielectric solvent like it is the case for the weak electrolyte solutions in electrochemistry. Through this dissociation process, the number of instantaneous charge carriers is an exponential function of the partial free energy of Na2O in the mixture, that is, in which is the free energy for the formation of a charge carrier in an arbitrary reference state, and the difference in partial free energy of in this reference state and the studied composition. The factor 2 in the exponent expresses that, similar to ionic crystals, the dissociation process generates two charged species, a charge carrier and its vacant site, both in the same amount.

Finally, the ionic conductivity variations with composition should follow the variations in the sodium partial free energy according to

This correlation has been directly evidenced in different glass-forming systems by simultaneous measurement of the ionic conductivity and alkali partial free energy using potentiometric or calorimetric techniques [13, 16]. Finally, in this approach, the large variations in ionic conductivity would be the result of large variation in with composition () as long as all other energetic terms may be considered as composition independent.

To interpret these large variations in with composition, Ravaine and Souquet [13, 14] and Pradel et al. [17] have developed a statistical model for binary oxide and sulfide systems similar to those already developed for molten silicates [18, 19]. Basically, the mixing enthalpy is deduced from the enthalpy variation in when a bridging oxygen is replaced by two nonbridging oxygens according the equilibrium: xy(13) and the mixing entropy () estimated from the different possible configurations for the oxygen, that is, ionic (O2−), non bridging (–O) and bridging (–O–).

The Gibbs-Duhem equation is then used to calculate the expression of the partial free energy of the alkali oxide as a function of the molar ratio in leading to

The enthalpic term is the sum of two terms, (<0) representing the variations in enthalpy associated to the bond reorganization according to (12)—of the order of some eV [20]—and a corrective term, (>0), expressing an eventual small decrease in this enthalpy as a function of the number of broken oxygen bridges. In other terms, the electronic reorganization implied in (12) modifies the bond strength of the remaining bridging oxygens as evidenced by X-ray structural analysis [21]. A positive value for indicates the progressive decrease in the energy balance of the exothermal reaction with the amount of modifier M2O. This decrease is also evidenced from the reported formation enthalpies of silicates on thermodynamic tables [20].

The dependence of on and, consequently, the number of charge carriers are due to both the entropic and enthalpic terms. In all cases, the important isothermal variations in ionic conductivity with are well represented by the relationship: using (14) to express as a function of . An assessment of the enthalpic term can be done by the research of the best fit with conductivities data as a function of . Such estimation can also be done by fitting variations in activation enthalpy with since the global activation enthalpy obeys the relation: including in the constant term the migration and dissociation enthalpies ( is a constant).

Note that in this model, the defect formation and migration energies are not estimated in absolute value; the model only proposes an interpretation of the variations of isothermal conductivity data with the amount of cation modifier dissolved in the GeO2 solvent.

4. Ionic Transport over the Glass Transition Temperature

The changes in transport properties above involve a cooperative mechanism of the neighboring atoms in addition to the low temperature process, where a new displacement mechanism is observed which may be described by the free volume approach. This mechanism can be understood by a local concentration of free volume distributed throughout the supercooled liquid which leads to the formation of a hole enabling the displacement of the diffusing species. The free volume concept was introduced by Doolittle [5] and further developed by Cohen and Turnbull [6].

At high temperature, that is, above , in addition to the low temperature process, another displacement mechanism is observed, which is schematized in [10]. Local deformations of the macromolecular chains enable the transfer of the defect to a neighboring position. This deformation requires local fluctuations of the free volume allocated to the chain segments [6]. Let us define * as the minimum free volume required in order that the chain movement results in a defect displacement. Hence, the probability, , that this minimum free volume is reached, is where is the mean free volume allocated to a chain segment. Its temperature dependence can be expressed by , where is the thermal expansion coefficient of the free volume, and corresponds to the volume at .

can then be rewritten as

This expression confirms that this second mechanism only appears above the ideal vitreous transition temperature, , at which the free volume disappears. Above this temperature, the cationic displacement may occur either by an activated jump or by an entropic free volume mechanism . As these mechanisms are exclusive, the total probability of a successful displacement reads

Under such conditions, the cationic mobility becomes and the cationic conductivity can be rewritten in the following manner:

The variation of this mobility as a function of temperature depends on the relative values of and . According to (18), a free volume migration mechanism is possible as soon as the temperature reaches . Nevertheless, it is detectable only for higher temperatures because the low temperature activated mechanism prevails below a particular temperature, , which can be seen as the “electrical vitreous temperature.” It is to say that, in the particular temperature range (), remains very small compared to , and the mobility expression (20) is reduced to (8). Hence, the conductivity expression is the same as the activated form proposed in (9), keeping its continuity from up to (in the mathematical sense). On the contrary, over , becomes preponderant and the conductivity expression now reads

This expression of is the product of two exponential terms where the first is representative of an activated mechanism and where the second describes a VFTH behavior. This kind of relationship has already been proposed to improve the viscosity-temperature fits for B2O3 or alkali silicate above their vitreous transition temperature [4]. We may also note that the preexponential term in this relationship contains the same microscopic characteristic parameters as in (9). It means that the limit of the conductivity-temperature product when the temperature tends to infinity does not depend on the nature of the displacement mechanism. Please note that above we assume that all activation enthalpy is due to the defect formation and will assume the same as below . The migration term is mainly due to the “chain segment” motion, that is, due to the formation of holes according to the free volume theory.

5. Numerical Fits and Discussion

For all glass compositions studied within the measured temperature range below , the conductivity was found to obey an Arrhenius behavior, according to (1). Examples are shown in Figure 1, where it is possible to note the increasing in when more sodium is added. For conductivity, the initial procedure was to obtain from the Arrhenius fit below the activation enthalpy . We have reported in Figure 3 the variations of calculated from experimental data using (16) in the range and observed a smooth decrease of between 1.13 to 0.72 eV (see Table 1). On the same graph are superimposed the expected calculated variations in deduced from the Ravaine-Souquet model, with  eV and  eV, giving reasonable results that can be comparable to other systems using the same approach [17]. Please note that this model does not allow to calculate the activation enthalpy in absolute value, but only its variations with composition through the parameter h. As far as the author knows, only recently binary silicate systems had activation enthalpy analysis compared with other theoretical models in near the same wide compositional and conduction ranges [2224].

On Figure 4 are reported the variations of as a function of for the compositions investigated. Despite some scattering, one observes a continuous and monotonous increase with , indicating the low influence of enthalpic terms, as expected. This trend was observed and explained in the same way in lithium silicate systems [25]. The full line superimposed represents the expected variations calculated from the weak electrolyte model. For instance, a variation of preexponential term with composition was experimentally observed in binary alkali silicate, borate, germanate,and tellurite systems [2629], ranging up to two orders of magnitude, with an average value of 104 K·S/cm.

Before analyzing the conductivity data above , we analyzed first the viscous flow behavior and assess reasonable values using DML (, , , ) parameter equation, as shown in Figure 2 and Table 1, using a least square analysis method. As DML and VFTH expressions give the same preexponential , we analyzed the viscous flow behavior using it freely, and the average value found was in agreement with experimental analysis done in dozens of silicate and borate systems [30, 31]. From this table, it was possible to calculate , and , resulting in reasonable values, giving one idea of the chemical bond reorganization above in this system, because ranged between . The DML equation uses four parameters and fits better than VFTH to experimental data, but this should be not surprising because the description of experimental results is most accurate by the equations that use the highest number of adjustable parameters.

About conductivity data above , it was possible to obtain the formation enthalpy term assuming that , as shown in Figure 1, thus easily finding the migration term . From this procedure we assumed that the defect formation does not depend on the nature of the displacement mechanism, that is, the knowledge of the formation and migration enthalpies was supposed to be the same below and above . From Table 1 it is possible to note that decreased with increasing , physically because less enthalpy is necessary to surpass the kinetic barrier, and near following the same trend observed below .

We may note that the respective orders of the magnitude of the two enthalpies ( and ) are reasonably in agreement with those generally observed in crystalline materials presenting extrinsic and intrinsic conductivity domains which allow, in that case, the separation of the two contributions [32]. It was not possible to obtain the entropy values and by adjusting (9); thus its contribution for preexponential term was assumed inserted in preexponential. Then, the Gibbs free energies and could not be determined from appropriated choice of respective enthalpies and entropies of formation and migration, respectively. Indeed, the preexponential values obtained from the Arrhenius fits was near to , in fact very reasonable results, as explained above.

The coherence between the low and high temperature conductivity parameters was obtained from experimental conductivity data using (23) and considering the same from viscosity data (5). We observed that, above , and it was possible to obtain easily assuming fixed, that is, because this term should have the same physical meaning in the full temperature range. All results are summarized in Table 1. The obtained ’s values are lower than considering DML equation, which means that the conduction process needs a lower free volume than the viscous flow. Unfortunately there are no viscosity measurements below to verify if the activation enthalpy for conductivity is lower than that for .

To be on the safe side we also tested (21), that involves the calculation of the and temperature dependent probabilities. As the ratio shows a sharp increasing above in the systems investigated, it confirms the use of the simplified form, (23). For instance, we found that the resulting values from (21) were similar to , which means that the migration enthalpy is related to the free volume available above , following similar procedure done by Caillot et al. [33].

We may now emphasize that (21) and (23) contain the usual VFTH “contribution” (3). The multiplying Arrhenius term takes into account the formation of defects; in fact this is an advantageous equation, that combines both the Arrhenius and VFTH equations with the free volume model. Such combinations of equations with two different exponential temperature dependencies may reflect different aspects of viscous flow mechanisms, for example, separation of the individual diffusing units from their neighbors followed by a translation over certain distances intro an appropriated vacancy [34].

5.1. An Assessment for Charge Carriers Concentration and Mobility above and below the Glass Transition Temperature

Table 1 allows a comparison between the enthalpies for charge carriers formation and migration deduced from the experimental values for activation enthalpy above and below . For the investigated compositions, the charge carriers formation require a formation enthalpy of about 1.10–1.90 eV, their migration enthalpies, a near fixed value close to 0.18 eV.

Using now these experimental values it is thus possible to estimate the relative concentration of the charge carriers in the glassy state using (7) and their mobility by means of (8). Variations of the relative concentration in charge carriers and their mobility are represented in Figures 5 and 6. These results are similar to analysis done in alkali disilicates [35]. Note that the combination of both figures gives the previous lines shown in Figure 1.

Mobility values at room temperature are quite similar to each other (), due to the same obtained values for migration enthalpy, . Such results are similar to the unique mobility measurement of silver done in glass systems, by means of the Hall effect on glasses [36]. The experimental value for mobility obtained from was at 25°C. Thus, our mobility results are in agreement with the hypothesis formulated by the weak electrolyte theory, that is, the independence of the mobility of the ionic carrier with glass composition. Analysis performed in such systems [37] showed similar results and conclusions as done in the present work.

In ionic crystals, the separation of the charge carriers formation enthalpy from the charge carriers migration enthalpy is done by the comparison of the conductivity activation enthalpies in the intrinsic and extrinsic domains [32]. The difference in activation enthalpy in the two domains corresponds to the defect formation enthalpy. In other words, for ionic crystals, it is a discontinuity in the origin of the charge carriers which allows to separate the two enthalpic terms. In the present case, it is not a discontinuity in the origin of the charge carriers formation but a discontinuity in the migration mechanism, over and below the glass transition temperature, which allows to separate the two contributions.

A possible error source on our calculations would be on fixing the jumping distance  Å for all systems under study, following (7). Other possible error concerns with few experimental data above [8]. Finally, the obtained values are quite acceptable, and it is possible to conclude that the proposal determination route to distinguish activation enthalpies of formation and migration is valid.

6. Conclusions

We verified two microscopic models developed to interpret both the ionic conductivity and viscosity variations in large temperature range surrounding the glass transition delimiting the supercooled liquid and vitreous solid domains. Below , the analyzed model shows that the conductivity obeys an Arrhenius law and that the relevant activation enthalpy is a simple function of the migration and formation enthalpies. Above the ideal glass transition temperature a cooperative displacement process, described by the free volume approach, is superposed onto the activated mechanism. Microscopically, this implies that the ionic displacement occurs either because the interstitial pair jumps over an energy barrier separating two adjacent sites or because this energy barrier vanishes consequently to a local rearrangement of surrounding atoms. The latter mechanism results in a conductivity increase appearing for temperatures higher than the transition temperature . A numerical fit of experimental data to viscosity-temperature model enables the determination of the characteristic parameters of both migration and formation processes, that were used to compare with similar model applied to conductivity data in the same temperature range. Moreover, by means of this identification, it becomes possible to calculate separately the migration and formation enthalpies. All determined preexponential parameters for viscosity and conductivity are consistent with the values expected from our current knowledge of ionic crystals and supercooled liquids.

About the conductivity mechanism below , the main interest of the Ravaine-Souquet model is to introduce, in addition, the concepts of solution and statistical thermodynamics to quantitatively justify the variations of the concentration in charge carriers with composition. This thermodynamic approach allows to justify the variations of the preexponential term and activation enthalpy with composition. Moreover, the weak electrolyte theory appears to be confirmed by this work, because the mobility of the ionic charge carriers is independent of their concentration, and the determination route to split activation enthalpies seems to work well, giving reasonable results. This work aims to contribute to future discussions to connect the full temperature dependence of the ionic conductivity from below to above .


The author acknowledges the Brazilian agencies FAPESP, Grants 04/10703-0 & 07/08179-9, and CNPq, contracts 305373/2009-9, 479799/2010-5, 300938/2012-8 and 471546/2012-7. The author offers special thanks to Professor Dr. Jean-Louis Souquet, INP-Grenoble, for explaining the weak electrolyte theory during his stay in Brazil and to Professor Dr. Vladimir Mikhailovich Fokin, Vavilov State Optical Institute/Russia, for constant encouragement.