<div dir="ltr">Hello,<div><br></div><div>I am a beginner in using Quantum Espresso. I have a 64 atom Si supercell for which I am trying to plot band diagrams. I have gone through the tutorial "PWSCF: First Steps" available at <a href="http://www.fisica.uniud.it/~giannozz/QE-Tutorial/tutorial_pwscf_ex.pdf">http://www.fisica.uniud.it/~giannozz/QE-Tutorial/tutorial_pwscf_ex.pdf</a> .</div>
<div>The problem I am facing is that when I run a scf calculation using 4*4*4 k-points, my output file has energies of 128 bands for every k-point. For example,</div><div><div><b>k = 0.0000 0.0000 0.0000 ( 12797 PWs) bands (ev):</b></div>
<div><b><br></b></div><div><b> -5.6509 -4.5410 -4.5410 -4.5410 -4.5410 -4.5410 -4.5410 -3.5502</b></div><div><b> -3.5502 -3.5502 -3.5502 -3.5502 -3.5502 -3.5502 -3.5502 -3.5502</b></div><div><b> -3.5502 -3.5502 -3.5502 -3.2434 -3.2434 -3.2434 -3.2434 -1.3979</b></div>
<div><b> -1.3979 -1.3979 -1.3978 -1.3978 -1.3978 -1.2206 -1.2206 -1.2206</b></div><div><b> -1.2206 -1.2206 -1.2206 -1.2205 -1.2205 -1.2205 -1.2205 -1.2205</b></div><div><b> -1.2205 -0.5951 -0.5951 -0.5951 -0.5951 0.9797 0.9797 0.9797</b></div>
<div><b> 0.9797 0.9797 0.9797 0.9797 0.9797 0.9797 0.9797 0.9797</b></div><div><b> 0.9797 2.5385 2.5385 2.5385 2.5385 2.5385 2.5385 2.5385</b></div><div><b> 2.5385 2.5385 2.5385 2.5385 2.5385 2.6504 2.6504 2.6504</b></div>
<div><b> 2.6504 2.6504 2.6504 2.6504 2.6504 2.6504 2.6504 2.6504</b></div><div><b> 2.6504 3.0031 3.0031 3.0031 3.0031 3.0031 3.0031 3.5849</b></div><div><b> 3.5849 3.5849 3.5849 3.5849 3.5849 4.5713 4.5713 4.5713</b></div>
<div><b> 4.5713 4.5713 4.5713 4.5713 4.5713 4.5713 4.5713 4.5713</b></div><div><b> 4.5713 5.1567 5.1567 5.1567 5.1567 5.1568 5.1568 5.1568</b></div><div><b> 5.1568 5.1568 5.1568 5.1568 5.1568 5.2932 5.2932 5.2932</b></div>
<div><b> 5.2932 5.2932 5.2932 5.2932 5.2932 6.5240 6.5240 6.5240</b></div></div><div><br></div><div>On doing calculation = 'bands' and nbnd = 8, the output file contains energies of only 8 bands for every k-point which are the lowermost bands. </div>
<div><div><b>k = 0.0000 0.0000 0.0000 band energies (ev):</b></div><div><b><br></b></div><div><b> -5.6509 -4.5410 -4.5410 -4.5410 -4.5410 -4.5410 -4.5410 -3.5502</b></div></div><div>Then I run a bands.x executable to genrate bands.dat. On using plotband.x, I get the following.</div>
<div><br></div><div><div><b>/opt/apps/espresso-5.0.1/bin/plotband.x</b></div><div><b>Input file > bands.dat</b></div><div><b>Reading 8 bands at 10 k-points</b></div><div><b>Range: -5.6510 -3.5500eV Emin, Emax > -6, -3</b></div>
</div><div><b>high-symmetry point: 0.0000 0.0000 0.0000 x coordinate 0.0000<br></b></div><div><b>...</b></div><div><b>high-symmetry point: -0.5000-0.5000-0.5000 x coordinate 0.0000<br></b></div><div><div><b>output file (xmgr) > si.bands.xmgr</b></div>
<div><b>bands in xmgr format written to file si.bands.xmgr </b></div>
<div><b>output file (ps) > <a href="http://si.bands.ps">si.bands.ps</a></b></div><div><b>Efermi > -3.987 </b></div><div><b>deltaE, reference E (for tics) 1.0, -3.987</b> </div></div><div><div><b>n= 2 3</b></div>
<div><b> 0.0000000E+00 0.0000000E+00</b></div><div><b> -5.651000 -5.581000 </b></div><div><b> NaN NaN </b> </div></div><div>....</div><div>....</div><div>....</div><div>and more NaNs and my <b><a href="http://si.bands.ps">si.bands.ps</a></b> file contains no plot.</div>
<div><br></div><div>I would be extremely thankful if someone could tell me where am I wrong and what different can I do to obtain the band plot (4 valence bands and 4 conduction bands).</div><div><br></div><div>Thanks,</div>
<div>Mitul Mundra</div><div>Final Year Dual Degree Student</div><div>Department of Chemical Engineering</div><div>IIT Kanpur, India.</div><div><br></div></div>