这里会显示出您选择的修订版和当前版本之间的差别。
后一修订版 | 前一修订版 | ||
adf:mosbauer [2016/09/27 09:39] – 创建 liu.jun | adf:mosbauer [2024/02/21 21:47] (当前版本) – [影响该计算精度的因素] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ====== 如何计算原子核处的电子密度(待续)====== | + | ====== 如何计算原子核处的电子密度、静电势、穆斯堡尔谱====== |
+ | |||
+ | 周期性体系的原子核处电子密度,可以用来计算Mössbauer谱。BAND使用STO基组与数值基组,精度、可信度高于GTO和平面波基组。 | ||
+ | |||
+ | =====参数设置===== | ||
+ | |||
+ | * Frozen Core务必选择None | ||
+ | * Relativity选择Scalar或Spin-Orbit | ||
+ | * Numerical Quality选择Good | ||
+ | * Details - Numerical Quality - K-Space选择Very Good或Excellent | ||
+ | |||
+ | 其他非必要选项参考[[adf: | ||
+ | |||
+ | * Vxc[rho(fit)]:在out文件中显示原子核处的交换相关势,但该势函数是由拟合出来的电子密度计算得到的。拟合电子密度是为了减小计算库伦积分的计算量,使用STO函数或者球谐函数来拟合出真实的电子密度,虽有误差,但几乎可忽略; | ||
+ | * rho(fit):在out文件中显示原子核处的电子密度,但该电子密度是拟合电子密度; | ||
+ | * rho(scf):在out文件中显示原子核处的电子密度,该电子密度是SCF计算得到的精确电子密度; | ||
+ | * V(coulomb/ | ||
+ | * rho(demormation/ | ||
+ | * rho(demormation/ | ||
+ | |||
+ | =====结果查看===== | ||
+ | 计算完毕之后,搜索:Properties at Nuclei,可以显示这些计算结果。例如下图是计算原子核电子密度的结果: | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | =====影响该计算精度的因素===== | ||
+ | - 该数值对基组很敏感,因此需要使用尽量大的基组; | ||
+ | - 原子核模型有点电荷模型和Gaussian型两种,对结果有影响,但一般不影响定性的大小关系,在Details > Relativity > Nuclear model里面可以设置; | ||
+ | - 是否采用相对论方法Scalar或者Spin-Orbit也有影响,后者更精确 | ||
+ | |||
+ | 对quadrupole splittings、Isomer shifts,参考ADF模块的处理方式:[[adf: |