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>     4 3 
2.5 4F 1( 6.00)         0.0000         0.0000         0.0000<br>     4 3 3.5 4F 1( 3.00)         0.0000         0.0000         0.0000<br>     6 0 0.5 6S 1( 2.00)        -0.2621        -0.1310        -3.5659<br>     5 2 1.5
 5D 1( 0.50)         0.0000         0.0000         0.0000<br>     5 2 2.5 5D 1( 0.00)         0.0000         0.0000         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           4           3<br>3<br>classical turning point too close to mesh           4           3<br>-4<br>classical turning point too close to mesh           5           2<br>
2<br>classical turning point too close to mesh           5           2<br>-3<br>...........<br>...........<br>      Wfc   4F  rcut= 1.400  Estimated cut-off energy=      156.83 Ry<br>      This function has    0 nodes for 0 < r <    
1.385<br>     ns=   4 l=   3 sum=    0.003498123 r(ikk)     2.754160645<br>     from gener_pseudo  : info #        -1<br>     chi too large beyond r_c<br><br>*********************************************************************************
<br><br>The input file is:<br><br> &input<br>   title='Tb',<br>   zed=65.0,<br>   rel=2,<br>   rmax = 100,<br>   iswitch =3,<br>   rlderiv=3.01,<br>   eminld=-6.0,<br>   emaxld=4.0,<br>   deld=0.02,<br>   nld=5,
<br>   config='[Kr] 4d10 5s2 5p5.5 4f9 6s2 5d0.5',<br>   dft='PBE'<br> &<br> &inputp<br>   pseudotype=3,<br>   lloc=0,<br>   file_pseudopw='Tb.pbe',<br>   nlcc =.true.,<br>   rho0=0.01,<br>
!   rcore= 1.1,<br>!   rcloc= 2.1<br> &<br>8<br>5S  1  0  2.00  0.00  2.20  2.20 0.50<br>5P  2  1  2.00  0.00  2.20  2.20 0.50<br>5P  2  1  3.50  0.00  2.20  2.20 1.50<br>4F  4  3  6.00  0.00  1.40  1.40 2.50<br>4F  4  3  
3.00  0.00  1.40  1.40 3.50<br>5D  3  2  0.50  0.00  2.10  2.10 1.50<br>5D  3  2  0.00  0.00  2.10  2.10 2.50<br>6S  1  0  2.00  0.00  2.60  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>               Beijing, China<br>Phone:    86-10-62582396
<br>Email:     <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>