VASPsol convergence issue for CeO2(111)

[neilsfantasy]

Dear Prof. Hennig,

I tried applying VASPsol to a cerium dioxide surface (CeO2(111); clean, without any adsorbate on it), and the wave function failed to converge despite my best efforts. To my knowledge, it works well for many other oxides such as titania, nickel oxide, and iridium oxide, and indeed I tried VASPsol with titanium using the same settings and did not have any convergence issues. Below I attach the relevant files for the ceria job. I would very much appreciate any help you can give me for getting this system converged.

POTCAR (using the 2012 PAW_PBE set)

KPOINTS:
Automatic mesh
0
Gamma
2 2 1
0 0 0

INCAR:
System = CeO2(111) 4x4x3 slab

LSOL = .TRUE.

NPAR = 1
GGA = 91
VOSKOWN= 1
ENCUT = 400.00 eV
LPLANE = .TRUE.
LREAL = Auto
ALGO = Fast
PREC = Normal
NELM = 70
EDIFF = 1e-5
IDIPOL = 3
ISMEAR = 0
SIGMA = 0.05
ISPIN = 1
ISIF = 2
NSW = 0
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 3 -1
LDAUU = 6 0
LDAUJ = 1 0

POSCAR:
CeO2(111)
3.884844655838890
2.0000000000000000 0.0000000000000000 0.0000000000000000
-1.0000000000000000 1.7320508075688772 0.0000000000000000
0.0000000000000000 0.0000000000000000 5.1775494819768904
Ce O
12 24
Selective dynamics
Direct
0.1666666666666643 0.3333333333333357 0.0394248564774671 F F F
0.6666666666666643 0.3333333333333357 0.0394248564774671 F F F
0.1666666666666643 0.8333333333333357 0.0394248564774671 F F F
0.6666666666666643 0.8333333333333357 0.0394248564774671 F F F
0.3333333333333357 0.1666666666666643 0.1971242823873425 F F F
0.8333333333333357 0.1666666666666643 0.1971242823873425 F F F
0.3333333333333357 0.6666666666666643 0.1971242823873425 F F F
0.8333333333333357 0.6666666666666643 0.1971242823873425 F F F
0.0000000000000000 0.0000000000000000 0.3532008151332775 T T T
0.5000000000000000 0.0000000000000000 0.3532008151332775 T T T
0.0000000000000000 0.5000000000000000 0.3532008151332775 T T T
0.5000000000000000 0.5000000000000000 0.3532008151332775 T T T
0.0000000000000000 0.0000000000000000 0.0000000000000000 F F F
0.5000000000000000 0.0000000000000000 0.0000000000000000 F F F
0.0000000000000000 0.5000000000000000 0.0000000000000000 F F F
0.5000000000000000 0.5000000000000000 0.0000000000000000 F F F
0.3333333333333357 0.1666666666666643 0.0788497129549413 F F F
0.8333333333333357 0.1666666666666643 0.0788497129549413 F F F
0.3333333333333357 0.6666666666666643 0.0788497129549413 F F F
0.8333333333333357 0.6666666666666643 0.0788497129549413 F F F
0.1666666666666643 0.3333333333333357 0.1576994259098754 F F F
0.6666666666666643 0.3333333333333357 0.1576994259098754 F F F
0.1666666666666643 0.8333333333333357 0.1576994259098754 F F F
0.6666666666666643 0.8333333333333357 0.1576994259098754 F F F
0.0000000000000000 0.0000000000000000 0.2365491388648095 F F F
0.5000000000000000 0.0000000000000000 0.2365491388648095 F F F
0.0000000000000000 0.5000000000000000 0.2365491388648095 F F F
0.5000000000000000 0.5000000000000000 0.2365491388648095 F F F
0.3333333333333357 0.1666666666666643 0.3137089924832708 T T T
0.8333333333333357 0.1666666666666643 0.3137089924832708 T T T
0.3333333333333357 0.6666666666666643 0.3137089924832708 T T T
0.8333333333333357 0.6666666666666643 0.3137089924832708 T T T
0.1666666666666643 0.3333333333333357 0.3918096184401207 T T T
0.6666666666666643 0.3333333333333357 0.3918096184401207 T T T
0.1666666666666643 0.8333333333333357 0.3918096184401207 T T T
0.6666666666666643 0.8333333333333357 0.3918096184401207 T T T

Excerpt of OSZICAR:
N E dE d eps ncg rms rms(c)
SOL: 1 -0.52812E+03 0.60794E-01 -0.52806E+03 103
DAV: 1 -0.820702626470E+03 -0.82070E+03 -0.32056E+01 1280 0.181E+01 0.988E+00
SOL: 2 -0.81573E+03 0.61764E-01 -0.81566E+03 119
RMM: 2 -0.110938402880E+04 -0.28868E+03 -0.54493E+01 1092 0.186E+01 0.556E+00
SOL: 3 -0.10752E+04 0.63630E-01 -0.10752E+04 121
RMM: 3 -0.377120855499E+04 -0.26618E+04 0.26545E+05 1096 0.392E+03 0.335E+02
SOL: 4 -0.10999E+04 0.22198E+00 -0.10997E+04 121
RMM: 4 -0.860774897451E+03 0.29104E+04 -0.32067E+03 1289 0.169E+02 0.147E+02
SOL: 5 -0.10494E+04 0.15546E+00 -0.10492E+04 118
RMM: 5 -0.112374524039E+04 -0.26297E+03 -0.26885E+03 1271 0.161E+02 0.347E+01
SOL: 6 -0.10422E+04 0.95755E-01 -0.10421E+04 112
RMM: 6 0.300259658379E+04 0.41263E+04 0.10527E+05 865 0.283E+03 0.261E+02
SOL: 7 -0.87971E+03 0.20732E+00 -0.87950E+03 117
RMM: 7 -0.918736494701E+03 -0.39213E+04 -0.66175E+03 1072 0.443E+02 0.186E+02
SOL: 8 -0.91527E+03 0.16960E+00 -0.91510E+03 116
RMM: 8 -0.103056369799E+04 -0.11183E+03 -0.40561E+03 1320 0.178E+02 0.508E+01
SOL: 9 -0.10156E+04 0.13218E+00 -0.10155E+04 119
RMM: 9 0.567275679760E+04 0.67033E+04 0.22001E+05 848 0.379E+03 0.271E+02
SOL: 10 -0.11002E+04 0.21664E+00 -0.10999E+04 120
RMM: 10 0.579347565339E+04 0.12072E+03 -0.31253E+04 782 0.884E+03 0.667E+02
SOL: 11 -0.12206E+04 0.27399E+00 -0.12203E+04 122
RMM: 11 -0.963696014479E+03 -0.67572E+04 -0.55210E+03 1195 0.252E+02 0.108E+02
SOL: 12 -0.10984E+04 0.24188E+00 -0.10981E+04 120
RMM: 12 -0.174837438657E+04 -0.78468E+03 -0.87259E+03 868 0.806E+02 0.116E+02
SOL: 13 -0.10492E+04 0.27924E+00 -0.10489E+04 117
RMM: 13 -0.104449713950E+04 0.70388E+03 -0.22427E+03 1292 0.149E+02 0.695E+01
SOL: 14 -0.10493E+04 0.23889E+00 -0.10491E+04 113
RMM: 14 -0.218489173678E+04 -0.11404E+04 -0.61078E+03 982 0.708E+02 0.120E+02
SOL: 15 -0.91098E+03 0.17805E+00 -0.91080E+03 104
RMM: 15 -0.100270769556E+04 0.11822E+04 -0.26748E+03 1279 0.149E+02 0.649E+01
SOL: 16 -0.80444E+03 0.16273E+00 -0.80427E+03 119
RMM: 16 0.566190980108E+04 0.66646E+04 0.25522E+05 1164 0.646E+03 0.680E+02
SOL: 17 -0.85792E+03 0.28412E+00 -0.85763E+03 123
RMM: 17 -0.483082404742E+03 -0.61450E+04 -0.61611E+03 1274 0.265E+02 0.127E+02
SOL: 18 -0.69718E+03 0.17535E+00 -0.69700E+03 125
RMM: 18 0.250364473675E+06 0.25085E+06 0.25895E+06 921 0.245E+04 0.221E+03
SOL: 19 -0.11238E+04 0.42367E+00 -0.11234E+04 124
RMM: 19 0.290791016175E+04 -0.24746E+06 -0.25034E+04 1298 0.151E+03 0.132E+03
SOL: 20 -0.92575E+03 0.32067E+00 -0.92543E+03 121
RMM: 20 0.121484422355E+04 -0.16931E+04 -0.27968E+04 867 0.338E+03 0.122E+03
SOL: 21 -0.61003E+03 0.33302E+00 -0.60970E+03 122
RMM: 21 -0.353509281836E+02 -0.12502E+04 -0.21922E+04 1348 0.132E+03 0.106E+03
...

I used previously converged WAVECAR and CHGCAR without VASPsol but that did not make any difference. I also tried various other things, such as turning off DFT+U, changing to gamma-point only, increasing the k-point density or cutoff energy, using an earlier set of PAW potentials, but nothing has helped.

Best,
Neil

[Crossedsticks]

Hello,

has this problem ever been solved? I have the same problem on my CeO2(111) surfaces and would like to know if there is anything one can do.