L-GALAXIES 2020: The formation and chemical evolution of stellar haloes in Milky Way Analogues and galaxy clusters

We present an analysis of the formation and chemical evolution of stellar haloes around (a) Milky Way Analogue (MWA) galaxies and (b) galaxy clusters in the L-Galaxies 2020 semi-analytic model of galaxy evolution. Observed stellar halo properties are better reproduced when assuming a gradual stripping model for the removal of cold gas and stars from satellites, compared to an instantaneous stripping model. The slope of the stellar mass – metallicity relation for MWA stellar haloes is in good agreement with that observed in the local Universe. This extends the good agreement between L-Galaxies 2020 and metallicity observations from the gas and stars inside galaxies to those outside. Halo stars contribute on average only ∼0.1 per cent of the total circumgalactic medium (CGM) enrichment by z = 0 in MWAs, ejecting predominantly carbon produced by AGB stars. Around a quarter of MWAs have a single ‘significant progenitor’ with a mean mass of ∼ 2.3 × 109M⊙, similar to that measured for Gaia Enceladus. For galaxy clusters, L-Galaxies 2020 shows good correspondence with observations of stellar halo mass fractions, but slightly over-predicts stellar masses. Assuming a Navarro-Frenk-White profile for the stellar halo mass distribution provides the best agreement. On average, the intracluster stellar component (ICS) is responsible for 5.4 per cent of the total intracluster medium (ICM) iron enrichment, exceeding the contribution from the brightest cluster galaxy (BCG) by z = 0. We find that considering gradual stripping of satellites and realistic radial profiles is crucial for accurately modelling stellar halo formation on all scales in semi-analytic models.

Rapid build-up of the stellar content in the protocluster core SPT2349−56 at z = 4.3

The protocluster SPT2349−56 at z = 4.3 contains one of the most actively star-forming cores known, yet constraints on the total stellar mass of this system are highly uncertain. We have therefore carried out deep optical and infrared observations of this system, probing rest-frame ultraviolet to infrared wavelengths. Using the positions of the spectroscopically-confirmed protocluster members, we identify counterparts and perform detailed source deblending, allowing us to fit spectral energy distributions in order to estimate stellar masses. We show that the galaxies in SPT2349−56 have stellar masses proportional to their high star-formation rates, consistent with other protocluster galaxies and field submillimetre galaxies (SMGs) around redshift 4. The galaxies in SPT2349−56 have on average lower molecular gas-to-stellar mass fractions and depletion timescales than field SMGs, although with considerable scatter. We construct the stellar-mass function for SPT2349−56 and compare it to the stellar-mass function of z = 1 galaxy clusters, finding consistent shapes between the two. We measure rest-frame galaxy ultraviolet half-light radii from our HST-F160W imaging, finding that on average the galaxies in our sample are similar in size to typical star-forming galaxies at these redshifts. However, the brightest HST-detected galaxy in our sample, found near the luminosity-weighted centre of the protocluster core, remains unresolved at this wavelength. Hydrodynamical simulations predict that the core galaxies will quickly merge into a brightest cluster galaxy, thus our observations provide a direct view of the early formation mechanisms of this class of object.

HST/WFC3 Grism Observations of z ∼ 1 Clusters: Evidence for Rapid Outside-in Environmental Quenching from Spatially Resolved Hα Maps

We present and publicly release (www.gclasshst.com) the first spatially resolved Hα maps of star-forming cluster galaxies at z ∼ 1, made possible with the Wide Field Camera 3 (WFC3) G141 grism on the Hubble Space Telescope (HST). Using a similar but updated method to 3D-HST in the field environment, we stack the Hα maps in bins of stellar mass, measure the half-light radius of the Hα distribution, and compare it to the stellar continuum. The ratio of the Hα to stellar continuum half-light radius, R [ H α / C ] = R eff , H α R eff , Cont , is smaller in the clusters by (6 ± 9)%, but statistically consistent within 1σ uncertainties. A negligible difference in R[Hα/C] with environment is surprising, given the higher quenched fractions in the clusters relative to the field. We postulate that the combination of high quenched fractions and no change in R[Hα/C] with environment can be reconciled if environmental quenching proceeds rapidly. We investigate this hypothesis by performing similar analysis on the spectroscopically confirmed, recently quenched cluster galaxies. 87% have Hα detections, with star formation rates 8 ± 1 times lower than star-forming cluster galaxies of similar stellar mass. Importantly, these galaxies have an R[Hα/C] that is (81 ± 8)% smaller than coeval star-forming field galaxies at fixed stellar mass. This suggests the environmental quenching process occurred outside-in. We conclude that disk truncation due to ram pressure stripping is occurring in cluster galaxies at z ∼ 1, but more rapidly and/or efficiently than in z ≲ 0.5 clusters, such that the effects on R[Hα/C] become observable just after the cluster galaxy has recently quenched.

RELICS: ICL Analysis of the z = 0.566 Merging Cluster WHL J013719.8–08284

We present a pilot study of the intracluster light (ICL) in massive clusters using imaging of the z = 0.566 cluster of galaxies WHL J013719.8–08284 observed by the RELICS project with the Hubble Space Telescope. We measure the ICL fraction in four optical ACS/WFC filters (F435W, F475W, F606W, and F814W) and five infrared WFC3/IR bands (F105W, F110W, F125W, F140W, and F160W). The ICL maps are calculated using the free-of-a-priori-assumptions algorithm CICLE, and the cluster membership is estimated from photometric properties. We find optical ICL fractions that range between ∼6% and 19%, in nice agreement with the values found in previous works for merging clusters. We also observe an ICL fraction excess between 3800 Å and 4800 Å, previously identified as a signature of merging clusters at 0.18 < z < 0.55. This excess suggests the presence of an enhanced population of young/low-metallicity stars in the ICL. All indicators thus point to WHL J013719.8–08284 as a disturbed cluster with a significant amount of recently injected stars, bluer than the average stars hosted by the cluster members and likely stripped out from infalling galaxies during the current merging event. Infrared ICL fractions are ∼50% higher than optical ones, which could be signatures of an older and/or higher-metallicity ICL population that can be associated with the buildup of the brightest cluster galaxy, passive evolution of previously injected young stars, or preprocessing in infalling groups. Finally, investigating the photometry of the cluster members, we tentatively conclude that WHL J013719.8–08284 fulfills the expected conditions for a fossil system progenitor.

Halo cluster shapes: insights from simulated galaxies and ICL with prospects for weak lensing applications

ABSTRACT We present a detailed study of the shapes and alignments of different galaxy cluster components using hydrodynamical simulations. We compute shape parameters from the dark matter (DM) distribution, the galaxy members and the intracluster light (ICL). We assess how well the DM cluster shape can be constrained by means of the identified galaxy member positions and the ICL. Further, we address the dilution factor introduced when estimating the cluster elongation using weak-lensing stacking techniques, which arises due to the misalignment between the total surface mass distribution and the distribution of luminous tracers. The dilution is computed considering the alignment between the DM and the brightest cluster galaxy, the galaxy members and the ICL. Our study shows that distributions of galaxy members and ICL are less spherical than the DM component, although both are well aligned with the semimajor axis of the latter. We find that the distribution of galaxy members hosted in more concentrated subhaloes is more elongated than the distribution of the DM. Moreover, these galaxies are better aligned with the DM component compared to the distribution of galaxies hosted in less concentrated subhaloes. We conclude that the positions of galaxy members can be used as suitable tracers to estimate the cluster surface density orientation, even when a low number of members is considered. Our results provide useful information for interpreting the constraints on the shapes of galaxy clusters in observational studies.

Probing the circumnuclear environment of NGC 1275 with High-Resolution X-ray spectroscopy

NGC 1275 is the Brightest Cluster Galaxy (BCG) in the Perseus cluster and hosts the active galactic nucleus (AGN) that is heating the central 100 kpc of the intracluster medium (ICM) atmosphere via a regulated feedback loop. Here we use a deep (490 ks) Cycle-19 Chandra High-Energy Transmission Grating (HETG) observation of NGC 1275 to study the anatomy of this AGN. The X-ray continuum is adequately described by an unabsorbed power-law with photon index Γ ≈ 1.9, creating strong tension with the detected column of molecular gas seen via HCN and HCO+ line absorption against the parsec-scale core/jet. This tension is resolved if we permit a composite X-ray source; allowing a column of $N_H\sim 8\times 10^{22}\hbox{${\rm \, cm}^{-2}\, $}$ to cover ∼15 per cent of the X-ray emitter does produce a significant improvement in the statistical quality of the spectral fit. We suggest that the dominant unabsorbed component corresponds to the accretion disk corona, and the sub-dominant X-ray component is the jet working surface and/or jet cocoon that is expanding into clumpy molecular gas. We suggest that this may be a common occurence in BCG-AGN. We conduct a search for photoionized absorbers/winds and fail to detect such a component, ruling out columns and ionization parameters often seen in many other Seyfert galaxies. We detect the 6.4 keV iron-Kα fluorescence line seen previously by XMM-Newton and Hitomi. We describe an analysis methodology which combines dispersive HETG spectra, non-dispersive microcalorimeter spectra, and sensitive XMM-Newton/EPIC spectra in order to constrain (sub)arcsec-scale extensions of the iron-Kα emission region.

Accurate modeling of the strong and weak lensing profiles for the galaxy clusters Abell 1689 and 1835

An accurate, spherically symmetric description of the mass distribution is presented for two quite virialized galaxy clusters, Abell 1689 and Abell 1835. A suitable regularization of the small eigenvalues of the covariance matrices is introduced. A stretched exponential profile is assumed for the brightest cluster galaxy. A similar stretched exponential for the dark matter and halo galaxies combined, functions well, as do thermal fermions for the dark matter and a standard profile for the halo galaxies. To discriminate between them, sensitive tests have been identified and applied. A definite verdict can follow from sharp data near the cluster centers and beyond 1 Mpc.

Discovery of an ancient massive merger event in the For- nax cluster galaxy NGC 1380

Driven by gravity, galaxies are expected to continuously grow through the merging of smaller systems. To derive their past merger history is challenging, as the accreted stars disperse quickly; yet, it is a needed step to test the theory of hierarchical evolution. The merger histories of the most massive Local Group spirals, the Milky Way and M31, have been re- cently uncovered by using the motion and chemistry of their individual stars. On the other hand, the details of the merger history of galaxies at larger distance have so far remained hidden. Here we report the discovery of an ancient, massive merger event in the lenticu- lar galaxy NGC 1380 in the Fornax cluster. By applying a recently developed population-orbital superposition model (Zhu at al 2020) to NGC 1380’s surface brightness as well as stellar kinematic, age, and metallicity maps from VLT/MUSE IFU data (Sarzi et al 2018), we obtain the stellar orbits, age and metallicity distributions of this galaxy. The highly radial orbits which make up an inner stellar halo are ∼ 13 Gyr old with metallicity Z/Z⊙ ∼ 1.2 and comprise a stellar mass of M∗,halo(r<2Re)∼3.4×10^10 M⊙. By comparing to analogues from the cosmological galaxy simulation TNG50 (Pillepich 2019), we find that the formation of the inner stellar halo of NGC 1380 requires a merger with a massive satellite galaxy with stellar mass of ∼ 3 × 10^10 M⊙ that occurred roughly ∼ 10 Gyr ago. Moreover, we infer the total accreted stellar mass of NGC 1380 to be ∼ 6 × 10^10 M⊙. The massive merger in NGC 1380 is the first major merger event found in a normal phase-mixed galaxy beyond the Local Volume, and it is the oldest and most massive one identified in nearby galaxies so far. Our chemo-dynamical method, when applied to extended deep IFU data and in combination with cosmological galaxy simulations, can quantitatively unravel the merger history of a large number of nearby galaxies.

WISDOM project – VIII. Multiscale feedback cycles in the brightest cluster galaxy NGC 0708

ABSTRACT We present high-resolution (synthesized beam size 0${_{.}^{\prime\prime}}$088 × 0${_{.}^{\prime\prime}}$083 or 25 × 23 pc2) Atacama Large Millimetre/submillimetre Array 12CO(2–1) line and 236 GHz continuum observations, as well as 5 GHz enhanced Multi-Element Radio Linked Interferometer Network (e-MERLIN) continuum observations, of NGC 0708; the brightest galaxy in the low-mass galaxy cluster Abell 262. The line observations reveal a turbulent, rotating disc of molecular gas in the core of the galaxy, and a high-velocity, blueshifted feature ≈0${_{.}^{\prime\prime}}$4 (≈113 pc) from its centre. The submillimetre continuum emission peaks at the nucleus, but extends towards this anomalous CO emission feature. No corresponding elongation is found on the same spatial scales at 5 GHz with e-MERLIN. We discuss potential causes for the anomalous blueshifted emission detected in this source, and conclude that it is most likely to be a low-mass in-falling filament of material condensing from the hot intracluster medium via chaotic cold accretion, but it is also possible that it is a jet-driven molecular outflow. We estimate the physical properties this structure has in these two scenarios, and show that either explanation is viable. We suggest future observations with integral field spectrographs will be able to determine the true cause of this anomalous emission, and provide further evidence for interaction between quenched cooling flows and mechanical feedback on both small and large scales in this source.

The WaZP galaxy cluster sample of the Dark Energy Survey Year 1

We present a new (2+1)D galaxy cluster finder based on photometric redshifts called Wavelet Z Photometric (wa zp) applied to DES first year (Y1A1) data. The results are compared to clusters detected by the South Pole Telescope (SPT) survey and the redMaPPer cluster finder, the latter based on the same photometric data. wa zp searches for clusters in wavelet-based density maps of galaxies selected in photometric redshift space without any assumption on the cluster galaxy populations. The comparison to other cluster samples was performed with a matching algorithm based on angular proximity and redshift difference of the clusters. It led to the development of a new approach to match two optical cluster samples, following an iterative approach to minimize incorrect associations. The wa zp cluster finder applied to DES Y1A1 galaxy survey (1,511.13 deg2 up to mi = 23 mag) led to the detection of 60,547 galaxy clusters with redshifts 0.05 &lt; z &lt; 0.9 and richness Ngals ≥ 5. Considering the overlapping regions and redshift ranges between the DES Y1A1 and SPT cluster surveys, all sz based SPT clusters are recovered by the wa zp sample. The comparison between wa zp and redMaPPer cluster samples showed an excellent overall agreement for clusters with richness Ngals (λ for redMaPPer) greater than 25 (20), with 95% recovery on both directions. Based on the cluster cross-match we explore the relative fragmentation of the two cluster samples and investigate the possible signatures of unmatched clusters.