Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602-5700, USA.
J Phys Chem B. 2012 Sep 13;116(36):11356-61. doi: 10.1021/jp305400z. Epub 2012 Aug 29.
Most theoretical studies of an electrical double layer, which is formed by an electrolyte in contact with a charged electrode, employ a primitive model in which the solvent is represented by a dielectric continuum. This implicit-solvent model is convenient because computations are comparatively simple. However, it suppresses oscillations in the density profiles of ionic species that result from the discreteness of the solvent molecules. Furthermore, the implicit-solvent model yields poor results for the capacitance. In comparison with experiment at fixed electrode charge density, it predicts a too small electrode potential, and the resultant capacitance is too large. This latter discrepancy can be compensated in part by postulating the existence of an often fictitious inner layer whose properties are parametrized to agree best with experiment. The use of an implicit solvent model and an inner layer helps in correlating experimental results but rests on a faulty microscopic picture. Unfortunately, explicit consideration of solvent molecules poses both theoretical and numerical difficulties and, as a result, studies using an explicit solvent model have been few and far between. In this study, we consider a simple nonprimitive or explicit solvent model in which each solvent molecule is represented by a dimer composed of touching positive and negative hard spheres, with a resulting dipole moment that is equal to that of a water molecule, and the ions are represented by charged hard spheres. The density profiles and charge-potential relationship of this model are examined using the classical density functional theory. We find that the introduction of an explicit solvent increases the electrode potential, at fixed electrode charge, without the need to postulate a parametrized inner layer. Because of the solvent polarity, the ion profiles become strong oscillatory and show local charge inversion near a highly charged electrode surface at all ion concentrations.
大多数关于双电层的理论研究都是在电解质与带电电极接触的情况下进行的,这些研究采用了一种原始模型,其中溶剂由介电连续体表示。这种隐溶剂模型很方便,因为计算相对简单。然而,它抑制了溶剂分子离散性导致的离子物种密度分布的振荡。此外,隐溶剂模型对电容的预测结果较差。与固定电极电荷密度下的实验相比,它预测的电极电位太小,因此电容太大。后一种差异可以通过假设存在一个通常是虚构的内层来部分补偿,该内层的性质被参数化,以与实验最佳匹配。隐溶剂模型和内层的使用有助于关联实验结果,但基于有缺陷的微观图像。不幸的是,溶剂分子的显式考虑既存在理论上的困难,也存在数值上的困难,因此,使用显式溶剂模型的研究很少。在这项研究中,我们考虑了一种简单的非原始或显式溶剂模型,其中每个溶剂分子由一个由接触的正、负硬球组成的二聚体表示,从而产生一个等于水分子的偶极矩,而离子则由带电硬球表示。使用经典密度泛函理论研究了该模型的密度分布和电荷-电位关系。我们发现,在固定电极电荷的情况下,引入显式溶剂会增加电极电位,而无需假设参数化的内层。由于溶剂的极性,离子分布变得强烈振荡,并在所有离子浓度下在带高电荷的电极表面附近显示局部电荷反转。