lennard-jones 流体固液相变的分子动力学模拟 - 南京工业大学学报 .pdf

第 19 卷第 2 期 南 京 化 工 大 学 学 报 Vol . 19 No. 2 1997 年 4 月 JOU RN AL O F NA NJ IN G U N IV ERSIT Y OF CH EM ICAL T ECHNOL OGY Apr . 1997 Lennard-Jones 周 健 陆小华 王延儒 时 钧 ( 南京化工大学化学工程系, 南京, 210009) 采用等温等容分子动力学模拟方法, 研究了 Lennard- J ones 流体从固相到液相体系结构 和性质的 化 。均方位移, 径向分布函数, 平动序参数等结构和性质的模拟表明体系经历了固液 相 。在对比密度为 0. 84 , 0 . 86, 0 . 88 和 0. 9 时, 体系熔点分别为 0 . 66~ 0 . 67 , 0. 71~ 0 . 72, 0. 82 ~ 0 . 83 和 0. 90 ~ 0. 91 。 L ennard- J ones 流体 固液相 正则系综( N VT ) 分子动力学( M D) O 642. 51 在许多工业生产中, 如冶金、聚合物的加工、熔融结晶分离混合物、非晶态合金材料制备、 半导体材料精制等都涉及到物质的熔解过程。对于熔解过程的机理, 目前尚无成熟的理论, 而 采用分子动力学( MD) 的分子模拟方法则有助于我们从分子水平上对上述过程加以理解。分 子模拟是从某一设定的位能模型出发, 在一定的条件下, 由体系粒子的微观性质来求算其宏观 性质, 它可以被用来检验理论的正确与否, 又可以将其结果与实验值相比较, 籍以检验和改进 模型, 同时它还可以模拟得到某些极限条件以及实验无法实现的情况下体系的一些信息, 因此 近二十年来计算机分子模拟得到了迅猛发展 1] 。Lennard-Jones 流体, 由于它具有真实流体之 间相互吸引和相互排斥的基本特征, 在计算机分子模拟中被广泛采用。 A lder 和 Wainw right 2] 首次对硬球流体进行了分子动力学模拟, 证明体系存在固体流体 相 。H ansen 等 3] 用 Monte Carlo 方法模拟了 L ennard-Jones 流体 的相 过程。Chokappa 等 4] 用等压分子动力学方法模拟了L ennard-Jones 流体的熔融过程。Anderson 5] 首先提出了 定压, 定温和定温定压三种新的平衡分子动力学方法, 其中定常温度 M D 方法最引入注目, 它 相当于模拟 N 、V 、T 不 的正则系综。因此采用定温 MD 方法模拟 Lennard-Jones 流体的熔 解过程是很有意义的。 1 M D 1984 年, Nose 6] 等提出了NVT 系综的扩展系统法。系统的H am ilton 函数可写为 国家 自然科学基金资助项 目 收稿 日期: 1996- 10- 18 第 2 期 周健等: Lennard- Jones 流体固液相 的分子动力学模拟 21 2 2 2 H = p / 2m s + ( q) + p / 2 Q + (f + 1) k T lns ( 1) i i s i 2 r i 2 2 s r i