Hydrodynamics and stability of galactic cooling flows

Using numerical techniques we study the global stability of cooling flows in X-ray luminous giant elliptical galaxies. As an unperturbed equilibrium state we choose the hydrostatic gas recycling model. Non-equilibrium radiative cooling, stellar mass loss, heating by type Ia supernovae, distributed mass deposition and thermal conductivity are included. Although the recycling model reproduces the basic X-ray observables, it appears to be unstable with respect to the development of inflow or outflow. In spherical symmetry the inflows are subject to a central cooling catastrophe, while the outflows saturate in a form of a subsonic galactic wind. Two-dimensional axisymmetric random velocity perturbations of the equilibrium model trigger the onset of a cooling catastrophe, which develops in an essentially non-spherical way. The simulations show a patchy pattern of mass deposition and the formation of hollow gas jets, which penetrate through the outflow down to the galaxy core. The X-ray observables of such a hybrid gas flow mimic those of the equilibrium recycling model, but the gas temperature exhibits a central depression. The mass deposition rate M ˙ consists of two contributions of similar size: (i) a hydrostatic one resembling that of the equilibrium model, and (ii) a dynamical one which is related to the jets and is more concentrated towards the centre. For a model galaxy, like NGC 4472, our 2D simulations predict M ˙ ≈ 2 M⊙ yr⁻¹ within the cooling radius for the advanced non-linear stage of the instability. We discuss the implications of these results to Hα nebulae and star formation in cooling flow galaxies and emphasize the need for high-resolution 3D simulations.     

Formation and evolution of elliptical galaxies and QSO activity

We present the results of a numerical code that combines multi-zone chemical evolution with 1D hydrodynamics to follow in detail the evolution and radial behaviour of gas and stars during the formation of elliptical galaxies. We use the model to explore the links between the evolution and formation of elliptical galaxies and QSO activity. The knowledge of the radial gas flows in the galaxy allows us to trace metallicity gradients, and, in particular, the formation of a high-metallicity core in ellipticals. The high-metallicity core is formed soon enough to explain the metal abundances inferred in high-redshift quasars. The star formation rate and the subsequent feedback regulate the episodes of wind, outflow and cooling flow, thus affecting the recycling of the gas and the chemical enrichment of the intergalactic medium. The evolution of the galaxy shows several stages, some of which are characterized by a complex flow pattern, with inflow in some regions and outflow in other regions. All models, however, exhibit during their late evolution a galactic wind at the outer boundary and, during their early evolution, an inflow towards the galactic nucleus. The characteristics of the inner inflow could explain the bolometric luminosity of a quasar lodged at the galactic centre as well as the evolution of the optical luminosity of quasars.     

Simulations of mixed-morphology supernova remnants with anisotropic thermal conduction

We explore the role of anisotropic thermal conduction on the evolution of supernova remnants (SNRs) through interstellar media with a range of densities via numerical simulations. We find that a remnant expanding in a dense environment can produce centre-bright hard X-ray emission within 20 kyr, and centre-bright soft X-ray emission within 60 kyr of the supernova event. In a more tenuous environment, the appearance of a centre-bright structure in hard X-rays is delayed until about 60 kyr. The soft X-ray emission from such a remnant may not become centre bright during its observable lifetime. This can explain the observations that show that mixed-morphology SNRs preferentially occur close to denser, molecular environments. Remnants expanding into denser environments tend to be smaller, making it easier for thermal conduction to make large changes in the temperatures of their hot gas bubbles. We show that the lower temperatures make it very favourable to use high-stage ions as diagnostics of the hot gas bubbles in SNRs. In particular, the distribution of O  viii transitions from shell bright at early epochs to centre bright at later epochs in the evolution of an SNR expanding in a dense interstellar medium when the physics of thermal conduction is included.     

The trouble with the Local Bubble

The model of a Local Hot Bubble has been widely accepted as providing a framework that can explain the ubiquitous presence of the soft X-ray background diffuse emission. We summarize the current knowledge on this local interstellar region, paying particular reference to observations that sample emission from the presumed local million degree K hot plasma. However, we have listed numerous observations that are seemingly in conflict with the concept of a hot Local Bubble. In particular, the discovery of solar wind charge exchange that can generate an appreciable soft X-ray background signal within the heliosphere, has led to a re-assessment of the generally accepted model that requires a hot local plasma. In order to explain the majority of observations of the local plasma, we forward two new speculative models that describe the physical state of the local interstellar gas. One possible scenario is similar to the present widely accepted model of the Local Hot Bubble, except that it accounts for only 50% of the soft X-ray emission currently detected in the galactic plane, has a lower thermal pressure than previously thought, and its hot plasma is not as hot as previously believed. Although such a model can solve several difficulties with the traditional hot Local Bubble model, a heating mechanism for the dimmer and cooler gas remains to be found. The second possible explanation is that of the ‘Hot Top’ model, in which the Local Cavity is an old supernova remnant in which no (or very little) million degree local plasma is presently required. Instead, the cavity is now thought to be filled with partially ionized cloudlets of temperature ∼7000 K that are surrounded by lower density envelopes of photo-ionized gas of temperature ∼20,000 K. Although this new scenario provides a natural explanation for many of the observations that were in conflict with the Local Hot Bubble model, we cannot (as yet) provide a satisfactory explanation or the emission levels observed in the B and Be ultra-soft X-ray bands.     

An XMM–Newton observation of the galaxy group MKW 4

We present an X-ray study of the galaxy group or poor cluster MKW 4. Working with XMM–Newton data we examine the distribution and properties of the hot gas which makes up the group halo. The inner halo shows some signs of structure, with circular or elliptical beta models providing a poor fit to the surface brightness profile. This may be evidence of large-scale motion in the inner halo, but we do not find evidence of sharp fronts or edges in the emission. The temperature of the halo declines in the core, with deprojected spectral fits showing a central temperature of ∼1.3 keV compared to ∼3 keV at 100 kpc. In the central ∼30 kpc of the group, multitemperature spectral models are required to fit the data, but they indicate a lack of gas at low temperatures. Steady-state cooling flow models provide poor fits to the inner regions of the group and the estimated cooling time of the gas is long except within the central dominant galaxy, NGC 4073. Abundance profiles show a sharp increase in the core of the group, with mean abundance rising by a factor of 2 in the centre of NGC 4073. Fitting individual elements shows the same trend, with high values of Fe, Si and S in the core. We estimate that ∼50 per cent of the Fe in the central 40 kpc was injected by Type Ia supernovae, in agreement with previous ASCA studies. Using our best-fitting surface brightness and temperature models, we calculate the mass, gas fraction, entropy and mass-to-light ratio of the group. At 100 kpc (∼0.1 virial radius) the total mass and gas entropy of the system (  ∼2 × 10¹³ M_⊙  and ∼300 keV cm²) are quite comparable to those of other systems of similar temperature, but the gas fraction is rather low (∼1 per cent). We conclude that MKW 4 is a fairly relaxed group, which has developed a strong central temperature gradient but not a large-scale cooling flow.     

Chandra X-ray observations of the 3C 295 cluster core

We examine the properties of the X-ray gas in the central regions of the distant ( z =0.46) , X-ray luminous cluster of galaxies surrounding the powerful radio source 3C 295, using observations made with the Chandra Observatory . Between radii of 50 and 500 kpc, the cluster gas is approximately isothermal with an emission-weighted temperature, kT ∼5 keV . Within the central 50-kpc radius this value drops to kT ∼3.7 keV . The spectral and imaging Chandra data indicate the presence of a cooling flow within the central 50-kpc radius of the cluster, with a mass deposition rate of approximately 280 M_⊙ yr⁻¹. We estimate an age for the cooling flow of 1–2 Gyr , which is approximately 1000 times older than the central radio source. We find no evidence in the X-ray spectra or images for significant heating of the X-ray gas by the radio source. We report the detection of an edge-like absorption feature in the spectrum for the central 50-kpc region, which may be caused by oxygen-enriched dust grains. The implied mass in metals seen in absorption could have been accumulated by the cooling flow over its lifetime. Combining the results on the X-ray gas density profile with radio measurements of the Faraday rotation measure in 3C 295, we estimate the magnetic field strength in the region of the cluster core to be B ∼12 μG .     

Molecular line profiles as diagnostics of protostellar collapse: modelling the 'blue asymmetry' in inside-out infall

The evolution of star-forming core analogues undergoing inside-out collapse is studied with a multipoint chemodynamical model which self-consistently computes the abundance distribution of chemical species in the core. For several collapse periods the output chemistry of infalling tracer species such as HCO⁺, CS and N₂H⁺ is then coupled to an accelerated Λ-iteration radiative transfer code, which predicts the emerging molecular line profiles using two different input gas/dust temperature distributions. We investigate the sensitivity of the predicted spectral line profiles and line asymmetry ratios to the core temperature distribution, the time-dependent model chemistry, as well as to ad hoc abundance distributions. The line asymmetry is found to be strongly dependent on the adopted chemical abundance distribution. In general, models with a warm central region show higher values of blue asymmetry in optically thick HCO⁺ and CS lines than models with a starless core temperature profile. We find that in the formal context of Shu-type inside-out infall, and in the absence of rotation or outflows, the relative blue asymmetry of certain HCO⁺ and CS transitions is a function of time and, subject to the foregoing caveats, can act as a collapse chronometer. The sensitivity of simulated HCO⁺ line profiles to linear radial variations, subsonic or supersonic, of the internal turbulence field is investigated in the separate case of static cores.     

Cooling flows in clusters of galaxies

Summary X-ray images and spectra of clusters of galaxies show strong evidence for cooling flows. In many clusters, the hot gas in the core is cooling at rates of 100Myr^–1 and greater. Few traces of the cooled gas have been observed, but it probably forms into low-mass stars (perhaps brown dwarf or even Jupiter-mass objects). X-ray surface-brightness profiles show that the cooling gas is highly inhomogeneous. Overdense gas cools rapidly to form cooled clumps distributed throughout the flow, with little of the gas ever reaching the cluster centre. Cooled and cooling clumps are disrupted because of their motion relative to the remainder of the gas, tending to produce small cooled fragments and, ultimately, low-mass stars. Large molecular clouds, which are the sites of massive star formation in our galaxy, do not occur in the outer parts of cooling flows. There is evidence of larger gas clumps and the formation of more massive stars in the central few kpc of some cooling flows. It is argued that cooling flows efficiently form dark matter. This has wider implications for the formation of dark matter in massive galaxies.     

On the influence of cooling and heating processes on the Parker instability – II. Numerical simulations

We investigate the Parker instability (PI) influenced by thermal processes in a non-adiabatic, gravitationally stratified interstellar medium and discuss a model including the photoionization heating together with the supplemental heating mechanisms postulated by Reynolds, Haffner and Tufte. A cooling rate due to radiative losses is described by an approximation to the realistic cooling function of Dalgarno and McCray for ionized interstellar gas. An unperturbed initial state of the system simultaneously represents both a magnetohydrostatic and thermal equilibrium, and is thermally stable. We perform a set of 3D numerical magnetohydrodynamic simulations using the zeusmp code. We find that PI developing in the presence of non-adiabatic effects promotes a transition of gas in magnetic valleys to a thermally unstable regime. We find that the region of initially enhanced density due to PI starts to condense more as the result of thermal instability action. The density in this region rises above the classical isothermal limit of two times the equilibrium value at the mid-plane. The maximum density in an evolved system reaches 10–40 times the equilibrium value at the mid-plane, and the structures so formed attain oval shapes. These results lead to the conclusion that PI, operating in the presence of realistic cooling and heating processes, can trigger the formation of dense clouds, which may give rise to giant molecular complexes.     

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 initial conditions of isolated star formation – V. ISOPHOT imaging and the temperature and energy balance of pre-stellar cores

ISO data taken with the long-wavelength imaging photo-polarimeter ISOPHOT are presented of 18 pre-stellar cores at three far-infrared wavelengths, 90, 170 and 200 μm. Most of the cores are detected clearly at 170 and 200 μm, but only one is detected strongly at 90 μm, indicating that mostly they are very cold, with typical temperatures of only ∼     . Colour temperature images are constructed for each of the cores. Most of the cores are seen either to be isothermal, or to have associated temperature gradients from the core centres to their edges, with all except one being cooler at the centre. We compare the data with previous ISOCAM absorption data, and calculate the energy balance for those cores in common between the two samples. We find that the energy radiated by each core in the far-infrared is similar to that absorbed at shorter wavelengths. Hence there is no evidence for a central heating source in any of the cores – even those for which previous evidence for core contraction exists. This is all consistent with external heating of the cores by the local interstellar radiation field, confirming their pre-stellar nature.     

Towards a holistic view of the heating and cooling of the intracluster medium

X-ray clusters are conventionally divided into two classes: 'cool core' (CC) clusters and 'non-cool core' (NCC) clusters. Yet relatively little attention has been given to the origins of this apparent dichotomy and, in particular, to the energetics and thermal histories of the two classes. We develop a model for the entropy profiles of clusters starting from the configuration established by gravitational shock heating and radiative cooling. At large radii, gravitational heating accounts for the observed profiles and their scalings well. However, at small and intermediate radii, radiative cooling and gravitational heating cannot be combined to explain the observed profiles of either CC or NCC clusters. The inferred entropy profiles of NCC clusters require that material is 'pre-heated' prior to cluster collapse in order to explain the absence of low-entropy (cool) material in these systems. We show that a similar modification is also required in CC clusters in order to match their entropy profiles at intermediate radii. In CC clusters, this modification is unstable, and an additional process is required to prevent cooling below a temperature of a few keV. We show that this can be achieved by adding a self-consistent active galactic nuclei (AGN) feedback loop in which the lowest entropy, most rapidly cooling material is heated and rises buoyantly to mix with material at larger radii. The resulting model does not require fine-tuning and is in excellent agreement with a wide variety of observational data from Chandra and XMM–Newton , including entropy and gas density profiles, the luminosity–temperature relation and high-resolution spectra. The spread in cluster core morphologies is seen to arise because of the steep dependence of the central cooling time on the initial level of pre-heating. Some of the other implications of this model are briefly discussed.     

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.     

On X-ray variability in narrow-line and broad-line active galactic nuclei

We propose a novel mathematical method to construct an exact polytropic sphere in self-gravitating hydrostatic equilibrium, improving the non-linear Poisson equation. The central boundary condition for the present equation requires a ratio of gas pressure to total one at the centre, which is uniquely identified by the whole mass and molecular weight of the system. The special solution derived from the Lane–Emden equation can be reproduced. This scheme is now available for modelling the molecular cloud cores in interstellar media. The mass–radius relation of the first core is found to be consistent with the recent results of radiation hydrodynamic simulations.     

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.     

An elementary model for the dust cycle in galaxies

We give an elementary model for the evolution of dust in galaxies, based on abundance arguments. The model takes account of grain core production in both supernovae and giant stars, and includes mantle growth in the interstellar medium. Destruction of grain cores does not appear to be a dominant effect. We show that a self-consistent picture can be made in which the interstellar dust mass is an approximately constant fraction of the heavy element mass in the interstellar medium. This result is demonstrated to be essentially independent of outflow or inflow of interstellar material.     

Cooling flows induced by Compton cooling resulting from luminous quasars in clusters of galaxies

We have studied the effects of Compton cooling on cooling flow by performing numerical hydrodynamic calculations of the time evolution of hot gas in clusters of galaxies with luminous quasars. We assumed various temperatures for the hot intracluster gas. We have shown that the Compton cooling resulting from very luminous quasars is effective in inducing cooling flow, before radiative cooling flow is realized. The mass flux resulting from the cooling flow increases with time, as the Compton-cooled region expands. However, the mass flux of the Compton cooling flow is not large, less than 1 M⊙ yr⁻¹ in our model, because the Compton-cooled region is limited to an inner galactic region around a quasar. Even after the quasar active phase has ceased, the cooling flow will continue for at least 10⁹ yr. The accreted mass is enough to explain X-ray absorption lines in high-redshift quasars, if the Compton-cooled gas is compressed by high-pressure intracluster gas.     

Radiative feedback from quasars and the growth of massive black holes in stellar spheroids

We discuss the importance of feedback via photoionization and Compton heating on the co-evolution of massive black holes (MBHs) at the centre of spheroidal galaxies, and their stellar and gaseous components. We first assess the energetics of the radiative feedback from a typical quasar on the ambient interstellar medium (ISM). We then demonstrate that the observed   M _BH–σ  relation could be established following the conversion of most of the gas of an elliptical progenitor into stars, specifically when the gas-to-stars mass ratio in the central regions has dropped to a low level  ∼0.01  or less, so that gas cooling is no longer able to keep up with the radiative heating by the growing central massive black hole (MBH). A considerable amount of the remaining gas will be expelled and both MBH accretion and star formation will proceed at significantly reduced rates thereafter, in agreement with observations of present-day ellipticals. We find further support for this scenario by evolving over an equivalent Hubble time a simple, physically based toy model that additionally takes into account the mass and energy return for the spheroid evolving stellar population, a physical ingredient often neglected in similar approaches.     

Fuelling quasars with hot gas

We consider a model for quasar formation in which massive black holes are formed and fuelled largely by the accretion of hot gas during the process of galaxy formation. In standard hierarchical collapse models, objects about the size of normal galaxies and larger form a dense hot atmosphere when they collapse. We show that if such an atmosphere forms a nearly 'maximal' cooling flow, then a central black hole can accrete at close to its Eddington limit. This leads to exponential growth of a seed black hole, resulting in a quasar in some cases. In this model, the first quasars form soon after the first collapses to produce hot gas. The hot gas is depleted as time progresses, mostly by cooling, so that the accretion rate eventually falls below the threshold for advection-dominated accretion, at which stage radiative efficiency plummets and any quasar turns off. A simple implementation of this model, incorporated into a semi-analytical model for galaxy formation, overproduces quasars when compared with observed luminosity functions, but is consistent with models of the X-ray background, which indicate that most accretion is obscured. It produces few quasars at high redshift owing to the lack of time needed to grow massive black holes. Quasar fuelling by hot gas provides a minimum level, sufficient to power most quasars at redshifts between one and two, to which other sources of fuel can be added. The results are sensitive to feedback effects, such as might result from radio jets and other outflows.