where the sums run over the N fragments in the environment and S n sites of a given n'th fragment, q s is the charge at site s, µ α s is the α Cartesian component of the dipole at site s, Θ αβ s is the αβ Cartesian component of the quadrupole at site s, and