Just want to check that I have got understanding correct.
To compute molecular polarizability need to use method such as Hartree fock(SCF), moller plesset, configuration interaction, coupled cluster etc. to get the wavefunctions and energy then just use sum over states formula.