24小时热门版块排行榜    

查看: 1236  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 十年后又回来了,论文投稿求助 +5 哈哈114477 2026-08-01 5/250 2026-08-08 15:14 by 仁砚薪传
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] fileCode有新解读? +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见