求助Fe、Mn替换掺杂时掺杂的形成能计算公式中Fe、Mn单原子能量的计算

Cuuuuu

计算公式如图
在CASTEP中完成计算
一开始采取构建晶胞，在中心添加单原子，计算出单点能，代入到公式中得出掺杂形成能
现在采用计算出单胞能量再除原子个数，获得单个原子能量
问题：
1：构建单胞后采用几何优化后的Final Enthalpy还是单点能计算的Final energy来除原子数获得单原子能量？
Mn单原子能量，采用几何优化时为-2863.840183eV; 单点能计算为-2858.086168eV（在计算时均+U了）。在进一步计算Mn掺杂形成能时，如果采用单点能结果，掺杂能为-0.0901eV，采用几何优化后的能量，掺杂能为5.05eV
2：在计算Fe、Mn单原子能量时的计算参数需不需要+U
在计算Fe单原子能量时，如果用几何优化的能量，+U时为-860.99245，不+U为-863.762012；如果有单点能的能量，+U为-860.2020145369，不＋U为-862.8916154260

谢谢大家！不胜感激！

guotong

为啥你的几何优化后的能量和单点能的能量相差这么大？ 单点能计算是额外提高了计算精度？