While the Band function or other phenomenological functions are commonly used to fit GRB prompt emission spectra, we propose a new parametric function that is inspired by an internal shock physical model. We use this function as a proxy of the model to confront it easily to GRB observations. We built a parametric function that represents the spectral form of the synthetic bursts provided by our internal shock synchrotron model (ISSM). We simulated the response of the Fermi instruments to the synthetic bursts and fitted the obtained count spectra to validate the ISSM function. Then, we applied this function to a sample of 74 bright GRBs detected by the Fermi/GBM, and we computed the width of their spectral energy distributions around their peak energy. For comparison, we fitted also the phenomenological functions that are commonly used in the literature. Finally, we performed a time-resolved analysis of the broadband spectrum of GRB 090926A, which was jointly detected by the Fermi GBM and LAT. The ISSM function reproduces 81% of the spectra in the GBM bright GRB sample, versus 59% for the Band function, for the same number of parameters. It gives also relatively good fits to the GRB 090926A spectra. The width of the MeV spectral component that is obtained from the fits of the ISSM function is slightly larger than the width from the Band fits, but it is smaller when observed over a wider energy range. Moreover, all of the 74 analysed spectra are found to be significantly wider than the synthetic synchrotron spectra. We discuss possible solutions to reconcile the observations with the internal shock synchrotron model, such as an improved modeling of the shock micro-physics or more accurate spectral measurements at MeV energies.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

GRB 131108A is a bright long Gamma-Ray Burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the \textit{Fermi Gamma-ray Space Telescope}. Dedicated temporal and spectral analyses reveal three $\gamma$-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the GBM band (10 keV--10 MeV). The high-energy light curve of GRB 131108A (100 MeV -- 10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the $\gamma$-ray flares finds that the three flares are 6 -- 20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 $\sim$ 2.0) $\times$ 10$^{-6}$ erg cm$^{-2}$, is comparable to that of emission during the steep decay phase, 1.7 $\times$ 10$^{-6}$ erg cm$^{-2}$. The total fluence from three $\gamma$-ray flares is 5.3 $\times$ 10$^{-6}$ erg cm$^{-2}$. The three $\gamma$-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in $\sim$ 50\% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect ($\hat{\alpha}$ = 2 + $\hat{\beta}$), which is the first observational evidence of the high-latitude emission in the GeV energy band.

GRB 160709A is one of the few bright short gamma-ray bursts detected by both the Gamma-ray Burst Monitor and the Large Area Telescope on board the $Fermi$ $Gamma$-$ray$ $Space$ $Telescope$. The $\gamma$-ray prompt emission of GRB 160709A is adequately fitted by combinations of three distinct components: (i) a nonthermal component described by a power law (PL) with a high-energy exponential cutoff, (ii) a thermal component modeled with a Planck function, and (iii) a second nonthermal component shaped by an additional PL crossing the whole $\gamma$-ray spectrum. While the thermal component dominates during $\sim$ 0.12 s of the main emission episode of GRB 160709A with an unusually high temperature of $\sim$ 340 keV, the nonthermal components dominate in the early and late time. The thermal component is consistent with the photospheric emission resulting in the following parameters: the size of the central engine, $R_{0}$ = $3.8 \substack{+5.9 \\ -1.8}$ $\times 10^{8}$ cm, the size of the photosphere, R$_{ph}$ = $7.4 \substack{+0.8 \\ -1.2}$ $\times 10^{10}$ cm, and a bulk Lorentz factor, $\Gamma$ = $728 \substack{+75 \\ -93}$ assuming a redshift of 1. The slope of the additional PL spectrum stays unchanged throughout the burst duration; however, its flux decreases continuously as a function of time. A standard external shock model has been tested for the additional PL component using the relation between the temporal and spectral indices (the closure relation). Each set of spectral and temporal indices from two energy bands (200 keV--40 MeV and 100 MeV--10 GeV) satisfies a distinct closure relation. From the closure relation test we derived the index for the electron spectral distribution, $p$ = 2.5 $\pm$ 0.1. The interaction of the jet with the interstellar environment is preferred over the interaction with the wind medium.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$\sigma$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15.

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

We report on a new time-resolved analysis of the GRB 090926A broad-band spectrum during its prompt phase, and on its interpretation in the framework of prompt emission models. We characterize the high-energy extra component of GRB 090926A using Pass 8 data from the Fermi Large Area Telescope (LAT), which offer a greater sensitivity than any data set used in previous studies of this burst, particularly in the 30-100 MeV energy band. Then, we combine the LAT data with the Fermi Gamma-ray Burst Monitor (GBM) in joint spectral fits in order to characterize the time evolution of the broad-band spectrum from keV to GeV energies. Careful attention is paid to the LAT systematic effects. Finally, we perform a temporal analysis of the light curves and we compute the variability time scales during and after the bright spike at ~10 s post-trigger. Our analysis confirms and better constrains the spectral break in the high-energy extra component, which has been previously reported during the bright spike. Furthermore, it reveals that the spectral attenuation persists at later times, with an increase of the break characteristic energy up to the GeV domain until the end of the prompt phase. We interpret the high-energy spectral break as caused by photon opacity to pair creation. Requiring that all emissions are produced above the photosphere of GRB 090926A, we compute the bulk Lorentz factor of the outflow, $\Gamma$. The latter decreases from 230 during the spike to 100 at the end of the prompt emission. Assuming, instead, that the spectral break reflects the natural curvature of the inverse Compton spectrum, lower limits corresponding to larger values of $\Gamma$ are also derived. Combined with the extreme temporal variability of GRB 090926A, these Lorentz factors lead to emission radii $R\sim10^{14}$ cm which are consistent with an internal origin of both the keV-MeV and GeV prompt emissions.

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($\Delta f / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.