这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 | ||
adf:dipolemoment [2016/07/07 14:23] – [组分2] liu.jun | adf:dipolemoment [2020/12/01 14:19] – liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======如何计算分子的电偶极矩====== | + | ======分子的电偶极矩的计算====== |
- | ADF模块任何一个计算,都会自动给出分子的电偶极矩,以及分子中每个原子的电偶极矩、单极电荷。整个分子的电偶极矩可以由分子中每个原子的电偶极矩、单极电荷、原子的坐标得到。如果是多个分子形成的弱相互作用体系,也可以通过这种方式,把其中每个分子的偶极矩计算出来。以下以一个H< | + | ADF模块任何一个计算,都会自动给出分子的电偶极矩,以及分子中每个原子的电偶极矩、单极电荷。整个分子的电偶极矩可以由分子中每个原子的电偶极矩、单极电荷、原子的坐标得到。如果是多个分子形成的弱相互作用体系,也可以通过这种方式,把其中每个分子的偶极矩计算出来。以下以一个H< |
- | =====找到所需数据===== | + | |
- | **找到原子的坐标**(ADF实际上把ADFinput中输入的分子进行了转动、平移,得到偏向对称化的坐标): | + | =====结果查看===== |
+ | |||
+ | **这个值在*.out文件中是可以直接找到的,程序已经自动帮我们计算好了。** | ||
+ | |||
+ | {{: | ||
+ | |||
+ | 一般情况,我们不需要自己手动计算。而且这个数字是来自拟合电子密度,因此有一点点误差,上图中第一个方框中的数值即是我们上面计算的值。第二个方框中是精确电子密度计算出来的值,二者有一点误差,但基本上可以忽略。 | ||
+ | |||
+ | 什么时候需要手动计算呢?也就是计算多个分子之间的诱导偶极、诱导力的时候,体系是多个分子,或者多个片段,需要手动的计算出各个分子片段的偶极矩,以此计算诱导偶极、诱导力。 | ||
+ | |||
+ | **总结:这里计算出了电偶极矩矢量、模,以及电偶极子的中心坐标。** | ||
+ | =====数据背后的计算方法===== | ||
+ | ====一、找到所需数据==== | ||
+ | **找到原子的坐标**(ADF实际上把AMSinput中输入的分子进行了转动、平移,得到偏向对称化的坐标): | ||
FRAGMENTS | FRAGMENTS | ||
行 23: | 行 35: | ||
3 H 0.090160 | 3 H 0.090160 | ||
- | =====计算分子偶极矩的两个组分===== | + | ====二、计算分子偶极矩的两个组分==== |
| | ||
- | ====组分1==== | + | 计算过程可以考虑使用Excel表格的计算功能。 |
+ | |||
+ | ===1,组分1=== | ||
**计算所有原子的偶极矩叠加**,原子dip总和: | **计算所有原子的偶极矩叠加**,原子dip总和: | ||
行 36: | 行 50: | ||
这是原子偶极总和对整个分子偶极矩的贡献。 | 这是原子偶极总和对整个分子偶极矩的贡献。 | ||
- | + | ===2,组分2=== | |
- | ====组分2==== | + | |
**计算原子的单极电荷对整个分子偶极矩的贡献**,也就是原子的**单极电荷**乘以各自坐标: | **计算原子的单极电荷对整个分子偶极矩的贡献**,也就是原子的**单极电荷**乘以各自坐标: | ||
行 63: | 行 76: | ||
正电中心坐标和负电中心坐标,二者的均值为**电偶极子的中心**。 | 正电中心坐标和负电中心坐标,二者的均值为**电偶极子的中心**。 | ||
- | =====计算分子的偶极矩===== | + | ====三、计算分子的偶极矩==== |
分子的偶极矩=dip总和+单极电荷对整个分子偶极矩的贡献: | 分子的偶极矩=dip总和+单极电荷对整个分子偶极矩的贡献: | ||
行 72: | 行 85: | ||
单位换算为Debye:0.6677344 a.u.=0.6677344*2.542=1.6973808 Debye | 单位换算为Debye:0.6677344 a.u.=0.6677344*2.542=1.6973808 Debye | ||
- | =====说明===== | ||
- | **这个值在*.out文件中是可以直接照到的,程序已经自动帮我们计算好了。** | ||
- | |||
- | {{: | ||
- | |||
- | 一般情况,我们不需要自己手动计算。而且这个数字是来自拟合电子密度,因此有一点点误差,上图中第一个方框中的数值即是我们上面计算的值。第二个方框中是精确电子密度计算出来的值,二者有一点误差,但基本上可以忽略。 | ||
- | |||
- | 什么时候需要手动计算呢?也就是计算多个分子之间的诱导偶极、诱导力的时候,体系是多个分子,或者多个片段,需要手动的计算出各个分子片段的偶极矩,以此计算诱导偶极、诱导力。 | ||
- | |||
- | **总结:这里计算出了电偶极矩矢量、模,以及电偶极子的中心坐标。** |