计算化学公社

 找回密码 Forget password
 注册 Register

Gaussian中如下计算的到的esp电荷变化是正确的吗以及为什么有一个很奇怪的点?

查看数: 1878 | 评论数: 7 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-10-31 10:09

正文摘要:

利用mk方法计算一价一氧化碳离子随键长变化时的原子电荷变化,输入如下 #p scan b3lyp/def2TZVP pop=(always,allorbitals,mk) geom=connectivity iop(6/42=6) scf=xqc Title Card Required 1 2 O  &n ...

回复 Reply

sobereva 发表于 Post on 2022-11-3 08:52:29
柯柯爱喝热可可 发表于 2022-11-2 15:21
谢谢社长!我对那个突变点进行了稳定性测试结果是显示波函数具有内部不稳定性的,进行了一些尝试但没解决 ...

仔细看《原子电荷计算方法的对比》(http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml

算Mulliken电荷用的基组不是越大越好。而且即便是都是比如6-31G*,不同原子电荷的差异也很显著,因为计算原理截然不同,上文说了
柯柯爱喝热可可 发表于 Post on 2022-11-2 16:03:52
wzkchem5 发表于 2022-11-2 15:49
本来Mulliken电荷就是不能适用于大基组的,基组越大越没用。single zeta基组的结果比较靠谱(但因为电子 ...

好的好的我明白了,谢谢您的解释!我是修物理的这学期刚跨选了量子化学的课程,了解的还比较少,之后会好好学习量化基础知识的,再次感谢
wzkchem5 发表于 Post on 2022-11-2 15:49:55
柯柯爱喝热可可 发表于 2022-11-2 08:21
谢谢社长!我对那个突变点进行了稳定性测试结果是显示波函数具有内部不稳定性的,进行了一些尝试但没解决 ...

本来Mulliken电荷就是不能适用于大基组的,基组越大越没用。single zeta基组的结果比较靠谱(但因为电子结构本身不靠谱,所以总的来说还是不可靠),double zeta基组的结果就只能看个定性趋势了,triple zeta基组可能连定性趋势都不可靠了。
一般量化教材都会讲为什么大基组Mulliken电荷结果完全不可信
柯柯爱喝热可可 发表于 Post on 2022-11-2 15:21:02
sobereva 发表于 2022-11-1 04:45
有可能超过某个点的时候电子结构就是有显著变化,是事实。也可能中途出现了波函数不稳定,可以对一些点单 ...

谢谢社长!我对那个突变点进行了稳定性测试结果是显示波函数具有内部不稳定性的,进行了一些尝试但没解决这个问题。
除此之外我还是比较好奇的是为什么def2tzvp基组计算得到的Mulliken电荷和esp的电荷差别这么大,Mulliken电荷不是可以作为原子电荷计算的粗略参考嘛,如果这两种走势都不同,是不是意味着这样的计算结果是错误的呢?
sobereva 发表于 Post on 2022-11-1 04:45:19
柯柯爱喝热可可 发表于 2022-10-31 11:25
谢谢sob老师!我看了Mulliken电荷,相同的地方也有一个很小的突变,而且和esp电荷差别很大,走势都是相反 ...

有可能超过某个点的时候电子结构就是有显著变化,是事实。也可能中途出现了波函数不稳定,可以对一些点单独算单点并检查波函数稳定性,和之前扫描的能量、原子电荷进行对比
sobereva 发表于 Post on 2022-10-31 11:03:41
看看Mulliken电荷的变化情况,如果是平滑的,MK电荷的突变可能是拟合点变化带来的巧合。可以再尝试CHELPG

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 11:21 , Processed in 0.170702 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list