24小时热门版块排行榜    

查看: 3269  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:41 by hZiFeudZyoGR
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:40 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:46 by hZiFeudZyoGR
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 jCd0dEvKHShX 2026-08-29 7/350 2026-08-30 19:45 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:13 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 17:47 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 17:45 by hZiFeudZyoGR
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 17:34 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 17:25 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:14 by hZiFeudZyoGR
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
信息提示
请填处理意见