BUG: excited configuration with nqfd not working

Running Abinit Ground state
1

Dear all,
I am trying to do a single point calculation with an excited configuration in which one electron is moved from the top of the valence band to the bottom of the conduction band by means of the keywords
occopt 9
ivalence 184
nqfd 1

After starting the first iteration of the SCF, abinit exits with the following error:

ab.abi (2.7 KB)
out.log (90.6 KB)

— !BUG
src_file: m_occ.F90
src_line: 718
mpi_rank: 0
message: |

 newocc: 
  The calling routine gives nelect=   0.36800000D+03
  The lowest bound is  -0.22772488D+01, with nelect=  0.00000000D+00
  The highest bound is   0.64408137D+00, with nelect=  0.12000000D+02

…

Is it a real bug or it is due to a mistake that I did in the input?
Please, attached you find the input and the standard output. If more information is needed, please, let me know.
Thanks a lot in advance for your help!

Antonio

2

Hi @Antonio,
You need to increase the number of bands in your input file. I got the same problem a while ago and increasing the number of bands (nband) fixed it.

best.

3

Hi @Carl_D ,
I did several tests; I increased it up to the double and I still get the same error. Can you send me your input that worked for you? Thanks.
Best

4

hi @antonio,
I computed the total band in your structure and found 512, I then set nband to 512, I did not check the ivalence value (I guess you got by considering the total valence electrons). I am running the calculation using the input that you gave in the post now. So far it’s running. Try that, once it’s completed I will keep you posted.

5

Hi @antonio,
The calculation went through smoothly. I tried to attach the output files but the forum is saying new user cannot attach files. I paste it here anway.

natom 24
typat 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3
ntypat 3
znucl 42 34 52
acell 1 1 1

rprim 13.3476882827743459 0.0000000000000000 0.0000000000000000
-6.6738441413871730 11.5594371346784737 0.0000000000000000
0.0000000000000000 0.0000000000000000 75.5890453150000070

xred 0.3330727758094342 0.1666666666666667 0.9068393005261008
0.8333333333333334 0.1666666666666667 0.9070998580500000
0.3333333333333333 0.6666666666666666 0.9070998580500000
0.8333333333333333 0.6666666666666666 0.9070998580500000
0.1666666666666667 0.3333333333333333 0.7426337735400002
0.6666666666666666 0.3333333333333333 0.7426337735400002
0.1666666666666667 0.8333333333333334 0.7426337735400002
0.6666666666666666 0.8333333333333334 0.7426337735400002
0.1666666666666667 0.3333333333333333 0.8677638392700000
0.6666666666666666 0.3333333333333333 0.8677638392700000
0.1666666666666667 0.8333333333333334 0.8677638392700000
0.6666666666666666 0.8333333333333334 0.8677638392700000
0.1666666666666667 0.3333333333333333 0.9464008108400001
0.6666666666666666 0.3333333333333333 0.9464008108400001
0.1666666666666667 0.8333333333333334 0.9464008108400001
0.6666666666666666 0.8333333333333334 0.9464008108400001
0.3333333333333333 0.1666666666666667 0.6973940219800001
0.8333333333333334 0.1666666666666667 0.6973940219800001
0.3333333333333333 0.6666666666666666 0.6973940219800001
0.8333333333333333 0.6666666666666666 0.6973940219800001
0.3333333333333333 0.1666666666666667 0.7878200549700000
0.8333333333333334 0.1666666666666667 0.7878200549700000
0.3333333333333333 0.6666666666666666 0.7878200549700000
0.8333333333333333 0.6666666666666666 0.7878200549700000
#Common settings for phonon calculations
#*****************************************
nsym 1

#Parameters of the SCF cycles
#****************************
iscf 17
nstep 10
ixc 11 # PBE
pp_dirpath “/storage/”
pseudos “Mo.psp8, Se.psp8, Te.psp8”
toldfe 1.0d-11
rmm_diis 1

#Plane wave basis and k-point grid
#*********************************
ecut 500 eV
ecutsm 0.5
ngkpt 7 7 1
nshiftk 1
shiftk 0.0 0.0 0.0
diemac 9.0
nband 512
kptopt 1

occopt 9
ivalence 184
nqfd 1

notice that I changed nstep to 10 to get it finish quickly.

6

@Carl_D thanks a lot for your test! I will use that value that seems to work for other similar systems too.
To avoid an overuse of the computational resources, it would be important to know what is the minimum number of empty bands that should be considered in order to avoid such bug. I will try to figure it out and I’ll post it here if I find it.

Thanks again for your help.

Best