Hi all,<br><br>When I use ld1.x to generate pseudopotential of Tb, I meet some problems below. <br>Can someone help me? Thank you very much.<br><br>*********************************************************************************
<br>classical turning point too close to mesh&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3<br>classical turning point too close to mesh&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -4<br>classical turning point too close to mesh&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2
<br>classical turning point too close to mesh&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3<br>...........<br>...........<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Wfc&nbsp;&nbsp; 4F&nbsp; rcut= 1.400&nbsp; Estimated cut-off energy=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 156.83 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; This function has&nbsp;&nbsp;&nbsp; 0 nodes for 0 &lt; r &lt;&nbsp;&nbsp;&nbsp; 
1.385<br>&nbsp;&nbsp;&nbsp;&nbsp; ns=&nbsp;&nbsp; 4 l=&nbsp;&nbsp; 3 sum=&nbsp;&nbsp;&nbsp; 0.003498123 r(ikk)&nbsp;&nbsp;&nbsp;&nbsp; 2.754160645<br>&nbsp;&nbsp;&nbsp;&nbsp; from gener_pseudo&nbsp; : info #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1<br>&nbsp;&nbsp;&nbsp;&nbsp; chi too large beyond r_c<br><br>*********************************************************************************
<br><br>The input file is:<br><br>&nbsp;&amp;input<br>&nbsp;&nbsp; title=&#39;Tb&#39;,<br>&nbsp;&nbsp; zed=65.0,<br>&nbsp;&nbsp; rel=2,<br>&nbsp;&nbsp; rmax = 100,<br>&nbsp;&nbsp; iswitch =3,<br>&nbsp;&nbsp; rlderiv=3.01,<br>&nbsp;&nbsp; eminld=-6.0,<br>&nbsp;&nbsp; emaxld=4.0,<br>&nbsp;&nbsp; deld=0.02,<br>&nbsp;&nbsp; nld=5,
<br>&nbsp;&nbsp; config=&#39;[Kr] 4d10 5s2 5p5.5 4f9 6s2 5d0.5&#39;,<br>&nbsp;&nbsp; dft=&#39;PBE&#39;<br>&nbsp;&amp;<br>&nbsp;&amp;inputp<br>&nbsp;&nbsp; pseudotype=3,<br>&nbsp;&nbsp; lloc=0,<br>&nbsp;&nbsp; file_pseudopw=&#39;Tb.pbe&#39;,<br>&nbsp;&nbsp; nlcc =.true.,<br>&nbsp;&nbsp; rho0=0.01,<br>
!&nbsp;&nbsp; rcore= 1.1,<br>!&nbsp;&nbsp; rcloc= 2.1<br>&nbsp;&amp;<br>8<br>5S&nbsp; 1&nbsp; 0&nbsp; 2.00&nbsp; 0.00&nbsp; 2.20&nbsp; 2.20 0.50<br>5P&nbsp; 2&nbsp; 1&nbsp; 2.00&nbsp; 0.00&nbsp; 2.20&nbsp; 2.20 0.50<br>5P&nbsp; 2&nbsp; 1&nbsp; 3.50&nbsp; 0.00&nbsp; 2.20&nbsp; 2.20 1.50<br>4F&nbsp; 4&nbsp; 3&nbsp; 6.00&nbsp; 0.00&nbsp; 1.40&nbsp; 1.40 2.50<br>4F&nbsp; 4&nbsp; 3&nbsp; 
3.00&nbsp; 0.00&nbsp; 1.40&nbsp; 1.40 3.50<br>5D&nbsp; 3&nbsp; 2&nbsp; 0.50&nbsp; 0.00&nbsp; 2.10&nbsp; 2.10 1.50<br>5D&nbsp; 3&nbsp; 2&nbsp; 0.00&nbsp; 0.00&nbsp; 2.10&nbsp; 2.10 2.50<br>6S&nbsp; 1&nbsp; 0&nbsp; 2.00&nbsp; 0.00&nbsp; 2.60&nbsp; 2.60 0.50<br><br>**************************************************************************************
<br><br>Yours,<br>Ma<br><br clear="all"><br>-- <br>========================================<br>Institute of Theoretical Physics, CAS<br>Address: No.55, Zhong-Guan-Cun East Road<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Beijing, China<br>Phone:&nbsp;&nbsp;&nbsp;&nbsp;86-10-62582396
<br>Email:&nbsp;&nbsp;&nbsp;&nbsp; <a href="mailto:fengjie.ma@gmail.com">fengjie.ma@gmail.com</a>