Henderson–Hasselbalch Equation

The Henderson–Hasselbalch equation is a fundamental logarithmic equation in analytical chemistry, biochemistry, and human physiology that relates the pH of an aqueous solution to the negative logarithm of the acid dissociation constant (p*K*_a) of a weak acid and the base-10 logarithm of the concentration ratio of its deprotonated conjugate base ([A⁻]) to its protonated weak acid ([HA]):

`` pH = pKₐ + log₁₀( [A⁻] / [HA] ) ``

First formulated in non-logarithmic mass-action form in 1908 by the American physiologist Lawrence Joseph Henderson to describe how carbonic acid and bicarbonate salts stabilize blood neutrality, the relation was recast into logarithmic pH notation in 1916 by the Danish physician and chemist Karl Albert Hasselbalch (following Søren Sørensen's 1909 introduction of the pH scale). Grounded in Brønsted–Lowry acid–base theory, the Henderson–Hasselbalch equation provides the quantitative framework for designing laboratory buffer solutions, predicting the ionization charge state of amino-acid side chains in proteins and enzymes, modeling drug transmembrane diffusion across lipid bilayers, and diagnosing clinical respiratory and metabolic acid–base disorders.

Derivation from the Equilibrium Dissociation Constant

Consider a weak Brønsted–Lowry monoprotic acid HA (such as acetic acid, CH₃COOH, or an ammonium ion, R–NH₃⁺) undergoing reversible proton transfer in water:

`` HA (aq) + H₂O (l) ⇌ H₃O⁺ (aq) + A⁻ (aq) ``

At thermodynamic equilibrium, the apparent molar acid dissociation constant (K_a) is defined by the law of mass action (incorporating the constant activity of bulk liquid water into K_a and writing H₃O⁺ simply as H⁺):

`` Kₐ = ( [H⁺] · [A⁻] ) / [HA] ``

Solving this algebraic expression for the hydronium ion concentration yields [H⁺] = K_a · ([HA] / [A⁻]). Taking the negative base-10 logarithm (−log₁₀) of both sides gives:

`` -log₁₀[H⁺] = -log₁₀(Kₐ) - log₁₀( [HA] / [A⁻] ) ``

Substituting the definitions pH = −log₁₀[H⁺] and p*K*_a = −log₁₀(K_a)—and inverting the numerator and denominator inside the logarithm to flip its sign—produces the canonical Henderson–Hasselbalch equation:

`` pH = pKₐ + log₁₀( [Conjugate Base] / [Weak Acid] ) = pKₐ + log₁₀( [A⁻] / [HA] ) ``

From the relationship between standard Gibbs free energy and the equilibrium constant (Δ*G*° = −*RT* ln K_a = 2.303 RT p*K*_a), every 1.0 unit increase in p*K*_a corresponds to a +5.71 kJ/mol (+1.36 kcal/mol at 25 °C) increase in the free energy required for deprotonation.

Ionic Strength and Debye–Hückel Activity Coefficients

In rigorous chemical thermodynamics, equilibrium constants are defined in terms of dimensionless chemical activities (a_i = γ_i[i]/c°) rather than raw molar concentrations, where γ_i is the activity coefficient:

`` pH = pKₐ° + log₁₀( [A⁻] / [HA] ) + log₁₀( γ_A⁻ / γ_HA ) ``

In dilute solutions with ionic strength I > 0.01 M (such as intracellular fluid or blood plasma, where I ≈ 0.15 M), the extended Debye–Hückel equation shows that charged ions experience electrostatic shielding from the surrounding ionic atmosphere (γ_i < 1), which lowers the apparent operational p*K*_a' relative to the infinite-dilution thermodynamic p*K*_a° (for example, shifting the second dissociation constant of phosphoric acid from p*K*_a2° = 7.20 at zero ionic strength to an apparent p*K*_a2' = 6.80–6.86 at physiological ionic strength and 37 °C).

The Van Slyke Buffer Capacity Equation

When a strong acid (HCl) or strong base (NaOH) is added to an equimolar mixture of a weak acid HA and its conjugate base A⁻, the added H⁺ ions are consumed by A⁻ to form HA, and added OH⁻ ions react with HA to form A⁻ + H₂O. Consequently, the ratio [A⁻]/[HA] changes gradually, resisting large shifts in pH. At the midpoint of a titration curve where [A⁻] = [HA], log₁₀(1) = 0 and pH = p*K*_a. A 10-fold ratio in either direction ([A⁻]/[HA] = 10:1 or 1:10) shifts the pH by ±1.0 unit (log₁₀(10) = +1; log₁₀(0.1) = −1), defining the effective buffering range of any monoprotic buffer as pH = p*K*_a ± 1.0 (spanning 9.1% to 90.9% ionization).

In 1922, Donald Van Slyke quantified the resistance of a buffer solution of total molar concentration C_buf = [HA] + [A⁻] to pH changes via the differential buffer capacity (β = dC_b / *d*pH = −*dC*_a / *d*pH, the moles of strong base or acid per liter required to change pH by 1.0 unit):

`` β = 2.303 · ( (K_w / [H⁺]) + [H⁺] + ( C_buf · Kₐ · [H⁺] ) / ( Kₐ + [H⁺] )² ) ``

Differentiating with respect to [H⁺] proves that for a given total buffer concentration C_buf, maximum buffer capacity occurs when [H⁺] = K_a (pH = p*K*_a), where the buffer term simplifies to β_max = 0.576 · C_buf.

The Open Bicarbonate–Carbon Dioxide Buffer System

Human arterial blood plasma is maintained within a narrow pH window of 7.35 to 7.45 (average pH 7.40, corresponding to [H⁺] = 40 nmol/L), primarily through the synergistic action of the extracellular bicarbonate buffer system and intracellular hemoglobin/protein histidine residues. Dissolved carbon dioxide (CO₂(d)) undergoes reversible hydration—accelerated >10⁶-fold in erythrocytes and renal tubules by the Zn²⁺ metalloenzyme carbonic anhydrase—to form carbonic acid (H₂CO₃), which rapidly dissociates into bicarbonate (HCO₃⁻) and H⁺:

`` CO₂ (g) ⇌ CO₂ (d) + H₂O ⇌ H₂CO₃ ⇌ H⁺ + HCO₃⁻ ``

Because true undissociated H₂CO₃ exists at <0.3% of dissolved CO₂ concentration, clinical physiology combines [CO₂(d)] + [H₂CO₃] into an effective carbonic acid pool whose concentration is given by Henry's Law: [CO₂(d)] = α · P(CO₂), where α = 0.0301 mmol·L⁻¹·mmHg⁻¹ is the solubility coefficient of CO₂ in plasma at 37 °C and P(CO₂) is the arterial partial pressure of carbon dioxide (in mmHg). With an apparent combined p*K*_a' = 6.10 at 37 °C, the physiological Henderson–Hasselbalch equation for arterial blood is:

`` pH = 6.10 + log₁₀( [HCO₃⁻] / (0.03 × P_CO₂) ) ``

At normal arterial values of [HCO₃⁻] = 24 mmol/L and P(CO₂) = 40 mmHg (so that 0.03 × 40 = 1.2 mmol/L dissolved CO₂), the base-to-acid ratio is 24 / 1.2 = 20 : 1. Because log₁₀(20) = 1.30, arterial blood pH equals 6.10 + 1.30 = 7.40.

Respiratory and Metabolic Acid–Base Disturbances

Although a closed buffer with p*K*_a = 6.10 would exhibit poor buffer capacity at pH 7.40 (where the ratio is 20:1 rather than 1:1), the physiological bicarbonate buffer operates as an open thermodynamic system coupled to two organs: the lungs regulate the denominator (P(CO₂)) within minutes via alveolar ventilation, while the kidneys regulate the numerator ([HCO₃⁻]) over hours to days via proximal tubule HCO₃⁻ reabsorption and distal H⁺/NH₄⁺ excretion:

• Respiratory Acidosis (↑ P(CO₂) > 45 mmHg due to hypoventilation or COPD) lowers pH; kidneys compensate by retaining and generating HCO₃⁻.
• Respiratory Alkalosis (↓ P(CO₂) < 35 mmHg due to hyperventilation or high altitude) raises pH; kidneys compensate by excreting HCO₃⁻.
• Metabolic Acidosis (↓ [HCO₃⁻] < 22 mM due to diabetic ketoacidosis, lactic acidosis, or renal failure) lowers pH; carotid and central chemoreceptors trigger compensatory hyperventilation (Kussmaul breathing) to blow off CO₂ and lower P(CO₂).
• Metabolic Alkalosis (↑ [HCO₃⁻] > 26 mM due to gastric HCl vomiting or loop diuretics) raises pH; respiratory centers modestly decrease ventilation to raise P(CO₂).

Why is the bicarbonate buffer so effective at blood pH 7.40 even though its pKa is 6.10?

In a sealed test tube (a closed system), a buffer with p*K*_a = 6.10 loses most of its buffering capacity at pH 7.40 because adding strong acid converts HCO₃⁻ into dissolved CO₂, rapidly lowering the [HCO₃⁻]/[CO₂] ratio. In the human body (an open system), any excess dissolved CO₂ formed when HCO₃⁻ neutralizes metabolic H⁺ is immediately exhaled by the lungs, holding arterial P(CO₂) fixed near 40 mmHg (1.2 mM). This physiological clamping of the denominator multiplies the effective buffer capacity of the bicarbonate system several-fold at pH 7.40.

How does the Henderson–Hasselbalch equation predict drug absorption and ion trapping across cell membranes?

Most pharmaceutical drugs are weak organic acids (HA ⇌ H⁺ + A⁻, such as aspirin with p*K*_a ≈ 3.5) or weak organic bases (BH⁺ ⇌ H⁺ + B, such as morphine with p*K*_a ≈ 8.0). Because only the uncharged species (HA or B) readily permeates hydrophobic lipid bilayers via simple diffusion, the Henderson–Hasselbalch equation determines the uncharged fraction in each compartment. In the acidic stomach lumen (pH ≈ 1.5), aspirin exists almost exclusively in its protonated, neutral HA form ([HA]/[A⁻] = 10^3.5−1.5 = 100:1) and diffuses rapidly into mucosal cells (pH 7.4), where it deprotonates to charged A⁻ ([A⁻]/[HA] ≈ 7,900:1) and becomes "ion-trapped."

When does the Henderson–Hasselbalch approximation fail in analytical chemistry?

In deriving the standard laboratory form of the Henderson–Hasselbalch equation, one assumes that the initial analytical concentrations placed into solution (C_HA and C_A−) remain unchanged at equilibrium ([HA] ≈ C_HA and [A⁻] ≈ C_A−). This approximation breaks down when (1) the acid is very strong or moderately strong (K_a > 10⁻², p*K*_a < 2) or very weak (p*K*_a > 12), or (2) the total buffer concentration is highly dilute (C_buf < 10⁻⁴ M or comparable to [H⁺] and [OH⁻]), requiring exact solution of the cubic charge-balance and water-autoionization (K_w) equations.

1. Henderson, L. J. (1908). "Concerning the relationship between the strength of acids and their capacity to preserve neutrality." American Journal of Physiology, 21(2), 173–179. https://doi.org/10.1152/ajplegacy.1908.21.2.173
2. Hasselbalch, K. A. (1916). "Die Berechnung der Wasserstoffzahl des Blutes aus der freien und gebundenen Kohlensäure desselben, und die Sauerstoffbindung des Blutes als Funktion der Wasserstoffzahl." Biochemische Zeitschrift, 78, 112–144. https://doi.org/10.1152/physrev.1924.4.3.347
3. Van Slyke, D. D. (1922). "On the measurement of buffer values and on the relationship of buffer value to the dissociation constant of the buffer and the concentration and reaction of the buffer solution." Journal of Biological Chemistry, 52(2), 525–570. https://doi.org/10.1016/S0021-9258(18)85842-285842-2)
4. Po, H. N., & Senozan, N. M. (2001). "The Henderson–Hasselbalch equation: Its history and limitations." Journal of Chemical Education, 78(11), 1499–1503. https://doi.org/10.1021/ed078p1499