<div dir="ltr"><br><div><span style="font-size:13px">Juerg Hutter,</span><br></div><div><span style="font-size:13px"><br></span></div><div><span style="font-size:13px">thank you very much for your answer.</span></div><div>It helped a lot. We still did not solve completely the problem. But, after we get the solution, we will comment on the forum about it.</div><div><span style="font-size:13px"><br></span></div><div>All the best,</div><div><br></div><div>Maicon P. Lourenço</div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-02-24 6:52 GMT-03:00  <span dir="ltr"><<a href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<br>
<br>
CP2K reads the second line of the SKF files in the following format:<br>
<br>
d,p,s (orb energy), energy, d,p,s (Hubbard), d,p,s, (occupation)<br>
<br>
the rest of the line is ignored. The energy value is not used.<br>
The two examples have both no d orbitals and therefore we can also<br>
ignore the d Hubbard value.<br>
<br>
The mio style files should work also in CP2K, but<br>
make sure that you only use the same basis (sp vs. spd),<br>
that you only use a single Hubbard value (no multiple values allowed in CP2K)<br>
that you are not using additional options like DFTB3.<br>
<br>
Another possible source of errors is the definition of the repulsion potential.<br>
CP2K defines the potential from the SKF file, while other programs do this<br>
through the input.<br>
Check for the value of "Repulsive potential energy:"<br>
<br>
regards<br>
<br>
Juerg Hutter<br>
<br>
--------------------------------------------------------------<br>
Juerg Hutter                         Phone : <a href="tel:%2B%2B41%2044%20635%204491" value="+41446354491">++41 44 635 4491</a><br>
Institut für Chemie C                FAX   : <a href="tel:%2B%2B41%2044%20635%206838" value="+41446356838">++41 44 635 6838</a><br>
Universität Zürich                   E-mail: <a href="mailto:hut...@chem.uzh.ch">hut...@chem.uzh.ch</a><br>
Winterthurerstrasse 190<br>
CH-8057 Zürich, Switzerland<br>
---------------------------------------------------------------<br>
<br>
-----<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a> wrote: -----To: cp2k <<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>><br>
From: Maicon Lourenco<br>
Sent by: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a><br>
Date: 02/23/2016 06:01PM<br>
Subject: [CP2K:7451] DFTB SKF format in CP2K<br>
<span class=""><br>
<br>
Dear all,<br>
Since we are trying to use CP2K to perform SCC-DFTB simulations, I've got some doubts about the second line of the SKF file as it is used in CP2K. For instance, we have the CP2K format and mio format for the pair O-O:<br>
<br>
CP2K style:<br>
<br>
0.02, 500,2<br>
0.0 -0.33213167 -0.87883246, -0.05414, 0.0  0.5222   0.5564 0.0 4.0 2.0 *T<br>
0.467495  0.523305*<br>
16.01,  19*0.0,<br>
<br>
DFTB (mio) style:<br>
<br>
0.02, 500,2<br>
0.0 -0.33213167 -0.87883246, -0.05414, 0.467495  0.523305   0.4954 0.0 4.0<br>
2.0<br>
16.01,  19*1.0,<br>
20*1.0,<br>
<br>
As far as we know, the meaning of the second line of the SKF is:<br>
d, p, s (valence atomic orbital energies), "spin polarization error", d, p, s (Hubbard parameter), d, p, s (valence orbital occupation).<br>
However, we do not get the order of the values in the second line of the CP2K O-O parameter and mainly the meaning of the "*T 0.467495  0.523305*".<br>
<br>
Another thing is, why is the "0.0" after the 0.05414 since it is not there in the mio format?<br>
<br>
As we need to use the mio format, how can I do that to use it properly in CP2K? We have observed some considerable total energy differences when the "T ..." is not present in the second line. How can I set that correctly?<br>
<br>
Thank you very much!<br>
<br>
Maicon Pierre Lourenço<br>
<br>
<br>
<br>
<br>
--<br>
<br>
</span>You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
<div class="HOEnZb"><div class="h5"><br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="noreferrer" target="_blank">https://groups.google.com/group/cp2k</a>.<br>
<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="noreferrer" target="_blank">https://groups.google.com/d/optout</a>.<br>
<br>
--<br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/SmxNCPgp-kI/unsubscribe" rel="noreferrer" target="_blank">https://groups.google.com/d/topic/cp2k/SmxNCPgp-kI/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="noreferrer" target="_blank">https://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="noreferrer" target="_blank">https://groups.google.com/d/optout</a>.<br>
</div></div></blockquote></div><br></div>