计算化学公社

标题: 多中心键级这一描述符应用的一些疑问 [打印本页]

作者
Author: radical 时间: 2022-4-19 11:47
标题: 多中心键级这一描述符应用的一些疑问
本帖最后由 radical 于 2022-4-19 14:36 编辑

(, 下载次数 Times of downloads: 5) 1.对于下图H4-四氮杂卟啉（加氢，并不是质子化,电荷与多重度：0 1，没有虚频），AICD的感应电流（磁场从屏幕里朝向屏幕外）和NICS(1)ZZ(数值为26.07）都显示出了明显的反芳香性，但通过考察它的MCI-pi(因为是非平面体系，已用multiwfn对占据轨道定域化），发现内环0.343，外环0.559，我之后，又单纯只算了MCI,，内环0.447，外环0.559，无论是哪个环，都应是是20个pi电子(4n),MCI的数值不应该这么大，对于类似的体系H0-四氮杂卟啉（去2个氢），它的AICD和NICS都显示它是非芳香性的，从MCI-pi(内环0.058，外环0.083)上看也是吻合的。
H4-四氮杂卟啉

(, 下载次数 Times of downloads: 13) (, 下载次数 Times of downloads: 14)
H0-四氮杂卟啉

(, 下载次数 Times of downloads: 16) (, 下载次数 Times of downloads: 12)

2.对四氮杂卟啉质子化（H3P1,电荷多重度： +1 1，没有虚频），同样是非平面体系，对轨道定域化处理后，MCI-PI值为负值（-0.428），MCi为0.538，请问负值是什么原因，即便是反芳香性不也应该是接近0吗
(, 下载次数 Times of downloads: 19)
fchk文件不知道为啥上传不了。。。显示过大

作者
Author: snljty 时间: 2022-4-19 13:50
fch压缩后上传，看首页的提醒 (, 下载次数 Times of downloads: 10)

作者
Author: radical 时间: 2022-4-19 14:27

snljty 发表于 2022-4-19 13:50
fch压缩后上传，看首页的提醒

好的，fch文件已上传

作者
Author: sobereva 时间: 2022-4-20 02:53
对这些体系没必要专门计算pi电子，本身sigma也不会对离域有贡献
对于原子数非常多的环，归一化的多中心键级的数值区分度不大。这种时候适合用原始的多中心键级对比芳香性的相对高低，也可以用一个芳香性已知较高的环的值作为参照来判断其它的环。

例如对于卟啉，在B3LYP/6-31G*级别下，我图中所示的路径的多中心键级是0.0000568579
(, 下载次数 Times of downloads: 10)

而你的体系，按照下图路径计算多中心键级（序号输入15,13,22,3,1,2,4,24,8,10,9,23,19,17,16,18,21,14），数值为0.0000076732，相比之下可知芳香性极低，要么非芳香性要么反芳香性
(, 下载次数 Times of downloads: 10)

作者
Author: radical 时间: 2022-4-20 09:32
本帖最后由 radical 于 2022-4-20 10:29 编辑

sobereva 发表于 2022-4-20 02:53
对这些体系没必要专门计算pi电子，本身sigma也不会对离域有贡献
对于原子数非常多的环，归一化的多中心键 ...

谢谢老师，那对于原子数不同的共轭路径，这个时候用原始的MCI数据去比较不同环的芳香性还有意义吗，我用AV1245去比较的不同环的芳香性得到的差值还是比较明显的，但问题是但对于个别物质，比如H4-四氮杂卟啉的芳香性描述与AICD和NICS并不是很符合

作者
Author: sobereva 时间: 2022-4-21 06:34

radical 发表于 2022-4-20 09:32
谢谢老师，那对于原子数不同的共轭路径，这个时候用原始的MCI数据去比较不同环的芳香性还有意义吗，我用A ...

原始的MCI不可能对不同原子数的环之间进行对比
对于原子数很多的路径，AV1245比归一化的MCI在数值上对芳香性强弱的区分度往往更高一些

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