Polarization of MeV gamma-rays and 511 keV line shape as probes of SNIa asymmetry and magnetic field

We discuss gamma-ray signatures associated with an asymmetric explosion and transport of positrons in SN Ia ejecta. In particular, Compton scattering of gamma-ray line photons can induce polarization in the continuum, which would be a direct probe of the asymmetries in the distribution of radioactive isotopes and/or of the scattering medium. Even more interesting would be a comparison of the shapes of $\gamma$-ray lines and that of the electron-positron annihilation line at 511 keV. The shapes of $\gamma$-ray lines associated with the decay of Co56 (e.g., lines at 847 and 1238 keV) directly reflect the velocity distribution of Co56. On the other hand, the 511 keV line arises from the annihilation of positions, which are also produced by the Co56 decay but can propagate through the ejecta before they slow down and annihilate. Therefore, the shape of the annihilation line might differ from other gamma-ray lines, providing constraints on the efficiency of positrons propagation through the ejecta and, as consequence, on the topology of magnetic fields in the ejecta and on the fraction of positrons that escape to the interstellar medium. We illustrate the above effects with two models aimed at capturing the main predicted signatures.Comment: 10 pages, 10 figures; replaced with accepted version (MNRAS

Sound wave generation by a spherically symmetric outburst and AGN Feedback in Galaxy Clusters

We consider the evolution of an outburst in a uniform medium under spherical symmetry, having in mind AGN feedback in the intra cluster medium (ICM). For a given density and pressure of the medium, the spatial structure and energy partition at a given time $t_{age}$ (since the onset of the outburst) are fully determined by the total injected energy $E_{inj}$ and the duration $t_b$ of the outburst. We are particularly interested in the late phase evolution when the strong shock transforms into a sound wave. We studied the energy partition during such transition with different combinations of $E_{inj}$ and $t_b$. For an instantaneous outburst with $t_b\rightarrow 0$, which corresponds to the extension of classic Sedov-Taylor solution with counter-pressure, the fraction of energy that can be carried away by sound waves is $\lesssim$12% of $E_{inj}$. As $t_b$ increases, the solution approaches the "slow piston" limit, with the fraction of energy in sound waves approaching zero. We then repeat the simulations using radial density and temperature profiles measured in Perseus and M87/Virgo clusters. We find that the results with a uniform medium broadly reproduce an outburst in more realistic conditions once proper scaling is applied. We also develop techniques to map intrinsic properties of an outburst $(E_{inj}, t_b$ and $t_{age})$ to the observables like the Mach number of the shock and radii of the shock and ejecta. For the Perseus cluster and M87, the estimated $(E_{inj}, t_b$ and $t_{age})$ agree with numerical simulations tailored for these objects with $20-30\%$ accuracy.Comment: Accepted by MNRAS, add one figure in appendix and minor changes in text based on referee's commen

Mass density slope of elliptical galaxies from strong lensing and resolved stellar kinematics

We discuss constraints on the mass density distribution (parameterized as $\rho\propto r^{-\gamma}$) in early-type galaxies provided by strong lensing and stellar kinematics data. The constraints come from mass measurements at two `pinch' radii. One `pinch' radius $r_1=2.2 R_{Einst}$ is defined such that the Einstein (i.e. aperture) mass can be converted to the spherical mass almost independently of the mass-model. Another `pinch' radius $r_2=R_{opt}$ is chosen so that the dynamical mass, derived from the line-of-sight velocity dispersion, is least sensitive to the anisotropy of stellar orbits. We verified the performance of this approach on a sample of simulated elliptical galaxies and on a sample of 15 SLACS lens galaxies at $0.01 \leq z \leq 0.35$, which have already been analysed in Barnabe et al. (2011) by the self-consistent joint lensing and kinematic code. For massive simulated galaxies the density slope $\gamma$ is recovered with an accuracy of $\sim 13\%$, unless $r_1$ and $r_2$ happen to be close to each other. For SLACS galaxies, we found good overall agreement with the results of Barnabe et al. (2011) with a sample-averaged slope $\gamma=2.1\pm0.05$. While the two-pinch-radii approach has larger statistical uncertainties, it is much simpler and uses only few arithmetic operations with directly observable quantities.Comment: accepted for publication in MNRA

Does the obscured AGN fraction really depend on luminosity?

We use a sample of 151 local non-blazar AGN selected from the INTEGRAL all-sky hard X-ray survey to investigate if the observed declining trend of the fraction of obscured (i.e. showing X-ray absorption) AGN with increasing luminosity is mostly an intrinsic or selection effect. Using a torus-obscuration model, we demonstrate that in addition to negative bias, due to absorption in the torus, in finding obscured AGN in hard X-ray flux limited surveys, there is also positive bias in finding unobscured AGN, due to Compton reflection in the torus. These biases can be even stronger taking into account plausible intrinsic collimation of hard X-ray emission along the axis of the obscuring torus. Given the AGN luminosity function, which steepens at high luminosities, these observational biases lead to a decreasing observed fraction of obscured AGN with increasing luminosity even if this fraction has no intrinsic luminosity dependence. We find that if the central hard X-ray source in AGN is isotropic, the intrinsic (i.e. corrected for biases) obscured AGN fraction still shows a declining trend with luminosity, although the intrinsic obscured fraction is significantly larger than the observed one: the actual fraction is larger than $\sim 85$% at $L\lesssim 10^{42.5}$ erg/s (17--60 keV), and decreases to $\lesssim 60$% at $L\gtrsim 10^{44}$ erg/s. In terms of the half-opening angle, $\theta$, of an obscuring torus, this implies that $\theta\lesssim 30$ deg in lower-luminosity AGN, and $\theta\gtrsim 45$ deg in higher-luminosity ones. If, however, the emission from the central SMBH is collimated as $dL/d\Omega\propto\cos\alpha$, the intrinsic dependence of the obscured AGN fraction is consistent with a luminosity-independent torus half-opening angle $\theta\sim 30$ deg.Comment: 20 pages, 18 figures, accepted for publication in MNRA

Non-Detection of X-Ray Emission From Sterile Neutrinos in Stacked Galaxy Spectra

We conduct a comprehensive search for X-ray emission lines from sterile neutrino dark matter, motivated by recent claims of unidentified emission lines in the stacked X-ray spectra of galaxy clusters and the centers of the Milky Way and M31. Since the claimed emission lines lie around 3.5 keV, we focus on galaxies and galaxy groups (masking the central regions), since these objects emit very little radiation above $\sim 2$ keV and offer a clean background against which to detect emission lines. We develop a formalism for maximizing the signal-to-noise of decaying dark matter emission lines by weighing each X-ray event according to the expected dark matter profile. In total, we examine 81 and 89 galaxies with Chandra and XMM-Newton respectively, totaling 15.0 and 14.6 Ms of integration time. We find no significant evidence of any emission lines, placing strong constraints on the mixing angle of sterile neutrinos with masses between 4.8-12.4 keV. In particular, if the 3.57 keV feature from Bulbul et al. (2014) were due to 7.1 keV sterile neutrino emission, we would have detected it at $4.4\sigma$ and $11.8\sigma$ in our two samples. The most conservative estimates of the systematic uncertainties reduce these constraints to $4.4\sigma$ and 7.8$\sigma$, or letting the line energy vary between 3.50 and 3.60 keV reduces these constraints to $2.7\sigma$ and $11.0\sigma$ respectively. Unlike previous constraints, our measurements do not depend on the model of the X-ray background or on the assumed logarithmic slope of the center of the dark matter profile.Comment: accepted to MNRA

Sound wave generation by a spherically symmetric outburst and AGN feedback in galaxy clusters II: impact of thermal conduction

We analyze the impact of thermal conduction on the appearance of a shock-heated gas shell which is produced when a spherically symmetric outburst of a supermassive black hole inflates bubbles of relativistic plasma at the center of a galaxy cluster. The presence of the hot and low-density shell can be used as an ancillary indicator for a high rate of energy release during the outburst, which is required to drive strong shocks into the gas. Here we show that conduction can effectively erase such shell, unless the diffusion of electrons is heavily suppressed. We conclude that a more robust proxy to the energy release rate is the ratio between the shock radius and bubble radius. We also revisited the issue of sound waves dissipation induced by thermal conduction in a scenario, where characteristic wavelength of the sound wave is set by the total energy of the outburst. For a fiducial short outburst model, the dissipation length does not exceed the cooling radius in a typical cluster, provided that the conduction is suppressed by a factor not larger than $\sim$100. For quasi-continuous energy injection neither the shock-heated shell nor the outgoing sound wave are important and the role of conduction is subdominant.Comment: 12 pages, 10 figures, MNRAS in pres

X-ray emission from warm-hot intergalactic medium: the role of resonantly scattered cosmic X-ray background

We revisit calculations of the X-ray emission from warm-hot intergalactic medium (WHIM) with particular focus on contribution from the resonantly scattered cosmic X-ray background (CXB). If the significant part of the CXB emission is resolved into point sources, the properties of the WHIM along the line of sight are recorded in the absorption lines in the stacked spectrum of resolved sources and in the emission lines in the remaining diffuse signal. For the strongest resonant lines, this implies a factor of $\sim30$ boost in emissivity compared to the intrinsic emissivity over the major part of the density-temperature parameter space region relevant for WHIM. The overall boost for the 0.5-1 keV band is $\sim4$, declining steeply at temperatures above $10^{6}$ K and over-densities $\delta\gtrsim100$. In addition to the emissivity boost, contribution of the resonant scattering changes relative intensities of the lines, so it should be taken into account when line-ratio-diagnostics from high resolution spectra or redshift determination from low resolution spectra are considered. Comparison between WHIM signatures in X-ray absorption and emission should allow differentiating truly diffuse gas of small overdensity from denser clumps having small filling factor by future X-ray missions.Comment: 15 pages, 13 figures, accepted for publication in MNRA

Standoff Distance of Bow Shocks in Galaxy Clusters as Proxy for Mach Number

X-ray observations of merging clusters provide many examples of bow shocks leading merging subclusters. While the Mach number of a shock can be estimated from the observed density jump using Rankine-Hugoniot condition, it reflects only the velocity of the shock itself and is generally not equal to the velocity of the infalling subcluster dark matter halo or to the velocity of the contact discontinuity separating gaseous atmospheres of the two subclusters. Here we systematically analyze additional information that can be obtained by measuring the standoff distance, i.e. the distance between the leading edge of the shock and the contact discontinuity that drives this shock. The standoff distance is influenced by a number of additional effects, e.g. (1) the gravitational pull of the main cluster (causing acceleration/deceleration of the infalling subcluster), (2) the density and pressure gradients of the atmosphere in the main cluster, (3) the non-spherical shape of the subcluster, and (4) projection effects. The first two effects tend to bias the standoff distance in the same direction, pushing the bow shock closer to (farther away from) the subcluster during the pre- (post-)merger stages. Particularly, in the post-merger stage, the shock could be much farther away from the subcluster than predicted by a model of a body moving at a constant speed in a uniform medium. This implies that a combination of the standoff distance with measurements of the Mach number from density/temperature jumps can provide important information on the merger, e.g. differentiating between the pre- and post-merger stages.Comment: 11 pages, 12 figures. Including major revision and matched to accepted version in MNRA

Runaway Merger Shocks in Galaxy Cluster Outskirts and Radio Relics

Moderately strong shocks arise naturally when two subclusters merge. For instance, when a smaller subcluster falls into the gravitational potential of a more massive cluster, a bow shock is formed and moves together with the subcluster. After pericenter passage, however, the subcluster is decelerated by the gravity of the main cluster, while the shock continues moving away from the cluster center. These shocks are considered as promising candidates for powering radio relics found in many clusters. The aim of this paper is to explore the fate of such shocks when they travel to the cluster outskirts, far from the place where the shocks were initiated. In a uniform medium, such a "runaway" shock should weaken with distance. However, as shocks move to large radii in galaxy clusters, the shock is moving down a steep density gradient that helps the shock to maintain its strength over a large distance. Observations and numerical simulations show that, beyond $R_{500}$, gas density profiles are as steep as, or steeper than, $\sim r^{-3}$, suggesting that there exists a "Habitable zone" for moderately strong shocks in cluster outskirts where the shock strength can be maintained or even amplified. A characteristic feature of runaway shocks is that the strong compression, relative to the initial state, is confined to a narrow region just behind the shock. Therefore, if such a shock runs over a region with a pre-existing population of relativistic particles, then the boost in radio emissivity, due to pure adiabatic compression, will also be confined to a narrow radial shell.Comment: 9 pages, 8 figures; published in MNRA