Oscillating Energy for Ionic Relaxation - HSE functional

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: Moderator, Global Moderator

Post Reply
Message
Author
shaun_tay
Newbie
Newbie
Posts: 2
Joined: Wed Oct 18, 2023 10:35 pm

Oscillating Energy for Ionic Relaxation - HSE functional

#1 Post by shaun_tay » Mon May 27, 2024 1:21 am

Dear VASP community,

I am currently running HSE functional calculations for NaFeFe(CN)6 using Fe_pv, with the specific defined spin states of 3.8 and 5.0 for the first and second iron sites, respectively.

Initially, I optimized my structure using the GGA functional with the Gamma point, and this step was successful. However, when I switched to the HSE functional (without using the previous WAVECAR), the SCF calculations began to oscillate back and forth during the first ionic step and failed to converge.

Despite troubleshooting the issue by setting LREAL to Auto as suggested in the OUTCAR file, the problem still persists. I have attached my INCAR and input files for your reference.

I would greatly appreciate any advice or guidance on how to resolve this issue. Thank you all :)

INCAR
ALGO = Normal
EDIFF = 1e-05
EDIFFG = -0.02
ENAUG = 1360
ENCUT = 680
GGA = Pe
HFSCREEN = 0.2
IBRION = 2
ISIF = 3
ISMEAR = 0
ISPIN = 2
LAECHG = True
LASPH = True
LCHARG = True
LDAU = False
LELF = False
LHFCALC = True
LMAXMIX = 4
LMIXTAU = True
LORBIT = 11
LREAL = False
LVTOT = False
LWAVE = False
MAGMOM = 4*0.6 4*3.8 4*5.0 48*0.6
NELM = 200
NSW = 99
PREC = Accurate
PRECFOCK = Fast
SIGMA = 0.2
You do not have the required permissions to view the files attached to this post.

henrique_miranda
Global Moderator
Global Moderator
Posts: 474
Joined: Mon Nov 04, 2019 12:41 pm
Contact:

Re: Oscillating Energy for Ionic Relaxation - HSE functional

#2 Post by henrique_miranda » Tue May 28, 2024 3:57 pm

Hard to say what is the reason for this.
I don't see anything obviously wrong in your INCAR.

You say you relaxed the cell with GGA for the Gamma point, but in your HSE calculation you have a 2x2x2 k-point sampling.
Would running a Gamma point only calculation with HSE make sense as well?

Maybe you can try setting `PRECFOCK=Normal` instead of `PRECFOCK=Fast`
The advice to use LREAL=Auto is probably related to the relatively large size of your cell, but it is unlikely to affect the electronic convergence.

shaun_tay
Newbie
Newbie
Posts: 2
Joined: Wed Oct 18, 2023 10:35 pm

Re: Oscillating Energy for Ionic Relaxation - HSE functional

#3 Post by shaun_tay » Tue Jun 04, 2024 1:07 pm

Hi Henrique,

Thank you for the tip!
You say you relaxed the cell with GGA for the Gamma point, but in your HSE calculation you have a 2x2x2 k-point sampling.
Would running a Gamma point only calculation with HSE make sense as well?
That makes sense. However, I was able to run the structure with a 2 x 2 x 2 k-point sampling so I suspect that k-points is not the issue.
Maybe you can try setting `PRECFOCK=Normal` instead of `PRECFOCK=Fast`
The calculation did not seem to converge and is slower than `PRECFOCK=Fast`.

However, upon setting ALGO = Normal, IMIX = 1.0, and AMIX = 0.4, the energy values no longer oscillate but instead increase to a very large positive value. Do you have any advice on how I can approach this problem? Thank you!

henrique_miranda
Global Moderator
Global Moderator
Posts: 474
Joined: Mon Nov 04, 2019 12:41 pm
Contact:

Re: Oscillating Energy for Ionic Relaxation - HSE functional

#4 Post by henrique_miranda » Thu Jun 06, 2024 2:59 pm

You are using the Davidson to diagonalize the Hamiltonian (ALGO=Normal) and then mix the charge density, you might consider using conjugate gradient (ALGO=All) with TIME=0.05.
https://www.vasp.at/wiki/index.php/Trou ... onvergence

It is not guaranteed that it will work, but worth a try.

Post Reply