Multiwfn forum / Question about spin contamination in Gaussian16

Re: Question about spin contamination in Gaussian16

dummy@example.com (sobereva) — Sun, 30 Jul 2023 14:21:16 +0000

The larger the deviation of printed to ideal , the larger the spin contamination. However, there is no widely accepted definition of percentage spin contamination.

Re: Question about spin contamination in Gaussian16

dummy@example.com (rikaaardoss) — Sun, 30 Jul 2023 13:50:15 +0000

Thanks for the clarification, I just have one last question: The values "0.6077" and "0.2087" what exactly do they refer to? What I thought was that 0.2087 was the spin contamination obtained after the annihilation process ( being the ideal value 0 for a singlet system, i.e. it is 20.87% contaminated).

Re: Question about spin contamination in Gaussian16

dummy@example.com (sobereva) — Sun, 30 Jul 2023 04:59:44 +0000

For ideal biradical singlet system, an alpha and a beta electron occupy an alpha spin orbital and a beta spin orbital, respectively, and the overlap between the two orbitals is negligible. While other occupied MOs are basically paired. In this case, according to the equation of evaluating expectation of S^2 operator in Section 3.100.5 of Multiwfn manual, will be 1.0.

Re: Question about spin contamination in Gaussian16

dummy@example.com (rikaaardoss) — Sat, 29 Jul 2023 19:36:16 +0000

Thanks for the answer, but it is not very clear to me. Why does a biradical singlet have a value of S = 1? Regarding the example I mentioned previously, the results show a high spin contamination?

Re: Question about spin contamination in Gaussian16

dummy@example.com (sobereva) — Sat, 29 Jul 2023 04:23:53 +0000

Ideal spin multiplicity of singlet is 0, therefore spin contamination currently is 0.9770. For an ideal singlet biradical, it is 1.0.

S is solved according to: S(S+1)=

Question about spin contamination in Gaussian16

dummy@example.com (rikaaardoss) — Fri, 28 Jul 2023 20:27:22 +0000

I have a small question regarding the assignment of spin contamination in a biradical (singlet) system. My output shows the following:

= 0.0000 = 0.0000 = 0.0000 = 0.9770 S= 0.6077
= 0.00000000000
Annihilation of the first spin contaminant:
S**2 before annihilation 0.9770, after 0.2087

As far as I understand, in this case, the spin contaminant taken is 0.2087 and the value 0.9770 represents the multiplicity of my compound, my question is S=0.6077 exactly what does it mean?

Thank you in advance