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