结构优化

首先利用建好结构，转化为分数坐标格式，命名为POSCAR
准备好四个文件：INCAR，POSCAR，POTCAR，KPOINTS
提交任务

#====================================================
!................................................
ISTART = 0              !存在wavecar时是1，否则是0
ICHARG = 2              !0，从初始波函数计算电荷密度，1，从ACHRGCAR文件读入，2，用原子电荷，密度组装
NSW = 300               !原子弛豫的步数，优化时设置，自洽时原子不动，所以为0
!.................................................
ENCUT = 400             !平面波基的截断能
ISMEAR = 0              !-5：bloch正四面体法，计算大块总能时使用，-4：不采用Bloch正四面体法，-1：费米，0：高斯展开        
SIGMA = 0.05
ISPIN = 2               !指电子体系的总自旋
IBRION = 2              !0，分子动力学模拟 1，准牛顿离子弛豫，适合初始结构和平衡结构接近的情况 2，CG算法的离子弛豫 
ISIF = 4                !改变元胞形状，不改变体积
!...................................................
EDIFF = 1E-8            !总能计算中的允许误差，默认值是1E-4
EDIFFG = -0.001         !离子弛豫运动的结束调节，默认是EDIFF*10,如果EDIFFG值为正数，是指各原子上的力之和；如果为负数表示每个原子上的力。和上一个参数是且的关系，两个都满足才行
LORBIT = 11             !控制是否输出投影波函数
POTIM = 0.4             !IBRION=0时，为分子动力学离子运动时间步长，IBRION=1，2，3时，为作用在力上的比例系数
LREAL=Auto              !投影操作在实空间进行，对投影操作进行了全优化
!.................................................
PREC = Accurate         !比较长，见中文说明书12页

LPLANE = .TRUE.         !与并行计算有关
LSCALU = .FALSE.        !没找到
ALGO = N                !选择电子优化算法，与IALGO=38相同，选择Davidson优化算法
ISYM = 2                !是否具有对称性，1具有，0不具有，2内存对电荷密度的处理更有效
!..................................................
LSORBIT = T    !     soc
GGA_COMPAT = F !    !!!!!
SAXIS = 0 0 1  !    !!!!!
!.................................................
#IVDW=12                #范德华修正
##.................................................
#LDIPOL=.TRUE.
#IDIPOL=3
#EFIELD=0.5
#DIPOL=0.5 0.5 0.5    #电场
#...................................................
LMAXMIX = 4             !d电子为4，f电子为6
LDAU=T                  !这两项不要改
LDAUTYPE=2              !
LDAUL= 2 -1 -1
LDAUU= 4.4 0.0 0.0      !加u的那个值对应的不为0
LDAUJ= 0.0 0.0 0.0
LASPH=T
#===================================================

在计算二维体系时，ISIF设置为4，不能让真空层消失，如果ISIF设置为3，需要加OPTCELL锁住c方向的晶格矢量。
OPTCELL设置如下：

1 1 0
1 1 0
0 0 0

常见错误
（1）ERROR: charge density could not be read from file CHGCAR for ICHARG>10：
解决：是由于与自洽的prec不一致造成的，修改PREC
（2）LAPACK: Routine ZPOTRF failed! 1 1 1
解决：ISIF=3时就出现报错，需要加上OPTCELL或者ISIF=4
（3）Sub-Space-Matrix is not hermitian in DAV
解决：减小POTIM 让IBRION=0，标准分子动力学模拟，通过POTIM控制步长
（4）WARNING: Sub-Space-Matrix is not hermitian in DAV
解决：此时换一种弛豫，加上IALGO＝48以及IALG=Fast ，减小potim
	  如果在自洽时出现，将PREC改为Accurate
（5）VERY BAD NEWS! internal error in subroutine IBZKPT:
 Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful...
方法一：ISYM=0
方法二：SYMPREC=1E-8
注：系统有一定的对称性但不太好，但程序的判定条件较低 (SYMPREC=1E-5，缺省？)而认为有对称，
如果设置上面两项，则需要的k点增多！在INCAR中添加了SYMPREC=1E-8问题得到解决。
（6）VERY BAD NEWS! internal error in subroutine INVGRP:
 inverse of rotation matrix was not found (increase SYMPREC)       9
 增大SYMPREC