In the era of exascale supercomputers, largeāscale, and longātime molecular dynamics (MD) calculations are expected to make breakthroughs in various fields of science and technology. Here, we propose a new algorithm to improve the parallelization performance of message passing interface (MPI)ācommunication in the MPIāparallelized fast multipole method (FMM) combined with MD calculations under threeādimensional periodic boundary conditions. Our approach enables a drastic reduction in the amount of communication data, including the atomic coordinates and multipole coefficients, both of which are required to calculate the electrostatic interaction by using the FMM. In communications of multipole coefficients, the reduction rate of communication data in the new algorithm relative to the amount of data in the conventional one increases as both the number of FMM levels and the number of MPI processes increase. The aforementioned rate increase could exceed 50% as the number of MPI processes becomes larger for very large systems. The proposed algorithm, named the minimumātransferred data (MTD) method, should enable largeāscale and longātime MD calculations to be calculated efficiently, under the condition of massive MPIāparallelization on exascale supercomputers.