<div dir="ltr"><div>Dear users of DFTB+<br></div><div><br></div><div>I just noticed that the SCC does not always converge to the same energy when optimizing iron clusters with a high number of unpaired electrons, depending on the initial charges or the initial number of unpaired electrons (if spin relaxation is allowed). I need to know if it is a normal behavior or if it is a bug in the code. I also noticed that spin relaxation only works when the initial number of unpaired electrons is slightly different than zero (for example for Fe6 it works with 0.1 initial unpaired electrons but not with 0, and for Fe11 it does not work even with 0.5 initial unpaired electrons) despite the fact that the most stable multiplicity is larger than zero. I would like to understand why it does not optimize the number of unpaired  when their initial value is too close to zero. Thanks in advance for your help.</div><div><br></div><div>This is an example of my input:</div><div><br></div><div>******************************************************************************************</div><div>Geometry = GenFormat {<br>6 C<br>Fe<br>     1   1  -1.5133154160E+00  -3.1163699090E+00   4.8395525190E+00<br>     2   1  -2.2869331200E+00  -1.1577302380E+00   7.7144216990E+00<br>     3   1  -5.9171879890E-01  -2.5667263510E+00   6.9211489160E+00<br>     4   1  -1.0719373520E+00  -7.3916752260E-01   5.5466727030E+00<br>     5   1  -3.2012161060E+00  -1.7091401130E+00   5.6341844250E+00<br>     6   1  -2.7174110570E+00  -3.5378755450E+00   7.0109896540E+00<br>}<br><br>Driver = ConjugateGradient {<br>  MovedAtoms = 1:-1<br>  ConvergentForcesOnly = No<br>  AppendGeometries = Yes<br>  MaxSteps = 300<br>}<br><br>Hamiltonian = DFTB {<br>  SCC = Yes<br>  MaxSCCIterations = 50<br>  Eigensolver = RelativelyRobust {}<br>  SlaterKosterFiles = Type2FileNames {<br>    Prefix = "./"<br>    Separator = "-"<br>    Suffix = ".skf"<br>    LowerCaseTypeName = No<br>  }<br>  MaxAngularMomentum {<br>    Fe = "d"<br>  }<br>  SpinPolarisation = Colinear {<br>    UnpairedElectrons = 0<br>    RelaxTotalSpin = Yes<br>  }<br>  SpinConstants = {<br>    Fe = {<br>      #  ss      sp      sd     ps       pp      pd     ds      dp      dd<br>      -0.016  -0.012  -0.003  -0.012  -0.029  -0.001  -0.003  -0.001  -0.015<br>    }<br>    ShellResolvedSpin = Yes<br>  }<br>  Mixer = Broyden {}<br>  Filling = Fermi {<br>    Temperature [K] = 1.0<br>  }<br>}<br><br>ParserOptions {<br>  ParserVersion = 5<br>}<br></div><div>******************************************************************************************</div><div><br></div><div>José Manuel<br></div><br></div>