Localized-orbital computation of linear and nonlinear susceptibilities

Andrea Dal Corso, Francesco Mauri

We present a method to compute high-order derivatives of the total energy which can be used in the framework of density functional theory. We provide a proof of the $2n+1$ theorem…