求助：vasp计算怎么得到吸附过程中的ΔG和熵值

yuzhang

  各位老师好，我用VASP做的水中碳纳米管的吸附，比较吸附强度大小的主要参数就是吸附能（∆Eads=Etotal−ECNTs−Eadsorbent）。审稿人提出这样的建议“To enhance thermodynamic rigor, the analysis should incorporate GibbsFree Energy (ΔG) rather than relying solely on internal energy to account forentropic penalties.”
      我不太明白他为什么要我考虑熵，我的论文里没有算能变和能垒，而且也不好实现。我想请问回复审稿意见能不能通过引用文献说明一下吸附能已经可以满足。如果不行的话VASP计算怎么能得到熵吗？
      麻烦各位老师答疑。

UW_0728.

计算热力学量需要做振动分析，VASP也不例外。可以简单参考一下VASP Wiki上做振动分析的例子
比较关键的词是IBRION。如果要做解析Hessian的振动分析设IBRION=7或8（8考虑对称性）；要做基于有限差分计算受力得到Hessian的则IBRION设5或6（6考虑），同时需要设定NFREE（一般为2，此时会在±x、±y、±z六个位移下计算受力，相当于N个原子下做6N次受力计算）、进一步设严EDIFF（比如1e-7甚至1e-8）等