分子动力学模拟中的变截断半径算法.pdfVIP

第27卷第4期 工程热物理学报 V01．27．No．4 2006年7月 JOURNALoFENGINEERINGTHERMOPHYSICSJul．．2006 分子动力学模拟中的变截断半径算法 蔡治勇 曾丹苓 刘娟芳 (重庆大学动力工程学院， 重庆 400044) 摘要 在分子动力学模拟中，人们往往根据所研究的对象选取一个固定的截断半径，而不同的作者所推荐的有效截断 半径值甚为分散．这一方面影响了模拟结果的精度，同时也造成了截断半径选取的困难。本文提出了一种变截断半径的新 算法，可根据对象所处的状态灵活地改变截断半径。文中以两相平衡系的非均匀各向异性区的模拟为例，采用所提出的方 法得到了更接近实验的结果，同时也节约了计算时间。 关键词，变截断半径；分子动力学模拟；算法 中图分类号：TKl23文献标识码：A 文章编号：0253—231x(2006)04一0556一03 ACoMPU7IlArIIoNALMETHoDWITHVARIABLECUToFF RADIUMINMoLECULARDYNAMICSSIMULAIoN CAI ZENG LIU Zhi—Ybng Dan—LingJuaI卜Fang ofP0wer Unjyersj亡y， (cd垤e En咖e刑ng，chongqjngChongqjng400044，∞ina) AbstractInm01ecular simulationa丘xedcutoffradiusis chosenfbrthe dynamics usuaUy calculation， andtherecommendedvalueofitis itcausestroubleofitsselectionandalsoinfluences scattering，thus the oftheresults．Inthe anew methodwithV甜iablecutoff precision presentpaper computational radiumis forthesimulationofatwo is to such suggested．Anexample phasesystemproVidedsupport asuggestion． method Key、Ⅳords、，ariablecuto圩radium；m01eculardynamicssimulation；computational 1前 言 匀相中表面张力的计算结果。结果表明：在均匀相 中取固定截断半径R=3．o一4．o，在非均匀相中， 在分子动力学模拟中，截断半径的选取对模拟 根据对象所处的状态灵活取截断半径时，各相的热 结果的影响甚大。传统的方法在模拟过程中通常选 物性计算能够同时得到满意的结果。与传统的定截 取一个固定的截断半径，由于气一液相平衡体系中， 断半径方法相比，变截断半径法计算精度更高，计 均匀相、非均匀相的分布特性不同，用统一的截断