24小时热门版块排行榜    

查看: 3256  |  回复: 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一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 03:31 by xxniao123
[基金申请] 基金不中,共勉 +5 eulota 2026-08-26 5/250 2026-08-26 23:33 by young90
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统查不到 +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
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见