应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps中相同方法前后两次建立的模型在nvt下对应的压力差别对结果有没有影响？
【调剂】北京石油化工学院2024年16个专业接受调剂
101/1返回列表
查看: 1427  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

木棉和我

木虫 (著名写手)

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

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

» 猜你喜欢
1楼2018-06-10 23:05:53
已阅   回复此楼   关注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 ]
赞一下(1人) 回复此楼
10楼2018-06-12 17:44:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 木棉和我 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 24年 申博 化学/材料 一作6篇sci +3 wangyp123 2024-04-23 4/200 2024-04-23 23:38 by 莹亮的星空
[论文投稿] 期刊推荐 20+3 木颜尘ip 2024-04-22 5/250 2024-04-23 23:03 by topedit
[考博] 真的好想读博！ +7 wangzhe_bs 2024-04-22 9/450 2024-04-23 22:20 by wangzhe_bs
[找工作] 工作问题求助 +3 向前?@ 2024-04-18 5/250 2024-04-23 21:25 by 向前?@
[教师之家] 期末给学生划重点都是什么话术啊 +13 luokereng 2024-04-20 15/750 2024-04-23 21:07 by WORLD0256
[论文投稿] 最近遇到这样一个问题 +4 asd123gfa689 2024-04-18 11/550 2024-04-23 20:14 by asd123gfa689
[论文投稿] 无under review，直接DIP，咋么办。 +10 lizhengke06 2024-04-19 16/800 2024-04-23 19:15 by 化学程序员
[考研] 0854－0855调剂 +5 shangannum1 2024-04-21 8/400 2024-04-23 17:08 by Dam026
[教师之家] 国内科研环境差导致国内科研水平差的根源 +8 zju2000 2024-04-18 8/400 2024-04-23 10:40 by Quakerbird
[论文投稿] 编辑返稿让改格式，这个时候能大修内容吗？ +6 双倍好运锦鲤 2024-04-17 9/450 2024-04-23 09:49 by fslaser
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
[教师之家] 河海大学常州校区新进的讲师待遇如何？考核要求如何？ +3 德国小球迷 2024-04-19 18/900 2024-04-22 19:23 by ZHONGWU_U
[论文投稿] 国家重点研发计划论文可以发MDPI里的期刊吗 +6 lin_2018 2024-04-17 6/300 2024-04-22 17:32 by 20094958
[考研] 博士生招生！！ +3 2458827884 2024-04-19 3/150 2024-04-22 12:37 by 钰波love
[有机交流] 兄弟们帮我看看这两个结构怎么合成 +6 xl2088131 2024-04-17 8/400 2024-04-21 22:39 by zyqchem
[硕博家园] 材料工程，二本硕士，值得读吗 +3 简风23 2024-04-20 6/300 2024-04-21 12:19 by 简风23
[考研] 298求调剂 +10 孙大大@ 2024-04-17 11/550 2024-04-20 18:33 by 371122gyd
[考研] 312求调剂 +3 Leroic000 2024-04-19 3/150 2024-04-20 18:19 by xinaishan
[论文投稿] 一审一个审稿人，小修，会怎么样呀？ +10 林师妹 2024-04-18 10/500 2024-04-20 18:08 by 肖虫家
[有机交流] 紫外光谱 50+3 54胡 2024-04-17 3/150 2024-04-19 10:48 by Nanamiwww
信息提示
关闭
请填处理意见
关闭