Add computation of <S^2> for Wavefunction types

[wavefunction91]

<S^2> is a common diagnostic for Wavefunction types, particularly single reference wavefunctions. It would be useful to compute this in the API, including fall throughs for exact / trivial cases (e.g. RHF).