Even though n Z is not including core compensation charges, the KS matrix contributions of Eq. ( ) are based on the total ground-state density n tot , thus leading to the mixed core compensation gradient terms of Eq. ( ). The right-hand side of the CPKS equations, R, is determined by the TDA eigenvalue equations and included GAPW expressions therefore do not include core compensation charges.