Source author record

Ruo-Yu Liu

Ruo-Yu Liu 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

42works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

42 published item(s)

preprint2025arXiv

Pulsar+escape: the recipe for LHAASO sources ?

Observations of our Galaxy at very-high and ultra-high photon energies have revealed a rich population of sources, many of which have significant angular extension and/or are positionally coincident with energetic pulsars. We assessed whether the properties of these gamma-ray emitters are consistent with the predictions of a model for the non-thermal leptonic emission of a pulsar wind nebula expanding within a supernova remnant, including the possibility of particle escape across the different components of the system up to the ambient medium. We used a multi-zone model framework that describes the transport and radiation of a population of ultra-relativistic electron-positron pairs. The dynamics of the system is retrieved from numerical hydrodynamics simulations, which allows to follow the evolution through advanced stages. Model predictions for typical parameter setups are compared to source catalogs issued by the H.E.S.S. and LHAASO collaborations. If particles remain confined in the nebula, 1-100TeV emission is only partially consistent with the properties of observed sources in terms of surface brightness, angular extent, and photon index. In particular, the systems spend a long time in a bright and soft state that has no equivalent in observations. Conversely, including the possibility of energy-dependent particle escape across the object results in a much better agreement, both in the 1-10TeV range explored with H.E.S.S. or LHAASO-WCDA, and in the 20-100TeV range explored by LHAASO-KM2A. Emission components with intermediate to low brightnesses and large to very large sizes result from particles trapped in the remnant or escaped in the surrounding medium. Particle escape, clearly seen in the form of misaligned jets or halos around middle-aged pulsars, is a very important process in much younger systems with 1-100kyr ages, and shapes their appearance as TeV/PeV sources.

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.

preprint2022arXiv

A PeVatron Candidate: Modelling the Boomerang Nebula in X-ray Band

Pulsar wind nebula (PWN) Boomerang and the associated supernova remnant (SNR) G106.3+2.7 are among candidates for the ultra-high-energy (UHE) gamma-ray counterparts published by LHAASO. Although the centroid of the extended source, LHAASO J2226+6057, deviates from the pulsar's position by about $0.3^\circ$, the source partially covers the PWN. Therefore, we cannot totally exclude the possibility that a part of the UHE emission comes from the PWN. Previous studies mainly focus on whether the SNR is a PeVatron, while neglecting the energetic PWN. Here, we explore the possibility of the Boomerang Nebula being a PeVatron candidate by studying its X-ray radiation. By modelling the diffusion of relativistic electrons injected in the PWN, we fit the radial profiles of the X-ray surface brightness and the photon index. The solution with a magnetic field $B=140μ$G can well reproduce the observed profiles and implies a severe suppression of IC scattering of electrons. Therefore, a proton component need be introduced to account for the UHE emission, in light of recent LHAASO's measurement on Crab Nebula, if future observations reveal part of the UHE emission originating from the PWN. In this sense, Boomerang Nebula would be a hadronic PeVatron.

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

On the Hadronic Origin of High Energy Emission of $γ$-ray Loud Narrow-Line Seyfert 1 PKS 1502+036

The radiation mechanism of Radio-Loud Narrow-Line Seyfert 1 (RL-NLS1) from X-ray to $γ$-ray bands remains an open question. While the leptonic model has been employed to explain the spectral energy distribution (SED), the hadronic process may potentially account for the high energy radiation of some $γ$-ray loud Narrow-Line Seyfert 1 (NLS1) as well. We study one of such RL-NLS1, PKS 1502+036, comparing the theoretical SEDs predicted by the leptonic model and the lepto-hadronic model to the observed one. For the hadronic processes, we take into account the proton synchrotron radiation and proton-photon interactions (including the Bethe-Heitler process and the photopion process) including the emission of pairs generated in the electromagnetic cascade initiated by these processes. Our results show that the leptonic model can reproduce the SED of this source, in which the X-ray to $γ$-ray radiation can be interpreted as the inverse Compton (IC) scattering. On the other hand, the proton synchrotron radiation can also explain the high energy component of SED although extreme parameters are needed. We also demonstrate that the $pγ$ interactions as well as the cascade process cannot explain SED. Our results imply that a leptonic origin is favored for the multi-wavelength emission of PKS 1502+036.

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

Chandra and XMM-Newton observations of A2256: cold fronts, merger shocks, and constraint on the IC emission

We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic (RR). The temperature and metallicity maps show clear evidence of a merger between the western subcluster (SC) and the primary cluster (PC). We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1 is associated with the infalling SC; CF2 is located in the east of the PC; and CF3 is to the west of the PC core. The other two edges at cluster outskirts are shock fronts (SFs): SF1 near the RR in the NW has Mach numbers derived from the temperature and the density jumps, respectively, of $M_T=1.62\pm0.12$ and $M_ρ=1.23\pm0.06$; SF2 in the SE has $M_T=1.54\pm0.05$ and $M_ρ=1.16\pm0.13$. In the region of the RR, there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse-Compton emission, and therefore set a lower limit on the magnetic field ($\sim$ 450 kpc from PC center) of $B>1.0\ μ$G for a single power-law electron spectrum or $B>0.4\ μ$G for a broken power-law electron spectrum. We propose a merger scenario including a PC, an SC, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects.

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

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

Gamma-ray heartbeat powered by the microquasar SS 433

Microquasars, the local siblings of extragalactic quasars, are binary systems comprising a compact object and a companion star. By accreting matter from their companions, microquasars launch powerful winds and jets, influencing the interstellar environment around them. Steady gamma-ray emission is expected to rise from their central objects, or from interactions between their outflows and the surrounding medium. The latter prediction was recently confirmed with the detection of SS 433 at high (TeV) energies. In this report, we analyze more than ten years of GeV gamma-ray data from the Fermi Gamma-ray Space Telescope on this source. Detailed scrutiny of the data reveal emission in the SS 433 vicinity, co-spatial with a gas enhancement, and hints for emission possibly associated with a terminal lobe of one of the jets. Both gamma-ray excesses are relatively far from the central binary, and the former shows evidence for a periodic variation at the precessional period of SS 433, linking it with the microquasar. This result challenges obvious interpretations and is unexpected from any previously published theoretical models. It provides us with a chance to unveil the particle transport from SS 433 and to probe the structure of the local magnetic field in its vicinity.

preprint2020arXiv

On the Unusually Large Spatial Extent of the TeV nebula HESS J1825-137: Implication from the Energy-Dependent Morphology

Deep observation of the High Energy Stereoscopic System (HESS) on the most extended pulsar wind nebula HESS J1825-137 reveals an enhanced energy-dependent morphology, providing useful information on the particle transport mechanism in the nebula. We find that the energy-dependent morphology is consistent with a diffusion-dominated transport of electrons/positrons. It provides an alternative possible interpretation for the unusually large spatial extent (i.e., $\gtrsim 100$pc) of the nebula, which could then be attributed to the diffusion of escaping electrons/positrons from a compact plerion. The influence of various model parameters on the energy-dependent extent of the nebula is studied in the diffusion-dominated scenario. We also show that the energy-dependent morphology of the nebula may also be used to study the spin-down history of the pulsar.

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

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

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

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.

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 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

Spin Evolution of Millisecond Magnetars with Hyperaccreting Fallback Disks: Implications for Early Afterglows of Gamma-Ray Bursts

The shallow decay phase or plateau phase of early afterglows of gamma-ray bursts (GRBs), discovered by Swift, is currently understood as being due to energy injection to a relativistic blast wave. One natural scenario for energy injection invokes a millisecond magnetar as the central engine of GRBs, because the conventional model of a pulsar predicts a nearly constant magnetic-dipole-radiation luminosity within the spin-down timescale. However, we note that significant brightening occurs in some early afterglows, which apparently conflicts with the above scenario. Here we propose a new model to explain this significant brightening phenomena by considering a hyperaccreting fallback disk around a newborn millisecond magnetar. We show that for typical values of the model parameters, sufficient angular momentum of the accreted matter is transferred to the magnetar and spins it up. It is this spin-up that leads to a dramatic increase of the magnetic dipole radiation luminosity with time and thus significant brightening of an early afterglow. Based on this model, we carry out numerical calculations and fit well early afterglows of 12 GRBs assuming sufficiently strong fallback accretion. If the accretion is very weak, our model turns out to be the conventional energy-injection scenario of a pulsar. Therefore, our model can provide a unified explanation for the shallow decay phase, plateaus, and significant brightening of early afterglows.

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

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

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.