Quantifying the Impact of Incompleteness on Identifying and Interpreting Galaxy Protocluster Populations with the TNG-Cluster Simulation

Abstract: We use the TNG-Cluster simulation to examine how stellar mass and star formation rate (SFR) incompleteness affect the identification of density peaks within galaxy protoclusters at different redshifts. We focus on 352 protoclusters, defined as the set of galaxies that will reside within the virialized region of a $z=0$ cluster halo with total mass $\sim10^{14.3-15.5}~{\rm M}{\odot}$, and consider only galaxies with ${M}{\star}> 10^{8.5}~{\rm M}{\odot}$ (our baseline) at any redshift. We find that ${M}{\star}$-limited (${M}{\star} > 10^{9.5}~{\rm M}{\odot}$) and SFR-limited ($\rm{SFR} > 10~{\rm M}_{\odot} \mathrm{yr}^{-1}$) subpopulations only recover the true highest galaxy density peak in $\lesssim40\%$ of cases within an accuracy of $1.0$ pMpc ($\sim 2-2.5$ arcmin) at $z > 2$. We find that at $z>2$ the highest galaxy density peaks do not generally coincide with the highest dark matter or stellar mass density peaks, and this separation is larger than 0.5 pMpc in $\sim80\%$ of cases. Consequently, the region surrounding the true highest galaxy density peaks are not generally sites of enhanced star formation or accelerated mass growth relative to the remainder of the protocluster. Lastly, we find that the $\sim4$ arcmin apertures typically used to define spectroscopically-confirmed protoclusters are generally much smaller than the $8$ arcmin apertures needed to study the progenitors of the most massive galaxy clusters at $z > 2$.