Pufferkapazität und Alkalinität

Puffer und alkalimetrische Titration

Ein Puffer ist ein Gemisch aus einer schwachen Säure und einer Base. Im Folgenden betrachten wir eine 2-protonige Säure H2A, zu welcher eine 1-wertige starke Base BOH hinzugeben bzw. hinzutitriert wird (mit B+ = Na+, K+ oder NH4+):1

(1) H2A + n BOH  =  BnH2-n A + n H2O

Beispiel: Für H2CO3 und NaOH entspricht dies im Spezialfall n = 0, 1, und 2 jeweils einer reinen H2CO3, NaHCO3 und Na2CO3-Lösung (siehe Äquivalenzpunkte).

Die Größe n agiert hier als stöchiometrischer Koeffizient. Sie hängt mit der Menge an zugegebener Base CB wie folgt zusammen:

(2)   n = CB /CT     oder     CB = n CT

Puffersystem: Säure + Base im chemischen Gleichgewicht

wobei CT die Säuremenge bezeichnet. Beide Größen CB und CT besitzen die Maßeinheit mol/L, während n dimensionslos ist.

Die Größe BnH2-n A in Reaktionsgleichung (1) ist im Wasser instabil; sie dissoziert in mehrere aquatische Spezies – siehe Bild rechts.

Die Stoffmenge der starken Base entspricht zudem der Alkalinität, Alk = CB. So gilt:

(3)   n = Alk /CT     oder     Alk = n CT

Variiert man nun den Wert von n (oder CB) in der Reaktionsgleichung (1), so spricht man von der alkalimetrischen Titration.

Anmerkung: Der Grund, weshalb wir unser Augenmerk auf 2-protonige Säuren richten, ist, dass zu dieser Gruppe die Kohlensäure H2CO3 als deren prominentester Vertreter gehört (also A-2 = CO3-2).

Notation. Im Weiteren kürzen wir die H+-Konzentration mit x ab. Der Zusammenhang mit dem pH-Wert ist dann gegeben durch:2

(4)   x = [H+] = 10-pH pH = –log [H+] = –log x

Ebenso führen wir die Größe “Alkalinität von reinem Wasser” wie folgt ein:

(5)   w  =  [OH-] – [H+]  =  Kw /x – x

wobei Kw = 10-14 die Gleichgewichtskonstante der Selbstionisation des Wassers ist (Autoprotolyse). Wie wir bald sehen werden, repräsentiert w den (additiven) Anteil der Autoprotolyse zur Alkalinität des Puffersystems. Für reines Wasser – definiert durch x = (Kw)1/2 bzw. pH = 7 – ist dieser Beitrag exakt Null: w = 0.

Im Spezialfall von reinem Wasser – definiert durch x = (Kw)1/2 bzw. pH = 7 – ist diese “Alkalinität” exakt Null: w = 0.

Ausgangspunkt: Das Gleichungssystem

Die mathematischen Gleichungen zur Beschreibung der alkalimetrischen Titration sind exakt dieselben, wie wir sie für 2-protonige Säuren im Zusammenhang mit ihren konjugierten Basen diskutiert haben:3

(6a) K1 =  {H+} {HA-} / {H2A} (1. Dissoziation)
(6b) K2 =  {H+} {A-2} / {HA-} (2. Dissoziation)
(6c) Kw =  {H+} {OH-} (Selbstdissoziation)
(6d) CT =  [H2A] + [HA-] + [A-2] (Massenbalance)
(6e) 0 =  [H+] + n [H2A] + (n-1) [HA-] + (n-2) [A-2] – [OH-] (Protonenbalance)

Im Unterschied zur bisherigen Betrachtung stellt n nun eine kontinuierliche Variable dar (anstelle eines festen, diskreten Werts von n = 0, 1, 2).

Der Spezialfall n = 0 beschreibt das reine Säure-Systems, bei dem [H+] (oder pH) durch Vorgabe von CT exakt bestimmt ist, und umgekehrt. Nun, durch Hinzunahme der neuen, kontinuierlichen Variablen n bekommt unser System einen neuen Freiheitsgrad.

Das obere Gleichungssystem kann in ein Polynom 4. Ordnung in x (= 10-pH) umgestellt werden:

(7)   x4 + {K1 + nCT} x3 + {K1K2 + (n–1)CT K1 – Kw} x2 + K1 {(n–2)CT K2 – Kw} x – K1K2Kw  =  0

Diese quartische Gleichung (nicht zu verwechseln mit “quadratische Gleichung”) liefert den pH-Wert für jedes Wertepaar von CT und n.

Kennt man den pH (bzw. x), dann lassen sich die Konzentrationen der drei gelösten Spezies H2 A, HA- und A-2 mittels der Verteilungskoeffizienten ai berechnen:

(8)   [H2-i A-i]  =  CT ai (x) für   i = 0, 1, 2

Das Gleichungssystem (6) einerseits, und 7 zusammen mit 8 andererseits, besitzen beide den selben Informationsgehalt. Im weiteren Verlauf werden wir noch andere analytische Formeln vorstellen, die exakt das Gleiche tun (und uns das nicht-triviale Lösen einer quartischen Gleichung ersparen).

Das Gleichungssystem mit Einbeziehung der Alkalinität

Wir wollen nun die kontinuierliche, maßeinheitslose Variable n durch die Alkalinität ersetzen. Dazu wird 6d in 6e eingesetzt, und wir erhalten

(9a)   n CT  =  ([OH-] – [H+]) + [HA-] + 2 [A-2]

Die rechte Seite ist nichts anderes als die Definition der Alkalinität (in Carbonatsystemen auch als m-Alkalinität bekannt):

(9b) Alkalinität: Alk  ≡  ([OH-] – [H+]) + [HA-] + 2 [A-2]

Ein Blick auf die letzten beiden Gleichungen bestätigt vollends 3: Alk = n CT.

Gleichungssystem. Ersetzt man die Protonenbalance in 6e durch die Alkalinität in 9b, dann folgt:

(10a) K1 =  {H+} {HA-} / {H2A} (1. Dissoziation)
(10b) K2 =  {H+} {A-2} / {HA-} (2. Dissoziation)
(10c) Kw =  {H+} {OH-} (Selbstdissoziation)
(10d) CT =  [H2A] + [HA-] + [A-2] (Massenbalance)
(10e) Alk =  ([OH-] – [H+]) + [HA-] + 2 [A-2] (Alkalinität)

Damit ändert sich auch das Polynom in 7 zu:

(11)   x4 + {K1 + Alk} x3 + {K1K2 + (Alk – CT) K1 – Kw} x2 + K1 {(Alk – 2CT) K2 – Kw} x – K1K2Kw  =  0

Die Spezies-Konzentrationen bleiben von dieser Transformation unangetastet. Für sie gilt nach wie vor 8.

Verteilungskoeffizienten – Ionization Fractions (a0, a1, a2)

Aus rein mathematischer Sicht sind die Verteilungskoeffizienten eine clevere Art und Weise, 5a und 5b miteinander zu verknüpfen. Die drei Verteilungskoeffizienten lauten:

(12a)   a0   =   [ 1 + K1/x + K1K2/x2 ]-1 = (x/K1) a1
(12b)   a1   =   [ x/K1 + 1 + K2/x ]-1    
(12c)   a2   =   [ x2/(K1K2) + x/K2 + 1 ]-1 = (K2/x) a1

Das sind allein Funktionen des pH (bzw. x); sonst enthalten sie nur die beiden Gleichgewichtskonstanten K1 und K2. Ihre pH-Abhängigkeit ist z.B. hier dargestellt. Andere Eigenschaften sind:

Alle drei Koeffizienten ergeben zusammen exakt 1:

(13a) a0 + a1 + a2 = 1 (Massenbalance)

Die Werte aller ai sind weder negativ noch größer 1 (selbst 0 und 1 sind, streng genommen, ausgeschlossen):

(13b) 0 < ai < 1 für   i = 0, 1, 2

Es gibt drei besondere pH-Werte (bzw. Werte von x) für die gilt:4

(14a)   pH = pK1 x = K1 a1 = a0 ≈ ½ a2 ≈ 0
(14b)   pH = ½ (pK1+pK2) x = (K1K2)1/2 a0 = a2 ≈ 0 a1 ≈ 1
(14c)   pH = pK2 x = K2 a1 = a2 ≈ ½ a0 ≈ 0

Asymptotische Werte. An den beiden entgegengesetzten Enden der pH-Skala werden folgende Grenzwerte erreicht:

(14d) stark sauer: pH → 0 (oder x → ∞): a0 = 1 a1 = 0 a2 = 0
(14e) stark basisch: pH → 14 (oder x → 0): a0 = 0 a1 = 0 a2 = 1

Die drei Verteilungskoeffizienten (a0, a1 und a2) sind Bausteine, aus denen wir alle nachfolgenden Beziehungen und Kennwerte aufbauen werden. Hier sei nur daran erinnert, dass im Gegensatz zur kontinuierlichen Variablen n, der Index i = 0, 1, 2 nur für ganzzahlige Werte reserviert ist. Das gilt sowohl für die Koeffizienten ai als auch für die drei aquatischen Spezies H2A, HA-, und A-2.

Exakte Beziehungen zwischen pH, CT und n (bzw. Alk = nCT)

Das Gleichungssystem (6a) bis (6e) – oder alternativ (10a) bis (10e) – beschreiben das System “2-protonige Säure + starke Base” (alkalimetrische Titration). Der Gleichgewichtszustand (d.h. die Spezies-Konzentrationen von H2A, HA- und A-2) ist damit eindeutig festgelegt, wenn man zwei von den drei Parametern CT, n und pH vorgibt (bzw. alternativ von CT, Alk und pH). Der dritte Parameter aus dieser Dreiergruppe ist dann auch festgelegt:

(16a)   pH (CT,n) = –log x1  mit  x1 = positive Lösung von 7   [i] = CT ai
(16b)   pH (CT,Alk) = –log x1  mit  x1 = positive Lösung von 11   [i] = CT ai
(17a)   n (CT,pH) = a1 + 2a2 + w/CT [i] = CT ai
(17b)   Alk (CT,pH) = CT (a1 + 2a2) + w [i] = CT ai
         
(18a)   CT (n,pH) = w/(n – a1 – 2a2) [i] = [w/(n – a1 – 2a2)] ai
(18b)   CT (Alk,pH) = (Alk – w)/(a1 + 2a2) [i] = [(Alk – w)/(a1 + 2a2)] ai

wobei anstelle von pH meist x (= 10-pH) verwendet wird. Mit [i] ist hier die molare Konzentration der drei aquatischen Spezies abgekürzt:

(19)   [i]  ≡  [H2-i A-i] für   i = 0, 1, 2

Alle hier gelisteten Gleichungen sind unterschiedliche Ausdrucksformen ein und desselben Sachverhalts, nämlich des Gleichungssystems (6a) bis (6e).

Diagramme. Die Nichtlinearität dieser Gleichungen veranschaulicht man sich am besten graphisch. Gegeben sei das reine Carbonat-Puffer-System. Die nachfolgenden Diagramme zeigen alle möglichen Kombinationen jeweils einer abhängigen Variable von zwei unabhängigen Parametern:

abhängige und unabhängige Variablen (CT, n, pH) des Carbonat-Puffersystems

Etwas andere Diagramme erhält man, wenn der “stöchiometrische Koeffizient” n durch die Alkalinität ersetzt wird (Alk = n CT):

abhängige und unabhängige Variablen (CT, Alkalinität, pH) des Carbonat-Puffersystems

Weitere Beispiele. Die Anwendung von 17a auf ein Carbonatsystem liefert Titrationskurven, wie sie auch im Diagramm unten als blaue Kurven eingetragen sind. Ein vollkommen anderes Bild ergibt sich, wenn man 18a für den Spezialfall n = 0, 1, 2 darstellt – siehe Äquivalenzpunkte.

Pufferintensität

Durch Ableitung der Titrationskurve n(pH) nach pH erhält man die Pufferintensität β. Im folgenden unterscheiden wir zwei Arten der Pufferintensität:

(20a) normiert: β  = dn / d pH [–]
(20b) nicht normiert: βc = dCB / d pH = β CT [mol/L]

mit CB = nCT. Beide Größen unterscheiden sich nur in der Maßeinheit, wobei die normierte Pufferintensität maßeinheitslos ist.

pH-gepuffertes und ungepuffertes System

Je steiler der Anstieg der Titrationskurve, desto größer ist die Pufferintensität β = dn/d(pH), d.h. desto größer ist der Widerstand des Systems gegenüber pH-Änderungen. Dort wo β sein Maximum hat, befindet sich auch der optimale Pufferbereich (pH±1 um den Maximalwert) – wie im unteren Diagramm:

optimaler Pufferbereich Carbonatsystem

Gleichungen. Den Ausgangspunkt der Berechnungen bildet die mit 17a beschriebene Titrationskurve:

(21)   n (CT,x)  = 

durch Ableitung nach dem pH erhält man die Pufferintensität (siehe Anhang):

(22a)   β =   ln 10 {(a2 + a0) – (a2 – a0)2} + ln 10 (w+2x) /CT

oder, in ausgeschriebener Form,

(22b)  

Aus 22b mit ihren durchweg positiven Summanden wird ersichtlich, dass β immer eine positive Größe ist (sie wird niemals negativ oder Null).

Extremalpunkte. Die lokalen Maxima (und Minima) der Pufferintensität befinden sich an folgenden Stellen:

x = K1 Maximum optimaler Pufferbereich
x = (K1K2)1/2 Minimum    
x = K2 Maximum optimaler Pufferbereich (nicht für kleine CT)

Man beachte, dies gilt nicht exakt; es ist eine Näherung, aber eine sehr gute. Die Extremwerte befinden sich an der Stelle, wo die erste Ableitung von β ihre Nullstellen hat.

Beispiel: Die Pufferintensität des einfachen CO2-System ist als grüne Kurve im Diagramm unten dargestellt.

Pufferkapazität. Die Pufferkapazität erhält man durch Integration der Pufferintensität:

(23) Pufferkapazität =

Ableitung der Pufferintensität nach pH

Mann kann nun einen Schritt weiter gehen, und die Pufferintensität β erneut nach dem pH-Wert ableiten. Das ergibt:

(24)  

mit der Abkürzung:

(25)   f(x) = (x/K1 – K2/x) (x/K1 + K2/x – 1 + 8 K2/K1)

Beispiel: Die Ableitung der Pufferintensität für das einfache CO2-System ist als rote Kurve im Diagramm unten dargestellt. Die Nullstellen von dβ/d(pH) sind als blaue Kreise eingezeichnet.

Anwendung auf das Carbonatsystem

Gegeben sei ein einfaches CO2-System mit CT = 100 mM, 10 mM und 1 mM H2CO3.5 Basierend auf den oberen Gleichungen werden drei Größen dargestellt:

•  n (normierte) Zugabemenge an Base6 (blaue Kurve)
•  β = dn/d(pH) (normierte) Pufferintensität (grüne Kurve)
•  dβ/d(pH) Ableitung der Pufferintensität (rote Kurve)

Alle drei Größen sind per Definition maßeinheitslos. Die kleinen blauen Kreise markieren die Nullstellen von dβ/d(pH).

aqion - Pufferintensität des Carbonatsystems

Da die (blaue) Titrationskurve immer ansteigend ist, muss die Pufferintensität β positiv sein (siehe grüne Kurve). Das entspricht dem Prinzip von Le Châtelier: Jede Lösung widersetzt sich einer pH-Änderung.

Anhang – Ableitungen nach pH

Die erste und zweite Ableitung von x = [H+] nach dem pH sind:

(A1)   dx/d(pH)  =  (–ln 10) x und   d2x/d(pH)2  =  (–ln 10)2 x

Dieses Ergebnis hilft uns, jede beliebige Funktion f(x) nach dem pH abzuleiten (mittels Kettenregel):

(A2)  

So folgt für die “Alkalinität des reinen Wassers” in 5, w(x) = Kw/x – x, das interessante Resultat:

(A3a)   dw/d(pH) =   ln 10 (Kw/x + x) =   ln 10 (w + 2x)
(A3b)   d2w/d(pH)2 =   (ln 10)2 (w + 2x – 2x) =   (ln 10)2 w

Alle höheren Ableitungen wiederholen sich nach dem gleichen Muster:

(A3c)  

Ionization Fractions. Die in 12a bis 12c eingeführten Koeffizienten erfüllen die folgenden Beziehungen (wobei die Summation von i = 0 bis 2 läuft):

(A4a)   Y0  ≡  Σi ai =  a0 + a1 + a2 =  1  
(A4b)   Y1  ≡  Σi i ai =  0⋅a0 + 1⋅a1 + 2⋅a2 =  a1 + 2 a2  
(A4c)   Yk  ≡  Σi ik ai =  0⋅a0 + 1⋅a1 + 2k⋅a2 =  a1 + 2k a2 (für alle k > 0)

Größen dieser Art werden auch als k-tes Moment bezeichnet. In unserem Fall verbirgt sich hinter diesen folgende Bedeutung: Y0 repräsentiert die Massenbalance, Y1 ist der Hauptbestandteil der Alkalinität und Y2 zusammen mit Y3 werden für die Pufferintensität und deren Ableitung gebraucht. Alle Yk (mit Ausnahme von Y0) sind positive Größen mit Wertebereich 0 < Yk < 2k.

Für die 1. Ableitung der Koeffizienten nach pH erhält man

(A5)   dai /d(pH)  =  ln 10 {i – (a1 + 2a2)} ai  =  ln 10 (i – Y1) ai für i = 0, 1, 2

Wendet man dies auf die Summen in A4a und A4b, so folgt

(A6a)   (d/d pH) Σi ai =   ln 10 Σi (i – Y1) ai =   ln 10 (Y1 – Y1)  =  0
(A6b)   (d/d pH) Σi i ai =   ln 10 Σi i (i – Y1) ai =   ln 10 (Y2 – Y12)

Die erste Beziehung ist offensichtlich, da es sich hier um die Ableitung einer Konstanten handelt, nämlich d1/d(pH) = 0. Im Allgemeinen gilt:

(A7)   d Yk /d pH   =   ln 10 Σi ik (i – Y1) ai   =   ln 10 (Yk+1 – Y1Yk)

Die 2. Ableitung von Y1 liefert

(A8a)   =  
      =   (ln 10)2 (Y3 – Y1Y2 – 2Y1 (Y2 – Y12))
      =   (ln 10)2 (Y3 – 3Y1Y2 + 2Y13)

und die 3. Ableitung von Y1 liefert

(A8b)   =  
        =   (ln 10)3 {(Y4–Y1Y3) – 3Y2 (Y2–Y12) – 3Y1 (Y3–Y1Y2) + 6Y12 (Y2–Y12)}

Pufferintensität & Co. Die oben hergeleiteten Hilfsgleichungen lassen sich nun auf die Pufferintensität und deren Ableitung – definiert in 22 und 24 – anwenden. Konkret gilt:

(A9a)   n = (Y1 + w/CT)
(A9b)   β =   d /d pH (Y1 + w/CT)
(A9c)   dβ/d pH =   d2/d (pH)2 (Y1 + w/CT)

Das ergibt:

(A10a)   n =   (Y1 + w/CT)
(A10b)   β =   (ln 10)  {(Y2 – Y12) + (w+2x) /CT}
(A10c)   dβ/d pH =   (ln 10)2 {(Y3 – 3Y1Y2 + 2Y13) + w/CT}

Man kann nun die Y’s durch die a’s ersetzen. Für den “Hauptbestandteil” von A10b folgt dann beispielsweise:

(A11)   Y2 – Y12   =   (a1 + 4a2) – (a1 + 2a2)2   =   (a2 + a0) – (a2 – a0)2

wobei in der letzten Beziehung A4a zur Anwendung kam.

Anmerkungen

  1. Zu diesem Thema gibt es sowohl eine PowerPoint-Präsentation als auch ein pdf

  2. Genau genommen basiert der pH-Wert nicht auf der Konzentration, sondern auf der Aktivität von H+: pH = –log {H+}. 

  3. Hier und im Folgendem wird angenommen, dass alle Aktivitäten (symbolisiert durch geschweifte Klammern) durch Stoffmengenkonzentrationen (symbolisiert durch eckige Klammern) ersetzt werden können. Dies ist statthaft für verdünnte Wässer und/oder beim Übergang zu konditionellen Gleichgewichtskonstanten. 

  4. Anstelle der Gleichgewichtskonstanten K verwendet man auch pK = –log K

  5. Die beiden Gleichgewichtskonstanten der Kohlensäure sind: K1 = 10-6.35 und K2 = 10-10.33 (siehe auch hier). 

  6. Auch negative Werte von n sind erlaubt. Damit wird der Entzug von Base simuliert, was rein formal das Gleiche ist wie die Zugabe einer starken Säure (z.B. HCl). 

[last modified: 2017-04-16]