SMAUG: a new technique for the deprojection of galaxy clusters

This paper presents a new technique for reconstructing the spatial distributions of hydrogen, temperature and metal abundance of a galaxy cluster. These quantities are worked out from the X-ray spectrum, modeled starting from few analytical functions describing their spatial distributions. These functions depend upon some parameters, determined by fitting the model to the observed spectrum. We have implemented this technique as a new model in the XSPEC software analysis package. We describe the details of the method, and apply it to work out the structure of the cluster A1795. We combine the observation of three satellites, exploiting the high spatial resolution of Chandra for the cluster core, the wide collecting area of XMM-Newton for the intermediate regions and the large field of view of Beppo-SAX for the outer regions. We also test the validity and precision of our method by i) comparing its results with those from a geometrical deprojection, ii) examining the spectral residuals at different radii of the cluster and iii) reprojecting the unfolded profiles and comparing them directly to the measured quantities. Our analytical method yields the parameters defining the spatial functions directly from the spectra. Their explicit knowledge allows a straightforward derivation of other indirect physical quantities like the gravitating mass, as well as a fast and easy estimate of the profiles uncertainties.Comment: 24 pages, 11 figures, 3 tables; emulateapj; accepted for publication in the Astrophysical Journa

The luminosity function of cluster galaxies. II. Data reduction procedures applied to the cluster Abell 496

We initiated a large project aimed to estimate the Luminosity Function of galaxies in clusters and to evaluate its relation to cluster morphology. With this paper we deem necessary to outline the general procedures of the data reduction and details of the data analysis. The cluster sample includes the brightest southern ROSAT all-sky survey clusters with z < 0.1. These have been observed in three colours g, r, i, and mapped up to a few core radii using a mosaic of CCD frames. E/S0 galaxies in the cluster core are singled out both by morphology (for the brightest galaxies), and by colour. The details of the data reduction procedure are illustrated via the analysis of the cluster Abell 496, which has been used as a pilot cluster for the whole program. The related photometric catalogue consists of 2355 objects. The limiting magnitudes (the reference Surface Brightness is given in parenthesis) in the various colours are respectively g(25.5) = 24.14, r(25.5) = 24.46, i(25.0) = 23.75$. These correspond to the limiting absolute magnitudes -12.28, -11.96 and -12.67 (H_0=50 km/sec/Mpc).Comment: 17 pages, 19 ps figures, aa.cl

Slow quench dynamics of the Kitaev model: anisotropic critical point and effect of disorder

We study the non-equilibrium slow dynamics for the Kitaev model both in the presence and the absence of disorder. For the case without disorder, we demonstrate, via an exact solution, that the model provides an example of a system with an anisotropic critical point and exhibits unusual scaling of defect density $n$ and residual energy $Q$ for a slow linear quench. We provide a general expression for the scaling of $n$ ($Q$) generated during a slow power-law dynamics, characterized by a rate $\tau^{-1}$ and exponent $\alpha$, from a gapped phase to an anisotropic quantum critical point in $d$ dimensions, for which the energy gap $\Delta_{\vec k} \sim k_i^z$ for $m$ momentum components ($i=1..m$) and $\sim k_i^{z'}$ for the rest $d-m$ components ($i=m+1..d$) with $z\le z'$: $n \sim \tau^{-[m + (d-m)z/z']\nu \alpha/(z\nu \alpha +1)}$ ($Q \sim \tau^{-[(m+z)+ (d-m)z/z']\nu \alpha/(z\nu \alpha +1)}$). These general expressions reproduce both the corresponding results for the Kitaev model as a special case for $d=z'=2$ and $m=z=\nu=1$ and the well-known scaling laws of $n$ and $Q$ for isotropic critical points for $z=z'$. We also present an exact computation of all non-zero, independent, multispin correlation functions of the Kitaev model for such a quench and discuss their spatial dependence. For the disordered Kitaev model, where the disorder is introduced via random choice of the link variables $D_n$ in the model's Fermionic representation, we find that $n \sim \tau^{-1/2}$ and $Q\sim \tau^{-1}$ ($Q\sim \tau^{-1/2}$) for a slow linear quench ending in the gapless (gapped) phase. We provide a qualitative explanation of such scaling.Comment: 10 pages, 11 Figs. v

Cluster mergers, core oscillations, and cold fronts

We use numerical simulations with hydrodynamics to demonstrate that a class of cold fronts in galaxy clusters can be attributed to oscillations of the dark matter distribution. The oscillations are initiated by the off-axis passage of a low-mass substructure. From the simulations, we derive three observable morphological features indicative of oscillations: 1) The existence of compressed isophotes; 2) The regions of compression must be alternate (opposite and staggered) and lie on an axis passing through the center of the cluster; 3) The gradient of each compression region must pass through the center of the cluster. Four of six clusters reported in the literature to have cold fronts have morphologies consistent with the presence of oscillations. The clusters with oscillations are A496, A1795, A2142, and RX J1720.1+2638. Galaxy clusters A2256 and A3667 are not consistent so the cold fronts are interpreted as group remnants. The oscillations may be able to provide sufficient energy to solve the cooling-flow problem and, importantly, provide it over an extended duration.Comment: Submitted to ApJ. 11 pages, 9 figure

Radiative cooling, heating and thermal conduction in M87

The crisis of the standard cooling flow model brought about by Chandra and XMM-Newton observations of galaxy clusters, has led to the development of several models which explore different heating processes in order to assess if they can quench the cooling flow. Among the most appealing mechanisms are thermal conduction and heating through buoyant gas deposited in the ICM by AGNs. We combine Virgo/M87 observations of three satellites (Chandra, XMM-Newton and Beppo-SAX) to inspect the dynamics of the ICM in the center of the cluster. Using the spectral deprojection technique, we derive the physical quantities describing the ICM and determine the extra-heating needed to balance the cooling flow assuming that thermal conduction operates at a fixed fraction of the Spitzer value. We assume that the extra-heating is due to buoyant gas and we fit the data using the model developed by Ruszkowski and Begelman (2002). We derive a scale radius for the model of $\sim 5$ kpc, which is comparable with the M87 AGN jet extension, and a required luminosity of the AGN of a $few \times 10^{42}$ erg s$^{-1}$, which is comparable to the observed AGN luminosity. We discuss a scenario where the buoyant bubbles are filled of relativistic particles and magnetic field responsible for the radio emission in M87. The AGN is supposed to be intermittent and to inject populations of buoyant bubbles through a succession of outbursts. We also study the X-ray cool component detected in the radio lobes and suggest that it is structured in blobs which are tied to the radio buoyant bubbles.Comment: 25 pages, 10 figures and 2 tables. Accepted for publication in Ap

The structure of Abell 1351: a bimodal galaxy cluster with peculiar diffuse radio emission

We aim to review the internal structure and dynamics of the Abell 1351 cluster, shown to host a radio halo with a quite irregular shape. Our analysis is based on radial velocity data for 135 galaxies obtained at the Telescopio Nazionale Galileo. We combine galaxy velocities and positions to select 95 cluster galaxy members and analyse the internal dynamics of the whole cluster. We also examine X-ray data retrieved from Chandra and XMM archives. We measure the cluster redshift, =0.325, the line-of-sight (LOS) velocity dispersion, \sigma_v~1500 km/s, and the X-ray temperature, kT~9 keV. From both X-ray and optical data independently, we estimate a large cluster mass, in the 1--4 $10^{15}$ M$_\odot$ range. We attribute the extremely high value of \sigma_v to the bimodality in the velocity distribution. We find evidence of a significant velocity gradient and optical 3D substructure. The X-ray analysis also shows many features in favour of a complex cluster structure, probably supporting an ongoing merger of substructures in Abell 1351. The observational scenario agrees with the presence of two main subclusters in the northern region, each with its brightest galaxy (BCG1 and BCG2), detected as the two most important X-ray substructures with a rest-frame LOS velocity difference of \Delta v~2500 km/s (in the rest frame) and probably being in large part aligned with the LOS. We conclude that Abell 1351 is a massive merging cluster. The details of the cluster structure allow us to interpret the quite asymmetric radio halo as a `normal' halo plus a southern relic, strongly supporting a previous suggestion based only on inspection of radio and preliminary X-ray data.Comment: 13 pages, 13 figures, 1 tabl

A Chandra archival study of the temperature and metal abundance profiles in hot Galaxy Clusters at 0.1 < z < 0.3

We present the analysis of the temperature and metallicity profiles of 12 galaxy clusters in the redshift range 0.1--0.3 selected from the Chandra archive with at least ~20,000 net ACIS counts and kT>6 keV. We divide the sample between 7 Cooling-Core (CC) and 5 Non-Cooling-Core (NCC) clusters according to their central cooling time. We find that single power-laws can describe properly both the temperature and metallicity profiles at radii larger than 0.1 r_180 in both CC and NCC systems, showing the NCC objects steeper profiles outwards. A significant deviation is only present in the inner 0.1 r_180. We perform a comparison of our sample with the De Grandi & Molendi BeppoSAX sample of local CC and NCC clusters, finding a complete agreement in the CC cluster profile and a marginally higher value (at ~1sigma) in the inner regions of the NCC clusters. The slope of the power-law describing kT(r) within 0.1 r_180 correlates strongly with the ratio between the cooling time and the age of the Universe at the cluster redshift, being the slope >0 and tau_c/tau_age<=0.6 in CC systems.Comment: 12 pages, 6 figures, Accepted for publication by the Astrophysical Journa

Where does the gas fueling star formation in BCGs originate?

We investigate the relationship between X-ray cooling and star formation in brightest cluster galaxies (BCGs). We present an X-ray spectral analysis of the inner regions, 10-40 kpc, of six nearby cool core clusters (z<0.35) observed with Chandra ACIS. This sample is selected on the basis of the high star formation rate (SFR) observed in the BCGs. We restrict our search for cooling gas to regions that are roughly cospatial with the starburst. We fit single- and multi-temperature mkcflow models to constrain the amount of isobarically cooling intracluster medium (ICM). We find that in all clusters, below a threshold temperature ranging between 0.9 and 3 keV, only upper limits can be obtained. In four out of six objects, the upper limits are significantly below the SFR and in two, namely A1835 and A1068, they are less than a tenth of the SFR. Our results suggests that a number of mechanisms conspire to hide the cooling signature in our spectra. In a few systems the lack of a cooling signature may be attributed to a relatively long delay time between the X-ray cooling and the star burst. However, for A1835 and A1068, where the X-ray cooling time is shorter than the timescale of the starburst, a possible explanation is that the region where gas cools out of the X-ray phase extends to very large radii, likely beyond the core of these systems.Comment: to appear in A&