[Next] [Previous] [Up] [Top]

EXAMPLES

Example 8.-- Surface Complexation

PHREEQC contains three surface-complexation models: (1) By default, an electrostatic double layer model is used with no explicit calculation of the diffuse-layer composition. (2) Alternatively, an electrostatic double layer model with explicit calculation of the diffuse-layer composition may be used (-diffuse_layer). (3) Finally, a non-electrostatic model may be selected (-no_edl). The electrostatic model is the diffuse double-layer model described in Dzombak and Morel (1990) with the following modifications: (1) surfaces may have more than two types of binding sites, (2) surface precipitation is not included, and (3) optionally, an alternative formulation for the charge-potential relationship, modified from Borkovec and Westall (1983), that explicitly calculates the composition of the diffuse layer can be employed (-diffuse_layer). The non-electrostatic model does not consider the effects of the development of surface charge on the formation of surface complexes, with the result that surface complexes are treated mathematically very much like aqueous complexes without activity coefficient terms.

The following example of the diffuse double-layer model is taken from Dzombak and Morel (1990, chapter 8) with no explicit calculation of the diffuse-layer composition. Zinc sorption on hydrous ferric oxide is simulated assuming two types of sites, weak and strong, are available on the oxide surface. Protons and zinc ions compete for the two types of binding sites, and equilibrium is described by mass-action equations. Activities of the surface species depend on the potential at the surface, which is due to the development of surface charge. The example considers the variation in sorption of zinc on hydrous ferric oxides as a function of pH for low zinc concentration (10-7 m) and high zinc concentration (10-4 m) in 0.1 m sodium nitrate electrolyte.

Surface-complexation reactions derived from the summary of Dzombak and Morel (1990) are contained in the default database files for PHREEQC. However, many of the intrinsic stability constants used in this example differ from the values in the default database files and definitions are thus included in the input file (table 16). Three keyword data blocks are required to define surface-complexation data for a simulation: SURFACE_MASTER_SPECIES, SURFACE_SPECIES, and SURFACE.The SURFACE_MASTER_SPECIES data block in the default database files selects surface species to be the master species for the binding sites of "Hfo" (hydrous ferric oxides). The name of a binding site is composed of a name for the surface, "Hfo" in the default database files, optionally followed by an underscore and a lowercase binding site designation,"Hfo_w" and "Hfo_s" for "weak" and "strong" in the database files. The underscore notation is necessary only if two or more binding sites exist for a single surface. The notation allows a mole-balance equation to be derived for each of the binding sites (Hfo_w and Hfo_s, in this example) and a single charge-potential or charge-balance equation for the surface (Hfo, in this example). Thus, the charge that develops on each binding site will enter into a single charge-potential or charge-balance equation for the surface.

Table 16. Input data set for example 8

TITLE Example 8.--Sorption of zinc on hydrous 
                  iron oxides
SURFACE_SPECIES
        Hfo_sOH  + H+ = Hfo_sOH2+
        log_k  7.18

        Hfo_sOH = Hfo_sO- + H+
        log_k  -8.82

        Hfo_sOH + Zn+2 = Hfo_sOZn+ + H+
        log_k  0.66

        Hfo_wOH  + H+ = Hfo_wOH2+
        log_k  7.18

        Hfo_wOH = Hfo_wO- + H+
        log_k  -8.82

        Hfo_wOH + Zn+2 = Hfo_wOZn+ + H+
        log_k  -2.32
SURFACE 1
        Hfo_sOH        5e-6    600.    0.09
        Hfo_wOH        2e-4    
SOLUTION 1
        -units  mmol/kgw
        pH      8.0
        Zn      0.0001 
        Na      100.    charge 
        N(5)    100.
SOLUTION 2
        -units  mmol/kgw
        pH      8.0
        Zn      0.1 
        Na      100.    charge 
        N(5)    100.
USE surface none
#
# Model definitions
#
PHASES
        Fix_H+
        H+ = H+
        log_k  0.0
END
#
#   Zn = 1e-7
#
SELECTED_OUTPUT
        -file ex8.pun
        -molalities     Zn+2    Hfo_wOZn+      Hfo_sOZn+
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.0   NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.25  NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.5   NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.75  NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.0   NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.25  NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.5   NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.75  NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.0   NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.25  NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.5   NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.75  NaOH    10.0
END
USE solution 1
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -8.0   NaOH    10.0
END
#
#   Zn = 1e-4
#
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.0   NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.25  NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.5   NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -5.75  NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.0   NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.25  NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.5   NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -6.75  NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.0   NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.25  NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.5   NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -7.75  NaOH    10.0
END
USE solution 2
USE surface 1
EQUILIBRIUM_PHASES 1
        Fix_H+   -8.0   NaOH    10.0
END
The chemical reactions and thermodynamic constants for all surface species, including the surface master species, are defined with the SURFACE_SPECIES data block. The mass-action equations taken from Dzombak and Morel (1990, p. 259) are given in the input data set (table 16), under keyword SURFACE_SPECIES. Note the activity coefficient or potential term is not included as part of the mass-action expression; the potential term is added internally by the program.

The composition and other characteristics of an assemblage of surfaces is defined with the SURFACE data block. The composition of multiple surfaces, each with multiple binding sites, may be defined within the data block. For each surface, the number of moles of sites, the initial composition of the surface, and the surface area must be defined. Although the composition of the surfaces may change due to reactions, the number of surfaces, moles of binding sites, and surface areas remain fixed until the end of the run or until the entire assemblage is redefined. In this example, one surface (Hfo) with two binding sites (Hfo_w and Hfo_s) is defined. The number of moles of strong binding sites, Hfo_s, is 5x10-6 sites and the number of moles of weak binding sites, Hfo_w, is 2x10-4. Initially, all surface sites are in the uncharged, protonated form. The surface area for the entire surface, Hfo, must be defined with two numbers, the area per mass of surface material (600 m2/g, in this example) and the total mass of surface material (0.09 g, in this example). The use of these two numbers to define surface area is traditional, but only the product of these numbers is used in the model to obtain the surface area; the individual numbers are not used separately. Surface area may be entered with the data for any of the binding sites for a surface; in this example, the surface area is entered with Hfo_s.

To complete the definition of the initial conditions for the simulations, two sodium nitrate solutions are defined with differing concentrations of zinc (SOLUTION 1 and 2 data blocks). A pseudo-phase, "Fix_H+" is defined with the PHASES data block. This phase is used in each of the reaction simulations to adjust pH to fixed values. Finally, the line "USE surface none" eliminates an implicitly defined reaction calculation for the first simulation. By default, if a SOLUTION and SURFACE data block are defined in a simulation, then the first solution defined in the simulation (SOLUTION 1 in this example) and the first surface defined in the simulation are put together (possibly with other assemblages and a gas phase) and allowed to equilibrate. The USE keyword with "surface none" removes the surface from any reaction calculated for the simulation, with the effect that no reaction calculation is performed because nothing is defined with which the solution may react. (The same logic applies to the EXCHANGE, GAS_PHASE, EQUILIBRIUM_PHASES, REACTION, REACTION_TEMPERATURE keywords that are defined within the input for a simulation. A reaction step is implicitly defined whenever a solution or mixture is defined in the simulation and any one of these keyword data blocks also is defined in the same simulation.).

The remaining simulations in the input data set equilibrate the surface assemblage with either solution 1 or solution 2 for pH values that range from 5 to 8. Each of the simulations uses the phase "Fix_H+" in an EQUILIBRIUM_PHASES keyword data block with varying saturation indices to adjust pH. The reaction NaOH is added or removed from each solution to produce a specified saturation index which, by the definition of the reaction for "Fix_H+" is numerically equal to the log of the hydrogen activity, or negative pH. Note that, although it is possible to attain the desired pH in all of these simulations, a pH that is sufficiently low will cause the program to fail because a very low pH can not be reached even by removing all of the sodium in solution.

The results of the simulation are plotted on figure 4 and are consistent with the results shown in Dzombak and Morel (1990, figure 8.9). Zinc is more strongly sorbed at high pH values than at low pH values. In addition, at low concentrations of zinc, the strong binding sites outcompete the weak binding sites for zinc over the entire pH range, and at high pH most of the zinc resides at the strong binding sites. At larger zinc concentrations, the strong binding sites predominate only at low pH. Because all the strong binding sites become filled at higher pH, most of the zinc resides at the more numerous weak binding sites at high pH and large zinc concentrations.

Table 16. Input data set for example 8

User's Guide to PHREEQC - 07 MAY 96
[Next] [Previous] [Up] [Top]

Generated with CERN WebMaker