电荷密度差分图判断氢键求助

yuexiatx

最近计算了一个NH3吸附在TiO2表面羟基上的结构，顺便做了一个电荷密度差分，
请问左边的N-H…O属于氢键吗？
能否用电荷密度差分判断氢键？

sobereva

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“电荷密度差分相关”改了，以后务必注意

密度差通常不用于判断氢键

根据几何关系通常就可以近似判断氢键存在性。
更好、更直观、更严格的做法是用Multiwfn做NCI或IGM或IGMH分析，看
使用Multiwfn结合CP2K通过NCI和IGM方法图形化考察固体和表面的弱相互作用
http://sobereva.com/588（http://bbs.keinsci.com/thread-21742-1-1.html）

虽然Multiwfn目前不支持VASP的波函数，但IGM分析只要提供几何结构就可以计算，Multiwfn可以直接读入POSCAR（用优化后的结构），然后做IGM分析
亦可使用CP2K算个单点得到波函数文件，然后用Multiwfn做效果更好的NCI或IGMH分析。