Selectivity Coefficients and ISE Errors

A calcium ISE dropped into seawater will lie to you with total confidence. It reports a calcium activity, producing a clean stable millivolt reading. That number is simply wrong. Magnesium and sodium ions, present at far higher concentration, are climbing partway onto the same ion-exchange sites. These are the sites the electrode uses to sense calcium. The electrode cannot distinguish "a little of what it wants" from "a lot of something it merely tolerates." Quantifying exactly how much an interfering ion corrupts the reading is the entire content of the Nikolsky-Eisenman equation. Learning to compute with it is what separates a technician who reads a meter from an analyst who trusts — or correctly distrusts — a result. This matters because ion-selective electrodes are deployed exactly where interferences are worst. Clinical serum has Na⁺, K⁺, Ca²⁺, Cl⁻ all coexisting at high and variable concentrations. The same is true of industrial brines, soil extracts, and natural waters loaded with competing cations. A K⁺ electrode used in a hospital blood-gas analyzer must report physiologically decisive potassium concentrations. The difference between 3.5 and 6.5 mM is the difference between a normal patient and one in cardiac danger. Meanwhile, sodium sits at roughly twenty-five times the potassium concentration in the same sample. If the selectivity coefficient for sodium interference were not tiny, and not known precisely, the instrument would be useless. This article builds on Ion-Selective Electrode Theory. It also builds on the membrane-potential ideas from The Glass pH Electrode. It turns "the electrode responds mostly to the ion of interest" into an actual quantitative error budget. The big picture is this: every ISE membrane exchange site has some finite affinity for ions other than the target. The Nikolsky-Eisenman equation packages that affinity into a single number per interferent. This number is the selectivity coefficient Ki,jpotK^{\text{pot}}_{i,j}. It converts an interferent's activity into an equivalent amount of the primary ion. Once selectivity coefficients are known (from IUPAC-standardized separate-solution or fixed-interference measurements), several things become possible. You can predict the electrode's detection limit. You can compute the systematic error a given interferent will cause. You can decide whether a method is fit for a real sample matrix at all. This article builds the Nikolsky-Eisenman equation from the mixed-membrane-potential picture. It derives the two standard experimental protocols for measuring Ki,jpotK^{\text{pot}}_{i,j}. It works three realistic examples. These are a foundational selectivity calculation, an applied clinical-serum error budget, and an exam-level detection-limit derivation. It closes with the traps that catch real analysts. These include nonlinear response near the detection limit, and activity-versus-concentration confusion. It also flags the wrong assumption that selectivity coefficients are matrix-independent constants.

Key takeaways
  • The Nikolsky-Eisenman equation extends the Nernst equation to include interfering ions weighted by selectivity coefficients Ki,jpotK^{\text{pot}}_{i,j}.
  • Ki,jpot≪1K^{\text{pot}}_{i,j} \ll 1 means the electrode strongly prefers the primary ion ii over interferent jj; Ki,jpot>1K^{\text{pot}}_{i,j} > 1 means it actually prefers the interferent.
  • Selectivity coefficients are measured by the separate-solution method (SSM) or the fixed-interference method (FIM), and the two can disagree for a nonideal membrane.
  • The detection limit of an ISE is set by the crossover point where the interference term equals the primary-ion term, not by electronic noise.
  • Systematic error from an interferent scales with Ki,jpot ajzi/zjK^{\text{pot}}_{i,j}\,a_j^{z_i/z_j} relative to aia_i — charge mismatch between ions makes this exponent matter enormously.
  • Key idea: a selectivity coefficient converts an interfering ion's activity into an "equivalent activity" of the primary ion, and the size of the resulting error depends on both the coefficient and the ratio of charges through a non-obvious exponent.
The core insight

The membrane cannot tell ions apart perfectly; it reports a Nernstian response to a weighted sum of activities, and the weight for each interferent — the selectivity coefficient — must be measured, not assumed, because it depends on the exchange equilibrium at the membrane surface, not on electronics.

Where this lives today

Clinical blood-gas and electrolyte analyzers (measuring Na⁺, K⁺, Ca²⁺, Cl⁻, and pH simultaneously from a single small blood sample) are built entirely around selectivity-coefficient engineering: valinomycin-based K⁺ membranes are specifically chosen because their KK,NapotK^{\text{pot}}_{\text{K,Na}} is small enough (typically 10−410^{-4}–10−510^{-5}) to survive the roughly 30-to-1 Na⁺-to-K⁺ activity ratio in serum, a calculation directly analogous to Example 2 above but with a membrane engineered to pass the clinical tolerance rather than fail it. Environmental and agricultural laboratories use fluoride, nitrate, and ammonium ISEs in soil extracts and natural waters where sulfate, chloride, and other major ions are present at far higher concentration than the trace analyte — precisely the charge-mismatch detection-limit calculation of Example 3, which is why total ionic strength adjustment buffers (TISAB) are added: they suppress activity-coefficient variability and, in the fluoride case, complex interfering polyvalent cations that would otherwise distort the free fluoride activity itself, on top of any pure selectivity-coefficient interference. Industrial process control uses ISEs for real-time monitoring of electroplating baths, boiler water, and fermentation broths, where the sample matrix is often far more concentrated and far less "clean" than any calibration standard; selectivity-coefficient error budgets, computed exactly as in this article, are used to decide whether an in-line ISE reading can be trusted directly or must be periodically cross-checked against a reference method such as ion chromatography. Biomedical and materials research groups developing new solid-contact and all-solid-state ISEs (replacing the classical liquid internal filling solution with a conducting polymer or ionic-liquid layer) report selectivity coefficients using the same SSM/FIM protocols as a primary performance metric, because a new membrane chemistry is not considered viable for real-sample deployment — for instance wearable sweat-electrolyte sensors — until its selectivity coefficients against the dominant physiological interferents are quantified and shown to meet the same kind of error budget worked through above. This same selectivity-and-error-budget logic carries forward directly into Potentiometric Titrations, where an ISE used as an endpoint-detecting indicator electrode must likewise be checked for interference from other ions present in the titration matrix before its inflection point can be trusted as a true stoichiometric endpoint.

Ask the tutor: If my Ca²⁺ ISE has KCa,Mgpot=0.015K^{\text{pot}}_{\text{Ca,Mg}} = 0.015 and my sample has both Mg²⁺ and Na⁺ present, how do I combine both interference terms correctly in the same detection-limit calculation?

The physical picture

An ion-selective membrane — glass, crystalline (LaF₃, AgCl/Ag₂S), or liquid ion-exchanger in PVC — works by establishing a phase-boundary potential that depends on the activity ratio of an ion between the bulk sample and the membrane surface, exactly as developed for the glass pH electrode. The Nernstian response assumes that only one ion, the "primary" ion ii, participates in this exchange equilibrium. Real membranes are never perfectly selective: exchange sites (fixed anionic groups on a glass surface, or a charged ionophore in a PVC liquid membrane) have some finite affinity for other ions of the same or opposite sign, denoted jj. When both ii and jj exchange onto the membrane simultaneously, the phase-boundary potential is no longer a function of aia_i alone — it becomes a function of an effective activity that blends aia_i and a scaled version of aja_j. The scaling factor is the potentiometric selectivity coefficient Ki,jpotK^{\text{pot}}_{i,j}. Physically it is a ratio of exchange-equilibrium constants and mobilities inside the membrane: how easily jj replaces ii at an exchange site, multiplied by how much more or less mobile jj is once inside the membrane phase, relative to ii. It is not a thermodynamic constant in the strict sense — it depends on membrane composition, ionic strength, and (for imperfect membranes) even on the concentration ratio used to measure it, which is why IUPAC prescribes standardized measurement protocols rather than leaving it undefined. The consequence for the analyst is that every ISE reading is really a reading of an "apparent activity" aiapp=ai+∑jKi,jpot ajzi/zja_i^{\text{app}} = a_i + \sum_j K^{\text{pot}}_{i,j}\,a_j^{z_i/z_j}, and everything downstream — calibration curves, detection limits, systematic bias in real samples — follows from taking this apparent activity seriously.

Figure
Mixed potential at an ion-exchange membrane: primary ion i and interferent j both partition onto exchange sites, and the phase-boundary potential responds to the weighted sum of their activities
membrane sample solution internal filling soln. i i j i j exchange sites i E = const − (RT/z_iF) ln[ a_i + Σ K(i,j) a_j^(z_i/z_j) ]

The equation at the bottom of that figure is the Nikolsky-Eisenman relation itself, and it is worth reading slowly before deriving it: the electrode measures a single effective activity that is a sum, not a product, of the primary and (scaled) interferent activities — meaning interference is additive on the activity scale, not multiplicative on the potential scale.

Figure
Selectivity-error map: reading error (mV and % activity) as K_pot and the interferent-to-primary activity ratio vary, for z_i = z_j = 1
Error from interferent j (same charge as i, z=1, 25°C) K_pot · a_j / a_i (log scale, decades) % error in reported a_i 0.001 0.01 0.1 1.0 ~11% at ratio 0.1 100% at ratio 1.0 error grows steeply once K·a_j ≳ 0.1·a_i

The core ideas

The full Nikolsky-Eisenman equation for a primary ion ii of charge ziz_i in the presence of interferents jj of charge zjz_j:

E=L−RTziFln⁡ ⁣(ai+∑jKi,jpot aj zi/zj)E = L - \frac{RT}{z_i F}\ln\!\left(a_i + \sum_j K^{\text{pot}}_{i,j}\, a_j^{\,z_i/z_j}\right)

where:

  • EE is the measured cell potential (V), LL a constant absorbing the reference electrode, asymmetry, and standard potentials,
  • R=8.314 J mol−1K−1R = 8.314\ \text{J mol}^{-1}\text{K}^{-1}, TT the absolute temperature (K), F=96485 C mol−1F = 96485\ \text{C mol}^{-1},
  • aia_i the activity of the primary ion, aja_j the activity of each interferent,
  • zi,zjz_i, z_j the signed charges of primary ion and interferent (same sign required for interference to occur at all — a cation-selective membrane is not perturbed by anions in this term),
  • Ki,jpotK^{\text{pot}}_{i,j} the (dimensionless, ion-pair-specific) potentiometric selectivity coefficient.

Reading it: the electrode behaves as though it were seeing a single ion at "apparent activity" ai+∑jKi,jpotajzi/zja_i + \sum_j K^{\text{pot}}_{i,j} a_j^{z_i/z_j}; the exponent zi/zjz_i/z_j is required so that the interferent term has the same dimensional weight in the exponential response as the primary-ion term — a doubly charged interferent competing with a singly charged primary ion must be raised to the ½ power. The systematic relative error caused by an interferent, expressed as a fractional error in the reported activity of ii (assuming a single interferent for clarity):

Δaiai=Ki,jpot aj zi/zjai\frac{\Delta a_i}{a_i} = \frac{K^{\text{pot}}_{i,j}\, a_j^{\,z_i/z_j}}{a_i}

where:

  • Δai\Delta a_i is the difference between the apparent (reported) activity and the true activity of ii,
  • all other symbols as above.

Reading it: this is the ratio of the two terms inside the logarithm's argument. It says the relative error is what matters, not the absolute size of Ki,jpotK^{\text{pot}}_{i,j} alone — a selectivity coefficient of 10−310^{-3} is harmless against an interferent at the same activity as the analyte, but ruinous against an interferent present at 10410^4 times the analyte's activity. The detection limit aiDLa_i^{\text{DL}} — the lowest activity of ii at which the calibration curve is still Nernstian rather than dominated by interference — is IUPAC-defined as the activity at the intersection of the two asymptotes of the calibration curve, which in the single-interferent, equal-charge case reduces to:

aiDL=Ki,jpot aj zi/zja_i^{\text{DL}} = K^{\text{pot}}_{i,j}\, a_j^{\,z_i/z_j}

where:

  • aiDLa_i^{\text{DL}} is the primary-ion activity at which the true term equals the interference term (the intersection of the low-concentration and high-concentration asymptotes on a semilog plot of EE vs. log⁡ai\log a_i).

Reading it: you cannot measure aia_i reliably below the level at which the interferent's contribution to the "apparent activity" equals the primary ion's own contribution — below this point, the electrode's response flattens and becomes governed entirely by the constant background of aja_j.

Derivation AFrom membrane exchange equilibria to the Nikolsky-Eisenman sum

Goal: show that a membrane exchanging two ions of possibly different charge produces exactly the additive, exponent-weighted form above, starting from the boundary potential and the ion-exchange equilibrium at the membrane surface.

Step 1 — Boundary potential for a single ion. For a membrane selective to ion ii alone, the phase-boundary (Donnan-type) potential is Nernstian:

Ei=const−RTziFln⁡aiE_i = \text{const} - \frac{RT}{z_i F}\ln a_i

This follows from equating the electrochemical potential of ii in the sample and in the membrane surface layer at equilibrium, exactly as in the derivation of the Nernstian glass electrode response.

Step 2 — Two ions competing for the same sites. Suppose the membrane also allows jj to occupy the same exchange sites via the ion-exchange reaction

i(aq)+j(mem)⇌i(mem)+j(aq),Kij=ai,mem aj,aqai,aq aj,memi(\text{aq}) + j(\text{mem}) \rightleftharpoons i(\text{mem}) + j(\text{aq}), \qquad K_{ij} = \frac{a_{i,\text{mem}}\,a_{j,\text{aq}}}{a_{i,\text{aq}}\,a_{j,\text{mem}}}

At the membrane surface, both ions establish local equilibrium with the sample simultaneously, and both contribute charge-carrying flux across the boundary layer. The total boundary potential can no longer depend on aia_i alone; it must depend on the effective "activity" that the membrane surface presents to the bulk of the membrane, which is a linear combination of the two ions' membrane-phase activities weighted by their mobilities uiu_i, uju_j inside the membrane phase.

Step 3 — Combining exchange equilibrium with mobility weighting. The mixed boundary potential (the full Nikolsky-Eisenman derivation, following the treatment of ion-exchanger membranes) takes the form

E=const−RTziFln⁡ ⁣(ai+ujKijui aj zi/zj)E = \text{const} - \frac{RT}{z_i F}\ln\!\left(a_i + \frac{u_j K_{ij}}{u_i}\, a_j^{\,z_i/z_j}\right)

The exponent zi/zjz_i/z_j arises because the exchange equilibrium constant KijK_{ij} relates activities raised to powers matching the stoichiometry of charge balance at the exchange site: exchanging one izi+i^{z_i+} ion requires zi/zjz_i/z_j ions of jzj+j^{z_j+} to maintain electroneutrality at the site, so the mass-action expression for KijK_{ij} naturally produces aja_j raised to that fractional power once solved for the membrane-phase activity of jj in terms of aia_i and aja_j.

Step 4 — Defining the selectivity coefficient. Everything that is membrane-specific and not easily separated experimentally — the true exchange constant KijK_{ij} and the mobility ratio uj/uiu_j/u_i — is collected into the single empirical quantity

Ki,jpot≡ujui KijK^{\text{pot}}_{i,j} \equiv \frac{u_j}{u_i}\,K_{ij}

so that the boundary potential is exactly the Nikolsky-Eisenman form written above. This is why Ki,jpotK^{\text{pot}}_{i,j} is called potentiometric: it is defined operationally through the potential response, not derived independently from a pure thermodynamic exchange constant, and so must be measured by electrochemical experiment.

Step 5 — Generalizing to many interferents. Because each interferent occupies exchange sites largely independently in the dilute-interference limit, the single-interferent term generalizes to a sum over all species sharing the membrane's charge sign, giving the full equation used throughout this article.

Go deeper · Derivation B · 1 minDetection limit as the intersection of two Nernstian asymptotes
Derivation BDetection limit as the intersection of two Nernstian asymptotes

Goal: derive the IUPAC detection-limit formula from the calibration curve directly, and show it collapses correctly to the pure-Nernstian limit when interference is negligible.

Step 1 — Set up the calibration curve. Hold the interferent activity fixed at some background level aja_j (as in a real sample matrix) and vary aia_i over many decades while recording EE. The Nikolsky-Eisenman equation predicts

E(ai)=L−RTziFln⁡ ⁣(ai+Ki,jpot aj zi/zj)E(a_i) = L - \frac{RT}{z_i F}\ln\!\left(a_i + K^{\text{pot}}_{i,j}\,a_j^{\,z_i/z_j}\right)

with LL and Ki,jpot ajzi/zjK^{\text{pot}}_{i,j}\,a_j^{z_i/z_j} both constant as aia_i is varied (the interferent's contribution is fixed because its own activity is held fixed).

Step 2 — High-activity asymptote. When ai≫Ki,jpot ajzi/zja_i \gg K^{\text{pot}}_{i,j}\,a_j^{z_i/z_j}, the constant term inside the logarithm is negligible and

E≈L−RTziFln⁡aiE \approx L - \frac{RT}{z_i F}\ln a_i

This is the ordinary Nernstian line: a straight segment of slope −RT/(ziF)⋅ln⁡(10)-RT/(z_iF)\cdot\ln(10) per decade on a plot of EE vs. log⁡ai\log a_i (about 59.2/zi59.2/z_i mV per decade at 25 °C).

Step 3 — Low-activity asymptote. When ai≪Ki,jpot ajzi/zja_i \ll K^{\text{pot}}_{i,j}\,a_j^{z_i/z_j}, the aia_i term is negligible and

E≈L−RTziFln⁡ ⁣(Ki,jpot aj zi/zj)E \approx L - \frac{RT}{z_i F}\ln\!\left(K^{\text{pot}}_{i,j}\,a_j^{\,z_i/z_j}\right)

which does not depend on aia_i at all — a horizontal plateau. The electrode has become, in effect, a (poor) electrode for the interferent, or simply insensitive to further dilution of the target.

Step 4 — Intersection defines the detection limit. The two asymptotic lines (Step 2's sloped line, Step 3's flat line) are both simple functions of log⁡ai\log a_i; IUPAC defines the detection limit as their point of intersection. Setting the two expressions for EE equal:

−RTziFln⁡ai=−RTziFln⁡ ⁣(Ki,jpot aj zi/zj)-\frac{RT}{z_i F}\ln a_i = -\frac{RT}{z_i F}\ln\!\left(K^{\text{pot}}_{i,j}\,a_j^{\,z_i/z_j}\right)

The prefactors and the constant LL cancel from both sides (they were common to Steps 2 and 3), leaving simply

ai=Ki,jpot aj zi/zja_i = K^{\text{pot}}_{i,j}\, a_j^{\,z_i/z_j}

Relabeling this crossover value of aia_i as aiDLa_i^{\text{DL}} gives exactly the detection-limit formula stated earlier.

Step 5 — Limiting-case check. As Ki,jpot→0K^{\text{pot}}_{i,j} \to 0 (a perfectly selective membrane), aiDL→0a_i^{\text{DL}} \to 0: an ideal electrode could in principle measure arbitrarily dilute primary ion, and the calibration curve would be Nernstian all the way down, consistent with there being no interference term left to compete with aia_i. This confirms the formula behaves correctly in the limit where it must.

Worked examples

Example 1 — Foundational: computing a selectivity coefficient from the separate-solution method (SSM). A calcium ISE (divalent primary ion, zCa=+2z_{\text{Ca}} = +2) is calibrated in a pure CaCl₂ solution and gives E1=−18.4 mVE_1 = -18.4\ \text{mV} (vs. a fixed reference) at aCa=1.00×10−2a_{\text{Ca}} = 1.00\times10^{-2}. In a separate pure NaCl solution (zNa=+1z_{\text{Na}} = +1) at aNa=1.00×10−1a_{\text{Na}} = 1.00\times10^{-1}, the same electrode gives E2=+42.6 mVE_2 = +42.6\ \text{mV}. Find KCa,NapotK^{\text{pot}}_{\text{Ca,Na}} at 25.0 °C. The SSM protocol: measure EE in a solution of pure ii at activity aia_i, then separately in a solution of pure jj at activity aja_j, and use the Nikolsky-Eisenman equation applied to each solution separately (in the second solution, ai=0a_i = 0, so the electrode responds purely as if to Ki,jpot ajzi/zjK^{\text{pot}}_{i,j}\,a_j^{z_i/z_j}). For the pure-Ca solution:

E1=L−RT2Fln⁡aCaE_1 = L - \frac{RT}{2F}\ln a_{\text{Ca}}

For the pure-Na solution, treating the electrode as responding to the "equivalent Ca activity" produced by Na:

E2=L−RT2Fln⁡ ⁣(KCa,Napot aNa 2/1)E_2 = L - \frac{RT}{2F}\ln\!\left(K^{\text{pot}}_{\text{Ca,Na}}\, a_{\text{Na}}^{\,2/1}\right)

Subtracting the second equation from the first eliminates LL:

E1−E2=−RT2Fln⁡ ⁣(aCaKCa,Napot aNa2)E_1 - E_2 = -\frac{RT}{2F}\ln\!\left(\frac{a_{\text{Ca}}}{K^{\text{pot}}_{\text{Ca,Na}}\,a_{\text{Na}}^{2}}\right)

At 25.0 °C, RT/F=0.025693 VRT/F = 0.025693\ \text{V}, so RT/(2F)=0.012846 V=12.846 mVRT/(2F) = 0.012846\ \text{V} = 12.846\ \text{mV}. Compute the left side:

E1−E2=−18.4−42.6=−61.0 mVE_1 - E_2 = -18.4 - 42.6 = -61.0\ \text{mV}

Solve for the log term:

ln⁡ ⁣(aCaKCa,Napot aNa2)=−(E1−E2)12.846 mV=61.012.846=4.749\ln\!\left(\frac{a_{\text{Ca}}}{K^{\text{pot}}_{\text{Ca,Na}}\,a_{\text{Na}}^{2}}\right) = \frac{-(E_1-E_2)}{12.846\ \text{mV}} = \frac{61.0}{12.846} = 4.749

Exponentiate:

aCaKCa,Napot aNa2=e4.749=115.4\frac{a_{\text{Ca}}}{K^{\text{pot}}_{\text{Ca,Na}}\,a_{\text{Na}}^{2}} = e^{4.749} = 115.4

Now insert aCa=1.00×10−2a_{\text{Ca}} = 1.00\times10^{-2} and aNa=1.00×10−1a_{\text{Na}} = 1.00\times10^{-1}, so aNa2=1.00×10−2a_{\text{Na}}^2 = 1.00\times10^{-2}:

KCa,Napot=aCa115.4×aNa2=1.00×10−2115.4×1.00×10−2=8.66×10−3K^{\text{pot}}_{\text{Ca,Na}} = \frac{a_{\text{Ca}}}{115.4 \times a_{\text{Na}}^2} = \frac{1.00\times10^{-2}}{115.4 \times 1.00\times10^{-2}} = 8.66\times10^{-3}

So KCa,Napot≈8.7×10−3K^{\text{pot}}_{\text{Ca,Na}} \approx 8.7\times10^{-3} (two sig figs, matching the precision of the millivolt readings). This is a physically reasonable value for a Ca-ISE with a phosphate-ester liquid ion-exchanger membrane, whose real published selectivity coefficients for Na⁺ interference typically fall in the 10−210^{-2}–10−310^{-3} range.

Limiting-case check: if the electrode were perfectly selective (KCa,Napot→0K^{\text{pot}}_{\text{Ca,Na}} \to 0), E2E_2 would have to approach −∞-\infty for a fixed finite aNaa_{\text{Na}} — the pure-interferent solution would generate no measurable signal at all, which matches the physical intuition that a perfectly selective membrane is utterly blind to the interferent.

Example 2 — Applied: error budget for a K⁺ ISE in serum. A potassium ISE has manufacturer-quoted selectivity coefficient KK,Napot=1.0×10−3K^{\text{pot}}_{\text{K,Na}} = 1.0\times10^{-3} (both monovalent, zK=zNa=+1z_{\text{K}} = z_{\text{Na}} = +1, so the exponent zi/zj=1z_i/z_j = 1). A patient's serum has aK=4.20×10−3a_{\text{K}} = 4.20\times10^{-3} (physiological ~4.2 mM, assuming activity ≈ concentration for this estimate at physiological ionic strength with a fixed ionic-strength background already accounted for in calibration) and aNa=1.40×10−1a_{\text{Na}} = 1.40\times10^{-1} (140 mM sodium, a typical serum value). Find the absolute and relative error in the reported potassium activity, and check whether this error is clinically tolerable (clinical labs typically require better than ±2% accuracy on serum K⁺). The apparent activity the electrode reports:

aKapp=aK+KK,Napot aNa 1a_{\text{K}}^{\text{app}} = a_{\text{K}} + K^{\text{pot}}_{\text{K,Na}}\,a_{\text{Na}}^{\,1}

Compute the interference term:

KK,Napot aNa=(1.0×10−3)(1.40×10−1)=1.40×10−4K^{\text{pot}}_{\text{K,Na}}\,a_{\text{Na}} = (1.0\times10^{-3})(1.40\times10^{-1}) = 1.40\times10^{-4}

So:

aKapp=4.20×10−3+1.40×10−4=4.34×10−3a_{\text{K}}^{\text{app}} = 4.20\times10^{-3} + 1.40\times10^{-4} = 4.34\times10^{-3}

The absolute error is ΔaK=1.40×10−4\Delta a_{\text{K}} = 1.40\times10^{-4}, and the relative error is:

ΔaKaK=1.40×10−44.20×10−3=3.33×10−2=3.33%\frac{\Delta a_{\text{K}}}{a_{\text{K}}} = \frac{1.40\times10^{-4}}{4.20\times10^{-3}} = 3.33\times10^{-2} = 3.33\%

This 3.3% positive bias exceeds the ±2% clinical tolerance — meaning that a K⁺ ISE with this selectivity coefficient, used naively on undiluted serum, would systematically over-report potassium by a clinically significant margin at the low end of the normal range, which is exactly the scenario clinical instruments must guard against (real clinical analyzers use ionophores, typically valinomycin-based membranes, engineered for KK,Napot∼10−4K^{\text{pot}}_{\text{K,Na}} \sim 10^{-4} to 10−510^{-5}, specifically to survive this test).

Limiting-case check: if aNaa_{\text{Na}} were reduced to zero (a hypothetical sodium-free sample), the interference term vanishes exactly and aKapp→aKa_{\text{K}}^{\text{app}} \to a_{\text{K}}, recovering the ideal Nernstian reading — confirming the error scales linearly with the interferent activity as it must for two ions of equal charge.

Example 3 — Exam-level: detection limit with a charge mismatch, and back-calculating a required selectivity coefficient. A fluoride ISE (LaF₃ crystal membrane, zF=−1z_{\text{F}} = -1) is to be used in groundwater containing a fixed background of sulfate interference, aSO4=5.0×10−4a_{\text{SO}_4} = 5.0\times10^{-4} (zSO4=−2z_{\text{SO}_4} = -2). A groundwater standards program requires the electrode to reliably measure fluoride down to aF=1.0×10−6a_{\text{F}} = 1.0\times10^{-6} (roughly 0.019 mg/L, a low but environmentally relevant level). What is the maximum tolerable KF,SO4potK^{\text{pot}}_{\text{F,SO}_4} for the detection limit to sit at or below this required level? Then verify what detection limit results if the electrode actually has KF,SO4pot=2.0×10−4K^{\text{pot}}_{\text{F,SO}_4} = 2.0\times10^{-4} (a realistic literature value for fluoride electrodes' generally excellent rejection of sulfate). First, note the exponent required: with zi=−1z_i = -1 (fluoride) and zj=−2z_j = -2 (sulfate), the Nikolsky-Eisenman exponent is

zizj=−1−2=12\frac{z_i}{z_j} = \frac{-1}{-2} = \frac{1}{2}

so the detection-limit formula becomes:

aFDL=KF,SO4pot aSO4 1/2a_{\text{F}}^{\text{DL}} = K^{\text{pot}}_{\text{F,SO}_4}\, a_{\text{SO}_4}^{\,1/2}

Part A — required selectivity coefficient. Require aFDL≤1.0×10−6a_{\text{F}}^{\text{DL}} \le 1.0\times10^{-6}:

KF,SO4pot≤aFDLaSO41/2=1.0×10−6(5.0×10−4)1/2K^{\text{pot}}_{\text{F,SO}_4} \le \frac{a_{\text{F}}^{\text{DL}}}{a_{\text{SO}_4}^{1/2}} = \frac{1.0\times10^{-6}}{(5.0\times10^{-4})^{1/2}}

Compute the square root:

(5.0×10−4)1/2=2.236×10−2(5.0\times10^{-4})^{1/2} = 2.236\times10^{-2}

So:

KF,SO4pot≤1.0×10−62.236×10−2=4.47×10−5K^{\text{pot}}_{\text{F,SO}_4} \le \frac{1.0\times10^{-6}}{2.236\times10^{-2}} = 4.47\times10^{-5}

The membrane must reject sulfate to better than KF,SO4pot=4.5×10−5K^{\text{pot}}_{\text{F,SO}_4} = 4.5\times10^{-5} for the stated detection limit to be achievable in this matrix.

Part B — actual detection limit with the literature value. With KF,SO4pot=2.0×10−4K^{\text{pot}}_{\text{F,SO}_4} = 2.0\times10^{-4}:

aFDL=(2.0×10−4)(2.236×10−2)=4.47×10−6a_{\text{F}}^{\text{DL}} = (2.0\times10^{-4})(2.236\times10^{-2}) = 4.47\times10^{-6}

This is roughly 4.5 times higher (worse) than the required 1.0×10−61.0\times10^{-6} — so the standard literature LaF₃ electrode, with this selectivity coefficient, would fail to meet the required low-level detection limit in this particular sulfate-rich groundwater; a fluoride-selective preconcentration step or an ionic-strength-adjustment buffer that also complexes or dilutes sulfate would be needed, or a lower-interference membrane batch must be sourced.

Limiting-case check: notice the half-power exponent means the detection limit depends only on the square root of the sulfate activity — a tenfold increase in sulfate activity worsens the detection limit by only about 10≈3.16×\sqrt{10} \approx 3.16\times, not tenfold, precisely because of the charge mismatch built into the Nikolsky-Eisenman exponent. If sulfate activity were instead raised to a full power (as it would be for an equal-charge interferent), the detection limit would degrade proportionally faster — confirming that the fractional exponent is doing real, checkable work in tempering the impact of a highly charged interferent.

Limiting cases, scope & common traps

Regime / assumptionWhat happensConsequence for the analyst
Ki,jpot→0K^{\text{pot}}_{i,j} \to 0Interference term vanishes; Nikolsky-Eisenman reduces to pure Nernst equationIdeal, fully selective electrode; detection limit set only by membrane solubility/leaching, not by this interferent
aj→0a_j \to 0 (interferent absent)Same reduction regardless of Ki,jpotK^{\text{pot}}_{i,j}Calibration in "clean" standards can hide an interference problem that appears only in real samples
zi=zjz_i = z_j (equal charge)Exponent zi/zj=1z_i/z_j = 1; error scales linearly with aja_jSimplest, most common textbook case (e.g., K⁺ vs. Na⁺)
zi≠zjz_i \ne z_j (charge mismatch)Fractional or multiple-valued exponent; error scales sublinearly or superlinearly with aja_jEasy to mis-apply the linear formula by forgetting the exponent — a frequent exam and lab error
ai≫Ki,jpotajzi/zja_i \gg K^{\text{pot}}_{i,j}a_j^{z_i/z_j}High end of calibration curve, fully NernstianNormal operating region; errors from interference negligible
ai≪Ki,jpotajzi/zja_i \ll K^{\text{pot}}_{i,j}a_j^{z_i/z_j}Plateau region below detection limitReading is meaningless — governed by interferent background, not analyte
SSM vs. FIM disagreementSelectivity coefficients from separate-solution and fixed-interference methods differ for nonideal membranesAlways state which protocol produced a quoted Ki,jpotK^{\text{pot}}_{i,j}; do not mix protocols when predicting real-sample error
Activity vs. concentrationNikolsky-Eisenman is strictly an activity equationIgnoring activity coefficients in high-ionic-strength samples (serum, brine) silently injects extra "error" that looks like bad selectivity but is really an uncorrected activity coefficient

The chemist behind it

The Nikolsky-Eisenman equation carries two names because it was built in two historical stages, separated by more than two decades and by the Cold War's division of electrochemical research communities. Boris Petrovich Nikolsky was a Soviet physical chemist working at Leningrad State University in the 1930s, at a time when the glass electrode's Nernstian pH response was already an established laboratory tool but its systematic deviations at high pH and in the presence of alkali-metal ions were not understood mechanistically. Nikolsky proposed, based on an ion-exchange model of the glass surface, that the glass membrane's potential should respond not to hydrogen ion activity alone but to a weighted sum of hydrogen and interfering alkali-metal-ion activities — introducing what is now recognized as the first selectivity-coefficient formalism, originally derived specifically for the H⁺/Na⁺ "alkaline error" of glass pH electrodes. His derivation, grounded in the same ion-exchange equilibrium logic used above, correctly predicted that the alkaline error should worsen at high pH (low H⁺ activity) and high Na⁺ activity — exactly the detection-limit crossover behavior this article derives in general form. The equation's generalization to arbitrary ion-selective membranes — crystalline, liquid-ion-exchanger, and neutral-carrier types alike, with arbitrary charges ziz_i and zjz_j — came from the American electrochemist James W. Ross and collaborators, alongside the broader systematization efforts of Garry A. Rechnitz and others through the 1960s, as the ISE field exploded beyond glass pH electrodes into calcium-, fluoride-, and later ammonia- and CO₂-gas-sensing electrodes. Ross's 1966 work on a practical liquid-membrane calcium electrode, and the accompanying formal extension of Nikolsky's exchange-equilibrium treatment to non-equal charges, produced the exponent zi/zjz_i/z_j in essentially the modern form. Eisenman's name entered through parallel theoretical work on glass and liquid-membrane ion selectivity mechanisms, particularly his systematic classification of alkali-metal cation selectivity sequences in glass electrodes, which gave the physical-chemical foundation for why particular membranes prefer particular ions — connecting the empirical selectivity coefficient to the underlying ion-exchange energetics.

Origins · 1 minBoris Petrovich Nikolsky and how this was discovered
BP
Boris Petrovich Nikolsky
1900–1990 · Leningrad State University, Soviet physical and radiochemistry

Nikolsky's career spanned physical chemistry, radiochemistry, and later nuclear chemistry within the Soviet academic system; his 1930s treatment of the glass-electrode alkaline error, built from an ion-exchange model of the silicate glass surface, is the direct theoretical ancestor of the general selectivity-coefficient formalism used across all modern ion-selective electrodes, from pH glass membranes to solid-state fluoride and calcium liquid-membrane sensors.

History
- 1930s — Nikolsky develops the ion-exchange model of the glass electrode and derives the first selectivity-weighted potential equation for the H⁺/Na⁺ alkaline error.
- 1950s — Eisenman's systematic studies of alkali-metal selectivity in glass membranes connect measured selectivity sequences to ion-exchange site energetics.
- 1966 — Ross reports a practical liquid ion-exchanger calcium-selective electrode, requiring the selectivity formalism to be generalized beyond equal-charge interferents.
- Late 1960s–1970s — IUPAC and the electroanalytical community standardize the separate-solution method (SSM) and fixed-interference method (FIM) for measuring Ki,jpotK^{\text{pot}}_{i,j} reproducibly across laboratories.
- 1976 onward — IUPAC recommendations formalize the detection-limit definition as the intersection of Nernstian and interference-plateau asymptotes, used industry-wide today.
What exactly is a "potentiometric selectivity coefficient," in plain terms?

It is a single number that tells you how much of an interfering ion's activity "looks like" the primary ion to the electrode. If Ki,jpot=10−3K^{\text{pot}}_{i,j} = 10^{-3}, then roughly a thousand-fold excess of interferent jj is needed to produce the same signal contribution as the primary ion ii at unit activity ratio. It is measured, not calculated from first principles, because it depends on the specific chemistry of the membrane's exchange sites and on ion mobilities inside the membrane phase — quantities that are not easily predicted from bulk solution thermodynamics alone.

Why does the exponent zi/zjz_i/z_j appear instead of a simple linear weighting?

The exponent comes from the stoichiometry of charge balance at the membrane's exchange sites: replacing one izi+i^{z_i+} ion at a site requires exactly zi/zjz_i/z_j ions of jzj+j^{z_j+} to keep the site electrically neutral during the exchange. When the ion-exchange equilibrium expression is solved for the interferent's effective contribution to the membrane's surface activity, this stoichiometric ratio surfaces directly as a fractional or multiple-valued power on aja_j. It is not an approximation — it is a direct consequence of charge conservation at the exchange equilibrium.

Are selectivity coefficients ever greater than 1?

Yes, and when this happens it means the membrane actually prefers the "interferent" over the nominal primary ion. This is common in poorly designed or poorly chosen membranes, and it is also exploited deliberately: some anion-selective electrodes are marketed as "nitrate electrodes" but have KNO3,ClO4pot>1K^{\text{pot}}_{\text{NO}_3,\text{ClO}_4} > 1, meaning they respond more strongly to perchlorate than to nitrate — useful information for anyone attempting to use such an electrode in a perchlorate-contaminated sample, and a caution against assuming a marketed ion name guarantees dominant selectivity for that ion.

How does temperature affect the Nikolsky-Eisenman equation and selectivity coefficients?

Temperature enters explicitly through the RT/ziFRT/z_iF prefactor (about 59.2 mV/decade for a monovalent ion at 25 °C, scaling linearly with absolute temperature), so the Nernstian slope itself is temperature-dependent and must be recalibrated or thermostatted for precise work. The selectivity coefficient Ki,jpotK^{\text{pot}}_{i,j} itself can also be weakly temperature-dependent, since it reflects an exchange equilibrium constant and a mobility ratio, both of which are functions of temperature, though this secondary effect is usually smaller than the direct Nernstian-slope temperature dependence for well-behaved membranes.

Does ionic strength affect selectivity coefficients, separately from the activity coefficients of the ions themselves?

Ionic strength affects the measured potential through the activity coefficients of both ii and jj (via something like the Debye-Hückel or Davies equation), which must be corrected for independently of the selectivity coefficient itself. However, very high ionic strength can also change the membrane's own exchange-site environment (swelling, site availability), which can shift the "true" Ki,jpotK^{\text{pot}}_{i,j} measured under those conditions — this is one reason calibration standards should be prepared at an ionic strength resembling the real sample, using an ionic strength adjustment buffer, rather than assuming a selectivity coefficient measured in dilute standards transfers unchanged to a concentrated real matrix.

Why does the detection limit depend on the interferent's activity rather than being a fixed property of the electrode alone?

Because the detection limit is defined as the crossover point where the interference term equals the primary-ion term, and the interference term is Ki,jpot ajzi/zjK^{\text{pot}}_{i,j}\,a_j^{z_i/z_j} — which obviously depends on how much interferent is actually present in the sample. The same electrode will have an excellent (low) detection limit in a clean matrix and a much worse (higher) detection limit in a matrix loaded with interferents, which is why manufacturer-quoted detection limits (usually measured in idealized low-interference standards) can be optimistic relative to performance in a real, dirty sample.

How is this connected to potentiometric titrations, where an ISE is used only to find an endpoint rather than an absolute concentration?

In a titration, the ISE tracks how activity changes as titrant is added, and the steep inflection at the equivalence point is usually far more robust to a modest selectivity-coefficient error than an absolute single-point activity measurement would be, because the interferent background is typically constant throughout the titration and mostly cancels in the derivative or inflection-point analysis. However, if the interferent's activity itself changes substantially during the titration (for example, if the titrant reacts with both the analyte and an interfering ion), Nikolsky-Eisenman-type errors can distort the shape of the titration curve near the endpoint, which is why interference is still checked even in electrodes used only for qualitative endpoint detection.

Before this

Where this leads