Hello again. I would like some advice please. I am successfully using Multiwfn to calculate localized orbitals and energies. However I am not sure if I am using the correct process for molecules that have an unpaired electron (unrestricted wavefunction) e.g. NO2. I have done localized energy calculations in such cases by combining the A Fock matrix with the B Fock matrix. Is this correct? The results show the A and B energies to be similar, but not the same. In some cases the A orbital is different from the corresponding B orbital. Can you give me any advice on the interpretation of the (different) A/B results? Thanks. Regards Donald.
This is very normal situation. For open-shell case, spin polarization effect makes energies of alpha orbitals different to that of beta orbitals.
]]>Please pay attention to this statement in Section 3.22 of Multiwfn manual:
The matrix elements should be provided in lower-triangular form, namely in this sequence: F(1,1) F(2,1) F(2,2) F(3,1) F(3,2) F(3,3) ... F(nbasis,nbasis), where nbasis is the total number of basis functions, the format is free.
Obviously, you don't need to carry out any mathematical transformation on the origin Fock matrix, you just need to ignore all F(i,j) elements whose j>i.
Best regards,
Tian
]]>