UT theoretical chemistry code forum

Hi all,

Now I am doing Bader analysis on a 100-atom crystal unit. I think I ran into some kind of trouble.
1. run "bader CHGCAR -ref CHGCAR_sum" I got "ERROR: should be no new maxima in edge refinement"
2. from OUTCAR I know my unit cell should have 660 electrons.
Bader calculation gave me "NUMBER OF ELECTRONS: 634.73046"
If I continue to use a higher NGX/NGY/NGZ, the results are:
=========================================================
NGX = 100
NGY = 100
NGZ = 240
-----------------------------------------------------------------------------------------------------
NUMBER OF BADER MAXIMA FOUND: 90
SIGNIFICANT MAXIMA FOUND: 90
VACUUM CHARGE: -0.1817
NUMBER OF ELECTRONS: 635.76616
==========================================================
NGX =165
NGY =165
NGZ =395
-----------------------------------------------------------------------------------------------------
NUMBER OF BADER MAXIMA FOUND: 90
SIGNIFICANT MAXIMA FOUND: 90
VACUUM CHARGE: -0.1817
NUMBER OF ELECTRONS: 636.57041
==========================================================

It seems like something is not right. If you post the files I'll take a look at the analysis.

Dear Dr. Graeme,

Just would like to let you know that I have solved my problem by running a previous version bader (Version 0.28a 07/12/12).

I can get output 'NUMBER OF ELECTRONS: 660.00000', which corresponds OUTCAR. Also ACF.dat is slightly different from previous result.

Thanks for the note; there was a bug in the new version which did not properly deal with charge density files containing negative values of the charge density. I'll update things right now.