[Wien] Formation energy of impurity charge states
Marcelo Barbosa
marcelo.b.barbosa at gmail.com
Thu Apr 11 12:34:30 CEST 2019
Thank you very much for your help!
Best regards,
Marcelo
> On 11 Apr 2019, at 16:16, SM Alay-e-Abbas <alayabbas at gmail.com> wrote:
>
> --> Since I need to calculate the electrostatic potential, is following method the correct approach? Is it also correct for spin polarized calculations?
> 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).
>
> --> Create the case.in5 file with the correct plane and number of points to be calculated
> and set ny = 1!
>
> --> Use “x lapw5 -d” to create the file lapw5.def and then change the unit 9 from “case.clmval” to “case.vtotal”
> for electrostatic potentials, it should be case.vcoul.
>
> --> Run lapw5 and get the values in case.rho
> The results you should be looking for are in case.rho_onedim
>
> On Thu, Apr 11, 2019 at 12:09 PM Marcelo Barbosa <marcelo.b.barbosa at gmail.com <mailto:marcelo.b.barbosa at gmail.com>> wrote:
> Dear Alay,
>
> Thank you very much for your help!
>
> 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.
> I just not entirely sure that I correctly understood the rest of the procedure by reading the user’s guide.
>
> Since I need to calculate the electrostatic potential, is following method the correct approach? Is it also correct for spin polarized calculations?
>
> - Create the case.in5 file with the correct plane and number of points to be calculated
> - Set “iuntits” to ATU in case.in5 to get the values in Ry
> - Set “cnorm” to VAL in case.in5
> - Use “x lapw5 -d” to create the file lapw5.def and then change the unit 9 from “case.clmval” to “case.vtotal”
> - Run lapw5 and get the values in case.rho
>
> Best regards,
> Marcelo
>
>
>> On 4 Apr 2019, at 21:30, SM Alay-e-Abbas <alayabbas at gmail.com <mailto:alayabbas at gmail.com>> wrote:
>>
>> Hello Marcelo,
>>
>> 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.
>>
>> Best Regards,
>>
>> Alay
>>
>> On Thu, Apr 4, 2019 at 5:40 PM Marcelo Barbosa <marcelo.b.barbosa at gmail.com <mailto:marcelo.b.barbosa at gmail.com>> wrote:
>> Dear Sirs,
>>
>> 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).
>>
>> The electrostatic potential is calculated in LAPW0, but how can I find the value corresponding to a specific position in my supercell?
>>
>>
>> Best regards,
>> Marcelo
>> _______________________________________________
>> Wien mailing list
>> Wien at zeus.theochem.tuwien.ac.at <mailto:Wien at 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>
>> _______________________________________________
>> Wien mailing list
>> Wien at zeus.theochem.tuwien.ac.at <mailto:Wien at 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>
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at <mailto:Wien at 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>
> _______________________________________________
> Wien mailing list
> Wien at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20190411/007ec5c8/attachment.html>
More information about the Wien
mailing list