<div dir="ltr"><div>--> Since I need to calculate the electrostatic potential, is following method the correct approach? Is it also correct for spin polarized calculations?</div><div><i>For electrostatic potentials you should be using case.vcoul. I don't think that the magnetic order matters here since there is only one output for electrostatic potential (*.vcoul). </i></div><div><br></div><div>--> Create the case.in5 file with the correct plane and number of points to be calculated</div><div><i>and s</i><i>et ny = 1!</i></div><div><i><br></i></div><div>--> Use “x lapw5 -d” to create the file lapw5.def and then change the unit 9 from “case.clmval” to “case.vtotal”<br></div><div><i>for electrostatic potentials, it should be case.vcoul</i>.  </div><div><br></div><div>--> Run lapw5 and get the values in case.rho</div><i>The results you should be looking for are in case.rho_onedim</i><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 11, 2019 at 12:09 PM Marcelo Barbosa <<a href="mailto:marcelo.b.barbosa@gmail.com">marcelo.b.barbosa@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 style="overflow-wrap: break-word;">Dear Alay,<div><br></div><div>Thank you very much for your help!</div><div><br></div><div>I am able to use your method to select the appropriate plane and the number of points in that plane to be calculated using lapw5.</div><div>I just not entirely sure that I correctly understood the rest of the procedure by reading the user’s guide.</div><div><br></div><div>Since I need to calculate the electrostatic potential, is following method the correct approach? Is it also correct for spin polarized calculations?</div><div><br></div><div>- Create the case.in5 file with the correct plane and number of points to be calculated</div><div>- Set “iuntits” to ATU in case.in5 to get the values in Ry</div><div>- Set “cnorm” to VAL in case.in5</div><div>- Use “x lapw5 -d” to create the file lapw5.def and then change the unit 9 from “case.clmval” to “case.vtotal”</div><div>- Run lapw5 and get the values in case.rho</div><div><br></div><div>Best regards,</div><div>Marcelo</div><div><br><div><div><br><blockquote type="cite"><div>On 4 Apr 2019, at 21:30, SM Alay-e-Abbas <<a href="mailto:alayabbas@gmail.com" target="_blank">alayabbas@gmail.com</a>> wrote:</div><br class="gmail-m_-2220089208851169109Apple-interchange-newline"><div><div dir="ltr">Hello Marcelo,<div><br></div><div>You may do this by selecting an appropriate atom centered plane (with reasonable width) and then setting npy = 1 in case.in5 before running lapw5. See section 8.13.3 of the userguide for more details.</div><div><br></div><div>Best Regards,</div><div><br></div><div>Alay </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 4, 2019 at 5:40 PM Marcelo Barbosa <<a href="mailto:marcelo.b.barbosa@gmail.com" target="_blank">marcelo.b.barbosa@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">Dear Sirs,<br>
<br>
When comparing the formation energies of two charge states of an impurity in a semiconductor as a function of the Fermi level, it is said that a correction term must be added to align the electrostatic potential from the supercell with the impurity (using a position far from the impurity) and from a supercell of the pure bulk material with the same size (see e.g. DOI: 10.1063/1.1682673).<br>
<br>
The electrostatic potential is calculated in LAPW0, but how can I find the value corresponding to a specific position in my supercell?<br>
<br>
<br>
Best regards,<br>
Marcelo<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.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" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>
_______________________________________________<br>Wien mailing list<br><a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.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/wien@zeus.theochem.tuwien.ac.at/index.html</a><br></div></blockquote></div><br></div></div></div>_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.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" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>