24小时热门版块排行榜    

查看: 11058  |  回复: 4

jiaozihao

新虫 (小有名气)

[求助] ms分子动力学模拟NVT及NPT应用问题? 已有1人参与

大家好!
       我最近在做水的分子动力学模拟,网上文献在对模拟结果性质进行分析时,会提到具体温度T和压力P。但是在实际模拟过程中,用NPT系综是很难准确达到我们所设定的压力,P一般都会在设定值附近涨落,那么我们是不是就可以认为此时(虽然压力有波动但能量曲线已经相对稳定时)的性质就是对应的我们所设定的T和P下的性质呢?
       还有一个问题就是做分子动力学模拟时,如果我们要使体系达到需要的密度和温度状态,是不是只要NVT一个系综就行了?或者是NVT+NVE(因为一旦有NPT加入,密度就和初始设定值不一样了)?但是这样的话最终的压力波动会很大
       如果我们要要使体系达到需要的温度和压力状态,是用NPT算平衡再用NVT计算性质对吗?
我是新书,上面这些只是我自己根据所学的一些猜测,希望大家不吝赐教!!!
衷心感谢!!!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jiaozihao: 金币+50, ★★★★★最佳答案 2018-03-09 00:02:32
总的评论:你的问题不局限于任何MD软件。

我最近在做水的分子动力学模拟,网上文献在对模拟结果性质进行分析时,会提到具体温度T和压力P。
但是在实际模拟过程中,用NPT系综是很难准确达到我们所设定的压力,P一般都会在设定值附近涨落,
那么我们是不是就可以认为此时(虽然压力有波动但能量曲线已经相对稳定时)的性质就是对应的我们所设定的T和P下的性质呢?

回答:当瞬时压强在目标压强上下涨落,而且瞬时压强的平均值等于目标压强时,就可以认为系统的压强等于目标压强。
温度得情形也类似。这是MD模拟的最基本假设,即任何物理量的测量值等于其时间平均值。


还有一个问题就是做分子动力学模拟时,如果我们要使体系达到需要的密度和温度状态,是不是只要NVT一个系综就行了?
或者是NVT+NVE(因为一旦有NPT加入,密度就和初始设定值不一样了)?但是这样的话最终的压力波动会很大

回答,如果事先就知道要达到怎样的密度,那么事先就会知道盒子的大小(如果原子数确定了的话)。此时,用NVT系综
模拟即可。没有NVT+NVE这样的说法。不要管压力波动大还是小,因为它肯定是有波动的,而且它的波动与体系粒子数
有很大关系。越大的体系波动越小。

如果我们要要使体系达到需要的温度和压力状态,是用NPT算平衡再用NVT计算性质对吗?

回答:如果要达到特定的温度和压强,就得用NPT系综模拟。当达到目标状态(平衡态)后,用什么系综
模拟来计算性质,取决于研究目标。
2楼2018-03-08 23:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaozihao

新虫 (小有名气)

嗯嗯,我懂了!真的是非常感谢您!!谢谢^ω^

发自小木虫Android客户端
3楼2018-03-09 00:04:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Quire_z

新虫 (初入文坛)

LZ你好,窝做的是液氮的分子动力学模拟,想问一下你用的MS软件是什么版本
4楼2018-03-12 21:04:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaozihao

新虫 (小有名气)

引用回帖:
4楼: Originally posted by Quire_z at 2018-03-12 21:04:07
LZ你好,窝做的是液氮的分子动力学模拟,想问一下你用的MS软件是什么版本

好像是2016

发自小木虫Android客户端
5楼2018-03-16 00:32:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiaozihao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:23 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见