We also use periodic boundary conditions in all 3 fundamental axes, to prevent the dominance of surface effects while effectively simulating bulk phenomena. The simulations were run on two nodes of the Univetsité de Picardie Jules Verne MatriCS computational platform with 500 GB of RAM and 2 AMD EPYC 7513 32-Core Processors.