Sound waves excitation by jet-inflated bubbles in clusters of galaxies

We show that repeated sound waves in the intracluster medium (ICM) can be excited by a single inflation episode of an opposite bubble pair. To reproduce this behaviour in numerical simulations, the bubbles should be inflated by jets, rather than being injected artificially as already full-blown bubbles. The multiple sound waves are excited by the motion of the bubble–ICM boundary that is caused by vortices inside the inflated bubbles and the backflow ('cocoon') of the ICM around the bubble. These sound waves form a structure that can account for the ripples observed in the Perseus cooling flow cluster. We inflate the bubbles using slow massive jets either with a very wide opening angle or that are narrow and precessing. The wide jets (or collimated fast winds) are slow in the sense that they are highly subrelativistic,   v _j∼ 0.01 c – 0.1 c   , and they are massive in the sense that the pair of bubbles carries back to the ICM a large fraction of the cooling mass, i.e.  ∼1–50 M_⊙ yr⁻¹  . We use a two-dimensional axisymmetric (referred to as 2.5D) hydrodynamical numerical code ( vh-1 ).     

The hydrodynamics of dead radio galaxies

We present a numerical investigation of dead, or relic, radio galaxies and the environmental impact that radio galaxy activity has on the host galaxy or galaxy cluster. We perform axisymmetric hydrodynamical calculations of light, supersonic, back-to-back jets propagating in a β -model galaxy/cluster atmosphere. We then shut down the jet activity and let the resulting structure evolve passively. The dead source undergoes an initial phase of pressure driven expansion until it achieves pressure equilibrium with its surroundings. Thereafter, buoyancy forces drive the evolution and lead to the formation of two oppositely directed plumes that float high into the galaxy/cluster atmosphere. These plumes entrain a significant amount of low entropy material from the galaxy/cluster core and lift it high into the atmosphere. An important result is that a large fraction (at least half) of the energy injected by the jet activity is thermalized in the interstellar medium (ISM)/intracluster medium (ICM) core. The whole ISM/ICM atmosphere inflates in order to regain hydrostatic equilibrium. This inflation is mediated by an approximately spherical disturbance which propagates into the atmosphere at the sound speed. The fact that such a large fraction of the injected energy is thermalized suggests that radio galaxies may have an important role in the overall energy budget of rich ISM/ICM atmospheres. In particular, they may act as a strong and highly time-dependent source of negative feedback for galaxy/cluster cooling flows.     

Stability of magnetic equilibria in radio bubbles

Current-carrying flows, in the laboratory and in astrophysical jets, can form remarkably stable magnetic structures. Decades of experience show that such flows often build equilibria that reverse field directions, evolving to a magnetohydrodynamic (MHD) Taylor state, which has remarkable stability properties. We model jets and the magnetic bubbles they build as reversed-field pinch equilibria by assuming the driver current to be stiff in the MHD sense. Taking the jet current as rigid and a fixed function of position, we prove a theorem: that the same, simple MHD stability conditions guarantee stability, even after the jet turns off. This means that magnetic structures harbouring a massive inventory of magnetic energy can persist long after the building jet current has died away. These may be the relic radio 'fossils', 'ghost bubbles' or 'magnetic balloons' found in clusters. These equilibria, which are under magnetic tension, will evolve, retaining the stability properties from that state. The remaining fossil is not a disordered ball of magnetic fields, but a stable structure under tension, able to respond to the slings and arrows of outside forces. Typically their Alfvén speeds greatly exceed the cluster sound speed, and so they can keep out hot cluster plasma, leading to X-ray ghosts. Passing shocks cannot easily destroy them, but can energize and light them up anew at radio frequencies. Bubbles can rise in the hot cluster plasma, perhaps detaching from the parent radio galaxy but stable against Rayleigh–Taylor and other modes.     

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.     

Buoyant radio plasma in clusters of galaxies

Radio galaxies are known to inflate lobes of hot relativistic plasmas into the intergalactic medium. Here we present hydrodynamical and magnetohydrodynamical simulations of these hot plasma bubbles in FR II objects. We focus on the later stages of their evolution after the jet has died down and after the bow shock that surrounded the lobes at earlier stages has vanished. We investigate the evolution of the plasma bubbles as they become subject to Rayleigh–Taylor instabilities. From our simulations we calculate the radio and X-ray emissivities of the bubbles and discuss their appearance in observations. Finally, we investigate the influence of large-scale magnetic fields on the evolution of the bubbles. The issues of re-acceleration and diffusion of relativistic particles are briefly discussed.     

The optical/X-ray connection: intra-cluster medium iron content and galaxy optical luminosity in 20 galaxy clusters

X-ray observations of galaxy clusters have shown that the intra-cluster gas has iron abundances of about one-third of the solar value. These observations also show that part (if not all) of the intra-cluster gas metals was produced within the member galaxies. We present a systematic analysis of 20 galaxy clusters to explore the connection between the iron mass and the total luminosity of early- and late-type galaxies, and of the brightest cluster galaxies (BCGs). From our results, the intra-cluster medium (ICM) iron mass seems to correlate better with the luminosity of the BCGs than with that of the red and blue galaxy populations. As the BCGs cannot produce alone the observed amount of iron, we suggest that ram-pressure plus tidal stripping acts together to enhance, at the same time, the BCG luminosities and the iron mass in the ICM. Through the analysis of the iron yield, we have also estimated that SN Ia are responsible for more than 50 per cent of the total iron in the ICM. This result corroborates the fact that ram-pressure contributes to the gas removal from galaxies to the ICM, being very efficient for clusters in the temperature range  2 < kT (keV) < 10  .     

Galaxies in clusters: the observational characteristics of bow shocks, wakes and tails

The dynamical signatures of the interaction between galaxies in clusters and the intracluster medium (ICM) can potentially yield significant information about the structure and dynamical history of clusters. To develop our understanding of this phenomenon we present results from numerical modelling of the galaxy–ICM interaction, as the galaxy moves through the cluster. The simulations have been performed for a broad range of ICM temperatures ( kT _cl=1, 4 and 8 keV), representative of poor clusters or groups through to rich clusters.
There are several dynamical features that can be identified in these simulations. For supersonic galaxy motion, a leading bow shock is present, and also a weak gravitationally focused wake or tail behind the galaxy (analogous to Bondi–Hoyle accretion). For galaxies with higher mass replenishment rates and a denser interstellar medium (ISM), the dominant feature is a dense ram-pressure stripped tail. In line with other simulations, we find that the ICM/galaxy–ISM interaction can result in complex time-dependent dynamics, with ram-pressure stripping occurring in an episodic manner.
In order to facilitate this comparison between the observational consequences of dynamical studies and X-ray observations we have calculated synthetic X-ray flux and hardness maps from these simulations. These calculations predict that the ram-pressure stripped tail will usually be the most visible feature, though in nearby galaxies the bow shock preceding the galaxy should also be apparent in deeper X-ray observations. We briefly discuss these results and compare them with X-ray observations of galaxies where there is evidence of such interactions.     

Morphology of flows and buoyant bubbles in the Virgo cluster

There is growing evidence that the active galactic nuclei (AGN) associated with the central elliptical galaxy in clusters of galaxies are playing an important role in the evolution of the intracluster medium (ICM) and clusters themselves. We use high-resolution three-dimensional simulations to study the interaction of the cavities created by AGN outflows (bubbles) with the ambient ICM. The gravitational potential of the cluster is modelled using the observed temperature and density profiles of the Virgo cluster. We demonstrate the importance of the hydrodynamical Kutta–Zhukovsky forces associated with the vortex ring structure of the bubbles, and discuss possible effects of diffusive processes on their evolution.     

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.     

Temperature and abundance profiles of hot gas in galaxy groups – II. Implications for feedback and ICM enrichment

We investigate the history of galactic feedback and chemical enrichment within a sample of 15 X-ray bright groups of galaxies, on the basis of the inferred Fe and Si distributions in the hot gas and the associated metal masses produced by core-collapse and Type Ia supernovae (SNe). Most of these cool-core groups show a central Fe and Si excess, which can be explained by prolonged enrichment by SN Ia and stellar winds in the central early-type galaxy alone, but with tentative evidence for additional processes contributing to core enrichment in hotter groups. Inferred metal mass-to-light ratios inside r ₅₀₀ show a positive correlation with total group mass but are generally significantly lower than in clusters, due to a combination of lower global intracluster medium (ICM) abundances and gas-to-light ratios in groups. This metal deficiency is present for products from both SN Ia and SN II, and suggests that metals were either synthesized, released from galaxies or retained within the ICM less efficiently in lower mass systems. We explore possible causes, including variations in galaxy formation and metal release efficiency, cooling out of metals, and gas and metal loss via active galactic nuclei (AGN) – or starburst-driven galactic winds from groups or their precursor filaments. Loss of enriched material from filaments coupled with post-collapse AGN feedback emerges as viable explanations, but we also find evidence for metals to have been released less efficiently from galaxies in cooler groups and for the ICM in these to appear chemically less evolved, possibly reflecting more extended star formation histories in less massive systems. Some implications for the hierarchical growth of clusters from groups are briefly discussed.     

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  .     

Interface modes and their instabilities in accretion disc boundary layers

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the innermost regions of accretion discs or boundary layers. Two simple models are considered in this paper. The magnetosphere–disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star–disc model involves a Keplerian disc terminated at the stellar atmosphere with high density and small density scaleheight. We find that the interface modes at the magnetosphere–disc boundary are generally unstable due to Rayleigh–Taylor and/or Kelvin–Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh–Taylor instability, and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin–Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh–Taylor or Kelvin–Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.     

Interactions of radio galaxies and the intracluster medium in Abell 160 and Abell 2462

We present Chandra and Very Large Array observations of two galaxy clusters, Abell 160 and Abell 2462, whose brightest cluster galaxies (BCGs) host wide angle tailed radio galaxies (WATs). We search for evidence of interactions between the radio emission and the hot, X-ray emitting gas, and we test various jet termination models. We find that both clusters have cool BCGs at the cluster centre, and that the scale of these cores (∼30–40 kpc for both sources) is of approximately the same scale as the length of the radio jets. For both sources, the jet flaring point is coincident with a steepening in the host cluster's temperature gradient, and similar results are found for 3C 465 and Hydra A. However, none of the published models of WAT formation offers a satisfactory explanation as to why this may be the case. Therefore, it is unclear what causes the sudden transition between the jet and the plume. Without accurate modelling, we cannot ascertain whether the steepening of the temperature gradient is the main cause of the transition, or merely a tracer of an underlying process.     

Hydrodynamical simulations of the Sunyaev–Zel'dovich effect: the kinetic effect

We use hydrodynamical N -body simulations to study the kinetic Sunyaev–Zel'dovich effect. We construct sets of maps, one square degree in size, in three different cosmological models. We confirm earlier calculations that on the scales studied the kinetic effect is much smaller than the thermal (except close to the thermal null point), with an rms dispersion smaller by about a factor of 5 in the Rayleigh–Jeans region. We study the redshift dependence of the rms distortion and the pixel distribution at the present epoch. We compute the angular power spectra of the maps, including their redshift dependence, and compare them with the thermal Sunyaev–Zel'dovich effect and with the expected cosmic microwave background anisotropy spectrum as well as with determinations by other authors. We correlate the kinetic effect with the thermal effect both pixel-by-pixel and for identified thermal sources in the maps to assess the extent to which the kinetic effect is enhanced in locations of strong thermal signal.     

Cluster–subcluster mergers and the formation of narrow-angle tailed radio sources

We examine the ROSAT PSPC X-ray properties of a sample of 15 Abell clusters containing 23 narrow-angle tailed (NAT) radio galaxies. We find that clusters with NATs show a significantly higher level of substructure than a similar sample of radio-quiet clusters, indicating that NAT radio sources are preferentially located in dynamically complex systems. Also, the velocity distribution of the NAT galaxies is similar to that of other cluster members; these velocities are inadequate for producing the ram pressure necessary to bend the radio jets. We therefore propose a new model for NAT formation, in which NATs are associated with dynamically complex clusters undergoing merger events. The U -shaped NAT morphology is produced in part by the merger-induced bulk motion of the ICM bending the jets.     

Chemical evolution of Elliptical Galaxies and the ICM

We present a new model for the chemical evolution of elliptical galaxies taking into account SN feedback, detailed nucleosynthesis and galactic winds. We discuss the effect of galactic winds on the chemical enrichment of the ICM and compute the energy per particle injected by the galaxies into the ICM. This revised version was published online in August 2006 with corrections to the Cover Date.     

Galaxy cluster masses without non-baryonic dark matter

We apply the modified acceleration law obtained from Einstein gravity coupled to a massive skew symmetric field,   F _μνλ  , to the problem of explaining X-ray galaxy cluster masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot intracluster medium (ICM) with King 'β-models', we show that the dynamical masses of the galaxy clusters resulting from our modified acceleration law fit the cluster gas masses for our sample of 106 clusters without the need of introducing a non-baryonic dark matter component. We are further able to show for our sample of 106 clusters that the distribution of gas in the ICM as a function of radial distance is well fitted by the dynamical mass distribution arising from our modified acceleration law without any additional dark matter component. In a previous work, we applied this theory to galaxy rotation curves and demonstrated good fits to our sample of 101 low surface brightness, high surface brightness and dwarf galaxies including 58 galaxies that were fitted photometrically with the single-parameter mass-to-light ratio ( M / L )_stars. The results obtained there were qualitatively similar to those obtained using Milgrom's phenomenological Modified Newtonian Dynamics (MOND) model, although the determined galaxy masses were quantitatively different, and MOND does not show a return to Keplerian behaviour at extragalactic distances. The results obtained here are compared to those obtained using Milgrom's phenomenological MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary dark matter component is included.     

Wakes of ram-pressure-stripped disc galaxies

We present a sample of 30 wide-angle tailed radio galaxies (WATs) that we use to constrain the jet speeds in these sources. We measure the distribution of the jet-sidedness ratios for the sample, and assuming that the jets are beamed, the jet speeds in the range (0.3–0.7) c are obtained. Whilst the core prominence of the sample, which ought to be a reliable indicator of beaming, shows little correlation with the jet sidedness, we argue that due to the peculiar nature of the WATs, core prominence is unlikely to be a good indicator of beaming in these sources. We further show that if the jets are fast and light, then the galaxy speeds required to bend the jets into C-shapes such as those seen in 0647+693 are reasonable for a galaxy in a merging or recently merged cluster.