title: Poisson-Boltzmann-Nernst-Planck model.
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.
authors: Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, USA.
authors: Qiong
authors: Q
authors: Zheng
authors: Guo-Wei
authors: GW
authors: Wei
references: Bull Math Biol. 2010 Aug;72(6):1562-622
references: 20169418
references: J Chem Phys. 2007 Sep 21;127(11):114106
references: 17887827
references: J Comput Phys. 2010 Jun 20;229(12):4431-4460
references: 20396650
references: J Chem Phys. 2009 Dec 7;131(21):215101
references: 19968368
references: Nat Nanotechnol. 2010 Apr;5(4):297-301
references: 20305644
references: Biophys J. 2003 Jun;84(6):3594-606
references: 12770869
references: Commun Comput Phys. 2013 Jan 1;13(1):285-324
references: 23550030
references: Math Biosci. 1992 Dec;112(2):189-206
references: 1283350
references: J Comput Chem. 2011 Mar;32(4):756-70
references: 20845420
references: J Chem Phys. 2007 Jun 28;126(24):244108
references: 17614538
references: Biopolymers. 1996 Mar;38(3):305-20
references: 8906967
references: Nucleic Acids Res. 2004 Jul 1;32(Web Server issue):W665-7
references: 15215472
references: Nucleic Acids Res. 2007 Jul;35(Web Server issue):W522-5
references: 17488841
references: J Comput Phys. 2010 Nov 1;229(22):8231-8258
references: 20938489
references: J Gen Physiol. 1999 Jun;113(6):789-94
references: 10352030
references: J Colloid Interface Sci. 2007 Nov 15;315(2):731-9
references: 17681522
references: Annu Rev Biophys Biophys Chem. 1987;16:205-26
references: 2439095
references: J Phys Chem A. 2007 Dec 13;111(49):12506-12
references: 17960920
references: Biophys J. 2000 Jul;79(1):80-93
references: 10866939
references: J Phys Chem B. 1998 Apr 30;102(18):3586-616
references: 24889800
references: J Comput Phys. 2011 Jun;230(13):5239-5262
references: 21552336
references: Biochim Biophys Acta. 2002 Oct 11;1565(2):267-86
references: 12409200
references: Lab Chip. 2008 Jan;8(1):23-33
references: 18094759
references: Biophys J. 1998 Dec;75(6):2830-44
references: 9826605
references: Phys Rev E Stat Nonlin Soft Matter Phys. 2007 Feb;75(2 Pt 1):021503
references: 17358344
references: Int J Heat Mass Transf. 2009 Jul;52(15-16):3504-3509
references: 20585410
references: J Comput Chem. 1996 Aug;17(11):1344-51
references: 25400153
references: J Math Biol. 2011 Dec;63(6):1139-200
references: 21279359
references: Phys Rev E Stat Nonlin Soft Matter Phys. 2007 Feb;75(2 Pt 1):021502
references: 17358343
references: Q Rev Biophys. 2001 Nov;34(4):473-561
references: 11852594
references: Q Rev Biophys. 2004 Feb;37(1):15-103
references: 17390604
references: Physica A. 2010 Apr 1;389(7):1329-1345
references: 20228879
references: Annu Rev Biophys Biophys Chem. 1986;15:29-57
references: 2424476