关注：

1） 电子结构如带隙的实验测量及理论计算方法  UPS、ARPES、GW

2） 理论计算过程中各参数的物理意义

3） 计算错误及原因分析

网络摘录：

http://blog.163.com/jjf_sxnu/blog/static/65117057201211994517175/

一、ERROR: charge density could not be read from file CHGCAR for ICHARG>10
And
ERROR: while reading WAVECAR, plane wave coefficients changed 23962, 11981

可能方案：sometimes an incorrect ISMEAR will cause this error；

please check whether the FFT meshes have changed；

the CHGCAR really is in the working directory at runtime；

the fft meshes of CHGCAR are compatible ；

More verbosely, grep the OUTCAR file from the SCF calculation and find the values of NGXF, NGYF, and NGZF. Set these parameters explicitly in the INCAR file that you use when calculating the band structure (the ICHARG=11 run).

2.  P4VASP导出并保存数据时须加后缀.dat

3. 如何绘制自旋向下和向上的能带图？  可从OUTCAR中的电子占据数判断。

实例：band计算OUTCAR文件分析

E-fermi :   6.3038     XC(G=0): -10.2027     alpha+bet :-14.3412

spin component 1

 k-point   1 :       0.0000    0.0000    0.0000   【gamma点？】
 band No.  band energies     occupation
     1     -39.9038      1.00000
     2     -39.8886      1.00000

    31      10.1345      0.00000
    32      10.2656      0.00000

k-point   2 :       0.0000    0.0000    0.0263  【到达下一个对称点之前经历的其他中间点，共18个中间点】
 band No.  band energies     occupation
     1     -39.9038      1.00000
     2     -39.8887      1.00000
     3     -13.1248      1.00000

k-point 140 :      -0.3330    0.6670    0.5000   【路径上的最后一个高对称点】
 band No.  band energies     occupation
     1     -39.8667      1.00000
     2     -39.8667      1.00000
     3     -13.0988      1.00000
     4     -13.0988      1.00000

...........

 spin component 2  【以下为自旋向下的情形】

k-point   1 :       0.0000    0.0000    0.0000
 band No.  band energies     occupation
     1     -37.2253      1.00000
     2     -37.2063      1.00000
     3     -11.2111      1.00000
     4     -10.5448      1.00000
     5     -10.5448      1.00000

......

k-point 140 :      -0.3330    0.6670    0.5000
 band No.  band energies     occupation
     1     -37.1794      1.00000
     2     -37.1794      1.00000
     3     -11.2095      1.00000
     4     -11.2095      1.00000
     5     -11.1972      1.00000

4.  关注固定原子位置优化情形

5.    CHECK_OCCUPATIONS: ERROR, system is not insulating.

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

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

附：

摘录：

如何判断 一个体系是否是 自旋极化体系，此外 很明显 上面的图不是自旋极化体系【从向上和向下的能带图是否重合判断？】

摘录2：

http://emuch.net/html/201305/5890125.html

如何用Vasp的结果画自旋向上和自旋向下分开的能带图

摘录3：

http://emuch.net/html/201002/1835984.html

是的，GW方法由于考虑了多体作用，不是单电子近似，是目前公认的计算能带宽度较准的方法之一。

6.PDOS