24小时热门版块排行榜    

查看: 3267  |  回复: 4

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客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

lujunjing

新虫 (初入文坛)

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

试试别的力场,看看那个接近?
4楼2016-07-23 12:51:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

studentone

新虫 (初入文坛)

您好,我现在进行气体的扩散模拟,在模拟过程中也遇到模拟密度大幅度下降的结果,请问楼主,怎么解决
学习
5楼2017-03-02 16:51:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员jEOp3v 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:43 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:51 by ZPa0EcMwuECS
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 03:26 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-30 01:58 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 01:14 by ZPa0EcMwuECS
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见