<div dir="ltr">Hi Máximo,<br><br>In my experience, especially when dealing with spinel LiMn2O4, you also need to address the Jahn-Teller instability induced by Mn3+. If you slightly distort the spinel unit cell then the two different Mn sites (formally 3+ and 4+) will naturally emerge during relaxation (of course, you need the right level of theory as Nicola said): Mn3+ sites will be in a tetragonally distorted octahedral environment and Mn4+ sites will be in a more octahedral environment. If you start from a high-symmetry unit cell, you can also try to break the initial high symmetry of the charge/spin density by manually distinguishing between the two Mn sites: just name one type of Mn let's say Mn1 and the other one Mn2 and then use the same pseudopotential for both. Then, you can initialize the magnetization to a different value for the two different species.<br><div><br></div><div>HTH</div><div>Marton Voros</div><div><br></div><div>--</div><div>Materials Science Division</div><div>Argonne National Laboratory</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 17, 2017 at 8:45 AM, Nicola Marzari <span dir="ltr"><<a href="mailto:nicola.marzari@epfl.ch" target="_blank">nicola.marzari@epfl.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
Indeed! At the very least, use DFT+U, and/or explore<br>
hybrid functionals. DFT+U+V seems to be very promising,<br>
but it's not in the code there.<br>
<br>
                                nicola<br>
<span class=""><br>
<br>
<br>
On 17/04/2017 15:36, Ary Junior wrote:<br>
> Hi,<br>
><br>
> Search for "constrained DFT" and you will find some documents that would<br>
> be useful for you, like this one<br>
><br>
> <a href="http://theossrv1.epfl.ch/Main/OxidationStates" rel="noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/<wbr>OxidationStates</a><br>
><br>
> Ary Ferreira<br>
><br>
> FAPESP postdoctoral fellow<br>
> UFSCar - Brazil<br>
><br>
</span><span class="">> On Sun, Apr 16, 2017 at 1:06 PM, Máximo Ramírez <<a href="mailto:aquiles011@gmail.com">aquiles011@gmail.com</a><br>
</span><span class="">> <mailto:<a href="mailto:aquiles011@gmail.com">aquiles011@gmail.com</a>>> wrote:<br>
><br>
>     Dear,<br>
><br>
>     I'm interested in simulate metal oxide that has two different<br>
>     valence state (oxidation number) in the same compound.<br>
><br>
>     For example:<br>
><br>
>     Is well-known from experimental results that spinel compounds have<br>
>     Mn ions with valence +3 and +4, LiMn3+Mn4+O4.<br>
><br>
>     How can I simulate this behavior in QE?<br>
><br>
>     And also, How can I post-process the two atoms like were different?<br>
><br>
>     Thank in Advance!<br>
><br>
><br>
</span>>     <<a href="https://www.researchgate.net/publication/245108648_Factors_affecting_Li_mobility_in_spinel_LiMn_2O_4-A_first-principles_study_by_GGA_and_GGA_U_methods" rel="noreferrer" target="_blank">https://www.researchgate.net/<wbr>publication/245108648_Factors_<wbr>affecting_Li_mobility_in_<wbr>spinel_LiMn_2O_4-A_first-<wbr>principles_study_by_GGA_and_<wbr>GGA_U_methods</a>><br>
<span class="">>     --<br>
>     Sincerely,<br>
><br>
>     Máximo Ramírez<br>
>     Research Assistant<br>
>     Nanotechnology Laboratory, INTEC<br>
>     Santo Domingo, Dominican Republic<br>
><br>
><br>
>     ______________________________<wbr>_________________<br>
>     Pw_forum mailing list<br>
</span>>     <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a> <mailto:<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>><br>
>     <a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
<span class="">>     <<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a>><br>
><br>
><br>
><br>
><br>
> --<br>
> <a href="http://lattes.cnpq.br/8221674673413336" rel="noreferrer" target="_blank">http://lattes.cnpq.br/<wbr>8221674673413336</a><br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
><br>
<br>
--<br>
</span>------------------------------<wbr>------------------------------<wbr>----------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
Director, National Centre for Competence in Research NCCR MARVEL, EPFL<br>
<a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/<wbr>Contact</a> <a href="http://nccr-marvel.ch/en/project" rel="noreferrer" target="_blank">http://nccr-marvel.ch/en/<wbr>project</a><br>
<div class="HOEnZb"><div class="h5">______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</div></div></blockquote></div><br></div>