Source author record

E. Churazov

E. Churazov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

86works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

86 published item(s)

preprint2021arXiv

The Galactic center chimneys: The base of the multiphase outflow of the Milky Way

Outflows and feedback are key ingredients of galaxy evolution. Evidence for an outflow arising from the Galactic center (GC) has recently been discovered at different wavelength. We show that the X-ray, radio, and infrared emissions are deeply interconnected, affecting one another and forming coherent features on scales of hundreds of parsecs, therefore indicating a common physical link associated with the GC outflow. We debate the location of the northern chimney and suggest that it might be located on the front side of the GC because of a significant tilt of the chimneys toward us. We report the presence of strong shocks at the interface between the chimneys and the interstellar medium, which are traced by radio and warm dust emission. We observe entrained molecular gas outflowing within the chimneys, revealing the multiphase nature of the outflow. In particular, the molecular outflow produces a long, strong, and structured shock along the northwestern wall of the chimney. Because of the different dynamical times of the various components of the outflow, the chimneys appear to be shaped by directed large-scale winds launched at different epochs. The data support the idea that the chimneys are embedded in an (often dominant) vertical magnetic field, which likely diverges with increasing latitude. We observe that the thermal pressure associated with the hot plasma appears to be smaller than the ram pressure of the molecular outflow and the magnetic pressure. This leaves open the possibility that either the main driver of the outflow is more powerful than the observed hot plasma, or the chimneys represent a "relic" of past and more powerful activity. These multiwavelength observations corroborate the idea that the chimneys represent the channel connecting the quasi-continuous, but intermittent, activity at the GC with the base of the Fermi bubbles.

preprint2020arXiv

Propagation of weak shocks in cool-core galaxy clusters in two-temperature magnetohydrodynamics with anisotropic thermal conduction

We investigate how different magnetohydrodynamic models of propagation of a weak (Mach number ~1.2) shock in the core of a galaxy cluster affect its observational appearance, using the Perseus cluster as our fiducial model. In particular, we study how thermal conduction, both isotropic and anisotropic, and ion-electron temperature equilibration modify a weak shock. Strong thermal conduction is expected to produce an electron temperature precursor. Less prominent pressure and density precursors are formed as well. A longer equilibration time largely reduces the density precursor, but does not change the electron temperature precursor much. When thermal conduction becomes anisotropic, the intracluster magnetic field imprints its characteristic spatial scale on the distortions of the shock induced by heat fluxes.

preprint2020arXiv

The discovery of radio halos in the Frontier Fields clusters Abell S1063 and Abell 370

Massive merging galaxy clusters often host diffuse Mpc-scale radio synchrotron emission. This emission originates from relativistic electrons in the ionized intracluster medium (ICM). An important question is how these synchrotron emitting relativistic electrons are accelerated. Our aim is to search for diffuse emission in the Frontier Fields clusters Abell S1063 and Abell 370 and characterize its properties. While these clusters are very massive and well studied at some other wavelengths, no diffuse emission has been reported for these clusters so far. We obtained 325 MHz Giant Metrewave Radio Telescope (GMRT) and 1--4 GHz Jansky Very Large Array (VLA) observations of Abell S1063 and Abell 370. We complement these data with Chandra and XMM-Newton X-ray observations. In our sensitive images, we discover radio halos in both clusters. In Abell S1063, a giant radio halo is found with a size of $\sim 1.2$ Mpc. The integrated spectral index between 325 MHz and 1.5 GHz is $-0.94\pm0.08$ and it steepens to $-1.77 \pm 0.20$ between 1.5 and 3.0 GHz. This spectral steepening provides support for the turbulent re-acceleration model for radio halo formation. Abell 370 hosts a faint radio halo mostly centred on the southern part of this binary merging cluster, with a size of $\sim 500-700$ kpc. The spectral index between 325 MHz and 1.5 GHz is $-1.10\pm0.09$. Both radio halos follow the known scaling relation between the cluster mass proxy $Y_{500}$ and radio power, consistent with the idea that they are related to ongoing cluster merger events.

preprint2020arXiv

The eROSITA X-ray telescope on SRG

eROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the primary instrument on the Spectrum-Roentgen-Gamma (SRG) mission, which was successfully launched on July 13, 2019, from the Baikonour cosmodrome. After the commissioning of the instrument and a subsequent calibration and performance verification phase, eROSITA started a survey of the entire sky on December 13, 2019. By the end of 2023, eight complete scans of the celestial sphere will have been performed, each lasting six months. At the end of this program, the eROSITA all-sky survey in the soft X-ray band (0.2--2.3\,keV) will be about 25 times more sensitive than the ROSAT All-Sky Survey, while in the hard band (2.3--8\,keV) it will provide the first ever true imaging survey of the sky. The eROSITA design driving science is the detection of large samples of galaxy clusters up to redshifts $z>1$ in order to study the large-scale structure of the universe and test cosmological models including Dark Energy. In addition, eROSITA is expected to yield a sample of a few million AGNs, including obscured objects, revolutionizing our view of the evolution of supermassive black holes. The survey will also provide new insights into a wide range of astrophysical phenomena, including X-ray binaries, active stars, and diffuse emission within the Galaxy. Results from early observations, some of which are presented here, confirm that the performance of the instrument is able to fulfil its scientific promise. With this paper, we aim to give a concise description of the instrument, its performance as measured on ground, its operation in space, and also the first results from in-orbit measurements.

preprint2020arXiv

X-raying molecular clouds with a short flare: probing statistics of gas density and velocity fields

We take advantage of a set of molecular cloud simulations to demonstrate a possibility to uncover statistical properties of the gas density and velocity fields using reflected emission of a short (with duration much less than the cloud's light-crossing time) X-ray flare. Such situation is relevant for the Central Molecular Zone of our Galaxy where several clouds get illuminated by a $\sim110$ yr-old flare from the supermassive black hole Sgr A*. Due to shortness of the flare ($Δt\lesssim1.6$ yrs), only a thin slice ($Δz\lesssim0.5$ pc) of the molecular gas contributes to the X-ray reflection signal at any given moment, and its surface brightness effectively probes the local gas density. This allows reconstructing the density probability distribution function over a broad range of scales with virtually no influence of attenuation, chemo-dynamical biases and projection effects. Such measurement is key to understanding the structure and star-formation potential of the clouds evolving under extreme conditions in the CMZ. For cloud parameters similar to the currently brightest in X-ray reflection molecular complex Sgr A, the sensitivity level of the best available data is sufficient only for marginal distinction between solenoidal and compressive forcing of turbulence. Future-generation X-ray observatories with large effective area and high spectral resolution will dramatically improve on that by minimising systematic uncertainties due to contaminating signals. Furthermore, measurement of the iron fluorescent line centroid with sub-eV accuracy in combination with the data on molecular line emission will allow direct investigation of the gas velocity field.

preprint2016arXiv

Arithmetic with X-ray images of galaxy clusters: effective equation of state for small-scale perturbations in the ICM

We discuss a novel technique of manipulating X-ray images of galaxy clusters to reveal the nature of small-scale density/temperature perturbations in the intra cluster medium (ICM). As we show, this technique can be used to differentiate between sound waves and isobaric perturbations in Chandra images of the Perseus and M87/Virgo clusters. The comparison of the manipulated images with the radio data and with the results of detailed spectral analysis shows that this approach successfully classifies the types of perturbations and helps to reveal their nature. For the central regions (5-100 kpc) of the M87 and Perseus clusters this analysis suggests that observed images are dominated by isobaric perturbations, followed by perturbations caused by bubbles of relativistic plasma and weak shocks. Such a hierarchy is best explained in a "slow" AGN feedback scenario, when much of the mechanical energy output of a central black hole is captured by the bubble enthalpy that is gradually released during buoyant rise of the bubbles. The "image arithmetic" works best for prominent structure and for datasets with excellent statistics, visualizing the perturbations with a given effective equation of state. The same approach can be extended to faint perturbations via cross-spectrum analysis of surface brightness fluctuations in X-ray images in different energy bands.

preprint2016arXiv

Frontier Fields Clusters: Deep Chandra Observations of the Complex Merger MACS J1149.6+2223

The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster.

preprint2016arXiv

Gamma-Ray emission from SN2014J near maximum optical light

The optical light curve of Type Ia supernovae (SNIa) is powered by thermalized gamma-rays produced by the decay of 56Ni and 56Co, the main radioactive isotopes synthesized by the thermonuclear explosion of a C/O white dwarf. Gamma-rays escaping the ejecta can be used as a diagnostic tool for studying the characteristics of the explosion. In particular, it is expected that the analysis of the early gamma emission, near the maximum of the optical light curve, could provide information about the distribution of the radioactive elements in the debris. In this paper, the gamma data obtained from SN2014J in M82 by the instruments on board of INTEGRAL are analyzed taking special care of the impact that the detailed spectral response has on the measurements of the intensity of the lines. The 158 keV emission of 56Ni has been detected in SN2014J at ~5 sigma at low energy with both ISGRI and SPI around the maximum of the optical light curve. After correcting the spectral response of the detector, the fluxes in the lines suggest that, in addition to the bulk of radioactive elements buried in the central layers of the debris, there is a plume of 56Ni, with a significance of ~3 sigma, moving at high velocity and receding from the observer. The mass of the plume is in the range of ~ 0.03-0.08 solar masses. No SNIa explosion model had predicted the mass and geometrical distribution of 56Ni suggested here. According to its optical properties, SN2014J looks as a normal SNIa. So it is extremely important to discern if it is also representative in the gamma-ray band.

preprint2016arXiv

Planck intermediate results. XLIII. The spectral energy distribution of dust in clusters of galaxies

Although infrared (IR) overall dust emission from clusters of galaxies has been statistically detected using data from the Infrared Astronomical Satellite (IRAS), it has not been possible to sample the spectral energy distribution (SED) of this emission over its peak, and thus to break the degeneracy between dust temperature and mass. By complementing the IRAS spectral coverage with Planck satellite data from 100 to 857 GHz, we provide new constraints on the IR spectrum of thermal dust emission in clusters of galaxies. We achieve this by using a stacking approach for a sample of several hundred objects from the Planck cluster sample; this procedure averages out fluctuations from the IR sky, allowing us to reach a significant detection of the faint cluster contribution. We also use the large frequency range probed by Planck, together with component-separation techniques, to remove the contamination from both cosmic microwave background anisotropies and the thermal Sunyaev-Zeldovich effect (tSZ) signal, which dominate below 353 GHz. By excluding dominant spurious signals or systematic effects, averaged detections are reported at frequencies between 353 and 5000 GHz. We confirm the presence of dust in clusters of galaxies at low and intermediate redshifts, yielding an SED with a shape similar to that of the Milky Way. Planck's beam does not allow us to investigate the detailed spatial distribution of this emission (e.g., whether it comes from intergalactic dust or simply the dust content of the cluster galaxies), but the radial distribution of the emission appears to follow that of the stacked SZ signal, and thus the extent of the clusters. The recovered SED allows us to constrain the dust mass responsible for the signal, as well as its temperature. We additionally explore the evolution of the IR emission as a function of cluster mass and redshift.

preprint2016arXiv

Polarization of thermal bremsstrahlung emission due to electron pressure anisotropy

Astrophysical plasmas are typically magnetized, with the Larmor radii of the charged particles many orders of magnitude smaller than their collisional mean free paths. The fundamental properties of such plasmas, e.g., conduction and viscosity, may depend on the instabilities driven by the anisotropy of the particle distribution functions and operating at scales comparable to the Larmor scales. We discuss a possibility that the pressure anisotropy of thermal electrons could produce polarization of thermal bremsstrahlung emission. In particular, we consider coherent large-scale motions in galaxy clusters to estimate the level of anisotropy driven by stretching of the magnetic-field lines by plasma flow and by heat fluxes associated with thermal gradients. Our estimate of the degree of polarization is $\sim 0.1 \%$ at energies $\gtrsim kT$. While this value is too low for the forthcoming generation of X-ray polarimeters, it is potentially an important proxy for the processes taking place at extremely small scales, which are impossible to resolve spatially. The absence of the effect at the predicted level may set a lower limit on the electron collisionality in the ICM. At the same time, the small value of the effect implies that it does not preclude the use of clusters as (unpolarized) calibration sources for X-ray polarimeters at this level of accuracy.

preprint2016arXiv

The discovery of lensed radio and X-ray sources behind the Frontier Fields cluster MACS J0717.5+3745 with the JVLA and Chandra

We report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

preprint2016arXiv

The nature and energetics of AGN-driven perturbations in the hot gas in the Perseus Cluster

Cores of relaxed galaxy clusters are often disturbed by AGN. Their Chandra observations revealed a wealth of structures induced by shocks, subsonic gas motions, bubbles of relativistic plasma, etc. In this paper, we determine the nature and energy content of gas fluctuations in the Perseus core by probing statistical properties of emissivity fluctuations imprinted in the soft- and hard-band X-ray images. About 80 per cent of the total variance of perturbations on ~ 8-70 kpc scales in the inner region have an isobaric nature, i.e., are consistent with slow displacements of the gas in pressure equilibrium with ambient medium. Observed variance translates to the ratio of non-thermal to thermal energy of ~13 per cent. In the region dominated by weak "ripples", about half of the total variance is also associated with isobaric perturbations on scales ~ a few tens of kpc. If these isobaric perturbations are induced by buoyantly rising bubbles, then these results suggest that most of the AGN-injected energy should first go into bubbles rather than into shocks. Using simulations of a shock propagating through the Perseus atmosphere, we found that models reproducing the observed features of a central shock have more than 50 per cent of the AGN-injected energy associated with the bubble enthalpy and only about 20 per cent is carried away with the shock. Such energy partition is consistent with the AGN-feedback model, mediated by bubbles of relativistic plasma, and supports the importance of turbulence in the balance between gas heating and radiative cooling.

preprint2015arXiv

(No) dimming of X-ray clusters beyond z~1 at fixed mass: crude redhshifts and masses from raw X-ray and SZ data

Scaling relations in the LCDM Cosmology predict that for a given mass the clusters formed at larger redshift are hotter, denser and therefore more luminous in X-rays than their local z~0 counterparts. This effect overturns the decrease in the observable X-ray flux so that it does not decrease at z > 1, similar to the SZ signal. Provided that scaling relations remain valid at larger redshifts, X-ray surveys will not miss massive clusters at any redshift, no matter how far they are. At the same time, the difference in scaling with mass and distance of the observable SZ and X-ray signals from galaxy clusters at redshifts $z\lesssim 2$ offers a possibility to crudely estimate the redshift and the mass of a cluster. This might be especially useful for preselection of massive high-redshift clusters and planning of optical follow-up for overlapping surveys in X-ray (e.g., by SRG/eRosita) and SZ (e.g. Planck, SPT and ACT).

preprint2015arXiv

Comparison of simple mass estimators for slowly rotating elliptical galaxies

We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $σ_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{σ_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.

preprint2015arXiv

Frontier Fields Clusters: Chandra and JVLA View of the Pre-Merging Cluster MACS J0416.1-2403

Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster constituents, including shock fronts, cold fronts, X-ray substructure, radio halos, and offsets between the dark matter and the gas components. Using observations from Chandra, the Jansky Very Large Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a multiwavelength analysis of the merging Frontier Fields cluster MACS J0416.1-2403 (z=0.396), which consists of a NE and a SW subclusters whose cores are separated on the sky by ~250 kpc. We find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity that corresponds to a compression factor of ~1.5. The discontinuity was most likely caused by the interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of the subcluster's center. For both the NE and the SW subclusters, the dark matter and the gas components are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of (1.06 +/- 0.09) x 10^{24} W Hz^{-1}, which is somewhat lower than expected based on the X-ray luminosity of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered a rare ultra-steep spectrum halo.

preprint2015arXiv

Gamma-rays from Type Ia supernova SN2014J

The whole set of INTEGRAL observations of type Ia supernova SN2014J, covering the period 19-162 days after the explosion has being analyzed. For spectral fitting the data are split into "early" and "late" periods covering days 19-35 and 50-162, respectively, optimized for $^{56}$Ni and $^{56}$Co lines. As expected for the early period much of the gamma-ray signal is confined to energies below $\sim$200 keV, while for the late period it is most strong above 400 keV. In particular, in the late period $^{56}$Co lines at 847 and 1248 keV are detected at 4.7 and 4.3 $σ$ respectively. The lightcurves in several representative energy bands are calculated for the entire period. The resulting spectra and lightcurves are compared with a subset of models. We confirm our previous finding that the gamma-ray data are broadly consistent with the expectations for canonical 1D models, such as delayed detonation or deflagration models for a near-Chandrasekhar mass WD. Late optical spectra (day 136 after the explosion) show rather symmetric Co and Fe lines profiles, suggesting that unless the viewing angle is special, the distribution of radioactive elements is symmetric in the ejecta.

preprint2015arXiv

On the Nature of X-ray Surface Brightness Fluctuations in M87

X-ray images of galaxy clusters and gas-rich elliptical galaxies show a wealth of small-scale features which reflect fluctuations in density and/or temperature of the intra-cluster medium. In this paper we study these fluctuations in M87/Virgo, to establish whether sound waves/shocks, bubbles or uplifted cold gas dominate the structure. We exploit the strong dependence of the emissivity on density and temperature in different energy bands to distinguish between these processes. Using simulations we demonstrate that our analysis recovers the leading type of fluctuation even in the presence of projection effects and temperature gradients. We confirm the isobaric nature of cool filaments of gas entrained by buoyantly rising bubbles, extending to 7' to the east and south-west, and the adiabatic nature of the weak shocks at 40" and 3' from the center. For features of 5--10 kpc, we show that the central 4'x 4' region is dominated by cool structures in pressure equilibrium with the ambient hotter gas while up to 30 percent of the variance in this region can be ascribed to adiabatic fluctuations. The remaining part of the central 14'x14' region, excluding the arms and shocks described above, is dominated by apparently isothermal fluctuations (bubbles) with a possible admixture (at the level of about 30 percent) of adiabatic (sound waves) and by isobaric structures. Larger features, of about 30 kpc, show a stronger contribution from isobaric fluctuations. The results broadly agree with an AGN feedback model mediated by bubbles of relativistic plasma.

preprint2015arXiv

Planck 2013 results. XXIX. The Planck catalogue of Sunyaev-Zeldovich sources: Addendum

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

preprint2015arXiv

Planck 2015 results. XXII. A map of the thermal Sunyaev-Zeldovich effect

We have constructed all-sky y-maps of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 30 to 857 GHz frequency channel maps from the Planck satellite survey. These reconstructed y-maps are delivered as part of the Planck 2015 release. The y-maps are characterised in terms of noise properties and residual foreground contamination, mainly thermal dust emission at large angular scales and CIB and extragalactic point sources at small angular scales. Specific masks are defined to minimize foreground residuals and systematics. Using these masks we compute the y-map angular power spectrum and higher order statistics. From these we conclude that the y-map is dominated by tSZ signal in the multipole range, 20-600. We compare the measured tSZ power spectrum and higher order statistics to various physically motivated models and discuss the implications of our results in terms of cluster physics and cosmology.

preprint2015arXiv

Planck 2015 results. XXIII. The thermal Sunyaev-Zeldovich effect--cosmic infrared background correlation

We use Planck data to detect the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and the infrared emission from the galaxies that make up the the cosmic infrared background (CIB). We first perform a stacking analysis towards Planck-confirmed galaxy clusters. We detect infrared emission produced by dusty galaxies inside these clusters and demonstrate that the infrared emission is about 50% more extended than the tSZ effect. Modelling the emission with a Navarro--Frenk--White profile, we find that the radial profile concentration parameter is $c_{500} = 1.00^{+0.18}_{-0.15}$. This indicates that infrared galaxies in the outskirts of clusters have higher infrared flux than cluster-core galaxies. We also study the cross-correlation between tSZ and CIB anisotropies, following three alternative approaches based on power spectrum analyses: (i) using a catalogue of confirmed clusters detected in Planck data; (ii) using an all-sky tSZ map built from Planck frequency maps; and (iii) using cross-spectra between Planck frequency maps. With the three different methods, we detect the tSZ-CIB cross-power spectrum at significance levels of (i) 6 $σ$, (ii) 3 $σ$, and (iii) 4 $σ$. We model the tSZ-CIB cross-correlation signature and compare predictions with the measurements. The amplitude of the cross-correlation relative to the fiducial model is $A_{\rm tSZ-CIB}= 1.2\pm0.3$. This result is consistent with predictions for the tSZ-CIB cross-correlation assuming the best-fit cosmological model from Planck 2015 results along with the tSZ and CIB scaling relations.

preprint2015arXiv

Planck intermediate results. XL. The Sunyaev-Zeldovich signal from the Virgo cluster

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck's wide angular scale and frequency coverage, together with its high sensitivity, allow a detailed study of this large object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure, which correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find a good agreement between the SZ signal (or Compton paranmeter, y_c) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Due to its proximity to us, the gas beyond the virial radius can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. We study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusive warm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions.

preprint2015arXiv

Stripped elliptical galaxies as probes of ICM physics: I. Tails, wakes, and flow patterns in and around stripped ellipticals

Elliptical cluster galaxies are progressively stripped of their atmospheres due to their motion through the intra-cluster medium (ICM). Deep X-ray observations reveal the fine-structure of the galaxy's remnant atmosphere and its gas tail and wake. This fine-structure depends on dynamic conditions (galaxy potential, initial gas contents, orbit through the host cluster), orbital stage (early infall, pre-/post-pericenter passage), and ICM plasma properties (thermal conductivity, viscosity, magnetic field structure). We aim to disentangle dynamic and plasma effects in order to use stripped ellipticals as probes of ICM plasma properties. This first paper of a series investigates the hydrodynamics of progressive gas stripping by means of inviscid hydrodynamical simulations. We distinguish a long-lasting initial relaxation phase and a quasi-steady stripping phase. During quasi-steady stripping, the ICM flow around the remnant atmosphere resembles the flow around solid bodies, including a `deadwater' region in the near wake. Gas is stripped from the remnant atmosphere predominantly at its sides via Kelvin-Helmholtz instabilities. The downstream atmosphere is largely shielded from the ICM wind and thus shaped into a tail. Observationally, both, this `remnant tail' and the stripped gas in the wake can appear as a `tail', but only in the wake can galactic gas mix with the ambient ICM. While the qualitative results are generic, the simulations presented here are tailored to the Virgo elliptical galaxy M89 (NGC 4552) for the most direct comparison to observations. Papers II and III of this series describe the effect of viscosity and compare to Chandra and XMM-Newton observations, respectively.

preprint2015arXiv

Stripped elliptical galaxies as probes of ICM physics: II. Stirred, but mixed? Viscous and inviscid gas stripping of the Virgo elliptical M89

Elliptical galaxies moving through the intra-cluster medium (ICM) are progressively stripped of their gaseous atmospheres. X-ray observations reveal the structure of galactic tails, wakes, and the interface between the galactic gas and the ICM. This fine-structure depends on dynamic conditions (galaxy potential, initial gas contents, orbit in the host cluster), orbital stage (early infall, pre-/post-pericenter passage), as well as on the still ill-constrained ICM plasma properties (thermal conductivity, viscosity, magnetic field structure). Paper I describes flow patterns and stages of inviscid gas stripping. Here we study the effect of a Spitzer-like temperature dependent viscosity corresponding to Reynolds numbers, Re, of 50 to 5000 with respect to the ICM flow around the remnant atmosphere. Global flow patterns are independent of viscosity in this Reynolds number range. Viscosity influences two aspects: In inviscid stripping, Kelvin-Helmholtz instabilities (KHIs) at the sides of the remnant atmosphere lead to observable horns or wings. Increasing viscosity suppresses KHIs of increasing length scale, and thus observable horns and wings. Furthermore, in inviscid stripping, stripped galactic gas can mix with the ambient ICM in the galaxy's wake. This mixing is suppressed increasingly with increasing viscosity, such that viscously stripped galaxies have long X-ray bright, cool wakes. We provide mock X-ray images for different stripping stages and conditions. While these qualitative results are generic, we tailor our simulations to the Virgo galaxy M89 (NGC 4552), where Re~ 50 corresponds to a viscosity of 10% of the Spitzer level. Paper III compares new deep Chandra and archival XMM-Newton data to our simulations.

preprint2014arXiv

Planck 2013 results. I. Overview of products and scientific results

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

preprint2014arXiv

Planck 2013 results. XXIX. Planck catalogue of Sunyaev-Zeldovich sources

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

preprint2014arXiv

Planck intermediate results. XXVI. Optical identification and redshifts of Planck clusters with the RTT150 telescope

We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with the Russian-Turkish 1.5-m telescope (RTT150), as a part of the optical follow-up programme undertaken by the Planck collaboration. During this time period approximately 20% of all dark and grey clear time available at the telescope was devoted to observations of Planck objects. Some observations of distant clusters were also done at the 6-m Bolshoy Telescope Azimutal'ny (BTA) of the Special Astrophysical Observatory of the Russian Academy of Sciences. In total, deep, direct images of more than one hundred fields were obtained in multiple filters. We identified 47 previously unknown galaxy clusters, 41 of which are included in the Planck catalogue of SZ sources. The redshifts of 65 Planck clusters were measured spectroscopically and 14 more were measured photometrically. We discuss the details of cluster optical identifications and redshift measurements. We also present new spectroscopic redhifts for 39 Planck clusters that were not included in the Planck SZ source catalogue and are published here for the first time.

preprint2014arXiv

Scattering of emission lines in galaxy cluster cores: measuring electron temperature

The central galaxies of some clusters can be strong emitters in the Ly$α$ and H$α$ lines. This emission may arise either from the cool/warm gas located in the cool core of the cluster or from the bright AGN within the central galaxy. The luminosities of such lines can be as high as $10^{42} - 10^{44}$ erg/s. This emission originating from the core of the cluster will get Thomson scattered by hot electrons of the intra-cluster medium (ICM) with an optical depth $\sim$ 0.01 giving rise to very broad ($Δλ/ λ\sim$ 15%) features in the scattered spectrum. We discuss the possibility of measuring the electron density and temperature using information on the flux and width of the highly broadened line features.

preprint2014arXiv

Stellar kinematics of X-ray bright massive elliptical galaxies

We discuss a simple and fast method for estimating masses of early-type galaxies from optical data and compare the results with X-ray derived masses. The optical method relies only on the most basic observables such as the surface brightness $I(R)$ and the line-of-sight velocity dispersion $σ_p(R)$ profiles and provides an anisotropy-independent estimate of the galaxy circular speed $V_c$. The mass-anisotropy degeneracy is effectively overcome by evaluating $V_c$ at a characteristic radius $R_{\rm sweet}$ defined from {\it local} properties of observed profiles. The sweet radius $R_{\rm sweet}$ is expected to lie close to $R_2$, where $I(R) \propto R^{-2}$, and not far from the effective radius $R_{\rm eff}$. We apply the method to a sample of five X-ray bright elliptical galaxies observed with the 6-m telescope BTA-6 in Russia. We then compare the optical $V_c$-estimate with the X-ray derived value, and discuss possible constraints on the non-thermal pressure in the hot gas and configuration of stellar orbits. We find that the average ratio of the optical $V_c$-estimate to the X-ray one is equal to $\approx 0.98$ with $11 \%$ scatter, i.e. there is no evidence for the large non-thermal pressure contribution in the gas at $\sim R_{\rm sweet}$. From analysis of the Lick indices H$β$, Mgb, Fe5270 and Fe5335, we calculate the mass of the stellar component within the sweet radius. We conclude that a typical dark matter fraction inside $R_{\rm sweet}$ in the sample galaxies is $\sim 60\%$ for the Salpeter IMF and $\sim 75 \%$ for the Kroupa IMF.

preprint2014arXiv

Suppression of local heat flux in a turbulent magnetized intracluster medium

X-ray observations of hot gas in galaxy clusters often show steeper temperature gradients across cold fronts -- contact discontinuities, driven by the differential gas motions. These sharp (a few kpc wide) surface brightness/temperature discontinuities would be quickly smeared out by the electron thermal conduction in unmagnetized plasma, suggesting significant suppression of the heat flow across the discontinuities. In fact, the character of the gas flow near cold fronts is favorable for suppression of conduction by aligning magnetic field lines along the discontinuities. We argue that a similar mechanism is operating in the bulk of the gas. Generic 3D random isotropic and incompressible motions increase the temperature gradients (in some places) and at the same time suppress the local conduction by aligning the magnetic field lines perpendicular to the local temperature gradient. We show that the suppression of the effective conductivity in the bulk of the gas can be linked to the increase of the frozen magnetic field energy density. On average the rate of decay of the temperature fluctuations $\mathrm{d}\langle δT^2 \rangle /\mathrm{d}t$ decreases as $\langle B^2 \rangle ^{-1/5}$.

preprint2014arXiv

The relation between gas density and velocity power spectra in galaxy clusters: high-resolution hydrodynamic simulations and the role of conduction

Exploring the ICM power spectrum can help us to probe the physics of galaxy clusters. Using high-resolution 3D plasma simulations, we study the statistics of the velocity field and its relation with the thermodynamic perturbations. The normalization of the ICM spectrum (density, entropy, or pressure) is linearly tied to the level of large-scale motions, which excite both gravity and sound waves due to stratification. For low 3D Mach number M~0.25, gravity waves mainly drive entropy perturbations, traced by preferentially tangential turbulence. For M>0.5, sound waves start to significantly contribute, passing the leading role to compressive pressure fluctuations, associated with isotropic (or slightly radial) turbulence. Density and temperature fluctuations are then characterized by the dominant process: isobaric (low M), adiabatic (high M), or isothermal (strong conduction). Most clusters reside in the intermediate regime, showing a mixture of gravity and sound waves, hence drifting towards isotropic velocities. Remarkably, regardless of the regime, the variance of density perturbations is comparable to the 1D Mach number. This linear relation allows to easily convert between gas motions and ICM perturbations, which can be exploited by Chandra, XMM data and by the forthcoming Astro-H. At intermediate and small scales (10-100 kpc), the turbulent velocities develop a Kolmogorov cascade. The thermodynamic perturbations act as effective tracers of the velocity field, broadly consistent with the Kolmogorov-Obukhov-Corrsin advection theory. Thermal conduction acts to damp the gas fluctuations, washing out the filamentary structures and steepening the spectrum, while leaving unaltered the velocity cascade. The ratio of the velocity and density spectrum thus inverts the downtrend shown by the non-diffusive models, allowing to probe the presence of significant conductivity in the ICM.

preprint2014arXiv

The X-ray coronae of the two brightest galaxies in the Coma cluster

We use deep Chandra X-ray Observatory observations to examine the coronae of the two brightest cluster galaxies in the Coma cluster of galaxies, NGC 4874 and NGC 4889. We find that NGC 4889 hosts a central depression in X-ray surface brightness consistent with a cavity or pair of cavities of radius 0.6 kpc. If the central cavity is associated with an AGN outburst and contains relativistic material, its enthalpy should be around 5x10^55 erg. The implied heating power of this cavity would be around an order of magnitude larger than the energy lost by X-ray emission. It would be the smallest and youngest known cavity in a brightest cluster galaxy and the lack of over pressuring implies heating is still gentle. In contrast, NGC 4874 does not show any evidence for cavities, although it hosts a well-known wide-angle-tail radio source which is visible outside the region occupied by the X-ray corona. These two galaxies show that AGN feedback can behave in varied ways in the same cluster environment.

preprint2014arXiv

Turbulent Heating in Galaxy Clusters Brightest in X-rays

The hot, X-ray-emitting intracluster medium (ICM) is the dominant baryonic constituent of clusters of galaxies. In the cores of many clusters, radiative energy losses from the ICM occur on timescales significantly shorter than the age of the system. Unchecked, this cooling would lead to massive accumulations of cold gas and vigorous star formation, in contradiction to observations. Various sources of energy capable of compensating these cooling losses have been proposed, the most promising being heating by the supermassive black holes in the central galaxies through inflation of bubbles of relativistic plasma. Regardless of the original source of energy, the question of how this energy is transferred to the ICM has remained open. Here we present a plausible solution to this question based on deep Chandra X-ray observatory data and a new data-analysis method that enables us to evaluate directly the ICM heating rate due to the dissipation of turbulence. We find that turbulent heating is sufficient to offset radiative cooling and indeed appears to balance it locally at each radius - it might therefore be the key element in resolving the gas cooling problem in cluster cores and, more universally, in atmospheres of X-ray gas-rich systems.

preprint2013arXiv

Bias from gas inhomogeneities in the pressure profiles as measured from X-ray and SZ observations

X-ray observations of galaxy clusters provide emission measure weighted spectra, arising from a range of density and temperature fluctuations in the intra-cluster medium (ICM). This is fitted to a single temperature plasma emission model to provide an estimate of the gas density and temperature, which are sensitive to the gas inhomogeneities. Therefore, X-ray observations yield a potentially biased estimate of the thermal gas pressure, P_X. At the same time Sunyaev-Zeldovich (SZ) observations directly measure the integrated gas pressure, P_SZ. If the X-ray pressure profiles are strongly biased with respect to the SZ, then one has the possibility to probe the gas inhomogeneities, even at scales unresolved by the current generation of telescopes. At the same time, a weak bias has implications for the interchangeable use of mass proxies like Y_SZ and Y_X as cosmological probes. In this paper we investigate the dependence of the bias, defined as b_P(r)=P_X(r)/P_SZ(r)-1, on the characteristics of fluctuations in the ICM taking into account the correlation between temperature and density fluctuations. We made a simple prediction of the irreducible bias in idealised X-ray vs SZ observations using multi-temperature plasma emission model. We also provide a simple fitting form to estimate the bias given the distribution of fluctuations. Analysing a sample of 16 clusters extracted from hydrodynamical simulations, we find that the median value of bias is within +/-3% within R_500, it decreases to -5% at R_500<r<1.5R_500 and then rises back to ~0% at r >~ 2R_500. The scatter of b_P(r) between individual relaxed clusters is small -- at the level of <0.03 within R_500, but turns significantly larger (0.25) and highly skewed at r>~ 1.5R_500. For any relaxed cluster we find |b_P(r)| < 15% within R_500, across different implementations of input physics in the simulations. [abridged]

preprint2013arXiv

Constraining turbulence and conduction in the hot ICM through density perturbations

Using 3D high-resolution hydrodynamic simulations, tracking both electrons and ions, we study the effects of turbulence and conduction in the hot intracluster medium. We show how the power spectrum of the gas density perturbations can accurately constrain both processes. The characteristic amplitude of density perturbations is linearly related to the strength of turbulence, i.e. the 3D Mach number, as A(k)_max = 0.25 M. The slope of A(k) in turn reflects the level of conduction. In a non-conductive medium, subsonic stirring motions advect density with a similar nearly Kolmogorov cascade. Increasing conduction (via magnetic suppression f = 0.001-1) progressively steepens the spectrum towards the Burgers-like regime. The turbulent Prandtl number defines the dynamic similarity of the flow; at scales where Pt < 100, the spectrum develops a significant decay. The transition is gentle for highly suppressed conduction, f < 0.001, while sharp in the opposite regime. For strong conductivity (f > 0.1), Pt ~ 100 occurs on spatial scales larger than the injection scale. The velocity spectrum is instead not much affected by conduction. The f > 0.1 regime affects the appearance of X-ray images, in which Kelvin-Helmholtz and Rayleigh-Taylor rolls and filaments are washed out. In a stratified system, weak/strong turbulence induces higher isobaric/adiabatic fluctuations, while conduction forces both modes towards the intermediate isothermal regime. We provide a general analytic fit which is applied to new deep Chandra observations of Coma. The observed spectrum is best consistent with strongly suppressed effective isotropic conduction, f = 0.001, and mild subsonic turbulence, M = 0.45 (with injection at ~250 kpc). The low conductivity corroborates the survival of sharp features in the ICM (cold fronts, filaments, bubbles), and indicates that cooling flows may not be balanced by conduction.

preprint2013arXiv

Linear Structures in the Core of the Coma Cluster of Galaxies

The hot X-ray emitting plasma in galaxy clusters is predicted to have turbulent motions which can contribute around ten percent of the cluster's central energy density. We report deep Chandra X-ray Observatory observations of the Coma cluster core, showing the presence of quasi-linear high-density arms spanning 150 kpc, consisting of low-entropy material likely stripped from merging subclusters. Two appear to be connected with a subgroup of galaxies at 650 kpc radius that is merging into the cluster, implying coherence over several hundred Myr. Such long lifetime implies that strong isotropic turbulence and conduction are suppressed in the core, despite the unrelaxed state of the cluster. Magnetic fields are presumably responsible. The structures seen in Coma present insight into the past Gyr of subcluster merger activity.

preprint2013arXiv

Powering of cool filaments in cluster cores by buoyant bubbles. I. Qualitative model

Cool-core clusters (e.g., Perseus or M87) often possess a network of bright gaseous filaments, observed in radio, infrared, optical and X-ray bands. We propose that these filaments are powered by the reconnection of the magnetic field in the wakes of buoyant bubbles. AGN-inflated bubbles of relativistic plasma rise buoyantly in the cluster atmosphere, stretching and amplifying the field in the wake to values of $\displaystyle β=8πP_{gas}/B^2\sim 1$. The field lines in the wake have opposite directions and are forced together as the bubble motion stretches the filament. This setup bears strong similarity to the coronal loops on the Sun or to the Earth's magneto-tail. The reconnection process naturally explains both the required level of local dissipation rate in filaments and the overall luminosity of filaments. The original source of power for the filaments is the potential energy of buoyant bubbles, inflated by the central AGN.

preprint2013arXiv

The Hot and Energetic Universe: AGN feedback in galaxy clusters and groups

Mechanical feedback via Active Galactic Nuclei (AGN) jets in the centres of galaxy groups and clusters is a crucial ingredient in current models of galaxy formation and cluster evolution. Jet feedback is believed to regulate gas cooling and thus star formation in the most massive galaxies, but a robust physical understanding of this feedback mode is currently lacking. The large collecting area, excellent spectral resolution and high spatial resolution of Athena+ will provide the breakthrough diagnostic ability necessary to develop this understanding, via: (1) the first kinematic measurements on relevant spatial scales of the hot gas in galaxy, group and cluster haloes as it absorbs the impact of AGN jets, and (2) vastly improved ability to map thermodynamic conditions on scales well-matched to the jets, lobes and gas disturbances produced by them. Athena+ will therefore determine for the first time how jet energy is dissipated and distributed in group and cluster gas, and how a feedback loop operates in group/cluster cores to regulate gas cooling and AGN fuelling. Athena+ will also establish firmly the cumulative impact of powerful radio galaxies on the evolution of baryons from the epoch of group/cluster formation to the present day.

preprint2012arXiv

A Mexican Hat with holes: calculating low resolution power spectra from data with gaps

A simple method for calculating a low-resolution power spectrum from data with gaps is described. The method is a modification of the $Δ$-variance method previously described by Stutzki and Ossenkopf. A Mexican Hat filter is used to single out fluctuations at a given spatial scale and the variance of the convolved image is calculated. The gaps in the image, defined by the mask, are corrected for by representing the Mexican Hat filter as a difference between two Gaussian filters with slightly different widths, convolving the image and mask with these filters and dividing the results before calculating the final filtered image. This method cleanly compensates for data gaps even if these have complicated shapes and cover a significant fraction of the data. The method was developed to deal with problematic 2D images, where irregular detector edges and masking of contaminating sources compromise the power spectrum estimates, but it can also be straightforwardly applied to 1D timing analysis or 3D data cubes from numerical simulations.

preprint2012arXiv

Constraints on the ICM velocity power spectrum from the X-ray lines width and shift

Future X-ray observations of galaxy clusters by high spectral resolution missions will provide spatially resolved measurements of the energy and width for the brightest emission lines in the intracluster medium (ICM) spectrum. In this paper we discuss various ways of using these high resolution data to constrain velocity power spectrum in galaxy clusters. We argue that variations of these quantities with the projected distance R in cool core clusters contain important information on the velocity field length scales in the ICM. The effective length $l_{\rm eff}$ along the line of sight, which provides dominant contribution to the line flux, increases with R, allowing one to probe the amplitude of the velocity variations at different spatial scales. In particular, we show that the width of the line as a function of R is closely linked to the structure function of the 3D velocity field. Yet another easily obtainable proxy of the velocity field length scales is the ratio of the amplitude of the projected velocity field (line energy) variations to the dispersion of the velocity along the line of sight (line width). Finally the projected velocity field can be easily converted into 3D velocity field, especially for clusters like Coma with an extended flat core in the surface brightness. Under assumption of a homogeneous isotropic Gaussian 3D velocity field we derived simple expressions relating the power spectrum of the 3D velocity field (or structure function) and the observables. The uncertainties in the observables, caused by stochastic nature of the velocity field, are estimated by making multiple realizations of the random Gaussian velocity field and evaluating the scatter in observables. If large scale motions are present in the ICM these uncertainties may dominate the statistical errors of line width and shift measurements.

preprint2012arXiv

Contribution of the accretion disk, hot corona, and obscuring torus to the luminosity of Seyfert galaxies: INTEGRAL and Spitzer observations

We estimate the relative contributions of the supermassive black hole (SMBH) accretion disk, corona, and obscuring torus to the bolometric luminosity of Seyfert galaxies, using Spizter mid-infrared (MIR) observations of a complete sample of 68 nearby active galactic nuclei from the INTEGRAL all-sky hard X-ray (HX) survey. This is the first HX-selected (above 15 keV) sample of AGNs with complementary high angular resolution, high signal to noise, MIR data. Correcting for the host galaxy contribution, we find a correlation between HX and MIR luminosities: L_MIR L_HX^(0.74+/-0.06). Assuming that the observed MIR emission is radiation from an accretion disk reprocessed in a surrounding dusty torus that subtends a solid angle decreasing with increasing luminosity (as inferred from the declining fraction of obscured AGNs), the intrinsic disk luminosity, L_D, is approximately proportional to the luminosity of the corona in the 2-300 keV energy band, L_C, with the L_D/L_C ratio varying by a factor of 2.1 around a mean value of 1.6. This ratio is a factor of ~2 smaller than for typical quasars producing the cosmic X-ray background (CXB). Therefore, over three orders of magnitude in luminosity, HX radiation carries a large, and roughly comparable, fraction of the bolometric output of AGNs. We estimate the cumulative bolometric luminosity density of local AGNs at ~(1-3)x10^40 erg/s/Mpc^3. Finally, the Compton temperature ranges between kT_c~2 and ~6 keV for nearby AGNs, compared to kT_c~2 keV for typical quasars, confirming that radiative heating of interstellar gas can play an important role in regulating SMBH growth.

preprint2012arXiv

Estimating turbulent velocities in the elliptical galaxies NGC 5044 and NGC 5813

The interstellar and intra-cluster medium in giant elliptical galaxies and clusters of galaxies is often assumed to be in hydrostatic equilibrium. Numerical simulations, however, show that about 5-30% of the pressure in a cluster is provided by turbulence induced by, for example, the central AGN and merger activity. We aim to put constraints on the turbulent velocities and turbulent pressure in the ICM of the giant elliptical galaxies NGC 5044 and NGC 5813 using XMM-Newton RGS observations. The magnitude of the turbulence is estimated using the Fe XVII lines at 15.01 A, 17.05 A, and 17.10 A in the RGS spectra. At low turbulent velocities, the gas becomes optically thick in the 15.01 A line due to resonant scattering, while the 17 A lines remain optically thin. By comparing the (I(17.05)+I(17.10))/I(15.01) line ratio from RGS with simulated line ratios for different Mach numbers, the level of turbulence is constrained. The measurement is limited by systematic uncertainties in the atomic data, which are at the 20-30% level. We find that the line ratio in NGC 5813 is significantly higher than in NGC 5044. This difference can be explained by a higher level of turbulence in NGC 5044. The high turbulent velocities and the fraction of the turbulent pressure support of >40% in NGC 5044, assuming isotropic turbulence, confirm that it is a highly disturbed system, probably due to an off-axis merger. The turbulent pressure support in NGC 5813 is more modest at 15-45%. The (I(17.05)+I(17.10))/I(15.01) line ratio in an optically thin plasma, calculated using AtomDB v2.0.1, is 2 sigma above the ratio measured in NGC 5044, which cannot be explained by resonant scattering. This shows that the discrepancies between theoretical, laboratory, and astrophysical data on Fe XVII lines need to be reduced to improve the accuracy of the determination of turbulent velocities using resonant scattering.

preprint2012arXiv

INTEGRAL constraints on the Galactic hard X-ray background from the Milky Way anticenter

We present results of a study of the Galactic ridge X-ray emission (GRXE) in hard X-rays with the IBIS telescope on board INTEGRAL in the region near the Galactic Anticenter (GA) at l=155 deg. We assumed a conservative 2 sigma upper limit on the flux from the GA in the 25-60 keV energy band of 1.25E-10 erg/s/cm^2 (12.8 mCrab) per IBIS field of view, or 6.6E-12 erg/s/cm^2 (0.7 mCrab) per degree longitude in the 135 deg. < l < 175 deg. region. This upper limit exceeds the expected GRXE intensity in the GA direction by an order of magnitude, given the near-infrared (NIR) surface brightness of the Milky Way in this region and the standard hard X-ray-to-NIR intensity ratio for the GRXE, assuming stellar origin. Based on the CGRO/EGRET surface brightness of the Galaxy above 100 MeV as a tracer of the cosmic-ray (CR) induced gamma-ray background, the expected GRXE flux in GA exceeds the measured 2 sigma upper limit by a factor of 4. Therefore, the non-detection of hard X-ray emission from the GA does not contradict the stellar nature of the GRXE, but is inconsistent with CR origin.

preprint2012arXiv

INTEGRAL/IBIS nine-year Galactic Hard X-Ray Survey

Context. The INTEGRAL observatory operating in a hard X-ray/gamma domain has gathered a large observational data set over nine years starting in 2003. Most of the observing time was dedicated to the Galactic source population study, making possible the deepest Galactic survey in hard X-rays ever compiled. Aims. We aim to perform a Galactic survey that can be used as the basis of Galactic source population studies, and perform mapping of the Milky Way in hard X-rays over the maximum exposure available at |b|<17.5 deg. Methods. We used sky reconstruction algorithms especially developed for the high quality imaging of INTEGRAL/IBIS data. Results. We present sky images, sensitivity maps, and catalogs of detected sources in the three energy bands 17-60, 17-35, and 35-80 keV in the Galactic plane at |b|<17.5 deg. The total number of sources in the reference 17-60 keV band includes 402 objects exceeding a 4.7 sigma detection threshold on the nine-year time-averaged map. Among the identified sources with known and tentatively identified natures, 253 are Galactic objects (108 low-mass X-ray binaries, 82 high-mass X-ray binaries, 36 cataclysmic variables, and 27 are of other types), and 115 are extragalactic objects, including 112 active galactic nuclei (AGNs) and 3 galaxy clusters. The sample of Galactic sources with S/N>4.7 sigma has an identification completeness of ~92%, which is valuable for population studies. Since the survey is based on the nine-year sky maps, it is optimized for persistent sources and may be biased against finding transients.

preprint2012arXiv

Kelvin-Helmholtz instabilities at the sloshing cold fronts in the Virgo cluster as a measure for the effective ICM viscosity

Sloshing cold fronts (CFs) arise from minor merger triggered gas sloshing. Their detailed structure depends on the properties of the intra-cluster medium (ICM): hydrodynamical simulations predict the CFs to be distorted by Kelvin-Helmholtz instabilities (KHIs), but aligned magnetic fields, viscosity, or thermal conduction can suppress the KHIs. Thus, observing the detailed structure of sloshing CFs can be used to constrain these ICM properties. Both smooth and distorted sloshing CFs have been observed, indicating that the KHI is suppressed in some clusters, but not in all. Consequently, we need to address at least some sloshing clusters individually before drawing general conclusions about the ICM properties. We present the first detailed attempt to constrain the ICM properties in a specific cluster from the structure of its sloshing CF. Proximity and brightness make the Virgo cluster an ideal target. We combine observations and Virgo-specific hydrodynamical sloshing simulations. Here we focus on a Spitzer-like temperature dependent viscosity as a mechanism to suppress the KHI, but discuss the alternative mechanisms in detail. We identify the CF at 90 kpc north and north-east of the Virgo center as the best location in the cluster to observe a possible KHI suppression. For viscosities $\gtrsim$ 10% of the Spitzer value KHIs at this CF are suppressed. We describe in detail the observable signatures at low and high viscosities, i.e. in the presence or absence of KHIs. We find indications for a low ICM viscosity in archival XMM-Newton data and demonstrate the detectability of the predicted features in deep Chandra observations.

preprint2012arXiv

Planck Intermediate Results. V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect

Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3 x R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from XMM-Newton in the region in which the profiles overlap (i.e., [0.1-1]R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is slightly flatter than most predictions from numerical simulations. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0, c500, gamma, alpha, beta] = [6.41, 1.81, 0.31, 1.33, 4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3 x R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

preprint2012arXiv

Planck intermediate results. X. Physics of the hot gas in the Coma cluster

We present an analysis of Planck satellite data on the Coma Cluster observed via the Sunyaev-Zeldovich effect. Planck is able, for the first time, to detect SZ emission up to r ~ 3 X R_500. We test previously proposed models for the pressure distribution in clusters against the azimuthally averaged data. We find that the Arnaud et al. universal pressure profile does not fit Coma, and that their pressure profile for merging systems provides a good fit of the data only at r<R_500: by r=2XR_500 it underestimates the observed y profile by a factor of ~2. This may indicate that at these larger radii either i) the cluster SZ emission is contaminated by unresolved SZ sources along the line of sight or ii) the pressure profile of Coma is higher at r>R_500 than the mean pressure profile predicted by the simulations. The Planck image shows significant local steepening of the y profile in two regions about half a degree to the west and to the south-east of the cluster centre. These features are consistent with the presence of shock fronts at these radii, and indeed the western feature was previously noticed in the ROSAT PSPC mosaic as well as in the radio. Using Planck y profiles extracted from corresponding sectors we find pressure jumps of 4.5+0.4-0.2 and 5.0+1.3-0.1 in the west and southeast, respectively. Assuming Rankine-Hugoniot pressure jump conditions, we deduce that the shock waves should propagate with Mach number M_w=2.03+0.09-0.04 and M_se=2.05+0.25-0.02 in the West and Southeast, respectively. Finally, we find that the y and radio-synchrotron signals are quasi-linearly correlated on Mpc scales with small intrinsic scatter. This implies either that the energy density of cosmic-ray electrons is relatively constant throughout the cluster, or that the magnetic fields fall off much more slowly with radius than previously thought.

preprint2012arXiv

Quantifying properties of ICM inhomogeneities

We present a new method to identify and characterize the structure of the intracluster medium (ICM) in simulated galaxy clusters. The method uses the median of gas properties, such as density and pressure, which we show to be very robust to the presence of gas inhomogeneities. In particular, we show that the radial profiles of median gas properties are smooth and do not exhibit fluctuations at locations of massive clumps in contrast to mean and mode properties. It is shown that distribution of gas properties in a given radial shell can be well described by a log-normal PDF and a tail. The former corresponds to a nearly hydrostatic bulk component, accounting for ~99% of the volume, while the tail corresponds to high density inhomogeneities. We show that this results in a simple and robust separation of the diffuse and clumpy components of the ICM. The FWHM of the density distribution grows with radius and varies from ~0.15 dex in cluster centre to ~0.5 dex at 2r_500 in relaxed clusters. The small scatter in the width between relaxed clusters suggests that the degree of inhomogeneity is a robust characteristic of the ICM. It broadly agrees with the amplitude of density perturbations in the Coma cluster. We discuss the origin of ICM density variations in spherical shells and show that less than 20% of the width can be attributed to the triaxiality of the cluster gravitational potential. As a link to X-ray observations of real clusters we evaluated the ICM clumping factor with and without high density inhomogeneities. We argue that these two cases represent upper and lower limits on the departure of the observed X-ray emissivity from the median value. We find that the typical value of the clumping factor in the bulk component of relaxed clusters varies from ~1.1-1.2 at r_500 up to ~1.3-1.4 at r_200, in broad agreement with recent observations.

preprint2012arXiv

Testing a simple recipe for estimating galaxy masses from minimal observational data

The accuracy and robustness of a simple method to estimate the total mass profile of a galaxy is tested using a sample of 65 cosmological zoom-simulations of individual galaxies. The method only requires information on the optical surface brightness and the projected velocity dispersion profiles and therefore can be applied even in case of poor observational data. In the simulated sample massive galaxies ($σ\simeq 200-400$ $\kms$) at redshift $z=0$ have almost isothermal rotation curves for broad range of radii (RMS $\simeq 5%$ for the circular speed deviations from a constant value over $0.5R_{\rm eff} < r < 3R_{\rm eff}$). For such galaxies the method recovers the unbiased value of the circular speed. The sample averaged deviation from the true circular speed is less than $\sim 1%$ with the scatter of $\simeq 5-8%$ (RMS) up to $R \simeq 5R_{\rm eff}$. Circular speed estimates of massive non-rotating simulated galaxies at higher redshifts ($z=1$ and $z=2$) are also almost unbiased and with the same scatter. For the least massive galaxies in the sample ($σ< 150$ $\kms$) at $z=0$ the RMS deviation is $\simeq 7-9%$ and the mean deviation is biased low by about $1-2%$. We also derive the circular velocity profile from the hydrostatic equilibrium (HE) equation for hot gas in the simulated galaxies. The accuracy of this estimate is about RMS $\simeq 4-5%$ for massive objects ($M > 6.5\times 10^{12} M_\odot$) and the HE estimate is biased low by $\simeq 3-4%$, which can be traced to the presence of gas motions. This implies that the simple mass estimate can be used to determine the mass of observed massive elliptical galaxies to an accuracy of $5-8 %$ and can be very useful for galaxy surveys.

preprint2012arXiv

X-Ray surface brightness and gas density fluctuations in the Coma cluster

X-ray surface brightness fluctuations in the core ($650 \times 650$ kpc) region of the Coma cluster observed with XMM-Newton and Chandra are analyzed using a 2D power spectrum approach. The resulting 2D spectra are converted to 3D power spectra of gas density fluctuations. Our independent analyses of the XMM-Newton and Chandra observations are in excellent agreement and provide the most sensitive measurements of surface brightness and density fluctuations for a hot cluster. We find that the characteristic amplitude of the volume filling density fluctuations relative to the smooth underlying density distribution varies from 7-10% on scales of $\sim$500 kpc down to $\sim$5% at scales $\sim$ 30 kpc. On smaller spatial scales, projection effects smear the density fluctuations by a large factor, precluding strong limits on the fluctuations in 3D. On the largest scales probed (hundreds of kpc), the dominant contributions to the observed fluctuations most likely arise from perturbations of the gravitational potential by the two most massive galaxies in Coma, NGC4874 and NGC4889, and the low entropy gas brought to the cluster by an infalling group. Other plausible sources of X-ray surface brightness fluctuations are discussed, including turbulence, metal abundance variations, and unresolved sources. Despite a variety of possible origins for density fluctuations, the gas in the Coma cluster core is remarkably homogeneous on scales from $\sim$ 500 to $\sim$30 kpc.

preprint2011arXiv

A wide Chandra view of the core of the Perseus cluster

We present new Chandra images of the X-ray emission from the core of the Perseus cluster of galaxies. The total observation time is now 1.4 Ms. New depressions in X-ray surface brightness are discovered to the north of NGC1275, which we interpret as old rising bubbles. They imply that bubbles are long-lived and do not readily breakup when rising in the hot cluster atmosphere. The existence of a 300 kpc long NNW-SSW bubble axis means there cannot be significant transverse large scale flows exceeding 100 km/s. Interesting spatial correlations are seen along that axis in early deep radio maps. A semi-circular cold front about 100 kpc west of the nucleus is seen. It separates an inner disturbed region dominated by the activity of the active nucleus of NGC1275 from the outer region where a subcluster merger dominates.

preprint2011arXiv

ICM Signatures from Multiple Outbursts in the Galaxy Group NGC 5813: A Window to AGN Feedback

We present results from Chandra observations of NGC 5813, the dominant central galaxy in a nearby galaxy group. We focus on three main results. 1) The diffuse gas shows clear signatures from three distinct outbursts of the central AGN, with three pairs of roughly collinear cavities. The inner two cavity pairs are associated with unambiguous elliptical shock fronts, with Mach numbers M~1.7 and M~1.5 for the inner and outer shocks, respectively. 2) The mean power of the two most recent outbursts differs by a factor of six, indicating that the mean jet power varies over long (~10^7 yr) time scales. 3) The heating from the shocks alone is sufficient to balance radiative cooling of the gas within at least the central 10 kpc, allowing feedback to operate isotropically at small radii.

preprint2011arXiv

Luminosity function of faint Galactic sources in the Chandra bulge field

We study the statistical properties of faint X-ray sources detected in the Chandra Bulge Field. The unprecedented sensitivity of the Chandra observations allows us to probe the population of faint Galactic X-ray sources down to luminosities L(2-10 keV)~1e30 erg/sec at the Galactic Center distance. We show that the luminosity function of these CBF sources agrees well with the luminosity function of sources in the Solar vicinity (Sazonov et al. 2006). The cumulative luminosity density of sources detected in the CBF in the luminosity range 1e30-1e32 erg/sec per unit stellar mass is L(2-10 keV)/M*=(1.7+/-0.3)e27 erg/sec/Msun. Taking into account sources in the luminosity range 1e32-1e34 erg/sec from Sazonov et al. (2006), the cumulative luminosity density in the broad luminosity range 1e30-1e34 erg/sec becomes L(2-10 keV)/M*=(2.4+/-0.4)e27 erg/sec/Msun. The majority of sources with the faintest luminosities should be active binary stars with hot coronae based on the available luminosity function of X-ray sources in the Solar environment.

preprint2011arXiv

ORIGIN: Metal Creation and Evolution from the Cosmic Dawn

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

preprint2011arXiv

Planck Early Results VIII: The all-sky Early Sunyaev-Zeldovich cluster sample

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

preprint2011arXiv

Planck Early Results XI: Calibration of the local galaxy cluster Sunyaev-Zeldovich scaling relations

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters.

preprint2011arXiv

Planck Early Results XXVI: Detection with Planck and confirmation by XMM-Newton of PLCK G266.6-27.3, an exceptionally X-ray luminous and massive galaxy cluster at z~1

We present first results on PLCK G266.6-27.3, a galaxy cluster candidate detected at a signal-to-noise ratio of 5 in the Planck All Sky survey. An XMM-Newton validation observation has allowed us to confirm that the candidate is a bona fide galaxy cluster. With these X-ray data we measure an accurate redshift, z = 0.94 +/- 0.02, and estimate the cluster mass to be M_500 = (7.8 +/- 0.8)e+14 solar masses. PLCK G266.6-27.3 is an exceptional system: its luminosity of L_X(0.5-2.0 keV)=(1.4 +/- 0.05)e+45 erg/s, equals that of the two most luminous known clusters in the z > 0.5 universe, and it is one of the most massive clusters at z~1. Moreover, unlike the majority of high-redshift clusters, PLCK G266.6-27.3 appears to be highly relaxed. This observation confirms Planck's capability of detecting high-redshift, high-mass clusters, and opens the way to the systematic study of population evolution in the exponential tail of the mass function.

preprint2011arXiv

Planck early results: Cluster Sunyaev-Zeldovich optical scaling relations

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters.

preprint2011arXiv

Planck early results. IX. XMM-Newton follow-up for validation of Planck cluster candidates

We present the XMM-Newton follow-up for confirmation of Planck cluster candidates. Twenty-five candidates have been observed to date using snapshot (~10 ksec) exposures, ten as part of a pilot programme to sample a low range of signal-to-noise ratios (4<S/N<6), and a further 15 in a programme to observe a sample of S/N>5 candidates. The sensitivity and spatial resolution of XMM-Newton allows unambiguous discrimination between clusters and false candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates are confirmed as extended X-ray sources. Seventeen are single clusters, the majority of which are found to have highly irregular and disturbed morphologies (about ~70%). The remaining four sources are multiple systems, including the unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to derive a redshift estimate from the X-ray Fe K line (albeit of variable quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a median redshift of z~0.37. A first determination is made of their X-ray properties including the characteristic size, which is used to improve the estimate of the SZ Compton parameter, Y_SZ. The follow-up validation programme has helped to optimise the Planck candidate selection process. It has also provided a preview of the X-ray properties of these newly-discovered clusters, allowing comparison with their SZ properties, and to the X-ray and SZ properties of known clusters observed in the Planck survey. Our results suggest that Planck may have started to reveal a non-negligible population of massive dynamically perturbed objects that is under-represented in X-ray surveys. However, despite their particular properties, these new clusters appear to follow the Y_SZ-Y_X relation established for X-ray selected objects, where Y_X is the product of the gas mass and temperature.

preprint2011arXiv

Planck Early Results. X. Statistical analysis of Sunyaev-Zeldovich scaling relations for X-ray galaxy clusters

All-sky data from the Planck survey and the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC) are combined to investigate the relationship between the thermal Sunyaev-Zeldovich (SZ) signal and X-ray luminosity. The sample comprises ~ 1600 X-ray clusters with redshifts up to ~ 1 and spans a wide range in X-ray luminosity. The SZ signal is extracted for each object individually, and the statistical significance of the measurement is maximised by averaging the SZ signal in bins of X-ray luminosity, total mass, or redshift. The SZ signal is detected at very high significance over more than two decades in X-ray luminosity (10^43 erg/s < L_500 E(z)^-7/3 < 2 X 10^45 erg/s). The relation between intrinsic SZ signal and X-ray luminosity is investigated and the measured SZ signal is compared to values predicted from X-ray data. Planck measurements and X-ray based predictions are found to be in excellent agreement over the whole explored luminosity range. No significant deviation from standard evolution of the scaling relations is detected. For the first time the intrinsic scatter in the scaling relation between SZ signal and X-ray luminosity is measured and found to be consistent with the one in the luminosity -- mass relation from X-ray studies. There is no evidence of any deficit in SZ signal strength in Planck data relative to expectations from the X-ray properties of clusters, underlining the robustness and consistency of our overall view of intra-cluster medium properties.

preprint2011arXiv

Resonant scattering in galaxy clusters for anisotropic gas motions on various spatial scales

The determination of characteristic amplitudes and anisotropy of hot gas motions in galaxy clusters from observations of the brightest resonance lines is discussed. Gas motions affect (i) the spectral line shape through the Doppler effect and (ii) the radial surface brightness profiles in lines during resonant scattering. Radiative transfer calculations have been performed by the Monte Carlo method in the FeXXV resonance line at 6.7 eV for the Perseus cluster (Abell426). It was shown that (a) radial motions reduce the scattering efficiency much more dramatically than purely tangential motions; (b) large-scale gas motions weakly affect the scattering efficiency. The uncertainty in measuring the characteristics of gas motions using resonant scattering has been estimated for existing and future observations of clusters.

preprint2010arXiv

Comparison of an approximately isothermal gravitational potentials of elliptical galaxies based on X-ray and optical data

We analyze six X-ray bright elliptical galaxies, observed with Chandra and XMM-Newton, and approximate their gravitational potentials by isothermal spheres phi(r)=v_c^2 ln(r) over a range of radii from ~0.5 to ~25 kpc. We then compare the circular speed v_c derived from X-ray data with the estimators available from optical data. In particular we discuss two simple and robust procedures for evaluating the circular speed of the galaxy using the observed optical surface brightness and the line-of-sight velocity dispersion profiles. The best fitting relation between the circular speeds derived from optical observations of stars and X-ray observations of hot gas is v_{c,opt}~ η* v_{c,X}, where η=1.10-1.15 (depending on the method), suggesting, albeit with large statistical and systematic uncertainties, that non-thermal pressure on average contributes ~20-30% of the gas thermal pressure.

preprint2010arXiv

Feedback under the microscope II: heating, gas uplift, and mixing in the nearest cluster core

Using a combination of deep 574ks Chandra data, XMM-Newton high-resolution spectra, and optical Halpha+NII images, we study the nature and spatial distribution of the multiphase plasma in M87. Our results provide direct observational evidence of `radio mode' AGN feedback in action, stripping the central galaxy of its lowest entropy gas and preventing star-formation. This low entropy gas was entrained with and uplifted by the buoyantly rising relativistic plasma, forming long "arms". These arms are likely oriented within 15-30 degrees of our line-of-sight. The mass of the uplifted gas in the arms is comparable to the gas mass in the approximately spherically symmetric 3.8 kpc core, demonstrating that the AGN has a profound effect on its immediate surroundings. The coolest X-ray emitting gas in M87 has a temperature of ~0.5 keV and is spatially coincident with Halpha+NII nebulae, forming a multiphase medium where the cooler gas phases are arranged in magnetized filaments. We place strong upper limits of 0.06 Msun/yr on the amount of plasma cooling radiatively from 0.5 keV and show that a uniform, volume-averaged heating mechanism could not be preventing the cool gas from further cooling. All of the bright Halpha filaments appear in the downstream region of the <3 Myr old shock front, at smaller radii than ~0.6'. We suggest that shocks induce shearing around the filaments, thereby promoting mixing of the cold gas with the ambient hot ICM via instabilities. By bringing hot thermal particles into contact with the cool, line-emitting gas, mixing can supply the power and ionizing particles needed to explain the observed optical spectra. Mixing of the coolest X-ray emitting plasma with the cold optical line emitting filamentary gas promotes efficient conduction between the two phases, allowing non-radiative cooling which could explain the lack of X-ray gas with temperatures under 0.5 keV.

preprint2010arXiv

INTEGRAL/IBIS 7-year All-Sky Hard X-Ray Survey. Part I: Image Reconstruction

This paper is the first in a series devoted to the hard X-ray whole sky survey performed by the INTEGRAL observatory over seven years. Here we present an improved method for image reconstruction with the IBIS coded mask telescope. The main improvements are related to the suppression of systematic effects which strongly limit sensitivity in the region of the Galactic Plane (GP), especially in the crowded field of the Galactic Center (GC). We extended the IBIS/ISGRI background model to take into account the Galactic Ridge X-ray Emission (GRXE). To suppress residual systematic artifacts on a reconstructed sky image we applied nonparametric sky image filtering based on wavelet decomposition. The implemented modifications of the sky reconstruction method decrease the systematic noise in the ~20 Ms deep field of GC by ~44%, and practically remove it from the high-latitude sky images. New observational data sets, along with an improved reconstruction algorithm, allow us to conduct the hard X-ray survey with the best currently available minimal sensitivity 3.7E-12 erg/s/cm2 ~0.26 mCrab in the 17-60 keV band at a 5 sigma detection level. The survey covers 90% of the sky down to the flux limit of 6.2E-11 erg/s/cm2 (~4.32 mCrab) and 10% of the sky area down to the flux limit of 8.6E-12 erg/s/cm2 (~0.60 mCrab).

preprint2010arXiv

INTEGRAL/IBIS 7-year All-Sky Hard X-Ray Survey. Part II: Catalog of Sources

This paper is the second in a series devoted to the hard X-ray (17-60 keV) whole sky survey performed by the INTEGRAL observatory over seven years. Here we present a catalog of detected sources which includes 521 objects, 449 of which exceed a 5 sigma detection threshold on the time-averaged map of the sky, and 53 were detected in various subsamples of exposures. Among the identified sources with known and suspected nature, 262 are Galactic (101 low-mass X-ray binaries, 95 high-mass X-ray binaries, 36 cataclysmic variables, and 30 of other types) and 219 are extragalactic, including 214 active galactic nuclei (AGNs), 4 galaxy clusters, and galaxy ESO 389-G 002. The extragalactic (|b|>5 deg) and Galactic (|b|<5 deg) persistently detected source samples are of high identification completeness (respectively ~96% and ~94%) and valuable for population studies.

preprint2010arXiv

Metal transport by gas sloshing in M87

We present the results of an XMM-Newton mosaic covering the central ~200 kpc of the nearby Virgo cluster. We focus on a strong surface brightness discontinuity in the outskirts of the brightest cluster galaxy, M87. Using both XMM-Newton and Suzaku, we derive accurate temperature and metallicity profiles across this feature and show that it is a cold front probably due to sloshing of the Virgo ICM. It is also associated with a discontinuity in the chemical composition. The gas in the inner, bright region of the front is ~40% more abundant in Fe than the gas outside the front, suggesting the important role of sloshing in transporting metals through the ICM. For the first time, we provide a quantitative estimate of the mass of Fe transported by a cold front. This amounts to ~6% of the total Fe mass within the radial range affected by sloshing, significantly more than the amount of metals transported by the AGN in the same cluster core. The very low Fe abundance of only ~0.2 solar immediately outside the cold front at a radius of 90 kpc suggests we are witnessing first-hand the transport of higher metallicity gas into a pristine region, whose abundance is typical of the cluster outskirts. The Mg/Fe and O/Fe abundance ratios remain approximately constant over the entire radial range between the centre of M87 and the faint side of the cold front, which requires the presence of a centrally peaked distribution not only for Fe but also for core-collapse type supernova products. This peak may stem from the star formation triggered as the BCG assembled during the protocluster phase.

preprint2010arXiv

Observational evidence for matter propagation in accretion flows

We study simultaneous X-ray and optical observations of three intermediate polars EX Hya, V1223 Sgr and TV Col with the aim to understand the propagation of matter in their accretion flows. We show that in all cases the power spectra of flux variability of binary systems in X-rays and in optical band are similar to each other and the majority of X-ray and optical fluxes are correlated with time lag <1 sec. These findings support the idea that optical emission of accretion disks, in these binary systems,largely originates as reprocessing of X-ray luminosity of their white dwarfs. In the best obtained dataset of EX Hya we see that the optical lightcurve unambiguously contains some component, which leads the X-ray emission by ~7 sec. We interpret this in the framework of the model of propagating fluctuations and thus deduce the time of travel of the matter from the innermost part of the truncated accretion disk to the white dwarf surface. This value agrees very well with the time expected for matter threaded onto the magnetosphere of the white dwarf to fall to its surface. The datasets of V1223 Sgr and TV Col in general confirm these findings,but have poorer quality.

preprint2010arXiv

Positron annihilation spectrum from the Galactic Centre region observed by SPI/INTEGRAL, revisited: annihilation in a cooling ISM?

We analyse SPI/INTEGRAL data on the 511 keV line from the Galactic Centre, accumulated over ~6 years of observations. We decompose the X-ray and soft gamma-ray emission of the central part of the Milky Way into a relatively compact "Bulge" and a more extended "Disk" components and report their spectral properties. The Bulge component shows a prominent 511 keV line and essentially no flux at 1.8 MeV, while the Disk component on the contrary contains a prominent 1.8 MeV line and a very weak annihilation line. We show that the spectral shape of the annihilation radiation (the narrow 511 keV line and the associated othro-positronium continuum) is surprisingly well described by a model of annihilation of hot positrons in a radiatively cooling interstellar medium (ISM). The model assumes that positrons are initially injected into a hot ($\sim 10^6$~K), volume filling ISM, which is allowed to freely cool via radiative losses. The annihilation time in such a medium is longer than the cooling time for temperatures higher than a few $10^{4}$~K. Thus, most of the positrons annihilate only after the gas has cooled down to $\sim 10^5$~K, giving rise to annihilation emission characteristic of a warm, ionized ISM.

preprint2010arXiv

Search for Outbursts in the Narrow 511-keV Line from Compact Sources Based on INTEGRAL Data

We present the results of a systematic search for outbursts in the narrow positron annihilation line on various time scales (5x10^4 - 10^6 s) based on the SPI/INTEGRAL data obtained from 2003 to 2008. We show that no outbursts were detected with a statistical significance higher than ~6 sigma for any of the time scales considered over the entire period of observations. We also show that, given the large number of independent trials, all of the observed spikes could be associated with purely statistical flux fluctuations and, in part, with a small systematic prediction error of the telescope's instrumental background. Based on the exposure achieved in ~6 yr of INTEGRAL operation, we provide conservative upper limits on the rate of outbursts with a given duration and flux in different parts of the sky.

preprint2010arXiv

Shocks and Cavities from Multiple Outbursts in the Galaxy Group NGC 5813: A Window to AGN Feedback

We present results from new Chandra, GMRT, and SOAR observations of NGC 5813, the dominant central galaxy in a nearby galaxy group. The system shows three pairs of collinear cavities at 1 kpc, 8 kpc, and 20 kpc from the central source, from three distinct outbursts of the central AGN, which occurred 3x10^6, 2x10^7, and 9x10^7 yr ago. The H-alpha and X-ray observations reveal filaments of cool gas that has been uplifted by the X-ray cavities. The inner two cavity pairs are filled with radio emitting plasma, and each pair is associated with an elliptical surface brightness edge, which we unambiguously identify as shocks (with measured temperature jumps) with Mach numbers of M~1.7 and M~1.5 for the inner and outer shocks, respectively. Such clear signatures from three distinct AGN outbursts in an otherwise dynamically relaxed system provide a unique opportunity to study AGN feedback and outburst history. The mean power of the two most recent outbursts differs by a factor of six, from 1.5--10x10^42 erg/s, indicating that the mean jet power changes significantly over long (~10^7 yr) timescales. The total energy output of the most recent outburst is also more than an order of magnitude less than the total energy of the previous outburst (1.5x10^56 erg versus 4x10^57 erg), which may be a result of the lower mean power, or may indicate that the most recent outburst is ongoing. The outburst interval implied by both the shock and cavity ages (~10^7 yr) indicates that, in this system, shock heating alone is sufficient to balance radiative cooling close to the central AGN, which is the relevant region for regulating feedback between the ICM and the central SMBH.

preprint2009arXiv

Constraints on turbulent pressure in the X-ray halos of giant elliptical galaxies from resonant scattering

The dense cores of X-ray emitting gaseous halos of large elliptical galaxies with temperatures below about 0.8 keV show two prominent Fe XVII emission features, which provide a sensitive diagnostic tool to measure the effects of resonant scattering. We present here high-resolution spectra of five bright nearby elliptical galaxies, obtained with the Reflection Grating Spectrometers (RGS) on the XMM-Newton satellite. The spectra for the cores of four of the galaxies show the Fe XVII line at 15.01 Angstrom being suppressed by resonant scattering. The data for NGC 4636 in particular allow the effects of resonant scattering to be studied in detail and to prove that the 15.01 Angstrom line is suppressed only in the dense core and not in the surrounding regions. Using deprojected density and temperature profiles for this galaxy obtained with the Chandra satellite, we model the radial intensity profiles of the strongest resonance lines, accounting for the effects of resonant scattering, for different values of the characteristic turbulent velocity. Comparing the model to the data, we find that the isotropic turbulent velocities on spatial scales smaller than about 1 kpc are less than 100 km/s and the turbulent pressure support in the galaxy core is smaller than 5% of the thermal pressure at the 90% confidence level, and less than 20% at 95% confidence. Neglecting the effects of resonant scattering in spectral fitting of the inner 2 kpc core of NGC 4636 will lead to underestimates of the chemical abundances of Fe and O by ~10-20%.

preprint2009arXiv

Correlation and time delays of the X-ray and optical emission of the Seyfert Galaxy NGC3783

We present simultaneous X-ray and optical B and V band light curves of the Seyfert Galaxy NGC3783 spanning 2 years. The flux in all bands is highly variable and the fluctuations are significantly correlated. As shown before by Stirpe et al. the optical bands vary simultaneously, with a delay of less than 1.5 days but both B and V bands lag the X-ray fluctuations by 3-9 days. This delay points at optical variability produced by X-ray reprocessing and the value of the lag places the reprocessor close to the broad line region. A power spectrum analysis of the light curve, however, shows that the X-ray variability has a power law shape bending to a steeper slope at a time-scale ~2.9 days while the variability amplitude in the optical bands continues to grow towards the longest time-scale covered, ~ 300 days. We show that the power spectra together with the small value of the time delay is inconsistent with a picture where all the optical variability is produced by X-ray reprocessing, though the small amplitude, rapid optical fluctuations might be produced in this way. We detect larger variability amplitudes on long time-scales in the optical bands than in the X-rays. This behaviour adds to similar results recently obtained for at least three other AGN and indicates a separate source of long term optical variability, possibly accretion rate or thermal fluctuations in the optically emitting accretion disc.

preprint2009arXiv

Quenching of the accretion disk strong aperiodic variability at the magnetospheric boundary

We study power density spectra (PDS) of X-ray flux variability in binary systems where the accretion flow is truncated by the magnetosphere. PDS of accreting X-ray pulsars where the neutron star is close to the corotation with the accretion disk at the magnetospheric boundary, have a distinct break/cutoff at the neutron star spin frequency. This break can naturally be explained in the "perturbation propagation" model, which assumes that at any given radius in the accretion disk stochastic perturbations are introduced to the flow with frequencies characteristic for this radius. These perturbations are then advected to the region of main energy release leading to a self-similar variability of X-ray flux P~f^{-1...-1.5}. The break in the PDS is then a natural manifestation of the transition from the disk to magnetospheric flow at the frequency characteristic for the accretion disk truncation radius (magnetospheric radius). The proximity of the PDS break frequency to the spin frequency in corotating pulsars strongly suggests that the typical variability time scale in accretion disks is close to the Keplerian one. In transient accreting X-ray pulsars characterized by large variations of the mass accretion rate during outbursts, the PDS break frequency follows the variations of the X-ray flux, reflecting the change of the magnetosphere size with the accretion rate. Above the break frequency the PDS steepens to ~f^{-2} law which holds over a broad frequency range. These results suggest that strong f^{-1...-1.5} aperiodic variability which is ubiquitous in accretion disks is not characteristic for magnetospheric flows.

preprint2009arXiv

The unusual X-ray morphology of NGC4636 revealed by deep Chandra observations: cavities and shocks created by past AGN outbursts

We present Chandra ACIS-I and ACIS-S observations ($\sim$200 ks in total) of the X-ray luminous elliptical galaxy NGC 4636, located in the outskirts of the Virgo cluster. A soft band (0.5-2 keV) image shows the presence of a bright core in the center surrounded by an extended X-ray corona and two pronounced quasi-symmetric, 8 kpc long, arm-like features. Each of this features defines the rimof an ellipsoidal bubble. An additional bubble-like feature, whose northern rim is located $\sim2$ kpc south of the north-eastern arm, is detected as well. We present surface brightness and temperature profiles across the rims of the bubbles, showing that their edges are sharp and characterized by temperature jumps of about 20-25%. Through a comparison of the observed profiles with theoretical shock models, we demonstrate that a scenario where the bubbles were produced by shocks, probably driven by energy deposited off-center by jets, is the most viable explanation to the X-ray morphology observed in the central part of NGC 4636.

preprint2008arXiv

Discovery of heavily-obscured AGN among 7 INTEGRAL hard X-ray sources observed by Chandra

We observed 7 INTEGRAL sources with the Chandra X-ray Observatory to refine their localization to ~2 arcsec and to study their X-ray spectra. Two sources are inferred to have a Galactic origin: IGR J08390-4833 is most likely a magnetic cataclysmic variable with a white dwarf spin period ~1,450 s; and IGR J21343+4738 is a high-mass X-ray binary. Five sources (IGR J02466-4222, IGR J09522-6231, IGR J14493-5534, IGR J14561-3738, and IGR J23523+5844) prove to be AGN with significant intrinsic X-ray absorption along the line of sight. Their redshifts and hard X-ray (17-60 keV) luminosities range from 0.025 to 0.25 and from ~2x10^43 to ~2x10^45 erg/s, respectively, with the distance to IGR J14493-5534 remaining unknown. The sources IGR J02466-4222 and IGR J14561-3738 are likely Compton-thick AGN with absorption column densities NH>10^24 cm^-2, and the former further appears to be one of the nearest X-ray bright, optically-normal galaxies. With the newly-identified sources, the number of heavily-obscured (NH>10^24 cm^-2) AGN detected by INTEGRAL has increased to ~10. Therefore, such objects constitute 10-15% of hard X-ray bright, non-blazar AGN in the local Universe. The small ratio (<<1%) of soft (0.5-8.0 keV) to hard (17-60 keV) band fluxes (Chandra to INTEGRAL) and the non-detection of optical narrow-line emission in some of the Compton-thick AGN in our sample suggests that there is a new class of objects in which the central massive black hole may be surrounded by a geometrically-thick dusty torus with a narrow ionization cone.

preprint2007arXiv

Filaments, Bubbles, and Weak Shocks in the Gaseous Atmosphere of M87

We present the first results from a 500 ksec Chandra ACIS-I observation of M87. At soft energies (0.5-1.0 keV), we detect filamentary structures associated with the eastern and southwestern X-ray and radio arms. Many filaments are spatially resolved with widths of ~300 pc. This filamentary structure is particularly striking in the eastern arm where we suggest the filaments are outer edges of a series of plasma-filled, buoyant bubbles whose ages differ by ~6 x 10^6 years. These X-ray structures may be influenced by magnetic filamentation. At hard energies (3.5-7.5 keV), we detect a nearly circular ring of outer radius 2.8' (13 kpc) which provides an unambiguous signature of a weak shock, driven by an outburst from the SMBH. The density rise in the shock is ~1.3 (Mach number, M~1.2). The observed spectral hardening in the ring corresponds to a temperature rise T_shock / T_0 ~ 1.2, or M~1.2, in agreement with the Mach number derived independently from the gas density. Thus, for the first time, we detect gas temperature and density jumps associated with a classical shock in the atmosphere around a supermassive black hole. We also detect two additional surface brightness edges and pressure enhancements at radii of ~0.6' and ~1'. The ~0.6' feature may be over-pressurized thermal gas surrounding the relativistic plasma in the radio cocoon, the ``piston'', produced by the current episode of AGN activity. The over-pressurized gas is surrounded by a cool gas shell. The ~1' feature may be an additional weak shock from a secondary outburst. In an earlier episode, the ``piston'' was responsible for driving the 2.8' shock.

preprint2006arXiv

Constraints on the late X-ray emission from the low-energy GRB 031203: INTEGRAL data

Comparison of the INTEGRAL upper limits on the hard X-ray flux before and after the low-energy GRB 031203 with the XMM measurements of the dust-scattered radiation at lower energies suggests that a significant fraction of the total burst energy could be released in the form of soft X-rays at an early afterglow stage with a characteristic duration of ~100-1000s. The overall time evolution of the GRB 031203 afterglow may have not differed qualitatively from the behavior of standard (i.e., more intense) bursts studied by the SWIFT observatory. The available data also admit the possibility that the dust-scattered radiation was associated with an additional soft component in the spectrum of the gamma-ray burst itself.

preprint2006arXiv

Hard X-ray luminosity function and absorption distribution of nearby AGN: INTEGRAL all-sky survey

We use the INTEGRAL all-sky hard X-ray survey to perform a statistical study of a representative sample of nearby AGN. Our entire all-sky sample consists of 127 AGN, of which 91 are confidently detected (>5 sigma) on the time-averaged map obtained with the IBIS/ISGRI instrument and 36 are detected only during single observations. Among the former there are 66 non-blazar AGN located at |b|>5 deg, where the survey's identification completeness is ~93%, which we use for calculating the AGN luminosity function and X-ray absorption distribution. In broad agreement with previous studies, we find that the fraction of obscured (log NH>22) objects is much higher (~70%) among the low-luminosity AGN (Lx<10^43.6 erg/s) than among the high-luminosity ones (Lx>10^43.6 erg/s), \~25%, where Lx is the luminosity in the 17-60 keV energy band. We also find that locally the fraction of Compton-thick AGN is less than 20% unless there is a significant population of AGN that are so strongly obscured that their observed hard X-ray luminosities fall below 10^40-10^41 erg/s, the effective limit of our survey. The constructed hard X-ray luminosity function has a canonical, smoothly broken power-law shape in the range 40<log Lx<45.5 with a characteristic luminosity of log L*=43.40+/-0.28. The estimated local luminosity density due to AGN with log Lx>40 is (1.4+/-0.3) 10^39 erg/s/Mpc^3 (17-60 keV). We demonstrate that the spectral shape and amplitude of the CXB are consistent with the simple scenario in which the NH distribution of AGN (for a given Lx/L*(z) ratio has not changed significantly since z~1.5, while the AGN luminosity function has experienced pure luminosity evolution.

preprint2003arXiv

An XMM-Newton view of the extended "filament" near the cluster of galaxies Abell 85

We have observed with XMM-Newton the extended 4 Mpc filament detected by the ROSAT PSPC in the neighbourhood of the cluster of galaxies Abell 85. We confirm that there is an extended feature, aligned at the same position angle as the major axis of the central cD, the bright cluster galaxies, and nearby groups and clusters. We find that the X-ray emission from the filament is best described by thermal emission with a temperature of about 2 keV, which is significantly lower than the ambient cluster medium, but is significantly higher than anticipated for a gas in a weakly bound extended filament. It is not clear whether this is a filament of diffuse emission, a chain of several groups of galaxies, or stripped gas from the infalling south blob.

preprint2003arXiv

XMM-Newton observations of the Perseus Cluster I: The temperature and surface brightness structure

We present preliminary results of the XMM-Newton 50 ksec observation of the Perseus cluster. The global east/west asymmetry of the gas temperature and surface brightness distributions, approximately aligned with the chain of bright galaxies, suggests an ongoing merger, although the modest degree of the observed asymmetry certainly excludes a major merger interpretation. The chain of galaxies probably traces the filament along which accretion has started some time ago and is continuing at the present time. A cold and dense (low entropy) cluster core like Perseus is probably well "protected" against the penetration of the gas of infalling groups and poor clusters whereas in non-cooling core clusters like Coma and A1367, infalling subclusters can penetrate deeply into the core region. In Perseus, gas associated with infalling groups may be stripped completely at the outskirts of the main cluster and only compression waves (shocks) may reach the central regions. We argue that the passage of such a wave(s) can qualitatively explain the overall horseshoe shaped appearance of the gas temperature map (the hot horseshoe surrounds the colder, low entropy core) as well as other features of the Perseus cluster core. As compression waves traverse the cluster core, they can induce oscillatory motion of the cluster gas which can generate multiple sharp "edges", on opposite sides or the central galaxy. Gas motions induced by mergers may be a natural way to explain the high frequency of "edges" seen in clusters with cooling cores.

preprint2003arXiv

XMM-Newton Observations of the Perseus Cluster II: Evidence for Gas Motions in the Core

The 5-9 keV spectrum of the inner ~100 kpc of the Perseus cluster measured by XMM-Newton can be well described by an optically thin plasma emission model as predicted by the APEC code, without any need for invoking a strong Ni overabundance or the effects of resonant scattering. For the strongest 6.7 keV line of He-like iron, the optical depth of the cluster, calculated using observed density, temperature and abundance profiles, is of order 3. The lack of evidence for resonant scattering effects implies gas motion in the core with a range in velocities of at least half of the sound velocity. If this motion has the character of small scale turbulence, then its dissipation would provide enough energy to compensate for radiative cooling of the gas. The activity of the supermassive black hole at the center of the cluster may be the driving force of the gas motion.

preprint2001arXiv

Evolution of Buoyant Bubbles in M87

The morphology of the X-ray and radio emitting features in the central $\sim$ 50 kpc region around the galaxy M87 strongly suggests that buoyant bubbles of cosmic rays (inflated by an earlier nuclear active phase of the galaxy) rise through the cooling gas at roughly half the sound speed. In the absence of strong surface tension, initially spherical bubbles will transform into tori as they rise through an external medium. Such structures can be identified in the radio images of the halo of M87. During their rise, bubbles will uplift relatively cool X-ray emitting gas from the central regions of the cooling flow to larger distances. This gas is colder than the ambient gas and has a higher volume emissivity. As a result, rising ``radio'' bubbles may be trailed by elongated X-ray features as indeed is observed in M87. We performed simple hydrodynamic simulations to qualitatively illustrate the evolution of buoyant bubbles in the M87 environment.

preprint2001arXiv

Merging Binary Clusters

We study three prominent bi-modal X-ray clusters: A3528, A1750 and A3395. Using observations taken with ROSAT and ASCA, we analyze the temperature and surface brightness distributions. We also analyze the velocity distributions of the three clusters using new measurements supplemented with previously published data. We examine both the overall cluster properties, as well as the two sub-cluster elements in each. These results are then applied to the determination of the overall cluster masses, and demonstrate excellent consistency between the various methods used. While the characteristic parameters of the sub-clusters are typical of isolated objects, our temperature results for the regions between the two sub-clusters clearly confirm the presence of merger activity. These three clusters represent a progression of equal-sized sub-cluster mergers, starting from initial contact to immediately before first core passage.

preprint2001arXiv

Reflection and noise in the low spectral state of GX339-4

We analyze RXTE/PCA observations of GX339-4 in the low spectral state from 1996--1997 and show that the pattern of its spectral and temporal variability is nearly identical to that of Cyg X-1. In particular, a tight correlation exists between the QPO centroid frequency and the spectral parameters. An increase of the QPO centroid frequency is accompanied with an increase of the amplitude of the reflected component and a steepening the slope of the underlying power law. Fourier frequency resolved spectral analysis showed, that the variability of the reflected component at frequencies higher than ~1-10 Hz is suppressed in comparison with that of the primary emission.

preprint1998arXiv

A Hot Spot in Coma

We study the temperature structure of the central part (r<18' ~0.7 h50**-1 Mpc) of the Coma cluster of galaxies using ASCA data. Two different analysis methods produce results in good agreement with each other and reveal the presence of interesting structures in the gas temperature distribution. Globally, the average temperature in the center of the cluster is 9.0 +/- 0.6 keV in good agreement with previous results. Superimposed on this, we find a cool area with temperatures of 4-6 keV associated with a filament of X-ray emission extending southeast from the cluster center detected by Vikhlinin and coworkers. We also find a hot spot with a temperature of around 13 keV displaced north from the central peak of emission. The distribution of the gas temperatures and relative specific entropies suggests that the cool features are most likely gas stripped from a galaxy group centered on NGC 4874 falling toward the core from outside, while the hot spot located ``ahead'' of this in-falling gas is due to shock heating. Thus our results suggest that we are observing Coma during a minor merger with a small group of galaxies associated with NGC 4874 shortly before the initial core passage.

preprint1998arXiv

The millisecond X-ray pulsar/burster SAX J1808.4-3658: the outburst light curve and the power law spectrum

The X-ray light curve and broad band spectral properties of the millisecond X-ray pulsar/burster SAX J1808.4-3658 are reported. In the course of RXTE observations in April--May 1998 the 3--150 keV luminosity of the source decreased by a factor of ~100 from the peak value of ~9E36 erg/s (for a 4 kpc distance). However, the spectrum was remarkably stable and maintained a roughly power law shape with a photon index of ~2 without a strong high energy cut-off below 100 keV, similar to that sometimes observed in the low spectral state of X-ray bursters. An approximation of the averaged spectrum with an exponentially cut-off power law with a superimposed reflected component yields the 90% lower limit on the e-folding energy ~270 keV. We speculate that Comptonization on the bulk motion in the radiation dominated shock might be a possible mechanism of spectral formation. The decaying part of the X-ray light curve was cut off abruptly at luminosity ~few 1E35 erg/s. Such behavior might be due to centrifugal inhibition of accretion (transition to the ``propeller'' regime) as well as to disk instability. In either case an upper limit on the neutron star magnetic field strength is B < few 1E7 Gauss.

preprint1997arXiv

Properties Of The Hard X-ray Radiation From The Black Hole Candidates: Cygnus X-1 And 1E1740.7-2942

The entire dataset of the GRANAT/SIGMA observations of Cyg X-1 and 1E1740.7-2942 in 1990-1994 was analyzed in order to search for correlations between primary observational characteristics of the hard X-ray (40-200 keV) emission - hard X-ray luminosity, hardness of the spectrum (quantified in terms of the best-fit thermal bremsstrahlung temperature kT) and the RMS of short-term flux variations. Although no strict point-to-point correlations were detected certain general tendencies are evident. It was found that for Cyg X-1 the spectral hardness is in general positively correlated with relative amplitude of short-term variability. The correlation of similar kind was found for X-ray transient GRO J0422+32 (X-ray Nova Persei 1992). For both sources an approximate correlation between kT and L_X was found. At low hard X-ray luminosity - below 10E37 erg/sec - kT increases with L_X. At higher luminosity the spectral hardness depends weaker or does not depend at all on the hard X-ray luminosity. The low luminosity end of these approximate correlations (low kT and low RMS) corresponds to extended episodes of very low hard X-ray flux occurred during SIGMA observations.