Nomenclature

Terminology

Species

Species are real things that exist in solution, for instance HCO3_{3}^{-}, CO2_{2}(aq), CaHCO3+_{3}^{+}. They never have negative mass.

Phase

A phase is a region of space that is physically distinct, mechanically separable, and homogeneous in its composition and properties. In geochemical systems, the phases are:

  • a fluid phase composed of water and its dissolved constituents;

  • potentially, a separate phase corresponding to each of the minerals in contact with the fluid phase;

  • potentially, separate gas phases. If different gases (CO2_{2}, O2_{2}, etc) are buffered separately with known fugacity, as is assumed in geochemistry, each gas counts as a distinct phase.

Components

Components are useful abstractions. They might or might not be real things, so can even have negative mass. For instance, consider an alkaline solution (more OH^{-} than H+^{+}). Choose H2_{2}O and H+^{+} to be the component basis set. Each hydroxyl ion, H2OH+OH\mathrm{H}_{2}\mathrm{O} - \mathrm{H}^{+} \rightarrow \mathrm{OH}^{-}, is made up of a water molecule less a hydrogen ion. Since this is an alkaline solution, the overall composition has a positive amount of water and a negative amount of H+^{+}. The molality of the H+^{+} species is positive, however.

A complete set of components forms a basis for the reaction system.

Basis

The set of components used in a geochemical model is the calculation's basis. There is no unique basis that describes a system, but the basis must satisfy these rules:

  1. each species and phase considered in the model must be a combination of the components in the basis;

  2. the number of components is minimal to satisfy (1)

  3. the components must be linearly independent of one another, that is, no reactions exist to form one component from the others.

That is, the basis is a linearly-independent complete set of components: every other chemical in the system is composed of these components, and none of these components can be composed from the others.

The number of components in the basis is 1+Ni+Nk+Nm+Np1 + N_{i} + N_{k} + N_{m} + N_{p} (the "1" comes from the water component), while the number of phases is 1+Nk+Nm+Np1 + N_{k} + N_{m} + N_{p}. The Gibbs phase rule, with fixed temperature and pressure, states that the number of degrees of freedom is the difference of these, or just NiN_{i}. Hence, given a constraint on the concentration or activity of each basis species, the system's equilibrium state can be determined. However, calculating the extent of the system — the amounts of fluid, minerals, gases and sorption sites — that are present, requires the extra 1+Nk+Nm+Np1 + N_{k} + N_{m} + N_{p} pieces of information.

Moles

A measure of the amount of a substance.

number of moles=number of particles (molecules, atoms, ions, etc)NA ,\mathrm{number\ of\ moles} = \frac{\mathrm{number\ of\ particles\ (molecules,\ atoms,\ ions,\ etc)}}{\mathcal{N}_{A}} \ , where NA\mathcal{N}_{A} is Avogadro's number.

Mole fraction

A measure of the amount of a constituent relative to the mixture.

mole fraction=number of moles of constituentnumber of moles of all constituents in the mixture\mathrm{mole\ fraction} = \frac{\mathrm{number\ of\ moles\ of\ constituent}}{\mathrm{number\ of\ moles\ of\ all\ constituents\ in\ the\ mixture}}

Molality

A measure of concentration of a constituent.

molality=number of moles of constituentmass (kg) of solvent (water)\mathrm{molality} = \frac{\mathrm{number\ of\ moles\ of\ constituent}}{\mathrm{mass\ (kg)\ of\ solvent\ (water)}}

Chemical potential

The partial derivative of the free energy, GG, with respect to the number of moles of that species: μ=G/n\mu = \partial G/\partial n, where temperature, pressure and the number of moles of other species are held fixed.

Activity

Defined in terms of the chemical potential, μ=μ0+RTloga ,\mu = \mu^{0} + RT\log a \ , where

  • μ\mu [J.mol1^{-1}] is the chemical potential

  • μ0\mu^{0} [J.mol1^{-1}] is the chemical potential in a standard state (it is independent of temperature, pressure, other constituents, etc)

  • R=8.314R = 8.314\ldots\,J.K1^{-1}.mol1^{-1} is the gas constant

  • TT (K) is temperature

  • log\log is the natural logarithm

  • aa [dimensionless] is the activity

For species equilibrium reactions, the activity is a dimensionless measure of effective concentration of a constituent in a solution a=γm/mref ,a = \gamma m/m_{\mathrm{ref}} \ , where γ\gamma is the activity coefficient, mm is the molality, and mref=1m_{\mathrm{ref}}=1\,mol.kg1^{-1}.

For minerals in AkA_{k}, a=1a=1.

For gases, the activity is proportional to the fugacity.

Mass action

This is an equation for the equilibrium constant of a reaction in terms of the activities of the reaction's species.

Saturation index

SI=log10(Q/K)\mathrm{SI} = \log_{10}(Q/K), where QQ is the activity product and KK is the equilibrium constant of the reaction. If SI>0\mathrm{SI}>0, the mineral is supersaturated, which means the system will try to precipitate the mineral.

Redox reaction

This involves the transfer of electrons. For instance Fe2++0.25O2+H+Fe3++0.5H2O\mathrm{Fe}^{2+} + 0.25\mathrm{O}_{2} + \mathrm{H}^{+} \rightleftharpoons \mathrm{Fe}^{3+} + 0.5\mathrm{H}_{2}\mathrm{O}, where charge is transferred between Fe2+\mathrm{Fe}^{2+} and O2\mathrm{O}_{2}. Fe2+\mathrm{Fe}^{2+} has lost an electron: it is oxidised. Often these reactions are split into half reactions that show the release/gain of an electron, e.g. Fe2+Fe3++e\mathrm{Fe}^{2+} \rightarrow \mathrm{Fe}^{3+} + \mathrm{e}^{-}. A redox couple is a reduced species and its corresponding oxidised form, e.g. Fe2+\mathrm{Fe}^{2+} / Fe3+\mathrm{Fe}^{3+}.

Nernst potential

Quantifies the redox potential for a redox half-reaction with activity product QQ and equilibrium constant KK: Eh=RTnFlog(Q/K)=Eh0RTnFlogQ .\mathrm{Eh} = -\frac{RT}{nF}\log(Q/K) = \mathrm{Eh}^{0} - \frac{RT}{nF}\log Q \ . Here RR is the gas constant, TT is the temperature (in Kelvin), FF is the Faraday constant, nn is the number of electrons consumed in the half reaction, Eh0^{0} is the standard potential, and log\log has base e.

Notation

Table 1: List of notation used in documentation

SymbolUnitsPhysical description
\rightleftharpoons-Denotes an equilibrium reaction, for instance B3C+2DB \rightleftharpoons 3C + 2D
AwA_{w}-Label for water.
AiA_{i}-Label for aqueous basis species. This is a label, not a quantity. For example, one of the AiA_{i} might be Na+^{+}. This includes decoupled non-kinetically-controlled alternative oxidation states as well as minerals that have zero precipitate.
AjA_{j}-Label for other aqueous species in equilibrium with the basis, which are called the "secondary" species. Included in this set are all coupled redox species and minerals that are not precipitated, but not species involved in sorption.
AjˉA_{\bar{j}}-Redox species that are governed by kinetics
AkA_{k}-Label for precipitated minerals in the equilibrium system of reactions. These will appear in the basis.
AkˉA_{\bar{k}}-Label for all components that are governed by a kinetic rate. This is a union of AjˉA_{\bar{j}}, AlˉA_{\bar{l}} and AqˉA_{\bar{q}}.
AlA_{l}-Label for all minerals, including those in AjA_{j} and AkA_{k} and even those that are kinetically controlled.
AlˉA_{\bar{l}}-Label for all minerals that are kinetically controlled.
AmA_{m}-Label for gases of known fugacity, so these will appear in the basis.
AnA_{n}-Label for all gases.
AqA_{q}-Label for sorbed species in equilibrium with the aqueous solution. These are combinations of basis species that have adsorbed onto sorbing sites on minerals, making a surface complex. The mass-action equations can be more complicated than for AjA_{j} and AlA_{l}. These do not get transported.
AqˉA_{\bar{q}}-Label for sorbed species that are kinetically controlled.
ApA_{p}-Label for sorbing sites, so these appear in the basis. In the Langmuir approach there is just one of these, and all sorbing species compete to sorb. In the surface-complexation approach there are Np1N_{p}\geq 1 of these.
aa-Activity. For species in equilibrium reactions a=γm/mrefa = \gamma m/m_{\mathrm{ref}}, where γ\gamma is the activity coefficient, mm is the molality, and mref=1m_{\mathrm{ref}}=1\,mol.kg1^{-1}.
CCmol.m3^{-3}Concentration of a component per volume of aqueous solution (the concentration of component per volume of porous skeleton is ϕC\phi C)
CT,C_{T,\ast}mol.m3^{-3}Concentration of mobile component per volume of aqueous solution
cwc_{w}kg.m3^{-3}Mass of solute water per volume of aqueous solution
ckc_{k}mol.m3^{-3}Concentration of mineral component AkA_{k} per volume of aqueous solution
ffPaGas fugacity
ϕ\phi-Porosity of porous skeleton. Its dynamics (if any) are not controlled by the geochemistry module
γ\gamma-Activity coefficient
IImol.kg1^{-1}Ionic strength: I=12imizi2I = \frac{1}{2}\sum_{i}m_{i}z_{i}^{2}, where the sum runs over the free ions in the solution, so neutral complexes are not counted.
IsI_{s}mol.kg1^{-1}Stoichiometric ionic strength: Is=12imizi2I_{s} = \frac{1}{2}\sum_{i}m_{i}z_{i}^{2}. Here, complete dissociation of complexes is assumed, and then the sum runs over all ions in the hypothetical solution.
KK-Equilibrium constant for a reaction
mmmol.kg1^{-1}Molality. The molality of unoccupied sorption sites is mpm_{p}.
MMmolMoles of a component
μ\muJ.mol1^{-1}Chemical potential: μ=G/n\mu = \partial G/\partial n — the derivative of the free energy, GG, with respect to the number of moles of that species, with all other things (temperature, number of moles of other species) held constant
nnkgMass of a component. For instance, nwn_{w} is the mass of solvent water
nkn_{k}molMole number of mineral kk.
NAN_{A}mol1^{-1}Avogadro constant: NA=6.02214076×1023N_{A} = 6.02214076\times 10^{23}\,mol1^{-1}
NA\mathcal{N}_{A}-Avogadro number: NA=6.02214076×1023\mathcal{N}_{A} = 6.02214076\times 10^{23}
NiN_{i}-Number of basis species AiA_{i}
NkN_{k}-Number of basis minerals AkA_{k}
NmN_{m}-Number of basis gases AmA_{m}
NpN_{p}-Number of basis sorbing sites ApA_{p}
ν\nu-Stoichiometric coefficient
QlQ_{l}-Activity product for mineral ll. This is used to construct the saturation index for the mineral
rkˉr_{\bar{k}}mol.s1^{-1}Reaction rate for kinetic species AkˉA_{\bar{k}}
RkˉR_{\bar{k}}mol.s1^{-1}.m3^{-3}Reaction rate per volume of aqueous solution for kinetic species AkˉA_{\bar{k}}
SIl_{l}-Saturation index for mineral ll. The mineral is supersaturated if this is positive, and undersaturated if it is negative
zz-Charge number of a substance