24小时热门版块排行榜    

查看: 1049  |  回复: 3

dreamfly_ykf

木虫 (正式写手)

[求助] nvt模拟要求表示压力,如何回答?

如题,使用nvt模拟了一个晶体,审稿人问题如下:
If the NVT runs were production run, please show the pressure variation during simulation run as supporting information
做模拟的都知道压力波动很大,尽管输出文件中有压力值,基本上都是在几百Mpa波动,和实验值是相差很大的
要做图获得模拟期间压力的变化,这样结果很不好看,说不定审稿人又提出,压力和实际不符和等问题
请问有人发文章时候,被问及压力波动吗,如何回答啊。
回复此楼
自古明王圣帝,犹须勤学,况凡庶乎!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luftf

木虫 (正式写手)

【答案】应助回帖

实际上压力为什么会波动那么大呢?比起加载的载荷GPa,几百MPa也没什么吧
2楼2016-02-02 16:46:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dreamfly_ykf: 金币+20, ★★★很有帮助, 谢谢 2016-02-24 09:09:18
瞬时压力波动大是NVE/NVT模拟中很常见的现象,为了减小波动最有效的方法是增加体系尺寸,因为统计力学中涨落是和粒子数的平方根成反比的。

http://www.gromacs.org/Documentation/Terminology/Pressure

这个问题之前在lammps mailing list里也有人问过,唯一的解决方法就是增加体系尺寸,你可以看看这个帖子及相关回帖
http://lammps.sandia.gov/threads/msg33152.html

如果你的体系比较小,那么瞬时压力波动大是再正常不过了。

你可以计算一下你的压力均值,看看跟设定值有多大差距。或者,可以像下面这张图,每隔一定步数(比如每100步)对压力取平均作图,如果这样计算出来的压力和你设定值差距不大就没问题。
http://www.bevanlab.biochem.vt.e ... zyme/07_equil2.html
3楼2016-02-04 07:27:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oxox6085

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dreamfly_ykf: 金币+10, 有帮助, 谢谢 2016-02-24 09:10:39
波动大正常,放到SI正文里代一笔就可以了
有些事情可能只有计算机模拟能告诉我们答案
4楼2016-02-06 08:38:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dreamfly_ykf 的主题更新
信息提示
请填处理意见