<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Is your question "Why do you get non-zero force for this fcc system with your given Cell_parameters and Atomic_positions?"</span></font></div>

<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">If you have a perfect fcc Cu structure, you should get 0 force. I think there are some numerical errors, it likes <<<2*0.5*<span class="Apple-style-span" style="font-size: 13px; ">6.81650937063832 .NE.6.81650937063832>>>. This kind of error "moves" the atoms off their ideal positions. <span class="Apple-style-span" style="font-size: small; ">I guess if you use crystal coordinate instead of Cartesian coordinate, you will get 0 force.</span></span></span></font></div>

<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">It's just my humble thought, I don't know for sure how pwscf read and process your coordinates.</span></font></div>

<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse; "><br></span></font></div>--------------------------------------------------<br>Duy Le<br>PhD Student<br>

Department of Physics<br>University of Central Florida.<br>
<br><br><div class="gmail_quote">On Wed, Oct 7, 2009 at 5:23 PM, Manoj Srivastava <span dir="ltr"><<a href="mailto:manoj@phys.ufl.edu">manoj@phys.ufl.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<div id=":471" class="ii gt">ATOMIC_POSITIONS<br>
  Cu  0.0   0.0   0.0<br>
  Cu  0.5   0.5   0.0<br>
  Cu  0.0   1.0   0.0<br>
  Cu  1.0   1.0   0.0<br>
  Cu  0.5   1.5   0.0<br>
<br>
K_POINTS (automatic)<br>
8 8 8 0 0 0<br>
CELL_PARAMETERS {cubic}<br>
 1.5   0.5  0.0<br>
-0.5   1.5  0.0<br>
 0.0   0.5  0.5</div></blockquote></div><br>