*NVT*, and *NPT* ensembles. For the numeric integration to solve the Newton's equations of motion, the program uses the rRESPA, a multiple time-step algorithm. The distance constraints between atoms are treated by the SHAKE, RATTLE, and ROLL algorithms. The program can handle all-atom force fields such as the CHARMM22 with CMAP, CHARMM36 with CMAP. In the near futre, the AMBER, and OPLSAA will be supported. Further, the program will also support free-energy calculation algorithms based on the thermodynamic integration method.

The program also equips several newly developed methods to execute highly parallelized molecular dynamics calculations. The program ensures excellent scalability by the use of algorithms that practically eliminate data copying for communications and arithmetic operations, algorithms with minimal communication latency, and a parallel bucket-relay communication algorithm for the upper-level multipole moments in the FMM. Moreover, the use of blocked arithmetic operations can avoid the need to reload data from memory to cache, ensuring very low cache-miss rates. A benchmark test on MODYLAS using 65,536 nodes of the K-computer showed that the overall calculation time per step including communications is 5 ms for a 10 million atom system. This means that the simulation of a 10 million atom system can progress by 35 ns per day. MODYLAS thus enables us to study large-scale real systems such as viruses, liposomes, protein aggregates, micelles, and polymers.