<div dir="ltr"><div>I was curious that whether similar issues can arise also for the case of Voronoi charge analysis.....Any comments?<br><br></div>Prasenjit<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On 31 May 2013 14:41, Giuseppe Mattioli <span dir="ltr"><<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Dear Sakhrawi Taoufek (and Ari, Giovanni, Axel)<br>
I've done some test with the Bader post processing tool few months ago. I have found/can suggest:<br>
<br>
1) It works well with PAW PPs only.<br>
<br>
2) After reaching "normal" convergence in the pw.x (vc-relax, relax, scf) calculation, a new scf calculation should be done by using a denser fft grid<br>
than the one reported in the pw.x output; I do not now if nscf are implemented in this case. The convergence of Bader charges should be checked<br>
against the denser grid, but a factor 2 seems to be always sufficient to ensure convergence.<br>
<br>
3) As pointed out by Ari, the pp.x calculations must be done by using plot_num= 17, that is, the all-electron reconstructed charge density.<br>
<br>
4) Atoms on the cell boundaries seem to lead to severe instabilities of the Bader charge estimates, but I did not perform in-depth tests in this case.<br>
<br>
HTH<br>
<br>
Giuseppe<br>
<div><div class="h5"><br>
On Friday 31 May 2013 09:55:09 Giovanni La Penna wrote:<br>
> Dear Dr Sakhrawi Taoufek,<br>
><br>
>   I used the Bader's analysis you mentioned on<br>
> cube files generated with PP:<br>
><br>
> <a href="http://dx.doi.org/10.1007/s00214-011-0955-3" target="_blank">http://dx.doi.org/10.1007/s00214-011-0955-3</a><br>
><br>
> There is a problem in the core contribution<br>
> to charge density that MUST be included. I made<br>
> this manually (see article above for some details).<br>
> Maybe, PP improvements nowadays allow this<br>
> operation more easily (and more accurately).<br>
><br>
> HTH,<br>
>                Giovanni<br>
><br>
> National research council of Italy (CNR)<br>
> Institute for chemistry of organometallic compounds (ICCOM)<br>
> Sesto Fiorentino (Firenze), Italy<br>
> tel.: <a href="tel:%2B39%200555225264" value="+390555225264">+39 0555225264</a> / skype: giovannilapenna<br>
> <a href="http://www.iccom.cnr.it/lapenna" target="_blank">http://www.iccom.cnr.it/lapenna</a><br>
><br>
> On Fri, 31 May 2013, Sakhrawi Taoufek wrote:<br>
> > Dear all,<br>
> > i want to study charge transfer in interface metal/oxide with bader<br>
> > analysis, i looked in the archive of the forum,  found this link<br>
> ><br>
> > <a href="http://theory.cm.utexas.edu/vtsttools/bader/" target="_blank">http://theory.cm.utexas.edu/vtsttools/bader/</a><br>
> > but corresponds to vasp code.<br>
> > how can i proceed after generating a cube with pp.x??<br>
> > thinks<br>
> ><br>
> ><br>
> > ====================<br>
> > Sakhraoui Taoufik<br>
> > Phd student, LMCN<br>
> > Monastir, Tunisia<br>
> > ====================<br>
<br>
<br>
</div></div>********************************************************<br>
- Article premier - Les hommes naissent et demeurent<br>
libres et ègaux en droits. Les distinctions sociales<br>
ne peuvent être fondèes que sur l'utilitè commune<br>
- Article 2 - Le but de toute association politique<br>
est la conservation des droits naturels et<br>
imprescriptibles de l'homme. Ces droits sont la libertè,<br>
la propriètè, la sùretè et la rèsistance à l'oppression.<br>
********************************************************<br>
<span class="HOEnZb"><font color="#888888"><br>
   Giuseppe Mattioli<br>
   CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
   v. Salaria Km 29,300 - C.P. 10<br>
   I 00015 - Monterotondo Stazione (RM)<br>
   Tel <a href="tel:%2B%2039%2006%2090672836" value="+390690672836">+ 39 06 90672836</a> - Fax <a href="tel:%2B39%2006%2090672316" value="+390690672316">+39 06 90672316</a><br>
   E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>><br>
   ResearcherID: F-6308-2012<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>PRASENJIT GHOSH,<br>IISER Pune,<br>First floor, Central Tower, Sai Trinity Building<br>Garware Circle, Sutarwadi, Pashan<br>Pune, Maharashtra 411021, India<br>
<br>Phone: +91 (20) 2590 8203<br>Fax: +91 (20) 2589 9790
</div>