<div dir="ltr">OK, I've noticed that the repulsive energy<br>of the initial structure in the above example<br>is different in CP2K and DFTB (5.2533030486<br>and 5.3189364428 Ha, respectively), so I've<br>now looked at a simpler structure (linear <br>TiO2 chains) to see where that difference<br>is coming from.<br><br>I've appended a slide showing the results<br>for various chain structures, which seem to<br>indicate that certain pairs are missing from the<br>repulsive energy sum in CP2K. Is this a bug<br>in the code, or must this be handled through<br>the input somehow? I've also attached the <br>CP2K in- and output files for the (N=1, L=5) <br>case of the overview.<br><br>/Maxime<br><br></div>