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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » MS模拟分子动力学NPT系综中压力的设置问题
51/1返回列表
查看: 3318  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

J伯爵

金虫 (小有名气)

[求助] MS模拟分子动力学NPT系综中压力的设置问题 已有1人参与

本人目前利用MS模拟计算聚合物的自由体积，主要是想解释关于气体渗透这一方面，但是目前在分子动力学这一步上有个疑问不知如何解决，就是我所计算的聚合物实测密度约1.4，到是我模拟出来的密度只有1.2几，模拟步骤首先是在1个大气压，298K下，进行1000ps的NPT，然后退火，接着再分别进行500ps的NPT和NVT。开始建了几个不通的模型，发现模拟出来的密度都比实测密度要小很多，然后我就尝试了在进行NPT分子动力学计算时，将压力由1个大气压分别提高到0.1Gpa，1Gpa和10Gpa进行尝试，发现模拟出来的密度终于大概接近了实测密度。现在比较纠结我这样的操作是否正确？NPT中这个压力的设置有没有什么硬性的要求，因为看了很多文献都是采用1个大气压进行模拟的。小弟基本是自学，没人带，现在是一脸懵逼，希望各位高手指导下到底如何操作，谢谢大家！

发自小木虫Android客户端
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-07-22 16:17:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lujunjing

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by J伯爵 at 2016-07-23 09:25:51
现在看别人的论文，他们最终模拟出来的密度和测试密度误差大部分是百分之一到百分之三，我现在算出来的误差有百分之十几，这样的话得最终得到的结果会不会不准确呢？
...

试试别的力场，看看那个接近？
赞一下 回复此楼
4楼2016-07-23 12:51:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

brucefan

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
有可能是你用的势能本身是不那么准确的。经典分子动力学中的势能都是经验性的，其准确性常常是值得怀疑的。你直接和实验比很难判断是你的模拟出了问题，还是所用势能本身就有问题。可以尝试与前人的模拟结果比较，也可以尝试换用别的势能。当然，你的模型也可能影响结果。至于压强，一个大气压和0.1GPa或者干脆零都不会有什么显著的区别，而10GPa显然是非常高的压强了。除非你研究高压情形，否则不要用那么大的压强。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
2楼2016-07-23 00:37:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

J伯爵

金虫 (小有名气)
引用回帖:
2楼: Originally posted by brucefan at 2016-07-23 00:37:50
有可能是你用的势能本身是不那么准确的。经典分子动力学中的势能都是经验性的，其准确性常常是值得怀疑的。你直接和实验比很难判断是你的模拟出了问题，还是所用势能本身就有问题。可以尝试与前人的模拟结果比较，也 ...

现在看别人的论文，他们最终模拟出来的密度和测试密度误差大部分是百分之一到百分之三，我现在算出来的误差有百分之十几，这样的话得最终得到的结果会不会不准确呢？

发自小木虫Android客户端
赞一下 回复此楼
3楼2016-07-23 09:25:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

studentone

新虫 (初入文坛)
您好，我现在进行气体的扩散模拟，在模拟过程中也遇到模拟密度大幅度下降的结果，请问楼主，怎么解决
赞一下 回复此楼
学习
5楼2017-03-02 16:51:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 22:52 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 22:16 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-16 20:50 by MHf40yA0Uhyb
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 20:50 by MHf40yA0Uhyb
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 19:15 by rRvxEj3zorCy
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:15 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:15 by i3KBRD19ERkG
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 16:12 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 14:39 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-16 06:48 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:24 by LH5NkK5Ud4bh
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 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.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭