瓦斯流动的分子动力模拟中势函数的研究.pdfVIP

第 22卷 第 5期 山 东 理 工 大 学 学 报 自然 科 学 版 VoL22NO．5 2008年 9月 JournalofShandongUniversityofTechnology NaturalScienceEdition Sep．2008 文章编号：1672～6197 2008 05—0019—03 瓦斯流动的分子动力模拟中势函数的研究 李希建，蔡立勇 贵州大学矿业学院，贵州 贵阳550003 摘 要：势能函数作为分子动力模拟的重要问题，它的选取及求解成为这一研究的热点．探讨 了分 子动力模拟中所涉及的一系列势能函数，选取了所研究的煤层瓦斯流动分子动力模拟所需要的势 函数，并根据临界常数法求取了此函数所涉及到的势能参数。最后计算出一个 甲烷分子与截断半径 处另一分子之间的L—J势能，为后期的模拟研究奠定了基础。 关键词：甲烷；分子动力学；势函数；势参数；临界常数 中图分类号 ：TD712；O313 文献标识码 ：A Potentialfunction inmoleculardynamics simulation offluid_。——flOW methane LIXi-jian，CAILi-yong MiningCollege，GuizhouUniversity，Guiyang550003，China · Abstract：Potentialfunciton isaimportantquestion inmoleculardynamicssimulation．Itsselec— tionandcalculation ishotinthisfield．Severalpotentialfunction，whichinvolvedinmoleculardy— namicssimulationwerediscussed．ThepotentialparameteractordingtocritiCalconstantmethod wascalculated．Atlast，amethaneandanothermethanesL—Jpotentialwasachieyed，whichlay thefoundation forourlaterresearch． Keywords：methane；moleculardynamics；potentialfunction；potentialparameter；criticalcon— stant 在分子动力学 中选择合适 的势能函数非常重 多体势进行描述． 要，它决定了计算工作量和计算模型与真实系统的 近似程度．由于物质系统的复杂性 以及原子间相互 1 势能函数 作用类型的不同，很难得到满足各种不同体系和物 质的一般性而又精确的势函数．针对不同的物质体 1．1 对势 系人们陆续发展了大量经验和半经验公式，分子 或 在分子动力学模拟初期，人们经常将对势用于模 原子 间势函数主要有对势和多体势两种．对势认为 拟中．1957年Alder和Wainwright首先在硬球模型下， 分子 或原子 之间的相互作用是两两之间的作用 ， 开创了利用分子动力学模拟方法研究物质宏观性质的 与其他原子的位置无关 ；在多原子体系中，一个原子 先例，当时采用的就是间断对势口]．1964年 Rahrnan应 的位置不同将影响空间一定范围内的电子云分布， 用连续对势对氩元素进行了研究_2]