VASP+Phonopy计算声子谱、偏态密度

1. 首先得到力FORCE_SETS/CONSTANTS

     VASP-INCAR文件-用于计算声子谱前的优化

       System = ***

          PREC = Accurate

          ENCUT = 400.0

          LREAL = .False

          ISMEAR = 0

          SIGMA = 0.05

          EDIFF = 1E-7

          EDIFFG = -1E-3

          NSW = 50

          ISIF = 3

          IBRION = 2

    VASP INCAR文件： DFPT in VASP- used to calculated phonon spectrum

          System = ***

          PREC = Accurate

          ENCUT = 400.0

          LREAL = .False

          ISMEAR = 0

          EDIFF = 1E-7

          NSW = 1

          IBRION = 8

          LWAVE = .FALSE

          LCHARG = .FALSE

          ADDGRID = .TRUE

      VASP POSCAR文件：cp SPOSCAR POSCAR（优化）

【应确保设置的KPOINTS产生的NKPTS个数达到20个左右。NKPT保持20个左右时，kmesh可取在0.04~0.05之间的某一数值，对计算结果影响不大，但相比0.03取值，可显著提高计算速度！！

优化完毕，check  OUTCAR文件中的频率，如果频率行包含i，则表明存在虚频；需修改参数继续优化，以检查是否是由于参数设置不当引起的虚频。如：grep 'cm-1' OUTCAR】

 对SPOSCR优化计算完成后，执行操作为： phonopy --fc vasprun.xml，得到FORCE_CONSTANTS

      VASP INCAR文件：冷冻原子法 used to calculated phonon spectrum

           ENCUT = 400.0

           LREAL = .False

           ISMEAR = 0

           EDIFF = 1E-7

           NSW = 0

           IBRION = 1

          LWAVE = .FALSE

           LCHARG = .FALSE

     VASP POSCAR文件：cp POSCAR-00i POSCAR（优化）

 对POSCR-00i静态计算完成后，执行操作为： phonopy --f $i/vasprun.xml得到FORCE_SETS

            （or phonopy-f disp-{001..003}/vasprun.xml）

2. 得到带结构所需文件及操作

    带结构所需文件：band.conf    FORCE_SETS/CONSTANTS POSCAR（扩胞前胞）    

         其中band.conf内容如下：

            ATOM_NAME = Ti S

           DIM = 2 3 2

           BAND =000 0.500 0.50.50 000 00.50 00.50.5 000 000.5

       BAND_LABELS= G Y C G  Z  D G B

       FORCE_CONSTANTS = WRITE

       BAND_CONNECTION = .TRUE.

    操作：1） phonopy -p --factor=521.471 band.conf  得到band.yaml

       >bandplot --gnuplot >Phononband.dat 将band.yaml转化为可用orgin绘图的数据文件

3. 得到声子谱态密度所需文件及操作

    所需文件： mesh.conf   FORCE_SETS    POSCAR（扩胞前的Unitcell）

         其中mesh.conf内容如下：

          ATOM_NAME = Ti S   （POSCAR原子分布 Ti Ti S S S S S S ）

         DIM = 2 3 2

      MP =20 20 20  （大一些）

      PDOS = 1 2, 3 4 5 6 7 8 （逗号隔开表示态密度按原子种类投影）

      FORCE_CONSTANTS = READ

   操作：phonopy -p --factor=521.471 mesh.conf，直接出来partial_dos.dat数据，用gnuplot作图

（1THz=4..413566meV http://halas.rice.edu/conversions）

plot "partial_dos.dat" u 1:($2+$3) w l lw 3 pt 4  ps 1.5 t 'Ti', "partial_dos.dat" u 1:($4+$5+$6+$7+$8+$9) w l lw 3 pt 8  ps 1.5 t 'S', "partial_dos.dat" u 1:($1+$2+$4+$5+$6+$7+$8+$9) w l lw 3 pt 6  ps 1.5 t 'Total'