24小时热门版块排行榜    

查看: 3257  |  回复: 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 09:45 by xxniao123
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +4 a089 2026-08-26 5/250 2026-08-27 09:24 by Howard28
[基金申请] 基金不中,共勉 +7 eulota 2026-08-26 7/350 2026-08-27 09:02 by hsding2026
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 为什么国自然不能直接公布 +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
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见