Hi Leonardo,<br><br>Thanks a lot for the detailed description of LDA-1/2 in the atomic code. I appreciate it.<br>Indeed, with CUT value of 3.67, I got the correct band-gap of Si. I went through one iteration of the paper,<br>
and I understand that the parameter CUT has to chosen such that it optimizes the excitation energy. However,<br>It seems that there is no general rule for specifying the fraction of electronic charge to be removed from the valence electronic configuration. For group IV element it is 0.25 (C is), on the other hand for II-VI or III-V compound it is 0.5 for anion p or cation d. What about the transition metal with incomplete d orbital (for example Ti, Zr, V, Cr, Mn)? The only transition metal which was discussed in the paper<br>
is Zn, which has completely filled d orbital. Also, what about the fractional charge to be removed from the valence of f-block elements?<br><br>Thank you, and best regards,<br><br>Rajan<br><br><div class="gmail_quote">On Thu, Mar 31, 2011 at 1:38 AM, Leonardo Matheus <span dir="ltr"><<a href="mailto:leommj@usp.br">leommj@usp.br</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi,<br>
<br>
I'm Leonardo M M Jorge, the one who put the LDA-1/2 method on QE. And<br>
yes, it is working!<br>
<br>
You indeed have to create a new pseudopotential with ld1.x and the<br>
flag "iswitch=4" (among other options, shown in the example file<br>
atomic_doc/pseudo-LDA-0.5/<a href="http://si.lda.0.5.in" target="_blank">si.lda.0.5.in</a> ), and then use it on a<br>
standard LDA calculation.<br>
<br>
You obtained results very similar to standard LDA because one of these<br>
other options is not well set. When generating the pseudo, besides the<br>
electronic configuration used (of a half-ionized atom, hence the name<br>
-0.5 ), there is also another parameter needed: rcutv, which is the<br>
radius where you cut the perturbation. The value set on the example<br>
file (1.0) is too small; for Silicon these values should be somewhere<br>
around 3.6.<br>
<br>
To obtain the right value for this rcutv parameter you need to do<br>
several calculations, changing it by small steps, and then get the<br>
rcutv where the gap of the system is an extremum (maximum or minimum).<br>
The explanation why this works is written on the reference you just<br>
cited ( L. G. Ferreira, M. Marques and L. K. Teles, Phys. Rev. B 78<br>
125116 (2008) ), but what's important is that it IS NOT a free<br>
parameter (one you keep changing until you get the result you want),<br>
as there is a procedure to obtain it.<br>
<br>
Also, as small variations of the rcutv parameters should give very<br>
small variations of the Electronic Gap (as it is an extremum for<br>
variations), we assume that you can obtain the right rcutv value for a<br>
small and simple system, and then use it on a bigger and more complex<br>
system (for a given atom; different atom kinds have different rcutv!).<br>
<br>
I hope I could explain the method for you, if you have any more<br>
questions feel free to write me.<br>
<br>
Leonardo Jorge<br>
PhD Student<br>
Universidade de São Paulo - Brazil<br>
<div><div></div><div class="h5">_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br>