Discrete and Continuous Dynamical Systems - Series B (DCDS-B)

Energy variational approach to study charge inversion (layering) near charged walls

Pages: 2725 - 2743, Volume 17, Issue 8, November 2012      doi:10.3934/dcdsb.2012.17.2725

       Abstract        References        Full Text (598.9K)       Related Articles

YunKyong Hyon - Department of Mechanical Engineering, University of Nevada, Reno, Reno, NV 89557, United States (email)
James E. Fonseca - Department of Molecular Biophysics & Physiology Rush Medical Center, 1653 West Congress, Parkway, Chicago, IL 60612, United States (email)
Bob Eisenberg - Department of Molecular Biophysics & Physiology Rush Medical Center, 1653 West Congress, Parkway, Chicago, IL 60612, United States (email)
Chun Liu - Department of Mathematics and Center for Materials Physics, Penn State University, University Park, PA 16802, United States (email)

Abstract: We introduce a mathematical model, which describes the charge inversion phenomena in systems with a charged wall or boundary. This model may prove helpful in understanding semiconductor devices, ion channels, and electrochemical systems like batteries that depend on complex distributions of charge for their function. The mathematical model is derived using the energy variational approach that takes into account ion diffusion, electrostatics, finite size effects, and specific boundary behavior. In ion dynamic theory, a well-known system of equations is the Poisson-Nernst-Planck (PNP) equation that includes entropic and electrostatic energy. The PNP type of equation can also be derived by the energy variational approach. However, the PNP equations have not produced the charge inversion/layering in charged wall situations presumably because the conventional PNP does not include the finite size of ions and other physical features needed to create the charge inversion. In this paper, we investigate the key features needed to produce the charge inversion phenomena using a mathematical model, the energy variational approach. One of the key features is a finite size (finite volume) effect, which is an unavoidable property of ions important for their dynamics on small scales. The other is an interfacial constraint to capture the spatial variation of electroneutrality in systems with charged walls. The interfacial constraint is established by the diffusive interface approach that approximately describes the boundary effect produced by the charged wall. The energy variational approach gives us a mathematically self-consistent way to introduce the interfacial constraint. We mainly discuss those two key features in this paper. Employing the energy variational approach, we derive a non-local partial differential equation with a total energy consisting of the entropic energy, electrostatic energy, repulsion energy representing the excluded volume effect, and the contribution of an interfacial constraint related to overall electroneutrality between bulk/bath and wall. The resulting mathematical model produces the charge inversion phenomena near charged walls. We compare the computational results of the mathematical model to those of Monte-Carlo computations.

Keywords:  Finite size effects, energetic variational approach, Poisson-Nernst-Planck equations, Monte-Carlo computations, hard sphere, Lennard-Jones repulsive potential, charge inversion, layering, numerical computations.
Mathematics Subject Classification:  Primary: 00A71, 49S05; Secondary: 65C05, 65N30.

Received: April 2011;      Revised: September 2011;      Available Online: July 2012.