一种周期库仑作用势优化法的改进.pdfVIP

第31卷 第 1期 计 算 物 理 Vo1．31．NO．1 2叭4年 1月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Jan．，2014 文章编号：1001—246X(2014)01—0121-06 一 种周期库仑作用势优化法的改进 李名锐 ， 周 刚 ， 初 哲 ， 马 坤 ， 钱秉文 (1．清华大学工程物理系 ，北京 100084；2．西北核技术研究所 ，西安 710024) 摘 要 ：利用简单多项式基组拟合变化平缓的长程作用部分，并增加多个限定条件 以提高优化势与库仑势的局部 重合度，对 Natoli优化法进行改进．在不增加计算复杂度的前提下，得到平均偏差明显降低的优化势；建议优化时划 分r空间与k空间的准则取 r．≥15．在稠密液氢分子模拟中，优化法获得的能量及压力均小于Natoli法与Ewald 法的结果． 关键词 ：长程作用 ；优化法 ；库仑势；Ewald求和 ；多体 中图分类号 ：0411．3 文献标识码 ：A 0 引言 随着计算机与分子模拟技术的不断发展，计算机分子模拟的发展趋势之一是既精确又高效快速地解决 系统内的Ⅳ-体相互作用问题…．采取周期性边界条件的分子模拟 ，不仅要计及模拟盒 内粒子间的相互作用 ， 还要考虑粒子与周期空间中镜像粒子 的相互作用 ，此时传统意义上的截断法 已不再适用于库仑势 的计 算 ，而库仑势的精度对模拟结果的真实性又起着关缝性影响作用 ，这就迫使人们对周期性边界条件下 的库仑作用进行大量研究 ，以解决因收敛速度缓慢而使计算变得异常复杂的难题．其中，应用最为广泛 的是 Ewald求和法 ，它的基本原理是将长程作用势分裂成在实空间与 Fourier空间内均为短程快速收敛 的两部分作用势之和，从而达到较高的收敛速度与计算精度． Natoli等 在极小化平均偏差的基础上 ，对 Ewald求和中实空间与 Fourier空间的划分提出了一种有效 可行的优化法 ，但利用3～5个非奇异的简单线性基组函数来拟合衰减快且原点处奇异的短程作用部分会存 有误差，且拟合出的短程作用势随意性也较大，很难同时严格满足短程作用部分在原点及截断距离处的限定 条件．本文在 Natoli法基础上，利用若干项简单基组拟合变化较为缓慢 的长程作用部分，从考虑周期空间内 粒子与其 自身镜像作用，以及增强优化势与库仑势间的重合度等角度 出发，分析设置 了多个限定条件 ，讨论 了优化法的局部改进性能；在此基础上，以稠密液氢为研究对象，探讨 了优化势对蒙卡模拟的影响． 1 优化法简介 在周期性边界条件下 ，模拟盒 ，J，，J，：，定义库仑作用势 (r)=1／r的镜像势为 Up(1，I)=∑ (1，+LI)=∑vkexp(ik’，)， (1) 式 中，Irl为两带电粒子间的距离 r， 为模拟盒 的布拉伐 晶格矢量 (rtL，nL，r,／L：)，n，n，z／．为整数 ，u为库 仑势t，(r)的Fourier分量 4*rg2 lkI～， =LxLL为模拟盒体积 ，k=21T(nx／L，n／L n／L：)．为克服上式求 和收敛缓慢的问题 ，采用类似于Ewald法 ，将 (1，1)分为在实空间和 Fourier空间均具较高收敛速度的 两部分之和 Oop(I，1)=∑Us(I，+LI)+∑ 2／1(I，．+Lf) ∑u，exP(ik’r)+∑．．u．exP(ik’r) (2) 收稿 151期 ：2013—03—10；修 回 日期 ：2013—06—24 基金项 目：国家部委基金(D413011301)资助项 目 作者简介 ：李名锐(1983一)，男，博士生，从事稠密液氢分子模拟研究，E—mail：limrui—nint@163．eom l22 计 算 物 理