24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3516  |  回复: 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的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

引用回帖:
11楼: Originally posted by mengyan142 at 2013-04-22 15:34:26
我就是用的最后的那个密度值,不过看文献里他的密度跟实验值差距特别小,我的还是比较大,是不在续一个NVT就会更好些呢?不过接着跑一段怕时间来不及,所以直接用npt的密度是不也可以?...

你需要看看你的 NPT 过程中盒子的边长变化情况是否最后趋向一个定值啊,文献中的密度值肯定是有时间平均的,如果你仅仅取NPT结果最后的那个值与实验结果对比,肯定会差异比较大的,后续一段NVT模拟,然后统计平均NVT结果的密度值再与文献对比,可能会好一些。

» 本帖已获得的红花(最新10朵)

13楼2013-04-22 15:52:52
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见