24小时热门版块排行榜    

查看: 2239  |  回复: 1

tzweir

金虫 (小有名气)

[求助] 分子模拟中压力和密度的差异问题

我们都知道水的密度在常温下大概是1左右,然而我用分子动力学模拟了一下却发现与这个常识有较大的差异。
使用软件:Material Studio 5.5,Forcite模块
1、水分子:在MS中画出水分子后Clean一下,然后电荷用QEq方法计算赋值
2、水盒子构成:采用AC模块,432个水分子构成一个盒子,设定密度为1,力场为COMPASS
3、结构优化:采用Forcite模块的Geo Opt优化盒子的结构
4、分子模拟Forcite Dynamic模块
系综NVT,温度:298K, Time step 1.0,设定模拟200 ps,Thermostat为Nose,但是整个体系在运行到70ps左右时,软件提示died unexpectedlly,如图(NVT_energy.jpg)能量曲线所示。
分子模拟中压力和密度的差异问题
NVT_energy.jpg
照理说现在水的密度为1,在298k下,他的压力应该是0.1 MPa(一个大气压)才对,但是查看Pressure发现此时体系内的压力大概有0.4 GPa,如图(NVT_pressure.jpg)所示。
分子模拟中压力和密度的差异问题-1
NVT_pressure.jpg
此外对所得水分子盒子作径向分布分析,如图(NVT_grwater.jpg)所示,怎么感觉第一个峰(2A左右)处的值有点偏小,那个峰的值是不是应该在2左右才对啊?
分子模拟中压力和密度的差异问题-2
NVT_grwater.jpg

此外,我也用NPT模式计算过,设定温度298K和压力0.0001 GPa时,水盒子的密度会最终稳定在0.5 g/mL左右的,而将压力提升至0.45 GPa时,用NPT跑水盒子的密度才会维持在1 g/mL,这个密度和压力的关系和我们的常识感觉有比较大的差异啊。

最后谢谢大家关注!
回复此楼

» 收录本帖的淘帖专辑推荐

MS计算模拟 分子模拟

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

多为成功找方法,少为失败找理由。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pigrass

专家顾问 (著名写手)

伊卡斯特的小猪

★
月只蓝: 金币+1, 感谢指导! 2014-07-22 18:49:48
我不清楚compass力场对纯水的适用性 这个你最好找找相关文献证明水是可以用compass做的

论坛也有人反映过compass计算水出错 http://muchong.com/html/201205/4458625.html
#有问题发论坛,不接受私聊#爱晶体,爱MaterialsStudio#MS专辑点击右边#http://emuch.net/taotie/41673.html#物性数据#http://webbook.nist.gov/chemistry/#
2楼2014-07-21 10:58:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员leB0zQ 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见