Source author record

Xiang-Yu Wang

Xiang-Yu Wang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

58works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

58 published item(s)

preprint2025arXiv

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Black holes (BH), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to tera-electronvolt (TeV) gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectrum extending to the ultrahigh-energy (UHE; photon energy $E>100$ TeV) band and one microquasar (Cygnus X-1) of spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory (LHAASO). Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles out of 1 peta-electronvolt (PeV), suggesting microquasars to be important contributors to Galactic cosmic rays especially around the `knee' region.

preprint2024arXiv

The BOAT GRB 221009A: a Poynting-Flux-Dominated Narrow Jet Surrounded by a Matter-Dominated Structured Jet Wing

We argue that the broad-band observations of the brightest-of-all-time (BOAT) GRB 221009A reveal a physical picture involving two jet components: a narrow ($\sim 0.6$ degree half opening angle) pencil-beam jet that has a Poynting-flux-dominated jet composition, and a broader matter-dominated jet with an angular structure. We discuss various observational evidence that supports such a picture. To treat the problem, we develop an analytical structured jet model for both forward and reverse shock emission from the matter dominated structured jet wing during the deceleration phase. We discuss the physical implications of such a two-component jet configuration for this particular burst and for GRBs in general. We argue that some bright X-ray flares could be similar narrow jets viewed slightly outside the narrow jet cone and that narrow jets may exist in many more GRBs without being detected.

preprint2022arXiv

A semi-analytical solution to the forward-reverse shock hydrodynamics of the gamma-ray burst afterglow

We extend the standard model of forward-reverse shock (FS-RS) for gamma-ray burst (GRB) afterglow to more general cases. On one hand, we derive the analytical solution to the hydrodynamics of the shocks in two limiting cases, i.e., an ultra-relativistic reverse shock case and a Newtonian reverse shock case. Based on the asymptotic solutions in these two limiting cases, we constitute a semi-analytical solution for the hydrodynamics of the shocks in the generic case, covering the mildly-relativistic reverse shock case. On the other hand, we derive the evolution of the system taking into account the condition of energy conservation which is not satisfied in the standard FS-RS model. A generic solution of semi-analytical expressions is also given. In both the extended standard FS-RS model (satisfying pressure balance condition) and the model satisfying energy conservation, we find that the results in the ultra-relativistic reverse shock case and in the early stage of the Newtonian reverse shock case are different from those in the standard FS-RS model by only a factor that close to one while the same initial conditions adopted. However, the asymptotic solutions in the limiting cases are not good approximations to those in the intermediate case. Our semi-analytical results agree well with the numerical results for a large range of model parameters, and hence can be easily employed to diagnose the physical quantities of the GRB shell and circumburst environment.

preprint2022arXiv

Constraining the baryon loading factor of AGN jets: implication from the gamma-ray emission of the Coma cluster

High-energy cosmic rays (CRs) can be accelerated in the relativistic jets of Active Galactic Nuclei (AGNs) powered by supermassive black holes. The baryon loading efficiency onto relativistic CR baryons from the accreting black holes is poorly constrained by observations so far. In this paper, we suggest that the $γ$-ray emission of galaxy clusters can be used to study the baryon loading factor of AGN jets, since CRs injected by AGN jets are completely confined in the galaxy clusters and sufficiently interact with intra-cluster medium via hadronic process, producing diffuse $γ$-rays. We study the propagation of CRs in the galaxy clusters and calculate the radial distribution of the gamma-rays in the galaxy cluster with different injection rates from AGNs. By comparison with the $γ$-ray flux and upper limits of the Coma cluster measured by $Fermi$-LAT and VERITAS, we find the upper limit of the average baryon loading factor {(defined as the efficiency with which the gravitational energy is converted into relativistic particles)} to be $η_{p, \mathrm{grav}} < 0.1$. The upper limit is much lower than that required to account for diffuse neutrino flux in the conventional blazar models.

preprint2022arXiv

Detection of diffuse gamma-ray emission towards a massive star forming region hosting Wolf-Rayet stars

Isotopic and elemental abundances seen in Galactic cosmic rays imply that $\sim20\%$ of the cosmic-ray (CR) nuclei are probably synthesized by massive Wolf-Rayet (WR) stars. Massive star clusters hosting WR and OB-type stars have been proposed as potential Galactic cosmic-ray accelerators for decades, in particular via diffusive shock acceleration at wind termination shocks. Here we report the analysis of {\em Fermi} Large Area Telescope's data towards the direction of Masgomas-6a, a young massive star cluster candidate hosting two WR stars. We detect an extended $γ$-ray source with $\rm{TS}=183$ in the vicinity of Masgomas-6a, spatially coincident with two unassociated {\em Fermi} 4FGL sources. We also present the CO observational results of molecular clouds in this region, using the data from the Milky Way Imaging Scroll Painting project. The $γ$-ray emission intensity correlates well with the distribution of molecular gas at the distance of Masgomas-6a, indicating that these gamma rays may be produced by CRs accelerated by massive stars in Masgomas-6a. At the distance of $3.9{\rm \ kpc}$ of Masgomas-6a, the luminosity of the extended source is $(1.81\pm0.02)\times 10^{35}{\rm \ erg \ s^{-1}}$. With a kinetic luminosity of $\sim 10^{37}{\rm erg \ s^{-1}}$ in the stellar winds, the WR stars are capable of powering the $γ$-ray emission via neutral pion decay resulted from cosmic ray $pp$ interactions. The size of the GeV source and the energetic requirement suggests a CR diffusion coefficient smaller than that in the Galactic interstellar medium, indicating strong suppression of CR diffusion in the molecular cloud.

preprint2022arXiv

Prospect of detecting TeV halos with LHAASO: in the framework of the anisotropic diffusion model

The particle diffusion coefficients of three TeV pulsar halos observed so far are inferred to be significantly smaller than the typical value of the interstellar medium (ISM). The anisotropic diffusion model ascribes the slow diffusion to the cross-field diffusion assuming sub-Alfv{é}nic turbulence in the ISM around the pulsar if the viewing angle between the observer's line-of-sight (LOS) to the pulsar and the local mean field direction is small. In general, the TeV halo's morphology under this model highly depends on the viewing angle, and an elongated, asymmetric morphology is predicted if the LOS is not approximately aligned with the local mean field direction. While the specific requirement of a small viewing angle is supposedly established only for a small fraction of TeV halos, TeV halos with apparent asymmetric morphology has not been detected. In this paper we will study the expectation of TeV halos measured by the TeV-PeV gamma-ray detector LHAASO in the framework of anisotropic diffusion model, with a particular focus on the influence of the viewing angle on the detectability. We show that a TeV halo is more detectable with a smaller viewing angle and this selection effect may explain why the morphologies of all three detected TeV halos so far are consistent being spherical. We also demonstrate that LHAASO is capable of detecting asymmetric TeV halos after several-year operation with reasonable source parameters. This can serve as a critical test of the anisotropic diffusion model.

preprint2022arXiv

The Large High Altitude Air Shower Observatory (LHAASO) Science Book (2021 Edition)

Since the science white paper of the Large High Altitude Air Shower Observatory (LHAASO) published on arXiv in 2019 [e-Print: 1905.02773 (astro-ph.HE)], LHAASO has completed the transition from a project to an operational gamma-ray astronomical observatory LHAASO is a new generation multi-component facility located in Daocheng, Sichuan province of China, at an altitude of 4410 meters. It aims at measuring with unprecedented sensitivity the spectrum, composition, and anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$~eV, and acting simultaneously as a wide aperture (one stereoradiant) continuously operating gamma-ray telescope in the energy range between 10$^{11}$ and $10^{15}$~eV with the designed sensitivity of 1.3\% of the Crab Unit (CU) above 100 TeV. LHAASO's capability of measuring simultaneously different shower components (electrons, muons, and Cherenkov/fluorescence light), will allow it to investigate the origin, acceleration, and propagation of CR through measurement of the energy spectrum, elemental composition, and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO will play a key role in CR physics and gamma-ray astronomy for a general and comprehensive exploration of the high energy universe and will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. The LHAASO Collaboration organized an editorial working group and finished all editorial work of this science book, to summarize the instrumental features and outline the prospects of scientific researches with the LHAASO experiment.

preprint2021arXiv

A unified model for orphan and multi-wavelength blazar flares

Blazars are a class of active galactic nuclei which host relativistic jets oriented close to the observer's line of sight. Blazars have very complex variability properties. Flares, namely flux variations around the mean value with a well-defined shape and duration, are one of the identifying properties of the blazar phenomenon. Blazars are known to exhibit multi-wavelength flares, but also "orphan" flares, namely flux changes that appear only in a specific energy range. Various models, sometimes at odds with each other, have been proposed to explain specific flares even for a single source, and cannot be synthesized into a coherent picture. In this paper, we propose a unified model for explaining orphan and multi-wavelength flares from blazars in a common framework. We assume that the blazar emission during a flare consists of two components: (i) a quasi-stable component that arises from the superposition of numerous but comparatively weak dissipation zones along the jet, forming the background (low-state) emission of the blazar, and (ii) a transient component, which is responsible for the sudden enhancement of the blazar flux, forming at a random distance along the jet by a strong energy dissipation event. Whether a multi-wavelength or orphan flare is emitted depends on the distance from the base of the jet where the dissipation occurs. Generally speaking, if the dissipation occurs at a small/large distance from the supermassive black hole, the inverse Compton/synchrotron radiation dominates and an orphan gamma-ray/optical flare tends to appear. On the other hand, we may expect a multi-wavelength flare if the dissipation occurs at an intermediate distance. We show that the model can successfully describe the spectral energy distribution of different flares from the flat spectrum radio quasar 3C 279 and the BL Lac object PKS 2155-304.

preprint2021arXiv

Fermi-LAT detection of extended gamma-ray emission in the vicinity of SNR G045.7$-$00.4: evidence for escaping cosmic rays interacting with the surrounding molecular clouds

We present the analysis of Fermi Large Area Telecope (LAT) data of the gamma-ray emission in the vicinity of a radio supernova remnant (SNR), G045.7$-$00.4. To study the origin of the gamma-ray emission, we also make use of the CO survey data of Milky Way Imaging Scroll Painting to study the massive molecular gas complex that surrounds the SNR. The whole size of the GeV emission is significantly larger than that of the radio morphology. Above 3 GeV, the GeV emission is resolved into two sources: one is spatially consistent with the position of the SNR with a size comparable to that of the radio emission, and the other is located outside of the western boundary of the SNR and spatially coincident with the densest region of the surrounding molecular cloud. We suggest that the GeV emission of the western source may arise from cosmic rays (CRs) which have escaped the SNR and illuminated the surrounding molecular cloud. We find that the gamma-ray spectra of the western source can be consistently explained by this scenario with a total energy of $\sim 10^{50}{\rm erg}$ in escaping CRs assuming the escape is isotropic.

preprint2021arXiv

Prospect of Detecting X-Ray Halos Around Middle-Aged Pulsars with eROSITA

The detection of extended TeV $γ$-ray emission (dubbed "TeV halos") from Geminga and Monogem pulsars by HAWC collaboration implies that the halo-like morphologies around middle-aged pulsars may be common. The $γ$-rays above 10 TeV are thought to arise from inverse Compton (IC) scattering of relativistic electrons/positrons in the pulsar halos off cosmic microwave background photons. In the meanwhile, these electrons and positrons can produce X-ray synchrotron emission in the interstellar magnetic field, resulting in a diffuse emission in the X-ray band (namely X-ray halos). Here, we study the prospect of detecting X-ray halos with eROSITA from 10 middle-aged pulsars with characteristic age larger than tens of thousands of years in the ATNF pulsar catalog. Assuming a benchmark value (i.e., $B = 3 \rm \, μG$) for the magnetic field, most of the X-ray halos are found to be bright enough to be detectable by eROSITA in the energy range of 0.5-2 keV during its four-year all-sky survey. Among these pulsar halos, three are supposed to produce X-ray fluxes above the eROSITA sensitivity of the first all-sky survey. Given the good angular resolution and the large field of view, eROSITA is expected to be able to measure the spatial distribution of the X-ray halos from sub-pc scale up to tens of pc scale. The intensity profiles of the X-ray halos are very useful to constrain the magnetic field and the energy-dependence of the diffusion coefficient in the pulsar halos.

preprint2020arXiv

A serendipitous discovery of GeV gamma-ray emission from supernova 2004dj in a survey of nearby star-forming galaxies with Fermi-LAT

The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time.

preprint2020arXiv

A two-zone blazar radiation model for "orphan" neutrino flares

In this work, we investigate the 2014-2015 neutrino flare associated with the blazar TXS 0506+056 and a recently discovered muon neutrino event IceCube-200107A in spatial coincidence with the blazar 4FGL J0955.1+3551, under the framework of a two-zone radiation model of blazars where an inner/outer blob close to/far from the supermassive black hole are invoked. An interesting feature that the two sources share in common is that no evidence of GeV gamma-ray activity is found during the neutrino detection period, probably implying a large opacity for GeV gamma rays in the neutrino production region. In our model, continuous particle acceleration/injection takes place in the inner blob at the jet base, where the hot X-ray corona of the supermassive black hole provides target photon fields for efficient neutrino production and strong GeV gamma-ray absorption. We show that this model can self-consistently interpret the neutrino emission from both two blazars in a large parameter space. In the meantime, the dissipation processes in outer blob are responsible for the simultaneous multi-wavelength emission of both sources. In agreement with previous studies of TXS 0506+056 and, an intense MeV emission from the induced electromagnetic cascade in the inner blob is robustly expected to accompany the neutrino flare in our model could be used to test the model with the next-generation MeV gamma-ray detector in the future.

preprint2020arXiv

Constraints on the intergalactic magnetic field from $γ$-ray observations of GRB 190114C

Very high energy photons from cosmological gamma-ray bursts (GRBs) are expected to interact with extragalactic background light (EBL) and produce electron-positron pairs when they propagate through intergalactic medium (IGM). These relativistic pairs will then up-scatter cosmic microwave background (CMB) photons and emit secondary GeV emission. Meanwhile, the motion of these pairs are deflected by intergalactic magnetic field (IGMF), so the secondary GeV photons arrive later than the primary emission. It has been suggested that the properties of the secondary GeV emission can be used to constrain IGMF. Recently, TeV gamma-ray emission has been detected, for the first time, from a GRB (GRB 190114C) by the MAGIC telescope and its steep ${\rm γ-ray}$ spectrum shows a clear evidence of absorption by EBL. We then constrain the IGMF with the GeV flux limit obtained from the $Fermi$-LAT observations. We find a limit of $>10^{-19.5}$ G for the coherence length of $λ\leq 1$ Mpc. Although this limit is weaker than that obtained by using blazars, it represents the first limit from ${\rm γ-ray}$ observations of GRBs, which provides an independent constraint on IGMF. We also find that, for transient ${\rm γ-ray}$ sources, one can choose a favorable time window to search for the echo emission at a particular energy.

preprint2020arXiv

Diffuse gamma-ray emission from the vicinity of young massive star cluster RSGC 1

We report the Fermi Large Area Telescope (Fermi-LAT) detection of the gamma-ray emission towards the young massive star cluster RSGC 1. Using the latest source catalog and diffuse background models, we found that the diffuse gamma-ray emission in this region can be resolved into three different components. The GeV gamma-ray emission from the region HESS J1837-069 has a photon index of 1.83 $\pm$ 0.08. Combining with the HESS and MAGIC data, we argue that the gamma-ray emission in this region likely originate from a pulsar wind nebula (PWN). The gamma-ray emission from the northwest part (region A) can be modelled by an ellipse with the semimajor and semiminor axis of 0.5$^{\circ}$ and 0.25$^{\circ}$, respectively. The GeV emission has a hard spectrum with a photon index of 2.05 $\pm$ 0.02 and is partially coincide with the TeV source MAGIC J1835-069. The possible origin of the gamma-ray emission in this region is the interaction of the cosmic rays (CRs) accelerated by SNR G24.7+0.6 or/and the OB cluster G25.18+0.26 with the surrounding gas clouds. The GeV gamma-ray emission from the southeast region (region B) can be modeled as an ellipse with the semimajor and semiminor axis of 0.9$^{\circ}$ and 0.5$^{\circ}$, respectively, and also reveals a hard gamma-ray spectrum. We argue that the most probable origin is the interaction of the accelerated protons in the young massive star cluster RSGC 1 with ambient gas clouds, and the total cosmic-ray (CR) proton energy is estimated to be as high as $\sim$ 1x10$^{50}$ erg.

preprint2020arXiv

Diffuse gamma-ray emission toward the massive star-forming region, W40

We report the detection of high-energy gamma-ray signal towards the young star-forming region, W40. Using 10-year Pass 8 data from the Fermi Large Area Telescope (Fermi-LAT), we extracted an extended gamma-ray excess region with a significance of about 18sigma. The radiation has a spectrum with a photon index of 2.49 +/- 0.01. The spatial correlation with the ionized gas content favors the hadronic origin of the gamma-ray emission. The total cosmic-ray (CR) proton energy in the gamma-ray production region is estimated to be the order of 10^47 erg. However, this could be a small fraction of the total energy released in cosmic rays (CRs) by local accelerators, presumably by massive stars, over the lifetime of the system. If so, W40, together with earlier detections of gamma-rays from Cygnus cocoon, Westerlund 1, Westerlund 2, NGC 3603, and 30 Dor C, supports the hypothesis that young star clusters are effective CR factories. The unique aspect of this result is that the gamma-ray emission is detected, for the first time, from a stellar cluster itself, rather than from the surrounding "cocoons".

preprint2020arXiv

Examining the secondary product origin of cosmic ray positrons with the latest AMS-02 data

Measurements of cosmic-ray (CR) positron fraction by PAMELA and other experiments have found an excess above 10 GeV relative to the standard predictions for secondary production in the interstellar medium (ISM). Although the excess has been mostly suggested to arise from some primary sources of positrons (such as pulsars and or annihilating dark matter particles), the almost constant flux ratio of $e^{+}/ \bar{p}$ argues for an alternative possibility that the excess positrons and antiprotons up to the highest energies are secondary products generated in hadronic interactions. Recently, Yang \& Aharonian (2019) revisit this possibility by assuming the presence of an additional population of CR nuclei sources. Here we examine this secondary product scenario using the \texttt{DRAGON} code, where the radiative loss of positrons is taken into account consistently. We confirm that the CR proton spectrum and the antiproton data can be explained by assuming the presence of an additional population of CR sources. However, the corresponding positron spectrum deviates from the measured data significantly above 100 GeV due to the strong radiative cooling. This suggests that, although hadronic interactions can explain the antiproton data, the corresponding secondary positron flux is still not enough to account for the AMS data. Hence contribution from some primary positron sources, such as pulsars or dark matter, is non-negligible.

preprint2020arXiv

Very High Energy Emission and Cascade Radiation of Gamma-Ray Burst Afterglows: Homogeneous Versus Wind External Media

Recent detection of sub-TeV emission from gamma-ray bursts (GRBs) represents a breakthrough in the GRB study. The multi-wavelength data of the afterglows of GRB 190114C support the synchrotron self-Compton (SSC) origin for its sub-TeV emission. We present a comparative analysis on the SSC emission of GRB afterglows in the homogeneous and wind environment in the framework of the forward shock model. The $γγ$ absorption of very high-energy photons due to pair production within the source and the Klein-Nishina effect on the inverse-Compton scattering are considered. Generally a higher SSC flux is expected for a larger circum-burst density due to a larger Compton parameter, but meanwhile the internal $γγ$ absorption is more severer for sub-TeV emission. The flux ratio between the SSC component and the synchrotron component decreases more quickly with time in the wind medium case than that in the homogenous-density medium case. The light curves of the SSC emission are also different for the two types of media. We also calculate the cascade emission resulted from the absorbed high-energy photons. In the ISM environment with $n> 1\,\rm cm^{-3}$, the cascade synchrotron emission could be comparable to the synchrotron emission of the primary electrons in the optical band, which may flatten the optical afterglow light curve at early time ($t<1$ h). In the wind medium with $A_{\ast}> 0.1$, the cascade emission in the eV-GeV band is comparable or even larger than the emission of the primary electrons at early time.

preprint2019arXiv

A two-zone model for blazar emission: implications for TXS 0506+056 and the neutrino event IceCube-170922A

A high-energy muon neutrino event, IceCube-170922A, was recently discovered in both spatial and temporal coincidence with a gamma-ray flare of the blazar TXS 0506+056. It has been shown, with standard one-zone models, that neutrinos can be produced in the blazar jet via hadronic interactions, but with a flux which is mostly limited by the X-ray data. In this work, we explore the neutrino production from TXS 0506+056 by invoking two physically distinct emission zones in the jet, separated by the broad line region (BLR). Using the Doppler-boosted radiation of the BLR as the target photon field, the inner zone accounts for the neutrino and gamma-ray emission via $pγ$ interactions and inverse Compton scattering respectively, while the outer zone produces the optical and X-ray emission via synchrotron and synchrotron self-Compton processes. The different conditions of the two zones allow us to suppress the X-ray emission from the electromagnetic cascade, and set a much higher upper limit on the muon neutrino flux (i.e., $\sim 10^{-11}\rm erg~cm^{-2}s^{-1}$) than in one-zone models. We compare, in detail, our scenario with one-zone models discussed in the literature, and argue that differentiating between such scenarios will become possible with next generation neutrino telescopes, such as IceCube-Gen2.

preprint2019arXiv

Discovery of a spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099: a TeV halo or supernova remnant ?

Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we search for possible GeV counterparts of three TeV halo candidates with Fermi Large Area Telescopes. We detect a new spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099, which is also detected by HESS (HESS J1912+101). We find that the size of the GeV source is significantly larger than that of the TeV emission measured by HESS, and a spatial template characteristic of a PWN can fit the GeV data. We suggest that the GeV source is an extended PWN powered by the central middle-aged pulsar PSR J1913+1011. This discovery favors the TeV halo scenario for the TeV source 2HWC J1912+099 (HESS J1912+101), although the possible shell-like morphology measured by HESS challenges this interpretation. Alternatively, the TeV emission could be dominated by a supernova remnant (SNR) via the hadronic process. Future multi-wavelength observations of the source and more precise measurements of the spatial profile of the TeV emission will be useful to distinguish between the two scenarios.

preprint2016arXiv

First detection of GeV emission from an ultraluminous infrared galaxy: Arp 220 as seen with the Fermi Large Area Telescope

Cosmic rays (CRs) in starburst galaxies produce high energy gamma-rays by colliding with the dense interstellar medium (ISM). Arp 220 is the nearest ultra luminous infrared galaxy (ULIRG) that has star-formation at extreme levels, so it has long been predicted to emit high-energy gamma-rays. However, no evidence of gamma-ray emission was found despite intense efforts of search. Here we report the discovery of high-energy gamma-ray emission above 200 MeV from Arp 220 at a confidence level of $\sim 6.3 σ$ using 7.5 years of \textsl {Fermi} Large Area Telescope observations. The gamma-ray emission shows no significant variability over the observation period and it is consistent with the quasi-linear scaling relation between the gamma-ray luminosity and total infrared luminosity for star-forming galaxies, suggesting that these gamma-rays arise from CR interactions. As the high density medium of Arp 220 makes it an ideal CR calorimeter, the gamma-ray luminosity can be used to measure the efficiency of powering CRs by supernova (SN) remnants given a known supernova rate in Arp 220. We find that this efficiency is about $4.2\pm2.6\%$ for CRs above 1 GeV.

preprint2016arXiv

How far are the sources of IceCube neutrinos? Constraints from the diffuse TeV gamma-ray background

The nearly isotropic distribution of the TeV-PeV neutrinos recently detected by IceCube suggests that they come from sources at distance beyond our Galaxy, but how far they are is largely unknown due to lack of any associations with known sources. In this paper, we propose that the cumulative TeV gamma-ray emission accompanying the production of neutrinos can be used to constrain the distance of these neutrino sources, since the opacity of TeV gamma rays due to absorption by the extragalactic background light (EBL) depends on the distance that these TeV gamma rays have travelled. As the diffuse extragalactic TeV background measured by \emph{Fermi} is much weaker than the expected cumulative flux associated with IceCube neutrinos, the majority of IceCube neutrinos, if their sources are transparent to TeV gamma rays, must come from distances larger than the horizon of TeV gamma rays. We find that above 80\% of the IceCube neutrinos should come from sources at redshift $z>0.5$. Thus, the chance for finding nearby sources correlated with IceCube neutrinos would be small. We also find that, to explain the flux of neutrinos under the TeV gamma-ray emission constraint, the redshift evolution of neutrino source density must be at least as fast as the the cosmic star-formation rate.

preprint2016arXiv

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

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.

preprint2016arXiv

Modeling the gamma-ray emission in the Galactic Center with a fading cosmic-ray accelerator

Recent HESS observations of the ~200 pc scale diffuse gamma-ray emission from the central molecular zone (CMZ) suggest the presence of a PeV cosmic-ray accelerator (PeVatron) located in the inner 10 pc region of the Galactic Center. Interestingly, the gamma-ray spectrum of the point-like source (HESS J1745-290) in the Galactic Center shows a cutoff at ~10 TeV, implying a cutoff around 100 TeV in the cosmic-ray proton spectrum. Here we propose that the gamma-ray emission from the inner and the outer regions may be explained self-consistently by run-away protons from a single, yet fading accelerator. In this model, gamma rays from the CMZ region are produced by protons injected in the past, while gamma rays from the inner region are produced by protons injected more recently. We suggest that the blast wave formed in a tidal disruption event (TDE) caused by the supermassive black hole (Sgr A*) could serve as such a fading accelerator. With typical parameters of the TDE blast wave, gamma-ray spectra of both the CMZ region and HESS J1745-290 can be reproduced simultaneously. Meanwhile, we find that the cosmic-ray energy density profile in the CMZ region may also be reproduced in the fading accelerator model when appropriate combinations of the particle injection history and the diffusion coefficient of cosmic rays are adopted.

preprint2016arXiv

Probing the environment of gravitational wave transient sources with TeV afterglow emission

Recently, Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) detected gravitational wave (GW) transients from mergers of binary black holes (BHs). The system may also produce a wide-angle, relativistic outflow if the claimed short GRB detected by GBM is in real association with GW 150914. It was suggested that mergers of double neutron stars (or neutron star-black hole binaries), another promising source of GW transients, also produce fast, wide-angle outflows. In this paper, we calculate the high-energy gamma-ray emission arising from the blast waves driven by these wide-angle outflows. We find that TeV emission arising from the inverse-Compton process in the relativistic outflow resulted from mergers of binary BHs similar to those in GW 150914 could be detectable by ground-based IACT telescopes such as Cherenkov Telescope Array (CTA) if the sources occur in { a dense medium with density $n > 0.3 cm^{-3}$}. For neutron star-neutron star (NS-NS) and NS-BH mergers, TeV emission from the wide-angle, mildly-relativistic outflow could be detected as well if they occur in a dense medium {with $n > 10-100 cm^{-3}$.} Thus TeV afterglow emission would be a useful probe of the environment of the GW transients, which could shed light on the evolution channels of the progenitors of GW transients.

preprint2016arXiv

Search for high energy gamma-ray emission from tidal disruption events with the Fermi Large Area Telescope

Massive black holes at galaxy center may tear apart a star when the star passes occasionally within the disruption radius, which is the so-called tidal disruption event(TDE). Most TDEs radiate with thermal emission resulted from the accretion disk, but three TDEs have been detected in bright non-thermal X-ray emission, which is interpreted as arising from the relativistic jets. Search for high-energy gamma-ray emission from one relativistic TDE (Swift J164449.3+573451) with the \textsl{Fermi} Large Area Telescope (LAT) has yielded non-detection. In this paper, we report the search for high energy emission from the other two relativistic TDEs (Swift J2058.4+0516 Swift J1112.2-8238) during the flare period. No significant GeV emission is found, with an upper limit fluence in LAT energy range being less than $1\%$ of that in X-rays. Compared with gamma-ray bursts (GRBs) and blazars, these TDEs have the lowest flux ratio between GeV emission and X-ray emission. The non-detection of high-energy emission from relativistic TDEs could be due to that the high-energy emission is absorbed by soft photons in the source. Based on this hypothesis, upper limits on the bulk Lorentz factors, $Γ\la 30$, are then obtained for the jets in these TDEs. We also search for high-energy gamma-ray emission from the nearest TDE discovered to date, ASASSN-14li. No significant GeV emission is found and an upper limit of $L(\rm 0.1-10 GeV)\le 4.4\times 10^{42}$ erg s$^{-1}$ (at $95\%$ confidence level) is obtained for the first $10^{7}$ s after the disruption.

preprint2016arXiv

Testing the equivalence principle and Lorentz invariance with PeV neutrinos from blazar flares

It was recently proposed that a giant flare of the blazar PKS B1424-418 at redshift $z=1.522$ is in association with a PeV-energy neutrino event detected by IceCube. Based on this association we here suggest that the flight time difference between the PeV neutrino and gamma-ray photons from blazar flares can be used to constrain the violations of equivalence principle (EP) and the Lorentz invariance for neutrinos. From the calculated Shapiro delay due to clusters or superclusters in the nearby universe, we find that violation of the equivalence principle for neutrinos and photons is constrained to an accuracy of at least $10^{-5}$, which is two orders of magnitude tighter than the constraint placed by MeV neutrinos from supernova 1987A. Lorentz invariance violation (LIV) arises in various quantum-gravity theories, which predicts an energy-dependent velocity of propagation in vacuum for photons and neutrinos. We find that the association of the PeV neutrino with the gamma-ray outburst set limits on the energy scale of possible LIV to $> 0.01 E_{pl}$ for linear LIV models and $>6\times 10^{-8} E_{pl}$ for quadratic order LIV models, where $E_{pl}$ is the Planck energy scale. These are the most stringent constraints on neutrino LIV for subluminal neutrinos.

preprint2015arXiv

Measuring the bulk Lorentz factors of gamma-ray bursts with Fermi

Gamma-ray bursts (GRBs) are powered by ultra-relativistic jets. Usually a minimum value of the Lorentz factor of the relativistic bulk motion is obtained based on the argument that the observed high energy photons ($\gg {\rm MeV}$) can escape without suffering from absorption due to pair production. The exact value, rather than a lower limit, of the Lorentz factor can be obtained if the spectral cutoff due to such absorption is detected. With the good spectral coverage of the Large Area Telescope (LAT) on \textsl{Fermi}, measurements of such cutoff become possible, and two cases (GRB~090926A and GRB~100724B) have been reported to have high-energy cutoffs or breaks. We systematically search for such high energy spectral cutoffs/breaks from the LAT and the Gamma-ray burst monitor (GBM) observations of the prompt emission of GRBs detected since August 2011. Six more GRBs are found to have cutoff-like spectral feature at energies of $\sim10-500$~MeV. Assuming that these cutoffs are caused by pair-production absorption within the source, the bulk Lorentz factors of these GRBs are obtained. We further find that the Lorentz factors are correlated with the isotropic gamma-ray luminosity of the bursts, indicating that more powerful GRB jets move faster.

preprint2015arXiv

Neutrinos from gamma-ray bursts: propagation of cosmic rays in their host galaxies

Gamma-ray bursts (GRBs) are proposed as candidate sources of ultra-high energy cosmic rays (UHECRs). We study the possibility that the PeV neutrinos recently observed by IceCube are produced by GRB cosmic rays interacting with the interstellar gas in the host galaxies. By studying the relation between the X-ray absorption column density N_H and the surface star-formation rate of GRB host galaxies, we find that N_H is a good indicator of the surface gas density of the host galaxies. Then we are able to calculate the neutrino production efficiency of CRs for GRBs with known N_H. We collect a sample of GRBs that have both measurements of N_H and accurate gamma-ray fluence, and attempt to calculate the accumulated neutrino flux based on the current knowledge about GRBs and their host galaxies. When the CR intensity produced by GRBs is normalized with the observed UHECR flux above $10^{19}{\rm eV}$, the accumulated neutrino flux at PeV energies is estimated to be about $(0.3\pm0.2)\times10^{-8} \rm{GeV\ cm^{-2}\ s^{-1}\ sr^{-1}} $ (per flavor) under the assumption that GRB energy production rate follows the cosmic star-formation rate and the favorable assumption about the CR diffusion coefficient. This flux is insufficient to account for the IceCube observations, but the estimate suffers from some assumptions in the calculation and thus we can not rule out this scenario at present.

preprint2015arXiv

Star-forming galaxies as the origin of the IceCube PeV neutrinos

Star-forming galaxies, due to their high star-formation rates and hence large number of supernova remnants therein, are huge reservoirs of cosmic rays (CRs). These CRs collide with gases in the galaxies and produce high-energy neutrinos through $pp$ collisions. In this paper, we calculate the neutrino production efficiency in star-forming galaxies by considering realistic galaxy properties, such as the gas density and galactic wind in star-forming galaxies. To calculate the accumulated neutrino flux, we use the infrared luminosity function of star-forming galaxies obtained by {\em Herschel} PEP/HerMES survey recently. The intensity of CRs producing PeV neutrinos in star-forming galaxies is normalized with the observed CR flux at EeV ({1\,EeV=$10^{18}\,$eV}), assuming that supernova remnants or hypernova remnants in star-forming galaxies can accelerate protons to EeV energies. Our calculations show that the accumulated neutrino emission produced by CRs in star-forming galaxies can account for the flux and spectrum of the sub-PeV/PeV neutrinos under reasonable assumptions on the CR confinement time in these galaxies.

preprint2015arXiv

Tidal disruption jets of supermassive black holes as hidden sources of cosmic rays: explaining the IceCube TeV-PeV neutrinos

Cosmic ray interactions that produce high-energy neutrinos also inevitably generate high-energy gamma rays, which finally contribute to the diffuse high-energy gamma-ray background after they escape the sources. It was recently found that, the high flux of neutrinos at $\sim30$ TeV detected by IceCube lead to a cumulative gamma-ray flux exceeding the Fermi isotropic gamma-ray background at 10-100 GeV, implying that the neutrinos are produced by hidden sources of cosmic rays, where GeV-TeV gamma-rays are not transparent. Here we suggest that relativistic jets in tidal disruption events (TDEs) of supermassive black holes are such hidden sources. We consider the jet propagation in an extended,optically thick envelope around the black hole, which is resulted from the ejected material during the disruption. While powerful jets can break free from the envelope, less powerful jets would be choked inside the envelope. The jets accelerate cosmic rays through internal shocks or reverse shocks and further produce neutrinos via interaction with the surrounding dense medium or photons. All three TDE jets discovered so far are not detected by Fermi/LAT, suggesting that GeV-TeV gamma-rays are absorbed in these jets. The cumulative neutrino flux from TDE jets can account for the neutrino flux observed by IceCube at PeV energies and may also account for the higher flux at $\sim30$ TeV if less powerful, choked jets are present in the majority of TDEs.

preprint2014arXiv

An inverse Compton origin for the 55 GeV photon in the late afterglow of GRB 130907A

The extended high-energy gamma-ray (>100 MeV) emission which occurs well after the prompt gamma-ray bursts (GRBs) is usually explained as the afterglow synchrotron radiation. Here we report the analysis of the Fermi Large Area Telescope observations of GRB 130907A. A 55 GeV photon compatible with the position of the burst was found at about 5 hours after the prompt phase. The probability that this photon is associated with GRB 130907A is higher than 99.96%. The energy of this photon exceeds the maximum synchrotron photon energy at this time and its occurrence thus challenges the synchrotron mechanism as the origin for the extended high-energy >10 GeV emission. Modeling of the broad-band spectral energy distribution suggests that such high energy photons can be produced by the synchrotron self-Compton emission of the afterglow.

preprint2014arXiv

Diffuse PeV neutrinos from EeV cosmic ray sources: semi-relativistic hypernova remnants in star-forming galaxies

We argue that the excess of sub-PeV/PeV neutrinos recently reported by IceCube could plausibly originate through pion-production processes in the same sources responsible for cosmic rays (CRs) with energy above the second knee around $10^{18}\,$eV. The pion production efficiency for escaping CRs that produce PeV neutrinos is required to be $\gtrsim 0.1$ in such sources. On the basis of current data, we identify semi-relativistic hypernova remants as possible sources that satisfy the requirements. By virtue of their fast ejecta, such objects can accelerate protons to EeV energies, which in turn can interact with the dense surrounding medium during propagation in their host galaxies to produce sufficient high-energy neutrinos via proton--proton ($pp$) collisions. Their accompanying gamma ray flux can remain below the diffuse isotropic gamma ray background observed by the {\it Fermi} Large Area Telescope (LAT). In order to test this scenario and discriminate from alternatives, the density of target protons/nuclei and the residence time of CRs in the interacting region are crucial uncertainties that need to be clarified. As long as the neutrinos and EeV CRs originate from the same source class, detection of $\gtrsim 10\,$PeV neutrinos may be expected within 5-10 years' operation of IceCube. Together with further observations in the PeV range, the neutrinos can help in revealing the currently unknown sources of EeV CRs.

preprint2014arXiv

Discovery of GeV emission from the direction of the luminous infrared galaxy NGC 2146

Recent detection of high-energy gamma-ray emission from starburst galaxies M82 and NGC 253 suggests that starburst galaxies are huge reservoirs of cosmic rays and these cosmic rays convert a significant fraction of their energy into gamma-rays by colliding with the dense interstellar medium. In this paper, we report the search for high-energy gamma-ray emission from several nearby star-forming and starburst galaxies using the 68 month data obtained with the Fermi Large Area Telescope. We found a $\sim5.5σ$ detection of gamma-ray emission above 200{\rm MeV} from a source spatially coincident with the location of the luminous infrared galaxy NGC~2146. Taking into account also the temporal and spectral properties of the gamma-ray emission, we suggest that the gamma-ray source is likely to be the counterpart of NGC~2146. The gamma-ray luminosity suggests that cosmic rays in NGC~2146 convert most of their energy into secondary pions, so NGC~2146 is a "proton calorimeter". It is also found that NGC~2146 obeys the quasi-linear scaling relation between the gamma-ray luminosity and total infrared luminosity for star-forming galaxies, strengthening the connection between massive star formation and gamma-ray emission of star-forming galaxies. Possible TeV emission from NGC~2146 is predicted and the implications for high-energy neutrino emission from starburst galaxies are discussed.

preprint2014arXiv

Photosphere emission in the X-Ray Flares of Swift Gamma-Ray Bursts and Implications for the Fireball Properties

X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitted with the absorbed single power-law or the Band function models. More interestingly, the joint spectra of the other 13 X-ray flares are fitted with the absorbed single power-law model plus a black body (BB) component. Phenomenally, the observed spectra of these 13 flares are analogous to several GRBs with a thermal component, but only with a much lower temperature of $kT=1\sim 3$ keV. Assuming that the thermal emission is the photosphere emission of the GRB fireball, we derive the fireball properties of the 13 flares that have redshift measurements, such as the bulk Lorentz factor $Γ_{\rm ph}$ of the outflow. The derived $Γ_{\rm ph}$ range from $50$ to $150$ and a relation of $Γ_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $Γ_0-L_{\rm iso}$ relation that are derived for the prompt gamma-ray emission. We discuss the physical implications of these results within the content of jet composition and radiation mechanism of GRBs and X-ray flares.

preprint2014arXiv

The diffuse gamma-ray flux associated with sub-PeV/PeV neutrinos from starburst galaxies

One attractive scenario for the excess of sub-PeV/PeV neutrinos recently reported by IceCube is that they are produced by cosmic rays in starburst galaxies colliding with the dense interstellar medium. These proton-proton ($pp$) collisions also produce high-energy gamma-rays, which finally contribute to the diffuse high-energy gamma-ray background. We calculate the diffuse gamma-ray flux with a semi-analytic approach and consider that the very high energy gamma-rays will be absorbed in the galaxies and converted into electron-position pairs, which then lose almost all their energy through synchrotron radiation in the strong magnetic fields in the starburst region. Since the synchrotron emission goes into energies below GeV, this synchrotron loss reduces the diffuse high-energy gamma-ray flux by a factor of about two, thus leaving more room for other sources to contribute to the gamma-ray background. For a $E_ν^{-2}$ neutrino spectrum, we find that the diffuse gamma-ray flux contributes about 20% of the observed diffuse gamma-ray background in the 100 GeV range. However, for a steeper neutrino spectrum, this synchrotron loss effect is less important, since the energy fraction in absorbed gamma-rays becomes lower.

preprint2013arXiv

Constraints on the source of ultra-high energy cosmic rays using anisotropy vs chemical composition

The joint analysis of anisotropy signals and chemical composition of ultra-high energy cosmic rays offers strong potential for shedding light on the sources of these particles. Following up on an earlier idea, this paper studies the anisotropies produced by protons of energy >E/Z, assuming that anisotropies at energy >E have been produced by nuclei of charge Z, which share the same magnetic rigidity. We calculate the number of secondary protons produced through photodisintegration of the primary heavy nuclei. Making the extreme assumption that the source does not inject any proton, we find that the source(s) responsible for anisotropies such as reported by the Pierre Auger Observatory should lie closer than ~20-30, 80-100 and 180-200 Mpc if the anisotropy signal is mainly composed of oxygen, silicon and iron nuclei respectively. A violation of this constraint would otherwise result in the secondary protons forming a more significant anisotropy signal at lower energies. Even if the source were located closer than this distance, it would require an extraordinary metallicity >120, 1600, 1100 times solar metallicity in the acceleration zone of the source, for oxygen, silicon and iron respectively, to ensure that the concomitantly injected protons not to produce a more significant low energy anisotropy. This offers interesting prospects for constraining the nature and the source of ultra-high energy cosmic rays with the increase in statistics expected from next generation detectors.

preprint2013arXiv

Discovery of an extra hard spectral component in the high-energy afterglow emission of GRB 130427A

The extended high-energy gamma-ray (>100 MeV) emission occurred after the prompt gamma-ray bursts (GRBs) is usually characterized by a single power-law spectrum, which has been explained as the afterglow synchrotron radiation. The afterglow inverse-Compton emission has long been predicted to be able to produce a high-energy component as well, but previous observations have not revealed such a signature clearly, probably due to the small number of >10 GeV photons even for the brightest GRBs known so far. In this Letter, we report on the Fermi Large Area Telescope (LAT) observations of the >100 MeV emission from the very bright and nearby GRB 130427A. We characterize the time-resolved spectra of the GeV emission from the GRB onset to the afterglow phase. By performing time-resolved spectral fits of GRB 130427A, we found a strong evidence of an extra hard spectral component that exists in the extended high-energy emission of this GRB. We argue that this hard component may arise from the afterglow inverse Compton emission.

preprint2013arXiv

Interpretation of the unprecedentedly long-lived high-energy emission of GRB 130427A

High energy photons (>100 MeV) are detected by the Fermi/LAT from GRB 130427A up to almost one day after the burst, with an extra hard spectral component being discovered in the high-energy afterglow. We show that this hard spectral component arises from afterglow synchrotron-self Compton emission. This scenario can explain the origin of >10 GeV photons detected up to ~30000s after the burst, which would be difficult to be explained by synchrotron radiation due to the limited maximum synchrotron photon energy. The lower energy multi-wavelength afterglow data can be fitted simultaneously by the afterglow synchrotron emission. The implication of detecting the SSC emission for the circumburst environment is discussed.

preprint2013arXiv

On the magnetization of gamma-ray burst blast waves

The origin of magnetic fields that permeate the blast waves of gamma-ray bursts (GRBs) is a long-standing problem. The present paper argues that in four GRBs revealing extended emission at >100 MeV, with follow-up in the radio, optical and X-ray domains at later times, this magnetization can be described as the partial decay of the micro-turbulence that is generated in the shock precursor. Assuming that the bulk of the extended emission >100 MeV can be interpreted as synchrotron emission of shock accelerated electrons, we model the multi-wavelength light curves of GRB 090902B, GRB 090323, GRB 090328 and GRB 110731A, using a simplified then a full synchrotron calculation with power-law-decaying microturbulence ε_B \propto t^{α_t} (t denotes the time since injection through the shock, in the comoving blast frame). We find that these models point to a consistent value of the decay exponent -0.5 < α_t < -0.4.

preprint2013arXiv

On the origin of >10 GeV photons in gamma-ray burst afterglows

Fermi/LAT has detected long-lasting high-energy photons (>100 MeV) from gamma-ray bursts (GRBs), with the highest energy photons reaching about 100 GeV. One proposed scenario is that they are produced by high-energy electrons accelerated in GRB forward shocks via synchrotron radiation. We study the maximum synchrotron photon energy in this scenario, considering the properties of the microturbluence magnetic fields behind the shock, as revealed by recent Particle-in-Cell simulations and theoretical analyses of relativistic collisionless shocks. Due to the small-scale nature of the micro-turbulent magnetic field, the Bohm acceleration approximation breaks down at such high energies. This effect leads to a typical maximum synchrotron photon of a few GeV at 100 s after the burst and this maximum synchrotron photon energy decreases quickly with time. We show that the fast decrease of the maximum synchrotron photon energy leads to a fast decay of the synchrotron flux. The 10-100 GeV photons detected after the prompt phase can not be produced by the synchrotron mechanism. They could originate from the synchrotron self-Compton emission of the early afterglow if the circum-burst density is sufficiently large, or from the external inverse-Compton process in the presence of central X-ray emission, such as X-ray flares and prompt high-latitude X-ray emission.

preprint2012arXiv

IceCube Non-detection of GRBs: Constraints on the Fireball Properties

The increasingly deep limit on the neutrino emission from gamma-ray bursts (GRBs) with IceCube observations has reached the level that could put useful constraints on the fireball properties. We first present a revised analytic calculation of the neutrino flux, which predicts a flux an order of magnitude lower than that obtained by the IceCube collaboration. For benchmark model parameters (e.g. the bulk Lorentz factor is Γ=10^{2.5}, the observed variability time for long GRBs is t_v=0.01 s and the ratio between the energy in accelerated protons and in radiation is η_p=10 for every burst) in the standard internal shock scenario, the predicted neutrino flux from 215 bursts during the period of the 40-string and 59-string configurations is found to be a factor of ~3 below the IceCube sensitivity. However, if we accept the recently found inherent relation between the bulk Lorentz factor and burst energy, the expected neutrino flux increases significantly and the spectral peak shifts to lower energy. In this case, the non-detection then implies that the baryon loading ratio should be η_p<10 if the variability time of long GRBs is fixed to t_v=0.01 s. Instead, if we relax the standard internal shock scenario but keep to assume η_p=10, the non-detection constrains the dissipation radius to be R>4x10^{12} cm assuming the same dissipation radius for every burst and benchmark parameters for fireballs. We also calculate the diffuse neutrino flux from GRBs for different luminosity functions existing in the literature. The expected flux exceeds the current IceCube limit for some luminosity functions, and thus the non-detection constrains η_p<10 in such cases when the variability time of long GRBs is fixed to t_v=0.01 s.

preprint2012arXiv

On the excess of ultra-high energy cosmic rays in the direction of Centaurus A

A posteriori anisotropy study of ultra-high energy cosmic rays (UHECRs) with the Pierre Auger Observatory (PAO) has shown evidence of excess of cosmic ray particles above 55 EeV within $18^{\circ}$ of the direction of the radio galaxy Centaurus A. However, the origin of the excess remains elusive. We simulate the propagation of different species of particles coming from the direction of Centaurus A in the Galactic magnetic fields, and find that only particles of nuclear charge $Z\la 10$ can avoid being deflected outside of the $18^{\circ}$ window of Centaurus A. On the other hand, considering the increasingly heavy composition of UHECRs at the highest energies measured by PAO, a plausible scenario for cosmic rays from the direction of Centaurus A can be found if they consist of intermediate-mass nuclei. The chemical composition of cosmic rays can be further constrained by lower-energy cosmic rays of the same rigidity. We find that cosmic ray acceleration in the lobes of Centaurus A is not favored, while acceleration in the stellar winds that are rich in intermediate-mass nuclei, could meet the requirement. This suggests that the observed excess may originate from cosmic ray accelerators induced by stellar explosions in the star-forming regions of Centaurus A and/or the Centaurus cluster located behind Centaurus A.

preprint2012arXiv

Origin of the GeV Emission During the X-ray Flaring Activity in GRB 100728A

Recently, Fermi-LAT detected GeV emission during the X-ray flaring activity in GRB 100728A. We study various scenarios for its origin. The hard spectrum of the GeV emission favors the external inverse-Compton origin in which X-ray flare photons are up-scattered by relativistic electrons in the external forward shock. This external IC scenario, with anisotropic scattering effect taken into account, can reproduce the temporal and spectral properties of the GeV emission in GRB 100728A.

preprint2012arXiv

Probing the structure of the outflow in the tidal disruption flare Sw J1644+57 with long-term radio emission

The recently discovered high-energy transient Sw J1644+57 is thought to arise from the tidal disruption of a passing star by a dormant massive black hole. The long-term, bright radio emission of Sw J1644+57 is believed to result from the synchrotron emission of the blast wave produced by an outflow expanding into the surrounding medium. Using the detailed multi-epoch radio spectral data, we are able to determine the total number of radiating electrons in the outflow at different times, and further the evolution of the cross section of the outflow with time. We find that the outflow gradually transits from a conical jet to a cylindrical one at later times. The transition may be due to collimation of the outflow by the pressure of the shocked jet cocoon that forms while the outflow is propagating in the ambient medium. Since cylindrical jets usually exist in AGNs and extragalactic jets, this may provide independent evidence that Sw J1644+57 signals the onset of an AGN.

preprint2012arXiv

Unusual Central Engine Activity in the Double Burst GRB 110709B

The double burst, GRB 110709B, triggered Swift/BAT twice at 21:32:39 UT and 21:43:45 UT, respectively, on 9 July 2011. This is the first time we observed a GRB with two BAT triggers. In this paper, we present simultaneous Swift and Konus-WIND observations of this unusual GRB and its afterglow. If the two events originated from the same physical progenitor, their different time-dependent spectral evolution suggests they must belong to different episodes of the central engine, which may be a magnetar-to-BH accretion system.

preprint2011arXiv

Constraining the Emissivity of Ultrahigh Energy Cosmic Rays in the Distant Universe with the Diffuse Gamma-ray Emission

Ultra-high cosmic rays (UHECRs) with energies >10^19 eV emitted at cosmological distances will be attenuated by cosmic microwave and infrared background radiation through photohadronic processes. Lower energy extra-galactic cosmic rays (~10^18-10^19 eV) can only travel a linear distance smaller than ~Gpc in a Hubble time due to the diffusion if the extra-galactic magnetic fields are as strong as nano Gauss. These prevent us from directly observing most of the UHECRs in the universe, and thus the observed UHECR intensity reflects only the emissivity in the nearby universe within hundreds of Mpc. However, UHECRs in the distant universe, through interactions with the cosmic background photons, produce UHE electrons and gamma-rays that in turn initiate electromagnetic cascades on cosmic background photons. This secondary cascade radiation forms part of the extragalactic diffuse GeV-TeV gamma-ray radiation and, unlike the original UHECRs, is observable. Motivated by new measurements of extragalactic diffuse gamma-ray background radiation by Fermi/LAT, we obtained upper limits placed on the UHECR emissivity in the distant universe by requiring that the cascade radiation they produce not exceed the observed levels. By comparison with the gamma-ray emissivity of candidate UHECR sources (such as GRBs and AGNs) at high-redshifts, we find that the obtained upper limit for a flat proton spectrum is ~10^1.5 times larger than the gamma-ray emissivity in GRBs and ~10 times smaller than the gamma-ray emissivity in BL Lac objects. In the case of iron nuclei composition, the derived upper limit of the UHECR emissivity is a factor of 3-5 times higher. Robust upper limit on the cosmogenic neutrino flux is further obtained, which is marginally reachable by the Icecube detector and the next-generation detector JEM-EUSO.

preprint2011arXiv

Nearby low-luminosity GRBs as the sources of ultra-high energy cosmic rays revisited

Low-luminosity gamma-ray bursts (GRBs) with luminosity . 10^49erg/s probably consititute a distinct population from the classic high-luminosity GRBs. They are the most luminous objects detected so far within ~ 100 Mpc, the horizon distance of ultra-high energy cosmic rays (UHECRs), so they are considered to be candidate sources of UHECRs. It was recently argued that the energy production rate in UHECRs is much larger than that in gamma-ray photons of long GRBs measured by the Fermi satellite, which, if true, would challenge the view that GRBs can be the sources of UHECRs. We here suggest that many of the low luminosity GRBs, due to their low luminosity, can not trigger the current GRB detectors and hence their contribution to the local gamma-ray energy production rate is missing. We find that the real local energy production rate by low-luminosity GRBs, taking into account the missing part, which constitutes a dominant fraction of the total amount, could be sufficient to account for the flux of UHECRs. Due to the low-luminosity, only intermediate-mass or heavy nuclei can be accelerated to ~ 10^20 eV. We discuss the acceleration and survival of these UHE nuclei in low-luminosity GRBs, especially in those missing low-luminosity GRBs. At last, the accompanying diffuse neutrino flux from the whole low-luminosity GRB population is calculated.

preprint2011arXiv

On the High Energy Emission of the Short GRB 090510

Long-lived high-energy (>100MeV) emission, a common feature of most Fermi-LAT detected gamma-ray burst, is detected up to \sim 10^2 s in the short GRB 090510. We study the origin of this long-lived high-energy emission, using broad-band observations including X-ray and optical data. We confirm that the late > 100 MeV, X-ray and optical emission can be naturally explained via synchrotron emission from an adiabatic forward shock propagating into a homogeneous ambient medium with low number density. The Klein-Nishina effects are found to be significant, and effects due to jet spreading and magnetic field amplification in the shock appear to be required. Under the constraints from the low-energy observations, the adiabatic forward shock synchrotron emission is consistent with the later-time (t>2s) high-energy emission, but falls below the early-time (t < 2s) high energy emission. Thus we argue that an extra high energy component is needed at early times. A standard reverse shock origin is found to be inconsistent with this extra component. Therefore, we attribute the early part of the high-energy emission (t< 2s) to the prompt component, and the long-lived high energy emission (t>2s) to the adiabatic forward shock synchrotron afterglow radiation. This avoids the requirement for an extremely high initial Lorentz factor.

preprint2011arXiv

Probing the tidal disruption flares of massive black holes with high-energy neutrinos

The recently discovered high-energy transient Swift J164449.3+573451 (Sw J1644+57) is thought to arise from the tidal disruption of a passing star by a dormant massive black hole. Modeling of the broadband emission suggests the presence of a powerful relativistic jet, which contributes dominantly to the observed X-ray emission. Here we suggest that protons can be accelerated to ultra-high energies by internal shocks occurring in the jets, but their flux is insufficient to account for the observed flux of ultra-high energy cosmic rays. High energy protons can produce ~0.1-10 PeV neutrinos through photomeson interactions with X-ray photons. The large X-ray fluence (7x10^-4 erg cm^-2) and high photopion efficiency, together with the insignificant cooling of secondary mesons, result in bright neutrino emission expected from Sw J1644+57 if the jet composition is matter-dominated. One to several neutrinos may be detected by a Km^3-scale detector from one tidal disruption event similar to Sw J1644+57, thereby providing a powerful probe of the composition of the jets.

preprint2011arXiv

Spectrum and Composition of Ultra-high Energy Cosmic Rays from Semi-relativistic Hypernovae

It has been suggested that hypernova remnants, with a substantial amount of energy in semi-relativistic ejecta, can accelerate intermediate mass or heavy nuclei to ultra-high energies and provide sufficient amount of energy in cosmic rays to account for the observed flux. We here calculate the expected energy spectrum and chemical composition of ultra-high energy cosmic rays from such semi-relativistic hypernovae. With a chemical composition equal to that of the hypernova ejecta and a flat or hard spectrum for cosmic rays at the sources, the spectrum and composition of the propagated cosmic rays observed at the Earth can be compatible with the measurements by the Pierre Auger Observatory.

preprint2010arXiv

Late-time detections of the X-ray afterglow of GRB 060729 with Chandra - the latest detections ever of an X-ray afterglow

We report on 5 Chandra observations of the X-ray afterglow of the Gamma-Ray Burst GRB 060729 performed between 2007 March and 2008 May. In all five observations the afterglow is clearly detected. The last Chandra pointing was performed on 2008-May-04, 642 days after the burst - the latest detection of a GRB X-ray afterglow ever. A reanalysis of the Swift XRT light curve together with the three detections by Chandra in 2007 reveals a break at about 1.0 Ms after the burst with a slight steepening of the decay slope from alpha = 1.32 to 1.61. This break coincides with a significant hardening of the X-ray spectrum, consistent with a cooling break in the wind medium scenario, in which the cooling frequency of the afterglow crosses the X-ray band. The last two Chandra observations in 2007 December and 2008 May provide evidence for another break at about one year after the burst. If interpreted as a jet break, this late-time break implies a jet half opening angle of about 14 degrees for a wind medium. Alternatively, this final break may have a spectral origin, in which case no jet break has been observed and the half-opening angle of the jet of GRB 060729 must be larger than about 15 degrees for a wind medium. We compare the X-ray afterglow of GRB 060729 in a wind environment with other bright X-ray afterglows, in particular GRBs 061121 and 080319B, and discuss why the X-ray afterglow of GRB 060729 is such an exceptionally long-lasting event.

preprint2010arXiv

Modeling the broadband emission of Fermi/LAT GRB 090902B

GRB 090902B, detected by Fermi Large Array Telescope (Fermi/LAT), shows extend high-energy emission (>100 MeV) up to 10^3 s after the burst, which decays with time in a power-law as t^{-1.5}. It has been also observed by several follow-up low-energy instruments, including an early optical detection around 5000 s after the burst. The optical emission at early time decays faster than t^{-1.6}, which has been suspected to originate from the reverse shock. We here explore the models that can possibly explain the the broadband afterglow emission of GRB 090902B. We find that the reverse shock model for the early optical emission would overpredict the radio afterglow flux that is inconsistent with observations. A partially radiative blast wave model, which though is able to produce a sufficiently steep decay slope, can not explain the broadband data of GRB 090902B. The two-component jet model, which consists of a narrow and bright jet component in the core and a surrounding wider and less energetic jet component, is shown to be able to explain the broadband afterglow data, including the LAT high-energy data after ~50 s and low-energy (radio, optical and X-ray) afterglow data. The early-time high-energy emission detected by LAT before ~50 s is likely due to internal origin as that of the sub-MeV emission. The highest energy (33 GeV) photon of GRB090902B detected at 80 s can be marginally accommodated within the forward shock emission under the optimistic condition that electrons are accelerated by the Bohm diffusive shock.

preprint2009arXiv

GRB 080916C: on the radiation origin of the prompt emission from KeV/MeV to GeV

Fermi observations of high-energy gamma-ray emission from GRB 080916C shows that its spectrum is consistent with the Band function from MeV to tens of GeV. Assuming one single emission mechanism dominates in the whole energy range, we show that this spectrum is consistent with synchrotron origin by shock-accelerated electrons. The simple electron inverse-Compton model and the hadronic model are found to be less viable. In the synchrotron scenario, the synchrotron self-Compton scattering is likely to be in the Klein-Nishina regime and therefore the resulting high-energy emission is subdominant, even though the magnetic field energy density is lower than that in relativistic electrons. The Klein-Nishina inverse-Compton cooling may also affect the low-energy electron number distribution and hence results in a low-energy synchrotron photon spectrum $n(ν)\proptoν^{-1}$ below the peak energy. Under the framework of the electron synchrotron interpretation, we constrain the shock microphysical parameters and derive a lower limit of the upstream magnetic fields. The detection of synchrotron emission extending to about 70 GeV in the source frame in GRB 080916C favors the Bohm diffusive shock acceleration if the bulk Lorentz factor of the relativistic outflow is not significantly greater than thousands.

preprint2009arXiv

High-energy gamma-ray afterglows from low-luminosity gamma-ray bursts

The observations of gamma-ray bursts (GRBs) such as 980425, 031203 and 060218, with luminosities much lower than those of other classic bursts, lead to the definition of a new class of GRBs -- low-luminosity GRBs. The nature of the outflow responsible for them is not clear yet. Two scenarios have been suggested: one is the conventional relativistic outflow with initial Lorentz factor of order of $Γ_0\ga 10$ and the other is a trans-relativistic outflow with $Γ_0\simeq 1-2$. Here we compare the high energy gamma-ray afterglow emission from these two different models, taking into account both synchrotron self inverse-Compton scattering (SSC) and the external inverse-Compton scattering due to photons from the cooling supernova or hypernova envelope (SNIC). We find that the conventional relativistic outflow model predicts a relatively high gamma-ray flux from SSC at early times ($<10^4 {\rm s}$ for typical parameters) with a rapidly decaying light curve, while in the trans-relativistic outflow model, one would expect a much flatter light curve of high-energy gamma-ray emission at early times, which could be dominated by both the SSC emission and SNIC emission, depending on the properties of the underlying supernova and the shock parameter $ε_e$ and $ε_B$. The Fermi Gamma-ray Space Telescope should be able to distinguish between the two models in the future.

preprint2009arXiv

Klein-Nishina effects on the high-energy afterglow emission of gamma-ray bursts

Extended high-energy(>100MeV) gamma-ray emission that lasts much longer than the prompt sub-MeV emission has been detected from quite a few gamma-ray bursts (GRBs) by Fermi Large Area Telescope (LAT) recently. A plausible scenario is that this emission is the afterglow synchrotron emission produced by electrons accelerated in the forward shocks. In this scenario, the electrons that produce synchrotron high-energy emission also undergo inverse-Compton (IC) loss and the IC scattering with the synchrotron photons should be in the Klein-Nishina regime. Here we study effects of the Klein-Nishina scattering on the high-energy synchrotron afterglow emission. We find that, at early times the Klein-Nishina suppression effect on those electrons that produce the high-energy emission is usually strong and therefore their inverse-Compton loss is small with a Compton parameter Y < a few for a wide range of parameter space. This leads to a relatively bright synchrotron afterglow at high energies that can be detected by Fermi LAT. As the Klein-Nishina suppression effect weakens with time, the inverse-Compton loss increases and could dominate over the synchrotron loss in some parameter space. This will lead to a faster temporal decay of the high-energy synchrotron emission than what is predicted by the standard synchrotron model, which may explain the observed rapid decay of the early high-energy gamma-ray emission in GRB090510 and GRB090902B.

preprint2009arXiv

The early high-energy afterglow emission from Short GRBs

We calculate the high energy afterglow emission from short Gamma-Ray Bursts (SGRBs) in the external shock model. There are two possible components contributing to the high energy afterglow: the electron synchrotron emission and the synchrotron self-Compton (SSC) emission. We find that for typical parameter values of SGRBs, the early high-energy afterglow emission in 10 MeV-10 GeV is dominated by the synchrotron emission. For a burst occurring at redshift z =0.1, the high-energy emission can be detectable by Fermi LAT if the blast wave has an energy E>=10^51 ergs and the fraction of energy in electrons is ε_e>=0.1 . This provides a possible explanation for the high energy tail of SGRB 081024B.

preprint2008arXiv

On the Origin and Survival of UHE Cosmic-Ray Nuclei in GRBs and Hypernovae

The chemical composition of the ultra-high-energy (UHE) cosmic rays serves as an important clue for their origin. Recent measurements of the elongation rates by the Pierre Auger Observatory hint at the possible presence of heavy or intermediate mass nuclei in the UHE cosmic rays. Gamma-ray bursts (GRBs) and hypernovae have been suggested as possible sources of the UHE cosmic rays. Here we derive the constraints on the physical conditions under which UHE heavy nuclei, if they are accelerated in these sources, can survive in their intense photon fields. We find that in the GRB external shock and in the hypernova scenarios, UHE nuclei can easily survive photo-disintegration. In the GRB internal shock scenario, UHE nuclei can also survive, provided the dissipation radius and/or the bulk Lorentz factor of the relativistic outflow are relatively large, or if the low-energy self-absorption break in the photon spectrum of the prompt emission occurs above several KeV. In internal shocks and in the other scenarios, intermediate-mass UHE nuclei have a higher probability of survival against photo-disintegration than UHE heavy nuclei such as Fe.

preprint2007arXiv

Swift and XMM-Newton Observations of the Extraordinary GRB 060729: An afterglow with a more than 100 days X-ray light curve

We report the results of the Swift and XMM observations of the Swift-discovered long Gamma-Ray Burst GRB 060729 ($T_{90}$=115s). The afterglow of this burst was exceptionally bright in X-rays as well as at UV/Optical wavelengths showing an unusually long slow decay phase ($α$=0.14\plm0.02) suggesting a larger energy injection phase at early times than in other bursts. The X-ray light curve displays a break at about 60 ks after the burst. The X-ray decay slope after the break is $α$=1.29\plm0.03. Up to 125 days after the burst we do not detect a jet break, suggesting that the jet opening angle is larger than 28 degrees. In the first 2 minutes after the burst (rest frame) the X-ray spectrum of the burst changed dramatically from a hard X-ray spectrum to a very soft one. We find that the X-ray spectra at this early phase can all be fitted by an absorbed single power law model or alternatively by a blackbody plus power law model. The power law fits show that the X-ray spectrum becomes steeper while the absorption column density decreases. In Swift's UV/Optical telescope the afterglow was clearly detected up to 9 days after the burst in all 6 filters and even longer in some of the UV filters with the latest detection in the UVW1 31 days after the burst. A break at about 50 ks is clearly detected in all 6 UVOT filters from a shallow decay slope of about 0.3 and a steeper decay slope of 1.3. In addition to the \swift observations we also present and discuss the results from a 61 ks ToO observation by XMM. (Abriviated)