Crossedsticks
commented
2 years ago 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.
Open neilsfantasy opened 6 years ago
Crossedsticks
commented
2 years ago 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.
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