| 查看: 1020 | 回复: 3 | ||||
dreamfly_ykf木虫 (正式写手)
|
[求助]
nvt模拟要求表示压力，如何回答？
|
|
如题，使用nvt模拟了一个晶体，审稿人问题如下： If the NVT runs were production run, please show the pressure variation during simulation run as supporting information 做模拟的都知道压力波动很大，尽管输出文件中有压力值，基本上都是在几百Mpa波动，和实验值是相差很大的 要做图获得模拟期间压力的变化，这样结果很不好看，说不定审稿人又提出，压力和实际不符和等问题 请问有人发文章时候，被问及压力波动吗，如何回答啊。 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
职称评审没过，求安慰
已经有49人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
三无产品还有机会吗
已经有6人回复
luftf
木虫 (正式写手)
- 应助: 51 (初中生)
- 金币: 3257.7
- 散金: 685
- 红花: 19
- 帖子: 541
- 在线: 200.2小时
- 虫号: 821742
- 注册: 2009-08-04
- 性别: GG
- 专业: 金属材料的力学行为
2楼2016-02-02 16:46:16
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dreamfly_ykf: 金币+20, ★★★很有帮助, 谢谢 2016-02-24 09:09:18
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
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 模拟EPI: 1
- 应助: 69 (初中生)
- 金币: 8044
- 散金: 222
- 红花: 39
- 帖子: 472
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
4楼2016-02-06 08:38:40