Thank you Davide for the clarification.<br><br><div class="gmail_quote">On Wed, Dec 16, 2009 at 6:30 PM, Davide Ceresoli <span dir="ltr"><<a href="mailto:ceresoli@mit.edu">ceresoli@mit.edu</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;">
<div class="im">On Wed, 16 Dec 2009, mohamed sabri majdoub wrote:<br>
<br>
> Dear Davide,<br>
><br>
> If I change celldim(1)=1.8897261 (in Bohr radius unit), this means that the crystal lattice<br>
> vectors are in Bohr too (same unit as alat).<br>
</div>Yes, the lattice vectors must be bohrradius unit.<br>
<div class="im"><br>
> Crystal lattice vectors:<br>
>     v1(1)  v1(2)  v1(3)    ... 1st lattice vector<br>
><br>
>     v2(1)  v2(2)  v2(3)    ... 2nd lattice vector<br>
>     v3(1)  v3(2)  v3(3)    ... 3rd lattice vector<br>
> In alat units if celldm(1) was specified or in a.u. otherwise.<br>
> Am I right? If yes then I need to change the following to Bohr units :<br>
><br>
> CELL_PARAMETERS<br>
> 8.51980 0.00000 0.00000<br>
> 0.00000 7.37600 0.00000<br>
> 0.00000 0.00000 22.6120<br>
</div>You can change them to bohr units and set celldm(1)=1.0, or<br>
you leave them unchanged but set celldm(1)=1.8897261.<br>
<br>
In any case, the product CELL_PARAMETERS*celldm(1) must always<br>
match your wanted lattice spacings in bohr units, not in angstrom.<br>
<br>
This is in contrast with ATOMIC_POSITIONS, where one can choose<br>
to enter them in angstrom/bohr/alat/crystal units.<br>
<br>
HTH.<br>
<font color="#888888"><br>
Davide<br>
</font><div><div></div><div class="h5"><br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">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>
</div></div></blockquote></div><br>