24小时热门版块排行榜    

查看: 10935  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +78 AI8RaGaPaSCR 2026-08-07 82/4100 2026-08-14 10:15 by 5mYCQrDYv1yS
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +14 医学老男孩 2026-08-13 36/1800 2026-08-14 09:41 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见