计算化学公社

标题: Multiwfn得到的偶极矩与ORCA和Gausian得到的差距较大 [打印本页]

作者
Author: qinsc 时间: 2024-1-27 19:49
标题: Multiwfn得到的偶极矩与ORCA和Gausian得到的差距较大
各位老师好，我使用orca在 wB97M-V/def2-TZVPD 级别下计算环己铵阳离子的偶极矩，参照的是sob老师的帖子 http://sobereva.com/558，计算得到的结果如图1，感觉不太对，这么明显的阳离子应该指向铵才对，随后用ORCA的%elprop 模块计算，得到了期望的结果如图2，高斯与orca同样的期望结果。
造成这样的结果是因为我在multiwfn中选错了什么吗？我使用主功能15，3 //基于内置原子密度的Hirshfeld权重函数，计算得到的结果。
请各位老师不吝赐教。

作者
Author: qinsc 时间: 2024-1-27 19:51
抱歉，不知道为什么图片顺序乱了，标注是对的，顺序错了。

作者
Author: ionexchangeC 时间: 2024-1-27 20:07
离子体系计算偶极矩没有意义

作者
Author: qinsc 时间: 2024-1-27 20:51

ionexchangeC 发表于 2024-1-27 20:07
离子体系计算偶极矩没有意义

确实，我去掉H离子两者结果就一致了，感谢指点。

作者
Author: sobereva 时间: 2024-1-28 01:00

qinsc 发表于 2024-1-27 19:51
抱歉，不知道为什么图片顺序乱了，标注是对的，顺序错了。

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

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

作者
Author: qinsc 时间: 2024-2-4 00:16

sobereva 发表于 2024-1-28 01:00
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

感谢sob老师提醒，我今后会注意。
静电荷不为0的体系偶极矩依赖原点是不是因为不同原点和离子的整体电荷间总有偶极，这一偶极依赖于原点位置，再加上分子的偶极才是整体的偶极，所以和原点位置相关呢？
对于这个离子，如果我原点和ORCA或者高斯默认的原点一致，结果是不是理论上就一样了呢？
还想问下sob老师，主功能15中计算偶极矩和300中5计算偶极矩的区别是不是只是分子偶极矩的两种计算方法，其计算结果本质是相同的，只是途径不同？这两种计算方式是否在不同方面有独特的优势呢？
谢谢sob老师。

作者
Author: sobereva 时间: 2024-2-4 04:10

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

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

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

3 主功能15是通过多中心数值积分得到的，主功能300的子功能5是解析积分得到的，前者耗时更高精度还更低

作者
Author: qinsc 时间: 2024-2-23 16:15

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

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

明白了，谢谢sob老师！

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