<div dir="ltr">Hello!<br><br>I've been trying out the DFTB module in CP2K on a simple <br>bulk TiO2 (rutile) structure. Strangely enough, it appears<br>impossible for the different optimizers to successfully<br>relax this structure. With BFGS, the atoms keep oscillating<br>along some kind of vibrational mode. With CG and LBFGS,<br>no such oscillations occur, but the forces cannot be brought<br>below circa 0.015 Hartree/Bohr.<br><br>I've attached a zip file with the input and output files for<br>the different optimizers. I've generated the Slater-Koster files<br>myself, and they work nicely with the DFTB+ code (without <br>any geometry optimization issues). I'm using a development <br>version of CP2K (6.0, git:8b033c3).<br><br>Could you help me understand what is going on? <br><br>Best,<br>Maxime<br></div>