We investigate the formation and evolution of stellar halos (SHs) around bright central galaxies (BCGs), focusing on their scaling relations, colors, and metallicities across cosmic time, and we compare model predictions with ultra-deep imaging data. We used the semianalytic model , applied to merger trees from high-resolution dark matter simulations, including an updated treatment of intracluster light (ICL) formation. SHs are defined as the stellar component within a physically motivated transition radius, linked to the structural properties of the host halo. Predictions are compared with observations from the VST Early-type GAlaxy Survey (VEGAS) and Fornax Deep Survey (FDS). FEGA25 The SH mass correlates well with the BCG and ICL masses, with tighter scatter in the SH--ICL relation. The transition radius peaks at 30--40 kpc nearly independent of redshift in the model predictions, but can reach sim400 kpc in the most massive halos, after z=0.5. SHs and ICL show nearly identical color distributions at all epochs, both reddening toward z=0. At z=2, SHs and the ICL are sim0.4 dex more metal--poor than BCGs, but the gap shrinks to sim0.1 dex by the present time. Observed colors are consistent with model predictions, while observed metallicities are lower, suggesting a larger contribution from disrupted dwarfs. SHs emerge as transition regions between BCGs and the ICL, dynamically and chemically coupled to both. Their properties depend on halo concentration, ICL formation efficiency, and the progenitor mass spectrum. Upcoming wide--field photometric and spectroscopic surveys (e.g. LSST, WEAVE, 4MOST) will provide crucial tests by mapping structure, metallicity, and kinematics in large galaxy samples.
Contini et al. (Tue,) studied this question.