Time-dependent Circulation Flows: Iron Enrichment in Cooling Flows with Heated Return Flows

We describe a new type of dynamical model for hot gas in galaxy groups and clusters in which gas moves simultaneously in both radial directions. Circulation flows are consistent with (1) the failure to observe cooling gas in X-ray spectra, (2) multiphase gas observed near the centers of these flows and (3) the accumulation of iron in the hot gas from Type Ia supernovae in the central galaxy. Dense inflowing gas cools, producing a positive central temperature gradient, as in normal cooling flows. Bubbles of hot, buoyant gas flow outward. Circulation flows eventually cool catastrophically if the outward flowing gas transports mass but no heat; to maintain the circulation both mass and energy must be supplied to the inflowing gas over a large volume, extending to the cooling radius. The rapid radial recirculation of gas produces a flat central core in the gas iron abundance, similar to many observations. We believe the circulation flows described here are the first gasdynamic, long-term evolutionary models that are in good agreement with all essential features observed in the hot gas: little or no gas cools as required by XMM spectra, the gas temperature increases outward near the center, and the gaseous iron abundance is about solar near the center and decreases outward.Comment: 17 pages (emulateapj5) with 6 figures; accepted by The Astrophysical Journa

Energetics of X-ray Cavities and Radio Lobes in Galaxy Clusters

We describe the formation and evolution of X-ray cavities in the hot gas of galaxy clusters. The cavities are formed only with relativistic cosmic rays that eventually diffuse into the surrounding gas. We explore the evolution of cavities formed with a wide range of cosmic ray diffusion rates. In previous numerical simulations cavities are formed by injecting ultra-hot but non-relativistic gas which increases the global thermal energy, offsetting radiative losses in the gas and helping to solve the cooling flow problem. Contrary to these results, we find that X-ray cavities formed solely by cosmic rays have a global cooling effect. As the cluster gas is displaced by cosmic rays, a global expansion of the cluster gas occurs with associated cooling that exceeds the heating by shock waves as the cavity forms. Most cosmic rays in our cavity evolutions do not move beyond the cooling radius even after 1 Gyr. The gas density is depressed by cosmic rays, becomes buoyant, and undergoes a significant outward mass transfer within the cooling radius, carrying cosmic rays and relatively low entropy gas to distant regions in the cluster where it remains for times exceeding the local cooling time in the hot gas. This post-cavity mass outflow due to cosmic ray buoyancy may contribute toward solving the cooling flow problem. We describe the energetics, size, stability and buoyant rise of X-ray cavities in detail, showing how each depends on the rate of cosmic ray diffusion.Comment: 17 pages, 8 figures, accepted by Ap

The X-Ray Concentration-Virial Mass Relation

We present the concentration (c)-virial mass (M) relation of 39 galaxy systems ranging in mass from individual early-type galaxies up to the most massive galaxy clusters, (0.06-20) x 10^{14} M_sun. We selected for analysis the most relaxed systems possessing the highest quality data currently available in the Chandra and XMM public data archives. A power-law model fitted to the X-ray c-M relation requires at high significance (6.6 sigma) that c decreases with increasing M, which is a general feature of CDM models. The median and scatter of the c-M relation produced by the flat, concordance LCDM model (Omega_m=0.3, sigma_8=0.9) agrees with the X-ray data provided the sample is comprised of the most relaxed, early forming systems, which is consistent with our selection criteria. Holding the rest of the cosmological parameters fixed to those in the concordance model the c-M relation requires 0.76< sigma_8 <1.07 (99% conf.), assuming a 10% upward bias in the concentrations for early forming systems. The tilted, low-sigma_8 model suggested by a new WMAP analysis is rejected at 99.99% confidence, but a model with the same tilt and normalization can be reconciled with the X-ray data by increasing the dark energy equation of state parameter to w ~ -0.8. When imposing the additional constraint of the tight relation between sigma_8 and Omega_m from studies of cluster abundances, the X-ray c-M relation excludes (>99% conf.) both open CDM models and flat CDM models with Omega_m ~1. This result provides novel evidence for a flat, low-Omega_m universe with dark energy using observations only in the local (z << 1) universe. Possible systematic errors in the X-ray mass measurements of a magnitude ~10% suggested by CDM simulations do not change our conclusions.Comment: Accepted for Publication in ApJ; 13 pages, 4 figures; minor clarifications and updates; correlation coefficients corrected in Table 1 (correct values were used in the analysis in previous versions); conclusions unchange

Heating Cooling Flows with Weak Shock Waves

The discovery of extended, approximately spherical weak shock waves in the hot intercluster gas in Perseus and Virgo has precipitated the notion that these waves may be the primary heating process that explains why so little gas cools to low temperatures. This type of heating has received additional support from recent gasdynamical models. We show here that outward propagating, dissipating waves deposit most of their energy near the center of the cluster atmosphere. Consequently, if the gas is heated by (intermittent) weak shocks for several Gyrs, the gas within 30-50 kpc is heated to temperatures that far exceed observed values. This heating can be avoided if dissipating shocks are sufficiently infrequent or weak so as not to be the primary source of global heating. Local PV and viscous heating associated with newly formed X-ray cavities are likely to be small, which is consistent with the low gas temperatures generally observed near the centers of groups and clusters where the cavities are located.Comment: 14 pages; Accepted by Ap

Self-Generated Magnetic Fields in Galactic Cooling Flows

Interstellar magnetic fields in elliptical galaxies are assumed to have their origin in stellar fields that accompany normal mass loss from an evolving population of old stars. The seed fields are amplified by interstellar turbulence driven by stellar mass loss and supernova events. These disordered fields are further amplified by time-dependent compression in the inward moving galactic cooling flow and are expected to dominate near the galactic core. Under favorable circumstances, fields similar in strength to those observed $B \sim 1-10~(r/10~kpc)^{-1.2}\mu$G can be generated solely from these natural galactic processes. In general the interstellar field throughout elliptical galaxies is determined by the outermost regions in the interstellar gas where the turbulent dynamo process can occur. Because of the long hydrodynamic flow times in galactic cooling flows, currently observed magnetic fields may result from periods of intense turbulent field amplification that occurred in the outer galaxy in the distant past. Particularly strong fields in ellipticals may result from ancient galactic mergers or shear turbulence introduced at the boundary between the interstellar gas and ambient cluster gas.Comment: 21 pages in AASTEX LaTeX with 2 figures; accepted by Astrophysical Journa

The Absence of Adiabatic Contraction of the Radial Dark Matter Profile in the Galaxy Cluster A2589

We present an X-ray analysis of the radial mass profile of the radio-quiet galaxy cluster A2589 between 0.015-0.25 r_vir using an XMM-Newton observation. Except for a ~16 kpc shift of the X-ray center of the R=45-60 kpc annulus, A2589 possesses a remarkably symmetrical X-ray image and is therefore an exceptional candidate for precision studies of its mass profile by applying hydrostatic equilibrium. The total gravitating matter profile is well described by the NFW model (fractional residuals <~10%) with c_vir=6.1 +/- 0.3 and M_vir = 3.3 +/- 0.3 x 10^{14} M_sun (r_vir = 1.74 +/- 0.05 Mpc) in excellent agreement with LCDM. When the mass of the hot ICM is subtracted from the gravitating matter profile, the NFW model fitted to the resulting dark matter (DM) profile produces essentially the same result. However, if a component accounting for the stellar mass (M_*) of the cD galaxy is included, then the NFW fit to the DM profile is substantially degraded in the central r ~50 kpc for reasonable M_*/L_V. Modifying the NFW DM halo by adiabatic contraction arising from the early condensation of stellar baryons in the cD galaxy further degrades the fit. The fit is improved substantially with a Sersic-like model recently suggested by high resolution N-body simulations but with an inverse Sersic index, alpha ~0.5, a factor of ~3 higher than predicted. We argue that neither random turbulent motions nor magnetic fields can provide sufficient non-thermal pressure support to reconcile the XMM mass profile with adiabatic contraction of a CDM halo assuming reasonable M_*/L_V. Our results support the scenario where, at least for galaxy clusters, processes during halo formation counteract adiabatic contraction so that the total gravitating mass in the core approximately follows the NFW profile.Comment: 15 pages, 11 figures, accepted for publication in ApJ. Minor changes to match published versio

Probing the Dark Matter and Gas Fraction in Relaxed Galaxy Groups with X-ray observations from Chandra and XMM

We present radial mass profiles within 0.3 r_vir for 16 relaxed galaxy groups-poor clusters (kT range 1-3 keV) selected for optimal mass constraints from the Chandra and XMM data archives. After accounting for the mass of hot gas, the resulting mass profiles are described well by a two-component model consisting of dark matter (DM), represented by an NFW model, and stars from the central galaxy. The stellar component is required only for 8 systems, for which reasonable stellar mass-to-light ratios (M/L_K) are obtained, assuming a Kroupa IMF. Modifying the NFW dark matter halo by adiabatic contraction does not improve the fit and yields systematically lower M/L_K. In contrast to previous results for massive clusters, we find that the NFW concentration parameter (c_vir) for groups decreases with increasing M_vir and is inconsistent with no variation at the 3 sigma level. The normalization and slope of the c_vir-M_vir relation are consistent with the standard LambdaCDM cosmological model with sigma_8 = 0.9. The small intrinsic scatter measured about the c_vir-M_vir relation implies the groups represent preferentially relaxed, early forming systems. The mean gas fraction (f =0.05 +/- 0.01) of the groups measured within an overdensity Delta=2500 is lower than for hot, massive clusters, but the fractional scatter (sigma_f/f=0.2) for groups is larger, implying a greater impact of feedback processes on groups, as expected.Comment: Accepted for Publication in ApJ; 30 pages, 9 figures. No changes from previous versio

Serendipitous XMM-Newton discovery of a cluster of galaxies at z=0.28

We report the discovery of a galaxy cluster serendipitously detected as an extended X-ray source in an offset observation of the group NGC 5044. The cluster redshift, z=0.281, determined from the optical spectrum of the brightest cluster galaxy, agrees with that inferred from the X-ray spectrum using the Fe K alpha complex of the hot ICM (z=0.27 +/- 0.01). Based on the 50 ks XMM observation, we find that within a radius of 383 kpc the cluster has an unabsorbed X-ray flux, f_X (0.5-2 keV) = 3.34 (+0.08, -0.13) x 10^{-13} erg/cm^2/s, a bolometric X-ray luminosity, L_X = 2.21 (+0.34, -0.19) x 10^{44} erg/s, kT = 3.57 +/- 0.12 keV, and metallicity, 0.60 +/- 0.09 solar. The cluster obeys the scaling relations for L_X and T observed at intermediate redshift. The mass derived from an isothermal NFW model fit is, M_vir = 3.89 +/- 0.35 x 10^{14} solar masses, with a concentration parameter, c = 6.7 +/- 0.4, consistent with the range of values expected in the concordance cosmological model for relaxed clusters. The optical properties suggest this could be a ``fossil cluster''.Comment: 5 pages, 4 colour figures, accepted for publication in Ap

Where Do Cooling Flows Cool?

Although only about 5 percent of the total baryonic mass in luminous elliptical galaxies is in the form of cooled interstellar gas, it is concentrated within the optical effective radius r_e where it influences the local dynamical mass. The mass of cooled gas must be spatially distributed since it greatly exceeds the masses of central black holes. We explore here the proposition that a population of low mass, optically dark stars is created from the cooled gas. We consider a wide variety of radial distributions for the interstellar cooling, but only a few are consistent with observed X-ray surface brightness profiles. In a region of concentrated interstellar cooling, the X-ray emission can exceed that observed, suggesting the presence of additional support by magnetic stresses or non-thermal pressure. In general we find that the mass of cooled gas contributes significantly to stellar dynamical mass to light ratios which vary with galactic radius. If the stars formed from cooled interstellar gas are optically luminous, their influence on the the mass to light ratio would be reduced. The mass of cooled gas inside r_e is sensitive to the rate that old stars lose mass, which is nearly independent of the initial mass function of the old stellar population.Comment: 18 pages with 6 figures; accepted by Astrophysical Journa