Implementation of spin polarisation

All data structures (density matrices, hamiltonians, etc.) need to include an extra dimension for spin. In addition, algorithms such as SPADE localisation and the associated routines used to construct the eigenvectors from the KS-eigenvalue problem need to be open shell.