We present the first results from ROMULUSC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. ROMULUSC, a zoom-in simulation of a halo with z = 0 mass 1014 ...M⊙, is run with the same sub-grid physics and resolution as ROMULUS25. With unprecedented mass and spatial resolution, ROMULUSC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that ROMULUSC results in an intracluster medium consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low-entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster centre. This enhanced quenched population extends beyond R200 and is in place at high redshift. Similarly, we predict that an SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest active galactic nucleus in the most massive galaxies.
Abstract
We analyse spatially resolved deep optical spectroscopy of brightestcluster galaxies (BCGs) located in 32 massive clusters with redshifts of 0.05 ≤ z ≤ 0.30 to investigate their velocity ...dispersion profiles. We compare these measurements to those of other massive early-type galaxies, as well as central group galaxies, where relevant. This unique, large sample extends to the most extreme of massive galaxies, spanning MK between –25.7 and –27.8 mag, and host cluster halo mass M500 up to 1.7 × 1015 M⊙. To compare the kinematic properties between brightest group and cluster members, we analyse similar spatially resolved long-slit spectroscopy for 23 nearby brightest group galaxies (BGGs) from the Complete Local-Volume Groups Sample. We find a surprisingly large variety in velocity dispersion slopes for BCGs, with a significantly larger fraction of positive slopes, unique compared to other (non-central) early-type galaxies as well as the majority of the brightest members of the groups. We find that the velocity dispersion slopes of the BCGs and BGGs correlate with the luminosity of the galaxies, and we quantify this correlation. It is not clear whether the full diversity in velocity dispersion slopes that we see is reproduced in simulations.
We present an analysis of the levels and evolution of star formation activity in a representative sample of 30 massive galaxy clusters at 0.15 < z < 0.30 from the Local Cluster Substructure Survey, ...combining wide-field Spitzer/MIPS 24 mum data with extensive spectroscopy of cluster members. The specific SFRs of massive (M > ~ 10 super(10) M sub(middot in circle)) star-forming cluster galaxies within r sub(200) are found to be systematically ~28% lower than their counterparts in the field at fixed stellar mass and redshift, a difference significant at the 8.7sigma level. This is the unambiguous signature of star formation in most (and possibly all) massive star-forming galaxies being slowly quenched upon accretion into massive clusters, their star formation rates (SFRs) declining exponentially on quenching timescales in the range 0.7-2.0 Gyr. We measure the mid-infrared Butcher-Oemler effect over the redshift range 0.0-0.4, finding rapid evolution in the fraction (functionof sub(SF)) of massive (M sub(K) < - 23.1) cluster galaxies within r sub(200) with SFRs > 3 M sub(middot in circle) yr super(-1), of the form functionof sub(SF) alpha (1 + z) super(7.6 + or - 1.1). We dissect the origins of the Butcher-Oemler effect, revealing it to be due to the combination of a ~3 x decline in the mean specific SFRs of star-forming cluster galaxies since z ~ 0.3 with a ~1.5 x decrease in number density. Two-thirds of this reduction in the specific SFRs of star-forming cluster galaxies is due to the steady cosmic decline in the specific SFRs among those field galaxies accreted into the clusters. The remaining one-third reflects an accelerated decline in the star formation activity of galaxies within clusters. The slow quenching of star formation in cluster galaxies is consistent with a gradual shut down of star formation in infalling spiral galaxies as they interact with the intracluster medium via ram-pressure stripping or starvation mechanisms. The observed sharp decline in star formation activity among cluster galaxies since z ~ 0.4 likely reflects the increased susceptibility of low-redshift spiral galaxies to gas removal mechanisms as their gas surface densities decrease with time. We find no evidence for the build-up of cluster S0 bulges via major nuclear starburst episodes.
ABSTRACT We present a study of the spatial distribution and kinematics of star-forming galaxies in 30 massive clusters at 0.15 < z < 0.30, combining wide-field Spitzer 24 m and GALEX near-ultraviolet ...imaging with highly complete spectroscopy of cluster members. The fraction (fSF) of star-forming cluster galaxies rises steadily with cluster-centric radius, increasing fivefold by 2r200, but remains well below field values even at 3r200. This suppression of star formation at large radii cannot be reproduced by models in which star formation is quenched in infalling field galaxies only once they pass within r200 of the cluster, but is consistent with some of them being first pre-processed within galaxy groups. Despite the increasing fSF-radius trend, the surface density of star-forming galaxies actually declines steadily with radius, falling ∼15× from the core to 2r200. This requires star formation to survive within recently accreted spirals for 2-3 Gyr to build up the apparent over-density of star-forming galaxies within clusters. The velocity dispersion profile of the star-forming galaxy population shows a sharp peak of 1.44 at 0.3r500, and is 10%-35% higher than that of the inactive cluster members at all cluster-centric radii, while their velocity distribution shows a flat, top-hat profile within r500. All of these results are consistent with star-forming cluster galaxies being an infalling population, but one that must also survive ∼0.5-2 Gyr beyond passing within r200. By comparing the observed distribution of star-forming galaxies in the stacked caustic diagram with predictions from the Millennium simulation, we obtain a best-fit model in which star formation rates decline exponentially on quenching timescales of 1.73 0.25 Gyr upon accretion into the cluster.
We present a weak-lensing study of the cluster A520 based on Advanced Camera for Surveys (ACS) data. The excellent data quality provides a mean source density of ~109 arcmin super(-2), which improves ...both resolution and significance of the mass reconstruction compared to a previous study based on Wide Field Planetary Camera 2 (WFPC2) images. We take care in removing instrumental effects such as the charge trailing due to radiation damage of the detector and the position-dependent point-spread function. This new ACS analysis confirms the previous claims that a substantial amount of dark mass is present between two luminous subclusters where we observe very little light. The centroid of the dark peak in the current ACS analysis is offset to the southwest by ~1' with respect to the centroid from the WFPC2 analysis. Interestingly, this new centroid is in better agreement with the location where the X-ray emission is strongest, and the mass-to-light ratio estimated with this centroid is much higher (813 + or - 78 M sub(middot in circle)/L sub(R) sub(middot in circle)) than the previous value; the aperture mass with the WFPC2 centroid provides a consistent mass. Although we cannot provide a definite explanation for the dark peak, we discuss a revised scenario, wherein dark matter with a more conventional range (sigma sub(DM)/m sub(DM) < 1 cm super(2) g super(-1)) of self-interacting cross-section can lead to the detection of this dark substructure. If supported by detailed numerical simulations, this hypothesis opens up the possibility that the A520 system can be used to establish a lower limit of the self-interacting cross-section of dark matter.
Virtually all massive galaxies, including our own, host central black holes ranging in mass from millions to billions of solar masses. The growth of these black holes releases vast amounts of energy ...that powers quasars and other weaker active galactic nuclei. A tiny fraction of this energy, if absorbed by the host galaxy, could halt star formation by heating and ejecting ambient gas. A central question in galaxy evolution is the degree to which this process has caused the decline of star formation in large elliptical galaxies, which typically have little cold gas and few young stars, unlike spiral galaxies.
We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at 0.15 < z < 0.3 from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing ...measurements of M
500 control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM–Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is
$\mathrel {\beta _{\rm X}}\;= 0.95\pm 0.05$
, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is
$\mathrel {\beta _{\rm P}}\;= 0.95\pm 0.04$
. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project, and Weighing the Giants agree on
$\mathrel {\beta _{\rm P}}\;\simeq 0.9{\rm -}0.95$
at 0.15 < z < 0.3. This small level of hydrostatic bias disagrees at ∼5σ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts.
ABSTRACT
We use stellar and dynamical mass profiles, combined with a stellar population analysis, of 32 brightest cluster galaxies (BCGs) at redshifts of 0.05 ≤$z$ ≤ 0.30, to place constraints on ...their stellar initial mass function (IMF). We measure the spatially resolved stellar population properties of the BCGs, and use it to derive their stellar mass-to-light ratios ($\Upsilon _{\star \rm POP}$). We find young stellar populations (<200 Myr) in the centres of 22 per cent of the sample, and constant $\Upsilon _{\star \rm POP}$ within 15 kpc for 60 per cent of the sample. We further use the stellar mass-to-light ratio from the dynamical mass profiles of the BCGs ($\Upsilon _{\star \rm DYN}$), modelled using a multi-Gaussian expansion and Jeans Anisotropic Method, with the dark matter contribution explicitly constrained from weak gravitational lensing measurements. We directly compare the stellar mass-to-light ratios derived from the two independent methods, $\Upsilon _{\star \rm POP}$ (assuming some IMF) to $\Upsilon _{\star \rm DYN}$ for the subsample of BCGs with no young stellar populations and constant $\Upsilon _{\star \rm POP}$. We find that for the majority of these BCGs, a Salpeter (or even more bottom-heavy) IMF is needed to reconcile the stellar population and dynamical modelling results although for a small number of BCGs, a Kroupa (or even lighter) IMF is preferred. For those BCGs better fit with a Salpeter IMF, we find that the mass-excess factor against velocity dispersion falls on an extrapolation (towards higher masses) of known literature correlations. We conclude that there is substantial scatter in the IMF amongst the highest mass galaxies.