add phi2deriv
Use Poisson's equation for the density
fix second derivatives of PowerSphericalPotentialwCutoff