The middle-aged supernova remnant (SNR) CTB 37A is known to interact with several dense molecular clouds through the detection of shocked ${\rm H_{2}}$ and OH 1720 MHz maser emission. In the present work, we use eight years of $\textit Fermi$-LAT Pass 8 data, with an improved point-spread function and an increased acceptance, to perform detailed morphological and spectral studies of the $\gamma$-ray emission toward CTB 37A from 200 MeV to 200 GeV. The best fit of the source extension is obtained for a very compact Gaussian model with a significance of 5.75$\sigma$ and a 68\% containment radius of $0.116^{\circ}$ $\pm$ $0.014^{\circ}_{\rm stat}$ $\pm$ $0.017^{\circ}_{\rm sys}$ above 1 GeV, which is larger than the TeV emission size. The energy spectrum is modeled as a LogParabola, resulting in a spectral index $\alpha$ = 1.92 $\pm$ 0.19 at 1 GeV and a curvature $\beta$ = 0.18 $\pm$ 0.05, which becomes softer than the TeV spectrum above 10 GeV. The SNR properties, including a dynamical age of 6000 yr, are derived assuming the Sedov phase. From the multiwavelength modeling of emission toward the remnant, we conclude that the nonthermal radio and GeV emission is mostly due to the reacceleration of preexisting cosmic rays (CRs) by radiative shocks in the adjacent clouds. Furthermore, the observational data allow us to constrain the total kinetic energy transferred to the trapped CRs in the clouds. Based on these facts, we infer a composite nature for CTB 37A to explain the broadband spectrum and to elucidate the nature of the observed $\gamma$-ray emission.

An accurate estimate of the interstellar gas density distribution is crucial to understanding the interstellar medium (ISM) and Galactic cosmic rays (CRs). To comprehend the ISM and CRs in a local environment, a study of the diffuse $\gamma$-ray emission in a mid-latitude region of the third quadrant was performed. The $\gamma$-ray data in the 0.1--25.6~GeV energy range of the Fermi Large Area Telescope (LAT) and other interstellar gas tracers such as the HI4PI survey data and the Planck dust thermal emission model were used, and the northern and southern regions were analyzed separately. The variation of the dust emission Dem with the total neutral gas column density NH was studied in high dust-temperature areas, and the NH/Dem ratio was calibrated using $\gamma$-ray data under the assumption of a uniform CR intensity in the studied regions. The measured integrated $\gamma$-ray emissivities above 100~MeV are $(1.58\pm0.04)\times10^{-26}~\mathrm{photons~s^{-1}~sr^{-1}~H\mbox{-}atom^{-1}}$ and $(1.59\pm0.02)\times10^{-26}~\mathrm{photons~s^{-1}~sr^{-1}~H\mbox{-}atom^{-1}}$ in the northern and southern regions, respectively, supporting the existence of a uniform CR intensity in the vicinity of the solar system. While most of the gas can be interpreted to be HI with a spin temperature of $T_\mathrm{S} = 125~\mathrm{K}$ or higher, an area dominated by optically thick HI with $T_\mathrm{S} \sim 40~\mathrm{K}$ was identified.

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 present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis on the first 7.4 years of \textit{Fermi} observations, and in two separate energy bands 0.1$-$0.8 GeV and 0.8$-$300 GeV, a total of 4547 flares has been detected with a significance greater than $6\sigma$ (before trials), on the time scale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources are identified. Likely counterparts, based on positional coincidence, have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of the freshly accelerated electrons is never harder than $p\sim$2.

We present a measurement of the cosmic-ray electron+positron spectrum between 7 GeV and 2 TeV performed with almost seven years of data collected with the Fermi Large Area Telescope. We find that the spectrum is well fit by a broken power law with a break energy at about 50 GeV. Above 50 GeV, the spectrum is well described by a single power law with a spectral index of $3.07 \pm 0.02 \; (\text{stat+syst}) \pm 0.04 \; (\text{energy measurement})$. An exponential cutoff lower than 1.8 TeV is excluded at 95\% CL.

A study of the interstellar medium (ISM) and cosmic rays (CRs) using Fermi Large Area Telescope (LAT) data, in a region encompassing the nearby molecular clouds MBM 53, 54, and 55 and a far-infrared loop-like structure in Pegasus, is reported. By comparing Planck dust thermal emission model with Fermi-LAT gamma-ray data, it was found that neither the dust radiance (R) nor the dust opacity at 353 GHz (tau353) were proportional to the total gas column density N(Htot) primarily because N(Htot)/R and N(Htot)/tau353 depend on the dust temperature (Td). The N(Htot) distribution was evaluated using gamma-ray data by assuming the regions of high Td} to be dominated by optically thin atomic hydrogen (HI) and by employing an empirical linear relation of N(Htot)/R to Td. It was determined that the mass of the gas not traced by the 21-cm or 2.6-mm surveys is ~25% of the mass of HI in the optically thin case and is larger than the mass of the molecular gas traced by carbon monoxide by a factor of up to 5. The measured gamma-ray emissivity spectrum is consistent with a model based on CR spectra measured at the Earth and the nuclear enhancement factor of <=1.5. It is, however, lower than local HI emissivities reported by previous Fermi-LAT studies employing different analysis methods and assumptions on ISM properties by 15%-20% in energies below a few GeV, even if we take account of the statistical and systematic uncertainties. The origin of the discrepancy is also discussed.

The first phase of the Alborz Observatory Array (Alborz-I) consists of 20 plastic scintillation detectors each one with surface area of 0.25 $m^{2}$ spread over an area of 40$\times$40 $m^{2}$ realized to the study of Extensive Air Showers around the $\it knee$ at the Sharif University of Technology campus. The first stage of the project including construction and operation of a prototype system has now been completed and the electronics that will be used in the array instrument has been tested under field conditions. In order to achieve a realistic estimate of the array performance, a large number of simulated CORSIKA~\cite{a} showers have been used. In the present work, theoretical results obtained in the study of different array layouts and trigger conditions are described. Using Monte Carlo simulations of showers the rate of detected events per day and the trigger probability functions, i.e., the probability for an extensive air shower to trigger a ground based array as a function of the shower core distance to the center of array are presented for energies above 1 TeV and zenith angles up to 60$^{\circ}$. Moreover, the angular resolution of the Alborz-I array is obtained.