<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear Bálint Aradi<div><br></div><div>Hi, I'm trying to reproduce the magnetic moments of iron clusters calculated by Köhler in his paper (C. Köhler et al., Chem. Phys., vol. 309, pp. 23, 2005) but I noticed that the values of spin constants and Hubbard Us used by him correspond to a different electronic configuration than the values used by default in DFTB+. When I use the parameters given in the manual which I think correspond to the configuration 3d74s1 I get very different magnetic moments than Köhler. When use spin and Hubbard constants for 3d64s2 (which I think is what Köhler used) my calculated magnetic moments are then in good agreement with his results except for a few cluster sizes.</div><div><br></div><div>In order to resolve this issues I have some questions:</div><div>¿Why are only spin and Hubbards U's constants for 3d74s1 configuration used in DFTB+, and not for 3d64s2 configuration?<br></div><div>In Köhler thesis there are Hubbards U's and spin constants for LDA and GGA, but I'm not sure if he is using LDA or GGA parameters for his magnetic moment calculations. I understand that Slater-Köster files for Iron in the trans3d set are for GGA, but then, do Slater-Köster files for LDA also exist? In that the case where can I find the LDA parameters?</div><div><br></div><div>As an example, for Fe15 the energies that I get  with different magnetic moments after relaxing the cluster (starting from the same symmetry than Köhler), have a minimum at 46 unpaired electrons, but according to the magnetic moment reported by him the minimum should be at 48.</div><div>The following are my DFTB+ energies:</div><div><br></div><div># Unpaired e / Energy (eV)</div><div><div>40<span style="white-space:pre-wrap">        </span>-593.2981</div><div>42<span style="white-space:pre-wrap">      </span>-593.6837</div><div>44<span style="white-space:pre-wrap">      </span>-594.0011</div><div>46<span style="white-space:pre-wrap">      </span>-594.0268</div><div>48<span style="white-space:pre-wrap">      </span>-593.9829</div><div>50<span style="white-space:pre-wrap">      </span>-592.7165</div></div><div><br></div><div>My hamiltonian is as follows, maybe there is some option that I missed to fix?</div><div><div><br></div><div>Hamiltonian = DFTB {</div><div>  SCC = Yes</div><div>  MaxSCCIterations = 50</div><div>  Eigensolver = RelativelyRobust {}</div><div>  SlaterKosterFiles = Type2FileNames {</div><div>    Prefix = "./"</div><div>    Separator = "-"</div><div>    Suffix = ".skf"</div><div>    LowerCaseTypeName = No</div><div>  }</div><div>  OrbitalResolvedSCC = Yes</div><div>  CustomisedHubbards {</div><div>#        U4s    U4p    U3d</div><div>    Fe = 0.237  0.214  0.447</div><div>  }</div><div>  MaxAngularMomentum {</div><div>    Fe = "d"</div><div>  }</div><div>  SpinPolarisation = Colinear {</div><div>    UnpairedElectrons = 6</div><div>    RelaxTotalSpin = Yes</div><div>  }</div><div>  SpinConstants = {</div><div>    Fe = {</div><div>#      Wss     Wsp     Wsd     Wps     Wpp     Wpd     Wds     Wdp     Wdd</div><div>      -0.015  -0.012  -0.003  -0.012  -0.019  -0.001  -0.003  -0.001  -0.017</div><div>    }</div><div>#    ShellResolvedSpin = Yes</div><div>  }</div><div>  Mixer = Broyden {}</div><div>}</div></div><div><br></div><div>Thank you in advance,</div><div><br></div><div>José Manuel</div></div></div></div></div></div></div></div>