Explaining the energetic AGN outburst of MS 0735+7421 with massive slow jets

By conducting axisymmetrical hydrodynamical numerical simulations (2.5 dimensional code) we show that slow, massive, wide jets can reproduce the morphology of the huge X-ray deficient bubble pair in the cluster of galaxies MS 0735+7421. The total energy of the jets, composed of the energy in the bubble pair and in the shock wave, is constrained by observations conducted by McNamara et al. to be  ∼10⁶² erg  . We show that two opposite jets that are active for ∼100 Myr, each with a launching half opening angle of  α≃ 70°  , an initial velocity of   v _j∼ 0.1 c   and a total mass loss rate of the two jets of     , can account for the observed morphology. Rapidly precessing narrow jets can be used instead of wide jets. In our model the cluster suffered from a cooling catastrophe ∼100 Myr ago. Most of the mass that cooled,  ∼10¹⁰ M_⊙  , was expelled back to the intracluster medium by the active galactic nuclei activity and is inside the bubbles now, ∼10 per cent formed stars and ∼10 per cent of the cold gas was accreted by the central black hole and was the source of the outburst energy. This type of activity is similar to that expected to occur in galaxy formation.     

A diffuse bubble-like radio-halo source MRC 0116+111: imprint of AGN feedback in a low-mass cluster of galaxies

We present detailed observations of MRC 0116+111, revealing a luminous, miniradio halo of ∼240-kpc diameter located at the centre of a cluster of galaxies at redshift   z = 0.131  . Our optical and multiwavelength Giant Metrewave Radio Telescope and Very Large Array radio observations reveal a highly unusual radio source: showing a pair of giant (∼100-kpc diameter) bubble-like diffuse structures, that are about three times larger than the analogous extended radio emission observed in M87 – the dominant central radio galaxy in the Virgo cluster. However, in MRC 0116+111 we do not detect any ongoing active galactic nucleus (AGN) activity, such as a compact core or active radio jets feeding the plasma bubbles. The radio emitting relativistic particles and magnetic fields were probably seeded in the past by a pair of radio jets originating in the AGN of the central cD galaxy. The extremely steep high-frequency radio spectrum of the north-western bubble, located ∼100 kpc from cluster centre, indicates radiation losses, possibly because having detached, it is rising buoyantly and moving away into the putative hot intracluster medium. The other bubble, closer to the cluster centre, shows signs of ongoing particle re-acceleration. We estimate that the radio jets which inflated these two bubbles might have also fed enough energy into the intracluster medium to create an enormous system of cavities and shock fronts, and to drive a massive outflow from the AGN, which could counter-balance and even quench a cooling flow. Therefore, this source presents an excellent opportunity to understand the energetics and the dynamical evolution of radio jet inflated plasma bubbles in the hot cluster atmosphere.     

On viscosity, conduction and sound waves in the intracluster medium

Recent X-ray and optical observations of the Perseus cluster indicate that a combination of weak shocks at small radii  (≳20  kpc)  and viscous and conductive dissipation of sound waves at larger radii is responsible for heating the intracluster medium and can balance radiative cooling of cluster cores. We discuss this mechanism more generally and show how the specific heating and cooling rates vary with temperature and radius. It appears that this heating mechanism is most effective above  10⁷  K  , which allows for radiative cooling to proceed within normal galaxy formation but stifles the growth of very massive galaxies. The scaling of the wavelength of sound waves with cluster temperature and feedback in the system are investigated.     

The weak shock in the core of the Perseus cluster

The dissipation of energy from sound waves and weak shocks is one of the most promising mechanisms for coupling active galactic nucleus (AGN) activity to the surrounding intracluster medium, and so offsetting cooling in cluster cores. We present a detailed analysis of the weak shock found in deep Chandra observations of the Perseus cluster core. A comparison of the spectra either side of the shock front shows that they are very similar. By performing a deprojection analysis of a sector containing the shock, we produce temperature and density profiles across the shock front. These show no evidence for a temperature jump coincident with the density jump. To understand this result, we model the shock formation using 1D hydrodynamic simulations including models with thermal conduction and  γ < 5/3  gas. These models do not agree well with the data, suggesting that further physics is needed to explain the shock structure. We suggest that an interaction between the shock and the Hα filaments could have a significant effect on cooling the post-shock gas.
We also calculate the thermal energy liberated by the weak shock. The total energy in the shocked region is about 3.5 times the work needed to inflate the bubbles adiabatically, and the power of the shock is around  6 × 10⁴⁴ erg s⁻¹  per bubble, just over  10⁴⁵ erg s⁻¹  in total.     

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.     

Exciting molecular hydrogen in the central galaxies of cooling flows

The origin of rovibrational H₂ emission in the central galaxies of cooling flow clusters is poorly understood. Here we address this issue using data from our near-infrared spectroscopic survey of 32 of the most line-luminous such systems, presented in the companion paper by Edge et al.
We consider excitation by X-rays from the surrounding intracluster medium (ICM), ultra-violet (UV) radiation from young stars, and shocks. The   v = 1–0  K -band lines with upper levels within  10⁴ K  of the ground state appear to be mostly thermalized (implying gas densities  ≳10⁵ cm⁻³  ), with the excitation temperature typically exceeding 2000 K, as found earlier by Jaffe, Bremer & van der Werf. Together with the lack of strong   v = 2–0  lines in the H -band, this rules out UV radiative fluorescence.
Using the cloudy photoionization code, we deduce that the H₂ lines can originate in a population of dense clouds, exposed to the same hot  ( T ∼ 50 000 K)  stellar continuum as the lower density gas which produces the bulk of the forbidden optical line emission in the Hα-luminous systems. This dense gas may be in the form of self-gravitating clouds deposited directly by the cooling flow, or may instead be produced in the high-pressure zones behind strong shocks. Furthermore, the shocked gas is likely to be gravitationally unstable, so collisions between the larger clouds may lead to the formation of globular clusters.     

Tracing gas motions in the Centaurus cluster

We apply the stochastic model of iron transport developed by Rebusco et al. to the Centaurus cluster. Using this model, we find that an effective diffusion coefficient D in the range  2 × 10²⁸–4 × 10²⁸ cm² s⁻¹  can approximately reproduce the observed abundance distribution. Reproducing the flat central profile and sharp drop around  30–70 kpc  , however, requires a diffusion coefficient that drops rapidly with radius so that   D > 4 × 10²⁸ cm² s⁻¹  only inside about  25 kpc  . Assuming that all transport is due to fully developed turbulence, which is also responsible for offsetting cooling in the cluster core, we calculate the length- and velocity-scales of energy injection. These length-scales are found to be up to a factor of ∼10 larger than expected if the turbulence is due to the inflation and rising of a bubble. We also calculate the turbulent thermal conductivity and find it is unlikely to be significant in preventing cooling.     

Global velocity field and bubbles in the blue compact dwarf galaxy Mrk 86

We have studied the velocity field of the blue compact dwarf galaxy Mrk 86 (NGC 2537) using data provided by 14 long-slit optical spectra obtained in 10 different orientations and positions. This kinematical information is complemented with narrow-band ([O  iii ]5007 Å and H α ) and broad-band ( B , V , Gunn r and K ) imaging. The analysis of the galaxy global velocity field suggests that the ionized gas could be distributed in a rotating inclined disc, with projected central angular velocity of Ω=34 km s⁻¹ kpc⁻¹. The comparison between the stellar, H  i and modelled dark matter density profile indicates that the total mass within its optical radius is dominated by the stellar component. Peculiarities observed in its velocity field can be explained by irregularities in the ionized gas distribution or local motions induced by star formation.
Kinematical evidences for two expanding bubbles, Mrk 86–B and Mrk 86–C, are given. They show expanding velocities of 34 and 17 km s⁻¹, H α luminosities of 3×10³⁸ and 1.7×10³⁹ erg s⁻¹, and physical radii of 374 and 120 pc, respectively. The change in the [S  ii ]/H α , [N  ii ]/H α , [O  ii ]/[O  iii ] and [O  iii ]/H β line ratios with the distance to the bubble precursor suggests a diminution in the ionization parameter and, in the case of Mrk 86–B, an enhancement of the shock-excited gas emission. The optical–near-infrared colours of the bubble precursors are characteristic of low‐metallicity star‐forming regions (∼0.2 Z_⊙) with burst strengths of about 1 per cent in mass.     

Starbursts and the triggering of the activity in nearby powerful radio galaxies

We present high-quality long-slit spectra for three nearby powerful radio galaxies – 3C 293, 3C 305 and PKS 1345+12. These were taken with the aim of characterizing the young stellar populations (YSP), and thereby investigating the evolution of the host galaxies, as well as the events that triggered the activity. Isochrone spectral synthesis modelling of the wide wavelength coverage spectra of nuclear and off-nuclear continuum-emitting regions have been used to estimate the ages, masses and luminosities of the YSP component, taking full account of reddening effects and potential contamination by activity-related components. We find that the YSP make a substantial contribution to the continuum flux in the off-nuclear regions on a radial scale of 1–20 kpc in all three objects. Moreover, in two objects we find evidence for reddened post-starburst stellar populations in the near-nuclear regions of the host galaxies. The YSP are relatively old (0.1–2 Gyr), massive  (10⁹ < M _YSP < 2 × 10¹⁰ M_⊙)  and make up a large proportion (∼1–50 per cent) of the total stellar mass in the regions of the galaxies sampled by the observations. Overall, these results are consistent with the idea that the nuclear activity of active galactic nuclei in some radio galaxies is triggered by major gas-rich mergers. Therefore, these radio galaxies form part of the subset of early-type galaxies that is evolving most rapidly in the local Universe. Intriguingly, the results also suggest that the radio jets are triggered relatively late in the merger sequence, and that there is an evolutionary link between radio galaxies and luminous/ultraluminous infrared galaxies.     

Impact of cosmic rays on Population III star formation

We explore the implications of a possible cosmic-ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of   z ∼ 15–20  , provided the CR flux was sufficient to yield an ionization rate greater than about 10⁻¹⁹ s⁻¹ near the centre of the minihalo. The presence of CRs with energies  ≲10⁷  eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as  ∼10 M_⊙  . CRs have a less-pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses  ≳10⁸ M_⊙  at the later stages of cosmological structure formation around   z ∼ 10–15  . In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.     

X-ray analysis of the cluster 3C 129

The cluster 3C 129 is classified as a rich cluster. An analysis of the properties of the cluster 3C 129 from ROSAT PSPC and HRI, Einstein IPC, and EXOSAT ME observations is presented. The mean temperature from a joint fit of the ROSAT PSPC and EXOSAT ME data is 5.5(±0.2) keV. The luminosity is 0.6×10⁴⁴ erg s⁻¹ in 0.2–2.4 keV and 2.7×10⁴⁴ erg s⁻¹ in 0.2–10 keV. We find a cooling flow with a rate of ∼84 M_⊙ yr⁻¹. The central gas density is 6×10⁻³ cm⁻³, and the ICM mass is 3.6×10¹³ M_⊙. The total cluster mass is ∼5×10¹⁴ M_⊙. The X-ray morphology shows an east–west elongation, which is evidence for a recent merger event. The radio source 3C 129.1 is located near the X-ray centre. Another cluster member galaxy (the radio galaxy 3C 129) is a prototype of head-tailed radio galaxies, and is located in the west part of the cluster. The tail points along the gradient of intracluster gas pressure. There are no significant point X-ray sources associated with the AGNs of the two radio galaxies.     

Active galactic nuclei jet mass loading and truncation by stellar winds

Active galactic nuclei can produce extremely powerful jets. While tightly collimated, the scale of these jets and the stellar density at galactic centres implies that there will be many jet/star interactions, which can mass load the jet through stellar winds. Previous work employed modest wind mass outflow rates, but this does not apply when mass loading is provided by a small number of high mass-loss stars. We construct a framework for jet mass loading by stellar winds for a broader spectrum of wind mass-loss rates than has previously been considered. Given the observed stellar mass distributions in galactic centres, we find that even highly efficient (0.1 Eddington luminosity) jets from supermassive black holes of masses M _BH≲ 10⁴ M_⊙ are rapidly mass loaded and quenched by stellar winds. For  10⁴ M_⊙ < M _BH < 10⁸ M_⊙  , the quenching length of highly efficient jets is independent of the jet's mechanical luminosity. Stellar wind mass loading is unable to quench efficient jets from more massive engines, but can account for the observed truncation of the inefficient M87 jet, and implies a baryon-dominated composition on scales ≳2 kpc therein even if the jet is initially pair plasma dominated.     

Further constraints on the evolution of Ks-selected galaxies in the GOODS/CDFS field

We have selected and analysed the properties of a sample of  2905 K_s < 21.5  galaxies in  ∼131 arcmin²  of the Great Observatories Origins Deep Survey (GOODS) Chandra Deep Field South (CDFS), to obtain further constraints on the evolution of K_s -selected galaxies with respect to the results already obtained in previous studies. We made use of the public deep multiwavelength imaging from the optical B through the infrared (IR) 4.5-μm bands, in conjunction with available spectroscopic and COMBO17 data in the CDFS, to construct an optimized redshift catalogue for our galaxy sample. We computed the K_s -band luminosity function and determined that its characteristic magnitude has a substantial brightening and a decreasing total density from   z = 0  to  〈 z 〉= 2.5  . We also analysed the colours and number density evolution of galaxies with different stellar masses. Within our sample, and in contrast to what is observed for less massive systems, the vast majority (∼85–90 per cent) of the most massive  ( M > 2.5 × 10¹¹ M_⊙)  local galaxies appear to be in place before redshift   z ∼ 1  . Around 65–70 per cent of the total assemble between redshifts   z = 1  and 3 and most of them display extremely red colours, suggesting that plausible star formation in these very massive systems should mainly proceed in obscured, short-time-scale bursts. The remaining fraction (up to ∼20 per cent) could be in place at even higher redshifts   z = 3–4  , pushing the first epoch of formation of massive galaxies beyond the limits of current near-IR surveys.     

Chandra measurements of the X-ray core and cluster of 3C 220.1

We report results of an 18-ks exposure with the ACIS instrument on Chandra of the powerful z =0.62 radio galaxy 3C 220.1. The X-ray emission separates into cluster gas of emission-weighted kT ∼5 keV , 0.7–12 keV luminosity (to a radius of 45 arcsec) of 5.6×10⁴⁴ erg s⁻¹ and unresolved emission (coincident with the radio core). While the extended X-ray emission is clearly thermal in nature, a straightforward cooling-flow model, even in conjunction with a point-source component, is a poor fit to the radial profile of the X-ray emission. This is despite the fact that the measured properties of the gas suggest a massive cooling flow of ∼130 M_⊙ yr⁻¹, and the data show weak evidence for a temperature gradient. The central unresolved X-ray emission has a power-law spectral energy index α ∼0.7 and 0.7–12 keV luminosity of 10⁴⁵ erg s⁻¹, and any intrinsic absorption is relatively small. The two-point spectrum of the core emission between radio and X-ray energies has α _rx=0.75 . Since this is a flatter spectrum than seen in other sources where the X-ray emission is presumed to be radio-related, regions close to the active galactic nucleus (AGN) in this source may dominate the central X-ray output, as is believed to be the case for lobe-dominated quasars. Simple unification models would be challenged if this were found to be the case for a large fraction of high-power radio galaxies.     

On the formation of massive galaxies: a simultaneous study of number density, size and intrinsic colour evolution in GOODS

The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the Hubble Space Telescope /Advanced Camera for Surveys images of the Great Observatories Origins Deep Survey (GOODS) North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin² with stellar masses above  3 × 10¹⁰ M_⊙  in the redshift range  0.4 < z < 1.2  . Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems  (≳3 × 10¹¹ M_⊙)  do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dF galaxy redshift survey, we find that the number density of massive early-type galaxies is consistent with no evolution between   z = 1.2  and 0, i.e. over an epoch spanning more than half of the current age of the Universe. We find large discrepancies between the predictions of semi-analytic models. Massive galaxies show very homogeneous intrinsic colour distributions, with nearly flat radial colour gradients, but with a significant negative correlation between stellar mass and colour gradient, such that red cores appear predominantly in massive galaxies. The distribution of half-light radii – when compared to   z ∼ 0  and   z > 1  samples – is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers.     

Bubbles, feedback and the intracluster medium: three-dimensional hydrodynamic simulations

We use a three-dimensional hydrodynamical code to simulate the effect of energy injection on cooling flows in the intracluster medium. Specifically, we compare a simulation of a 10¹⁵-M_⊙ cluster with radiative cooling only with a second simulation in which thermal energy is injected 31 kpc off-centre, over 64 kpc³ at a rate of     for 50 Myr. The heat injection forms a hot, low-density bubble which quickly rises, dragging behind it material from the cluster core. The rising bubble pushes with it a shell of gas which expands and cools. We find the appearance of the bubble in X-ray temperature and luminosity to be in good qualitative agreement with recent Chandra observations of cluster cores. Toward the end of the simulation, at 600 Myr, the displaced gas begins to fall back toward the core, and the subsequent turbulence is very efficient at mixing the low- and high-entropy gas. The result is that the cooling flow is disrupted for up to ∼ 50 Myr after the injection of energy ceases. Thus this mechanism provides a very efficient method for regulating cooling flows, if the injection events occur with a 1:1 duty cycle.     

The growth of central and satellite galaxies in cosmological smoothed particle hydrodynamics simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to  7 × 10⁹ M_⊙  . Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2–5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5–1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1–0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since   z = 1  , 27 per cent of central galaxies (above  3 × 10¹⁰ M_⊙  ) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain 'central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.     

Chandra finds that X-ray jets are common in low-power radio galaxies

We present results for the first three low-power radio galaxies from the B2 bright sample to have been observed with Chandra . Two have kiloparsec-scale radio jets, and in both Chandra resolves jet X-ray emission, and detects soft X-ray core emission and an X-ray-emitting galaxy-scale atmosphere of luminosity a few ×10⁴¹ erg s⁻¹. These are the first detections of X-ray jets in low-power radio galaxies more distant than Centaurus A and M87. The cooling time of the galaxy-scale gas implies mass infall rates of the order of 1 M_⊙ yr⁻¹. The gas pressure near the jets is comparable to the minimum pressure in the jets, implying that the X-ray-emitting gas may play an important role in jet dynamics. The third B2 radio galaxy has no kiloparsec-scale radio jet, and here only soft X-ray emission from the core is detected. The ratio of X-ray to radio flux is similar for the jets and cores, and the results favour a synchrotron origin for the emission. Kiloparsec-scale radio jets are detected in the X-ray in ∼7-ks exposures with Chandra more readily than in the optical via Hubble Space Telescope snapshot surveys.     

Dry mergers: a crucial test for galaxy formation

We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of  ∼6 × 10⁻⁵ Mpc⁻³ Gyr⁻¹  at   z ≤ 1  and a steep decline at larger z . Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of   M _{*, c} ∼ 6.3 × 10¹⁰ M_⊙  by mergers only. As a consequence of this transition, we find that around   M _{*, c}   , the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z . We however find that only 10–20 per cent of galaxies more massive than   M _{*, c}   experience dry major mergers within their last Gyr at any given redshift   z ≤ 1  .     

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  .