<style class="ke-style">
[list-style-type] {padding-left:20px;list-style-position:inside}
[list-style-type] li {margin:0}
[list-style-type] li:before, span.ke-list-item-matter {font-family:"sans serif",tahoma,verdana,helvetica}
[list-style-type] li p,[list-style-type] li h1,[list-style-type] li h2,[list-style-type] li h3,[list-style-type] li h4,[list-style-type] li h5,[list-style-type] li div,[list-style-type] li blockquote{display:inline;word-break:break-all}
[list-style-type] li table {display:inline-block;vertical-align:top}
p{margin:0}
td {word-break: break-word}
.default-font-1752467044080{
}
</style><div class="default-font-1752467044080" dir="ltr"><p style="font-size:14px;">Dear all</p><p style="font-size:14px;">  I am trying to do band unfolding calculation, but after I read the tutorial online, I only find the auxiliary software unfold-x provided by <span style="font-size:14px;">Dr.Pietro Bonfa’.I followed the tutorial and calculated my system step by step but I'm stopped by some problem.</span></p><p style="font-size:14px;"><span style="font-size:14px;">  </span></p><p style="font-size:14px;"><span style="font-size:14px;">  Firstly, I calculated the primitive cell in the following sequence</span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;"><em><strong>mpirun -np 32 pw.x < scf.in > scf.out</strong></em></span></span></p><p style="font-size:14px;"><strong><em><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;">mpirun -np 32 pw.x < bands.in > bands.out</span></span></em></strong></p><p style="font-size:14px;"><strong><em><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;">mpirun -np 1 bands.x < pp.bands.in > pp.bands.out</span></span></em><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;"></span></span></strong><span style="font-size:14px;"></span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;">Then I obtained <span style="text-wrap-mode:nowrap;">VSe2_bands.dat.gnu</span> which describe the band structure in its' two columns. After that I followed the tutorial for the supercell calculation. The supercell I choose is 2*2*1. I calculated by the following step:</span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;"><em><strong>mpirun -np 32 pw.x < scf.in > scf.out</strong></em></span></span></p><p style="font-size:14px;"><strong><em><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;">mpirun -np 1 unklist.x < unklist.in > unklist.out</span></span></em></strong></p><p style="font-size:14px;"><strong><em><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;">mpirun -np 32 pw.x < bands.in > bands.out</span></span></em></strong></p><p style="font-size:14px;"><strong><em><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;">mpirun -np 1 unfold.x < unfold.in > unfold.out</span></span></em><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;"></span></span></strong><span style="font-size:14px;"></span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;">After that, the file <span style="text-wrap-mode:nowrap;">VSe2_Unfolded_SPFN.dat</span> generated successfully. But when I use gnuplot to draw the picture, something went wrong. My gnuplot code is in the following lines </span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;"><span style="text-wrap-mode:nowrap;"><em><strong>set terminal pngcairo enhanced font "Verdana,10"</strong></em></span></span></p><p><strong><em><span style="text-wrap-mode:nowrap;">set output 'output.png'</span></em></strong></p><p><strong><em><span style="text-wrap-mode:nowrap;">TPIBA=1.1845361</span></em></strong></p><p><strong><em><span style="text-wrap-mode:nowrap;">set pm3d map interpolate 2,2</span></em></strong></p><p><strong><em><span style="text-wrap-mode:nowrap;">set yrange [-3.5:0.5]</span></em></strong></p><p><strong><em><span style="text-wrap-mode:nowrap;">#splot 'VSe2_Unfolded_SPFN.dat' binary record=(121,-1) format='%double' u ($1/TPIBA):2:3 notitle</span></em></strong></p><p><strong><em><span style="text-wrap-mode:nowrap;">splot 'VSe2_Unfolded_SPFN.dat' binary record=(121,-1) format='%double' u ($1/TPIBA):2:3 notitle, '../KGM/VSe2_bands.dat.gnu' u 1:2:2 w p ps 0.1 notitle</span></em></strong></p><p><strong><em><span style="text-wrap-mode:nowrap;">set output</span></em></strong><em><span style="text-wrap-mode:nowrap;"></span></em><span style="text-wrap-mode:nowrap;"></span></p><p><br></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span><br></p><p style="font-size:14px;"><span style="font-size:14px;">and the output png file is like this </span></p><p style="font-size:14px;"><span style="font-size:14px;"> </span></p><img src="cid:30def768$1$1980742b0f6$Coremail$wangzongyi$mail.ustc.edu.cn" alt="" _owidth="640" _oheight="480" width="640px" height="480px; width: 640px"><p><br></p><p style="font-size:14px;"><span style="font-size:14px;">It shows that the two data file are not compatible and the figure is really terrible.</span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;">What should I change? Could you please help me?</span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;">Thank you very much!</span></p><p style="font-size:14px;"><span style="font-size:14px;"><br></span></p><p style="font-size:14px;"><span style="font-size:14px;">Zongyi Wang</span></p></div>