Page 1 of 1

Question regarding MBJ functional?

Posted: Mon Jun 26, 2023 1:34 pm
by bikash.patra
Dear all,

In the case of MBJ calculation, let's say I have specified in my input file the following two lines.
GGA = PE
METAGGA= MBJ

Does this mean the exchange-correlation energy functional is taken from PBE?


With regards,
Bikash.

Re: Question regarding MBJ functional?

Posted: Mon Jun 26, 2023 2:08 pm
by fabien_tran1
Hi,

In general, it does not make sense to specify both GGA and METAGGA tags simultaneously in INCAR, because METAGGA will be used and GGA ignored. So, Either GGA or METAGGA should be used. Actually, since vasp.6.4.0 the simultaneous use of GGA and METAGGA leads to a stop of the calculation.

When METAGGA=MBJ is chosen, the xc energy functional is LDA and this can not be changed because this is hard coded. The LDA call is "CALL CORPBE ..." (only the LDA component of PBE is considered) in the block "ELSEIF (ID_METAGGA==30)" in metagga.F.