乱谈激发态的计算方法

小范范1989

sobereva 发表于 2015-1-14 03:10
问题没有那么复杂。激发态势能面形状对泛函的敏感度也明显没有激发能那么高，不必考虑得过于精细。
简单 ...

谢谢sob的指点，也就是说：
1：一般，优化单重激发态使用的是b3lyp  或者是pbe0 。优化三重激发态用的是m062x
2：如果说，想不区分，都用一个来做，就用m062x。
3：如果说更细分：激发态是CT态，需要更多的HF来描述比较准确，所以使用m062x。
                 激发态是LE态，就用b3lyp  pbe0。
4：总的来说，泛函的选择，需要根据计算的分子的特点来选定，没有说是一成不变的，但是，上面是我们总结的经验，一般情况是这样的，
不知道以上我的理解是否有错误，谢谢sob指点
ps：我贴出的那个jacs文章，激发态用的什么方法？

谢谢sob

ter20

Sob可否在提到GGA时把LDA和GGA型的泛函罗列一下，或者提示相关信息的来源，现在还是有些不能完全记清楚，百度或google到的感觉都不全，汗

sobereva

ter20 发表于 2015-3-10 16:03
Sob可否在提到GGA时把LDA和GGA型的泛函罗列一下，或者提示相关信息的来源，现在还是有些不能完全记清楚，百 ...

GGA或meta-GGA有BLYP、TPSS、 BPW91、BP86、PW91PW91、PBE、OLYP、mPW1PW91、M06L、XLYP、VSXC、B95、SSB、KT1/2/3、HCTH等等

LDA就那么几种，交换部分是Slater，相关部分是VWN3、VWN5或PW。

sobereva

小范范1989 发表于 2015-1-14 13:21
谢谢sob的指点，也就是说：
1：一般，优化单重激发态使用的是b3lyp  或者是pbe0 。优化三重激发态用的是 ...

基本就是这样。
那个文章属于非主流做法，自己调的。

sobereva

小范范1989 发表于 2015-1-14 13:21
谢谢sob的指点，也就是说：
1：一般，优化单重激发态使用的是b3lyp  或者是pbe0 。优化三重激发态用的是 ...

再补充一下，虽然小分子单重态价层激发适合用你说的PBE0，但对于有机大共轭分子，建议用HF成份高一些的，M062X是个不错的选择。GGA以及HF成分低的情况，包括B3LYP，对有机大共轭分子容易出现ghost态。

小范范1989

sobereva 发表于 2015-3-10 22:21
再补充一下，虽然小分子单重态价层激发适合用你说的PBE0，但对于有机大共轭分子，建议用HF成份高一些的， ...

嗯，谢谢sob的指点。very much

Gratitude_zy

谢谢sob老师！

mengyingbian

sob大神，急需求助！我是高斯计算菜鸟，已经被nto计算折磨的身心俱疲了，我进行了nto计算，终于计算过程中没有报错，但是不知道怎么证明已经完成了nto计算，还有当我用fchk文件看nto轨道的时候我应该点完cube以后选哪种类型啊，我可以选的类型只有：molecular orbital；total density；alpha density；ESP。
以下是我做的三次计算的输入文件：
1
%nprocshared=3
%mem=200MW
%chk=E:\g9new\bmy\tpesinglet1energyNTO2.chk
# td=(singlets,nstates=10,root=2)/6-31+g(d,p) geom=connectivity density=transition=2 m062x
2
%nprocshared=3
%mem=200MW
%chk=E:\g9new\bmy\tpesinglet1energyNTO2.chk
# td=(singlets,nstates=10,root=2)/6-31+g(d,p) pop=nto density=transition=2 geom=check guess=read m062x
3
%nprocshared=3
%mem=200MW
%chk=E:\g9new\bmy\tpesinglet1energyNTO2.chk
# density=(check,transition=2) geom=allcheck guess=(read,only) pop=savento

求大神指教！

sobereva

mengyingbian 发表于 2016-3-13 10:16
sob大神，急需求助！我是高斯计算菜鸟，已经被nto计算折磨的身心俱疲了，我进行了nto计算，终于计算过程中 ...

为何搞这么复杂，直接按照此文的做法弄就完了：
跃迁密度分析方法-自然跃迁轨道(NTO)简介
http://sobereva.com/91
要算S0->S2的NTO，就写
#P m062x/6-31g(d,p) td=nstates=10 density=transition=2 pop(NTO,saveNTO)
然后用Multiwfn载入fch文件，进主功能0就能看到NTO了。gview的话，直接在轨道观看界面，看到的也是NTO。

mengyingbian

太感动了，竟然这么快就得到了回复！！我已经折腾了很久了。。。
直接用gview在轨道观看界面是载入new cube的，然后会有很多选项，比如看homo lumo什么的，不知道选哪个，而且输出文件中关于激发态的概率算完nto后会不会也有专门关于它的概率

sobereva

mengyingbian 发表于 2016-3-13 10:26
太感动了，竟然这么快就得到了回复！！我已经折腾了很久了。。。
直接用gview在轨道观看界面是载入new cub ...

你那个是观看等值面界面。我说的是Edit-MOs那个界面，就和观看普通分子轨道一样，因为用了saveNTO后NTO已经代替MO存到了chk里了。

没有什么所谓的激发态概率，怎么解释NTO对激发的贡献在前面提到的博文里已经说了。

mengyingbian

我的edit mos界面是附件中的样子

mengyingbian

sobereva 发表于 2016-3-13 10:30
你那个是观看等值面界面。我说的是Edit-MOs那个界面，就和观看普通分子轨道一样，因为用了saveNTO后NTO已 ...

是不是我算的是错的，下面是您给的例子
Alpha  occ. eigenvalues --    0.00000   0.00000   0.00000   0.00000   0.00000
Alpha  occ. eigenvalues --    0.00000   0.00000   0.00001   0.00003   0.00005
Alpha  occ. eigenvalues --    0.00008   0.00018   0.00077   0.00611   0.99276
Alpha virt. eigenvalues --    0.99276   0.00611   0.00077   0.00018   0.00008
Alpha virt. eigenvalues --    0.00005   0.00003   0.00001   0.00000   0.00000
Alpha virt. eigenvalues --    0.00000   0.00000   0.00000   0.00000   0.00000
但是我算完后occ.部分是负值。。。。

sobereva

mengyingbian 发表于 2016-3-13 10:50
是不是我算的是错的，下面是您给的例子
Alpha  occ. eigenvalues --    0.00000   0.00000   0.00000   ...

不代表是错的。
你先把我文中的例子重复一遍在搞你的体系。

mengyingbian

sobereva 发表于 2016-3-13 12:56
不代表是错的。
你先把我文中的例子重复一遍在搞你的体系。

好的 谢谢大神！！遇到问题还得请教您！！！