<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div>Dear Lorenzo Paulatto</div><div><br></div><div>Thank you for your reply and sorry for write back so late.</div><div><br></div><div>Actually, in the tutorial "<span class="fontstyle0">Notes on pseudopotential generation</span>" by <span class="fontstyle0">Paolo Giannozzi, it is said that </span>"Heavy alkali metals (Rb, Cs, maybe also K) have a large polarizable core. PP¡¯s <span class="fontstyle0">with just one electron may not always give satisfactory results</span> ", so "<span class="fontstyle0">promoting the highest core states </span><span class="fontstyle2">ns </span><span class="fontstyle0">and </span><span class="fontstyle2">np</span><span class="fontstyle0">, or </span><span class="fontstyle2">nd</span><span class="fontstyle0">, into valence may </span><span class="fontstyle0">be a computationally expensive but obliged way to improve poor transferability</span> ". However, all the PBE PPs of Rb, Cs, even K are with just one electron, and none of them is ionized. This is the first reason I want to generate PP.</div><div><br></div><div>The second reason is that for K, the generated PP gives excellent phonon dispersion results compared with the usual one. </div><div><br></div><div>For the other two elements, Rb, Cs, the usual PBE PP could not give good results about phonon calculation, and it seems that I need to generate ionized PPs.</div><div><br></div><div>However, I got the error like below when generating ionized PBE PP for Rb and Cs, while everything is ok for the ionized PBE PP generation of K.</div><div><pre style="font-family: courier, "courier new", monospace; margin-top: 0px; margin-bottom: 0px;">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
      Error in routine run_pseudo (1):
      Errors in PS-KS equation
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</pre></div><div><pre style="font-family: courier, "courier new", monospace; margin-top: 0px; margin-bottom: 0px;">I tried using different channel, but it did not work.</pre></div><div><br></div><div>Any suggestion to solved this problem is welcome. </div><div><br></div><div>Best regards</div><div>Zhishuo Huang</div><div><br></div></div><br><br><span title="neteasefooter"><p> </p></span>