<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:宋体;
        panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:"\@宋体";
        panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:#0563C1;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:#954F72;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal><span style='color:#1F497D'>Thank you for your reply ! <o:p></o:p></span></p><p class=MsoNormal><span style='color:#1F497D'>I tried this strategy fingers crossed, but it didn’t work, as could be expected. On a single unit cell short of 2 electrons, they mostly left from the 3 I atoms (+0.5 each compared to the “normal” case), and just +0.25 for the Sn atom. Furthermore, as I want to model it as a defect, I will eventually want to have a supercell (maybe 3x3x3), with only one of the Sn being in the unusual oxidation state. So I guess that would make that approach even less realistic, as there is no reason that the code should decide to arbitrarily deplete one of the atoms rather than the others.<o:p></o:p></span></p><p class=MsoNormal><span style='color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='color:#1F497D'>I will try to have a closer look to the discussion link you sent in your other message and see if I can apply the idea to my case.<o:p></o:p></span></p><p class=MsoNormal><span style='color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='color:#1F497D'>Julien<o:p></o:p></span></p><p class=MsoNormal><span style='color:#1F497D'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b>From:</b> Marzari Nicola [mailto:nicola.marzari@epfl.ch] <br><b>Sent:</b> mercredi 24 avril 2019 12:38<br><b>To:</b> Quantum Espresso users Forum; julien_barbaud@sjtu.edu.cn<br><b>Subject:</b> Re: [QE-users] forcing Sn(iv) or Sn(0) charges in a MASnI3 crystal<o:p></o:p></p></div></div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal><span style='font-size:12.0pt'><o:p> </o:p></span></p></div><p class=MsoNormal style='margin-bottom:12.0pt'>Sorry, different issue - you can change the number of electrons in the unit cell, and hope they go or leave from the right place.<o:p></o:p></p><div id=AppleMailSignature><p class=MsoNormal>Sent from a tiny keyboard... Contact info: <o:p></o:p></p><div><div><p class=MsoNormal><a href="http://theossrv1.epfl.ch/Main/Contact">http://theossrv1.epfl.ch/Main/Contact</a><o:p></o:p></p></div></div></div><div><p class=MsoNormal style='margin-bottom:12.0pt'><br>On 24 Apr 2019, at 00:32, Marzari Nicola <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>> wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal style='margin-bottom:12.0pt'>Here is a discussion <a href="http://theossrv1.epfl.ch/Main/OxidationStates">http://theossrv1.epfl.ch/Main/OxidationStates</a><o:p></o:p></p><div id=AppleMailSignature><p class=MsoNormal>Sent from a tiny keyboard... Contact info: <o:p></o:p></p><div><div><p class=MsoNormal><a href="http://theossrv1.epfl.ch/Main/Contact">http://theossrv1.epfl.ch/Main/Contact</a><o:p></o:p></p></div></div></div><div><p class=MsoNormal style='margin-bottom:12.0pt'><br>On 24 Apr 2019, at 00:28, JULIEN, CLAUDE, PIERRE BARBAUD <<a href="mailto:julien_barbaud@sjtu.edu.cn">julien_barbaud@sjtu.edu.cn</a>> wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><p class=MsoNormal><span lang=FR>Dear users,</span><o:p></o:p></p><p class=MsoNormal><span lang=FR> </span><o:p></o:p></p><p class=MsoNormal>I am working on a MASnI3 crystal. In this structure, the Sn can usually be considered as a Sn2+ cation. I ran some calculations on the system, and performed a Bader charge analysis on an all-electron paw charge density. It seems to confirm that the tin is in Sn(2) configuration with bader charge 48.3 instead of 50 (this is, by the way, inconsistent with the results of Lowdin analysis as implemented in projwfc.x, which gives pretty much the full 4d10 5s2 5p2 orbitals of Sn(0) ).<o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>So everything is as expected so far (from the Bader point of view at least). However, I would like to model a MASnI3 with a “defect” consisting of an “oxidized or reduced” tin atom given by Sn(iv) or Sn(0) in this material. Indeed, it was reported <i>in J. Mater. Chem. A, 2018,6, 21389-21395</i> that some degradation mechanisms can lead to the presence of such states, and I want to explore their consequence on the material properties.<o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>However, I am not sure how to tackle this. My first idea was that I probably needed to create a pseudopotential with 2 missing or additional valence electrons. But on second thought, this method might be valid if we have missing core electrons, but not for valence. I highly doubt that  it would give me the expected result once I place it in the crystal lattice, given that there is really no reason for those additional electrons to gently “stay on their starting atom”, so to speak.<o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>So is there a reliable method to study such a system by “forcing” an oxidization state for an atom in a crystal ? This task seems to be made difficult by the very subjective definition of an “atomic charge” in the framework of quantum mechanics and DFT…<o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>Thanks in advance<o:p></o:p></p><p class=MsoNormal>Julien<o:p></o:p></p></div></blockquote><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman","serif"'>_______________________________________________<br>Quantum Espresso is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso">www.max-centre.eu/quantum-espresso</a>)<br>users mailing list <a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br><a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><o:p></o:p></span></p></div></blockquote></div></blockquote><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman","serif"'>_______________________________________________<br>Quantum Espresso is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso">www.max-centre.eu/quantum-espresso</a>)<br>users mailing list <a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br><a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><o:p></o:p></span></p></div></blockquote></div></body></html>