计算化学公社

标题: 用fchk和wfn文件计算的激发态ADCH电荷不同 [打印本页]

作者
Author: 小强 时间: 2017-12-5 22:45
标题: 用fchk和wfn文件计算的激发态ADCH电荷不同
请教各位老师，如题所述, 我用如下方法计算水分子激发态下的ADCH原子电荷,发现用fch和wfn文件得到的结果不同。1) 用#p opt m62x/6-31g* 优化了水分子的基态结构；
2) 用#p opt m62x/6-31g* td density out=wfn 优化了S1态结构,得到的chk文件按照Multiwfn手册操作，转化为fchk文件。
用所得到的fchk文件通过Multiwfn计算得到的ADCH原子电荷如下:
======= Summary of atomic dipole moment corrected (ADC) charges =======
Atom: 1O Corrected charge: -1.298593
Atom: 2H Corrected charge: 0.649286
Atom: 3H Corrected charge: 0.649286

用所得到的wfn文件通过Multiwfn计算得到的ADCH原子电荷如下:
======= Summary of atomic dipole moment corrected (ADC) charges =======
Atom: 1O Corrected charge: -0.123166
Atom: 2H Corrected charge: 0.061584
Atom: 3H Corrected charge: 0.061584
二者差的还挺大的,而且通过fchk文件得到的ADCH电荷与基态下的也不一样,基态ADCH电荷如下:

======= Summary of atomic dipole moment corrected (ADC) charges =======
Atom: 1O Corrected charge: -0.755662
Atom: 2H Corrected charge: 0.377808
Atom: 3H Corrected charge: 0.377808

请各位指正下问题出在哪了，而且不明白此时通过fchk文件得到的ADCH电荷是对应与什么状态下的，还是没有意义呢？先谢谢各位了!

作者
Author: sobereva 时间: 2017-12-5 23:31
注意看
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379
只有当你把TDDFT激发态自然轨道转存到chk文件里，基于fch计算的结果才是对应激发态的，否则算的是相应结构下基态的。

作者
Author: 小强 时间: 2017-12-6 10:58

sobereva 发表于 2017-12-5 23:31
注意看
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379

谢谢老师回复。帖子看过了，按照帖子中提示的操作保存fchk文件，对于垂直激发的激发态结构的ADCH电荷计算，fchk文件与wfn文件得到的结果一致。但是当我用优化后的激发态结构计算时，这两文件得到的结果还是不一样。请问老师, 对于优化的激发态结构，fchk文件处理有什么需要注意的地方吗？

作者
Author: sobereva 时间: 2017-12-6 13:15

小强发表于 2017-12-6 10:58
谢谢老师回复。帖子看过了，按照帖子中提示的操作保存fchk文件，对于垂直激发的激发态结构的ADCH电荷计算 ...

看不懂你的描述。

“垂直激发的激发态结构的ADCH电荷计算，fchk文件与wfn文件得到的结果一致”不可能一致。你用TD还out=wfn，则wfn里是激发态自然轨道（对较新版本来说，此时默认用density），而此时fch里还是基态DFT轨道，结果显然不同。

作者
Author: 小强 时间: 2017-12-6 14:18

sobereva 发表于 2017-12-6 13:15
看不懂你的描述。

“垂直激发的激发态结构的ADCH电荷计算，fchk文件与wfn文件得到的结果一致”不可 ...

首先感谢老师的耐心指导。我的问题简单来说就是: 我想计算激发态极小点结构的ADCH电荷；fchk文件是通过您给的博文里提到的方法得到的(此时,fchk文件记录的应该也是激发态轨道呢），然而通过该文件得到的ADCH电荷与wfn文件得到的ADCH电荷还是不一样。(ps: 之所以提到垂直激发的激发态结构的ADCH计算，并与wfn文件的结果进行比较，我是想验证转化fchk文件的操作没有出错。) 但愿表述清楚了，希望得到老师的指导。

作者
Author: sobereva 时间: 2017-12-6 17:54

小强发表于 2017-12-6 14:18
首先感谢老师的耐心指导。我的问题简单来说就是: 我想计算激发态极小点结构的ADCH电荷；fchk文件是通过您 ...

分别载入fch和wfn，在主功能0里选择显示所有轨道，看显示的轨道占据数。如果二者不同，肯定是激发态自然轨道没有存到fch里，或者没有写saveNO

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