Solution of the nonlinear Poisson-Boltzmann equation using pseudo-transient continuation and the finite element method.

The nonlinear Poisson-Boltzmann (PB) equation is solved using Newton-Krylov iterations coupled with pseudo-transient continuation. The PB potential is used to compute the electrostatic energy and evaluate the force on a user-specified contour. The PB solver is embedded in a existing, 3D, massively parallel, unstructured-grid, finite element code. Either Dirichlet or mixed boundary conditions are allowed. The latter specifies surface charges, approximates far-field conditions, or linearizes conditions "regulating" the surface charge. Stability and robustness are proved using results for backward Euler differencing of diffusion equations. Potentials and energies of charged spheres and plates are computed and results compared to analysis. An approximation to the potential of the nonlinear, spherical charge is derived by combining two analytic formulae. The potential and force due to a conical probe interacting with a flat plate are computed for two types of boundary conditions: constant potential and constant charge. The second case is compared with direct force measurements by chemical force microscopy. The problem is highly nonlinear-surface potentials of the linear and nonlinear PB equations differ by over an order of magnitude. Comparison of the simulated and experimentally measured forces shows that approximately half of the surface carboxylic acid groups, of density 1/(0.2 nm2), ionize in the electrolyte implying surface charges of 0.4 C/m2, surface potentials of 0.27 V, and a force of 0.6 nN when the probe and plate are 8.7 nm apart.

[1]  C. Kelley,et al.  Convergence Analysis of Pseudo-Transient Continuation , 1998 .

[2]  A. Lenhoff,et al.  Energetic and Entropic Contributions to the Interaction of Unequal Spherical Double Layers , 1994 .

[3]  Barry Honig,et al.  Calculating total electrostatic energies with the nonlinear Poisson-Boltzmann equation , 1990 .

[4]  V. Adrian Parsegian,et al.  DNA‐Inspired Electrostatics , 2000 .

[5]  Charles M. Lieber,et al.  DIRECT GROWTH OF SINGLE-WALLED CARBON NANOTUBE SCANNING PROBE MICROSCOPY TIPS , 1999 .

[6]  B. Montgomery Pettitt,et al.  Numerical Considerations in the Computation of the Electrostatic Free Energy of Interaction within the Poisson-Boltzmann Theory , 1997 .

[7]  Douglas A. Lauffenburger,et al.  NUMERICAL SOLUTION OF THE NONLINEAR POISSON-BOLTZMANN EQUATION FOR A MEMBRANE-ELECTROLYTE SYSTEM , 1994 .

[8]  M. Lozada-Cassou,et al.  The application of the hypernetted chain approximation to the electrical double layer: Comparison with Monte Carlo results for symmetric salts , 1982 .

[9]  Bowen,et al.  Adaptive Finite-Element Solution of the Nonlinear Poisson-Boltzmann Equation: A Charged Spherical Particle at Various Distances from a Charged Cylindrical Pore in a Charged Planar Surface , 1997, Journal of colloid and interface science.

[10]  Richard A. Friesner,et al.  Numerical solution of the Poisson–Boltzmann equation using tetrahedral finite‐element meshes , 1997 .

[11]  L. B. Bhuiyan,et al.  An improved modified Poisson–Boltzmann equation in electric-double-layer theory , 1983 .

[12]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[13]  A. I. Shestakov,et al.  Combining Cell- and Point-Centered Methods in 3D, Unstructured-Grid Radiation–Hydrodynamic Codes , 2001 .

[14]  Carlo Caccamo,et al.  An improved closure for the Born–Green–Yvon equation for the electric double layer , 1986 .

[15]  D. Chapman,et al.  LI. A contribution to the theory of electrocapillarity , 1913 .

[16]  J. Andrew McCammon,et al.  Computation of electrostatic forces on solvated molecules using the Poisson-Boltzmann equation , 1993 .

[17]  Huan-Xiang Zhou,et al.  Macromolecular electrostatic energy within the nonlinear Poisson–Boltzmann equation , 1994 .

[18]  G. Gouy.  Sur la fonction électrocapillaire , 1916 .

[19]  Derek Y. C. Chan,et al.  Computation of forces between spherical colloidal particles : nonlinear Poisson-Boltzmann theory , 1994 .

[20]  G. Torrie,et al.  Electrical double layers. II. Monte Carlo and HNC studies of image effects , 1982 .

[21]  Derek Y. C. Chan,et al.  Calculations of Electric Double-Layer Force and Interaction Free Energy between Dissimilar Surfaces , 1995 .

[22]  C. Radke,et al.  Double layer interactions between charge-regulated colloidal surfaces: pair potentials for spherical particles bearing ionogenic surface groups , 1993 .

[23]  J. Kirkwood,et al.  On the Theory of Strong Electrolyte Solutions , 1934 .

[24]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[25]  J. Bechhoefer,et al.  Calibration of atomic‐force microscope tips , 1993 .

[26]  A. I. Shestakov,et al.  The radiation-hydrodynamic ICF3D code , 2000 .

[27]  D. Henderson,et al.  Mixtures of hard ions and dipoles against a charged wall: The Ornstein–Zernike equation, some exact results, and the mean spherical approximation , 1981 .

[28]  Charles M. Lieber,et al.  Force Titrations and Ionization State Sensitive Imaging of Functional Groups in Aqueous Solutions by Chemical Force Microscopy , 1997 .

[29]  J. P. Valleau,et al.  Electrical double layers. I. Monte Carlo study of a uniformly charged surface , 1980 .

[30]  M. Plischke,et al.  Pair correlation functions and density profiles in the primitive model of the electric double layer , 1988 .

[31]  D. Chan,et al.  The structure of electrolytes at charged surfaces: Ion–dipole mixtures , 1980 .

[32]  H. T. Davis,et al.  A nonlocal free-energy density-functional approximation for the electrical double layer , 1990 .

[33]  M. Fujihira,et al.  Effect of tip shape on force-distance curves for AFM in aqueous electrolytes , 1994 .

[34]  Thomas W. Healy,et al.  Electrical double layer interactions under regulation by surface ionization equilibria–dissimilar amphoteric surfaces , 1976 .

[35]  Malcolm E. Davis,et al.  Electrostatics in biomolecular structure and dynamics , 1990 .

[36]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[37]  J. P. Valleau,et al.  Electrical double layers. 4. Limitations of the Gouy-Chapman theory , 1982 .

[38]  S. Levine,et al.  Theory of the electric double layer using a modified poisson–boltzman equation , 1980 .

[39]  J. Theodoor,et al.  The role of energy and entropy in the electrical double layer , 1990 .