Systematic Treatment of Equilibrium

Clinical laboratories computing blood pH and buffer capacity use the exact bicarbonate/carbonic-acid/phosphate simultaneous system — the Henderson–Hasselbalch equation used at the bedside is a simplification of precisely the kind of system built here, valid only because physiological pH sits far from extreme p$K_a$ values, a fact verified by the same systematic setup. Environmental chemists modeling ocean acidification solve carbonate-system equations with four to six coupled species (dissolved CO₂, carbonic acid, bicarbonate, carbonate, borate species, water) simultaneously, using software (CO2SYS and similar) that implements exactly the mass-balance/charge-balance/equilibrium-constant scheme described here, generalized to dozens of species for full seawater speciation. Pharmaceutical formulation scientists use systematic equilibrium setups to predict drug solubility and ionization state across the gastrointestinal pH range, where a single active pharmaceutical ingredient may have three or four ionizable groups whose simultaneous protonation states determine bioavailability. Complexometric titration methods (EDTA titrations for water hardness, metal assays) require simultaneous treatment of metal-ligand formation equilibria, competing proton equilibria on the ligand, and any competing hydroxide or auxiliary-complexing-agent equilibria — precisely the multi-equilibrium census demonstrated in the silver-ammonia example. Geochemists modeling mineral solubility and groundwater speciation, and electrochemists designing battery electrolytes with multiple redox-active and acid-base-active species, all rely on the same disciplined enumeration: list every species, write every independent equilibrium, write every mass balance, write the one charge balance, count degrees of freedom, and only then solve — whether by hand with justified approximations or numerically with iterative solvers, the topic that follows directly from this framework.

Ask the tutor: How do I decide, once the system is fully assembled, which terms in the charge balance are actually safe to drop before solving?

Consider dissolving $C_0$ mol of acetic acid in water to make one liter of solution. A beginner's instinct is to write

and then assume $[\text{H}^+] \approx [\text{OAc}^-]$, solving a quadratic. This works for acetic acid alone, but the assumption $[\text{H}^+] \approx [\text{OAc}^-]$ is not an axiom — it is a consequence of the charge balance, valid only because water's own ions are negligible at that particular concentration. The systematic treatment reveals precisely when that shortcut is legitimate and when it silently fails (e.g., extremely dilute acid, where water autoionization is not negligible). The physical picture is a census of every species truly present, followed by a matching census of every independent constraint on those species. Species arise from: (1) whatever solute(s) were dissolved, broken into all their acid/base/conjugate forms; (2) water itself, plus $\text{H}_3\text{O}^+$ and $\text{OH}^-$ from its autoionization; (3) any counter-ions from salts, which participate only in charge balance (they are "spectators" chemically but not electrically). Constraints arise from: equilibrium constants (one per independent reaction), mass balances (one per conserved atom-group, i.e., one per distinct "total" you introduced when dissolving reagents), and exactly one charge balance for the whole solution. A powerful way to see why the charge balance is not redundant with mass balance is to recognize that mass balance conserves atoms regardless of what species they are in, while charge balance conserves net electrostatic neutrality of the bulk solution — a completely independent physical fact (solutions are electroneutral because bulk charge separation costs enormous electrostatic energy). These are physically different constraints, and both are needed.

Notice in the figure that six chemically distinct dissolved unknowns appear: $\text{H}_2\text{CO}_3$, $\text{HCO}_3^-$, $\text{CO}_3^{2-}$, $\text{H}_3\text{O}^+$, $\text{OH}^-$ — that is five, plus if a counter-ion (say $\text{Na}^+$ from a $\text{Na}_2\text{CO}_3$ source) is present, six. The equation count must match exactly.

Mass balance for a conserved group introduced at formal (analytical) concentration $C_0$:

where:

• $C_0$ is the total formal concentration of the atom/group as originally weighed in,
• the sum runs over every species containing that atom/group,
• $n_i$ is the number of that atom/group per formula unit of species $i$ (usually 1, but e.g. 2 for a group appearing twice in one species).

Reading: mass balance says atoms do not disappear — they merely redistribute among all the protonation/complexation states available.

Charge balance (electroneutrality):

where:

• $z_j^+$ and $z_k^-$ are the magnitudes of charge on each cationic/anionic species,
• every ionic species in solution — including $\text{H}_3\text{O}^+$ and $\text{OH}^-$ — must appear,
• the equation is a single scalar equality, valid regardless of how many equilibria are occurring simultaneously.

Reading: bulk solution has zero net charge; multiply each concentration by its charge magnitude before summing, because a doubly charged ion contributes twice the electrostatic "weight" of a singly charged one at the same molarity.

Equilibrium-constant expressions (one per independent reaction):

where:

• $\nu_{i,m}$ is the stoichiometric coefficient (signed: positive for products, negative for reactants) of species $m$ in reaction $i$,
• $K_i$ is the thermodynamic (or, at fixed ionic strength, the conditional) equilibrium constant for reaction $i$,
• "independent" means reaction $i$ cannot be obtained by linear combination of the others already listed.

Degrees-of-freedom rule, the master bookkeeping check:

where:

• $N_{\text{species}}$ counts every distinct dissolved chemical species (ions and neutrals, including $\text{H}_3\text{O}^+$, $\text{OH}^-$, but not solid water/solvent activity, which is fixed at 1),
• $N_{\text{equilibria}}$ counts independent equilibrium reactions (including water autoionization),
• $N_{\text{mass balances}}$ counts distinct conserved totals introduced by dissolution (one per reagent "family," e.g., one for all carbonate-containing species, one for all sodium-containing species),
• the "+1" is because there is always exactly one charge balance for the whole solution, never more.

Reading: if the left and right sides of this counting equation do not match, either a species was missed, an equilibrium was double-counted, or a mass balance was forgotten — the algebra cannot proceed correctly until they match.

Example 1 (foundational) — full setup for a weak monoprotic acid, no solving. A chemist dissolves $0.0742\ \text{M}$ chloroacetic acid (ClCH$_2$COOH, abbreviated HA, $K_a = 1.36\times10^{-3}$) in water. Set up the complete simultaneous system. Species: HA, A⁻, H₃O⁺, OH⁻ → $N_{\text{species}} = 4$. Equilibria (2):

Mass balance (1):

Charge balance (1):

Degrees-of-freedom check: $N_{\text{species}} = 4$; RHS $= 2 + 1 + 1 = 4$. ✓ Exactly determined. This is the complete, exact system — nothing approximated yet. A later approximation step (recognizing $K_a = 1.36\times10^{-3}$ is not tiny relative to $C_0 = 0.0742$, so the "$x \ll C_0$" shortcut is questionable here, requiring the full quadratic $[\text{H}_3\text{O}^+]^2/(C_0-[\text{H}_3\text{O}^+]) \approx K_a$) belongs to the solving stage, not this setup stage — but note the setup already tells you which equations that quadratic came from and why $[\text{OH}^-]$ was dropped from the charge balance (because $[\text{H}_3\text{O}^+]$ here will be of order $\sqrt{K_aC_0} \approx \sqrt{(1.36\times10^{-3})(0.0742)} \approx 1.0\times10^{-2}\ \text{M}$, some $10^{10}$-fold larger than the $[\text{OH}^-] \approx K_w/[\text{H}_3\text{O}^+] \approx 10^{-12}\ \text{M}$ it would introduce).

Example 2 (applied) — full system for a buffer with an added common ion, EDTA-adjacent case. An analyst prepares a solution that is $0.0500\ \text{M}$ in NaH$_2$PO$_4$ and $0.0300\ \text{M}$ in Na$_2$HPO$_4$ (a phosphate buffer, relevant to controlling pH during a complexometric titration with EDTA, since many metal–EDTA formation constants are pH-dependent). Set up the complete system, treating phosphoric acid's three dissociations explicitly. Species: H₃PO₄, H₂PO₄⁻, HPO₄²⁻, PO₄³⁻, Na⁺, H₃O⁺, OH⁻ → $N_{\text{species}} = 7$. Equilibria (4):

Mass balances (2): Phosphate total is the sum of both salts' contributions, since both sources feed the same conjugate family:

Sodium comes from $\text{NaH}_2\text{PO}_4$ (1 Na per formula unit) and $\text{Na}_2\text{HPO}_4$ (2 Na per formula unit):

Charge balance (1):

Degrees-of-freedom check: $N_{\text{species}} = 7$; RHS $= 4\ (\text{equilibria}) + 2\ (\text{mass balances}) + 1 = 7$. ✓ This system is now fully posed for exact numerical solution (e.g., by Newton–Raphson on the resulting polynomial in $[\text{H}_3\text{O}^+]$), or for the standard Henderson–Hasselbalch-style approximation once it is verified — from this system, not assumed — that $[\text{H}_3\text{PO}_4]$ and $[\text{PO}_4^{3-}]$ are negligible compared to $[\text{H}_2\text{PO}_4^-]$ and $[\text{HPO}_4^{2-}]$ near pH $\approx \text{p}K_{a2} = 7.20$.

Example 3 (exam/competition level) — mixed-ligand system: silver, ammonia, and chloride. A solution is prepared by mixing AgNO$_3$ ($C_{\text{Ag}} = 2.00\times10^{-3}\ \text{M}$ formal), NH$_3$ ($C_{\text{N}} = 0.150\ \text{M}$ formal, added as free base), and NaCl ($C_{\text{Cl}} = 1.00\times10^{-3}\ \text{M}$ formal), forming silver–ammine complexes and possibly AgCl(s). Assume for this setup that all species remain dissolved (no precipitate — this is checked separately via $K_{sp}$ comparison, foreshadowing solubility-equilibrium problems). Set up the full system. Species: Ag⁺, Ag(NH₃)⁺, Ag(NH₃)₂⁺, NH₃, NH₄⁺, Cl⁻, Na⁺, H₃O⁺, OH⁻ → $N_{\text{species}} = 9$. Equilibria (5): stepwise silver–ammine formation, ammonia protonation, and water autoionization:

Chloride is here assumed to form no significant complex with silver under these dilute, ammonia-rich conditions (a claim that would itself be checked quantitatively in a full solubility treatment), so it contributes no equilibrium expression beyond appearing as a free ion — this is a modeling decision stated explicitly, not hidden. That is only 4 independent equilibria listed; a possible fifth (direct reaction of $\text{Ag}^+$ with $\text{OH}^-$ to precipitate $\text{Ag}_2\text{O}$/AgOH) is excluded by the problem's stated assumption of no precipitate, which is exactly the kind of assumption a systematic setup forces you to state rather than bury. Mass balances (3):

Notice the nitrogen mass balance: ammonia bound in the complexes must be counted in the nitrogen total too, with the correct stoichiometric weight (1 or 2 NH₃ per complex ion) — a subtlety easy to miss if mass balances are written carelessly. Charge balance (1):

Here $[\text{Na}^+] = 1.00\times10^{-3}\ \text{M}$ (from NaCl) and $[\text{NO}_3^-] = 2.00\times10^{-3}\ \text{M}$ (from AgNO₃) are fixed spectator concentrations, not unknowns requiring their own equilibrium — but they are still unknowns in the counting sense and must appear as terms determined by separate trivial mass balances ($[\text{Na}^+] = C_{\text{Cl}}$ since NaCl dissolves 1:1, and $[\text{NO}_3^-] = C_{\text{Ag}}$ since AgNO₃ dissolves 1:1). Full degrees-of-freedom check: including Na⁺ and NO₃⁻ as unknowns, $N_{\text{species}} = 11$ (Ag⁺, Ag(NH₃)⁺, Ag(NH₃)₂²⁺-labeled-as-1+, NH₃, NH₄⁺, Cl⁻, Na⁺, NO₃⁻, H₃O⁺, OH⁻ — that's 10; recount: Ag⁺(1), Ag(NH₃)⁺(2), Ag(NH₃)₂⁺(3), NH₃(4), NH₄⁺(5), Cl⁻(6), Na⁺(7), NO₃⁻(8), H₃O⁺(9), OH⁻(10) → $N_{\text{species}}=10$). Equations: 4 equilibria + 5 mass balances ($C_{\text{Ag}}$, $C_{\text{N}}$, $C_{\text{Cl}}$, $[\text{Na}^+]=C_{\text{Cl}}$, $[\text{NO}_3^-]=C_{\text{Ag}}$) + 1 charge balance $= 4+5+1 = 10$. ✓ The count closes exactly once every trivial spectator relation is included as its own mass balance — a common source of "off by one" errors in exam settings, resolved only by disciplined enumeration.

The systematic algebraic treatment of multi-equilibrium systems, as taught in modern analytical chemistry, is most closely associated with the work of James N. Butler, whose 1964 text Ionic Equilibrium: A Mathematical Approach formalized the mass-balance/charge-balance/equilibrium-constant framework into the general procedure used in essentially all modern analytical chemistry textbooks (Skoog, Harris, Christian). Butler's contribution was not a new equilibrium constant or a new reaction, but a method — insisting that every species be enumerated, every conservation law written explicitly, and the resulting system solved by systematic iteration (he was an early advocate of numerical, computer-assisted solution of the exact equations rather than case-by-case approximations), decades before this became routine. The deeper conceptual roots go back to Niels Bjerrum, who in the 1910s and 1920s developed stepwise formation constants for polyprotic acids and metal complexes and recognized that overlapping equilibria required simultaneous treatment rather than sequential, isolated calculations — the same insight this article formalizes.