<div dir="ltr">Dear DFTB users,<div>When I am running DFTB minimization it is diverging SCC error is reaching to 10^29. Please find below the attached input file. Please help me in resolving the issue.</div><div><br><div class="gmail_extra"><div style="font-size:12.8px">Geometry = GenFormat {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><<< "input.gen" </div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">}</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Driver = ConjugateGradient {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> MovedAtoms = 1:-1</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> MaxSteps = 5000</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> MaxForceComponent = 1e-4</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> AppendGeometries = Yes</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">}</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Hamiltonian = DFTB {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  SCC = Yes</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  SCCTolerance = 1e-6</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  MaxSCCIterations = 5000</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  SlaterKosterFiles = Type2Filenames {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   Prefix = "/home/trans3d-0-1/"</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   Separator = "-"</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   Suffix = ".skf"</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   LowerCaseTypeName = No</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   MaxAngularMomentum {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    Fe = "d"</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">     C = "p"</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">     O = "p"</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">     H = "s"</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   OrbitalResolvedSCC = Yes</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   SpinPolarisation = Colinear {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   UnpairedElectrons = 0.0</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    RelaxTotalSpin = Yes</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    InitialSpins = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">     AtomSpin = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">      Atoms = 1:100</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">      SpinPerAtom = 4.0</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">     }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    AtomSpin = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">      Atoms = 101:119</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">      SpinPerAtom = 0.0</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">     }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> </div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   </div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   SpinConstants = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    Fe = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    -0.016 -0.012 -0.003 -0.012 -0.029 -0.001 -0.003 -0.001 -0.015</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">      }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    C = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   -0.031 -0.025 -0.025 -0.023 </div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   O = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   -0.035 -0.030 -0.030 -0.028</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  H = {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   -0.072</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">   }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  Filling = MethfesselPaxton {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    Temperature [Kelvin] = 300.0</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">    Order = 2</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  ReadInitialCharges = No</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  KPointsAndWeights = SuperCellFolding {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  1 0 0</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  0 1 0</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  0 0 1</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  0.0 0.0 0.0</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">  }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">}</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Options {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> WriteDetailedOut = Yes</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">}</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">ParserOptions {</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> ParserVersion = 4</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">}</div><div><div class="gmail_signature"><div dir="ltr">Regards,<div>Dharmendr Kumar</div><div><br></div></div></div></div>
<br><div class="gmail_quote">On Fri, Nov 17, 2017 at 6:40 PM, Dharmendra Kumar <span dir="ltr"><<a href="mailto:dkumar.rav@gmail.com" target="_blank">dkumar.rav@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><span style="color:rgb(136,136,136);font-size:12.8px">Dear Dr. Bálint,</span><br><div><font color="#888888"><span style="font-size:12.8px">Thank you for your immediate response.</span></font></div><div><font color="#888888"><span style="font-size:12.8px">I will incorporate the suggested change and let you know if error persists.</span></font></div><div><font color="#888888"><span style="font-size:12.8px"><br></span></font></div><div><font color="#888888"><span style="font-size:12.8px">Regards,</span></font></div><div><font color="#888888"><span style="font-size:12.8px">Dharmendra</span></font></div><div><font color="#888888"><span style="font-size:12.8px"> </span></font></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail-m_7882145858044714609gmail_signature"><div dir="ltr">Regards,<div>Dharmendr Kumar</div><div>TCS-Pune</div></div></div></div>
<br><div class="gmail_quote"><div><div class="gmail-h5">On Fri, Nov 17, 2017 at 6:29 PM, Bálint Aradi <span dir="ltr"><<a href="mailto:aradi@uni-bremen.de" target="_blank">aradi@uni-bremen.de</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><div class="gmail-h5">Dear Dharmendra,<br>
<br>
> <<a href="https://www.dftb.org/parameters/download/trans3d/" rel="noreferrer" target="_blank">https://www.dftb.org/paramete<wbr>rs/download/trans3d/</a>> paramete<wbr>rs but the<br>
<br>
As the description of the parameter sets states, that set was created to<br>
describe organic molecules with single 3d transition metals in it, so I<br>
would not suggest to use it in general for slabs.<br>
<br>
On the other hand, the Fe-Fe interaction in the trans3d set happens to<br>
be exactly the same as the one in the pbc-set (which was meant to be<br>
used with periodic systems), so there is no reason, why it should not<br>
work. However, make sure you set OrbitalResolvedScc = Yes in the DFTB-block.<br>
<br>
  Best regards,<br>
<br>
  Bálint<br>
<span class="gmail-m_7882145858044714609HOEnZb"><font color="#888888"><br>
--<br>
Dr. Bálint Aradi<br>
Bremen Center for Computational Materials Science, University of Bremen<br>
<a href="http://www.bccms.uni-bremen.de/cms/people/b-aradi/" rel="noreferrer" target="_blank">http://www.bccms.uni-bremen.de<wbr>/cms/people/b-aradi/</a><br>
<br>
<br>
</font></span><br></div></div>______________________________<wbr>_________________<br>
DFTB-Plus-User mailing list<br>
<a href="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de" target="_blank">DFTB-Plus-User@mailman.zfn.uni<wbr>-bremen.de</a><br>
<a href="https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user" rel="noreferrer" target="_blank">https://mailman.zfn.uni-bremen<wbr>.de/cgi-bin/mailman/listinfo/<wbr>dftb-plus-user</a><br></blockquote></div><br></div>
</blockquote></div><br></div></div></div>