[Pw_forum] More pseudopotential troubles

Miguel Martínez Canales wmbmacam at lg.ehu.es
Wed Oct 11 12:15:53 CEST 2006


Dear people,

Although your previous comments on Ca pseudopotential were very helpful, 
now I'm finding another problem. I may be wrong to think that the input 
attached generates an US pseudopotential with semicore states, using l=2 as 
the pseudization channel and two reference energies for the s and p 
channels. Anyway, when I try to run this input

  &input
    atom='Ca',
    beta = 0.2,
    rel=1,
    config='[Ar] 3d0 4s2 4p0',
    iswitch=3,
    dft='PBE'
  /
  &inputp
    pseudotype=3,
    lloc=2,
    nlcc =.true.,
    file_pseudopw='Ca.pbe-nsp-1.UPF',
    rho0=0.01
/
5
3S  1  0  2.00  0.00  1.10  1.30
3S  2  0  0.00 -0.27  2.80  2.90
3P  2  1  6.00  0.00  1.20  1.30
3P  3  1  0.00 -0.10  2.80  2.90
3D  3  2  0.00  0.00  2.90  2.90
  &test
    nconf = 1,
    file_pseudo='Ca.pbe-nsp-1.UPF',
    ecutmin = 20.0,
    ecutmax = 60.0,
/
3
3S  2  0  2.00  0.00  2.80  2.90
3P  3  1  0.00  0.00  2.80  2.90
3D  3  2  0.00  0.00  2.90  2.90

I get the following error:
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
      from ascheqps : error #         1
      too many attempts
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

      stopping ...

Increasing beta doesn't help. Neither does increasing the maximum number of 
iterations in ascheqps.f90. I'm currently looking at that module to 
understand what's going wrong. Nevertheless, if you saw something stupid in 
my input, pleas tell me.

Thanks in advance,

Miguel

PS: In the excited state for the s channel, 3S 2 0, do I have to write the 
pseudization energy? Or is 0.00 the corresponding energy to the excited state?



More information about the Pw_forum mailing list