ABINIT Discussion Forums

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

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.