My Community

Dear all,

Previously I posted this question in the Bugreports section which is definitely wrong place for such a question

I try to calculate the total charge density by means of Band decomposed chargedensity calculation (i.e. PARCHG file) with SOC included (LSORBIT = .TRUE.). I follow the instructions for the PARCHG calculation in manual. The outcome is very puzzling - when I integrate my PARCHG I see that the number of electrons is smaller than NBANDS.

Does anybody has idea why it is so?

When I do similar calculation without SOC everything is OK.

=========================

Some input details.

1st run - static SCF calculation (only relevant tags)
ISTART = 0
ICHARG = 1
ISYM = 0
#LWAVE = .FALSE.
#LCHARG = .FALSE.
NSW = 0
IBRION = -1

2nd run - PARCHG calculation - EVERYTHING IS THE SAME EXCEPT FOR
ISTART = 1
ICHARG = 11
LPARD = .TRUE.
NBMOD = -2
EINT = -17 # only occupied bands
LSEPB = .FALSE.
LSEPK = .FALSE.
===========================

Any suggestions are greatly appreciated.

Cheers,

Misha

Hi,

We're sorry that we didn’t answer your question. This does not live up to the quality of support that we aim to provide. The team has since expanded. If we can still help with your problem, please ask again in a new post, linking to this one, and we will answer as quickly as possible.

Best wishes,

VASP