representing the reciprocal space vectors. More explicit expressions for E ovlp and E self and their derivatives are summarized in the supplementary information. The number of plane waves needed for an all-electron description would increase with the square of the nuclear charge and is therefore not feasible in routine applications. In order to avoid this bottleneck the external potential contributions of h uses pseudo potential for core electron contributions