<div dir="ltr"><div>Dear QE users,</div><div><br></div><div>I am working with a hexagonal symmetry crystal and I relaxed all internal coordinates and cell parameters with vc-relax-bfgs algorithm.</div><div><br></div><div>I have some concern about my results. So, I would like some advice.  </div>
<div><br></div><div>My initial cell parameter are:</div><div><br></div><div><br></div><div>45.48740530000000000<span class="" style="white-space:pre">        </span>0.00000000000000000<span class="" style="white-space:pre">       </span>0.00000000000000000</div>
<div>-22.74370265000000000<span class="" style="white-space:pre">       </span>39.39324854203890000<span class="" style="white-space:pre">      </span>0.00000000000000000</div><div>0.00000000000000000<span class="" style="white-space:pre">     </span>0.00000000000000000<span class="" style="white-space:pre">       </span>10.31035723491920000</div>
<div><br></div><div>with vectores  'a'  and 'b', as:</div><div><br></div><div>a = b = 45.48740530000000000</div><div>c = 10.31035723491920000</div><div><br></div><div>After optimization I got this value:</div>
<div><br></div><div>44.53851082397360000<span class="" style="white-space:pre">       </span>0.00000000000000000<span class="" style="white-space:pre">       </span>0.00000000000000000</div><div>-22.26925541198680000<span class="" style="white-space:pre">   </span>38.57148182127630000<span class="" style="white-space:pre">      </span>0.00000000000000000</div>
<div>0.00000000000000000<span class="" style="white-space:pre"> </span>0.00000000000000000<span class="" style="white-space:pre">       </span>10.12292974744850000</div><div><br></div><div>In my calculation all operations symmetry  are conserved. However if use the cell parameters as:</div>
<div><br></div><div>a = 44,53851082397360000</div><div>b = 44,53851082482830000</div><div>c = 10,12292974744850000</div><div><br></div><div>The 'a'  and 'b' parameter aren't exactly the same.</div><div>
<br></div><div>I did a test, where I use</div><div><br></div><div>ibrav = 4</div><div>celldm(1)=  44.53851082397360000,</div><div>celldm(3) = 0.227284872 ,</div><div><br></div><div>The component forces of atoms become high than  (>10-3). </div>
<div><br></div><div>To keep the component forces of atoms less than >10-3,  I have to use ibrav = 0 and specify the matriz of cell parameters. </div><div><br></div><div>With this process, Can I still consider my cell as hexagonal symmetry?</div>
<div><br></div><div>Cheers</div><div><br></div>-- <br><div dir="ltr">José Gadelha da Silva Filho<div><br></div><div>Aluno de Doutorado do Curso de Pós-Graduação em Física</div><div>da Universidade Federal do Ceará.</div></div>

</div>