计算化学公社

标题: ECDA计算片段A对HOMO的贡献与体系的HOMO在各个片段分布比例相差很大 [打印本页]

作者
Author: dlut211 时间: 2019-5-21 17:49
标题: ECDA计算片段A对HOMO的贡献与体系的HOMO在各个片段分布比例相差很大
老师，您好。1 我们采用Multiwfn的ECDA模块计算了片段A对HOMO的贡献，发现总贡献约为40%。

2 我们做了HOMO的分布，从HOMO分布图看，很少的分布在片段A上。图形比例约10%。
想问老师，怎么解释这两者差距很大。
计算细节：
G16 TDDFT cc-pVDZ

谢谢老师。

作者
Author: sobereva 时间: 2019-5-22 03:04
那个是CDA模块，不是ECDA模块。ECDA只是个给出片段间电荷转移量的分析方法，是CDA模块会顺带做的分析。

要算轨道成份，用下文的做法，用不着CDA模块
谈谈轨道成份的计算方法
http://sobereva.com/131

仅当你要考察片段间电荷转移细节、考察片段轨道是怎么组合成复合物轨道、绘制轨道相互作用图的时候才用得着CDA模块。

作者
Author: dlut211 时间: 2019-5-22 11:58

sobereva 发表于 2019-5-22 03:04
那个是CDA模块，不是ECDA模块。ECDA只是个给出片段间电荷转移量的分析方法，是CDA模块会顺带做的分析。

...

谢谢老师的回复，我没有描述清楚问题，我们是考察片段轨道是怎么组合成复合物轨道。片段1-H2O对复合物HOMO贡献为37.25%，而HOMO图中1-H2O所占比例太小。想问下老师这两者是否有可比性，怎么来解释这个差异。
(, 下载次数 Times of downloads: 39)

麻烦老师了。

作者
Author: sobereva 时间: 2019-5-22 12:31
看HOMO其余部分都是什么轨道构成的
注意序号别弄错了，不能带弥散函数

作者
Author: dlut211 时间: 2019-5-22 12:53

sobereva 发表于 2019-5-22 12:31
看HOMO其余部分都是什么轨道构成的
注意序号别弄错了，不能带弥散函数

谢谢sobereva老师。
1  HOMO其余部分主要分布在片段CO2上，约为43.52%。但从图形上看，绝大部分HOMO分布在片段CO2。
2  我们用的基组为cc-pVDZ (相当于6-31G**级别),HOMO作图也是该基组计算得到的fchk。
3  原子序号检查了，分析时没有出错。

作者
Author: sobereva 时间: 2019-5-22 14:00

dlut211 发表于 2019-5-22 12:53
谢谢sobereva老师。
1 HOMO其余部分主要分布在片段CO2上，约为43.52%。但从图形上看，绝大部分HOMO分布 ...

没相关细节没法回复
相关文件压缩后上传，如果压缩后超过3MB，传到网盘里

作者
Author: dlut211 时间: 2019-5-23 11:20

sobereva 发表于 2019-5-22 14:00
没相关细节没法回复
相关文件压缩后上传，如果压缩后超过3MB，传到网盘里

附件中包含：HOMO、轨道组成、体系和各个片段的fchk文件。
谢谢老师。

作者
Author: sobereva 时间: 2019-5-24 05:03

dlut211 发表于 2019-5-23 11:20
附件中包含：HOMO、轨道组成、体系和各个片段的fchk文件。
谢谢老师。

fragment-3里的是log文件不是fchk文件，我没法考察
scf=qc这种关键词尽量不要用，不仅慢还容易收敛到不稳定波函数
每个片段的gjf也一起上传

作者
Author: dlut211 时间: 2019-5-24 10:06

sobereva 发表于 2019-5-24 05:03
fragment-3里的是log文件不是fchk文件，我没法考察
scf=qc这种关键词尽量不要用，不仅慢还容易收敛到不 ...

感谢老师的耐心回复，重新准备了文件，请您查阅。

作者
Author: sobereva 时间: 2019-5-25 03:58

dlut211 发表于 2019-5-24 10:06
感谢老师的耐心回复，重新准备了文件，请您查阅。

你的片段里的原子顺序和整体不一样，此时没法做CDA
你的复合物里CO2的原子序号是1、2、4，但是CO2片段的文件里它们是挨着的，所以依次载入每个片段后对应的原子顺序和整体不符，这势必会导致分析结果错乱

作者
Author: dlut211 时间: 2019-5-25 16:03

sobereva 发表于 2019-5-25 03:58
你的片段里的原子顺序和整体不一样，此时没法做CDA
你的复合物里CO2的原子序号是1、2、4，但是CO2片段的 ...

谢谢老师，按照您说的，已解决。

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