1.自洽計(jì)算

復(fù)制計(jì)算后的文件到當(dāng)前文件夾
cp …/phonon/pb.save/charge-density.dat pb.save/
cp …/phonon/pb.save/data-file-schema.xml pb.save/

2.band.in

/pw.x -nk 4 -in pb.band.in > pb.band.out &control calculation='bands' restart_mode='from_scratch', prefix='pb', pseudo_dir = '../', outdir='./' / &system ibrav= 2, celldm(1) = 9.2225583816, nat= 1, ntyp= 1, ecutwfc = 30.0 occupations='smearing', smearing='marzari-vanderbilt', degauss=0.05 nbnd=10 / &electrons conv_thr = 1.0d-10 mixing_beta = 0.7 / ATOMIC_SPECIES Pb 207.2 pb_s.UPF ATOMIC_POSITIONS Pb 0.00 0.00 0.00 K_POINTS crystal_b #這里是高對(duì)稱路徑 7 0.000 0.000 0.000 30 0.500 0.000 0.500 30 0.500 0.250 0.750 30 0.500 0.500 0.500 30 0.375 0.375 0.750 30 0.000 0.000 0.000 30 0.500 0.500 0.500 30

注：nbnd選擇
以scf.out中這個(gè)量為基準(zhǔn)

number of electrons = 8.00

絕緣體取1/2,金屬在絕緣體的基礎(chǔ)上加20%，加空帶會(huì)拖慢計(jì)算速度，但有利于完整能帶結(jié)構(gòu)。
PS.這個(gè)量是怎么來(lái)的呢？
贗勢(shì)文件打開(kāi)搜索Z valence
所有原子的價(jià)電子加起來(lái)就是它。

3.bands.in

/bands.x -in pb.bands.in > pb.bands.out &bands prefix = ’pb’ outdir = ’out/’ filband = ’bands.dat’ lsym = .false. /

4.畫(huà)圖

qe/plotband.x

直接在主節(jié)點(diǎn)運(yùn)行，回答它的問(wèn)題

$ ~/soft/q-e-qe-6.7MaX-Release/bin/plotband.x Input file > bands.dat Reading 10 bands at 251 k-points file with representations not compatible with bands Range: -25.4690 21.7530eV Emin, Emax, [firstk, lastk] > -26,22 high-symmetry point: 0.0000 0.0000 0.0000 x coordinate 0.0000 high-symmetry point: 0.0000 0.2802 0.0000 x coordinate 0.2802 high-symmetry point: 0.1401 0.1401 0.0000 x coordinate 0.4784 high-symmetry point: 0.1401 0.1401 0.1401 x coordinate 0.6185 high-symmetry point: 0.0000 0.0000 0.0000 x coordinate 0.8612 high-symmetry point: 0.1401 0.1401 0.0000 x coordinate 1.0594 output file (gnuplot/xmgr) > plot.band.dat bands in gnuplot/xmgr format written to file CaH6.band.dat

最后把plot.band.dat拖去畫(huà)圖，注意，y列數(shù)據(jù)要減去費(fèi)米能的值
費(fèi)米能的值在scf.out里找

the Fermi energy is 14.5455 ev

大功告成

5. FatBand 分析

https://blog.csdn.net/weixin_45973682/article/details/108194472

projwfc.fat.in

&projwfcoutdir='./'prefix='pwscf'lsym=.false.filproj = 'fatband'/

運(yùn)行

projwfc.x < projwfc.fat.in > projwfc.fat.out

畫(huà)圖 復(fù)制下面文件到計(jì)算文件夾

DFTtoolbox-master/examples/qe/soc/qe_pp.py from DFTtoolbox.qe import postproc import os # Parameter ======================================================== run_task=[3,4] ！選擇計(jì)算內(nèi)容，畫(huà)胖帶分析就選3，4 wkdir=os.path.dirname(os.path.realpath(__file__)) # band_read & fatband_read Ef=17.2829 ！scf.out里有，這是Fermi能 #band_plot kdiv=[50,50,50,50,50,50] ！這是高對(duì)稱點(diǎn)插值，和band.in保持一致 klabel=['$\Gamma$','H','N','P','$\Gamma$','N'] Ebound=[-26,10] ！這是你想畫(huà)的圖的范圍，和文獻(xiàn)對(duì)應(yīng) #fatband_plot state_grp=[[1],[2]] 這個(gè)量的選擇看后面，詳細(xì)說(shuō) # Main ================================================================ pp=postproc(wkdir) for task in run_task:if task==1: #'band_read':pp.band_read(Ef=Ef,bandfile='bands.dat')elif task==2: #'band_plot':pp.band_plot(kdiv=kdiv,klabel=klabel,Ebound=Ebound)elif task==3: #'fatband_read':pp.fatband_read(Ef=Ef,projout='proj.fat.out',projprefix='fatband')elif task==4: #'fatband_plot':pp.fatband_plot(state_grp=state_grp,kdiv=kdiv,klabel=klabel,Ebound=Ebound)elif task==5: #pdos_read:pp.pdos_read(Ef=Ef)elif task==6:pp.pdos_plot(state_grp=state_grp,Ebound=Ebound)

注意！使用的贗勢(shì)要有<PP_PSWFC>部分才可以
**【state_grp】**的選擇
先隨便寫(xiě)運(yùn)行一下
會(huì)得到下面輸出

=> check fatband-state.dat for state_info: [' 0 => H-1 orb-1 ( 1 / 0 / 1 / +1)\n'' 1 => H-2 orb-1 ( 2 / 0 / 1 / +1)\n'' 2 => H-3 orb-1 ( 3 / 0 / 1 / +1)\n'' 3 => S-4 orb-1 ( 4 / 0 / 1 / +1)\n'' 4 => S-4 orb-2 ( 4 / 1 / 1 / +1)\n'' 5 => S-4 orb-2 ( 4 / 1 / 2 / +1)\n'' 6 => S-4 orb-2 ( 4 / 1 / 3 / +1)\n']

我這是H3S的例子，H的贗勢(shì)文件中PP_PSWFC里有1s軌道，S的贗勢(shì)文件中PP_PSWFC里有3s，3p軌道。
所以，你想看H對(duì)能帶的貢獻(xiàn)就state_grp=[[0,1,2]]

想看S的3s軌道就state_grp=[[3]]，想看S的3p軌道就state_grp=[[4,5,6]]
當(dāng)然，這三張圖可以寫(xiě)在一起state_grp=[[0,1,2],[3],[4,5,6]]

總結(jié)

以上是生活随笔為你收集整理的使用QE绘制Band Structure的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。