Stellar halos of bright central galaxies II: Scaling relations, colors and metallicity evolution with redshift

Using the updated FEGA25 semi-analytic model, this study reveals that stellar halos around bright central galaxies act as dynamically and chemically coupled transition regions to the intracluster light, with their scaling relations, colors, and metallicity evolution across cosmic time showing strong agreement with observed colors but suggesting a larger contribution from disrupted dwarf galaxies than currently predicted.

The Gist

Imagine a massive, ancient city (a Bright Central Galaxy, or BCG) sitting in the middle of a vast, foggy countryside. This city has a dense, bustling downtown (the Bulge/Inner Galaxy), a sprawling suburban ring (the Stellar Halo), and a thin, misty haze stretching out for hundreds of miles into the surrounding fields (the Intracluster Light or ICL).

For a long time, astronomers struggled to tell where the "suburbs" ended and the "fog" began. This paper is like a detective story where the authors use a super-advanced computer simulation (a "time machine") to figure out how these suburbs formed, what they are made of, and how they change over billions of years.

Here is the story of their findings, broken down into simple concepts:

1. The "Fog" and the "Suburbs" are actually the same family

In the past, scientists thought the "suburbs" (Stellar Halos) and the "fog" (Intracluster Light) were two different things. But this paper argues they are actually two sides of the same coin.

• The Analogy: Imagine a family reunion. The "downtown" is the main house where the family lives. The "suburbs" are the cousins who moved to the edge of the property. The "fog" is the distant relatives living in the next town over.
• The Discovery: The authors found that the "suburbs" are just the inner part of the "fog." They are made of the exact same stuff (stars stripped from smaller galaxies) and they move together. You can't really draw a hard line between them; it's a smooth transition.

2. How big is the "Suburb"? (The Transition Radius)

The authors wanted to know: How far out does the "suburb" go before it becomes just "fog"?

• The Finding: In their simulation, this boundary (called the Transition Radius) is usually about 30–40 kilometers (in cosmic terms, that's huge!). However, in the biggest, most massive galaxies, this "suburb" can stretch out to 400 kilometers.
• The Catch: When they compared their computer model to real photos taken by telescopes (like VEGAS and FDS), the real "suburbs" looked smaller.
• Why the difference? It's like trying to measure a cloud. In the computer, they defined the edge based on physics (where stars stop being bound to the galaxy). In real life, astronomers define the edge based on how the light looks in a photo. The definitions are slightly different, so the measurements don't match perfectly, but the general trend is the same.

3. The "Makeup" of the Stars (Colors and Metallicity)

Astronomers look at the color of stars to guess their age and metallicity (how many heavy elements like iron they have). Think of "metallicity" like the seasoning in a soup.

• Old Soup (Low Metallicity): Made of simple ingredients (hydrogen/helium).
• Rich Soup (High Metallicity): Made with lots of spices (heavy elements) created by previous generations of stars.

The Evolution of Flavor:

• In the Past (High Redshift): The "downtown" (BCG) was like a rich, spicy soup. The "suburbs" and "fog" were like bland, watery soup. There was a big difference in flavor (about 0.4 "units" of difference).
• Today (Low Redshift): Over time, the "suburbs" and "fog" got richer. They absorbed stars from bigger, more "seasoned" galaxies that were torn apart. Now, the difference between the downtown and the suburbs is tiny (only 0.1 "units"). They have almost the same flavor!

The Twist:
When the authors looked at real data from the Fornax galaxy cluster, the "suburbs" there were still quite bland (low metallicity). This suggests that in the real world, these halos are made mostly of tiny, weak galaxies (dwarfs) that were torn apart, rather than the big, rich ones the computer model predicted. It's like the real city's suburbs were built by small, poor families, while the simulation assumed they were built by wealthy families.

4. The "Time Travel" Aspect

The paper tracks these galaxies from 10 billion years ago (when the universe was young) to today.

• The Trend: Everything gets "redder" (older) as time goes on. The stars stop forming new blue babies and just age.
• The Surprise: The "suburbs" and the "fog" age at the same rate and look almost identical in color. You cannot tell them apart just by looking at a photo; you need to know their "family history" (kinematics) to tell them apart.

The Big Picture Conclusion

This paper tells us that Stellar Halos are not isolated islands. They are the "inner ring" of the vast cosmic fog surrounding a galaxy.

• They are coupled: The galaxy, its halo, and the surrounding fog are all chemically and dynamically linked.
• They are history books: The stars in the halo are fossils. By studying them, we can see which smaller galaxies were eaten by the big one in the past.
• The Future: The authors say that upcoming giant telescopes (like LSST) and spectroscopic surveys (like WEAVE) will act like high-resolution microscopes. They will finally let us taste the "soup" of thousands of galaxies, helping us understand exactly which "ingredients" (small vs. large galaxies) built the halos we see today.

In short: The paper confirms that the "suburbs" of giant galaxies are just the inner edge of the cosmic fog, formed by the slow, messy process of galaxies eating each other over billions of years. The only way to truly understand them is to look at their chemical history, not just their shape.

Technical Summary

Here is a detailed technical summary of the paper "Stellar halos of bright central galaxies II: Scaling relations, colours and metallicity evolution with redshift" by Contini et al.

1. Problem Statement

Stellar halos (SHs) are diffuse, low-surface-brightness envelopes surrounding galaxies that preserve the fossil record of hierarchical assembly and accretion events. While previous studies have established correlations between halo properties and host mass, there remains a need to understand:

• The precise physical definition and scaling relations of SHs relative to Brightest Central Galaxies (BCGs) and Intracluster Light (ICL).
• The evolutionary history of SH properties (mass, size, color, metallicity) across cosmic time ($z=2$ to $z=0$).
• The distinction between SHs and the surrounding ICL, particularly whether they are chemically and dynamically distinct components or part of a continuous transition.
• The ability of semi-analytic models (SAMs) to reproduce ultra-deep observational data from surveys like VEGAS and Fornax Deep Survey (FDS).

2. Methodology

Model Framework:

• FEGA25: The authors utilize the semi-analytic model (SAM) FEGA25 (Formation and Evolution of GAlaxies), applied to merger trees extracted from high-resolution dark matter simulations (GADGET-4).
• Physics: The model includes updated prescriptions for baryonic physics, specifically AGN feedback (AGNeject1 option), supernova feedback, and a sophisticated star formation law.
• ICL and SH Definition:
  • ICL Formation: Modeled via three channels: stellar stripping (dominant), mergers (minor/major), and pre-processing.
  • Transition Radius ($R_{trans}$): The SH is defined as the stellar component within a physically motivated transition radius, $R_{trans}$. This radius is derived from the concentration of the ICL ($c_{ICL}$), which is linked to the dark matter halo concentration ($c_{200}$) via a scaling factor $\gamma$ (typically 1–3).
  • SH Mass: Stars within $R_{trans}$ are classified as SH. Stability criteria ensure SH mass does not exceed the BCG mass; excess is redistributed.
• Simulations: Three cosmological simulations (YS300, YS200, YS50HR) with varying volumes and resolutions were used to cover halo masses from Milky Way-sized systems to massive clusters, following Planck 2018 cosmology.

Observational Data:

• Surveys: Data from the VST Early-type GAlaxy Survey (VEGAS) and Fornax Deep Survey (FDS) provided photometric colors ($g-r$, $r-i$).
• Metallicity: Complemented by spectroscopic data from Fornax3D (F3D) and M3G surveys.
• Comparison Strategy: Since observational definitions of SH boundaries differ from the model's $R_{trans}$, the authors used an empirical calibration to map observational transition radii ($R_2$) to model predictions, treating the comparison statistically rather than one-to-one.

3. Key Contributions

1. Physical Definition of SHs: The paper proposes a rigorous, physics-based definition of SHs as a transition region between the bound BCG and the unbound ICL, determined by halo concentration rather than arbitrary photometric cuts.
2. Redshift Evolution: It provides a comprehensive timeline of SH evolution from $z=2$ to $z=0$, tracking mass, size, color, and metallicity simultaneously.
3. Model-Data Synthesis: It bridges the gap between theoretical SAMs and ultra-deep imaging, demonstrating that SHs and ICL are chemically and color-wise indistinguishable in the model, challenging the utility of broadband photometry alone for separating these components.
4. Scaling Relations: It establishes tight linear scaling relations between SH mass and both BCG and ICL masses, identifying halo concentration and ICL formation efficiency as primary drivers.

4. Key Results

Scaling Relations and Structure:

• Mass Correlations: SH mass correlates tightly with both BCG and ICL masses. The SH–ICL relation shows significantly less scatter than the SH–BCG relation, confirming the direct origin of SH stars from the ICL.
• Drivers: SH mass increases in less concentrated (typically more massive) halos and in systems with higher ICL formation efficiency.
• Transition Radius ($R_{trans}$): The model predicts $R_{trans}$ peaks at 30–40 kpc, largely independent of redshift. However, in the most massive halos, $R_{trans}$ can reach $\sim400$ kpc at $z=0$.
• Observational Offset: Model-predicted $R_{trans}$ values are systematically larger than those inferred from VEGAS observations. This is attributed to different definitions (physical vs. photometric decomposition) and sample selection effects.

Colors and Metallicities:

• Color Evolution: All components (BCG, SH, ICL) redden from $z=2$ to $z=0$. At $z=0$, SHs and ICL have nearly identical color distributions ($g-r$ and $r-i$), making them indistinguishable via broadband photometry. BCGs remain slightly redder on average.
• Metallicity Evolution:
  • High Redshift ($z=2$): SHs and ICL are $\sim0.4$ dex more metal-poor than BCGs.
  • Low Redshift ($z=0$): The metallicity gap shrinks to $\sim0.1$ dex as SHs/ICL become enriched through the accretion of more massive, chemically evolved satellites.
  • Observational Discrepancy: Observed metallicities (FDS) are systematically lower than model predictions. This suggests that observed halos in the local universe may have a larger contribution from disrupted low-mass dwarf galaxies compared to the central-galaxy dominated sample in the simulation.

5. Significance and Conclusions

• Nature of Stellar Halos: The study concludes that SHs are not isolated components but represent the inner manifestation of the ICL, dynamically and chemically coupled to the outer galaxy environment. They serve as a transition zone between the bound BCG and the diffuse ICL.
• Future Surveys: The findings imply that broadband colors are insufficient to separate SHs from ICL. Future wide-field surveys (e.g., LSST, WEAVE, 4MOST) must rely on chemo-dynamical tracers (metallicity gradients, $\alpha$-enhancement, and kinematics) to disentangle these components.
• Progenitor Mass Spectrum: The evolution of metallicity gradients serves as a diagnostic for the mass spectrum of accreted satellites: early assembly is dominated by metal-poor dwarfs, while late-time growth is driven by the stripping of massive, metal-rich satellites.
• Model Validation: The FEGA25 model successfully reproduces the statistical trends of SH scaling relations and colors, validating the physical treatment of ICL formation and the transition radius definition, despite discrepancies in absolute metallicity values which point to environmental selection effects.