咨询文献中新指数的计算，谢谢指点。

小范范1989

如图，文献中的J的计算方式，但是我感觉这样计算是不是J为0吗？我的理解一定错了。
再就是为什么要-E（Homo）（N）必须要等于IP？
谢谢指点。

sobereva

对于精确的泛函，HOMO的能量和电离能精确相等。但对于我们实际用的泛函，显然不满足这个条件，于是让J（HOMO的能量和电离能的差值）最小化，就能优化参数。

小范范1989

sobereva 发表于 2015-8-29 19:07
对于精确的泛函，HOMO的能量和电离能精确相等。但对于我们实际用的泛函，显然不满足这个条件，于是让J（HOM ...

偶，这样子，谢谢sob老师的指点

小范范1989

sobereva 发表于 2015-8-29 19:07
对于精确的泛函，HOMO的能量和电离能精确相等。但对于我们实际用的泛函，显然不满足这个条件，于是让J（HOM ...

sob老师好，回顾之前的文献，我又不太理解一个地方。所以又来麻烦老师您。这个地方，和我原来咨询的这个J怎么计算不一样了？而且文献中所说：N体系的电子亲和能=N+1体系的电离能？按说这两个应该不一样吧？
没太看懂，谢谢老师指点。

ps：1：文章中说Tamm-Dancoff approximation （TDA）这和我的通常用的TD-DFT有什么区别吗？比如说，我用TD-DFT的方法计算，我可以写为b3lyp/6-31g* td....那要是用TDA的方式呢？泛函和基组不也是这样写？如果是这样的话，高斯09是怎样区别这两种方式的？换句话说，我怎么样确定采用TDA的方式？
谢谢sob老师的指点。thanks

sobereva

小范范1989 发表于 2015-9-28 08:16
sob老师好，回顾之前的文献，我又不太理解一个地方。所以又来麻烦老师您。这个地方，和我原来咨询的这个J ...

J2就是计算[E_HOMO(N)+IP(N)]^2 + [E_HOMO(N+1)+IP(N+1)]^2
原先计算J的时候只考虑了N没有考虑N+1态，而且没有像这样取平方而只是取绝对值。

N+1的电离能一般叫做electron detachment energy，和N的电子亲和能一样表达式都是E(N)-E(N+1)。如果是绝热的，那么二者等价，如果是垂直的，那么前者用的是N+1平衡结构的，后者用的是N电子平衡结构的，因此会不同。


把TD关键词改为TDA就行了，其它都不变。

小范范1989

sobereva 发表于 2015-9-28 22:36
J2就是计算[E_HOMO(N)+IP(N)]^2 + [E_HOMO(N+1)+IP(N+1)]^2
原先计算J的时候只考虑了N没有考虑N+1态，而 ...

偶，明白了，谢谢sob老师。

youshoubian2015

我想追问一下，TDA-DFT是不是和TDA是相同的？

小范范1989

youshoubian2015 发表于 2015-11-27 22:15
我想追问一下，TDA-DFT是不是和TDA是相同的？

1：我记得昨天我看了一个文献：Tamm-Dancoff approximation （TDA）忽略了B矩阵也不知道啥的，就是当时sob老师给的那篇初始文献：JCP, 115, 3540。有介绍，我要是没记错的话。
2：TDA好像是比TD快，但是，需要高斯新的版本才支持TDA。
3：具体的使用方法就是和sob老师说的这样：把TD关键词改为TDA就行了，其它都不变

sobereva

youshoubian2015 发表于 2015-11-27 22:15
我想追问一下，TDA-DFT是不是和TDA是相同的？

TDA是指Tamm-Dancoff approximation，基于它来做TDDFT（Time-dependent DFT），全称应当写为TDA-TDDFT。为了省事就简写为了TDA-DFT，甚至直接简写成TDA。

youshoubian2015

晓得了，谢谢啊！