Solubility equilibrium

Solubility equilibrium is a type of dynamic equilibrium that exists when a chemical compound in the solid state is in chemical equilibrium with a solution of that compound. The solid may dissolve unchanged, with dissociation, or with chemical reaction with another constituent of the solution, such as acid or alkali. Each solubility equilibrium is characterized by a temperature-dependent solubility product which functions like an equilibrium constant. Solubility equilibria are important in pharmaceutical, environmental and many other scenarios.

Definitions

A solubility equilibrium exists when a chemical compound in the solid state is in chemical equilibrium with a solution containing the compound. This type of equilibrium is an example of dynamic equilibrium in that some individual molecules migrate between the solid and solution phases such that the rates of dissolution and precipitation are equal to one another. When equilibrium is established and the solid has not all dissolved, the solution is said to be saturated. The concentration of the solute in a saturated solution is known as the solubility. Units of solubility may be molar (mol dm⁻³) or expressed as mass per unit volume, such as μg mL⁻¹. Solubility is temperature dependent. A solution containing a higher concentration of solute than the solubility is said to be supersaturated. A supersaturated solution may be induced to come to equilibrium by the addition of a "seed" which may be a tiny crystal of the solute, or a tiny solid particle, which initiates precipitation.^{[ citation needed ]}

There are three main types of solubility equilibria.

1. Simple dissolution.
2. Dissolution with dissociation reaction. This is characteristic of salts. The equilibrium constant is known in this case as a solubility product.
3. Dissolution with ionization reaction. This is characteristic of the dissolution of weak acids or weak bases in aqueous media of varying pH.

In each case an equilibrium constant can be specified as a quotient of activities. This equilibrium constant is dimensionless as activity is a dimensionless quantity. However, use of activities is very inconvenient, so the equilibrium constant is usually divided by the quotient of activity coefficients, to become a quotient of concentrations. See Equilibrium chemistry#Equilibrium constant for details. Moreover, the activity of a solid is, by definition, equal to 1 so it is omitted from the defining expression.

For a chemical equilibrium

$${\displaystyle \mathrm {A} _{p}\mathrm {B} _{q}\leftrightharpoons p\mathrm {A} +q\mathrm {B} }$$

the solubility product, K_sp for the compound A_pB_q is defined as follows

$${\displaystyle K_{\mathrm {sp} }=[\mathrm {A} ]^{p}[\mathrm {B} ]^{q}}$$

where [A] and [B] are the concentrations of A and B in a saturated solution. A solubility product has a similar functionality to an equilibrium constant though formally K_sp has the dimension of (concentration)^p+q.

Effects of conditions

Temperature effect

Solubility is sensitive to changes in temperature. For example, sugar is more soluble in hot water than cool water. It occurs because solubility products, like other types of equilibrium constants, are functions of temperature. In accordance with Le Chatelier's Principle, when the dissolution process is endothermic (heat is absorbed), solubility increases with rising temperature. This effect is the basis for the process of recrystallization, which can be used to purify a chemical compound. When dissolution is exothermic (heat is released) solubility decreases with rising temperature. ^[1] Sodium sulfate shows increasing solubility with temperature below about 32.4 °C, but a decreasing solubility at higher temperature. ^[2] This is because the solid phase is the decahydrate (Na
_²SO
_⁴·10H
_²O) below the transition temperature, but a different hydrate above that temperature.^{[ citation needed ]}

The dependence on temperature of solubility for an ideal solution (achieved for low solubility substances) is given by the following expression containing the enthalpy of melting, Δ_mH, and the mole fraction ${\displaystyle x_{i}}$ of the solute at saturation:

$${\displaystyle \left({\frac {\partial \ln x_{i}}{\partial T}}\right)_{P}={\frac {{\bar {H}}_{i,\mathrm {aq} }-H_{i,\mathrm {cr} }}{RT^{2}}}}$$

where ${\displaystyle {\bar {H}}_{i,\mathrm {aq} }}$ is the partial molar enthalpy of the solute at infinite dilution and ${\displaystyle H_{i,\mathrm {cr} }}$ the enthalpy per mole of the pure crystal. ^[3]

This differential expression for a non-electrolyte can be integrated on a temperature interval to give: ^[4]

$${\displaystyle \ln x_{i}={\frac {\Delta _{m}H_{i}}{R}}\left({\frac {1}{T_{f}}}-{\frac {1}{T}}\right)}$$

For nonideal solutions activity of the solute at saturation appears instead of mole fraction solubility in the derivative with respect to temperature:

$${\displaystyle \left({\frac {\partial \ln a_{i}}{\partial T}}\right)_{P}={\frac {H_{i,\mathrm {aq} }-H_{i,\mathrm {cr} }}{RT^{2}}}}$$

Common-ion effect

The common-ion effect is the effect of decreased solubility of one salt when another salt that has an ion in common with it is also present. For example, the solubility of silver chloride, AgCl, is lowered when sodium chloride, a source of the common ion chloride, is added to a suspension of AgCl in water. ^[5]

$${\displaystyle \mathrm {AgCl(s)\leftrightharpoons Ag^{+}(aq)+Cl^{-}(aq)} }$$

The solubility, S, in the absence of a common ion can be calculated as follows. The concentrations [Ag⁺] and [Cl⁻] are equal because one mole of AgCl would dissociate into one mole of Ag⁺ and one mole of Cl⁻. Let the concentration of [Ag⁺(aq)] be denoted by x. Then

$${\displaystyle K_{\mathrm {sp} }=\mathrm {[Ag^{+}][Cl^{-}]} =x^{2}}$$

$${\displaystyle {\text{Solubility}}=\mathrm {[Ag^{+}]=[Cl^{-}]} =x={\sqrt {K_{\mathrm {sp} }}}}$$

K_sp for AgCl is equal to 1.77×10⁻¹⁰ mol² dm⁻⁶ at 25 °C, so the solubility is 1.33×10⁻⁵ mol dm⁻³.

Now suppose that sodium chloride is also present, at a concentration of 0.01 mol dm⁻³ = 0.01 M. The solubility, ignoring any possible effect of the sodium ions, is now calculated by

$${\displaystyle K_{\mathrm {sp} }=\mathrm {[Ag^{+}][Cl^{-}]} =x(0.01\,{\text{M}}+x)}$$

This is a quadratic equation in x, which is also equal to the solubility.

$${\displaystyle x^{2}+0.01\,{\text{M}}\,x-K_{sp}=0}$$

In the case of silver chloride, x² is very much smaller than 0.01 M x, so the first term can be ignored. Therefore

$${\displaystyle {\text{Solubility}}=\mathrm {[Ag^{+}]} =x={\frac {K_{\mathrm {sp} }}{0.01\,{\text{M}}}}=\mathrm {1.77\times 10^{-8}\,mol\,dm^{-3}} }$$

a considerable reduction from 1.33×10⁻⁵ mol dm⁻³. In gravimetric analysis for silver, the reduction in solubility due to the common ion effect is used to ensure "complete" precipitation of AgCl.

Particle size effect

The thermodynamic solubility constant is defined for large monocrystals. Solubility will increase with decreasing size of solute particle (or droplet) because of the additional surface energy. This effect is generally small unless particles become very small, typically smaller than 1 μm. The effect of the particle size on solubility constant can be quantified as follows:

$${\displaystyle \log(^{*}K_{A})=\log(^{*}K_{A\to 0})+{\frac {\gamma A_{\mathrm {m} }}{3.454RT}}}$$

where *K_A is the solubility constant for the solute particles with the molar surface area A, *K_A→0 is the solubility constant for substance with molar surface area tending to zero (i.e., when the particles are large), γ is the surface tension of the solute particle in the solvent, A_m is the molar surface area of the solute (in m²/mol), R is the universal gas constant, and T is the absolute temperature. ^[6]

Salt effects

The salt effects ^[7] (salting in and salting-out) refers to the fact that the presence of a salt which has no ion in common with the solute, has an effect on the ionic strength of the solution and hence on activity coefficients, so that the equilibrium constant, expressed as a concentration quotient, changes.

Phase effect

Equilibria are defined for specific crystal phases. Therefore, the solubility product is expected to be different depending on the phase of the solid. For example, aragonite and calcite will have different solubility products even though they have both the same chemical identity (calcium carbonate). Under any given conditions one phase will be thermodynamically more stable than the other; therefore, this phase will form when thermodynamic equilibrium is established. However, kinetic factors may favor the formation the unfavorable precipitate (e.g. aragonite), which is then said to be in a metastable state.^{[ citation needed ]}

In pharmacology, the metastable state is sometimes referred to as amorphous state. Amorphous drugs have higher solubility than their crystalline counterparts due to the absence of long-distance interactions inherent in crystal lattice. Thus, it takes less energy to solvate the molecules in amorphous phase. The effect of amorphous phase on solubility is widely used to make drugs more soluble. ^{[8] [9]}

Pressure effect

For condensed phases (solids and liquids), the pressure dependence of solubility is typically weak and usually neglected in practice. Assuming an ideal solution, the dependence can be quantified as:

$${\displaystyle \left({\frac {\partial \ln x_{i}}{\partial P}}\right)_{T}=-{\frac {{\bar {V}}_{i,\mathrm {aq} }-V_{i,\mathrm {cr} }}{RT}}}$$

where ${\displaystyle x_{i}}$ is the mole fraction of the ${\displaystyle i}$-th component in the solution, ${\displaystyle P}$ is the pressure, ${\displaystyle T}$ is the absolute temperature, ${\displaystyle {\bar {V}}_{i,{\text{aq}}}}$ is the partial molar volume of the ${\displaystyle i}$th component in the solution, ${\displaystyle V_{i,{\text{cr}}}}$ is the partial molar volume of the ${\displaystyle i}$th component in the dissolving solid, and ${\displaystyle R}$ is the universal gas constant. ^[10]

The pressure dependence of solubility does occasionally have practical significance. For example, precipitation fouling of oil fields and wells by calcium sulfate (which decreases its solubility with decreasing pressure) can result in decreased productivity with time.

Quantitative aspects

Simple dissolution

Dissolution of an organic solid can be described as an equilibrium between the substance in its solid and dissolved forms. For example, when sucrose (table sugar) forms a saturated solution

$${\displaystyle \mathrm {C_{12}H_{22}O_{11}(s)\leftrightharpoons C_{12}H_{22}O_{11}(aq)} }$$

An equilibrium expression for this reaction can be written, as for any chemical reaction (products over reactants):

$${\displaystyle K^{\ominus }={\frac {\left\{\mathrm {{C}_{12}{H}_{22}{O}_{11}(aq)} \right\}}{\left\{\mathrm {{C}_{12}{H}_{22}{O}_{11}(s)} \right\}}}}$$

where K^o is called the thermodynamic solubility constant. The braces indicate activity. The activity of a pure solid is, by definition, unity. Therefore

$${\displaystyle K^{\ominus }=\left\{\mathrm {{C}_{12}{H}_{22}{O}_{11}(aq)} \right\}}$$

The activity of a substance, A, in solution can be expressed as the product of the concentration, [A], and an activity coefficient, γ. When K^o is divided by γ, the solubility constant, K_s,

$${\displaystyle K_{\mathrm {s} }=\left[\mathrm {{C}_{12}{H}_{22}{O}_{11}(aq)} \right]}$$

is obtained. This is equivalent to defining the standard state as the saturated solution so that the activity coefficient is equal to one. The solubility constant is a true constant only if the activity coefficient is not affected by the presence of any other solutes that may be present. The unit of the solubility constant is the same as the unit of the concentration of the solute. For sucrose K_s = 1.971 mol dm⁻³ at 25 °C. This shows that the solubility of sucrose at 25 °C is nearly 2 mol dm⁻³ (540 g/L). Sucrose is unusual in that it does not easily form a supersaturated solution at higher concentrations, as do most other carbohydrates.

Dissolution with dissociation

Ionic compounds normally dissociate into their constituent ions when they dissolve in water. For example, for silver chloride:

$${\displaystyle {\ce {AgCl_{(s)}<=> Ag^+_{(aq)}{}+ Cl^-_{(aq)}}}}$$

The expression for the equilibrium constant for this reaction is:

$${\displaystyle K^{\ominus }={\frac {\left\{{\ce {Ag+}}_{{\ce {(aq)}}}\right\}\left\{{\ce {Cl-}}_{{\ce {(aq)}}}\right\}}{\left\{{\ce {AgCl_{(s)}}}\right\}}}=\left\{{\ce {Ag+}}_{{\ce {(aq)}}}\right\}\left\{{\ce {Cl-}}_{{\ce {(aq)}}}\right\}}$$

where ${\displaystyle K^{\ominus }}$ is the thermodynamic equilibrium constant and braces indicate activity. The activity of a pure solid is, by definition, equal to one.

When the solubility of the salt is very low the activity coefficients of the ions in solution are nearly equal to one. By setting them to be actually equal to one this expression reduces to the solubility product expression:

$${\displaystyle K_{{\ce {sp}}}=[{\ce {Ag+}}][{\ce {Cl-}}]=[{\ce {Ag+}}]^{2}=[{\ce {Cl-}}]^{2}.}$$

For 2:2 and 3:3 salts, such as CaSO₄ and FePO₄, the general expression for the solubility product is the same as for a 1:1 electrolyte

$${\displaystyle \mathrm {AB} \leftrightharpoons \mathrm {A} ^{p+}+\mathrm {B} ^{p-}}$$

${\displaystyle K_{sp}=\mathrm {[A][B]} =\mathrm {[A]^{2}} =\mathrm {[B]^{2}} }$ (electrical charges are omitted in general expressions, for simplicity of notation)

With an unsymmetrical salt like Ca(OH)₂ the solubility expression is given by

$${\displaystyle {\ce {Ca(OH)_2 <=> {Ca}^{2+}+ 2OH^-}}}$$

$${\displaystyle K_{sp}={\ce {[Ca]}}{\ce {[OH]}}^{2}}$$

Since the concentration of hydroxide ions is twice the concentration of calcium ions this reduces to ${\displaystyle \mathrm {K_{sp}=4[Ca]^{3}} }$

In general, with the chemical equilibrium

$${\displaystyle {\ce {A}}_{p}{\ce {B}}_{q}~{\ce {\leftrightharpoons }}~p{\ce {A}}^{n+}+q{\ce {B}}^{m-}}$$

$${\displaystyle {\ce {[B]}}={\frac {q}{p}}{\ce {[A]}}}$$

and the following table, showing the relationship between the solubility of a compound and the value of its solubility product, can be derived. ^[11]

Salt	p	q	Solubility, S
AgCl Ca(SO4) Fe(PO4)	1	1	√Ksp
Na2(SO4) Ca(OH)2	2 1	1 2	${\displaystyle {\sqrt[{3}]{K_{{\ce {sp}}} \over 4}}}$
Na3(PO4) FeCl3	3 1	1 3	${\displaystyle {\sqrt[{4}]{K_{{\ce {sp}}} \over 27}}}$
Al2(SO4)3 Ca3(PO4)2	2 3	3 2	${\displaystyle {\sqrt[{5}]{K_{{\ce {sp}}} \over 108}}}$
Mp(An)q	p	q	${\displaystyle {\sqrt[{p+q}]{K_{{\ce {sp}}} \over p^{p}q^{q}}}}$

Solubility products are often expressed in logarithmic form. Thus, for calcium sulfate, with K_sp = 4.93×10⁻⁵ mol² dm⁻⁶, log K_sp = −4.32. The smaller the value of K_sp, or the more negative the log value, the lower the solubility.

Some salts are not fully dissociated in solution. Examples include MgSO₄, famously discovered by Manfred Eigen to be present in seawater as both an inner sphere complex and an outer sphere complex. ^[12] The solubility of such salts is calculated by the method outlined in dissolution with reaction.

Hydroxides

The solubility product for the hydroxide of a metal ion, Mⁿ⁺, is usually defined, as follows:

$${\displaystyle \mathrm {M(OH)_{n}\leftrightharpoons \mathrm {M^{n+}+nOH^{-}} } }$$

$${\displaystyle K_{sp}=\mathrm {[M^{n+}][OH^{-}]^{n}} }$$

However, general-purpose computer programs are designed to use hydrogen ion concentrations with the alternative definitions.

$${\displaystyle \mathrm {M(OH)_{n}+nH^{+}\leftrightharpoons M^{n+}+nH_{2}O} }$$

$${\displaystyle K_{\text{sp}}^{*}=\mathrm {[M^{n+}][H^{+}]^{-n}} }$$

For hydroxides, solubility products are often given in a modified form, K*_sp, using hydrogen ion concentration in place of hydroxide ion concentration. The two values are related by the self-ionization constant for water, K_w. ^[13]

$${\displaystyle K_{\mathrm {w} }=[\mathrm {H^{+}} ][\mathrm {OH^{-}} ]}$$

$${\displaystyle K_{\text{sp}}^{*}={\frac {K_{\text{sp}}}{(K_{\text{w}})^{n}}}}$$

$${\displaystyle \log K_{\text{sp}}^{*}=\log K_{\text{sp}}-n\log K_{\text{w}}}$$

For example, at ambient temperature, for calcium hydroxide, Ca(OH)₂, lg K_sp is ca. −5 and lg K*_sp ≈ −5 + 2 × 14 ≈ 23.

Dissolution with reaction

When a concentrated solution of ammonia is added to a suspension of silver chloride dissolution occurs because a complex of Ag is formed

A typical reaction with dissolution involves a weak base, B, dissolving in an acidic aqueous solution.

$${\displaystyle \mathrm {B} \mathrm {(s)} +\mathrm {H} ^{+}\mathrm {(aq)} \leftrightharpoons \mathrm {BH} ^{+}(\mathrm {aq)} }$$

This reaction is very important for pharmaceutical products. ^[14] Dissolution of weak acids in alkaline media is similarly important.

$${\displaystyle \mathrm {HA(s)+OH^{-}(aq)\leftrightharpoons A^{-}(aq)+H_{2}O} }$$

The uncharged molecule usually has lower solubility than the ionic form, so solubility depends on pH and the acid dissociation constant of the solute. The term "intrinsic solubility" is used to describe the solubility of the un-ionized form in the absence of acid or alkali.

Leaching of aluminium salts from rocks and soil by acid rain is another example of dissolution with reaction: alumino-silicates are bases which react with the acid to form soluble species, such as Al³⁺(aq).

Formation of a chemical complex may also change solubility. A well-known example is the addition of a concentrated solution of ammonia to a suspension of silver chloride, in which dissolution is favoured by the formation of an ammine complex.

$${\displaystyle \mathrm {AgCl(s)+2NH_{3}(aq)\leftrightharpoons [Ag(NH_{3})_{2}]^{+}(aq)+Cl^{-}(aq)} }$$

When sufficient ammonia is added to a suspension of silver chloride, the solid dissolves. The addition of water softeners to washing powders to inhibit the formation of soap scum provides an example of practical importance.

Experimental determination

The determination of solubility is fraught with difficulties. ^[6] First and foremost is the difficulty in establishing that the system is in equilibrium at the chosen temperature. This is because both precipitation and dissolution reactions may be extremely slow. If the process is very slow solvent evaporation may be an issue. Supersaturation may occur. With very insoluble substances, the concentrations in solution are very low and difficult to determine. The methods used fall broadly into two categories, static and dynamic.

Static methods

In static methods a mixture is brought to equilibrium and the concentration of a species in the solution phase is determined by chemical analysis. This usually requires separation of the solid and solution phases. In order to do this the equilibration and separation should be performed in a thermostatted room. ^[15] Very low concentrations can be measured if a radioactive tracer is incorporated in the solid phase.

A variation of the static method is to add a solution of the substance in a non-aqueous solvent, such as dimethyl sulfoxide, to an aqueous buffer mixture. ^[16] Immediate precipitation may occur giving a cloudy mixture. The solubility measured for such a mixture is known as "kinetic solubility". The cloudiness is due to the fact that the precipitate particles are very small resulting in Tyndall scattering. In fact the particles are so small that the particle size effect comes into play and kinetic solubility is often greater than equilibrium solubility. Over time the cloudiness will disappear as the size of the crystallites increases, and eventually equilibrium will be reached in a process known as precipitate ageing. ^[17]

Dynamic methods

Solubility values of organic acids, bases, and ampholytes of pharmaceutical interest may be obtained by a process called "Chasing equilibrium solubility". ^[18] In this procedure, a quantity of substance is first dissolved at a pH where it exists predominantly in its ionized form and then a precipitate of the neutral (un-ionized) species is formed by changing the pH. Subsequently, the rate of change of pH due to precipitation or dissolution is monitored and strong acid and base titrant are added to adjust the pH to discover the equilibrium conditions when the two rates are equal. The advantage of this method is that it is relatively fast as the quantity of precipitate formed is quite small. However, the performance of the method may be affected by the formation supersaturated solutions.

See also

• Solubility table: A table of solubilities of mostly inorganic salts at temperatures between 0 and 100 °C.
• Solvent models

References

1. Pauling, Linus (1970). General Chemistry. Dover Publishing. p. 450.
2. Linke, W.F.; Seidell, A. (1965). Solubilities of Inorganic and Metal Organic Compounds (4th ed.). Van Nostrand. ISBN   0-8412-0097-1.
3. Kenneth Denbigh, The Principles of Chemical Equilibrium, 1957, p. 257
4. Peter Atkins, Physical Chemistry, p. 153 (8th edition)
5. Housecroft, C. E.; Sharpe, A. G. (2008). Inorganic Chemistry (3rd ed.). Prentice Hall. ISBN   978-0-13-175553-6. Section 6.10.
6. Hefter, G. T.; Tomkins, R. P. T., eds. (2003). The Experimental Determination of Solubilities. Wiley-Blackwell. ISBN   0-471-49708-8.
7. Mendham, J.; Denney, R. C.; Barnes, J. D.; Thomas, M. J. K. (2000), Vogel's Quantitative Chemical Analysis (6th ed.), New York: Prentice Hall, ISBN   0-582-22628-7 Section 2.14
8. Hsieh, Yi-Ling; Ilevbare, Grace A.; Van Eerdenbrugh, Bernard; Box, Karl J.; Sanchez-Felix, Manuel Vincente; Taylor, Lynne S. (2012-05-12). "pH-Induced Precipitation Behavior of Weakly Basic Compounds: Determination of Extent and Duration of Supersaturation Using Potentiometric Titration and Correlation to Solid State Properties". Pharmaceutical Research. 29 (10): 2738–2753. doi:10.1007/s11095-012-0759-8. ISSN   0724-8741. PMID   22580905. S2CID   15502736.
9. Dengale, Swapnil Jayant; Grohganz, Holger; Rades, Thomas; Löbmann, Korbinian (May 2016). "Recent advances in co-amorphous drug formulations". Advanced Drug Delivery Reviews. 100: 116–125. doi:10.1016/j.addr.2015.12.009. ISSN   0169-409X. PMID   26805787.
10. Gutman, E. M. (1994). Mechanochemistry of Solid Surfaces. World Scientific Publishing.
11. Skoog, Douglas A; West, Donald M; Holler, F James (2004). "9B-5". Fundamentals of Analytical Chemistry (8th ed.). Brooks/Cole. pp. 238–242. ISBN   0030355230.
12. Eigen, Manfred (1967). "Nobel lecture" (PDF). Nobel Prize.
13. Baes, C. F.; Mesmer, R. E. (1976). The Hydrolysis of Cations. New York: Wiley.
14. Payghan, Santosh (2008). "Potential Of Solubility In Drug Discovery And development". Pharminfo.net. Archived from the original on March 30, 2010. Retrieved 5 July 2010.
15. Rossotti, F. J. C.; Rossotti, H. (1961). "Chapter 9: Solubility". The Determination of Stability Constants . McGraw-Hill.
16. Aqueous solubility measurement – kinetic vs. thermodynamic methods Archived July 11, 2009, at the Wayback Machine
17. Mendham, J.; Denney, R. C.; Barnes, J. D.; Thomas, M. J. K. (2000), Vogel's Quantitative Chemical Analysis (6th ed.), New York: Prentice Hall, ISBN   0-582-22628-7 Chapter 11: Gravimetric analysis
18. Stuart, M.; Box, K. (2005). "Chasing Equilibrium: Measuring the Intrinsic Solubility of Weak Acids and Bases". Analytical Chemistry. 77 (4): 983–990. doi:10.1021/ac048767n. PMID   15858976.

External links

• Housecroft, C. E.; Sharpe, A. G. (2008). Inorganic Chemistry (3rd ed.). Prentice Hall. ISBN   978-0-13-175553-6. Section 6.9: Solubilities of ionic salts. Includes a discussion of the thermodynamics of dissolution.
• IUPAC–NIST solubility database
• Solubility products of simple inorganic compounds
• Solvent activity along a saturation line and solubility
• Solubility challenge: Predict solubilities from a data base of 100 molecules. The database, of mostly compounds of pharmaceutical interest, is available at One hundred molecules with solubilities (Text file, tab separated).

A number of computer programs are available to do the calculations. They include:

• CHEMEQL: A comprehensive computer program for the calculation of thermodynamic equilibrium concentrations of species in homogeneous and heterogeneous systems. Many geochemical applications.
• JESS: All types of chemical equilibria can be modelled including protonation, complex formation, redox, solubility and adsorption interactions. Includes an extensive database.
• MINEQL+: A chemical equilibrium modeling system for aqueous systems. Handles a wide range of pH, redox, solubility and sorption scenarios.
• PHREEQC: USGS software designed to perform a wide variety of low-temperature aqueous geochemical calculations, including reactive transport in one dimension.
• MINTEQ: A chemical equilibrium model for the calculation of metal speciation, solubility equilibria etc. for natural waters.
• WinSGW: A Windows version of the SOLGASWATER computer program.