No Arabic abstract
We compute 3D gasdynamical models of jet outflows from the central AGN, that carry mass as well as energy to the hot gas in galaxy clusters and groups. These flows have many attractive attributes for solving the cooling flow problem: why the hot gas temperature and density profiles resemble cooling flows but show no spectral evidence of cooling to low temperatures. Subrelativistic jets, described by a few parameters, are assumed to be activated when gas flows toward or cools near a central SMBH. Using approximate models for a rich cluster (A1795), a poor cluster (2A 0336+096) and a group (NGC 5044), we show that mass-carrying jets with intermediate mechanical efficiencies ($sim10^{-3}$) can reduce for many Gyr the global cooling rate to or below the low values implied by X-spectra, while maintaining $T$ and $rho$ profiles similar to those observed, at least in clusters. Groups are much more sensitive to AGN heating and present extreme time variability in both profiles. Finally, the intermittency of the feedback generates multiple generations of X-ray cavities similar to those observed in Perseus cluster and elsewhere. Thus we also study the formation of buoyant bubbles and weak shocks in the ICM, along with the injection of metals by SNIa and stellar winds.
It is now widely accepted that heating processes play a fundamental role in galaxy clusters, struggling in an intricate but fascinating `dance with its antagonist, radiative cooling. Last generation observations, especially X-ray, are giving us tiny hints about the notes of this endless ballet. Cavities, shocks, turbulence and wide absorption-lines indicate the central active nucleus is injecting huge amount of energy in the intracluster medium. However, which is the real dominant engine of self-regulated heating? One of the model we propose are massive subrelativistic outflows, probably generated by a wind disc or just the result of the entrainment on kpc scale by the fast radio jet. Using a modified version of AMR code FLASH 3.2, we explored several feedback mechanisms which self-regulate the mechanical power. Two are the best schemes that answer our primary question, id est quenching cooling flow and at the same time preserving a cool core appearance for a long term evolution (7 Gyr): one more explosive (with efficiencies 0.005 - 0.01), triggered by central cooled gas, and the other gentler, ignited by hot gas Bondi accretion (with efficiency 0.1). These three-dimensional simulations show that the total energy injected is not the key aspect, but the results strongly depend on how energy is given to the ICM. We follow the dynamics of best model (temperature, density, SB maps and profiles) and produce many observable predictions: buoyant bubbles, ripples, turbulence, iron abundance maps and hydrostatic equilibrium deviation. We present a deep discussion of merits and flaws of all our models, with a critical eye towards observational concordance.
We assess the importance of AGN outflows with respect to the metal enrichment of the intracluster medium (ICM) in galaxy clusters. We use combined N-body and hydrodynamic simulations, along with a semi-numerical galaxy formation and evolution model. Using assumptions based on observations, we attribute outflows of metal-rich gas initiated by AGN activity to a certain fraction of our model galaxies. The gas is added to the model ICM, where the evolution of the metallicity distribution is calculated by the hydrodynamic simulations. For the parameters describing the AGN content of clusters and their outflow properties, we use the observationally most favorable values. We find that AGNs have the potential to contribute significantly to the metal content of the ICM or even explain the complete abundance, which is typically ~0.5 Z_sun in core regions. Furthermore, the metals end up being inhomogeneously distributed, in accordance with observations.
AGN heating, through massive subrelativistic outflows, might be the key to solve the long-lasting `cooling flow problem in cosmological systems. In a previous paper, we showed that cold accretion feedback and, to a lesser degree, Bondi self-regulated models are in fact able to quench cooling rates for several Gyr, at the same time preserving the mainc ool core features, like observed density and temperature profiles. Is it true also for lighter systems, such as galaxy groups? The answer is globally yes, although with remarkable differences. Adopting a modified version of the AMR code FLASH 3.2, we found that successful 3D simulations with cold and Bondi models are almost convergent in the galaxy group environment, with mechanical efficiencies in the range 5.e-4 - 1.e-3 and 5.e-2 - 1.e-1, respectively. The evolutionary storyline of galaxy groups is dominated by a quasi-continuous gentle injection with sub-Eddington outflows (with mechanical power and velocity around 1.e44 erg/s and 1.e4 km/s). The cold and hybrid accretion models present, in addition, very short quiescence periods, followed by moderate outbursts (10 times the previous phase), which generate a series of 10-20 kpc size cavities with high density contrast, temperatures similar to the ambient medium and cold rims. After shock heating, a phase of turbulence promotes gas mixing and diffusion of metals, which peak along jet-axis (up to 40 kpc) during active phases. At this stage the tunnel, produced by the enduring outflow (hard to detect in the mock SBx maps), is easily fragmented, producing tiny buoyant bubbles, typically a few kpc in size. In contrast to galaxy clusters, the AGN self-regulated feedback has to be persistent, with a `delicate touch, rather than rare and explosive strokes. This evolutionary difference dictates in the end that galaxy groups are not scaled-do
Previous numerical studies have shown that in protostellar outflows, the mass-velocity distribution $m(v)$ can be well described by a broken power law $propto v^{- gamma}$. On the other hand, recent observations of a sample of outflows show that the CO intensity-velocity distribution, closely related to $m(v)$, follows an exponential law $propto exp(-v/v_0)$. In the present work, we revisit the physical origin of the mass-velocity relationship $m(v)$ in jet-driven protostellar outflows. We investigate the respective contributions of the different regions of the outflow, from the swept-up ambient gas to the jet. We performed 3D numerical simulations of a protostellar jet propagating into a molecular cloud using the hydrodynamical code Yguazu-a. The code takes into account atomic and ionic species and was modified to include the H$_2$ gas. We find that by excluding the jet contribution, $m(v)$ is satisfyingly fitted with a single exponential law, with $v_0$ well in the range of observational values. The jet contribution results in additional components in the mass-velocity relationship. This empirical mass-velocity relationship is found to be valid locally in the outflow. The exponent $v_0$ is almost constant in time and for a given level of mixing between the ambient medium and the jet material. In general, $v_0$ displays only a weak spatial dependence. A simple modeling of the L1157 outflow successfully reproduces the various components of the observed CO intensity-velocity relationship. Our simulations indicate that these components trace the outflow cavity of swept-up gas and the material entrained along the jet, respectively. The CO intensity-velocity exponential law is naturally explained by the jet-driven outflow model. The entrained material plays an important role in shaping the mass-velocity profile.
Using Chandra data for a sample of 26 galaxy groups, we constrained the central cooling times (CCTs) of the ICM and classified the groups as strong cool-core (SCC), weak cool-core (WCC) and non-cool-core (NCC) based on their CCTs. The total radio luminosity of the brightest cluster galaxy (BCG) was obtained using radio catalog data and literature, which was compared to the CCT to understand the link between gas cooling and radio output. We determined K-band luminosities of the BCG with 2MASS data, and used it to constrain the masses of the SMBH, which were then compared to the radio output. We also tested for correlations between the BCG luminosity and the overall X-ray luminosity and mass of the group. The observed cool-core/non-cool-core fractions for groups are comparable to those of clusters. However, notable differences are seen. For clusters, all SCCs have a central temperature drop, but for groups, this is not the case as some SCCs have centrally rising temperature profiles. While for the cluster sample, all SCC clusters have a central radio source as opposed to only 45% of the NCCs, for the group sample, all NCC groups have a central radio source as opposed to 77% of the SCC groups. For clusters, there are indications of an anticorrelation trend between radio luminosity and CCT which is absent for the groups. Indications of a trend of radio luminosity with black hole mass observed in SCC clusters is absent for groups. The strong correlation observed between the BCG luminosity and the cluster X-ray luminosity/cluster mass weakens significantly for groups. We conclude that there are important differences between clusters and groups within the ICM cooling/AGN feedback paradigm.