Dear Gabriele Sclauzero,<div><br></div><div>Thank you so much for your help. I did an atomic calculation according to your suggestion using that pseudo-potential. Sadly, the occupation number is not correct.</div><div>first step </div>
<div>atom    1   Tr[ns(na)]=   7.0000000 </div><div>and after optimization</div><div>atom    1   Tr[ns(na)]=   0.0748899</div><div><br></div><div>The pseudo-potential I got from website of QE named <span class="Apple-style-span" style="font-family: Verdana, Arial, Helvetica, sans-serif; font-size: 12px; color: rgb(92, 106, 92); line-height: 15px; "><a href="http://www.pwscf.org/pseudo/1.3/UPF/Co.pbe-sp-mt_gipaw.UPF" class="aranciolink" style="color: rgb(255, 153, 0); font-size: 12px; line-height: 15px; font-weight: bold; text-decoration: underline; font-family: Verdana, Arial, Helvetica, sans-serif; ">Co.pbe-sp-mt_gipaw.UPF</a></span></div>
<div><span class="Apple-style-span" style="font-family: Verdana, Arial, Helvetica, sans-serif; font-size: 12px; color: rgb(92, 106, 92); line-height: 15px; "><a href="http://www.pwscf.org/pseudo/1.3/UPF/Co.pbe-sp-mt_gipaw.UPF" class="aranciolink" style="color: rgb(255, 153, 0); font-size: 12px; line-height: 15px; font-weight: bold; text-decoration: underline; font-family: Verdana, Arial, Helvetica, sans-serif; "></a></span>and header of this pseudo-potential is copied below</div>
<span class="Apple-style-span" style="font-family: 'Times New Roman'; font-size: medium; "><pre style="word-wrap: break-word; white-space: pre-wrap; "><PP_HEADER>
   0                   Version Number
  Co                   Element
   NC                  Norm - Conserving pseudopotential
    F                  Nonlinear Core Correction
 SLA  PW   PBX  PBC    PBE  Exchange-Correlation functional
   17.00000000000      Z valence
 -287.58146877626      Total energy
      0.000      0.000 Suggested cutoff for wfc and rho
    2                  Max angular momentum component
 1193                  Number of points in mesh
    6    2             Number of Wavefunctions, Number of Projectors
 Wavefunctions         nl  l   occ
                       3S  0  2.00
                       4S  0  0.00
                       3P  1  6.00
                       4P  1  0.00
                       3D  2  7.00
                       4D  2  0.00
</PP_HEADER>
</pre></span><div>According to your post, I don't think it can be any problem, but atomic calculation I did is not correct. Can you give me further help about what can go wrong? I looked inside offset_atom_wfc.f90, but I have no idea how to check if pw picks right atomic wavefunction.... Can you give me some hints about that? Thank you so much!</div>
<div><br></div><div>Bests</div><div>-- <br>Jia Chen<br><br>Dept of Chemistry<br>Princeton University<br>Princeton, NJ 08544<br>
</div>