计算化学公社

 找回密码 Forget password
 注册 Register

Multiwfn得到的偶极矩与ORCA和Gausian得到的差距较大

查看数: 1154 | 评论数: 7 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-1-27 19:49

正文摘要:

各位老师好,我使用orca在 wB97M-V/def2-TZVPD 级别下计算环己铵阳离子的偶极矩,参照的是sob老师的帖子 http://sobereva.com/558,计算得到的结果如图1,感觉不太对,这么明显的阳离子应该指向铵才对,随后用ORCA的 ...

回复 Reply

qinsc 发表于 Post on 2024-2-23 16:15:09
sobereva 发表于 2024-2-4 04:10
1 根据偶极矩表达式,自己推一下就知道带净电荷的体系的偶极矩对原点存在依赖性

2 是。注意高斯的标准 ...

明白了,谢谢sob老师!
sobereva 发表于 Post on 2024-2-4 04:10:00
qinsc 发表于 2024-2-4 00:16
感谢sob老师提醒,我今后会注意。
静电荷不为0的体系偶极矩依赖原点是不是因为不同原点和离子的整体电荷 ...

1 根据偶极矩表达式,自己推一下就知道带净电荷的体系的偶极矩对原点存在依赖性

2 是。注意高斯的标准朝向问题
谈谈Gaussian中的对称性与nosymm关键词的使用
http://sobereva.com/297http://bbs.keinsci.com/thread-1355-1-1.html

3 主功能15是通过多中心数值积分得到的,主功能300的子功能5是解析积分得到的,前者耗时更高精度还更低
qinsc 发表于 Post on 2024-2-4 00:16:54
sobereva 发表于 2024-1-28 01:00
有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员 ...

感谢sob老师提醒,我今后会注意。
静电荷不为0的体系偶极矩依赖原点是不是因为不同原点和离子的整体电荷间总有偶极,这一偶极依赖于原点位置,再加上分子的偶极才是整体的偶极,所以和原点位置相关呢?
对于这个离子,如果我原点和ORCA或者高斯默认的原点一致,结果是不是理论上就一样了呢?
还想问下sob老师,主功能15中计算偶极矩和300中5计算偶极矩的区别是不是只是分子偶极矩的两种计算方法,其计算结果本质是相同的,只是途径不同?这两种计算方式是否在不同方面有独特的优势呢?
谢谢sob老师。
sobereva 发表于 Post on 2024-1-28 01:00:40
qinsc 发表于 2024-1-27 19:51
抱歉,不知道为什么图片顺序乱了,标注是对的,顺序错了。

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。

净电荷不是0的体系,偶极矩是依赖于原点的选取的,所以计算离子的偶极矩没什么实际意义,而且就算要在不同程序间横向对比,也必须保证坐标是一致的。
而且计算偶极矩跟计算原子电荷没有直接联系,Multiwfn里你若要基于当前波函数文件里的信息计算偶极矩和多极矩,最好的做法是用主功能300的子功能5,会通过解析的积分计算,又精确又快。
qinsc 发表于 Post on 2024-1-27 20:51:33
ionexchangeC 发表于 2024-1-27 20:07
离子体系计算偶极矩没有意义

确实,我去掉H离子两者结果就一致了,感谢指点。
ionexchangeC 发表于 Post on 2024-1-27 20:07:12
离子体系计算偶极矩没有意义
qinsc 发表于 Post on 2024-1-27 19:51:27
抱歉,不知道为什么图片顺序乱了,标注是对的,顺序错了。

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

GMT+8, 2025-8-16 15:26 , Processed in 0.310991 second(s), 25 queries , Gzip On.

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