<div dir="ltr"><p style="margin-bottom:0in">Dear Wien2k users</p>
<p style="margin-bottom:0in">I have successfully optimized the
volume and various ratios of  lattice parameters of an orthorhombic
system (Bi2Fe4O9) without any problem (0.33% error in a, 0.17% error in b
 and 2.86% error in c lattice parameter from the corresponding Exp 
values). Now when I try to relax the atomic
positions of the system using MSR1a, it is, unexpectedly, very slow 
(around 300 iterations) in
charge convergence. The structure for this calculation  has optimized
lattice parameters  and ratios with experimental atomic positions that
I want to relax. I have used following command 
</p>
<p style="margin-bottom:0in">runsp_lapw -p -i 999 -ec 0.0001 -cc
0.001 -fc 0.50</p>
<p style="margin-bottom:0in">Just to check I started a new simple
scf calculation which also ran into trouble due to the slow
convergence of charge density. Interestingly all the calculations
with experimental parameters and Optimizer generated structures for
option 1, 3 and 4 converge nicely (around 80 iterations with the above 
mentioned command). I really don't know what is wrong
with the calculation. Any help in this regard would be highly
appreciated. I have tried to attach scf files but they do not go through
 due to large size. I have attached first iteration and log file
 for a quick examination so that, if possible, someone can spot a 
problem.  If you need any
further information I would be happy to provide it. 
</p>
<p style="margin-bottom:0in">I am using WIEN2k_13.1 on a 12 core
machine with opensuse 11.2. 
</p>
<p style="margin-bottom:0in">Thanks in advance,</p><br><div class="gmail_signature">Dr. Shafqat Hussain Shah<br></div>
</div>