Abstract
Upon the addition of protons to an aqueous solution containing multiple buffers, the final H+ concentration ([H+]) at equilibrium is determined by the partitioning of added H+ among the various buffer components. In the analysis of acid-base chemistry, the Henderson-Hasselbalch equation and the Stewart strong ion formulation can only describe (rather than predict) the equilibrium pH following a proton load since these formulas calculate the equilibrium pH only when the reactant concentrations at equilibrium(1) 1The term "equilibrium" refers to the steady state proton and reactant concentrations when the buffering of excess protons by the various buffers is complete. are already known. In this regard, it is simpler to directly measure the equilibrium pH rather than measure the equilibrium reactant concentrations to calculate the equilibrium pH. As these formulas cannot predict the final equilibrium [H+] following a proton load to a multiple-buffered aqueous solution, we developed a new quantitative approach for predicting the equilibrium [H+] that is based on the preequilibrium(2)2 The term "preequilibrium" refers to the initial proton and reactant concentrations immediately upon addition of protons and before the buffering of excess protons by the various buffers. concentrations of all buffers in an aqueous solution. The mathematical model used to derive our equation is based on proton transfer buffer equilibria without requiring the incorporation of electroneutrality considerations. The model consists of a quartic polynomial equation that is derived based solely on the partitioning of H+ among the various buffer components. We tested the accuracy of the model using aqueous solutions with various buffers and measured the equilibrium pH values following the addition of HCl. Our results confirmed the accuracy of our new equation (r2 = 1; measured pH vs. predicted pH), indicating that it quantitatively accounts for the underlying acid-base phenomenology.