Trouble with hybrid functional calculations using Ferwe=-2

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Locked
Message
Author
Kousei
Newbie
Newbie
Posts: 10
Joined: Wed Oct 27, 2021 6:19 am

Trouble with hybrid functional calculations using Ferwe=-2

#1 Post by Kousei » Wed Oct 27, 2021 8:28 am

Dear All,

I am currently trying to perform a structural optimization of the excited state of CaF2:Eu2+ using the FERWE tag.
Here is the INCAR of the calculation.

Code: Select all

System = CaF2_Eu2+
ENCUT = 500
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = All
TIME = 0.4
FERWE = 390*1.0 1*0.0 1*1.0 58*0.0 390*1.0 1*0.0 1*1.0 58*0.0 390*1.0 1*0.0 1*1.0 58*0.0 390*1.0 1*0.0 1*1.0 58*0.0
FERDO = 384*1.0 66*0.0 384*1.0 66*0.0 384*1.0 66*0.0 384*1.0 66*0.0
IBRION = 2
ICHARG = 2
ISMEAR = -2
ISTART = 1
ISPIN = 2
ISIF = 2
LORBIT = 11
MAGMOM = 1*7.0 323*0
NUPDOWN = 7
NBANDS = 450
NSW = 50
NPAR = 1
SIGMA = 0.1
The calculation completes without any error messages, but I have identified a problem where the occupancy that should have been specified in the FERWE tag does not hold until the end of the calculation. I would like to know if there is a possible cause for this.

Many thanks in advance,
Kousei

martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: Trouble with hybrid functional calculations using Ferwe=-2

#2 Post by martin.schlipf » Wed Oct 27, 2021 12:44 pm

Did you try whether this problem is specific to your system or holds in general? FERWE should not be updated differently for normal DFT calculations, as far as I know, so it seems you should be able to reproduce this with a simpler system. It is generally a good idea to try the feature that you want in a small system to avoid wasting computational resources if something goes wrong.

For the future, please also remember that we need a complete set of input files to reproduce the issue.

Kousei
Newbie
Newbie
Posts: 10
Joined: Wed Oct 27, 2021 6:19 am

Re: Trouble with hybrid functional calculations using Ferwe=-2

#3 Post by Kousei » Thu Oct 28, 2021 9:31 am

Thank you for your comment.
I followed your suggestion and did the calculation in the unit cell, but still occured the same problem.
I have attached the files and would appreciate your advice.
You do not have the required permissions to view the files attached to this post.

martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: Trouble with hybrid functional calculations using Ferwe=-2

#4 Post by martin.schlipf » Thu Oct 28, 2021 9:22 pm

Thank you for this input. I could reproduce the issue. It seems that there is a workaround by using ALGO = Normal instead, if that works for your system.

Kousei
Newbie
Newbie
Posts: 10
Joined: Wed Oct 27, 2021 6:19 am

Re: Trouble with hybrid functional calculations using Ferwe=-2

#5 Post by Kousei » Mon Nov 01, 2021 4:37 am

Thank you very much.
By changing the ALGO to Normal, the calculation converged well.

martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: Trouble with hybrid functional calculations using Ferwe=-2

#6 Post by martin.schlipf » Wed Nov 10, 2021 8:34 am

We have analyzed this in a bit more detail and updated the wiki entries for ISMEAR, LDIAG, and ALGO.

The bottom line is that you should use ALGO = VeryFast or ALGO = Damped to avoid bands changing the order. In either case you should set LDIAG = F to avoid a diagonalization step that can also change the order of the bands.

Locked