(Abridged) The direct precise measurements of spectra of cosmic rays (CR) species in the wide energy range form a basis for propagation models, for interpretation of gamma-ray and microwave observations, and for disentangling possible signatures of new phenomena. Composition and spectra of CR species are vital for studies of galactic nucleosynthesis and high-energy processes in variety of environments and on different scales. Since the beginning of its operation, AMS-02 has delivered outstanding quality measurements of the spectra of CR antiprotons, electrons, positrons, and nuclei: H-O, Ne, Mg, Si. However, the spectra of heavier nuclei are expected to be available only later in the mission, which is currently extended through the end of the International Space Station (ISS) in 2030. Meanwhile, a comparison of published AMS-02 results with data from earlier experiments, such as HEAO-3-C2, indicate that earlier data may be affected by significant undocumented systematic errors. This creates an awkward situation: the new data are still incomplete, while using the old data to make up for the lack of AMS-02 measurements could result in significant errors. In this paper we are using the existing AMS-02 data to test HEAO-3-C2 results. We show that a fraction of HEAO-3-C2 data match available AMS-02 measurements quite well and can be used together with Voyager 1 and ACE-CRIS data to make predictions for the local interstellar spectra (LIS) of nuclei that are not yet released by the AMS-02 collaboration. We are also updating our already published LIS to provide a complete set of LIS from H-Ni in the energy range from 1 MeV/nucleon to ~100-500 TeV/nucleon thus covering 8-9 orders of magnitude in energy. Our calculations employ the GalProp-HelMod framework that is proved to be a reliable tool in deriving the LIS of CR antiprotons, electrons, and nuclei from H to O.

Cosmic rays (CRs) in the Galaxy are an important dynamical component of the interstellar medium (ISM) that interact with the other major components (interstellar magnetic and radiation fields, and gas) to produce broadband interstellar emissions that span the electromagnetic spectrum. The standard modelling of CR propagation and production of the associated emissions is based on a steady-state assumption, where the CR source spatial density is described using a smoothly varying function of position that does not evolve with time. While this is a convenient approximation, reality is otherwise where primary CRs are produced in and about highly localised regions, e.g., supernova remnants, which have finite lifetimes. In this paper we use the latest version of the GALPROP CR propagation code to model time-dependent CR injection and propagation through the ISM from a realistic three-dimensional discretised CR source density distribution, together with full three-dimensional models for the other major ISM components, and make predictions of the associated broadband non-thermal emissions. We compare the predictions for the discretised and equivalent steady-state model, finding that the former predicts novel features in the broadband non-thermal emissions that are absent for the steady-state case. Some of features predicted by the discretised model may be observable in all-sky observations made by WMAP and Planck, the recently launched eROSITA, the Fermi-LAT, and ground-based observations by HESS, HAWC, and the forthcoming CTA. The non-thermal emissions predicted by the discretised model may also provide explanations of puzzling anomalies in high-energy gamma ray data, such as the Fermi-LAT north/south asymmetry and residuals like the so-called "Fermi bubbles".

Local interstellar spectra (LIS) of secondary cosmic ray (CR) nuclei, lithium, beryllium, boron, and partially secondary nitrogen, are derived in the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. The lithium spectrum appears somewhat flatter at high energies compared to other secondary species that may imply a primary lithium component. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for the model-data comparison. The proposed LIS accommodate the low-energy interstellar spectra measured by Voyager 1, HEAO-3, and ACE/CRIS as well as the high-energy observations by PAMELA, AMS-02, and earlier experiments that are made deep in the heliosphere. The interstellar and heliospheric propagation parameters derived in this study are consistent with our earlier results for propagation of CR protons, helium, carbon, oxygen, antiprotons, and electrons.

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 High Altitude Water Cherenkov (HAWC) telescope recently observed extended emission around the Geminga and PSR~B0656+14 pulsar wind nebulae (PWNe). These observations have been used to estimate cosmic-ray (CR) diffusion coefficients near the PWNe that appear to be more than two orders of magnitude smaller than that typically derived for the interstellar medium from the measured abundances of secondary species in CRs. Two-zone diffusion models have been proposed as a solution to this discrepancy, where the slower diffusion zone (SDZ) is confined to a small region around the PWN. Such models are shown to successfully reproduce the HAWC observations of the Geminga PWN while retaining consistency with other CR data. It is found that the size of the SDZ influences the predicted positron flux and the spectral shape of the extended $\gamma$-ray emission at lower energies that can be observed with the {\it Fermi} Large Area Telescope ({\it Fermi} LAT). If the two observed PWNe are not unique, then it is likely that there are similar pockets of slow diffusion around many CR sources elsewhere in the Milky Way. The consequences of such picture for Galactic CR propagation is explored.

The current version of the Fermi Large Area Telescope data (P8R2) has been publicly available since June 2015, with the caveat that the residual background of all event classes, except ULTRACLEANVETO, was not fully isotropic: it was enhanced by a factor ~2 at 1-3 GeV within ~20 deg of the Ecliptic compared to the poles. By investigating the residual background using data only, we were able to find two sources of residual background: one due to non-interacting heavy ions and one due to cosmic-ray electrons leaking through the ribbons of the Anti-Coincidence Detector, the latter source being responsible for the background anisotropy. A set of simple cuts allows us to reject these events while losing less than 1% of the SOURCE class acceptance. This new selection has been used to produce a new version of the LAT data (P8R3).

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.

Identification of the cosmic-ray (CR) `PeVatrons', which are sources capable of accelerating particles to $\sim10^{15}$ eV energies and higher, may lead to resolving the long-standing question of the origin of the spectral feature in the all-particle CR spectrum known as the `knee'. Because CRs with these energies are deflected by interstellar magnetic fields identification of individual sources and determination of their spectral characteristics is more likely via very high energy $\gamma$-ray emissions, which provide the necessary directional information. However, pair production on the interstellar radiation field (ISRF) and cosmic microwave background leads to steepening of the high-energy tails of $\gamma$-ray spectra, and should be corrected for to enable true properties of the spectrum at source to be recovered. Employing recently developed three-dimensional ISRF models this paper quantifies the pair-absorption effect on spectra for sources in the Galactic centre direction at 8.5 kpc and 23.5 kpc distance, with the latter corresponding to the far side of the Galactic stellar disc where it is expected that discrimination of spectral features $>10$ TeV will be possible by the forthcoming Cherenkov Telescope Array. The estimates made suggest spectral cutoffs could be underestimated by factors of a few in the energy range so far sampled by TeV $\gamma$-ray telescopes. As an example to illustrate this, the recent HESS measurements of diffuse $\gamma$-ray emissions possibly associated with injection of CRs nearby Sgr A$^*$ are ISRF-corrected, and estimates of the spectral cutoff are re-evaluated. It is found that it could be higher by up to a factor $\sim 2$, indicating that these emissions may be consistent with a CR accelerator with a spectral cutoff of at least 1 PeV at the 95% confidence level.

Local interstellar spectra (LIS) of primary cosmic ray (CR) nuclei, such as helium, oxygen, and mostly primary carbon are derived for the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined into a single framework that is used to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. The developed iterative maximum-likelihood method uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The interstellar and heliospheric propagation parameters derived in this study are consistent with our prior analyses using the same methodology for propagation of CR protons, helium, antiprotons, and electrons. The resulting LIS accommodate a variety of measurements made in the local interstellar space (Voyager 1) and deep inside the heliosphere at low (ACE/CRIS, HEAO-3) and high energies (PAMELA, AMS-02).

Direct measurements of cosmic ray (CR) species combined with observations of their associated gamma-ray emissions can be used to constrain models of CR propagation, trace the structure of the Galaxy, and search for signatures of new physics. The spatial density distribution of the interstellar gas is a vital element for all these studies. So far models have employed the 2D cylindrically symmetric geometry, but their accuracy is well behind that of the available data. In this paper, 3D spatial density models for the neutral and molecular hydrogen are constructed based on empirical model fitting to gas line-survey data. The developed density models incorporate spiral arms and account for the warping of the disk, and the increasing gas scale height with radial distance from the Galactic center. They are employed together with the GALPROP CR propagation code to investigate how the new 3D gas models affect calculations of CR propagation and high-energy gamma-ray intensity maps. The calculations made reveal non-trivial features that are directly related to the new gas models. The best-fit values for propagation model parameters employing 3D gas models are presented and they differ significantly from the values derived with the 2D gas density models that have been widely used. The combination of 3D CR and gas density models provide a more realistic basis for the interpretation of non-thermal emissions from the Galaxy.

The local interstellar spectrum (LIS) of cosmic-ray (CR) electrons for the energy range 1 MeV to 1 TeV is derived using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The optimized HelMod parameters are then used to adjust GALPROP parameters to predict a refined LIS with the procedure repeated subject to a convergence criterion. The parameter optimization uses an extensive data set of proton spectra from 1997-2015. The proposed CR electron LIS accommodates both the low-energy interstellar spectra measured by Voyager 1 as well as the high-energy observations by PAMELA and AMS-02 that are made deep in the heliosphere; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The interstellar and heliospheric propagation parameters derived in this study agree well with our earlier results for CR protons, helium nuclei, and anti-protons propagation and LIS obtained in the same framework.

Information on cosmic-ray (CR) composition comes from direct CR measurements while their distribution in the Galaxy is evaluated from observations of their associated diffuse emission in the range from radio to gamma rays. Even though the main interaction processes are identified, more and more precise observations provide an opportunity to study more subtle effects and pose a challenge to the propagation models. GALPROP is a sophisticated CR propagation code that is being developed for about 20 years. It provides a unified framework for interpretations of data from many different types of experiments. It is used for a description of direct CR measurements and associated interstellar emissions (radio to gamma rays), thereby providing important information about CR injection and propagation in the interstellar medium. By accounting for all relevant observables at a time, the GALPROP code brings together theoretical predictions, interpretation of the most recent observations, and helps to reveal the signatures of new phenomena. In this paper we review latest applications of GALPROP and address ongoing and near future improvements. We are discussing effects of different propagation models, and of the transition from cylindrically symmetrical models to a proper 3D description of the components of the interstellar medium and the source distribution.

High-energy gamma rays of interstellar origin are produced by the interaction of cosmic-ray (CR) particles with the diffuse gas and radiation fields in the Galaxy. The main features of this emission are well understood and are reproduced by existing CR propagation models employing 2D Galactocentric cylindrically symmetrical geometry. However, the high-quality data from instruments like the Fermi Large Area Telescope reveal significant deviations from the model predictions on few to tens of degree scales indicating the need to include the details of the Galactic spiral structure and thus require 3D spatial modelling. In this paper the high-energy interstellar emissions from the Galaxy are calculated using the new release of the GALPROP code employing 3D spatial models for the CR source and interstellar radiation field (ISRF) densities. Three models for the spatial distribution of CR sources are used that are differentiated by their relative proportion of input luminosity attributed to the smooth disc or spiral arms. Two ISRF models are developed based on stellar and dust spatial density distributions taken from the literature that reproduce local near- to far-infrared observations. The interstellar emission models that include arms and bulges for the CR source and ISRF densities provide plausible physical interpretations for features found in the residual maps from high-energy gamma-ray data analysis. The 3D models for CR and ISRF densities provide a more realistic basis that can be used for the interpretation of the non-thermal interstellar emissions from the Galaxy.

Local interstellar spectra (LIS) for protons, helium and antiprotons are built using the most recent experimental results combined with the state-of-the-art models for propagation in the Galaxy and heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod that provides modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameters optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997-2015. The proposed LIS accommodate both the low energy interstellar CR spectra measured by Voyager 1 and the high energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.

Propagation of cosmic rays (CRs) from their sources to the observer is described mainly as plain diffusion at high energies, while at lower energies there are other physical processes involved, both in the interstellar space and in the heliosphere. The latter was a subject of considerable uncertainty until recently. New data obtained by several CR missions can be used to and the local interstellar spectra (LIS) of CR species that would significantly reduce the uncertainties associated with the heliospheric propagation. In this paper we present the LIS of CR protons and helium outside the heliospheric boundary. The proposed LIS are tuned to accommodate both, the low energy CR spectra measured by Voyager 1, and the high energy observations publicly released by BESS, Pamela, AMS-01 and AMS-02. The proton and helium LIS are derived by combining CR propagation in the Galaxy, as described by GALPROP, with the heliospheric modulation computed using the HelMod Monte Carlo Tool. The proposed LIS are tuned to reproduce the modulated spectra for both, high and low, levels of solar activity.

HelMod is a code evaluating the transport of Galactic cosmic rays through the inner heliosphere down to Earth. It is based on a 2-D Monte Carlo approach and includes a general description of the symmetric and antisymmetric parts of the diffusion tensor, thus, properly treating the particle drift effects. The model has been tuned in order to fit the data observed outside the ecliptic plane at several distances from the Earth and the spectra observed near the Earth for both, high and low solar activity levels. A stand-alone python module, fully compatible with GalProp, was developed for a comprehensive calculation of solar modulation effects, resulting in a newly suggested set of local interstellar spectra.

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.

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.

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.

While supernova remnants (SNRs) are widely thought to be powerful cosmic-ray accelerators, indirect evidence comes from a small number of well-studied cases. Here we systematically determine the gamma-ray emission detected by the Fermi Large Area Telescope (LAT) from all known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions containing SNRs, accounting for systematic uncertainties caused by source confusion and instrumental response. We have also developed a method to explore some systematic effects on SNR properties caused by the modeling of the interstellar emission (IEM). The IEM contributes substantially to gamma-ray emission in the regions where SNRs are located. To explore the systematics we consider different model construction methods, different model input parameters, and independently fit the model components to the gamma-ray data. We will describe this analysis method in detail. In the First Fermi LAT SNR Catalog there are 30 sources classified as SNRs, using spatial overlap with the radio position. For all the remaining regions we evaluated upper limits on SNRs' emission. In this work we will present a study of the aggregate characteristics of SNRs, such as comparisons between GeV and radio sizes as well as fluxes and spectral indexes and with TeV.

It is widely accepted that cosmic rays (CRs) up to at least PeV energies are Galactic in origin. Accelerated particles are injected into the interstellar medium where they propagate to the farthest reaches of the Milky Way, including a surrounding halo. The composition of CRs coming to the solar system can be measured directly and has been used to infer the details of CR propagation that are extrapolated to the whole Galaxy. In contrast, indirect methods, such as observations of gamma-ray emission from CR interactions with interstellar gas, have been employed to directly probe the CR densities in distant locations throughout the Galactic plane. In this article we use 73 months of data from the Fermi Large Area Telescope in the energy range between 300 MeV and 10 GeV to search for gamma-ray emission produced by CR interactions in several high- and intermediate-velocity clouds located at up to ~ 7 kpc above the Galactic plane. We achieve the first detection of intermediate-velocity clouds in gamma rays and set upper limits on the emission from the remaining targets, thereby tracing the distribution of CR nuclei in the halo for the first time. We find that the gamma-ray emissivity per H atom decreases with increasing distance from the plane at 97.5% confidence level. This corroborates the notion that CRs at the relevant energies originate in the Galactic disk. The emissivity of the upper intermediate-velocity Arch hints at a 50% decline of CR densities within 2 kpc from the plane. We compare our results to predictions of CR propagation models.

We present the analysis of Fermi Large Area Telescope (LAT) $\gamma$-ray observations of HB~21 (G89.0+4.7). We detect significant $\gamma$-ray emission associated with the remnant: the flux >100 MeV is $9.4\pm0.8(stat)\pm1.6(syst)\times10^{-11}$ erg cm$^{-2}$ s$^{-1}$. HB 21 is well modeled by a uniform disk centered at $l= 88{\deg}.75\pm 0{\deg}.04$, $b = +4{\deg}.65 \pm 0{\deg}.06$ with a radius of $1{\deg}.19 \pm 0{\deg}.06$. The $\gamma$-ray spectrum shows clear evidence of curvature, suggesting a cutoff or break in the underlying particle population at an energy of a few GeV. We complement $\gamma$-ray observations with the analysis of the WMAP 7-year data from 23 to 93 GHz, achieving the first detection of HB 21 at these frequencies. In combination with archival radio data, the radio spectrum shows a spectral break which helps to constrain the relativistic electron spectrum, hence parameters of simple non-thermal radiation models. In one-zone models multiwavelength data favor the origin of $\gamma$ rays from nucleon-nucleon collisions. A single population of electrons cannot produce both $\gamma$ rays through bremsstrahlung and radio emission through synchrotron radiation. A predominantly inverse-Compton origin of the $\gamma$-ray emission is disfavored because it requires lower interstellar densities than are inferred for HB 21. In the hadronic-dominated scenarios accelerated nuclei contribute a total energy of $\sim 3 \times10^{49}$ ergs, while in a two-zone bremsstrahlung-dominated scenario the total energy in accelerated particles is $\sim1\times10^{49}$ ergs.

Galactic interstellar emission contributes substantially to Fermi LAT observations in the Galactic plane, the location of the majority of supernova remnants (SNRs). To explore some systematic effects on SNRs' properties caused by interstellar emission modeling, we have developed a method comparing the official LAT interstellar emission model results to eight alternative models. We created the eight alternative Galactic interstellar models by varying a few input parameters to GALPROP, namely the height of the cosmic ray propagation halo, cosmic ray source distribution in the Galaxy, and atomic hydrogen spin temperature. We have analyzed eight representative SNRs chosen to encompass a range of Galactic locations, extensions, and spectral properties using the eight different interstellar emission models. We will present the results and method in detail and discuss the implications for studies such as the 1st Fermi LAT SNR Catalog.

The Sun and the Moon are quiescent gamma-ray sources that are clearly detectable in Fermi-LAT data. While moving through the sky, the Sun and the Moon can be a significant background in the analysis of Fermi-LAT data if they pass through the region of interest. Accurate modeling of their intensity is needed in this case, accounting for the correct exposure of their track along the sky. We present the Solar System Tools (SST) which calculate the templates of the Sun's and the Moon's intensity in the sky for a given observing period and a model of their steady emission. These tools are available in the standard Fermi-LAT Science Tools.

Galactic interstellar emission contributes substantially to Fermi LAT observations in the Galactic plane, where the majority of Supernova Remnants (SNRs) are located. We have developed a method to explore some systematic effects on SNRs' properties caused by interstellar emission modeling. We created eight alternative Galactic interstellar models by varying a few input parameters to GALPROP, namely the height of the cosmic ray propagation halo, cosmic ray source distribution in the Galaxy, and atomic hydrogen spin temperature. We have analyzed eight representative candidate SNRs chosen to encompass a range of Galactic locations, extensions, and spectral properties using the eight different interstellar emission models. The models were fitted to the LAT data with free independent normalization coefficients for the various components of the model along the line of sight in each region of interest. We will discuss the results and compare them with those obtained with the official LAT interstellar emission model.

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

Research in many areas of modern physics such as, e.g., indirect searches for dark matter and particle acceleration in SNR shocks, rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions, The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, that uses astrophysical information, nuclear and particle data as input to self-consistently predict CRs, gamma rays, synchrotron and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

We use the GALPROP code for cosmic-ray (CR) propagation to calculate the broad-band luminosity spectrum of the Milky Way related to CR propagation and interactions in the interstellar medium. This includes gamma-ray emission from the production and subsequent decay of neutral pions, bremsstrahlung, and inverse Compton scattering, and synchrotron radiation. The Galaxy is found to be nearly a CR electron calorimeter, but {\it only} if gamma ray emitting processes are taken into account. Synchrotron radiation alone accounts for only one third of the total electron energy losses with ~10-20% of the total synchrotron emission from secondary CR electrons and positrons. The relationship between far-infrared and radio luminosity that we find from our models is consistent with that found for galaxies in general. The results will be useful for understanding the connection between diffuse emissions from radio through gamma rays in ``normal'' (non-AGN dominated) galaxies, as well as for estimating the broad-band extragalactic diffuse background from these kinds of galaxies.

The diffuse high-energy gamma-ray emission of the Milky Way arises from interactions of cosmic-rays (CRs) with interstellar gas and radiation field in the Galaxy. The neutral hydrogen (H I) gas component is by far the most massive and broadly distributed component of the interstellar medium. Using the 21-cm emission line from the hyperfine structure transition of atomic hydrogen it is possible to determine the column density of H I if the spin temperature (Ts) of the emitting gas is known. Studies of diffuse gamma-ray emission have generally relied on the assumption of a fixed, constant spin temperature for all H I in the Milky Way. Unfortunately, observations of H I in absorption against bright background sources has shown it to vary greatly with location in the Milky Way. We will discuss methods for better handling of spin temperatures for Galactic diffuse emission modeling using the Fermi-LAT data and direct observation of the spin temperature using H I absorption.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

Aims: We present multiwavelength observations of one of the most distant gamma-ray bursts detected so far, GRB 080913. Based on these observations, we consider whether it could be classified as a short-duration GRB and discuss the implications for the progenitor nature and energy extraction mechanisms. Methods: Multiwavelength X-ray, near IR and millimetre observations were made between 20.7 hours and 16.8 days after the event. Results: Whereas a very faint afterglow was seen at the 3.5m CAHA telescope in the nIR, the X-ray afterglow was clearly detected in both Swift and XMM-Newton observations. An upper limit is reported in the mm range. We have modeled the data assuming a collimated $\theta_0$ $\gtrsim$ 3$^\circ$ blast wave with an energy injection at 0.5 days carrying $5\sim 10^{52}$ erg or approximately 12 times the initial energy of the blast wave. We find that GRB 080913 shares many of the gamma-ray diagnostics with the more recent burst GRB 090423 for being classified as short had they ocurred at low redshift. If the progenitor were a compact binary merger, it is likely composed by a NS and BH. The Blandford-Znajek (BZ) mechanism is the preferred one to extract energy from the central, maximally-rotating BH. Both the magnetic field close to the event horizon (B) and the BH mass ($M_{bh}$) are restricted within a relatively narrow range, such that $(B / 3\times 10^{16} \rm{G}) (M_{bh} / 7 M_\odot) \sim 1$. Similar constraints on the central BH hold for collapsar-like progenitor systems if the BZ-mechanism works for the system at hand.

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

The cosmic-ray propagation code GALPROP has recently been enhanced in the areas of Galactic magnetic field and synchrotron radiation. A full 3D treatment is now implemented. We describe these enhancements and show applications to the interpretation of non-thermal radio surveys and gamma rays, to probe the Galactic cosmic-ray electron spectrum and magnetic fields.

Aims. Understand the shape and implications of the multiband light curve of GRB 050408, an X-ray rich (XRR) burst. Methods. We present a multiband optical light curve, covering the time from the onset of the gamma-ray event to several months after, when we only detect the host galaxy. Together with X-ray, millimetre and radio observations we compile what, to our knowledge, is the most complete multiband coverage of an XRR burst afterglow to date. Results. The optical and X-ray light curve is characterised by an early flattening and an intense bump peaking around 6 days after the burst onset. We explain the former by an off-axis viewed jet, in agreement with the predictions made for XRR by some models, and the latter with an energy injection equivalent in intensity to the initial shock. The analysis of the spectral flux distribution reveals an extinction compatible with a low chemical enrichment surrounding the burst. Together with the detection of an underlying starburst host galaxy we can strengthen the link between XRR and classical long-duration bursts.

Since the discovery of the first short-hard gamma-ray burst afterglows in 2005, the handful of observed events have been found to be embedded in nearby (z < 1), bright underlying galaxies. We present multiwavelength observations of the short-duration burst GRB 060121, which is the first observed to clearly outshine its host galaxy (by a factor >10^2). A photometric redshift for this event places the progenitor at a most probable redshift of z = 4.6, with a less probable scenario of z = 1.7. In either case, GRB 060121 could be the farthermost short-duration GRB detected to date and implies an isotropic-equivalent energy release in gamma-rays comparable to that seen in long-duration bursts. We discuss the implications of the released energy on the nature of the progenitor. These results suggest that GRB 060121 may belong to a family of energetic short-duration events, lying at z > 1 and whose optical afterglows would outshine their host galaxies, unlike the first short-duration GRBs observed in 2005. The possibility of GRB 060121 being an intermediate duration burst is also discussed.

We extend the standard fireball model, widely used to interpret gamma-ray burst (GRB) afterglow light curves, to include energy injections, and apply the model to the afterglow light curves of GRB 990510, GRB 000301C and GRB 010222. We show that discrete energy injections can cause temporal variations in the optical light curves and present fits to the light curves of GRB 000301C as an example. A continuous injection may be required to interpret other bursts such as GRB 010222. The extended model accounts reasonably well for the observations in all bands ranging from X-rays to radio wavelengths. In some cases, the radio light curves indicate that additional model ingredients may be needed.

In this study optical/near-infrared(NIR) broad band photometry and optical spectroscopic observations of the GRB 030329 host galaxy are presented. The Spectral Energy Distribution (SED) of the host is consistent with a starburst galaxy template with a dominant stellar population age of ~150 Myr and an extinction Av ~0.6. Analysis of the spectral emission lines shows that the host is likely a low metallicity galaxy. Two independent diagnostics, based on the restframe UV continuum and the [OII] line flux, provide a consistent unextincted star formation rate of SFR ~0.6 Mo yr^-1. The low absolute magnitude of the host (M_B ~ -16.5) implies a high specific star formation rate value, SSFR = ~34 Mo yr^-1 (L/L*)^-1.

We report the result of a search for Ly-alpha emission from the host galaxies of the gamma-ray bursts (GRBs) 030226 (z = 1.986), 021004 (z = 2.335) and 020124 (z = 3.198). We find that the host galaxy of GRB 021004 is an extended (around 8 kpc) strong Ly-alpha emitter with a restframe equivalent width (EW) of 68^{+12}_{-11} AA, and a star-formation rate of 10.6 +/- 2.0 M_sun/yr. We do not detect the hosts of GRB 030226 and GRB 020124, but the upper limits on their Ly-alpha fluxes do not rule out large restframe EWs. In the fields of GRB 021004 and GRB 030226 we find seven and five other galaxies, respectively, with excess emission in the narrow-band filter. These galaxies are candidate Ly-alpha emitting galaxies in the environment of the host galaxies. We have also compiled a list of all z > 2 GRB hosts, and demonstrate that a scenario where they trace star formation in an unbiased way is compatible with current observational constraints. Fitting the z=3 luminosity function (LF) under this assumption, results in a characteristic luminosity of R* = 24.6 and a faint end slope of alpha = -1.55, consistent with the LF measured for Lyman-break galaxies.

Several GRB afterglow light curves deviate strongly from the power law decay observed in most bursts. We show that these variations can be accounted for by including refreshed shocks in the standard fireball model previously used to interpret the overall afterglow behavior. As an example we consider GRB 021004 that exhibited strong light curve variations and has a reasonably well time-resolved polarimetry. We show that the light curves in the R-band, X-rays and in the radio can be accounted for by four energy injection episodes in addition to the initial event. The polarization variations are shown to be a consequence of the injections.