| 查看: 2194 | 回复: 1 | ||||
[求助]
分子模拟中压力和密度的差异问题
|
|
我们都知道水的密度在常温下大概是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）所示。  NVT_pressure.jpg 此外对所得水分子盒子作径向分布分析，如图（NVT_grwater.jpg）所示，怎么感觉第一个峰（2A左右）处的值有点偏小，那个峰的值是不是应该在2左右才对啊？  NVT_grwater.jpg 此外，我也用NPT模式计算过，设定温度298K和压力0.0001 GPa时，水盒子的密度会最终稳定在0.5 g/mL左右的，而将压力提升至0.45 GPa时，用NPT跑水盒子的密度才会维持在1 g/mL，这个密度和压力的关系和我们的常识感觉有比较大的差异啊。 最后谢谢大家关注！ |
回复此楼» 收录本帖的淘帖专辑推荐
 MS计算模拟 | 分子模拟 |
» 猜你喜欢
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
明天放榜？
已经有4人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于正压力的求解的实现问题
已经有4人回复
- compass力场模拟苯酚密度
已经有3人回复
- 分子模拟中MSD的问题讨论
已经有25人回复
- 模拟中的密度计算
已经有3人回复
- 如何计算溶剂盒子密度
已经有15人回复
- 周期性和非周期性态密度的区别
已经有3人回复
- MS里的gulp能计算合金中某一种元素的声子态密度吗？
已经有10人回复
- 跑动力学密度越跑越大了？？？？？
已经有8人回复
- 用什么软件把三张态密度图画在一个图中？
已经有11人回复
- Discover dyanimics 模拟醇的水溶液时的氢键能为什么总是零？
已经有7人回复
- DPD 关于粒子数密度的取值
已经有3人回复
- 采用MS分子模拟计算双马来酰亚胺单体的内聚能密度和溶解度参数
已经有5人回复
- MS-dmol3电荷密度分析
已经有14人回复
- 请教分子力学与分子动力学的区别
已经有14人回复
- 【求助】溶剂盒子的密度？
已经有5人回复
- 【求助】计算固体热膨胀系数
已经有8人回复
- 【求助】Discover模块下的压力单位GPa和平常的压力单位一样吗？
已经有3人回复
pigrass
专家顾问 (著名写手)
伊卡斯特的小猪
-
专家经验: +503 - 模拟EPI: 4
- 应助: 377 (硕士)
- 金币: 10372
- 散金: 3286
- 红花: 167
- 帖子: 2101
- 在线: 1968.4小时
- 虫号: 1463681
- 注册: 2011-10-27
- 性别: GG
- 专业: 高分子材料结构与性能
- 管辖: 分子模拟
★
月只蓝: 金币+1, 感谢指导！ 2014-07-22 18:49:48
月只蓝: 金币+1, 感谢指导！ 2014-07-22 18:49:48
|
我不清楚compass力场对纯水的适用性 这个你最好找找相关文献证明水是可以用compass做的 论坛也有人反映过compass计算水出错 http://muchong.com/html/201205/4458625.html |
2楼2014-07-21 10:58:47
