Charge Vacancy Calculation Error

Problems running VASP: crashes, internal errors, "wrong" results.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
tommy91779

Charge Vacancy Calculation Error

#1 Post by tommy91779 » Wed Jun 24, 2009 5:08 pm

Hello all,

I am doing a set of vacancy calculations for an AX system where A has 12 electrons and X has 6 electrons. Considing a 3x3x2 supercell, I have 36 AX pairs which gives 648 electrons. For a neutral X vacancy I of course have 642 electrons. I would now like to perform a 1+ charged X vacancy calculation. I set the flag NELECT=641. UNfortunately when I perform the calculation, I obtain the following error:

LDA part: xc-table for Ceperly-Alder, standard interpolation
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: wrap around errors must be expected
FFT: planning ... 4
reading WAVECAR
WARNING: random wavefunctions but no delay for mixing, default for NELMDL
entering main loop
N E dE d eps ncg rms rms(c)
DAV: 1 0.322421237694E+04 0.32242E+04 -0.21631E+05 18008 0.916E+02
DAV: 2 -0.335114282016E+02 -0.32577E+04 -0.31628E+04 19704 0.281E+02
DAV: 3 -0.246281923418E+03 -0.21277E+03 -0.21077E+03 20008 0.679E+01
DAV: 4 -0.249997819101E+03 -0.37159E+01 -0.36964E+01 19648 0.816E+00
DAV: 5 -0.250072145117E+03 -0.74326E-01 -0.74188E-01 19912 0.946E-01BRMIX: very serious problems
the old and the new charge density differ
old charge density: 642.00000 new 641.00001


BElow is my INCAR file

Electronic minimisation
IBRION = 2
ISIF = 2
NSW = 300

LREAL = Auto

NELECT = 641



Does anybody have any insight into my particular problem??


Thanks in advance

-Tom
Last edited by tommy91779 on Wed Jun 24, 2009 5:08 pm, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Charge Vacancy Calculation Error

#2 Post by admin » Fri Jul 17, 2009 11:58 am

this actually is no error: for the initial 5 steps, vasp just 'initializes' the Hamiltonian and the wavefunctions from the POTCAR file (without any charge update, i.e ignoring the NELECT tag).
After the 5th electronic step, NELECT is recognized to be different from the sum of valence electrons as calculated from POTCAR.
The written 'warning' just tells you about this (cell is not charge-neutral), vasp will continue to run for 641 valence electrons.
please keep in mind that for charged systems, it may improve convergence to add set IDIPOL (and even apply multipole corrections, please have a look at the respective chapter of the online manual)
Last edited by admin on Fri Jul 17, 2009 11:58 am, edited 1 time in total.

Post Reply