计算化学公社

标题: SAPT与AIM IGMH分析结果对应问题 [打印本页]

作者
Author: 黑色桃花 时间: 2023-10-14 14:49
标题: SAPT与AIM IGMH分析结果对应问题
各位老师，大家好！
我选用了同一个阳离子，采用三个不同的阴离子，研究其性能。现在，在AIM分析里，1号离子体系氢键数量多（虽然都是弱氢键），而且在IGMH里可以看到1号离子体系范德华相互作用的面积是最大的。而2号离子体系，尽管也有氢键（仅有一个)，还有一个的电子密度较大0.06左右，它的范德华作用面积排第二。三号离子体系的氢键数量比2号的多，没1号的多，范德华作用面积最小。
但是，我在计算SAPT能量分解的时候，2号的静电相互作用最大，其次是1号，最后是3号。就是我一直不理解，为啥1号的氢键数量多，且范德华作用面积大，他的静电相互作用为啥不是最大的，而是2号最大，因为我觉得相互吸引作用中，无非就是氢键等还有就是范德华相互作用。
请各位老师赐教，谢谢！

作者
Author: sobereva 时间: 2023-10-14 17:40
文字描述过于抽象，把相关的图摆出来

作者
Author: 牧生 时间: 2023-10-14 18:42
不完全一定IGMH看到蓝色等值面，且VMD看到氢键，就认为是吸引作用，也不一定氢键数量多，吸引作用就强。
这个水二聚体的坐标，你可以拿去做IGMH分析和SAPT能量分解。IGMH图像上看到就是蓝色等值面，看起来是吸引；但是SAPT能量分解就是正值，应该就是排斥作用。我更相信SAPT能量分解的结果。

O    48.650  10.280  14.060
H    48.160  11.040  13.760
H    49.080  10.570  14.870
O    47.840  12.710  13.450
H    47.700  13.090  14.310
H    46.960  12.540  13.110

作者
Author: 黑色桃花 时间: 2023-10-14 19:02

牧生发表于 2023-10-14 18:42
不完全一定IGMH看到蓝色等值面，且VMD看到氢键，就认为是吸引作用，也不一定氢键数量多，吸引作用就强。
...

哇，您真的懂我的意思了，雷姆，那我想请教您一下，阴阳离子之间的极性是否代表着一把双刃剑，就是极性大了可能导致交换互斥能增加，但是也正是极性大，两者的静电相互作用也会增加，这个就得看谁的作用更大了，您说我说的对吗

作者
Author: sobereva 时间: 2023-10-14 19:55
IGMH等值面的颜色判断吸引还是互斥一般是对于平衡结构而言的，不是平衡结构颜色说明不了什么。诸如即便强行把原子间距离拉近乃至相互作用能为正了，sign(lambda2)仍然可能是负的，导致等值面颜色是蓝的。

作者
Author: 黑色桃花 时间: 2023-10-14 20:07
老师，我现在明白您的意思了，我现在把图上传了，静电相互作用能大的那个是2号的，我现在不理解的是为撒他静电大，同时交换互斥也大，我觉得是不是因为2号体系阴离子的极性太强导致的，二号阴离子是N3-，一号是N5-，请您批评指教

作者
Author: 黑色桃花 时间: 2023-10-14 20:08

sobereva 发表于 2023-10-14 19:55
IGMH等值面的颜色判断吸引还是互斥一般是对于平衡结构而言的，不是平衡结构颜色说明不了什么。诸如即便强行 ...

老师，具体图我上传了，请您批评指教，谢谢卢老师

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