Reference articles on history, science, culture and more
Encyclopedia

Bjerrum plot

Graph of polyprotic acid concentration as a function of pH

Image credit is listed at the end of this article.

A Bjerrum plot (named after Niels Bjerrum), sometimes also known as a Sillén diagram (after Lars Gunnar Sillén), or a Hägg diagram (after Gunnar Hägg) is a graph of the concentrations of the different species of a polyprotic acid in a solution, as a function of pH, when the solution is at equilibrium. Due to the many orders of magnitude spanned by the concentrations, they are commonly plotted on a logarithmic scale. Sometimes the ratios of the concentrations are plotted rather than the actual concentrations. Occasionally H+ and OH are also plotted.

Most often, the carbonate system is plotted, where the polyprotic acid is carbonic acid (a diprotic acid), and the different species are dissolved carbon dioxide, carbonic acid, bicarbonate, and carbonate. In acidic conditions, the dominant form is CO2; in basic (alkaline) conditions, the dominant form is CO2−
3
; and in between, the dominant form is HCO
3
. At every pH, the concentration of carbonic acid is assumed to be negligible compared to the concentration of dissolved CO
2
, and so is often omitted from Bjerrum plots. These plots are very helpful in solution chemistry and natural water chemistry. In the example given here, it illustrates the response of seawater pH and carbonate speciation due to the input of man-made CO
2
emission by the fossil fuel combustion.

The Bjerrum plots for other polyprotic acids, including silicic, boric, sulfuric and phosphoric acids, are other commonly used examples.

01Bjerrum plot equations for carbonate system

If carbon dioxide, carbonic acid, hydrogen ions, bicarbonate and carbonate are all dissolved in water, and at chemical equilibrium, their equilibrium concentrations are often assumed to be given by:

{\begin{aligned}[]\left[{\textrm {CO}}_{2}\right]_{\text{eq}}&={\frac {\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}+K_{1}\left[{\textrm {H}}^{+}\right]_{\text{eq}}+K_{1}K_{2}}}\times {\textrm {DIC}},\\[3pt]\left[{\textrm {HCO}}_{3}^{-}\right]_{\text{eq}}&={\frac {K_{1}\left[{\textrm {H}}^{+}\right]_{\text{eq}}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}+K_{1}\left[{\textrm {H}}^{+}\right]_{\text{eq}}+K_{1}K_{2}}}\times {\textrm {DIC}},\\[3pt]\left[{\textrm {CO}}_{3}^{2-}\right]_{\text{eq}}&={\frac {K_{1}K_{2}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}+K_{1}\left[{\textrm {H}}^{+}\right]_{\text{eq}}+K_{1}K_{2}}}\times {\textrm {DIC}},\end{aligned}}

where the subscript 'eq' denotes that these are equilibrium concentrations, K1 is the equilibrium constant for the reaction CO
2
+ H
2
O
H+ + HCO
3
(i.e. the first acid dissociation constant for carbonic acid), K2 is the equilibrium constant for the reaction HCO
3
H+ + CO2−
3
(i.e. the second acid dissociation constant for carbonic acid), and DIC is the (unchanging) total concentration of dissolved inorganic carbon in the system, i.e. [CO2] + [HCO
3
] + [CO2−
3
]. K1, K2 and DIC each have units of a concentration, e.g. mol/L.

A Bjerrum plot is obtained by using these three equations to plot these three species against pH = −log10 [H+]eq, for given K1, K2 and DIC. The fractions in these equations give the three species' relative proportions, and so if DIC is unknown, or the actual concentrations are unimportant, these proportions may be plotted instead.

These three equations show that the curves for CO2 and HCO
3
intersect at [H+]eq = K1, and the curves for HCO
3
and CO2−
3
intersect at [H+]eq = K2. Therefore, the values of K1 and K2 that were used to create a given Bjerrum plot can easily be found from that plot, by reading off the concentrations at these points of intersection. An example with linear Y axis is shown in the accompanying graph. The values of K1 and K2, and therefore the curves in the Bjerrum plot, vary substantially with temperature and salinity.

02Chemical and mathematical derivation of Bjerrum plot equations for carbonate system

Suppose that the reactions between carbon dioxide, hydrogen ions, bicarbonate and carbonate ions, all dissolved in water, are as follows:

CO
2
+ H
2
O
⇌ H+ + HCO
3
1
HCO
3
⇌ H+ + CO2−
3
2

Note that reaction 1 is actually the combination of two elementary reactions:

CO
2
+ H
2
O
H
2
CO
3
H+ + HCO
3

Assuming the mass action law applies to these two reactions, that water is abundant, and that the different chemical species are always well-mixed, their rate equations are

{\begin{aligned}{\frac {{\textrm {d}}\left[{\textrm {CO}}_{2}\right]}{{\textrm {d}}t}}&=-k_{1}\left[{\textrm {CO}}_{2}\right]+k_{-1}\left[{\textrm {H}}^{+}\right]\left[{\textrm {HCO}}_{3}^{-}\right],\\{\frac {{\textrm {d}}\left[{\textrm {H}}^{+}\right]}{{\textrm {d}}t}}&=k_{1}\left[{\textrm {CO}}_{2}\right]-k_{-1}\left[{\textrm {H}}^{+}\right]\left[{\textrm {HCO}}_{3}^{-}\right]+k_{2}\left[{\textrm {HCO}}_{3}^{-}\right]-k_{-2}\left[{\textrm {H}}^{+}\right]\left[{\textrm {CO}}_{3}^{2-}\right],\\{\frac {{\textrm {d}}\left[{\textrm {HCO}}_{3}^{-}\right]}{{\textrm {d}}t}}&=k_{1}\left[{\textrm {CO}}_{2}\right]-k_{-1}\left[{\textrm {H}}^{+}\right]\left[{\textrm {HCO}}_{3}^{-}\right]-k_{2}\left[{\textrm {HCO}}_{3}^{-}\right]+k_{-2}\left[{\textrm {H}}^{+}\right]\left[{\textrm {CO}}_{3}^{2-}\right],\\{\frac {{\textrm {d}}\left[{\textrm {CO}}_{3}^{2-}\right]}{{\textrm {d}}t}}&=k_{2}\left[{\textrm {HCO}}_{3}^{-}\right]-k_{-2}\left[{\textrm {H}}^{+}\right]\left[{\textrm {CO}}_{3}^{2-}\right]\end{aligned}}

where [] denotes concentration, t is time, and K1 and k−1 are appropriate proportionality constants for reaction 1, called respectively the forwards and reverse rate constants for this reaction. (Similarly K2 and k−2 for reaction 2.)

At any equilibrium, the concentrations are unchanging, hence the left hand sides of these equations are zero. Then, from the first of these four equations, the ratio of reaction 1's rate constants equals the ratio of its equilibrium concentrations, and this ratio, called K1, is called the equilibrium constant for reaction 1, i.e.

K_{1}={\frac {k_{1}}{k_{-1}}}={\frac {[{\textrm {H}}^{+}]_{\text{eq}}[{\textrm {HCO}}_{3}^{-}]_{\text{eq}}}{[{\textrm {CO}}_{2}]_{\text{eq}}}} 3

where the subscript 'eq' denotes that these are equilibrium concentrations.

Similarly, from the fourth equation for the equilibrium constant K2 for reaction 2,

K_{2}={\frac {k_{2}}{k_{-2}}}={\frac {\left[{\textrm {H}}^{+}\right]_{\text{eq}}\left[{\textrm {CO}}_{3}^{2-}\right]_{\text{eq}}}{\left[{\textrm {HCO}}_{3}^{-}\right]_{\text{eq}}}} 4

Rearranging 3 gives

\left[{\textrm {HCO}}_{3}^{-}\right]_{\text{eq}}={\frac {K_{1}\left[{\textrm {CO}}_{2}\right]_{\text{eq}}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}}}         5

and rearranging 4, then substituting in 5, gives

\left[{\textrm {CO}}_{3}^{2-}\right]_{\text{eq}}={\frac {K_{2}\left[{\textrm {HCO}}_{3}^{-}\right]_{\text{eq}}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}}}={\frac {K_{1}K_{2}\left[{\textrm {CO}}_{2}\right]_{\text{eq}}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}}}         6

The total concentration of dissolved inorganic carbon in the system is given by substituting in 5 and 6:

{\begin{aligned}{\textrm {DIC}}&=\left[{\textrm {CO}}_{2}\right]+\left[{\textrm {HCO}}_{3}^{-}\right]+\left[{\textrm {CO}}_{3}^{2-}\right]\\&=\left[{\textrm {CO}}_{2}\right]_{\text{eq}}\left(1+{\frac {K_{1}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}}}+{\frac {K_{1}K_{2}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}}}\right)\\&=\left[{\textrm {CO}}_{2}\right]_{\text{eq}}\left({\frac {\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}+K_{1}\left[{\textrm {H}}^{+}\right]_{\text{eq}}+K_{1}K_{2}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}}}\right)\end{aligned}}

Re-arranging this gives the equation for CO
2
:

\left[{\textrm {CO}}_{2}\right]_{\text{eq}}={\frac {\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}}{\left[{\textrm {H}}^{+}\right]_{\text{eq}}^{2}+K_{1}\left[{\textrm {H}}^{+}\right]_{\text{eq}}+K_{1}K_{2}}}\times {\textrm {DIC}}

The equations for HCO
3
and CO2−
3
are obtained by substituting this into 5 and 6.

Watch videos about Bjerrum plotExplainers and documentaries on YouTube (opens in a new tab)

Sources and credits

This article is adapted from the Wikipedia article Bjerrum plot, written by its contributors and licensed under CC BY-SA 4.0. Fathomly has changed the layout, removed citation markers, navigation and maintenance notices, and adjusted punctuation. This adapted version is shared under the same license. For references, see the original article.

Images, from Wikimedia Commons:

Fathomly is not affiliated with or endorsed by the Wikimedia Foundation. Spotted a problem? Tell us.

Continue exploring

Related topics

Charlot equation

The Charlot equation, named after Gaston Charlot, is used in analytical chemistry to relate the hydrogen ion concentration, and therefore the pH, with the formal analytical concentration of an acid and its conjugate base. It can be used for computing the pH of buffer solutions when the approximations of the Henderson-Hasselbalch equation break down.

Gran plot

A Gran plot is a common means of standardizing a titrate or titrant by estimating the equivalence volume or end point in a strong acid-strong base titration or in a potentiometric titration. Such plots have been also used to calibrate glass electrodes, to estimate the carbonate content of aqueous solutions, and to estimate the Ka values (acid dissociation constants) of weak acids and bases from titration data.

Henderson-Hasselbalch equation

In chemistry and biochemistry, the pH of weakly acidic chemical solutions can be estimated using the Henderson-Hasselbalch Equation: pH = p K a + log 10 ⁡ {\displaystyle {\ce {pH}}={\ce {p}}K_{{\ce {a}}}+\log _{10}\left({\frac {[{\ce {Base}}]}{[{\ce {Acid}}]}}\right)} The equation relates the pH of the weak acid to the numerical value of the acid dissociation constant, Ka, of the acid, and the ratio of the concentrations of the acid and its conjugate base. Acid-base Equilibrium Reaction H A ( a c i d ) ⇌ A − ( b a s e ) + H + {\displaystyle \mathrm {{\underset {(acid)}{HA}}\rightleftharpoons {\underset {(base)}{A^{-}}}+H^{+}} } The Henderson-Hasselbalch equation is often used for estimating the pH of buffer solutions by approximating the actual concentration ratio as the ratio of the analytical concentrations of the acid and of a salt, MA. It is also useful for determining the volumes of the reagents needed before preparing buffer solutions, which prevents unnecessary waste of chemical reagents that may need to be further neutralized by even more reagents before they are safe to expose.