LJ势氩系统分子动力学模拟中截断半径的选择  

Application of cutoff distance selection in molecular dynamics simulation of LJ argon system

作　　者：孙晨阳 侯超峰[1] 葛蔚[1,2] Chenyang SUN;Chaofeng HOU;Wei GE(Institute of Process Engineering,Chinese Academy of Sciences,Beijing 100190,China;School of Chemical Engineering,University of Chinese Academy of Sciences,Beijing 100049,China)

机构地区：[1]中国科学院过程工程研究所,北京100190 [2]中国科学院大学化学工程学院,北京100049

出　　处：《过程工程学报》2021年第3期259-264,共6页The Chinese Journal of Process Engineering

基　　金：国家自然科学基金(编号:21776280,91834303,91934302);中国科学院资助项目(编号:XDC01040100,QYZDJ-SSW-JSC029);北京市自然科学基金委员会-北京市教育委员会联合资助项目(编号:KZ201910017019);多相复杂系统国家重点实验室自主研究课题(编号:MPCS-2019-A-10)。

摘　　要：近年来在分子动力学方法研究LJ(Lennard-Jones)势氩系统时,越来越多的计算建议采用4.5σ甚至更大的截断半径(σ为氩原子直径)。本工作研究了不同截断半径对等温等压NPT系综氩系统熔沸点相图的影响,分析了不同截断半径下在氩系统的熔点及液相区域不同热力学状态点的径向分布函数和速度自相关函数。结果表明,以熔点为基准,在距离其相同液态温区分率的热力学温度点能获得相同的热力学性质,2.5σ的截断半径在模拟的准确性和计算性能上均能满足模拟的要求,本工作为液氩模拟中截断半径的选择指明了一个新的思考方向。In Lennard-Jones(LJ)potential argon system investigated by molecular dynamics simulation,the different cutoff distances are frequently applied to calculate the interactive force between atoms,and some reports have given out the influence of the cutoff distances on the simulation systems.More and more calculations suggest to use 4.5?or even larger truncation distances(?is the diameter of argon atom)to obtain the more accurate thermodynamic properties of the systems.In this work,a simple method was proposed to solve the problem of superheating encountered in the calculation of melting point by direct heating,where an independent track and ensemble at each temperature point are run.And then,the effect of different cutoff distances on the phase diagrams of the melting and boiling points of argon system in the NPT ensemble was studied.The melting point was in good agreement with the experimental and theoretically calculated values when the cutoff distance of 2.5?was used.However,the deviation from the experimental melting point became more evident when the larger cutoff distances were employed.In order to find out the underlying mechanism behind the deviation,the radial distribution functions and velocity autocorrelation function at the melting points and different thermodynamic states of the liquid argon with different cutoff distances were analyzed.It was found that the same thermodynamic properties can be obtained at the corresponding thermodynamic state points under different truncation distances.The mapping between the different thermodynamic state points was understandable due to the different thermodynamic states at the same temperatures under the varied truncation distances,and was beneficial to significantly reduce the computational workload at the smaller cutoff distance.This work proposed an exploratory way for the selection of the cutoff distance in the simulation of liquid argon,where the truncation distance of 2.5?can meet the requirements of computational accuracy and performance in the simulations.

关 键 词：分子动力学模拟 LJ势 截断半径 液氩 

分 类 号：TQ460.64[化学工程—制药化工]