这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:chargepartition [2017/04/13 22:41] – liu.jun | adf:chargepartition [2020/11/16 19:56] (当前版本) – 创建 liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ====== | + | ======使用FDE方法强制体系的某个区域带电====== |
- | + | 该方法,在2020.101版中,存在bug,建议使用2019版,参考:[[adf:chargepartition2019]] | |
- | 对于某些特殊问题,我们希望人为地干预电荷分布,例如水分子解离为< | + | |
- | + | ||
- | 本例较为极端,但实际上对于两个配体互相靠近的情况,都存在本例中的问题,也就是电子分布是不合理的,尤其是HOMO、LUMO轨道的空间分布,往往是不正确的。 | + | |
- | + | ||
- | 这种情况下,我们可以使用ADF中的FDE方法来实现。此处我们用H-OH来作为示范的例子,其他大的分子体系,计算流程是一样的。 | + | |
- | + | ||
- | + | ||
- | 1,对一个< | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | **注意,Integration那一项,在新版本(2014版)中,增加了Spline Zlm fit选项。这个选项是新版本的默认方法,但对于本文的FDE功能,则应该取消,即:该选项选 None (use STO fit)。后面的几个ADFinput文件,也应该这样设置。** | + | |
- | + | ||
- | 这一步计算得到H+.t21文件,保留此文件,下一步计算会使用; | + | |
- | + | ||
- | 2,计算质子逐步远离分子,例如该质子离羟基3.96埃的轨道与能量: | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | 将质子和剩余部分分成两个区域(如图所示,表示为红色和绿色)——使用ADF的分区功能: | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | 将质子设置为FDE区域(点击下图中+,然后选择H_作为FDE区域,也就是将要限制H的带电量保持不变,勾选relax表示允许两个片段的电子进行弛豫,达到平衡): | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | <color grey> | + | |
- | + | ||
- | 同时使用上一步中得到的< | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | 保存,并运行。其中保存的时候,弹出H+片段的结果是否覆盖,选No,运行; | + | |
- | + | ||
- | 3,结果查看: | + | |
- | 在SCM LOGO > View窗口查看Mulliken电荷(注意此处< | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | 显示这样的计算结果,所有的电子确实分布在羟基上,质子上没有电子; | + | |
- | 分别显示所有的电子轨道: | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | 分别在orbitals(occupied)选中各个轨道,进而分别显示这些轨道,发现所有的电子轨道都分布在羟基上: | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | 显示总的电子密度空间分布,也可以看见所有电子都分布在< | + | |
- | + | ||
- | {{ : | + | |
- | 到此完成。 | + | |
- | + | ||
- | ADF软件提供免费试用(一般为一个月),试用申请方式参见**费米科技维基百科:[[adf:trial|]]** | + |