In theoretical chemistry, Specific ion Interaction Theory (SIT theory) is a theory used to estimate single-ion activity coefficients in electrolyte solutions at relatively high concentrations.^[1][2] It does so by taking into consideration interaction coefficients between the various ions present in solution. Interaction coefficients are determined from equilibrium constant values obtained with solutions at various ionic strengths. The determination of SIT interaction coefficients also yields the value of the equilibrium constant at infinite dilution.

Background

The need for this theory arises from the need to derive activity coefficients of solutes when their concentrations are too high to be predicted accurately by Debye–Hückel theory. These activity coefficients are needed because an equilibrium constant is defined in thermodynamics as a quotient of activities but is usually measured using concentrations. The protonation of a monobasic acid will be used to simplify the exposition. The equilibrium for protonation of the conjugate base, A⁻ of the acid, may be written as

${\displaystyle {\ce {H+ + A- <=> HA))}$

for which

${\displaystyle K={\frac ((\ce {\{HA\))))((\ce {\{A^-\}\{H^+\))))))$

where {HA} signifies an activity of the chemical species HA etc.. The role of water in the equilibrium has been ignored as in all but the most concentrated solutions the activity of water is a constant. K is defined here as an association constant, the reciprocal of an acid dissociation constant.

Each activity term can be expressed as the product of a concentration and an activity coefficient. For example,

${\displaystyle \{HA\}=[HA]\times \gamma _{HA))$

where the square brackets signify a concentration and γ is an activity coefficient. Thus the equilibrium constant can be expressed as a product of a concentration quotient and an activity coefficient quotient.

${\displaystyle K={\frac ((\ce {[HA]))}((\ce {[H^+][A^-]))))\times {\frac {\gamma _((\ce {HA)))){\gamma _((\ce {H^+))}\gamma _((\ce {A^-))))))$

Taking logarithms.

${\displaystyle \log K=\log K^{0}+\log \gamma _((\ce {HA))}-\log \gamma _((\ce {H^+))}-\log \gamma _((\ce {A^-))))$

K⁰ is the hypothetical value that the equilibrium constant would have if the solution of the acid were so dilute that the activity coefficients were all equal to one.

It is common practise to determine equilibrium constants in solutions containing an electrolyte at high ionic strength such that the activity coefficients are effectively constant. However, when the ionic strength is changed the measured equilibrium constant will also change, so there is a need to estimate individual (single ion) activity coefficients. Debye–Huckel theory provides a means to do this, but it is accurate only at very low concentrations. Hence the need for an extension to Debye–Hückel theory. Two main approaches have been used. SIT theory, discussed here and Pitzer equations.^[3][4]

Development

SIT theory was first proposed by Brønsted^[5] and was further developed by Guggenheim.^[1] Scatchard^[6] extended the theory to allow the interaction coefficients to vary with ionic strength. The theory was mainly of theoretical interest until 1945 because of the difficulty of determining equilibrium constants before the glass electrode was invented. Subsequently, Ciavatta^[2] developed the theory further.

The activity coefficient of the jth ion in solution is written as γ_j when concentrations are on the molal concentration scale and as y_j when concentrations are on the molar concentration scale. (The molality scale is preferred in thermodynamics because molal concentrations are independent of temperature). The basic idea of SIT theory is that the activity coefficient can be expressed as

${\displaystyle \log \gamma _{j}=-z_{j}^{2}{\frac {0.51{\sqrt {I))}{1+1.5{\sqrt {I))))+\sum _{k}\epsilon _{jk}m_{k))$ (molalities)

or

${\displaystyle \log y_{j}=-z_{j}^{2}{\frac {0.51{\sqrt {I))}{1+1.5{\sqrt {I))))+\sum _{k}b_{jk}c_{k))$ (molar concentrations)

where z is the electrical charge on the ion, I is the ionic strength, ε and b are interaction coefficients and m and c are concentrations. The summation extends over the other ions present in solution, which includes the ions produced by the background electrolyte. The first term in these expressions comes from Debye–Hückel theory. The second term shows how the contributions from "interaction" are dependent on concentration. Thus, the interaction coefficients are used as corrections to Debye–Hückel theory when concentrations are higher than the region of validity of that theory.

The activity coefficient of a neutral species can be assumed to depend linearly on ionic strength, as in

${\displaystyle \log \gamma =k_{m}I\,}$

where k_m is a Sechenov coefficient.^[7]

In the example of a monobasic acid HA, assuming that the background electrolyte is the salt NaNO₃, the interaction coefficients will be for interaction between H⁺ and NO₃⁻, and between A⁻ and Na⁺.

Determination and application

Firstly, equilibrium constants are determined at a number of different ionic strengths, at a chosen temperature and particular background electrolyte. The interaction coefficients are then determined by fitting to the observed equilibrium constant values. The procedure also provides the value of K at infinite dilution. It is not limited to monobasic acids.^[8] and can also be applied to metal complexes.^[9] The SIT and Pitzer approaches have been compared recently.^[10] The Bromley equation^[11] has also been compared to both SIT and Pitzer equations.^[12]