13 PitzerKS JSC 13 265 1998
13 PitzerKS JSC 13 265 1998
13 PitzerKS JSC 13 265 1998
4, 1999
While the original ion-interaction (Pitzer) equations of 1973 were adequate for
many electrolytes to the limit of solubility, additional terms are needed for some
systems of large solubility. A simple pattern of ionic-strength dependence is
proposed for third, fourth, and higher virial coefficients. It is found to be very
effective in representing the complex behavior of CaCl2(aq) as well as that of
MgCl2(aq) at 25°C. Equations without ion association are presented for each
system that are valid for the full range to 11.0 and 5.9 mol-kg -1 , respectively,
as well as simpler equations for limited molality ranges.
1. INTRODUCTION
The primary advance in papers 1 and 2 of this series(1,2) arose from the
recognition and accurate representation of the ionic-strength dependence of
1
Department of Chemistry and Lawrence Berkeley National Laboratory, University of Califor-
nia, Berkeley, California 94720-1460. Current address for Peiming Wang: OLI Systems, Inc.,
108 American Road, Morris Plains, New Jersey 07950.
2
Geosciences and Environmental Technologies, Earth and Environmental Sciences Directorate,
Lawrence Livermore National Laboratory, University of California, Livermore, California 94550.
3
School of Environmental Sciences, University of East Anglia, Norwich NR4 7TJ, U.K.
4
Deceased December 26, 1997.
f
We dedicate this paper to the memory of Kenneth S. Pitzer in recognition of his many
invaluable contributions to solution chemistry.
265
0095-9782/99/0400-0265$16.00/0 C 1999 Plenum Publishing Corporation
266 Pitzer, Wang, Rard, and Clegg
the second virial coefficient in the low molality range. Basic theory for
electrolyte solutions(3) provides for ionic-strength dependence of the various
virial coefficients. However, for most systems and with accurate second
virials, the third virial coefficient is so small that any ionic-strength depen-
dence is negligible. Thus, the ion-interaction (Pitzer) equations were presented
with third virial coefficients independent of ionic strength.(1,2,4,5) They are,
of course, dependent on pressure and temperature. This was found satisfactory
for mixed as well as pure electrolytes(4) and even extending to solid saturation
molalities, in most cases, in the extensive investigations of Harvie and Weare(6)
and of Filippov and associates.(7) Subsequently, fourth and higher virial coeffi-
cients were added for electrolytes extending to very high molality. Examples
are CdCl2, ZnCl2, LiCl, and CaCl2. For the CdCl2, ZnCl2, and LiCl, these
treatments were quite satisfactory,(8,9) but for the important case of CaCl2,
there are now very accurate data(10) and the 1985 fit of Ananthaswamy and
Atkinson (11) with virial coefficients through the sixth order, although good, can
now be improved. This paper presents and tests for CaC!2(aq) and MgCI2(aq) a
simple formulation of ionic-strength dependent terms extending to any order.
Recently, Archer(12) introduced an ionic-strength term at the third virial
level for NaBr, which is soluble to 9 mol-kg - 1 at 25°C. He also found that
this addition improved the fit to the very accurate data for NaCl(aq).(13)
Archer's extension with an ionic-strength dependent third virial term was
generalized by Clegg, Rard, and Pitzer(14) to self-associated electrolytes and to
electrolyte mixtures. They applied these equations to solutions of H2SO4(aq),
which contain equilibrium mixtures of H+, SO2-, and HSO4-. Their model provides
a very accurate representation of the available osmotic and activity coefficients,
enthalpies of dilution, and heat capacities to 6.1 mol-kg -1 . The degrees of dissocia-
tion of the HSO4 ion calculated with this model are in good agreement with
experimental values from spectroscopic measurements.
Rard and Clegg(10) used similar self-association models for CaCl2(aq), which
include CaCl+ ion pairs. There is ample evidence for ion pair formation above
about 5 mol-kg-1 in these solutions,(15) but neutron diffraction data(16) show no
significant association at 4.5 mol-kg -1 . Rard and Clegg(10) obtained an accurate
model for CaCl2(aq) valid to 8.0 mol-kg -1 , and, by including additional virial
terms, a second model valid over the full molality range to 11.0 mol-kg -1 .
However, the association constant K°(CaCl+) is poorly known at 25°C, in contrast
to the more precisely known value of K(HSO4). This uncertainty for K°(CaCl+)
gives rise to large uncertainties in the predicted degree of dissociation of CaCl+
and the corresponding ionic activity coefficients for CaCl2(aq) solutions, which
then depend strongly on the assumed value of K°(CaCl+) and on whether higher-
order virial terms are included in the model.
Both models presented by Rard and Clegg yield very accurate values
for the stoichiometric quantities c and G±, and thus their equations are
Ionic Strength Dependence 267
with possible additional terms in B(3), a3, etc.; in addition, B(2) is often zero.
If we define
Then the second virial expression for the excess Gibbs energy is
For the natural logarithm of the mean activity coefficient, one has
Thus, for the fourth virial, if there are three terms, one has for the excess
Gibbs energy and with xDO = 0
with
A complete expression through the sixth virial for the osmotic coefficient
is then
At the third virial level, the only difference from the Archer term for
the osmotic coefficient is the use of the first power of I instead of I1/2 in xCL.
However, for the Gibbs energy and the activity coefficient, the expressions
corresponding to Eqs. (5) and (17b) are more complex for the Archer term
and progressively much more complex for the fourth or higher virials.
At each level, the ionic-strength-independent term, B(0), C(0), D(0), etc.,
is obtained with a = 0, x = 0. It is only at the second virial level that there
is useful guidance from basic theory for the A's. The other nonzero a values
are chosen for optimum fit to experimental data. While purely empirical, the
higher order terms with nonzero a are very useful in fitting experimental
data since their contribution is relatively localized in a limited molality range.
Comparison with Eqs. (7), (8), and (10) then yields, for the D and E terms
Then, for mixtures, terms of high order such as O(Na, Cd, 3Cl) are added if
Ionic Strength Dependence 271
Since the basic objective was equations that are not only accurate, but
also convenient, for use for either pure CaCl2 or for mixed solutions, several
equations of increasing complexity are presented for increasing ranges of
molality. We consider first the complex equation, but still without ion associa-
tion, for the full range from 0 to 11 mol-kg -1 (i.e., 3.5 mol-kg -1 above
saturation). It extends to the fifth virial coefficient with nine terms in all.
The parameters are given in Table I, along with those for the simpler equations.
The Debye-Huckel parameter AP = 0.391476 kg l / 2 -mol - l / 2 is based on the
dielectric constant equation of Archer and Wang.(l9) There are five ionic
strength-dependent terms with a values selected after various trial calcula-
tions. The B, C, D, and E values are obtained by the least-squares method.
For mixtures the molality of CaCl2 seldom exceeds 7 mol-kg -1 and the
ionic strength is less than 21 mol-kg -1 . In addition, the interpretation of the
fifth virial is ambiguous with respect to CaCl4 and Ca2Cl3 interactions.
Hence, an equation limited to the fourth virial is desirable; we present one
valid to 7 mol-kg -1 with seven terms. The fourth virial is certainly dominated
by CaCl3 interactions and can be treated as a w(Ca, 3Cl) term in the Filippov
system described above. This equation, valid to 7 mol-kg - 1 , also suffices for
Table I. Parameter Values for CaCl(aq) in Eqs. (14-17) with 3, 5, 7, and 9 Ion-
Interaction Parameters
Parameter Units np = 9a np = 7a np = 5a np = 3b
B(0) kg-mol-1 0.569081 0.464899 0.446867 0.31588
B(1)
(2)
kg-mol -1 1.59385 1.55840 1.57819 1.6140
B kg-mol -1 0 0 0 0
ABI kg1/2-mol1/2 3.0 2.5 2.45 2.0
C(0) kg 2 -mol -2 -0.0556802 0.00459539 -0.00663855 -0.00017103
C(2)(1) kg2- mol - 2 0.215755 -0.00941176 -0.0151049
C kg 2 -mol -2 -0.317714 -0.0729924 -0.0421769
ACI kg-mol -1 0.15 0.15 0.15
AC2 kg-mol-1 0.25 0.25 0.25
D(()) kg-mol-3 0.00282839 -0.000745381
D(1) kg3-mol-3 0.0170624 -0.00241580
D(2) kg 3 -mol -3 -0.0244910
ADI kg3/2-mol-3/2 0.030 0.030
AD2 kg3/2-mol-3/2 0.035
(0)
E kg4-mol-4 -0.0000509632
mmax mol-kg -1 11.0 7.0 5.2 2.5
Sf 0.0033 0.0026 0.0022
Number of points 430 355 293
a
This study.
b
Pitzer and Mayorga (Ref. 2).
c
Standard deviation of the fit.
Ionic Strength Dependence 273
Fig. 1. Comparison of (a) osmotic coefficients and (b) activity coefficients of CaCI2(aq) from
Rard and Clegg (Ref. 10) with those calculated using equations from this study (Eqs. 14 and
16; np = 5, 7, 9) and from Pitzer and Mayorga (Ref. 2; np = 3). np = 9; np =
7; np = 5; np = 3.
274 Pitzer, Wang, Rard, and Clegg
Table II. Comparison of the Osmotic Coefficients of CaCl2(aq) from Table 11 of Rard
and Clegg with Those Calculated Using Eq. (14) and Ion-Interaction Parameters in
Table I
ma Rard and Clegg np = 9b np = 7b np = 5b np = 3c
0.001 0.9623 0.9623 0.9622 0.9622 0.9621
0.002 0.9493 0.9494 0.9492 0.9492 0.9490
0.005 0.9275 0.9275 0.9271 0.9272 0.9269
0.01 0.9078 0.9079 0.9074 0.9074 0.9072
0.02 0.8871 0.8871 0.8866 0.8867 0.8866
0.05 0.8628 0.8627 0.8627 0.8629 0.8636
0.10 0.8527 0.8527 0.8533 0.8535 0.8553
0.20 0.8588 0.8589 0.8593 0.8594 0.8613
0.30 0.8748 0.8753 0.8751 0.8750 0.8758
0.40 0.8944 0.8957 0.8950 0.8948 0.8941
0.50 0.9163 0.9181 0.9172 0.9170 0.9150
0.60 0.9401 0.9419 0.9411 0.9410 0.9382
0.70 0.9655 0.9668 0.9664 0.9664 0.9633
0.80 0.9924 0.9928 0.9929 0.9929 0.9900
0.90 1.0205 1.0200 1.0204 1.0205 1.0180
1.00 1.0495 1.0482 1.0489 1.0490 1.0474
1.20 1.1100 1.1078 1.1087 1.1088 1.1092
1.40 1.1733 1.1715 1.1721 1.1721 1.1745
1.60 1.2393 1.2387 1.2387 1.2386 1.2426
1.80 1.3079 1.3089 1.3084 1.3082 1.3129
2.00 1.3793 1.3816 1.3807 1.3805 1.3850
2.25 1.4725 1.4754 1.4744 1.4743 1.4772
2.50 1.5699 1.5717 1.5712 1.5712 1.5711
2.75 1.6705 1.6702 1.6703 1.6706
3.00 1.7731 1.7705 1.7713 1.7718
3.25 1.8759 1.8723 1.8736 1.8740
3.50 1.9778 1.9751 1.9765 1.9767
3.75 2.0788 2.0785 2.0795 2.0793
4.00 2.1796 2.1818 2.1820 2.1812
4.25 2.2805 2.2841 2.2834 2.2819
4.50 2.3808 2.3845 2.3828 2.3811
4.75 2.4792 2.4820 2.4795 2.4784
5.00 2.5743 2.5753 2.5727 2.5735
5.25 2.6643 2.6635 2.6615 2.6661
5.50 2.7479 2.7457 2.7450
5.75 2.8241 2.8209 2.8220
6.00 2.8920 2.8887 2.8918
6.25 2.9514 2.9486 2.9531
6.50 3.0022 3.0005 3.0049
6.75 3.0447 3.0443 3.0463
7.00 3.0796 3.0804 3.0762
7.25 3.1074 3.1092
7.50 3.1291 3.1314
7.75 3.1453 3.1476
Ionic Strength Dependence 275
Table III. Comparison of the Activity Coefficients of CaCl2(aq) from Table 11 of Rard
and Clegg with Those Calculated Using Eq. (16) and Ion-Interaction Parameters in
Table I
cients from these data are essentially the same and need not be discussed in
detail here. Only the values given large weights in these reviews were included
in the present calculations. A few more recent measurements of Kuschel and
Seidel(30) were included.
For the very dilute range, the values from the freezing point measure-
ments of Gibbard and Gossmann(31) can be converted using thermal data to
yield accurate osmotic coefficients at 25°C; these were included up to 1.02
mol-kg -1 . In addition, in the extremely dilute range below 0.01 mol-kg-1,
activity coefficient values from diffusion data of Harned(32) were included.
For an accurate representation over the full range to 5.9 mol-kg -1 , an
equation extending to the fourth virial with six terms was required. Since
there is no ambiguity of interpretation of the fourth virial as an MgCl3 term,
the equation is fully satisfactory for use in mixed electrolytes. Nevertheless,
it is convenient to have a simpler equation if it is valid over a substantial
range. Thus, we also present an equation with only B(0), B(1), C(0), C(1) terms
that is valid to 4.5 mol-kg -1 . This is essentially the same pattern as the
Archer equations(12,13) for NaBr(aq) and NaCl(aq), the only difference being
in the simpler ionic strength dependence of the C(1) term in the present
equation. The parameters for these equations are in Table IV.
It is interesting to compare the equation first presented in 1973(2) and evaluate
its validity for further use. As expected, the MgCl2 np = 3 equation deviates in
the range above 2.0 mol-kg-1 where the Rard-Miller measurements(21) differ
significantly from the earlier values. Below 2.0 mol-kg -1 , the agreement shown
in Fig. 2 is acceptable for many purposes. However, the new np = 4 equation
278 Pitzer, Wang, Rard, and Clegg
Table IV. Parameter Values for MgCl2(aq) in Eqs. (14-17) with 3, 4, and 6 Ion-
Interaction Parameters
Parameter Units np = 6a np = 4a np = 3b
B(0) kg-mol -1 0.563113 0.512290 0.35235
B(1) kg-mol -1 1 .95628 1 .97367 1.6815
kg-mol -1
(2)
B 0 0 0
ABI kg1/2-mol-1/2 3.0 2.8 2.0
C(0) kg 2 -mol -2 -0.0279311 -0.00676278 0.0025959
C(1) kg 2 -mol -2 0.0468272 -0.0665035
C(2) kg 2 -mol -2 -0.129805
ACI kg-mol -1 0.15 0.23
AC2 kg-mol -1 0.25
D(0) kg3-mol-3 0.00131506
mmax mol-kg -1 5.9 4.5 2
of 0.0028 0.0030
Number of weighted points 176 150
a
This study.
b
Pitzer and Mayorga (Ref. 2).
c
Standard deviation of the fit.
Fig. 2. Differences of the experimental osmotic coefficients (symbols) of MgCI2(aq) from the
model equation, Eq. (14), and parameters given in Table IV, fit to 5.9 mol-kg -1 (np = 6): +
Ref. 21; • Ref. 27; A Ref. 29; X Ref. 28; * Ref. 22; * Ref. 24; - Ref. 25; - Ref. 26; a Ref.
30; 0 Ref. 31; o Ref. 23. Differences are also shown for the model fit to 4.5 mol-kg-1 (np
= 4, dotted curve) and the model by Pitzer and Mayorga (Ref. 2; np = 3, dashed curve) from
those fit to 5.9 mol-kg-1 (np = 6).
Ionic Strength Dependence 279
is no more complex than the Archer equations(12,13) that are coming into general
use; hence, the new np = 4 equation should be used for new research.
Tables V and VI give osmotic and activity coefficient values for these
equations; Fig. 2 shows comparisons with the array of measurements for the
appropriate range of molality.
6. DISCUSSION
If there is significant ion association at very low concentration, it is
unambiguous and can be evaluated from thermodynamic, conductance, or
structural measurements. Aqueous sulfuric acid is an example. In addition,
for cases where there is partial association at low concentration with redissoci-
ation at higher molality, as in MgSO4(aq) and other 2-2 electrolytes, the
method with B(2) = -K/2 is satisfactory(5,33) At higher concentrations, the
population of positive and negative ion near neighbors is substantial on a
Debye-Huckel basis and is subject to large variation for differences in short-
range forces. Thus, any assignment to ion association is ambiguous and
alternate descriptions may be both simpler and preferable for various reasons.
The method described in this paper is shown to give an accurate descrip-
tion of the CaCl2(aq) system at 25°C where any significant ion association
occurs only at high molality. This system has important advantages over a
representation that includes an assumed ion-association constant. One is
convenience of calculation, even for the pure electrolyte. The calculation of
the fractional association at each point and the corresponding ionic strength is
avoided. For mixed electrolytes, it is more important to avoid the unnecessary
complexity introduced by the assumption of additional species of ambiguous
molality. Thus, we strongly recommend the present system when higher-
order terms are needed, in addition to the long-established formulation (1,2,4,5)
with B(0), B(1), B(2), and C terms.
ACKNOWLEDGMENTS
The work at Berkeley was supported by the Director, Office of Energy
Research, Office of Basic Energy Sciences, Division of Chemical Sciences
of the U. S. Department of Energy under Contract No. DE-AC03-76SF00098.
The contribution of J. A. R. was performed under the auspices of the Office
of Basic Energy Sciences (Geosciences) of the U. S. Department of Energy
by the Lawrence Livermore National Laboratory under contract No. W-7405-
ENG-48. That of S.L.C. was supported by an Advanced Research Fellowship
(GT5/93/AAPS/2) from the Natural Environment Research Council of the
U.K.
280 Pitzer, Wang, Rard, and Clegg
Table V. Calculated Osmotic Coefficients of MgCl2(aq) Using Eq. (14) and Ion-
Interaction Parameters in Table IV
ma np = 6b np = 4b np = 3c
Table VI. Calculated Activity Coefficients of MgCl2 (aq) Using Eq. (16) and Ion-
Interaction Parameters in Table IV
ma np = 6b np = 4b np = 3c
REFERENCES
1. K. S. Pitzer, J. Phys. Chem. 77, 268 (1973).
2. K. S. Pitzer and G. Mayorga, 7. Phys. Chem. 77, 2300 (1973).
3. H. L. Friedman, Ionic Solution Theory, Chap. 3 (Wiley-Interscience, New York, 1962).
4. K. S. Pitzer and J. J. Kim, J. Am. Chem. Soc. 96, 5701 (1974).
5. K. S. Pitzer, in Activity Coefficients in Electrolyte Solutions, 2nd edn., Chap. 3, K. S.
Pitzer, ed. (CRC Press, Boca Raton, FL, 1991).
6. C. E. Harvie and J. H. Weare, Geochim. Cosmochim. Acta 44, 981 (1980).
7. V. K. Filippov and V. I. Nokhrin, Russ. J. Inorg. Chem. 30, 282, 289 (1985).
8. V. K. Filippov, N. A. Charykov, and V. Rumyantsev, Dokl. Akad. Nauk, SSSR, Fiz. Khim.
296, 665 (1987) (Engl. Transl. p. 936).
9. R. G. Anstiss and K. S. Pitzer, J. Solution Chem. 20, 849 (1991).
10. J. A. Rard and S. L. Clegg, J. Chem. Eng. Data 42, 819 (1997).
11. J. Ananthaswamy and G. Atkinson, 7. Chem. Eng. Data 30, 120 (1985).
12. D. G. Archer, J. Phys. Chem. Ref. Data 20, 509 (1991).
13. D. G. Archer, 7. Phys. Chem. Ref. Data 21, 793 (1992).
14. S. L. Clegg, J. A. Rard, and K. S. Pitzer, J. Chem. Soc. Faraday Trans. 90, 1875 (1994).
15. R. C. Phutela and K. S. Pitzer, J. Solution Chem. 12, 201 (1983).
16. N. A. Hewish, G. W. Neilson, and J. E. Enderby, Nature 297, 138 (1982).
17. J. A. Rard, A. Habenschuss, and F. H. Spedding, J. Chem. Eng. Data 22, 180 (1977).
18. R. A. Robinson and R. H. Stokes, Electrolyte Solutions, 2nd edn. revised, (Butterworths,
London, 1965).
19. D. G. Archer and P. Wang, J. Phys. Chem. Ref. Data 19, 371 (1990).
20. R. N. Goldberg and R. L. Nuttall, J. Phys. Chem. Ref. Data 7, 263 (1978).
21. J. A. Rard and D. G. Miller, J. Chem. Eng. Data 26, 33, 38 (1981).
22. R. H. Stokes, Trans. Faraday Soc. 41, 642 (1945).
23. R. A. Robinson and R. H. Stokes, Trans. Faraday Soc. 36, 733 (1940).
24. R. A. Robinson and V. E. Bower, 7. Res. Natl. Bur. Stand. U.S. A70, 305 (1966).
25. R. F. Platford, 7. Phys. Chem. 72, 4053 (1968).
26. Y. C. Wu, R. M. Rush, and G. Scatchard, J. Phys. Chem. 72, 4048 (1968); 73, 2047 (1969).
27. Yu. G. Frolov, V. P. Nikolaev, M. Kh. Karapet'yants, and K. K. Vlasenko, Russ. J. Phys.
Chem. 45, 1054 (1971).
28. J. Padova and D. Saad, 7. Solution Chem. 6, 57 (1977).
29. D. Saad, J. Padova, and Y. Marcus, 7. Solution Chem. 4, 983 (1975).
30. F. Kuschel and J. Seidel, 7. Chem. Eng. Data 30, 440 (1985).
31. H. F. Gibbard, Jr., and A. F. Gossmann, 7. Solution Chem. 3, 385 (1974).
32. H. S. Harned, in The Structure of Electrolytic Solutions, Chap. 10, W. J. Hamer, ed., (Wiley,
New York, 1959).
33. K. S. Pitzer and G. Mayorga, 7. Solution Chem. 3, 539 (1974).