Numerical Solution of 3D PoissonNernstPlanck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment
Abstract
We have developed efficient numerical algorithms for the solution of 3D steadystate PoissonNernstPlanck equations (PNP) with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by finite difference scheme and solved iteratively by Gummel method with relaxation. The NernstPlanck equations are transformed into Laplace equations through the Slotboom transformation. Algebraic multigrid method is then applied to efficiently solve the Poisson equation and the transformed NernstPlanck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed which reduces computational complexity from O(N2) to O(NlogN) where N is the number of grid points. Integrals involving Dirac delta function are evaluated directly by coordinate transformation which yields more accurate result compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for Li ion batteries are shown to be in good agreement with the experimental data and the results from previous studies.
 Authors:
 Publication Date:
 Research Org.:
 Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 1170083
 Report Number(s):
 PNNLSA98014
 DOE Contract Number:
 AC0576RL01830
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Communications in Computational Physics, 16(5):12981322
 Country of Publication:
 United States
 Language:
 English
Citation Formats
Meng, Da, Zheng, Bin, Lin, Guang, and Sushko, Maria L. Numerical Solution of 3D PoissonNernstPlanck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment. United States: N. p., 2014.
Web.
Meng, Da, Zheng, Bin, Lin, Guang, & Sushko, Maria L. Numerical Solution of 3D PoissonNernstPlanck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment. United States.
Meng, Da, Zheng, Bin, Lin, Guang, and Sushko, Maria L. Fri .
"Numerical Solution of 3D PoissonNernstPlanck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment". United States.
doi:.
@article{osti_1170083,
title = {Numerical Solution of 3D PoissonNernstPlanck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment},
author = {Meng, Da and Zheng, Bin and Lin, Guang and Sushko, Maria L.},
abstractNote = {We have developed efficient numerical algorithms for the solution of 3D steadystate PoissonNernstPlanck equations (PNP) with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by finite difference scheme and solved iteratively by Gummel method with relaxation. The NernstPlanck equations are transformed into Laplace equations through the Slotboom transformation. Algebraic multigrid method is then applied to efficiently solve the Poisson equation and the transformed NernstPlanck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed which reduces computational complexity from O(N2) to O(NlogN) where N is the number of grid points. Integrals involving Dirac delta function are evaluated directly by coordinate transformation which yields more accurate result compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for Li ion batteries are shown to be in good agreement with the experimental data and the results from previous studies.},
doi = {},
journal = {Communications in Computational Physics, 16(5):12981322},
number = ,
volume = ,
place = {United States},
year = {Fri Aug 29 00:00:00 EDT 2014},
month = {Fri Aug 29 00:00:00 EDT 2014}
}

A PoissonNernstPlanckFermi (PNPF) theory is developed for studying ionic transport through biological ion channels. Our goal is to deal with the finite size of particle using a Fermi like distribution without calculating the forces between the particles, because they are both expensive and tricky to compute. We include the steric effect of ions and water molecules with nonuniform sizes and interstitial voids, the correlation effect of crowded ions with different valences, and the screening effect of water molecules in an inhomogeneous aqueous electrolyte. Including the finite volume of water and the voids between particles is an important new part ofmore »

Analytical solution of the PoissonNernstPlanck equations for an electrochemical system close to electroneutrality
Single charge densities and the potential are used to describe models of electrochemical systems. These quantities can be calculated by solving a system of time dependent nonlinear coupled partial differential equations, the PoissonNernstPlanck equations. Assuming small deviations from the electroneutral equilibrium, the linearized and decoupled equations are solved for a radial symmetric geometry, which represents the interface between a cell and a sensor device. The densities and the potential are expressed by FourierBessels series. The system considered has a ratio between the Debyelength and its geometric dimension on the order of 10{sup −4} so the FourierBessel series can be approximatedmore » 
Relativistic theory of classical collisional transport in a weakly coupled plasma: Solution of Fokker{endash}Planck equation
A generalized Fokker{endash}Planck equation is developed in the pattern of a wellknown collision operator [Beliaev and Budker, Sov. Phys. Dokl. {bold 1}, 218 (1957)] and solved for a relativistic singly charged electron{endash}ion plasma under the impact of crossed electric and magnetic fields. Using the differential form of the collision operator and the modified Rosenbluth potentials appropriate for the relativistic plasma various crossfield transport coefficients are derived analytically in the limit of nonrelativistic, moderately relativistic and ultrarelativistic temperatures. {copyright} {ital 1997 American Institute of Physics.} 
Error analysis of finite element method for Poisson–Nernst–Planck equations
A priori error estimates of finite element method for timedependent PoissonNernstPlanck equations are studied in this work. We obtain the optimal error estimates in L∞(H1) and L2(H1) norms, and suboptimal error estimates in L∞(L2) norm, with linear element, and optimal error estimates in L∞(L2) norm with quadratic or higherorder element, for both semi and fully discrete finite element approximations. Numerical experiments are also given to validate the theoretical results. 
PoissonNernstPlanck equations for simulating biomolecular diffusionreaction processes I: Finite element solutions
In this paper we developed accurate finite element methods for solving 3D PoissonNernstPlanck (PNP) equations with singular permanent charges for simulating electrodiffusion in solvated biomolecular systems. The electrostatic Poisson equation was defined in the biomolecules and in the solvent, while the NernstPlanck equation was defined only in the solvent. We applied a stable regularization scheme to remove the singular component of the electrostatic potential induced by the permanent charges inside biomolecules, and formulated regular, wellposed PNP equations. An inexactNewton method was used to solve the coupled nonlinear elliptic equations for the steady problems; while an AdamsBashforthCrankNicolson method was devised formore »