We adopt the two-point correlation function (2PCF) as a statistical tool to quantify the spatial clustering of halo stars for galaxy systems spanning a wide range in host halo virial mass ( 11.25<log10M200c/M⊙<15 ) and redshifts (0 < z < 1.5) from the IllustrisTNG simulations. Consistent with a previous study, we identify clear correlations between the strength of the 2PCF signals and galaxy formation redshifts, but over a much wider mass range. We find that such correlations are slightly stronger at higher redshifts and get weakened with the increase of host halo mass. We demonstrate that the spatial clustering of halo stars is affected by two factors: (1) the clustering gets gradually weakened as time passes (phase mixing); (2) newly accreted stars at more recent times would increase the clustering. For more massive galaxy systems, they assemble late, and the newly accreted stars would increase the clustering. The late assembly of massive systems may also help to explain the weaker correlations between the 2PCF signals and the galaxy formation redshifts in massive halos, as their 2PCFs are affected more by recently accreted stars, while formation redshift characterizes mass accretion on a much longer timescale. We find that the orbits of satellite galaxies in more massive halos maintain larger radial anisotropy, reflecting the more active accretion state of their hosts while also contributing to their stronger mass loss rates.
Read full abstract