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