<div style="inset: 0px auto auto 0px; position: absolute; overflow: hidden; z-index: 1; display: flex; flex-direction: column; width: 1920px; height: 1007px;"><span role="main" style="contain: style; display: block; overflow-y: auto; z-index: 1; height: 943px; margin-left: 280px; transition: margin-left 0.25s cubic-bezier(0.4, 0, 0.2, 1) 0s, visibility 0s linear 0s; visibility: visible; opacity: 1;"><span style="contain: style; display: block; height: 943px; margin-inline-start: 16px; overflow: hidden; z-index: 1;"><div style="display: flex; flex-direction: column; height: 943px; position: relative;"><div style="overflow-y: auto; outline: none;"><div role="list" aria-label="Convergence problem with cobalt oxide Cu3O4"><span tabindex="0" role="listitem" aria-expanded="true" style="border-bottom: none; padding-top: 8px; padding-left: 0px; border-left: 2px solid rgb(77, 144, 240);"><div tabindex="-1" style="outline: none;"><div style="display: flex;"><div style="flex-grow: 1; min-width: 0px;"><div role="region" aria-labelledby="c243" style="margin: 12px 0px; overflow: auto; padding-right: 20px;">Dear CP2K users,<div><br /></div><div>I am currently working with a series of transition metal oxides and CP2K is working fine but when I try to converge a calculation on the mixed valence Cu oxide (Cu3O4) my calculations are not converging. I have read that this oxide is a tricky system for DFT, so I would be really grateful if you could give me some suggestions to solve this problem. I attatch the input file that I am using. I tried to use different mixing, diagonalization and smearing methods with the same problem.</div><div><br /></div><div>&GLOBAL<br />  PROJECT Co3O4<br />  RUN_TYPE ENERGY<br />  PRINT_LEVEL MEDIUM<br />&END GLOBAL<br /><br />&FORCE_EVAL<br />  METHOD Quickstep<br />  &DFT<br />    BASIS_SET_FILE_NAME  BASIS_MOLOPT_UZH<br />    POTENTIAL_FILE_NAME  POTENTIAL_UZH<br />    UKS<br />    MULTIPLICITY 1<br />    &POISSON<br />      PERIODIC XYZ<br />    &END POISSON<br />    &XC<br />      &XC_FUNCTIONAL PBE<br />       &PBE<br />        PARAMETRIZATION PBESOL<br />       &END<br />      &END XC_FUNCTIONAL<br />    &END XC<br />    &MGRID<br />      NGRIDS 4<br />      CUTOFF 800<br />      REL_CUTOFF 50<br />    &END MGRID<br />    &QS<br />      METHOD GPW<br />      EXTRAPOLATION ASPC<br />      EPS_DEFAULT 1E-014<br />    &END QS<br />    &SCF<br />      SCF_GUESS ATOMIC<br />      EPS_SCF 1.0E-7<br />      ADDED_MOS 20<br />      MAX_SCF 20<br />      &OUTER_SCF<br />       MAX_SCF 300<br />       EPS_SCF 1.0E-7<br />      &END<br />      &SMEAR ON<br />        METHOD FERMI_DIRAC<br />        ELECTRONIC_TEMPERATURE [K] 3000<br />      &END SMEAR<br />      &DIAGONALIZATION<br />        ALGORITHM STANDARD<br />      &END DIAGONALIZATION<br />      &MIXING<br />        METHOD BROYDEN_MIXING<br />        ALPHA 0.02<br />        BETA 0.01<br />      &END MIXING<br />     &END SCF<br />  &END DFT<br />  &SUBSYS<br />    &CELL<br />      ABC 5.938 5.945 5.938<br />      ALPHA_BETA_GAMMA 60 60 60<br />      PERIODIC XYZ<br />    &END CELL<br />    &TOPOLOGY<br />      COORD_FILE_FORMAT xyz<br />      COORD_FILE_NAME  ./Co3O4-2.xyz<br />    &END<br />    &KIND O<br />      ELEMENT O<br />      BASIS_SET DZVP-MOLOPT-PBE-GTH-q6<br />      POTENTIAL GTH-PBE-q6<br />      &BS<br />      &ALPHA<br />      N     2<br />      L     1<br />      NEL  +2<br />      &END ALPHA<br />      &BETA<br />      N     2<br />      L     1<br />      NEL  +2<br />     &END BETA<br />     &END BS<br />    &END KIND<br />    &KIND Co<br />      ELEMENT Co<br />      BASIS_SET DZVP-MOLOPT-PBE-GTH-q17<br />      POTENTIAL GTH-PBE-q17<br />     &BS<br />      &ALPHA<br />       N    4  3<br />       L    0  2<br />       NEL -2 -1<br />      &END ALPHA<br />      &BETA<br />       N    4  3<br />       L    0  2<br />       NEL -2 -1<br />      &END BETA<br />      &END BS<br />      &DFT_PLUS_U<br />       L 2<br />       U_MINUS_J [eV] 2.0<br />     &END DFT_PLUS_U<br />    &END KIND<br />    &KIND Co1<br />      ELEMENT Co<br />      BASIS_SET DZVP-MOLOPT-PBE-GTH-q17<br />      POTENTIAL GTH-PBE-q17<br />     &BS<br />      &ALPHA<br />       N    4  3<br />       L    0  2<br />       NEL -2 -3<br />      &END ALPHA<br />      &BETA<br />       N    4  3<br />       L    0  2<br />       NEL -2 +3<br />      &END BETA<br />      &END BS<br />      &DFT_PLUS_U<br />       L 2<br />       U_MINUS_J [eV] 2.0<br />     &END DFT_PLUS_U<br />    &END KIND<br />    &KIND Co2<br />      ELEMENT Co<br />      BASIS_SET DZVP-MOLOPT-PBE-GTH-q17<br />      POTENTIAL GTH-PBE-q17<br />     &BS<br />      &BETA<br />       N    4  3<br />       L    0  2<br />       NEL -2 -3<br />      &END BETA<br />      &ALPHA<br />       N    4  3<br />       L    0  2<br />       NEL -2 +3<br />      &END ALPHA<br />      &END BS<br />      &DFT_PLUS_U<br />       L 2<br />       U_MINUS_J [eV] 2.0<br />     &END DFT_PLUS_U<br />    &END KIND<br />  &END SUBSYS<br />&END FORCE_EVAL<br /></div><div><br /></div><div>Thank you in advanced.</div><div>Best regards,</div><div>Diego.</div></div></div></div></div></span></div></div></div></span></span></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/a81091b7-fc7e-4089-939a-7a6e1823a25an%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/a81091b7-fc7e-4089-939a-7a6e1823a25an%40googlegroups.com</a>.<br />