No Arabic abstract
We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the `stellar splashback feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (mu >> 32 mag/arcsec^2) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid.
We present an analysis of the intracluster light in the Frontier Field Cluster MACS J1149.5+2223 (z=0.544), which combines new and archival Hubble WFC3/IR imaging to provide continuous radial coverage out to 2.8 Mpc from the brightest cluster galaxy. Employing careful treatment of potential systematic biases and using data at the largest radii to determine the background sky level, we reconstruct the surface brightness profile out to a radius of 2 Mpc. This radius is the largest to which the intracluster light (ICL) has been measured for an individual cluster. Within this radius, we measure a total luminosity of 1.5e13 Lsun for the brightest cluster galaxy plus ICL light. From the profile and its logarithmic slope, we identify the transition from the brightest cluster galaxy to ICL at r~70 kpc. Remarkably, we also detect an inflection in the profile centered in the 1.2-1.7 Mpc (0.37-0.52 r200m) radial bin, a signature of an infall caustic in the stellar distribution. Based upon the shape and strength of the feature, we interpret it as potentially being at the splashback radius, although the radius is smaller than theoretical predictions. If this is the splashback radius, then it is the first such detection in the ICL and the first detection of the splashback radius for an individual cluster. Similar analyses should be possible with the other Frontier Field clusters, and eventually with clusters from the Euclid and Roman missions.
We have explored the outskirts of dark matter haloes out to 2.5 times the virial radius using a large sample of halos drawn from Illustris, along with a set of zoom simulations (MUGS). Using these, we make a systematic exploration of the shape profile beyond R$_{vir}$. In the mean sphericity profile of Illustris halos we identify a dip close to the virial radius, which is robust across a broad range of masses and infall rates. The inner edge of this feature may be related to the virial radius and the outer edge with the splashback radius. Due to the high halo-to-halo variation this result is visible only on average. However, in four individual halos in the MUGS sample, a decrease in the sphericity and a subsequent recovery is evident close to the splashback radius. We find that this feature persists for several Gyr, growing with the halo. This feature appears at the interface between the spherical halo density distribution and the filamentary structure in the environment. The shape feature is strongest when there is a high rate of infall, implying that the effect is due to the mixing of accreting and virializing material. The filamentary velocity field becomes rapidly mixed in the halo region inside the virial radius, with the area between this and the splashback radius serving as the transition region. We also identify a long-lasting and smoothly evolving splashback region in the radial density gradient in many of the MUGS halos.
We study the Intra-Halo Stellar Component (IHSC) of Milky Way-mass systems up to galaxy clusters in the Horizon-AGN cosmological hydrodynamical simulation. We identify the IHSC using an improved phase-space galaxy finder algorithm which provides an adaptive, physically motivated and shape-independent definition of this stellar component, that can be applied to halos of arbitrary masses. We explore the IHSC mass fraction-total halos stellar mass, $f_{M*,IHSC}-M*$, relation and the physical drivers of its scatter. We find that on average the $f_{M*,IHSC}$ increases with $M_{*,tot}$, with the scatter decreasing strongly with mass from 2 dex at $M_{*,tot}sim10^{11}M_odot$ to 0.3 dex at group masses. At high masses, $M_{*,tot}>10^{11.5}M_odot$, $f_{M*,IHSC}$ increases with the number of substructures, and with the mass ratio between the central galaxy and largest satellite, at fixed $M_{*,tot}$. From mid-size groups and systems below $M_{*,tot}<10^{12}M_odot$, we find that the central galaxys stellar rotation-to-dispersion velocity ratio, V/{sigma}, displays the strongest (anti)-correlation with $f_{M*,IHSC}$ at fixed $M_{*,tot}$ of all the galaxy and halo properties explored, transitioning from $f_{M*,IHSC}$<0.1% for high V/{sigma}, to $f_{M*,IHSC}sim5$% for low V/{sigma} galaxies. By studying the $f_{M*,IHSC}$ temporal evolution, we find that, in the former, mergers not always take place, but if they did, they happened early (z>1), while the high $f_{M*,IHSC}$ population displays a much more active merger history. In the case of massive groups and galaxy clusters, $M_{*,tot}>10^{12}M_odot$, a fraction $f_{M*,IHSC}sim$10-20% is reached at $zsim1$ and then they evolve across lines of constant $f_{M*,IHSC}$ modulo some small perturbations. Because of the limited simulations volume, the latter is only tentative and requires a larger sample of simulated galaxy clusters to confirm.
The largest stellar halos in the universe are found in massive galaxy clusters, where interactions and mergers of galaxies, along with the cluster tidal field, all act to strip stars from their host galaxies and feed the diffuse intracluster light (ICL) and extended halos of brightest cluster galaxies (BCGs). Studies of the nearby Virgo Cluster reveal a variety of accretion signatures imprinted in the morphology and stellar populations of its ICL. While simulations suggest the ICL should grow with time, attempts to track this evolution across clusters spanning a range of mass and redshift have proved difficult due to a variety of observational and definitional issues. Meanwhile, studies of nearby galaxy groups reveal the earliest stages of ICL formation: the extremely diffuse tidal streams formed during interactions in the group environment.
The intracluster light (ICL) is a luminous component of galaxy clusters composed of stars that are gravitationally bound to the cluster potential but do not belong to the individual galaxies. Previous studies of the ICL have shown that its formation and evolution are intimately linked to the evolutionary stage of the cluster. Thus, the analysis of the ICL in the Coma cluster will give insights into the main processes driving the dynamics in this highly complex system. Using a recently developed technique, we measure the ICL fraction in Coma at several wavelengths, using the J-PLUS unique filter system. The combination of narrow- and broadband filters provides valuable information on the dynamical state of the cluster, the ICL stellar types, and the morphology of the diffuse light. We use the Chebyshev-Fourier Intracluster Light Estimator (CICLE) to disentangle the ICL from the light of the galaxies, and to robustly measure the ICL fraction in seven J-PLUS filters. We obtain the ICL fraction distribution of the Coma cluster at different optical wavelengths, which varies from $sim 7%-21%$, showing the highest values in the narrowband filters J0395, J0410, and J0430. This ICL fraction excess is distinctive pattern recently observed in dynamically active clusters (mergers), indicating a higher amount of bluer stars in the ICL compared to the cluster galaxies. Both the high ICL fractions and the excess in the bluer filters are indicative of a merging state. The presence of younger/lower-metallicity stars the ICL suggests that the main mechanism of ICL formation for the Coma cluster is the stripping of the stars in the outskirts of infalling galaxies and, possibly, the disruption of dwarf galaxies during past/ongoing mergers.