“…Furthermore, the cluster backsplash galaxies [19,20] and shock radius [21] have also been well addressed. The advanced baryon models in hydrodynamic simulations allow us to perform a detailed investigation on the cluster properties, such as profiles [22,23], substructure and its baryonic content [24][25][26][27], dynamical state and morphology [28,29], cluster (non-)thermalization [30,31], the fundamental plane [32], and the cluster mass biases [33, for hydrostatic-equilibrium assumption], [34, for sigma-mass relation] and [35, for weaklensing]. Lastly, comparing to the void/field region runs in this project allows us to study the effect of environment [36]; including the self-interacting dark matter run allow us to constrain the dark matter cross-section [37]; It further help us to examine the chameleon gravity [38].…”