Star-forming galaxies are huge reservoirs of cosmic rays (CRs) and these CRs convert a significant fraction of their energy into $\gamma$-rays by colliding with the interstellar medium (ISM). Several nearby star-forming galaxies have been detected in GeV-TeV $\gamma$-rays. It is also found that the $\gamma$-ray luminosities in 0.1-100 GeV correlate well with indicators of star formation rates of the galaxies, such as the total infrared (IR) luminosity. In this paper, we report a systematic search for possible $\gamma$-ray emission from galaxies in the IRAS Revised Bright Galaxies Sample, using 11.4 years of $\gamma$-ray data taken by the Fermi Large Area Telescope (LAT). Two new galaxies, M33 and Arp 299, are detected significantly. The two galaxies are consistent with the empirical correlation between the $\gamma$-ray luminosity and total infrared luminosity, suggesting that their $\gamma$-ray emissions should mainly originate from CRs interacting with ISM. Nevertheless, there is a tentative evidence that the flux of the $\gamma$-ray emission from Arp~299 is variable. If the variability is true, part of the emission from Arp 299 should originate from the obscured AGN in this interacting galaxy system. In addition, we find that the $\gamma$-ray excess from M33 is located at the northeast region of the galaxy, where a supergiant H II region, NGC604, resides. This indicates that some bright star-forming regions in spiral galaxies could play a dominant role in the galaxy in producing $\gamma$-ray emission.

In this work, we collect quasi-simultaneous infrared, optical, X-ray and $\gamma$-ray data of 60 $Fermi$-4LAC flat spectrum radio quasars (FSRQs). In the framework of the conventional one-zone leptonic model, we investigate the physical properties of $Fermi$-4LAC FSRQs' jets by modeling their quasi-simultaneous spectral energy distributions (SEDs). Our main results are summarized as follows. (1) There is a linear correlation between synchrotron peak frequency and curvature of the electron energy distribution. As suggested by previous works, the slope of the best linear fitting equation of this correlation is consistent with statistic acceleration which needs a fluctuation of fractional acceleration gain. (2) The gamma-ray dissipation regions are located at the range from 0.1 to 10 pc away from the super-massive black hole, and located outside the broad-line region (BLR) and within the dusty torus (DT). (3) A size relation $P_{\rm e}$ (the kinetic power carried in relativistic electrons) $\sim$ $P_{\rm B}$ (Poynting flux) $\leq$ $P_{\rm r}$ (the radiative power ) $<$ $P_{\rm p}$ (the kinetic power in cold protons) is found in our modeling. Among them, $P_{\rm e}\sim P_{\rm B}$ suggests that SEDs of almost all FSRQs with parameters are close to equipartition between the magnetic field and the relativistic electrons. The $P_{\rm e} < P_{\rm r}$ suggest that the most energy of the relativistic electrons are dissipated by EC radiation for FSRQs. (4) There is an anti-correlation between the peak energy of SEDs ($\gamma_{\rm peak}$) and the jet power ($P_{\rm jet}$), which is consistent with the blazar sequence.

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

The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time.

Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we search for possible GeV counterparts of three TeV halo candidates with Fermi Large Area Telescopes. We detect a new spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099, which is also detected by HESS (HESS J1912+101). We find that the size of the GeV source is significantly larger than that of the TeV emission measured by HESS, and a spatial template characteristic of a PWN can fit the GeV data. We suggest that the GeV source is an extended PWN powered by the central middle-aged pulsar PSR J1913+1011. This discovery favors the TeV halo scenario for the TeV source 2HWC J1912+099 (HESS J1912+101), although the possible shell-like morphology measured by HESS challenges this interpretation. Alternatively, the TeV emission could be dominated by a supernova remnant (SNR) via the hadronic process. Future multi-wavelength observations of the source and more precise measurements of the spatial profile of the TeV emission will be useful to distinguish between the two scenarios.

Recently, very high-energy photons above 100 GeV were reported to be detected from GRB~190114C and GRB~180720B at, respectively, 100-1000 s and 10 hours after the burst. We model the available broad-band data of both GRBs with the synchrotron plus synchrotron self-Compton (SSC) emission of the afterglow shocks. We find that the sub-TeV emission of GRB~180720B can be interpreted as the SSC emission from afterglow shocks expanding in a constant density circum-burst medium. The SSC emission of GRB~190114C dominates over the synchrotron component from GeV energies at $\sim100$ s, which can explain the possible hard spectrum of the GeV emission at this time. The extrapolated flux of this SSC component to sub-TeV energies can explain the high-significance detection of GRB~190114C by the MAGIC telescope. The parameter values (such as the circum-burst density and shock microphysical parameters) in the modeling are not unusual for both GRBs, implying that the detection of sub-TeV photons from these two bursts should be attributed to their large burst energies and low redshifts.

We analyze 10 years of Fermi-LAT data towards the SS433/W50 region. With the latest source catalog and diffuse background models, the gamma-ray excess from SS433/W50 is detected with a significance of 6{\sigma} in the photon energy range of 500 MeV - 10 GeV. Our analysis indicates that an extended flat disk morphology is preferred over a point-source description, suggesting that the GeV emission region is much larger than that of the TeV emission detected by HAWC. The size of the GeV emission is instead consistent with the extent of the radio nebula W50, a supernova remnant being distorted by the jets, so we suggest that the GeV emission may originate from this supernova remnant. The spectral result of the GeV emission is also consistent with an supernova remnant origin. We also derive the GeV flux upper limits on the TeV emission region, which put moderate constrains on the leptonic models to explain the multiwavelength data.

It has long been suggested that nearby pulsars within $\sim 1 \,{\rm kpc}$ are the leading candidate of the 10-500 GeV cosmic-ray positron excess measured by PAMELA and other experiments. The recent measurement of surface brightness profile of TeV nebulae surrounding Geminga and PSR~B0656+14 by the High-Altitude Water Cherenkov Observatory (HAWC) suggests inefficient diffusion of particles from the sources, giving rise to a debate on the pulsar interpretation of the cosmic-ray positron excess. Here we argue that GeV observations provide more direct constraints on the positron density in the TeV nebulae in the energy range of 10-500 GeV and hence on the origin of the observed positron excess. Motivated by this, we search for GeV emission from the TeV nebulae with the \textsl{Fermi} Large Area Telescope (LAT). No spatially-extended GeV emission is detected from these two TeV nebulae in the framework of two-zone diffusion spatial templates, suggesting a relatively low density of GeV electrons/positrons in the TeV nebulae. A joint modelling of the data from HAWC and \textsl{Fermi}-LAT disfavors Geminga and PSR~B0656+14 as the dominant source of the positron excess at $\sim 50-500$ GeV for the usual Kolmogorov-type diffusion, while for an energy-independent diffusion, a dominant part of the positron excess contributed by them cannot be ruled out by the current data.

Star-forming regions on different scales, such as giant molecular clouds in our Galaxy and star-forming galaxies, emit GeV gamma-rays. These are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. It has recently been shown that the gamma-ray luminosity ($L_\gamma$) of star-forming galaxies is well correlated with their star formation rates (SFR). We investigated \textsl{Fermi} data of eight Galactic molecular clouds in the Gould belt and found that molecular clouds do not follow the $L_\gamma-{\rm SFR}$ correlation of star-forming galaxies. We also compared the scaling relations of gamma-ray luminosity, SFR, and the gas mass for molecular clouds and star-forming galaxies. Using a multiple-variable regression analysis, we found different dependences of gamma-ray emission on SFR or mass for molecular clouds and star-forming galaxies. This suggests that different mechanisms may govern the production of gamma-rays in these two types of sources. Specifically, the strong dependence on mass supports that gamma-ray emission of molecular clouds primarily comes from {\em \textup{passive}} interaction by diffuse Galactic CRs, whereas the strong dependence on SFR supports that gamma-ray emission of star-forming galaxies originates from CRs that are accelerated by local {\em \textup{active}} sources.

Many galaxy clusters have giant halos of non-thermal radio emission, indicating the presence of relativistic electrons in the clusters. Relativistic protons may also be accelerated by merger and/or accretion shocks in galaxy clusters. These cosmic-ray (CR) electrons and/or protons are expected to produce gamma-rays through inverse-Compton scatterings or inelastic $pp$ collisions respectively. Despite of intense efforts in searching for high-energy gamma-ray emission from galaxy clusters, conclusive evidence is still missing so far. Here we report the discovery of $\ge 200$ MeV gamma-ray emission from the Coma cluster direction with an unbinned likelihood analysis of the 9 years of {\it Fermi}-LAT Pass 8 data. The gamma-ray emission shows a spatial morphology roughly coincident with the giant radio halo, with an apparent excess at the southwest of the cluster. Using the test statistic analysis, we further find tentative evidence that the gamma-ray emission at the Coma center is spatially extended. The extended component has an integral energy flux of $\sim 2\times 10^{-12}{\rm \ erg\ cm^{-2}\ s^{-1}}$ in the energy range of 0.2 - 300 GeV and the spectrum is soft with a photon index of $\simeq-2.7$. Interpreting the gamma-ray emission as arising from CR proton interaction, we find that the volume-averaged value of the CR to thermal pressure ratio in the Coma cluster is about $\sim 2\%$. Our results show that galaxy clusters are likely a new type of GeV gamma-ray sources, and they are probably also giant reservoirs of CR protons.

The Einstein's weak equivalence principle (WEP) is one of the foundational assumptions of general relativity and some other gravity theories. In the theory of parametrized post-Newtonian (PPN), the difference between the PPN parameters $\gamma$ of different particles or the same type of particle with different energies, $\Delta \gamma$, represents the violation of WEP. Current constraints on $\Delta \gamma$ are derived from the observed time delay between correlated particles of astronomical sources. However, the observed time delay is contaminated by other effects, such as the time delays due to different particle emission times, the potential Lorentz invariance violation and none-zero photon rest mass. Therefore, current constraints are only upper limits. Here we propose a new method to test WEP based on the fact that the gravitational time delay is direction-dependent while others are not. This is the first method which can naturally correct other time delay effects. Using the time delay measurements of BASTE gamma-ray burst (GRB) sample and the gravitational potential of local super galaxy cluster Laniakea, we find that the constraint on $\Delta \gamma$ of different energy photons can be as low as $10^{-14}$. In the future, if more gravitational wave events and fast radio bursts with much more precise time delay measurements are observed, this method can give a reliable and tight constraint on WEP.

Recently, gamma-ray emission in the direction of Coma, with a TS value of $\sim 40$, has been reported. In this work we will discuss the possibility of such a residual emission coming from dark matter annihilation. Our results show that the gamma-ray emission within the Coma region is also spatially correlated to the mass distribution derived from weak gravitational lensing measurements very well. However the dark matter models are not supported by the spectral analysis results and constraints by observations of other targets. Thus we derive the upper limits of the dark matter annihilation cross section according to the observation of the Coma region.

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 present time-resolved spectral analysis of the steep decay segments of 29 bright X-ray flares of gamma-ray bursts (GRBs) observed with the Swift/X-ray telescope, and model their lightcurves and spectral index evolution behaviors with the curvature effect model. Our results show that the observed rapid flux decay and strong spectral index evolution with time can be well fit with this model, and the derived characteristic timescales ($t_c$) are in the range of $33\sim 264$ seconds. Using an empirical relation between the peak luminosity and the Lorentz factor derived from the prompt gamma-rays, we estimate the Lorentz factors of the flares ($\Gamma_{\rm X}$). We obtain $\Gamma_{\rm X}=17\sim 87$ with a median value of $52$, which is smaller than the initial Lorentz factors of prompt gamma-ray fireballs. With the derived $t_c$ and $\Gamma_{\rm X}$, we constrain the radiating regions of 13 X-ray flares, yielding $R_{\rm X}=(0.2\sim 1.1)\times 10^{16}$ cm, which are smaller than the radii of the afterglow fireballs at the peak times of the flares. A long evolution feature from prompt gamma-ray phase to the X-ray epoch is found by incorporating our results with a sample of GRBs whose initial Lorentz factors are available in literatures, i.e., $\Gamma\propto [t_{p}/(1+z)]^{-0.69\pm 0.06}$. These results may shed lights on the long term evolution of GRB central engines.

X-ray flares are generally supposed to be produced by the later central engine activities, and may share the similar physical origin with prompt emission of gamma-ray bursts (GRBs). In this paper, we have analyzed all significant X-ray flares from the GRBs observed by {\em Swift} from April 2005 to March 2015. The catalog contains 468 bright X-ray flares, including 200 flares with redshifts. We obtain the fitting results of X-ray flares, such as start time, peak time, duration, peak flux, fluence, peak luminosity, and mean luminosity. The peak luminosity decreases with peak time, following a power-law behavior $L_p \propto T_{peak,z}^{-1.27}$. The flare duration increases with peak time. The 0.3-10 keV isotropic energy of X-ray flares distribution is a lognormal peaked at $10^{51.2}$ erg. We also study the frequency distributions of flare parameters, including energies, durations, peak fluxes, rise times, decay times and waiting times. Power-law distributions of energies, durations, peak fluxes, and waiting times are found in GRB X-ray flares and solar flares. These distributions could be well explained by a fractal-diffusive, self-organized criticality model. Some theoretical models basing on magnetic reconnection have been proposed to explain X-ray flares. Our result shows that the relativistic jets of GRBs may be Poynting-flux dominated.

X-ray flares have been discovered in black hole systems, such as gamma-ray bursts, the tidal disruption event Swift J1644+57, the supermassive black hole Sagittarius A$^*$ at the center of our Galaxy, and some active galactic nuclei. Their occurrences are always companied by relativistic jets. However, it is still unknown whether there is a physical analogy among such X-ray flares produced in black hole systems spanning nine orders of magnitude in mass. Here we report the observed data of X-ray flares, and show that they have three statistical properties similar to solar flares, including power-law distributions of energies, durations, and waiting times, which both can be explained by a fractal-diffusive self-organized criticality model. These statistical similarities, together with the fact that solar flares are triggered by a magnetic reconnection process, suggest that all of the X-ray flares are consistent with magnetic reconnection events, implying that their concomitant relativistic jets may be magnetically dominated.

X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitted with the absorbed single power-law or the Band function models. More interestingly, the joint spectra of the other 13 X-ray flares are fitted with the absorbed single power-law model plus a black body (BB) component. Phenomenally, the observed spectra of these 13 flares are analogous to several GRBs with a thermal component, but only with a much lower temperature of $kT=1\sim 3$ keV. Assuming that the thermal emission is the photosphere emission of the GRB fireball, we derive the fireball properties of the 13 flares that have redshift measurements, such as the bulk Lorentz factor $\Gamma_{\rm ph}$ of the outflow. The derived $\Gamma_{\rm ph}$ range from $50$ to $150$ and a relation of $\Gamma_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $\Gamma_0-L_{\rm iso}$ relation that are derived for the prompt gamma-ray emission. We discuss the physical implications of these results within the content of jet composition and radiation mechanism of GRBs and X-ray flares.

We jointly analyze the gamma-ray burst (GRB) data observed with BAT and XRT on board the Swift mission to present a global view on the internal energy dissipation processes in GRBs, including precursors, prompt gamma-ray emission, extended soft gamma-ray emission, and late X-ray flares. The Bayesian block method is utilized to analyze the BAT lightcurves to identify various emission episodes. Our results suggest that these emission components likely share a same physical origin, which is repeated activation of the GRB central engine. What we observe in the gamma-ray band may be the tip-of-iceberg of more extended underlying activities. The precursor emission, which is detected in about 10% of {\em Swift} GRBs, is preferably detected in those GRBs that have a massive star core-collapse origin. The soft extended emission (EE) tail, on the other hand, is preferably detected in those GRBs that have a compact star merger origin. Bright X-ray emission is detected during the BAT quiescent phases prior to subsequent gamma-ray peaks, implying that X-ray emission may be detectable prior the BAT trigger time. Future GRB alert instruments with soft X-ray capability would be essential to reveal the early stage of GRB central engine activities, sheding light into jet composition and jet launching mechanism in GRBs.

GRB 121027A is un-usual with its extremely long-lasting, energetic X-ray flares. The total energy release in X-ray flares is about one order of magnitude higher than prompt gamma-rays, making it special from most long GRBs. We show that while the prompt gamma-ray emission satisfies the empirical E_{iso}-E_{p} relation of typical long GRBs, the X-ray flares, whose spectra can be fit with a cutoff-power-law model with well-constrained E_p, significantly deviate from such a relation. Nonetheless, a time-resolved spectral analysis of X-ray flares suggest that the X-ray emission is consistent with the L_{iso}-E_{p} relation of long GRBs. We constrain the minimum Lorentz factor of the X-ray flares to be ~14, which is consistent with the Gamma-L_{iso} relation. Our results imply that prompt gamma-ray emission and late X-ray flares share the similar radiation mechanism, but originate from the outflows with different Lorentz factors. We search for similar GRBs from the Swift GRB archives, and find that the z=6.29 GRB 050904 is a carbon copy of GRB 121027A. The long-lasting, energetic X-ray flares in these GRBs demand significant accretion at late times, which point towards a large-radius progenitor star.

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{\gamma, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{{\gamma}, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $\gamma$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors.