请问，Ni/CeO2结构优化，能量收敛了，但是力一直不收敛怎么调整？

LeePanda

Incar如下：
SYSTEM = PROJECT_NAME






# ========== 基础计算模式设置 ==========
ISTART = 0        # 续算模式，读取CHGCAR/WAVECAR加速收敛
ICHARG = 2        # 读取电荷密度，续算核心参数
ISPIN  = 2        # 自旋极化（磁性计算必开）
LORBIT = 11       # 输出磁矩/态密度，不影响收敛


# ========== 基础计算 ==========
PREC   = Normal   # 精度适中，兼顾速度与稳定性
ENCUT  = 450      # 保持不变，适配赝势即可
LREAL  = Auto     # 团簇/表面体系适配，保持不变
ALGO   = Fast     # 替换All：Fast算法对磁性+U体系收敛更快
NELM   = 60      # SCF最大迭代步数60




# ========== 磁性初始磁矩 MAGMOM ==========
# 顺序：O 53个 → Ce 27个 → Ni 13个
MAGMOM = 53*0.0   27*0.5   13*0.5  # 保持不变，初始磁矩合理


# ========== DFT+U 设置 ==========
LDAU=.True.
LDAUTYPE  = 2     # 保持不变，常用的DFT+U类型
LDAUL     = -1  3  2    # O(p=0), Ce(f=3), Ni(d=2)
LDAUU     = 0.0 4.5 3.5 # U值合理，保持不变
LDAUJ     = 0.0 0.0 0.0
LMAXMIX   = 6        # 替换4：Ce(f轨道l=3)，设为6更稳定电子结构
LASPH     = .TRUE.   # 保持不变，球形谐波修正必开


# ========== 展宽（氧化物+金属团簇标准） ==========
ISMEAR   = 2      # 高斯展宽，保持不变
SIGMA    = 0.05    # 展宽值合理，保持不变




# ========== 结构优化 ==========
IBRION   = 1        # 替换2：准牛顿法（RMM-DIIS），力收敛比共轭梯度更高效
POTIM    = 0.2      # 替换0.3：降低步长，减少Ce原子震荡
NSW      = 60      # 核心调整：新任务跑100步优化
ISIF     = 2        # 只优化原子，保持不变


# ========== 收敛阈值（适配续算） ==========
EDIFF    = 1E-5     # 保持不变，能量收敛阈值合理
EDIFFG   = -0.04    # 替换-0.04：适度调松力收敛阈值，先让体系稳定


# ========== 电子混合参数（核心修改：防震荡）==========
AMIX     = 0.02      # 替换0.05：降低电荷混合系数，减少震荡
BMIX     = 0.001     # 保持不变
AMIN     = 0.001     # 保持不变
BMIX_MIX = 0.001     # 保持不变
MAXMIX   = 100       # 替换80：增加最大混合步数，给足收敛时间
IMIX     = 4         # Anderson混合，鲁棒性最强（适配磁性+U体系）


# ========== 泛函与弱相互作用修正 ==========
GGA    = PE       # PBE泛函，保持不变
IVDW   = 11       # DFT-D3色散修正，保持不变


# ========== 输出 ==========
LCHARG = .TRUE.   # 保存电荷密度，续算用
LWAVE  = .TRUE.   # 保存波函数，续算核心（当前80步的波函数会被读取）
LVTOT  = .FALSE.
NWRITE = 2        # 输出精简，保持不变


# ========== 并行效率设置 ==========
KPAR     = 1
NCORE  = 8        # 96核最优参数，保持不变
LMAXFOCK = 6      # 替换4：适配Ce的f轨道，稳定Fock矩阵计算





尝试了好久，就是不收敛。如图：



请问，下一步应该怎么调整呀？