Computing Second-Order Functional Derivatives with Respect to the External Potential

Following the increasing interest in the higher-order (functional) derivatives of conceptual density functional theory, we developed and implemented a method for calculating second-order functional derivatives with respect to the external potential. Our method is theoretically exact but involves two numerical approximations: the functional derivatives are expanded in a basis set, and the values of the corresponding expansion coefficients are …