问题：

vasp check_occupations: error system is not insulating basement

对于INCAR

cat > INCAR <<!
ISMEAR =  -5 !0
SIGMA  =  0.1 !0.05

GGA    = PE
EDIFF  = 1E-8

ALGO = scGW
NELM = 4
LSPECTRAL = .TRUE.
LOPTICS=.TRUE.
NOMEGA = 64
LRPA = .FALSE.
#PRECFOCK = Fast
#ENCUTGW = 700

NBANDS = 48
#NEDOS = 2000
#EMIN=-16.0
#EMAX=16.0

LOPTICS = .TRUE. ; LPEAD = .TRUE.

#VASP2WANNIER90
LWANNIER90=.TRUE.
!

网络摘录：

http://emuch.net/html/201406/7568893.html

计算极化强度的各位同行，我想请教一个用VASP计算极化强度的问题，铁电相是半导体可以算出极化强度，但是找到的参考相在计算的时候，有报错不是半导体或绝缘体，OUTCAR中得不到极化值，这时候该怎么办呢？而且看手册上说VASP是可以计算insulating system的，如果不是呢？需要改什么参数才能计算，ISMEAR我改过，但是没有用。如果大家遇到同样的问题，希望能交流

我前面用ISMEAR=0做静态和极化，但是铁电相也会报错有金属性，我就在做极化的时候改成-5，结果好了， 但是参考相没有什么作用，怎么改都是金属性。都想了两个月了也没有解决，有什么想法我们一起讨论吧

好的，这个是算极化的时候用的最基本 的INCAR：
SYSTEM = PKT
ENCUT=500
PREC=Accurate
ISMEAR =-5
SIGMA=0.05
EDIFF=1E-6
EDIFFG=1E-4
ISTART=0
ICHARG=2
IBRION=2
LCALCPOL=.TRUE.
DIPOL=0.3 0.3 0.3
说这个最基本是因为这个是算两层超晶格极化时用的，10个原子，如果再多一层或者N层的话，里边ISMEAR我需要改成0，不然经常出现错误；DIPOL我只是避开原子的坐标以及保证极化前后原子在0.3一侧就可以了，其他的我只设过DIPOL=0.125，目前还只是这两种情况；之于DIPOL不一样，电子和离子偶极矩差别很大这个我没专门对比过结果是否有很大出入，只是猜想 用那个公式加减处理后应该也差不多；还有就是多层的INCAR里加个LREAL；
下面这个是算LaAlO3/PbTiO3极化时INCAR，反复试了很多次，在网上搜的侯老师举得一个算Born的例子；能跑起来，只是我不知道怎么提取BERRY算法的结果，这个你知道应该grep 什么吗？
SYSTEM=PTOLAO
ENCUT=450
ISMEAR=0
SIGMA=0.002
EDIFF=1E-3
ISTART=1
ICHARG=11
IBRION=1
PREC=N
LBERRY=.TRUE.
NPPSTR=4
IGPAR=3
DIPOL=0.3 0.3 0.3
LDAU=.TRUE.
LDAUTYPE=2
LDAUL=-1 -1 3 -1 -1
LDAUU=0 0 6 0 0
LDAUJ=0 0 0 0 0
LDAUPRINT=2
静态计算的INCAR一般都没什么吧，把优化结构用的INCAR东西删一删，我一般留个PREC  ISMEAR  EDIFF ISYSTEM ,可能还加个LREAL就可以了，反正算后边的能带和DOS更重要吧