Dear CP2K community,<div><br /></div><div>I am trying to optimize geometries for the reactant and product states of an electron transfer reaction, where the movement of electron is being tracked by constraining charges using CDFT. The coordinates and wavefunction were restarted from a previous DFT PBE-D3 geometry optimization. However,  I have noticed for all of the test cases that reaching optimization criteria is very difficult. The charges take forever to optimize at each geometry step and sometimes start to give a drastically wrong number of valence electrons. </div><div>Do you have any idea why this might be happening and if there is any way to speed up this process? I have tried these calculations with Becke/Hirshfeld partitioning schemes and BFGS/CG optimizers and ran into the same problem every time.</div><div><br /></div><div>I have attached the zipped folder containing all the necessary files to run the calculation for your reference (cp2k.psmp -i energy_run.inp). </div><div><br /></div><div>Thank you all!</div><div><br /></div><div>Regards,</div><div>Akanksha</div><div><br /></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/35771128-b809-4e0b-85af-4874f9c4f7f5n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/35771128-b809-4e0b-85af-4874f9c4f7f5n%40googlegroups.com</a>.<br />