The GaMD implementation is such that only a small computational and communication (in parallel) overhead is added compared to cMD. The GaMD-dih and GaMD-dualwat have been considered on the STMV system (1 066 624 atoms) with the AMOEBA PFF and the 10 fs outer time-step HMR BAOAB-RESPA1 multiple-time-step integrator. 17 V100 GPUs from the national Jean Zay supercalculator have been used for all the benchmark computations. Similar scalability studies have been performed on the Jean Zay multi-CPUs (Figure ). The AMOEBA GPU simulations were performed on a single node as the multi-node extension of the AMOEBA PFF within the Tinker-HP package is still under development.On 1 and 2 GPUs (Figure ), the GaMD data communications are negligible, 1%. On 4 GPUs, the communications are increasing and the performance decreases by 7%. Overall, the use of GaMD only slightly alters the performance. This high scalability opens the door to simulate at a high-accuracy large complex biomolecular systems with PFFs.