“…For many extended systems such a description is insufficient to account for optical excitations because the electron-hole attraction is not properly accounted for. However, dielectric properties, in particular collective plasmon excitations, are generally accurately reproduced by this approach 7,8 , and quantitative agreement with electron energy loss experiments have been reported for bulk metals 9,10 , surfaces 11,12 , graphene-based systems 13,14 , semiconductors 15,16 and even supercondutors 17 . Furthermore, the accurate evaluation of the density response function at the RPA or ALDA level is a prerequisite for implementation of most post-DFT schemes, such as RPA correlation energy 18 , exact-exchange optimized-effective-potential methods 19 , the GW approximation for quasi-particle excitations 20,21 , and the Bethe-Salpeter equation 21,22 for optical excitations.…”