Dear Matt and other CP2K gurus,<div><br></div><div><div>I hope you are well! We are developing a periodic embedding method in ChemShell to perform QM/MM MD simulations, typically on adsorbate/metal systems, using CP2K.</div><div><br></div><div>The attached O2/Pd calculation (PBE) converges extremely slowly. And it doesn't converge unless the density mixing ratio (&MIXING ALPHA) is only 2%-5%. Any higher ratio results in strong SCF fluctuation and it would never converge (I tried 2000 cycles). Does this imply the initial guess is too far away from a stationary point?</div><div><br></div><div>Various other options I tried but didn't seem to help:<br></div><ul><li>smearing</li><li>level shift</li><li>initial guess</li><li>other mixing methods</li><li>multi-grid cutoff and real cutoff</li><li>smooth XC grid methods</li></ul><div>I hope the supercell parameters are correct (because CP2K doesn't accept redundant atoms in the cell). Any idea how to speed up the SCF convergence to make QM/MM MD simulation more useful, please?</div><div><br></div><div>Thank you very much.<br></div></div><div><br></div><div>Kind regards,</div><div><br></div><div>You</div>