RMM convergence issues with LEPSILON

Problems running VASP, runtime errors, technical questions.
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

RMM convergence issues with LEPSILON

Postby jcg_forum » Fri May 18, 2018 4:07 pm

Hi all,

I'm having some troubles properly converging linear response calculations using the LEPSILON tag. I have a slightly complex system (layered perovskite), which is elongated in one direction. In this direction, the RMM algo used to solve the Sternheimer equation including local field effects at the DFT level seems not to converge. Here a short sample of an OSZICAR:

Code: Select all
 Linear response to external field, progress :
  Direction:   1
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.850584656635E+01   -0.85058E+01   -0.26103E+01 13340   0.423E+00
RMM:   2    -0.874603943914E+01   -0.24019E+00   -0.19829E-01 11342   0.968E-01    0.597E+00
RMM:   3    -0.748967270442E+01    0.12564E+01   -0.10579E+00 12067   0.138E+00    0.292E+00
RMM:   4    -0.730654704832E+01    0.18313E+00   -0.12068E+00 12239   0.130E+00    0.111E+00
RMM:   5    -0.731457232913E+01   -0.80253E-02   -0.13850E+00 12843   0.685E-01    0.524E-01
RMM:   6    -0.726168940721E+01    0.52883E-01   -0.12933E+00 12661   0.258E-01    0.365E-01
RMM:   7    -0.711357711134E+01    0.14811E+00   -0.24243E+00 11848   0.345E-01    0.164E-01
RMM:   8    -0.710307277916E+01    0.10504E-01   -0.32406E+00 13050   0.162E-01    0.102E-01
RMM:   9    -0.701358586939E+01    0.89487E-01   -0.27760E+00 12605   0.941E-02    0.457E-02
RMM:  10    -0.705063689592E+01   -0.37051E-01   -0.32854E+00 12669   0.475E-02    0.220E-02
RMM:  11    -0.700366510845E+01    0.46972E-01   -0.28143E+00 12758   0.276E-02    0.563E-03
RMM:  12    -0.705342455625E+01   -0.49759E-01   -0.33081E+00 13589   0.162E-02    0.343E-03
RMM:  13    -0.700526629954E+01    0.48158E-01   -0.28231E+00 13420   0.120E-02    0.941E-04
RMM:  14    -0.705390675078E+01   -0.48640E-01   -0.33095E+00 15276   0.951E-03    0.666E-04
RMM:  15    -0.700535516486E+01    0.48552E-01   -0.28229E+00 15425   0.836E-03    0.237E-04
RMM:  16    -0.705400285065E+01   -0.48648E-01   -0.33092E+00 15764   0.778E-03    0.148E-04
RMM:  17    -0.700534540025E+01    0.48657E-01   -0.28228E+00 15958   0.749E-03    0.135E-04
RMM:  18    -0.705399439894E+01   -0.48649E-01   -0.33093E+00 16055   0.736E-03    0.883E-05
...
RMM: 147    -0.700538452600E+01    0.48607E-01   -0.28232E+00 16132   0.721E-03    0.676E-05
RMM: 148    -0.705398501495E+01   -0.48600E-01   -0.33092E+00 16089   0.720E-03    0.535E-05
RMM: 149    -0.700536610283E+01    0.48619E-01   -0.28230E+00 16136   0.720E-03    0.677E-05
RMM: 150    -0.705398678068E+01   -0.48621E-01   -0.33093E+00 16144   0.721E-03    0.535E-05


In this case, the algo did not converged in the maximum number of iterations allowed, with a charge-sloshing-like behavior.

Sometimes this behavior is visible, but the algorithm "converges", without reaching the required EDIFF (I'm not sure EDIFF makes sense for this solver, actually). I've checked the code, there's another alternative criterion of "convergence" for the residue minimization algo of LEPSILON, based on the difference in residue between the current iteration and the two previous one. I don't know if that is a reliable convergence criterion - I think that it's more of a way to see if the solver just won't improve things with more iterations. Actually in the OUTCAR, the same message is always printed ("aborting loop because EDIFF is reached"), which is incorrect.

I therefore have three questions:

  • How sensitive are results of linear response computations to the RMM convergence? Can we trust poorly converged results (I could try that on my own in fact)?
  • Is there a way to improve the RMM convergence? Tightening the kmesh seems to help a bit, but not that much. What about mixing parameters - do these apply to this very specific algo?
  • I have absolutely no issues converging a standard scf on these systems. I also didn't have issues with RMM convergence in LEPSILON runs on other systems (bulk Si, etc). Is this convergence issue something often observed?

many thanks,
Arthur
jcg_forum
Newbie
Newbie
 
Posts: 1
Joined: Thu May 17, 2018 10:49 pm

Return to Bugreports

Who is online

Users browsing this forum: No registered users and 0 guests