where it should be understood that all matrices on the left hand side are evaluated at the reference coordinates R 0 . This result is analogous to what we would obtain if we were instead expanding about some reference density obtained from a Kohn-Sham calculation. In that case, the resulting Kohn-Sham Euler equation allows us to cancel out the electrostatic components with the functional derivative of the kinetic and XC components, less the constant