In our task group approach, the short-and long-range MM atoms are divided among the groups, and the grid planes are distributed between the MPI processes of each group (see illustration in Fig. ). Thus, each MPI process receives a subset of short-and long-range MM atoms and a subset of grid planes, and computes the corresponding quantities of eqs 3-7.