“…[10]. In fact, by constructing the localized perturbation (typically requiring costly calculations in supercells) as a series of independent monochromatic perturbations in the primitive unit cell, it improves significantly the computational efficiency, accuracy, user-friendliness, and automation [2,3], as also demonstrated by several recent applications [14,41,15,42,43,44,45,46,47,48,49]. Key to this successful implementation of the LR-cDFT is indeed the capability to express perturbation theory in reciprocal space as in the calculation of phonons using DFPT [50,51,52].…”