计算化学公社

标题: 求助 VASP二氧化钛掺杂铜电子步无法收敛恳请指教 [打印本页]

作者
Author: 2955778771 时间: 2025-11-24 21:15
标题: 求助 VASP二氧化钛掺杂铜电子步无法收敛恳请指教
本帖最后由 2955778771 于 2025-11-25 11:08 编辑

(, 下载次数 Times of downloads: 2)
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

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

作者
Author: LIJINWEI 时间: 2026-1-5 10:30

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)