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.

On the other hand: I do not know if the VASP calculation is correct either ...



On 03/10/2014 08:06 AM, shamik chakrabarti wrote:
Dear Prof. Blaha,

*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.*

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.

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 & we
actually will get an increment by this method?.....

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).

Looking forward to your further discussion which may shed some light
over our problems.

with regards,


On Mon, Mar 10, 2014 at 12:16 PM, Peter Blaha
<pbl...@theochem.tuwien.ac.at <mailto:pbl...@theochem.tuwien.ac.at>> wrote:

    The options "optimize vol, c/a or b/a" are "special" options for
    certain properties, but for sure a sequential optimization along
    these lines does NOT give the minimum a,b,c.

    In your case it seems you should use option [6], vary a,b,c.

    Furthermore: it is vital, that you optimize internal coordinates
    (MSR1a or min_lapw) for each of these steps (modify the resulting
    script).

    When you are talking about Fe2+,3+; you have to think whether GGA+U
    is necessary ??


    On 03/10/2014 06:32 AM, shamik chakrabarti wrote:

        Dear wien2k users,

               We have simulated structural, elcetronic properties and
        cathode
        characteristics (in terms of voltage) of a Li based material.
        When we
        extract some Li ions from the parent compounds and optimize its
        volume
        and lattice parameters, we have obtained a decrements in volume.
        This
        decrements has been followed from the concern that Fe+2 ion get
        oxidized
        in Fe+3 ions in the delithiated compounds & hence there is an
        decrements
        in Fe-O bond length in comparison to the parent (lithiated) one.

        However, simulation report on the same compound, computed in
        VASP, has
        shown that there is actually an expansion of volume after Li
        extraction.
        This expansion occurs as when Li ions are taken out from the parent
        materials there is a decrements of attraction between two layers
        connected with Li ions previously. At least this is the explanation
        given by earlier researchers.

        we optimize the material and its Li extracted counterpart by
        following
        the methods:
        (1) volume optimization by keeping a:b:c constant, (2)
        optimizing c/a
        ratio and (3) optimizing b/a ratio

        we have taken lattice parameters after (3) rd operation
        (optimization)
        and calculated  DOS and voltage.

        we have also seen....that after 3rd optimization

        (1) there is an increment in b/a ratio in the delithiated part
        than in
        the lithiated part in agreement with VASP
        (2) However, unit cell volume of the de-lithiated part is less
        than the
        lithiated compound due to more dominant reduction of Fe-O bond
        lengths
        in comparison to the increment of layer-layer separation due to
        Li ion
        extractions in our calculation.

        Now my queries are:

        (1) Is the reduction of Fe-O bond length after Li extraction is more
        than the decrements of layer-layer separation in FPLAPW
        calculation?...&
        that's why we have achieved decrements of volume instead of
        increments
        as obtained in VASP?

        (2) Whether this difference came out from some inherent fact of
        FPLAPW
        and PAW (VASP) approach?

        (3) Whether this is the limitation of the wien2k code itself?

        (4) Or, whether our methodology of calculation by first volume
        optimization & then further optimizing c/a and b/a ration....is
        actually
        a wrong path...& we should choose some other way for lattice
        parameter
        optimization after Li extraction.

        Sorry for such a long email. However, the answer of these
        queries are
        very necessary for our further research and I believe, this
        discussion
        will help the community as a whole.

        Thanks in advance,

        with regards,

        --
        Shamik Chakrabarti
        Senior Research Fellow
        Dept. of Physics & Meteorology
        Material Processing & Solid State Ionics Lab
        IIT Kharagpur
        Kharagpur 721302
        INDIA


        _________________________________________________
        Wien mailing list
        w...@zeus.theochem.tuwien.ac.__at
        <mailto:Wien@zeus.theochem.tuwien.ac.at>
        http://zeus.theochem.tuwien.__ac.at/mailman/listinfo/wien
        <http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien>
        SEARCH the MAILING-LIST at:
        
http://www.mail-archive.com/__wien@zeus.theochem.tuwien.ac.__at/index.html
        <http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html>


    --

                                           P.Blaha
    
------------------------------__------------------------------__--------------
    Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
    Phone: +43-1-58801-165300             FAX: +43-1-58801-165982
    Email: bl...@theochem.tuwien.ac.at
    <mailto:bl...@theochem.tuwien.ac.at>    WWW:
    http://info.tuwien.ac.at/__theochem/
    <http://info.tuwien.ac.at/theochem/>
    
------------------------------__------------------------------__--------------
    _________________________________________________
    Wien mailing list
    w...@zeus.theochem.tuwien.ac.__at
    <mailto:Wien@zeus.theochem.tuwien.ac.at>
    http://zeus.theochem.tuwien.__ac.at/mailman/listinfo/wien
    <http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien>
    SEARCH the MAILING-LIST at:
    http://www.mail-archive.com/__wien@zeus.theochem.tuwien.ac.__at/index.html
    <http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html>




--
Shamik Chakrabarti
Senior Research Fellow
Dept. of Physics & Meteorology
Material Processing & Solid State Ionics Lab
IIT Kharagpur
Kharagpur 721302
INDIA


_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html


--

                                      P.Blaha
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982
Email: bl...@theochem.tuwien.ac.at WWW: http://info.tuwien.ac.at/theochem/
--------------------------------------------------------------------------
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Reply via email to