<div dir="ltr">Dear Prof. Blaha,<div><br></div><div>                        <b>we have also used the option 6 for optimization of delithiated compound...but in that case also we have obtained decrements in unit cell volume with Li ion extraction.</b></div>
<div><br></div><div>We have also used GGA, and GGA+U for lattice parameter optimization. As, usual with GGA+U approach we have obtained higher lattice parameters than the case with GGA. However, decrements of unit cell volume with Li extraction has been apparent for both the cases.</div>
<div><br></div><div>Yes it could be the fact that we may need to modify the script for simultaneous optimization of coordinates with lattice parameter optimization....whether it could lead to the solution?....whether the decrements of unit cell volume with Li delithiation can be avoided &amp; we actually will get an increment by this method?.....</div>
<div><br></div><div>Most importantly, although we are getting a volume reduction with Li de-intercalation (opposite to VASP)...we are getting very well matched values of density of states and voltage (from total energy consideration).</div>
<div><br></div><div>Looking forward to your further discussion which may shed some light over our problems.</div><div><br></div><div>with regards,</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Mar 10, 2014 at 12:16 PM, Peter Blaha <span dir="ltr">&lt;<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">The options &quot;optimize vol, c/a or b/a&quot; are &quot;special&quot; options for certain properties, but for sure a sequential optimization along these lines does NOT give the minimum a,b,c.<br>

<br>
In your case it seems you should use option [6], vary a,b,c.<br>
<br>
Furthermore: it is vital, that you optimize internal coordinates (MSR1a or min_lapw) for each of these steps (modify the resulting script).<br>
<br>
When you are talking about Fe2+,3+; you have to think whether GGA+U is necessary ??<div><div class="h5"><br>
<br>
On 03/10/2014 06:32 AM, shamik chakrabarti wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Dear wien2k users,<br>
<br>
      We have simulated structural, elcetronic properties and cathode<br>
characteristics (in terms of voltage) of a Li based material. When we<br>
extract some Li ions from the parent compounds and optimize its volume<br>
and lattice parameters, we have obtained a decrements in volume. This<br>
decrements has been followed from the concern that Fe+2 ion get oxidized<br>
in Fe+3 ions in the delithiated compounds &amp; hence there is an decrements<br>
in Fe-O bond length in comparison to the parent (lithiated) one.<br>
<br>
However, simulation report on the same compound, computed in VASP, has<br>
shown that there is actually an expansion of volume after Li extraction.<br>
This expansion occurs as when Li ions are taken out from the parent<br>
materials there is a decrements of attraction between two layers<br>
connected with Li ions previously. At least this is the explanation<br>
given by earlier researchers.<br>
<br>
we optimize the material and its Li extracted counterpart by following<br>
the methods:<br>
(1) volume optimization by keeping a:b:c constant, (2) optimizing c/a<br>
ratio and (3) optimizing b/a ratio<br>
<br>
we have taken lattice parameters after (3) rd operation (optimization)<br>
and calculated  DOS and voltage.<br>
<br>
we have also seen....that after 3rd optimization<br>
<br>
(1) there is an increment in b/a ratio in the delithiated part than in<br>
the lithiated part in agreement with VASP<br>
(2) However, unit cell volume of the de-lithiated part is less than the<br>
lithiated compound due to more dominant reduction of Fe-O bond lengths<br>
in comparison to the increment of layer-layer separation due to Li ion<br>
extractions in our calculation.<br>
<br>
Now my queries are:<br>
<br>
(1) Is the reduction of Fe-O bond length after Li extraction is more<br>
than the decrements of layer-layer separation in FPLAPW calculation?...&amp;<br>
that&#39;s why we have achieved decrements of volume instead of increments<br>
as obtained in VASP?<br>
<br>
(2) Whether this difference came out from some inherent fact of FPLAPW<br>
and PAW (VASP) approach?<br>
<br>
(3) Whether this is the limitation of the wien2k code itself?<br>
<br>
(4) Or, whether our methodology of calculation by first volume<br>
optimization &amp; then further optimizing c/a and b/a ration....is actually<br>
a wrong path...&amp; we should choose some other way for lattice parameter<br>
optimization after Li extraction.<br>
<br>
Sorry for such a long email. However, the answer of these queries are<br>
very necessary for our further research and I believe, this discussion<br>
will help the community as a whole.<br>
<br>
Thanks in advance,<br>
<br>
with regards,<br>
<br>
--<br>
Shamik Chakrabarti<br>
Senior Research Fellow<br>
Dept. of Physics &amp; Meteorology<br>
Material Processing &amp; Solid State Ionics Lab<br>
IIT Kharagpur<br>
Kharagpur 721302<br>
INDIA<br>
<br>
<br></div></div>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
<br>
</blockquote>
<br>
-- <br>
<br>
                                      P.Blaha<br>
------------------------------<u></u>------------------------------<u></u>--------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a>    WWW: <a href="http://info.tuwien.ac.at/theochem/" target="_blank">http://info.tuwien.ac.at/<u></u>theochem/</a><br>
------------------------------<u></u>------------------------------<u></u>--------------<br>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Shamik Chakrabarti<br>Senior Research Fellow<br>Dept. of Physics &amp; Meteorology<br>Material Processing &amp; Solid State Ionics Lab<br>IIT Kharagpur<br>Kharagpur 721302<br>
INDIA
</div>