声子谱计算实例一（计算Si的声子谱） (1).pdf

例子一 ：计算Si的声子色散谱 1 ）有限位移法 （又称直接法 ，冻声子法frozen-phonon法 ） 准备好INCAR ，POSCAR-UNITCELL （经过弛豫之后得到的稳定结构的 POSCAR ，也就是说已经完成了弛豫优化 ），KPOINTS ，POTCAR ， band.conf 这里给出我所准备的4个输入文件 ： INCAR START = 0 ICHARG = 2 ISIF = 2 NSW = 0 IBRION = -1 LWAVE = F LCHARG = F ENCUT = 500 EDIFF = 1E-8 #EDIFFG = -1E-3 不能设定EDIFFG参数为负值 ，一定要注意 ISMEAR = 0 LREAL = F SIGMA = 0.05 PREC = A POTIM = 0.35 #NPAR = 6 NELMIN = 4 POSCAR-UNITCELL Si                                         1.00000000000000        3.8669111948285479    0.0000000000000000     0.0000000000000000    1.9334555974142740    3.3488433288773667     0.0000000000002515    1.9334555974142740    1.1162811096595966     3.1573197693163806    Si     2 Direct     0.0000000000000000  0.0000000000000000  0.0000000000000000     0.2500000000000000  0.2500000000000000  0.2500000000000000 KPOINTS K-Mesh Generated with KP-Resolved Value  0 Gamma    3   3   3 0.0  0.0  0.0 POTCAR文件并不复杂 ，这里就不给出了 band.conf ATOM_NAME = Si DIM = 4 4 4 #FORCE_CONSTANTS = READ BAND = 0 0 0   0.5 0 0.5   0.5 0.25 0.75   0.375 0.375 0.75   0 0 0   0.5  0.5 0.5   0.625 0.25 0.625   0.5 0.25 0.75   0.5 0.5 0.5   0.375 0.375 0.75  #上面给出的是高对称点路径 ，注意上面从BAND开始到最后是同一行的内 容 ，中间不需要加enter换行 a)建立超胞 （需要多大的超胞 ，不同的材料不同的结构其实很难说一致 ，但 是基本上超胞越大计算结果得到虚频的可能性更小 ） phonopy -d --dim="4 4 4" -c POSCAR-UNITCELL 执行上述指令之后 ，可以得到SPOSCAR超胞文件 ，disp.yaml和POSCAR- {number}文件。 b ）计算力常数 mv POSCAR-001 POSCAR （这样必要输入文件INCAR、POSCAR、KPOINTS、POTCAR就已经准备 完毕了 ，通过静态计算 ，就可以得到vasprun.xml文件 ） phonopy -f vasprun.xml-* （得到FORCE_SETS文件 ） c ）得到声子谱 有了文件FORCE_SETS文件之后 ，结合前面准备的band.conf文件 ： phonopy --dim="4 4 4" -c POSCAR-UNITCELL -p band.conf （这样就可以得到band.yaml文件 ，从中可以提取得到声子谱能带数据 ） bandplot --gnuplot band.yaml > band.dat