24小时热门版块排行榜    

查看: 1743  |  回复: 9

木棉和我

木虫 (职业作家)

[求助] lammps中相同方法前后两次建立的模型在nvt下对应的压力差别对结果有没有影响? 已有1人参与

各位虫友好。用lammps做高温熔体的结构和相变。发现即便采用同样的方法建立非晶态模型,在同样的语句下(如fix 1 all nvt temp 4000 4000 0.1),前后两次建立的模型对应显示的压力差别很大,有1 000 000bar左右波动的,也有在10 000bar左右波动的。我想知道两次建模的模型下,对应压力差别这么大,其对模拟结果有没有影响?如果有,怎样消除这样影响?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
用 fix npt 控温控压。否则你无法控制压强。

[ 发自手机版 http://muchong.com/3g ]
2楼2018-06-11 02:25:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木棉和我

木虫 (职业作家)

针对计算结构信息的情况,只用nvt是可以的,但是就会出现这种情况。

发自小木虫Android客户端
3楼2018-06-11 02:33:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木棉和我

木虫 (职业作家)

引用回帖:
2楼: Originally posted by brucefan at 2018-06-11 02:25:38
用 fix npt 控温控压。否则你无法控制压强。

问题是npt和nvt不能同时联用啊。

发自小木虫Android客户端
4楼2018-06-11 02:34:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 木棉和我 at 2018-06-11 02:34:45
问题是npt和nvt不能同时联用啊。
...

fix npt 是即控温也控压,实现等温等压系综。

[ 发自手机版 http://muchong.com/3g ]
5楼2018-06-11 02:36:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

另外,你不要太关注压强和温度的涨落是多大(如果你不是故意研究这些涨落),而且要关注分块(例如每 ps平均一次)平均后的温度和压强变化趋势。画图也有助于自己和别人理解。

[ 发自手机版 http://muchong.com/3g ]
6楼2018-06-11 02:40:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木棉和我

木虫 (职业作家)

引用回帖:
5楼: Originally posted by brucefan at 2018-06-11 02:36:51
fix npt 是即控温也控压,实现等温等压系综。
...

关注结构信息,很少看到在npt下统计的,一般都是在nvt或nve下统计。另外,如果上来就用nvt,那么显示的压力的大小好像与建模时设定的密度有关。如果设置的密度偏大,那么压力会特别大,就会出现out of range的情况。

发自小木虫Android客户端
7楼2018-06-11 11:01:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木棉和我

木虫 (职业作家)

引用回帖:
6楼: Originally posted by brucefan at 2018-06-11 02:40:55
另外,你不要太关注压强和温度的涨落是多大(如果你不是故意研究这些涨落),而且要关注分块(例如每 ps平均一次)平均后的温度和压强变化趋势。画图也有助于自己和别人理解。
...

我用的步长0.001fs,热力学信息每1000步输出一次,正好是1ps输出一次。如果算1ps下的热力学评论,可以设置每100步输出一次,导出后再进行评论?

发自小木虫Android客户端
8楼2018-06-11 11:04:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木棉和我

木虫 (职业作家)

引用回帖:
6楼: Originally posted by brucefan at 2018-06-11 02:40:55
另外,你不要太关注压强和温度的涨落是多大(如果你不是故意研究这些涨落),而且要关注分块(例如每 ps平均一次)平均后的温度和压强变化趋势。画图也有助于自己和别人理解。
...

樊老师。分块平均是数据输出后的后处理中完成,还是lammpsz中有命令可以直接输出分块的结果?
9楼2018-06-12 16:25:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
木棉和我: 金币+20, ★★★★★最佳答案, 非常感谢 2018-06-12 18:02:22
引用回帖:
9楼: Originally posted by 木棉和我 at 2018-06-12 16:25:32
樊老师。分块平均是数据输出后的后处理中完成,还是lammpsz中有命令可以直接输出分块的结果?...

自己算就行。lammps只要输出数据即可。我是说你问问题的时候可以用图展示你的结果,不然问题太抽象。你最近问的几个问题都是简单的控温控压问题,那么你需要的是观察温度和压强的时间演化如何依赖于你的模拟参数。自这样你就会通过测试获得经验。

[ 发自手机版 http://muchong.com/3g ]
10楼2018-06-12 17:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 木棉和我 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 化学工程321分求调剂 +9 大米饭! 2026-03-15 9/450 2026-03-16 16:41 by 我的船我的海
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 考研化学学硕调剂,一志愿985 +3 张vvvv 2026-03-15 3/150 2026-03-16 16:36 by houyaoxu
[考研] 东南大学364求调剂 +4 JasonYuiui 2026-03-15 4/200 2026-03-16 08:36 by Linda Hu
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见