Computing polarizabilities without a Hamiltonian matrix

► We calculate polarizabilities and hyper-polarizabilities without computing or storing a Hamiltonian matrix. ► Reduced memory cost. ► Matrix-vector products evaluated by exploiting the structure of the basis.