index previous next

Triprotic Titration Model and Data Analysis

OpenOffice Calc Versions


Screen shot. Click to see larger graphic

A model of the titration curve of a weak triprotic acid titrated by a strong base. For a triprotic acid, that is a quintic equation - too complex to evaluate by hand, but easy for a spreadsheet to handle. By using the sliders. you can change the acid concentration, acid volume, base concentrations, and the three acid constants of the acid (K1a, K2a, and K3a). (You can even change the Kw of water in cell J2). Assumptions: Activity effects are ignored. The pH is measured correctly without error at the ends of the pH range.

View Equations (.txt)
Download spreadsheet in OpenOffice format (.ods)

Note: to run these spreadsheets, you have to first download the OpenOffice installer (download from OpenOffice), then install it (by double-clicking on the installer file that you just downloaded), and then download my spreadsheets from this page.  Once OpenOffice is installed, you can run my spreadsheets just by double-clicking on them.  Note 1: Don't use version 3.1. There is a bug in OpenOffice 3.1 that causes bad x-axis scaling on some of my graphs.  The problems does not occur in version 3.0 or in the most recent version 3.2.  Note 2: Downloading these files with Interent Explorer will change the file types from ".ods" to ".zip"; you will have to edit the file names and change the extensions back to ".ods" for them to work properly.  This problem does not occur in Firefox or in Chrome.

The data entry version ("TriproticTitrationData") is similar to above except that it has a space to type in some experimental titration data (volume and pH). After performing a pH titration of a weak triprotic acid, students type their pH/titrant volume data. By adjusting the parameters of the model and observing graphically the fit between the experimental data (circles) and the calculated model (line) they can estimate the unknown parameters, such as the pKs of the acid. The spreadsheet plots the experimental data on the same plot (as blue dots). The sliders can then be adjusted until the theoretical curve (red line) can be made to match the data (blue dots). The RMS error between the data and the model is also displayed.

../SimpleModels/TriproticTitrationData.GIF
Data entry version. Click to see larger graphic

Download data entry version with blank data table

For demonstration purposes, you can also download versions with sample data entered for the following titration: Acid concentration=0.1M, acid volume=10mL; titrant (NaOH) concentration=0.1M ; acid dissociation constants pK1=2.0; pK2=7.0;pK3=12.0. Titrant volume measured to the nearest 0.1 mL; pH measured to the nearest 0.1 pH unit.  You can try to measure the experimental acid concentration by adjusting the acid concentration slider, or the pKa's by adjusting the three pKa sliders, to minimize the RMS error of the fit.  
Download data entry version with sample data already entered

This version is like the above except that the adjustable variables have already been adjusted for the best possible fit.
Download data entry version with variables adjusted for best fit


WingZ Version


Click to see larger view
Screen shot. Click to see larger view
[OpenOffice and Excel Versions]

Download links: TriproticTitration.wkz;

Wingz player application and basic set of simulation modules, for windows PCs or Macintosh

Other related simulations:
Monoprotic Titration Curve model

[Return to Index]

Interpretation of the Titration Curve of a Triprotic Acid

The titration of a triprotic weak acid H3A with a strong base (e.g. NaOH) involves six solution species: H3O+, OH, H3A, H2A-1, HA-2, and the base cation, e.g. Na+. Other variables are the initial acid concentration (Ainit), initial base concentration (Binit), volume of acid titrated (Va), and volume of base added (Vb).

The algebraic description of a triprotic titration (neglecting activities, as usual) is completely specified by three equations for the three ionization steps (K1, K2, and K3), the equation for water ionization (Kw), equations for the mass and charge balance, and in addition two equations to account for dilution of the acid and base concentrations during the titration.

Using the usual "computer algebra" notation, with H representing the hydronium ion and CA representing the total concentration of acid in all forms:

H*H2A/H3A==K1	expression for K1
H*HA/H2A==K2 expression for K2
H*A/HA==K3 expression for K3
H*OH==Kw water
CA==H3A+H2A+HA+A mass balance for total concentration of acid
H+B==H2A+2*HA+3*A+OH charge balance (B = the base cation)
B==Vb Binit/(Va+Vb) concentration of base during titration
CA==Va Ainit/(Va+Vb) concentration of acid during titration
Eliminating the variables H2A, HA, and OH, B, and CA between these equations and solving for Vb using a computer algebra program (e.g. Mathematica or Maple) yields:
Vb = -((Va*(H^5 - Ainit*H^3*K1 + H^4*K1 - 2*Ainit*H^2*K1*K2 + H^3*K1*K2 - 
3*Ainit*H*K1*K2*K3 + H^2*K1*K2*K3 - H^3*Kw - H^2*K1*Kw - H*K1*K2*Kw -
K1*K2*K3*Kw)) / ((Binit*H + H^2 - Kw)*(H^3 + H^2*K1 + H*K1*K2 + K1*K2*K3)))
This expression gives the volume of base Vb as a function of hydrogen ion concentration H, the three Ks, Kw, the volume of acid Va, and the initial concentrations of acid and base Ainit and Binit. It is completely general and works for all concentrations and for Ks. The titration curve is obtained by plotting Vb on the x-axis and pH (=-log(H)) on the y-axis. This is referred to as an inverse solution, because we usually think of Vb as the independent variable and H as the dependent variable. In fact, it is in principle possible to solve this expression directly for H as a function of Vb, but the solution is extremely complex. Fortunately, this is not necessary for our purposes because we simply want to know how closely the theoretical expression above describes the experimental data. Our titration data will consist of pairs of experimentally measured volumes and pHs, and we could use either an expression for H as a function of Vb (to predict the volume Vb at each measured H) or an expression for Vb as a function of H (to predict the volume Vb at each measured H). We will do the latter.

To facilitate evaluating the above expression, the spreadsheets on this page already contain this equation, so you won't have to type it in. Follow the instructions in the scrolling text field at the bottom of the screen. Adjust the variable parameters to get the best possible fit between your experimental data and the theoretical curve (line). Get a print-out of the screen and submit with your lab report.



(c) 1991, 2015. This page is part of Interactive Computer Models for Analytical Chemistry Instruction, created and maintained by Prof. Tom O'Haver , Professor Emeritus, The University of Maryland at College Park. Comments, suggestions and questions should be directed to Prof. O'Haver at toh@umd.edu. Revised June 2009.  Number of unique visits since May 17, 2008: