High-redshift blazars are valuable tools to study the early Universe. So far only a handful of $\gamma$-ray blazars have been found at redshifts above 3. $\gamma$-ray signals are detected in the direction of PMN J2219-2719 ($z=3.63$) and PMN J2321-0827 ($z=3.16$) by analyzing the 10-year $Fermi$-LAT Pass 8 data. PMN J2219-2719 is not distinguished from the background in the global analysis. During the 5-month epoch, the TS value is 47.8 and the flux is more than 10 times of the 10-year averaged flux. In addition, the angular distance between the $\gamma$-ray position and the radio position of PMN J2219-2719 is only ${0.04}^{\circ}$. Moreover, the $\gamma$-ray and infrared light curves of long time scale are very similar, which support the association between the $\gamma$-ray source and PMN J2219-2719. The global analysis of PMN J2321-0827 suggest a new $\gamma$-ray source, during the flare phase, the TS value is 61.4 and the $\gamma$-ray flux increased significantly. The association probability suggests that PMN J2321-0827 may be the counterpart of the new $\gamma$-ray source. In the future, the number of high-redshift $\gamma$-ray sources will increase by combining $Fermi$-LAT and the upcoming Large Synoptic Survey Telescope.

GB 1508+5714 is a high-redshift blazar (z = 4.3) and a spectrally soft gamma-ray source has been found in its direction. By analyzing 11.4-yr Fermi-LAT data, significant long-term variability of the gamma-ray source is confirmed. More importantly, a gamma-ray emission flare appeared in an epoch of several tens of days in year 2018, when the flux is about four times of the value from the global fit. Meanwhile, optical flares displayed in both i- and r-bands from the Zwicky Transient Facility light curves. Detections of the simultaneous gamma-ray and optical brightening provide a decisive proof to pin down the association between the gamma-ray source and GB 1508+5714, which makes it the first identified gamma-ray blazar beyond redshift 4. Future multiwavelength campaigns are needed to further investigate its jet properties, and, upcoming wide-deep-fast optical telescopes together with the gamma-ray sky surveyors will shed lights on the role that the AGN jets play in the early cosmic time.

For on-axis typical short gamma-ray bursts (sGRBs), the forward shock emission is usually so bright that renders the identification of kilonovae (also known as macronovae) in the early afterglow ($t<0.5$ d) phase rather challenging. This is why previously no thermal-like kilonova component has been identified at such early time except in the off-axis dim GRB 170817A associated with GW170817. Here we report the identification of an unusual optical radiation component in GRB 070809 at $t\sim 0.47$ d, thanks plausibly to the very-weak/subdominant forward shock emission. The optical emission with a very red spectrum is well in excess of the extrapolation of the X-ray emission that is distinguished by an unusually hard spectrum, which is at odds with the forward shock afterglow prediction but can be naturally interpreted as a kilonova. Our finding supports the speculation that kilonovae are ubiquitous , and demonstrates the possibility of revealing the neutron star merger origin with the early afterglow data of some typical sGRBs that take place well beyond the sensitive radius of the advanced gravitational wave detectors and hence the opportunity of organizing dedicated follow-up observations for events of interest.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

Circinus galaxy is a nearby composite starburst/AGN system. In this work we re-analyze the GeV emission from Circinus with 10 years of {\it Fermi}-LAT Pass 8 data. In the energy range of 1-500 GeV, the spectrum can be well fitted by a power-law model with a photon index of $\Gamma$ = $2.20\pm0.14$, and its photon flux is $(5.90\pm1.04) \times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$. Our 0.1-500 GeV flux is several times lower than that reported in the previous literature, which is roughly in compliance with the empirical relation for star-forming and local group galaxies and might be reproduced by the interaction between cosmic rays and the interstellar medium. The ratio between the $\gamma$-ray luminosity and the total infrared luminosity is near the proton calorimetric limit, indicating that Circinus may be a proton calorimeter. However, marginal evidence for variability of the $\gamma$-ray emission is found in the timing analysis, which may indicate the activity of AGN jet. More {\it Fermi}-LAT data and future observation of CTA are required to fully reveal the origin of its $\gamma$-ray emission.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

We report on OVRO, WISE, Swift and Fermi-LAT observations of the high redshift blazar CGRaBS J0733+0456, from which significant flux variations in radio, infrared (IR) as well as $\gamma$-ray domains are detected. Particularly, the amplitude of long-term IR variation is over one order of magnitude and the IR variability timescale can be constrained as short as a few hours in the source frame. The IR and $\gamma$-ray light curves are found to be rather similar, and the strong quasi-simultaneous infrared and $\gamma$-ray flares are proved to be unique among the nearby $\gamma$-ray sources. This is the first time to identify a $\gamma-$ray blazar at redshift $z\geq 3$ with multi-wavelength flux variations (flares). Broadband spectral energy distributions in different flux states are constructed and theoretically described. The $\gamma-$ray flares from some blazars as distant as redshift $\sim$ 5 are expected to be detectable for Fermi-LAT.

HESS J1640-465 is an extended TeV $\gamma$-ray source and its $\gamma$-ray emission whether from the shell of a supernova remnant (SNR) or a pulsar wind nebula (PWN) is still under debate. We reanalyze the GeV $\gamma$-ray data in the field of HESS J1640-465 using eight years of Pass 8 data recorded by the Fermi Large Area Telescope. An extended GeV $\gamma$-ray source positionally coincident with HESS J1640-465 is found. Its photon spectrum can be described by a power-law with an index of $1.42\pm0.19$ in the energy range of 10-500 GeV, and smoothly connects with the TeV spectrum of HESS J1640-465. The broadband spectrum of HESS J1640-465 can be well fit by a leptonic model with a broken power-law spectrum of electrons with an exponential cut-off at $\sim$ 300 TeV. The spectral properties of HESS J1640-465 are broadly consistent with the characteristics of other sources identified as PWNe, such as the correlations between high-energy luminosity ratios and the physical parameters of pulsar, including spin-down luminosity $\dot{E}$ and characteristic age $\tau_c$. All these pieces of evidence support that the $\gamma$-ray emission of HESS J1640-465 may originate from the PWN powered by PSR J1640-4631 rather than the shell of the SNR G338.3-0.0.

We report the detection of a transient $\gamma$-ray source in the direction of B3 1428+422 ($z=4.72$) by analyzing the 110-month {\it Fermi}-LAT Pass 8 data. The new transient $\gamma-$ray source is far away from the Galactic plane and has a rather soft spectrum, in agreement with being a high redshift blazar. We suggest that the newly discovered transient is the $\gamma$-ray counterpart of B3 1428+422, which could be the {\it most distant} GeV source detected so far. The detection of a group of such distant $\gamma-$ray blazars will be helpful to reconstruct the evolution of the luminosity function and to study the extragalactic background light at such high redshifts.

Recently the IceCube collaboration and 15 other collaborations reported the spatial and temporal coincidence between the neutrino event IceCube-170922A and the radio-TeV activity of the blazar TXS 0506+056. Their further analysis on 9.5 years of IceCube data discovered neutrino flare between September 2014 and March 2015, when TXS 0506+056 is however in "quiescent" state. We analyze the Fermi-LAT data in that direction, and find another bright GeV source PKS 0502+049, which is at an angle of $1.2^{\circ}$ from TXS 0506+056, with strong activties during the neutrino flare. No other bright GeV source was detected in the region of interest. Though PKS 0502+049 is $1.2^\circ$ separated from TXS 0506+056, it locates within the directional reconstruction uncertainties of 7 neutrinos, out of the 13 neutrino events during the neutrino flare. Together with the observed high flux of the $\gamma$-ray flare, it may be unreasonable to fully discard the (partial) contribution of PKS 0502+049 to the neutrino flare. The single source assumption used in the neutrino data analysis might need to be revisited.

On 2017 September 22 the IceCube neutrino observatory detected a track-like, very-high-energy event (IceCube-170922A) that is spatially associated with TXS 0506+056, a quasar at a redshift of $z=0.3365$. This source is characterized by the increased acitivies in a very wide energy range (from radio to TeV) during these days. To investigate the possible connection of the PeV neutrino emission with the GeV activity of blazars, in this work we select 116 bright sources and analyze their lightcurves and spectra. We focus on the sources displaying GeV activities. Among these blazars, TXS 0506+056 seems to be typical in many aspects but is distinguished by the very strong GeV activties. We suggest to search for neutrino outburst in the historical data of IceCube, as recently done for TXS 0506+056, from the directions of these more energetic and harder blazars with strong GeV activities.

The composition and energy dissipation in jets are two of the fundamental questions of jet physics that are not fully understood. In this paper, we attempt to constrain the composition, magnetization as well as radiative efficiency for blazar with the recently released low-frequency radio catalog of the TIFR GMRT Sky Survey at 150 MHz. The jet power estimated from the low-frequency radio emission is much lower than that derived from spectral energy distribution fitting assuming one proton per electron. Assuming the jet power estimated from low-frequency radio emission is physical, the fraction of electron/positron pairs can be constrained with $n_{\rm pairs}/n_{\rm p} \sim 10$. By comparing the power carried by magnetic field and radiation with the jet power estimated from the low-frequency radio emission, we find both relatively high magnetization parameter of $\sigma \sim 0.5$ and radiative efficiency of $\eta \sim 0.4$ in the dissipation region of blazars. These results suggest that the magnetic reconnection processes may play an important role in the energy dissipation of blazars. We also explore the connection between these three parameters ($n_{\rm pairs}/n_{\rm p}$, $\sigma$, and $\eta$) and the black hole mass, disk luminosity as well as Eddington ratio. No significant correlation is found, except that $\sigma$ shows possible correlation with disk luminosity.

In this work, we search for $\gamma$-ray emission from the directions of some nearby Milky Way dwarf spheroidal galaxies (dSphs) and candidates with the publicly-available Pass 8 data of Fermi-LAT. Our sample includes 12 sources with the distances $<50$ kpc. Very weak $\gamma$-ray excesses ($\sim 2\sigma$) are found in some dSphs/candidates, consistent with that reported in the previous literature. Intriguingly, the peak test statistic (TS) value of the weak emission from the direction of Reticulum II rises continually. If interpreted as dark matter (DM) annihilation, the peak TS value is 13.5 for the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ and the DM mass of $m_\chi \sim 16$ GeV. The combination of all these nearby sources yields a more significant (with local significance $> 4\sigma$) $\gamma$-ray signal.

We report simultaneous observations at different energy bands in radio, optical, UV, X-rays and $\gamma$ rays of the flat-spectrum radio-quasar 3C 345. We built the light curve of the source at different frequencies from 2008, the beginning of the \textit{Fermi} all-sky survey, to 2016, using new data and public archives. In particular we obtained several optical spectra, to study the behavior of emission lines and the continuum in different activity states and to derive the black hole mass. 3C 345 showed two flaring episodes in 2009, which occurred simultaneously in $\gamma$ ray, optical/UV and X-rays, and were later followed in radio. The source shows an inverse Compton dominated spectral energy distribution, which moved from higher to lower frequencies from the high to the low state. The reverberation of emission lines during one outburst event allowed us to constrain the location of production of $\gamma$ rays very close to the broad-line region, and possibly in the jet-base. We report the observation of an increased accretion after the outburst, possibly induced by the decrease of magnetic field intensity with respect to the low state.

The axion-like particle (ALP)-photon mixing in the magnetic field around $\gamma$-ray sources or along the line-of-sight could induce oscillation between photons and ALPs, which then causes irregularities in the $\gamma$-ray spectra. In this work we try to search for such spectral irregularities in the spectrum of PKS $2155-304$ using 8.6 years of the Fermi Large Area Telescope (Fermi-LAT) data. No significant evidence for the presence of ALP-photon oscillation is obtained, and the parameter space of ALPs is constrained. The exclusion region sensitively depends on the poorly known magnetic field of host galaxy cluster of PKS $2155-304$. If the magnetic field is as high as $\sim 10~{\rm \mu G}$, the "hole"-like parameter region allowed in Ref.~\cite{Fermi} can be ruled out.

We present our photometric monitoring of a flat spectrum radio quasar (FSRQ) 3C 454.3 at Yunnan observatories from 2006 to 2011. We find that the optical color of 3C 454.3 shows obvious redder-when-brighter trend, which reaches a saturation stage when the source is brighter than 15.15 mag at V band. We perform a simulation with multiple values of disk luminosity and spectral index to reproduce the magnitude-color diagram. The results show that the contamination caused by the disk radiation alone is difficult to produce the observed color variability. The variability properties during the outburst in December 2009 are also compared with $\gamma$-ray data derived from Fermi $\gamma$-ray space telescope. The flux variation of these two bands follow a linear relation with $F_{\gamma} \propto F_R^{1.14\pm0.07}$, which provides an observational evidence for external Compton process in 3C 454.3. Meanwhile, this flux correlation indicates that electron injection is the main mechanism for variability origin. We also explore the variation of the flux ratio $F_{\gamma}/F_R$ and the detailed structures in the lightcurves, and discuss some possible origins for the detailed variability behaviors.

The remnant of supernova explosion is widely believed to be the acceleration site of high-energy cosmic ray particles. The acceleration timescale is, however, typically very long. Here we report the detection of a variable $\gamma$-ray source with the Fermi Large Area Telescope, which is positionally and temporally consistent with a peculiar supernova, iPTF14hls. A quasi-stellar object SDSS J092054.04+504251.5, which is probably a blazar according to the infrared data, is found in the error circle of the $\gamma$-ray source. More data about the $\gamma$-ray source and SDSS J092054.04+504251.5 are needed to confirm their association. On the other hand, if the association between the $\gamma$-ray source and the supernova is confirmed, this would be the first time to detect high-energy $\gamma$-ray emission from a supernova, suggesting very fast particle acceleration by supernova explosions.

We report the detection of GeV $\gamma$-ray emission from supernova remnant HESS J1731-347 using 9 years of {\it Fermi} Large Area Telescope data. We find a slightly extended GeV source in the direction of HESS J1731-347. The spectrum above 1 GeV can be fitted by a power-law with an index of $\Gamma = 1.77\pm0.14$, and the GeV spectrum connects smoothly with the TeV spectrum of HESS J1731-347. Either a hadronic-leptonic or a pure leptonic model can fit the multi-wavelength spectral energy distribution of the source. However, the hard GeV $\gamma$-ray spectrum is more naturally produced in a leptonic (inverse Compton scattering) scenario, under the framework of diffusive shock acceleration. We also searched for the GeV $\gamma$-ray emission from the nearby TeV source HESS J1729-345. No significant GeV $\gamma$-ray emission is found, and upper limits are derived.

High-redshift blazars are one of the most powerful monsters in the universe and $\gamma$-ray variability carries crucial information of their relativistic jets. In this work we present results of the first systematical temporal analysis of {\it Fermi}-LAT data of all known seven $\gamma$-ray blazars beyond redshift 3. Significant long-term $\gamma$-ray variability is found from five sources in monthly $\gamma$-ray light curves, in which three of them are reported for the first time. Furthermore, intraday $\gamma$-ray variations are detected from NVSS J053954$-$283956 and NVSS J080518$+$614423. Doubling variability timescale of the former source is limited as short as $\lesssim$ 1 hour (at the source frame). Together with variability amplitude over one order of magnitude, NVSS J053954$-$283956 is the most distant $\gamma$-ray flaring blazar so far. Meanwhile, intraday optical variability of NVSS J163547$+$362930 is found based on archival PTF/iPTF light curve. Benefited from multiwavelength activity of these sources, constraints of their Doppler factors as well as locations of $\gamma$-ray radiation region and indications for the SDSS high redshift jetted active galactic nuclei deficit are discussed.

In the narrow-line Seyfert 1 galaxy 1H 0707-495, recently a transient quasi-periodic oscillation (QPO) signal with a frequency of $\sim 2.6\times 10^{-4}$ Hz has been detected at a high statistical significance. Here, we reanalyze the same set of XMM-Newton data observed on 2008 February 4 with the Weighted-Wavelet Z-transform (WWZ) method. In addition to confirming the previous findings, we also find another QPO signal with a frequency of $\sim 1.2\times 10^{-4}$ Hz in a separated X-ray emission phase at the significance level of $\sim 3.7\sigma$. The signal is also found fitting an auto-regressive model though at a lower significance. The frequency ratio between these two signals is $\sim 2:1$. The analysis of other XMM-Newton measurements of 1H 0707-495 also reveals the presence of the $\sim 2.6\times 10^{-4}$ Hz ($\sim 1.2\times 10^{-4}$ Hz) QPO signal on 2007 May 14 (2010 September 17) at the significance level of $\sim 4.2\sigma$ ($\sim 3.5\sigma$). The QPO frequency found in this work follows the $f_{QPO}-M_{BH}$ relation reported in previous works spanning from stellar-mass to supermassive black holes. This is the first time to observe two separated transient X-ray QPO signals in active galactic nuclei (AGNs), which sheds new light on the physics of accreting supermassive black holes.

On 17 August 2017, a gravitational wave event (GW170817) and an associated short gamma-ray burst (GRB 170817A) from a binary neutron star merger had been detected. The followup optical/infrared observations also identified the macronova/kilonova emission (AT2017gfo). In this work we discuss some implications of the remarkable GW170817/GRB 170817A/AT2017gfo association. We show that the $\sim 1.7$s time delay between the gravitational wave (GW) and GRB signals imposes very tight constraint on the superluminal movement of gravitational waves (i.e., the relative departure of GW velocity from the speed of light is $\leq 4.3\times 10^{-16}$) or the possible violation of weak equivalence principle (i.e., the difference of the gamma-ray and GW trajectories in the gravitational field of the galaxy and the local universe should be within a factor of $\sim 3.4\times 10^{-9}$). The so-called Dark Matter Emulators and a class of contender models for cosmic acceleration ("Covariant Galileon") are ruled out, too. The successful identification of Lanthanide elements in the macronova/kilonova spectrum also excludes the possibility that the progenitors of GRB 170817A are a binary strange star system. The high neutron star merger rate (inferred from both the local sGRB data and the gravitational wave data) together with the significant ejected mass strongly suggest that such mergers are the prime sites of heavy r-process nucleosynthesis.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Narrow line Seyfert 1 galaxies (NLS1s) usually do not host relativistic jets and the detection of their $\gamma$-ray emissions is rarely reported. Evidences for strong relativistic beaming effect have been found in all $\gamma$-ray NLS1s to date and hence their jets are believed to be closely aligned. In this work we perform a systematic analysis of 7-year {\it Fermi} Large Area Telescope (LAT) data of 18 steep radio spectrum NLS1s (S-NLS1s) and report the {\it first} two detections of $\gamma$-ray emissions in this class, increasing the total number of known $\gamma$-ray NLS1s to 9. A strong and stable $\gamma$-ray signal (over 8$\sigma$) is produced by SDSS J144318.56+472556.7, and a significant $\gamma$-ray emission (4.5$\sigma$) of SDSS J130522.74+511640.2 has been detected in a monthly flare in November 2013. A thorough investigation of their multiwavelength properties from radio to X rays has been performed, finding no significant variability. Along with their spectral energy distributions and the radio morphology, these facts strongly suggest that these two sources harbor mis-aligned and possibly underdeveloped relativistic jets, which provides valuable targets to reveal the formation and evolution of relativistic jets under extreme circumstances.

Puppis A is a very famous and extensively studied supernova remnant (SNR) that shows strong evidence of shock-cloud interaction. We re-analyze the GeV $\gamma$-ray emission of it using seven years Pass 8 data recorded by the Fermi Large Area Telescope (Fermi-LAT). The morphology of the $\gamma$-ray emission is more compatible with that of the thermal X-ray and IR emissions than the radio image, which suggests a possible correlation between the gamma-ray emitting region and dense clouds. The $\gamma$-ray spectrum in the energy range of 1-500 GeV shows a break at $7.92\pm1.91$ GeV with the photon indices of $1.81\pm0.08$ below the break and $2.53\pm0.12$ above the break, which can naturally explain the lack of TeV $\gamma$-ray emission from Puppis A. The multi-wavelength observations favor a hadronic origin for the $\gamma$-ray emission.

We search for $\gamma$-ray and optical periodic modulations in a distant flat spectrum radio quasar (FSRQ) PKS 0426-380 (the redshift $z=1.1$). Using two techniques (i.e., the maximum likelihood optimization and the exposure-weighted aperture photometry), we obtain $\gamma$-ray light curves from \emph{Fermi}-LAT Pass 8 data covering from 2008 August to 2016 December. We then analyze the light curves with the Lomb-Scargle Periodogram (LSP) and the Weighted Wavelet Z-transform (WWZ). A $\gamma$-ray quasi-periodicity with a period of 3.35 $\pm$ 0.68 years is found at the significance-level of $\simeq3.6\ \sigma$. The optical-UV flux covering from 2005 August to 2013 April provided by ASI SCIENCE DATA CENTER is also analyzed, but no significant quasi-periodicity is found. It should be pointed out that the result of the optical-UV data could be tentative because of the incomplete of the data. Further long-term multiwavelength monitoring of this FSRQ is needed to confirm its quasi-periodicity.

GRB 160625B is an extremely-bright outburst with well-monitored afterglow emission. The geometry-corrected energy is high up to $\sim 5.2\times10^{52}$ erg or even $\sim 8\times 10^{52}$ erg, rendering it the most energetic GRB prompt emission recorded so far. We analyzed the time-resolved spectra of the prompt emission and found that in some intervals there were likely thermal-radiation components and the high energy emission were characterized by significant cutoff. The bulk Lorentz factors of the outflow material are estimated accordingly. We found out that the Lorentz factors derived in the thermal-radiation model are consistent with the luminosity-Lorentz factor correlation found in other bursts as well as in GRB 090902B for the time-resolved thermal-radiation components. While the spectral cutoff model yields much lower Lorentz factors that are in tension with the constraints set by the electron pair Compoton scattering process. We then suggest that these spectral cutoffs are more likely related to the particle acceleration process and that one should be careful in estimating the Lorentz factors if the spectrum cuts at a rather low energy (e.g., $\sim$ tens MeV). The nature of the central engine has also been discussed and a stellar-mass black hole is likely favored.

We examine the gamma-ray quasi-periodic variability of PKS 2155-304 with the latest publicly available \emph{Fermi}-LAT Pass 8 data which covers the years from 2008 August to 2016 October. We produce the light curves in two ways, i.e., the exposure-weighted aperture photometry and the maximum likelihood optimization. The light curves are then analyzed by using Lomb-Scargle Periodogram (LSP) and Weighted Wavelet Z-transform (WWZ); and the results reveal a significant quasi-periodicity with a period of $1.74\pm0.13$ years and a significance of $\sim$ 4.9 $\sigma$. The constraint of multifrequencies quasi-periodic variabilities on blazar emission model is discussed.

We report the detection of a GeV $\gamma$-ray source which is spatially overlapping and thus very likely associated with the unidentified very-high-energy (VHE) $\gamma$-ray source HESS J1427-608 with the Pass 8 data recorded by the Fermi Large Area Telescope. The photon spectrum of this source is best described by a power-law with an index of $1.85\pm0.17$ in the energy range of $3-500$ GeV, and the measured flux connects smoothly with that of HESS J1427-608 at a few hundred GeV. This source shows no significant extension and time variation. The broadband GeV-TeV emission over four decades of energies can be well fitted by a single power-law function with an index of 2.0, without obvious indication of spectral cutoff toward high energies. Such a result implies that HESS J1427-608 may be a PeV particle accelerator. We discuss possible nature of HESS J1427-608 according to the multi-wavelength spectral fittings. Given the relatively large errors, either a leptonic or a hadronic model can explain the multi-wavelength data from radio to VHE $\gamma$-rays. The inferred magnetic field strength is a few $\mu$G, which is smaller than typical values of supernova remnants (SNRs), and is consistent with some pulsar wind nebulae (PWNe). On the other hand, the flat $\gamma$-ray spectrum is slightly different from typical PWNe but similar to that of some known SNRs.

III Zw 2 is the prototype of radio-intermediate quasars. Although there is the evidence of possessing strong jet, significant $\gamma$-ray emission has not been reported before. In this work, we carry out a detailed analysis of the latest {\it Fermi}-LAT {\it Pass} 8 data. No significant $\gamma$-ray signal has been detected in the time-averaged 7-year {\it Fermi}-LAT data of III Zw 2, we however have identified two distinct $\gamma$-ray flares with isotropic luminosities of $\sim 10^{45}$ erg $\rm s^{-1}$. Multiwavelength data analysis (including also the optical photometric observations from Yunnan Observatories) are presented and the main finding is the simultaneous optical and $\gamma$-ray flares of III Zw 2 appearing in Nov. 2009. Violent $\gamma$-ray variability with doubling timescale of 2.5 hours was detected in another $\gamma$-ray flare in May 2010, for which the 3-hour $\gamma$-ray peak flux is $\sim 250$ times of the average flux in 7 years. Rather similar behaviors are observed in blazars and the blazar model can reasonably reproduce the spectral energy distribution of III Zw 2 in a wide energy range, strongly suggesting that its central engine resembles that of blazars. In view of its core which shares radio similarities with young radio sources together with the weak extended radio lobe emission, we suggest that III Zw 2 harbors a recurrent activity core and serves as a valuable target for investigating the fuelling and triggering of the activity in radio loud AGNs.

Very recently the Dark Energy Survey (DES) Collaboration has released their second group of Dwarf spheroidal (dSph) galaxy candidates. With the publicly-available Pass 8 data of Fermi-LAT we search for $\gamma-$ray emissions from the directions of these eight newly discovered dSph galaxy candidates. No statistically significant $\gamma-$ray signal has been found in the combined analysis of these sources. With the empirically estimated J-factors of these sources, the constraint on the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ is comparable to the bound set by the joint analysis of fifteen previously known dSphs with kinematically constrained J-factors for the dark matter mass $m_\chi>250$ GeV. In the direction of Tucana III (DES J2356-5935), one of the nearest dSph galaxy candidates that is $\sim 25$ kpc away, there is a weak $\gamma-$ray signal and its peak test statistic (TS) value for the dark matter annihilation channel $\chi\chi\rightarrow \tau^{+}\tau^{-1}$ is $\approx 6.7$ at $m_\chi \sim 15$ GeV. The significance of the possible signal likely increases with time. More data is highly needed to pin down the physical origin of such a GeV excess.

We systematically analyze the 6-year {\it Fermi}/LAT data of the lobe-dominated quasars (LDQs) in the complete LDQ sample from 3CRR survey and report the discovery of high-energy $\gamma$-ray emission from 3C 275.1. The $\gamma$-ray emission of 3C 207 is confirmed and significant variability of the lightcurve is identified. We do not find statistically significant $\gamma$-ray emission from other LDQs. 3C 275.1 is the known $\gamma$-ray quasar with the lowest core dominance parameter (i.e., $R=0.11$). We also show that both the northern radio hotspot and parsec jet models can reasonably reproduce the $\gamma$-ray data. The parsec jet model, however, is favored by the potential $\gamma$-ray variability at the timescale of months. We suggest that some dimmer $\gamma$-ray LDQs will be detected in the future and LDQs could contribute non-negligibly to the extragalactic $\gamma$-ray background.

We report a detailed analysis of $\gamma$-ray data of 2FGL J0509.9+1802 observed by the Large Area Telescope on board {\it Fermi} satellite, especially focusing on April 2013 when extraordinary $\gamma$-ray variability has been detected. Localization of $\gamma$-ray emission during this epoch suggests that the $\gamma$-ray source only associates with PKS 0507+17. The $\gamma$-ray emission of PKS 0507+17 is identified at the first time. The daily peak flux is over two orders of magnitude higher than the first two-year average flux, giving an isotropic $\gamma$-ray luminosity of $\simeq4\times10^{48}$ erg $\rm s^{-1}$. Rapid $\gamma$-ray variability with doubling time of 2-3 hours has been detected. Such a short doubling timescale has been detected for only a few bright blazars and indicates a location of $\gamma$-ray emission inside the broad line region. Together with the bluer-when-brighter $\gamma$-ray spectra, the variability phenomena could be well explained by the classic flat-spectrum radio quasar variability model that includes a fast injection of accelerated electrons and the external Compton cooling process.

BL Lacertae (BL Lac) objects at high redshifts ($z\geq 2$) are rarely detected. Through careful analysis of the SDSS spectrum, BZB J1450+5201 is confirmed to be a high-$z$ BL Lac object with $z\geq$ 2.471 by identifying the Ly$\alpha$ 1216 and CIV 1548/1550 absorption lines. This indicates that BZB J1450+5201 is the most distant BL Lac object discovered to date. Careful analysis of the five-year \fermi data of 2FGL J1451.0+5159 shows that its $\gamma$-ray emission is robust with confidence level of 6.2$\sigma$ at 1-3 GeV and 6.7$\sigma$ at 3-10 GeV, and that the confusion of bright neighbor is negligible, which can not be fixed in the analysis of the two-year data. Meanwhile, 2FGL J1451.0+5159 is confirmed to be associated with BZB J1450+5201 using the five-year data. The analysis of multiwavelength data, from radio to $\gamma$-ray energies, indicates BZB J1450+5201 is an intermediate synchrotron peaked (ISP) source and consistent with distributions of other ISP sources at lower redshifts in the second LAT AGN catalog. The pure SSC model seems to be disfavoured, while scattering of weak external emission plus SSC process can provide a satisfactory description of the broadband emission.

S5 0716+714 is a typical BL Lacertae object. In this paper we present the analysis and results of long term simultaneous observations in the radio, near-infrared, optical, X-ray and $\gamma$-ray bands, together with our own photometric observations for this source. The light curves show that the variability amplitudes in $\gamma$-ray and optical bands are larger than those in the hard X-ray and radio bands and that the spectral energy distribution (SED) peaks move to shorter wavelengths when the source becomes brighter, which are similar to other blazars, i.e., more variable at wavelengths shorter than the SED peak frequencies. Analysis shows that the characteristic variability timescales in the 14.5 GHz, the optical, the X-ray, and the $\gamma$-ray bands are comparable to each other. The variations of the hard X-ray and 14.5 GHz emissions are correlated with zero-lag, so are the V band and $\gamma$-ray variations, which are consistent with the leptonic models. Coincidences of $\gamma$-ray and optical flares with a dramatic change of the optical polarization are detected. Hadronic models do not have the same nature explanation for these observations as the leptonic models. A strong optical flare correlating a $\gamma$-ray flare whose peak flux is lower than the average flux is detected. Leptonic model can explain this variability phenomenon through simultaneous SED modeling. Different leptonic models are distinguished by average SED modeling. The synchrotron plus synchrotron self-Compton (SSC) model is ruled out due to the extreme input parameters. Scattering of external seed photons, such as the hot dust or broad line region emission, and the SSC process are probably both needed to explain the $\gamma$-ray emission of S5 0716+714.