<div dir="ltr">Hi Lorenzo,<div><br></div><div>the two tutorial should be quite equivalent, although the one Federico found is more recent and probably I updated some details. </div><div><br></div><div>Regarding your questions: as far as I understand you are only interested in U on one specific atomic species of your materials. You could do that by doing two perturbations: one on the atom you are interested in, the second on all the other atoms (which you treat as a "macro-atom") at the same time. Of course the occupation of this second "extended atom" is the sum of the occupations of the atoms it is composed by. you then construct a 2x2 response matrix that, once inverted, should give you the U you are looking for (throwing away the one for the "macro-atom"). </div><div>This only limits the number of perturbation but still requires a supercell of sufficiently large size.</div><div><br></div><div>We are working on a way to make the whole procedure more efficient, through automatizing the linear-response, but the code is not ready yet.</div><div><br></div><div>Best,<br></div><div><br></div><div>Matteo</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 2, 2016 at 5:55 PM, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:lorenzo.paulatto@impmc.upmc.fr" target="_blank">lorenzo.paulatto@impmc.upmc.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
I have a question about how to compute the abinitio U, I'm following the<br>
online Santa Barbara 2009 tutorial (and trying to use the script therein) and<br>
PRB 71, 35105 (2005).<br>
<br>
I have a system that is not very large but quite heavy to compute, and have<br>
several atomic species (3 at first, then 4) but only one is problematic, as it<br>
is could be Cerium or other Lanthanoids.<br>
<br>
It would be quite anoying and very cpu time consuming for me to properly*<br>
compute U by finite response for all the species, but I'm not sure that<br>
1. it is possible with the available toolchain<br>
2. it makes sense<br>
<br>
Also, if somebody (e.g. Matteo) has some more recent and/or more user friendly<br>
code to compute dalpha/dn I'm a taker.<br>
<br>
thank you for your help!<br>
<br>
*) by properly I mean take care that it is converged with the supercell size<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Dr. Lorenzo Paulatto<br>
IdR @ IMPMC -- CNRS & Université Paris 6<br>
<a href="tel:%2B33%20%280%291%2044%20275%20084" value="+33144275084">+33 (0)1 44 275 084</a> / skype: paulatz<br>
<a href="http://www.impmc.upmc.fr/~paulatto/" rel="noreferrer" target="_blank">http://www.impmc.upmc.fr/~<wbr>paulatto/</a><br>
23-24/4é16 Boîte courrier 115,<br>
4 place Jussieu 75252 Paris Cédex 05<br>
<br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</font></span></blockquote></div><br></div>