These simulations are performed for a set of 324 Lagrangian regions of radius \(15 \ h^{-1}{\mathrm {Mpc}}\) , each centred on a galaxy cluster with \(M_{500} > 4.6\times 10^{14} \ h^{-1} {\mathrm {M}_{\odot }}\) . These clusters are the ones originally recognised as the most massive galaxy clusters (virial mass \(M > 8\times 10^{14} \ h^{-1}{\mathrm {M}_{\odot }}\) ) identified at \(z=0\)  in the MultiDarkPlanck2 box of the simulation MultiDarkPublically available at: https://www.cosmosim.org/ [1]}. This is a dark–matter (DM)–only simulation and the MultiDarkPlanck2 box is a periodically defined cube with size \(\sim 1.5 \ \mathrm {Gpc}\)  comoving length. This cube contains \(3840^3\)  DM particles, each of mass \(1.5\times 10^{9} \ h^{-1}{\mathrm {M}_{\odot }}\) . As suggested by its name, this simulation assumes a cosmology based on the best–fit model and cosmological parameter estimates from the Planck 2015 cosmological results [2]}.
