{{Short description|Thermodynamic extension of Debye–Hückel theory}} '''Pitzer equations'''<ref name=pitzer>{{cite book|last1=Pitzer|first1=Kenneth S.|title=Activity coefficients in electrolyte solutions|date=1991|publisher=CRC Press|location=Boca Raton|isbn=0849354153|edition=2nd|chapter=Chapter 3: Ion Interaction Approach: Theory and Data Correlation}}</ref> are important for the understanding of the behaviour of ions dissolved in natural waters such as rivers, lakes and sea-water.<ref name=stumm_morgan>{{cite book |title=Water Chemistry |last=Stumm |first=W. |author2=Morgan, J.J. |year=1996 |publisher=Wiley |location=New York |isbn=0-471-05196-9 |url-access=registration |url=https://archive.org/details/waterchemistry00snoerich }}</ref><ref name=aquatic>{{cite book |title=Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters |last=Snoeyink |first=V.L. |author2=Jenkins, D. |year=1980 |publisher=Wiley |location=New York |isbn=0-471-51185-4}}</ref><ref>{{cite book |title=Chemical Oceanography |last=Millero |first=F.J. |edition=3rd |year=2006 |publisher=Taylor and Francis |location=London |isbn=0-8493-2280-4}}</ref> They were first described by physical chemist Kenneth Pitzer.<ref name="Connick">{{cite journal|last1=E. Connick|first1=Robert E. Connick|title=Kenneth Pitzer, 6 January 1914 · 26 December 1997|journal=Proceedings of the American Philosophical Society|date=December 2000|volume=14|issue=4|pages=479–483|jstor=1515624}}</ref> The parameters of the Pitzer equations are linear combinations of parameters, of a virial expansion of the excess Gibbs free energy, which characterise interactions amongst ions and solvent. The derivation is thermodynamically rigorous at a given level of expansion. The parameters may be derived from various experimental data such as the osmotic coefficient, mixed ion activity coefficients, and salt solubility. They can be used to calculate mixed ion activity coefficients and water activities in solutions of high ionic strength for which the Debye–Hückel theory is no longer adequate. They are more rigorous than the equations of specific ion interaction theory (SIT theory), but Pitzer parameters are more difficult to determine experimentally than SIT parameters.

== Historical development ==

A starting point for the development can be taken as the virial equation of state for a gas. :<math> PV= R T + B P + C P^2 + D P^3 \dots </math> where <math>P</math> is the pressure, <math>V</math> is the volume, <math>T</math> is the temperature and <math>B, C, D</math> ... are known as virial coefficients. The first term on the right-hand side is for an ideal gas. The remaining terms quantify the departure from the ideal gas law with changing pressure, <math>P</math>. It can be shown by statistical mechanics that the second virial coefficient arises from the intermolecular forces between ''pairs'' of molecules, the third virial coefficient involves interactions between three molecules, etc. This theory was developed by McMillan and Mayer.<ref>{{cite journal|last=McMillan|first=W.G.|author2=Mayer, J.E. |year=1945|title=The statistical thermodynamics of multicomponent systems |journal =J. Chem. Phys.|volume=13|pages=276|doi=10.1063/1.1724036|issue=7 |bibcode = 1945JChPh..13..276M }}</ref>

Solutions of uncharged molecules can be treated by a modification of the McMillan-Mayer theory. However, when a solution contains electrolytes, electrostatic interactions must also be taken into account. The Debye–Hückel theory<ref>{{cite journal|last=Debye|first=P.|author2=Hückel, E.|year=1923|title=Zur Theorie der Electrolyte|journal=Phys. Z.|volume=24|pages=185}}</ref> was based on the assumption that each ion was surrounded by a spherical "cloud" or ionic atmosphere made up of ions of the opposite charge. Expressions were derived for the variation of single-ion activity coefficients as a function of ionic strength. This theory was very successful for dilute solutions of 1:1 electrolytes and, as discussed below, the Debye–Hückel expressions are still valid at sufficiently low concentrations. The values calculated with Debye–Hückel theory diverge more and more from observed values as the concentrations and/or ionic charges increases. Moreover, Debye–Hückel theory takes no account of the specific properties of ions such as size or shape.

Brønsted had independently proposed an empirical equation,<ref>{{cite journal| last=Brønsted |first=J.N. |year=1922 |title=Studies on solubility IV. The principle of the specific interaction of ions|journal=J. Am. Chem. Soc.|volume=44|pages=877–898|doi=10.1021/ja01426a001| issue=5|url=https://zenodo.org/record/1428790 }}</ref> :<math> \ln{\gamma} = - \alpha m^{1/2} - 2 \beta m </math> :<math> 1-\varphi = (\alpha/3) m^{1/2} + \beta m </math> in which the activity coefficient depended not only on ionic strength, but also on the concentration, ''m'', of the specific ion through the parameter ''β''. This is the basis of SIT theory. It was further developed by Guggenheim.<ref name=G>{{cite journal|last=Guggenheim|first=E.A.|author2=Turgeon, J.C. |year=1955|title=Specific interaction of ions|journal=Trans. Faraday Soc. |volume=51|pages=747–761|doi=10.1039/TF9555100747}}</ref> Scatchard<ref>{{cite journal| last=Scatchard |first=G. |year=1936 |title=Concentrated solutions of strong electrolytes|journal=Chem. Rev.|volume=19|pages=309–327|doi=10.1021/cr60064a008| issue=3}}</ref> extended the theory to allow the interaction coefficients to vary with ionic strength. Note that the second form of Brønsted's equation is an expression for the osmotic coefficient. Measurement of osmotic coefficients provides one means for determining mean activity coefficients.

== The Pitzer parameters ==

The exposition begins with a virial expansion of the excess Gibbs free energy<ref name=pitzer/>{{rp|p=84}} :<math>\frac{G^{ex}}{W_wRT} = f(I) +\sum_i \sum_j b_ib_j\lambda_{ij}(I)+\sum_i \sum_j \sum_kb_ib_jb_k\mu_{ijk}+\cdots</math> ''W<sub>w</sub>'' is the mass of the water in kilograms,'' b<sub>i</sub>, b<sub>j</sub>'' ... are the molalities of the ions and <math display="inline">I = \tfrac{1}{2} \sum_i b_i {z_i}^2 </math> is the molal ionic strength. The first term, ''f(I)'' represents a Debye–Hückel extended law (see below). The quantities ''λ<sub>ij</sub>(I)'' represent the short-range interactions in the presence of solvent between solute particles ''i'' and ''j''. This binary interaction parameter or second virial coefficient depends on ionic strength, on the particular species ''i'' and ''j'' and the temperature and pressure. The quantities ''μ''<sub>''ijk''</sub> represent the interactions between three particles. Higher terms may also be included in the virial expansion.

Next, the free energy is expressed as the sum of chemical potentials, or partial molal free energy, :<math>G= \sum_i \mu_i\cdot N_i = \sum_i \left (\mu^0_i +RT \ln b_i\gamma_i \right )\cdot N_i</math> and an expression for the activity coefficient is obtained by differentiating the virial expansion with respect to a molality b. :<math>\ln \gamma_i = \frac{\partial(\frac{G^{ex}}{W_wRT})}{\partial b_i} =\frac{z_i^2}{2}f' +2\sum_j \lambda_{ij}b_j +\frac{z_i^2}{2}\sum_j\sum_k \lambda'_{jk} b_jb_k + 3\sum_j\sum_k \mu_{ijk} b_jb_k+ \cdots </math> And molal osmotic coefficient: :<math>\phi-1=\left(\sum_ib_i\right)^{-1}\left[If'-f + \sum_i\sum_j\left(\lambda_{ij}+I\lambda'_{ij} \right)b_ib_j +2\sum_i\sum_j\sum_k \mu_{ijk} b_ib_jb_k + \cdots\right]</math> However, these forms are not used directly, because it is not possible to determine the coefficients independently due to charge neutrality constraints.

Instead the sums are re-worked in terms of observable (charge-neutral) combinations and some choices are made about their functional forms, which we will see below:<ref name=pitzer/> * The <math>\lambda_{ij}</math>-related terms for salt pairs are gathered together into observable combinations (<math>B_{ca}</math> values) and then a specific exponential functional form is imposed on them in terms of <math>\beta</math> parameters, see below. * The <math>\mu_{ijk}</math>-related terms for salt pairs are likewise gathered into observable combinations (<math>C_{ca}</math> values). * All ternary interactions <math>\mu_{ijk}</math> involving three ions of same sign are set to 0. * For mixed electrolytes, new combinations appear: <math>\Phi_{cc'}</math>, <math>\Phi_{aa'}</math>, and ternary interactions <math>\psi_{cc'a}</math>, <math>\psi_{caa'}</math>. A special functional form is applied to <math>\Phi</math> in cases of unsymmetrical mixing (two ions present with same sign of charge but different magnitude). * For electrolytes combined with neutral solutes, <math>\lambda_{ij}</math> remain as-is, but they are assumed to be independent of ionic strength.

These observable combinations then provide a set of free parameters which are then empirically fit to experimental data.

=== Pure electrolyte case ===

Consider a simple electrolyte ''M''<sub>''p''</sub>''X''<sub>''q''</sub> with molal concentration ''m'', dissolved to ions ''M''<sup>''z''<sup>+</sup></sup> and ''X''<sup>''z''<sup>−</sup></sup>, with ionic molal concentrations <math>b_M = pm</math> and <math>b_X = qm</math>.

The Pitzer parameters <math>f^\phi</math>, <math>B^\phi_{MX}</math> and <math>C^\phi_{MX}</math> are defined as :<math>f^\phi=\frac{f'-\frac{f}{I}}{2}</math> :<math>B^\phi_{MX}=\lambda_{MX}+I\lambda'_{MX} +\left(\frac{p}{2q}\right)\left(\lambda_{MM}+I\lambda'_{MM}\right)+\left(\frac{q}{2p}\right)\left(\lambda_{XX}+I\lambda'_{XX}\right)</math> :<math>C^\phi_{MX} =\left[\frac{3}{\sqrt{pq}}\right] \left(p\mu_{MMX}+q\mu_{MXX}\right). </math> (Terms involving <math>\mu_{MMM}</math> and <math>\mu_{XXX}</math> are not included as interactions between three ions of the same charge are unlikely to occur except in very concentrated solutions.)

With these definitions, the expression for the (molal-basis) osmotic coefficient becomes<ref name=pitzer/>{{rp|p=87}} :<math>\phi-1=|z^+z^-|f^\phi+m\left(\frac{2pq}{p+q}\right)B^\phi_{MX} +m^2\left[2\frac{(pq)^{3/2}}{p+q}\right]C^\phi_{MX}. </math> A similar expression is obtained for the (molal-basis) mean activity coefficient:<ref name=pitzer/>{{rp|p=88}} :<math>\ln \gamma_\pm =|z^+z^-|f^\gamma+m\left(\frac{2pq}{p+q}\right)B^\gamma_{MX} +m^2\left[2\frac{(pq)^{3/2}}{p+q}\right]C^\gamma_{MX}, </math> where <math>f^\gamma</math>, <math>B^\gamma_{MX}</math> and <math>C^\gamma_{MX}</math> are related to <math>f^\phi</math>, <math>B^\phi_{MX}</math> and <math>C^\phi_{MX}</math>, but distinct.

Finally, some forms are imposed on the coefficients based on a mixture of theoretical and empirical observations: <ul> <li> The term ''f''<sup>''φ''</sup> is defined to be an extended Debye–Hückel term:<ref name=pitzer/>{{rp|p=87}} :<math>f^\phi = -A_\phi \frac{I^{1/2}}{1 + b I^{1/2}}</math> with <math>A_{\phi}</math> being calculated in terms of the solvent dielectric constant, and <math>b=1.2~\mathrm{kg}^{1/2}\mathrm{mol}^{-1/2}</math> is defined as a universal empirical parameter (note this <math>b</math> should not be confused with the molality <math>b_i</math>). </li> <li> The ''B'' parameter was found empirically to show an ionic strength dependence which could be expressed as :<math>B^\phi_{MX}=\beta^{(0)}_{MX} + \beta^{(1)}_{MX} e^{-\alpha \sqrt I},</math> or sometimes with a second term which can often capture ion pairing effects without requiring explicit ion association accounting:<ref name=pitzer/>{{rp|pp=87,93–95}} :<math>B^\phi_{MX}=\beta^{(0)}_{MX} + \beta^{(1)}_{MX} e^{-\alpha_1 \sqrt I} + \beta^{(2)}_{MX} e^{-\alpha_2 \sqrt I}.</math> (with specific values of <math>\alpha_1</math> and <math>\alpha_2</math> being chosen depending on the ion charges). The empirical Pitzer data tables therefore list <math>\beta^{(0)}_{MX}</math>, <math>\beta^{(1)}_{MX}</math>, and sometimes <math>\beta^{(2)}_{MX}</math> (usually 0), whereas <math>C^\phi_{MX}</math> is directly tabulated. </li> </ul>

''Note on ion association:'' If ion pairing is included as an ion association equilibrium with an explicit separate solute species (with its own separately-accounted concentration), then the empirical values <math>\beta^{(n)}</math> and <math>C^\phi</math> will change completely. Moreover, this choice fundamentally redefines the meaning and numerical values of ionic molalities, ionic strength, mean activity coefficients, and even the osmotic coefficient. Ionic mean activities and solvent activity are, however, thermodynamically independent of this accounting choice.

=== General case: mixed electrolytes, neutral solutes, and single-ion activities ===

Pitzer defines the above pure electrolyte case (yielding osmotic coefficient and mean activity) to be mathematically simple, then 'works backwards' to deduce the general case in a way that is consistent with the pure electrolyte case. Pitzer thus arrives at the following Gibbs energy:<ref name=pitzer/>{{rp|p=88}}

<math display="block"> \begin{aligned} \frac{G^{ex}}{W_w RT} &= f(I) \\ &\quad + 2 \sum_{c} \sum_{a} b_c b_a \left[ B_{ca} + \left( \sum_{c} b_c z_c \right) C_{ca} \right] \\ &\quad + \mathop{\sum\sum}_{c < c'} b_c b_{c'} \left[ 2\Phi_{cc'} + \sum_{a} b_a \psi_{cc'a} \right] \\ &\quad + \mathop{\sum\sum}_{a < a'} b_a b_{a'} \left[ 2\Phi_{aa'} + \sum_{c} b_c \psi_{caa'} \right] \\ &\quad + 2 \sum_{n} \sum_{c} b_n b_c \lambda_{nc} + 2 \sum_{n} \sum_{a} b_n b_a \lambda_{na} \\ &\quad + 2 \mathop{\sum\sum}_{n < n'} b_n b_{n'} \lambda_{nn'} + \sum_{n} {b_n}^2 \lambda_{nn} + \dots. \end{aligned} </math>

where <math>c</math> indices are positive ions (cations), <math>a</math> are negative ions (anions), and <math>n</math> are neutral solutes. Note that <math>B_{ca}</math>, <math>\Phi_{cc'}</math>, <math>\Phi_{aa'}</math> are also functions of ionic strength.

This form (with the <math>\dots</math> truncated) is then the actual master thermodynamic equation underlying Pitzer theory, and by differentiation it yields all other observable quantities (activity coefficients, osmotic coefficients).<ref name=pitzer/>{{rp|p=89}}

This expression omits various terms proportional to the total charge <math display="inline">\sum_{i} z_i b_i</math>, and therefore it yields different single-ion activities compared to the original <math>G^{ex}</math>, but only in an unobservable way.<ref name=pitzer/>{{rp|p=91}} The raw Pitzer single-ion activities are not experimentally observable on their own, and in practice they are combined into observables (like the mean activity above), or they are at least transformed to obey well-known activity conventions (such as the MacInnes convention).

The various terms are all consistent with the pure electrolyte case. For example, the full Debye-Huckel term is chosen to be<ref name=pitzer/>{{rp|pp=88,120–122}} <math display="block">f(I) = -(4 I A_{\phi}/b)\ln(1 + b I^{1/2}) </math> which is consistent with the <math>f^\phi</math> above.

Expressions for the interaction coefficients <math>B_{ca}</math>, <math>C_{ca}</math>, <math>\Phi</math>, and <math>\psi</math> can be found in standard references.<ref name=pitzer/>{{rp|pp=87–88}} It is crucial to note that the parameters appearing in the Gibbs energy differ from the tabulated parameters often labeled with <math>\phi</math> or <math>\gamma</math> superscripts. For example, <math>B_{ca}</math> is the fundamental interaction parameter, while <math>B_{ca}^\phi</math> and <math>B_{ca}^\gamma</math> are derivatives used for osmotic and activity coefficient calculations respectively.

== Commentary ==

These equations were applied to an extensive range of experimental data at 25&nbsp;°C with excellent agreement to about 6&nbsp;mol&nbsp;kg<sup>−1</sup> for various types of electrolyte.<ref>{{cite journal|last=Pitzer|first=K.S.|author2=Mayorga, G. |year=1973|title=Thermodynamics of Electrolytes, II. Activity and osmotic coefficients with one or both ions univalent|journal=J. Phys. Chem.|volume=77|pages=2300–2308|doi=10.1021/j100638a009|issue=19|url=https://escholarship.org/uc/item/1rc6r4vm}}</ref><ref>{{cite journal|last=Pitzer|first=K.S.|author2=Mayorga, G. |year=1974|title=Thermodynamics of Electrolytes. III. Activity and osmotic coefficients for 2–2 electrolytes |journal=J. Solution. Chem.|volume=3|pages=539–546|doi=10.1007/BF00648138|issue=7|url=http://www.escholarship.org/uc/item/0q16265x}}</ref> The treatment can be extended to mixed electrolytes<ref>{{cite journal|last=Pitzer|first=K.S.|author2=Kim, J.J. |year=1974|title=Thermodynamics of electrolytes. IV. Activity and osmotic coefficients for mixed electrolytes|journal=J. Am. Chem. Soc.|volume=96|issue=18|pages=5701–5707|doi=10.1021/ja00825a004|url=http://www.escholarship.org/uc/item/7rq7244j}}</ref> and to include association equilibria.<ref name=pitzer/>{{rp|p=93}} Values for the parameters ''&beta;''<sup>(0)</sup>, ''&beta;''<sup>(1)</sup> and ''C'' for inorganic and organic acids, bases and salts have been tabulated.<ref>Pitzer (1991), Tables 2-11</ref> Temperature and pressure variation is also discussed.

One area of application of Pitzer parameters is to describe the ionic strength variation of equilibrium constants measured as concentration quotients. Both SIT and Pitzer parameters have been used in this context, For example, both sets of parameters were calculated for some uranium complexes and were found to account equally well for the ionic strength dependence of the stability constants.<ref>{{cite journal|last=Crea|first=F.|author2=Foti, C. |author3=Sammartano, S. |year=2008|title=Sequestering ability of polycarboxylic acids towards dioxouranium(V)|journal=Talanta|volume=28|pages=775–778|doi=10.1016/j.talanta.2007.12.009|issue=3}}</ref>

Pitzer parameters and SIT theory have been extensively compared. There are more parameters in the Pitzer equations than in the SIT equations. Because of this the Pitzer equations provide for more precise modelling of mean activity coefficient data and equilibrium constants. However, the determination of the greater number of Pitzer parameters means that they are more difficult to determine.<ref name=grenthe>{{cite book|last=Grenthe|first=I.|author2=Puigdomenech, I. |title=Modelling in aquatic chemistry|publisher=Nuclear Energy Agency, O.E.C.D|year=1997|isbn=92-64-15569-4}} Chapter 9, Estimation of medium effects on thermodynamic data</ref>

== Compilation of Pitzer parameters ==

Besides the set of parameters obtained by Pitzer et al. in the 1970s mentioned in the previous section. Kim and Frederick<ref>{{cite journal|last1=Kim|first1=Hee Taik|last2=Frederick|first2=William J.|title=Evaluation of Pitzer ion interaction parameters of aqueous electrolytes at 25.degree.C. 1. Single salt parameters|journal=Journal of Chemical & Engineering Data|date=April 1988|volume=33|issue=2|pages=177–184|doi=10.1021/je00052a035}}</ref><ref>{{cite journal|last1=Kim|first1=Hee Taik|last2=Frederick|first2=William J.|title=Evaluation of Pitzer ion interaction parameters of aqueous mixed electrolyte solutions at 25.degree.C. 2. Ternary mixing parameters|journal=Journal of Chemical & Engineering Data|date=July 1988|volume=33|issue=3|pages=278–283|doi=10.1021/je00053a017}}</ref> published the Pitzer parameters for 304 single salts in aqueous solutions at 298.15 K, extended the model to the concentration range up to the saturation point. Those parameters are widely used, however, many complex electrolytes including ones with organic anions or cations, which are very significant in some related fields, were not summarized in their paper.

For some complex electrolytes, Ge et al.<ref>{{cite journal|last1=Ge|first1=Xinlei|last2=Zhang|first2=Mei|last3=Guo|first3=Min|last4=Wang|first4=Xidong|title=Correlation and Prediction of Thermodynamic Properties of Some Complex Aqueous Electrolytes by the Modified Three-Characteristic-Parameter Correlation Model|journal=Journal of Chemical & Engineering Data|date=April 2008|volume=53|issue=4|pages=950–958|doi=10.1021/je7006499}}</ref> obtained the new set of Pitzer parameters using up-to-date measured or critically reviewed osmotic coefficient or activity coefficient data.

== Comparable activity coefficient models == Besides the well-known Pitzer-like equations, there is a simple and easy-to-use semi-empirical model, which is called the three-characteristic-parameter correlation (TCPC) model. It was first proposed by Lin et al.<ref>{{cite journal|last1=Lin|first1=Cheng-Long|last2=Lee|first2=Liang-Sun|last3=Tseng|first3=Hsieng-Cheng|title=Thermodynamic behavior of electrolyte solutions|journal=Fluid Phase Equilibria|date=September 1993|volume=90|issue=1|pages=57–79|doi=10.1016/0378-3812(93)85004-6}}</ref> It is a combination of the Pitzer long-range interaction and short-range solvation effect:

:ln ''&gamma;''&nbsp;= ln ''&gamma;''<sup>PDH</sup> + ln ''&gamma;''<sup>SV</sup>

Ge et al.<ref>{{cite journal|last1=Ge|first1=Xinlei|last2=Wang|first2=Xidong|last3=Zhang|first3=Mei|last4=Seetharaman|first4=Seshadri|title=Correlation and Prediction of Activity and Osmotic Coefficients of Aqueous Electrolytes at 298.15 K by the Modified TCPC Model|journal=Journal of Chemical & Engineering Data|date=March 2007|volume=52|issue=2|pages=538–547|doi=10.1021/je060451k}}</ref> modified this model, and obtained the TCPC parameters for a larger number of single salt aqueous solutions. This model was also extended for a number of electrolytes dissolved in methanol, ethanol, 2-propanol, and so on.<ref>{{cite journal|last1=Ge|first1=Xinlei|last2=Zhang|first2=Mei|last3=Guo|first3=Min|last4=Wang|first4=Xidong|title=Correlation and Prediction of Thermodynamic Properties of Nonaqueous Electrolytes by the Modified TCPC Model|journal=Journal of Chemical & Engineering Data|date=January 2008|volume=53|issue=1|pages=149–159|doi=10.1021/je700446q}}</ref> Temperature dependent parameters for a number of common single salts were also compiled, available at.<ref>{{cite journal|last1=Ge|first1=Xinlei|last2=Wang|first2=Xidong|title=A Simple Two-Parameter Correlation Model for Aqueous Electrolyte Solutions across a Wide Range of Temperatures|journal=Journal of Chemical & Engineering Data|date=12 February 2009|volume=54|issue=2|pages=179–186|doi=10.1021/je800483q}}</ref>

The performance of the TCPC model in correlation with the measured activity coefficient or osmotic coefficients is found to be comparable with Pitzer-like models.

Due to its empirical aspects, the Pitzer modelling framework has a number of well-known limitations.<ref>{{Cite journal |last=Rowland |first=Darren |last2=Königsberger |first2=E. |last3=Hefter |first3=G. |last4=May |first4=P.M. |date=2015 |title=Aqueous electrolyte solution modelling: some limitations of the Pitzer equations |url=https://www.sciencedirect.com/science/article/abs/pii/S0883292714002327 |journal=Applied Geochemistry |volume=55 |pages=170–183 |doi=10.1016/j.apgeochem.2014.09.021 |via=|url-access=subscription }}</ref> Most importantly, to improve the fits to experimental data, different variations of the equations have been described. Extrapolations, especially in the temperature and pressure domain, are generally problematic. One alternative modelling approach<ref>{{Cite journal |last=May |first=Peter M. |last2=May |first2=Eric |date=2024 |title=Ion Trios: Cause of Ion Specific Interactions in Aqueous Solutions and Path to a Better pH Definition |url=https://pubs.acs.org/doi/10.1021/acsomega.4c07525 |journal=ACS Omega |volume=9 |issue=46 |pages=46373–46386 |doi=10.1021/acsomega.4c07525 |doi-access=free|pmc=11579776 }}</ref> has been specifically designed to address this extrapolation issue by reducing the number of equation parameters while maintaining similar predictive precision and accuracy.

== See also == * Bromley equation * Davies equation * Osmotic coefficient

== References == {{Reflist}}

Category:Thermodynamic equations Category:Chemical thermodynamics Category:Equilibrium chemistry Category:Electrochemical equations