Dear cp2k users,<div><br></div><div>I built a Mg slab with Fe attachment using vesta and tried to optimize it with cp2k. It works well for spin restricted calculation but when it comes to spin unrestricted calculation, it didn't converge.  I set UKS .True. , multipilicity =5( spin moment is roughly 2.3 for iron and 2 iron in my slab), maginezation=2.3. </div><div><br></div><div>Here is my input.  I appreciated all your suggestions.</div><div><br></div><div>Best regards,</div><div>Mengxuan</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT_NAME Mg_Feslab</div><div>  PRINT_LEVEL MEDIUM</div><div>  RUN_TYPE GEO_OPT</div><div>&END GLOBAL</div><div>&MOTION</div><div>  &GEO_OPT</div><div>  &END GEO_OPT</div><div>&END MOTION</div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  STRESS_TENSOR ANALYTICAL</div><div>  &PRINT</div><div>    &FORCES On</div><div>    &END FORCES</div><div>  &END PRINT</div><div>  &SUBSYS</div><div>    &COORD</div><div>Mg    0.000001    1.860046   11.675848</div><div>Mg    0.000001    1.860046   16.865114</div><div>Mg    0.000001    1.860046   22.054379</div><div>Mg    0.000001    1.860046   27.243645</div><div>Fe    0.000001    1.860046   32.432911</div><div>Fe    1.610850    0.930024    9.081215</div><div>Mg    1.610850    0.930024   14.270481</div><div>Mg    1.610850    0.930024   19.459747</div><div>Mg    1.610850    0.930024   24.649012</div><div>Mg    1.610850    0.930024   29.838278</div><div>    &END COORD</div><div>    &CELL</div><div>      A             3.2216992378         0.0000000000         0.0000000000</div><div>      B  -1.6108503481         2.7900737882         0.0000000000</div><div>      C  0.0000000000         0.0000000000        41.5141258240</div><div>    &END CELL</div><div>    &KIND Mg</div><div>      POTENTIAL GTH-PBE-q2</div><div>      ELEMENT Mg</div><div>      BASIS_SET spd-1z-s-2z-rc7pt5-r06pt0-plus-diffuse-s-8pt0</div><div>    &END KIND</div><div> &KIND Fe</div><div>      POTENTIAL GTH-PBE-q16</div><div>      ELEMENT Fe</div><div>      BASIS_SET TZV2P-MOLOPT-SR-GTH-q16</div><div>      MAGNETIZATION 2.5</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>  &DFT</div><div>    POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div>    BASIS_SET_FILE_NAME PLATO_BASIS</div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT_UCL</div><div>    UKS .TRUE.</div><div>    MULTIPLICITY 5</div><div>    RELAX_MULTIPLICITY 0.05</div><div><br></div><div>    &MGRID</div><div>      CUTOFF [eV] 10000</div><div>      REL_CUTOFF [eV] 700</div><div>      NGRIDS 4</div><div>    &END MGRID</div><div>    &KPOINTS</div><div>      SCHEME MONKHORST-PACK 10 10 1</div><div>    &END KPOINTS</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>      ADDED_MOS 8</div><div>      MAX_SCF 300</div><div>      EPS_SCF 1.0E-7</div><div>      &MIXING T</div><div>        METHOD BROYDEN_MIXING</div><div>        ALPHA 0.1</div><div>        NBUFFER 8</div><div>      &END MIXING</div><div>      &DIAGONALIZATION ON</div><div>        ALGORITHM Standard</div><div>      &END DIAGONALIZATION</div><div>      &SMEAR ON</div><div>        METHOD FERMI_DIRAC</div><div>        ELECTRONIC_TEMPERATURE [K] 300</div><div>      &END SMEAR</div><div>    &END SCF</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-10</div><div>    &END QS</div><div>  &END DFT</div><div>&END FORCE_EVAL</div></div><div><br></div>