The first measurements from the CIBER experiment of extragalactic background light (EBL) in near-infrared (NIR) band exhibit a higher intensity than those inferred through $\gamma$-ray observations. Recent theoretical-EBL intensities are typically consistent with the very high energy (VHE) $\gamma$-ray observations. Yet, it is possible that the excess NIR radiation is a new component of EBL and not in tension with the TeV spectra of distant blazars, since the hypothetical axion-like particle (ALP) may lead to a reduced opacity of the Universe for VHE $\gamma$-rays. In order to probe whether the excess component arises mainly from EBL, thirteen observed spectra in high energy and VHE ranges from ten distant TeV BL Lac objects are fitted by four theoretical spectra which involve theoretical EBL (Gilmore ), Gilmore's EBL model including photon/ALP coupling, Gilmore's EBL with CIBER excess and the latter including photon/ALP coupling respectively. We find the goodness of fit for the model with CIBER excess can be improved with a significance of $7.6~\sigma$ after including the photon/ALP coupling; Thus, the ALP/photon mixing mechanism can effectively alleviate the tension; However, the Gilmore EBL model, on the whole, is more compatible with the observed spectra compared to those with ALP, although individual blazars such as PKS 1424+240 and 1ES 1101-232 prefer the ALP-model. Our results suggest that the recent EBL models can solely explain the VHE $\gamma$-ray observations, and assuming the existence of the ALP to alleviate the tension is not required in a statistical sense, thus the excess over the EBL models is less likely to be a new EBL component.

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.

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.

We analyzed 199 GRBs detected by Fermi-LAT during the years 2008-2019. We found 67 photons at energies >=10 GeV, which come from 34 GRBs. Out of these 34 GRBs, Fermi-LAT detects significant (TS>=4) afterglow 0.1-200 GeV photons from 25 GRBs. We present time-integrated 0.1-200 GeV spectra of these 25 GRBs. The spectra of a significant fraction (9/25) of these GRBs revealed a harder component above an energy break of 0.3-2 GeV. While shock synchrotron may account for the photons at the lower energy end, high energy photons above the break is naturally explained by synchrotron self-Compton (SSC) emission. We perform broadband model fit to the X-ray-LAT emission of GRB 131231A. Comparing the afterglow spectra of these 25 GRBs with other Fermi-LAT detected GRBs, we found that the power-law index distribution is similar for the two populations. This may indicate that the additional high-energy component may also exist in Fermi-LAT GRBs in general.

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 ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\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}.

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.

We investigate the nature of the unidentified very-high-energy (VHE) gamma-ray object, HESS J1832-093, in a multi-wavelength context. Based on X-ray variability and spectral index ($\Gamma_X\sim\,1.5$), and its broad-band spectrum (which was remarkably similar to HESS J0632+057, a confirmed gamma-ray binary), HESS J1832-093 has been considered to be a strong gamma-ray binary candidate in previous works. In this work, we provide further evidence for this scenario. We obtained a spectrum of its IR counterpart using Gemini/Flamingo, finding absorption lines that are usually seen in massive stars, in particular O stars. We also obtained a rather steep ATCA spectrum ($\alpha=-1.18^{+1.04}_{-0.88}$) which prefers a gamma-ray binary over an AGN scenario. Based on spatial-spectral analysis and variability search, we found that 4FGL J1832.9-0913 is possible to be associated with SNR G22.7-0.2 rather than with HESS J1832-093 only.

The supernova remnant (SNR) HESS J1731-347 is a young SNR which displays a non-thermal X-ray and TeV shell structure. A molecular cloud at a distance of 3.2 kpc is spatially coincident with the western part of the SNR, and it is likely hit by the SNR. The X-ray emission from this part of the shell is much lower than from the rest of the SNR. Moreover, a compact GeV emission region coincident with the cloud has been detected with a soft spectrum. These observations seem to imply a shock-cloud collision scenario at this area, where the stalled shock can no longer accelerate super-TeV electrons or maintain strong magnetic turbulence downstream, while the GeV cosmic rays (CRs) are released through this stalled shock. To test this hypothesis, we have performed a detailed Fermi-LAT reanalysis of the HESS J1731-347 region with over 9 years of data. We find that the compact GeV emission region displays a spectral power-law index of -2.4, whereas the GeV emission from the rest of the SNR (excluding the cloud region) has an index of -1.8. A hadronic model involving a shock-cloud collision scenario is built to explain the -ray emission from this area. It consists of three CR sources: run-away super-TeV CRs that have escaped from the fast shock, leaked GeV CRs from the stalled shock, and the local CR sea. The X-ray and -ray emission of the SNR excluding the shock-cloud interaction region is explained in a one-zone leptonic model. Our shock-cloud collision model explains well the GeV-TeV observations from both cloud regions around HESS J1731-347, i.e. from the cloud in contact with the SNR and from the more distant cloud which is coincident with the nearby TeV source HESS J1729-345. We find however that the leaked GeV CRs from the shock-cloud collision do not necessarily dominate the GeV emission from the clouds, due to a comparable contribution from the local CR sea.

We report X-ray data analysis results obtained from Chandra, XMM-Newton, NuSTAR and Swift observations of PSR J2032+4127 taken before, during, and after the periastron on 2017 November 13. We found the first clear evidence of a change in the X-ray spectral index over the passage period, thanks to a broad and sensitive spectral coverage by XMM-Newton and NuSTAR. We analysed the joint XMM-Newton and NuSTAR observation epochs with power-law and broken power-law model. We have obtained change in spectral parameters before and after the periastron passage for both models. The spectra get softened after the passage. The evolution of the spectral index and break energy before and after the periastron may indicate a change in the physical state of shock-accelerated electrons.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

GRB 190114C is a famous Gamma-Ray Burst (GRB) due to its detection at sub-TeV energies by MAGIC, seen at redshift z = 0.42. This burst is one of the brightest GRB detected by \fermi. We present a detailed analysis of GRB 190114C prompt emission, using the two \fermi detectors: GBM and LAT. The \emph{LAT low energy events} (LLE) data is also considered. A joint GBM-LAT analysis reveals a sub-$GeV$ spectral cutoff. A similar high energy cutoff was reported in GRB 160509A and GRB 100724B earlier, and a handful of other sources. The cutoff can be explained by the intrinsic opacity due to pair production within the emitting region. Such morphology in these GRBs suggests that they belong to one specific class having a similar source of the radiation mechanism. GRB 190114C shows a transition from non-thermal to a quasi-thermal-like spectrum along with radiation due to external shock. From spectrum analysis and $Lorentz$ factor evolution from the trigger time to late emission, considering the fact that sub-TeV photons are detected in MAGIC, we are able to draw an emission mechanism picture, where the prompt emission spectrum is more consistent with spectrum via photospheric dissipation with presence of external shock emission simultaneously.

The three major large-scale, diffuse $\gamma$-ray structures of the Milky Way are the Galactic disk, a bulge-like GeV excess towards the Galactic center, and the Fermi bubble. Whether such structures can also be present in other normal galaxies remains an open question. M31, as the nearest massive normal galaxy, holds promise for spatially-resolving the $\gamma$-ray emission. Based on more than 8 years of Fermi-LAT observations, we use (1) disk, (2) bulge, and (3) disk-plus-bulge templates to model the spatial distribution of the $\gamma$-ray emission from M31. Among these, the disk-plus-bulge template delivers the best-fit, in which the bulge component has a TS value 25.7 and a photon-index of $2.57\pm0.17$, providing strong evidence for a centrally-concentrated gamma-ray emission from M31, that is analogous to the Galactic center excess. The total 0.2-300 GeV gamma-ray luminosity from this bulge component is $(1.16\pm0.14)\times10^{38} erg s^{-1}$, which would require $\sim1.5\times10^5$ millisecond pulsars, if they were the dominant source. We also search for a Fermi bubble-like structure in M31 using the full dataset (pass8), but no significant evidence is found. In addition, a likelihood analysis using only photons with the most accurate reconstructed direction (i.e., PSF3-only data) reveals a 4.8 $\sigma$ point-like source located at $\sim$ 10kpc to the northwest of the M31 disk, with a luminosity of $(0.97\pm0.27)\times10^{38} erg s^{-1}$ and a photon-index of $2.31\pm0.18$. Lacking of a counterpart on the southeast side of the disk, the relation between this point-like source and a bubble-like structure remains elusive.

PSR J1119-6127 is a radio pulsar that behaved magnetar like bursts, and we performed a comprehensive investigation of this pulsar using the archival high-energy observations obtained after its outburst in 2016 July. After the 2016 outburst, specific regions on the neutron star surface were heated up to > 0.3 and > 1 keV from about 0.2 keV. A hard non-thermal spectral component with a photon index < 0.5 related to the magnetospheric emission can be resolved from the NuSTAR spectra above 10 keV. We find that the thermal emitting regions did not cool down and gradually shrank by about (20-35)% four months after the outburst. Hard X-ray pulsations were detected with NuSTAR immediately after the outburst, at 5 sigma confidence level and with a background-subtracted pulsed fraction of (40+-10)%. However, the signal became undetectable after a few days. Using Fermi data, we found that the gamma-ray emission in 0.5-300 GeV was suppressed along with the disappearance of the radio pulsations. This is likely caused by a reconfiguration of the magnetic field. We also discovered that the timing noise evolved dramatically, and the spin-down rate significantly increased after the 2016 glitch. We proposed that post-outburst temporal and spectral behaviors from radio to gamma-ray bands were caused by changes of the magnetosphere structure, pair plasma injection and the shrinking emission sites on the neutron star.

We report a detailed analysis of the orbital properties of binary millisecond pulsar (MSP) with a white dwarf (WD) companion. Positive correlations between the orbital period $P_{\rm b}$ and eccentricity $\epsilon$ are found in two classes of MSP binaries with a He WD and with a CO/ONeMg WD, though their trends are different. The distribution of $P_{\rm b}$ is not uniform. Deficiency of sources at $P_{\rm b}\sim35-50$~days (Gap 1) have been mentioned in previous studies. On the other hand, another gap at $P_{\rm b}\sim2.5-4.5$~days (Gap 2) is identified for the first time. Inspection of the relation between $P_{\rm b}$ and the companion masses $M_{\rm c}$ revealed the subpopulations of MSP binaries with a He WD separated by Gap 1, above which $P_{\rm b}$ is independent of $M_{\rm c}$ (horizontal branch) but below which $P_{\rm b}$ correlates strongly with $M_{\rm c}$ (lower branch). Distinctive horizontal branch and lower branch separated by Gap 2 were identified for the MSP binaries with a CO/ONeMg WD at shorter $P_{\rm b}$ and higher $M_{\rm c}$. Generally, $M_{\rm c}$ are higher in the horizontal branch than in the lower branch for the MSP binaries with a He WD. These properties can be explained in terms of a binary orbital evolution scenario in which the WD companion was ablated by a pulsar wind in the post mass-transfer phase.

We have conducted a systematic survey for the X-ray properties of millisecond pulsars (MSPs). Currently, there are 47 MSPs with confirmed X-ray detections. We have also placed the upper limits for the X-ray emission from the other 36 MSPs by using the archival data. We have normalized their X-ray luminosities $L_{x}$ and their effective photon indices $\Gamma$ into a homogeneous data set, which enable us to carry out a detailed statistical analysis. Based on our censored sample, we report a relation of $L_{x}\simeq10^{31.05}\left(\dot{E}/10^{35}\right)^{1.31}$ erg/s (2-10 keV) for the MSPs. The inferred X-ray conversion efficiency is found to be lower than previously reported estimate that could be affected by selection bias. $L_{x}$ also correlates/anti-correlates with the magnetic field strength at the light cylinder $B_{LC}$/characteristic age $\tau$. On the other hand, there is no correlation between $L_{x}$ and their surface magnetic field strength $B_{s}$. We have further divided the sample into four classes: (i) black-widows, (ii) redbacks, (iii) isolated MSPs and (iv) other MSP binaries, and compare the properties among them. We noted that while the rotational parameters and the orbital periods of redbacks and black-widow are similar, $L_{x}$ of redbacks are significantly higher than those of black-widows in the 2-10 keV band. Also the $\Gamma$ of redbacks are apparently smaller than those of black-widows, which indicates the X-ray emission of redbacks are harder than that of black-widows. This can be explained by the different contribution of intrabinary shocks in the X-ray emission of these two classes.

GeV flares from PSR B1259-63/LS 2883 were seen starting around 30 days after the two periastron passages in 2010 and 2014. The flares are clearly delayed compared to the occurrence of the X-ray and TeV flux peaks during the post-periastron disk crossing. While several attempts have been put forward to explain this phenomenon, the origin of these GeV flares remains a puzzle. Here we present a detailed analysis of the observational data taken by the Fermi and Swift observatories over the 2017 September periastron passage. For the first time, we find short-lived but powerful GeV flares on time scales of down to three hours. The onset of the GeV flaring period in 2017 is also delayed compared to those seen in 2011 and 2014. Supplemented by a re-analysis of previous data, we compare the Fermi/LAT, Swift/XRT and Swift/UVOT light curves in 2017 with those taken over the 2010 and 2014 periastrons, and difference in UVOT light curves are noted.

We present the first results from a joint XMM-Newton and NuSTAR observation of the gamma-ray emitting millisecond pulsar compact binary PSR J2129-0429. X-ray emission up to about 40 keV is detected and the joint spectrum can be modeled with a power-law plus a neutron star atmosphere model. At a distance of 1.4 kpc, the 0.3-79 keV luminosity is $3.56\times10^{32}$ erg s$^{-1}$. We also detected the 0.64-day binary orbital period with a double-peaked structure across the wavebands. By combining the updated Fermi GeV data, we modeled the broadband spectral energy distribution as well as the X-ray modulation with an intrabinary model involving shock interaction between pulsar wind and outflow from the companion star. Lastly, we reported a high-resolution X-ray image provided by Chandra to rule out the proposed pulsar wind nebula associated with PSR J2129-0429.

Supernova remnant (SNR)\,W28 is well known for its classic hadronic scenario, in which the TeV cosmic rays (CRs) released at the early stage of this mid-aged SNR are illuminating nearby molecular clouds (MCs). Overwhelming evidences have shown that the northeast of the SNR (W28-North) has already encountered with the MC clumps. Through this broken shell -- W28-North, we believe the CRs with energy down to $<$1\,GeV to be able to be injected into nearby MCs. To further testify this hadronic scenario, we first analyse the 9 years Fermi-LAT data in/around W28 with energy down to 0.3\,GeV. Our Fermi-LAT analysis display a 10-200 GeV skymap which spatially matches well with the known TeV sources -- HESS\,J1801-233 (W28-North), HESS\,J1800-240\,A,\,B\,\&\,C (240\,A\,B\,\&\,C). At low energy band, we has discovered a 0.5-1\,GeV blob located to the south of 240\,B\,\&\,C, and a low flux of 0.3-1\,GeV at 240\,A. A hadronic model is build to explain our analysis results and previous multi-wavelength observations of W28. Our model consists of three CR sources: The run-away CRs escaped from a strong shock; The leaked GeV CRs from the broken shell -- W28-North; And the local CR sea. Through modelling the SNR evolution, CR acceleration\,\&\,releasing, we have explained the GeV-TeV emission in/around SNR\,W28 (except for 240\,A) in one model. Both the damping of the magnetic waves by the neutrals and the decreased acceleration efficiency are taken into account in our model due to the mid-age of SNR W28.

We present the Neil Gehrels Swift Observatory (Swift), Fermi Large Area Telescope (Fermi-LAT), and Karl G. Jansky Very Large Array (VLA) observations of the gamma-ray binary PSR J2032+4127/MT91 213, of which the periastron passage has just occurred in November 2017. In the Swift X-ray light curve, the flux was steadily increasing before mid-October 2017, however, a sharp X-ray dip on a weekly time-scale is seen during the periastron passage, followed by a post-periastron X-ray flare lasting for ~20 days. We suggest that the X-ray dip is caused by (i) an increase of the magnetization parameter at the shock, and (ii) the suppression due to the Doppler boosting effect. The 20-day post-periastron flare could be a consequence of the Be stellar disk passage by the pulsar. An orbital GeV modulation is also expected in our model, however, no significant variability is seen in the Fermi-LAT light curve. We suspect that the GeV emission resulted from the interaction between the binary's members is hidden behind the bright magnetospheric emission of the pulsar. Pulsar gating technique would be useful to remove the magnetospheric emission and recover the predicted GeV modulation, if an accurate radio timing solution over the periastron passage is provided in the future.

The angular distribution of GRB jets is not yet clear. The observed luminosity of GRB 170817A is the lowest among all known short GRBs, which is best explained by the fact that our line of sight is outside of the jet opening angle, $\theta_{obs}$ $>$ $\theta_j$, where $\theta_{obs}$ is the angle between our line of sight and jet axis. Inferred by gravitational wave observations, as well as radio and X-ray afterglow modeling of GRB 170817A, it is likely that $\theta_{obs}$ $\sim$ 20$^{\circ}$ -- 28$^{\circ}$. In this work, we quantitatively consider two scenarios of angular energy distribution of GRB ejecta: top-hat jet, and structured jet with a power law index $s$. For top-hat jet model, we get a large $\theta_j$ (e.g., $\theta_{j}$$>$$10^{\circ}$), a rather high local (i.e., z$<$0.01) short GRB rate $\sim$8--15$\times10^3$ $Gpc^{-3}$$yr^{-1}$ (estimated 90$\sim$1850 $ Gpc^{-3}$$yr^{-1}$ in Fong et al. 2015) and an extremely high $E_{peak,0}$(on-axis,rest-frame)$>$7.5$\times10^4$keV($\sim$500keV for typical short GRB). For structured jet model, we use $\theta_{obs}$ to give limits on $s$ and $\theta_{j}$ for a typical on-axis luminosity short GRB (e.g., $10^{49}$ erg s$^{-1}$$\sim$$10^{51}$ erg s$^{-1}$), and a low on-axis luminosity case (e.g., $10^{49}$ erg s$^{-1}$) gives more reasonable values of $s$. A structured jet model is more feasible for GRB 170817A than a top-hat jet model, due to the rather high local short GRB rate and the extremely high on-axis $E_{peak,0}$ almost rules out the top-hat jet model. GRB 170817A is likely a low on-axis luminosity GRB ($10^{49}$ erg s$^{-1}$) with a structured jet.

We report the analysis result of UV/X-ray emission from AR~Scorpii, which is an intermediate polar (IP) composed of a magnetic white dwarf and a M-type star, with the XMM-Newton data. The X-ray/UV emission clearly shows a large variation over the orbit, and their intensity maximum (or minimum) is located at the superior conjunction (or inferior conjunction) of the M-type star orbit. The hardness ratio of the X-ray emission shows a small variation over the orbital phase, and shows no indication of the absorption by an accretion column. These properties are naturally explained by the emission from the M-type star surface rather than from the accretion column on the WD's star similar to the usual IPs. Beside, the observed X-ray emission also modulates with WD's spin with a pulse fraction of $\sim 14\%$. The peak position is aligned in the optical/UV/X-ray band. This supports the hypothesis that the electrons in AR~Scorpii are accelerated to a relativistic speed, and emit non-thermal photons via the synchrotron radiation. In the X-ray bands, the evidence of the power-law spectrum is found in the pulsed component, although the observed emission is dominated by the optically thin thermal plasma emissions with several different temperatures. It is considered that the magnetic dissipation/reconnection process on the M-type star surface heats up the plasma to a temperature of several keV, and also accelerates the electrons to the relativistic speed. The relativistic electrons are trapped in the WD's closed magnetic field lines by the magnetic mirror effect. In this model, the observed pulsed component is explained by the emissions from the first magnetic mirror point.

The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the "Band" function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each "sub-burst" being $\sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{\rm p, iso}\sim 4\times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.

The High Energy Stereoscopic System (H.E.S.S.) has been searching for counterparts of Gamma Ray Bursts (GRBs) for many years. In 2012 the system was upgraded with a fifth $28$ m diameter telescope (CT5) which is equipped with faster motors for rapid repointing, marking the start of the second phase of H.E.S.S. operation (H.E.S.S. II). CT5s large light collection area of $600\,{\rm m}^{2}$ improves the sensitivity to low-energy gamma-rays and even extends the energy range below $100$ GeV. The search for counterparts continues now in the energy range of tens of GeV to tens of TeV. A detection in this energy range would open a new window to the part of the spectrum of these highly energetic explosions which Fermi-LAT has only successfully detected in a reduced subset of events, with rather limited statistics. In the past years, H.E.S.S. has performed followup observations based on GRB detections by Swift-BAT and Fermi-GBM/-LAT. This Target of Opportunity observation program was carried out with a generalised Target of Opportunity Alert system. This contribution will highlight key features of the Target of Opportunity Alert system, present follow-up statistics of GRBs as well as detailed results of promising follow-up observations.

A recent study by Posselt et al. (2017) reported the deepest X-ray investigation of the Geminga pulsar wind nebula (PWN) by using \emph{Chandra X-ray Observatory}. In comparison with previous studies of this system, a number of new findings have been reported and we found these suggest the possible variabilities in various components of this PWN. This motivates us to carry out a dedicated search for the morphological and spectral variations of this complex nebula. We have discovered variabilities on timescales from a few days to a few months from different components of the nebula. The fastest change occurred in the circumstellar environment at a rate of 80 per cent of the speed of light. One of the most spectacular results is the wiggling of a half light-year long tail as an extension of the jet, which is significantly bent by the ram pressure. The jet wiggling occurred at a rate of about 20 per cent of the speed of light. This twisted structure can possibly be a result of a propagating torsional Alf\`{v}en wave. We have also found evidence of spectral hardening along this tail for a period of about nine months.

We search for the gamma-ray counterparts of stellar-mass black holes using long-term Fermi archive to investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, by applying the pulsar outer-gap model to their magnetosphere. When a black hole transient (BHT) is in a low-hard or quiescent state, the radiatively inefficient accretion flow cannot emit enough MeV photons that are required to sustain the force-free magnetosphere in the polar funnel via two-photon collisions. In this charge-starved gap region, an electric field arises along the magnetic field lines to accelerate electrons and positrons into ultra-relativistic energies. These relativistic leptons emit copious gamma-rays via the curvature and inverse-Compton (IC) processes. It is found that these gamma-ray emissions exhibit a flaring activity when the plasma accretion rate stays typically between 0.01 and 0.005 percent of the Eddington value for rapidly rotating, stellar-mass black holes. By analyzing the detection limit determined from archival Fermi/LAT data, we find that the 7-year averaged duty cycle of such flaring activities should be less than 5% and 10% for XTE J1118+480 and 1A 0620-00, respectively, and that the detection limit is comparable to the theoretical prediction for V404 Cyg. It is predicted that the gap emission can be discriminated from the jet emission, if we investigate the high-energy spectral behaviour or observe nearby BHTs during deep quiescence simultaneously in infrared wavelength and very-high energies.

We present the high-energy emission properties of GRB 160509A, from its prompt mission to late afterglow phase. GRB 160509A contains two emission episodes: 0-40s and 280-420s after the burst onset (t0). The relatively high fluence of GRB 160509A allows us to establish an evolving spectrum above 100 MeV. During the first emission episode, the >100 MeV spectrum is soft with \Gamma=>3.0, which can be smoothly connected to keV energies with a Band function with a high-energy cutoff. The >100 MeV spectrum rapidly changes to a hard spectrum with \Gamma<=1.5 after t0+40s. The existence of very energetic photons, e.g., a 52 GeV that arrives t0+77 seconds, and a 29 GeV that arrives t0+70 ks, is hard to reconcile by the synchrotron emission from forward-shock electrons, but likely due to inverse Compton mechanism (e.g., synchrotron self-Compton emission). A soft spectrum (\Gamma~2) between 300s and 1000s after the burst onset is also found at a significance of about 2 standard deviations, which suggests a different emission mechanism at work for this short period of time. GRB 160509A represents the latest example where inverse Compton emission has to be taken into account in explaining the afterglow GeV emission, which had been suggested long before the launch of Fermi LAT.

A glitch of a pulsar is known as a sudden increase in the spin frequency and spin-down rate (frequency time derivative), and it can be caused by a sudden rel\ ease of the stress built up in the solid crust of the star or pinned vortices in the superfluid interior. PSR J2021+4026 is the first pulsar that shows a significant change in the gamma-ray flux and pulse profile at the glitch that occurred around 2011 October 16. We report the results of timing and spectral analysis of PSR~J2021+4026 using $\sim$ 8~yr Fermi-LAT data. We find that the pulsar stayed at a high spin-down rate ($\sim 4\%$ higher than the pre-glitch value) and a low gamma-ray state ($\sim 18\%$ lower) for about 3~yr after the glitch. Around 2014 December, the spin-down rate and gamma-ray flux gradually returned to pre-glitch values within a time scale of a few months. The phase-resolved spectra and pulse profiles after the relaxation are also consistent with those before the glitch. The observed long-term evolution of the spin-down rate and the gamma-ray flux indicates that the glitch triggered a mode change in the global magnetosphere. We speculate that the glitch changed the local magnetic field structure around the polar cap and/or the inclination angle of the dipole axis, leading to a change in the electric current circulating in the magnetosphere.

Around a rapidly rotating black hole (BH), when the plasma accretion rate is much less than the Eddington rate, the radiatively inefficient accretion flow (RIAF) cannot supply enough MeV photons that are capable of materializing as pairs. In such a charge-starved BH magnetosphere, the force-free condition breaks down in the polar funnels. Applying the pulsar outer-magnetospheric lepton accelerator theory to super-massive BHs, we demonstrate that a strong electric field arises along the magnetic field lines in the direct vicinity of the event horizon in the funnels, that the electrons and positrons are accelerated up to 100~TeV in this vacuum gap, and that these leptons emit copious photons via inverse-Compton (IC) process between 0.1~TeV and 30~TeV for a distant observer. It is found that these IC fluxes will be detectable with Imaging Atmospheric Cherenkov Telescopes, provided that a low-luminosity active galactic nucleus is located within 1~Mpc for a million-solar-mass central BH or within 30~Mpc for a billion-solar-mass central BH. These very-high-energy fluxes are beamed in a relatively small solid angle around the rotation axis because of the inhomogeneous and anisotropic distribution of the RIAF photon field, and show an anti-correlation with the RIAF submillimeter fluxes. The gap luminosity little depends on the three-dimensional magnetic-field configuration, because the Goldreich-Julian charge density, and hence the exerted electric field is essentially governed by the frame-dragging effect, not by the magnetic field configuration.

We report our recent Swift, NuSTAR, and XMM-Newton X-ray and Lijiang optical observations on PSR J2032+4127/MT91 213, the gamma-ray binary candidate with a period of 45-50 years. The coming periastron of the system was predicted to be in November 2017, around which high-energy flares from keV to TeV are expected. Recent studies with Chandra and Swift X-ray observations taken in 2015/16 showed that its X-ray emission has been brighter by a factors of ~10 than that before 2013, probably revealing some on-going activities between the pulsar wind and the stellar wind. Our new Swift/XRT lightcurve shows no strong evidence of a single vigorous brightening trend, but rather several strong X-ray flares on weekly to monthly timescales with a slowly brightening baseline, namely the low state. The NuSTAR and XMM-Newton observations taken during the flaring and the low states, respectively, show a denser environment and a softer power-law index during the flaring state, implying that the pulsar wind interacted with stronger stellar winds of the companion to produce the flares. These precursors would be crucial in studying the predicted giant outburst from this extreme gamma-ray binary during the periastron passage in late 2017.

It has been suggested that high-energy gamma-ray emission ($>100{\rm MeV}$) of nearby star-forming galaxies may be produced predominantly by cosmic rays colliding with the interstellar medium through neutral pion decay. The pion-decay mechanism predicts a unique spectral signature in the gamma-ray spectrum, characterized by a fast rising spectrum and a spectral break below a few hundreds of MeV. We here report the evidence of a spectral break around 500 MeV in the disk emission of Large Magellanic Cloud (LMC), which is found in the analysis of the gamma-ray data extending down to 60 MeV observed by {\it Fermi}-Large Area Telescope. The break is well consistent with the pion-decay model for the gamma-ray emission, although leptonic models, such as the electron bremsstrahlung emission, cannot be ruled out completely.

The non-repeating fast radio bursts (FRBs) could arise from catastrophic stellar explosions or magnetar giant flares, so relativistic blast waves might be produced in these events. Motivated by this, we here search for GeV counterparts to all non-repeating FRBs with Fermi Large Area Telescope (LAT), including FRB 131104 that is claimed to be possibly associated with a gamma-ray transient candidate detected by Swift Burst Alert Telescope (BAT). FRB 131104 enters the field of view (FoV) of LAT ~5000s after the burst time, so we are only able to search for the GeV afterglow emission during this period, but no significant GeV emission is detected. We also perform a search for GeV emissions from other FRBs, but no significant GeV emissions are detected either. Upper limit fluences in the range of (4.7-29.2)x10^(-7) erg/cm2 are obtained, and then the upper limits of the isotropic blast wave kinetic energy of about (1-200)x10^(53) erg are inferred under certain assumptions. Although the current limits on the isotropic blast wave energy are not sufficiently stringent to rule out the connection between FRBs and GRB-like transients, future more sensitive observations with Fermi or Imaging Atmospheric Cherenkov Telescopes might be able to constrain the connection.

We report on the first NuSTAR observation of the gamma-ray emitting millisecond pulsar binary PSR J1723-2837. X-ray radiation up to 79 keV is clearly detected and the simultaneous NuSTAR and Swift spectrum is well described by an absorbed power-law with a photon index of ~1.3. We also find X-ray modulations in the 3-10 keV, 10-20 keV, 20-79 keV, and 3-79 keV bands at the 14.8-hr binary orbital period. All these are entirely consistent with previous X-ray observations below 10 keV. This new hard X-ray observation of PSR J1723-2837 provides strong evidence that the X-rays are from the intrabinary shock via an interaction between the pulsar wind and the outflow from the companion star. We discuss how the NuSTAR observation constrains the physical parameters of the intrabinary shock model.

The supernova remnant Kes 73 and/or the magnetar 1E 1841-045 at its center can deposit a large amount of energy to the surroundings and is potentially responsible for particle acceleration. Using the data taken with the Fermi Large Area Telescope (LAT), we confirmed the presence of an extended source whose centroid position is highly consistent with this magnetar/supernova-remnant pair. Its emission is intense from 100 MeV to $>$100 GeV. Its LAT spectrum can be decoupled into two components which are respectively governed by two different mechanisms. According to the young age of this system, the magnetar is seemingly a necessary and sufficient source for the downward-curved spectrum below 10 GeV, as the observed $<$10 GeV flux is too high for the supernova remnant to account for. On the other hand, the supernova remnant is reasonably responsible for the hard spectrum above 10 GeV. Further studies of this region in the TeV regime is required, so that we can perform physically meaningful comparisons of the $>$10 GeV spectrum and the TeV spectrum.

PSR J2032+4127 is a radio-loud gamma-ray-emitting pulsar; it is orbiting around a high-mass Be type star with a very long orbital period of 25-50years, and is approaching periastron, which will occur in late 2017/early 2018. This system comprises with a young pulsar and a Be type star, which is similar to the so-called gamma-ray binary PSR~B1259-63/LS2883. It is expected therefore that PSR J2032+4127 shows an enhancement of high-energy emission caused by the interaction between the pulsar wind and Be wind/disk around periastron. Ho et al. recently reported a rapid increase in the X-ray flux from this system. In this paper, we also confirm a rapid increase in the X-ray flux along the orbit, while the GeV flux shows no significant change. We discuss the high-energy emissions from the shock caused by the pulsar wind and stellar wind interaction and examine the properties of the pulsar wind in this binary system. We argue that the rate of increase of the X-ray flux observed by Swift indicates (1) a variation of the momentum ratio of the two-wind interaction region along the orbit, or (2) an evolution of the magnetization parameter of the pulsar wind with the radial distance from the pulsar. We also discuss the pulsar wind/Be disk interaction at the periastron passage, and propose the possibility of formation of an accretion disk around the pulsar. We model high-energy emissions through the inverse-Compton scattering process of the cold-relativistic pulsar wind off soft photons from the accretion disk.

We investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, applying the pulsar outer-gap model to black hole magnetospheres. During a low accretion phase, the radiatively inefficient accretion flow (RIAF) cannot emit enough MeV photons that are needed to sustain the force-free magnetosphere via two-photon collisions. In such a charge-starved region (or a gap), an electric field arises along the magnetic field lines to accelerate electrons and positrons into ultra-relativistic energies. These relativistic leptons emit copious gamma-rays via curvature and inverse-Compton (IC) processes. Some of such gamma-rays collide with the submillimeter-IR photons emitted from the RIAF to materialize as pairs, which polarize to partially screen the original acceleration electric field. It is found that the gap gamma-ray luminosity increases with decreasing accretion rate. However, if the accretion rate decreases too much, the diminished RIAF soft photon field can no longer sustain a stationary pair production within the gap. As long as a stationary gap is formed, the magnetosphere becomes force-free outside the gap by the cascaded pairs, irrespective of the BH mass. If a nearby stellar-mass black hole (BH) is in quiescence, or if a galactic intermediate-mass BH is in a very low accretion state, its curvature and IC emissions are found to be detectable with Fermi/LAT and imaging atmospheric Cherenkov telescopes (IACT). If a low-luminosity active galactic nucleus is located within a few tens of Mpc, the IC emission from its super-massive BH is marginally detectable with IACT.

The region around SGR 1806-20 and its host stellar cluster Cl* 1806-20 is a potentially important site of particle acceleration. The soft $\gamma-$ray repeater and Cl* 1806-20, which also contains several very massive stars including a luminous blue variable hypergiant LBV 1806-20, are capable of depositing a large amount of energy to the surroundings. Using the data taken with the \emph{Fermi} Large Area Telescope (LAT), we identified an extended LAT source to the south-west of Cl* 1806-20. The centroid of the 1-50~GeV emission is consistent with that of HESS J1808-204 (until now unidentified). The LAT spectrum is best-fit by a broken power-law with the break energy $E_\mathrm{b}=297\pm15$ MeV. The index above $E_\mathrm{b}$ is $2.60\pm0.04$, and is consistent with the flux and spectral index above 100 GeV for HESS J1808-204, suggesting an association between the two sources. Meanwhile, the interacting supernova remnant SNR G9.7-0.0 is also a potential contributor to the LAT flux. A tentative flux enhancement at the MeV band during a 45-day interval (2011 Jan 21 - 2011 Mar 7) is also reported. We discuss possible origins of the extended LAT source in the context of both leptonic and hadronic scenarios.

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 $\gamma$-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 $\gamma$-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.

Over the last few years, the data obtained using the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope has provided new insights on high-energy processes in globular clusters, particularly those involving compact objects such as Millisecond Pulsars (MSPs). Gamma-ray emission in the 100 MeV to 10 GeV range has been detected from more than a dozen globular clusters in our galaxy, including 47 Tucanae and Terzan 5. Based on a sample of known gamma-ray globular clusters, the empirical relations between gamma-ray luminosity and properties of globular clusters such as their stellar encounter rate, metallicity, and possible optical and infrared photon energy densities, have been derived. The measured gamma-ray spectra are generally described by a power law with a cut-off at a few gigaelectronvolts. Together with the detection of pulsed gamma-rays from two MSPs in two different globular clusters, such spectral signature lends support to the hypothesis that gamma-rays from globular clusters represent collective curvature emission from magnetospheres of MSPs in the clusters. Alternative models, involving Inverse-Compton (IC) emission of relativistic electrons that are accelerated close to MSPs or pulsar wind nebula shocks, have also been suggested. Observations at >100 GeV by using Fermi/LAT and atmospheric Cherenkov telescopes such as H.E.S.S.-II, MAGIC-II, VERITAS, and CTA will help to settle some questions unanswered by current data.

We report the results from a detailed $\gamma-$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 $\gamma-$ray spectrum can be modeled with a single power-law with a photon index of $\Gamma\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 $\gamma-$ray pulsar. This makes it possibly associated with PSR B1737-20 or PSR J1739-3023.

Gamma-ray bursts (GRBs) are some of the most energetic and exotic events in the Universe, however their behaviour at the highest energies (>10 GeV) is largely unknown. Although the Fermi-LAT space telescope has detected several GRBs in this energy range, it is limited by the relatively small collection area of the instrument. The H.E.S.S. experiment has now entered its second phase by adding a fifth telescope of 600 m$^{2}$ mirror area to the centre of the array. This new telescope increases the energy range of the array, allowing it to probe the sub-100 GeV range while maintaining the large collection area of ground based gamma-ray observatories, essential to probing short-term variability at these energies. We will present a description of the GRB observation scheme used by the H.E.S.S. experiment, summarising the behaviour and performance of the rapid GRB repointing system, the conditions under which potential GRB repointings are made and the data analysis scheme used for these observations.

Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $\gamma-$ray sources in this catalog are identified as promising MSP candidates based on their $\gamma$-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 $\gamma$-ray fluxes. The X-ray to $\gamma$-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.

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

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.

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.

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.

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\sigma$). The detection of pulsation is significant above $4\sigma$ at $>79$ GeV and above $3\sigma$ 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.

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.

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.

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

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.

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 $\Gamma\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.

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

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.

GRB 131231A was detected by the Large Area Telescope onboard Fermi Space Gamma-ray Telescope. The high energy gamma-ray ($> 100$ MeV) afterglow emission spectrum is $F_\nu \propto \nu^{-0.54\pm0.15}$ in the first $\sim 1300$ s after the trigger and the most energetic photon has an energy $\sim 62$ GeV arriving at $t\sim 520$ s. With reasonable parameters of the GRB outflow as well as the density of the circum-burst medium, the synchrotron radiation of electrons or protons accelerated at an external forward shock have difficulty accounting for the data. The synchrotron self-Compton radiation of the forward shock-accelerated electrons, instead, can account for both the spectrum and temporal behavior of the GeV afterglow emission. We also show that the prospect for detecting GRB 131231A$-$like GRBs with Cherenkov Telescope Array (CTA) is promising.

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.

The extended high-energy gamma-ray (>100 MeV) emission occurring after the prompt gamma-ray bursts (GRBs) is usually characterized by a single power-law spectrum, which has been explained as the afterglow synchrotron radiation. We report on the Fermi Large Area Telescope (LAT) observations of the >100 MeV emission from the very bright and nearby GRB 130427A, up to ~100 GeV. By performing time-resolved spectral fits of GRB 130427A, we found a strong evidence of an extra hard spectral component above a few GeV that exists in the extended high-energy emission of this GRB. This extra spectral component may represent the first clear evidence of the long sought-after afterglow inverse Compton emission. Prospects for observations at the very high-energy gamma-rays, i.e., above 100 GeV, are described.

We have investigated the X-ray and $\gamma$-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 $\gamma-$ray emission in $0.1-300$ GeV from this system at a significance of $\sim6\sigma$ for the first time. The $\gamma-$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.

GRB 130821A was detected by Fermi-GBM/LAT, Konus-Wind, SPI-ACS/INTEGRAL, RHESSI and Mars Odyssey-HEND. Although the data of GRB 130821A are very limited, we show in this work that the high energy gamma-ray emission (i.e., above 100 MeV) alone imposes tight constraint on the density profile of the circum-burst medium as well as the initial Lorentz factor of the outflow. The temporal behavior of the high energy gamma-ray emission is consistent with the forward shock synchrotron radiation model and the circum-burst medium likely has a constant-density profile. The Lorentz factor is about a few hundred, similar to other bright GRBs.

The very high energy (VHE; >100 GeV) source HESS J0632+057 has been recently confirmed as a \gamma-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 \gamma-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 \mu 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 \gamma-ray emission.

A nearby super-luminous burst GRB 130427A was simultaneously detected by six $\gamma$-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 $\gamma-$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 $\gamma-$rays at energies above 10 GeV, while five were detected in the same time interval. A possible unified model for the prompt soft $\gamma-$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.

H.E.S.S. is an array of five Imaging Atmospheric Cherenkov Telescopes (IACTs) located 1800 m above sea level in the Khomas Highland of Namibia and is sensitive to very-high-energy (VHE) gamma rays between tens of GeV to tens of TeV. The very-high background rejection capabilities of IACTs provide excellent sensitivity of H.E.S.S. to GRBs. In this contribution the status of the H.E.S.S. GRB programme, already started in 2003, is reviewed. A highlight is the recent addition of the fifth telescope, which is the world's largest IACT. Its 600 square metre mirror lowers the energy threshold to tens of GeV and provides an effective area that is ten thousands of times larger than Fermi-LAT at these energies. The higher performance drive system will reduce the response time to a GRB alert, which will significantly enhance the chances of a H.E.S.S. GRB detection. Recent results on selected GRBs will be shown.

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

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.

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.

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 {\gamma}-ray pulsations above 200 MeV from the AXP, 1E 2259+586, which puts the current theoretical models of {\gamma}-ray emission mechanisms of magnetars into challenge. We speculate that the high-energy {\gamma}-rays originate from the outer magnetosphere of the magnetar.

Gamma-ray bursts (GRBs) that emit photons at GeV energies form a small but significant population of GRBs. However, the number of GRBs whose GeV-emitting period is simultaneously observed in X-rays remains small. We report gamma-ray observations of GRB 110625A using Fermi's Large Area Telescope (LAT) in the energy range 100 MeV to 20 GeV. Gamma-ray emission at these energies was clearly detected using data taken between 180s and 580s after the burst, an epoch after the prompt emission phase. The GeV light curve differs from a simple power-law decay, and probably consists of two emission periods. Simultaneous Swift/XRT observations did not show flaring behaviors as in the case of GRB 100728A. We discuss the possibility that the GeV emission is the synchrotron self-Compton radiation of underlying ultraviolet flares.

We report the discovery of pulsed {\gamma}-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 {\gamma}-ray with a signicance of \sim 8{\sigma} in the 0.2 - 20 GeV band. The 200 MeV - 20 GeV {\gamma}-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{\sigma}) {\gamma}-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 {\gamma}-rays and X-rays with theoretical models, we suggest that the {\gamma}-ray emission is from the outer gap while the X-rays can be from intra-binary shock and pulsar magnetospheric synchrotron emissions.

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.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

We report on the discovery of gamma-ray emission from several globular clusters (GCs), including Terzan 5, the second known gamma-ray GCs. By now, more than a dozen GCs are known to emit gamma-rays of energies above 100 MeV, thus enabling us to carry out the first detailed correlation study with several cluster properties. We found strong correlations between the observed gamma-ray luminosities and four cluster parameters: stellar encounter rate, metallicity [Fe/H], and energy densities of the soft photons at the cluster locations. These "fundamental planes" of gamma-ray GCs put an intimate relation of the observed gamma-rays to the underlying millisecond pulsar population and have important implications on the origin of the gamma-ray emission of GCs.

To understand Galactic objects that emits GeV-TeV emission, a spatial correlation study between the Fermi bright source catalog and TeV source population was carried out in Tam et al. (2010), finding that a significant number of very high-energy (VHE; E >100 GeV) sources are also emitting at GeV energies. We extended our previous study utilizing the first Fermi catalog (1FGL) sources. A cross-correlation comparison of the 1FGL sources was carried out with the VHE gamma-ray sources in the literature as of May 2011. While it is found that a significant number of VHE gamma-ray sources were also detected in the GeV band, the GeV-TeV spectra of some of these spatially coincident sources cannot be described by a single spectral component. While some of these cases are gamma-ray pulsars accompanied by VHE gamma-ray emitting nebulae, we present cases where the 100 MeV to multi-TeV spectra of coincident 1FGL/VHE source pairs do not seem to be well fit by a single spectral component.

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.

The binary system PSR B1259-63/LS 2883 consists of a 47.8 ms radio pulsar that orbits the companion Be star with a period of 3.4 years in a highly eccentric orbit. The system is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital-phase-dependent variability in all observed frequencies. Here we report on the discovery of >100 MeV gamma-rays from PSR B1259-63/LS 2883 through the 2010 pariastron passage. Using data collected with the Large Area Telescope aboard Fermi from 33 days before pariastron to 75 days after pariastron, PSR B1259-63/LS 2883 is detected at a significance of 13.6 standard deviations. The gamma-ray light curve is highly variable over the above period, with changing photon index that correlates with gamma-ray flux. In particular, two major flares that occur after the pariastron passage were observed. The onset of gamma-ray emission occurs close to, but not at the same orbital phases as, the two disk passages that occur ~1 month before and ~1 month after the pariastron passage. The fact that the GeV orbital light curve is different from that of the X-ray and TeV light curves strongly suggests that GeV gamma-ray emission originates from a different component. We speculate that the observed GeV flares may be resulting from Doppler boosting effects.

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

The gamma-ray pulsar PSR B1706-44 and the adjacent supernova remnant (SNR) candidate G343.1-2.3 were observed by H.E.S.S. during a dedicated observation campaign in 2007. As a result of this observation campaign, a new source of very-high-energy (VHE; E > 100 GeV) gamma-ray emission, HESS J1708-443, was detected with a statistical significance of 7 sigma, although no significant point-like emission was detected at the position of the energetic pulsar itself. In this paper, the morphological and spectral analyses of the newly-discovered TeV source are presented. The centroid of HESS J1708-443 is considerably offset from the pulsar and located near the apparent center of the SNR, at RA(J2000) = 17h08m11s +/- 17s and Dec(J2000) = -44d20' +/- 4'. The source is found to be significantly more extended than the H.E.S.S. point spread function (~0.1 deg), with an intrinsic Gaussian width of 0.29 deg +/- 0.04 deg. Its integral flux between 1 and 10 TeV is ~ 3.8 x 10^-12 ph cm^-2 s^-1, equivalent to 17% of the Crab Nebula flux in the same energy range. The measured energy spectrum is well-fit by a power law with a relatively hard photon index Gamma = 2.0 +/- 0.1 (stat) +/- 0.2 (sys). Additional multi-wavelength data, including 330 MHz VLA observations, were used to investigate the VHE gamma-ray source's possible associations with the pulsar wind nebula of PSR B1706-44 and/or with the complex radio structure of the partial shell-type SNR G343.1-2.3.

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.

We have investigated the properties of a group of $\gamma$-ray emitting globular clusters (GCs) which have recently been uncovered in our Galaxy. By correlating the observed $\gamma$-ray luminosities $L_{\gamma}$ with various cluster properties, we probe the origin of the high energy photons from these GCs. We report $L_{\gamma}$ is positively correlated with the encounter rate $\Gamma_{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_{\gamma}$ increase with the energy densities of the soft photon at the cluster location. Furthermore, the two-dimensional regression analysis suggests that $L_{\gamma}$, soft photon densities, and $\Gamma_{c}$/$\left[{\rm Fe/H}\right]$ possibly span fundamental planes which potentially provide better predictions for the $\gamma$-ray properties of GCs.

The low-mass X-ray binary (LMXB) system FIRST J102347.6+003841 hosts a newly born millisecond pulsar (MSP) PSR J1023+0038 that was revealed as the first and only known rotation-powered MSP in a quiescent LMXB. While the system is shown to have an accretion disk before 2002, it remains unclear how the accretion disk has been removed in order to reveal the radio pulsation in 2007. In this Letter, we report the discovery of gamma-rays spatially consistent with FIRST J102347.6+003841, at a significance of 7 standard deviations, using data obtained by the Fermi Gamma-ray Space Telescope. The gamma-ray spectrum can be described by a power law (PL) with a photon index of 2.9+-0.2, resulting in an energy flux above 200 MeV of (5.5+-0.9)x10^{-12} erg cm^{-2}s^{-1}. The gamma-rays likely originate from the MSP PSR J1023+0038, but also possibly from an intrabinary shock between the pulsar and its companion star. To complement the gamma-ray study, we also re-investigate the XMM-Newton data taken in 2004 and 2008. Our X-ray spectral analysis suggests that a broken PL with two distinct photon indices describes the X-ray data significantly better than a single PL. This indicates that there exists two components and that both components appear to vary with the orbital phase. The evidence for gamma-ray emission conforms with a recent suggestion that gamma-rays from PSR J1023+0038 may be responsible for ejecting the disk material out of the system.

Context: The detection of gamma-rays in the very-high-energy (VHE) range (100 GeV-100 TeV) offers the possibility of studying the parent population of ultrarelativistic particles found in astrophysical sources, so it is useful for understanding the underlying astrophysical processes in nonthermal sources. Aim: The discovery of the VHE gamma-ray source HESS J1507-622 is reported and possibilities regarding its nature are investigated. Methods: The H.E.S.S. array of imaging atmospheric Cherenkov telescopes (IACTs) has a high sensitivity compared with previous instruments (~1% of the Crab flux in 25 hours observation time for a 5 sigma point-source detection) and has a large field of view (~5 deg in diameter). HESS J1507-622 was discovered within the ongoing H.E.S.S. survey of the inner Galaxy, and the source was also studied by means of dedicated multiwavelength observations. Results: A Galactic gamma-ray source, HESS J1507-622, located ~3.5 deg from the Galactic plane was detected with a statistical significance > 9 sigma. Its energy spectrum is well fitted by a power law with spectral index \Gamma = 2.24 +/- 0.16_{stat} +/- 0.20_{sys} and a flux above 1 TeV of (1.5 +/- 0.4_{stat} +/- 0.3_{sys}) X 10^{-12} cm^{-2} s^{-1}. Possible interpretations (considering both hadronic and leptonic models) of the VHE gamma-ray emission are discussed in the absence of an obvious counterpart.

Aims. Previous observations with the H.E.S.S. telescope array revealed the existence of extended very-high-energy (VHE; E>100 GeV) {\gamma}-ray emission, HESS J1023-575, coincident with the young stellar cluster Westerlund 2. At the time of discovery, the origin of the observed emission was not unambiguously identified, and follow-up observations have been performed to further investigate the nature of this {\gamma}-ray source. Methods. The Carina region towards the open cluster Westerlund 2 has been re-observed, increasing the total exposure to 45.9 h. The combined dataset includes 33 h of new data and now permits a search for energy-dependent morphology and detailed spectroscopy. Results. A new, hard spectrum VHE {\gamma}-ray source, HESSJ1026-582, was discovered with a statistical significance of 7{\sigma}. It is positionally coincident with the Fermi LAT pulsar PSR J1028-5819. The positional coincidence and radio/{\gamma}-ray characteristics of the LAT pulsar favors a scenario where the TeV emission originates from a pulsar wind nebula. The nature of HESS J1023-575 is discussed in light of the deep H.E.S.S. observations and recent multi-wavelength discoveries, including the Fermi LAT pulsar PSRJ1022-5746 and giant molecular clouds in the region. Despite the improved VHE dataset, a clear identification of the object responsible for the VHE emission from HESS J1023-575 is not yet possible, and contribution from the nearby high-energy pulsar and/or the open cluster remains a possibility.

PKS 0548-322 (z=0.069) is a ``high-frequency-peaked'' BL Lac object and a candidate very high energy (VHE, E>100 GeV) gamma-ray emitter, due to its high X-ray and radio flux. Observations at the VHE band provide insights into the origin of very energetic particles present in this source and the radiation processes at work. We report observations made between October 2004 and January 2008 with the H.E.S.S. array, a four imaging atmospheric-Cherenkov telescopes. Contemporaneous UV and X-ray observations with the Swift satellite in November 2006 are also reported. PKS 0548-322 is detected for the first time in the VHE band with H.E.S.S. We measure an excess of 216 gamma-rays corresponding to a significance of 5.6 standard deviations. The photon spectrum of the source is described by a power-law, with a photon index of Gamma=2.86 +/- 0.34 (stat) +/- 0.10 (sys). The integral flux above 200 GeV is 1.3 % of the flux of the Crab Nebula, and is consistent with being constant in time. Contemporaneous Swift/XRT observations reveal an X-ray flux between 2 and 10 keV of F_{2-10 keV}=2.3 +/- 0.2 x 10^{-11} erg.cm^{-2}. s^{-1}, an intermediate intensity state with respect to previous observations. The spectral energy distribution can be reproduced using a simple one-zone synchrotron self Compton model, with parameters similar those observed for other sources of this type.

Observations of very high energy gamma-rays from blazars provide information about acceleration mechanisms occurring in their innermost regions. Studies of variability in these objects allow a better understanding of the mechanisms at play. To investigate the spectral and temporal variability of VHE (>100 GeV) gamma-rays of the well-known high-frequency-peaked BL Lac object PKS 2155-304 with the H.E.S.S. imaging atmospheric Cherenkov telescopes over a wide range of flux states. Data collected from 2005 to 2007 are analyzed. Spectra are derived on time scales ranging from 3 years to 4 minutes. Light curve variability is studied through doubling timescales and structure functions, and is compared with red noise process simulations. The source is found to be in a low state from 2005 to 2007, except for a set of exceptional flares which occurred in July 2006. The quiescent state of the source is characterized by an associated mean flux level of 4.32 +/-0.09 x 10^-11 cm^-2 s^-1 above 200 GeV, or approximately 15% of the Crab Nebula, and a power law photon index of 3.53 +/-0.06. During the flares of July 2006, doubling timescales of ~2 min are found. The spectral index variation is examined over two orders of magnitude in flux, yielding different behaviour at low and high fluxes,which is a new phenomenon in VHE gamma-ray emitting blazars. The variability amplitude characterized by the fractional r.m.s. is strongly energy-dependent and is proportional to E^(0.19 +/- 0.01). The light curve r.m.s. correlates with the flux. This is the signature of a multiplicative process which can be accounted for as a red noise with a Fourier index of ~2. This unique data set shows evidence for a low level gamma-ray emission state from PKS 2155-304, which possibly has a different origin than the outbursts. The discovery of the light curve lognormal behaviour might be an indicator ..

AIMS: The properties of the broad-band emission from the high-frequency peaked BL Lac H 2356-309 (z=0.165) are investigated. METHODS: Very High Energy (VHE; E > 100 GeV) observations of H 2356-309 were performed with the High Energy Stereoscopic System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM optical monitoring observations were also made in 2007. RESULTS: A strong VHE signal, ~13 sigma total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The integral flux above 240 GeV is I(>240 GeV) = (3.06 +- 0.26 {stat} +- 0.61 {syst}) x 10^{-12} cm^{-2} s^{-1}, corresponding to ~1.6% of the flux observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a power law (photon index of Gamma = 3.06 +- 0.15 {stat} +- 0.10 {syst}). Significant small-amplitude variations in the VHE flux from H 2356-309 are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by a hard, broken-power-law spectrum on both nights. CONCLUSIONS: The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed X-ray flux is at a historically low level.

Recent theoretical predictions of the lowest very high energy (VHE) luminosity of SN 1006 are only a factor 5 below the previously published H.E.S.S. upper limit, thus motivating further in-depth observations of this source. Deep observations at VHE energies (above 100 GeV) were carried out with the High Energy Stereoscopic System (H.E.S.S.) of Cherenkov Telescopes from 2003 to 2008. More than 100 hours of data have been collected and subjected to an improved analysis procedure. Observations resulted in the detection of VHE gamma-rays from SN 1006. The measured gamma-ray spectrum is compatible with a power-law, the flux is of the order of 1% of that detected from the Crab Nebula, and is thus consistent with the previously established H.E.S.S. upper limit. The source exhibits a bipolar morphology, which is strongly correlated with non-thermal X-rays. Because the thickness of the VHE-shell is compatible with emission from a thin rim, particle acceleration in shock waves is likely to be the origin of the gamma-ray signal. The measured flux level can be accounted for by inverse Compton emission, but a mixed scenario that includes leptonic and hadronic components and takes into account the ambient matter density inferred from observations also leads to a satisfactory description of the multi-wavelength spectrum.

Very high-energy (VHE; E>100 GeV) gamma-rays have been detected from a wide range of astronomical objects, such as pulsar wind nebulae (PWNe), supernova remnants (SNRs), giant molecular clouds, gamma-ray binaries, the Galactic Center, active galactic nuclei (AGN), radio galaxies, starburst galaxies, and possibly star-forming regions as well. At lower energies, observations using the Large Area Telescope (LAT) onboard Fermi provide a rich set of data which can be used to study the behavior of cosmic accelerators in the MeV to TeV energy bands. In particular, the improved angular resolution of current telescopes in both bands compared to previous instruments significantly reduces source confusion and facilitates the identification of associated counterparts at lower energies. In this paper, a comprehensive search for VHE gamma-ray sources which are spatially coincident with Galactic Fermi/LAT bright sources is performed, and the available MeV to TeV spectra of coincident sources are compared. It is found that bright LAT GeV sources are correlated with TeV sources, in contrast to previous studies using EGRET data. Moreover, a single spectral component seems unable to describe the MeV to TeV spectra of many coincident GeV/TeV sources. It has been suggested that gamma-ray pulsars may be accompanied by VHE gamma-ray emitting nebulae, a hypothesis that can be tested with VHE observations of these pulsars.

Very high-energy (VHE; E > 100 GeV) gamma-rays have been detected from a wide range of astronomical objects, such as SNRs, pulsars and pulsar wind nebulae, AGN, gamma-ray binaries, molecular clouds, and possibly star-forming regions as well. At lower energies, sources detected using Large Area Telescope (LAT) aboard Fermi provide a rich set of data which can be used to study the behavior of cosmic accelerators in the GeV to TeV energy bands. In particular, the improved angular resolution in both bands compared to previous instruments significantly reduces source confusion and facilitates the identification of associated counterparts at lower energies. In this talk, a comprehensive search for VHE gamma-ray sources which are spatially coincident with Galactic Fermi/LAT bright sources is performed, and the GeV to TeV spectra of selected coincident sources are shown. It is found that LAT bright GeV sources are correlated to TeV sources, in contrast with previous studies using EGRET data.

The synchrotron self-Compton (SSC) emission from Gamma-ray Burst (GRB) forward shock can extend to the very-high-energy (VHE; $E_\gamma > $100 GeV) range. Such high energy photons are rare and are attenuated by the cosmic infrared background before reaching us. In this work, we discuss the prospect to detect these VHE photons using the current ground-based Cherenkov detectors. Our calculated results are consistent with the upper limits obtained with several Cherenkov detectors for GRB 030329, GRB 050509B, and GRB 060505 during the afterglow phase. For 5 bursts in our nearby GRB sample (except for GRB 030329), current ground-based Cherenkov detectors would not be expected to detect the modeled VHE signal. Only for those very bright and nearby bursts like GRB 030329, detection of VHE photons is possible under favorable observing conditions and a delayed observation time of $\la$10 hours.

The High Energy Stereoscopic System (H.E.S.S.) consists of four Imaging Atmospheric Cherenkov Telescopes (IACTs) in Namibia for the detection of cosmic very-high-energy (VHE) gamma-rays. Gamma-ray bursts (GRBs) are among the potential VHE gamma-ray sources. VHE gamma-emission from GRBs is predicted by many GRB models. Because of its generally fast-fading nature in many wavebands, the time evolution of any VHE gamma-radiation is still unknown. In order to probe the largely unexplored VHE gamma-ray spectra of GRBs, a GRB observing program has been set up by the H.E.S.S. collaboration. With the high sensitivity of the H.E.S.S. array, VHE gamma-ray flux levels predicted by GRB models are well within reach. We report the H.E.S.S. observations of and results from some of the reported GRB positions during March 2003 - May 2006.

Gamma-ray Bursts (GRBs) are among the potential extragalactic sources of very-high-energy (VHE) gamma-rays. We discuss the prospects of detecting VHE gamma-rays with current ground-based Cherenkov instruments during the afterglow phase. Using the fireball model, we calculate the synchrotron self-Compton (SSC) emission from forward-shock electrons. The modeled results are compared with the observational afterglow data taken with and/or the sensitivity level of ground-based VHE instruments (e.g. STACEE, H.E.S.S., MAGIC, VERITAS, and Whipple). We find that modeled SSC emission from bright and nearby bursts such as GRB 030329 are detectable by these instruments even with a delayed observation time of ~10 hours.

H.E.S.S. (High Energy Stereoscopic System), which is designed to detect TeV gamma-rays, is a system of four Imaging Atmospheric Cherenkov Telescopes situated in Namibia. The system has been shown to be very successful in detecting and observing galactic and extra-galactic TeV sources. In order to explore the highest energy end of GRB spectra, a GRB observing program has been established in the H.E.S.S. collaboration. Here we introduce our GRB observing program and report on its current status.

The association of long gamma-ray bursts (GRBs) with star forming regions and the idea of massive stars as progenitors of GRBs are widely accepted. Because of their short lifetimes, it is very likely that massive stars are still embedded in dense molecular clouds when they give birth to GRBs. Stellar winds from GRB progenitors can create low-density bubbles with sizes and densities strongly depending on the initial ambient density. A boundary between the bubble and the dense molecular cloud must exist with the density at the boundary increasing from that of the bubble to that of the outer cloud. We have calculated the lightcurves of the afterglows in such environments with three regions: the stellar wind region, the boundary, and the molecular cloud. We show that the interaction between the cylindrical jet and the density boundary can result in a re-brightening of the afterglow occurring as early as ~1 day after the GRB. We compare our models with the optical afterglows of GRB970508, GRB000301C, and GRB030226. We find that the values of our model parameters, including the radius of the wind bubble, the densities in the bubble and in the outer molecular cloud are within typical ranges.