| 查看: 1237 | 回复: 3 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
dreamfly_ykf木虫 (正式写手)
|
[求助]
nvt模拟要求表示压力，如何回答？
|
||||
|
如题，使用nvt模拟了一个晶体，审稿人问题如下： If the NVT runs were production run, please show the pressure variation during simulation run as supporting information 做模拟的都知道压力波动很大，尽管输出文件中有压力值，基本上都是在几百Mpa波动，和实验值是相差很大的 要做图获得模拟期间压力的变化，这样结果很不好看，说不定审稿人又提出，压力和实际不符和等问题 请问有人发文章时候，被问及压力波动吗，如何回答啊。 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有10人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
关于代码变化问题，想知道的进来
已经有17人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
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
luftf
木虫 (正式写手)
- 应助: 51 (初中生)
- 金币: 3257.7
- 散金: 685
- 红花: 19
- 帖子: 541
- 在线: 200.2小时
- 虫号: 821742
- 注册: 2009-08-04
- 性别: GG
- 专业: 金属材料的力学行为
2楼2016-02-02 16:46:16
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 模拟EPI: 1
- 应助: 69 (初中生)
- 金币: 8077
- 散金: 222
- 红花: 39
- 帖子: 473
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
4楼2016-02-06 08:38:40
