Page 360 - Physical Chemistry
P. 360

lev38627_ch11.qxd  3/14/08  1:10 PM  Page 341





                                                                                                                 341
                      The programs MINTEQA2 (www.epa.gov/ceampubl/mmedia/minteq/index.htm;                   Section 11.7
                  see also www.lwr.kth.se/English/oursoftware/vminteq), PHREEQC (wwwbrr.cr.usgs.  Temperature and Pressure Dependences
                                                                                                    of the Equilibrium Constant
                  gov/projects/GWC_coupled/phreeqc), and MINEQL  (www.mineql.com) are com-
                  putationally similar programs for calculations on natural and laboratory aqueous sys-
                  tems of ionic strength less than 0.5 mol/kg. Ion activity coefficients are estimated
                  from the Debye–Hückel formula (10.57) with a term proportional to I added, using
                                                                              m
                  parameters fitted for each species. If these parameters are not available, the Davies
                  equation is used. Activity coefficients for uncharged solute species are estimated from
                  the formula
                                             log   g   0.11I >m°2
                                                          m
                                                10
                                                   i
                  The experimentally determined activity coefficients of several uncharged solutes in
                  aqueous solution have been found to fit the equation log g   b (I /m°), where b i
                                                                             m
                                                                           i
                                                                   10  i
                  differs for different species. For several species in water at 25°C, b is roughly 0.1.
                                                                             i
                  PHREEQC can also use the Pitzer method for activity coefficients.
                      The program EQS4WIN (www.mathtrek.com) solves multiple-equilibria prob-
                  lems that can involve several phases but no activity or fugacity coefficients are used,
                  the phases being assumed to be ideal. A free demonstration version is available.
                      ChemSage [www.esm-software.com/chemsage/; G. Eriksson and K. Hack, Metal-
                  lurg. Trans. B, 21B, 1013 (1990)] is a general-purpose DOS program that does nonideal
                  chemical equilibrium calculations in multiphase systems at temperatures up to 6000 K
                  and pressures to 1 Mbar. For nonideal gases, fugacity coefficients are found using the
                  virial equation. For aqueous electrolyte solutions, the Pitzer equations are used. For
                  nonelectrolyte solutions, the user can choose from several models for estimating activ-
                  ity coefficients. FactSage (www.factsage.com) is a Windows successor to ChemSage.
                      HSC Chemistry (www.esm-software.com/hsc/ and www.hsc-chemistry.net) is a
                  Windows program for multiphase chemical equilibrium calculations. Gases are as-
                  sumed ideal and the user must specify activity coefficients for species in solution.
                      The programs Equilib-Web (www.crct.polymtl.ca/equiweb.html) and Aqualib-
                  Web (www.crct.polymtl.ca/aquaweb.html) do online calculations for ideal-gas sys-
                  tems or ideally dilute aqueous solutions and are limited to systems with no more than
                  three reactants and five elements.



                    11.7         TEMPERATURE AND PRESSURE DEPENDENCES
                                 OF THE EQUILIBRIUM CONSTANT

                  From  G°   RT ln K° [Eq. (11.4)], we have
                                              ln K°   ¢G°>RT                        (11.31)
                  where  G°     n m° is the standard change in Gibbs energy (all species in their stan-
                                  i
                                i
                                    i
                  dard states). For gases and for pure liquids and solids, we chose a fixed-pressure stan-
                  dard state (P°   1 bar), so here  G° and hence K° are independent of pressure and
                  depend only on T. For liquid and solid solutions, we chose a variable-pressure stan-
                  dard state, with standard-state pressure equal to the actual pressure of the solution, so
                  here  G° and K° are functions of both T and P.
                      Differentiation of (11.31) with respect to T gives

                           0 ln K°    ¢G°    10¢G°>0T 2 P  ¢G°   ¢S°    ¢G°   T ¢S°
                         a       b
                            0T    P   RT  2      RT        RT  2  RT        RT  2
                                               0 ln K°     ¢H°
                                               a      b                            (11.32)*
                                                 0T    P   RT 2
   355   356   357   358   359   360   361   362   363   364   365