8
$\begingroup$

Since I have started working on monolayers I was using the ISIF=3 tag for relaxation, this was after recompiling VASP with the modified "constr_cell_relax.F" file. There was no problem in relaxing pure materials, but with the doping of 4x4x1 supercells I don't get convergence in relaxation and the calculations are taking too too much time. On the other hand my friend is using the tag ISIF=4 and he is getting always converged results.

Now I am a little bit confused with which tag should be used, since I think that because of vacuum along the c axis there should be no relaxation along it, so ISIF=4 tag shouldn't be used. I want to know your opinions, especially if you have worked on relaxing doped monolayers.

$\endgroup$
8
  • 1
    $\begingroup$ The convergence problems may be related to ALGO and IBRION, rather than ISIF. $\endgroup$
    – Jack
    Feb 17, 2021 at 4:31
  • 2
    $\begingroup$ Take a look at this post: mattermodeling.stackexchange.com/questions/2022/…. By the way, the ISIF=3 tag with recompiled VASP [ modified constr_cell_relax.F] is suitable for monolayer relaxation. Also this post: mattermodeling.stackexchange.com/questions/2032/… $\endgroup$
    – Jack
    Feb 17, 2021 at 10:46
  • 2
    $\begingroup$ @NikeDattani Thanks for your suggestions. Later I will post an answer. $\endgroup$
    – Jack
    Sep 5, 2021 at 23:41
  • 1
    $\begingroup$ @Jack any luck coming up with an answer. Also, I was able to edit your 3 comments into just 1. $\endgroup$
    – Tyberius
    Nov 19, 2021 at 2:17
  • 1
    $\begingroup$ @ChiKou were you able to figure out any of this? Is there any update you can give us? $\endgroup$ Jan 13 at 0:53

1 Answer 1

4
$\begingroup$

When dealing with monolayers, I have found that using the ISIF=4 tag is better and is less time-consuming. As mentioned in Vaspwiki, it allows the relaxation of both lattice constants and atomic coordinates with a fixed volume. However, what is important here is uncommenting the FCELL(3,3)=SAVE line in constr_cell_relax.F file and recompiling VASP. This trick prevents the relaxation along the z-direction, which is negligible in monolayers.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.