求助从电荷密度角度对不同分子的吸附能的分析

Daisy_zju

论坛的各位老师好，我最近使用vasp计算了对苯二酚和二甲基苯酚两个小分子的吸附能，发现对苯二酚虽然有两个氧原子，理论上吸附作用可能更强，实际计算出的吸附能却高于二甲基苯酚。为了分析原因，我尝试了通过计算bader电荷得到两个分子上的电荷转移结果（如下图），但是仍然无法得到有效的结果，目前正在尝试做吸附末态结构的电荷密度差图，想请教各位老师是否还可以通过何种方法分析是什么作用或者什么基团对吸附的贡献导致两个分子出现了反常的现象？谢谢！

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“从电荷密度角度对吸附能的分析”改了，以后务必注意

这种分子体系极为不适合用VASP，效率很低，功能上限制极大，还没法结合Multiwfn做丰富的波函数分析。应当用Gaussian、ORCA等量子化学程序来算。如果你是指某个固体表面吸附小分子，必须说清楚是什么固体，把截图发出来（最好发二者吸附上去优化后的结构图），否则谁也不知道吸附的机制是什么样。

不要拿Bader电荷（正规名字是AIM电荷）说事，这是非常烂的原子电荷，仔细看
《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）。讨论静电作用时，应当用对范德华表面静电势描述较好的原子电荷，而AIM电荷在这方面是知名的原子电荷里最烂的。

不要什么时候都拿静电相互作用说事，在非极性表面上色散作用绝对不容忽视，甚至色散作用才是主导，此时拿原子电荷讨论毫无意义，参看下文
18碳环（cyclo[18]carbon）与石墨烯的相互作用：基于簇模型的研究一例
http://sobereva.com/615（http://bbs.keinsci.com/thread-25306-1-1.html）
谈谈范德华势以及在Multiwfn中的计算、分析和绘制
http://sobereva.com/551（http://bbs.keinsci.com/thread-17271-1-1.html）