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