Gran plot

A Gran plot (also known as Gran titration or the Gran method) 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 K_a values (acid dissociation constants) of weak acids and bases from titration data.

Gran plots use linear approximations of the a priori non-linear relationships between the measured quantity, pH or electromotive potential (emf), and the titrant volume. Other types of concentration measures, such as spectrophotometric absorbances or NMR chemical shifts, can in principle be similarly treated. These approximations are only valid near, but not at, the end point, and so the method differs from end point estimations by way of first- and second-derivative plots, which require data at the end point. Gran plots were originally devised for graphical determinations in pre-computer times, wherein an x-y plot on paper would be manually extrapolated to estimate the x-intercept. The graphing and visual estimation of the end point have been replaced by more accurate least-squares analyses since the advent of modern computers and enabling software packages, especially spreadsheet programs with built-in least-squares functionality.

Basis of the calculations

The Gran plot is based on the Nernst equation which can be written as

${\displaystyle E=E^{0}+s\log\{H^{+}\}}$

where E is a measured electrode potential, E⁰ is a standard electrode potential, s is the slope, ideally equal to RT/nF, and {H⁺} is the activity of the hydrogen ion. The expression rearranges to

${\displaystyle [H^{+}]=10^{\frac {E-E^{0}}{s}}\ or\ [H^{+}]=10^{-pH}}$

depending on whether the electrode is calibrated in millivolts or pH. For convenience the concentration, [H⁺], is used in place of activity. In a titration of strong acid with strong alkali, the analytical concentration of the hydrogen ion is obtained from the initial concentration of acid, C_i and the amount of alkali added during titration.

${\displaystyle [H^{+}]={\frac {C_{i}v_{i}-c_{OH}v}{v_{i}+v}}}$

where v_i is the initial volume of solution, c_OH is the concentration of alkali in the burette and v is the titre volume. Equating the two expressions for [H⁺] and simplifying, the following expression is obtained

${\displaystyle C_{i}v_{i}-c_{OH}v=(v_{i}+v)10^{\frac {E-E^{0}}{s}}\ or\ =(v_{i}+v)10^{-pH}}$

A plot of ${\displaystyle (v_{i}+v)10^{\frac {E-E^{0}}{s}}\ or\ (v_{i}+v)10^{-pH}}$ against v will be a straight line. If E⁰ and s are known from electrode calibration, where the line crosses the x-axis indicates the volume at the equivalence point, ${\displaystyle C_{i}v_{i}=c_{OH}v}$. Alternatively, this plot can be used for electrode calibration by finding the values of E⁰ and s that give the best straight line.

Titrating strong acid with strong base

For a strong acid-strong base titration monitored by pH, we have at any i'th point in the titration

${\displaystyle K_{w^{}}=[H^{+}]_{i}[OH^{-}]_{i}}$

where K_w is the water autoprotolysis constant.

If titrating an acid of initial volume ${\displaystyle v_{0^{}}}$ and concentration ${\displaystyle [H^{+}]_{0^{}}}$ with base of concentration ${\displaystyle [OH^{-}]_{0^{}}}$, then at any i'th point in the titration with titrant volume ${\displaystyle v_{i^{}}}$,

${\displaystyle {\frac {v_{0}[H^{+}]_{0}-v_{i}[OH^{-}]_{0}}{v_{0}+v_{i}}}{\begin{cases}\approx [H^{+}]_{i}{\text{ or }}10^{-pH_{i}}&{\text{ when }}v_{0^{}}[H^{+}]_{0}>v_{i}[OH^{-}]_{0}{\text{ (acidic region)}}\\=0&{\text{ when }}v_{0^{}}[H^{+}]_{0}=v_{i}[OH^{-}]_{0}{\text{ (equivalence point)}}\\\approx -[OH^{-}]_{i}{\text{ or }}-K_{w}10^{pH_{i}}&{\text{ when }}v_{0^{}}[H^{+}]_{0}<v_{i}[OH^{-}]_{0}{\text{ (alkaline region)}}\end{cases}}}$

At the equivalence point, the equivalence volume ${\displaystyle v_{e^{}}=v_{i^{}}}$.

Thus,

• a plot of ${\displaystyle ({v_{0}+v_{i}})10^{-pH_{i}}{\text{ vs. }}v_{i^{}}}$ will have a linear region before equivalence, with slope ${\displaystyle -[OH^{-}]_{0^{}}}$
• and a plot of ${\displaystyle ({v_{0}+v_{i}})10^{pH_{i}}{\text{ vs. }}v_{i^{}}}$ will have a linear region after equivalence, with slope ${\displaystyle [OH^{-}]_{0}/K_{w^{}}}$
• both plots will have ${\displaystyle v_{e}=v_{0}[H^{+}]_{0}/[OH^{-}]_{0^{}}}$ as intercept

The equivalence volume is used to compute whichever of ${\displaystyle [H^{+}]_{0^{}}}$ or ${\displaystyle [OH^{-}]_{0^{}}}$ is unknown.

The pH meter is usually calibrated with buffer solutions at known pH values before starting the titration. The ionic strength can be kept constant by judicious choice of acid and base. For instance, HCl titrated with NaOH of approximately the same concentration will replace H⁺ with an ion (Na⁺) of the same charge at the same concentration, to keep the ionic strength fairly constant. Otherwise, a relatively high concentration of background electrolyte can be used, or the activity quotient can be computed.[1]

Titrating strong base with strong acid

Mirror-image plots are obtained if titrating the base with the acid, and the signs of the slopes are reversed.

${\displaystyle {\frac {v_{0}[OH^{-}]_{0}-v_{i}[H^{+}]_{0}}{v_{0}+v_{i}}}{\begin{cases}\approx [OH^{-}]_{i}{\text{ or }}K_{w}10^{pH_{i}}&{\text{ when }}v_{0^{}}[OH^{-}]_{0}>v_{i}[H^{+}]_{0}{\text{ (alkaline region)}}\\=0&{\text{ when }}v_{0^{}}[OH^{-}]_{0}=v_{i}[H^{+}]_{0}{\text{ (equivalence point)}}\\\approx -[H^{+}]_{i}{\text{ or }}-10^{-pH_{i}}&{\text{ when }}v_{0^{}}[OH^{-}]_{0}<v_{i}[H^{+}]_{0}{\text{ (acidic region)}}\end{cases}}}$

Click on the image to view in full size.
Figure 1. Sample Gran plots using data from Butler (1998). Note that the filled circles indicate the data points included in the least-squares computations to give the fitted dashed lines.

Click on the image to view in full size.
Figure 2. Sample Gran plots using data from "an on-line source". Retrieved 2008-02-18. Only the region near equivalence is shown, as data far from equivalence deviate strongly from linearity. Note that the filled circles indicate the data points included in the least-squares computations to give the fitted dashed lines.

Hence,

• a plot of ${\displaystyle ({v_{0}+v_{i}})10^{pH_{i}}{\text{ vs. }}v_{i^{}}}$ will have a linear region before equivalence with slope ${\displaystyle -[H^{+}]_{0}/K_{w^{}}}$
• and a plot of ${\displaystyle ({v_{0}+v_{i}})10^{-pH_{i}}{\text{ vs. }}v_{i^{}}}$ will have a linear region after equivalence with slope ${\displaystyle [H^{+}]_{0^{}}}$
• both plots will have ${\displaystyle v_{e}=v_{0}[OH^{-}]_{0^{}}/[H^{+}]_{0}}$ as x-intercept

Figure 1 gives sample Gran plots of a strong base-strong acid titration.

Concentrations and Dissociation Constants of Weak Acids

The method can be used to estimate the dissociation constants of weak acids, as well as their concentrations (Gran, 1952). With an acid represented by HA, where

${\displaystyle K_{a}={\frac {[H^{+}]_{i}[A^{-}]_{i}}{[HA]_{i}}}}$ ,

we have at any i'th point in the titration of a volume ${\displaystyle v_{0}}$ of acid at a concentration ${\displaystyle [HA]_{0}}$ by base of concentration ${\displaystyle [OH^{-}]_{0}}$. In the linear regions away from equivalence,

${\displaystyle [HA]_{i}\approx {\frac {v_{0}[HA]_{0}-v_{i}[OH^{-}]_{0}}{v_{0}+v_{i}}}}$  and

${\displaystyle [A^{-}]_{i}\approx {\frac {v_{i}[OH^{-}]_{0}}{v_{0}+v_{i}}}}$ 

are valid approximations, whence

${\displaystyle K_{a}\approx {\frac {10^{-pH_{i}}v_{i}[OH^{-}]_{0}}{v_{0}[HA]_{0}-v_{i}[OH^{-}]_{0}}}}$ , or

${\displaystyle K_{a}(v_{0}{\frac {[HA]_{0}}{[OH^{-}]_{0}}}-v_{i})\approx 10^{-pH_{i}}v_{i}}$ or, because ${\displaystyle v_{e}=v_{0}{\frac {[HA]_{0}}{[OH^{-}]_{0}}}}$ ,

${\displaystyle K_{a}(v_{e}-v_{i})\approx 10^{-pH_{i}}v_{i}}$.

A plot of ${\displaystyle 10^{-pH_{i}}v_{i}}$ versus ${\displaystyle v_{i^{}}}$ will have a slope ${\displaystyle -K_{a^{}}}$ over the linear acidic region and an extrapolated x-intercept ${\displaystyle v_{e^{}}}$, from which either ${\displaystyle [HA]_{0^{}}}$ or ${\displaystyle [OH^{-}]_{0^{}}}$ can be computed.[1] The alkaline region is treated in the same manner as for a titration of strong acid. Figure 2 gives an example; in this example, the two x-intercepts differ by about 0.2 mL but this is a small discrepancy, given the large equivalence volume (0.5% error).

Similar equations can be written for the titration of a weak base by strong acid (Gran, 1952; Harris, 1998).

Carbonate content

Martell and Motekaitis (1992) use the most linear regions and exploit the difference in equivalence volumes between acid-side and base-side plots during an acid-base titration to estimate the adventitious CO₂ content in the base solution. This is illustrated in the sample Gran plots of Figure 1. In that situation, the extra acid used to neutralize the carbonate, by double protonation, in volume ${\displaystyle v_{0^{}}}$ of titrate is ${\displaystyle (v_{e}-v_{e}^{\prime })[H^{+}]_{0^{}}=2v_{0}[CO_{2}]_{0}}$. In the opposite case of a titration of acid by base, the carbonate content is similarly computed from ${\displaystyle (v_{e}^{\prime }-v_{e})[OH^{-}]_{0^{}}=2v_{e}^{\prime }[CO_{2}]_{0}}$, where ${\displaystyle v_{e}^{\prime }}$ is the base-side equivalence volume (from Martell and Motekaitis).

When the total CO₂ content is significant, as in natural waters and alkaline effluents, two or three inflections can be seen in the pH-volume curves owing to buffering by higher concentrations of bicarbonate and carbonate. As discussed by Stumm and Morgan (1981), the analysis of such waters can use up to six Gran plots from a single titration to estimate the multiple end points and measure the total alkalinity and the carbonate and/or bicarbonate contents.

Potentiometric monitoring of H⁺

To use potentiometric (e.m.f.) measurements ${\displaystyle E_{i^{}}}$ in monitoring the ${\displaystyle H^{+_{}}}$ concentration in place of ${\displaystyle pH_{i^{}}}$ readings, one can trivially set ${\displaystyle -log_{10}[H^{+}]_{i}=b_{0}-b_{1}E_{i^{}}}$ and apply the same equations as above, where ${\displaystyle b_{0^{}}}$ is the offset correction ${\displaystyle nFE_{0^{}}/RT}$, and ${\displaystyle b_{1^{}}}$ is a slope correction ${\displaystyle nF^{_{}}/RT}$ (1/59.2 pH units/mV at 25°C), such that ${\displaystyle -b_{1}E_{i^{}}}$ replaces ${\displaystyle pH_{i^{}}}$.

Thus, as before for a titration of strong acid by strong base,

• a plot of ${\displaystyle ({v_{0}+v_{i}})10^{b_{1}E_{i}}}$ vs. ${\displaystyle v_{i^{}}}$ will have a linear region before equivalence, with slope ${\displaystyle -10^{b_{0}}[OH^{-}]_{0^{}}}$
• and a plot of ${\displaystyle ({v_{0}+v_{i}})10^{-b_{1}E_{i}}}$ vs. ${\displaystyle v_{i^{}}}$ will have a linear region after equivalence, with slope ${\displaystyle 10^{-b_{0}}[OH^{-}]_{0}/K_{w^{}}}$
• both plots will have ${\displaystyle v_{e}=v_{0}[H^{+}]_{0}/[OH^{-}]_{0^{}}}$ as intercept and, as before, the acid-side equivalence volume can be used to standardize whichever concentration is unknown, and the difference between acid-side and base-side equivalence volumes can be used to estimate the carbonate content

Analogous plots can be drawn using data from a titration of base by acid.

Electrode calibration

Note that the above analysis requires prior knowledge of ${\displaystyle b_{0^{}}}$ and ${\displaystyle b_{1^{}}}$.

If a pH electrode is not well calibrated, an offset correction can be computed in situ from the acid-side Gran slope:

• For a titration of acid by base, the acid-side slope (${\displaystyle -[OH^{-}]_{0^{}}10^{b_{0}}}$) can serve to compute ${\displaystyle b_{0^{}}}$ using a known value of ${\displaystyle [OH^{-}]_{0^{}}}$ or using the value given by the equivalence volume. ${\displaystyle K_{w^{}}}$ can then be computed from the base-side slope.
• For a titration of base by acid, as illustrated in the sample plots, the acid-side slope (${\displaystyle [H^{+}]_{0^{}}10^{b_{0}}}$) is similarly used to compute ${\displaystyle b_{0^{}}}$ and the base-side slope (${\displaystyle -[H^{+}]_{0^{}}10^{-b_{0}}/K_{w}}$) is used to compute ${\displaystyle K_{w^{}}}$ using a known value of ${\displaystyle [H^{+}]_{0^{}}}$ or using the value given by the acid-side equivalence volume.

In the sample data illustrated in Figure 1, this offset correction was not insignificant, at -0.054 pH units.

The value of ${\displaystyle b_{1^{}}}$, however, may deviate from its theoretical value and can only be assessed by a proper calibration of the electrode. Calibration of an electrode is often performed using buffers of known pH, or by performing a titration of strong acid with strong base. In that case, a constant ionic strength can be maintained, and ${\displaystyle [H^{+_{}}]_{i}}$ is known at all titration points if both ${\displaystyle [H_{}^{+}]_{0}}$ and ${\displaystyle [OH_{}^{-}]_{0}}$ are known (and should be directly related to primary standards). For instance, Martell and Motekaitis (1992) calculated the pH value expected at the start of the titration, having earlier titrated the acid and base solutions against primary standards, then adjusted the pH electrode reading accordingly, but this does not afford a slope correction if one is needed.

Based on earlier work by McBryde (1969), Gans and O'Sullivan (2000) describe an iterative approach to arrive at both ${\displaystyle b_{0^{}}}$ and ${\displaystyle b_{1^{}}}$ values in the relation ${\displaystyle -log_{10}[H^{+}]_{i}=b_{0}-b_{1}E_{i^{}}}$, from a titration of strong acid by strong base:

1. ${\displaystyle b_{0^{}}}$ is first estimated from the acidic data according to Rossotti and Rossotti (1965), and ${\displaystyle b_{1^{}}}$ is initially taken to have its theoretical value;
2. modified Gran function plots are drawn, using ${\displaystyle ({v_{0}+v_{i}})10^{b_{1}E_{i}-b_{0}}}$ vs. ${\displaystyle v_{i^{}}}$ on the acidic side of equivalence and ${\displaystyle ({v_{0}+v_{i}})K_{w}10^{-b_{1}E_{i}+b_{0}}}$ vs. ${\displaystyle v_{i^{}}}$ on the alkaline side, and the equivalence volumes ${\displaystyle v_{e^{}}}$ and ${\displaystyle v_{e}^{\prime }}$ are computed therefrom, as before;
3. as before, the difference in the equivalence volumes is used to compute the carbonate content but also to calculate an 'effective base concentration' ${\displaystyle v_{e}[OH^{-}]_{0^{}}/v_{e}^{\prime }}$ for the alkaline side of equivalence;
4. approximate ${\displaystyle [H^{+}]_{i^{}}}$ values are computed on the acid side as ${\displaystyle (v_{0^{}}[H^{+}]_{0}-v_{i}[OH^{-}]_{0})/(v_{0}+v_{i})}$ and on the alkaline side as ${\displaystyle (v_{e}[H^{+}]_{0}-v_{i}v_{e}[OH^{-}]_{0}/v_{e}^{\prime })/(v_{0}+v_{i})}$;
5. the initial definition ${\displaystyle -log_{10}[H^{+}]_{i}=b_{0}-b_{1}E_{i^{}}}$ is rewritten as ${\displaystyle E_{i^{}}=b_{0}/b_{1}+(1/b_{1})log_{10}[H^{+}]_{i}}$, and the ${\displaystyle E_{i^{}}}$ data are plotted against ${\displaystyle log_{10}[H^{+}]_{i^{}}}$, using those ${\displaystyle [H^{+}]_{i^{}}}$ values corresponding to pH values in the 2.5-4.5 and 10.7-11.5 ranges (the linear response range for a glass electrode that avoids variation of junction potentials and/or alkaline error at the pH extrema, and that additionally avoids measurement uncertainties near the equivalence point as well as computational errors from the neglect of ${\displaystyle [OH^{-}]_{i^{}}}$ on the acid side and the neglect of the carbonate/bicarbonate equilibrium on the alkaline side); a linear least-squares treatment provides ${\displaystyle 1/b_{1^{}}}$ as slope and ${\displaystyle b_{0}/b_{1^{}}}$ as intercept;
6. steps 2 and 3 are repeated with the new ${\displaystyle b_{0^{}}}$ and ${\displaystyle b_{1^{}}}$ values for greater precision in the equivalence volumes and the CO₂ content.

The procedure could in principle be modified for titrations of base by acid. A computer program named GLEE (for GLass Electrode Evaluation) implements this approach on titrations of acid by base for electrode calibration. This program additionally can compute (by a separate, non-linear least-squares process) a 'correction' for the base concentration. An advantage of this method of electrode calibration is that it can be performed in the same medium of constant ionic strength which may later be used for the determination of equilibrium constants.
Note that the regular Gran functions will provide the required equivalence volumes and, as ${\displaystyle b_{1^{}}}$ is initially set at its theoretical value, the initial estimate for ${\displaystyle b_{0^{}}}$ in step 1 can be had from the slope of the regular acid-side Gran function as detailed earlier. Note too that this procedure computes the CO₂ content and can indeed be combined with a complete standardization of the base, using the definition of ${\displaystyle v_{e^{}}}$ to compute ${\displaystyle [OH^{-}]_{0^{}}}$. Finally, the usable pH range could be extended by solving the quadratic ${\displaystyle (v_{0{^{}}}[H^{+}]_{0}-v_{i}[OH^{-}]_{0})/(v_{0}+v_{i})=[H^{+}]_{i}-K_{w}/[H^{+}]_{i}}$ for ${\displaystyle [H^{+}]_{i^{}}}$.

Potentiometric monitoring of other species

Click on the image to view in full size.
Figure 3. Sample Gran plots using data from a titration of Cl⁻ by Ag⁺ monitored potentiometrically. The potentials were converted to [Ag⁺] values for plotting. Note that the filled circles indicate the data points included in the least-squares computations to give the fitted dashed lines.

Potentiometric data are also used to monitor species other than ${\displaystyle H^{+_{}}}$. When monitoring any species ${\displaystyle S^{1_{}}}$ by potentiometry, one can apply the same formalism with ${\displaystyle -log_{10}[S^{1}]_{i}=b_{0}-b_{1}E_{i^{}}}$. Thus, a titration of a solution of another species ${\displaystyle S^{2_{}}}$ by species ${\displaystyle S^{1_{}}}$ is analogous to a pH-monitored titration of base by acid, whence either ${\displaystyle ({v_{0}+v_{i}})10^{b_{1}E_{i}}}$ or ${\displaystyle ({v_{0}+v_{i}})10^{-b_{1}E_{i}}}$ plotted versus ${\displaystyle v_{i^{}}}$ will have an x-intercept ${\displaystyle v_{0}[S^{0}]_{0}/[S^{1}]_{0^{}}}$. In the opposite titration of ${\displaystyle S^{1_{}}}$ by ${\displaystyle S^{2_{}}}$, the equivalence volume will be ${\displaystyle v_{0}[S^{1}]_{0}/[S^{0}]_{0^{}}}$. The significance of the slopes will depend on the interactions between the two species, whether associating in solution or precipitating together (Gran, 1952). Usually, the only result of interest is the equivalence point. However, the before-equivalence slope could in principle be used to assess the solubility product ${\displaystyle K_{sp^{}}}$ in the same way as ${\displaystyle K_{w^{}}}$ can be determined from acid-base titrations, although other ion-pair association interactions may be occurring as well.[2]

To illustrate, consider a titration of Cl⁻ by Ag⁺ monitored potentiometrically: ${\displaystyle {\frac {v_{0}[Cl^{-}]_{0}-v_{i}[Ag^{+}]_{0}}{v_{0}+v_{i}}}{\begin{cases}\approx [Cl^{-}]_{i}{\text{ or }}K_{sp}10^{-b_{1}E_{i}+b_{0}}&{\text{ when }}v_{0^{}}[Cl^{-}]_{0}>v_{i}[Ag^{+}]_{0}{\text{ (before equivalence)}}\\=0&{\text{ when }}v_{0^{}}[Cl^{-}]_{0}=v_{i}[Ag^{+}]_{0}{\text{ (equivalence point)}}\\\approx -[Ag^{+}]_{i}{\text{ or }}-10^{b_{1}E_{i}-b_{0}}&{\text{ when }}v_{0^{}}[Cl^{-}]_{0}<v_{i}[Ag^{+}]_{0}{\text{ (after equivalence)}}\end{cases}}}$

Hence,

• a plot of ${\displaystyle ({v_{0}+v_{i}})10^{-b_{1}E_{i}}{\text{ vs. }}v_{i^{}}}$ will have a linear region before equivalence, with slope ${\displaystyle -[Ag^{+}]_{0}10^{-b_{0}}/K_{sp^{}}}$
• and a plot of ${\displaystyle ({v_{0}+v_{i}})10^{b_{1}E_{i}}{\text{ vs. }}v_{i^{}}}$ will have a linear region after equivalence, with slope ${\displaystyle [Ag^{+}]_{0^{}}10^{b_{0}}}$
• in both plots, the x-intercept is ${\displaystyle v_{e}=v_{0}[Cl^{-}]_{0^{}}/[Ag^{+}]_{0}}$

Figure 3 gives sample plots of potentiometric titration data.

Non-ideal behaviour

In any titration lacking buffering components, both before-equivalence and beyond-equivalence plots should ideally cross the x axis at the same point. Non-ideal behaviour can result from measurement errors (e.g. a poorly calibrated electrode, an insufficient equilibration time before recording the electrode reading, drifts in ionic strength), sampling errors (e.g. low data densities in the linear regions) or an incomplete chemical model (e.g. the presence of titratable impurities such as carbonate in the base, or incomplete precipitation in potentiometric titrations of dilute solutions, for which Gran et al. (1981) propose alternate approaches). Buffle et al. (1972) discuss a number of error sources.

Because the ${\displaystyle 10^{pH_{i}}}$ or ${\displaystyle 10^{-pH_{i}}}$ terms in the Gran functions only asymptotically tend toward, and never reach, the x axis, curvature approaching the equivalence point is to be expected in all cases. However, there is disagreement among practitioners as to which data to plot, whether using only data on one side of equivalence or on both sides, and whether to select data nearest equivalence or in the most linear portions:[3][4] using the data nearest the equivalence point will enable the two x-intercepts to be more coincident with each other and to better coincide with estimates from derivative plots, while using acid-side data in an acid-base titration presumably minimizes interference from titratable (buffering) impurities, such as bicarbonate/carbonate in the base (see Carbonate content), and the effect of a drifting ionic strength. In the sample plots displayed in the Figures, the most linear regions (the data represented by filled circles) were selected for the least-squares computations of slopes and intercepts. Data selection is always subjective.

References

• Buffle, J., Parthasarathy, N. and Monnier, D. (1972): Errors in the Gran addition method. Part I. Theoretical Calculation of Statistical Errors; Anal. Chim. Acta 59, 427-438; Buffle, J. (1972): Anal. Chim. Acta 59, 439.
• Butler, J. N. (1991): Carbon Dioxide Equilibria and Their Applications; CRC Press: Boca Raton, FL.
• Butler, J. N. (1998): Ionic Equilibrium: Solubility and pH Calculations; Wiley-Interscience. Chap. 3.
• Gans, P. and O'Sullivan, B. (2000): GLEE, a new computer program for glass electrode calibration; Talanta, 51, 33–37.
• Gran, G. (1950): Determination of the equivalence point in potentiometric titrations, Acta Chemica Scandinavica, 4, 559-577.
• Gran, G. (1952): Determination of the equivalence point in potentiometric titrations—Part II, Analyst, 77, 661-671.
• Gran, G., Johansson, A. and Johansson, S. (1981): Automatic Titration by Stepwise Addition of Equal Volumes of Titrant Part VII. Potentiometric Precipitation Titrations, Analyst, 106, 1109-1118.
• Harris, D. C.: Quantitative Chemical Analysis, 5th Ed.; W.H. Freeman & Co., New. York, NY, 1998.
• Martell, A. E. and Motekaitis, R. J.: The determination and use of stability constants, Wiley-VCH, 1992.
• McBryde, W. A. E. (1969): Analyst, 94, 337.
• Rossotti, F. J. C. and Rossotti, H. (1965): J. Chem. Ed., 42, 375
• Skoog, D. A., West, D. M., Holler, F. J. and Crouch, S. R. (2003): Fundamentals of Analytical Chemistry: An Introduction, 8th Ed., Brooks and Cole, Chap. 37.
• Stumm, W. and Morgan, J. J. (1981): Aquatic chemistry, 2nd Ed.; John Wiley & Sons, New York.

Notes

1. If the ionic strength is not constant, one can correct this expression for the changing activity quotient. See Harris (1998) or this "on-line document" (PDF). Retrieved 2008-02-17.
2. Gran et al. (1981) give a more detailed treatment that takes into account other complex species in a titration of Cl⁻ by Ag⁺ (Ag₂Cl⁺ and AgCl₂⁻, notably) and in other precipitation titrations, in order to compute equivalence volumes of dilute solutions, when precipitation is incomplete.
3. "On-line document" (PDF). Archived from the original (PDF) on December 17, 2008. Retrieved 2008-02-17.) and W. Wolbach of Depaul University ("On-line document" (PDF). Retrieved 2008-02-17.) recommend using last 10-20% volume data before the equivalence point, and the latter recognizes that base-side Gran plots from titrations of acid by base (i.e after the equivalence point) can be used to assess the CO₂ content in the base. Similarly, W. E. Brewer and J. L. Ferry of the University of South Carolina recommend using those data within 10% before equivalence ("On-line document" (PDF). Retrieved 2008-02-17.). K. Kuwata of Macalester College recommends that students choose whichever data region gives the straightest line before equivalence ("On-line document" (PDF). Retrieved 2008-02-17.). D. L. Zellmer of the California State University at Fresno asks students to plot data from both sides of equivalence, using data furthest from equivalence, and to assess the errors in order to determine whether or not the two estimates of the equivalence volumes are significantly different (pH data: "On-line document". Retrieved 2008-02-17.; potentiometric titration of chloride ion with silver ion: "On-line document". Retrieved 2008-02-17.).
4. Butler (1991) discusses the issue of data selection, and also examines interferences from titratable impurities such as borate and phosphate.