<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div style="width: 1076.34px;">Dr. Lorenzo Paulatto,</div><div style="width: 1076.34px;"><br></div><div style="width: 1076.34px;">Thank you so much for the detailed information.</div><div style="width: 1076.34px;"><br></div><div style="width: 1076.34px;">However, I did not find the suitable PP in <span style="font-family: arial; white-space: pre-wrap;">GBRV library or</span><span style="font-family: arial; white-space: pre-wrap;"> the </span><span style="font-family: arial; white-space: pre-wrap;">pslibrary.</span></div><div style="width: 1076.34px;"><br></div><div style="width: 1076.34px;">Actually, I want to generate Troullier-Martins norm-conserving PBE PP with the ionized configurations (3p)6(4s)0(3d)0, (4p)6(5s)0(4d)0, and (5p)6(6s)0(5d)0 for K, Rb, and Cs, respectively. The nonlinear core correction is applied to all the three elements and the scalar relativistic effects are taken into account in Rb and Cs.</div><div style="width: 1076.34px;"><br></div><div style="width: 1076.34px;">The input file for Rb is shown below:</div><div style="width: 1076.34px;">###############################</div><div style="width: 1076.34px;"><div style="width: 1076.34px;"> &input</div><div style="width: 1076.34px;">        iswitch=3,</div><div style="width: 1076.34px;"><br></div><div style="width: 1076.34px;">        rlderiv=1.5,</div><div style="width: 1076.34px;">        eminld=-8.0,</div><div style="width: 1076.34px;">        emaxld=0.0,</div><div style="width: 1076.34px;">        deld=0.01d0,</div><div style="width: 1076.34px;">        nld=4,</div><div style="width: 1076.34px;"><br></div><div style="width: 1076.34px;">        rel=1,</div><div style="width: 1076.34px;">        zed=37.0,</div><div style="width: 1076.34px;">        config="[Ar] 3d10 4s2 4p6 5s0 4d0"</div><div style="width: 1076.34px;">        dft='PBE',</div><div style="width: 1076.34px;"> /</div><div style="width: 1076.34px;"> &inputp</div><div style="width: 1076.34px;">   tm=.true.</div><div style="width: 1076.34px;">   lloc=-1,</div><div style="width: 1076.34px;">   nlcc=.true.</div><div style="width: 1076.34px;">   rcloc=1.5</div><div style="width: 1076.34px;">   rcore=1.2,</div><div style="width: 1076.34px;">   pseudotype=1,</div><div style="width: 1076.34px;">   file_pseudopw='Rb.pbe-mt_scalar.UPF',</div><div style="width: 1076.34px;"> /</div><div style="width: 1076.34px;">3</div><div style="width: 1076.34px;">5S  5  0  0.00  0.00  1.50  1.50</div><div style="width: 1076.34px;">4P  4  1  6.00  0.00  1.50  1.50</div><div style="width: 1076.34px;">4D  4  2  0.00  0.00  1.50  1.50</div></div><div style="width: 1076.34px;">###############################</div><div style="width: 1076.34px;">I tried to change some parameters like nlcc, rcloc, etc. it doesn't work, always with the error information:</div><div style="width: 1076.34px;"><br></div><div style="width: 1076.34px;"><div style="width: 1076.34px;"> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div style="width: 1076.34px;">     Error in routine run_pseudo (1):</div><div style="width: 1076.34px;">     Errors in PS-KS equation</div><div style="width: 1076.34px;"> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div style="width: 1076.34px;"><br></div></div><div style="width: 1076.34px;">Best regards</div><div style="width: 1076.34px;">Zhishuo Huang</div></div><br><br><span title="neteasefooter"><p> </p></span>