<br>Dear All Pw-users !<br><br>I&#39;m a new pw user. <br>I want to calculate total energy of B-GaSe  crystal (hexagonal structure). <br>Altough I tried so hard, I couldn&#39;t set up unit cell of B-GaSe in Pwscf.<br>I couldn t see correct atomic structure of GaSe. I don t know where ư am wrong.<br>
Can anybody help me ! please ! <br>I used parameters in the attached files. Also, my input file follow as <br>
<br>(lattice parameters a=3.73A  , c=15.919A,  z1=0.1003 and z2=0.1701<br>  <br><br>control<br>    calculation=&#39;scf&#39;<br>    restart_mode=&#39;from_scratch&#39;,<br>    tstress = .true.<br>    tprnfor = .true.<br>
    prefix=&#39;810-28&#39;,<br>    pseudo_dir = &#39;/home/projet/pseudo/&#39;,<br>
    outdir=&#39;/home/projet/tp/&#39;<br> /<br> &amp;system<br>    ibrav=0, <br>    celldm(1)=7.051,  <br>    nat=8, <br>    ntyp=2,<br>    ecutwfc =50.0, <br>    ecutrho = 600.0,<br> /<br> &amp;electrons<br>    diagonalization=&#39;david&#39;<br>

    conv_thr = 1.0e-6<br>    mixing_beta = 0.7<br> /<br>ATOMIC_SPECIES<br> Se  1.0 Se_pz.upf<br> Ga  1.0 Ga_pz.upf<br>ATOMIC_POSITIONS {alat}<br>Se -0.5     0.28867  0.1003<br>Ga  0.0     0.0      0.1701<br>Ga  0.0     0.0      0.3299<br>

Se -0.5     0.28867  0.3997<br>Se  0.0     0.0      0.6003<br>Ga -0.5     0.28867  0.6701<br>Ga -0.5     0.28867  0.8299<br>Se  0.0     0.0      0.8997 <br>K_POINTS {automatic}<br>6 6 6 0 0 0<br><br>CELL_PARAMETERS {hexagonal}<br>

 1.0   0.0   0.0<br> 0.5   0.866 0.0<br> 0.0   0.0   4.26<br>