Page 1 of 1
Positive energy for the case of METAGGA calculations
Posted: Fri May 27, 2022 11:55 am
by chandan_kumarvishwakarma1
Dear Users,
I am trying to run scf calculations using METAGGA functional, but in the output file code is printing positive energy.
I am not able to rectify the problem. Please help me regarding the same.
I am attaching all the inputs and few outputs (OSZICAR and OUTCAR), Please find the attachment.
---
Regards
Chandan
IIT Delhi, India
Re: Positive energy for the case of METAGGA calculations
Posted: Fri May 27, 2022 12:47 pm
by martin.schlipf
You are using the MBJ functional which is not suitable for total energy calculations.
VaspWiki wrote:The MBJ functional is a potential-only functional, i.e., there is no corresponding MBJ exchange-correlation energy, instead E_xc is taken from LDA. This means that MBJ calculations can never be self-consistent with respect to the total energy, and thus we cannot compute Hellmann-Feynman forces (i.e., no ionic relaxation, etc.). Actually, MBJ calculations aim solely at a description of the electronic properties, primarily band gaps. (
source)
Re: Positive energy for the case of METAGGA calculations
Posted: Fri May 27, 2022 1:41 pm
by chandan_kumarvishwakarma1
Dear Martin,
Thanks for very prompt reply, I want to mention one point that I am not using MBJ for energy calculations.
I am using this for electronic properties but in between I came to know that in OSZICAR file positive energy is printing.
This was strange for me so I raise this question.
If you can comment more on that then it will be helpful for me.
---
Thanks
Chandan
Re: Positive energy for the case of METAGGA calculations
Posted: Sat May 28, 2022 8:36 am
by martin.schlipf
Use the MBJ functional to get better band structures. It modifies the potential to open the gap compared to local DFT.
However the change to the potential also affects the total energy so it cannot be used to compare structures or for relaxations anymore. The only reason to still print it is to judge whether your electronic structure is converged.
Re: Positive energy for the case of METAGGA calculations
Posted: Thu Jun 02, 2022 1:59 pm
by chandan_kumarvishwakarma1
Dear Martin,
Thanks for your suggestion and clarifications.
---
Regards
Chandan