|
本帖最后由 Mars-allen 于 2024-12-25 19:49 编辑
由于材料中O在高电压会有氧化还原反应,Mn的d轨道会影响O的荷电状态,为了相对准确的计算结构和能量考虑使用GGA+U的来计算,但是目前计算之后发现两个问题
1计算报错结构优化过程中出现了:ZBRENT: fatal error in bracketing please rerun with smaller EDIFF, or copy CONTCAR to POSCAR and continue 的错误,查阅后发现说是能量和力震荡造成的有三种说法1是修改IBRION=1 而后添加参数;2是POTIM = 0.5数值太大造成震荡;3说我结构优化步骤有问题应该把晶胞优化和原子优化分开先让ISIF=2再ISIF= 3那种相对正确?
2计算速度及其慢我第二个模型144个原子结构优化了一天多还没出来结果我的参数否有问题造成计算缓慢?附INCAR如下:(计算硬件参数1个节点128核)
INCAR
#Starting parameters
SYSTEM = Li2MnO3-P-Lirich
ISTART = 0
ICHARG = 2
LWAVE = .F.
LCHARG = .F.
NWRITE = 1
#MAGMOM=1*1.22*0.2
LMAXPAW = 4
LMAXMIX = 4
ADDGRID = .TRUE.
LASPH = .TRUE.
ISMEAR = 0
SIGMA = 0.1
#LORBIT=11
#Exchange correlation functional
GGA=PE#PBE;91:PW91;AM05;PS:PBEsO1;
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = -1 2 -1
LDAUU = 0.0 4.9 0.0
LDAUJ = 0.0 0.0 0.0
#Parameters for Ionic relaxations
IBRION = 2
POTIM = 0.5
NSW = 300
ISIF = 3
ISYM = 0
EDIFFG = -0.01
#other parameters
PREC = A
NCORE = 8
|
|