LAMMPS里的溫度是靠原子動(dòng)能算出來(lái)的，具體可以看官方手冊(cè)里compute temp那部分。
我們平時(shí)用dump命令輸出的lammpstrj文件，里面存了每個(gè)原子在各個(gè)時(shí)刻的位置、速度、動(dòng)能、勢(shì)能這些信息（當(dāng)然得在dump里設(shè)好了才會(huì)存）。只要你這個(gè)文件里有動(dòng)能數(shù)據(jù)，畫(huà)溫度分布圖就沒(méi)問(wèn)題。
推薦一個(gè)方法：用OVITO導(dǎo)入lammpstrj然后上色。
OVITO是專門處理LAMMPS結(jié)果的免費(fèi)軟件，Win和Linux都能用，搜一下官網(wǎng)就能下。它能按原子的某個(gè)屬性來(lái)上色，比如x坐標(biāo)、動(dòng)能啥的，操作特別簡(jiǎn)單。不過(guò)缺點(diǎn)是圖出來(lái)是球狀原子模型，屬于離散的點(diǎn)，不是那種連續(xù)平滑的分布圖。比如有張圖顯示石墨被C60撞了之后表面溫度分布（側(cè)面視角），每個(gè)原子的溫度都標(biāo)出來(lái)了，看得清清楚楚。
我們平時(shí)用dump命令輸出的lammpstrj文件，里面存了每個(gè)原子在各個(gè)時(shí)刻的位置、速度、動(dòng)能、勢(shì)能這些信息（當(dāng)然得在dump里設(shè)好了才會(huì)存）。只要你這個(gè)文件里有動(dòng)能數(shù)據(jù)，畫(huà)溫度分布圖就沒(méi)問(wèn)題。
推薦一個(gè)方法：用OVITO導(dǎo)入lammpstrj然后上色。
OVITO是專門處理LAMMPS結(jié)果的免費(fèi)軟件，Win和Linux都能用，搜一下官網(wǎng)就能下。它能按原子的某個(gè)屬性來(lái)上色，比如x坐標(biāo)、動(dòng)能啥的，操作特別簡(jiǎn)單。不過(guò)缺點(diǎn)是圖出來(lái)是球狀原子模型，屬于離散的點(diǎn)，不是那種連續(xù)平滑的分布圖。比如有張圖顯示石墨被C60撞了之后表面溫度分布（側(cè)面視角），每個(gè)原子的溫度都標(biāo)出來(lái)了，看得清清楚楚。