Researcher profile

Ruo-Yu Liu

Ruo-Yu Liu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 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&#39;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&#39;s line-of-sight (LOS) to the pulsar and the local mean field direction is small. In general, the TeV halo&#39;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&#39;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&#39;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 &#34;orphan&#34; 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 &#34;TeV halos&#34;) 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&#39;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 &#34;orphan&#34; 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 &#34;TeV Halos&#34;, 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.