We present a new approach for numerical solutions of ab initio quantum chemistry systems. The main idea of the approach, which we call canonical diagonalization, is to diagonalize directly the second-quantized Hamiltonian by a sequence of numerical canonical transformations.