Hi all,<br><br>In the unfinshed output-file, I find that the energies of 4F and 5D orbitals are always <br>zero no matter which electronic configuration / matching radii / pseudization energies <br>I choose.<br><br>&nbsp;&nbsp;&nbsp;&nbsp; 4 3 
2.5 4F 1( 6.00)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000<br>&nbsp;&nbsp;&nbsp;&nbsp; 4 3 3.5 4F 1( 3.00)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000<br>&nbsp;&nbsp;&nbsp;&nbsp; 6 0 0.5 6S 1( 2.00)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.2621&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.1310&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3.5659<br>&nbsp;&nbsp;&nbsp;&nbsp; 5 2 1.5
 5D 1( 0.50)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000<br>&nbsp;&nbsp;&nbsp;&nbsp; 5 2 2.5 5D 1( 0.00)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000<br><br>Yours,<br>Ma<br><br><div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>Hi all,<br><br>When I use ld1.x to generate pseudopotential of Tb, I meet some problems<br>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<br>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<br>-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<br>
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<br>-3<br>...........<br>...........<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Wfc&nbsp;&nbsp; 4F&nbsp;&nbsp;rcut= 1.400&nbsp;&nbsp;Estimated cut-off energy=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;156.83 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;This function has&nbsp;&nbsp;&nbsp;&nbsp;0 nodes for 0 &lt; r &lt;&nbsp;&nbsp;&nbsp;&nbsp;
1.385<br>&nbsp;&nbsp;&nbsp;&nbsp; ns=&nbsp;&nbsp; 4 l=&nbsp;&nbsp; 3 sum=&nbsp;&nbsp;&nbsp;&nbsp;0.003498123 r(ikk)&nbsp;&nbsp;&nbsp;&nbsp; 2.754160645<br>&nbsp;&nbsp;&nbsp;&nbsp; from gener_pseudo&nbsp;&nbsp;: info #&nbsp;&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> &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> &amp;<br> &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> &amp;<br>8<br>5S&nbsp;&nbsp;1&nbsp;&nbsp;0&nbsp;&nbsp;2.00&nbsp;&nbsp;0.00&nbsp;&nbsp;2.20&nbsp;&nbsp;2.20 0.50<br>5P&nbsp;&nbsp;2&nbsp;&nbsp;1&nbsp;&nbsp;2.00&nbsp;&nbsp;0.00&nbsp;&nbsp;2.20&nbsp;&nbsp;2.20 0.50<br>5P&nbsp;&nbsp;2&nbsp;&nbsp;1&nbsp;&nbsp;3.50&nbsp;&nbsp;0.00&nbsp;&nbsp;2.20&nbsp;&nbsp;2.20 1.50<br>4F&nbsp;&nbsp;4&nbsp;&nbsp;3&nbsp;&nbsp;6.00&nbsp;&nbsp;0.00&nbsp;&nbsp;1.40&nbsp;&nbsp;1.40 2.50<br>4F&nbsp;&nbsp;4&nbsp;&nbsp;3&nbsp;&nbsp;
3.00&nbsp;&nbsp;0.00&nbsp;&nbsp;1.40&nbsp;&nbsp;1.40 3.50<br>5D&nbsp;&nbsp;3&nbsp;&nbsp;2&nbsp;&nbsp;0.50&nbsp;&nbsp;0.00&nbsp;&nbsp;2.10&nbsp;&nbsp;2.10 1.50<br>5D&nbsp;&nbsp;3&nbsp;&nbsp;2&nbsp;&nbsp;0.00&nbsp;&nbsp;0.00&nbsp;&nbsp;2.10&nbsp;&nbsp;2.10 2.50<br>6S&nbsp;&nbsp;1&nbsp;&nbsp;0&nbsp;&nbsp;2.00&nbsp;&nbsp;0.00&nbsp;&nbsp;2.60&nbsp;&nbsp;2.60 0.50<br><br>**************************************************************************************
<br><br>Yours,<br>Ma<br><br><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><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20080107/b4389d7f/attachment-0001.htm">
http://www.democritos.it/pipermail/pw_forum/attachments/20080107/b4389d7f/attachment-0001.htm</a><br><br><br></blockquote></div>