论文部分内容阅读
溶液中离子的关联效应可以描述离子间的相互作用,包括尺寸效应和静电关联作用等。在溶液体系模拟中,它可以对离子浓度、静电势能以及其他一些相关物理性质现象等给出合理有效的实验预测。在本文中,我们介绍溶液中的离子关联效应,将其加入到传统连续模型中,通过模型分析和数值模拟研究其对溶剂化能、电流、浓度分布以及离子选择性的影响与作用。尺寸修正Poisson-Nernst-Planck(SMPNP)或尺寸修正Poisson-Boltzmann(SMPB)模型在传统连续模型中加入了尺寸效应。分析表明在平衡态下,尺寸修正模型中根据自由能泛函推导出的溶剂化能存在一个简单的近似计算公式,与Poisson-Boltzmann(PB)模型中广泛应用的近似公式完全相同。模拟结果发现,加入尺寸效应后溶剂化能的绝对值和电流都会减小(实际溶剂化能值升高)。溶剂化能的近似计算与准确计算结果相差不大,因此尺寸效应对溶剂化能的影响主要体现在近似计算部分;对电流各部分的数值结果表示其降低主要体现在电场驱动项。 本研究统计力学出发新导出了一个局域近似的三维LHSPNP模型。在自由能泛函中增加超额Helmholtz自由能,结合密度泛函理论(DFT),运用Rosenfeld的FMT表示硬球之间的排斥作用,并对积分形式的硬球超额化学势进行局域近似就得到简化的LHSPNP模型。在均匀分布的单组分溶液中,通过状态方程的研究发现,FMT可以成功预测前三项维里系数,推导得出的LHS模型能够精确给出前两项,而SMPNP模型则只能精确确定第一项。对带电球形粒子的模拟结果表明,LHSPNP和SMPNP的结果在很多情形下都比较接近,但是与传统连续模型的结果相差很大。在尺寸修正模型和LHSPNP模型的数值计算方面,我们将电势分解为三部分加和处理与固定电荷对应的Dirac-delta函数项,用有限元方法和Gummel方法对耦合偏微分方程迭代求解,结合并行自适应有限元程序包PHG对大规模计算进行并行加快计算速度。同时为了保证收敛性,在迭代过程中加入对电势和浓度的松弛算法。在平板、圆柱和球形对称构型中,FMT及MFMT中相应的加权密度和硬球超额化学势表达式中的积分都可以进行简化计算,对此我们给出了严格的推导。在平衡态下圆柱构型中,我们结合Poisson方程,MFMT和静电关联作用模拟离子分布,研究其中的离子选择性。根据对称性,电势和浓度都只与一个自变量即到中轴的径向距离有关,此时Poisson方程可以约化为常微分方程,结合平衡态下浓度的表达式进行耦合求解。我们发现对于只有尺寸差异的两种反离子,在圆柱半径小的时候尺寸大的离子选择性通过的机会更大,而圆柱半径大时则恰恰相反。