| 4 | 1/1 | 返回列表 |
| 查看: 1090 | 回復(fù): 3 | ||||
dreamfly_ykf木蟲 (正式寫手)
|
[求助]
nvt模擬要求表示壓力,如何回答?
|
|
如題,使用nvt模擬了一個(gè)晶體,審稿人問題如下: If the NVT runs were production run, please show the pressure variation during simulation run as supporting information 做模擬的都知道壓力波動(dòng)很大,盡管輸出文件中有壓力值,基本上都是在幾百M(fèi)pa波動(dòng),和實(shí)驗(yàn)值是相差很大的 要做圖獲得模擬期間壓力的變化,這樣結(jié)果很不好看,說不定審稿人又提出,壓力和實(shí)際不符和等問題 請(qǐng)問有人發(fā)文章時(shí)候,被問及壓力波動(dòng)嗎,如何回答啊。 |
分子動(dòng)力學(xué)模擬 |

木蟲 (正式寫手)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
瞬時(shí)壓力波動(dòng)大是NVE/NVT模擬中很常見的現(xiàn)象,為了減小波動(dòng)最有效的方法是增加體系尺寸,因?yàn)榻y(tǒng)計(jì)力學(xué)中漲落是和粒子數(shù)的平方根成反比的。 http://www.gromacs.org/Documentation/Terminology/Pressure 這個(gè)問題之前在lammps mailing list里也有人問過,唯一的解決方法就是增加體系尺寸,你可以看看這個(gè)帖子及相關(guān)回帖 http://lammps.sandia.gov/threads/msg33152.html 如果你的體系比較小,那么瞬時(shí)壓力波動(dòng)大是再正常不過了。 你可以計(jì)算一下你的壓力均值,看看跟設(shè)定值有多大差距;蛘,可以像下面這張圖,每隔一定步數(shù)(比如每100步)對(duì)壓力取平均作圖,如果這樣計(jì)算出來的壓力和你設(shè)定值差距不大就沒問題。 http://www.bevanlab.biochem.vt.e ... zyme/07_equil2.html |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +129 |

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|