in the auxiliary fitting basis functions { A }. Several strategies to determine the expansion coe cients can be employed, but in our case this is done by minimising the least squares error of the electric field, commonly referred as robust density fitting. Thereby, the coe cients are obtained by the minimization of the positive definite functional