Analytical Hartree–Fock calculation of the dynamical polarizabilities α, β, and γ of molecules in solution