Dear all<div>I am using DFTtoolbox to draw fat band for Ni3SiTe8. However, I am puzzled. </div><div>The crystral structure is obtained from experiment, so I omit the procedure of doing structure relax. Than I submit the commands one by one </div><div><b>srun -n 128 pw.x <scf.in >scf.out</b></div><div><b>srun -n 128 pw.x <bands.in >bands.out</b></div><div><b>srun -n 1 bands.x <pp.bands.in >pp.bands.out  </b>(<i>Nb3SiTe8_bands.dat Nb3SiTe8_bands.dat.gnu Nb3SiTe8_bands.dat.rap</i> is obtained after this step while <i>Nb3SiTe8_bands.dat.rap</i> is an empty file)</div><div><b>srun -n 2 projwfc.x <projwfc.in >projwfc.out   </b>(<i>Nb3SiTe8_proj.dat.projwfc_up</i> is obtained after this step)</div><div>After that, I used the code provided by DFTtoolbox to draw my fatband plot. I used the postproc.py file given in the webside <a href="https://github.com/pipidog/DFTtoolbox/blob/master/build/lib/DFTtoolbox/postproc.py" target="_blank">DFTtoolbox/build/lib/DFTtoolbox/postproc.py at master · pipidog/DFTtoolbox · GitHub</a> I haven't made any change to the code. I changed some parameter of the file qe_pp.py this is my code (Nb3SiTe8_proj.dat.projwfc_up is the file I obtained in the fouth step)</div><div><div><b>import postproc</b></div><div><b>import os</b></div><div><b># Parameter ========================================================</b></div><div><b>run_task=[1,2,3,4]</b></div><div><b>wkdir=os.path.dirname(os.path.realpath(Nb3SiTe8_proj.dat.projwfc_up))</b></div><div><b># band_read & fatband_read</b></div><div><b>Ef=8.65</b></div><div><b>#band_plot</b></div><div><b>kdiv=[15,7,5,15,13,9,5,10,9,1]</b></div><div><b>klabel=['$\Gamma$','X','W','K','$\Gamma$','L','U','W','L','K']</b></div><div><b>Ebound=[-5,5]</b></div><div><b>#fatband_plot</b></div><div><b>state_grp=[['1:1/2/a/a'],['2:2/1/a/a']]</b></div><div><b># Main ================================================================</b></div><div><b>pp=postproc(wkdir)</b></div><div><b>for task in run_task:</b></div><div><b>    if task==1: #'band_read':</b></div><div><b>        pp.band_read(Ef=Ef,bandfile='pw.bands.out')</b></div><div><b>    elif task==2: #'band_plot':</b></div><div><b>        pp.band_plot(kdiv=kdiv,klabel=klabel,Ebound=Ebound)</b></div><div><b>    elif task==3: #'fatband_read':</b></div><div><b>        pp.fatband_read(Ef=Ef,projout='projwfc.fat.out',projprefix='fatband')</b></div><div><b>    elif task==4: #'fatband_plot':</b></div><div><b>        pp.fatband_plot(state_grp=state_grp,kdiv=kdiv,klabel=klabel,Ebound=Ebound)</b></div><div><b>    elif task==5: #pdos_read:</b></div><div><b>    <span style="white-space:pre">      </span>pp.pdos_read(Ef=Ef)</b></div><div><b>    elif task==6:</b></div><div><b>        pp.pdos_plot(state_grp=state_grp,Ebound=Ebound)</b></div></div><div><br></div><div>However, the fatband haven't been draw. I only obtained another python file <b>postproc.cpython-37.pyc</b> what's wrong? Where should I change? Or what should I write to obtain the fatband plot?</div><div>Could you please help me? Thank you very much!</div><div><br></div><div>Zongyi Wang </div>