计算化学公社

标题: 请问通过量化计算可以获得分子的正负电荷中心位置吗 [打印本页]

作者
Author: fntg 时间: 2021-1-14 15:09
标题: 请问通过量化计算可以获得分子的正负电荷中心位置吗
或者可以获得正负电荷中心之间的距离吗
谢谢

作者
Author: wzkchem5 时间: 2021-1-14 17:11
这里取决于你怎么定义正负电荷
你既可以把原子核算作正电荷，把电子算作负电荷，也可以把带正电的原子算作正电荷，把带负电的原子算作负电荷，两种算法是不等价的。

作者
Author: fntg 时间: 2021-1-14 17:20

wzkchem5 发表于 2021-1-14 17:11
这里取决于你怎么定义正负电荷
你既可以把原子核算作正电荷，把电子算作负电荷，也可以把带正电的原子算作 ...

谢谢回复。
想知道和偶极矩计算相对应的正负电荷中心和它们之间的距离。

作者
Author: 北大-陶豫 时间: 2021-1-14 17:40

wzkchem5 发表于 2021-1-14 17:11
这里取决于你怎么定义正负电荷
你既可以把原子核算作正电荷，把电子算作负电荷，也可以把带正电的原子算作 ...

楼上已经说得很清楚了……
你既可以把原子核算作正电荷，把电子算作负电荷，这时候总电荷数大，正负电荷中心距离小；
你也可以把带正电的原子算作正电荷，把带负电的原子算作负电荷，这时候总电荷数小，正负电荷中心距离大；

无论如何计算，总的偶极矩都是一样的：
μ=∫rρ(r)dr（这里ρ(r)包括原子核和电子的贡献，要求∫ρ(r)dr=0）

作者
Author: wzkchem5 时间: 2021-1-14 18:32

fntg 发表于 2021-1-14 17:20
谢谢回复。
想知道和偶极矩计算相对应的正负电荷中心和它们之间的距离。

不影响我的问题，我说的那两种算法算出来的电荷和距离都是不相同的，但是偶极矩是相同的。
所以你必须结合你的实际需求，从我说的两种定义里面选一种，两种都有道理，所以我不能替你选。

作者
Author: 北大-陶豫 时间: 2021-1-14 23:53

北大-陶豫发表于 2021-1-14 17:40
楼上已经说得很清楚了……
你既可以把原子核算作正电荷，把电子算作负电荷，这时候总电荷数大，正负电荷 ...

我这么说也许你能明白：
偶极矩公式 μ=∫rρ(r)dr，我把某个原子核的一部分正电荷和它周围一部分球对称的电子密度单独拿出来，那么单独拿出来的部分可以做到 ∫rρ(r)dr=0，因此我们只要考虑剩余部分的 ∫rρ(r)dr 即可。所以，考虑原子核和所有电子形成的偶极矩，可以转化为考虑带一定电荷（有正有负）的原子们形成的偶极矩。

作者
Author: sobereva 时间: 2021-1-15 05:44
下载目前Multiwfn官网（http://sobereva.com/multiwfn）上最新版Multiwfn，载入fch、wfn等含有波函数信息的文件，进入主功能300的子功能5，屏幕上就会显示根据下式计算的正电荷中心和负电荷中心
(, 下载次数 Times of downloads: 58)
其中q_A指A的原子核电荷，X_A指A原子的X坐标，rho是电子密度函数，x是电子坐标x分量

作者
Author: fntg 时间: 2021-1-15 08:48

sobereva 发表于 2021-1-15 05:44
下载目前Multiwfn官网（http://sobereva.com/multiwfn）上最新版Multiwfn，载入fch、wfn等含有波函数信息的 ...

谢谢!

作者
Author: fntg 时间: 2021-1-15 08:48

北大-陶豫发表于 2021-1-14 23:53
我这么说也许你能明白：
偶极矩公式 μ=∫rρ(r)dr，我把某个原子核的一部分正电荷和它周围一部分球对称 ...

谢谢！

作者
Author: fntg 时间: 2021-1-15 08:49

wzkchem5 发表于 2021-1-14 18:32
不影响我的问题，我说的那两种算法算出来的电荷和距离都是不相同的，但是偶极矩是相同的。
所以你必须结 ...

谢谢！

作者
Author: fntg 时间: 2021-1-15 09:41

sobereva 发表于 2021-1-15 05:44
下载目前Multiwfn官网（http://sobereva.com/multiwfn）上最新版Multiwfn，载入fch、wfn等含有波函数信息的 ...

sob老师，我用Multiwfn试了试，输出了偶极矩和四极矩等，但是没有输出正负电荷中心。
如果只需要知道正负电荷中心的距离的话，是不是用计算出的偶极矩除以分子总的（核/电子）电量就可以了？
谢谢。

作者
Author: 北大-陶豫 时间: 2021-1-15 10:47

fntg 发表于 2021-1-15 09:41
sob老师，我用Multiwfn试了试，输出了偶极矩和四极矩等，但是没有输出正负电荷中心。
如果只需要知道正 ...

如果你把原子核算作正电荷、把电子算作负电荷，而不是把带正电的原子算作正电荷、把带负电的原子算作负电荷，那么就确实可以这么做。

作者
Author: fntg 时间: 2021-1-15 10:52

北大-陶豫发表于 2021-1-15 10:47
如果你把原子核算作正电荷、把电子算作负电荷，而不是把带正电的原子算作正电荷、把带负电的原子算作负电 ...

谢谢！

作者
Author: sobereva 时间: 2021-1-16 01:29

fntg 发表于 2021-1-15 09:41
sob老师，我用Multiwfn试了试，输出了偶极矩和四极矩等，但是没有输出正负电荷中心。
如果只需要知道正 ...

没输出绝对是因为版本不是最新的
我已经明确强调了必须用最新版本

(, 下载次数 Times of downloads: 46)

作者
Author: fntg 时间: 2021-1-16 09:55

sobereva 发表于 2021-1-16 01:29
没输出绝对是因为版本不是最新的
我已经明确强调了必须用最新版本

谢谢！

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