<html><head></head><body><div class="ydp91e973acyahoo-style-wrap" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 13px;"><div></div>
        <div dir="ltr" data-setdir="false">Hello,</div><div dir="ltr" data-setdir="false">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 dir="ltr" data-setdir="false">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 dir="ltr" data-setdir="false">As for your TS simulations:<br></div><div dir="ltr" data-setdir="false">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 dir="ltr" data-setdir="false">Yours,</div><div dir="ltr" data-setdir="false">Mat<br></div><div><br></div>
        
        </div><div id="yahoo_quoted_6322974839" class="yahoo_quoted">
            <div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
                
                <div>
                    On Tuesday, 3 December 2019, 14:49:47 CET, Martin Richter <m-richter93@gmx.de> wrote:
                </div>
                <div><br></div>
                <div><br></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 ymailto="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de" href="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de">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><br></div>
            </div>
        </div></body></html>