<div dir="ltr"><div><div><div><div><div><div><div><div><div>Dear all<br><br></div>Hi,<br><br></div>This is my calculation to reproduce 3D bands structure of graphene.<br><br></div>1) I calculate calculation = 'scf', and then calculation = 'bands', as following calculation at end of this message <br>
</div>2) I read bands by <a href="http://bands.in">bands.in</a> as following calculation at end of this message </div><br>However; <br><br></div>I face three problems:<br></div>1) the results that I plotted by GNU plot is totally different by calculation of tight-binding and it is meaningless for me. However I do not have any problem two plot 2D band structure<br>
<br></div>2) I would like to produce the mesh grid of Kx and Ky by myself; but whenever I try to do this I input Kx and Ky by myself QE considers it as path and not as mesh grid . <br><br></div><div>3) I try two reproduce Kx and Ky by kpoint.x in this case the QE considers it as path in BZ again not as mesh grid.(I changed the automatic to crystal the only way that QE works )<br>
<br></div><div>I think my problems is I cannot reproduce a proper mesh grid for 3D calculation by QE.<br></div><div>I deeply appreciate you if you help me to understand and plot the 3D bands structure in QE. Thanks<br></div>
<br><div><div><div><div><div><div>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br><div> &CONTROL<br>       calculation = 'scf',<br>      restart_mode = 'from_scratch',<br>        pseudo_dir = './pseudo/',<br>
            outdir = './tmp/',<br>            prefix = 'slg',<br>          <br> /<br> &SYSTEM<br>             ibrav = 4,<br>                 a = 2.439,<br>                 c = 10<br>               nat = 2,<br>
              ntyp = 1,<br>           ecutwfc = 40,<br>              nosym=.TRUE.<br>              noinv=.TRUE.<br><br><br> /<br> &ELECTRONS<br><br> /<br>ATOMIC_SPECIES<br>C 12.0107 C.pz-rrkjus.UPF<br>ATOMIC_POSITIONS {crystal}<br>
C  0.333333333  0.666666666  0.500000000<br>C  0.666666666  0.333333333  0.500000000<br>K_POINTS {automatic}<br>  42 42 1   0 0 0<br><br>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br><br> &bands<br>    prefix  = 'slg'<br>
    outdir='./tmp'<br>    filband = 'bands.dat'<br>    plot_2d = .TRUE.<br> /<br>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br><br></div><div>Thank you !<br><br></div><div>Pourya<br></div></div>
</div></div></div></div></div></div>