用户工具

站点工具


adf:mosbauer

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

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

参数设置

Frozen Core务必选择None,Relativity选择Scalar,其他选项参考单点计算与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也有影响,后者更精确
adf/mosbauer.txt · 最后更改: 2019/12/18 17:29 由 liu.jun

© 2014-2020 费米科技(京ICP备14023855号