Biophysical Journal 15: 143-162 (1975)
© 1975 the Biophysical Society

This Article Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Olivares, W Articles by McQuarrie, D A Search for Related Content PubMed PubMed Citation Articles by Olivares, W Articles by McQuarrie, D A

On the theory of ionic solutions.

ABSTRACT

One of the purposes of this paper is to assess the degree of applicability of the nonlinear Poisson-Boltzmann equation. In order to do this we compare the thermodynamic properties calculated through this equation with Monte Carlo data on 1-1 and 2-2 electrolytes described by the restricted primitive model, in which the ions are modeled by hard spheres with a coulombic potential and the solvent is modeled as a continuum dielectric medium of uniform dielectric constant epsilon. We choose Monte Carlo data rather than real experimental data since all parameters are completely specified and there is no liberty for "adjustment." Thus this serves as a definitive test. In addition, we present a simple but numerically accurate alternative approximation scheme which is not only numerically superior to the Poisson-Boltzmann equation but avoids the necessity of solving a nonlinear partial differential equation which is approximate in the first place. The new approximation scheme that is presented here is suggested by recent developments in the statistical mechanical theories of ionic solutions which are reviewed in the Introduction. Although these theories themselves yield exceedingly good comparison with experimental (Monte Carlo) data, they involve fairly advanced theoretical and mathematical techniques and do not appear to be readily solvable for other than very simple geometries. The two approximations suggested here require only the solution of the linear Debye-Hückel equation, which has been solved for a variety of systems. These two approximations are simple to apply and yield good thermodynamic properties up to concentrations of 2 M for the restricted primitive model. In addition, they have a sound theoretical foundation and are offered as a substitute for the difficult-to-solve nonlinear Poisson-Boltzmann equation.