Qian-Hui Chen

Stellar halos are the only major stellar component of disk galaxies that lack systematic observational characterization, yet they encode critical information about galaxy merger histories. We present the first systematic census of stellar halos in a large, flux-limited sample of 169 high-inclination central galaxies with stellar masses 7.3 ≤ log M ⋆ / M ⊙ ≤ 11.0 and redshift z < 0.1, using Hyper Suprime-Cam Subaru Strategic Survey Deep optical images. Stellar halos are detected in 93 galaxies, primarily through their low isophotal ellipticities in the outskirts, improving upon conventional methods of stellar halo identification. The halo detection rate reaches ∼50% at log M ⋆ / M ⊙ > 9.9 and ≳70% for Milky Way (MW)-mass galaxies. We derive halo surface brightness profiles, colors, and masses, finding that stellar halos generally follow power-law radial profiles. Higher-mass galaxies, on average, exhibit smaller power-law indices and larger halo mass fractions, indicating more extended halos and more active merger histories. A significant stellar halo color–mass correlation, driven mainly by the mass–metallicity relation, suggests dominance by a few massive accretion events. MW-mass galaxies have a median stellar halo fraction of 10% ± 5%. Among nearby galaxies with halo measurements within 25 Mpc, two-thirds (including the MW) lie below the mean stellar halo fraction–galaxy mass relation. Overall, the nearby galaxies show a median halo deficit of ∼0.3 dex, implying unusually quiescent merger histories. We show that this deficit follows a broader trend in which typical halo fractions increase with heliocentric distance, tracking the gradual rise in matter density toward the cosmic average by z ≲ 0.07.

Stellar halos are the only major stellar component of disk galaxies that lack systematic observational characterization, yet they encode critical information about galaxy merger histories. We present the first systematic census of stellar halos in a large, flux-limited sample of 169 high-inclination central galaxies with stellar masses 7.3 <= log Mstar/Msun <= 11.0 and redshift z < 0.1, using HSC-SSP Deep optical images. Stellar halos are detected in 93 galaxies, primarily through their low isophotal ellipticities in the outskirts, improving upon conventional methods of stellar halo identification. The halo detection rate reaches ~ 50% at log Mstar/Msun > 9.9 and >= 70% for Milky Way (MW)-mass galaxies. We derive halo surface brightness profiles, colors, and masses, finding that stellar halos generally follow power-law radial profiles. Higher-mass galaxies, on average, exhibit smaller power-law indices and larger halo mass fractions, indicating more extended halos and more active merger histories. A significant stellar halo color-mass correlation, driven mainly by the mass-metallicity relation, suggests dominance by a few massive accretion events. MW-mass galaxies have a median stellar halo fraction of 10% +/- 5%. Among nearby galaxies with halo measurements within 25 Mpc, two thirds (including the MW) lie below the mean stellar halo fraction-galaxy mass relation. Overall, the nearby galaxies show a median halo deficit of ~ 0.3 dex, implying unusually quiescent merger histories. We show that this deficit follows a broader trend in which typical halo fractions increase with heliocentric distance, tracking the gradual rise in matter density toward the cosmic average by z <= 0.07.

The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star-forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across broad mass (10^8.0 M _⊙ ≤ M _⋆ ≲ 10^12.0 M _⊙) and redshift (0 ≤ z ≤ 8) ranges. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of 10^10.0–10^11.0 M _⊙, while galaxies with stellar mass < 10¹⁰ M _⊙ and >10¹¹ M _⊙ exhibit weaker redshift evolution. Our result of negative gradients at high redshift contrast with the many positive and flat gradients in the 1 < z < 4 observational literature. At z > 6, the negative gradients observed with JWST and the Atacama Large Millimeter/submillimeter Array are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients.

Radial metallicity gradients are fundamental to understanding galaxy formation and evolution. In our high-resolution simulation of a NIHAO-UHD Milky Way analogue, we analyze the linearity, scatter, spatial coherence, and age-related variations of metallicity gradients using young stars and gas. While a global linear model generally captures the gradient, it ever so slightly overestimates metallicity in the inner galaxy and underestimates it in the outer regions of our simulated galaxy. Both a quadratic model, showing an initially steeper gradient that smoothly flattens outward, and a piecewise linear model with a break radius around 9.3-11.5~kpc (2.4-3.0 effective radii) fit the data equally better. The spread of [Fe/H] of young stars in the simulation increases by tenfold from the innermost to the outer galaxy at a radius of 20 kpc. We find that stars born at similar times along radial spirals drive this spread in the outer galaxy, with a chemical under- and over-enhancement of up to 0.1 dex at leading and trailing regions of such spirals, respectively. This localised chemical variance highlights the need to examine radial and azimuthal selection effects for both Galactic and extragalactic observational studies. The arguably idealised but volume-complete simulations suggest that future studies should not only test linear and piecewise linear gradients, but also non-linear functions such as quadratic ones to test for a smooth gradient rather than one with a break radius. Either finding would help to determine the importance of different enrichment or mixing pathways and thus our understanding of galaxy formation and evolution scenarios.

The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy’s stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift (z ≲ 4) observational studies suggest that the FMR is redshift-invariant, recent high-z JWST data indicate deviations from the FMR established at low-z. In this study, we utilize the FMR to predict the evolution of the normalisation of the mass-metallicity relation (MZR) using the cosmological simulations Illustris, IllustrisTNG, EAGLE, and SIMBA. Our findings demonstrate that a z = 0 calibrated FMR struggles to predict the evolution in the MZR of each simulation. To quantify the divergence of the predictions, we introduce the concepts of a “static” FMR, where the role of the SFR in setting the normalization of the MZR does not change with redshift, and a “dynamic” FMR, where the role of SFR evolves over time. We find static FMRs in SIMBA and dynamic FMRs in Illustris, IllustrisTNG and EAGLE. We suggest that the differences between these models likely points to the subtle differences in the implementation of the baryon cycle. Moreover, we echo recent JWST results at z > 4 by finding significant offsets from the FMR in IllustrisTNG and EAGLE, suggesting that the observed FMR may have a similar dynamic trend as these simulations. Overall, our findings imply that the current FMR framework neglects important time variations of these simulations’ baryon cycles.