An algorithm for massively parallel computers is developed for energy calculations of second-order Møller–Plesset (MP2) perturbation theory with numerical quadratures. Message Passing Interface (MPI) and Open Multi-Processing (OpenMP) technologies are utilized for inter-node and intra-node parallelization, respectively. Computational tasks and intermediates are distributed across nodes by dividing quadrature points, and the distributed data are stored in memory. Benchmark calculations were performed on 256–8,192 CPU cores, and we observed the speed-ups 4,534–6,266 for 8,192 cores. A large calculation for fullerene (C60) with aug-cc-pCVTZ (3,540 basis functions) was completed in ca. 4.8 h on 8,192 cores without invoking molecular symmetry.