<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>Hello,</div>

<div>thank you for suggesting ASE, looks like something I can use.</div>

<div>And thanks for pointing out that (embarrassing) error. When I correct for that I get dE = -17 kJ/mol so now the TS is lower in energy, which is non-physical in this case. </div>

<div>I would go ahead and try other TS geometries to test. Do you think the input specification is sufficient for that or do I need shell resolved spin, specify initial spins,  ... ?</div>

<div> </div>

<div>
<div>Thanks and best regards,</div>

<div>Martin</div>
</div>

<div> 
<div> 
<div name="quote" style="margin:10px 5px 5px 10px; padding: 10px 0 10px 10px; border-left:2px solid #C3D9E5; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div style="margin:0 0 10px 0;"><b>Gesendet:</b> Dienstag, 03. Dezember 2019 um 16:33 Uhr<br/>
<b>Von:</b> "Mat Toliday via DFTB-Plus-User" <dftb-plus-user@mailman.zfn.uni-bremen.de><br/>
<b>An:</b> "User list for DFTB+ related questions" <dftb-plus-user@mailman.zfn.uni-bremen.de><br/>
<b>Cc:</b> "Mat Toliday" <djlets2004@yahoo.co.uk><br/>
<b>Betreff:</b> Re: [DFTB-Plus-User] Spin polarization and other</div>

<div name="quoted-content">
<div class="ydp91e973acyahoo-style-wrap" style="font-family: Helvetica Neue , Helvetica , Arial , sans-serif;font-size: 13.0px;">
<div> </div>

<div>Hello,</div>

<div>I don't think you can use DFTB+ inside of Gaussian but it claims to be capable of DFTB calculations. I spent a  couple of hours trying to modify the Slater-Koster files and then lost interest...</div>

<div>I don't think DFTB+ has the geometry optimisers needed for TS searching but these can be found using another package alongside DFTB+, namely the Atomic simulation Environment which has lots of useful geometry/optimiser tools and will interface with DFTB+.</div>

<div>As for your TS simulations:</div>

<div>I notice between your two input files there is a change in the total charge. This is most likely the the reason the dE value is too high.</div>

<div>Yours,</div>

<div>Mat</div>

<div> </div>
</div>

<div class="yahoo_quoted" id="yahoo_quoted_6322974839">
<div style="font-family: "Helvetica Neue" , Helvetica , Arial , sans-serif;font-size: 13.0px;color: rgb(38,40,42);">
<div>On Tuesday, 3 December 2019, 14:49:47 CET, Martin Richter <m-richter93@gmx.de> wrote:</div>

<div> </div>

<div> </div>

<div>
<div id="yiv9323306855">
<div>
<div style="font-family: Verdana;font-size: 12.0px;">
<div>Hello dear DFTB+ community,</div>

<div> </div>

<div>I have a few general questions and one specific problem. Id be very thankflul if you could help me. </div>

<div> </div>

<div>- I assume you cannot use DFTB+ via gaussian? I have seen options for DFTB and DFTBA though. </div>

<div>- I saw a question on 'spin polarization calculation' in the achrive, but no answer. Does this simply mean no one responded yet or did I miss something?</div>

<div>- Are there anywhere some more examples? I went through the tutorials of the homepage, but cannot find more useful examples. </div>

<div>- Can you perform transition state optimization ?</div>

<div> </div>

<div>And now my question regarding spin polarization. Im fairly new to QM chemistry so pardon any harsh mistakes. I want to study a hydrogen atom transfer reaction. To begin I just looked at CH4 + CH3* --> CH3* +CH4 (star * means radical). For the energy difference between precomplex and transition state (TS) I get from gaussian calculation  dE = 58 kJ/mol and literature tells me  dE = 39 kJ/mol. With my attempt in DFTB+ I get dE = 506 kJ/mol, which is way too much.</div>

<div> </div>

<div>Im wondering if I just made some mistakes or if maybe DFTB+ is not suited for this.</div>

<div> </div>

<div>Best regards and thanks,</div>

<div>Martin</div>

<div> </div>

<div>My input files look like this: </div>

<div>For TS:</div>

<div> </div>

<div>
<div>Geometry = GenFormat {<br/>
9 C<br/>
  C H<br/>
    1    1     -1.3492310000      0.0000000000     -0.0000040000                 <br/>
    2    2     -1.6387210000     -1.0390570000     -0.2023440000                 <br/>
    3    2     -1.6388130000      0.3443150000      1.0009990000                 <br/>
    4    2     -1.6386560000      0.6947620000     -0.7986980000                 <br/>
    5    1      1.3492320000      0.0000020000      0.0000070000                 <br/>
    6    2      1.6387010000     -0.6957490000      0.7978260000                 <br/>
    7    2      1.6387040000     -0.3431010000     -1.0014410000                 <br/>
    8    2      1.6387850000      1.0388000000      0.2035920000                 <br/>
    9    2     -0.0000060000      0.0000200000      0.0000490000                                <br/>
}</div>

<div>Driver = {}</div>

<div>Hamiltonian = DFTB {<br/>
  Scc = Yes<br/>
  SlaterKosterFiles {<br/>
    H-H = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/H-H.skf"<br/>
    C-H = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/C-H.skf"<br/>
    H-C = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/H-C.skf"<br/>
    C-C = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/C-C.skf"<br/>
  }<br/>
  MaxAngularMomentum {<br/>
    C = "p"    <br/>
    H = "s"<br/>
  }<br/>
  Charge = -1.0000000000000000<br/>
SpinPolarisation = Colinear {<br/>
UnpairedElectrons = 1.0<br/>
}<br/>
SpinConstants = {<br/>
C={-0.023}<br/>
H={-0.072}<br/>
}<br/>
 <br/>
}</div>

<div>Options {}</div>

<div>Analysis {<br/>
  CalculateForces = Yes<br/>
}</div>

<div>ParserOptions {<br/>
  ParserVersion = 7<br/>
}</div>

<div> </div>

<div> </div>

<div>and for the precomplex:</div>

<div> </div>

<div>
<div>Geometry = GenFormat {<br/>
9 C<br/>
  C H<br/>
 1           1            -4.673     1.29559  -0.1981 <br/>
 2           2            -4.68347   1.67679   0.82361 <br/>
 3           2            -4.47712   1.97238  -1.03062 <br/>
 4           2            -4.8825    0.24263  -0.38913 <br/>
 5           1            -0.56698   0.5267    0.17636 <br/>
 6           2            -0.15029   0.1488   -0.76893 <br/>
 7           2            -0.42145  -0.22437   0.96669 <br/>
 8           2            -1.64234   0.72488   0.05129 <br/>
 9           2            -0.05174   1.45716   0.45691                                          <br/>
}</div>

<div>Driver = {}</div>

<div>Hamiltonian = DFTB {<br/>
  Scc = Yes<br/>
  SlaterKosterFiles {<br/>
    H-H = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/H-H.skf"<br/>
    C-H = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/C-H.skf"<br/>
    H-C = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/H-C.skf"<br/>
    C-C = "/home/richtemn/Documents/DFTB+/recipes/slakos/mio-ext/C-C.skf"<br/>
  }<br/>
  MaxAngularMomentum {<br/>
    C = "p"    <br/>
    H = "s"<br/>
  }<br/>
  Charge = 0.0000000000000000<br/>
 SpinPolarisation = Colinear {<br/>
  UnpairedElectrons = 1.0<br/>
 }</div>

<div> SpinConstants = {<br/>
  C={-0.023}<br/>
  H={-0.072}<br/>
 }</div>

<div>}</div>

<div>Options {}</div>

<div>Analysis {<br/>
  CalculateForces = Yes<br/>
}</div>

<div>ParserOptions {<br/>
  ParserVersion = 7<br/>
}</div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br/>
DFTB-Plus-User mailing list<br/>
<a href="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de" onclick="parent.window.location.href='mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de'; return false;" target="_blank">DFTB-Plus-User@mailman.zfn.uni-bremen.de</a><br/>
<a href="https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user" target="_blank">https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user</a></div>
</div>
</div>
_______________________________________________ DFTB-Plus-User mailing list DFTB-Plus-User@mailman.zfn.uni-bremen.de <a href="https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user" target="_blank">https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user</a></div>
</div>
</div>
</div></div></body></html>