如何计算原子核处的电子密度、静电势、穆斯堡尔谱

周期性体系的原子核处电子密度，可以用来计算Mössbauer谱。BAND使用STO基组与数值基组，精度、可信度高于GTO和平面波基组。

• Frozen Core务必选择None
• Relativity选择Scalar或Spin-Orbit
• Numerical Quality选择Good
• Details - Numerical Quality - K-Space选择Very Good或Excellent

其他非必要选项参考【入门基础教程】单点计算与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，可以显示这些计算结果。例如下图是计算原子核电子密度的结果：

1. 该数值对基组很敏感，因此需要使用尽量大的基组；
2. 原子核模型有点电荷模型和Gaussian型两种，对结果有影响，但一般不影响定性的大小关系，在Details > Relativity > Nuclear model里面可以设置；
3. 是否采用相对论方法Scalar或者Spin-Orbit也有影响，后者更精确

对quadrupole splittings、Isomer shifts，参考ADF模块的处理方式：原子核处电子密度、Isomer shifts、Mössbauer 四极分裂与 Mössbauer 谱