Thank you so much dear Lorenzo...<br><br><div class="gmail_quote">On Tue, Oct 27, 2009 at 9:32 PM, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:paulatto@sissa.it" target="_blank">paulatto@sissa.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
In data 27 ottobre 2009 alle ore 15:00:25, mohnish pandey<br>
<<a href="mailto:mohnish.iitk@gmail.com" target="_blank">mohnish.iitk@gmail.com</a>> ha scritto:<br>
<div>> I can fix this problem by preserving the symmetry and<br>
> getting eqb. lattice constan.?<br>
<br>
</div>Dear Mohnish,<br>
the code actually do preserve the symmetry, it is enforced actively when<br>
computing the forces and cannot be broken during the relax procedure.<br>
What's happening is that the system has a lower symmetry than you think<br>
since the very beginning.<br>
<br>
There are two possible reasons:<br>
1. you have not specified the initial atomic positions with sufficient<br>
precision, the code consider an atom to be in a high symmetry position if<br>
it is less than 1.d-6 away from it (I think in units of alat=celldm(0),<br>
but I would have to check)<br>
2. the symmetry of your system has a free parameter you where not aware<br>
of, e.g. in a Wurzite structure the distance between alternate planes can<br>
change in a specific way without breaking symmetry.<br>
<br>
These two problems have two corresponding solutions:<br>
1. specify more digits, or (with QE 4.1 and next) specify the position as<br>
a fraction, e.g.<br>
   Ga 1/3 2/3 0<br>
note that xcrysden cannot yet interpret an input file written like this<br>
and it is still a bit experimental<br>
2. lock the atoms in position, you can do that adding three zeroes after<br>
the atom's coordinate in the input, e.g.:<br>
  Ga  0.3333333 0.6666667 0.000000  0 0 0<br>
have a look at the manual (Doc/INPUT_PW.txt) for details on details on<br>
this notation<br>
<br>
best regards<br>
<br>
<br>
<br>
--<br>
Lorenzo Paulatto<br>
SISSA  &  DEMOCRITOS (Trieste)<br>
phone: +39 040 3787 511<br>
skype: paulatz<br>
www:   <a href="http://people.sissa.it/%7Epaulatto/" target="_blank">http://people.sissa.it/~paulatto/</a><br>
<br>
     *** save italian brains ***<br>
  <a href="http://saveitalianbrains.wordpress.com/" target="_blank">http://saveitalianbrains.wordpress.com/</a><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Mohnish Pandey<br>Y6262,4th Year Undergraduate,<br>Department of Chemical Engineering,<br>IIT KANPUR<br>