Convergence problems with dipole correction

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

Post Reply
Message
Author
weronika_walkosz3
Newbie
Newbie
Posts: 6
Joined: Wed Jun 02, 2021 2:45 pm

Convergence problems with dipole correction

#1 Post by weronika_walkosz3 » Tue Jun 21, 2022 6:29 pm

Hello,
I am trying to model Rh(111) in the presence of oxygen adsorbed at one side only. Hence, I am using dipole corrections (LDIPOL = T, IDIPOL = 3). With increasing O coverage, I face energy convergence problems. Can someone comment on how to overcome this issue? Some of my jobs converged when dipole correction was not used.

fabien_tran1
Global Moderator
Global Moderator
Posts: 316
Joined: Mon Sep 13, 2021 11:02 am

Re: Convergence problems with dipole correction

#2 Post by fabien_tran1 » Tue Jun 21, 2022 7:08 pm

Hi,

Could you please provide the INCAR, POSCAR, KPOINTS and OUTCAR files.

weronika_walkosz3
Newbie
Newbie
Posts: 6
Joined: Wed Jun 02, 2021 2:45 pm

Re: Convergence problems with dipole correction

#3 Post by weronika_walkosz3 » Tue Jun 21, 2022 9:53 pm

Thanks for the reply; I appreciate it.
How do I attach the files (is it possible) to the response? OUTCAR is too big to share otherwise.

fabien_tran1
Global Moderator
Global Moderator
Posts: 316
Joined: Mon Sep 13, 2021 11:02 am

Re: Convergence problems with dipole correction

#4 Post by fabien_tran1 » Wed Jun 22, 2022 7:30 am

The most convenient way is to create a zip file and then attach it by clicking on "Attachments". If the OUTCAR is too large, then just grep for TOTEN in OUTCAR ("grep TOTEN OUTCAR") and display us the results (I just want to see how the total energy behaves along the iterations).

weronika_walkosz3
Newbie
Newbie
Posts: 6
Joined: Wed Jun 02, 2021 2:45 pm

Re: Convergence problems with dipole correction

#5 Post by weronika_walkosz3 » Wed Jun 22, 2022 3:48 pm

Attached are the files. Thanks so much in advance for your help.
You do not have the required permissions to view the files attached to this post.

fabien_tran1
Global Moderator
Global Moderator
Posts: 316
Joined: Mon Sep 13, 2021 11:02 am

Re: Convergence problems with dipole correction

#6 Post by fabien_tran1 » Thu Jun 23, 2022 8:30 am

Hi,

Your calculation with dipole correction stopped because of "Error EDDDAV: Call to ZHEGV failed ..." (last line in OUTCAR). At iteration 37, the density and orbitals were so strange that the diagonalization of the Hamiltonian failed. To some extent this may be a question of luck, since this error could have happen (or not happen) at any early stage of the calculation when the density/orbitals have still not acquired a reasonable shape. For the moment, a suggestion is to restart the calculation with a different mixing scheme (see IALGO) or changing AMIX or BMIX for instance.

By the way, did you restart the calculation with dipole correction using the WAVECAR from the calculation without dipole correction as suggested at LDIPOL (where the DIPOL tag is also discussed)?

There is also this related (but old) topic here: https://www.vasp.at/forum/viewtopic.php?f=4&t=13121

weronika_walkosz3
Newbie
Newbie
Posts: 6
Joined: Wed Jun 02, 2021 2:45 pm

Re: Convergence problems with dipole correction

#7 Post by weronika_walkosz3 » Mon Jun 27, 2022 6:52 pm

Thanks so much.

Post Reply