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.

Simulating and analysing detailed observations of astrophysical sources for very high energy (VHE) experiments, like the Cherenkov Telescope Array (CTA), can be a demanding task especially in terms of CPU consumption and required storage. In this context, we propose an innovative cloud computing architecture based on Amazon Web Services (AWS) aiming to decrease the amount of time required to simulate and analyse a given field by distributing the workload and exploiting the large computational power offered by AWS. We detail how the various services offered by the Amazon online platform are jointly used in our architecture and we report a comparison of the execution times required for simulating observations of a test source with the CTA, by a single machine and the cloud-based approach. We find that, by using AWS, we can run our simulations more than 2 orders of magnitude faster than by using a general purpose workstation for the same cost. We suggest to consider this method when observations need to be simulated, analysed, and concluded within short timescales.

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.

e-ASTROGAM (`enhanced ASTROGAM') is a breakthrough Observatory mission dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. In the largely unexplored MeV-GeV domain, 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 Galactic ecosystems. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, 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 the promise of eLISA. Keywords: High-energy gamma-ray astronomy, High-energy astrophysics, Nuclear Astrophysics, Compton and Pair creation telescope, Gamma-ray bursts, Active Galactic Nuclei, Jets, Outflows, Multiwavelength observations of the Universe, Counterparts of gravitational waves, Fermi, Dark Matter, Nucleosynthesis, Early Universe, Supernovae, Cosmic Rays, Cosmic antimatter.

We report the Fermi Large Area Telescope (LAT) detections of high-energy (>100 MeV) gamma-ray emission from two recent optically bright classical novae, V1369 Centauri 2013 and V5668 Sagittarii 2015. At early times, Fermi target-of-opportunity observations prompted by their optical discoveries provided enhanced LAT exposure that enabled the detections of gamma-ray onsets beginning ~2 days after their first optical peaks. Significant gamma-ray emission was found extending to 39-55 days after their initial LAT detections, with systematically fainter and longer duration emission compared to previous gamma-ray detected classical novae. These novae were distinguished by multiple bright optical peaks that encompassed the timespans of the observed gamma rays. The gamma-ray light curves and spectra of the two novae are presented along with representative hadronic and leptonic models, and comparisons to other novae detected by the LAT are discussed.

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.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

The Cherenkov Telescope Array (CTA) will be the largest cosmic gamma ray detector ever built in the world. It will be installed at two different sites in the North and South hemispheres and should be operational for about 30 years. In order to cover the desired energy range, the CTA is composed of typically 50-100 collecting telescopes of various sizes (from 6 to 24-m diameters). Most of them are equipped with a focal plane camera consisting of 1500 to 2000 Photomultipliers (PM) equipped with light concentrating optics, whose double function is to maximize the amount of Cherenkov light detected by the photo-sensors, and to block any stray light originating from the terrestrial environment. Two different optical solutions have been designed, respectively based on a Compound Parabolic Concentrator (CPC), and on a purely dioptric concentrating lens. In this communication are described the technical specifications, optical designs and performance of the different solutions envisioned for all these light concentrators. The current status of their prototyping activities is also given.

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

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

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

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

Laue lenses are an emerging technology based on diffraction in crystals that allows the concentration of soft gamma rays. This kind of optics that works in the 100 keV - 1.5 MeV band can be used to realize an high-sensitivity and high-angular resolution telescope (in a narrow field of view). This paper reviews the recent progresses that have been done in the development of efficient crystals, in the design study and in the modelisation of the answer of Laue lenses. Through the example of a new concept of 20 m focal length lens focusing in the 100 keV - 600 keV band, the performance of a telescope based on a Laue lens is presented. This lens uses the most efficient mosaic crystals in each sub-energy range in order to yield the maximum reflectivity. Imaging capabilities are investigated and shows promising results.

Apart from the Milky Way, the Large Magellanic Cloud (LMC) is the only other normal star-forming galaxy that was conclusively detected in high energy (> 100 MeV) gamma rays by the Energetic Gamma Ray Telescope (EGRET) on the Compton Gamma-Ray Observatory. However, the sensitivity of EGRET was sufficient only to marginally resolve the LMC. We report on measurements of the gamma-ray emission from the LMC by the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. For the first time an externally viewed star-forming galaxy is well resolved in gamma rays. We discuss the distribution of the LMC diffuse emission as seen by the LAT and implications for cosmic-ray physics.

The concept of a gamma-ray telescope based on a Laue lens offers the possibility to increase the sensitivity by more than an order of magnitude with respect to existing instruments. Laue lenses have been developed by our collaboration for several years : the main achievement of this R&D program was the CLAIRE lens prototype. Since then, the endeavour has been oriented towards the development of efficient diffracting elements (crystal slabs), the aim being to step from a technological Laue lens to a scientifically exploitable lens. The latest mission concept featuring a gamma-ray lens is the European Gamma-Ray Imager (GRI) which intends to make use of the Laue lens to cover energies from 200 keV to 1300 keV. Investigations of two promising materials, low mosaicity copper and gradient concentration silicon-germanium are presented in this paper. The measurements have been performed during three runs on beamline ID15A of the European Synchrotron Radiation Facility, and on the GAMS 4 instrument of the Institute Laue-Langevin (both in Grenoble, France) using highly monochromatic beam of energy close to 500 keV. Despite it was not perfectly homogeneous, the presented copper crystal exhibits peak reflectivity of 25% in accordance with theoretical predictions, and a mosaicity around 26 arcsec, the ideal range for the realization of a Laue lens such as GRI. Silicon-germanium featuring a constant gradient have been measured for the very first time at 500 keV. Two samples showed a quite homogeneous reflectivity reaching 26%, which is far from the 48% already observed in experimental crystals but a very encouraging beginning. This results have been used to estimate the performance of the GRI Laue lens design.

We report the detection of both the 67.9 and 78.4 keV 44Sc gamma-ray lines in Cassiopeia A with the INTEGRAL IBIS/ISGRI instrument. Besides the robustness provided by spectro-imaging observations, the main improvements compared to previous measurements are a clear separation of the two 44Sc lines together with an improved significance of the detection of the hard X-ray continuum up to 100 keV. These allow us to refine the determination of the 44Ti yield and to constrain the nature of the nonthermal continuum emission. By combining COMPTEL, BeppoSAX/PDS and ISGRI measurements, we find a line flux of (2.5 +/- 0.3)*10(-5) cm(-2) s(-1) leading to a synthesized 44Ti mass of 1.6 (+0.6-0.3)*10(-4) solar mass. This high value suggests that Cas A is peculiar in comparison to other young supernova remnants, from which so far no line emission from 44Ti decay has been unambiguously detected.

We present a spectral analysis of the e+e- annihilation emission from the Galactic Centre region based on the first year of measurements made with the spectrometer SPI of the INTEGRAL mission. We have found that the annihilation spectrum can be modelled by the sum of a narrow and a broad 511 keV line plus an ortho-Ps continuum. The broad line is detected with a flux of (0.35+/-0.11)e-3 s-1 cm-2. The measured width of 5.4+/-1.2 keV FWHM is in agreement with the expected broadening of 511 keV photons emitted in the annihilation of Ps that are formed by the charge exchange process of slowing down positrons with H atoms. The flux of the narrow line is (0.72+/-0.12)e-3 s-1 cm-2 and its width is 1.3+/-0.4 keV FWHM. The measured ortho-Ps continuum flux yields a fraction of Ps of (96.7+/-2.2)%. To derive in what phase of the interstellar medium positrons annihilate, we have fitted annihilation models calculated for each phase to the data. We have found that 49(+2,-23)% of the annihilation emission comes from the warm neutral phase and 51(+3,-2)% from the warm ionized phase. While we may not exclude that less than 23% of the emission might come from cold gas, we have constrained the fraction of annihilation emission from molecular clouds and hot gas to be less than 8% and 0.5%, respectively. We have compared our knowledge of the interstellar medium in the bulge and the propagation of positrons with our results and found that they are in good agreement if the sources are diffusively distributed and if the initial kinetic energy of positrons is lower than a few MeV. Despite its large filling factor, the lack of annihilation emission from the hot gas is due to its low density, which allows positrons to escape this phase.

The history of X-ray astronomy has shown that any advancement in our knowledge of the X-ray sky is strictly related to an increase in instrument sensitivity. At energies above 60 keV, there are interesting prospects for greatly improving the limiting sensitivity of the current generation of direct viewing telescopes (with or without coded masks), offered by the use of Laue lenses. We will discuss below the development status of a Hard X-Ray focusing Telescope (HAXTEL) based on Laue lenses with a broad bandpass (from 60 to 600 keV) for the study of the X-ray continuum of celestial sources. We show two examplesof multi-lens configurations with expected sensitivity orders of magnitude better ($\sim 1 \times 10^{-8}$ photons cm$^{-2}$ s$^{-1}$ keV$^{-1}$ at 200 keV) than that achieved so far. With this unprecedented sensitivity, very exciting astrophysical prospects are opened.

The status of the SPI energy calibration after the first year of INTEGRAL operations is reported. We have studied the gain variations and we have demonstrated that the most important parameter is the germanium detector temperature. This study permits us to determine the limits of our calibration method and the frequency of calibrations needed.

It is believed that core-collapse supernovae (CCSN), occurring at a rate about once per century, have seeded the interstellar medium with long-lived radioisotopes such as 60Fe (half-life 1.5 Myr), which can be detected by the gamma rays emitted when they beta-decay. Here we report the detection of the 60Fe decay lines at 1173 keV and 1333 keV with fluxes 3.7 +/- 1.1 x 10(-5) ph cm(-2) s(-1) per line, in spectra taken by the SPI spectrometer on board INTEGRAL during its first year. The same analysis applied to the 1809 keV line of 26Al yielded a line flux ratio 60Fe/26Al = 0.11 +/- 0.03. This supports the hypothesis that there is an extra source of 26Al in addition to CCSN.

The status of the analysis of the electron/positron annihilation radiation performed by INTEGRAL/SPI in the 4th GALACTIC QUADRANT is reported. We use data from the first two Galactic Center Deep Exposures (GCDE) and from the Galactic Plane Scans (GPS). The analysis presented here is focused on the spectroscopic aspects of the electron/positron annihilation radiation. Background substraction and model fitting methods are described, and the parameters of the 511 keV line (flux, energy, and width) are deduced.

The center of our Galaxy is a known strong source of electron-positron 511-keV annihilation radiation. Thus far, however, there have been no reliable detections of annihilation radiation outside of the central radian of our Galaxy. One of the primary objectives of the INTEGRAL (INTErnational Gamma-RAy Astrophysics Laboratory) mission, launched in Oct. 2002, is the detailed study of this radiation. The Spectrometer on INTEGRAL (SPI) is a high resolution coded-aperture gamma-ray telescope with an unprecedented combination of sensitivity, angular resolution and energy resolution. We report results from the first 10 months of observation. During this period a significant fraction of the observing time was spent in or near the Galactic Plane. No positive annihilation flux was detected outside of the central region (|l| > 40 deg) of our Galaxy. In this paper we describe the observations and data analysis methods and give limits on the 511-keV flux.

In order to allow for an efficient and flexible scientific analysis of data from the SPI imaging spectrometer aboard INTEGRAL, I developed a set of analysis executables that are publicly available through the internet. The software is fully compatible with the ISDC data format. It complements software that is actually available through ISDC, and will be included in the next ISDC software release. This paper describes the design of my software system and provides a brief introduction to the executables.

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.

The search for gamma-ray line emission from the radioactive decay of 60Fe figures among the prime scientific objectives of the INTEGRAL mission. 60Fe is believed to be primarily produced in core-collapse supernovae, although other sites, such as carbon deflagration supernovae or intermediate mass AGB stars have also been suggested. We present first results of our search for the 1173 and 1332 keV gamma-ray lines of 60Fe decay in the SPI data of the first mission year. So far we can only report upper flux limits whose levels are determined by systematic uncertainties in the treatment of the instrumental background. Once we understand the background uncertainties better, the SPI data should provide stringent constraints on 60Fe nucleosynthesis in our Galaxy.

We provide first constraints on the morphology of the 511 keV line emission from the galactic centre region on basis of data taken with the spectrometer SPI on the INTEGRAL gamma-ray observatory. The data suggest an azimuthally symmetric galactic bulge component with FWHM of ~9 deg with a 2 sigma uncertainty range covering 6-18 deg. The 511 keV line flux in the bulge component amounts to (9.9+4.7-2.1) 10e-4 ph cm-2 s-1. No evidence for a galactic disk component has been found so far; upper 2 sigma flux limits in the range (1.4-3.4) 10e-3 ph cm-2 s-1 have been obtained that depend on the assumed disk morphology. These limits correspond to lower limits on the bulge-to-disk ratio of 0.3-0.6.

With the launch of ESA's INTEGRAL satellite in october 2002, a gamma-ray observatory will be placed in orbit providing a multiwavelength coverage from a few keV up to 10 MeV for the study of high energy phenomena in the universe. Among the major scientific objectives of INTEGRAL are the study of compact objects, and in particular of microquasars, stellar nucleosynthesis, high energy transients, particle acceleration and interaction processes, and galactic diffuse emission. In this review I will expose the science prospects of INTEGRAL in the above mentioned fields.

With the launch of ESA's INTEGRAL satellite in october 2002, a gamma-ray observatory will become available to the scientific community that combines imaging and spectroscopic capacities in the 20 keV to 10 MeV energy range. In this paper, we summarise the science prospects of the SPI spectrometer aboard INTEGRAL which provides unprecedented spectral resolution with good imaging capabilities and high sensitivity. Emphasise will be given to the key objectives of SPI which are the determination of the galactic positron origin, the study of galactic star formation and nucleosynthesis activity via the radioactive trace isotopes 26Al and 60Fe, and the study of nucleosynthesis in supernova explosions and the gamma-ray line emission from their remnants.

Gamma-ray and microwave observations of the Cygnus region reveal an intense signal of 1.809 MeV line emission, attributed to radioactive decay of 26Al, that is closely correlated with 53 GHz free-free emission, originating from the ionised interstellar medium. We modelled both emissions using a multi-wavelength evolutionary synthesis code for massive star associations that we applied to the known massive star populations in Cygnus. For all OB associations and young open clusters in the field, we determined the population age, distance, and richness as well as the uncertainties in all these quantities from published photometric and spectroscopic data. We propagate the population uncertainties in model uncertainties by means of a Bayesian method. The young globular cluster Cyg OB2 turns out to be the dominant 26Al nucleosynthesis and ionisation source in Cygnus. Our model reproduces the ionising luminosity of the Cygnus region very well, yet it underestimates 26Al production by about a factor of 2. We attribute this underestimation to shortcomings of current nucleosynthesis models, and suggest the inclusion of stellar rotation as possible mechanism to enhance 26Al production. We also modelled 60Fe nucleosynthesis in the Cygnus region, yet the small number of recent supernova events suggests only little 60Fe production. Consequently, a detection of the 1.137 MeV and 1.332 MeV decay lines of 60Fe from Cygnus by the upcoming INTEGRAL observatory is not expected.

Evolutionary synthesis models are a fundamental tool to interpret the properties of observed stellar systems. In order to achieve a meaningful comparison between models and real data, it is necessary to calibrate the models themselves, i.e. to evaluate the dispersion due to the discreteness of star formation as well as the possible model errors. In this paper we show that linear interpolations in the log M - log t_k plane, that are customary in the evaluation of isochrones in evolutionary synthesis codes, produce unphysical results. We also show that some of the methods used in the calculation of time-integrated quantities (kinetic energy, and total ejected masses of different elements) may produce unrealistic results. We propose alternative solutions to solve both problems. Moreover, we have quantified the expected dispersion of these quantities due to stochastic effects in stellar populations. As a particular result, we show that the dispersion in the 14N/12C ratio increases with time.

The detection of 1.809 MeV gamma-ray line emission from the Cygnus X complex by the COMPTEL telescope is one of the most convincing proves of massive star nucleosynthesis in our Galaxy. The Cygnus X complex is an extremely active nearby region, containing several star forming regions, OB associations and young open star clusters. It houses some of the most massive stars known in our Galaxy and concentrates a large number of Wolf-Rayet stars. Thermal radio continuum emission and intense [C II] line emission reveals widespread ionisation, and at least 60 individual H II regions have been identified. In order to understand the 1.809 MeV gamma-ray line emission from the Cygnus X complex, and to compare the observations to theoretical nucleosynthesis calculations, we modelled the multi-wavelength spectrum of the region by means of an evolutionary synthesis model. Our investigation leads us to the following conclusions: 1. Stellar wind ejection is the dominant mechanism for the observed 26Al enrichment in the Cygnus X region. 2. Cyg OB2 is by far the dominating massive star association in Cygnus X and 1.809 MeV line emission from 26Al produced in this association should be detectable by the spectrometer SPI on INTEGRAL. 3. There is only low supernova activity in the Cygnus X complex and the 60Fe lines should be below INTEGRAL's detection sensitivity. The detectability of an individual massive star cluster (Cyg OB2) by INTEGRAL provides a fantastic opportunity for future nucleosynthesis studies using gamma-ray line spectroscopy. We will explore the scientific potential of such observations, and provide estimates for gamma-ray line intensity distributions based on the massive star census of the Cygnus X region.

The BATSE mission aboard CGRO can be used to observe hard X-ray sources by using the Earth occultation method. This method relies on measuring a step in the count rate profile in each BATSE detector as a source rises above or sets below the Earth's limb. A major problem in determining the step sizes (and hence the flux) is in extracting the steps from the varying background. A technique for flat-fielding the response of gamma ray detectors has been developed at Southampton. The technique uses a dynamic Monte-Carlo model to simulate the dominant components of the gamma-ray background encountered by the experiment at any point in its orbit. A maximum likelihood imaging method is also being developed that will be used to make a sky survey with all 9 years of the BATSE Continuous data set in the 20 - 500 keV range. An all sky map of 25 - 35 keV emission has been made using 60 days of data and has a 3 sigma flux sensitivity of 14 mCrab.

After the recent beautiful results on gamma-ray lines obtained with CGRO, the INTEGRAL mission with the imaging-spectrometer SPI will set the next milestone, combining improved sensitivity and angular resolution with a considerable increase in spectral resolution. SPI is expected to provide significant new information on galactic nucleosynthesis processes and star formation activity, as traced by the distributions of annihilation radiation and radioactive isotopes such as 26Al and 60Fe. The unprecedented spectral resolution will allow the study of dynamic processes in stellar mass ejections and will provide access to kinematic distance estimates for gamma-ray line sources. The study of supernovae and their remnants will be prime objectives for SPI observations. Nearby type Ia SN, within 15 Mpc or so, are in reach of the instrument and a few such events are expected during the lifetime of INTEGRAL. Young galactic supernova remnants, possibly hidden by interstellar dust, may be unveiled by their characteristic gamma-ray line signature from the radioactive decay of 44Ti, as has been demonstrated by COMPTEL for Cas A and possibly RX J0852.0-4622. Classical novae are also among the SPI targets, which may observe the gamma-ray lines from radioactive 7Be and 22Na. Such observations can constrain the physics of the nova explosions and will allow to evaluate their role as nucleosynthesis sites. The interaction of cosmic rays with the dense matter in molecular clouds may be another source of gamma-ray lines that is potentially accessible to SPI. Finally after the SIGMA results on Nova Muscae and 1E1740.7-2942, and a possible 2.223 MeV line detection by COMPTEL, the search for lines from X novae is another way to participate in the understanding of the physical conditions in these close binary systems.

We have developed a new diagnostic tool for the study of gamma-ray emission lines from radioactive isotopes (such as 26Al and 60Fe in conjunction with other multi-wavelength observables of Galactic clusters, associations, and alike objects. Our evolutionary synthesis models are based on the code of Cervino & Mas-Hesse (1994), which has been updated to include recent stellar evolution tracks, new stellar atmospheres for OB and WR stars, and nucleosynthetic yields from massive stars during hydrostatic burning phases and explosive SN II and SN Ib events. The temporal evolution of 26Al and 60Fe production, the equivalent yield of 26Al per ionising O7 V star (yal), and other observables are predicted for a coeval population. We also employed our model to estimate the steady state nucleosynthesis of a population of solar metallicity. To allow for a fully quantitative analysis of existing and future multi-wavelength observations, we propose a Bayesian approach that allows the inclusion of IMF richness effects and observational uncertainties in the analysis. In particular, a Monte Carlo technique is adopted to estimate probability distributions for all observables of interest. We outline the procedure of exploiting these distributions by applying our model to a fictive massive star association. Applications to existing observations of the Cygnus and Vela regions will be discussed in companion papers.

The morphology and stellar content of the Cygnus OB2 association has been determined using 2MASS infrared observations in the J, H, and K bands. The analysis reveals a spherically symmetric association of 2 deg in diameter with a half light radius of 13', corresponding to Rh = 6.4 pc at an assumed distance of 1.7 kpc. The interstellar extinction for member stars ranges from Av = 5m to 20m, which led to a considerable underestimation of the association size and richness in former optical studies. From the infrared colour-magnitude diagram, the number of OB member stars is estimated to 2600 +/- 400, while the number of O stars amounts to 120 +/- 20. This is the largest number of O stars ever found in a galactic massive star association. The slope of the initial mass function has been determined from the colour-magnitude diagram to Gamma=-1.6 +/- 0.1. The total mass of Cyg OB2 is estimated to (4-10)*10^4 Msol, where the primary uncertainty comes from the unknown lower mass cut-off. Using the radial density profile of the association, the central mass density is determined to rho0=40-150 Msol pc^-3. Considering the mass, density, and size of Cyg OB2 it seems untenable to classify this object still as OB association. Cygnus OB2 more closely resembles a young globular cluster like those observed in the Large Magellanic Cloud or in extragalactic star forming regions. It is therefore suggested to re-classify Cygnus OB2 as young globular cluster - an idea which goes back to Reddish et al. (1966). Cygnus OB2 would then be the first object of this class in the Milky Way.

Gamma-ray line observations provide a versatile tool for studies of nucleosynthesis processes and supernova physics. In particular, the observation of radioactive species in the interstellar medium probes recent nucleosynthesis activity on various time-scales for different kinds of sources. Considerable progress in gamma-ray instrumentation during the last decades has led to the discovery of several cosmic gamma-ray lines. In this review, recent observational results are presented and their astrophysical implications are discussed. Prospects of gamma-ray line astronomy will be explored in view of the future INTEGRAL mission.

Gamma-ray line observations provide a versatile tool for studies of nucleosynthesis processes and supernova physics. In particular, the observation of radioactive species in the interstellar medium probes recent nucleosynthesis activity on various time-scales for different kinds of sources. Considerable progress in gamma-ray instrumentation during the last decades has led to the discovery of several cosmic gamma-ray lines. The best studied of these lines is today the 1.809 MeV line attributed to the decay of radioactive 26Al within the interstellar medium. In this review, recent observational results are presented and their astrophysical implications are discussed.

OB associations and young open clusters constitute the most prolific nucleosynthesis sites in our Galaxy. The combined activity of stellar winds and core-collapse supernovae ejects significant amounts of freshly synthesised nuclei into the interstellar medium. Radioactive isotopes, such as 26Al or 60Fe, that have been co-produced in such events may eventually be observed by gamma-ray instruments through their characteristic decay-line signatures. In order to study such signatures and to derive constraints on the involved nucleosynthesis processes, we developed a multi-wavelength evolutionary synthesis model for massive star associations. We apply this model to associations and clusters in the well-studied Cygnus region for which we re-determined the stellar census based on photometric and spectroscopic data. In particular we study the relation between 1.809 MeV gamma-ray line emission and ionising flux, since the latter has turned out to provide an excellent tracer of the global galactic 1.809 MeV emission. We compare our model to COMPTEL 1.8 MeV gamma-ray line observations from which we derive limits on the relative contributions from massive stars and core-collapse supernovae to the actual 26Al content in this region. Based on our model we make predictions about the expected 26Al and 60Fe line signatures in the Cygnus region. These predictions make the Cygnus region a prime target for the future INTEGRAL mission.

We present a new algorithm, called Multiresolution Regularized Expectation Maximization (MREM), for the reconstruction of gamma-ray intensity maps from COMPTEL data. The algorithm is based on the iterative Richardson-Lucy scheme to which we added a wavelet thresholding step in order to eliminate image-noise in the reconstruction. The wavelet thresholding explicitly accounts for spatial correlations in the data, and adapts the angular resolution locally, depending on the significance of the signal in the data. We compare the performance of MREM to that of the maximum entropy and the Richardson-Lucy algorithms by means of Monte-Carlo simulations of COMPTEL 1.809 MeV gamma-ray line observations. The simulations demonstrate that the maximum entropy and Richardson-Lucy algorithms provide virtually identical reconstructions which are heavily disturbed by image noise. MREM largely suppresses this noise in the reconstructions, showing only the significant structures that are present in the data. Application of MREM to COMPTEL 1.8 MeV gamma-ray line data results in a 1.809 MeV sky map that is much smoother than the maximum entropy or Richardson-Lucy reconstructions presented previously. The essential features of this map are (1) an asymmetric galactic ridge emission reaching from l=45 deg to l=240 deg, (2) a bright localised emission feature in the Cygnus region around (l,b)=(80 deg,0 deg), (3) two emission spots at l=317 deg and l=332 deg situated in the galactic plane, and (4) an extended emission region around (l,b)=(160 deg, 0 deg). Comparison of the MREM map to the simulated reconstructions demonstrates that the 1.809 MeV emission is confined to the galactic plane.

We present an evolutionary model for gamma-ray line emission from superbubbles based on shock acceleration of metal-rich stellar ejecta. Application to the Orion OB1 association shows that gamma-ray lines at the detection threshold of the SPI spectrometer aboard of INTEGRAL are expected, making this region an interesting target for studies of the interaction of supernova shocks with the interstellar medium.

The comparison of COMPTEL 1.8 MeV gamma-ray line data to an extended data base of tracer maps revealed an extremely close correlation between 1.8 MeV gamma-ray line emission on the one side and 53 GHz microwave free-free emission and 158 um far-infrared line emission on the other side (Knodlseder et al., these proceedings). We demonstrate that the two tracer maps are correlated to the extreme Population I, tracing the distribution of very massive stars in the Galaxy. This strongly suggests that these stars are also the dominant source of galactic 26Al, making Wolf-Rayet stars and/or core collapse supernovae the primary candidates.

Deconvolving COMPTEL gamma-ray data into images presents a major methodological challenge. We developed a new algorithm called Multiresolution Regularized Expectation Maximization (MREM), which explicitly accounts for spatial correlations in the image by using wavelets (Dixon et al. 1998). We demonstrate that MREM largely suppresses image noise in the reconstruction, showing only significant structures that are present in the data. Application to 1.8 MeV data results in a sky map that is much smoother than the maximum entropy reconstructions presented previously, but it shows the same characteristic emission features which have been established earlier.

We examined the question of a possible correlation between 1.8 MeV emission and intensity distributions observed at other wavelengths by comparing the 1.8 MeV data to an extended set of all-sky maps, covering the entire explored wavelength range from the radio band up to high-energy gamma-rays. This analysis revealed that tracers of the old stellar population or the local interstellar medium provide only a poor description of the data. Tracers of the young stellar population considerably improve the fit. Residuals are minimized for the free-free emission map obtained by the DMR at microwave wavelengths and the 158 um fine-structure line map of C+ observed by FIRAS. Within the statistics of the present data, both maps provide an overall satisfactory fit. We therefore conclude that adequate tracers of 26Al sources have now been identified. Implications for 26Al source types are discussed (see also Knodlseder, these proceedings).