请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 338|回复: 6

[量化理论] TDDFT的求解过程疑问

[复制链接]

279

帖子

2

威望

1361

eV
积分
1680

Level 5 (御坂)

发表于 2018-9-13 14:10:42 | 显示全部楼层 |阅读模式
我在阅读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及各位知道的老师解惑,谢谢!

1万

帖子

25

威望

1万

eV
积分
35741

管理员

公社社长

发表于 2018-9-13 16:36:30 | 显示全部楼层
1 是
2 是
3 后HF、TDDFT等牵扯到多组态/多行列式的方法,肯定都需要做积分变换,从AO双电子积分变换成MO双电子积分
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

279

帖子

2

威望

1361

eV
积分
1680

Level 5 (御坂)

 楼主| 发表于 2018-9-13 19:27:37 | 显示全部楼层
sobereva 发表于 2018-9-13 16:36
1 是
2 是
3 后HF、TDDFT等牵扯到多组态/多行列式的方法,肯定都需要做积分变换,从AO双电子积分变换成MO ...

明白了~ 谢谢sob老师~

279

帖子

2

威望

1361

eV
积分
1680

Level 5 (御坂)

 楼主| 发表于 2018-9-13 19:45:39 | 显示全部楼层
本帖最后由 让你变成回忆 于 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老师了~

1万

帖子

25

威望

1万

eV
积分
35741

管理员

公社社长

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

X是激发组态系数(对应TDDFT输出里的->),Y是去激发组态系数(对应<-)
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

279

帖子

2

威望

1361

eV
积分
1680

Level 5 (御坂)

 楼主| 发表于 2018-9-14 08:11:26 | 显示全部楼层
sobereva 发表于 2018-9-14 05:31
X是激发组态系数(对应TDDFT输出里的->),Y是去激发组态系数(对应

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

1万

帖子

25

威望

1万

eV
积分
35741

管理员

公社社长

发表于 2018-9-15 05:57:42 | 显示全部楼层
让你变成回忆 发表于 2018-9-14 08:11
嗯,我是这样理解的。  但是从rwf文件里提取出来的数据的个数仍然远多于输出文件里的组态系数,请问sob老 ...

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

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

更多的就是去看Gaussian程序员手册或者源码了
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-11-17 07:19 , Processed in 0.107448 second(s), 23 queries .

快速回复 返回顶部 返回列表