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.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

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.

The first simultaneous detection of a short gamma-ray burst (SGRB) with a gravitational-wave (GW) signal ushered in a new era of multi-messenger astronomy. In order to increase the number of SGRB-GW simultaneous detections, we need full sky coverage in the gamma-ray regime. BurstCube, a CubeSat for Gravitational Wave Counterparts, aims to expand sky coverage in order to detect and localize gamma-ray bursts (GRBs). BurstCube will be comprised of 4 Cesium Iodide scintillators coupled to arrays of Silicon photo-multipliers on a 6U CubeSat bus (a single U corresponds to cubic unit $\sim$10 cm $\times$ 10 cm $\times$ 10 cm) and will be sensitive to gamma-rays between 50 keV and 1 MeV, the ideal energy range for GRB prompt emission. BurstCube will assist current observatories, such as $Swift$ and $Fermi$, in the detection of GRBs as well as provide astronomical context to gravitational wave events detected by Advanced LIGO, Advanced Virgo, and KAGRA. BurstCube is currently in its development and testing phase to prepare for launch readiness in the fall of 2021. We present the mission concept, preliminary performance, and status.

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.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Indirect searches for dark matter through Standard Model products of its annihilation generally assume a cross-section which is dominated by a term independent of velocity ($s$-wave annihilation). However, in many DM models an $s$-wave annihilation cross-section is absent or helicity suppressed. To reproduce the correct DM relic density in these models, the leading term in the cross section is proportional to the DM velocity squared ($p$-wave annihilation). Indirect detection of such $p$-wave DM is difficult because the average velocities of DM in galaxies today are orders of magnitude slower than the DM velocity at the time of decoupling from the primordial thermal plasma, suppressing the annihilation cross-section today by some five orders of magnitude relative to its value at freeze out. Thus $p$-wave DM is out of reach of traditional searches for DM annihilations in the Galactic halo. Near the region of influence of a central supermassive black hole, such as Sgr A$^*$, however, DM can form a localized over-density known as a `spike'. In such spikes the DM is predicted to be both concentrated in space and accelerated to higher velocities, allowing the $\gamma$-ray signature from its annihilation to potentially be detectable above the background. We use the $Fermi$ Large Area Telescope to search for the $\gamma$-ray signature of $p$-wave annihilating DM from a spike around Sgr A$^*$ in the energy range 10 GeV-600 GeV. Such a signal would appear as a point source and would have a sharp line or box-like spectral features difficult to mimic with standard astrophysical processes, indicating a DM origin. We find no significant excess of $\gamma$ rays in this range, and we place upper limits on the flux in $\gamma$-ray boxes originating from the Galactic Center. This result, the first of its kind, is interpreted in the context of different models of the DM density near Sgr A$^*$.

The era of precision cosmology has revealed that about 85% of the matter in the universe is dark matter. Two well-motivated candidates are weakly interacting massive particles (WIMPs) and weakly interacting sub-eV particles (WISPs) (e.g. axions). Both WIMPs and WISPs possess distinct {\gamma}-ray signatures. Over the last decade, data taken between 50 MeV to >300 GeV by the Fermi Large Area Telescope (Fermi-LAT) have provided stringent constraints on both classes of dark matter models. Thus far, there are no conclusive detections. However, there is an intriguing {\gamma}-ray excess associated with the Galactic center that could be explained by WIMP annihilation. At lower energies, the poor angular resolution of the Fermi-LAT makes source identification challenging, inhibiting our ability to more sensitively probe both the Galactic center excess, as well as lower-mass WIMP and WISP models. Additionally, targeted WISP searches (e.g., those probing supernovae and blazars) would greatly benefit from enhanced energy resolution and polarization measurements in the MeV range. To address these issues, a new telescope that is optimized for MeV observations is needed. Such an instrument would allow us to explore new areas of dark matter parameter space and provide unprecedented access to its particle nature.

Latest precise cosmic-ray (CR) measurements and present gamma-ray observations have started challenging our understanding of CR transport and interaction in the Galaxy. Moreover, because the density of CRs is similar to the density of the magnetic field, gas, and starlight in the interstellar medium (ISM), CRs are expected to affect the ISM dynamics, including the physical and chemical processes that determine transport and star formation. In this context, observations of gamma-ray emission at MeV energies produced by the low-energy CRs are very important and urgent. A telescope covering the energy range between ~0.1 MeV and a few GeV with a sensitivity more than an order of magnitude better than previous instruments would allow for the first time to study in detail the low-energy CRs, providing information on their sources, their spectra throughout the Galaxy, their abundances, transport properties, and their role on the evolution of the Galaxy and star formation. Here we discuss the scientific prospects for studies of CRs, ISM (gas, interstellar photons, and magnetic fields) and associated gamma-ray emissions with such an instrument.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

The 511 keV line from positron annihilation in the Galaxy was the first $\gamma$-ray line detected to originate from outside our solar system. Going into the fifth decade since the discovery, the source of positrons is still unconfirmed and remains one of the enduring mysteries in $\gamma$-ray astronomy. With a large flux of $\sim$10$^{-3}$ $\gamma$/cm$^{2}$/s, after 15 years in operation INTEGRAL/SPI has detected the 511 keV line at $>50\sigma$ and has performed high-resolution spectral studies which conclude that Galactic positrons predominantly annihilate at low energies in warm phases of the interstellar medium. The results from imaging are less certain, but show a spatial distribution with a strong concentration in the center of the Galaxy. The observed emission from the Galactic disk has low surface brightness and the scale height is poorly constrained, therefore, the shear number of annihilating positrons in our Galaxy is still not well know. Positrons produced in $\beta^+$-decay of nucleosynthesis products, such as $^{26}$Al, can account for some of the annihilation emission in the disk, but the observed spatial distribution, in particular the excess in the Galactic bulge, remains difficult to explain. Additionally, one of the largest uncertainties in these studies is the unknown distance that positrons propagate before annihilation. In this paper, we will summarize the current knowledge base of Galactic positrons, and discuss how next-generation instruments could finally provide the answers.

Most violent and energetic processes in our universe, including mergers of compact objects, explosions of massive stars and extreme accretion events, produce copious amounts of X-rays. X-ray follow-up is an efficient tool for identifying transients because (1) X-rays can quickly localize transients with large error circles, and (2) X-rays reveal the nature of transients that may not have unique signatures at other wavelengths. In this white paper, we identify key science questions about several extragalactic multi-messenger and multi-wavelength transients, and demonstrate how X-ray follow-up helps answer these questions

Electromagnetic observations of the sky have been the basis for our study of the Universe for millennia, cosmic ray studies are now entering their second century, the first neutrinos from an astrophysical source were identified three decades ago, and gravitational waves were directly detected only four years ago. Detections of these messengers are now common. Astrophysics will undergo a revolution in the 2020s as multimessenger detections become routine. The 8th Astro2020 Thematic Area is Multimessenger Astronomy and Astrophysics, which includes the identification of the sources of gravitational waves, astrophysical and cosmogenic neutrinos, cosmic rays, and gamma-rays, and the coordinated multimessenger and multiwavelength follow-ups. Identifying and characterizing multimessenger sources enables science throughout and beyond astrophysics. Success in the multimessenger era requires: (i) sensitive coverage of the non-electromagnetic messengers, (ii) full coverage of the electromagnetic spectrum, with either fast-response observations or broad and deep high-cadence surveys, and (iii) improved collaboration, communication, and notification platforms.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

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 preliminary results from the Fermi High-Latitude Extended Sources Catalog (FHES), a comprehensive search for spatially extended gamma-ray sources at high Galactic latitudes ($|b|>5^\circ$) based on data from the Fermi Large Area Telescope (LAT). While the majority of high-latitude LAT sources are extragalactic blazars that appear point-like within the LAT angular resolution, there are several physics scenarios that predict the existence of populations of spatially extended sources. If Dark Matter consists of Weakly Interacting Massive Particles, the annihilation or decay of these particles in subhalos of the Milky Way would appear as a population of unassociated gamma-ray sources with finite angular extent. Gamma-ray emission from blazars could also be extended (so-called pair halos) due to the deflection of electron-positron pairs in the intergalactic magnetic field (IGMF). The pairs are produced in the absorption of gamma rays in the intergalactic medium and subsequently up-scatter photons of background radiation fields to gamma-ray energies. Measurement of pair halos could provide constraints on the strength and coherence length scale of the IGMF. In a dedicated search, we find 21 extended sources and 16 sources not previously characterized as extended. Limits on the flux of the extended source components are used to derive constraints on the strength of the IGMF using spectral and spatial templates derived from Monte Carlo simulations of electromagnetic cascades. This allows us to constrain the IGMF to be stronger than $3\times10^{-16}\,$G for a coherence length $\lambda \gtrsim 10\,$kpc.

BurstCube will detect long GRBs, attributed to the collapse of massive stars, short GRBs (sGRBs), resulting from binary neutron star mergers, as well as other gamma-ray transients in the energy range 10-1000 keV. sGRBs are of particular interest because they are predicted to be the counterparts of gravitational wave (GW) sources soon to be detectable by LIGO/Virgo. BurstCube contains 4 CsI scintillators coupled with arrays of compact low-power Silicon photomultipliers (SiPMs) on a 6U Dellingr bus, a flagship modular platform that is easily modifiable for a variety of 6U CubeSat architectures. BurstCube will complement existing facilities such as Swift and Fermi in the short term, and provide a means for GRB detection, localization, and characterization in the interim time before the next generation future gamma-ray mission flies, as well as space-qualify SiPMs and test technologies for future use on larger gamma-ray missions. The ultimate configuration of BurstCube is to have a set of $\sim10$ BurstCubes to provide all-sky coverage to GRBs for substantially lower cost than a full-scale mission.

Fermipy is an open-source python framework that facilitates analysis of data collected by the Fermi Large Area Telescope (LAT). Fermipy is built on the Fermi Science Tools, the publicly available software suite provided by NASA for the LAT mission. Fermipy provides a high-level interface for analyzing LAT data in a simple and reproducible way. The current feature set includes methods for extracting spectral energy distributions and lightcurves, generating test statistic maps, finding new source candidates, and fitting source position and extension. Fermipy leverages functionality from other scientific python packages including NumPy, SciPy, Matplotlib, and Astropy and is organized as a community-developed package following an open-source development model. We review the current functionality of Fermipy and plans for future development.

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.

We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for char- acterizing the upper limits across a large area of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, dif- ferences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914.

The nature of dark matter is a longstanding enigma of physics; it may consist of particles beyond the Standard Model that are still elusive to experiments. Among indirect search techniques, which look for stable products from the annihilation or decay of dark matter particles, or from axions coupling to high-energy photons, observations of the $\gamma$-ray sky have come to prominence over the last few years, because of the excellent sensitivity of the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope mission. The LAT energy range from 20 MeV to above 300 GeV is particularly well suited for searching for products of the interactions of dark matter particles. In this report we describe methods used to search for evidence of dark matter with the LAT, and review the status of searches performed with up to six years of LAT data. We also discuss the factors that determine the sensitivities of these searches, including the magnitudes of the signals and the relevant backgrounds, considering both statistical and systematic uncertainties. We project the expected sensitivities of each search method for 10 and 15 years of LAT data taking. In particular, we find that the sensitivity of searches targeting dwarf galaxies, which provide the best limits currently, will improve faster than the square root of observing time. Current LAT limits for dwarf galaxies using six years of data reach the thermal relic level for masses up to 120 GeV for the $b\bar{b}$ annihilation channel for reasonable dark matter density profiles. With projected discoveries of additional dwarfs, these limits could extend to about 250 GeV. With as much as 15 years of LAT data these searches would be sensitive to dark matter annihilations at the thermal relic cross section for masses to greater than 400 GeV (200 GeV) in the $b\bar{b}$ ($\tau^+ \tau^-$) annihilation channels.

The Small Magellanic Cloud (SMC) is the second-largest satellite galaxy of the Milky Way and is only 60 kpc away. As a nearby, massive, and dense object with relatively low astrophysical backgrounds, it is a natural target for dark matter indirect detection searches. In this work, we use six years of Pass 8 data from the Fermi Large Area Telescope to search for gamma-ray signals of dark matter annihilation in the SMC. Using data-driven fits to the gamma-ray backgrounds, and a combination of N-body simulations and direct measurements of rotation curves to estimate the SMC DM density profile, we found that the SMC was well described by standard astrophysical sources, and no signal from dark matter annihilation was detected. We set conservative upper limits on the dark matter annihilation cross section. These constraints are in agreement with stronger constraints set by searches in the Large Magellanic Cloud and approach the canonical thermal relic cross section at dark matter masses lower than 10 GeV in the $b\bar{b}$ and $\tau^+\tau^-$ channels.

We report on a search for monochromatic $\gamma$-ray features in the spectra of galaxy clusters observed by the \emph{Fermi} Large Area Telescope. Galaxy clusters are the largest structures in the Universe that are bound by dark matter (DM), making them an important testing ground for possible self-interactions or decays of the DM particles. Monochromatic $\gamma$-ray lines provide a unique signature due to the absence of astrophysical backgrounds and are as such considered a smoking-gun signature for new physics. An unbinned joint likelihood analysis of the sixteen most promising clusters using five years of data at energies between 10 and 400 GeV revealed no significant features. For the case of self-annihilation, we set upper limits on the monochromatic velocity-averaged interaction cross section. These limits are compatible with those obtained from observations of the Galactic Center, albeit weaker due to the larger distance to the studied clusters.

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 for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $\gamma$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3\deg that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $\gamma$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle \sigma v \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $\gamma$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$.

The gamma-ray energy range from a few hundred keV to a few hundred MeV has remained largely unexplored, mainly due to the challenging nature of the measurements, since the pi- oneering, but limited, observations by COMPTEL on the Compton Gamma-Ray Observatory (1991-2000). This energy range is a transition region between thermal and nonthermal processes, and accurate measurements are critical for answering a broad range of astrophysical questions. We are developing a MIDEX-scale wide-aperture discovery mission, ComPair (Compton-Pair Production Space Telescope), to investigate the energy range from 200 keV to > 500 MeV with high energy and angular resolution and with sensitivity approaching a factor of 20-50 better than COMPTEL. This instrument will be equally capable to detect both Compton-scattering events at lower energy and pair-production events at higher energy. ComPair will build on the her- itage of successful space missions including Fermi LAT, AGILE, AMS and PAMELA, and will utilize well-developed space-qualified detector technologies including Si-strip and CdZnTe-strip detectors, heavy inorganic scintillators, and plastic scintillators.

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

Above tens of GeV, gamma-ray observations with the Fermi Large Area Telescope (LAT) can be dominated by statistical uncertainties due to the low flux of sources and the limited acceptance. We are developing a new event class which can improve the acceptance: the "Calorimeter-only (CalOnly)" event class. The LAT has three detectors: the tracker, the calorimeter, and the anti-coincidence detector. While the conventional event classes require information from the tracker, the CalOnly event class is meant to be used when there is no usable tracker information. Although CalOnly events have poor angular resolution and a worse signal/background separation compared to those LAT events with usable tracker information, they can increase the instrument acceptance above few tens of GeV, where the performance of Fermi-LAT is limited by low photon statistics. In these proceedings we explain the concept and report some preliminary characteristics of this novel analysis.