Charge neutrality and mass action

So far we have worked in one direction: given the concentrations, we read off the species voltages. Real problems usually run the other way around. Some of the voltages are pinned for us, whether by a contact, a reservoir, or a reaction with a solid phase, and what we are left to determine is where the standard-state ladder ViV^\circ_i sits. Pinning the ViV_i does not pin the ladder, and pinning the ladder does not pin the ViV_i; the two are tied together only through a constraint that we have so far set aside: bulk charge neutrality.

Handled algebraically this quickly becomes tedious, since every concentration depends exponentially on the gap between a carrier's voltage and its standard state. On a ViV_i diagram, though, the whole business reduces to a single motion: sliding the rigid ViV^\circ_i ladder up or down until the net charge balances.

As a reminder, the ideal-dilute relation from the solutions topic can be rearranged to express a concentration in terms of that gap:

ci=cexp ⁣(ziF(ViVi)RT).c_i = c^\circ \exp\!\left(\frac{z_i F\,(V_i - V^\circ_i)}{RT}\right).

A carrier's concentration is fixed entirely by how far its ViV_i floats away from its own rung on the ladder, or as we will see it here: where the ViV^\circ_i moves relative to the fixed ViV_i.

Neutrality sets the ladder

Suppose we fix a single carrier's voltage and then slide the ladder. Because only the gap ViViV_i - V^\circ_i enters the formula, moving the ladder down widens the gap and enriches a pinned cation, while moving it up depletes it; for an anion the directions reverse. Every mobile carrier responds to the slide this way at once, exponentially, each with the sign of its own charge.

This is exactly the handle that bulk charge neutrality needs. Away from interfaces, a conductor carries essentially no net charge, so

iziFci+ρbg=0,\sum_i z_i F c_i + \rho_{\mathrm{bg}} = 0,

where ρbg\rho_{\mathrm{bg}} is any fixed, immobile background charge. As the ladder slides down, every cation in this sum enriches and every anion depletes, so the net charge climbs monotonically and passes through zero at exactly one position. Once the ViV_i are given, the ladder floats to precisely the offset that neutrality requires, and to no other.

VX+V_{X^{+}}VX+V^\circ_{X^{+}}VY2+V_{Y^{2+}}VY2+V^\circ_{Y^{2+}}VZV_{Z^{-}}VZV^\circ_{Z^{-}}Solution−0.15−0.10−0.050.000.050.100.150.20Voltage (V) — arbitrary offset
Net chargeladder offset00ρbg\rho_{\mathrm{bg}}ρbg+iziFci\textstyle \rho_{\mathrm{bg}} + \sum_i z_i F c_i

A multi-ion solution with three ViV_i fixed, and some fixed negative background charge. Top: species voltage levels and the ViV^\circ_i ladder. Bottom: the total charge density as a function of the ViV^\circ_i ladder offset, crossing zero at the single neutral position. Sliding the ladder is how the system "finds" neutrality; the tick on the slider marks the neutral offset.

Electrostatic mass action: a pair invariant

Neutrality moves the ladder around: dissolve something new, or change the background charge, and the ladder must slide to a new neutral offset, carrying every concentration with it. Certain combinations of concentration, however, come through untouched. Pin two carriers, a species AA (VAV_{\mathrm{A}}) and a species BB (VBV_{\mathrm{B}}), and let the ladder go wherever neutrality sends it: one particular combination of the two concentrations stays fixed. To find it, expand the pinned difference VAVBV_A - V_B with the ideal-dilute relation:

VAVB=(VAVB)+RTzAFlncAcRTzBFlncBc.V_A - V_B = \left(V^\circ_A - V^\circ_B\right) + \frac{RT}{z_A F}\ln\frac{c_A}{c^\circ} - \frac{RT}{z_B F}\ln\frac{c_B}{c^\circ}.

The left side is pinned, and the rung spacing VAVBV^\circ_A - V^\circ_B is a rigid property of the medium, blind to the ladder's offset. The two concentration terms are therefore left holding a constant between them, and exponentiating turns that into

(cAc)1/zA ⁣(cBc)1/zB=K,\left(\frac{c_A}{c^\circ}\right)^{1/z_A}\!\left(\frac{c_B}{c^\circ}\right)^{-1/z_B} = K,

with

K=exp ⁣(F[(VAVB)(VAVB)]RT).K = \exp\!\left(\frac{F\big[(V_A - V_B) - (V^\circ_A - V^\circ_B)\big]}{RT}\right).

KK depends only on the pinned difference VAVBV_A - V_B and a material constant, and so it holds fixed no matter where the ladder happens to float. (More rigorously, activities can be used[1].)

This is a mass action law, though with one important difference from its chemical namesake: nothing here is actually reacting.[2] The two concentrations are tied together purely electrostatically, both responding to the same shift of the ladder. For a +1+1 cation paired with a 1-1 anion both exponents reduce to +1+1, and the law takes the familiar product form,

cAcB=const,c_A\, c_B = \text{const},

so the cation can be enriched only at the anion's expense, and the reverse.

VAV_AVAV^\circ_AVBV^\circ_BVBV_BcA=100.85cc_A = 10^{-0.85}c^\circcB=100.85cc_B = 10^{-0.85}c^\circSolution−0.30−0.25−0.20−0.15−0.10−0.050.000.050.100.150.200.250.30Voltage (V) — arbitrary offset

Two carriers pinned (VAV_A, VBV_B fixed), with a slider for the ViV^\circ_i ladder offset. As the ladder slides, cAc_A and cBc_B trade off against each other, but their mass-action product holds constant. Other ions are not shown; the slider stands in for whatever unseen charges the ladder is neutralizing.

Note the same combination holds just as well for two carriers of the same sign. Pairing Na+\mathrm{Na}^+ with Al3+\mathrm{Al}^{3+}, for instance, gives cNa+/(cAl3+)1/3=constc_{\mathrm{Na}^+}/(c_{\mathrm{Al}^{3+}})^{1/3} = \text{const}. The fractional power[3] records that a slide of the ladder moves each log-concentration in proportion to ziz_i, three times as far for the Al3+\mathrm{Al}^{3+} as for the Na+\mathrm{Na}^+: the combination above is locked, yet the plain ratio between the two ions remains free to range widely. This goes for every pair: if all three of VNa+V_{\mathrm{Na}^+}, VAl3+V_{\mathrm{Al}^{3+}}, VClV_{\mathrm{Cl}^-} were fixed, then we would have a mass-action law pinning all three pairwise products: cNa+/(cAl3+)1/3c_{\mathrm{Na}^+}/(c_{\mathrm{Al}^{3+}})^{1/3} and (cAl3+)1/3cCl(c_{\mathrm{Al}^{3+}})^{1/3}c_{\mathrm{Cl}^-} and cNa+cClc_{\mathrm{Na}^+}c_{\mathrm{Cl}^-}.

This pair of tools, neutrality to place the ladder and the invariants that survive its motion, accounts for a whole family of effects that, in their conventional presentations, appear quite unrelated. Each is a disturbance to the neutrality balance, met by a shift of the ladder.

Dopants

The semiconductor case is the same pair of statements in the physicist's notation, and a simple one. The ionized donors (ND+N_{\mathrm{D}}^+) and acceptors (NAN_{\mathrm{A}}^-) cannot move, so they have no ViV_i of their own; they enter the balance only as background charge:

ρbg=+eND+eNA,\rho_{\mathrm{bg}} = +eN_{\mathrm{D}}^+ -eN_{\mathrm{A}}^- ,

which the mobile charges must compensate for:

ene+enh++ρbg=0.-e\, n_{\mathrm{e}^-} + e\, n_{\mathrm{h}^+} + \rho_{\mathrm{bg}} = 0.

Equilibrium pins the electron and hole rails together, Ve=Vh+V_{\mathrm{e}^-} = V_{\mathrm{h}^+} (the Fermi level), leaving the band-edge ladder free to seek neutrality. For the two-carrier case, the electrostatic mass-action law above gives us a shortcut to help find neutrality:

nenh+=K,n_{\mathrm{e}^-}\, n_{\mathrm{h}^+} = K,

which one usually sees as the semiconductor mass action law, np=ni2np = {n_i}^2, with essentially the same meaning.[4]

The two equations with two unknowns close the system: together they combine into a quadratic equation that places nen_{\mathrm{e}^-} and nh+n_{\mathrm{h}^+}.

Semiconductor Vh+V^\circ_{\mathrm{h}^+}Ve=Vh+V_{\mathrm{e}^-}{=}V_{\mathrm{h}^+}VeV^\circ_{\mathrm{e}^-}Species voltage
Net chargeladder offset00ρbg\rho_{\mathrm{bg}}ρ\rho

← acceptors ($p$)  ·  donors ($n$) →

The same solve for a doped semiconductor. Top: the two carriers share a single pinned rail Ve=Vh+V_{\mathrm{e}^-} = V_{\mathrm{h}^+} (the Fermi level), and it is the band-edge ladder VeV^\circ_{\mathrm{e}^-}, Vh+V^\circ_{\mathrm{h}^+} that floats. Bottom: net charge versus ladder offset. This time the slider sets the fixed dopant charge ρbg\rho_{\mathrm{bg}}, not the offset: raising it lifts the ρ\rho curve bodily, sliding its zero-crossing (the neutral offset) to a new place, and the ladder follows. Donors (ρbg>0\rho_{\mathrm{bg}} > 0) drive it nn-type, lifting the conduction edge toward the rail; acceptors do the reverse.

(More detailed treatments model the dopants as incompletely ionized sites which, while immobile, are thermodynamically active, their occupation equilibrating with VeV_{\mathrm{e}^-}. Finding neutrality in this case no longer reduces to a simple quadratic equation, yet the ladder slide remains the correct way to find it.)

The common-ion effect

The same reasoning recovers a classic result of solution chemistry. Consider a solution saturated with AgCl\mathrm{AgCl}, so that the solid pins the difference VAg+VClV_{\mathrm{Ag}^+} - V_{\mathrm{Cl}^-} at its saturation value: a two-carrier mass-action constraint, with the product cAg+cClc_{\mathrm{Ag}^+} c_{\mathrm{Cl}^-} locked at the solubility product. Now dissolve a little KCl\mathrm{KCl}. The K+\mathrm{K}^+ has no salt of its own to pin it, and so it behaves as a pure positive disturbance to neutrality, a mobile dopant. The ladder shifts to restore the balance, and because the product is locked, that shift enriches Cl\mathrm{Cl}^- while depleting Ag+\mathrm{Ag}^+, which precipitates back onto the solid. (The solid serves as the reservoir that holds the VAg+VClV_{\mathrm{Ag}^+} - V_{\mathrm{Cl}^-} gap rigid throughout.) The result is the textbook common-ion effect, in which added chloride suppresses the solubility of silver, and in the ViV_i picture it appears directly as a shift of the ladder.

Silver wireAgCl coatingSolution drawn +2.30 V from its IUPAC-referenced position (per-species display offset)drawn +2.30 V from its IUPAC-referenced position (per-species display offset)VeV_{\mathrm{e}^{-}}VClV_{\mathrm{Cl}^{-}}VClV_{\mathrm{Cl}^{-}}^\circVAg+V_{\mathrm{Ag}^{+}}VAg+V_{\mathrm{Ag}^{+}}^\circVK+V_{\mathrm{K}^{+}}VK+V_{\mathrm{K}^{+}}^\circVH+V_{\mathrm{H}^{+}}^\circ−1.5−1.0−0.50.00.51.01.5Species Voltage (V) — per-species offsets ⌇

Concentrations (mol/L):
[K⁺]  [Cl⁻]
[Ag⁺]

The common-ion effect in action. The solid pins the VAg+V_{\mathrm{Ag}^+}VClV_{\mathrm{Cl}^-} gap while added KCl supplies a mobile-dopant K+\mathrm{K}^+, shifting the whole ladder: chloride enriched, silver depleted, product locked (readouts below). The VH+V^\circ_{\mathrm{H}^+} rung rides along as a bystander to show the shift is shared by every rung. The checkbox runs the counter-experiment, dosing AgNO3\mathrm{AgNO_3} instead: the same ladder shift in reverse, now suppressing chloride.

Trapped charge and the Donnan potential

If the disturbance is confined to one side of a membrane, the effect becomes spatial. Let a membrane pass the small mobile ions (Na+\mathrm{Na}^+ and Cl\mathrm{Cl}^-) while blocking a large polyion, say a protein Pr\mathrm{Pr}^-, on the right. The mobile ions equilibrate across the membrane, so VNa+V_{\mathrm{Na}^+} and VClV_{\mathrm{Cl}^-} run perfectly flat through it. The trapped Pr\mathrm{Pr}^-, however, is a background charge present on the right alone, and so the ladder must sit at a different offset on each side in order to keep both sides neutral. That step in the ViV^\circ_i ladder across the membrane is the Donnan potential. With the mobile ViV_i flat and the ladder stepped, the concentrations step as well: the protein side is left richer in Na+\mathrm{Na}^+ and poorer in Cl\mathrm{Cl}^-, which is precisely Donnan exclusion.

SalineSaline + Pr⁻ drawn -2.87 V from its IUPAC-referenced position (per-species display offset)drawn -2.87 V from its IUPAC-referenced position (per-species display offset)VClV_{\mathrm{Cl}^{-}}VClV_{\mathrm{Cl}^{-}}^\circVNa+V_{\mathrm{Na}^{+}}^\circVNa+V_{\mathrm{Na}^{+}}0.00.20.40.60.81.01.21.41.6Species Voltage (V) — per-species offsets ⌇

Right side: [Na⁺]  [Cl⁻] mol/L · Donnan step mV

Gibbs–Donnan equilibrium. The mobile-ion traces VNa+V_{\mathrm{Na}^+} and VClV_{\mathrm{Cl}^-} run flat across the membrane; the ViV^\circ_i ladder steps between the two sides to neutralize the trapped protein, and that step is the Donnan potential. A slider sets the blocked-protein concentration on the right. (Saline at 1 mmol/L1\ \mathrm{mmol/L}; the drawn rung spacing is compressed for display, but the gaps and the step are to scale.)

A different solvent on each side

Let the two sides instead be genuinely different solvents, water against an immiscible organic phase such as nitrobenzene. The flat-ViV_i equilibrium remains the same, but now each side has a distinct internal ViV^\circ_i ladder structure. This liquid–liquid interface is known by the acronym ITIES, interface between two immiscible electrolyte solutions. The step in each ion's rung, ΔVi\Delta V^\circ_i, folds together an electrostatic and a solvation contribution; a uniform Donnan shift Δϕ\Delta \phi, which would parse the two apart, is neither well defined here nor pertinent to the thermodynamics.[5] The total steps do carry meaning, measurable in charge-neutral combinations, subject beyond that to the same single-ion ambiguity as the activities of the nonideality topic.

WaterOrganic VC+V_{\mathrm{C}^{+}}^\circVC+V_{\mathrm{C}^{+}}VAV_{\mathrm{A}^{-}}VAV_{\mathrm{A}^{-}}^\circSpecies Voltage (V)

A liquid–liquid interface (ITIES): water against an immiscible organic solvent. A shared ion runs flat across the boundary (ViV_i continuous), but the two solvents' ViV^\circ_i ladders step by a different amount for each ion, set by how each prefers one solvent over the other. The (schematically drawn) electric double layers appear as bends near the junction, which we will cover more in the basic electrostatics topic.

A semiconductor reader will recognize this deforming ladder: a junction between two different media is a heterojunction. There too each species' edge moves by a different amount at the contact (the conduction and valence band offsets differ), and no single potential step can describe what is happening.

Takeaways

Everything in this topic comes down to a single motion: sliding the rigid ViV^\circ_i ladder to the one offset at which the bulk is neutral. The pinned ViV_i values and any fixed background charge together decide where that offset falls. Mass-action products emerge independent of the ladder, since the offset cancels, and dopants, the common-ion effect, and the Donnan potential turn out to be the same neutrality-restoring shift in three different guises. A boundary between two solvents generalizes the picture one step further, to a ladder that deforms as well as shifts.

The Donnan case, in which a fixed background charge on one side of a boundary forces a step in the ladder, has an exact counterpart in the solid state. The next topic sets the two side by side: the electrochemical bipolar membrane and its semiconductor twin, the pn junction.

NEXT TOPIC: Bipolar membranes and pn junctions


  1. In a nonideal solution, replace each ci/cc_i/c^\circ with the activity aia_i. The left-hand combination (aA)1/zA(aB)1/zB(a_A)^{1/z_A}(a_B)^{-1/z_B} is then a mean-activity-type product, which is exactly the charge-balanced, ladder-independent quantity that remains measurable even when single-ion activities are not. In fact this ladder-shifting is exactly the source of the ion-activity ambiguity; see the nonideality topic. ↩︎

  2. The concentrations do change as the ladder slides, so ions must be coming and going somewhere: to and from whatever reservoirs pin the ViV_i. Those exchanges may themselves be reactions (an ion plating onto an electrode, a salt precipitating out), but each ties a carrier to its own reservoir; nothing converts AA into BB. ↩︎

  3. In chemistry one usually clears the fractional exponents by raising the relation to the power zAzBz_A z_B, giving the tidy integer form cNa+3/cAl3+=constc_{\mathrm{Na}^+}^{3}/c_{\mathrm{Al}^{3+}} = \text{const}. It is the same law; the fractional version is just the one written per unit of charge, which is the normalization the ViV_i picture invites. ↩︎

  4. Actually, the way we have defined KK means it would vary out of equilibrium, i.e. in the case of split quasi-Fermi levels VeVh+V_{\mathrm{e}^-} \neq V_{\mathrm{h}^+} as on the previous page. The value ni2{n_i}^2 is strictly the equilibrium value of KK. ↩︎

  5. The electrochemical literature does parse the step, as ΔVi=Δϕ+ΔGtr,i/(ziF)\Delta V^\circ_i = \Delta\phi + \Delta G^\circ_{\mathrm{tr},i}/(z_i F) with the ion's standard Gibbs energy of transfer ΔGtr,i\Delta G^\circ_{\mathrm{tr},i} carrying the solvation part. This requires an extrathermodynamic ϕ\phi convention fixed for both solvents; chemists have generally adopted TATB, which (in effect) anchors each solvent's ϕ\phi midway between the ViV^\circ_i rungs of two bulky ions built to be twins. The TATB twins get their audit in ϕ\phi under the microscope, and the per-solvent ϕ\phi choice is illustrated in Offsets galore. ↩︎