<p>Dear wien_2k users</p>
<p>I want to do ion relaxation for PrBa2Cu3O7 system. I used NEWT methods and this command: min -j'run_lapw –I –i 99 –fc 1'.</p>
<p>In case_n.scf (n=1,2,3,…) the forces have very large fluctuation in final iteration. For example the forces of 4.ATOM in case_3.scf are :</p>
<p>Iteration: 12&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;13&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14</p>
<p>FOR004: 33.467&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 33.079&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;3.214</p>
<p>(14 is the final iteration)</p>
<p>I have one important question .When –fc 1 was used, why after a large fluctuation of forces in final iteration of scf n (n=1,2,3,…) (about 30 mRy/au &gt;&gt;1)the program converged (and went to n+1 step)?</p>
<p>Thanks</p>
<p>Vahid Ghanbarian.</p>
<p>University of Tehran. IRAN. </p>