计算化学公社

标题: ORCA输出的偶极矩与Multiwfn输出的差异问题 [打印本页]

作者
Author: sdjgdu 时间: 2020-8-16 11:57
标题: ORCA输出的偶极矩与Multiwfn输出的差异问题
对同一个体系，ORCA输出的偶极矩只有 6.1 Debye，而使用multiwfn （15->2）得到的偶极矩竟然有 57.0 Debye，整整差一个数量级。为什么有这么大的差距，哪个准确呢？

作者
Author: sobereva 时间: 2020-8-16 12:03
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的严重不恰当标题“偶极矩”改了，以后务必注意，下次将扣分处理。

作者
Author: sobereva 时间: 2020-8-16 12:04
说清楚什么版本的Multiwfn（精确到更新日期），上传压缩后的给Multiwfn用的文件

作者
Author: sdjgdu 时间: 2020-8-16 12:51
谢谢提醒! 使用 Version 3.7(dev), release date: 2020-Aug-2，nultiwfn分析利用.molden文件，文件压缩后仍然有38mb无法上传

作者
Author: sdjgdu 时间: 2020-8-16 13:03
文件以上传到百度网盘，链接: https://pan.baidu.com/s/1kssU32FKSJYCaBn7dcvM5A 提取码: ci6f

作者
Author: sobereva 时间: 2020-8-16 21:25
你这个体系有净电荷，讨论偶极矩没实际意义
把ORCA输出文件也压缩后上传，我感觉你看到的ORCA的输出数据和这个molden文件肯定不对应，毕竟体系又有净电荷，体系中心又偏离坐标原点非常显著，偶极矩注定很大。注意看看到底输出文件的坐标和molden里的坐标是否一致。

作者
Author: sdjgdu 时间: 2020-8-17 11:26
比较输出与molden，结构没问题，只不过单位不一样而已（一个是埃，一个是bohr）。输出文件见附件，是否带电荷的系统不适合用偶极矩表述？另有没有定性讨论不同复合物其水中溶解能力的描述量（极性表面比率可否？）

作者
Author: snljty 时间: 2020-8-17 11:58

sdjgdu 发表于 2020-8-17 11:26
比较输出与molden，结构没问题，只不过单位不一样而已（一个是埃，一个是bohr）。输出文件见附件，是否带电 ...

电场展开后有静电荷项、偶极矩项、四极矩项、八极矩项，...可以证明，仅当前n项全为0的时候，第n+1项如果是非0项，才和原点的选取无关。进而可知此时第n+2项和原点的选取有直接关系。对于有静电荷的体系，选不同的坐标原点得到的偶极矩不一样。

作者
Author: sobereva 时间: 2020-8-17 12:08
你最好拿个单点任务进行对比
当前输出信息也太臃肿了，还牵扯到有限差分之类会改变结构的事情

作者
Author: sdjgdu 时间: 2020-8-17 12:17

snljty 发表于 2020-8-17 11:58
电场展开后有静电荷项、偶极矩项、四极矩项、八极矩项，...可以证明，仅当前n项全为0的时候，第n+1项如果 ...

谢谢！

作者
Author: sdjgdu 时间: 2020-8-17 12:18

sobereva 发表于 2020-8-17 12:08
你最好拿个单点任务进行对比
当前输出信息也太臃肿了，还牵扯到有限差分之类会改变结构的事情

谢谢！我算下试试

作者
Author: hubo 时间: 2023-8-10 16:09

sdjgdu 发表于 2020-8-17 12:18
谢谢！我算下试试

楼主你好，请问一下你后来弄明白原因了吗

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