本來(lái)一開(kāi)始打算用QE來(lái)算的，結(jié)果QE新版本bandup都不支持，banduppy例子都報(bào)錯(cuò)無(wú)語(yǔ)子，還是VASP吧。
他這里給了四步作為參考，但是實(shí)際上主要就是跑一個(gè)均勻網(wǎng)格點(diǎn)的自洽，和一個(gè)高對(duì)稱(chēng)點(diǎn)上的非自洽，然后做反折疊
首先第一步得到超胞CHGCAR
輸入文件參考BANDUP的example的第一步，提交任務(wù)之后得到CHGCAR。
然后第二步主要是為了得到超胞非自洽計(jì)算的K點(diǎn)路徑，這一步需要分別給出單胞和超胞的晶格坐標(biāo)文件prim_cell_lattice.in與supercell_lattice.in。還有單胞的高對(duì)稱(chēng)K點(diǎn)路徑
Kpoints along some high symmetry directions on the pcbz of Si (diamond structure)
23 27 9 29
Line-mode
Reciprocal
0.500 0.500 0.500 1 ! L
0.000 0.000 0.000 1 ! G
0.000 0.000 0.000 1 ! G
0.500 0.000 0.500 1 ! X
0.500 0.000 0.500 1 ! X
0.625 0.250 0.625 1 ! U
0.375 0.375 0.750 1 ! K
0.000 0.000 0.000 1 ! G
命令是bandup kpts-sc-get -no_symm_avg
這一步會(huì)生成KPOINTS_supercell.out文件后面會(huì)用到。
-no_symm_avg這個(gè)標(biāo)簽是為了減少需要計(jì)算的k點(diǎn)個(gè)數(shù)，減小計(jì)算量和WAVECAR大小。
然后是第三步
這一步利用前面第一步生成的電荷密度來(lái)計(jì)算高對(duì)稱(chēng)點(diǎn)的波函數(shù)
INCAR為
SYSTEM = Graphene
Flags which affect the runtime for parallel jobs (and don't change the results)
ALGO = Fast
LPLANE = .TRUE.
Start parameter for this Run:
ISTART = 0 ! job : 0-new 1-cont 2-samecut
ICHARG = 11 ! charge: 1-file 2-atom 10-const
LWAVE = T ! Write the WAVECAR file
LCHARG = F ! Do not write CHGCAR
Electronic Relaxation
EDIFF = 1E-7
Ionic Relaxation
IBRION = -1
ISIF = 2
NSW=0
Band Structure related values
ISMEAR = 0 ! Gaussian smearing
SIGMA = 0.025
NBANDS = 48
Precision
PREC = Normal
for consistency with static (total_energy) run:
NGX = 24; NGY = 24; NGZ = 24
NGXF = 48; NGYF = 48; NGZF = 48
ln -s ../step_1/CHGCAR CHGCAR
POSCAR POTCAR與第一步的一樣，K點(diǎn)從第二步那個(gè)復(fù)制過(guò)來(lái)
cp ../step_2_*/KPOINTS_supercell.out KPOINTS
然后提交任務(wù)得到波函數(shù)WAVECAR
第四步就要開(kāi)始做能帶反折疊了，這一步其實(shí)就是一句命令，但是例子沒(méi)有給它是怎么設(shè)置前邊得到的波函數(shù)的信息的。這句命令是
bandup unfold -emin -13 -emax 6 -dE 0.050 -no_symm_avg
會(huì)生成一個(gè)unfolded_EBS_not-symmetry_averaged.dat文件，畫(huà)出來(lái)就行了。
Quantum Espresso其實(shí)步驟也是差不多的，不同的地方主要在于首先bandup讀取波函數(shù)的方式還是QE老的格式，QE在6.4開(kāi)始采用了全新的存儲(chǔ)方式，在此之前的版本也必須編譯的時(shí)候加上一個(gè)標(biāo)簽
configure --disable-xml
而之后的版本有一個(gè)人開(kāi)發(fā)了banduppy聲稱(chēng)能夠解決問(wèn)題，我還沒(méi)試過(guò)。
另外不同的地方就在于第四步需要加一個(gè)指定prefix和outdir的標(biāo)簽
bandup unfold -qe -prefix TiSe2 -outdir ../step_3*/outdir -emin -2 -emax 2 -dE 0.005 -no_symm_avg
