LREAL=.FALSE. LPLANE=.TRUE. NBANDS=500 NPAR=4 ISYM=-1这些参数主要是为了进行静态计算,并得到WAVECAR文件。接下来,我们需要定义能量窗口和基组函数:bash COHPstartEnergy=-30 COHPendEnergy=10 basisfunctions Fe 3d 4s basisfunctions C 2p 2s basisfunctions N 2p 2s basisfunctions O 2p 2s这里,我们定义了...
ENCUT = 520 LREAL = .FALSE. IBRION = 2 EDIFFG =-0.01 NELM = 100 ISIF = 3 NSW = 100 ISMEAR = 1 SIGMA = 0.05 ISPIN = 2 MAGMOM = 4*0.6 LMAXMIX= 4 LASPH = True INCAR.SCF(ELF) PREC = Accurate ALGO = Fast EDIFF = 1E-6 ENCUT = 520 LREAL = .FALSE. IBRION = -1 EDIFFG =-...
SYSTEM = InSeISTART = 1NWRITE = 2PREC = AccurateENCUT = 500GGA = PENSW = 0ISIF = 2ISYM = 2IBRION = -1NELM = 80EDIFF = 1E-05EDIFFG = -0.01ALGO = NormalLDIAG = .TRUE.LREAL = .FALSE.ISMEAR = 0SIGMA = 0.05ICHARG = 2LORBIT...
Global ParametersISTART=0(Read existing wavefunction;ifthere)ICHARG=2(Non-self-consistent: GGA/LDA band structures)LREAL=F(Projection operators: automatic)ENCUT=500(Cut-off energyforplane wave basis set, in eV)PREC=Accurate(Precision level)LWAVE=.FALSE.(Write WAVECAR or not)LCHARG=.FALSE.(Write...
NELMIN=8一般用6到8,最小的电子scf数.太少的话,收敛的不好.LREAL=A APACO=10径向分布函数距离,单位是埃.NPACO=200径向分布函数插的点数. LCHARG=F尽量不写电荷密度,否则CHG文件太大.TEBEG=300初始温度. TEEND=300终态温度。不设的话,等于TEBEG. ...
3、事实上一般需要LREAL=F的情况不是很多，除非对能量/电子结构计算要求特别精确或体系这些性质对数值比较...
1)根据手册描述,若体系原子大于 20 个,可以考虑使用 LREAL=Auto。 2)若关注的性质需要误差很小,比如极小的带隙 (meV 范围),磁各向异性能 (多在 μeV 范围内),建议使用 LREAL=.FALSE. 。 3)同样,若做分步优化,在第一步粗精度优化时,若体系原子大于 20 个,个人会考虑使用 LREAL=Auto。
LREAL: 决定投射是在实空间还是倒易空间进行。取值:.TRUE.(实空间)/.FALSE.(倒易空间),缺省值:.FALSE. 用于求解赝势的非局域部分用到的一个积分,在倒格空间里采用平面波基组求解,在实空间里则采用积分球求解。 其他还有两个选项:O or On,A or Auto。
LREAL=AUTO是在实空间中计算，速度会比.FALSE.快非常多。如果从结构优化开始就用AUTO，后面就一直用AUTO...
(一)LREAL= AUTO ;SIGMA=0.05 第一次得到的能量为-.16588630E+01 在OUTCAR中用/E-fermi搜索 可以看到电子占据是有问题的,不符合真实情况 (二)LREAL= .FALSE. ;SIGMA=0.05 第二次得到的能量为-.15453431E+01 在OUTCAR中用/E-fermi搜索 可以看到电子占据还是有问题的,但和(一)差别很大,这也一定程度上提...