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.

A critical review of the standard paradigm for the origin of Galactic cosmic rays is presented. Recent measurements of local and far-away cosmic rays reveal unexpected behaviours, which challenge the commonly accepted scenario. These recent findings are discussed, together with long-standing open issues. Despite the progress made thanks to ever-improving observational techniques and theoretical investigations, at present our understanding of the origin and of the behaviour of cosmic rays remains incomplete. We believe it is still unclear whether a modification of the standard paradigm, or rather a radical change of the paradigm itself is needed in order to interpret all the available data on cosmic rays within a self-consistent scenario.

COMPTEL was the Compton telescope on NASA's Compton Gamma Ray Observatory CGRO launched in April 1991 and which was re-entered in June 2000. COMPTEL covered the energy range 0.75 to 30 MeV, and performed a full-sky survey which is still unique in this range, with no followup mission yet approved. This remains a major uncharted region, and the heritage data from COMPTEL are still our main source of information. Data analysis has continued at MPE however, since the data were never fully analysed during the mission or in the period following, and improvements in analysis techniques and computer power make this possible.

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.

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.

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.

We investigate the impact of the recently published ISRF model of the Milky Way by Popescu et al. on the CR electrons and positrons in the context of cosmic-ray transport modelling. We also study predictions for diffuse Galactic gamma radiation and underline the importance of the increased ISRF intensities.We use the PICARD code for solving the CR transport equation to obtain predictions based on Galactic cosmic-ray electron fluxes. We show that the new ISRF yields gamma-ray intensity increases, most particular at the Galactic Center and in the Galactic Plane. The impact is largest at energies around 220 GeV.

We analyze the processes governing cosmic-ray (CR) penetration into molecular clouds and the resulting generation of gamma-ray emission. The density of CRs inside a cloud is depleted at lower energies due to the self-excited MHD turbulence. The depletion depends on the effective gas column density ("size") of the cloud. We consider two different environments where the depletion effect is expected to be observed. For the Central Molecular Zone, the expected range of CR energy depletion is $E\lesssim 10$ GeV, leading to the depletion of gamma-ray flux below $E_\gamma\approx 2$ GeV. This effect can be important for the interpretation of the GeV gamma-ray excess in the Galactic Center, which has been revealed from the standard model of CR propagation (assuming the CR spectrum inside a cloud to be equal to the interstellar spectrum). Furthermore, recent observations of some local molecular clouds suggest the depletion of the gamma-ray emission, indicating possible self-modulation of the penetrating low-energy CRs.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

As the next step toward an improved large scale Galactic magnetic field model, we present a simple comparison of polarised synchrotron and thermal dust emission on the Galactic plane. We find that the field configuration in our previous model that reproduces the polarised synchrotron is not compatible with the WMAP 94 GHz polarised emission data. In particular, the high degree of dust polarisation in the outer Galaxy (90deg < l < 270deg) implies that the fields in the dust-emitting regions are more ordered than the average of synchrotron-emitting regions. This new dust information allows us to constrain the spatial mixing of the coherent and random magnetic field components in the outer Galaxy. The inner Galaxy differs in polarisation degree and apparently requires a more complicated scenario than our current model. In the scenario that each interstellar component (including fields and now dust) follows a spiral arm modulation, as observed in external galaxies, the changing degree of ordering of the fields in dust-emitting regions may imply that the dust arms and the field component arms are shifted as a varying function of Galacto-centric radius. We discuss the implications for how the spiral arm compression affects the various components of the magnetised interstellar medium but conclude that improved data such as that expected from the Planck satellite will be required for a thorough analysis.

We extend previous work modeling the Galactic magnetic field in the plane using synchrotron emission in total and polarised intensity. In this work, we include a more realistic treatment of the cosmic-ray electrons using the GALPROP propagation code optimized to match the existing high-energy data. This addition reduces the degeneracies in our previous analysis and when combined with an additional observed synchrotron frequency allows us to study the low-energy end of the cosmic-ray electron spectrum in a way that has not previously been done. For a pure diffusion propagation, we find a low-energy injection spectrum slightly harder than generally assumed; for J(E) \propto E^{\alpha}, we find {\alpha} = -1.34 \pm 0.12, implying a very sharp break with the spectrum above a few GeV. This then predicts a synchrotron brightness temperature spectral index, {\beta}, on the Galactic plane that is -2.8 < {\beta} < -2.74 below a few GHz and -2.98 < {\beta} < -2.91 up to 23 GHz. We find that models including cosmic-ray re-acceleration processes appear to be incompatible with the synchrotron data.

Cosmic rays (CRs) propagate in the Milky Way and interact with the interstellar medium and magnetic fields. These interactions produce emissions that span the electromagnetic spectrum, and are an invaluable tool for understanding the intensities and spectra of CRs in distant regions, far beyond those probed by direct CR measurements. We present updates on the study of CR properties by combining multi-frequency observations of the interstellar emission and latest CR direct measurements with propagation models.

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.

The first evidence of the gamma-ray emission from the quiescent Sun was found in the archival EGRET data that was later confirmed by Fermi-LAT observations with high significance. This emission is produced by Galactic cosmic rays (CRs) penetrating the inner heliosphere and inter- acting with the solar atmosphere and optical photons. The solar emission is characterized by two spatially and spectrally distinct components: (i) disk emission due to the CR cascades in the solar atmosphere, and (ii) spatially extended inverse Compton (IC) emission due to the CR electrons scattering off of solar photons. The intensity of both components associated with Galactic CRs anti-correlate with the level of the solar activity being the brightest during solar minimum. In this paper we discuss updates of the models of the IC component of the emission based on CR measurements made at different levels of solar activity, and we make predictions for e- ASTROGAM and AMEGO, proposed low-energy gamma-ray missions.

INTEGRAL/SPI is a space-based coded mask telescope featuring a 19-element Germanium detector array for high-resolution gamma-ray spectroscopy, encapsulated in a scintillation detector assembly that provide a veto for background from charged particles. In space, cosmic rays irradiate spacecraft and instrument, which results in a large instrumental background from activation of those materials, and leads to deterioration of the charge collection properties of the Ge detectors. We use 13.5 years of INTEGRAL/SPI spectra for each detector and for each pointing of the satellite. Spectral fits provide us with details about separated background components. From the strongest background lines, we first determine how the spectral response changes with time. We then determine how the instrumental background components change in intensities and other characteristics, most-importantly their relative distribution among detectors. Spectral resolution of Ge detectors in space degrades with time, up to 15% within half a year, consistently for all detectors, and across the SPI energy range. Semi-annual annealing operations recover these losses, yet there is a small long-term degradation. The intensity of instrumental background varies anti-correlated to solar activity, in general. There are significant differences among different lines and with respect to continuum. Background lines are found to have a characteristic, well-defined and long-term consistent intensity ratio among detectors. We use this to categorise lines in groups of similar behaviour. The dataset of spectral-response and background parameters as fitted across the INTEGRAL mission allows studies of SPI spectral response and background behaviour in a broad perspective, and efficiently supports precision modelling of instrumental background.

The Picard code for the numerical solution of the Galactic cosmic ray propagation problem allows for high-resolution models that acknowledge the 3D structure of our Galaxy. Picard was used to determine diffuse gamma-ray emission of the Galaxy over the energy range from 100 MeV to 100 TeV. We discuss the impact of a cosmic-ray source distribution aligned with the Galactic spiral arms for a range of such spiral-arm models. As expected, the impact on the gamma-ray emission is most distinct in the inverse-Compton channel, where imprints of the spiral arms are visible and yield predictions that are no longer symmetric to the rotational axis of the Milkyway. We will illustrate these differences by a direct comparison to results from previous axially symmetric Galactic propagation models: we find differences in the gamma-ray flux both on global scales and on local scales related to the spiral arm tangents. We compare gamma-ray flux and spectra at on-arm vs. off-arm projections and characterize the differences to axially symmetric models.

We make quantitative estimates of the power supplied to the Galactic cosmic ray population by second-order Fermi acceleration in the interstellar medium, or as it is usually termed in cosmic ray propagation studies, diffusive reacceleration. Using recent results on the local interstellar spectrum from the Voyager missions we show that for parameter values, in particular the Alfv\'en speed, typically used in propagation codes such as Galprop to fit the B/C ratio, the power contributed by diffusive reacceleration is significant and can be of order 50\% of the total Galactic cosmic ray power. The implications for the damping of interstellar turbulence are briefly considered.

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

We extend the description of the isotropic and anisotropic random component of the small-scale magnetic field within the existing magnetic field model of the Milky Way from Jansson&Farrar, by including random realizations of the small-scale component. Using a magnetic-field power spectrum with Gaussian random fields, the NE2001 model for the thermal electrons and the Galactic cosmic-ray electron distribution from the current GALPROP model we derive full-sky maps for the total and polarized synchrotron intensity as well as the Faraday rotation-measure distribution. While previous work assumed that small-scale fluctuations average out along the line-of-sight or which only computed ensemble averages of random fields, we show that these fluctuations need to be carefully taken into account. Comparing with observational data we obtain not only good agreement with 408 MHz total and WMAP7 22 GHz polarized intensity emission maps, but also an improved agreement with Galactic foreground rotation-measure maps and power spectra, whose amplitude and shape strongly depend on the parameters of the random field. We demonstrate that a correlation length of ~220pc (50pc being a 5\sigma lower limit) is needed to match the slope of the observed power spectrum of Galactic foreground rotation-measure maps. Using multiple realizations allows us also to infer errors on individual observables. We find that previously-used amplitudes for random and anisotropic random magnetic field components need to be rescaled by factors of ~0.3 and 0.6 to account for the new small-scale contributions. Our model predicts a rotation measure of -2.8\pm7.1 rad/m^2 and 4.4\pm11.0 rad/m^2 for the north and south Galactic poles respectively, in good agreement with observations. Applying our model to deflections of ultra-high-energy cosmic rays we infer a mean deflection of ~3.5\pm1.1 degree for 60 EeV protons arriving from CenA.

Recent models for the large-scale Galactic magnetic fields in the literature have been largely constrained by synchrotron emission and Faraday rotation measures. We use three different but representative models to compare their predicted polarized synchrotron and dust emission with that measured by the Planck satellite. We first update these models to match the Planck synchrotron products using a common model for the cosmic-ray leptons. We discuss the impact on this analysis of the ongoing problems of component separation in the Planck microwave bands and of the uncertain cosmic-ray spectrum. In particular, the inferred degree of ordering in the magnetic fields is sensitive to these systematic uncertainties, and we further show the importance of considering the expected variations in the observables in addition to their mean morphology. We then compare the resulting simulated emission to the observed dust polarization and find that the dust predictions do not match the morphology in the Planck data but underpredict the dust polarization away from the plane. We modify one of the models to roughly match both observables at high latitudes by increasing the field ordering in the thin disc near the observer. Though this specific analysis is dependent on the component separation issues, we present the improved model as a proof of concept for how these studies can be advanced in future using complementary information from ongoing and planned observational projects.

Microquasars are stellar-mass black holes accreting matter from a companion star and ejecting plasma jets at almost the speed of light. They are analogues of quasars that contain supermassive black holes of $10^6$ to $10^{10}$ solar masses. Accretion in microquasars varies on much shorter timescales than in quasars and occasionally produces exceptionally bright X-ray flares. How the flares are produced is unclear, as is the mechanism for launching the relativistic jets and their composition. An emission line near 511 kiloelectronvolts has long been sought in the emission spectrum of microquasars as evidence for the expected electron-positron plasma. Transient high-energy spectral features have been reported in two objects, but their positron interpretation remains contentious. Here we report observations of $\gamma$-ray emission from the microquasar V404 Cygni during a recent period of strong flaring activity. The emission spectrum around 511 kiloelectronvolts shows clear signatures of variable positron annihilation, which implies a high rate of positron production. This supports the earlier conjecture that microquasars may be the main sources of the electron-positron plasma responsible for the bright diffuse emission of annihilation $\gamma$-rays in the bulge region of our Galaxy. Additionally, microquasars could be the origin of the observed megaelectronvolt continuum excess in the inner Galaxy.

We present the results of the most complete ever scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine learning package. This is the first such study to separate out low-mass isotopes ($p$, $\bar p$ and He) from the usual light elements (Be, B, C, N, O). We find that the propagation parameters that best fit $p$, $\bar p$, He data are significantly different from those that fit light elements, including the B/C and $^{10}$Be/$^9$Be secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center.

The annihilation of positrons in the Galaxy's interstellar medium produces characteristic gamma-rays with a line at 511 keV. This emission has been observed with the spectrometer SPI on INTEGRAL, confirming a puzzling morphology with bright emission from an extended bulge-like region, and faint disk emission. Most plausible sources of positrons are believed to be distributed throughout the disk of the Galaxy. We aim to constrain characteristic spectral shapes for different spatial components in the disk and bulge with the high-resolution gamma-ray spectrometer SPI, based on a new instrumental background method and detailed multi-component sky model fitting. We confirm the detection of the main extended components of characteristic annihilation gamma-ray signatures at 58$\sigma$ significance in the line. The total Galactic line intensity amounts to $(2.7\pm0.3)\times10^{-3}~\mathrm{ph~cm^{-2}~s^{-1}}$ for our assumed spatial model. We derive spectra for the bulge and disk, and a central point-like and at the position of Sgr A*, and discuss spectral differences. The bulge shows a 511 keV line intensity of $(0.96\pm0.07)\times10^{-3}~\mathrm{ph~cm^{-2}~s^{-1}}$ and an ortho-positronium continuum equivalent to a positronium fraction of ($1.08\pm0.03$). The 2D-Gaussian representing the disk has an extent of $60^{+10}_{-5}$ deg in longitude and a latitudinal extent of $10.5^{+2.5}_{-1.5}$ deg; the line intensity is $(1.7\pm0.4)\times10^{-3}~\mathrm{ph~cm^{-2}~s^{-1}}$ with a marginal detection of ortho-positronium and an overall diffuse Galactic continuum of $(5.85\pm1.05)\times10^{-5}~\mathrm{ph~cm^{-2}~s^{-1}~keV^{-1}}$ at 511 keV. The disk shows no flux asymmetry between positive and negative longitudes, although spectral details differ. The flux ratio between bulge and disk is ($0.6\pm0.1$). The central source has an intensity of $(0.8\pm0.2)\times10^{-4}~\mathrm{ph~cm^{-2}~s^{-1}}$.

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.

How a galaxy regulates its SNe energy into different interstellar/circumgalactic medium components strongly affects galaxy evolution. Based on the JVLA D-configuration C- (6 GHz) and L-band (1.6 GHz) continuum observations, we perform statistical analysis comparing multi-wavelength properties of the CHANG-ES galaxies. The high-quality JVLA data and edge-on orientation enable us for the first time to include the halo into the energy budget for a complete radio-flux-limited sample. We find tight correlations of $L_{\rm radio}$ with the mid-IR-based SFR. The normalization of our $I_{\rm 1.6GHz}/{\rm W~Hz^{-1}}-{\rm SFR}$ relation is $\sim$2-3 times of those obtained for face-on galaxies, probably a result of enhanced IR extinction at high inclination. We also find tight correlations between $L_{\rm radio}$ and the SNe energy injection rate $\dot{E}_{\rm SN(Ia+CC)}$, indicating the energy loss via synchrotron radio continuum accounts for $\sim0.1\%$ of $\dot{E}_{\rm SN}$, comparable to the energy contained in CR electrons. The integrated C-to-L-band spectral index is $\alpha\sim0.5-1.1$ for non-AGN galaxies, indicating a dominance by the diffuse synchrotron component. The low-scatter $L_{\rm radio}-{\rm SFR}$/$L_{\rm radio}-\dot{E}_{\rm SN (Ia+CC)}$ relationships have super-linear logarithmic slopes at $\sim2~\sigma$ in L-band ($1.132\pm0.067$/$1.175\pm0.102$) while consistent with linear in C-band ($1.057\pm0.075$/$1.100\pm0.123$). The super-linearity could be naturally reproduced with non-calorimeter models for galaxy disks. Using Chandra halo X-ray measurements, we find sub-linear $L_{\rm X}-L_{\rm radio}$ relations. These results indicate that the observed radio halo of a starburst galaxy is close to electron calorimeter, and a galaxy with higher SFR tends to distribute an increased fraction of SNe energy into radio emission (than X-ray).

We present results obtained from our newly developed Galactic cosmic-ray transport code PICARD, that solves the cosmic-ray transport equation. This code allows for the computation of cosmic-ray spectra and the resulting gamma-ray emission. Relying on contemporary numerical solvers allows for efficient computation of models with deca-parsec resolution. PICARD can handle locally anisotropic spatial diffusion acknowledging a full diffusion tensor. We used this framework to investigate the transition from axisymmetric to spiral-arm cosmic-ray source distributions. Wherever possible we compare model predictions with constraining observables in cosmic-ray astrophysics.

We search for evidence of dark matter (DM) annihilation in the isotropic gamma-ray background (IGRB) measured with 50 months of Fermi Large Area Telescope (LAT) observations. An improved theoretical description of the cosmological DM annihilation signal, based on two complementary techniques and assuming generic weakly interacting massive particle (WIMP) properties, renders more precise predictions compared to previous work. More specifically, we estimate the cosmologically-induced gamma-ray intensity to have an uncertainty of a factor ~20 in canonical setups. We consistently include both the Galactic and extragalactic signals under the same theoretical framework, and study the impact of the former on the IGRB spectrum derivation. We find no evidence for a DM signal and we set limits on the DM-induced isotropic gamma-ray signal. Our limits are competitive for DM particle masses up to tens of TeV and, indeed, are the strongest limits derived from Fermi LAT data at TeV energies. This is possible thanks to the new Fermi LAT IGRB measurement, which now extends up to an energy of 820 GeV. We quantify uncertainties in detail and show the potential this type of search offers for testing the WIMP paradigm with a complementary and truly cosmological probe of DM particle signals.

We present the first part of the observations made for the Continuum Halos in Nearby Galaxies, an EVLA Survey (CHANG-ES) project. The aim of the CHANG-ES project is to study and characterize the nature of radio halos, their prevalence as well as their magnetic fields, and the cosmic rays illuminating these fields. This paper reports observations with the compact D configuration of the Karl G. Jansky Very Large Array (VLA) for the sample of 35 nearby edge-on galaxies of CHANG-ES. With the new wide bandwidth capabilities of the VLA, an unprecedented sensitivity was achieved for all polarization products. The beam resolution is an average of 9.6" and 36" with noise levels reaching approximately 6 and 30 microJy per beam for C- and L-bands, respectively (robust weighting). We present intensity maps in these two frequency bands (C and L), with different weightings, as well as spectral index maps, polarization maps, and new measurements of star formation rates (SFRs). The data products described herein are available to the public in the CHANG-ES data release available at www.queensu.ca/changes. We also present evidence of a trend among galaxies with larger halos having higher SFR surface density, and we show, for the first time, a radio continuum image of the median galaxy, taking advantage of the collective signal-to-noise ratio of 30 of our galaxies. This image shows clearly that a typical spiral galaxy is surrounded by a halo of magnetic fields and cosmic rays.

Cosmic-ray scattering on magnetic turbulence leads to spatial diffusive propagation; if the scattering medium is moving, this will inevitably also cause changes in the momentum of the particles, so-called diffusive reacceleration. This can be described as diffusion in momentum space. Diffusive reacceleration has often been invoked to explain the peak observed in secondary-to-primary ratios at a few GeV, in particular Boron-to-Carbon. This avoids the necessity to postulate an ad-hoc break in the spatial diffusive coefficient, and has become almost a standard in modelling cosmic-ray spectra. However, at the levels invoked, the process implies a significant input of energy from the interstellar medium into cosmic rays, so that in such models interstellar space competes with the usual accelerators like supernova remnants. The questions arise: is reacceleration really occurring at the high level required to explain secondary-to-primary ratios? and are the energy requirements physically plausible? We address this issue using both analytical and numerical models of cosmic-ray propagation.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

Some recent extensions to the GALPROP cosmic-ray propagation package are described. The enhancements include: an accurate solution option, improved convection formulation, alternative spatial boundary conditions, polarized synchrotron emission, new magnetic field models, updated gamma-ray production cross-sections, free-free radio emission and absorption, primary positrons, additional injection spectral breaks, deuterium production by pp fusion, hadronic energy losses, improved HEALPix skymap format, compatibility with latest HEALPix release, and various bug fixes. The Explanatory Supplement has been extensively updated, including details of these enhancements. A compatible plot package GALPLOT for GALPROP output is also provided, as well as other related software.

Cosmic rays (CRs) interact with the gas, the radiation field and the magnetic field in the Milky Way, producing diffuse emission from radio to gamma rays. Observations of this diffuse emission and comparison with detailed predictions are powerful tools to unveil the CR properties and to study CR propagation. We present various GALPROP CR propagation scenarios based on current CR measurements. The predicted synchrotron emission is compared to radio surveys, and synchrotron temperature maps from WMAP and Planck, while the predicted interstellar gamma-ray emission is compared to Fermi-LAT observations. We show how multi-wavelength observations of the Galactic diffuse emission can be used to help constrain the CR lepton spectrum and propagation. Finally we discuss how radio and microwave data could be used in understanding the diffuse Galactic gamma-ray emission observed with Fermi-LAT, especially at low energies.

Precise gamma-ray emissivities from cosmic-ray interactions with interstellar gas have been recently derived using Fermi-LAT data, and used to constrain the local interstellar spectra of protons and leptons. We report on a continuing effort to exploit these emissivities combined with the latest hadronic gamma-ray production cross-sections and other constraints such as synchrotron emission for the leptonic component. The interstellar spectra provide important information for heliospheric modulation, and cosmic-ray origin and propagation.

Planck data when combined with ancillary data provide a unique opportunity to separate the diffuse emission components of the inner Galaxy. The purpose of the paper is to elucidate the morphology of the various emission components in the strong star-formation region lying inside the solar radius and to clarify the relationship between the various components. The region of the Galactic plane covered is l=300-0-60deg where star-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1deg and 2deg, which correspond to FWHM z-widths of 100-200pc at a typical distance of 6kpc. The four emission components studied here are synchrotron, free-free, anomalous microwave emission (AME), and thermal (vibrational) dust emission. These components are identified by constructing spectral energy distributions (SEDs) at positions along the Galactic plane using the wide frequency coverage of Planck (28.4-857GHz) in combination with low-frequency radio data at 0.408-2.3GHz plus WMAP data at 23-94GHz, along with far-infrared (FIR) data from DIRBE and IRAS. The free-free component is determined from radio recombination line (RRL) data. AME is found to be comparable in brightness to the free-free emission on the Galactic plane in the frequency range 20-40GHz with a width in latitude similar to that of the thermal dust; it comprises 45+/-1% of the total 28.4GHz emission in the longitude range l=300-0-60deg. The free-free component is the narrowest, reflecting the fact that it is produced by current star-formation as traced by the narrow distribution of OB stars. It is the dominant emission on the plane between 60 and 100GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM ~100pc...(abridged)

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 study the impact of possible spiral-arm distributions of Galactic cosmic-ray sources on the flux of various cosmic-ray nuclei throughout our Galaxy. We investigate model cosmic-ray spectra at the nominal position of the sun and at different positions within the Galaxy. The modelling is performed using the recently introduced numerical cosmic ray propagation code \textsc{Picard}. Assuming non-axisymmetric cosmic ray source distributions yields new insights on the behaviour of primary versus secondary nuclei. We find that primary cosmic rays are more strongly confined to the vicinity of the sources, while the distribution of secondary cosmic rays is much more homogeneous compared to the primaries. This leads to stronger spatial variation in secondary to primary ratios when compared to axisymmetric source distribution models. A good fit to the cosmic-ray data at Earth can be accomplished in different spiral-arm models, although leading to decisively different spatial distributions of the cosmic-ray flux. This results in very different cosmic ray anisotropies, where even a good fit to the data becomes possible. Consequently, we advocate directions to seek best fit propagation parameters that take into account the higher complexity introduced by the spiral-arm structure on the cosmic-ray distribution. We specifically investigate whether the flux at Earth is representative for a large fraction of the Galaxy. The variance among possible spiral-arm models allows us to quantify the spatial variation of the cosmic-ray flux within the Galaxy in presence of non-axisymmetric source distributions.

Shortened abstract: Observations of the nearby Chamaeleon clouds in gamma rays with the Fermi Large Area Telescope and in thermal dust emission with Planck and IRAS have been used with the HI and CO radio data to (i) map the gas column densities in the different phases and at the dark neutral medium (DNM) transition between the HI-bright and CO-bright media; (ii) constrain the CO-to-$H_2$ conversion factor, $X_{CO}$; (iii) probe the dust properties per gas nucleon in each gas phase and spatially across the clouds. We have separated clouds in velocity in HI and CO emission and modelled the 0.4-100 GeV intensity, the dust optical depth at 353 GHz, the thermal radiance of the large grains, and an estimate of the dust extinction empirically corrected for the starlight intensity, $A_{VQ}$. The gamma-ray emissivity spectra confirm that the GeV-TeV cosmic rays uniformly permeate all gas phases up to the CO cores. The dust and cosmic rays reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the HI-DNM-CO transitions and CO-dark $H_2$ fractions for separate clouds. The corrected extinction provides the best fit to the total gas traced by the gamma rays, but we find evidence for a rise in $A_{VQ}/N_H$ and a steep rise in opacity, with increasing $N_H$ and $H_2$ fraction, and with decreasing dust temperature. We observe less variations for the specific power of the grains, except for a decline by half in the CO cores. This combined information suggests grain evolution. The gamma rays and dust radiance yield consistent $X_{CO}$ estimates near $0.7\times10^{20}$ cm$^{-2}$ (K km/s)$^{-1}$. The other dust tracers yield biased values because of the grain opacity rise in the CO clouds. These results also confirm a factor of 2 difference between $X_{CO}$ estimates at pc and kpc scales.

Planck has mapped the microwave sky in nine frequency bands between 30 and 857 GHz in temperature and seven bands between 30 and 353 GHz in polarization. In this paper we consider the problem of diffuse astrophysical component separation, and process these maps within a Bayesian framework to derive a consistent set of full-sky astrophysical component maps. For the temperature analysis, we combine the Planck observations with the 9-year WMAP sky maps and the Haslam et al. 408 MHz map to derive a joint model of CMB, synchrotron, free-free, spinning dust, CO, line emission in the 94 and 100 GHz channels, and thermal dust emission. Full-sky maps are provided with angular resolutions varying between 7.5 arcmin and 1 deg. Global parameters (monopoles, dipoles, relative calibration, and bandpass errors) are fitted jointly with the sky model, and best-fit values are tabulated. For polarization, the model includes CMB, synchrotron, and thermal dust emission. These models provide excellent fits to the observed data, with rms temperature residuals smaller than 4 uK over 93% of the sky for all Planck frequencies up to 353 GHz, and fractional errors smaller than 1% in the remaining 7% of the sky. The main limitations of the temperature model at the lower frequencies are degeneracies among the spinning dust, free-free, and synchrotron components; additional observations from external low-frequency experiments will be essential to break these. The main limitations of the temperature model at the higher frequencies are uncertainties in the 545 and 857 GHz calibration and zero-points. For polarization, the main outstanding issues are instrumental systematics in the 100-353 GHz bands on large angular scales in the form of temperature-to-polarization leakage, uncertainties in the analog-to-digital conversion, and very long time constant corrections, all of which are expected to improve in the near future.

The origin of the extragalactic $\gamma$-ray background (EGB) has been debated for some time. { The EGB comprises the $\gamma$-ray emission from resolved and unresolved extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, as well as radiation from truly diffuse processes.} This letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the {\it Fermi} Large Area Telescope (LAT). We show that blazars account for 50$^{+12}_{-11}$\,\% of the EGB photons ($>$0.1\,GeV), and that {\it Fermi}-LAT has already resolved $\sim$70\,\% of this contribution. Blazars, and in particular low-luminosity hard-spectrum nearby sources like BL Lacs, are responsible for most of the EGB emission above 100\,GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cut-off observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1--820\,GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark-matter annihilation cross section.

Observations of the properties of dense molecular clouds are critical in understanding the process of star-formation. One of the most important, but least understood, is the role of the magnetic fields. We discuss the possibility of using high-resolution, high-sensitivity radio observations with the SKA to measure for the first time the in-situ synchrotron radiation from these molecular clouds. If the cosmic-ray (CR) particles penetrate clouds as expected, then we can measure the B-field strength directly using radio data. So far, this signature has never been detected from the collapsing clouds themselves and would be a unique probe of the magnetic field. Dense cores are typically ~0.05 pc in size, corresponding to ~arcsec at ~kpc distances, and flux density estimates are ~mJy at 1 GHz. The SKA should be able to readily detect directly, for the first time, along lines-of-sight that are not contaminated by thermal emission or complex foreground/background synchrotron emission. Polarised synchrotron may also be detectable providing additional information about the regular/turbulent fields.

Observations of the properties of dense molecular clouds are critical in understanding the process of star-formation. One of the most important, but least understood, is the role of the magnetic fields. We discuss the possibility of using high-resolution, high-sensitivity radio observations to measure the in-situ synchrotron radiation from these molecular clouds. If the cosmic-ray particles penetrate clouds as expected, then we can measure the B-field strength directly using radio data. So far, this signature has never been detected from the collapsing clouds themselves and would be a unique probe of the magnetic field. Dense cores are typically ~0.05 pc in size, corresponding to ~arcsec at ~kpc distances, and flux density estimates are ~mJy at 1 GHz. They should be detectable, for example with the Square Kilometre Array.

There is evidence that the distribution of suspected cosmic ray sources are associated with the spiral arm structure of galaxies. It is yet not clearly understood what effect such a cosmic ray source distribution has on the particle transport in our Galaxy. We use the PICARD code to perform high-resolution 3D simulations of electrons and protons in galactic propagation scenarios that include four-arm and two-arm logarithmic spiral cosmic ray source distributions with and without a central bar structure as well as the spiral arm configuration of the NE2001 model for the distribution of free electrons in the Milky Way. Results of these simulation are compared to an axisymmetric radial source distribution. Also, effects on the cosmic ray flux and spectra due to different positions of the Earth relative to the spiral structure are studied. We find that high energy electrons are strongly confined to their sources and the obtained spectra largely depend on the Earth's position relative to the spiral arms. Similar finding have been obtained for low energy protons and electrons albeit at smaller magnitude. We find that even fractional contributions of a spiral arm component to the total cosmic ray source distribution influences the spectra on the Earth. This is apparent when compared to an axisymmetric radial source distribution as well as with respect to the Earth's position relative to the spiral arm structure. We demonstrate that the presence of a Galactic bar manifests itself as an overall excess of low energy electrons at the Earth. Using a spiral arm geometry as a cosmic ray source distributions offers a genuine new quality of modelling and is used to explain features in cosmic ray spectra at the Earth that are else-wise attributed to other propagation effects. We show that realistic cosmic ray propagation scenarios have to acknowledge non-axisymmetric source distributions.[Abbreviated]

The {\gamma}-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse {\gamma}-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

We present a detailed investigation of the $\gamma$-ray emission in the vicinity of the supernova remnant (SNR) W28 (G6.4$-$0.1) observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. We detected significant $\gamma$-ray emission spatially coincident with TeV sources HESS J1800$-$240A, B, and C, located outside the radio boundary of the SNR. Their spectra in the 2-100 GeV band are consistent with the extrapolation of the power-law spectra of the TeV sources. We also identified a new source of GeV emission, dubbed Source W, which lies outside the boundary of TeV sources and coincides with radio emission from the western part of W28. All of the GeV $\gamma$-ray sources overlap with molecular clouds in the velocity range from 0 to 20 km s$^{-1}$. Under the assumption that the $\gamma$-ray emission towards HESS J1800-240A, B, and C comes from $\pi^0$ decay due to the interaction between the molecular clouds and cosmic rays (CRs) escaping from W28, they can be naturally explained by a single model in which the CR diffusion coefficient is smaller than the theoretical expectation in the interstellar space. The total energy of the CRs escaping from W28 is constrained through the same modeling to be larger than $\sim$ 2 $\times$ 10$^{49}$ erg. The emission from Source W can also be explained with the same CR escape scenario.

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $\gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $\gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7\sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $\gamma$-ray flux from individual clusters in our sample.

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $\gamma$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

Cosmic-ray (CR) leptons produce radio synchrotron radiation by gyrating in interstellar magnetic fields (B-field). Details of B-fields, CR electron distributions and propagation are still uncertain. We present developments in our modelling of Galactic radio emission with the GALPROP code. It now includes calculations of radio polarization, absorption, and free-free emission. Total and polarized synchrotron emission are investigated in the context of physical model of CR propagation. Predictions are compared with radio data from 22 MHz to 2.3 GHz, and Wilkinson Microwave Anisotropy Probe data at 23 GHz. Spatial and spectral effects on the synchrotron modelling with different CR distribution, propagation halo size and CR propagation models are presented. We find that all-sky total intensity and polarization maps are reasonably reproduced by including an anisotropic B-field, with comparable intensity to the regular one defined by rotation measures. A halo size of 10 kpc, which is larger than usually assumed, is favoured. This work provides a basis for further studies on foreground emission with the Planck satellite and on interstellar gamma-ray emission with Fermi-Large Area Telescope.

We report the discovery of gamma-ray emission from the Circinus galaxy using the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. Circinus is a nearby (~4 Mpc) starburst with a heavily obscured Seyfert-type active nucleus, bipolar radio lobes perpendicular to the spiral disk, and kpc-scale jet-like structures. Our analysis of 0.1-100 GeV events collected during 4 years of LAT observations reveals a significant (~ 7.3 sigma) excess above the background. We find no indications of variability or spatial extension beyond the LAT point-spread function. A power-law model used to describe the 0.1-100 GeV gamma-ray spectrum yields a flux of (18.8+/-5.8)x10^{-9} ph cm^{-2} s^{-1} and photon index 2.19+/-0.12, corresponding to an isotropic gamma-ray luminosity of 3 x 10^{40} erg s^{-1}. This observed gamma-ray luminosity exceeds the luminosity expected from cosmic-ray interactions in the interstellar medium and inverse Compton radiation from the radio lobes. Thus the origin of the GeV excess requires further investigation.

We show the necessity for a new approach towards comprehensive and consistent simulations of the propagation of galactic cosmic rays. Our developments are optimised for addressing the spatially 3-dimensional inhomogeneous diffusion problem and utilise contemporary numerical methods. We aim to address the transport problem in a full 3-dimensional environment. For that, we test the transition from 2D to 3D simulation results within an existing propagation code. We present sub-kpc scale simulations that allow the investigation of small- scale structures regarding different model conditions such as variety regarding non-axisymmetric cosmic ray source distributions. These results are discussed critically and motivate our development of a new transport code for galactic cosmic rays. The capabilities of this code are outlined.

StellarICs is software to compute gamma-ray emission from inverse-Compton scattering by cosmic-ray leptons in the heliosphere and in the photospheres of stars. It includes a set of cosmic-ray spectra and a formulation of their modulation, but it can be used for any user-defined modulation model and lepton spectra. Inverse-Compton emission provides a unique probe of cosmic-ray leptons in the heliosphere between the earth and the sun. The software is publicly available and it is under continuing development.

Interstellar synchrotron emission depends on Galactic magnetic fields and on cosmic-ray leptons. Observations of radio emission are an important tool for studying cosmic-ray propagation models and interstellar electron spectrum and distribution in the Galaxy. We present the latest developments in our modeling of Galactic synchrotron emission with the GALPROP code, including polarization, absorption, and free-free emission. Using surveys over a wide range of radio frequencies and polarization measurements, we derive constraints on the low-energy interstellar cosmic-ray electron spectrum, magnetic fields and cosmic-ray propagation models. This work is of interest for studies of interstellar gamma-ray emission with Fermi-LAT, and synchrotron for the Planck mission.

More than 90% of the Galactic gas-related gamma-ray emissivity above 1 GeV is attributed to the decay of neutral pions formed in collisions between cosmic rays and interstellar matter, with lepton-induced processes becoming increasingly important below 1 GeV. Given the high-quality measurements of the gamma-ray emissivity of local interstellar gas between ~50 MeV and ~4 GeV obtained with the Large Area Telescope on board the Fermi space observatory, it is timely to re-investigate this topic in detail, including the hadronic production mechanisms. The emissivity spectrum will allow the interstellar cosmic-ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input to solar modulation models. A method for such an analysis and illustrative results are presented.

Secondary nuclear production physics is receiving increased attention given the high-quality measurements of the gamma-ray emissivity of local interstellar gas between ~50 MeV and ~40 GeV, obtained with the Large Area Telescope on board the Fermi space observatory. More than 90% of the gas-related emissivity above 1 GeV is attributed to gamma-rays from the decay of neutral pions formed in collisions between cosmic rays and interstellar matter, with lepton-induced processes becoming increasingly important below 1 GeV. The elementary kinematics of neutral pion production and decay are re-examined in light of two physics questions: does isobaric production follow a scaling behavior? and what is the minimum proton kinetic energy needed to make a gamma-ray of a certain energy formed through intermediate pi0 production? The emissivity spectrum will allow the interstellar cosmic-ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input to solar modulation models. A method for such an analysis and illustrative results are presented.

We discuss the relation between cosmic-ray induced interstellar gamma-ray and radio emission from our Galaxy, and emphasize the importance of their parallel study. We give an overview of results from Fermi-LAT on the interstellar gamma-ray emission, and then focus on complementary studies of the radio emission from the Galaxy.

We present our software to compute gamma-ray emission from inverse-Compton scattering by cosmic-ray leptons in the heliosphere, as well as in the photospheres of stars. It includes a formulation of modulation in the heliosphere, but it can be used for any user-defined modulation model. Profiles and spectra are output to FITS files in a variety of forms for convenient use. Also included are general-purpose inverse-Compton routines with other features like energy loss rates and emissivity for any user-defined target photon and lepton spectra. The software is publicly available and it is under continuing development.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

Recent observations of high ionization rates of molecular hydrogen in diffuse interstellar clouds point to a distinct low-energy cosmic-ray component. Supposing that this component is made of nuclei, two models for the origin of such particles are explored and low-energy cosmic-ray spectra are calculated which, added to the standard cosmic ray spectra, produce the observed ionization rates. The clearest evidence of the presence of such low-energy nuclei between a few MeV per nucleon and several hundred MeV per nucleon in the interstellar medium would be a detection of nuclear \gamma-ray line emission in the range E_ 0.1 - 10 MeV, which is strongly produced in their collisions with the interstellar gas and dust. Using a recent \gamma-ray cross section compilation for nuclear collisions, \gamma-ray line emission spectra are calculated alongside with the high-energy \gamma-ray emission due to {\pi} 0 decay, the latter providing normalization of the absolute fluxes by comparison with Fermi-LAT observations of the diffuse emission above E \gamma = 0.1 GeV. Our predicted fluxes of strong nuclear \gamma-ray lines from the inner Galaxy are well below the detection sensitivies of INTEGRAL, but a detection, especially of the 4.4-MeV line, seems possible with new-generation \gamma-ray telescopes based on available technology. We predict also strong \gamma-ray continuum emission in the 1-8 MeV range, which in a large part of our model space for low-energy cosmic rays exceeds considerably estimated instrument sensitivities of future telescopes.

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 \mu\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

We present the first results from the CHANG-ES survey, a new survey of 35 edge-on galaxies to search for both in-disk as well as extra-planar radio continuum emission. The motivation and science case for the survey are presented in a companion paper (Paper I). In this paper (Paper II), we outline the observations and data reduction steps required for wide-band calibration and mapping of EVLA data, including polarization, based on C-array test observations of NGC 4631. With modest on-source observing times (30 minutes at 1.5 GHz and 75 minutes at 6 GHz for the test data) we have achieved best rms noise levels of 22 and 3.5 $\mu$Jy beam$^{-1}$ at 1.5 GHz and 6 GHz, respectively. New disk-halo features have been detected, among them two at 1.5 GHz that appear as loops in projection. We present the first 1.5 GHz spectral index map of NGC 4631 to be formed from a single wide-band observation in a single array configuration. This map represents tangent slopes to the intensities within the band centered at 1.5 GHz, rather than fits across widely separated frequencies as has been done in the past and is also the highest spatial resolution spectral index map yet presented for this galaxy. The average spectral index in the disk is $\bar\alpha_{1.5 GHz}\,=\,-0.84\,\pm\,0.05$ indicating that the emission is largely non-thermal, but a small global thermal contribution is sufficient to explain a positive curvature term in the spectral index over the band. Two specific star forming regions have spectral indices that are consistent with thermal emission. Polarization results (uncorrected for internal Faraday rotation) are consistent with previous observations and also reveal some new features. On broad scales, we find strong support for the notion that magnetic fields constrain the X-ray emitting hot gas.

We introduce a new survey to map the radio continuum halos of a sample of 35 edge-on spiral galaxies at 1.5 GHz and 6 GHz in all polarization products. The survey is exploiting the new wide bandwidth capabilities of the Karl G. Jansky Very Large Array (i.e. the Expanded Very Large Array, or EVLA) in a variety of array configurations (B, C, and D) in order to compile the most comprehensive data set yet obtained for the study of radio halo properties. This is the first survey of radio halos to include all polarization products. In this first paper, we outline the scientific motivation of the survey, the specific science goals, and the expected improvements in noise levels and spatial coverage from the survey. Our goals include investigating the physical conditions and origin of halos, characterizing cosmic ray transport and wind speed, measuring Faraday rotation and mapping the magnetic field, probing the in-disk and extraplanar far-infrared - radio continuum relation, and reconciling non-thermal radio emission with high-energy gamma-ray models. The sample size allows us to search for correlations between radio halos and other properties, including environment, star formation rate, and the presence of AGNs. In a companion paper (Paper II) we outline the data reduction steps and present the first results of the survey for the galaxy, NGC 4631.

The Galaxy hosts a widespread population of low-energy positrons revealed by successive generations of gamma-ray telescopes through a bright annihilation emission from the bulge region, with a fainter contribution from the inner disk. The exact origin of these particles remains currently unknown. We estimate the contribution to the annihilation signal of positrons generated in the decay of radioactive 26Al, 56Ni and 44Ti. We adapted the GALPROP propagation code to simulate the transport and annihilation of radioactivity positrons in a model of our Galaxy. Using plausible source spatial distributions, we explored several possible propagation scenarios to account for the large uncertainties on the transport of ~1MeV positrons in the interstellar medium. We then compared the predicted intensity distributions to the INTEGRAL/SPI observations. We obtain similar intensity distributions with small bulge-to-disk ratios, even for extreme large-scale transport prescriptions. At least half of the positrons annihilate close to their sources, even when they are allowed to travel far away. In the high-diffusion, ballistic case, up to 40% of them escape the Galaxy. In proportion, this affects bulge positrons more than disk positrons because they are injected further off the plane in a tenuous medium, while disk positrons are mostly injected in the dense molecular ring. The predicted intensity distributions are fully consistent with the observed longitudinally-extended disk-like emission, but the transport scenario cannot be strongly constrained by the current data. Nucleosynthesis positrons alone cannot account for the observed annihilation emission in the frame of our model. An additional component is needed to explain the strong bulge contribution, and the latter is very likely concentrated in the central regions if positrons have initial energies in the 100keV-1MeV range.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

Using infrared, radio, and gamma-ray data,we investigate the propagation characteristics of cosmic-ray (CR) electrons and nuclei in the 30 Doradus (30\,Dor) star-forming region in the Large Magellanic Cloud (LMC) using a phenomenological model based on the radio-far-infrared correlation within galaxies. Employing a correlation analysis, we derive an average propagation length of \sim 100-140 pc for \sim 3 GeV CR electrons resident in 30 Dor from consideration of the radio and infrared data. Assuming that the observed gamma-ray emission towards 30 Dor is associated with the star-forming region, and applying the same methodology to the infrared and gamma-ray data, we estimate a \sim 20 GeV propagation length of 200-320 pc for the CR nuclei. This is approximately twice as large as for\sim 3 GeV CR electrons, corresponding to a spatial diffusion coefficient that is \sim 4 times higher, scaling as (R/GV)^{\delta} with \delta \approx 0.7-0.8 depending on the smearing kernel used in the correlation analysis. This value is in agreement with the results found by extending the correlation analysis to include \sim 70 GeV CR nuclei traced by the 3-10 GeV gamma-ray data (\delta \approx 0.66+/-0.23). Using the mean age of the stellar populations in 30 Dor and the results from our correlation analysis, we estimate a diffusion coefficient D_{R} \approx 0.9-1.0 \times10^27 (R/GV)^{0.7} cm^2 s^-1. We compare the values of the CR electron propagation length and surface brightness for 30 Dor and the LMC as a whole with those of entire disk galaxies. We find that the trend of decreasing average CR propagation distance with increasing disk-averaged star formation activity holds for the LMC, and extends down to single star-forming regions, at least for the case of 30 Dor.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

Context: The relation between Galactic cosmic-ray electrons, magnetic fields and synchrotron radiation. Aims: We exploit synchrotron radiation to constrain the low-energy interstellar electron spectrum, using various radio surveys and connecting with electron data from Fermi-LAT and other experiments. Methods: The GALPROP programme for cosmic-ray propagation, gamma-ray and synchrotron radiation is used. Secondary electrons and positrons are included. Propagation models based on cosmic-ray and gamma-ray data are tested against synchrotron data from 22 MHz to 94 GHz. Results: The synchrotron data confirm the need for a low-energy break in the cosmic-ray electron injection spectrum. The interstellar spectrum below a few GeV has to be lower than standard models predict, and this suggests less solar modulation than usually assumed. Reacceleration models are more difficult to reconcile with the synchrotron constraints. We show that secondary leptons are important for the interpretation of synchrotron emission. We also consider a cosmic-ray propagation origin for the low-energy break. Conclusions: Exploiting the complementary information on cosmic rays and synchrotron gives unique and essential constraints on electrons, and has implications for gamma rays. This connection is especially relevant now in view of the ongoing PLANCK and Fermi missions.

GALPROP is a numerical code for calculating the galactic propagation of relativistic charged particles and the diffuse emissions produced during their propagation. The code incorporates as much realistic astrophysical input as possible together with latest theoretical developments and has become a de facto standard in astrophysics of cosmic rays. We present GALPROP WebRun, a service to the scientific community enabling easy use of the freely available GALPROP code via web browsers. In addition, we introduce the latest GALPROP version 54, available through this service.

Significant advances have been made in the understanding of the diffuse Galactic hard X-ray continuum emission using data from the INTEGRAL observatory. The diffuse hard power-law component seen with the INTEGRAL/SPI spectrometer has been identified with inverse-Compton emission from relativistic (GeV) electrons on the cosmic microwave background and Galactic interstellar radiation field. In the present analysis, SPI data from 2003 to 2009, with a total exposure time of ~ 10^8 s, are used to derive the Galactic ridge hard X-ray spatial distribution and spectrum between 20 keV and 2.4 MeV. Both are consistent with predictions from the GALPROP code. The good agreement between measured and predicted emission from keV to GeV energies suggests that the correct production mechanisms have been identified. We discuss the potential of the SPI data to provide an indirect probe of the interstellar cosmic-ray electron distribution, in particular for energies below a few GeV.

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $\gamma$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

In over two years of operation Fermi-LAT has revolutionized our knowledge of the gamma-ray sky. Interstellar gamma rays are part of this new era and allow unprecedented tests for models of cosmic rays in the Galaxy. The extension to lower energies with INTEGRAL/SPI data is also evolving. The global multiwavelength luminosity of the Milky Way has been derived, with implications for the Galactic energy balance and the radio-FIR correlation.

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.

The first gamma-ray line originating from outside the solar system that was ever detected is the 511 keV emission from positron annihilation in the Galaxy. Despite 30 years of intense theoretical and observational investigation, the main sources of positrons have not been identified up to now. Observations in the 1990's with OSSE/CGRO showed that the emission is strongly concentrated towards the Galactic bulge. In the 2000's, the SPI instrument aboard ESA's INTEGRAL gamma-ray observatory allowed scientists to measure that emission across the entire Galaxy, revealing that the bulge/disk luminosity ratio is larger than observed in any other wavelength. This mapping prompted a number of novel explanations, including rather "exotic ones (e.g. dark matter annihilation). However, conventional astrophysical sources, like type Ia supernovae, microquasars or X-ray binaries, are still plausible candidates for a large fraction of the observed total 511 keV emission of the bulge. A closer study of the subject reveals new layers of complexity, since positrons may propagate far away from their production sites, making it difficult to infer the underlying source distribution from the observed map of 511 keV emission. However, contrary to the rather well understood propagation of high energy (>GeV) particles of Galactic cosmic rays, understanding the propagation of low energy (~MeV) positrons in the turbulent, magnetized interstellar medium, still remains a formidable challenge. We review the spectral and imaging properties of the observed 511 keV emission and we critically discuss candidate positron sources and models of positron propagation in the Galaxy.

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.

We aim to exploit the available INTEGRAL/SPI data to provide time-averaged spectra of the brightest soft gamma-ray sources. Employing a maximum-likelihood fit technique for our SPI data analysis, we take as input to our source model the source catalog derived by Bouchet et al. (2008) from a SPI all-sky study. We use the first four years of public SPI data and extract spectra between 25 keV and 1 MeV for the 20 catalog sources detected by Bouchet et al. at 200 - 600 keV with >= 2.5 sigma. In order to verify our analysis, we also extract spectra for the same sources from the corresponding INTEGRAL/ISGRI data. We fit adequate spectral models to the energy range 25-1000 keV for SPI and 25-600 keV for ISGRI. We use our spectra from the Crab (which is among the 20 sources studied here) and an empty location in a crowded field to derive an estimation of the systematic errors. The agreement between our SPI and ISGRI measurements is good if we normalise them on the Crab spectrum. Our SPI flux measurements also agree well with those by Bouchet et al. (2008). All 20 sources in our sample are detected independently in the bands 25-100 keV and 100-200 keV. At 200-600 keV we detect eight sources, at 600-1000 keV we detect two sources. Our spectra agree well with the results from previous publications where available. For six of the 14 XRBs in our sample we find evidence for a hard powerlaw-component which becomes dominant above the cutoff energy of the thermal Comptonization component. In two of these cases, our study provides the first indication of such emission. For the others, our results confirm previous studies. Our Crab spectrum (from 1.3 Ms exposure), shows a significant flux in all points and is well described by a powerlaw with a break near 100 keV and spectral indices 2.11 and 2.20.

The Fermi-LAT experiment recently reported high precision measurements of the spectrum of cosmic-ray electrons-plus-positrons (CRE) between 20 GeV and 1 TeV. The spectrum shows no prominent spectral features, and is significantly harder than that inferred from several previous experiments. Here we discuss several interpretations of the Fermi results based either on a single large scale Galactic CRE component or by invoking additional electron-positron primary sources, e.g. nearby pulsars or particle Dark Matter annihilation. We show that while the reported Fermi-LAT data alone can be interpreted in terms of a single component scenario, when combined with other complementary experimental results, specifically the CRE spectrum measured by H.E.S.S. and especially the positron fraction reported by PAMELA between 1 and 100 GeV, that class of models fails to provide a consistent interpretation. Rather, we find that several combinations of parameters, involving both the pulsar and dark matter scenarios, allow a consistent description of those results. We also briefly discuss the possibility of discriminating between the pulsar and dark matter interpretations by looking for a possible anisotropy in the CRE flux.

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.

The Galactic cosmic-ray propagation code GALPROP is designed to make predictions of many kinds of data self-consistently, including direct cosmic-ray measurements, gamma rays and synchrotron radiation. In the decade since its conception it has undergone considerable development, which is continuing. A public version was made available a few years ago, supported by a website. We describe the new features of the current version which will become part of the next public release. Plans for future developments are also mentioned.

Despite a century of cosmic-ray measurements there seems to have been no attempt to collect these data systematically for use by the CR community. The result is that everyone makes their own collection as required, a large duplication of effort. Therefore we have started a project to place published Galactic CR measurements in a database available online. It currently addresses energies up to 100 TeV, including elemental spectra, secondary/primary nuclei ratios, antiprotons, electrons and positrons. It is updated regularly as data appears in the literature. It is supported by access software. The community is encouraged to participate by providing data, pointing out errors or omissions, and suggestions for improvements.

The LAT instrument on the Fermi Gamma-Ray Space Telescope is performing an all-sky survey from 20 MeV to 300 GeV with unprecedented statistics and angular resolution. This is providing a wealth of new information on the non-thermal emission from the Galactic interstellar medium with implications for cosmic rays and Galactic structure. First results at intermediate latitudes have already shown good agreement with predictions based on direct measurements of cosmic rays, suggesting that at least the local (within about 1 kpc from the Sun) gamma-ray emission is understood. We will present the first spectra from regions over the sky using the LAT data, and profiles for selected energies. The aim here is to evaluate the agreement with the models and assess what we can expect to learn as this analysis matures.

Gamma-ray line emission from the radioactive decay of 26Al reflects nucleosynthesis in massive stars and supernovae. We use INTEGRAL 26Al measurements to characterize the distribution and characteristics of 26Al source regions throughout the Galaxy. We detect the 26Al line from the inner Galaxy at 28\sigma significance. The line appears narrow, and we constrain broadening in the source regions to <1.3 keV (2\sigma). Different sky distribution models do not significantly affect those large-scale results. The 26Al intensity for the inner Galaxy is derived as (2.9\pm 0.2)\times 10^{-4}ph\ cm^{-2} s^{-1} rad^{-1}, consistent with earlier results from COMPTEL and SPI data. This can be translated to an 26Al mass of 2.7\pm 0.7 solar mass in the Galaxy as a whole. The 26Al intensity is also confirmed to be somewhat brighter in the 4th than in the 1st quadrant (ratio 1.3\pm 0.2). 26Al spectra separately derived for regions along the Galactic plane show clear line centroid shifts, attributed largely to the Galaxy's large-scale rotation. The 26Al line toward the direction of the Aquila region (20deg < l < 40deg) appears somewhat broadened. Latitudinal variations of 26Al emission towards the inner Galaxy are studied, finding a latitudinal scale height of 130^{+120}_{-70} pc (1\sigma) for 26Al in the inner Galaxy and a hint of peculiar 26Al emission towards the region l<0deg, b>5deg.

This is a report on the findings of the SNR/cosmic-ray working group for the white paper on the status and future of ground-based gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe shell-type supernova remnants and diffuse emission from cosmic rays at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to study the acceleration of relativistic charged particles which is one of the main unsolved, yet fundamental, problems in modern astrophysics. The acceleration of particles relies on interactions between energetic particles and magnetic turbulence. In the case of SNRs we can perform spatially resolved studies in systems with known geometry, and the plasma physics deduced from these observations will help us to understand other systems where rapid particle acceleration is believed to occur and where observations as detailed as those of SNRs are not possible.

We use Swift/BAT Earth occultation data at different geomagnetic latitudes to derive a sensitive measurement of the Cosmic X-ray background (CXB) and of the Earth albedo emission in the 15--200 keV band. We compare our CXB spectrum with recent (INTEGRAL, BeppoSAX) and past results (HEAO-1) and find good agreement. Using an independent measurement of the CXB spectrum we are able to confirm our results. This study shows that the BAT CXB spectrum has a normalization ~8(+/-3)% larger than the HEAO-1 measurement. The BAT accurate Earth albedo spectrum can be used to predict the level of photon background for satellites in low Earth and mid inclination orbits.

Diffuse $\gamma$-ray emission produced by the interaction of cosmic-ray particles with matter and radiation in the Galaxy can be used to probe the distribution of cosmic rays and their sources in different regions of the Galaxy. With its large field of view and long observation time, the Milagro Gamma Ray Observatory is an ideal instrument for surveying large regions of the Northern Hemisphere sky and for detecting diffuse $\gamma$-ray emission at very high energies. Here, the spatial distribution and the flux of the diffuse $\gamma$-ray emission in the TeV energy range with a median energy of 15 TeV for Galactic longitudes between 30$^\circ$ and 110$^\circ$ and between 136$^\circ$ and 216$^\circ$ and for Galactic latitudes between -10$^\circ$ and 10$^\circ$ are determined. The measured fluxes are consistent with predictions of the GALPROP model everywhere except for the Cygnus region ($l\in[65^\circ,85^\circ]$). For the Cygnus region, the flux is twice the predicted value. This excess can be explained by the presence of active cosmic ray sources accelerating hadrons which interact with the local dense interstellar medium and produce gamma rays through pion decay.

We investigate the sensitivity of the Gamma-ray Large Area Space Telescope (GLAST) to indirectly detect weakly interacting massive particles (WIMPs) through the $\gamma$-ray signal that their pair annihilation produces. WIMPs are among the favorite candidates to explain the compelling evidence that about 80% of the mass in the Universe is non-baryonic dark matter (DM). They are serendipitously motivated by various extensions of the standard model of particle physics such as Supersymmetry and Universal Extra Dimensions (UED). With its unprecedented sensitivity and its very large energy range (20 MeV to more than 300 GeV) the main instrument on board the GLAST satellite, the Large Area Telescope (LAT), will open a new window of discovery. As our estimates show, the LAT will be able to detect an indirect DM signature for a large class of WIMP models given a cuspy profile for the DM distribution. Using the current state of the art Monte Carlo and event reconstruction software developed within the LAT collaboration, we present preliminary sensitivity studies for several possible sources inside and outside the Galaxy. We also discuss the potential of the LAT to detect UED via the electron/positron channel. Diffuse background modeling and other background issues that will be important in setting limits or seeing a signal are presented.

A recent re-determination of the non-thermal component of the hard X-ray to soft gamma-ray emission from the Galactic ridge, using the SPI instrument on the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) Observatory, is shown to be well reproduced as inverse-Compton emission from the interstellar medium. Both cosmic-ray primary electrons and secondary electrons and positrons contribute to the emission. The prediction uses the GALPROP model and includes a new calculation of the interstellar radiation field. This may solve a long-standing mystery of the origin of this emission, and potentially opens a new window on Galactic cosmic rays.

Context: The Sun has recently been predicted to be an extended source of gamma-ray emission, produced by inverse-Compton (IC) scattering of cosmic-ray (CR) electrons on the solar radiation field. The emission was predicted to be extended and a confusing foreground for the diffuse extragalactic background even at large angular distances from the Sun. The solar disk is also expected to be a steady gamma-ray source. While these emissions are expected to be readily detectable in the future by GLAST, the situation for available EGRET data is more challenging. Aims: The theory of gamma-ray emission from IC scattering on the solar radiation field by Galactic CR electrons is given in detail. This is used as the basis for detection and model verification using EGRET data. Methods: We present a detailed study of the solar emission using the EGRET database, accounting for the effect of the emission from 3C 279, the moon, and other sources, which interfere with the solar emission. The analysis was performed for 2 energy ranges, above 300 MeV and for 100-300 MeV, as well as for the combination to improve the detection statistics. The technique was tested on the moon signal, with our results consistent with previous work. Results: Analyzing the EGRET database, we find evidence of emission from the solar disk and its halo. The observations are compared with our model for the extended emission. The spectrum of the solar disk emission and the spectrum of the extended emission have been obtained. The spectrum of the moon is also given. Conclusions: The observed intensity distribution and the flux are consistent with the predicted model of IC gamma-rays from the halo around the Sun.

We have processed the data accumulated with INTEGRAL SPI during 4 years (~ 51 Ms) to study the Galactic ``diffuse'' emission morphology in the 20 keV to 8 MeV energy range. To achieve this objective, we have derived simultaneously an all-sky census of emitting sources and images of the Galactic Ridge (GR) emission. In the central radian, the resolved point source emission amounts to 88%, 91% and 68% of the total emission in the 25-50, 50-100 and 100-300 keV domains respectively. We have compared the GR emission spatial distribution to those obtained from CO and NIR maps, and quantified our results through latitude and longitude profiles. Below 50 keV, the SPI data are better traced by the latter, supporting a stellar origin for this emission. Furthermore, we found that the GR emission spectrum follows a power law with a photon index ~ 1.55 above 50 keV while an additional component is required below this energy. This component shows a cutoff around 30 keV, reinforcing a stellar origin, as proposed by Krivonos et al. (2007). The annihilation diffuse emission component is extracted simultaneously, leading to the determination of the related parameters (positronium flux and fraction). A specific discussion is devoted to the annihilation line distribution since a significant emission is detected over a region as large as ~80 degrees by ~10 degrees potentially associated with the disk or halo surrounding the central regions of our Galaxy.

We applied the Maximum Likelihood method, as an image reconstruction algorithm, to the BAT X-ray Survey (BXS). This method was specifically designed to preserve the full statistical information in the data and to avoid mosaicking of many exposures with different pointing directions, thus reducing systematic errors when co-adding images. We reconstructed, in the 14-170 keV energy band, the image of a 90x90 deg$^2$ sky region, centered on (RA,DEC)=105$^{\circ}$,-25$^{\circ}$, which BAT surveyed with an exposure time of $\sim1$ Ms (in Nov. 2005). The best sensitivity in our image is $\sim0.85$ mCrab or $2.0\times 10^{-11}$ erg cm$^{-2}$. We detect 49 hard X-ray sources above the 4.5 $\sigma$ level; of these, only 12 were previously known as hard X-ray sources ($>$15 keV). Swift/XRT observations allowed us to firmly identify the counterparts for 15 objects, while 2 objects have Einstein IPC counterparts \citep{harris90}; in addition to those, we found a likely counterpart for 13 objects by correlating our sample with the ROSAT All-Sky Survey Bright Source Catalog \citep{voges99}. 7 objects remain unidentified. Analysis of the noise properties of our image shows that $\sim75$% of the area is surveyed to a flux limit of $\sim$1 mCrab. This study shows that the coupling of the Maximum Likelihood method to the most sensitive, all-sky surveying, hard X-ray instrument, BAT, is able to probe for the first time the hard X-ray sky to the mCrab flux level. The successful application of this method to BAT demonstrates that it could also be applied with advantage to similar instruments like INTEGRAL-IBIS.

The spectrometer SPI on INTEGRAL allows for the first time simultaneous imaging of diffuse and point-like emission in the hard X-ray and soft gamma-ray regime. To fully exploit the capabilities of the instrument, we implemented the MREM image deconvolution algorithm, initially developed for COMPTEL data analysis, to SPI data analysis. We present the performances of the algorithm by means of simulations and apply it to data accumulated during the first 2 mission years of INTEGRAL. Skymaps are presented for the 1809 keV gamma-ray line, attributed to the radioactive decay of 26Al, and for continuum energy bands, covering the range 20 keV - 3 MeV. The 1809 keV map indicates that emission is clearly detected by SPI from the inner Galactic radian and from the Cygnus region. The continuum maps reveal the transition between a point-source dominated hard X-ray sky to a diffuse emission dominated soft gamma-ray sky. From the skymaps, we extract a Galactic ridge emission spectrum that matches well SPI results obtained by model fitting. By comparing our spectrum with the cumulative flux measured by IBIS from point sources, we find indications for the existence of an unresolved or diffuse emission component above ~100 keV.

Inverse Compton scattering by relativistic electrons produces a major component of the diffuse emission from the Galaxy. The photon fields involved are the cosmic microwave background and the interstellar radiation field from stars and dust. Calculations of the inverse Compton distribution have usually assumed a smooth ISRF, but in fact a large part of the Galactic luminosity comes from the most luminous stars which are rare. Therefore we expect the ISRF, and hence the inverse Compton emission, to be clumpy. We also show that some of the most luminous stars may be visible to GLAST. In this paper we give an update on our previous work including examples of the intensity distribution around stars, and the predicted spectrum of Cygnus OB2.

In this concluding part of the series of three papers dedicated to the Swift/BAT hard X-ray survey (BXS), we focus on the X-ray spectral analysis and statistical properties of the source sample. Using a dedicated method to extract time-averaged spectra of BAT sources we show that Galactic sources have, generally, softer spectra than extragalactic objects and that Seyfert 2 galaxies are harder than Seyfert 1s. The averaged spectrum of all Seyfert galaxies is consistent with a power-law with photon index of 2.00 (+/-0.07). The cumulative flux-number relation for the extragalactic sources in the 14-170 keV band is best described by a power-law with a slope alpha=1.55 (+/-0.20) and a normalization of 9.6$\pm1.9 \times 10^{-3}$ AGN deg$^{-2}$ (or 396(+/-80) AGN all-sky) above a flux level of 2$\times 10^{-11}$erg cm$^{-2}$ s$^{-1}$ (~0.85 mCrab). The integration of the cumulative flux per unit area indicates that BAT resolves 1-2% of the X-ray background emission in the 14-170 keV band. A sub-sample of 24 extragalactic sources above the 4.5 sigma detection limit is used to study the statistical properties of AGN. This sample comprises local Seyfert galaxies (z=0.026, median value) and ~10% blazars. We find that 55% of the Seyfert galaxies are absorbed by column densities of Log(N_H)>22, but that none is a bona fide Compton-thick. This study shows the capabilities of BAT to probe the hard X-ray sky to the mCrab level.

Diffuse emission from the Milky Way dominates the gamma-ray sky. About 80% of the high-energy luminosity of the Milky Way comes from processes in the interstellar medium. The Galactic diffuse emission traces interactions of energetic particles, primarily protons and electrons, with the interstellar gas and radiation field, thus delivering information about cosmic-ray spectra and interstellar mass in distant locations. Additionally, the Galactic diffuse emission is the celestial foreground for the study of gamma-ray point sources and the extragalactic diffuse gamma-ray emission. We report on the latest developments in the modelling of the Galactic diffuse emission, which will be used for the Gamma Ray Large Area Space Telescope (GLAST) investigations.

We use the GALPROP code and the Advanced Composition Explorer (ACE) data to derive the cosmic ray (CR) isotopic composition at the sources. The composition is derived for two propagation models, diffusive reacceleration and plain diffusion. We show that the compositions derived assuming these two propagation models are different. We also compare the isotopic composition at the sources with the latest solar composition.

The origin of the soft gamma-ray (200 keV - 1 MeV) galactic ridge emission is one of the long-standing mysteries in the field of high-energy astrophysics. Population studies at lower energies have shown that emission from accreting compact objects gradually recedes in this domain, leaving place to another source of gamma-ray emission that is characterised by a hard power-law spectrum extending from 100 keV up to 100 MeV The nature of this hard component has remained so far elusive, partly due to the lack of sufficiently sensitive imaging telescopes that would be able to unveil the spatial distribution of the emission. The SPI telescope aboard INTEGRAL allows now for the first time the simultaneous imaging of diffuse and point-like emission in the soft gamma-ray regime. We present here all-sky images of the soft gamma-ray continuum emission that clearly reveal the morphology of the different emission components. We discuss the implications of our results on the nature of underlying emission processes and we put our results in perspective of GLAST studies of diffuse galactic continuum emission.

Gamma-ray line emission from radioactive decay of 60Fe provides constraints on nucleosynthesis in massive stars and supernovae. The spectrometer SPI on board INTEGRAL has accumulated nearly three years of data on gamma-ray emission from the Galactic plane. We have analyzed these data with suitable instrumental-background models and sky distributions to produce high-resolution spectra of Galactic emission. We detect the gamma-ray lines from 60Fe decay at 1173 and 1333 keV, obtaining an improvement over our earlier measurement of both lines with now 4.9 sigma significance for the combination of the two lines. The average flux per line is (4.4 \pm 0.9) \times 10^{-5} ph cm^{-2} s^{-1} rad^{-1} for the inner Galaxy region. Deriving the Galactic 26Al gamma-ray line flux with using the same set of observations and analysis method, we determine the flux ratio of 60Fe/26Al gamma-rays as 0.148 \pm 0.06. The current theoretical predictions are still consistent with our result.

The large-scale Galactic interstellar radiation field (ISRF) is the result of stellar emission and dust re-processing of starlight. Where the energy density of the ISRF is high (e.g., the Galactic Centre), the dominant gamma-ray emission in individual supernova remnants (SNRs), such as G0.9+0.1, may come from inverse Compton (IC) scattering of the ISRF. Several models of the ISRF exist. The most recent one, which has been calculated by us, predicts a significantly higher ISRF than the well-used model of Mathis, Mezger, and Panagia. However,comparison with data is limited to local observations. Based on our current estimate of the ISRF we predict the gamma-ray emission in the SNRs G0.9+0.1 and RXJ1713, and pair-production absorption features above 20 TeV in the spectra of G0.9+0.1, J1713-381, and J1634-472. We discuss how GLAST, along with current and future very high energy instruments, may be able to provide upper bounds on the large-scale ISRF.

Diffuse emission is produced in energetic cosmic ray (CR) interactions, mainly protons and electrons, with the interstellar gas and radiation field and contains the information about particle spectra in distant regions of the Galaxy. It may also contain information about exotic processes such as dark matter annihilation, black hole evaporation etc. A model of the diffuse emission is important for determination of the source positions and spectra. Calculation of the Galactic diffuse continuum gamma-ray emission requires a model for CR propagation as the first step. Such a model is based on theory of particle transport in the interstellar medium as well as on many kinds of data provided by different experiments in Astrophysics and Particle and Nuclear Physics. Such data include: secondary particle and isotopic production cross sections, total interaction nuclear cross sections and lifetimes of radioactive species, gas mass calibrations and gas distribution in the Galaxy (H_2, H I, H II), interstellar radiation field, CR source distribution and particle spectra at the sources, magnetic field, energy losses, gamma-ray and synchrotron production mechanisms, and many other issues. We are continuously improving the GALPROP model and the code to keep up with a flow of new data. Improvement in any field may affect the Galactic diffuse continuum gamma-ray emission model used as a background model by the GLAST LAT instrument. Here we report about the latest improvements of the GALPROP and the diffuse emission model.

The Sun was recently predicted to be an extended source of gamma-ray emission, produced by inverse-Compton scattering of cosmic-ray electrons with the solar radiation. The emission was predicted to contribute to the diffuse extragalactic background even at large angular distances from the Sun. While this emission is expected to be readily detectable in future by GLAST, the situation for available EGRET data is more challenging. We present a detailed study of the EGRET database, using a time dependent analysis, accounting for the effect of the emission from 3C 279, the moon, and other sources, which interfere with the solar signal. The technique has been tested on the moon signal, with results consistent with previous work. We find clear evidence for emission from the Sun and its vicinity. The observations are compared with our model for the extended emission.

The imaging spectrometer SPI on board ESA's INTEGRAL observatory provides us with an unprecedented view of positron annihilation in our Galaxy. The first sky maps in the 511 keV annihilation line and in the positronium continuum from SPI showed a puzzling concentration of annihilation radiation in the Galactic bulge region. By now, more than twice as many INTEGRAL observations are available, offering new clues to the origin of Galactic positrons. We present the current status of our analyses of this augmented data set. We now detect significant emission from outside the Galactic bulge region. The 511 keV line is clearly detected from the Galactic disk; in addition, there is a tantalizing hint at possible halo-like emission. The available data do not yet permit to discern whether the emission around the bulge region originates from a halo-like component or from a disk component that is very extended in latitude.

We survey the theory and experimental tests for the propagation of cosmic rays in the Galaxy up to energies of 10^15 eV. A guide to the previous reviews and essential literature is given, followed by an exposition of basic principles. The basic ideas of cosmic-ray propagation are described, and the physical origin of its processes are explained. The various techniques for computing the observational consequences of the theory are described and contrasted. These include analytical and numerical techniques. We present the comparison of models with data including direct and indirect - especially gamma-ray - observations, and indicate what we can learn about cosmic-ray propagation. Some particular important topics including electrons and antiparticles are chosen for discussion.

The diffuse gamma radiation arising from the interaction of cosmic ray particles with matter and radiation in the Galaxy is one of the few probes available to study the origin of the cosmic rays. Milagro is a water Cherenkov detector that continuously views the entire overhead sky. The large field-of-view combined with the long observation time makes Milagro the most sensitive instrument available for the study of large, low surface brightness sources such as the diffuse gamma radiation arising from interactions of cosmic radiation with interstellar matter. In this paper we present spatial and flux measurements of TeV gamma-ray emission from the Cygnus Region. The TeV image shows at least one new source MGRO J2019+37 as well as correlations with the matter density in the region as would be expected from cosmic-ray proton interactions. However, the TeV gamma-ray flux as measured at ~12 TeV from the Cygnus region (after excluding MGRO J2019+37) exceeds that predicted from a conventional model of cosmic ray production and propagation. This observation indicates the existence of either hard-spectrum cosmic-ray sources and/or other sources of TeV gamma rays in the region.

The abundance of Li, Be, and B isotopes in galactic cosmic rays (GCR) between E=50-200 MeV/nucleon has been observed by the Cosmic Ray Isotope Spectrometer (CRIS) on NASA's ACE mission since 1997 with high statistical accuracy. Precise observations of Li, Be, B can be used to constrain GCR propagation models. \iffalse Precise observations of Li, Be, and B in addition to well-measured production cross-sections are used to further constrain GCR propagation models. \fi We find that a diffusive reacceleration model with parameters that best match CRIS results (e.g. B/C, Li/C, etc) are also consistent with other GCR observations. A $\sim$15--20% overproduction of Li and Be in the model predictions is attributed to uncertainties in the production cross-section data. The latter becomes a significant limitation to the study of rare GCR species that are generated predominantly via spallation.

Population synthesis is used to study the contribution from undetected sources to the Galactic ridge emission measured by EGRET. Synthesized source counts are compared with the 3rd EGRET catalogue at low and high latitudes. For pulsar-like populations, 5-10% of the emission >100 MeV comes from sources below the EGRET threshold. A steeper luminosity function can increase this to 20% without violating EGRET source statistics. Less luminous populations can produce much higher values without being detected. Since the unresolved source spectrum is different from the interstellar spectrum, it could provide an explanation of the observed MeV and GeV excesses above the predictions, and we give an explicit example of how this could work.

Inverse Compton scattering by relativistic electrons produces a major component of the diffuse emission from the Galaxy. The photon fields involved are the cosmic microwave background and the interstellar radiation field from stars and dust. Calculations of the inverse Compton distribution have usually assumed a smooth ISRF, but in fact a large part of the Galactic luminosity comes from the most luminous stars which are rare. Therefore we expect the ISRF, and hence the inverse Compton emission, to be clumpy at some level which could be detectable by instruments such as GLAST. Even individual nearby luminous stars could be detectable assuming just the normal cosmic-ray electron spectrum. We present the basic formalism required and give possible candidate stars to be detected, and make predictions for GLAST. Then we apply the formalism to OB associations and the Sun, showing that the IC emission produced is not negligible compared to the sensitivity of current or coming detectors. We estimate that the gamma-ray flux from the halo around the Sun contributes to the diffuse background emission at about the few percent level.

We discuss uncertainties and possible sources of errors associated with the determination of the diffuse Galactic gamma-ray emission using the EGRET data. Most of the issues will be relevant also in the GLAST era. The focus here is on issues that impact evaluation of dark matter annihilation signals against the diffuse gamma-ray emission of the Milky Way.

The evidence for particle acceleration in supernova shells comes from electrons whose synchrotron emission is observed in radio and X-rays. Recent observations by the HESS instrument reveal that supernova remnants also emit TeV gamma-rays; long awaited experimental evidence that supernova remnants can accelerate cosmic rays up to the ``knee'' energies. Still, uncertainty exists whether these gamma-rays are produced by electrons via inverse Compton scattering or by protons via neutral pion decay. The multi-wavelength spectra of supernova remnants can be fitted with both mechanisms, although a preference is often given to neutral pion decay due to the spectral shape at very high energies. A recent study of the interstellar radiation field indicates that its energy density, especially in the inner Galaxy, is higher than previously thought. In this paper we evaluate the effect of the interstellar radiation field on the inverse Compton emission of electrons accelerated in a supernova remnant located at different distances from the Galactic Centre. We show that contribution of optical and infra-red photons to the inverse Compton emission may exceed the contribution of cosmic microwave background and in some cases broaden the resulted gamma-ray spectrum. Additionally, we show that if a supernova remnant is located close to the Galactic Centre its gamma-ray spectrum will exhibit a ``universal'' cutoff at very high energies due to the Klein-Nishina effect and not due to the cut-off of the electron spectrum. As an example, we apply our calculations to the supernova remnants RX J1713.7-3946 and G0.9+0.1 recently observed by HESS.

The attenuation of very high energy gamma rays by pair production on the Galactic interstellar radiation field has long been thought of as negligible. However, a new calculation of the interstellar radiation field consistent with multi-wavelength observations by DIRBE and FIRAS indicates that the energy density of the Galactic interstellar radiation field is higher, particularly in the Galactic center, than previously thought. We have made a calculation of the attenuation of very high energy gamma rays in the Galaxy using this new interstellar radiation field which takes into account its nonuniform spatial and angular distributions. We find that the maximum attenuation occurs around 100 TeV at the level of about 25% for sources located at the Galactic center, and is important for both Galactic and extragalactic sources.

The physical processes involved in diffusion of Galactic cosmic rays in the interstellar medium are addressed. We study the possibility that the nonlinear MHD cascade sets the power-law spectrum of turbulence which scatters charged energetic particles. We find that the dissipation of waves due to the resonant interaction with cosmic ray particles may terminate the Kraichnan-type cascade below wavelengths 10^13 cm. The effect of this wave dissipation has been incorporated in the GALPROP numerical propagation code in order to asses the impact on measurable astrophysical data. The energy-dependence of the cosmic-ray diffusion coefficient found in the resulting self-consistent model may explain the peaks in the secondary to primary nuclei ratios observed at about 1 GeV/nucleon.

Gamma-rays from radioactive 26Al (half life ~7.2 10^5 yr) provide a 'snapshot' view of ongoing nucleosynthesis in the Galaxy. The Galaxy is relatively transparent to such gamma-rays, and emission has been found concentrated along the plane of the Galaxy. This led to the conclusion1 that massive stars throughout the Galaxy dominate the production of 26Al. On the other hand, meteoritic data show locally-produced 26Al, perhaps from spallation reactions in the protosolar disk. Furthermore, prominent gamma-ray emission from the Cygnus region suggests that a substantial fraction of Galactic 26Al could originate in localized star-forming regions. Here we report high spectral resolution measurements of 26Al emission at 1808.65 keV, which demonstrate that the 26Al source regions corotate with the Galaxy, supporting its Galaxy-wide origin. We determine a present-day equilibrium mass of 2.8 (+/-0.8) M_sol of 26Al. We use this to estimate that the frequency of core collapse (i.e. type Ib/c and type II) supernovae to be 1.9(+/- 1.1) events per century.

We performed a spectroscopic study of the 1809 keV gamma-ray line from 26Al decay in the Galaxy using the SPI imaging spectrometer with its high-resolution Ge detector camera on the INTEGRAL observatory. We analyzed observations of the first two mission years, fitting spectra from all 7130 telescope pointings in narrow energy bins to models of instrumental background and the 26Al sky. Instrumental background is estimated from independent tracers of cosmic-ray activation. The shape of the 26Al signal is compared to the instrumental response to extract the width of the celestial line. We detect the 26Al line at \~16sigma significance. The line is broadened only slightly, if at all; we constrain the width to be below 2.8 keV (FWHM, 2 sigma). The average Doppler velocities of 26Al at the time of its decay in the interstellar medium (decay time~1.04 My) therefore are probably around 100 km/s, in agreement with expectations from Galactic rotation and interstellar turbulence. The flux and spatial distribution of the emission are found consistent with previous observations. The derived amount of 26Al in the Galaxy is 2.8 (+/-0.8) M_solar.

The INTEGRAL observatory has been performing a deep survey of the Galactic central radian since 2003, with the goal of both extracting a catalog of sources and gaining insight into the Galactic diffuse emission. This paper concentrates on the estimation of the total point sources emission contribution. It is now clear that unresolved point sources contribute to the observed diffuse emission; the increasing sensitivity of instruments with time has lead to a steady decrease in estimates of this ``diffuse emission''. We have analysed the first year data obtained with the spectrometer and imager SPI on board INTEGRAL. First, a catalog of 63 hard X-ray sources detected, time-averaged, during our 2003 Galactic plane survey, is derived. Second, after extracting the spectra of the sources detected by SPI, their combined contribution is compared to the total (resolved and unresolved) emission from the Galactic ridge. The data analysis is complex: it requires us to split the total emission into several components, as discrete sources and diffuse emission are superimposed in SPI data. The main result is that point source emission dominates in the hard X-ray/soft $\gamma$-ray domain, and contributes around 90 % of the total emission around 100 keV, while above 250 keV, diffuse electron-positron annihilation, through its three-photon positronium continuum with a positronium fraction $\sim$ 0.97 and the 511 keV electron-positron line, dominates over the sources.

The Galactic diffuse emission is potentially able to reveal much about the sources and propagation of cosmic rays (CR), their spectra and intensities in distant locations. It can possibly unveil WIMP dark matter (DM) through its annihilation signatures. The extragalactic background may provide vital information about the early stages of the universe, neutralino annihilation, and unresolved sources (blazars?) and their cosmological evolution. The gamma-ray instrument EGRET on the CGRO contributed much to the exploration of the Galactic diffuse emission. The new NASA Gamma-ray Large Area Space Telescope (GLAST) is scheduled for launch in 2007; study of the diffuse gamma-ray emission is one of the priority goals. We describe current understanding of the diffuse emission and its potential for future discoveries.

The diffuse continuum emission from the Galactic plane in the energy range 18-1000 keV has been studied using 16 Ms of data from the SPI instrument on INTEGRAL. With such an exposure we can exploit the imaging properties of SPI to achieve a good separation of point sources from the various diffuse components. Using a candidate-source catalogue derived with IBIS on INTEGRAL and a number of sky distribution models we obtained spectra resolved in Galactic longitude. We can identify spectral components of a diffuse continuum of power law shape with index about 1.7, a positron annihilation component with a continuum from positronium and the line at 511 keV, and a second, roughly power-law component from detected point sources. Our analysis confirms the concentration of positron annihilation emission in the inner region (|l|<10), the disk (10<|l|<30) being at least a factor 7 weaker in this emission. The power-law component in contrast drops by only a factor 2, showing a quite different longitude distribution and spatial origin. Detectable sources constitute about 90% of the total Galactic emission between 20 and 60 keV, but have a steeper spectrum than the diffuse emission, their contribution to the total emission dropping rapidly to a small fraction at higher energies. In the SPI energy range the flux is lower than found by OSSE, probably due to the more complete accounting for sources by SPI. The power-law emission is difficult to explain as of interstellar origin, inverse Compton giving at most 10%, and instead a population of unresolved point sources is proposed as a possible origin, AXPs with their spectra hardening above 100 keV being plausible candidates. We present a broadband spectrum of the emission.

Cosmic-ray electrons and positrons propagating in the Galaxy produce diffuse gamma-rays via the inverse Compton (IC) process. The low energy target photon populations with which the cosmic-rays interact during propagation are produced by stars, this stellar light being reprocessed by Galactic dust. Detailed modelling of the Galactic stellar distribution, dust distribution, and treatment of the absorption and scattering of light is therefore required to obtain accurate models for the low energy Galactic photon distribution and spectrum. Using a realistic Galactic stellar distribution model, and dust distribution, we calculate the diffuse radiation field from stars in the Galaxy (the `optical' radiation field), including absorption and scattering. Using a dust heating code, we self-consistently calculate the infra-red radiation field for the same dust model used for the optical calculation; both transient and equilibrium heating are included. We present the calculated radiation field spectra and distributions, and will use these to calculate the expected Galactic diffuse IC gamma-ray spectrum.

The GALPROP model for cosmic-ray propagation is able to make explicit predictions for the distribution of galactic diffuse gamma-rays. We compare different propagation models with gamma-ray spectra measured by EGRET for various regions of the sky. This allows sensitive tests of alternative explanations for the apparent excess emission observed at GeV gamma-rays. We find that a population of hard-spectrum gamma-ray sources cannot be solely responsible for the excess since it also appears at high latitudes; on the other hand a hard cosmic-ray electron spectrum cannot explain the gamma-ray excess in the inner Galaxy. By normalizing the cosmic ray spectra within reasonable bounds under preservation of their shape we are able to obtain our best prediction of the Galactic component of diffuse gamma rays, and show that away from the Galactic plane it gives an accurate prediction of the observed gamma-ray intensities. On this basis we reevaluate the extragalactic gamma-ray background. We find that for some energies previous work underestimated the Galactic contribution and hence overestimated the background. The new EGRB spectrum shows a positive curvature similar to that expected for models of the extragalactic gamma-ray emission based on contributions from unresolved blazars.

The nuclei fraction in cosmic rays (CR) far exceeds the fraction of other CR species, such as antiprotons, electrons, and positrons. Thus the majority of information obtained from CR studies is based on interpretation of isotopic abundances using CR propagation models where the nuclear data and isotopic production cross sections in p- and alpha-induced reactions are the key elements. This paper presents an introduction to the astrophysics of CR and diffuse gamma rays and discusses some of the puzzles that have emerged recently due to more precise data and improved propagation models. Merging with cosmology and particle physics, astrophysics of CR has become a very dynamic field with a large potential of breakthrough and discoveries in the near future. Exploiting the data collected by the CR experiments to the fullest requires accurate nuclear cross sections.

The gamma-ray source GRO J0852-4642, discovered by COMPTEL, is a possible counterpart of the supernova remnant RX J0852-4622. Detection of radioactive decay from 44Ti nuclei would proof it to be the youngest and nearest supernova remnant known so far. During the first year of INTEGRAL core program, the Vela region was observed twice in all for more than 2000 ks. Among other nucleosynthesis studies, one of the most important scientific goals of this observation is the detection of 44Ti gamma-ray lines expected at 68 keV, 78 keV and 1157 keV. For this purpose the INTEGRAL Spectrometer (SPI), with its very high energy resolution thanks to its germanium detector camera is the key instrument, permitting a precise determination of gamma-ray line intensities and profiles. The upper limit for the 78.4 keV 44Ti gamma-ray line emission derived from the first analysis is 1.1 x 10^{-4} gamma/cm2s. This value is mainly dominated by systematic uncertainties in the treatment of the instrumental background. By accumulating more observation time in the next years of the mission and by improving the background understanding, a reliable 44Ti flux for GRO J0852-4642 or an upper limit which constrains the COMPTEL flux can be expected.

The spectrometer SPI, one of the two main instruments of the INTEGRAL spacecraft, has strong capabilities in the field of Gamma-Ray Burst (GRB) detections. In its 16 deg field of view (FoV) SPI is able to trigger and to localize GRBs with an accuracy for strong bursts better than 1 deg. The expected GRB detection rate is about one per month.

We present a study of the compatibility of some current models of the diffuse Galactic continuum gamma rays with EGRET data. A set of regions sampling the whole sky is chosen to provide a comprehensive range of tests. The range of EGRET data used is extended to 100 GeV. The models are computed with our GALPROP cosmic-ray propagation and gamma-ray production code. We confirm that the "conventional model" based on the locally observed electron and nucleon spectra is inadequate, for all sky regions. A conventional model plus hard sources in the inner Galaxy is also inadequate, since this cannot explain the GeV excess away from the Galactic plane. Models with a hard electron injection spectrum are inconsistent with the local spectrum even considering the expected fluctuations; they are also inconsistent with the EGRET data above 10 GeV. We present a new model which fits the spectrum in all sky regions adequately. Secondary antiproton data were used to fix the Galactic average proton spectrum, while the electron spectrum is adjusted using the spectrum of diffuse emission itself. The derived electron and proton spectra are compatible with those measured locally considering fluctuations due to energy losses, propagation, or possibly details of Galactic structure. This model requires a much less dramatic variation in the electron spectrum than models with a hard electron injection spectrum, and moreover it fits the gamma-ray spectrum better and to the highest EGRET energies. It gives a good representation of the latitude distribution of the gamma-ray emission from the plane to the poles, and of the longitude distribution. We show that secondary positrons and electrons make an essential contribution to Galactic diffuse gamma-ray emission.

We use the GALPROP model for cosmic-ray propagation to obtain a new estimate of the Galactic component of gamma rays, and show that away from the Galactic plane it gives an accurate prediction of the observed EGRET intensities in the energy range 30 MeV - 50 GeV. On this basis we re-evaluate the extragalactic gamma-ray background. We find that for some energies previous work underestimated the Galactic contribution at high latitudes and hence overestimated the background. Our new background spectrum shows a positive curvature similar to that expected for models of the extragalactic emission based on the blazar population.

We present a solution to the apparent discrepancy between the radial gradient in the diffuse Galactic gamma-ray emissivity and the distribution of supernova remnants, believed to be the sources of cosmic rays. Recent determinations of the pulsar distribution have made the discrepancy even more apparent. The problem is shown to be plausibly solved by a variation in the CO-to-N(H2) scaling factor. If this factor increases by a factor of 5-10 from the inner to the outer Galaxy, as expected from the Galactic metallicity gradient, we show that the source distribution required to match the radial gradient of gamma-ray can be reconciled with the distribution of supernova remnants as traced by current studies of pulsars. The resulting model fits the EGRET gamma-ray profiles extremely well in longitude, and reproduces the mid-latitude inner Galaxy intensities better than previous models.

INTEGRAL Spectrometer SPI data from the first year of the Galactic Centre Deep Exposure has been analysed for the diffuse continuum from the Galactic ridge. A new catalogue of sources from the INTEGRAL Imager IBIS has been used to account for their contribution to the celestial signal. Apparently diffuse emission is detected at a level ~10% of the total source flux. A comparison of the spectrum of diffuse emission with that from an analysis of IBIS data alone shows that they are consistent. The question of the contribution of unresolved sources to this ridge emission is still open.

The nature of the soft gamma-ray (20-200 keV) Galactic emission has been a matter of debate for a long time. Previous experiments have tried to separate the point source contribution from the real interstellar emission, but with a rather poor spatial resolution, they concluded that the interstellar emission could be a large fraction of the total Galactic emission. INTEGRAL, having both high resolution and high sensitivity, is well suited to reassess more precisely this problem. Using the INTEGRAL core program Galactic Center Deep Exposure (GCDE), we estimate the contribution of detected point sources to the total Galactic flux.

Al(26) radioactivity traces recent nucleosynthesis throughout the Galaxy, and is known to be produced in massive stars and novae. The map from its decay gamma-ray line suggests massive stars to dominate, but high-resolution line spectroscopy is expected to supplement imaging of Al(26) source regions and thus to help decide about the Al(26) injection process and interstellar environment, hence about the specific massive-star subgroup and phase which produces interstellar Al(26). The INTEGRAL Spectrometer SPI has observed Galactic Al(26) radioactivity in its 1809 keV gamma-ray line during its first inner-Galaxy survey. Instrumental background lines make analysis difficult; yet, a clear signal from the inner Galaxy agrees with expectations. In particular, SPI has constrained the line width to exclude previously-reported line broadenings corresponding to velocities >500 km/s. The signal-to-background ratio of percent implies that detector response and background modeling need to be fine-tuned to eventually enable line shape deconvolution in order to extract source location information along the line of sight.

We present the preliminary results of the INTEGRAL survey of the Large Magellanic Cloud. The observations have been carried out in January 2003 (about 10^6 s) and January 2004 (about 4x10^5 s). Here we concentrate on the bright sources LMC X-1, LMC X-2, LMC X-3 located in our satellite galaxy, and on the serendipitous detections of the Galactic Low Mass X-ray Binary EXO 0748-676 and of the Seyfert 2 galaxy IRAS 04575-7537.

"Diffuse" gamma rays consist of several components: truly diffuse emission from the interstellar medium, the extragalactic background, whose origin is not firmly established yet, and the contribution from unresolved and faint Galactic point sources. One approach to unravel these components is to study the diffuse emission from the interstellar medium, which traces the interactions of high energy particles with interstellar gas and radiation fields. Because of its origin such emission is potentially able to reveal much about the sources and propagation of cosmic rays. The extragalactic background, if reliably determined, can be used in cosmological and blazar studies. Studying the derived "average" spectrum of faint Galactic sources may be able to give a clue to the nature of the emitting objects.

The SPI instrument has been launched on-board the INTEGRAL observatory on October 17, 2002. SPI is a spectrometer devoted to the sky observation in the 20 keV-8 MeV energy range using 19 germanium detectors. The performance of the cryogenic system is nominal and allows to cool the 19 kg of germanium down to 85 K with a comfortable margin. The energy resolution of the whole camera is 2.5 keV at 1.1 MeV. This resolution degrades with time due to particle irradiation in space. We show that the annealing process allows the recovery of the initial performance. The anticoincidence shield works as expected, with a low threshold at 75 keV, reducing the GeD background by a factor of 20. The digital front-end electronics system allows the perfect alignement in time of all the signals as well as the optimisation of the dead time (12%). We demonstrate that SPI is able to map regions as complex as the galactic plane. The obtained spectrum of the Crab nebula validates the present version of our response matrix. The 3 $\sigma$ sensitivity of the instrument at 1 MeV is 8 10$^{-7}$ph$\cdot$cm$^{-2}\cdot$s$^{-1}\cdot$keV$^{-1}$ for the continuum and 3 10$^{-5}$ph$\cdot$cm$^{-2}\cdot$s$^{-1}$ for narrow lines.

SPI, the Spectrometer on board the ESA INTEGRAL satellite, to be launched in October 2002, will study the gamma-ray sky in the 20 keV to 8 MeV energy band with a spectral resolution of 2 keV for photons of 1 MeV, thanks to its 19 germanium detectors spanning an active area of 500 cm2. A coded mask imaging technique provides a 2 deg angular resolution. The 16 deg field of view is defined by an active BGO veto shield, furthermore used for background rejection. In April 2001 the flight model of SPI underwent a one-month calibration campaign at CEA in Bruy\`eres le Ch\^atel using low intensity radioactive sources and the CEA accelerator for homogeneity measurements and high intensity radioactive sources for imaging performance measurements. After integration of all scientific payloads (the spectrometer SPI, the imager IBIS and the monitors JEM-X and OMC) on the INTEGRAL satellite, a cross-calibration campaign has been performed at the ESA center in Noordwijk. A set of sources has been placed in the field of view of the different instruments in order to compare their performances and determine their mutual influence. Some of those sources had already been used in Bruy\`eres during the SPI standalone test. For the lowest energy band calibration an X-ray generator has been used. We report on the scientific goals of this calibration activity, and present the measurements performed as well as some preliminary results.

We present the analysis of the first observations of the Cygnus region by the SPI spectrometer onboard the Integral Gamma Ray Observatory, encompassing ${\sim}$ 600 ks of data. Three sources namely Cyg X-1, Cyg X-3 and EXO 2030+375 were clearly detected. Our data illustrate the temporal variability of Cyg X-1 in the energy range from 20 keV to 300 keV. The spectral analysis shows a remarkable stability of the Cyg X-1 spectra when averaged over one day timescale. The other goal of these observations is SPI inflight calibration and performance verification. The latest objective has been achieved as demonstrated by the results presented in this paper.

We report the first measurements of the 511 keV line emission from the Galactic Centre (GC) region performed with the spectrometer SPI on the space observatory INTEGRAL (International Gamma-Ray Astrophysics Laboratory). Taking into account the range of spatial distribution models which are consistent with the data, we derive a flux of $9.9^{+4.7}_{-2.1} \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ and an intrinsic line width of $2.95^{+0.45}_{-0.51}$ keV (FWHM). The results are consistent with other high-spectroscopy measurements, though the width is found to be at the upper bound of previously reported values.

Three calibration campaigns of the spectrometer SPI have been performed before launch in order to determine the instrument characteristics, such as the effective detection area, the spectral resolution and the angular resolution. Absolute determination of the effective area has been obtained from simulations and measurements. At 1 MeV, the effective area is 65 cm^2 for a point source on the optical axis, the spectral resolution ~2.3 keV. The angular resolution is better than 2.5 deg and the source separation capability about 1 deg. Some temperature dependant parameters will require permanent in-flight calibration.

The spectrometer SPI, one of the two main instruments of the INTEGRAL spacecraft, offers significant gamma-ray burst detection capabilities. In its 35 deg (full width) field of view SPI is able to localise gamma-ray bursts at a mean rate of ~ 0.8/month. With its large anticoincidence shield of 512 kg of BGO crystals SPI is able to detect gamma-ray bursts quasi omni-directionally with a very high sensitivity. Burst alerts of the anticoincidence shield are distributed by the INTEGRAL Burst Alert System. In the first 8 months of the mission about 0.8/day gamma-ray burst candidates and 0.3/day gamma-ray burst positions were obtained with the anticoincidence shield by interplanetary network triangulations with other spacecrafts.

The INTEGRAL Science Data Centre (ISDC) provides the INTEGRAL data and means to analyse them to the scientific community. The ISDC runs a gamma ray burst alert system that provides the position of gamma ray bursts on the sky within seconds to the community. It operates a quick-look analysis of the data within few hours that detects new and unexpected sources as well as it monitors the instruments. The ISDC processes the data through a standard analysis the results of which are provided to the observers together with their data.

Scans of the Galactic plane performed at regular intervals constitute a key element of the guaranteed time observations of the INTEGRAL observing programme. These scans are done for two reasons: frequent monitoring of the Galactic plane in order to detect transient sources, and time resolved mapping of the Galactic plane in continuum and diffuse line emission. This paper describes first results obtained from the Galactic plane scans executed so far during the early phase (Dec 2002 - May 2003) of the nominal mission.

Recent measurements of the CR antiproton flux have been shown to pose a problem for conventional propagation models. In particular, models consistent with secondary/primary nuclei ratio in CR produce too few antiprotons, while matching the ratio and the antiproton flux requires ad hoc assumptions. This may indicate an additional local CR component or new phenomena in CR propagation in the Galaxy. We discuss several possibilities which may cause this problem.

The origin and evolution of isotopes of the lightest elements H2, He3, Li, Be, B in the universe is a key problem in such fields as astrophysics of CR, Galactic evolution, non-thermal nucleosynthesis, and cosmological studies. One of the major sources of these species is spallation by CR nuclei in the interstellar medium. On the other hand, it is the B/C ratio in CR and Be10 abundance which are used to fix the propagation parameters and thus the spallation rate. We study the production and Galactic propagation of isotopes of elements Z<6 using the numerical propagation code GALPROP and updated production cross sections.

The GALPROP model for cosmic-ray (CR) propagation produces explicit predictions for the angular distribution of Galactic diffuse gamma rays. We compare our current models with EGRET spectra for various regions of the sky. This allows a critical test of alternative hypotheses for the observed GeV excess. We show that a population of hard-spectrum gamma-ray sources cannot be solely responsible for the excess since it also appears at high latitudes; on the other hand a hard CR electron spectrum model cannot explain the gamma-ray excess in the inner Galaxy. Hence some combination of these explanations is suggested.

We use the GALPROP model for cosmic-ray (CR) propagation to obtain our best prediction of the Galactic component of gamma rays, and show that away from the Galactic plane it gives an accurate prediction of the observed EGRET intensities. On this basis we re-evaluate the isotropic extragalactic gamma-ray background (EGRB). We find that for some energies previous work underestimated the Galactic contribution and hence overestimated the background. Our new EGRB spectrum shows a positive curvature similar to that expected for models of the extragalactic emission based on the blazar population.

The spectrometer SPI, one of the two main instruments of the INTEGRAL spacecraft, has strong capabilities in the Field of Gamma-Ray Burst (GRB) detections. In its 16 degree Field of view (FoV) SPI is able to trigger and to localize GRBs. With its large anticoincidence shield (ACS) of 512 kg of BGO crystals SPI is able to detect GRBs quasi omnidirectionally with a very high sensitivity. The ACS GRB alerts will provide GRB arrival times with high accuracy but with no or very rough positional information. The expected GRB detection rate in SPI's FoV will be one per month and for the ACS around 300 per year. At MPE two SPI software contributions to the real-time INTEGRAL burst-alert system (IBAS) at the INTEGRAL science data centre ISDC have been developed. The SPI-ACS branch of IBAS will produce burst alerts and light-curves with 50 ms resolution. It is planned to use ACS burst alerts in the 3rd interplanetary network. The SPI-FoV branch of IBAS is currently under development at MPE. The system is using the energy and timing information of single and multiple events detected by the Germanium-camera of SPI. Using the imaging algorithm developed at the University of Birmingham the system is expected to locate strong bursts with an accuracy of better than 1 degree.

In a recent paper (Moskalenko et al., 2002), it has been shown that the flux of secondary cosmic ray (CR) antiprotons appears to be contradictory to measurements of secondary to primary nuclei ratios in cosmic rays when calculated in the same Galactic propagation model. The contradiction appears as a value of the diffusion coefficient necessary to match the secondary ratios antiproton/proton and B/C. In particlular, it was shown that the reacceleration models designed to match secondary to primary nuclei ratios produce too few antiprotons. It is, however, clear that some reacceleration is unavoidable in the turbulent interstellar medium. Here we discuss an idea of how to improve reacceleration model by allowing for the damping of interstellar turbulences on the small scale by cosmic rays, mostly protons. This would lead to increase in the mean free path lenghts at low energies, the well-known phenomena empirically discovered in the Leaky Box models, thus producing less secondary nuclei. Antiprotons will remain almost non-affected due to their high energy threshold of production cross section.

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) produce too few antiprotons. In the present paper we discuss one possibility to overcome these difficulties. Using the measured antiproton flux AND B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local "unprocessed" component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z<=28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three "radioactive clock" isotopes in CR, 10Be, 26Al, 36Cl, are all consistent and indicate a halo size z_h~4 kpc based on the most accurate data taken by the ACE spacecraft.

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., boron/carbon) produce too few antiprotons. Matching both the secondary to primary nuclei ratio and the antiproton flux requires artificial breaks in the diffusion coefficient and the primary injection spectrum suggesting the need for other approaches. In the present paper we discuss one possibility to overcome these difficulties. Using the measured antiproton flux AND B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local "unprocessed" component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z<=28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three "radioactive clock" isotopes in CR, 10Be, 26Al, 36Cl, are all consistent and indicate a halo size z_h~4 kpc based on the most accurate data taken by the ACE spacecraft.

We have used observations of Cygnus X-1 from the Compton Gamma-Ray Observatory (CGRO) and BeppoSAX to study the variation in the MeV gamma-ray emission between the hard and soft spectral states, using spectra that cover the energy range 20 keV up to 10 MeV. These data provide evidence for significant spectral variability at energies above 1 MeV. In particular, whereas the hard X-ray flux decreases during the soft state, the flux at energies above 1 MeV increases, resulting in a significantly harder gamma-ray spectrum at energies above 1 MeV. This behavior is consistent with the general picture of galactic black hole candidates having two distinct spectral forms at soft gamma-ray energies. These data extend this picture, for the first time, to energies above 1 MeV. We have used two different hybrid thermal/non-thermal Comptonization models to fit broad band spectral data obtained in both the hard and soft spectral states. These fits provide a quantitative estimate of the electron distribution and allow us to probe the physical changes that take place during transitions between the low and high X-ray states. We find that there is a significant increase (by a factor of ~4) in the bolometric luminosity as the source moves from the hard state to the soft state. Furthermore, the presence of a non-thermal tail in the Comptonizing electron distribution provides significant constraints on the magnetic field in the source region.

We investigate the origin of the nonthermal X-ray emission from the Galactic ridge in the range 10-200 keV. We consider bremsstrahlung of subrelativistic cosmic ray protons and electrons as production processes. From the solution of the kinetic equations describing the processes of particle in-situ acceleration and spatial propagation we derive parameters of the spectra for protons and electrons. For in-situ acceleration the flux of accelerated particles consists mainly of protons since the ratio of the accelerated protons to electrons is large and the flux of nuclei with charges Z>1 is strongly suppressed. Our analysis has shown that the proton bremsstrahlung model is free from the previously discussed problems of the high gamma-ray flux at MeV or at hundreds of MeV energy regions. However, the model with proton bremsstrahlung requires specific restrictions on the proton spectrum.The real problem of the proton bremsstrahlung model is the large pressure of the accelerated protons in the disk. For the electron origin of the ridge flux we do not have the problems of the gamma-ray line and continuous fluxes. On the other hand the electron bremsstrahlung model requires more effective acceleration of background particles.

High-energy collisions of cosmic-ray nuclei with interstellar gas are believed to be the mechanism producing the majority of cosmic ray antiprotons. Due to the kinematics of the process they are created with a nonzero momentum; the characteristic spectral shape with a maximum at ~2 GeV and a sharp decrease towards lower energies makes antiprotons a unique probe of models for particle propagation in the Galaxy and modulation in the heliosphere. On the other hand, accurate calculation of the secondary antiproton flux provides a ``background'' for searches for exotic signals from the annihilation of supersymmetric particles and primordial black hole evaporation. Recently new data with large statistics on both low and high energy antiproton fluxes have become available which allow such tests to be performed. We use our propagation code GALPROP to calculate interstellar cosmic-ray propagation for a variety of models. We show that there is no simple model capable of accurately describing the whole variety of data: boron/carbon and sub-iron/iron ratios, spectra of protons, helium, antiprotons, positrons, electrons, and diffuse gamma rays. We find that only a model with a break in the diffusion coefficient plus convection can reproduce measurements of cosmic-ray species, and the reproduction of primaries (p, He) can be further improved by introducing a break in the primary injection spectra. For our best-fit model we make predictions of proton and antiproton fluxes near the Earth for different modulation levels and magnetic polarity using a steady-state drift model of propagation in the heliosphere.

We present the time-averaged characteristics of the Crab pulsar in the 0.75-30 MeV energy window using data from the imaging Compton Telescope COMPTEL aboard the Compton Gamma-Ray Observatory (CGRO) collected over its 9 year mission. Exploiting the exceptionally long COMPTEL exposure on the Crab allowed us to derive significantly improved COMPTEL spectra for the Crab nebula and pulsar emissions, and for the first time to accurately determine at low-energy gamma-rays the pulse profile as a function of energy. These timing data, showing the well-known main pulse and second pulse at a phase separation of 0.4 with strong bridge emission, are studied together with data obtained at soft/hard X-ray energies from the ROSAT HRI, BeppoSAX LECS, MECS and PDS, at soft gamma-rays from CGRO BATSE and at high-energy gamma-rays from CGRO EGRET in order to obtain a coherent high-energy picture of the Crab pulsar from 0.1 keV up to 10 GeV.

We use a new version of our numerical model for particle propagation in the Galaxy to study radioactive secondaries. For evaluation of the production cross sections we use the Los Alamos compilation of all available experimental cross sections together with calculations using the improved Cascade-Exciton Model code CEM2k. Using the radioactive secondary ratios 26Al/27Al, 36Cl/Cl, 54Mn/Mn, we show how the improved cross-section calculations together with the new propagation code allow us to better constrain the size of the CR halo.

The GALPROP cosmic-ray (CR) propagation model has been extended to three dimensions including the effects of stochastic SNR sources, a comprehensive cross-section database, and nuclear reaction networks. A brief description of the new code is given and some illustrative results for the distribution of CR protons presented. Results for electrons and gamma rays are given in an accompanying paper.

High energy collisions of cosmic ray (CR) nuclei with interstellar gas are believed to be the mechanism producing the majority of CR antiprotons. The distinguishing spectral shape with a maximum at 2 GeV and a sharp decrease towards lower energies makes antiprotons a unique probe of the models of particle propagation in the Galaxy and modulation in the heliosphere. Besides, accurate calculation of the secondary antiproton flux provides a ``background'' for searches for exotic signals from the annihilation of supersymmetric particles and primordial black hole evaporation. Recently new data with large statistics on the antiproton flux have become available which allow for such tests to be performed. We use our 3D Galactic cosmic ray propagation code GALPROP to calculate interstellar propagation in several models. For our best model we make predictions of proton and antiproton fluxes near the Earth for different modulation levels and polarity using a steady-state drift model for heliospheric modulation.

In studies of cosmic-ray (CR) propagation and diffuse continuum gamma-ray emission from the Galaxy it has usually been assumed that the source function can be taken as smooth and time-independent. However, especially for electrons at high energies where energy losses are rapid, the effect of the stochastic nature of the sources becomes apparent and indeed has been invoked to explain the GeV excess in the diffuse emission observed by EGRET. In order to address this problem in detail a model with explicit time-dependence and a stochastic SNR population has been developed, which follows the propagation in three dimensions. The results indicate that although the inhomogeneities are large they are insufficient to easily explain the GeV excess. However the fluctuations should show up in the gamma-ray distribution at high energies and this should be observable with GLAST. Estimates of the TeV continuum emission from the plane are consistent with the Whipple upper limit.

We have used observations from CGRO to study the variation in the MeV emission of Cygnus X-1 between its low and high X-ray states. These data provide a measurement of the spectral variability above 1 MeV. The high state MeV spectrum is found to be much harder than that of the low state MeV spectrum. In particular, the power-law emission seen at hard X-ray energies in the high state spectrum (with a photon spectral index of 2.6) is found to extend out to at least 5 MeV, with no evidence for any cutoff. Here we present the data and describe our efforts to model both the low state and high state spectra using a hybrid thermal/nonthermal model in which the emission results from the Comptonization of an electron population that consists of both a thermal and nonthermal component.

Most cosmic ray antiprotons observed near the Earth are secondaries produced in collisions of energetic cosmic ray (CR) particles with interstellar gas. The spectrum of secondary antiprotons is expected to peak at ~2 GeV and decrease sharply at lower energies. This leaves a low energy window in which to look for signatures of exotic processes such as evaporation of primordial black holes or dark matter annihilation. In the inner heliosphere, however, modulation of CRs by the solar wind makes analysis difficult. Detecting these antiprotons outside the heliosphere on an interstellar probe removes most of the complications of modulation. We present a new calculation of the expected secondary antiproton flux (the background) as well as a preliminary design of a light-weight, low-power instrument for the interstellar probe to make such measurements.

A new numerical model of particle propagation in the Galaxy has been developed, which allows the study of cosmic-ray and gamma-ray production and propagation in 2D or 3D, including a full reaction network. This is a further development of the code which has been used for studies of cosmic ray reacceleration, Galactic halo size, antiprotons and positrons in cosmic rays, the interpretation of diffuse continuum gamma rays, and dark matter. In this paper we illustrate recent results focussing on B/C, sub-Fe/Fe, ACE radioactive isotope data, source abundances and antiprotons. From the radioactive nuclei we derive a range of 3-7 kpc for the height of the cosmic-ray halo.

The Compton Gamma-Ray Observatory (CGRO) has observed the Cygnus region on several occasions since its launch in 1991. The data collected by the COMPTEL experiment on CGRO represent the most sensitive observations to date of Cygnus X-1 in the 0.75-30 MeV range. A spectrum accumulated by COMPTEL over 10 weeks of observation time shows significant evidence for emission extending out to several MeV. We have combined these data with contemporaneous data from both BATSE and OSSE to produce a broad-band gamma-ray spectrum, corresponding to the low X-ray state of Cygnus X-1, extending from 50 keV up to approximately 5 MeV. Although there is no evidence for any broad line-like emissions in the MeV region, these data further confirm the presence of a hard tail at energies above several hundred keV. In particular, the spectrum at MeV energies can be described as a power-law with a photon spectral index of $\alpha$ = -3.2, with no evidence for a cutoff at high energies. For the 200 keV to 5 MeV spectrum, we provide a quantitative description of the underlying electron spectrum, in the context of a hybrid thermal/non-thermal model for the emission. The electron spectrum can be described by a thermal Maxwellian with a temperature of $kT_e$ = 86 keV and a non-thermal power-law component with a spectral index of $p_e$ = 4.5. The spectral data presented here should provide a useful basis for further theoretical modeling.

A new study of the diffuse Galactic gamma-ray continuum radiation is presented, using a cosmic-ray propagation model which includes nucleons, antiprotons, electrons, positrons, and synchrotron radiation. Our treatment of the inverse Compton (IC) scattering includes the effect of anisotropic scattering in the Galactic interstellar radiation field (ISRF) and a new evaluation of the ISRF itself. Models based on locally measured electron and nucleon spectra and synchrotron constraints are consistent with gamma-ray measurements in the 30-500 MeV range, but outside this range excesses are apparent. A harder nucleon spectrum is considered but fitting to gamma rays causes it to violate limits from positrons and antiprotons. A harder interstellar electron spectrum allows the gamma-ray spectrum to be fitted above 1 GeV as well, and this can be further improved when combined with a modified nucleon spectrum which still respects the limits imposed by antiprotons and positrons. A large electron/IC halo is proposed which reproduces well the high-latitude variation of gamma-ray emission. The halo contribution of Galactic emission to the high-latitude gamma-ray intensity is large, with implications for the study of the diffuse extragalactic component and signatures of dark matter. The constraints provided by the radio synchrotron spectral index do not allow all of the <30 MeV gamma-ray emission to be explained in terms of a steep electron spectrum unless this takes the form of a sharp upturn below 200 MeV. This leads us to prefer a source population as the origin of the excess low-energy gamma rays.

In previous work, we have used data from the first three years of the CGRO mission to assemble a broad-band $\gamma$-ray spectrum of the galactic black hole candidate Cygnus X-1. Contemporaneous data from the COMPTEL, OSSE and BATSE experiments on CGRO were selected on the basis of the hard X-ray flux (45--140 keV) as measured by BATSE. This provided a spectrum of Cygnus X-1 in its canonical low X-ray state (as measured at energies below 10 keV), covering the energy range from 50 keV to 5 MeV. Here we report on a comparison of this spectrum to a COMPTEL-OSSE spectrum collected during a high X-ray state of Cygnus X-1 (May, 1996). These data provide evidence for significant spectral variability at energies above 1 MeV. In particular, whereas the hard X-ray flux {\it decreases} during the high X-ray state, the flux at energies above 1 MeV {\it increases}, resulting in a significantly harder high energy spectrum. This behavior is consistent with the general picture of galactic black hole candidates having two distinct spectral forms at soft $\gamma$-ray energies. These data extend this picture, for the first time, to energies above 1 MeV.

We use all COMPTEL data from the beginning of the CGRO mission (April '91) up to the end of CGRO Cycle 6 (November '97) to carry out all-sky point source analyses in the four standard COMPTEL energy bands for different time periods. We apply our standard maximum-likelihood method to generate all-sky significance and flux maps for point sources by subtracting off the diffuse emission components via model fitting. In addition, fluxes of known sources have been determined for individual CGRO Phases/Cycles to generate lightcurves with a time resolution of the order of one year. The goal of the analysis is to derive quantitative results -- significances, fluxes, light curves -- of our brightest and most significant sources such as 3C 273, and to search for additional new COMPTEL sources, showing up in time-averaged maps only.

Models for the diffuse Galactic continuum emission and synchrotron radiation show that it is difficult to reproduce observations of both of these from the same population of cosmic-ray electrons. This indicates that an important contributor to the emission below 10 MeV could be an unresolved point-source population. We suggest that these could be Crab-like sources in the inner Galaxy. Alternatively a sharp upturn in the electron spectrum below 200 MeV is required.

Galactic diffuse continuum gamma-ray emission is intricately related to cosmic-ray physics and radio astronomy. We describe recent results from an approach which endeavours to take advantage of this. Information from cosmic-ray composition constrains the propagation of cosmic rays; this in turn can be used as input for gamma-ray models. The GeV gamma-ray excess cannot be explained as neutral pion decay resulting from a hard nucleon spectrum without violating antiproton and positron data; the best explanation at present appears to be inverse-Compton emission from a hard interstellar electron spectrum. One consequence is an increased importance of Galactic inverse Compton for estimates of the extragalactic background. At low energies, an additional point-source component of gamma-rays seems to be necessary.

The inverse Compton scattering of interstellar photons off cosmic-ray electrons seems to play a more important role in the generation of diffuse emission from the Galaxy than thought before. The background radiation field of the Galaxy is highly anisotropic since it is dominated by the radiation from the Galactic plane. An observer in the Galactic plane thus sees mostly head-on scatterings even if the distribution of the cosmic-ray electrons is isotropic. This is especially evident when considering inverse Compton scattering by electrons in the halo, i.e. the diffuse emission at high Galactic latitudes. We derive formulas for this process and show that the anisotropy of the interstellar radiation field has a significant effect on the intensity and angular distribution of the Galactic diffuse gamma rays, which can increase the high-latitude Galactic gamma-ray flux up to 40%. This effect should be taken into account when calculating the Galactic emission for extragalactic background estimates.

We have developed a model which aims to reproduce observational data of many kinds related to cosmic-ray (CR) origin and propagation: direct measurements of nuclei, antiprotons, electrons and positrons, gamma-rays, and synchrotron radiation. Our main results include evaluation of diffusion/convection and reacceleration models, estimates of the halo size, calculations of the interstellar positron and antiproton spectra, evaluation of alternative hypotheses of nucleon and electron interstellar spectra, and computation of the Galactic diffuse gamma-ray emission. Recently our CR propagation code has been generalized to include fragmentation networks of arbitrary complexity. The code can now provide an alternative to leaky-box calculations for full isotopic abundance calculations and has the advantage of including the spatial dimension which is essential for radioactive nuclei. Preliminary predictions for sub-Fe/Fe, 10Be/9Be and 26Al/27Al are presented in anticipation of new experimental isotopic data. We show that combining information from classical CR studies with gamma-ray and other data leads to tighter constraints on CR origin and propagation.

We have made a calculation of the propagation of positrons from dark-matter particle annihilation in the Galactic halo for different models of the dark matter halo distribution using our 3D code. We show that the Green's functions are not very sensitive to the dark matter distribution for the same local dark matter energy density. We compare our predictions with computed cosmic ray positron spectra ("background") for the "conventional" cosmic-ray nucleon spectrum which matches the local measurements, and a modified spectrum which respects the limits imposed by measurements of diffuse Galactic gamma-rays, antiprotons, and positrons. We conclude that significant detection of a dark matter signal requires favourable conditions and precise measurements unless the dark matter is clumpy which would produce a stronger signal. Although our conclusion qualitatively agrees with that of previous authors, it is based on a more realistic model of particle propagation and thus reduces the scope for future speculations. Reliable background evaluation requires new accurate positron measurements and further developments in modelling production and propagation of cosmic ray species in the Galaxy.

The cosmic-ray propagation code GALPROP has been generalized to include fragmentation networks of arbitrary complexity. The code can now provide an alternative to leaky-box calculations for full isotopic abundance calculations and has the advantage of including the spatial dimension which is essential for radioactive nuclei. Preliminary predictions for sub-Fe/Fe, 10Be/9Be, 26Al/27Al, 36Cl/Cl, and 54Mn/Mn are presented in anticipation of new experimental isotopic data.

Recent evidence for a large Galactic halo, based on cosmic-ray radioactive nuclei, implies a significant contribution from inverse Compton emission at high Galactic latitudes. We present predictions for the expected intensity distribution, and show that the EGRET gamma-ray latitude distribution is well reproduced from the plane to the poles. We show that the Galactic component at high latitudes may be comparable to the extragalactic emission in some energy ranges.

We have made a calculation of the propagation of positrons from dark-matter particle annihilation in the Galactic halo in different models of the dark matter halo distribution using our 3D code, and present fits to our numerical propagation Green's functions. We show that the Green's functions are not very sensitive to the dark matter distribution for the same local dark matter energy density. We compare our predictions with computed cosmic ray positron spectra (``background'') for the ``conventional'' CR nucleon spectrum which matches the local measurements, and a modified spectrum which respects the limits imposed by measurements of diffuse Galactic gamma-rays, antiprotons, and positrons. We conclude that significant detection of a dark matter signal requires favourable conditions and precise measurements unless the dark matter is clumpy which would produce a stronger signal. Although our conclusion qualitatively agrees with that of previous authors, it is based on a more realistic model of particle propagation and thus reduces the scope for future speculations. Reliable background evaluation requires new accurate positron measurements and further developments in modelling production and propagation of cosmic ray species in the Galaxy.

We have developed a model which aims to reproduce observational data of many kinds related to cosmic-ray origin and propagation: direct measurements of nuclei, antiprotons, electrons and positrons, gamma-rays, and synchrotron radiation. Our main results include evaluation of diffusion/convection and reacceleration models, estimates of the halo size, calculations of the interstellar positron and antiproton spectra, evaluation of alternative hypotheses of nucleon and electron interstellar spectra, and computation of the Galactic diffuse gamma-ray emission. We show that combining information from classical cosmic-ray studies with gamma-ray and other data leads to tighter constraints on cosmic-ray origin and propagation.

We are constructing a model which aims to reproduce observational data of many kinds related to cosmic-ray origin and propagation: direct measurements of nuclei, antiprotons, electrons and positrons, gamma-rays, and synchrotron radiation. These data provide many independent constraints on any model. Propagation of primary and secondary nucleons, primary and secondary electrons and positrons are calculated self-consistently. Fragmentation and energy losses are computed using realistic distributions for the interstellar gas and radiation fields, and diffusive reacceleration is also incorporated. The models are adjusted to agree with the observed cosmic-ray B/C and 10Be/9Be ratios. Our main results include evaluation of diffusion/convection and reacceleration models, estimates of the halo size, calculations of the interstellar positron and antiproton spectra, evaluation of alternative hypotheses of hard nucleon and hard electron interstellar spectra, and computation of the Galactic diffuse gamma-ray emission.

Inverse Compton scattering appears to play a more important role in the diffuse Galactic continuum emission than previously thought, from MeV to GeV energies. We compare models having a large inverse Compton component with EGRET data, and find good agreement in the longitude and latitude distributions at low and high energies. We test an alternative explanation for the >1 GeV gamma-ray excess, the hard nucleon spectrum, using secondary antiprotons and positrons. At lower energies to fit the COMPTEL and OSSE data as diffuse emission requires either a steep upturn in the electron spectrum below 200 MeV or a population of discrete sources.

Large-scale skymapping with COMPTEL using the full survey database presents challenging problems on account of the complex response and time-variable background. A new approach which attempts to address some of these problems is described, in which the information about each observation is preserved throughout the analysis. In this method, a maximum-entropy algorithm is used to determine image and background simultaneously. Because of the extreme computing requirements, the method has been implemented on a parallel computer, which brings large gains since the response computation is fully parallelizable. The zero level is left undetermined in this method. Results using data from 7 years of COMPTEL data are presented.

We present quantitative statistical evidence for a $\gamma$-ray emission halo surrounding the Galaxy. Maps of the emission are derived. EGRET data were analyzed in a wavelet-based non-parametric hypothesis testing framework, using a model of expected diffuse (Galactic + isotropic) emission as a null hypothesis. The results show a statistically significant large scale halo surrounding the center of the Milky Way as seen from Earth. The halo flux at high latitudes is somewhat smaller than the isotropic gamma-ray flux at the same energy, though of the same order (O(10^(-7)--10^(-6)) ph/cm^2/s/sr above 1 GeV).

The excess of continuum gamma-ray emission from the Galaxy above 1 GeV is an unsolved puzzle. It may indicate that the interstellar nucleon or electron spectra are harder than local direct measurements, as could be the case if a local source of cosmic rays were to dominate the nearby flux. It is however difficult to distinguish between the two cases. Cosmic-ray secondary antiprotons provide a way to resolve this issue. We have made a calculation of the cosmic-ray secondary antiproton spectrum in our model, which computes self-consistently propagation of primary and secondary nucleons, and electrons. Fragmentation and energy losses are computed using realistic distributions for the interstellar gas and radiation fields, and diffusive reacceleration is also incorporated. Our study shows that accurate measurements of the antiproton flux, especially at high energies, could provide a diagnostic of the interstellar nucleon spectrum allowing us to test the hard nucleon spectrum hypothesis. Present antiproton data above 3 GeV indicate that it can already be excluded at the few sigma level.

We describe a method for the numerical computation of the propagation of primary and secondary nucleons, primary electrons, and secondary positrons and electrons. Fragmentation and energy losses are computed using realistic distributions for the interstellar gas and radiation fields, and diffusive reacceleration is also incorporated. The models are adjusted to agree with the observed cosmic-ray B/C and 10Be/9Be ratios. Models with diffusion and convection do not account well for the observed energy dependence of B/C, while models with reacceleration reproduce this easily. The height of the halo propagation region is determined, using recent 10Be/9Be measurements, as >4 kpc for diffusion/convection models and 4-12 kpc for reacceleration models. For convection models we set an upper limit on the velocity gradient of dV/dz < 7 km/s/kpc. The radial distribution of cosmic-ray sources required is broader than current estimates of the SNR distribution for all halo sizes. Full details of the numerical method used to solve the cosmic-ray propagation equation are given.

An extensive program for the calculation of galactic cosmic-ray propagation has been developed. Primary and secondary nucleons, primary and secondary electrons, secondary positrons and antiprotons are included. Fragmentation and energy losses are computed using realistic distributions for the interstellar gas and radiation fields. Models with diffusion and convection only do not account naturally for the observed energy dependence of B/C, while models with reacceleration reproduce this easily. The height of the halo propagation region is determined, using recent 10Be/9Be measurements, as greater than 4 kpc. The radial distribution of cosmic-ray sources required is broader than current estimates of the SNR distribution for all halo sizes. Our results include an estimate of cosmic-ray antiproton and positron spectra, and the Galactic diffuse gamma-ray emission (see accompanying paper: Moskalenko 1998b).

Inverse Compton scattering appears to play a more important role in the diffuse Galactic continuum emission than previously thought, from MeV to GeV energies. We compare models having a large inverse Compton component with EGRET data, and find good agreement in the longitude and latitude distributions at low and high energies. We test an alternative explanation for the >1 GeV gamma-ray excess, the hard nucleon spectrum, using secondary antiprotons and positrons.

We report on a study of the gamma-ray continuum emission from the Crab supernova nebula and on a search for nuclear de-excitation gamma-ray lines. Crab is the brightest continuum source in the 1-10 MeV gamma-ray sky, and its continuum radiation is most likely of synchrotron origin. It is a likely source of cosmic rays through shock acceleration and thus a potential candidate for gamma-ray line emission from nuclear interactions. Five years of COMPTEL observations enable a fine spectral binning to investigate the behaviour of the 0.75-30 MeV emission in detail and to search for nuclear de-excitation lines on top of the continuum. The nebular spectrum shows a break at the edge of the COMPTEL energy range and connects well to the EGRET spectrum, probably reflecting electron energy losses in the synchrotron emission scenario. Such a smooth continuum model alone may not be sufficient to explain the observations. A weak bump in the spectrum at 1-2 MeV may be present. No significant evidence for distinct line emission is seen, but the presence of a blend of line features or another synchrotron component cannot be excluded.

We have made a new calculation of the cosmic-ray secondary positron spectrum using a diffusive halo model for Galactic cosmic-ray propagation. The code computes self-consistently the spectra of primary and secondary nucleons, primary electrons, and secondary positrons and electrons. The models are first adjusted to agree with the observed cosmic-ray Boron/Carbon ratio, and the interstellar proton and Helium spectra are then computed; these spectra are used to obtain the source function for the secondary positrons/electrons which are finally propagated with the same model parameters. The primary electron spectrum is evaluated, again using the same model. Fragmentation and energy losses are computed using realistic distributions for the interstellar gas and radiation fields, and diffusive reacceleration is also incorporated. Our study includes a critical re-evaluation of the secondary decay calculation for positrons. The predicted positron fraction is in good agreement with the measurements up to 10 GeV, beyond which the observed flux is higher than that calculated. Since the positron fraction is now accurately measured in the 1-10 GeV range our primary electron spectrum should be a good estimate of the true interstellar spectrum in this range, of interest for gamma ray and solar modulation studies. We further show that a harder interstellar nucleon spectrum, similar to that suggested to explain EGRET diffuse Galactic gamma ray observations above 1 GeV, can reproduce the positron observations above 10 GeV without requiring a primary positron component.

An extensive program for the calculation of galactic cosmic-ray propagation has been developed. This is a continuation of the work described in Strong & Youssefi (1995). The main motivation for developing this code is the prediction of diffuse Galactic gamma rays for comparison with data from the CGRO instruments EGRET, COMPTEL, and OSSE. The basic spatial propagation mechanisms are (momentum-dependent) diffusion, convection, while in momentum space energy loss and diffusive reacceleration are treated. Primary and secondary nucleons, primary and secondary electrons, and secondary positrons are included. Fragmentation and energy losses are computed using realistic distributions for the interstellar gas and radiation fields. This study indicates that it is possible to construct a model satisfying a wide range of observational constraints and provides a basis for future developments.

An extensive program for the calculation of galactic cosmic-ray propagation has been developed. Primary and secondary nucleons, primary and secondary electrons, and secondary positrons are included. The basic spatial propagation mechanisms are (momentum-dependent) diffusion, convection, while in momentum space energy loss and diffusive reacceleration are treated. Fragmentation and energy losses are computed using realistic distributions for the interstellar gas and radiation fields. Predictions of the diffuse Galactic gamma ray and synchrotron emission are made. Our new estimate of the cosmic-ray positron spectrum is compared with the latest experimental data.

We present an analysis of COMPTEL observations made between November 1991 and May 1994 of 2CG 135+01, a bright gamma-ray source located near the Galactic plane. At energies above 1 MeV, an excess consistent with the position of 2CG 135+01 is detected in the sum of the observations, at flux levels which are a factor of 10-100 below those published in the past. The detection significance of this excess, when the possible presence of underlying Galactic diffuse emission is neglected, is 6.6 sigma for 3 degrees of freedom. The differential photon spectrum in the 1-30 MeV energy range can be described by a power law with a spectral index of $1.95^{+0.2}_{-0.3}$. Due to the uncertainties involved in modelling the Galactic-disk diffuse emission underneath the source, the absolute flux levels must be considered uncertain by a factor of two. They are consistent with the extrapolation of the time-averaged spectrum of 2CG 135+01 measured with EGRET, thereby strengthening the identification. No significant temporal correlation between the gamma-ray emission and the monitored radio emission of the possible counterpart radio source GT 0236+610 (showing a 26.5 day modulation) is found.