如何计算原子核处的电子密度、静电势、穆斯堡尔谱
周期性体系的原子核处电子密度,可以用来计算Mössbauer谱。BAND使用STO基组与数值基组,精度、可信度高于GTO和平面波基组。
参数设置
其他非必要选项参考【入门基础教程】单点计算与BAND的基本参数设置、计算效率即可。在Properties > At nuclei勾选:
Vxc[rho(fit)]:在out文件中显示原子核处的交换相关势,但该势函数是由拟合出来的电子密度计算得到的。拟合电子密度是为了减小计算库伦积分的计算量,使用STO函数或者球谐函数来拟合出真实的电子密度,虽有误差,但几乎可忽略;
rho(fit):在out文件中显示原子核处的电子密度,但该电子密度是拟合电子密度;
rho(scf):在out文件中显示原子核处的电子密度,该电子密度是SCF计算得到的精确电子密度;
V(coulomb/scf):在out文件中显示原子核处的静电势,该静电势由SCF自洽迭代的电子密度精确地得到;
rho(demormation/fit):在out文件中显示原子核处的形变电子密度,也就是由片段形成晶体的过程中,原子核处的电子密度变化量。该密度是拟合出来的电子密度;
rho(demormation/scf):在out文件中显示原子核处的形变电子密度,也就是由片段形成晶体的过程中,原子核处的电子密度变化量。该密度是SCF计算得到的精确电子密度;
结果查看
计算完毕之后,搜索:Properties at Nuclei,可以显示这些计算结果。例如下图是计算原子核电子密度的结果:
影响该计算精度的因素
该数值对基组很敏感,因此需要使用尽量大的基组;
原子核模型有点电荷模型和Gaussian型两种,对结果有影响,但一般不影响定性的大小关系,在Details > Relativity > Nuclear model里面可以设置;
是否采用相对论方法Scalar或者Spin-Orbit也有影响,后者更精确
对quadrupole splittings、Isomer shifts,参考ADF模块的处理方式:原子核处电子密度、Isomer shifts、Mössbauer 四极分裂与 Mössbauer 谱