熔融KCl中微观空洞及其分布的计算机模拟研究.pdfVIP

2008年 12月 应用数学与计算数学学报 第22卷 第2期 Dec．，2008 COMM ．ON APPL．MATH．AND COMPUT V01．22 No．2 熔融 KC1中微观空洞及其分布的计算机模拟研究 苏会娟 汤正诠 摘要 以分子动力学计算机模拟所得KCI瞬态构型为原始数据，本文计算了不同温 度的系统中存在的微观空洞并用径向分布函数描述微观空洞的统计分布，进而讨论了微观 空洞的分布随温度变化的特点和不同半径的微观空洞的分布特点．我们也讨论了微观空洞 的体积与系统 自扩散系数、粘滞性之间的关系． 关键词 计算机模拟，微观空洞，统计分布 TheCom puterSimulation Study ofM icroscopic Holesand TheirD istribution in M olten KCl SuHuijuan TangZhengquan Abstract TakingtheKCIinstantaneousconfigurationsobtainedbycomputersim— ulationwiththemoleculardynamicsmethodassourcedata，wecalculatethemicroscopic holesinKC1system atdifferenttemperaturesanddescribetheirstatisticaldistribution usingtheradialdistributionfunction．Furthermore．wediscussthecharacteristicthatthe distribution ofmicroscopicholeschangeswithtemperatureandthedistribution chra ac— teristicofmicroscopicholesofdifferentradii．W ealsodiscusstherelationship amongthe volumeofmicroscopicholes，selhdiffusioncoefficientandshera viscosity． Keywordscomputersimulation，microscopichole，statisticaldistribution l 问题的提出 分子动力学计算机模拟能从时间和空间两方面再现系统的宏观性质，同时又储藏 着大量的微观信息，是结构 ．功能研究的有效方法，如Woodcock对简单离子系统所 做的工作 l【】_微观空洞的分布规律不仅对理解扩散、电导、粘度等性质，而且对描述 系统的结构特征和状态变化都是十分重要的 [213】． 本文的计算表 明在熔融 KC1系中存在与离子尺寸相当的空洞，它们在熔化过程中 迅速地增长． 2 计算方法 根据分子动力学计算机模拟所获的相轨道，在周期性边界条件下，若能在模拟立 方元胞中的某位置加入一个不比组分粒子小的 虚“粒子 ”，而且该 虚“粒子 ”又不与 本文 2008年 7月2日收到． 1．上海大学理学院数学系，上海，200444；DepartmentofMathematics，ShanghaiUniversity，Shanghai 200444，China 102 应用数学与计算数学学报 22卷 任何其它粒子 (包括组分粒子和先前已经加入的 虚“粒子 ”)相重叠 (包括部分重叠)， 则称在该处存在一个微观空洞。可容纳的最大虚粒子半径就是该位置存在的微观空洞 的半径． 为嵌入虚粒子，必须对 已有粒子 (包括组分粒子和 已经嵌入的 虚“粒子 ”)进行搜 索．搜索 以四个粒子组成的粒子组为单元．在模拟立方元胞中依次选择一个粒子为母 粒子，每选定一个母粒子，在它周 围充分大的范围内任意选择三个粒子，这四个粒子 便构成了一个单元组．为了使四粒子间嵌入的虚粒子最大，应当将虚粒子的中心选得 与四个粒子的表面等距．求得