利用gaussian进行NTO计算的相关问题

seantan521

大家好，事情是这样的，今天听同事说了一种计算NTO的方法，结合公社上其他人的帖子，有一些疑问：

（方法一）同事的方法：

优化完S1和T1后，先做一次energy计算：# td=(singlet，nstates=3，root=1） b3lyp/6-311g(d)

然后再这样算一次：# td=(singlet，nstates=3，root1） b3lyp/6-311g(d)  pop=(nto,minimal,savento) density=(check,transition=1) geom=connectivity chkbas

同事的意思是说，因为S1和T1优化完之后，log文件里有很多组excited state 1，pop后的key words无法识别应该提取哪一组excited state 1，所以她又算一次energy。

（方法二）公社帖子里的方法：

直接进行以下计算：# td=(singlet，nstates=3，root1） b3lyp/6-311g(d)  pop(saveNTO,NTO) density=transition=1


问题一：怎么评价同事的这种方法？
问题二：是不是如果之前的优化用了一种基组，比如说b3lyp/6-311g(d)，这里计算NTO的时候最好也用同一种呢？
问题三：方法二里，到底应该是基于优化后的S0还是优化后的S1和T1做计算啊？

感谢大家解惑答疑

sobereva

按照此文方法做，是绝对正确的，别管别的
使用Multiwfn做自然跃迁轨道(NTO)分析
http://sobereva.com/377

用S0的结构做相当于分析垂直激发的情况，用S1、T1优化的结构分别分析S0-S1和S0-T1的跃迁相当于分析荧光和磷光发射的情况。

电子激发计算不要求和几何优化用相同级别。

让你变成回忆

sobereva 发表于 2018-4-6 00:37
按照此文方法做，是绝对正确的，别管别的
使用Multiwfn做自然跃迁轨道(NTO)分析
http://sobereva.com/377 ...

完全赞同sob老师的推荐。Multiwfn的批量做NTO的确方便~