Real space pseudopotentials have a number of advantages in solving for the electronic structure of materials. These advantages include ease of implementation, implementation on highly parallel systems, and great flexibility for describing partially periodic systems. One limitation of this approach, shared by other electronic structure methods, is the slow convergence of interatomic forces when compared to total energies. For real space methods, this requires a fine grid to converge a solution of the Kohn–Sham problem, which is accompanied by concurrent increase in memory and additional matrix-vector multiplications. Here we introduce a method to expedite the computation of interatomic forces by employing a high order integration technique. We demonstrate the usefulness of this technique by calculating accurate bond lengths and vibrational frequencies for molecules and nanocrystals without using fine real space grids.