计算化学公社

标题: TDDFT的求解过程疑问 [打印本页]

作者
Author: 让你变成回忆 时间: 2018-9-13 14:10
标题: TDDFT的求解过程疑问
我在阅读Gaussian做TDDFT计算的输出文件发现，在NFC上面有这样一行：Range of M.O.s used for correlation: 2 35，而此时的NFC数值为1。而其他的TDDFT计算输出文件也有类似规律：即Range of M.O.s used for correlation后的一个数总是NFC后面的数字加1。为此，我有如下几个疑问：

（1）NFC的全称是否是：Number of Freeze Core？意思是对内层的轨道不考虑？
（2）在TDDFT计算中，首先是做基态计算得到分子轨道，然后根据占据轨道和空轨道构造组态函数，那么这里的意思是否是在考虑占据轨道的时候没有考虑最内层的轨道？
（3）看了下这是Gaussian的L801做的计算，我查阅手册发现L801的功能是：Initializes transformation of 2-electron integrals，意思是否是对AO积分->MO积分的变换？那与这里TD的计算有什么关系呢？

恳请sob及各位知道的老师解惑，谢谢！

作者
Author: sobereva 时间: 2018-9-13 16:36
1 是
2 是
3 后HF、TDDFT等牵扯到多组态/多行列式的方法，肯定都需要做积分变换，从AO双电子积分变换成MO双电子积分

作者
Author: 让你变成回忆 时间: 2018-9-13 19:27

sobereva 发表于 2018-9-13 16:36
1 是
2 是
3 后HF、TDDFT等牵扯到多组态/多行列式的方法，肯定都需要做积分变换，从AO双电子积分变换成MO ...

明白了~ 谢谢sob老师~

作者
Author: 让你变成回忆 时间: 2018-9-13 19:45
本帖最后由让你变成回忆于 2018-9-13 19:47 编辑

sobereva 发表于 2018-9-13 16:36
1 是
2 是
3 后HF、TDDFT等牵扯到多组态/多行列式的方法，肯定都需要做积分变换，从AO双电子积分变换成MO ...

sob老师您好，这里再向您请教一个问题：
最近我在看pysoc的代码，其中一个是从rwf文件从提取出CIS的系数（对应于635R槽位），但是提取出来的数据要比log文件里输出的多得多。
为了能够很好的对比，我把log文件中输出的阈值改为了IOP(9/40=400)，然后在激发态信息输出部分，所有的占据轨道->空轨道的组态均能输出（包括所有的去激发），统计个数发现为：占据轨道数目（不包括内层轨道）*空轨道数目*2（说明IOP设置得足够小），即便我乘上计算的激发态的总数目，仍然比从rwf文件中提出出来的数据个数少得多？我看pysoc的源代码中是利用XY这样的词语来描述的，然后从Gaussian官网看到有<X+Y|X-Y>=1这样的条件。
所以想请教您：rwf中的X、Y与log文件中的跃迁组态系数之间的关系是什么？
麻烦sob老师了~

作者
Author: sobereva 时间: 2018-9-14 05:31

让你变成回忆发表于 2018-9-13 19:45
sob老师您好，这里再向您请教一个问题：
最近我在看pysoc的代码，其中一个是从rwf文件从提取出CIS的系数 ...

X是激发组态系数（对应TDDFT输出里的->），Y是去激发组态系数（对应<-）

作者
Author: 让你变成回忆 时间: 2018-9-14 08:11

sobereva 发表于 2018-9-14 05:31
X是激发组态系数（对应TDDFT输出里的->），Y是去激发组态系数（对应

嗯，我是这样理解的。但是从rwf文件里提取出来的数据的个数仍然远多于输出文件里的组态系数，请问sob老师，您知道这是为啥吗？我想彻底了解从rwf提取出的信息，具体可以看看哪些资料呢？

作者
Author: sobereva 时间: 2018-9-15 05:57

让你变成回忆发表于 2018-9-14 08:11
嗯，我是这样理解的。但是从rwf文件里提取出来的数据的个数仍然远多于输出文件里的组态系数，请问sob老 ...

找个很小的体系用很小基组计算，对比一下输出应该能判断出来原因

这个没看过的话可以看看
提取Gaussian的rwf文件信息的工具rwfdump简介
http://sobereva.com/112

更多的就是去看Gaussian程序员手册或者源码了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)