求助 VASP二氧化钛掺杂铜电子步无法收敛 恳请指教

2955778771

1.含钛空位的二氧化钛纳米片掺杂的铜，这个材料我之前做过很多关于掺杂镍的计算

2.现在想拓展一下体系，发现一模一样的模型和位点，掺杂Mn Fe Co Ni Zn都可以正常收敛
   如：RMM:  29    -0.919194498842E+03   -0.53133E-05   -0.68064E-05  3456   0.121E-02
          2 F= -.91919450E+03 E0= -.91918199E+03  d E =-.556719E+00  mag=     0.0000

3.但是换成Cu电子步无法收敛，INCAR KPOINT完全一致，POSCAR只是把Ni换成Cu
如:    RMM: 200    -0.777448739943E+03   -0.67072E+00   -0.73460E+01  7377   0.200E+01
在第5步离子步后，电子步在200时仍无法收敛  
非常不理解是为什么 明明只改了POSCAR中Ni变为Cu,而其他Mn Fe Co Ni Zn都可以正常收敛
恳请大佬赐教！！

SYSTEM = CuO(111) #用户可以自由输入注释
ISTART = 0    #随机生成波函数
ICHARG = 2   #通过原子电荷密度叠加生成电荷密度
PREC = Normal    #总体计算精度
LREAL = A   #确定投影算法在倒空间计算
IBRION = 2    #共轭梯度(GC)算法的离子弛豫
ISIF = 3    #进行原子弛豫同时优化晶格
NSW = 200  #离子运动的最大步数
POTIM = 0.2  #分子动力学运行的时间步长
EDIFFG = -0.02   #原子弛豫收敛的标准
ENCUT = 450 eV    #确定平面波截断能
NELM = 200    #允许电子迭代得到最大步数
EDIFF = 1E-05     #电子自洽收敛标准参数
LCHARG =T   #确定是否在CHGCAR和CHG中输出电荷密度
LWAVE = F    #确定是否在WAVECAR文件中输出波函数
LAECHG = T    #确定是否输出AGE
ISMEAR = 0      #Gaussian smearing方法（氧化物，半导体）
SIGMA = 0.05     #展宽
ALGO = Fast     #确定电子优化的算法 （变为N/V也不收敛）
NCORE= 4


Gamma
   2   1   1


PAW_PBE Cu 22Jun2005                  
   11.0000000000000     
parameters from PSCTR are:
   SHA256 =  0e71e558f37a4e8407882b4609ceaae59938afaad9ced71dfe8490308a0ab630 Cu/POTCAR
   COPYR  = (c) Copyright 22Jun2005 Georg Kresse
   COPYR  = This file is part of the software VASP. Any use, copying, and all other rights are regulated by the VASP license agreement.
   COPYR  = If you do not have a valid VASP license, you may not use, copy or distribute this file.
   VRHFIN =Cu: d10 p1
   LEXCH  = PE
   EATOM  =  1390.9808 eV,  102.2342 Ry
   TITEL  = PAW_PBE Cu 22Jun2005

LIJINWEI

2955778771 发表于 2025-11-25 18:45
已解决，分别使用DFT+U或者ALGO =Conjugate都可以收敛

您好，加U是给Cu加U吗?如果是给Cu加U，你是通过什么确定的？例如是看CuO的文献，看Cu设置的U多少，在您这个体系就设置多少？

2955778771

已解决，分别使用DFT+U或者ALGO =Conjugate都可以收敛