24小时热门版块排行榜    

查看: 3583  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

mengyan142

木虫 (小有名气)

[求助] 如何计算溶剂盒子密度

现在我手头上有一个大分子,然后在它的周围加上了10A的水分子作为溶剂盒子,请问我该怎么去计算这个盒子的密度呢?
而且我看到有文献里面是这样提到的:“A periodic solvent box was created and equilibrated by MD simulations. The calculated densities of 0.692 g/cm3 (heptane) and 1.004 g/cm3 (( R)-1-phenyl ethanol) of the box in equilibrium were in good agreement with the experimental values, 0.684 g/cm3 and 1.012 g/cm3".
就是想知道为什么我那个图片上的是0.8,与水的实际密度1.0差这么多,而文献中却几乎一样?

rongjihezi.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengyan142

木虫 (小有名气)

送红花一朵
引用回帖:
13楼: Originally posted by jiaoyixiong at 2013-04-22 15:52:52
你需要看看你的 NPT 过程中盒子的边长变化情况是否最后趋向一个定值啊,文献中的密度值肯定是有时间平均的,如果你仅仅取NPT结果最后的那个值与实验结果对比,肯定会差异比较大的,后续一段NVT模拟,然后统计平均N ...

您好,NPT 过程中盒子的边长变化情况在哪个文件看呢?我在输出文件里只看到了体积的变化,我把out文件上传了,闲了您帮我看看吧,
还有个问题,这里的密度只是溶剂密度么?还是包括了蛋白质,最后平均下来的密度呢?
谢谢您了....

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : npt.out
  • 2013-04-23 15:40:17, 31.81 M
14楼2013-04-23 15:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

mengyan142

木虫 (小有名气)

上面那个图不对,截错了,应该为:solvatebox myy TIP3PBOX 10
                                                solute vdw bounding box :                 49.327  41.467  53.940
                                                total bounding box for atom centers:   69.327  61.467 73.940
                                                solvent unit box :                             18.774  18.774  18.774
                                                total vdw box size :                           72.535  64.737  77.076
                                                volume :361923.915 A3
                                                total mass 174342.130 amu, density 0.800 g/cc
                                                added 8290 residues
2楼2013-04-19 20:08:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

death_snake

金虫 (初入文坛)

同学,这个貌似在你的mdcrd文件中应该可以看到的。应该是1没错,至于它这个是怎么算的。。我也不是特别清楚哎。。
sky king cover the land tiger
3楼2013-04-19 20:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

death_snake

金虫 (初入文坛)


ben_ladeng: 金币+1, 鼓励交流 2013-04-22 21:27:30
引用回帖:
3楼: Originally posted by death_snake at 2013-04-19 20:19:11
同学,这个貌似在你的mdcrd文件中应该可以看到的。应该是1没错,至于它这个是怎么算的。。我也不是特别清楚哎。。

不对,说错了,是Out文件。。。
sky king cover the land tiger
4楼2013-04-19 20:22:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见