<p>Dear all,<br></p><p>I'm trying DFTB-calculation of CP2K for molecules consisted of more than three heavy atoms (e.g. CO2).</p><p>I have one problem that the atomic-forces depend on the sequence of heavy atoms in the input-file.</p><p>For example, the following input coordinates-(A) and (B) showed different results at the calculated forces,</p><p>(A),,,{at angstrom-unit}</p><p>O  0.0  0.0  +1.2</p><p>C  0.0   0.0    0.0</p><p>O   0.0  0.0  -1.2</p><p>(B),,,</p><p>C  0.0  0.0   0.0</p><p>O  0.0  0.0  +1.2</p><p>O  0.0  0.0  -1.2</p><p>Could anyone tell me whether the this problem has been occurred or reported?</p><p>To be more precise, the above dependency was not resolved by the following approaches</p><p>{Periodicity, SCC-calculation and changing compile condition (e.g. gfortran to ifort within Linux_x86)}.</p><p>The only solution was to copy&paste the Slater-Koster file of "co.spl" to "oc.spl",</p><p>or to focus on hydrogen compounds (e.g. CH4, H2O, etc.).</p>