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.

Diffusion Limit of Kinetic Equations for Multiple Species Charged Particles

Diffusive dynamics of charge regulated macro-ion solutions

The diffusive limits of two species Vlasov-Maxwell-Boltzmann equations

Global strong solutions of the Vlasov-Poisson-Boltzmann system in bounded domains

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

Charged Fluids in Porous Media

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

Global solutions to the Nernst-Planck-Euler system on bounded domain

Poisson-Boltzmann-Nernst-Planck model

Quasineutral limit of the electro-diffusion model arising in electrohydrodynamics

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

Diffusiophoresis of charge-regulated nanoparticles comprising finite ion size and electrostatic correlation effects

Charge Transport in one Dimension:Dissipative and Non-Dissipative Space-Charge Limited Currents

Homogenization of Poisson-Nernst-Planck equations for multiple species in a porous medium

On the Three-dimensional Nernst-Planck-Boussinesq System

Global weak solutions for the compressible Poisson-Nernst-Planck-Navier-Stokes System

Global Regularity for Nernst-Planck-Navier-Stokes Systems with Mixed Boundary Conditions

Interior Electroneutrality in Nernst–Planck–Navier–Stokes Systems

A Probabilistic Mean Field Limit for the Vlasov-Poisson System for Ions

Diffusion Limit with Optimal Convergence Rate of Classical Solutions to the Vlasov-Maxwell-Boltzmann System

Boltzmann Diffusive Limit with Maxwell Boundary Condition