<div>For anyone who reads this thread! It seems that it is possible to calculate the spin contamination value manually. One needs to include the eigenvectors into the CP2K output by adding the PRINT statement to the DFT block and print out the eigenvectors (https://manual.cp2k.org/cp2k-6_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/MO.html#EIGENVECTORS). I also specified the multiplicity of the calculation explicitly in the input.<br></div><div><br></div><div>The equation for the calculation of spin contamination is described here: https://en.wikipedia.org/wiki/Spin_contamination . And the eigenvectors are used to calculate the last term in the equation - I believe it is the sum of the overlap between the occupied alpha and beta orbitals. The overlap can be found by multiplying each alpha eigenvector with the corresponding beta eigenvector and then adding the results together.</div><div><br></div><div>I added an example of how to do it as an attachment.</div><div><br></div><div>Best Regards,</div><div>Mario<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, 3 December 2020 at 16:38:38 UTC ma...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>It seems that Juerg was right, it is not implemented for semi-empirical methods. The information about spin contamination became available right after I changed AM1 to the XC functional PBE.</div><div><br></div><div>Thank you both for your input. Would you think this is something the developers would be interested in adding in?</div><div><br></div><div>Best Regards,</div><div>Mario<br></div><div><br></div><div><br></div></blockquote></div>