<html><head></head><body><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></body></html>