<div dir="ltr"><div dir="ltr">Dear Kevin<div>Hi again</div><div>Sorry i bothered you many times for asking many questions, but i am really confused about my case</div><div>After using the paw-pp that you suggested, i obtained E(f-vacancy)~2.0 eV!! here in attachment, i gave my input. i used a 128-atom supercell for a perfect crystal (E=32199.76950786 Ry) and by removing a Fe atom, i had E=31870.37251766 Ry and isolated Fe E=329.18792433 Ry. It seems that i cannot obtain a better result than 2 eV. I appreciated if you give me any suggestions. the paper that you mentioned talked about "Ecorr". what exactly is this parameter?</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jun 1, 2019 at 11:27 PM Kevin May <<a href="mailto:kmay@mit.edu" target="_blank">kmay@mit.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div dir="ltr" class="gmail_attr"><span class="gmail_default" style="font-family:arial,helvetica,sans-serif">Hi Mehrdad,</span></div><div dir="ltr" class="gmail_attr"><span class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></span></div><div dir="ltr" class="gmail_attr"><span class="gmail_default" style="font-family:arial,helvetica,sans-serif"></span>On Sat, Jun 1, 2019 at 1:58 PM mehrdad zamzamian <<a href="mailto:mehrdad.zamzamian@gmail.com" target="_blank">mehrdad.zamzamian@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">With regard<div>Dear Kevin</div><div>1) Actually, my reference is computational material science 44(2008) 690-694 (although they reported delta(E)= -0.456 eV/atom that i don't know what it is. because i said that it must be 1.6 eV according to <a href="https://doi.org/10.1080/09506608.2018.1560984" target="_blank">https://doi.org/10.1080/09506608.2018.1560984</a>). i also used spin-polarized (0.2 for Fe) but the same result was obtained (not better than 2.2 eV). I also used constant volume (by using relax calculation not vc_relax), but i had the same results.</div></div></blockquote><div><br></div><div style="font-family:arial,helvetica,sans-serif">The first paper you mention does not calculate vacancy formation energy, the -0.456 eV/atom is the bulk formation energy of orthorhombic Fe3C from GGA. The second paper is a review article that cites the paper by Jiang et al. that I mentioned in my last email for their ~1.6 eV formation energy of a Fe vacancy.<br></div><div style="font-family:arial,helvetica,sans-serif"></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>2) thanks for the proposed links, but in that reference they offered two non-consistent pps (C.pbe-n-kjpaw_psl.1.0.0.UPF and Fe.pbe-spn-kjpaw_psl.0.2.1.UPF) that i cannot use for Fe3C. i used these pp:</div><div>C.pbe-n-rrkjus_psl.1.0.0.UPF<br></div><div><span class="gmail_default" style="font-family:arial,helvetica,sans-serif"></span>Fe.pbe-spn-rrkjus_psl.1.0.0.UPF<br></div></div></blockquote><div><br></div><div style="font-family:arial,helvetica,sans-serif">I'm not sure what you mean by "non-consistent", or why you can't use those PAW datasets for your calculation. If you need to use ultrasoft for another reason, I'd still recommend using Fe.pbe-spn-rrkjus_psl.0.2.1.UPF instead of the psl.1.0.0 version. <br></div><div style="font-family:arial,helvetica,sans-serif"><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div></div><div>I am really confused about how can i calculate this rather simple parameter. i should mention that i calculate the Fe vacancy with:</div><div>E(Fe-vacancy)=E(perfect Fe3C)- E(Fe3C with lack of one Fe atom)- E(single Fe)</div></div></blockquote><div><br></div><div style="font-family:arial,helvetica,sans-serif">By E(single Fe) you mean 1/2 the energy of ferromagnetic BCC Fe (2 Fe per unit cell), rather than an isolated Fe atom, right? That is the chemical potential reference for Fe that was used in calculating the 1.6 eV formation energy value. And for what it's worth I don't think point defect calculations are that trivial.<br></div><div style="font-family:arial,helvetica,sans-serif"></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>3) thanks very much for the reference Rev. Mod. Phys. 86, 253 (2014).</div></div></blockquote><div><br></div><div style="font-family:arial,helvetica,sans-serif">Regarding chemical potential reference, see section II.B.2 of this article.<br></div><div style="font-family:arial,helvetica,sans-serif"><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div><div>regards</div><div><br></div><div><div style="color:rgb(80,0,80)">Seyed Mehrdad Zamzamian</div><div style="color:rgb(80,0,80)">Sharif University of Technology, Tehran, Iran</div><div style="color:rgb(80,0,80)">Energy engineering department</div><div>E-mail: <a href="mailto:mehrdad.zamzamian@gmail.com" target="_blank">mehrdad.zamzamian@gmail.com</a><span class="gmail_default" style="font-family:arial,helvetica,sans-serif"></span></div></div></div></div></blockquote><div><span class="gmail_default" style="font-family:arial,helvetica,sans-serif"></span></div><div><span class="gmail_default" style="font-family:arial,helvetica,sans-serif">Good luck,</span></div><div><span class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br>
</span></div><div><div style="font-family:arial,helvetica,sans-serif">Kevin May, PhD<br></div><div style="font-family:arial,helvetica,sans-serif">Postdoctoral Associate</div><div style="font-family:arial,helvetica,sans-serif">Department of Materials Science and Engineering</div><div style="font-family:arial,helvetica,sans-serif">Massachusetts Institute of Technology</div><span class="gmail_default" style="font-family:arial,helvetica,sans-serif"></span><br></div></div></div>
</blockquote></div>
</div>