<div dir="ltr">Sir,<div><br></div><div> The forces in our fully optimized structure is less than 3mRy/a.u. per atom...also we have done coordinate optimization by taking fully optimized lattice parameters......</div><div>
<br></div><div>Yes...we have not done the simultaneous optimization of lattice parameters and atomic coordinates...while we have done all the optimizations in steps: volume, c/a, b/a, atomic coordinates.</div><div><br></div>
<div>Also, there are several peer reviewed papers on calculation on the same material with VASP and they have all achieved volume expansion.......</div><div> </div><div>Sir...is this mean ...from this discussion..it is coming out that we should try with simultaneous optimization of lattice parameters and atomic coordinates for final verification....?</div>
<div><br></div><div>but as the energy is mainly depends on volume...& not on the lattice parameters...I really have doubt whether this method would change the trend of volume contraction towards expansion.</div></div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Mar 10, 2014 at 12:46 PM, Peter Blaha <span dir="ltr"><<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">You cannot expect that a,b,c is correct when you fix atomic positions (which are allowed to vary according to symmetry). But without doing these calculations (or at least checking how large are the forces in your "best" structure) one cannot say anything.<br>
<br>
On the other hand: I do not know if the VASP calculation is correct either ...<br>
<br>
<br>
<br>
On 03/10/2014 08:06 AM, shamik chakrabarti wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Prof. Blaha,<br>
<br>
*we have also used the option 6 for optimization of delithiated<div class=""><br>
compound...but in that case also we have obtained decrements in unit<br></div>
cell volume with Li ion extraction.*<div class=""><br>
<br>
We have also used GGA, and GGA+U for lattice parameter optimization. As,<br>
usual with GGA+U approach we have obtained higher lattice parameters<br>
than the case with GGA. However, decrements of unit cell volume with Li<br>
extraction has been apparent for both the cases.<br>
<br>
Yes it could be the fact that we may need to modify the script for<br>
simultaneous optimization of coordinates with lattice parameter<br>
optimization....whether it could lead to the solution?....whether the<br>
decrements of unit cell volume with Li delithiation can be avoided & we<br>
actually will get an increment by this method?.....<br>
<br>
Most importantly, although we are getting a volume reduction with Li<br>
de-intercalation (opposite to VASP)...we are getting very well matched<br>
values of density of states and voltage (from total energy consideration).<br>
<br>
Looking forward to your further discussion which may shed some light<br>
over our problems.<br>
<br>
with regards,<br>
<br>
<br>
On Mon, Mar 10, 2014 at 12:16 PM, Peter Blaha<br></div><div><div class="h5">
<<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a> <mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.<u></u>tuwien.ac.at</a>>> wrote:<br>
<br>
The options "optimize vol, c/a or b/a" are "special" options for<br>
certain properties, but for sure a sequential optimization along<br>
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<br>
(MSR1a or min_lapw) for each of these steps (modify the resulting<br>
script).<br>
<br>
When you are talking about Fe2+,3+; you have to think whether GGA+U<br>
is necessary ??<br>
<br>
<br>
On 03/10/2014 06:32 AM, shamik chakrabarti wrote:<br>
<br>
Dear wien2k users,<br>
<br>
We have simulated structural, elcetronic properties and<br>
cathode<br>
characteristics (in terms of voltage) of a Li based material.<br>
When we<br>
extract some Li ions from the parent compounds and optimize its<br>
volume<br>
and lattice parameters, we have obtained a decrements in volume.<br>
This<br>
decrements has been followed from the concern that Fe+2 ion get<br>
oxidized<br>
in Fe+3 ions in the delithiated compounds & hence there is an<br>
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<br>
VASP, has<br>
shown that there is actually an expansion of volume after Li<br>
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<br>
following<br>
the methods:<br>
(1) volume optimization by keeping a:b:c constant, (2)<br>
optimizing c/a<br>
ratio and (3) optimizing b/a ratio<br>
<br>
we have taken lattice parameters after (3) rd operation<br>
(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<br>
than in<br>
the lithiated part in agreement with VASP<br>
(2) However, unit cell volume of the de-lithiated part is less<br>
than the<br>
lithiated compound due to more dominant reduction of Fe-O bond<br>
lengths<br>
in comparison to the increment of layer-layer separation due to<br>
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<br>
calculation?...&<br>
that's why we have achieved decrements of volume instead of<br>
increments<br>
as obtained in VASP?<br>
<br>
(2) Whether this difference came out from some inherent fact of<br>
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 & then further optimizing c/a and b/a ration....is<br>
actually<br>
a wrong path...& we should choose some other way for lattice<br>
parameter<br>
optimization after Li extraction.<br>
<br>
Sorry for such a long email. However, the answer of these<br>
queries are<br>
very necessary for our further research and I believe, this<br>
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 & Meteorology<br>
Material Processing & Solid State Ionics Lab<br>
IIT Kharagpur<br>
Kharagpur 721302<br>
INDIA<br>
<br>
<br></div></div>
______________________________<u></u>___________________<br>
Wien mailing list<br>
Wien@zeus.theochem.tuwien.ac._<u></u>_at<br>
<mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.<u></u>tuwien.ac.at</a>><br>
<a href="http://zeus.theochem.tuwien." target="_blank">http://zeus.theochem.tuwien.</a>__<a href="http://ac.at/mailman/listinfo/wien" target="_blank"><u></u>ac.at/mailman/listinfo/wien</a><div class=""><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:<br></div>
<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>
<<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>
<br>
--<br>
<br>
P.Blaha<br>
------------------------------<u></u>__----------------------------<u></u>--__--------------<div class=""><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><br></div>
<mailto:<a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.<u></u>ac.at</a>> WWW:<br>
<a href="http://info.tuwien.ac.at/__theochem/" target="_blank">http://info.tuwien.ac.at/__<u></u>theochem/</a><br>
<<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>
Wien@zeus.theochem.tuwien.ac._<u></u>_at<br>
<mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.<u></u>tuwien.ac.at</a>><br>
<a href="http://zeus.theochem.tuwien." target="_blank">http://zeus.theochem.tuwien.</a>__<a href="http://ac.at/mailman/listinfo/wien" target="_blank"><u></u>ac.at/mailman/listinfo/wien</a><div class=""><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:<br></div>
<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><div><div class="h5"><br>
<<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>
<br>
<br>
<br>
--<br>
Shamik Chakrabarti<br>
Senior Research Fellow<br>
Dept. of Physics & Meteorology<br>
Material Processing & Solid State Ionics Lab<br>
IIT Kharagpur<br>
Kharagpur 721302<br>
INDIA<br>
<br>
<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>
<br>
</div></div></blockquote><div class="HOEnZb"><div class="h5">
<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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Shamik Chakrabarti<br>Senior Research Fellow<br>Dept. of Physics & Meteorology<br>Material Processing & Solid State Ionics Lab<br>IIT Kharagpur<br>
Kharagpur 721302<br>INDIA
</div>