Researcher profile

Xiang-Yu Wang

Xiang-Yu Wang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
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

19 published item(s)

preprint2025arXiv

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

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

preprint2024arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

Isotopic and elemental abundances seen in Galactic cosmic rays imply that $\sim20\%$ of the cosmic-ray (CR) nuclei are probably synthesized by massive Wolf-Rayet (WR) stars. Massive star clusters hosting WR and OB-type stars have been proposed as potential Galactic cosmic-ray accelerators for decades, in particular via diffusive shock acceleration at wind termination shocks. Here we report the analysis of {\em Fermi} Large Area Telescope&#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

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

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2019arXiv

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

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

preprint2019arXiv

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

Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named &#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.