<div dir="ltr">Dear all,<br><br>I am trying to calculate the atomic charge of BaMoO3 with Bader analysis. I am using Quantum ESPRESSO for post-processing and Critic2 for Bader analyses. the exchange-correlation function is PBEsol and the type of pseudopotentials are PAW. so I used the plot_num=21 for pp.x calculation. <br>I have read rest of the Pw_forum for bader analysis but this statement has a lot of ambiguity. as I found out, for Bader analysis I have to use denser FFT (but how much?). for this purpose I have tried different nr1, nr2, nr3:<br><br>====================================<br>SCF.in:<br>nr1=600, nr2=600, nr3=600<br><br>SCF.out:<br>     Dense  grid:   307531 G-vectors     FFT dimensions: ( 600, 600, 600)<br>     Smooth grid:    59085 G-vectors     FFT dimensions: ( 576, 576, 576)<br><br>Critic2 Out:<br>* Critical point list, final report (non-equivalent cps)<br>  Topological class (n|b|r|c): 3(5) 3(42) 0(0) 2(11) <br>  Morse sum: -48<br><br>* Integrated atomic properties<br># Integrable properties 1 to 4<br># Id   cp   ncp   Name  Z   mult     Volume            Pop             Lap            Charge       <br>  1    1    1      Ba   56  --   1.37246087E+02  6.70321481E+01  7.09104580E-02  6.70321481E+01<br>  2    2    2      Mo   42  --   6.97318032E+01  4.13474024E+01  2.91573989E+00  4.13474024E+01<br>  3    3    3      O    8   --   7.70440931E+01  9.21734416E+00 -9.95550129E-01  9.21734416E+00<br>  4    4    3      O    8   --   7.70440931E+01  9.21734416E+00 -9.95550129E-01  9.21734416E+00<br>  5    5    3      O    8   --   7.70440931E+01  9.21734416E+00 -9.95550128E-01  9.21734416E+00<br>------------------------------------------------------------------------------------------------<br>  Sum                            4.38110170E+02  1.36031583E+02 -4.21448737E-08  1.36031583E+02<br>====================================<br><br><br>====================================<br>SCF.in:<br>nr1=800, nr2=800, nr3=800<br><br>SCF.out:<br>     Dense  grid:   429697 G-vectors     FFT dimensions: ( 800, 800, 800)<br>     Smooth grid:   108767 G-vectors     FFT dimensions: ( 750, 750, 750)<br><br>Critic2 Out:<br>* Critical point list, final report (non-equivalent cps)<br>  Topological class (n|b|r|c): 3(5) 0(0) 3(60) 3(17) <br>  Morse sum: 48<br><br><br>* Integrated atomic properties<br># Integrable properties 1 to 4<br># Id   cp   ncp   Name  Z   mult     Volume            Pop             Lap             Charge         <br>  1    1    1      Ba   56  --   1.37215953E+02  5.97476664E+01  2.80099625E-02  5.97476664E+01<br>  2    2    2      Mo   42  --   6.97266839E+01  4.04985130E+01  1.20777162E+00  4.04985130E+01<br>  3    3    3      O    8   --   7.70558444E+01  9.21662377E+00 -4.11927139E-01  9.21662377E+00<br>  4    4    3      O    8   --   7.70558444E+01  9.21662377E+00 -4.11927139E-01  9.21662377E+00<br>  5    5    3      O    8   --   7.70558444E+01  9.21662377E+00 -4.11927139E-01  9.21662377E+00<br>------------------------------------------------------------------------------------------------<br>  Sum                            4.38110170E+02  1.27896051E+02  1.70404130E-07  1.27896051E+02<br>====================================<br><br><br><br>as you can see, for different "nr" I found out different atomic charge. even with a Dense grid around 429697 G-vectors, the charge does not get converged (Ba = 59.7).<br><br> <br><br><br>Please if somebody has a clue about this let me know.<br><br><br><br><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span style="color:rgb(153,153,153)"><span style="font-family:verdana,sans-serif"><b>Mahdi FaghihNasiri</b></span></span><span><font color="#888888"><br style="font-family:verdana,sans-serif"></font></span><span><font color="#888888">Department of Physics<br>Shahrood University of Technology<br>Shahrood,<span style="font-family:verdana,sans-serif"> Iran</span><font color="#888888"><b><span style="font-family:verdana,sans-serif"></span></b></font></font></span></div></div></div></div></div>
</div>