24小时热门版块排行榜    

查看: 4525  |  回复: 12

lucaca

铁杆木虫 (著名写手)


[交流] NPT系综下判断是否平衡的标志

跑了一段NPT,600K,现在温度变化大概是600±40K,这个浮动可以判断体系在这一系综下达到平衡了吗?
回复此楼

» 收录本帖的淘帖专辑推荐

晶体模型计算分析 MDs-Gromacs

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)


★
lucaca(金币+1): 谢谢参与
判断是否平衡的话可以用能量啊,RMSD之类的信息
2楼2013-07-05 23:54:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbn

银虫 (小有名气)


★
lucaca(金币+1): 谢谢参与
NPT判断是否平衡难道不是看密度的波动么?
3楼2013-07-06 03:53:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lucaca

铁杆木虫 (著名写手)


引用回帖:
2楼: Originally posted by shengxiang at 2013-07-05 23:54:53
判断是否平衡的话可以用能量啊,RMSD之类的信息

那么比如说RMSD,算出来的数值为多少的时候就算是平衡了?
4楼2013-07-06 05:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lucaca

铁杆木虫 (著名写手)


引用回帖:
3楼: Originally posted by wbn at 2013-07-06 03:53:40
NPT判断是否平衡难道不是看密度的波动么?

可能是我表述的不清楚。我先做了200ps的NPT,密度最后基本上是趋于稳定,大概是1.00±0.05,但是温度的波动一直在减小,于是我就又做了200ps,密度基本上不变,但温度从600±60K减小到600±40K,如果再往下跑的话密度基本上还是不变的,但温度就不太确定了,所以就问一下,±40k 是否可以算作误差范围内的波动。
5楼2013-07-06 05:21:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbn

银虫 (小有名气)


★ ★ ★ ★ ★
lucaca: 金币+5 2013-07-08 01:58:29
引用回帖:
5楼: Originally posted by lucaca at 2013-07-06 05:21:03
可能是我表述的不清楚。我先做了200ps的NPT,密度最后基本上是趋于稳定,大概是1.00±0.05,但是温度的波动一直在减小,于是我就又做了200ps,密度基本上不变,但温度从600±60K减小到600±40K,如果再往下跑的话密 ...

关注温度的波动这个意义不大吧..因为分子模拟是小体系,而温度是宏观量,体系的温度是由分子的速度和动量算出来的,所以波动一定很大。所以我觉得判断体系是否平衡主要看能量,密度等等是不是收敛,rmsd是不是一条直线,以及测得的各项物理量是不是与实际接近。温度,压力什么的平均值对就行了

» 本帖已获得的红花(最新10朵)

6楼2013-07-06 10:04:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyulinup

金虫 (正式写手)


★
lucaca(金币+1): 谢谢参与
送红花一朵
引用回帖:
6楼: Originally posted by wbn at 2013-07-06 10:04:39
关注温度的波动这个意义不大吧..因为分子模拟是小体系,而温度是宏观量,体系的温度是由分子的速度和动量算出来的,所以波动一定很大。所以我觉得判断体系是否平衡主要看能量,密度等等是不是收敛,rmsd是不是一条 ...

前辈,如何得到能量,密度呢 ?如何用gui 提取STATIS能量数据吧?求指点,感激不尽!!
8楼2014-06-26 15:19:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beeMD

新虫 (初入文坛)


★
lucaca(金币+1): 谢谢参与
如果是做常规MD,数据收集,跑得时间越长越好,ns级的MD基本上不能说是一定平衡了,扩散过程绝大多数不可能是ns级的,还要注意用对了控温控压的方法

BeeMD.org开发团队
9楼2014-06-26 18:15:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖仅楼主可见
10楼2014-08-13 17:00:13
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

CQU_wjq

铁杆木虫 (著名写手)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by wbn at 2013-07-06 03:53:40
NPT判断是否平衡难道不是看密度的波动么?

请问一下NPT模拟的时候压力该如何设置,(我的对象是水溶液);如果用MS的话,在哪里看密度的波动呢,谢谢~
11楼2014-09-11 18:19:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wbn

银虫 (小有名气)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
11楼: Originally posted by CQU_wjq at 2014-09-11 18:19:38
请问一下NPT模拟的时候压力该如何设置,(我的对象是水溶液);如果用MS的话,在哪里看密度的波动呢,谢谢~...

没用过MS,不会,不好意思...
12楼2014-09-12 22:59:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CQU_wjq

铁杆木虫 (著名写手)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
12楼: Originally posted by wbn at 2014-09-12 22:59:00
没用过MS,不会,不好意思......

Thank you all the same
13楼2014-09-13 08:43:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2013-07-07 17:50   回复  
lucaca(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 lucaca 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 04:21 by Equinoxhua
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 16:55 by zzy933268
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:41 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:27 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:14 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:41 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 5/250 2026-11-12 08:15 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 05:00 by 9fOTqj074tIg
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 04:57 by 9fOTqj074tIg
信息提示
请填处理意见