计算化学公社

标题: 如何计算一部分分子的偶极矩来进行向量叠加? [打印本页]

作者
Author:
落花飘雨679    时间: 2018-11-27 20:36
标题: 如何计算一部分分子的偶极矩来进行向量叠加?
老师,您好,在计算偶极矩的过程中,看到文章中由计算一个分子的两部分然后通过矢量叠加得到整个分子的偶极矩,相关文字与图示如下图所示,我知道整个分子的偶极矩如何计算,但是如何计算一部分分子的偶极矩来进行向量叠加?感谢老师的回复。

作者
Author:
adia    时间: 2018-11-27 21:20
可以提供一下文章信息吗,搜文字和图片居然没搜到,想了解一下相关的东西,谢谢!
作者
Author:
sobereva    时间: 2018-11-27 21:48
用Multiwfn计算ADCH电荷,产生.chg文件,此文件里将含有原子电荷信息。将.chg载入Multiwfn时,屏幕上就会提示基于原子电荷计算的偶极矩。手动将.chg中另一个片段部分删除,则载入Multiwfn显示的就是剩下的片段的偶极矩。因此可以容易地计算两个片段各自的偶极矩。
作者
Author:
落花飘雨679    时间: 2018-11-27 22:19
adia 发表于 2018-11-27 21:20
可以提供一下文章信息吗,搜文字和图片居然没搜到,想了解一下相关的东西,谢谢!

你好,文章为Adv. Mater. 2018, 30, 1800052
作者
Author:
落花飘雨679    时间: 2018-11-27 22:24
sobereva 发表于 2018-11-27 21:48
用Multiwfn计算ADCH电荷,产生.chg文件,此文件里将含有原子电荷信息。将.chg载入Multiwfn时,屏幕上就会提 ...

谢谢老师的回复,我按照您说的方法计算一下,还有一个问题需要麻烦老师帮忙,在优化分子的过程中,如何在考虑分子中的非共价相互作用的情况下得到最优构象?即优化分子的过程中如何将非共价相互作用考虑在内?
作者
Author:
sobereva    时间: 2018-11-27 22:37
落花飘雨679 发表于 2018-11-27 22:24
谢谢老师的回复,我按照您说的方法计算一下,还有一个问题需要麻烦老师帮忙,在优化分子的过程中,如何在 ...

你的计算级别若能合理描述弱相互作用,自然而然就考虑了
作者
Author:
落花飘雨679    时间: 2018-11-27 22:44
sobereva 发表于 2018-11-27 22:37
你的计算级别若能合理描述弱相互作用,自然而然就考虑了

老师,我懂了,如果这个体系中存在弱相互作用,我用考虑弱相互作用的泛函来计算就可以了。谢谢您的耐心解答。
作者
Author:
落花飘雨679    时间: 2018-11-29 19:27
sobereva 发表于 2018-11-27 21:48
用Multiwfn计算ADCH电荷,产生.chg文件,此文件里将含有原子电荷信息。将.chg载入Multiwfn时,屏幕上就会提 ...

老师您好,我用您说的方法将.chg文件导入Multiwfn中,显示总的偶极矩为 0.094789 Debye,而用高斯得到的结果为 0.1728Debye,请问老师我要以哪个结果为准?谢谢老师。
作者
Author:
sobereva    时间: 2018-11-29 23:32
落花飘雨679 发表于 2018-11-29 19:27
老师您好,我用您说的方法将.chg文件导入Multiwfn中,显示总的偶极矩为 0.094789 Debye,而用高斯得到的 ...

是bug,计算的时候单位转换疏忽了,我已经修复了并在官网上更新了程序。谢谢发现这个bug
作者
Author:
落花飘雨679    时间: 2018-11-30 08:14
sobereva 发表于 2018-11-29 23:32
是bug,计算的时候单位转换疏忽了,我已经修复了并在官网上更新了程序。谢谢发现这个bug

现在一致了,谢谢老师,老师我如果手动删掉了原子,计算其中一部分的偶极矩,这部分偶极矩的方向通过什么方法确定呢?
作者
Author:
sobereva    时间: 2018-11-30 11:04
落花飘雨679 发表于 2018-11-30 08:14
现在一致了,谢谢老师,老师我如果手动删掉了原子,计算其中一部分的偶极矩,这部分偶极矩的方向通过什么 ...

就是∑R(i)q(i)方式计算,R是原子坐标矢量,q是原子电荷,加和是对于chg里所有原子而言的




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3