ASR violation in AlAs fermi smearing

Phonons, DFPT, electron-phonon, electric-field response, mechanical response…

Moderators: mverstra, joaocarloscabreu

Locked
KonradG
Posts: 2
Joined: Tue Feb 16, 2021 10:13 pm

ASR violation in AlAs fermi smearing

Post by KonradG » Tue Feb 16, 2021 10:34 pm

Hello

I am having trouble running phonon calculations with metallic compounds and I went back to the tutorial and tried running AlAs as a metal (occopt 3) and ended up with huge differences in the DFPT calculations.

Am I missing something obvious? I am using the standard pseudopotentials from the abinit tutorial and all I changed was the occupation scheme between the two runs.

Thank You
Attachments
AlAsOccopt3.out
(67.34 KiB) Downloaded 301 times
AlAsOccopt1.out
(55.76 KiB) Downloaded 295 times
AlAsOccopt3.in
(725 Bytes) Downloaded 296 times
AlAsOccopt1.in
(727 Bytes) Downloaded 292 times

KonradG
Posts: 2
Joined: Tue Feb 16, 2021 10:13 pm

Re: ASR violation in AlAs fermi smearing

Post by KonradG » Thu Apr 15, 2021 4:36 pm

Hello, I have since solved my issue.

The problem lies with the default value of nbdbuf. It automatically assumes that there are 2 or more unoccupied bands in the ground state but for small systems this is not true. For AlAs, the code was freezing an occupied band when doing RF calculations thus giving incorrect results.

The solution is to set nbdbuf to 0.

Locked