| 查看: 1209 | 回复: 2 | ||
[求助]
使用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，我觉得不太合理，不太理解为什么会出现这种情况。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有65人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- towhee中计算临界温度和密度
已经有3人回复
- 关于气液平衡的简单问题
已经有24人回复
- 【讨论】GCMC系综中化学势如何确定
已经有10人回复
2楼2017-05-14 15:57:19
3楼2017-10-26 13:39:41
