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 4 3 3<br>classical turning point too close to mesh 4 3 -4<br>classical turning point too close to mesh 5 2 2
<br>classical turning point too close to mesh 5 2 -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 clear="all"><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>