This talk will discuss multigrid methods for solving the nonlinear equations occuring in quantum chemistry. High order finite difference representations are used for the Laplacian, and all particles are represented discretely on a grid in real space. The equations are solved on a heirarchy of levels. The eigenvalue solver is coupled with a Poisson solver to handle the self consistency to locate the ground state of the electrons. A new method for high order conservative composite domain multigrid calculations will be discussed which will allow for increased resolution around the atomic nuclei. This is a generalization of the previous method of Bai and Brandt.