中国物理B ›› 2010, Vol. 19 ›› Issue (4): 48702-048702.doi: 10.1088/1674-1056/19/4/048702

• • 上一篇    下一篇

Finite size effect of ions and dipoles close to charged interfaces

童朝晖, 诸跃进   

  1. Department of Physics, Ningbo University, Ningbo 315211, China
  • 收稿日期:2009-06-20 修回日期:2009-09-27 出版日期:2010-04-15 发布日期:2010-04-15
  • 基金资助:
    Project supported by the National Natural Science Foundation of China (Grant Nos.~20954001 and 10774079), the Natural Science Foundation of Zhejiang Province of China (Grant No.~Y7080401), the Natural Science Foundation of Ningbo City (Grant No.~2009A6100

Finite size effect of ions and dipoles close to charged interfaces

Tong Chao-Hui(童朝晖) and Zhu Yue-Jin(诸跃进)   

  1. Department of Physics, Ningbo University, Ningbo 315211, China
  • Received:2009-06-20 Revised:2009-09-27 Online:2010-04-15 Published:2010-04-15
  • Supported by:
    Project supported by the National Natural Science Foundation of China (Grant Nos.~20954001 and 10774079), the Natural Science Foundation of Zhejiang Province of China (Grant No.~Y7080401), the Natural Science Foundation of Ningbo City (Grant No.~2009A6100

摘要: The modified dipolar Poisson--Boltzmann (MDPB) equation, where the electrostatics of the dipolar interactions of solvent molecules and also the finite size effects of ions and dipolar solvent molecules are explicitly taken into account on a mean-field level, is studied numerically for a two-plate system with oppositely charged surfaces. The MDPB equation is solved numerically, using the nonlinear Multigrid method, for one-dimensional finite volume meshes. For a high enough surface charge density, numerical results of the MDPB equation reveal that the effective dielectric constant decreases with the increase of the surface charge density. Furthermore, increasing the salt concentration leads to the decrease of the effective dielectric constant close to the charged surfaces. This decrease of the effective dielectric constant with the surface charge density is opposite to the trend from the dipolar Poisson--Boltzmann (DPB) equation. This seemingly inconsistent result is due to the fact that the mean-field approach breaks down in such highly charged systems where the counterions and dipoles are strongly attracted to the charged surfaces and form a quasi two-dimensional layer. In the weak-coupling regime with the electrostatic coupling parameter (the ratio of Bjerrum length to Gouy--Chapman length) $\varXi < 1$, where the MDPB equation works, the effective dielectric constant is independent of the distance from the charged surfaces and there is no accumulation of dipoles near the charged surfaces. Therefore, there are no physical and computational advantages for the MDPB equation over the modified Poisson--Boltzmann (MPB) equation where the effect of dipolar interactions of solvent dipoles is implicitly taken into account in the renormalised dielectric constant.

Abstract: The modified dipolar Poisson--Boltzmann (MDPB) equation, where the electrostatics of the dipolar interactions of solvent molecules and also the finite size effects of ions and dipolar solvent molecules are explicitly taken into account on a mean-field level, is studied numerically for a two-plate system with oppositely charged surfaces. The MDPB equation is solved numerically, using the nonlinear Multigrid method, for one-dimensional finite volume meshes. For a high enough surface charge density, numerical results of the MDPB equation reveal that the effective dielectric constant decreases with the increase of the surface charge density. Furthermore, increasing the salt concentration leads to the decrease of the effective dielectric constant close to the charged surfaces. This decrease of the effective dielectric constant with the surface charge density is opposite to the trend from the dipolar Poisson--Boltzmann (DPB) equation. This seemingly inconsistent result is due to the fact that the mean-field approach breaks down in such highly charged systems where the counterions and dipoles are strongly attracted to the charged surfaces and form a quasi two-dimensional layer. In the weak-coupling regime with the electrostatic coupling parameter (the ratio of Bjerrum length to Gouy--Chapman length) $\varXi < 1$, where the MDPB equation works, the effective dielectric constant is independent of the distance from the charged surfaces and there is no accumulation of dipoles near the charged surfaces. Therefore, there are no physical and computational advantages for the MDPB equation over the modified Poisson--Boltzmann (MPB) equation where the effect of dipolar interactions of solvent dipoles is implicitly taken into account in the renormalised dielectric constant.

Key words: dipole, Poisson--Boltzmann equation, charged interface, polarisation

中图分类号:  (Permittivity (dielectric function))

  • 77.22.Ch
73.25.+i (Surface conductivity and carrier phenomena) 41.20.Cv (Electrostatics; Poisson and Laplace equations, boundary-value problems)