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