Dear CP2K users and developers<br><br>I'm a relatively new CP2K user, and I'm trying to do MD simulations of a dimer of a small molecule in water with the GFN-xTB method. To prepare the simulations I first minimize and then equilibrate the boxes in a NPT ensemble (300 K, 1 bar) with the AMBER force field, followed by a short NVT run (also with AMBER FF). I then start the xTB simulation from the final point of the classical MD run, in the NVT ensemble (300 K).  However, I noticed that voids appear in the xTB simulation boxes, probably indicative of a wrong system density.<br><br>To reproduce the issue in a simpler system, I ran a NPT simulation with xTB (300 K, 1 bar) on a box of 126 water molecules in CP2K. The box was pre-equilibrated in a NPT run (also 300 K, 1 bar), with the AMBER FF TIP3P parameters, yielding a density of 1.023 g/cm3. This box was then subjected to a minimization and a short NVT equilibration with xTB, prior to the NPT equilibration. Indeed, I observed that the xTB NPT run compressed the box, to a final density of 1.481 g/cm3.<br><br>I'm attaching here the input for the xTB NPT simulation of water. The sections that are commented are settings which I also tried, without success, namely:<br>- Adding a NONBONDED potential to xTB (taken from the test file xTB/regtest-debug/h2o-atprop_nonbonded.inp)<br><div>- Changing the settings in the POISSON/EWALD section to non default values (up to RCUT = 10.0 Angstrom)</div><div><br></div>Note that I also ran the test xTB/regtest-debug/h2o-atprop_nonbonded.inp for 20000 steps and obtained an abnormally high density (1.722 g/cm3).<br>All runs were performed in CP2K v8.2, but I also noticed the same issues when running in version 7.1.<br><br>Could you please take a look at the input file and tell me if I'm missing something?<br><br>Best regards<br>Pedro Maximiano