This keyword data block is used to define the amounts of an assemblage of pure phases that can react reversibly with the aqueous phase. When the phases included in this keyword data block are brought into contact with an aqueous solution, each phase will dissolve or precipitate to achieve equilibrium or will dissolve completely. Pure phases include minerals with fixed composition and gases with fixed partial pressures. Two types of input are available: in one type, the phase itself reacts to equilibrium (or a specified saturation index or gas partial pressure); in the other type, an alternative reaction occurs to the extent necessary to reach equilibrium (or a specified saturation index or gas partial pressure) with the specified pure phase.
Line 0: EQUILIBRIUM_PHASES 1 Define amounts of phases in assemblage. Line 1a: Chalcedony 0.0 0.0 Line 1b: CO2(g) -3.5 1.0 Line 1c: Gibbsite(c) 0.0 KAlSi3O8 1.0 Line 1d: Calcite 1.0 Gypsum 1.0 Line 1e: pH_Fix -5.0 HCl 10.0
Line 0: EQUILIBRIUM_PHASES [ number ] [ description ]
EQUILIBRIUM_PHASES is the keyword for the data block. Optionally, EQUILIBRIUM, EQUILIBRIA, PURE_PHASES, PURE.
number --Positive number to designate the following phase assemblage and its composition. A range of numbers may also be given in the form m-n , where m and n are positive integers, m is less than n , and the two numbers are separated by a hyphen without intervening spaces. Default is 1.
description --Optional comment that describes the phase assemblage.
Line 1: phase name [ saturation index [( alternative formula or alternative phase ) [ amount ]]]
phase name --Name of a phase. The phase must be defined with PHASES input, either in the database file or in the current or previous simulations of the run. The name must be spelled identically to the name used in PHASES input (except for case).
saturation index --Target saturation index for the pure phase in the aqueous phase (line 1a); for gases, this number is the log of the partial pressure (line 1b). The target saturation index (partial pressure) may not be attained if the amount of the phase in the assemblage is insufficient. Default is 0.0.
alternative formula --Chemical formula that is added (or removed) to attain the target saturation index (or log partial pressure). By default, the mineral defined by phase name dissolves or precipitates to attain the target saturation index. If alternative formula is entered, phase name does not react; the stoichiometry of alternative formula is added or removed from the aqueous phase to attain the target saturation index. Alternative formula must be a legitimate chemical formula composed of elements defined to the program. Line 1c indicates that the stoichiometry given by alternative formula , KAlSi 3 O 8 (potassium feldspar), will be added or removed from the aqueous phase until gibbsite equilibrium is attained. The alternative formula and alternative phase are mutually exclusive fields.
alternative phase --The chemical formula defined for alternative phase is added (or removed) to attain the target saturation index (or log partial pressure). By default, the mineral defined by phase name dissolves or precipitates to attain the target saturation index. If alternative phase is entered, phase name does not react; the stoichiometry of the alternative phase is added or removed from the aqueous phase to attain the target saturation index. Alternative phase must be defined through PHASES input (either in the database file or in the present or previous simulations). Line 1d indicates that the phase gypsum will be added to or removed from the aqueous phase until calcite equilibrium is attained. The alternative formula and alternative phase are mutually exclusive fields.
amount --Moles of the phase in the phase assemblage or moles of the alternative reaction. This number of moles defines the maximum amount of the mineral or gas that can dissolve. It may be possible to dissolve the entire amount without reaching the target saturation index, in which case the solution will have a smaller saturation index for this phase than the target saturation index. If amount is equal to zero, then the phase can not dissolve, but will precipitate if the solution becomes supersaturated with the phase. Default is 10.0 moles.
If just one number is included on line 1, it is assumed to be the target saturation index (or log partial pressure) and the amount of the phase defaults to 10.0 mol. If two numbers are included on the line, the first is the target saturation index and the second is the amount of the phase present. Line 1 may be repeated to define all pure phases that are assumed to react reversibly. It is possible to include a pure phase that has an amount of zero (line 1a). In this case, chalcedony can not dissolve, but can precipitate if the solution is supersaturated with chalcedony, either by initial conditions, or through dissolution of pure phases or other specified reactions (mixing, stoichiometric or kinetic reactions). It is possible to maintain constant pH conditions by specification of an alternative formula and a special phase ( PHASES input). Line 1e would maintain a pH of 5.0 by adding HCl, provided a phase named "pH_Fix" were defined with reaction H + = H + and log K = 0.0 (see example 8, in "Examples"). (Note: If the acid, HCl, is specified and, in fact, a base is needed to attain pH 5.0, it is possible the program will fail to find a solution to the algebraic equations.)
The number of exchange sites can be related to the moles of a phase that are present in an EQUILIBRIUM_PHASES phase assemblage (see EXCHANGE). As the moles of the phase increase or decrease, the number of exchange sites will increase or decrease. Likewise, the number of surface sites can be related to the moles of a phase that are present in an EQUILIBRIUM_PHASES phase assemblage (see SURFACE).
For batch reactions, after a pure-phase assemblage has reacted with the solution, it is possible to save the resulting assemblage composition (that is, the identity, target saturation index, and moles of each phase) with the SAVE keyword. If the new composition is not saved, the assemblage composition will remain the same as it was before the batch reaction. After it has been defined or saved, the assemblage may be used in subsequent simulations by the USE keyword. TRANSPORT and ADVECTION calculations automatically update the pure-phase assemblage and SAVE has no effect during these calculations.
The keyword EQUILIBRIUM_PHASES is used in example problems 2, 3, 5, 6, 7, 8, 9, 10, and 14.
ADVECTION, EXCHANGE, PHASES, SAVE equilibrium_phases, SURFACE, TRANSPORT, and USE equilibrium_phases.