Deep Chandra observation and numerical studies of the nearest cluster cold front in the sky

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M 87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4–8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin–Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are ∼10 per cent brighter than the surrounding gas and are separated by ∼15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches ∼5–10 per cent of the thermal pressure, reducing the gas density between the bright features

Galaxy Clusters as Reservoirs of Heavy Dark Matter and High-Energy Cosmic Rays: Constraints from Neutrino Observations

Galaxy Clusters (GCs) are the largest reservoirs of both dark matter and cosmic rays (CRs). Dark matter self-annihilation can lead to a high luminosity in gamma rays and neutrinos, enhanced by a strong degree of clustering in dark matter substructures. Hadronic CR interactions can also lead to a high luminosity in gamma rays and neutrinos, enhanced by the confinement of CRs from cluster accretion/merger shocks and active galactic nuclei. We show that IceCube/KM3Net observations of high-energy neutrinos can probe the nature of GCs and the separate dark matter and CR emission processes, taking into account how the results depend on the still-substantial uncertainties. Neutrino observations are relevant at high energies, especially at >10 TeV. Our results should be useful for improving experimental searches for high-energy neutrino emission. Neutrino telescopes are sensitive to extended sources formed by dark matter substructures and CRs distributed over large scales. Recent observations by Fermi and imaging atmospheric Cherenkov telescopes have placed interesting constraints on the gamma-ray emission from GCs. We also provide calculations of the gamma-ray fluxes, taking into account electromagnetic cascades inside GCs, which can be important for injections at sufficiently high energies. This also allows us to extend previous gamma-ray constraints to very high dark matter masses and significant CR injections at very high energies. Using both neutrinos and gamma rays, which can lead to comparable constraints, will allow more complete understandings of GCs. Neutrinos are essential for some dark matter annihilation channels, and for hadronic instead of electronic CRs. Our results suggest that the multi-messenger observations of GCs will be able to give useful constraints on specific models of dark matter and CRs. [Abstract abridged.]Comment: 31 pages, 20 figures, 1 table, accepted for publication in JCAP, references and discussions adde

Unidentified EGRET Sources and the Extragalactic Gamma-Ray Background

The large majority of EGRET point sources remain to this day without an identified low-energy counterpart. Whatever the nature of the EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse gamma-ray background: if most unidentified objects are extragalactic, faint unresolved sources of the same class contribute to the background, as a distinct extragalactic population; on the other hand, if most unidentified sources are Galactic, their counterparts in external galaxies will contribute to the unresolved emission from these systems. Understanding this component of the gamma-ray background, along with other guaranteed contributions from known sources, is essential in any attempt to use gamma-ray observations to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the extragalactic gamma-ray background is likely to be important, and we find that it is. Additionally, we comment on how the anticipated GLAST measurement of the diffuse gamma-ray background will change, depending on the nature of the majority of these sources.Comment: 6 pages, 3 figures, to appear in proceedings of "The Multi-Messenger Approach to High Energy Gamma-Ray Sources", Barcelona, 4-7 July 2006; comments welcom

The upstream magnetic field of collisionless GRB shocks: constraint by Fermi-LAT observations

Long-lived >100 MeV emission has been a common feature of most Fermi-LAT detected gamma-ray bursts (GRBs), e.g., detected up to ~10^3s in long GRBs 080916C and 090902B and ~10^2s in short GRB 090510. This emission is consistent with being produced by synchrotron emission of electrons accelerated to high energy by the relativistic collisionless shock propagating into the weakly magnetized medium. Here we show that this high-energy afterglow emission constrains the preshock magnetic field to satisfy 1(n/1cc)^{9/8} mG<B<10^2(n/1cc)^{3/8}mG, where n is the preshock density, more stringent than the previous constraint by X-ray afterglow observations on day scale. This suggests that the preshock magnetic field is strongly amplified, most likely by the streaming of high energy shock accelerated particles.Comment: 9 pages, JCAP accepte

Probing mSUGRA via the Extreme Universe Space Observatory

An analysis is carried out within mSUGRA of the estimated number of events originating from upward moving ultra-high energy neutralinos that could be detected by the Extreme Universe Space Observatory (EUSO). The analysis exploits a recently proposed technique that differentiates ultra-high energy neutralinos from ultra-high energy neutrinos using their different absorption lengths in the Earth's crust. It is shown that for a significant part of the parameter space, where the neutralino is mostly a Bino and with squark mass $\sim 1$ TeV, EUSO could see ultra-high energy neutralino events with essentially no background. In the energy range 10^9 GeV < E < 10^11 GeV, the unprecedented aperture of EUSO makes the telescope sensitive to neutralino fluxes as low as 1.1 \times 10^{-6} (E/GeV)^{-1.3} GeV^{-1} cm^{-2} yr^{-1} sr^{-1}, at the 95% CL. Such a hard spectrum is characteristic of supermassive particles' $N$-body hadronic decay. The case in which the flux of ultra-high energy neutralinos is produced via decay of metastable heavy particles with uniform distribution throughout the universe is analyzed in detail. The normalization of the ratio of the relics' density to their lifetime has been fixed so that the baryon flux produced in the supermassive particle decays contributes to about 1/3 of the events reported by the AGASA Collaboration below 10^{11} GeV, and hence the associated GeV gamma-ray flux is in complete agreement with EGRET data. For this particular case, EUSO will collect between 4 and 5 neutralino events (with 0.3 of background) in ~ 3 yr of running. NASA's planned mission, the Orbiting Wide-angle Light-collectors (OWL), is also briefly discussed in this context.Comment: Some discussion added, final version to be published in Physical Review

Squeeze-and-Breathe Evolutionary Monte Carlo Optimisation with Local Search Acceleration and its application to parameter fitting

Motivation: Estimating parameters from data is a key stage of the modelling process, particularly in biological systems where many parameters need to be estimated from sparse and noisy data sets. Over the years, a variety of heuristics have been proposed to solve this complex optimisation problem, with good results in some cases yet with limitations in the biological setting. Results: In this work, we develop an algorithm for model parameter fitting that combines ideas from evolutionary algorithms, sequential Monte Carlo and direct search optimisation. Our method performs well even when the order of magnitude and/or the range of the parameters is unknown. The method refines iteratively a sequence of parameter distributions through local optimisation combined with partial resampling from a historical prior defined over the support of all previous iterations. We exemplify our method with biological models using both simulated and real experimental data and estimate the parameters efficiently even in the absence of a priori knowledge about the parameters.Comment: 15 Pages, 3 Figures, 6 Tables; Availability: Matlab code available from the authors upon reques

High-Energy Cosmology: gamma rays and neutrinos from beyond the galaxy

Our knowledge of the high-energy universe is undergoing a period of rapid change as new astronomical detectors of high-energy radiation start to operate at their design sensitivities. Now is a boomtime for high-energy astrophysics, with new discoveries from Swift and HESS, results from MAGIC and VERITAS starting to be reported, the upcoming launches of the gamma-ray space telescopes GLAST and AGILE, and anticipated data releases from IceCube and Auger. A formalism for calculating statistical properties of cosmological gamma-ray sources is presented. Application is made to model calculations of the statistical distributions of gamma-ray and neutrino emission from (i) beamed sources, specifically, long-duration GRBs, blazars, and extagalactic microquasars, and (ii) unbeamed sources, including normal galaxies, starburst galaxies and clusters. Expressions for the integrated intensities of faint beamed and unbeamed high-energy radiation sources are also derived. A toy model for the background intensity of radiation from dark-matter annihilation taking place in the early universe is constructed. Estimates for the gamma-ray fluxes of local group galaxies, starburst, and infrared luminous galaxies are briefly reviewed. Because the brightest extragalactic gamma-ray sources are flaring sources, and these are the best targets for sources of PeV -- EeV neutrinos and ultra-high energy cosmic rays, rapidly slewing all-sky telescopes like MAGIC and an all-sky gamma-ray observatory beyond Milagro will be crucial for optimal science return in the multi-messenger age.Comment: 10 pages, 3 figs, accepted for publication in the Barcelona Conference on Multimessenger Astronomy; corrected eq. 27, revised Fig. 3, added 2 ref

Gamma-Ray Luminosity Function of Blazars and the Cosmic Gamma-Ray Background: Evidence for the Luminosity Dependent Density Evolution

We present a comprehensive study for the gamma-ray luminosity function (GLF) of blazars and their contribution to the extragalactic diffuse gamma-ray background (EGRB). Radio and gamma-ray luminosity correlation is introduced to take into account the radio detectability which is important for the blazar identification. Previous studies considered only pure luminosity evolution (PLE) or pure density evolution, but we introduce the luminosity dependent density evolution (LDDE) model, which is favored from the evolution of X-ray luminosity function (XLF) of AGNs. The model parameters are constrained by likelihood analyses about the observed redshift and gamma-ray flux distributions of the EGRET blazars. We find that the LDDE model gives a better fit to the observed distributions than the PLE model, indicating that the LDDE model is also appropriate for gamma-ray blazars, and that the jet activity is universally correlated with the accretion history of AGNs. The normalization between the GLF and XLF is consistent with the unified picture of AGNs, when the beaming and a reasonable duty cycle of jet activity are taken into account. We then find that only 25--50% of the EGRB can be explained by unresolved blazars with the best-fit LDDE parameters. Unresolved blazars can account for all the EGRB only with a steeper index of the faint-end slope of the GLF, which is marginally consistent with the EGRET data but inconsistent with that of the XLF. Therefore unresolved AGNs cannot be the dominant source of the EGRB, unless there is a new population of gamma-ray emitting AGNs that evolves differently from the XLF of AGNs. Predictions for the GLAST mission are made, and we find that the best-fit LDDE model predicts about 3000 blazars in the entire sky, which is considerably fewer than a previous estimate.Comment: 13 pages, 12 figures, accepted by ApJ; minor typos corrected and some figures revised, main conclusions essentially unchange

Statistical Mechanics of Canonical-Dissipative Systems and Applications to Swarm Dynamics

We develop the theory of canonical-dissipative systems, based on the assumption that both the conservative and the dissipative elements of the dynamics are determined by invariants of motion. In this case, known solutions for conservative systems can be used for an extension of the dynamics, which also includes elements such as the take-up/dissipation of energy. This way, a rather complex dynamics can be mapped to an analytically tractable model, while still covering important features of non-equilibrium systems. In our paper, this approach is used to derive a rather general swarm model that considers (a) the energetic conditions of swarming, i.e. for active motion, (b) interactions between the particles based on global couplings. We derive analytical expressions for the non-equilibrium velocity distribution and the mean squared displacement of the swarm. Further, we investigate the influence of different global couplings on the overall behavior of the swarm by means of particle-based computer simulations and compare them with the analytical estimations.Comment: 14 pages incl. 13 figures. v2: misprints in Eq. (40) corrected, ref. updated. For related work see also: http://summa.physik.hu-berlin.de/~frank/active.htm