Source author record

P. H. T. Tam

P. H. T. Tam 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

49works
3topics
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

49 published item(s)

preprint2026arXiv

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Identifying Galactic PeVatrons (PeV particle accelerators) from the ultra-high-energy (UHE, >100 TeV) $γ$-ray sources plays a crucial role in revealing the origin of Galactic cosmic rays. The UHE source 1LHAASO J1857+0203u is suggested to be associated with HESS J1858+020, which may be attributed to the possible PeVatron candidate supernova remnant (SNR) G35.6$-$0.4 or HII region G35.6$-$0.5. We perform detailed analysis on the very-high-energy and UHE $γ$-ray emissions towards this region with data from the Large High Altitude Air Shower Observatory (LHAASO). 1LHAASO J1857+0203u is detected with a significance of 11.6$σ$ above 100 TeV, indicating the presence of a PeVatron. It has an extension of $\sim 0.18^\circ$ with a power-law (PL) spectral index of $\sim$2.5 in 1-25 TeV and a point-like emission with a PL spectral index of $\sim$3.2 above 25 TeV. Using the archival CO and HI data, we identify some molecular and atomic clouds that may be associated with the TeV $γ$-ray emissions. Our modelling indicates that the TeV $γ$-ray emissions are unlikely to arise from the clouds illuminated by the protons that escaped from SNR G35.6$-$0.4. In the scenario that HII region G35.6$-$0.5 could accelerate particles to the UHE band, the observed GeV-TeV $γ$-ray emission could be well explained by a hadronic model with a PL spectral index of $\sim$2.0 and cutoff energy of $\sim$450 TeV. However, an evolved pulsar wind nebula origin cannot be ruled out.

preprint2026arXiv

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

preprint2026arXiv

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $γ$-rays detectable by arrays such as LHAASO. These $γ$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $pγ$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $γ$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $γ$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $γ$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

preprint2026arXiv

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

preprint2026arXiv

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude $b\approx 10.5^\circ$. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21$σ$ and 17$σ$ at 8$-$100 TeV and $>$100 TeV, respectively. The corresponding extensions are determined to be 0.23$^{\circ}\pm$0.03$^{\circ}$ and 0.17$^{\circ}\pm$0.03$^{\circ}$. The emission is proposed to originate from the relativistic electrons and positrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function $dN/dE = (42.4\pm4.1)(\frac{E}{20\rm\ TeV})^{-2.31\pm0.11}\exp(-\frac{E}{110\pm25\rm\ TeV})$ $\rm\ TeV^{-1}\ cm^{-2}\ s^{-1}$in the energy range from 8 TeV to 300 TeV, implying a steady-state parent electron spectrum $dN_e/dE_e\propto (\frac{E_e}{100\rm\ TeV})^{-3.13\pm0.16}\exp[(\frac{-E_e}{373\pm70\rm\ TeV})^2]$ at energies above $\approx 50 \rm\ TeV$. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to $\approx 4.5\rm\ μG$. To satisfy the multi-wavelength spectrum and the $γ$-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

preprint2026arXiv

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

We report the discovery of a new unidentified extended $γ$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $α=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($α\leq 1.75$) spectrum with a cutoff at $\approx 50$ TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

preprint2026arXiv

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

preprint2026arXiv

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range.

preprint2026arXiv

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories.

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

preprint2026arXiv

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Recently the LHAASO Collaboration published the detection of 12 ultra-high-energy gamma-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV gamma rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy gamma rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^5 times the Planck scale M_{pl} and that the second-order LIV scale is >10^{-3}M_{pl}. Both limits improve by at least one order of magnitude the previous results.

preprint2026arXiv

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of $8.2σ$ for $E>25$~TeV assuming a Gaussian template. The best-fit location is (R.A., Dec.)$=(95^{\circ}\!.47\pm0^{\circ}\!.11,\,37^{\circ}\!.92 \pm0^{\circ}\!.09)$, and the extension is $0^{\circ}\!.40\pm0^{\circ}\!.07$. The energy spectrum can be described by a power-law spectrum with an index of ${-2.92 \pm 0.17_{\rm stat} \pm 0.02_{\rm sys} }$. No clear extended multi-wavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of $\sim160$~TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

preprint2026arXiv

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments.

preprint2026arXiv

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

preprint2026arXiv

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

The diffuse Galactic gamma-ray emission is a very important tool used to study the propagation and interaction of cosmic rays in the Milky Way. In this work, we report the measurements of the diffuse emission from the Galactic plane, covering Galactic longitudes from $15^{\circ}$ to $235^{\circ}$ and latitudes from $-5^{\circ}$ to $+5^{\circ}$, in an energy range of 1 TeV to 25 TeV, with the Water Cherenkov Detector Array (WCDA) of the Large High Altitude Air Shower Observatory (LHAASO). After masking the sky regions of known sources, the diffuse emission is detected with $24.6σ$ and $9.1σ$ significance in the inner Galactic plane and outer Galactic plane, respectively. The WCDA spectra in both regions can be well described by a power-law function, with spectral indices of $-2.67\pm0.05_{\rm stat}$ in the inner region and $-2.83\pm0.19_{\rm stat}$ in the outer region, respectively. Combined with the Square Kilometer Array (KM2A) measurements at higher energies, a clear softening of the spectrum is found in the inner region, with change of spectral indices by $\sim0.5$ at a break energy around $30$ TeV. The fluxes of the diffuse emission are higher by a factor of $1.5-2.7$ than the model prediction assuming local CR spectra and the gas column density, which are consistent with those measured by the KM2A. Along Galactic longitude, the spatial distribution of the diffuse emission shows deviation from that of the gas column density. The spectral shape of the diffuse emission are possibly variation in different longitude region. The WCDA measurements bridge the gap between the low-energy measurements by space detectors and the ultra-high-energy observations by LHAASO-KM2A and other experiments. These results suggest that improved modeling of the wide-band diffuse emission is required.

preprint2026arXiv

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties.

preprint2025arXiv

Observation of the $γ$-ray Emission from W43 with LHAASO

In this paper, we report the detection of the very-high-energy (VHE, $ 100{\rm\ GeV} < E < 100{\rm\ TeV} $) and ultra-high-energy (UHE, $E > 100\rm\ TeV$) $γ$-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended $γ$-ray source was detected with a significance of ${\sim}16\,σ$ by KM2A and ${\sim}17\,σ$ by WCDA, respectively. The angular extension of this $γ$-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the $γ$-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

preprint2025arXiv

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4.

preprint2021arXiv

HESS J1858+020: A GeV-TeV source possibly powered by CRs from SNR G35.6-0.4

Context: The supernova remnant (SNR) G35.6-0.4 shows a non-thermal radio shell, however, no γ-ray or X-ray counterparts have been found for it thus far. One TeV source, HESS J1858+020, was found near the SNR and this source is spatially associated with some clouds at 3.6 kpc. Aims: To attain a better understanding of the origin of HESS J1858+020, we further investigate the association between SNR cosmic rays (CRs) and the clouds through the Fermi-LAT analysis and hadronic modeling. Methods: We performed the Fermi-LAT analysis to explore the GeV emission in and around the SNR. We explored the SNR physics with previously observed multi-wavelength data. We built a hadronic model using runaway CRs of the SNR to explain the GeV-TeV observation. Results: We found a hard GeV source (SrcX2) that is spatially coincident with both HESS J1858+020 and a molecular cloud complex at 3.6 kpc. In addition, a soft GeV source (SrcX1) was found at the northern edge of the SNR. The GeV spectrum of SrcX2 connects well with the TeV spectrum of HESS J1858+020. The entire γ-ray spectrum ranges from several GeV up to tens of TeV and it follows a power-law with an index of ~2.15. We discuss several pieces of observational evidence to support the middle-aged SNR argument. Using runaway CRs from the SNR, our hadronic model explains the GeV-TeV emission at HESS J1858+020, with a diffusion coefficient that is much lower than the Galactic value.

preprint2020arXiv

Multi-wavelength observations of the BL Lac object Fermi J1544-0649: one year after its awakening

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($Γ_γ\sim$1.7) in the GeV band and at times also very hard (($Γ_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

preprint2020arXiv

Peculiar prompt emission and afterglow in H.E.S.S. detected GRB 190829A

We present the results of a detailed investigation of the prompt and afterglow emission in the HESS detected GRB 190829A. Swift and Fermi observations of the prompt phase of this GRB reveal two isolated sub-bursts or episodes, separated by a quiescent phase. The energetic and the spectral properties of the first episode are in stark contrast to the second. The first episode, which has a higher spectral peak of $\sim 120\:\text{keV}$ and a low isotropic energy $\sim 10^{50}\:\text{erg}$ is an outlier to the Amati correlation and marginally satisfies the Yonetoku correlation. However, the energetically dominant second episode has lower peak energy and is consistent with the above correlations. We compared this GRB to other low luminosity GRBs (LLGRBs). Prompt emission of LLGRBs also indicates a relativistic shock breakout origin of the radiation. For GRB 190829A, some of the properties of a shock breakout origin are satisfied. However, the absence of an accompanying thermal component and energy above the shock breakout critical limit precludes a shock breakout origin. In the afterglow, an unusual long-lasting late time flare of duration $\sim 10^4\:\text{s}$ is observed. We also analyzed the late-time \fermi-LAT emission that encapsulates the H.E.S.S. detection. Some of the LAT photons are likely to be associated with the source. All the above observational facts suggest GRB 190829A is a peculiar low luminosity GRB that is not powered by a shock breakout, and with an unusual rebrightening due to a patchy emission or a refreshed shock during the afterglow. Furthermore, our results show that TeV energy photons seem common in both high luminosity GRBs and LLGRBs.

preprint2020arXiv

Periodicity search for Pulsar Binaries with TESS

Pulsar binaries, in particular redback systems, provide good sources to study the pulsar wind flow and its interaction with the companion stars. {\it Fermi}-LAT have proposed probable pulsar binary candidates in its catalogs. To identify pulsar binary sources from the catalog, orbital modulation search of binary candidates is an effective way. TESS observes in survey mode for a large part of the sky and thus provide an excellent data set to periodicity search of pulsar binary candidates by observing the flux variation, thought to mainly come from the stellar companion. Using TESS data we look for flux modulation of five pulsar binaries (or candidates) with reported orbital periods, including PSR J1023+0038, 3FGL J0523.3-2528, 3FGL J0212.1+5320, 3FGL J0744.1-2523 and PSR J1417-4402, demonstrating that TESS photometric data are very useful in identifying periodicities of redback-like systems. This method can be effective in searches for new pulsar binaries or similar binary systems in the future.

preprint2020arXiv

Repeated state change of variable gamma-ray pulsar, PSR J2021+4026

PSR J2021+4026 is a radio-quiet gamma-ray pulsar and the first pulsar that shows state change of the gamma-ray emission and spin-down rate. The state change of PSR J2021+4026 was first observed at 2011 October, at which the pulsar changes the state from high gamma-ray flux/low spin-down rate state to low gamma-ray flux/high spin-down rate st\ ate. In December 2014, PSR J2021+4026 recovered the state before the 2011 state change over a timescale of a few months. We report that the long term evolution of the gamma-ray flux and timing behavior suggests that PSR J2021+4026 changed the state near 2018 February 1st and entered a new low gamma-ray flux/high spin-down rate state. At the 2018 state change, the averaged flux dropped from $(1.29\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1}}$ to $(1.12\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1\ }}$, which has the similar behavior to the case of 2011 event. The spin-down rate has increased by $\sim 3\%$ in the new state since the 2018 state change. The shapes of pulse profile and spectrum in GeV bands also changed at the 2018 event, and they are consistent with behavior at the 2011 state change. Our results probably suggest that PSR J2021+4026 is switching between different states with a timescale of several years, like some radio pulsars (e.g. PSR~B1828-11). PSR J2021+4026 will provide a unique opportunity to study the mechanism of the state switching.

preprint2020arXiv

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered.

preprint2016arXiv

Discovery of an X-ray Emitting Contact Binary System 2MASS J11201034$-$2201340

We report the detection of orbital modulation, a model solution, and X-ray properties of a newly discovered contact binary, 2MASS J11201034$-$2201340. We serendipitously found this X-ray point source outside the error ellipse when searching for possible X-ray counterparts of $γ$-ray millisecond pulsars among the unidentified objects detected by the {\it Fermi Gamma-ray Space Telescope}. The optical counterpart of the X-ray source (unrelated to the $γ$-ray source) was then identified using archival databases. The long-term CRTS survey detected a precise signal with a period of $P=0.28876208(56)$ days. A follow-up observation made by the SLT telescope of Lulin Observatory revealed the binary nature of the object. Utilizing archived photometric data of multi-band surveys, we construct the spectral energy distribution, which is well fitted by a K2V spectral template. The fitting result of the orbital profile using the Wilson--Devinney code suggests that 2MASS J11201034-2201340 is a short-period A-type contact binary and the more massive component has a cool spot. The X-ray emission was first noted in observations made by \swift\, then further confirmed and characterized by an \xmm\ observation. The X-ray spectrum can be described by a power law or thermal Bremsstrahlung. Unfortunately, we could not observe significant X-ray orbital modulation. Finally, according to the spectral energy distribution, this system is estimated to be 690 pc from Earth with a calculated X-ray intensity of $(0.7-1.5)\times 10^{30}$ erg s$^{-1}$, which is in the expected range of an X-ray emitting contact binary.

preprint2016arXiv

Observing two dark accelerators around the Galactic Centre with Fermi Large Area Telescope

We report the results from a detailed $γ-$ray investigation in the field of two "dark accelerators", HESS J1745-303 and HESS J1741-302, with $6.9$ years of data obtained by the Fermi Large Area Telescope. For HESS J1745-303, we found that its MeV-GeV emission is mainly originated from the "Region A" of the TeV feature. Its $γ-$ray spectrum can be modeled with a single power-law with a photon index of $Γ\sim2.5$ from few hundreds MeV to TeV. Moreover, an elongated feature, which extends from "Region A" toward northwest for $\sim1.3^{\circ}$, is discovered for the first time. The orientation of this feature is similar to that of a large scale atomic/molecular gas distribution. For HESS J1741-302, our analysis does not yield any MeV-GeV counterpart for this unidentified TeV source. On the other hand, we have detected a new point source, Fermi J1740.1-3013, serendipitously. Its spectrum is apparently curved which resembles that of a $γ-$ray pulsar. This makes it possibly associated with PSR B1737-20 or PSR J1739-3023.

preprint2015arXiv

Exploring the intrabinary shock from the redback millisecond pulsar PSR J2129-0429

We have investigated the intrabinary shock emission from the redback millisecond pulsar PSR J2129-0429 with XMM-Newton and Fermi. Orbital modulation in X-ray and UV can be clearly seen. Its X-ray modulation has a double-peak structure with a dip in between. The observed X-rays are non-thermal dominant which can be modeled by a power-law with a photon index of ~1.2. Intrabinary shock can be the origin of the observed X-rays. The UV light curve is resulted from the ellipsoidal modulation of the companion. Modeling the UV light curve prefers a large viewing angle. The heating effect of the UV light curve is found to be negligible which suggests the high energy radiation beam of PSR J2129-0429 does not direct toward its companion. On the other hand, no significant orbital modulation can be found in gamma-ray which suggests the majority of the gamma-rays come from the pulsar.

preprint2015arXiv

Searches for millisecond pulsar candidates among the unidentified Fermi objects

Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $γ-$ray sources in this catalog are identified as promising MSP candidates based on their $γ$-ray properties. Through the X-ray analysis, we have detected possible X-ray counterparts, localized to an arcsecond accuracy. We have systematically estimated their X-ray fluxes and compared with the corresponding $γ$-ray fluxes. The X-ray to $γ$-ray flux ratios for 2FGL J1653.6-0159 and 2FGL J1946.4-5402 are comparable with the typical value for pulsars. For 2FGL J1625.2-0020, 2FGL J1653.6-0159 and 2FGL J1946.4-5402, their candidate X-ray counterparts are bright enough for performing a detailed spectral and temporal analysis to discriminate their thermal/non thermal nature and search for the periodic signal. We have also searched for possible optical/IR counterparts at the X-ray positions. For the optical/IR source coincident with the brightest X-ray object that associated with 2FGL J1120.0-2204, its spectral energy distribution is comparable with a late-type star. Evidence for the variability has also been found by examining its optical light curve. All the aforementioned 2FGL sources resemble a pulsar in one or more aspects, which make them as the promising targets for follow-up investigations.

preprint2014arXiv

Chandra detection of a new diffuse X-ray component from the globular cluster 47 Tucanae

In re-analyzing the archival Chandra data of the globular cluster 47 Tucanae, we have detected a new diffuse X-ray emission feature within the half-mass radius of the cluster. The spectrum of the diffuse emission can be described by a power-law model plus a plasma component with photon index $Γ\sim1.0$ and plasma temperature $kT\sim0.2$ keV. While the thermal component is apparently uniform, the non-thermal contribution falls off exponentially from the core. The observed properties could possibly be explained in the context of multiple shocks resulted from the collisions among the stellar wind in the cluster and the inverse Compton scattering between the pulsar wind and the relic photons.

preprint2014arXiv

Discovery of an ultracompact gamma-ray millisecond pulsar binary candidate

We report multi-wavelength observations of the unidentified Fermi object 2FGL J1653.6-0159. With the help of high-resolution X-ray observation, we have identified an X-ray and optical counterpart of 2FGL J1653.6-0159. The source exhibits a periodic modulation of 75 min in optical and possibly also in X-ray. We suggest that 2FGL J1653.6-0159 is a compact binary system with an orbital period of 75 min. Combining the gamma-ray and X-ray properties, 2FGL J1653.6-0159 is potentially a black widow/redback type gamma-ray millisecond pulsar (MSP). The optical and X-ray lightcurve profiles show that the companion is mildly heated by the high-energy emission and the X-rays are from intrabinary shock. Although no radio pulsation has been detected yet, we estimated that the spin period of the MSP is ~2ms based on a theoretical model. If pulsation can be confirmed in the future, 2FGL J1653.6-0159 will become the first ultracompact rotation-powered MSP.

preprint2014arXiv

Fermi-LAT Detection of Pulsed Gamma-rays Above 50 GeV from the Vela Pulsar

The First \fermi-LAT Catalog of Sources Above 10 GeV reported evidence of pulsed emission above 25 GeV from 12 pulsars, including the Vela pulsar, which showed evidence of pulsation at $>37$ GeV energy bands. Using 62 months of \fermi-LAT data, we analyzed the gamma-ray emission from the Vela pulsar and searched for pulsed emission above 50 GeV. Having confirmed the significance of the pulsation in 30-50 GeV with the H-test (p-value $\sim10^{-77}$), we extracted its pulse profile using the Bayesian block algorithm and compared it with the distribution of the 5 observed photons above 50 GeV using the likelihood ratio test. Pulsation was significantly detected for photons above 50 GeV with p-value $=3\times10^{-5}$ ($4.2σ$). The detection of pulsation is significant above $4σ$ at $>79$ GeV and above $3σ$ at $>90$ GeV energy bands, making this the highest energy pulsation significantly detected by the LAT. We explore non-stationary outer gap scenario of the very high-energy emissions from the Vela pulsar.

preprint2014arXiv

High-energy emissions from the gamma-ray binary LS 5039

We study mechanisms of multi-wavelength emissions (X-ray, GeV and TeV gamma-rays) from the gamma-ray binary LS~5039. This paper is composed of two parts. In the first part, we report on results of observational analysis using four year data of \fermi\ Large Area Telescope. Due to the improvement of instrumental response function and increase of the statistics, the observational uncertainties of the spectrum in $\sim$100-300 MeV bands and $>10$GeV bands are significantly improved. The present data analysis suggests that the 0.1-100GeV emissions from LS~5039 contain three different components; (i) the first component contributes to $<$1GeV emissions around superior conjunction, (ii) the second component dominates in 1-10GeV energy bands and (iii) the third component is compatible to lower energy tail of the TeV emissions. In the second part, we develop an emission model to explain the properties of the phase-resolved emissions in multi-wavelength observations. Assuming that LS~5039 includes a pulsar, we argue that both emissions from magnetospheric outer gap and inverse-Compton scattering process of cold-relativistic pulsar wind contribute to the observed GeV emissions. We assume that the pulsar is wrapped by two kinds of termination shock; Shock-I due to the interaction between the pulsar wind and the stellar wind and Shock-II due to the effect of the orbital motion. We propose that the X-rays are produced by the synchrotron radiation at Shock-I region and the TeV gamma-rays are produced by the inverse-Compton scattering process at Shock-II region.

preprint2014arXiv

High-energy observations of PSR B1259-63/LS 2883 through the 2014 periastron passage: connecting X-rays to the GeV flare

The binary system PSR B1259-63/LS 2883 is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital phase-dependent variability in these frequencies. The first detection of GeV gamma-rays from the system was made around the 2010 periastron passage. In this Letter, we present an analysis of X-ray and gamma-ray data obtained by the Swift/XRT, NuSTAR/FPM, and Fermi/LAT, through the recent periastron passage which occurred on 2014 May 4. While PSR B1259-63/LS 2883 was not detected by the LAT before and during this passage, we show that the GeV flares occurred at a similar orbital phase as in early 2011, thus establishing the repetitive nature of the post-periastron GeV flares. Multiple flares each lasting for a few days have been observed and short-term variability is seen as well. We also found X-ray flux variation contemporaneous with the GeV flare for the first time. A strong evidence of the keV-to-GeV connection came from the broadband high-energy spectra, which we interpret as synchrotron radiation from the shocked pulsar wind.

preprint2014arXiv

Multi-wavelength emissions from the millisecond pulsar binary PSR J1023+0038 during an accretion active state

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that 1) its gamma-rays suddenly brightened within a few days in June/July 2013 and has remained at a high gamma-ray state for several months; 2) both UV and X-ray fluxes have increased by roughly an order of magnitude, and 3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100-500 seconds and 50-100 seconds, respectively. Our model suggests that a newly formed accretion disk due to the sudden increase of the stellar wind could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk and the pulsar is still powered by rotation.

preprint2014arXiv

NuSTAR observations and broadband spectral energy distribution modeling of the millisecond pulsar binary PSR J1023+0038

We report the first hard X-ray (3-79 keV) observations of the millisecond pulsar (MSP) binary PSR J1023+0038 using NuSTAR. This system has been shown transiting between a low-mass X-ray binary (LMXB) state and a rotation-powered MSP state. The NuSTAR observations were taken in both LMXB state and rotation-powered state. The source is clearly seen in both states up to ~79 keV. During the LMXB state, the 3-79 keV flux is about a factor of 10 higher that in the rotation-powered state. The hard X-rays show clear orbital modulation during the X-ray faint rotation-powered state but the X-ray orbital period is not detected in the X-ray bright LMXB state. In addition, the X-ray spectrum changes from a flat power-law spectrum during the rotation-powered state to a steeper power-law spectrum in the LMXB state. We suggest that the hard X-rays are due to the intra-binary shock from the interaction between the pulsar wind and the injected material from the low-mass companion star. During the rotation-powered MSP state, the X-ray orbital modulation is due to Doppler boosting of the shocked pulsar wind. At the LMXB state, the evaporating matter of the accretion disk due to the gamma-ray irradiation from the pulsar stops almost all the pulsar wind, resulting the disappearance of the X-ray orbital modulation.

preprint2014arXiv

Search for TeV Gamma-ray Emission from GRB 100621A, an extremely bright GRB in X-rays, with H.E.S.S

The long gamma-ray burst (GRB) 100621A, at the time the brightest X-ray transient ever detected by Swift-XRT in the $0.3\textrm{--}10$ keV range, has been observed with the H.E.S.S. imaging air Cherenkov telescope array, sensitive to gamma radiation in the very-high-energy (VHE, $>100$ GeV) regime. Due to its relatively small redshift of $z\sim0.5$, the favourable position in the southern sky and the relatively short follow-up time ($<700 \rm{s}$ after the satellite trigger) of the H.E.S.S. observations, this GRB could be within the sensitivity reach of the H.E.S.S. instrument. The analysis of the H.E.S.S. data shows no indication of emission and yields an integral flux upper limit above $\sim$380 GeV of $4.2\times10^{-12} \rm cm^{-2}s^{-1}$ (95 % confidence level), assuming a simple Band function extension model. A comparison to a spectral-temporal model, normalised to the prompt flux at sub-MeV energies, constraints the existence of a temporally extended and strong additional hard power law, as has been observed in the other bright X-ray GRB 130427A. A comparison between the H.E.S.S. upper limit and the contemporaneous energy output in X-rays constrains the ratio between the X-ray and VHE gamma-ray fluxes to be greater than 0.4. This value is an important quantity for modelling the afterglow and can constrain leptonic emission scenarios, where leptons are responsible for the X-ray emission and might produce VHE gamma rays.

preprint2014arXiv

The birth of radio millisecond pulsars and their high-energy signature

Millisecond pulsars (MSPs) are thought to born in low-mass X-ray binaries when the neutron star has gained enough angular momentum from the accreting materials of its companion star. It is generally believed that a radio MSP is born when the neutron star stops accreting and enters a rotation-powered state. Exactly what happens during the transition time was poorly understood until a year ago. In the past year, observations have revealed a few objects that not only switched from one state to the other (as predicted in the above picture), but also have swung between the two states within weeks to years. In this work, we present observations of two of these transition objects (PSR J1023+0038 and XSS J12270-4859) and a theoretical framework that tries to explain their high-energy radiation.

preprint2013arXiv

Exploring the X-ray and gamma-ray properties of the redback millisecond pulsar PSR J1723-2837

We have investigated the X-ray and $γ$-ray properties of the redback millisecond pulsar PSR J1723-2837 with XMM-Newton, Chandra and Fermi. We have discovered the X-ray orbital modulation of this binary system with the minimum that coincides with the phases of radio eclipse. The X-ray emission is clearly non-thermal in nature which can be well described by a simple power-law with a photon index of $\sim1.2$. The phase-averaged luminosity is $\sim9\times10^{31}$ erg/s in 0.3-10 keV which consumes $\sim0.2\%$ of the spin-down power. We have detected the $γ-$ray emission in $0.1-300$ GeV from this system at a significance of $\sim6σ$ for the first time. The $γ-$rays in this energy range consumes $\sim2\%$ of the spin-down power and can be modeled by a power-law with a photon index of $\sim2.6$. We discuss the high energy properties of the new redback in the context of a intrabinary shock model.

preprint2013arXiv

High energy emission of GRB 130427A: evidence for inverse Compton radiation

A nearby super-luminous burst GRB 130427A was simultaneously detected by six $γ$-ray space telescopes ({\it Swift}, Fermi-GBM/LAT, Konus-Wind, SPI-ACS/INTEGRAL, AGILE and RHESSI) and by three RAPTOR full-sky persistent monitors. The isotropic $γ-$ray energy release is of $\sim 10^{54}$ erg, rendering it the most powerful explosion among the GRBs with a redshift $z\leq 0.5$. The emission above 100 MeV lasted about one day and four photons are at energies greater than 40 GeV. We show that the count rate of 100 MeV-100 GeV emission may be mainly accounted for by the forward shock synchrotron radiation and the inverse Compton radiation likely dominates at GeV-TeV energies. In particular, an inverse Compton radiation origin is established for the $\sim (95.3,~47.3,~41.4,~38.5,~32)$ GeV photons arriving at $t\sim (243,~256.3,~610.6,~3409.8,~34366.2)$ s after the trigger of Fermi-GBM. Interestingly, the external-inverse-Compton-scattering of the prompt emission (the second episode, i.e., $t\sim 120-260$ s) by the forward-shock-accelerated electrons is expected to produce a few $γ-$rays at energies above 10 GeV, while five were detected in the same time interval. A possible unified model for the prompt soft $γ-$ray, optical and GeV emission of GRB 130427A, GRB 080319B and GRB 090902B is outlined. Implication of the null detection of $>1$ TeV neutrinos from GRB 130427A by IceCube is discussed.

preprint2013arXiv

Search for pulsed gamma-ray emission from globular cluster M28

Using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, we have searched for the gamma-ray pulsations from the direction of globular cluster M28 (NGC 6626). We report the discovery of a signal with the frequency consistent with that of the energetic millisecond pulsar (MSP) PSR B1821-24 in M28. A weighted H-test test statisic (TS) of 28.8 is attained which corresponds to a chance probability of ~1e-5 (4.3-sigma detection). With a phase-resolved analysis, the pulsed component is found to contribute ~25% of the total observed gamma-ray emission from the cluster. On the other hand, the unpulsed level provides a constraint for the underlying MSP population and the fundamental plane relations for the scenario of inverse Compton scattering. Follow-up timing observations in radio/X-ray are encouraged for further investigating this periodic signal candidate.

preprint2013arXiv

The missing GeV γ-ray binary: Searching for HESS J0632+057 with Fermi-LAT

The very high energy (VHE; >100 GeV) source HESS J0632+057 has been recently confirmed as a γ-ray binary, a subclass of the high mass X-ray binary (HMXB) population, through the detection of an orbital period of 321 days. We performed a deep search for the emission of HESS J0632+057 in the GeV energy range using data from the Fermi Large Area Telescope (LAT). The analysis was challenging due to the source being located in close proximity to the bright γ-ray pulsar PSR J0633+0632 and lying in a crowded region of the Galactic plane where there is prominent diffuse emission. We formulated a Bayesian block algorithm adapted to work with weighted photon counts, in order to define the off-pulse phases of PSR J0633+0632. A detailed spectral-spatial model of a 5 deg circular region centred on the known location of HESS J0632+057 was generated to accurately model the LAT data. No significant emission from the location of HESS J0632+057 was detected in the 0.1-100 GeV energy range integrating over ~3.5 years of data; with a 95% flux upper limit of F_{0.1-100 GeV} < 3 x 10-8 ph cm-2 s-1. A search for emission over different phases of the orbit also yielded no significant detection. A search for source emission on shorter timescales (days--months) did not yield any significant detections. We also report the results of a search for radio pulsations using the 100-m Green Bank Telescope (GBT). No periodic signals or individual dispersed bursts of a likely astronomical origin were detected. We estimated the flux density limit of < 90/40 μJy at 2/9 GHz. The LAT flux upper limits combined with the detection of HESS J0632+057 in the 136-400 TeV energy band by the MAGIC collaboration imply that the VHE spectrum must turn over at energies <136 GeV placing constraints on any theoretical models invoked to explain the γ-ray emission.

preprint2012arXiv

Discovery of an unidentified Fermi object as a black widow-like millisecond pulsar

The Fermi Gamma-ray Space Telescope has revolutionized our knowledge of the gamma-ray pulsar population, leading to the discovery of almost 100 gamma-ray pulsars and dozens of gamma-ray millisecond pulsars (MSPs). Although the outer-gap model predicts different sites of emission for the radio and gamma-ray pulsars, until now all of the known gamma-ray MSPs have been visible in the radio. Here we report the discovery of a "radio-quiet" gamma-ray emitting MSP candidate by using Fermi, Chandra, Swift, and optical observations. The X-ray and gamma-ray properties of the source are consistent with known gamma-ray pulsars. We also found a 4.63-hr orbital period in optical and X-ray data. We suggest that the source is a black widow-like MSP with a ~0.1 solar-mass late-type companion star. Based on the profile of the optical and X-ray light-curves, the companion star is believed to be heated by the pulsar while the X-ray emissions originate from pulsar magnetosphere and/or from intra-binary shock. No radio detection of the source has been reported yet and although no gamma-ray/radio pulsation has been found, we estimated that the spin period of the MSP is ~3-5 ms based on the inferred gamma-ray luminosity.

preprint2012arXiv

Discovery of γ-ray pulsation and X-ray emission from the black widow pulsar PSR J2051-0827

We report the discovery of pulsed γ-ray emission and X-ray emission from the black widow millisecond pulsar PSR J2051-0827 by using the data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope and the Advanced CCD Imaging Spectrometer array (ACIS-S) on the Chandra X-ray Observatory. Using 3 years of LAT data, PSR J2051-0827 is clearly detected in γ-ray with a signicance of \sim 8σ in the 0.2 - 20 GeV band. The 200 MeV - 20 GeV γ-ray spectrum of PSR J2051-0827 can be modeled by a simple power- law with a photon index of 2.46 \pm 0.15. Significant (\sim 5σ) γ-ray pulsations at the radio period were detected. PSR J2051-0827 was also detected in soft (0.3-7 keV) X-ray with Chandra. By comparing the observed γ-rays and X-rays with theoretical models, we suggest that the γ-ray emission is from the outer gap while the X-rays can be from intra-binary shock and pulsar magnetospheric synchrotron emissions.

preprint2012arXiv

Orbital Phase Dependent Gamma-ray emissions from the Black Widow Pulsar

We report on evidence for orbital phase-dependence of the gamma-ray emission from PSR B1957+20 black widow system by using the data of the Fermi Large Area Telescope. We divide an orbital cycle into two regions: a region containing the inferior conjunction, and the other region containing rest of the orbital cycle. We show that the observed spectra for the different orbital regions are fitted by different functional forms. The spectrum of the orbital region containing inferior conjunction can be described by a power-law with an exponential cutoff (PLE) model, which gives the best-fit model for the orbital phase that does not contain the inferior conjunction, plus an extra component above ~2.7 GeV. The emission above 3 GeV in this region is detected with a ~7-sigma confidence level. The gamma-ray data above ~2.7 GeV are observed to be modulated at the orbital period at the ~2.3-sigma level. We anticipate that the PLE component dominating below ~2.7 GeV originates from the pulsar magnetosphere. We also show that the inverse-Compton scattering of the thermal radiation of the companion star off a "cold" ultra-relativistic pulsar wind can explain the extra component above ~2.7 GeV. The black widow pulsar PSR B1957+20 may be the member of a new class of object, in the sense that the system is showing gamma-ray emission with both magnetospheric and pulsar wind origins.

preprint2012arXiv

Pulsed γ-ray emission from magnetar 1E 2259+586

Anomalous x-ray pulsars (AXPs) are thought to be magnetars which are young isolated neutron stars with extremely strong magnetic fields of > 10^14Gauss. Their tremendous magnetic fields inferred from the spin parameters provide a huge energy reservoir to power the observed x-ray emission. High-energy emission above 0.3 MeV has never been detected despite intensive search. Here, we present the possible Fermi Large Area Telescope (LAT) detection of γ-ray pulsations above 200 MeV from the AXP, 1E 2259+586, which puts the current theoretical models of γ-ray emission mechanisms of magnetars into challenge. We speculate that the high-energy γ-rays originate from the outer magnetosphere of the magnetar.

preprint2011arXiv

Discovery of gamma-ray emission from the supernova remnant Kes 17 with Fermi Large Area Telescope

We report the discovery of GeV emission at the position of supernova remnant Kes 17 by using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Kes 17 can be clearly detected with a significance of ~12 sigma in the 1 - 20 GeV range. Moreover, a number of gamma-ray sources were detected in its vicinity. The gamma-ray spectrum of Kes 17 can be well described by a simple power-law with a photon index of ~ 2.4. Together with the multi-wavelength evidence for its interactions with the nearby molecular cloud, the gamma-ray detection suggests that Kes 17 is a candidate acceleration site for cosmic-rays.

preprint2011arXiv

Exploring the dark accelerator HESS J1745-303 with Fermi Large Area Telescope

We present a detailed analysis of the gamma-ray emission from HESS J1745-303 with the data obtained by the Fermi Gamma-ray Space Telescope in the first ~29 months observation.The source can be clearly detected at the level of ~18-sigma and ~6-sigma in 1-20 GeV and 10-20 GeV respectively. Different from the results obtained by the Compton Gamma-ray Observatory, we do not find any evidence of variability. Most of emission in 10-20 GeV is found to coincide with the region C of HESS J1745-303. A simple power-law is sufficient to describe the GeV spectrum with a photon index of ~2.6. The power-law spectrum inferred in the GeV regime can be connected to that of a particular spatial component of HESS J1745-303 in 1-10 TeV without any spectral break. These properties impose independent constraints for understanding the nature of this "dark particle accelerator".

preprint2011arXiv

Gamma-ray emission from the globular clusters Liller 1, M80, NGC 6139, NGC 6541, NGC 6624, and NGC 6752

Globular clusters (GCs) are emerging as a new class of gamma-ray emitters, thanks to the data obtained from the Fermi Gamma-ray Space Telescope. By now, eight GCs are known to emit gamma-rays at energies >100~MeV. Based on the stellar encounter rate of the GCs, we identify potential gamma-ray emitting GCs out of all known GCs that have not been studied in details before. In this paper, we report the discovery of a number of new gamma-ray GCs: Liller 1, NGC 6624, and NGC 6752, and evidence for gamma-ray emission from M80, NGC 6139, and NGC 6541, in which gamma-rays were found within the GC tidal radius. With one of the highest metallicity among all GCs in the Milky Way, the gamma-ray luminosity of Liller 1 is found to be the highest of all known gamma-ray GCs. In addition, we confirm a previous report of significant gamma-ray emitting region next to NGC 6441. We briefly discuss the observed offset of gamma-rays from some GC cores. The increasing number of known gamma-ray GCs at distances out to ~10 kpc is important for us to understand the gamma-ray emitting mechanism and provides an alternative probe to the underlying millisecond pulsar populations of the GCs.

preprint2011arXiv

The Fundamental Plane of Gamma-ray Globular Clusters

We have investigated the properties of a group of $γ$-ray emitting globular clusters (GCs) which have recently been uncovered in our Galaxy. By correlating the observed $γ$-ray luminosities $L_γ$ with various cluster properties, we probe the origin of the high energy photons from these GCs. We report $L_γ$ is positively correlated with the encounter rate $Γ_{c}$ and the metalicity $\left[{\rm Fe/H}\right]$ which place an intimate link between the gamma-ray emission and the millisecond pulsar population. We also find a tendency that $L_γ$ increase with the energy densities of the soft photon at the cluster location. Furthermore, the two-dimensional regression analysis suggests that $L_γ$, soft photon densities, and $Γ_{c}$/$\left[{\rm Fe/H}\right]$ possibly span fundamental planes which potentially provide better predictions for the $γ$-ray properties of GCs.