本文利用PYthon編寫了分子動力學模擬程序,并利用該程序?qū)e分子體系進行了詳細研究。分別研究了不同初始條件,不同邊界條件,截斷點位置等的研究。在分子數(shù)為500,分子初速度為500的初始條件下,體系最終經(jīng)過2.25e-11s發(fā)展成為麥克斯韋平衡體系。分別研究周期性邊界條件與剛性邊界對系統(tǒng)發(fā)展的影響,研究發(fā)現(xiàn)周期性邊界條件與剛性邊界條件最終對系統(tǒng)的發(fā)展沒有明顯影響。對截斷點進行研究后發(fā)現(xiàn),必須選取大于1.1*sigma(sigma為L-J勢能中的常數(shù),由實驗可測定)作為截斷,如果小于該值,體系由于數(shù)值誤差將偏離物理實際。同時研究發(fā)現(xiàn)勢函數(shù)對于體系的發(fā)展也具有重要影響。
標簽：
PYthon
上傳時間：
2022-06-19
上傳用戶：shjgzh
