Feedback heating by cosmic rays in clusters of galaxies

Recent observations show that the cooling flows in the central regions of galaxy clusters are highly suppressed. Observed active galactic nuclei (AGN)-induced cavities/bubbles are a leading candidate for suppressing cooling, usually via some form of mechanical heating. At the same time, observed X-ray cavities and synchrotron emission point towards a significant non-thermal particle population. Previous studies have focused on the dynamical effects of cosmic ray pressure support, but none has built successful models in which cosmic ray heating is significant. Here, we investigate a new model of AGN heating, in which the intracluster medium is efficiently heated by cosmic rays, which are injected into the intra-cluster medium (ICM) through diffusion or the shredding of the bubbles by Rayleigh–Taylor or Kelvin–Helmholtz instabilities. We include thermal conduction as well. Using numerical simulations, we show that the cooling catastrophe is efficiently suppressed. The cluster quickly relaxes to a quasi-equilibrium state with a highly reduced accretion rate and temperature and density profiles which match observations. Unlike the conduction-only case, no fine-tuning of the Spitzer conduction suppression factor f is needed. The cosmic ray pressure, P _c/ P _g ≲ 0.1 and ∇ P _c≲ 0.1ρ g , is well within observational bounds. Cosmic ray heating is a very attractive alternative to mechanical heating, and may become particularly compelling if Gamma-ray Large Array Space Telescope ( GLAST ) detects the γ-ray signature of cosmic rays in clusters.     

Observable consequences of kinetic and thermal AGN feedback in elliptical galaxies and galaxy clusters

We have constructed an analytical model of active galactic nuclei (AGN) feedback and studied its implications for elliptical galaxies and galaxy clusters. The results show that momentum injection above a critical value will eject material from low-mass elliptical galaxies, and leads to an X-ray luminosity, L _X, that is  ∝σ⁸⁻¹⁰  , depending on the AGN fuelling mechanism, where σ is the velocity dispersion of the hot gas. This result agrees well with both observations and semi-analytic models. In more massive ellipticals and clusters, AGN outflows quickly become buoyancy dominated. This necessarily means that heating by a central cluster AGN redistributes the intracluster medium (ICM) such that the mass of hot gas, within the cooling radius, should be  ∝ L _X(< r _cool)/[ g ( r _cool)σ]  , where   g ( r _cool)  is the gravitational acceleration at the cooling radius. This prediction is confirmed using observations of seven clusters. The same mechanism also defines a critical ICM cooling time of  ∼0.5 Gyr  , which is in reasonable agreement with recent observations showing that star formation and AGN activity are triggered below a universal cooling time threshold.     

A statistically selected Chandra sample of 20 galaxy clusters – II. Gas properties and cool core/non-cool core bimodality

We investigate the thermodynamic and chemical structure of the intracluster medium (ICM) across a statistical sample of 20 galaxy clusters analysed with the Chandra X-ray satellite. In particular, we focus on the scaling properties of the gas density, metallicity and entropy and the comparison between clusters with and without cool cores (CCs). We find marked differences between the two categories except for the gas metallicity, which declines strongly with radius for all clusters  ( Z ∝ r ^−0.31)  , outside  ∼0.02 r ₅₀₀  . The scaling of gas entropy is non-self-similar and we find clear evidence of bimodality in the distribution of logarithmic slopes of the entropy profiles. With only one exception, the steeper sloped entropy profiles are found in CC clusters whereas the flatter slope population are all non-CC clusters. We explore the role of thermal conduction in stabilizing the ICM and conclude that this mechanism alone is sufficient to balance cooling in non-CC clusters. However, CC clusters appear to form a distinct population in which heating from feedback is required in addition to conduction. Under the assumption that non-CC clusters are thermally stabilized by conduction alone, we find the distribution of Spitzer conduction suppression factors, f _c, to be lognormal, with a log (base 10) mean of  −1.50 ± 0.03  (i.e.   f _c= 0.032  ) and log standard deviation  0.39 ± 0.02  .     

Environmental dependence of active galactic nuclei activity in the supercluster A901/2

We present XMM data for the supercluster A901/2, at   z ∼ 0.17  , which is combined with deep imaging and 17-band photometric redshifts (from the COMBO-17 survey), two degree field (2dF) spectra and Spitzer 24 μm data, to identify active galactic nuclei (AGN) in the supercluster. The 90 ksec XMM image contains 139 point sources, of which 11 are identified as supercluster AGN with   L _{X(0.5−7.5 keV)} > 1.7 × 10⁴¹ erg cm⁻² s⁻¹  . The host galaxies have   M _R < −20  and only two of eight sources with spectra could have been identified as AGN by the detected optical emission lines. Using a large sample of 795 supercluster galaxies, we define control samples of massive galaxies with no detected AGN. The local environments of the AGN and control samples differ at ≳98 per cent significance. The AGN host galaxies lie predominantly in areas of moderate projected galaxy density and with more local blue galaxies than the control sample, with the exception of one very bright type I AGN very near the centre of a cluster. These environments are similar to, but not limited to, cluster outskirts and blue groups. Despite the large number of potential host galaxies, no AGN are found in regions with the highest galaxy density (excluding some cluster cores where emission from the intra-cluster medium obscures moderate luminosity AGN). AGN are also absent from the areas with lowest galaxy density. We conclude that the prevalence of cluster AGN is linked to their environment.     

The distribution of active galactic nuclei in a large sample of galaxy clusters

We present an analysis of the X-ray point source populations in 182 Chandra images of galaxy clusters at   z > 0.1  with exposure time >10 ks, as well as 44 non-cluster fields. The analysis of the number and flux of these sources, using a detailed pipeline to predict the distribution of non-cluster sources in each field, reveals an excess of X-ray point sources associated with the galaxy clusters. A sample of 148 galaxy clusters at  0.1 < z < 0.9  , with no other nearby clusters, shows an excess of 230 cluster sources in total, an average of ∼1.5 sources per cluster. The lack of optical data for these clusters limits the physical interpretation of this result, as we cannot calculate the fraction of cluster galaxies hosting X-ray sources. However, the fluxes of the excess sources indicate that over half of them are very likely to be active galactic nuclei (AGN), and the radial distribution shows that they are quite evenly distributed over the central 1 Mpc of the cluster, with almost no sources found beyond this radius. We also use this pipeline to successfully reproduce the results of previous studies, particularly the higher density of sources in the central 0.5 Mpc of a few cluster fields, but show that these conclusions are not generally valid for this larger sample of clusters. We conclude that some of these differences may be due to the sample properties, such as the size and redshift of the clusters studied, or a lack of publications for cluster fields with no excess sources. This paper also presents the basic X-ray properties of the galaxy clusters, and in subsequent papers in this series the dependence of the AGN population on these cluster properties will be evaluated.
In addition the properties of over 9500 X-ray point sources in the fields of galaxy clusters are tabulated in a separate catalogue available online or at http://www.sc.eso.org~rgilmour .     

Hydrodynamical simulations of cluster formation with central AGN heating

We analyse a hydrodynamical simulation model for the recurrent heating of the central intra-cluster medium (ICM) by active galactic nuclei (AGN). Besides the self-gravity of the dark matter and gas components, our approach includes the radiative cooling and photoheating of the gas, as well as a subresolution multiphase model for star formation and supernova feedback. Additionally, we incorporate a periodic heating mechanism in the form of hot, buoyant bubbles, injected into the intragalactic medium (IGM) during the active phases of the accreting central AGN. We use simulations of isolated cluster haloes of different masses to study the bubble dynamics and heat transport into the IGM. We also apply our model to self-consistent cosmological simulations of the formation of galaxy clusters with a range of masses. Our numerical schemes explore a variety of different assumptions for the spatial configuration of AGN-driven bubbles, for their duty cycles and for the energy injection mechanism, in order to obtain better constraints on the underlying physical picture. We argue that AGN heating can substantially affect the properties of both the stellar and gaseous components of clusters of galaxies. Most importantly, it alters the properties of the central dominant (cD) galaxy by reducing the mass deposition rate of freshly cooled gas out of the ICM, thereby offering an energetically plausible solution to the cooling-flow problem. At the same time, this leads to reduced or eliminated star formation in the central cD galaxy, giving it red stellar colours as observed.     

ASCA spectroscopy of the luminous infrared galaxy NGC 6240: X-ray emission from a starburst and a buried active nucleus

We present an X-ray spectroscopic study of the prototype far-infrared galaxy NGC 6240 from ASCA . The soft X-ray spectrum (below 2 keV) shows clear signatures of thermal emission well described by a multitemperature optically thin plasma, which probably originates in a powerful starburst. Strong hard X-ray emission is also detected with ASCA and its spectrum above 3 keV is extremely flat with a prominent iron K line complex, very similar to that seen in the Seyfert 2 galaxy NGC 1068 but about an order of magnitude more luminous ( L _3−10keV ≈ 1.4 × 10⁴² erg s⁻¹). The hard X-ray spectrum indicates that only reflected X-rays of an active galactic nucleus (AGN) buried in a heavy obscuration ( N _H > 2 × 10²⁴ cm⁻²) are visible. This is evidence for an AGN in NGC 6240, emitting possibly at a quasar luminosity (∼ 10⁴⁵ erg s⁻¹), and suggests its significant contribution to the far-infrared luminosity.     

Are fossil groups a challenge of the cold dark matter paradigm?

We study six groups and clusters of galaxies suggested in the literature to be 'fossil' systems (i.e. to have luminous diffuse X-ray emission and a magnitude gap of at least 2 mag R between the first and the second ranked member within half of the virial radius), each having good quality X-ray data and Sloan Digital Sky Survey (SDSS) spectroscopic or photometric coverage out to the virial radius. The poor cluster AWM 4 is clearly established as a fossil system, and we confirm the fossil nature of four other systems (RX J1331.5+1108, RX J1340.6+4018, RX J1256.0+2556 and RX J1416.4+2315), while the cluster RX J1552.2+2013 is disqualified as fossil system. For all systems, we present the luminosity functions within 0.5 and 1 virial radius that are consistent, within the uncertainties, with the universal luminosity function of clusters. For the five bona fide fossil systems, having a mass range  2 × 10¹³–3 × 10¹⁴ M_⊙  , we compute accurate cumulative substructure distribution functions (CSDFs) and compare them with the CSDFs of observed and simulated groups/clusters available in the literature. We demonstrate that the CSDFs of fossil systems are consistent with those of normal observed clusters and do not lack any substructure with respect to simulated galaxy systems in the cosmological Λ cold dark matter (ΛCDM) framework. In particular, this holds for the archetype fossil group RX J1340.6+4018 as well, contrary to earlier claims.     

Evidence of a link between the evolution of clusters and their AGN fraction

We discuss the optical properties, X-ray detections and active galactic nucleus (AGN) populations of four clusters at   z ∼ 1  in the Subaru–XMM Deep Field (SXDF). The velocity distribution and plausible extended X-ray detections are examined, as well as the number of X-ray point sources and radio sources associated with the clusters. We find that the two clusters that appear virialized and have an extended X-ray detection contain few, if any, AGN, whereas the two pre-virialized clusters have a large AGN population. This constitutes evidence that the AGN fraction in clusters is linked to the clusters' evolutionary stage. The number of X-ray AGN in the pre-virialized clusters is consistent with an overdensity of a factor of ∼200; the radio AGN appear to be clustered with a factor of 3 to 6 higher. The median K -band luminosities of   L_K = 1.7 ± 0.7 L *  for the X-ray sources and   L_K = 2.3 ± 0.1 L *  for the radio sources support the theory that these AGN are triggered by galaxy interaction and merging events in sub-groups with low internal velocity distributions, which make up the cluster environment in a pre-virialization evolutionary stage.     

Where do 'red and dead' early-type void galaxies come from?

Void regions of the Universe offer a special environment for studying cosmology and galaxy formation, which may expose weaknesses in our understanding of these phenomena. Although galaxies in voids are observed to be predominately gas rich, star forming and blue, a subpopulation of bright-red void galaxies can also be found, whose star formation was shutdown long ago. Are the same processes that quench star formation in denser regions of the Universe also at work in voids?
We compare the luminosity function of void galaxies in the 2dF Galaxy Redshift Survey, to those from a galaxy formation model built on the Millennium simulation. We show that a global star formation suppression mechanism in the form of low-luminosity 'radio-mode' active galactic nuclei (AGN) heating is sufficient to reproduce the observed population of void early types. Radio-mode heating is environment independent other than its dependence on dark matter halo mass, where, above a critical mass threshold of approximately   M _vir∼ 10^12.5 M_⊙  , gas cooling on to the galaxy is suppressed and star formation subsequently fades. In the Millennium simulation, the void halo mass function is shifted with respect to denser environments, but still maintains a high-mass tail above this critical threshold. In such void haloes, radio-mode heating remains efficient and red galaxies are found; collectively these galaxies match the observed space density without any modification to the model. Consequently, galaxies living in vastly different large-scale environments but hosted by haloes of similar mass are predicted to have similar properties, consistent with observations.     

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.     

Supersonic motions of galaxies in clusters

We study motions of galaxies in galaxy clusters formed in the concordance Λ cold dark matter cosmology. We use high-resolution cosmological simulations that follow the dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing the motions of galaxies and the properties of intracluster gas in a sample of eight simulated clusters at z = 0, we study the velocity dispersion profiles of the dark matter, gas and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ≈1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move slightly faster than the dark matter particles. The magnitude of the velocity bias,   b _v ≈ 1.1  , is, however, smaller than the bias estimated for subhaloes in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ≈20–30 per cent of the gas sound speed. These random motions provide about 10 per cent of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.     

The far-infrared–submillimetre spectral energy distribution of high-redshift quasars

We combine photometric observations of high-redshift     quasars, obtained at submillimetre to millimetre wavelengths, to obtain a mean far-infrared (rest-frame) spectral energy distribution (SED) of the thermal emission from dust, parametrized by a single temperature ( T ) and power-law emissivity index ( β ). The best-fitting values are     and     . Our method exploits the redshift spread of this set of quasars, which allows us to sample the SED at a larger number of rest wavelengths than is possible for a single object: the wavelength range extends down to ∼60 μm, and therefore samples the turnover in the greybody curve for these temperatures. This parametrization is of use to any studies that extrapolate from a flux at a single wavelength, for example to infer dust masses and far-infrared luminosities.
We interpret the cool, submillimetre component as arising from dust heated by star formation in the host galaxy of the quasar, although we do not exclude the presence of dust heated directly by the active galactic nucleus (AGN). Applying the mean SED to the data, we derive consistent star formation rates ∼1000 M_⊙ yr⁻¹ and dust masses ∼10⁹ M_⊙, and investigate a simple scheme of AGN and host galaxy co-evolution to account for these quantities. The time-scale for formation of the host galaxy is     , and the luminous quasar phase occurs towards the end of this period, just before the reservoir of cold gas is depleted. Given the youth of the Universe at     (1.6 Gyr), the coexistence of a massive black hole and a luminous starburst at high redshifts is a powerful constraint on models of quasar host galaxy formation.     

Enhanced mergers of galaxies in low-redshift clusters

An ensemble cluster has been formed from a data set comprising a complete magnitude-limited sample of 680 giant galaxies  ( M ⁰ _B ≲−19)  in eight low-redshift clusters, normalized by the velocity dispersions and virial radii for the early-type cluster populations. Distinct galaxy populations have been identified, including an infall population. A majority (50–70 per cent or greater) of the infall population are found to be in interacting or merging systems characterized by slow gravitational encounters. The observed enhancement of galaxy–galaxy encounters in the infall population compared to the field can be explained by gravitational shocking. It is shown that disc galaxy mergers in the infall population integrated over the estimated lifetime of the cluster (∼10 Gyr) can readily account for the present cluster S0 population.     

On the nature of EXO 0423.4−0840

We present a Chandra observation of the candidate BL Lac object EXO 0423.4−0840. The X-ray emission from EXO 0423.4−0840 is clearly extended, and is associated with an optical early-type galaxy, MCG-01-12-005, at the centre of cluster ClG 0422-09. We do not detect a point source that can be associated with a BL Lac, but we found a small radio source in the centre of MCG-01-12-005. The cluster gas temperature mapped by the Chandra observation drops continuously from 80 kpc towards the centre, and is locally single phase. We measure a metallicity profile that declines outwards with a value  0.8 Z_⊙  in the centre, dropping to  0.35 Z_⊙  at larger radius, which we interpret as a superposition of cluster gas and a dense interstellar medium (ISM) in the central galaxy. Although the temperature profile suggests that conduction is not efficient, the ISM and intra-cluster medium seem not to have mixed. The entropy profile declines continuously towards the centre, in agreement with recent results on groups and clusters. The radio source appears to have had some effect in terms of gas heating, as seen in the small-scale (∼10 kpc) entropy core, and the asymmetric hard emission on the same scale.     

The cluster environments of powerful radio-loud and radio-quiet active galactic nuclei

The spatial clustering amplitude ( B _gq) is determined for a sample of 44 powerful active galactic nuclei (AGN) at z ≃0.2. No significant difference is detected in the richness of the cluster environments of the radio-loud and radio-quiet subsamples, both of which typically inhabit environments as rich as Abell class ≃0. Comparison with radio luminosity-matched samples from Hill & Lilly and Wold et al. suggests that there is no epoch-dependent change in environment richness out to at least z ≥0.5 for either radio galaxies or radio quasars. Comparison with the APM cluster survey shows that, contrary to current folklore, powerful AGN do not avoid rich clusters, but rather display a spread in cluster environment, which is perfectly consistent with being drawn at random from the massive elliptical population. Finally, we argue that virtually all Abell class ≃0 clusters contained an active galaxy during the epoch of peak quasar activity at z ∼2.5.     

Radio study of X-ray clusters of galaxies — I. A3528: a pre-merging cluster in the core of the Shapley Supercluster

As part of an extensive radio–IR–optical–X-ray study of ROSAT clusters of galaxies in the Hydra region we have observed the bimodal Abell cluster A3528, located in the core of the Shapley Supercluster ( z  ≃ 0.053), with the Molonglo Observatory Synthesis Telescope at 843 MHz and the Australia Telescope Compact Array at 1.4 and 2.4 GHz. This is part I in a series of papers which looks at the relationship between the radio and X-ray emission in samples of ROSAT selected clusters.   The radio source characteristics — tailed morphologies and steep spectra — are consistent with the effects of a dense intracluster medium and the pre-merging environment of A3528. In particular, we present evidence that the minor member of the radio-loud dumbbell galaxy located at the centre of the northern component of A3528 is on a plunging orbit. We speculate that this orbit may have been induced by the tidal interactions between the merging components of A3528. In addition, the radio source associated with the dominant member of the dumbbell galaxy exhibits many of the characteristics of compact steep spectrum sources. We argue that the radio emission from this source was triggered ∼ 10⁶ yr ago by tidal interactions between the two members of the dumbbell galaxy, strengthening the argument that compact steep spectrum (CSS) sources are young.   Re-analysis of archive pointed Position Sensitive Proportional Counter (PSPC) data using multiresolution filtering suggests the presence of an AGN and/or a cooling flow in the southern component of A3528.     

Structural parameters for globular clusters in NGC 5128 – III. ACS surface brightness profiles and model fits

We present internal surface brightness profiles, based on Hubble Space Telescope /ACS imaging in the F 606 W bandpass, for 131 globular cluster (GC) candidates with luminosities   L ≃ 10⁴–3 × 10⁶ L_⊙  in the giant elliptical galaxy NGC 5128. Several structural models are fitted to the profile of each cluster and combined with mass-to-light ratios ( M / L values) from population-synthesis models, to derive a catalogue of fundamental structural and dynamical parameters parallel in form to the catalogues recently produced by McLaughlin & van der Marel and by Barmby et al. for GCs and massive young star clusters in Local Group galaxies. As part of this, we provide corrected and extended parameter estimates for another 18 clusters in NGC 5128, which we observed previously. We show that, like GCs in the Milky Way and some of its satellites, the majority of globulars in NGC 5128 are well fitted by isotropic Wilson models, which have intrinsically more distended envelope structures than the standard King lowered isothermal spheres. We use our models to predict internal velocity dispersions for every cluster in our sample. These predictions agree well in general with the observed dispersions in a small number of clusters for which spectroscopic data are available. In a subsequent paper, we use these results to investigate scaling relations for GCs in NGC 5128.