计算化学公社

标题: 求助:NiO(311)优化弛豫时结构变化很大。 [打印本页]

作者
Author:
loovfnd    时间: 2021-6-24 09:19
标题: 求助:NiO(311)优化弛豫时结构变化很大。
NiO优化弛豫时结构变化很大(底层固定)。请问初始结构和输入文件有什么问题呢。
初始构型如下:
(, 下载次数 Times of downloads: 20)
首次优化时INCAR文件如下:
SYSTEM =  qvasp

INCAR-basic
   PREC   = Normal
   ENCUT  = 400 eV
   NELMIN= 5
   LREAL = Auto

#   ALGO = Fast
   EDIFF  = 1E-4
   ISMEAR = 0
   SIGMA = 0.05
   ISYM = 0

#   LDAU=.TRUE.  ## For LDA + U set
#   LDAUTYPE=2
#   LDAUL=2 -1
#   LDAUU=7     0
#   LDAUJ=1     0

   ISPIN = 2    ## For magnetism set
#  MAGMOM   = 2.0 -2.0 2*0
   NPAR = 4     ## Important to reduce calculate time, equal sqrt(cpu cores)

Ionic Relaxation
   ICHARG = 1
   EDIFFG = -0.05        force
   NSW    = 300          max steps for IOM
   IBRION = 2
   POTIM = 0.5
   ISIF   = 2            suface:2 crystal:3
   PSTRESS= 0
   LCHARG = .FALSE.
   LWAVE  =  F

# Use 'qvasp -cif' to transer format (VESTA or MS)

(, 下载次数 Times of downloads: 12)


再次优化时INCAR文件加上了 LDIPOL = .TRUE. IDIPOL = 3
构型变化还是很大
(, 下载次数 Times of downloads: 11)

作者
Author:
dmh1998dmh    时间: 2021-7-3 16:21
结构有保持化学计量比吗?而且上下表面不对称呐
另外有优化过晶格参数吗?
作者
Author:
loovfnd    时间: 2021-7-5 10:37
dmh1998dmh 发表于 2021-7-3 16:21
结构有保持化学计量比吗?而且上下表面不对称呐
另外有优化过晶格参数吗?

保持了化学计量比了。优化过晶格参数了。上下不对称要怎么处理?
作者
Author:
dmh1998dmh    时间: 2021-7-7 15:35
loovfnd 发表于 2021-7-5 10:37
保持了化学计量比了。优化过晶格参数了。上下不对称要怎么处理?

重新切slab,使得上下表面对称,做成非极性的表面
作者
Author:
jiangning198511    时间: 2021-7-7 16:27
dmh1998dmh 发表于 2021-7-7 15:35
重新切slab,使得上下表面对称,做成非极性的表面

估计上下表面对称就不一定保持化学计量比了
作者
Author:
cokoy    时间: 2021-7-7 18:07
设置一下磁矩初猜




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3