<div dir="ltr"><div>Dear cp2k community,</div><div><br></div><div>I am trying to perform a SE MD simulation with 64 water-molecules with PBC using the SCP-NDDO method, as described in Murdachaew, Mundy, Schenter, Laino and Hutter, J. Phys. Chem. A 115, 6048 (2011). Unfortunately, I have not been able to reproduce the reported results yet. </div><div><br></div><div>I know that the usage of PBC in SE simulations needs a lot of insight into the algorithms. Hence I tried to optimize the parameters in the SE and POISSON/EWALD sections, but after a couple of MD steps the calculation ceases to convergence. </div><div><br></div><div>Therefore I tried using PERIODIC EWALD_GKS in the SE-section instead of PERIODIC EWALD. According to the manual and some testings this method disables the usage of multipoles (the manual rather says that PERIODIC EWALD enables multipoles) and is independent from the definitions inside the COULOMB and EXCHANGE methods. </div><div>However, the revised version works well till 600 steps, then the kinetic energy of the system starts rising and a few H-O bonds seem to break. </div><div><br></div><div>What is the difference between EWALD and EWALD_GKS?</div><div>Is an outer SCF loop of TYPE SCP necessary? </div><div><br></div><div>Many thanks in advance!</div><div>Bassirian</div></div>