<br><br><div id="htc_header" style="">----- Reply message -----<br>From: "Ary Junior" <aryjunior@gmail.com><br>To: "PWSCF Forum" <pw_forum@pwscf.org><br>Subject: [Pw_forum] Mo pseudopotential<br>Date: Tue, Mar 26, 2013 2:38 pm<br><br></div><br>Dear all,<br><br>I was looking for a PBE USPP for molybdenum and checked that the .UPF files available in the QE website are the ones provided in pslibrary 0.2... So, based in the input file I created a test for the Mo.pbe-spn-rrkjus_psl.0.2.UPF file like this:<br>




<br> &input<br>   title='Mo',<br>   zed=42.0,<br>   rel=1,<br>   config='[Kr] 4d5.0 5s1 5p0',<br>   iswitch=2,<br>   dft='PBE',<br>   rlderiv=2.93, eminld=-3.0, emaxld=3.0, deld=0.01, nld=3<br>




 /<br> &test<br>   file_pseudo='Mo.pbe-spn-rrkjus_psl.0.2.UPF',<br>   nconf=1, <br>   configts(1)='4d5.0 5s1 5p0',<br>   ecutmin=50, ecutmax=200, decut=50<br> /<br><br>Besides some warnings after the "All-electron run" section, the difference between the AE and the PS eigenvalues are zero and it is ok...<br>




<br>     n l     nl             e AE (Ry)        e PS (Ry)    De AE-PS (Ry)<br>     3 2     4D   1( 5.00)       -<a href="tel:0.27582" value="+39027582" target="_blank">0.27582</a>       -<a href="tel:0.27582" value="+39027582" target="_blank">0.27582</a>        0.00000<br>



     1 0     5S   1( 1.00)       -<a href="tel:0.30033" value="+39030033" target="_blank">0.30033</a>       -<a href="tel:0.30033" value="+39030033" target="_blank">0.30033</a>        0.00000<br>
     2 1     5P   1( 0.00)       -0.07669       -0.07669        0.00000<br><br>However, when compare these eigenvalues with an atomic calculation using the spherical Bessel functions I got a weird result, even with a high cutoff energy:<br>


<br>     Cutoff (Ry) :  200.0<br>                           N = 1       N = 2       N = 3<br>     E(L=0) =       -23.1157 Ry  -13.3019 Ry   -9.0812 Ry<br>

     E(L=1) =       -20.7569 Ry  -11.9179 Ry   -7.7609 Ry<br>     E(L=2) =       -17.0812 Ry   -9.9089 Ry   -6.4930 Ry<br><br>Moreover, the AE vs PS logarithmic derivatives are <span lang="en"><span>quite</span> <span>discrepant, even at lower r values</span></span>... The point is that after analyzing the shape of the AE orbitals I couldn't understand the so low rcut values in the pslibrary input file in both versions, 0.2 and 0.3... I have done the same test for other PBE USPPs available in the QE website like Ni.pbe-n-rrkjus_psl.0.1.UPF, O.pbe-n-rrkjus_psl.0.1.UPF and S.pbe-n-rrkjus_psl.0.1.UPF and they are fine... So, I tried to reset the input file with radii more compatible with the AE wavefunctions in my conception, besides change the core corrections like this:<br>


<br> &input<br>   title='Mo',<br>   prefix='Mo',<br>

   zed=42.0,<br>   config='[Kr] 5s1 4d5 5p0',<br>   iswitch=3,<br>   dft='PBE',<br>   rlderiv=2.93, eminld=-3.0, emaxld=3.0, deld=0.01, nld=3,<br> /<br> &inputp<br>   pseudotype=3,<br>   file_pseudopw='Mo.PBE-spn-rrkjus_aryjr.UPF',<br>




   author='ADC',<br>   lloc=-1,<br>   rcloc=3.2,<br>   which_augfun='PSQ',<br>   rmatch_augfun=1.3,<br>   nlcc=.true.,<br>   new_core_ps=.true.,<br>   rcore=1.0,<br>   tm=.true.<br> /<br>6<br>5S  1  0  1.00  0.00  2.70  3.30  0.0<br>




5S  1  0  0.00  7.00  2.50  3.30  0.0<br>5P  2  1  0.00  0.00  3.20  3.30  0.0<br>5P  2  1  0.00  7.00  3.20  3.30  0.0<br>4D  3  2  5.00  0.00  1.50  1.70  0.0<br>4D  3  2  0.00 -0.30  0.70  1.70  0.0<br><br>Then, I think that everything is fine with the  AE and the PS eigenvalues:<br>




<br>     n l     nl             e AE (Ry)        e PS (Ry)    De AE-PS (Ry)<br>     1 0     5S   1( 1.00)       -<a href="tel:0.30033" value="+39030033" target="_blank">0.30033</a>       -<a href="tel:0.30033" value="+39030033" target="_blank">0.30033</a>        0.00000<br>



     3 2     4D   1( 5.00)       -<a href="tel:0.27582" value="+39027582" target="_blank">0.27582</a>       -<a href="tel:0.27582" value="+39027582" target="_blank">0.27582</a>        0.00000<br>
     2 1     5P   1( 0.00)       -0.07669       -0.07669       -0.00000<br><br>And after compare these eigenvalues with an atomic calculation using the spherical Bessel functions again, things look better now:<br><br>     Cutoff (Ry) :   50.0<br>

                           N = 1       N = 2       N = 3<br>     E(L=0) =        -0.3002 Ry   -0.0021 Ry    0.0380 Ry<br>


     E(L=1) =        -0.0766 Ry    0.0213 Ry    0.0681 Ry<br>     E(L=2) =        -0.2681 Ry    0.0313 Ry    0.0692 Ry<br><br>     Cutoff (Ry) :  200.0<br>                           N = 1       N = 2       N = 3<br>     E(L=0) =        -0.3002 Ry   -0.0021 Ry    0.0380 Ry<br>




     E(L=1) =        -0.0766 Ry    0.0213 Ry    0.0681 Ry<br>     E(L=2) =        -0.2683 Ry    0.0313 Ry    0.0692 Ry<br><br>The AE vs PS logarithmic derivatives are being well reproduced in energies lower than 1 Ry.... I also tested chemical states like the ionized Mo4+ and Mo6+, and the results are listed below:<br>

<br>     1 0     5S   1( 1.00)       -<a href="tel:0.30033" value="+39030033" target="_blank">0.30033</a>       -<a href="tel:0.30033" value="+39030033" target="_blank">0.30033</a>        0.00000<br>

     3 2     4D   1( 5.00)       -<a href="tel:0.27582" value="+39027582" target="_blank">0.27582</a>       -<a href="tel:0.27582" value="+39027582" target="_blank">0.27582</a>        0.00000<br>     2 1     5P   1( 0.00)       -0.07669       -0.07669       -0.00000<br>


     Etot =   -8097.243530 Ry,   -4048.621765 Ha, -110168.600949 eV<br>
     Etotps =   -32.753518 Ry,     -16.376759 Ha,    -445.634282 eV<br><br>     1 0     5S   1( 0.00)       -2.69065       -2.78173        0.09108<br>     3 2     4D   1( 2.00)       -3.40343       -3.51010        <a href="tel:0.10667" value="+39010667" target="_blank">0.10667</a><br>



     2 1     5P   1( 0.00)       -2.17762       -2.27343        0.09581<br>     dEtot_ae =       6.511981 Ry<br>     dEtot_ps =       6.625031 Ry,   Delta E=      -<a href="tel:0.113050" value="+390113050" target="_blank">0.113050</a> Ry<br>


<br>     1 0     5S   1( 0.00)       -4.30796       -4.58873        <a href="tel:0.28077" value="+39028077" target="_blank">0.28077</a><br>
     3 2     4D   1( 0.00)       -5.60739       -5.87530        <a href="tel:0.26791" value="+39026791" target="_blank">0.26791</a><br>     2 1     5P   1( 0.00)       -3.66423       -3.99180        0.32756<br>     dEtot_ae =      15.474602 Ry<br>


     dEtot_ps =      15.954369 Ry,   Delta E=      -0.479767 Ry<br>
<br>My question: <span lang="en"><span>is my</span> <span>procedure described</span> <span>above</span> <span>wrong</span><span>?</span></span> Study and understand the generation of PSPs is not easy for me and I'm trying to learn it grounded in some textbooks as well as some documents spread on the internet and of course the input files <span lang="en"><span>provided in the pslibrary project. But this PBE USPP for </span></span>molybdenum made me a little bit confused. I tested both in a MoS2 unit cell with PBE-D and the cell volume as well as the interatomic distances are fine if compared with the experiments, with differences lower than 3%... One more question: how can I improve the result for the Mo4+ and Mo6+ ionized states?<br>

<br>Thank you very much!<br><br>Ary Junior<br><br>Doctoral Student<br>
Chemistry Department<br>Universidade Federal de Juiz de Fora - MG - Brasil<br clear="all"><br>-- <br><a href="http://lattes.cnpq.br/8221674673413336" target="_blank">http://lattes.cnpq.br/8221674673413336</a>