24小时热门版块排行榜

返回列表

维尼MM

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8.5
帖子: 2
在线: 4.1小时
虫号: 2294265
注册: 2013-02-22
性别: MM
专业: 凝聚态物性I:结构、力学和

[求助] 使用Towhee计算气体吸附遇到的化学势问题及其他

我最近在学习用Towhee做气体吸附，遇到了一些问题，towhee的说明文件又写得非常简略，我自己也想了很久但是百思不得其解，所以非常希望在这里能得到大家的帮助。
我在做的是N2在77K条件下在mesoporous silica MCM-41中的吸附，silica是固定不动的，用的是uvt条件。这是我出现的一些问题：
1) 最后输出的文件里面，average和block average是什么区别呢？我发现两个给出的各项chemical potential值挺不一样，而且都和我最初设定的值有出入，比如我设了－800K的化学势（对应相对气压在0.01), average 给出的结果：
                     units type  average
u (Density)             K 3    -689.409
u (NVT Insertion)       K 3    -896.226
u (NpT Insertion)       K 3    -896.226
u (Den. + NVT Insert)    K 3 -1585.634
u (Den. + NpT Insert)    K 3 -1585.634
u (Gibbs Total)          K 3 -1586.374

blocks average 给出的结果：
                           units type box average    standard deviation
u (Gibbs Total)             K 3 1    -846.397    501.838
u (NpT Insertion)             K 3 1    -157.476    501.718
u (NVT Insertion)             K 3 1    -157.476    501.718
2）block average最后一列的deviation值是不是太大了？我不太清楚怎么判定吸附有没有平衡，我现在的方法是看block average的化学势变化趋势，如果变化缓慢差不多在一个值附近波动，我就概认为平衡了。
3）有时候某些block的化学势会等于0，我觉得不太合理，不太理解为什么会出现这种情况。

回复此楼