<br>Dear All Pw-users !<br><br>I'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'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='scf'<br>    restart_mode='from_scratch',<br>    tstress = .true.<br>    tprnfor = .true.<br>
    prefix='810-28',<br>    pseudo_dir = '/home/projet/pseudo/',<br>
    outdir='/home/projet/tp/'<br> /<br> &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> &electrons<br>    diagonalization='david'<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>