Abstract:The Poisson-Nernst-Planck (PNP) model is based on a mean-field approximation of ion interactions and continuum descriptions of concentration and electrostatic potential. It provides qualitative explanation and increasingly quantitative predictions of experimental measurements for the ion transport problems in many areas such as semiconductor devices, nanofluidic systems, and biological systems, despite many limitations. While the PNP model gives a good prediction of the ion transport phenomenon for chemical, physical, and biological systems, the number of equations to be solved and the number of diffusion coefficient profiles to be determined for the calculation directly depend on the number of ion species in the system, since each ion species corresponds to one Nernst-Planck equation and one position-dependent diffusion coefficient profile. In a complex system with multiple ion species, the PNP can be computationally expensive and parameter demanding, as experimental measurements of diffusion coefficient profiles are generally quite limited for most confined regions such as ion channels, nanostructures and nanopores. We propose an alternative model to reduce number of Nernst-Planck equations to be solved in complex chemical and biological systems with multiple ion species by substituting Nernst-Planck equations with Boltzmann distributions of ion concentrations. As such, we solve the coupled Poisson-Boltzmann and Nernst-Planck (PBNP) equations, instead of the PNP equations. The proposed PBNP equations are derived from a total energy functional by using the variational principle. We design a number of computational techniques, including the Dirichlet to Neumann mapping, the matched interface and boundary, and relaxation based iterative procedure, to ensure efficient solution of the proposed PBNP equations. Two protein molecules, cytochrome c551 and Gramicidin A, are employed to validate the proposed model under a wide range of bulk ion concentrations and external voltages. Extensive numerical experiments show that there is an excellent consistency between the results predicted from the present PBNP model and those obtained from the PNP model in terms of the electrostatic potentials, ion concentration profiles, and current-voltage (I-V) curves. The present PBNP model is further validated by a comparison with experimental measurements of I-V curves under various ion bulk concentrations. Numerical experiments indicate that the proposed PBNP model is more efficient than the original PNP model in terms of simulation time.

A local approximation of fundamental measure theory incorporated into three dimensional Poisson-Nernst-Planck equations to account for hard sphere repulsion among ions

Poisson-Nernst-Planck Systems for Ion Flow with a Local Hard-Sphere Potential for Ion Size Effects.

A Modified Poisson--Nernst--Planck Model with Excluded Volume Effect: Theory and Numerical Implementation

Ion Size Effects on Individual Fluxes Via Poisson-Nernst-Planck Systems with Bikerman's Local Hard-Sphere Potential: Analysis Without Electroneutrality Boundary Conditions

Numerical Methods for a Poisson-Nernst-Planck-Fermi Model

Poisson-Boltzmann-Nernst-Planck model

Finite ion size effects on I-V relations via Poisson-Nernst-Planck systems with two cations: A case study

Ion Transport in Dipolar Medium I: A Local Dielectric Poisson–Nernst–Planck/Poisson–Boltzmann Model

Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

Poisson-Nernst-Planck-Fermi Theory for Ion Channels

A Fundamental-Measure Theory for Inhomogeneous Associating Fluids

Molecular Mean-Field Theory of Ionic Solutions: a Poisson-Nernst-Planck-Bikerman Model

Scaled Particle Theory for Multicomponent Hard Sphere Fluids Confined in Random Porous Media

A Novel Method to Describe the Interaction Pressure Between Charged Plates with Application of the Weighted Correlation Approach

Poisson-Boltzmann theory with non-linear ion correlations

Density-functional Theory of Spherical Electric Double Layers and Zeta Potentials of Colloidal Particles in Restricted-Primitive-model Electrolyte Solutions.

Fundamental measure theory for non-spherical hard particles: predicting liquid crystal properties from the particle shape

Finite-Difference Approximations and Local Algorithm for the Poisson and Poisson-Boltzmann Electrostatics

Finite element method for the numerical simulation of modified Poisson-Nernst-Planck/Navier-Stokes model

Mean-Field Description Of Ionic Size Effects With Nonuniform Ionic Sizes: A Numerical Approach

Studies on Ionic Flows via Poisson–Nernst–Planck Systems with Bikerman's Local Hard-Sphere Potentials under Relaxed Neutral Boundary Conditions