Showing 1–50 of 149 results for author: Aharonian, F A

Particle acceleration in shearing flows: the self-generation of turbulent spine-sheath structures in relativistic MHD jet simulations

Authors: Jie-Shuang Wang, Brian Reville, Yosuke Mizuno, Frank M. Rieger, Felix A. Aharonian

Abstract: X-ray observations of several kiloparsec-scale extragalactic jets favour a synchrotron origin. The short cooling times of the emitting electrons requires distributed acceleration of electrons up to sub-PeV energies. In a previous paper, we found that this can be self-consistently explained by a shear acceleration model, where particles are accelerated to produce power-law spectra with a spectral i… ▽ More X-ray observations of several kiloparsec-scale extragalactic jets favour a synchrotron origin. The short cooling times of the emitting electrons requires distributed acceleration of electrons up to sub-PeV energies. In a previous paper, we found that this can be self-consistently explained by a shear acceleration model, where particles are accelerated to produce power-law spectra with a spectral index being determined mainly by the velocity profile and turbulence spectrum. In this paper, we perform 3D relativistic magneto-hydrodynamic simulations to investigate the formation of a spine-sheath structure and the development of turbulence for a relativistic jet propagating into a static cocoon. We explore different spine velocities and magnetic field profiles with values being chosen to match typical Fanaroff-Riley type I/II jets. We find that in all cases a sheath is generated on the interface of the spine and the cocoon mainly due to the Kelvin-Helmholtz instability. The large scale velocity profile in the sheath is close to linear. Turbulence develops in both the spine and the sheath, with a turbulent velocity spectrum consistent with Kolmogorov-scaling. The implications for shear particle acceleration are explored, with a focus on the particle spectral index. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 9 pages, 11 figures, 1 table, accepted for publication in MNRAS

Does the Geminga, Monogem and PSR J0622+3749 $γ$-ray halos imply slow diffusion around pulsars?

Authors: Sarah Recchia, Mattia Di Mauro, Felix A. Aharonian, Luca Orusa, Fiorenza Donato, Stefano Gabici, Silvia Manconi

Abstract: The HAWC Collaboration has reported the detection of an extended $γ$-ray emission around the Geminga and Monogem pulsars of a few degree extension. Very recently, the LHAASO Collaboration released also the data for an extended $γ$-ray emission around the pulsar PSR J0622+3749. This flux can be explained with electrons and positrons injected from these sources and their inverse Compton Scattering o… ▽ More The HAWC Collaboration has reported the detection of an extended $γ$-ray emission around the Geminga and Monogem pulsars of a few degree extension. Very recently, the LHAASO Collaboration released also the data for an extended $γ$-ray emission around the pulsar PSR J0622+3749. This flux can be explained with electrons and positrons injected from these sources and their inverse Compton Scattering on the interstellar radiation fields. So far the size of such $γ-$ray halos has been interpreted as the result of the diffusion coefficient around the sources being about two orders of magnitude smaller than the average in the Galaxy. However, this conclusion is driven by the assumption that particles propagate diffusively right away after the injection without taking into account the ballistic propagation. The propagation of cosmic-ray leptons in the proximity of the Geminga, Monogem and PSR J0622+3749 pulsars is examined here considering the transition from the quasi-ballistic, valid for the most recently injected particles, to the diffusive transport regime. For typical interstellar values of the diffusion coefficient, the quasi-ballistic regime dominates the lepton distribution up to distances of a few tens of parsec from the pulsar for particle energies above $\sim 10$ TeV. In this regime the resulting $γ-$ray source tends to be rather compact, despite particles travel a long distance. Indeed, for larger values of the diffusion coefficient, particles propagate ballistically up to larger distances with the result of a more point-like $γ-$ray source. When such transition is taken into account, a good fit to the HAWC and LHAASO $γ-$ray data around Geminga, Monogem and PSR J0622+3749 is obtained without the need to invoke a strong suppression of the diffusion coefficient. △ Less

Submitted 22 October, 2021; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 10 pages, 4 figures. This version includes other two sources (Monogem and PSR J0622+3749) with respect to the first version

Particle acceleration in shearing flows: the case for large-scale jets

Authors: Jie-Shuang Wang, Brian Reville, Ruo-Yu Liu, Frank M. Rieger, Felix A. Aharonian

Abstract: X-ray observations of kilo-parsec scale jets indicate that a synchrotron origin of the sustained non-thermal emission is likely. This requires distributed acceleration of electrons up to near PeV energies along the jet. The underlying acceleration mechanism is still unclear. Shear acceleration is a promising candidate, as velocity-shear stratification is a natural consequence of the collimated flo… ▽ More X-ray observations of kilo-parsec scale jets indicate that a synchrotron origin of the sustained non-thermal emission is likely. This requires distributed acceleration of electrons up to near PeV energies along the jet. The underlying acceleration mechanism is still unclear. Shear acceleration is a promising candidate, as velocity-shear stratification is a natural consequence of the collimated flow of a jet. We study the details of shear acceleration by solving the steady-state Fokker-Planck-type equation and provide a simple general solution for trans-relativistic jets for a range of magnetohydrodynamic turbulent power-law spectra. In general, the accelerated particle population is a power-law spectrum with an exponential-like cut-off, where the power-law index is determined by the turbulence spectrum and the balance of escape and acceleration of particles. Adopting a simple linearly decreasing velocity profile in the boundary of large-scale jets, we find that the multi-wavelength spectral energy distribution of X-ray jets, such as Centaurus A and 3C 273, can be reproduced with electrons that are accelerated up to $\sim$ PeV. In kpc-scale jets, protons may be accelerated up to $\sim$ EeV, supporting the hypothesis that large-scale jets are strong candidates for ultra-high-energy-cosmic-ray sources within the framework of shear acceleration. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Comments: 9 pages, 3 figures, 1 table, accepted to MNRAS

Giant cosmic ray halos around M31 and the Milky Way

Authors: S. Recchia, S. Gabici, F. A. Aharonian, V. Niro

Abstract: Recently, a diffuse emission of 1-100 GeV $γ$-rays has been detected from the direction of Andromeda. The emission is centered on the galaxy, and extends for $\sim 100-200$ kpc away from its center. Explaining the extended $γ-$ray emission within the framework of standard scenarios for the escape of cosmic rays injected in the galactic disk or in the galactic center is problematic. In this paper,… ▽ More Recently, a diffuse emission of 1-100 GeV $γ$-rays has been detected from the direction of Andromeda. The emission is centered on the galaxy, and extends for $\sim 100-200$ kpc away from its center. Explaining the extended $γ-$ray emission within the framework of standard scenarios for the escape of cosmic rays injected in the galactic disk or in the galactic center is problematic. In this paper, we argue that a cosmic ray origin (either leptonic or hadronic) of the $γ$-ray emission is possible in the framework of non standard cosmic ray propagation scenarios or in the case of particle acceleration taking place in the galaxy's halo. It would imply the existence of a giant cosmic ray halo surrounding M31, possibly powered by the galaxy nuclear activity, or by accretion of intergalactic gas. Remarkably, if cosmic ray halos, as the one observed around M31, are a common feature of galaxies, including our own, the interactions between cosmic ray protons and the Milky Way circumgalactic gas could also explain the isotropic diffuse flux of neutrinos observed by Icecube. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: Submitted to ApJ, 13 pages, 3 figures

Time Variability of Nonthermal X-ray Stripes in Tycho's Supernova Remnant with Chandra

Authors: Tomoyuki Okuno, Takaaki Tanaka, Hiroyuki Uchida, Felix A. Aharonian, Yasunobu Uchiyama, Takeshi Go Tsuru, Masamune Matsuda

Abstract: Analyzing Chandra data of Tycho's supernova remnant (SNR) taken in 2000, 2003, 2007, 2009, and 2015, we search for time variable features of synchrotron X-rays in the southwestern part of the SNR, where stripe structures of hard X-ray emission were previous found. By comparing X-ray images obtained at each epoch, we discover a knot-like structure in the northernmost part of the stripe region becam… ▽ More Analyzing Chandra data of Tycho's supernova remnant (SNR) taken in 2000, 2003, 2007, 2009, and 2015, we search for time variable features of synchrotron X-rays in the southwestern part of the SNR, where stripe structures of hard X-ray emission were previous found. By comparing X-ray images obtained at each epoch, we discover a knot-like structure in the northernmost part of the stripe region became brighter particularly in 2015. We also find a bright filamentary structure gradually became fainter and narrower as it moved outward. Our spectral analysis reveal that not only the nonthermal X-ray flux but also the photon indices of the knot-like structure change from year to year. During the period from 2000 to 2015, the small knot shows brightening of $\sim 70\%$ and hardening of $ΔΓ\sim 0.45$. The time variability can be explained if the magnetic field is amplified to $\sim 100~\mathrm{μG}$ and/or if magnetic turbulence significantly changes with time. △ Less

Submitted 15 April, 2020; v1 submitted 22 March, 2020; originally announced March 2020.

Comments: 8 pages, 3 figures, 2 tables, accepted for publication in ApJ

Journal ref: 2020, ApJ, 894, 50

Magnetic absorption of VHE photons in the magnetosphere of the Crab pulsar

Authors: S. V. Bogovalov, I. Contopoulos, A. Prosekin, I. Tronin, F. A. Aharonian

Abstract: The detection of the pulsed $\sim 1 $~TeV gamma-ray emission from the Crab pulsar reported by MAGIC and VERITAS collaborations demands a substantial revision of existing models of particle acceleration in the pulsar magnetosphere. In this regard model independent restrictions on the possible production site of the VHE photons become an important issue. In this paper, we consider limitations impose… ▽ More The detection of the pulsed $\sim 1 $~TeV gamma-ray emission from the Crab pulsar reported by MAGIC and VERITAS collaborations demands a substantial revision of existing models of particle acceleration in the pulsar magnetosphere. In this regard model independent restrictions on the possible production site of the VHE photons become an important issue. In this paper, we consider limitations imposed by the process of conversion of VHE gamma rays into $e^{\pm}$ pairs in the magnetic field of the pulsar magnetosphere. Photons with energies exceeding 1~TeV are effectively absorbed even at large distances from the surface of the neutron star. Our calculations of magnetic absorption in the force-free magnetosphere show that the twisting of the magnetic field due to the pulsar rotation makes the magnetosphere more transparent compared to the dipole magnetosphere. The gamma-ray absorption appears stronger for photons emitted in the direction of rotation than in the opposite direction. There is a small angular cone inside which the magnetosphere is relatively transparent and photons with energy $1.5$~TeV can escape from distances beyond $0.1$~light cylinder radius ($R_{\rm{lc}}$). The emission surface from where photons can be emitted in the observer's direction further restricts the sites of VHE gamma-ray production. For the observation angle $57^{\circ}$ relative to the Crab pulsar axis of rotation and the orthogonal rotation, the emission surface in the open field line region is located as close as $0.4\,R_{\rm{lc}}$ from the stellar surface for a dipole magnetic field, and $0.1\,R_{\rm{lc}}$ for a force-free magnetic field. △ Less

Submitted 4 February, 2019; originally announced February 2019.

Comments: 22 pages, 9 figures

Journal ref: MNRAS, 2018

A local fading accelerator and the origin of TeV cosmic ray electrons

Authors: S. Recchia, S. Gabici, F. A. Aharonian, J. Vink

Abstract: The total cosmic ray electron spectrum (electrons plus positrons) exhibits a break at a particle energy of $\sim 1\rm~TeV$ and extends without any attenuation up to $\rm \sim 20~ TeV $. Synchrotron and inverse Compton energy losses strongly constrain both the age and the distance of the potential sources of TeV and multi-TeV electrons to $\rm\approx 10^5~yr$ and $\rm \approx 100-500~pc$, depending… ▽ More The total cosmic ray electron spectrum (electrons plus positrons) exhibits a break at a particle energy of $\sim 1\rm~TeV$ and extends without any attenuation up to $\rm \sim 20~ TeV $. Synchrotron and inverse Compton energy losses strongly constrain both the age and the distance of the potential sources of TeV and multi-TeV electrons to $\rm\approx 10^5~yr$ and $\rm \approx 100-500~pc$, depending on both the absolute value and energy dependence of the cosmic ray diffusion coefficient. This suggests that only a few, or just one nearby discrete source may explain the observed spectrum of high energy electrons. On the other hand the measured positron fraction, after initially increasing with particle energy, saturates at a level well below 0.5 and likely drops above $\sim 400-500$ GeV. This means that the local source(s) of TeV electrons should not produce positrons in equal amount, ruling out scenarios involving pulsars/pulsar winds as the main sources of high energy leptons. In this paper we show that a single, local, and fading source can naturally account for the entire spectrum of cosmic ray electrons in the TeV domain. Even though the nature of such source remains unclear, we discuss known cosmic ray accelerators, such as supernova remnant and stellar wind shocks, which are believed to accelerate preferentially electrons rather than positrons. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Journal ref: Phys. Rev. D 99, 103022 (2019)

NuSTAR Detection of Nonthermal Bremsstrahlung from the Supernova Remnant W49B

Authors: Takaaki Tanaka, Hiroya Yamaguchi, Daniel R. Wik, Felix A. Aharonian, Aya Bamba, Daniel Castro, Adam R. Foster, Robert Petre, Jeonghee Rho, Randall K. Smith, Hiroyuki Uchida, Yasunobu Uchiyama, Brian J. Williams

Abstract: We report on NuSTAR observations of the mixed morphology supernova remnant (SNR) W49B, focusing on its nonthermal emission. Whereas radio observations as well as recent gamma-ray observations evidenced particle acceleration in this SNR, nonthermal X-ray emission has not been reported so far. With the unprecedented sensitivity of NuSTAR in the hard X-ray band, we detect a significant power-law-like… ▽ More We report on NuSTAR observations of the mixed morphology supernova remnant (SNR) W49B, focusing on its nonthermal emission. Whereas radio observations as well as recent gamma-ray observations evidenced particle acceleration in this SNR, nonthermal X-ray emission has not been reported so far. With the unprecedented sensitivity of NuSTAR in the hard X-ray band, we detect a significant power-law-like component extending up to $\sim 20~{\rm keV}$, most probably of nonthermal origin. The newly discovered component has a photon index of $Γ=1.4^{+1.0}_{-1.1}$ with an energy flux between 10 and 20 keV of $(3.3 \pm 0.7) \times 10^{-13}~{\rm erg}~{\rm cm}^{-2}~{\rm s}^{-1}$. The emission mechanism is discussed based on the NuSTAR data combined with those in other wavelengths in the literature. The NuSTAR data, in terms both of the spectral slope and of the flux, are best interpreted as nonthermal electron bremsstrahlung. If this scenario is the case, then the NuSTAR emission provides a new probe to sub-relativistic particles accelerated in the SNR. △ Less

Submitted 16 October, 2018; originally announced October 2018.

Comments: 7 pages, 3 figures, 2 tables, accepted for publication in ApJ Letters

Journal ref: 2018, ApJL, 866, L26

Particle Acceleration in Mildly Relativistic Shearing Flows: the Interplay of Systematic and Stochastic Effects, and the Origin of the Extended High-energy Emission in AGN Jets

Authors: Ruo-Yu Liu, F. M. Rieger, F. A. Aharonian

Abstract: The origin of the extended X-ray emission in the large-scale jets of active galactic nuclei (AGNs) poses challenges to conventional models of acceleration and emission. Although the electron synchrotron radiation is considered the most feasible radiation mechanism, the formation of the continuous large-scale X-ray structure remains an open issue. As astrophysical jets are expected to exhibit some… ▽ More The origin of the extended X-ray emission in the large-scale jets of active galactic nuclei (AGNs) poses challenges to conventional models of acceleration and emission. Although the electron synchrotron radiation is considered the most feasible radiation mechanism, the formation of the continuous large-scale X-ray structure remains an open issue. As astrophysical jets are expected to exhibit some turbulence and shearing motion, we here investigate the potential of shearing flows to facilitate an extended acceleration of particles and evaluate its impact on the resultant particle distribution. Our treatment incorporates systematic shear and stochastic second-order Fermi effects. We show that for typical parameters applicable to large-scale AGN jets, stochastic second-order Fermi acceleration, which always accompanies shear particle acceleration, can play an important role in facilitating the whole process of particle energization. We study the time-dependent evolution of the resultant particle distribution in the presence of second-order Fermi acceleration, shear acceleration, and synchrotron losses using a simple Fokker--Planck approach and provide illustrations for the possible emergence of a complex (multicomponent) particle energy distribution with different spectral branches. We present examples for typical parameters applicable to large-scale AGN jets, indicating the relevance of the underlying processes for understanding the extended X-ray emission and the origin of ultrahigh-energy cosmic rays. △ Less

Submitted 4 June, 2017; originally announced June 2017.

Comments: 26 pages, 8 figures; to appear in ApJ

Scenarios for ultrafast gamma-ray variability in AGN

Authors: F. A. Aharonian, M. V. Barkov, D. Khangulyan

Abstract: We analyze three scenarios to address the challenge of ultrafast gamma-ray variability reported from active galactic nuclei. We focus on the energy requirements imposed by these scenarios: (i) external cloud in the jet, (ii) relativistic blob propagating through the jet material, and (iii) production of high-energy gamma rays in the magnetosphere gaps. We show that while the first two scenarios ar… ▽ More We analyze three scenarios to address the challenge of ultrafast gamma-ray variability reported from active galactic nuclei. We focus on the energy requirements imposed by these scenarios: (i) external cloud in the jet, (ii) relativistic blob propagating through the jet material, and (iii) production of high-energy gamma rays in the magnetosphere gaps. We show that while the first two scenarios are not constrained by the flare luminosity, there is a robust upper limit on the luminosity of flares generated in the black hole magnetosphere. This limit depends weakly on the mass of the central black hole and is determined by the accretion disk magnetization, viewing angle, and the pair multiplicity. For the most favorable values of these parameters, the luminosity for 5-minute flares is limited by $2\times10^{43}\rm\,erg\,s^{-1}$, which excludes a black hole magnetosphere origin of the flare detected from IC310. In the scopes of scenarios (i) and (ii), the jet power, which is required to explain the IC310 flare, exceeds the jet power estimated based on the radio data. To resolve this discrepancy in the framework of the scenario (ii), it is sufficient to assume that the relativistic blobs are not distributed isotropically in the jet reference frame. A realization of scenario (i) demands that the jet power during the flare exceeds by a factor $10^2$ the power of the radio jet relevant to a timescale of $10^8$ years. △ Less

Submitted 16 May, 2017; v1 submitted 26 April, 2017; originally announced April 2017.

Comments: 15 pages, accepted by ApJ

Interstellar gas towards the TeV gamma-ray sources HESS J1640-465 and HESS J1641-463

Authors: J. C. Lau, G. Rowell, M. G. Burton, Y. Fukui, F. A. Aharonian, I. Oya, J. Vink, S. Ohm, S. Casanova

Abstract: We present a detailed analysis of the interstellar medium towards the TeV $γ$-ray sources HESS J1640$-$465 and HESS J1641$-$463 using results from the Mopra Southern Galactic Plane CO Survey and from a Mopra 7 mm-wavelength study. The $γ$-ray sources are positionally coincident with two supernova remnants G338.3$-$0.0 and G338.5+0.1 respectively. A bright complex of HII regions connect the two SNR… ▽ More We present a detailed analysis of the interstellar medium towards the TeV $γ$-ray sources HESS J1640$-$465 and HESS J1641$-$463 using results from the Mopra Southern Galactic Plane CO Survey and from a Mopra 7 mm-wavelength study. The $γ$-ray sources are positionally coincident with two supernova remnants G338.3$-$0.0 and G338.5+0.1 respectively. A bright complex of HII regions connect the two SNRs and TeV objects. Observations in the CO(1-0) transition lines reveal substantial amounts of diffuse gas positionally coincident with the $γ$-ray sources at multiple velocities along the line of sight, while 7 mm observations in CS, SiO, HC$_{3}$N and CH$_{3}$OH transition lines reveal regions of dense, shocked gas. Archival HI data from the Southern Galactic Plane Survey was used to account for the diffuse atomic gas. Physical parameters of the gas towards the TeV sources were calculated from the data. We find that for a hadronic origin for the $γ$-ray emission, the cosmic-ray enhancement rates are $\sim 10^{3}$ and $10^{2}$ times the local solar value for HESS J1640$-$465 and HESS J1641$-$463 respectively. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Comments: 20 pages, 13 figures, 6 tables, accepted for publication by MNRAS

Hitomi constraints on the 3.5 keV line in the Perseus galaxy cluster

Authors: Hitomi Collaboration, Felix A. Aharonian, Hiroki Akamatsu, Fumie Akimoto, Steven W. Allen, Lorella Angelini, Keith A. Arnaud, Marc Audard, Hisamitsu Awaki, Magnus Axelsson, Aya Bamba, Marshall W. Bautz, Roger D. Blandford, Laura W. Brenneman, Gregory V. Brown, Esra Bulbul, Edward M. Cackett, Maria Chernyakova, Meng P. Chiao, Paolo Coppi, Elisa Costantini, Jelle de Plaa, Jan-Willem den Herder, Chris Done, Tadayasu Dotani , et al. (193 additional authors not shown)

Abstract: High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-New… ▽ More High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment. △ Less

Submitted 27 February, 2017; v1 submitted 25 July, 2016; originally announced July 2016.

Comments: Discussion of systematics significantly expanded. 9 pages, 5 figures; ApJ Lett. in press

Journal ref: ApJ, 837, L15 (2017)

Polarization of radiation of electrons in highly turbulent magnetic fields

Authors: A. Yu. Prosekin, S. R. Kelner, F. A. Aharonian

Abstract: We study the polarization properties of the jitter and synchrotron radiation produced by electrons in highly turbulent anisotropic magnetic fields. The net polarization is provided by the geometry of the magnetic field the directions of which are parallel to a certain plane. Such conditions may appear in the relativistic shocks during the amplification of the magnetic field through the so-called W… ▽ More We study the polarization properties of the jitter and synchrotron radiation produced by electrons in highly turbulent anisotropic magnetic fields. The net polarization is provided by the geometry of the magnetic field the directions of which are parallel to a certain plane. Such conditions may appear in the relativistic shocks during the amplification of the magnetic field through the so-called Weibel instability. While the polarization properties of the jitter radiation allows extraction of direct information on the turbulence spectrum as well as the geometry of magnetic field, the polarization of the synchrotron radiation reflects the distribution of the magnetic field over its strength. For the isotropic distribution of monoenergetic electrons, we found that the degree of polarization of the synchrotron radiation is larger than the polarization of the jitter radiation. For the power-law energy distribution of electrons the relation between the degree of polarization of synchrotron and jitter radiation depends on the spectral index of the distribution. △ Less

Submitted 6 July, 2016; originally announced July 2016.

Journal ref: Phys. Rev. D 94, 063010 (2016)

Indication of a Local "Fog" of Sub-Ankle UHECR

Authors: Ruo-Yu Liu, Andrew M. Taylor, Xiang-Yu Wang, Felix A. Aharonian

Abstract: During their propagation through intergalactic space, ultrahigh energy cosmic rays (UHECRs) interact with the background radiation fields. These interactions give rise to energetic electron/positron pairs and photons which in turn feed electromagnetic cascades, contributing to the isotropic gamma-ray background (IGRB). The gamma-ray flux level generated in this way highly depends upon the UHECR pr… ▽ More During their propagation through intergalactic space, ultrahigh energy cosmic rays (UHECRs) interact with the background radiation fields. These interactions give rise to energetic electron/positron pairs and photons which in turn feed electromagnetic cascades, contributing to the isotropic gamma-ray background (IGRB). The gamma-ray flux level generated in this way highly depends upon the UHECR propagation distance, as well as the evolution of their sources with redshift. Recently, the Fermi-LAT collaboration reported that the majority of the total extragalactic gamma-ray flux originates from extragalactic point sources. This posits a stringent upper limit on the IGRB generated via UHECR propagation, and subsequently constrains their abundance in the distant Universe. Focusing on the contribution of UHECR at energies below the ankle within a narrow energy band ($(1-4)\times 10^{18}$eV), we calculate the diffuse gamma-ray flux generated through UHECR propagation, normalizing the total cosmic ray energy budget in this band to that measured. We find that in order to not over-produce the new IGRB limit, a local "fog" of UHECR produced by nearby sources may exist, with a possible non-negligible contribution from our Galaxy. Following the assumption that a given fraction of the observed IGRB at 820GeV originates from UHECR, we obtain a constraint on the maximum distance for the majority of their sources. With other unresolved source populations still contaminating the new IGRB limit, and UHECR above the ankle invariably contributing also to this background, the results presented here are rather conservative. △ Less

Submitted 4 September, 2016; v1 submitted 10 March, 2016; originally announced March 2016.

Comments: 6 pages, 2 figures, comments are welcome; accepted by PRD

Journal ref: Phys. Rev. D 94, 043008 (2016)

On transition of propagation of relativistic particles from the ballistic to the diffusion regime

Authors: A. Y. Prosekin, S. R. Kelner, F. A. Aharonian

Abstract: A stationary distribution function that describes the entire processes of propagation of relativistic particles, including the transition between the ballistic and diffusion regimes, is obtained. The spacial component of the constructed function satisfies to the first two moments of the Boltzmann equation. The angular part of the distribution provides accurate values for the angular moments derive… ▽ More A stationary distribution function that describes the entire processes of propagation of relativistic particles, including the transition between the ballistic and diffusion regimes, is obtained. The spacial component of the constructed function satisfies to the first two moments of the Boltzmann equation. The angular part of the distribution provides accurate values for the angular moments derived from the Boltzmann equation, and gives a correct expression in the limit of small-angle approximation. Using the derived function, we studied the gamma-ray images produced through the $pp$ interaction of relativistic particles with gas clouds in the proximity of the accelerator. In general, the morphology and the energy spectra of gamma-rays significantly deviate from the "standard" results corresponding to the propagation of relativistic particles strictly in the diffusion regime. △ Less

Submitted 22 June, 2015; originally announced June 2015.

Journal ref: Phys. Rev. D 92, 083003 (2015)

Star-Jet Interactions and Gamma-Ray Outbursts from 3c454.3

Authors: M. V. Barkov, D. V. Khangulyan, V. Bosch-Ramon, F. A. Aharonian, A. V. Dorodnitsyn

Abstract: We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar \c454. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calcula… ▽ More We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar \c454. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calculations of the related high-energy radiation. The model readily explains the day-long, variable on timescales of hours, GeV gamma-ray flare from \c454, observed during November 2010 on top of a weeks-long plateau. In the proposed scenario, the plateau state is caused by a strong wind generated by the heating of the star atmosphere by nonthermal particles accelerated at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key model parameters of the source, including the mass of the central black hole: $M_{\rm BH}\simeq 10^9 M_{\odot}$, the total jet power: $L_{\rm j}\simeq 10^{48}\,\rm erg\,s^{-1}$, and the Doppler factor of the gamma-ray emitting clouds, $δ\simeq 20$. Whereas we do not specify the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in the context of the proposed model. We argue that synchrotron radiation of protons has certain advantages compared to other radiation channels of directly accelerated electrons. △ Less

Submitted 9 February, 2015; originally announced February 2015.

Comments: 2014 Fermi Symposium proceedings - eConf C14102.1, 7 pages, 3 figures

Gamma-ray emission from young supernova remnants: hadronic or leptonic?

Authors: S. Gabici, F. A. Aharonian

Abstract: The debate on the nature of the gamma-ray emission from young supernova remnants is still open. Ascribing such emission to hadronic rather than leptonic processes would provide an evidence for the acceleration of protons and nuclei, and this fact would fit with the very popular (but not proven) paradigm that supernova remnants are the sources of Galactic cosmic rays. Here, we discuss this issue wi… ▽ More The debate on the nature of the gamma-ray emission from young supernova remnants is still open. Ascribing such emission to hadronic rather than leptonic processes would provide an evidence for the acceleration of protons and nuclei, and this fact would fit with the very popular (but not proven) paradigm that supernova remnants are the sources of Galactic cosmic rays. Here, we discuss this issue with a particular focus on the best studied gamma-ray-bright supernova remnant: RX~J1713.7-3946. △ Less

Submitted 2 February, 2015; originally announced February 2015.

Comments: invited talk at RICAP 2014

Synchro-curvature radiation of charged particles in the strong curved magnetic fields

Authors: S. R. Kelner, A. Yu. Prosekin, F. A. Aharonian

Abstract: It is generally believed that the radiation of relativistic particles in a curved magnetic field proceeds in either the synchrotron or the curvature radiation modes. In this paper we show that in strong curved magnetic fields a significant fraction of the energy of relativistic electrons can be radiated away in the intermediate, the so-called synchro-curvature regime. Because of the persistent cha… ▽ More It is generally believed that the radiation of relativistic particles in a curved magnetic field proceeds in either the synchrotron or the curvature radiation modes. In this paper we show that in strong curved magnetic fields a significant fraction of the energy of relativistic electrons can be radiated away in the intermediate, the so-called synchro-curvature regime. Because of the persistent change of the trajectory curvature, the radiation varies with the frequency of particle gyration. While this effect can be ignored in the synchrotron and curvature regimes, the variability plays a key role in the formation of the synchro-curvature radiation. Using the Hamiltonian formalism, we find that the particle trajectory has the form of a helix wound around the drift trajectory. This allows us to calculate analytically the intensity and energy distribution of prompt radiation in the general case of magnetic bremsstrahlung in the curved magnetic field. We show that the transition to the limit of the synchrotron and curvature radiation regimes is determined by the relation between the drift velocity and the component of the particle velocity perpendicular to the drift trajectory. The detailed numerical calculations, which take into account the energy losses of particles, confirm the principal conclusions based on the simplified analytical treatment of the problem, and allow us to analyze quantitatively the transition between different radiation regimes for a broad range of initial pitch angles. We argue that in the case of realization of specific configurations of the electric and magnetic fields, the gamma-ray emission of the pulsar magnetospheres can be dominated by the component radiated in the synchro-curvature regime. △ Less

Submitted 19 January, 2015; originally announced January 2015.

Comments: this article supersedes arXiv:1207.6903 and arXiv:1305.0783

Journal ref: The Astronomical Journal 149, 33 (2015)

Hadronic gamma-rays from RX J1713.7-3946?

Authors: Stefano Gabici, Felix A. Aharonian

Abstract: RXJ1713.7-3946 is a key object to check the supernova remnant paradigm of the origin of Galactic cosmic rays. While the origin of its gamma-ray emission (hadronic versus leptonic) is still debated, the hard spectrum at GeV energies reported by the Fermi collaboration is generally interpreted as a strong argument in favor of a leptonic scenario. On the contrary, we show that hadronic interactions c… ▽ More RXJ1713.7-3946 is a key object to check the supernova remnant paradigm of the origin of Galactic cosmic rays. While the origin of its gamma-ray emission (hadronic versus leptonic) is still debated, the hard spectrum at GeV energies reported by the Fermi collaboration is generally interpreted as a strong argument in favor of a leptonic scenario. On the contrary, we show that hadronic interactions can naturally explain the gamma-ray spectrum if gas clumps are present in the supernova remnant shell. The absence of thermal X-rays from the remnant fits well within this scenario. △ Less

Submitted 9 June, 2014; originally announced June 2014.

Comments: 4 pages, 2 figures. submitted

Simple analytical approximations for treatment of inverse Compton scattering of relativistic electrons in the black-body radiation field

Authors: Dmitry Khangulyan, Felix A. Aharonian, Stanislav R. Kelner

Abstract: The inverse Compton (IC) scattering of relativistic electrons is one of the major gamma-ray production mechanisms in different environments. Often the target photons for the IC scattering are dominated by black (or grey) body radiation. In this case, the precise treatment of the characteristics of IC radiation requires numerical integrations over the Planckian distribution. Formally, analytical in… ▽ More The inverse Compton (IC) scattering of relativistic electrons is one of the major gamma-ray production mechanisms in different environments. Often the target photons for the IC scattering are dominated by black (or grey) body radiation. In this case, the precise treatment of the characteristics of IC radiation requires numerical integrations over the Planckian distribution. Formally, analytical integrations are also possible but they result in series of several special functions; this limits the efficiency of usage of these expressions. The aim of this work is the derivation of approximate analytical presentations which would provide adequate accuracy for the calculations of the energy spectra of up-scattered radiation, the rate of electron energy losses, and the mean energy of emitted photons. Such formulae have been obtained by merging the analytical asymptotic limits. The coefficients in these expressions are calculated via the least square fitting of the results of numerical integrations. The simple analytical presentations, obtained for both the isotropic and anisotropic target radiation fields, provide adequate (as good as $1\%$) accuracy for broad astrophysical applications. △ Less

Submitted 24 January, 2014; v1 submitted 29 October, 2013; originally announced October 2013.

Comments: 16 pages, 11 figures, accepted for publication in ApJ

The beaming pattern of External Compton Emission from relativistic outflows

Authors: S. R. Kelner, E. Lefa, F. M. Rieger, F. A. Aharonian

Abstract: The beaming pattern of radiation emitted by a relativistically moving source like jets in microquasars, AGN and GRBs, is a key issue for understanding of acceleration and radiation processes in these objects. In this paper we introduce a formalism based on a solution of the photon transfer equation to study the beaming patterns for emission produced by electrons accelerated in the jet and upscatte… ▽ More The beaming pattern of radiation emitted by a relativistically moving source like jets in microquasars, AGN and GRBs, is a key issue for understanding of acceleration and radiation processes in these objects. In this paper we introduce a formalism based on a solution of the photon transfer equation to study the beaming patterns for emission produced by electrons accelerated in the jet and upscattering photons of low-energy radiation fields of external origin (the so-called External Compton scenario). The formalism allows us to treat non-stationary, non-homogeneous and anisotropic distributions of electrons, but assuming homogeneous/isotropic and non-variable target photon fields. We demonstrate the non-negligible impact of the anisotropy in the electron distribution on angular and spectral characteristics of the EC radiation. △ Less

Submitted 23 August, 2013; originally announced August 2013.

Comments: 7 pages, 2 figures

Nonthermal radiation of young supernova remnants: the case of Cas A

Authors: V. N. Zirakashvili, F. A. Aharonian, R. Yang, E. Ona-Wilhelmi, R. J. Tuffs

Abstract: The processes responsible for the broad-band radiation of the young supernova remnant Cas A are explored using a new code which is designed for a detailed treatment of the diffusive shock acceleration of particles in nonlinear regime. The model is based on spherically symmetric hydrodynamic equations complemented with transport equations for relativistic particles. Electrons, protons and the oxyge… ▽ More The processes responsible for the broad-band radiation of the young supernova remnant Cas A are explored using a new code which is designed for a detailed treatment of the diffusive shock acceleration of particles in nonlinear regime. The model is based on spherically symmetric hydrodynamic equations complemented with transport equations for relativistic particles. Electrons, protons and the oxygen ions accelerated by forward and reverse shocks are included in the numerical calculations. We show that the available multi-wavelength observations in the radio, X-ray and gamma-ray bands can be best explained by invoking particle acceleration by both forward and reversed shocks. Although the TeV gamma-ray observations can be interpreted by interactions of both accelerated electrons and protons/ions, the measurements by Fermi LAT at energies below 1 GeV give a tentative preference to the hadronic origin of gamma-rays. Then, the acceleration efficiency in this source, despite the previous claims, should be very high; 25% of the explosion energy (or approximately $3\cdot 10^{50}$ erg) should already be converted to cosmic rays, mainly by the forward shock. At the same time, the model calculations do not provide extension of the maximum energy of accelerated protons beyond 100 TeV. In this model, the acceleration of electrons is dominated by the reverse shock; the required $10^{48}$ erg can be achieved under the assumption that the injection of electrons (positrons) is supported by the radioactive decay of $^{44}$Ti. △ Less

Submitted 16 August, 2013; originally announced August 2013.

Comments: 28 pages, 14 figures, 3 tables, submitted to ApJ

Star-Jet Interactions and Gamma-Ray Outbursts from 3C454.3

Authors: D. V. Khangulyan, M. V. Barkov, V. Bosch-Ramon, F. A. Aharonian, A. V. Dorodnitsyn

Abstract: We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar 3C454.3. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calcu… ▽ More We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar 3C454.3. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calculations of the related high-energy radiation. The model readily explains the day-long, variable on timescales of hours, GeV gamma-ray flare from 3C454.3, observed during November 2010 on top of a weeks-long plateau. In the proposed scenario, the plateau state is caused by a strong wind generated by the heating of the star atmosphere by nonthermal particles accelerated at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key model parameters of the source, including the mass of the central black hole: $M_{\rm BH}\simeq 10^9 M_{\odot}$, the total jet power: $L_{\rm j}\simeq 10^{48}\,\rm erg\,s^{-1}$, and the Doppler factor of the gamma-ray emitting clouds, $δ\simeq 20$. Whereas we do not specify the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in the context of the proposed model. We argue that synchrotron radiation of protons has certain advantages compared to other radiation channels of directly accelerated electrons. △ Less

Submitted 27 May, 2013; v1 submitted 22 May, 2013; originally announced May 2013.

Comments: 16 pages, 5 figures, submitted to ApJ

Synchrotron-to-curvature transition regime of radiation of charged particles in a dipole magnetic field

Authors: A. Yu. Prosekin, S. R. Kelner, F. A. Aharonian

Abstract: The details of trajectories of charged particles become increasingly important for proper understanding of processes of formation of radiation in strong and curved magnetic fields. Because of damping of the perpendicular component of motion, the particle's pitch angle could be decreased by many orders of magnitude leading to the change of the radiation regime -- from synchrotron to the curvature m… ▽ More The details of trajectories of charged particles become increasingly important for proper understanding of processes of formation of radiation in strong and curved magnetic fields. Because of damping of the perpendicular component of motion, the particle's pitch angle could be decreased by many orders of magnitude leading to the change of the radiation regime -- from synchrotron to the curvature mode. To explore the character of this transition, we solve numerically the equations of motion of a test particle in a dipole magnetic field, and calculate the energy spectrum of magnetic bremsstrahlung self-consistently, i.e. without a priori assumptions on the radiation regime. In this way we can trace the transitions between the synchrotron and curvature regimes, as well as study the third (intermediate or the so-called synchro-curvature) regime. We briefly discuss three interesting astrophysical scenarios, the radiation of electrons in the pulsar magnetosphere in the polar cap and outer gap models, as well as the radiation of ultrahigh energy protons in the magnetosphere of a massive black hole, and demonstrate that in these models the synchrotron, synchro-curvature and curvature regimes can be realized with quite different relative contributions to the total emission. △ Less

Submitted 22 January, 2015; v1 submitted 3 May, 2013; originally announced May 2013.

Comments: the reprocessed version of this article merged with arXiv:1207.6903 can be found in arXiv:1501.04994

On the Jitter Radiation

Authors: Stanislav R. Kelner, Felix A. Aharonian, Dmitry Khangulyan

Abstract: In a small scale turbulent medium, when the nonrelativistic Larmor radius $R_{\rm L}=mc^2/eB$ exceeds the correlation length $λ$ of the magnetic field, the magnetic bremsstrahlung of charged relativistic particles unavoidably proceeds in the so-called jitter radiation regime. The cooling timescale of parent particles is identical to the synchrotron cooling time, thus this radiation regime can be p… ▽ More In a small scale turbulent medium, when the nonrelativistic Larmor radius $R_{\rm L}=mc^2/eB$ exceeds the correlation length $λ$ of the magnetic field, the magnetic bremsstrahlung of charged relativistic particles unavoidably proceeds in the so-called jitter radiation regime. The cooling timescale of parent particles is identical to the synchrotron cooling time, thus this radiation regime can be produced with very high efficiency in different astrophysical sources characterized by high turbulence. The jitter radiation has distinct spectral features shifted, compared to synchrotron radiation, towards high energies. This makes the jitter mechanism an attractive broad-band gamma-ray production channel which in highly magnetized and turbulent environments can compete or even dominate over other high energy radiation mechanisms. In this paper we present a novel study on spectral properties of the jitter radiation performed within the framework of perturbation theory. The derived general expression for the spectral power of radiation is presented in a compact and convenient for numerical calculations form. △ Less

Submitted 30 May, 2013; v1 submitted 1 April, 2013; originally announced April 2013.

Comments: 17 pages, 6 figures, submitted to ApJ

Evidence for a second component in the high-energy core emission from Centaurus A?

Authors: N. Sahakyan, R. Yang, F. A. Aharonian, F. M. Rieger

Abstract: We report on an analysis of \fermi data from four year of observations of the nearby radio galaxy Centaurus A (Cen A). The increased photon statistics results in a detection of high-energy ($>\:100$ MeV) $γ$-rays up to 50 GeV from the core of Cen A, with a detection significance of about 44$σ$. The average gamma-ray spectrum of the core reveals evidence for a possible deviation from a simple power… ▽ More We report on an analysis of \fermi data from four year of observations of the nearby radio galaxy Centaurus A (Cen A). The increased photon statistics results in a detection of high-energy ($>\:100$ MeV) $γ$-rays up to 50 GeV from the core of Cen A, with a detection significance of about 44$σ$. The average gamma-ray spectrum of the core reveals evidence for a possible deviation from a simple power-law. A likelihood analysis with a broken power-law model shows that the photon index becomes harder above $E_b \simeq 4$ GeV, changing from $Γ_1=2.74\pm0.03$ below to $Γ_2=2.09\pm0.20$ above. This hardening could be caused by the contribution of an additional high-energy component beyond the common synchrotron-self Compton jet emission. A variability analysis of the light curve with 15-, 30-, and 60-day bins does not provide evidence for variability for any of the components. Indications for a possible variability of the observed flux are found on 45-day time scale, but the statistics do not allow us to make a definite conclusion in this regards. We compare our results with the spectrum reported by H.E.S.S. in the TeV energy range and discuss possible origins for the hardening observed. △ Less

Submitted 28 February, 2013; originally announced February 2013.

Comments: submitted to ApJ Letters

TeV Astronomy

Authors: Frank M. Rieger, Emma de Ona-Wilhelmi, Felix A. Aharonian

Abstract: With the successful realization of the current-generation of ground-based detectors, TeV Astronomy has entered into a new era. We review recent advances in VHE astronomy, focusing on the potential of Imaging Atmospheric Cherenkov Telescopes (IACTs), and highlight astrophysical implications of the results obtained within recent years. With the successful realization of the current-generation of ground-based detectors, TeV Astronomy has entered into a new era. We review recent advances in VHE astronomy, focusing on the potential of Imaging Atmospheric Cherenkov Telescopes (IACTs), and highlight astrophysical implications of the results obtained within recent years. △ Less

Submitted 22 February, 2013; originally announced February 2013.

Comments: Invited review article to appear in the special issue of "Frontiers of Physics" on "High Energy Astrophysics" (eds. B. Zhang & P. Meszaros); review, 33 pages, 23 figures, 305 references

Analytic Solution for Self-regulated Collective Escape of Cosmic Rays from their Acceleration Sites

Authors: M. A. Malkov, P. H. Diamond, R. Z. Sagdeev, F. A. Aharonian, I. V. Moskalenko

Abstract: Supernova remnants (SNRs), as the major contributors to the galactic cosmic rays (CR), are believed to maintain an average CR spectrum by diffusive shock acceleration (DSA) regardless of the way they release CRs into the interstellar medium (ISM). However, the interaction of the CRs with nearby gas clouds crucially depends on the release mechanism. We call into question two aspects of a popular pa… ▽ More Supernova remnants (SNRs), as the major contributors to the galactic cosmic rays (CR), are believed to maintain an average CR spectrum by diffusive shock acceleration (DSA) regardless of the way they release CRs into the interstellar medium (ISM). However, the interaction of the CRs with nearby gas clouds crucially depends on the release mechanism. We call into question two aspects of a popular paradigm of the CR injection into the ISM, according to which they passively and isotropically diffuse in the prescribed magnetic fluctuations as test particles. First, we treat the escaping CR and the Alfven waves excited by them on an equal footing. Second, we adopt field aligned CR escape outside the source, where the waves become weak. An exact analytic self-similar solution for a CR "cloud" released by a dimmed accelerator strongly deviates from the test-particle result. The CR diffusion coefficient $D_{NL}$ is strongly suppressed compared to its background ISM value $D_{ISM}$: $D_{NL}\sim D_{ISM}\exp(-Π)<< D_{ISM}$ for sufficiently high field-line-integrated CR partial pressure, $Π$. When $Π>>1$, the CRs drive Alfven waves efficiently enough to build a transport barrier that strongly reduces the leakage. The solution has a spectral break at $p=p_{br}$, where $p_{br}$ satisfies the following equation $D_{NL}(p_{br})\simeq z^{2}/t$. △ Less

Submitted 20 March, 2013; v1 submitted 19 July, 2012; originally announced July 2012.

Comments: 6 figures, aastex, 31 pages

Journal ref: The Astrophysical Journal, Volume 768, Issue 1, article id. 73, (2013)

On the spectral shape of radiation due to Inverse Compton Scattering close to the maximum cut-off

Authors: E. Lefa, S. R. Kelner, F. A. Aharonian

Abstract: The spectral shape of radiation due to Inverse Compton Scattering is analyzed, in the Thomson and the Klein-Nishina regime, for electron distributions with exponential cut-off. We derive analytical, asymptotic expressions for the spectrum close to the maximum cut-off region. We consider monoenergetic, Planckian and Synchrotron photons as target photon fields. These approximations provide a direct… ▽ More The spectral shape of radiation due to Inverse Compton Scattering is analyzed, in the Thomson and the Klein-Nishina regime, for electron distributions with exponential cut-off. We derive analytical, asymptotic expressions for the spectrum close to the maximum cut-off region. We consider monoenergetic, Planckian and Synchrotron photons as target photon fields. These approximations provide a direct link between the distribution of parent electrons and the up-scattered spectrum at the cut-off region. △ Less

Submitted 13 May, 2012; originally announced May 2012.

Comments: 27 pages, 15 figures, accepted for publication in ApJ

Interpretation of the flares of M87 at TeV energies in the cloud-jet interaction scenario

Authors: Maxim V. Barkov, Valenti Bosch-Ramon, Felix A. Aharonian

Abstract: Active galactic nuclei with misaligned jets have been recently established as a class of high-energy gamma-ray sources. M87, a nearby representative of this class, shows fast TeV variability on timescales less than one day. We present calculations performed in the framework of the scenario in which gamma-ray flares in non-blazar active galactic nuclei are produced by a red giant or a gas cloud int… ▽ More Active galactic nuclei with misaligned jets have been recently established as a class of high-energy gamma-ray sources. M87, a nearby representative of this class, shows fast TeV variability on timescales less than one day. We present calculations performed in the framework of the scenario in which gamma-ray flares in non-blazar active galactic nuclei are produced by a red giant or a gas cloud interacting with the jet. We show that both the light curve and energy spectrum of the spectacular April 2010 flare can be reproduced by this model, assuming that a relatively massive cloud of approx 1.e29 g penetrates into the jet at few tens of Schwarzschild radii from the super-massive black hole. △ Less

Submitted 27 June, 2012; v1 submitted 27 February, 2012; originally announced February 2012.

Comments: 8 pages, 8 figures, accepted by ApJ

Magnetic substructure in the northern Fermi Bubble revealed by polarized WMAP emission

Authors: David I. Jones, Roland M. Crocker, Wolfgang Reich, Jürgen Ott, Felix A. Aharonian

Abstract: We report a correspondence between giant, polarized microwave structures emerging north from the Galactic plane near the Galactic center and a number of GeV gamma-ray features, including the eastern edge of the recently-discovered northern Fermi Bubble. The polarized microwave features also correspond to structures seen in the all-sky 408 MHz total intensity data, including the Galactic center spu… ▽ More We report a correspondence between giant, polarized microwave structures emerging north from the Galactic plane near the Galactic center and a number of GeV gamma-ray features, including the eastern edge of the recently-discovered northern Fermi Bubble. The polarized microwave features also correspond to structures seen in the all-sky 408 MHz total intensity data, including the Galactic center spur. The magnetic field structure revealed by the polarization data at 23 GHz suggests that neither the emission coincident with the Bubble edge nor the Galactic center spur are likely to be features of the local ISM. On the basis of the observed morphological correspondences, similar inferred spectra, and the similar energetics of all sources, we suggest a direct connection between the Galactic center spur and the northern Fermi Bubble. △ Less

Submitted 2 February, 2012; v1 submitted 21 January, 2012; originally announced January 2012.

Comments: Accepted for publication in The Astrophysical Journal Letters after minor changes

Journal ref: 2012 ApJ 747 L12

"Leading blob" model in a stochastic acceleration scenario: the case of the 2009 flare of Mkn~501

Authors: E. Lefa, F. A. Aharonian, F. M. Rieger

Abstract: Evidence for very hard, intrinsic $γ$-ray source spectra, as inferred after correction for absorption in the extragalactic background light (EBL), has interesting implications for the acceleration and radiation mechanisms acting in blazars. A key issue so far has been the dependance of the hardness of the $γ$-ray spectrum on different existing EBL models. The recent {\it Fermi} observations of Mkn… ▽ More Evidence for very hard, intrinsic $γ$-ray source spectra, as inferred after correction for absorption in the extragalactic background light (EBL), has interesting implications for the acceleration and radiation mechanisms acting in blazars. A key issue so far has been the dependance of the hardness of the $γ$-ray spectrum on different existing EBL models. The recent {\it Fermi} observations of Mkn~501 now provide additional evidence for the presence of hard intrinsic $γ$-ray spectra independent of EBL uncertainties. Relativistic Maxwellian-type electron energy distributions that are formed in stochastic acceleration scenarios offer a plausible interpretation for such hard source spectra. Here we show that the combined emission from different components with Maxwellian-type distributions could in principle also account for more softer and broader, power law-like emission spectra. We introduce a "leading blob" scenario, applicable to active flaring episodes, when one (or few) of these components become distinct over the "background" emission, producing hard spectral features and/or hardening of the observed spectra. We show that this model can explain the peculiar high-energy characteristics of Mkn~501 in 2009, with evidence for flaring activity and strong spectral hardening at the highest $γ$-ray energies. △ Less

Submitted 23 August, 2011; originally announced August 2011.

Comments: 19 pages, 5 figures, submitted to ApJL

Post-Periastron Gamma Ray Flare from PSR B1259-63/LS 2883 as a Result of Comptonization of the Cold Pulsar Wind

Authors: D. Khangulyan, F. A. Aharonian, S. V. Bogovalov, M. Ribo

Abstract: We argue that the bright flare of the binary pulsar \object{PSR B1259$-$63/LS2883} detected by the {\it Fermi} Large Area Telescope (LAT), is due to the inverse Compton (IC) scattering of the unshocked electron-positron pulsar wind with a Lorentz factor $Γ_0 \approx 10^4$. The combination of two effects both linked to the circumstellar disk (CD), is a key element in the proposed model. The first e… ▽ More We argue that the bright flare of the binary pulsar \object{PSR B1259$-$63/LS2883} detected by the {\it Fermi} Large Area Telescope (LAT), is due to the inverse Compton (IC) scattering of the unshocked electron-positron pulsar wind with a Lorentz factor $Γ_0 \approx 10^4$. The combination of two effects both linked to the circumstellar disk (CD), is a key element in the proposed model. The first effect is related to the impact of the surrounding medium on the termination of the pulsar wind. Inside the disk, the "early" termination of the wind results in suppression of its gamma-ray luminosity. When the pulsar escapes the disk, the conditions for termination of the wind undergo significant changes. This would lead to a dramatic increase of the pulsar wind zone, and thus to the proportional increase of the gamma-ray flux. On the other hand, if the parts of the CD disturbed by the pulsar can supply infrared photons of density high enough for efficient Comptonization of the wind, almost the entire kinetic energy of the pulsar wind would be converted to radiation, thus the gamma-ray luminosity of the wind could approach to the level of the pulsar's spin-down luminosity as reported by the {\it Fermi} collaboration. △ Less

Submitted 8 May, 2012; v1 submitted 24 July, 2011; originally announced July 2011.

Comments: 14 pages, 4 figures

Journal ref: Khangulyan et al. 2012 ApJ 752 L17

Modelling interaction of relativistic and non-relativistic winds in binary system PSR B1259-63/SS2883 - II. Impact of magnetization and anisotropy of the pulsar wind

Authors: S. Bogovalov, D. Khangulyan, A. V. Koldoba, G. V. Ustyugova, F. A. Aharonian

Abstract: In this paper, we present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with its environment in binary system. We compare the impact of both the magnetic field and the wind anisotropy to the benchmark case of a purely hydrodynamical (HD) interaction of isotropic winds, which has been studied in detail by Bogovalov et al. (2008). W… ▽ More In this paper, we present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with its environment in binary system. We compare the impact of both the magnetic field and the wind anisotropy to the benchmark case of a purely hydrodynamical (HD) interaction of isotropic winds, which has been studied in detail by Bogovalov et al. (2008). We consider the interaction in axisymmetric approximation, i.e. the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star and the effects related to the pulsar orbiting are neglected. The impact of the magnetic field is studied for the case of weak magnetization (with magnetization parameter $σ<0.1$), which is consistent with conventional models of pulsar winds. The effects related to anisotropy in pulsar winds are modeled assuming that the kinetic energy flux in a non-magnetized pulsar wind is strongly anisotropic, with the minimum at the pulsar rotation axis and the maximum in the perpendicular direction. We show that, although both considered effects change the shape of the region occupied by the terminated pulsar wind, their impact appears to be small. In particular, for the magnetization of the pulsar wind below 0.1, the magnetic field pressure remains well below the plasma pressure in the post-shock region. Thus, in the case of interaction of a pulsar with the stellar wind environment (opposite to the case of plerions, i.e. the pulsar interaction with the interstellar medium, when the magnetic field becomes dynamically important independently on the wind magnetization) the HD approach represents a feasible approximation for numerical modelling. △ Less

Submitted 24 July, 2011; originally announced July 2011.

Comments: 14 pages, 5 figures, submitted to MNRAS

The need for a local source of UHE CR nuclei

Authors: Andrew M. Taylor, Markus Ahlers, Felix A. Aharonian

Abstract: Recent results of the Pierre Auger (Auger) fluorescence detectors indicate an increasingly heavy composition of ultra-high energy (UHE) cosmic rays (CRs). Assuming that this trend continues up to the highest energies observed by the Auger surface detectors we derive the constraints this places on the local source distribution of UHE CR nuclei. Utilizing an analytic description of UHE CR propagatio… ▽ More Recent results of the Pierre Auger (Auger) fluorescence detectors indicate an increasingly heavy composition of ultra-high energy (UHE) cosmic rays (CRs). Assuming that this trend continues up to the highest energies observed by the Auger surface detectors we derive the constraints this places on the local source distribution of UHE CR nuclei. Utilizing an analytic description of UHE CR propagation we derive the expected spectra and composition for a wide range of source emission spectra. We find that sources of intermediate-to-heavy nuclei are consistent with the observed spectra and composition data above the ankle. This consistency requires the presence of nearby sources within 60 Mpc and 80 Mpc for silicon and iron only sources, respectively. The necessity of these local sources becomes even more compelling in the presence nano-Gauss local extragalactic magnetic fields. △ Less

Submitted 13 February, 2012; v1 submitted 11 July, 2011; originally announced July 2011.

Comments: 10 pages, 8 figures

Journal ref: Phys.Rev.D84:105007,2011

Modeling the hard TeV spectra of blazars 1ES 0229+200 and 3C 66A with an internal absorption scenario

Authors: Olga Zacharopoulou, Dmitry Khangulyan, Felix A. Aharonian, Luigi Costamante

Abstract: We study the applicability of the idea of internal absorption of gamma-rays produced through synchrotron radiation of ultrarelativistic protons in highly magnetized blobs to 1ES 0229+200 and 3C 66A, the two TeV blazars which show unusually hard intrinsic gamma-ray spectra after being corrected for the intergalactic absorption. We show that for certain combinations of reasonable model parameters, e… ▽ More We study the applicability of the idea of internal absorption of gamma-rays produced through synchrotron radiation of ultrarelativistic protons in highly magnetized blobs to 1ES 0229+200 and 3C 66A, the two TeV blazars which show unusually hard intrinsic gamma-ray spectra after being corrected for the intergalactic absorption. We show that for certain combinations of reasonable model parameters, even with quite modest energy requirements, the scenario allows a self-consistent explanation of the non-thermal emission of these objects in the keV, GeV, and TeV energy bands. △ Less

Submitted 15 June, 2011; originally announced June 2011.

Comments: 34 pages, 5 figures, 1 table, accepted for publication in ApJ

Mechanics and kinetics in the Friedmann-Lemaitre-Robertson-Walker space-times

Authors: S. R. Kelner, A. Yu. Prosekin, F. A. Aharonian

Abstract: Using the standard canonical formalism, the equations of mechanics and kinetics in the Friedmann-Lemaitre-Robertson-Walker (FLRW) space-times in Cartesian coordinates have been obtained. The transformation law of the generalized momentum under the shift of the origin of the coordinate system has been found, and the form invariance of the Hamiltonian function relative to the shift transformation ha… ▽ More Using the standard canonical formalism, the equations of mechanics and kinetics in the Friedmann-Lemaitre-Robertson-Walker (FLRW) space-times in Cartesian coordinates have been obtained. The transformation law of the generalized momentum under the shift of the origin of the coordinate system has been found, and the form invariance of the Hamiltonian function relative to the shift transformation has been proved. The general solution of the collisionless Boltzmann equation has been found. In the case of the homogeneous distribution the solutions of the kinetic equation for several simple, but important for applications, cases have been obtained. △ Less

Submitted 11 May, 2011; originally announced May 2011.

Journal ref: Phys.Rev.D84:044016,2011

Non-variable cosmologically distant gamma-ray emitters as an imprint of propagation of ultra-high-energy protons

Authors: A. Yu. Prosekin, S. R. Kelner, F. A. Aharonian

Abstract: The acceleration cites of ultra-high-energy (UHE) protons can be traced by the footprint left by these particles propagating through cosmic microwave background (CMB) radiation. Secondary electrons produced in extended region of several tens of Mpc emit their energy via synchrotron radiation predominantly in the initial direction of the parent protons. It forms a non-variable and compact (almost p… ▽ More The acceleration cites of ultra-high-energy (UHE) protons can be traced by the footprint left by these particles propagating through cosmic microwave background (CMB) radiation. Secondary electrons produced in extended region of several tens of Mpc emit their energy via synchrotron radiation predominantly in the initial direction of the parent protons. It forms a non-variable and compact (almost point-like) source of high energy gamma rays. The importance of this effect is increased for cosmologically distant objects; because of severe energy losses, UHE protons cannot achieve us even in the case of extremely weak intergalactic magnetic fields. Moreover, at high redshifts the energy conversion from protons to secondary particles becomes significantly more effective due to the denser and more energetic CMB in the past. This increases the chances of UHE cosmic rays to be traced by the secondary synchrotron gamma radiation. We discuss the energy budget and the redshift dependence of the efficiency of energy transfer from UHE protons to synchrotron radiation. The angular and spectral distributions of radiation in the gamma- and X-ray energy bands are calculated and discussed in the context of their detectability by Fermi LAT and Chandra observatories. △ Less

Submitted 10 May, 2011; originally announced May 2011.

Molecular Clouds as Cosmic Ray Laboratories

Authors: Sabrina Casanova, Felix A. Aharonian, Yasuo Fukui, Stefano Gabici, David I. Jones, Akiko Kawamura, Toshikazu Onishi, Gavin Rowell, Hidetoshi Sano, Kazufumi Torii, Francesca Volpe, Hiroaki Yamamoto

Abstract: We will here discuss how the gamma-ray emission from molecular clouds can be used to probe the cosmic ray flux in distant regions of the Galaxy and to constrain the highly unknown cosmic ray diffusion coefficient. In particular we will discuss the GeV to TeV emission from runaway cosmic rays penetrating molecular clouds close to young and old supernova remnants and in molecular clouds illuminated… ▽ More We will here discuss how the gamma-ray emission from molecular clouds can be used to probe the cosmic ray flux in distant regions of the Galaxy and to constrain the highly unknown cosmic ray diffusion coefficient. In particular we will discuss the GeV to TeV emission from runaway cosmic rays penetrating molecular clouds close to young and old supernova remnants and in molecular clouds illuminated by the background cosmic ray flux. △ Less

Submitted 6 April, 2011; originally announced April 2011.

Comments: to appear on Proceedings of 25th Texas Symposium on Relativistic Astrophysics

Astrophysical parameters of LS2883 and implications for the PSR B1259-63 gamma-ray binary

Authors: Ignacio Negueruela, Marc Ribó, Artemio Herrero, Javier Lorenzo, Dmitry Khangulyan, Felix A. Aharonian

Abstract: Only a few binary systems with compact objects display TeV emission. The physical properties of the companion stars represent basic input to understand the physical mechanisms behind the particle acceleration, emission, and absorption processes in these so-called gamma-ray binaries. Here we present high-resolution and high signal-to-noise optical spectra of LS2883, the Be star forming a gamma-ray… ▽ More Only a few binary systems with compact objects display TeV emission. The physical properties of the companion stars represent basic input to understand the physical mechanisms behind the particle acceleration, emission, and absorption processes in these so-called gamma-ray binaries. Here we present high-resolution and high signal-to-noise optical spectra of LS2883, the Be star forming a gamma-ray binary with the young non-accreting pulsar PSR B1259-63, showing it to rotate faster and be significantly earlier and more luminous than previously thought. Analysis of the interstellar lines suggest that the system is located at the same distance as (and thus is likely a member of) CenOB1. Taking the distance to the association, d=2.3kpc, and a color excess of E(B-V)=0.85 for LS2883, results in Mv=-4.4. Because of fast rotation, LS2883 is oblate (R_eq=9.7R_sun and R_pole=8.1R_sun) and presents a temperature gradient (T_eq=27500K, log g_eq=3.7; T_pole=34000K, log g_pole=4.1). If the star did not rotate, it would have parameters corresponding to a late O-type star. We estimate its luminosity at log(L*/L_sun)=4.79, and its mass at M=30M_sun. The mass function then implies an inclination of the binary system i_orb=23deg,slightly smaller than previous estimates. We discuss the implications of these new astrophysical parameters of LS2883 for the production of high energy and very high energy gamma rays in the PSR B1259-63/LS2883 gamma-ray binary system. In particular, the stellar properties are very important for prediction of the line-like bulk Comptonization component from the unshocked ultra-relativistic pulsar wind. △ Less

Submitted 23 March, 2011; originally announced March 2011.

Comments: Accepted for publication in ApJL

Journal ref: ApJL 732, L11 (2011)

Rapid TeV variability in Blazars as result of Jet-Star Interaction

Authors: Maxim V. Barkov, Felix A. Aharonian, Sergey V. Bogovalov, Stanislav R. Kelner, Dmitriy V. Khangulyan

Abstract: We propose a new model for the description of ultra-short flares from TeV blazars by compact magnetized condensations (blobs), produced when red giant stars cross the jet close to the central black hole. Our study includes a simple dynamical model for the evolution of the envelope lost by the star in the jet, and its high energy nonthermal emission through different leptonic and hadronic radiation… ▽ More We propose a new model for the description of ultra-short flares from TeV blazars by compact magnetized condensations (blobs), produced when red giant stars cross the jet close to the central black hole. Our study includes a simple dynamical model for the evolution of the envelope lost by the star in the jet, and its high energy nonthermal emission through different leptonic and hadronic radiation mechanisms. We show that the fragmented envelope of the star can be accelerated to Lorentz factors up to 100 and radiate effectively the available energy in gamma-rays predominantly through proton synchrotron radiation or external inverse Compton scattering of electrons. The model can readily explain the minute-scale TeV flares on top of longer (typical time-scales of days) gamma-ray variability as observed from the blazar PKS 2155-304. In the framework of the proposed scenario, the key parameters of the source are robustly constrained. In the case of proton synchrotron origin of the emission a mass of the central black hole of $M_{\rm BH}\approx 10^8 M_{\odot}$, a total jet power of $L_{\rm j} \approx 2\times 10^{47} \, \rm erg\,s^{-1}$ and a Doppler factor, of the gamma-ray emitting blobs, of $δ\geq 40$ are required. Whilst for the external inverse Compton model, parameters of $M_{\rm BH}\approx 10^8 M_{\odot}$, $L_{\rm j} \approx 10^{46} \, \rm erg\,s^{-1}$ and the $δ\geq 150$ are required. △ Less

Submitted 23 December, 2011; v1 submitted 8 December, 2010; originally announced December 2010.

Comments: 25 pages, 11 figures, Submitted to ApJ

Radioactivity and Electron Acceleration in Supernova Remnants

Authors: V. N. Zirakashvili, F. A. Aharonian

Abstract: We argue that the decays of radioactive nuclei related to $^{44}$Ti and $^{56}$Ni ejected during supernova explosions can provide a vast pool of mildly relativistic positrons and electrons which are further accelerated to ultrarelativistic energies by reverse and forward shocks. This interesting link between two independent processes - the radioactivity and the particle acceleration - can be a clu… ▽ More We argue that the decays of radioactive nuclei related to $^{44}$Ti and $^{56}$Ni ejected during supernova explosions can provide a vast pool of mildly relativistic positrons and electrons which are further accelerated to ultrarelativistic energies by reverse and forward shocks. This interesting link between two independent processes - the radioactivity and the particle acceleration - can be a clue for solution of the well known theoretical problem of electron injection in supernova remnants. In the case of the brightest radio source Cas A, we demonstrate that the radioactivity can supply adequate number of energetic electrons and positrons for interpretation of observational data provided that they are stochastically pre-accelerated in the upstream regions of the forward and reverse shocks. △ Less

Submitted 7 October, 2011; v1 submitted 22 November, 2010; originally announced November 2010.

Comments: 6 pages, 1 figure, revised version accepted to Phys.Rev. D

Constraints on the cosmic ray diffusion coefficient in the W28 region from gamma-ray observations

Authors: S. Gabici, S. Casanova, F. A. Aharonian, G. Rowell

Abstract: GeV and TeV gamma rays have been detected from the supernova remnant W28 and its surroundings. Such emission correlates quite well with the position of dense and massive molecular clouds and thus it is often interpreted as the result of hadronic cosmic ray interactions in the dense gas. Constraints on the cosmic ray diffusion coefficient in the region can be obtained, under the assumption that the… ▽ More GeV and TeV gamma rays have been detected from the supernova remnant W28 and its surroundings. Such emission correlates quite well with the position of dense and massive molecular clouds and thus it is often interpreted as the result of hadronic cosmic ray interactions in the dense gas. Constraints on the cosmic ray diffusion coefficient in the region can be obtained, under the assumption that the cosmic rays responsible for the gamma ray emission have been accelerated in the past at the supernova remnant shock, and subsequently escaped in the surrounding medium. In this scenario, gamma ray observations can be explained only if the diffusion coefficient in the region surrounding the supernova remnant is significantly suppressed with respect to the average galactic one. △ Less

Submitted 27 September, 2010; originally announced September 2010.

Comments: To appear in the proceedings of "Journées de la SF2A 2010" Marseille 21-24 June 2010, 4 pages, 4 figures

Galactic center at very high-energies

Authors: M. Chernyakova, D. Malyshev, F. A. Aharonian, R. M. Crocker, D. I. Jones

Abstract: Employing data collected during the first 25 months' observations by the Fermi-LAT, we describe and subsequently seek to model the very high energy (>300 MeV) emission from the central few parsecs of our Galaxy. We analyze the morphological, spectral and temporal characteristics of the central source, 1FGL J1745.6-2900. Remarkably, the data show a clear, statistically significant signal at energie… ▽ More Employing data collected during the first 25 months' observations by the Fermi-LAT, we describe and subsequently seek to model the very high energy (>300 MeV) emission from the central few parsecs of our Galaxy. We analyze the morphological, spectral and temporal characteristics of the central source, 1FGL J1745.6-2900. Remarkably, the data show a clear, statistically significant signal at energies above 10 GeV, where the Fermi-LAT has an excellent angular resolution comparable to the angular resolution of HESS at TeV energies, which makes meaningful the joint analysis of the Fermi and HESS data. Our analysis does not show statistically significant variability of 1FGL J1745.6-2900. Using the combination of Fermi data on 1FGL J1745.6-2900 and HESS data on the coincident, TeV source HESS J1745-290, we show that the spectrum of the central gamma-ray source is inflected with a relatively steep spectral region matching between the flatter spectrum found at both low and high energies. We seek to model the gamma-ray production in the inner 10 pc of the Galaxy and examine, in particular, cosmic ray (CR) proton propagation scenarios that reproduce the observed spectrum of the central source. We show that a model that instantiates a transition from diffusive propagation of the CR protons at low energy to almost rectilinear propagation at high energies (given a reasonable energy-dependence of the assumed diffusion coefficient) can well explain the spectral phenomenology. In general, however, we find considerable degeneracy between different parameter choices which will only be broken with the addition of morphological information that gamma-ray telescopes cannot deliver given current angular resolution limits.We argue that a future analysis done in combination with higher-resolution radio continuum data holds out the promise of breaking this degeneracy. △ Less

Submitted 14 September, 2010; originally announced September 2010.

Comments: submitted to ApJ

Journal ref: The Astrophysical Journal, Volume 726, Issue 2, article id. 60 (2011)

A simple model for electron plasma heating in supernova remnants

Authors: D. Malyshev, S. Gabici, L. O'C. Drury, F. A. Aharonian

Abstract: Context: Multiwavelength observations of supernova remnants can be explained within the framework of diffusive shock acceleration theory, which allows effective conversion of the explosion energy into cosmic rays. Although the models of nonlinear shocks describe reasonably well the nonthermal component of emission, certain issues, including the heating of the thermal electron plasma and the relate… ▽ More Context: Multiwavelength observations of supernova remnants can be explained within the framework of diffusive shock acceleration theory, which allows effective conversion of the explosion energy into cosmic rays. Although the models of nonlinear shocks describe reasonably well the nonthermal component of emission, certain issues, including the heating of the thermal electron plasma and the related X-ray emission, still remain open. Methods: Numerical solution of the equations of the Chevalier model for supernova remnant evolution, coupled with Coulomb scattering heating of the electrons. Results: The electron temperature and the X-ray thermal Bremsstrahlung emission from supernova remnants have been calculated as functions of the relevant parameters. Since only the Coulomb mechanism was considered for electron heating, the values obtained for the electron temperatures should be treated as lower limits. Results from this work can be useful to constrain model parameters for observed SNRs. △ Less

Submitted 7 July, 2010; v1 submitted 6 July, 2010; originally announced July 2010.

Comments: Accepted to A&A as a research note

Angular, spectral, and time distributions of highest energy protons and associated secondary gamma-rays and neutrinos propagating through extragalactic magnetic and radiation fields

Authors: F. A. Aharonian, S. R. Kelner, A. Y. Prosekin

Abstract: The angular, spectral and temporal features of the highest energy protons and accompanying them secondary neutrinos and synchrotron gamma-rays propagating through the intergalactic magnetic and radiation fields are studied using the analytical solutions of the Boltzmann transport equation obtained in the limit of the small-angle and continuous-energy-loss approximation. The angular, spectral and temporal features of the highest energy protons and accompanying them secondary neutrinos and synchrotron gamma-rays propagating through the intergalactic magnetic and radiation fields are studied using the analytical solutions of the Boltzmann transport equation obtained in the limit of the small-angle and continuous-energy-loss approximation. △ Less

Submitted 5 June, 2010; originally announced June 2010.

Comments: 21 pages, 13 figures

Journal ref: Phys.Rev.D82:043002,2010

Gamma-ray flares from red giant/jet interactions in AGN

Authors: Maxim V. Barkov, Felix A. Aharonian, Valenti Bosch-Ramon

Abstract: Non-blazar AGN have been recently established as a class of gamma-ray sources. M87, a nearby representative of this class, show fast TeV variability on timescales of a few days. We suggest a scenario of flare gamma-ray emission in non-blazar AGN based on a red giant interacting with the jet at the base. We solve the hydrodynamical equations that describe the evolution of the envelope of a red gian… ▽ More Non-blazar AGN have been recently established as a class of gamma-ray sources. M87, a nearby representative of this class, show fast TeV variability on timescales of a few days. We suggest a scenario of flare gamma-ray emission in non-blazar AGN based on a red giant interacting with the jet at the base. We solve the hydrodynamical equations that describe the evolution of the envelope of a red giant blown by the impact of the jet. If the red giant is at least slightly tidally disrupted by the supermassive black hole, enough stellar material will be blown by the jet, expanding quickly until a significant part of the jet is shocked. This process can render suitable conditions for energy dissipation and proton acceleration, which could explain the detected day-scale TeV flares from M87 via proton-proton collisions. Since the produced radiation would be unbeamed, such an events should be mostly detected from non-blazar AGN. They may be frequent phenomena, detectable in the GeV-TeV range even up to distances of $\sim 1$ Gpc for the most powerful jets. The counterparts at lower energies are expected to be not too bright.} {M87, and nearby non-blazar AGN in general, can be fast variable sources of gamma-rays through red giant/jet interactions. △ Less

Submitted 3 September, 2010; v1 submitted 28 May, 2010; originally announced May 2010.

Comments: 8 pages, 4 figures

Journal ref: Astrophys.J.724:1517-1523,2010

AGILE detection of GeV gamma-ray emission from the SNR W28

Authors: A. Giuliani, M. Tavani, A. Bulgarelli, E. Striani, S. Sabatini, M. Cardillo, Y. Fukui, A. Kawamura, A. Ohama, N. Furukawa, K. Torii, F. A. Aharonian, F. Verrecchia, A. Argan, G. Barbiellini, P. A. Caraveo, P. W. Cattaneo, A. W. Chen, V. Cocco, E. Costa, F. D'Ammando, E. Del Monte, G. De Paris, G. Di Cocco, I. Donnarumma , et al. (37 additional authors not shown)

Abstract: Supernova remnants (SNRs) are believed to be the main sources of Galactic cosmic rays. Molecular clouds associated with SNRs can produce gamma-ray emission through the interaction of accelerated particles with the concentrated gas. The middle aged SNR W28, for its associated system of dense molecular clouds, provides an excellent opportunity to test this hypothesis. We present the AGILE/GRID obse… ▽ More Supernova remnants (SNRs) are believed to be the main sources of Galactic cosmic rays. Molecular clouds associated with SNRs can produce gamma-ray emission through the interaction of accelerated particles with the concentrated gas. The middle aged SNR W28, for its associated system of dense molecular clouds, provides an excellent opportunity to test this hypothesis. We present the AGILE/GRID observations of SNR W28, and compare them with observations at other wavelengths (TeV and 12CO J=1-->0 molecular line emission). The gamma-ray flux detected by AGILE from the dominant source associated with W28 is (14 +- 5) 10^-8 ph cm^-2 s^-1 for E > 400 MeV. This source is positionally well correlated with the TeV emission observed by the HESS telescope. The local variations of the GeV to TeV flux ratio suggest a difference between the CR spectra of the north-west and south molecular cloud complexes. A model based on a hadronic-induced interaction and diffusion with two molecular clouds at different distances from the W28 shell can explain both the morphological and spectral features observed by AGILE in the MeV-GeV energy range and by the HESS telescope in the TeV energy range. The combined set of AGILE and H.E.S.S. data strongly support a hadronic model for the gamma-ray production in W28. △ Less

Submitted 5 May, 2010; originally announced May 2010.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

Journal ref: Astron.Astrophys.516:L11,2010

Modeling the gamma-ray emission produced by runaway cosmic rays in the environment of RX J1713.7-3946

Authors: S. Casanova, D. I. Jones, F. A. Aharonian, Y. Fukui, S. Gabici, A. Kawamura, T. Onishi, G. Rowell, H. Sano, K. Torii, H. Yamamoto

Abstract: Diffusive shock acceleration in supernova remnants is the most widely invoked paradigm to explain the Galactic cosmic ray spectrum. Cosmic rays escaping supernova remnants diffuse in the interstellar medium and collide with the ambient atomic and molecular gas. From such collisions gamma-rays are created, which can possibly provide the first evidence of a parent population of runaway cosmic rays.… ▽ More Diffusive shock acceleration in supernova remnants is the most widely invoked paradigm to explain the Galactic cosmic ray spectrum. Cosmic rays escaping supernova remnants diffuse in the interstellar medium and collide with the ambient atomic and molecular gas. From such collisions gamma-rays are created, which can possibly provide the first evidence of a parent population of runaway cosmic rays. We present model predictions for the GeV to TeV gamma-ray emission produced by the collisions of runaway cosmic rays with the gas in the environment surrounding the shell-type supernova remnant RX J1713.7-3946. The spectral and spatial distributions of the emission, which depend upon the source age, the source injection history, the diffusion regime and the distribution of the ambient gas, as mapped by the LAB and NANTEN surveys, are studied in detail. In particular, we find for the region surrounding RX J1713-3946, that depending on the energy one is observing at, one may observe startlingly different spectra or may not detect any enhanced emission with respect to the diffuse emission contributed by background cosmic rays. This result has important implications for current and future gamma-ray experiments. △ Less

Submitted 20 October, 2010; v1 submitted 1 March, 2010; originally announced March 2010.

Comments: version published on PASJ

Journal ref: Publ. Astron. Soc. Japan 62, pp.1127-1134 (2010)

Acceleration and radiation of ultra-high energy protons in galaxy clusters

Authors: G. Vannoni, F. A. Aharonian, S. Gabici, S. R. Kelner, A. Prosekin

Abstract: Clusters of galaxies are believed to be capable to accelerate protons at accretion shocks to energies exceeding 10^18 eV. At these energies, the losses caused by interactions of cosmic rays with photons of the Cosmic Microwave Background Radiation (CMBR) become effective and determine the maximum energy of protons and the shape of the energy spectrum in the cutoff region. The aim of this work is… ▽ More Clusters of galaxies are believed to be capable to accelerate protons at accretion shocks to energies exceeding 10^18 eV. At these energies, the losses caused by interactions of cosmic rays with photons of the Cosmic Microwave Background Radiation (CMBR) become effective and determine the maximum energy of protons and the shape of the energy spectrum in the cutoff region. The aim of this work is the study of the formation of the energy spectrum of accelerated protons at accretion shocks of galaxy clusters and of the characteristics of their broad band emission. The proton energy distribution is calculated self-consistently via a time-dependent numerical treatment of the shock acceleration process which takes into account the proton energy losses due to interactions with the CMBR. We calculate the energy distribution of accelerated protons, as well as the flux of broad-band emission produced by secondary electrons and positrons via synchrotron and inverse Compton scattering processes. We find that the downstream and upstream regions contribute almost at the same level to the emission. For the typical parameters characterising galaxy clusters, the synchrotron and IC peaks in the spectral energy distributions appear at comparable flux levels. For an efficient acceleration, the expected emission components in the X-ray and gamma-ray band are close to the detection threshold of current generation instruments, and will be possibly detected with the future generation of detectors. △ Less

Submitted 29 October, 2009; originally announced October 2009.

Comments: 9 pages, 13 figures, submitted to Astronomy & Astrophysics