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.

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

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

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

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

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

Self-consistent field model for strong electrostatic correlations and inhomogeneous dielectric media

Asymptotic Analysis on Dielectric Boundary Effects of Modified Poisson--Nernst--Planck Equations

Structure-Preserving and Efficient Numerical Methods for Ion Transport

Analysis of a Poisson-Nernst-Planck-Fermi system for charge transport in ion channels

Analysis of improved Nernst–Planck–Poisson models of compressible isothermal electrolytes

Poisson-Boltzmann-Nernst-Planck model

Energy Dissipative and Positivity Preserving Schemes for Large-Convection Ion Transport with Steric and Solvation Effects.

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

Mean-field Theory and Computation of Electrostatics with Ionic Concentration Dependent Dielectrics

On Applicability Of Poisson-Boltzmann Equation For Micro- And Nanoscale Electroosmotic Flows

Incorporating Dipolar Solvents with Variable Density in Poisson-Boltzmann Electrostatics

The Calculus of Boundary Variations and the Dielectric Boundary Force in the Poisson-Boltzmann Theory for Molecular Solvation

The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients

Microscopic formulation of non-local electrostatics in polar liquids embedding polarizable ions

An Efficient Finite Element Iterative Method for Solving a Nonuniform Size Modified Poisson-Boltzmann Ion Channel Model