写在前面,刚上手一周左右,建立模型后跑起来就是比较机械化的流程了。在这里写一个简明教程算是备忘。如有不足之处欢迎指正。
---------------------------------------------
使用如图所示的sorption模块对所建立的模型进行计算。
TASK选择自己需要的模块,我这里选择的是Fixed pressure。
设置steps & Temperature。
选择method。
简单粗暴的选择依据:简单分子选metropolis,链状分子选configurational bias
选择吸附质sorbates,将浓度转化为压力fugacity并设置。
选择力场forcefield,一般选用COMPASS/COMPASS II。
取消C.A.前面的勾选后,点击calculate和assign,然后勾选C.A.。
设置Total charge tolerance,数值越小,我们后面观察到的精度越高。
设置静电力为Ewald & Group,设置范德华力为Atom based。
注意:红圈里的截断半径不能超过模型所在盒子短边边长的一半
然后就可以run了。
---------------------------------------------
几个小时后,这是successful的样子,吸附质被吸附在模型表面。
这是结束后输出的文件。点击红圈里的文件,可以看到里面有loading和Isosteric heats这两个我们需要的参数。
---------------------------------------------
题外,下面两张是百度到的结果分析技巧。