Fast radio bursts (FRBs) are short (millisecond) radio pulses originating from enigmatic sources at extragalactic distances so far lacking a detection in other energy bands. Magnetized neutron stars (magnetars) have been considered as the sources powering the FRBs, but the connection is controversial because of differing energetics and the lack of radio and X-ray detections with similar characteristics in the two classes. We report here the detection by the AGILE satellite on April 28, 2020 of an X-ray burst in coincidence with the very bright radio burst from the Galactic magnetar SGR 1935+2154. The burst detected by AGILE in the hard X-ray band (18-60 keV) lasts about 0.5 seconds, it is spectrally cutoff above 80 keV, and implies an isotropically emitted energy ~ $10^{40}$ erg. This event is remarkable in many ways: it shows for the first time that a magnetar can produce X-ray bursts in coincidence with FRB-like radio bursts; it also suggests that FRBs associated with magnetars may emit X-ray bursts of both magnetospheric and radio-pulse types that may be discovered in nearby sources. Guided by this detection, we discuss SGR 1935+2154 in the context of FRBs, and especially focus on the class of repeating-FRBs. Based on energetics, magnetars with fields B ~ $10^{15}$ G may power the majority of repeating-FRBs. Nearby repeating-FRBs offer a unique occasion to consolidate the FRB-magnetar connection, and we present new data on the X-ray monitoring of nearby FRBs. Our detection enlightens and constrains the physical process leading to FRBs: contrary to previous expectations, high-brightness temperature radio emission coexists with spectrally-cutoff X-ray radiation.

We focus on two repeating fast radio bursts (FRBs) recently detected by the CHIME/FRB experiment in 2018--2019 (Source 1: 180916.J0158+65, and Source 2: 181030.J1054+73). These sources have low excess dispersion measures (DMs) ($ < 100 \rm \, pc \, cm^{-3}$ and $ < 20 \rm \, pc \, cm^{-3}$, respectively), implying relatively small maximal distances. They were repeatedly observed by AGILE in the MeV--GeV energy range. We do not detect prompt emission simultaneously with these repeating events. This search is particularly significant for the submillisecond and millisecond integrations obtainable by AGILE. The sources are constrained to emit a MeV-fluence in the millisecond range below $F'_{MeV} = 10^{-8} \, \rm erg \, cm^{-2}$ corresponding to an isotropic energy near $E_{MeV,UL} \simeq 2 \times 10^{46}\,$erg for a distance of 150 Mpc (applicable to Source 1). We also searched for $\gamma$-ray emission for time intervals up to 100 days, obtaining 3$\,\sigma$ upper limits (ULs) for the average isotropic luminosity above 50 MeV, $L_{\gamma,UL} \simeq \,$(5-10)$\,\times 10^{43} \rm \, erg \, s^{-1}$. For a source distance near 100 kpc (possibly applicable to Source 2), our ULs imply $E_{MeV,UL}\simeq10^{40} \rm erg$, and $L_{\gamma,UL} \simeq \,$2$\,\times 10^{37} \rm \, erg \, s^{-1}$. Our results are significant in constraining the high-energy emission of underlying sources such as magnetars, or other phenomena related to extragalactic compact objects, and show the prompt emission to be lower than the peak of the 2004 magnetar outburst of SGR 1806-20 for source distances less than about 100 Mpc.

The Fermi-LAT telescope has measured the extragalactic gamma-ray background (EGB) generated by the ensemble of all extragalactic sources. The energy distribution of the EGB is well described as a power-law with a spectral index approximately equal to 2.3, and an exponential cutoff, that is consistent with being the effect of absorption of high energy photons assuming an emission that is an unbroken power-law spectrum. The observations of the Fermi telescope have also resolved the EGB, determining that most of it is formed by the sum of the contributions of Active Galactic Nuclei (AGN) of the blazar class. The spectra of the individual AGN sources have a broad range of spectral shapes, and the brightest and most precisely measured sources have curved spectra that can be fitted with the "log-parabola" (or log-normal) form. It might appear surprising that the sum of contributions with different spectral shapes generate a spectrum that has a simple power-law form. We suggest that this fact reveals some important properties for the ensemble of the extragalactic sources, and indicates that the blazar high energy emission can be considered as a "critical phenomenon". A curved log--parabola (log-normal) form is also required to fit the spectra of many Galactic gamma-ray sources, including the emission from the brightest supernova remnants (SNR), and an intriguing possibility is that also the Galactic cosmic rays are injected in interstellar space by their sources not with a single universal power-law spectrum, as predicted by the commonly accepted models, but instead with a broad range of spectral shapes.

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.

Recent observations of cosmic ray protons in the energy range $10^2$--$10^5$~GeV have revealed that the spectrum cannot be described by a simple power law. A hardening of the spectrum around an energy of order few hundred~GeV, first observed by the magnetic spectrometers PAMELA and AMS02, has now been confirmed by several calorimeter detectors (ATIC, CREAM, CALET, NUCLEON and DAMPE). These new measurements reach higher energy and indicate that the hardening corresponds to a larger step in spectral index than what estimated by the magnetic spectrometers. Data at still higher energy (by CREAM, NUCLEON and DAMPE) show that the proton spectrum undergoes a marked softening at $E \approx 10^4$~GeV. Understanding the origin of these unexpected spectral features is a significant challenge for models of the Galactic cosmic rays. An important open question is whether additional features are present in the proton spectrum between the softening and the "Knee". Extensive Air Shower detectors, using unfolding procedures that require the modeling of cosmic ray showers in the atmosphere, estimated the proton flux below and around the Knee (at $E \simeq 3$~PeV). These results however have large systematic uncertainties and are in poor agreement with each other. The measurement in the PeV energy range, recently presented by IceTop/IceCube, indicates a proton flux higher than extrapolations of the direct measurements calculated assuming a constant slope, and therefore requires the existence of an additional spectral hardening below the Knee. A clarification of this point is very important for an understanding of the origin of the Galactic cosmic rays, and is also essential for a precise calculation of the spectra of atmospheric neutrinos in the energy range ($E \gtrsim 10$~TeV) where they constitute the foreground for the emerging astrophysical $\nu$ signal.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Recent precision measurements of the flux of cosmic ray positrons by the Alpha Magnetic Spectrometer show that the spectrum has a marked softening feature for energies close to one TeV. A possible interpretation of this result is that the observed feature measures the maximum energy of a new hard source of positrons perhaps associated to dark matter self--annihilation or decay, or to positron accelerators. A gradual hardening of the positron flux centered at $E \simeq 25$~GeV can also be understood as the signature of the transition where the new source overtakes the conventional component due to secondary production. This interpretation is simple and attractive, but it is not unique. The alternative possibility, that the positron flux is entirely of secondary origin, remains viable. In such a scenario the spectral softening observed by AMS for positrons is generated by energy loss effects, and a feature of similar, but not identical structure should be also visible in the $e^-$ spectrum. Spectral features similar to both the hardening and softening of the positron flux are in fact observed for electrons and call for a consistent explanation. Precision measurements of the $e^+$ and $e^-$ spectra in the TeV and multi--TeV energy range are crucial to clarify the problem.

The origin of cosmic neutrinos is still largely unknown. Using data obtained by the gamma-ray imager on board of the AGILE satellite, we systematically searched for transient gamma-ray sources above 100 MeV that are temporally and spatially coincident with ten recent high-energy neutrino IceCube events. We find three AGILE candidate sources that can be considered possible counterparts to neutrino events. Detecting 3 gamma-ray/neutrino associations out of 10 IceCube events is shown to be unlikely due to a chance coincidence. One of the sources is related to the BL Lac source TXS 0506+056. For the other two AGILE gamma-ray sources there are no obvious known counterparts, and both Galactic and extragalactic origin should be considered.

The cosmic ray energy spectra encode very important information about the mechanisms that generate relativistic particles in the Milky Way, and about the properties of the Galaxy that control their propagation. Relativistic electrons and positrons traveling in interstellar space lose energy much more rapidly than more massive particles such as protons and nuclei, with a rate that grows quadratically with the particle energy $E$. One therefore expects that the effects of energy loss should leave observable signatures in the $e^\mp$ spectra, in the form of softenings centered at the critical energy $E^*$. This quantity is determined by the condition that the total energy loss suffered by particles during their residence time in the Galaxy is of the same order of the initial energy. If the electrons and positrons are accelerated in discrete (quasi) point--like astrophysical objects, such as Supernova explosions or Pulsars, the stochastic nature of the sources should also leave observable signatures in the $e^\mp$ energy spectra and angular distributions above a second (higher) critical energy $E^\dagger$, determined by the condition that particles with $E \gtrsim E^\dagger$ can propagate only for a maximum distance shorter than the average separation between sources. In this work we discuss the theoretically expectations for the signatures of the energy loss effects on the electron and positron spectra, and compare these predictions with the existing observations. Recent measurements of the $(e^- + e^+$) flux have discovered the existence of a prominent spectral break at $E \simeq 1$~TeV. This spectral feature can perhaps be identified with the critical energy $E^*$. An alternative hypothesis is to assmume that $E^*$ has a much lower value of order few GeV. Resolving this ambiguity is of great importance for our understanding of Galactic cosmic rays.

The observation of the diffuse Galactic gamma ray flux is the most powerful tool to study cosmic rays in different regions of the Galaxy, because the energy and angular distributions of the photons encode information about the density and spectral shape of relativistic particles in the entire Milky Way. An open problem of fundamental importance is whether cosmic rays in distant regions of the Milky Way have the same spectral shape observed at the Earth or not. If the spectral shape of protons and nuclei is equal in all the Galaxy, the dominant, hadronic component of the diffuse gamma ray flux must have an angular distribution that, after correcting for absorption effects, is energy independent. To study experimentally the validity of this factorization of the energy and angular dependence of the diffuse flux it is necessary to compare observations in a very broad energy range. The extension of the observations to energies $E_\gamma \simeq 0.1$--10~PeV is of great interest, because it allows the study of the cosmic ray spectra around the feature known as the "knee". The absorption probability for photons in this energy range is not negligible, and distorts the energy and angular distributions of the diffuse flux, therefore a precise calculation of the absorption effects is necessary for the interpretation of the data. In this work we present predictions of the diffuse gamma ray flux at very high energy, constructed under different hypothesis for the space dependence of the cosmic ray energy spectra, and discuss the potential of the observations for present and future detectors.

This paper discusses the relation between the study of the fluxes of cosmic rays, gamma rays and neutrinos, and the connection of these observations with the newly born field of gravitational wave astronomy.

The cosmic ray energy distributions contain spectral features, that is narrow energy regions where the slope of the spectrum changes rapidly. The identification and study of these features is of great importance to understand the astrophysical mechanisms of acceleration and propagation that form the spectra. In first approximation a spectral feature in often described as a discontinuous change in slope, however very valuable information is also contained in its width, that is the length of the interval in logarithm of energy where the change in spectral index develops. In this work we discuss the best way to define and parametrize the width a spectral feature, and for illustration discuss some of the most prominent known structures.

The spectral shape of cosmic ray positrons and antiprotons has been accurately measured in the broad kinetic energy range 1-350 GeV. In the higher part of this range (E > 30 GeV) the e+ and pbar are both well described by power laws with spectral indices gamma[e+] = 2.77 +-0.02 and gamma[pbar] = 2.78 +- 0.04 that are approximately equal to each other and to the spectral index of protons. In the same energy range the positron/antiproton flux ratio has the approximately constant value 2.04+-0.04, that is consistent with being equal to the ratio e_/pbar calculated for the conventional mechanism of production, where the antiparticles are created as secondaries in the inelastic interactions of primary cosmic rays with interstellar gas. The positron/antiproton ratio at lower energy is significantly higher (reaching the approximate value e+/pbar = 100 for E around 1 GeV), but in the entire energy range 1-350 GeV, the flux ratio is consistent with being equal to ratio of the production rates in the conventional mechanism, as the production of low energy antiprotons is kinematically suppressed in collisions with a target at rest. These results strongly suggest that cosmic ray positrons and antiprotons have a common origin as secondaries in hadronic interactions. This conclusion has broad implications for the astrophysics of cosmic rays in the Galaxy.

Galactic Gamma ray astronomy at very high energy (E > 30 TeV) is a vital tool in the study of the non-thermal universe. The interpretation of the observations in this energy region requires the precise modeling of the attenuation of photons due to pair production interactions, where the targets are the radiation fields present in interstellar space. For gamma rays with energy E > 300 TeV the attenuation is mostly due to the photons of the Cosmic Microwave Background Radiation (CMBR). At lower energy the most important target are infrared photons with wavelengths in the range 50-500 micron emitted by dust. The evaluation of the attenuation requires a good knowledge of the density, and energy and angular distributions of the target photons for all positions in the Galaxy. In this work we discuss a simple model for the infrared radiation that depends on only few parameters associated to the space and temperature distributions of the emitting dust. The model allows to compute with good accuracy the effects of absorption for any space and energy distribution of the diffuse Galactic gamma ray emission. The absorption probability due to the Galactic infrared radiation is maximum for E around 150 TeV, and can be as large as P_abs around 0.45 for distant sources on lines of sight that pass close to the Galactic Center. The systematic uncertainties on the absorption probability are estimated as Delta P_abs < 0.08.

The standard way to model the cosmic ray solar modulations is via the Parker equation, that is as the effect of diffusion in the turbulent magnetic of an expanding solar wind. Calculations performed with this method that do not include a description of the regular magnetic field in the heliosphere predict, in disagreement with the observations, equal modulations for particles and antiparticles. The effects of the regular heliospheric field, that break the the particle/anti-particle symmetry, have been included in the Parker equation adding convection terms associated to the magnetic drift velocity of charged particles moving in non--homogeneous magnetic field. In this work we take a completely different approach and study the propagation of charged particles in the heliosphere assuming only the presence of the regular magnetic field, and completely neglecting the random component. Assuming that the field is purely magnetic in the wind frame, one can deduce the existence of a large scale electric field, that results in important energy losses for charged particles that traverse the heliosphere and reach the Earth. The energy loss (Delta E) due to the large scale regular electric field depends on the particle electric charge, is proportional to the absolute value of the electric charge, and to a good approximation is independent from the particle energy and direction. We speculate that this deterministic energy loss is an important, or perhaps even the leading contribution to the solar modulation effects.

The most reliable method to estimate the residence time of cosmic rays in the Galaxy is based on the study of the suppression, due to decay, of the flux of unstable nuclei such as beryllium-10, that have lifetime of appropriate duration. The Cosmic Ray Isotope Spectrometer (CRIS) collaboration has measured the ratio between the fluxes of beryllium-10 and beryllium-9 in the energy range E_0 \simeq 70-145 MeV/nucleon, and has used the data to estimate an escape time tau_{ esc} = 15.0 +- 1.6 Myr. This widely quoted result has been obtained in the framework of a simple leaky-box model where the distributions of escape time and age for stable particles in the Galaxy are identical and have exponential form. In general, the escape time and age distributions do not coincide, they are not unique (because they depend on the injection or observation point), and do not have a simple exponential shape. It is therefore necessary to discuss the measurement of the beryllium ratio in a framework that is more general and more realistic than the leaky-box model. In this work we compute the escape time and age distributions of cosmic rays in the Galaxy in a model based on diffusion that is much more realistic than the simple leaky-box, but that remains sufficiently simple to have exact analytic solutions. Using the age distributions of the model to interpret the measurements of the beryllium-10 suppression, one obtains a cosmic ray residence time that is significantly longer (a factor 2 to 4 depending on the extension of the cosmic ray halo) than the leaky-box estimate. This revised residence time implies a proportional reduction of the power needed to generate the galactic cosmic rays.

High-energy phenomena in the cosmos, and in particular processes leading to the emission of gamma- rays in the energy range 10 MeV - 100 GeV, play a very special role in the understanding of our Universe. This energy range is indeed associated with non-thermal phenomena and challenging particle acceleration processes. The technology involved in detecting gamma-rays is challenging and drives our ability to develop improved instruments for a large variety of applications. GAMMA-LIGHT is a Small Mission which aims at an unprecedented advance of our knowledge in many sectors of astrophysical and Earth studies research. The Mission will open a new observational window in the low-energy gamma-ray range 10-50 MeV, and is configured to make substantial advances compared with the previous and current gamma-ray experiments (AGILE and Fermi). The improvement is based on an exquisite angular resolution achieved by GAMMA-LIGHT using state-of-the-art Silicon technology with innovative data acquisition. GAMMA-LIGHT will address all astrophysics issues left open by the current generation of instruments. In particular, the breakthrough angular resolution in the energy range 100 MeV - 1 GeV is crucial to resolve patchy and complex features of diffuse sources in the Galaxy as well as increasing the point source sensitivity. This proposal addresses scientific topics of great interest to the community, with particular emphasis on multifrequency correlation studies involving radio, optical, IR, X-ray, soft gamma-ray and TeV emission. At the end of this decade several new observatories will be operational including LOFAR, SKA, ALMA, HAWK, CTA. GAMMA-LIGHT will "fill the vacuum" in the 10 MeV-10 GeV band, and will provide invaluable data for the understanding of cosmic and terrestrial high-energy sources.

We present a variability study of a sample of bright gamma-ray (30 MeV -- 50 GeV) sources. This sample is an extension of the first AGILE catalogue of gamma-ray sources (1AGL), obtained using the complete set of AGILE observations in pointing mode performed during a 2.3 year period from July 9, 2007 until October 30, 2009. The dataset of AGILE pointed observations covers a long time interval and its gamma-ray data archive is useful for monitoring studies of medium-to-high brightness gamma-ray sources. In the analysis reported here, we used data obtained with an improved event filter that covers a wider field of view, on a much larger (about 27.5 months) dataset, integrating data on observation block time scales, which mostly range from a few days to thirty days. The data processing resulted in a better characterized source list than 1AGL was, and includes 54 sources, 7 of which are new high galactic latitude (|BII| >= 5) sources, 8 are new sources on the galactic plane, and 20 sources from the previous catalogue with revised positions. Eight 1AGL sources (2 high-latitude and 6 on the galactic plane) were not detected in the final processing either because of low OB exposure and/or due to their position in complex galactic regions. We report the results in a catalogue of all the detections obtained in each single OB, including the variability results for each of these sources. In particular, we found that 12 sources out of 42 or 11 out of 53 are variable, depending on the variability index used, where 42 and 53 are the number of sources for which these indices could be calculated. Seven of the 11 variable sources are blazars, the others are Crab pulsar+nebula, LS I +61{\deg}303, Cyg X-3, and 1AGLR J2021+4030.

Context: AGILE is a gamma-ray astrophysics mission which has been in orbit since 23 April 2007 and continues to operate reliably. The gamma-ray detector, AGILE-GRID, has observed Galactic and extragalactic sources, many of which were collected in the first AGILE Catalog. Aims: We present the calibration of the AGILE-GRID using in-flight data and Monte Carlo simulations, producing Instrument Response Functions (IRFs) for the effective area A_eff), Energy Dispersion Probability (EDP), and Point Spread Function (PSF), each as a function of incident direction in instrument coordinates and energy. Methods: We performed Monte Carlo simulations at different gamma-ray energies and incident angles, including background rejection filters and Kalman filter-based gamma-ray reconstruction. Long integrations of in-flight observations of the Vela, Crab and Geminga sources in broad and narrow energy bands were used to validate and improve the accuracy of the instrument response functions. Results: The weighted average PSFs as a function of spectra correspond well to the data for all sources and energy bands. Conclusions: Changes in the interpolation of the PSF from Monte Carlo data and in the procedure for construction of the energy-weighted effective areas have improved the correspondence between predicted and observed fluxes and spectra of celestial calibration sources, reducing false positives and obviating the need for post-hoc energy-dependent scaling factors. The new IRFs have been publicly available from the Agile Science Data Centre since November 25, 2011, while the changes in the analysis software will be distributed in an upcoming release.

Recently the IceCube collaboration has reported the observation of 28 contained events with a visible energy in the interval between 60 TeV and 1.5 PeV, and has argued that this detection is evidence, with a statistical significance of more than four standard deviations, for the existence of an astrophysical neutrino flux that accounts for a large fraction of the events. In this work we analyze the arguments that allow to identify a component of astrophysical origin in the high energy neutrino flux separating it from atmospheric neutrinos. An astrophysical origin for a large fraction of the IceCube contained events is the simplest and most natural explanation of the data but, conservatively, an atmospheric origin cannot yet be entirely ruled out. This ambiguity should soon be resolved.

In 2008 AGILE and Fermi detected gamma-ray flaring activity from the unidentified EGRET source 3EG J1236+0457, recently associated with a flat spectrum radio quasar GB6 J1239+0443 at z=1.762. The optical counterpart of the gamma-ray source underwent a flux enhancement of a factor 15-30 in 6 years, and of ~10 in six months. We interpret this flare-up in terms of a transition from an accretion-disk dominated emission to a synchrotron-jet dominated one. We analysed a Sloan Digital Sky Survey (SDSS) archival optical spectrum taken during a period of low radio and optical activity of the source. We estimated the mass of the central black hole using the width of the CIV emission line. In our work, we have also investigated SDSS archival optical photometric data and UV GALEX observations to estimate the thermal-disk emission contribution of GB6 J1239+0443. Our analysis of the gamma-ray data taken during the flaring episodes indicates a flat gamma-ray spectrum, with an extension of up to 15 GeV, with no statistically-relevant sign of absorption from the broad line region, suggesting that the blazar-zone is located beyond the broad line region. This result is confirmed by the modeling of the broad-band spectral energy distribution (well constrained by the available multiwavelength data) of the flaring activity periods and by the accretion disk luminosity and black hole mass estimated by us using archival data.

During its first 2 years of operation, the gamma-ray AGILE satellite accumulated an extensive dataset for the Galactic plane. The data have been monitored for transient sources and several gamma-ray sources were detected. Their variability and possible association were studied. In this talk we will focus on the results of extensive observations of the Carina Region during the time period 2007 July - 2009 January, for a total livetime of ~130 days. The region is extremely complex, hosting massive star formation, with the remarkable colliding wind binary Eta Carinae, massive star clusters and HII regions (e.g. NGC 3324, RCW49, Westerlund II) and a giant molecular cloud extending over 150 pc (between l=284.7 and l=289). The Carina Nebula itself is the largest and IR highest surface brightness nebula of the Southern emisphere. We monitored several gamma ray sources in the Carina Region. In particular we detect a gamma ray source (1AGL J1043-5931) consistent with the position of Eta Carinae and report a remarkable 2-days gamma-ray flaring episode from this source on 2008 Oct 11-13. If 1AGL J1043-5931 is associated with the Eta Car system, our data provides the long sought first detection above 100 MeV of a colliding wind binary.

At the core of the AGILE scientific instrument, designed to operate on a satellite, there is the Gamma Ray Imaging Detector (GRID) consisting of a Silicon Tracker (ST), a Cesium Iodide Mini-Calorimeter and an Anti-Coincidence system of plastic scintillator bars. The ST needs an on-ground calibration with a $\gamma$-ray beam to validate the simulation used to calculate the energy response function and the effective area versus the energy and the direction of the $\gamma$ rays. A tagged $\gamma$-ray beam line was designed at the Beam Test Facility (BTF) of the INFN Laboratori Nazionali of Frascati (LNF), based on an electron beam generating $\gamma$ rays through bremsstrahlung in a position-sensitive target. The $\gamma$-ray energy is deduced by difference with the post-bremsstrahlung electron energy \cite{prest}-\cite{hasan}. The electron energy is measured by a spectrometer consisting of a dipole magnet and an array of position sensitive silicon strip detectors, the Photon Tagging System (PTS). The use of the combined BTF-PTS system as tagged photon beam requires understanding the efficiency of $\gamma$-ray tagging, the probability of fake tagging, the energy resolution and the relation of the PTS hit position versus the $\gamma$-ray energy. This paper describes this study comparing data taken during the AGILE calibration occurred in 2005 with simulation.

We report the AGILE gamma-ray observations and the results of the multiwavelength campaigns on seven flaring blazars detected by the mission: During two multiwavelength campaigns, we observed gamma-ray activity from two Flat Spectrum Radio Quasars of the Virgo region, e.g. 3C 279 and 3C 273 (the latter being the first extragalactic source simultaneously observed with the gamma-ray telescope and the hard X ray imager of the mission). Due to the large FOV of the AGILE/GRID instrument, we achieved an almost continuous coverage of the FSRQ 3C 454.3. The source showed flux above 10E-6 photons/cm2/s (E > 100 MeV) and showed day by day variability during all the AGILE observing periods. In the EGRET era, the source was found in high gamma-ray activity only once. An other blazar, PKS 1510-089 was frequently found in high gamma-ray activity. S5 0716+71, an intermediate BL Lac object, exhibited a very high gamma-ray activity and fast gamma-ray variability during a period of intense optical activity. We observed high gamma-ray activity from W Comae, a BL Lac object, and Mrk 421, an high energy peaked BL Lac object. For this source, a multiwavelength campaign from optical to TeV has been performed.

Since its early phases of operation, the AGILE mission is successfully observing Gamma Ray Bursts (GRBs) in the hard X-ray band with the SuperAGILE imager and in the MeV range with the Mini-Calorimeter. Up to now, three firm GRB detections were obtained above 25 MeV and some bursts were detected with lower statistical confidence in the same energy band. When a GRB is localized, either by SuperAGILE or Swift/BAT or INTEGRAL/IBIS or Fermi/GBM or IPN, inside the field of view of the Gamma Ray Imager of AGILE, a detection is searched for in the gamma ray band or an upper limit is provided. A promising result of AGILE is the detection of very short gamma ray transients, a few ms in duration and possibly identified with Terrestrial Gamma-ray Flashes. In this paper we show the current status of the observation of Gamma Ray Bursts and Terrestrial Gamma-ray Flashes with AGILE.

The AGILE scientific instrument has been calibrated with a tagged $\gamma$-ray beam at the Beam Test Facility (BTF) of the INFN Laboratori Nazionali di Frascati (LNF). The goal of the calibration was the measure of the Point Spread Function (PSF) as a function of the photon energy and incident angle and the validation of the Monte Carlo (MC) simulation of the silicon tracker operation. The calibration setup is described and some preliminary results are presented.

Cygnus X-3 (Cyg X-3) is a well-known microquasar producing variable emission at all wavelengths. Cyg X-3 is a prominent X-ray binary producing relativistic jets, and studying its high energy emission is crucial for the understanding of the fundamental acceleration processes in accreting compact objects. Aims. Our goal is to study extreme particle acceleration and {\gamma}-ray production above 100 MeV during special spectral states of Cyg X- 3 usually characterized by a low hard X-ray flux and enhanced soft X-ray states. We observed Cyg X-3 with the AGILE satellite in extended time intervals from 2009 Jun.-Jul., and 2009 Nov.-2010 Jul. We report here the results of the AGILE {\gamma}-ray monitoring of Cyg X-3 as well as the results from extensive multiwavelength campaigns involving radio (RATAN-600, AMI-LA and Mets\"{a}hovi Radio Observatories) and X-ray monitoring data (XTE and Swift). We detect a series of repeated {\gamma}-ray flaring activity from Cyg X-3 that correlate with the soft X-ray states and episodes of decreasing or non-detectable hard X-ray emission. Furthermore, we detect {\gamma}-ray enhanced emission that tends to be associated with radio flares greater than 1 Jy at 15 GHz, confirming a trend already detected in previous observations. The source remained active above 100 MeV for an extended period of time (almost 1.5 months in 2009 Jun.-Jul. and 1 month in 2010 May). We study in detail the short timescale {\gamma}-ray flares that occurred before or near the radio peaks. Our results confirm the transient nature of the extreme particle acceleration from the microquasar Cyg X-3. A series of repeated {\gamma}-ray flares shows correlations with radio and X-ray emission confirming a well established trend of emission. We compare our results with Fermi-LAT and MAGIC TeV observations of Cyg X-3.

The observation of Gamma Ray Bursts (GRBs) in the gamma-ray band has been advanced by the AGILE and Fermi satellites after the era of the Compton Gamma-Ray Observatory. AGILE and Fermi are showing that the GeV-bright GRBs share a set of common features, particularly the high fluence from the keV up to the GeV energy bands, the high value of the minimum Lorentz factor, the presence of an extended emission of gamma-rays, often delayed with respect to lower energies, and finally the possible presence of multiple spectral components. GRB 100724B, localised in a joint effort by Fermi and the InterPlanetary Newtork, is the brightest burst detected in gamma-rays so far by AGILE. Characteristic features of GRB 100724B are the simultaneous emissions at MeV and GeV, without delayed onset nor time lag as shown by the analysis of the cross correlation function, and the significant spectral evolution in hard X-rays over the event duration. In this paper we show the analysis of the AGILE data of GRB 100724B and we discuss its features in the context of the bursts observed so far in gamma-rays and the recently proposed models.

The workshop on "Hadron-Hadron and Cosmic-Ray Interactions at multi-TeV Energies" held at the ECT* centre (Trento) in Nov.-Dec. 2010 gathered together both theorists and experimentalists to discuss issues of the physics of high-energy hadronic interactions of common interest for the particle, nuclear and cosmic-ray communities. QCD results from collider experiments -- mostly from the LHC but also from the Tevatron, RHIC and HERA -- were discussed and compared to various hadronic Monte Carlo generators, aiming at an improvement of our theoretical understanding of soft, semi-hard and hard parton dynamics. The latest cosmic-ray results from various ground-based observatories were also presented with an emphasis on the phenomenological modeling of the first hadronic interactions of the extended air-showers generated in the Earth atmosphere. These mini-proceedings consist of an introduction and short summaries of the talks presented at the meeting.

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.

The cosmic ray interaction event generator Sibyll is widely used in extensive air shower simulations for cosmic ray and neutrino experiments. Charmed particle production has been added to the Monte Carlo with a phenomenological, non-perturbative model that properly accounts for charm production in the forward direction. As prompt decays of charm can become a significant background for neutrino detection, proper simulation of charmed particles is very important. We compare charmed meson and baryon production to accelerator data.

The well known Crab Nebula is at the center of the SN1054 supernova remnant. It consists of a rotationally-powered pulsar interacting with a surrounding nebula through a relativistic particle wind. The emissions originating from the pulsar and nebula have been considered to be essentially stable. Here we report the detection of strong gamma-ray (100 MeV-10 GeV) flares observed by the AGILE satellite in September, 2010 and October, 2007. In both cases, the unpulsed flux increased by a factor of 3 compared to the non-flaring flux. The flare luminosity and short timescale favor an origin near the pulsar, and we discuss Chandra Observatory X-ray and HST optical follow-up observations of the nebula. Our observations challenge standard models of nebular emission and require power-law acceleration by shock-driven plasma wave turbulence within a ~1-day timescale.

Identification of gamma-ray-emitting Galactic sources is a long-standing problem in astrophysics. One such source, 1AGL J2022+4032, coincident with the interior of the radio shell of the supernova remnant Gamma Cygni (SNR G78.2+2.1) in the Cygnus Region, has recently been identified by Fermi as a gamma-ray pulsar, LAT PSR J2021+4026. We present long-term observations of 1AGL J2022+4032 with the AGILE gamma-ray telescope, measuring its flux and light curve. We compare the light curve of 1AGL J2022+4032 with that of 1AGL J2021+3652 (PSR J2021+3651), showing that the flux variability of 1AGL J2022+4032 appears to be greater than the level predicted from statistical and systematic effects and producing detailed simulations to estimate the probability of the apparent observed variability. We evaluate the possibility that the gamma-ray emission may be due to the superposition of two or more point sources, some of which may be variable, considering a number of possible counterparts. We consider the possibility of a nearby X-ray quiet microquasar contributing to the flux of 1AGL J2022+4032 to be more likely than the hypotheses of a background blazar or intrinsic gamma-ray variabilty of LAT PSR J2021+4026.

The observation of astrophysical neutrinos requires a detailed understanding of the atmospheric neutrino background. Since neutrinos are produced in meson decays together with a charged lepton, important constraints on this background can be obtained from the measurement of the atmospheric muon flux. Muons, however, can also be produced as mu+ mu- pairs by purely electromagnetic processes. We use the Z-moment method to study and compare the contributions to the atmospheric muon and neutrino fluxes from different sources (pi/K decay, charmed and unflavored hadron decay, and photon conversion into a muon pair). We pay special attention to the contribution from unflavored mesons (eta, eta', rho0, omega and phi). These mesons are abundant in air showers, their lifetimes are much shorter than those of charged pions or kaons, and they have decay branching ratios of order 10^-4 into final states containing a muon pair. We show that they may be the dominant source of muons at E_mu >10^3 TeV.

Terrestrial Gamma-Ray Flashes (TGFs) are very short bursts of high energy photons and electrons originating in Earth's atmosphere. We present here a localization study of TGFs carried out at gamma-ray energies above 20 MeV based on an innovative event selection method. We use the AGILE satellite Silicon Tracker data that for the first time have been correlated with TGFs detected by the AGILE Mini-Calorimeter. We detect 8 TGFs with gamma-ray photons of energies above 20 MeV localized by the AGILE gamma-ray imager with an accuracy of 5-10 degrees at 50 MeV. Remarkably, all TGF-associated gamma rays are compatible with a terrestrial production site closer to the sub-satellite point than 400 km. Considering that our gamma rays reach the AGILE satellite at 540 km altitude with limited scattering or attenuation, our measurements provide the first precise direct localization of TGFs from space.

We present the results of new Agile observations of PSR B1509-58 performed over a period of 2.5 years following the detection obtained with a subset of the present data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier by Comptel up to 30 MeV: a broad asymmetric first peak reaching its maximum 0.39 +/- 0.02 cycles after the radio peak plus a second peak at 0.94 +/- 0.03. The gamma-ray spectral energy distribution of the pulsed flux detected by Comptel and Agile is well described by a power-law (photon index alpha=1.87+/-0.09) with a remarkable cutoff at E_c=81 +/- 20 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The pulsar luminosity at E > 1 MeV is $L_{\gamma}=4.2^{+0.5}_{-0.2} \times10^{35}$ erg/s, assuming a distance of 5.2 kpc, which implies a spin-down conversion efficiency to gamma-rays of $\sim 0.03$. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In this interpretation our spectrum constrains the magnetic altitude of the emission point(s) at 3 km above the neutron star surface, implying that the attenuation may not be as strong as formerly suggested because pair production can substitute photon splitting in regions of the magnetosphere where the magnetic field becomes too low to sustain photon splitting. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce Agile lightcurves and spectra match the polarization measurements.

We present the gamma-ray data of the extraordinary flaring activity above 100 MeV from the flat spectrum radio quasar 3C 454.3 detected by AGILE during the month of December 2009. 3C 454.3, that has been among the most active blazars of the FSRQ type since 2007, was detected in the gamma-ray range with a progressively rising flux since November 10, 2009. The gamma-ray flux reached a value comparable with that of the Vela pulsar on December 2, 2009. Remarkably, between December 2 and 3, 2009 the source more than doubled its gamma-ray emission and became the brightest gamma-ray source in the sky with a peak flux of F_{\gamma,p} = (2000 \pm 400) x 10^-8 ph cm^-2 s^-1 for a 1-day integration above 100 MeV. The gamma-ray intensity decreased in the following days with the source flux remaining at large values near F \simeq (1000 \pm 200) x 10^-8 ph cm^-2 s^-1 for more than a week. This exceptional gamma-ray flare dissipated among the largest ever detected intrinsic radiated power in gamma-rays above 100 MeV (L_{\gamma, source, peak} \simeq 3 x 10^46 erg s^-1, for a relativistic Doppler factor of {\delta} \simeq 30). The total isotropic irradiated energy of the month-long episode in the range 100 MeV - 3 GeV is E_{\gamma,iso} \simeq 10^56 erg. We report the intensity and spectral evolution of the gamma-ray emission across the flaring episode. We briefly discuss the important theoretical implications of our detection.

During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.

Supernova remnants (SNRs) are believed to be the main sources of Galactic cosmic rays. Molecular clouds associated with SNRs can produce gamma-ray emission through the interaction of accelerated particles with the concentrated gas. The middle aged SNR W28, for its associated system of dense molecular clouds, provides an excellent opportunity to test this hypothesis. We present the AGILE/GRID observations of SNR W28, and compare them with observations at other wavelengths (TeV and 12CO J=1-->0 molecular line emission). The gamma-ray flux detected by AGILE from the dominant source associated with W28 is (14 +- 5) 10^-8 ph cm^-2 s^-1 for E > 400 MeV. This source is positionally well correlated with the TeV emission observed by the HESS telescope. The local variations of the GeV to TeV flux ratio suggest a difference between the CR spectra of the north-west and south molecular cloud complexes. A model based on a hadronic-induced interaction and diffusion with two molecular clouds at different distances from the W28 shell can explain both the morphological and spectral features observed by AGILE in the MeV-GeV energy range and by the HESS telescope in the TeV energy range. The combined set of AGILE and H.E.S.S. data strongly support a hadronic model for the gamma-ray production in W28.

We present the observation of Cyg X-1 in hard spectral state performed during the AGILE Science Verification Phase and Observing Cycle 1 in hard X-rays (with SuperAGILE) and gamma rays (with the GRID) and lasting for about 160 days with a live time of $\sim 6$ Ms. We investigate the variability of Cyg X-1 in hard X-rays at different timescales, from $\sim 300$ s up to one day, and we apply different tools of timing analysis, such as the autocorrelation function, the first order structure function and the Lomb-Scargle periodogram, to our data (from SuperAGILE) and to the simultaneous data in soft X-rays (from RXTE/ASM). We conclude our investigation with a search for emission in the energy range above 100 MeV with the maximum likelihood technique. In the hard X-ray band the flux of Cyg X-1 shows its typical erratic fluctuations at all timescales with variations of about a factor of two that do not affect significantly the shape of the energy spectrum. From the first order structure function we find that the X-ray emission of Cyg X-1 is characterized by \textit{antipersistence}, indication of a negative feedback mechanism at work. In the gamma ray data a statistically significant point-like source at the position of Cyg X-1 is not found and the upper limit on the flux is $\mathrm{5 \times 10^{-8} \; ph \; cm^{-2} \; s^{-1}}$, over the whole observation (160 days). Finally we compare our upper limit in gamma rays with the expectation of various models of the Cyg X-1 emission, of both hadronic and leptonic origin, in the GeV -- TeV band. The time history of Cyg X-1 in the hard X-ray band over 13 months (not continuous) is shown. Different tools of analysis do not provide fully converging results of the characteristic timescales in the system, suggesting that the timescales found in the structure function are not intrinsic to the physics of the source.

Cygnus X-1 is the archetypal black hole (BH) binary system in our Galaxy. We report the main results of an extensive search for transient gamma-ray emission from Cygnus X-1 carried out in the energy range 100 MeV - 3 GeV by the AGILE satellite, during the period 2007 July - 2009 October. The total exposure time is about 300 days, during which the source was in the "hard" X-ray spectral state. We divided the observing intervals in 2 or 4 week periods, and searched for transient and persistent emission. We report an episode of significant transient gamma-ray emission detected on 2009, October 16 in a position compatible with Cygnus X-1 optical position. This episode, occurred during a hard spectral state of Cygnus X-1, shows that a 1-2 day time variable emission above 100 MeV can be produced during hard spectral states, having important theoretical implications for current Comptonization models for Cygnus X-1 and other microquasars. Except for this one short timescale episode, no significant gamma-ray emission was detected by AGILE. By integrating all available data we obtain a 2$\sigma$ upper limit for the total integrated flux of $F_{\gamma,U.L.} = 3 \times 10^{-8} \rm ph cm^{-2} s^{-1}$ in the energy range 100 MeV - 3 GeV. We then clearly establish the existence of a spectral cutoff in the energy range 1-100 MeV that applies to the typical hard state outside the flaring period and that confirms the historically known spectral cutoff above 1 MeV.

We report on 18 months of multiwavelength observations of the blazar 3C 454.3 (Crazy Diamond) carried out in July 2007-January 2009. We show the results of the AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008-January 2009. During the May 2008-January 2009 period, the source average flux was highly variable, from an average gamma-ray flux F(E>100MeV) > 200E-8 ph/cm2/s in May-June 2008, to F(E>100MeV)~80E-8 ph/cm2/s in October 2008-January 2009. The average gamma-ray spectrum between 100 MeV and 1 GeV can be fit by a simple power law (Gamma_GRID ~ 2.0 to 2.2). Only 3-sigma upper limits can be derived in the 20-60 keV energy band with Super-AGILE. During July-August 2007 and May-June 2008, RXTE measured a flux of F(3-20 keV)= 8.4E-11 erg/cm2/s, and F(3-20 keV)=4.5E-11 erg/cm2/s, respectively and a constant photon index Gamma_PCA=1.65. Swift/XRT observations were carried out during all AGILE campaigns, obtaining a F(2-10 keV)=(0.9-7.5)E-11 erg/cm2/s and a photon index Gamma_XRT=1.33-2.04. BAT measured an average flux of ~5 mCrab. GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period from the radio to the optical. A correlation analysis between the optical and the gamma-ray fluxes shows a time lag of tau=-0.4 days. An analysis of 15 GHz and 43 GHz VLBI core radio flux observations shows an increasing trend of the core radio flux, anti- correlated with the higher frequency data. The modeling SEDs, and the behavior of the long-term light curves in different energy bands, allow us to compare the jet properties during different emission states, and to study the geometrical properties of the jet on a time-span longer than one year.

The Supernova Remnant (SNR) IC 443 is an intermediate-age remnant well known for its radio, optical, X-ray and gamma-ray energy emissions. In this Letter we study the gamma-ray emission above 100 MeV from IC 443 as obtained by the AGILE satellite. A distinct pattern of diffuse emission in the energy range 100 MeV-3 GeV is detected across the SNR with its prominent maximum (source "A") localized in the Northeastern shell with a flux F = (47 \pm 10) 10^{-8} photons cm^{-2} s^{-1} above 100 MeV. This location is the site of the strongest shock interaction between the SNR blast wave and the dense circumstellar medium. Source "A" is not coincident with the TeV source located 0.4 degree away and associated with a dense molecular cloud complex in the SNR central region. From our observations, and from the lack of detectable diffuse TeV emission from its Northeastern rim, we demonstrate that electrons cannot be the main emitters of gamma-rays in the range 0.1-10 GeV at the site of the strongest SNR shock. The intensity, spectral characteristics, and location of the most prominent gamma-ray emission together with the absence of co-spatial detectable TeV emission are consistent only with a hadronic model of cosmic-ray acceleration in the SNR. A high-density molecular cloud (cloud "E") provides a remarkable "target" for nucleonic interactions of accelerated hadrons: our results show enhanced gamma-ray production near the molecular cloud/shocked shell interaction site. IC 443 provides the first unambiguous evidence of cosmic-ray acceleration by SNRs.

Pulsars are known to power winds of relativistic particles that can produce bright nebulae by interacting with the surrounding medium. These pulsar wind nebulae (PWNe) are observed in the radio, optical, x-rays and, in some cases, also at TeV energies, but the lack of information in the gamma-ray band prevents from drawing a comprehensive multiwavelength picture of their phenomenology and emission mechanisms. Using data from the AGILE satellite, we detected the Vela pulsar wind nebula in the energy range from 100 MeV to 3 GeV. This result constrains the particle population responsible for the GeV emission, probing multivavelength PWN models, and establishes a class of gamma-ray emitters that could account for a fraction of the unidentified Galactic gamma-ray sources.

We present the long-term monitoring of the High Mass X-ray Binary GX 301-2 performed with the SuperAGILE instrument on-board the AGILE mission. The source was monitored in the 20-60 keV energy band during the first year of the mission from 2007 July 17 to 2008 August 31, covering about one whole orbital period and three more pre-periastron passages for a total net observation time of about 3.7 Ms. The SuperAGILE dataset represents one of the most continuous and complete monitoring at hard X-ray energies of the 41.5 day long binary period available to date. The source behavior was characterized at all orbital phases in terms of hard X-ray flux, spectral hardness, spin period history, pulsed fraction and pulse shape profile. We also complemented the SuperAGILE observations with the soft X-ray data of the RossiXTE/ASM. Our analysis shows a clear orbital modulation of the spectral hardness, with peaks in correspondence with the pre-periastron flare and near phase 0.25. The hardness peaks we found could be related with the wind-plus-stream accretion model proposed in order to explain the orbital light curve modulation of GX 301-2. Timing analysis of the pulsar spin period shows that the secular trend of the about 680 s pulse period is consistent with the previous observations, although there is evidence of a slight decrease in the spin-down rate. The analysis of the hard X-ray pulsed emission also showed a variable pulse shape profile as a function of the orbital phase, with substructures detected near the passage at the periastron, and a clear modulation of the pulsed fraction, which appears in turn strongly anti-correlated with the source intensity.

The study of relativistic particle acceleration is a major topic of high-energy astrophysics. It is well known that massive black holes in active galaxies can release a substantial fraction of their accretion power into energetic particles, producing gamma-rays and relativistic jets. Galactic microquasars (hosting a compact star of 1-10 solar masses which accretes matter from a binary companion) also produce relativistic jets. However, no direct evidence of particle acceleration above GeV energies has ever been obtained in microquasar ejections, leaving open the issue of the occurrence and timing of extreme matter energization during jet formation. Here we report the detection of transient gamma-ray emission above 100 MeV from the microquasar Cygnus X-3, an exceptional X-ray binary which sporadically produces powerful radio jets. Four gamma-ray flares (each lasting 1-2 days) were detected by the AGILE satellite simultaneously with special spectral states of Cygnus X-3 during the period mid-2007/mid-2009. Our observations show that very efficient particle acceleration and gamma-ray propagation out of the inner disk of a microquasar usually occur a few days before major relativistic jet ejections. Flaring particle energies can be thousands of times larger than previously detected maximum values (with Lorentz factors of 105 and 102 for electrons and protons, respectively). We show that the transitional nature of gamma-ray flares and particle acceleration above GeV energies in Cygnus X-3 is clearly linked to special radio/X-ray states preceding strong radio flares. Thus gamma-rays provide unique insight into the nature of physical processes in microquasars.

The cosmic ray interaction event generator Sibyll is widely used in extensive air shower simulations. We describe in detail the properties of Sibyll 2.1 and the differences with the original version 1.7. The major structural improvements are the possibility to have multiple soft interactions, introduction of new parton density functions, and an improved treatment of diffraction. Sibyll 2.1 gives better agreement with fixed target and collider data, especially for the inelastic cross sections and multiplicities of secondary particles. Shortcomings and suggestions for future improvements are also discussed.

SuperAGILE is the hard X-ray monitor of the AGILE gamma ray mission, in orbit since 23$^{rd}$ April 2007. It is an imaging experiment based on a set of four independent silicon strip detectors, equipped with one-dimensional coded masks, operating in the nominal energy range 18-60 keV. The main goal of SuperAGILE is the observation of cosmic sources simultaneously with the main gamma-ray AGILE experiment, the Gamma Ray Imaging Detector (GRID). Given its $\sim$steradian-wide field of view and its $\sim$15 mCrab day-sensitivity, SuperAGILE is also well suited for the long-term monitoring of Galactic compact objects and the detection of bright transients. The SuperAGILE detector properties and design allow for a 6 arcmin angular resolution in each of the two independent orthogonal projections of the celestial coordinates. Photon by photon data are continuously available by the experiment telemetry, and are used to derive images and fluxes of individual sources, with integration times depending on the source intensity and position in the field of view. In this paper we report on the main scientific results achieved by SuperAGILE over its first two years in orbit, until April 2009.

We report on the second AGILE multiwavelength campaign of the blazar 3C 454.3 during the first half of December 2007. This campaign involved AGILE, Spitzer, Swift,Suzaku,the WEBT consortium,the REM and MITSuME telescopes,offering a broad band coverage that allowed for a simultaneous sampling of the synchrotron and inverse Compton (IC) emissions.The 2-week AGILE monitoring was accompanied by radio to optical monitoring by WEBT and REM and by sparse observations in mid-Infrared and soft/hard X-ray energy bands performed by means of Target of Opportunity observations by Spitzer, Swift and Suzaku, respectively.The source was detected with an average flux of~250x10^{-8}ph cm^-2s^-1 above 100 MeV,typical of its flaring states.The simultaneous optical and gamma-ray monitoring allowed us to study the time-lag associated with the variability in the two energy bands, resulting in a possible ~1-day delay of the gamma-ray emission with respect to the optical one. From the simultaneous optical and gamma-ray fast flare detected on December 12, we can constrain the delay between the gamma-ray and optical emissions within 12 hours. Moreover, we obtain three Spectral Energy Distributions (SEDs) with simultaneous data for 2007 December 5, 13, 15, characterized by the widest multifrequency coverage. We found that a model with an external Compton on seed photons by a standard disk and reprocessed by the Broad Line Regions does not describe in a satisfactory way the SEDs of 2007 December 5, 13 and 15. An additional contribution, possibly from the hot corona with T=10^6 K surrounding the jet, is required to account simultaneously for the softness of the synchrotron and the hardness of the inverse Compton emissions during those epochs.

We report results from an intensive multiwavelength campaign on the intermediate-frequency-peaked BL Lacertae object W Com (z=0.102) during a strong outburst of very high energy gamma-ray emission in June 2008. The very high energy gamma-ray signal was detected by VERITAS on 2008 June 7-8 with a flux F(>200 GeV) = (5.7+-0.6)x10^-11 cm-2s-1, about three times brighter than during the discovery of gamma-ray emission from W Com by VERITAS in 2008 March. The initial detection of this flare by VERITAS at energies above 200 GeV was followed by observations in high energy gamma-rays (AGILE, E>100 MeV), and X-rays (Swift and XMM-Newton), and at UV, and ground-based optical and radio monitoring through the GASP-WEBT consortium and other observatories. Here we describe the multiwavelength data and derive the spectral energy distribution (SED) of the source from contemporaneous data taken throughout the flare.

JEM-EUSO, on board of the Japanese Exploration Module of the International Space Station, is being proposed as the first space observatory devoted to UHECR. Its privileged position at 430 km above the Earth surface, combined with a large field of view, innovative optics and a high efficiency focal surface, results in an unprecedented exposure which significantly surpasses that of the largest ground observatories. The large number of events expected above the GZK threshold for photo-pion production by protons will allow the directional identification of individual sources and the determination of their spectra, i.e., doing astronomy and astrophysics through the particle channel. Similar goals can be achieved in the case of light UHECR nuclei. Furthermore, the atmospheric target volume ($\sim 10^{12}$ ton) makes the possibility of neutrino observation a highlight of the mission. Other exploratory objectives include the detection of extreme energy gammas and the study of Galactic magnetic fields as well as global observations of the earth's atmosphere, including clouds, night-glows, plasma discharges, and meteors. In this contribution we will describe the scientific objectives of JEM-EUSO.

Short gamma-ray bursts (GRBs), typically lasting less than 2 s, are a special class of GRBs of great interest. We report the detection by the AGILE satellite of the short GRB 090510 which shows two clearly distinct emission phases: a prompt phase lasting ~ 200 msec and a second phase lasting tens of seconds. The prompt phase is relatively intense in the 0.3-10 MeV range with a spectrum characterized by a large peak/cutoff energy near 3 MeV, in this phase, no significant high-energy gamma-ray emission is detected. At the end of the prompt phase, intense gamma-ray emission above 30 MeV is detected showing a power-law time decay of the flux of the type t^-1.3 and a broad-band spectrum remarkably different from that of the prompt phase. It extends from sub-MeV to hundreds of MeV energies with a photon index alpha ~ 1.5. GRB 090510 provides the first case of a short GRB with delayed gamma-ray emission. We present the timing and spectral data of GRB 090510 and briefly discuss its remarkable properties within the current models of gamma-ray emission of short GRBs.

We report the detection by the AGILE satellite of a rapid gamma-ray flare from the powerful gamma-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been continuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). Moreover, the gamma-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 2008. In the period 1-16 March 2008, AGILE detected gamma-ray emission from PKS 1510-089 at a significance level of 6.2-sigma with an average flux over the entire period of (84 +/- 17) x 10^{-8} photons cm^{-2} s^{-1} for photon energies above 100 MeV. After a predefined satellite re-pointing, between 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-sigma, with an average flux (E > 100 MeV) of (134 +/- 29) x 10^{-8} photons cm^{-2} s^{-1} and a peak level of (281 +/- 68) x 10^{-8} photons cm^{-2} s^{-1} with daily integration. During the observing period January-April 2008, the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at millimetric frequencies. Moreover, in the X-ray band the Swift/XRT observations seem to show an harder-when-brighter behaviour of the source spectrum. The spectral energy distribution of mid-March 2008 is modelled with a homogeneous one-zone synchrotron self Compton emission plus contributions from inverse Compton scattering of external photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of Seyfert-like components, such as the little blue bump and the big blue bump.

The identification of very high energy photons is of great importance for the understanding of the origin of extreme energy cosmic rays (EECR). Several can be the sources of high energy photons at Earth. A guaranteed component is the flux of high energy photons expected as a consequence of the interaction of cosmic rays with the cosmic photon background. Another contribution may be expected as by-product at the acceleration sites of protons and nuclei, although such flux should be strongly suppressed for distant sources. On the other hand, top-down scenarios involving the decay of super heavy relic particles or topological defects, even if not currently favored, have as a characteristic signature an increasingly dominant flux of photons at the highest energies. In this work we study the statistical separation between hadron and photon showers at energies where both, LPM effect and magnetospheric interactions are important for the development of the cascades. We consider a detector with the same orbital characteristics as JEM-EUSO, but disregard trigger and reconstruction efficiencies, in order to define the maximum ideal discrimination power attainable.

Neutrinos are a very promising messenger at tens of EeV and above. They can be produced by several channels, namely as by products of hadronic interactions at the sources, as the main products of the decay of super massive particles and, in a guaranteed way, as the result of the propagation of UHECR through the bath of microwave relic photons. A new era of very large exposure space observatories, of which the JEM-EUSO mission is a prime example, is on the horizon and, with it, it is even larger the possibility of astrophysical neutrino detection at the highest energies. In the present work we use a combination of the PYTHIA interaction code with the CONEX shower simulation package in order to produce fast one-dimensional simulations of neutrino initiated showers in air. We make a detail study of the structure of the corresponding longitudinal profiles, but focus our physical analysis mainly on the development of showers at mid and high altitudes, where they can be an interesting target for space fluorescence observatories.

We present the first catalog of high-confidence gamma-ray sources detected by the AGILE satellite during observations performed from July 9, 2007 to June 30, 2008. Catalogued sources are detected by merging all the available data over the entire time period. AGILE, launched in April 2007, is an ASI mission devoted to gamma-ray observations in the 30 MeV - 50 GeV energy range, with simultaneous X-ray imaging capability in the 18-60 keV band. This catalog is based on Gamma-Ray Imaging Detector (GRID) data for energies greater than 100 MeV. For the first AGILE catalog we adopted a conservative analysis, with a high-quality event filter optimized to select gamma-ray events within the central zone of the instrument Field of View (radius of 40 degrees). This is a significance-limited (4 sigma) catalog, and it is not a complete flux-limited sample due to the non-uniform first year AGILE sky coverage. The catalog includes 47 sources, 21 of which are associated with confirmed or candidate pulsars, 13 with Blazars (7 FSRQ, 4 BL Lacs, 2 unknown type), 2 with HMXRBs, 2 with SNRs, 1 with a colliding-wind binary system, 8 with unidentified sources.

Since its early phases of operation, the AGILE satellite is observing Gamma Ray Bursts (GRBs) over an energy range potentially spanning six orders of magnitude. In the hard X-ray band the SuperAGILE imager provides localization of about one GRB/month plus the detection of 1-2 GRBs per month out of its field of view. The Mini-Calorimeter detects about one GRB/week in the 350keV - 100MeV energy range, plus several other transients at very short time scales. In fact, the on-board MCAL trigger logic, implemented for the first time on time windows as short as 300 microseconds, is particularly suitable for very short bursts detection. The Gamma-Ray Imaging Detector (GRID), sensitive in the 30MeV - 30GeV range, firmly detected three GRBs (GRB080514B, GRB090401B and GRB090510) plus some other candidates at a lower significance level. Moreover, all GRBs localized by other spacecrafts inside the GRID field of view are currently searched for possible detection, and upper limits are provided. In this paper we review the status of the observation of GRBs with AGILE, with particular emphasis on high energy and short time scales detections.

We present the results of extensive observations by the gamma-ray AGILE satellite of the Galactic region hosting the Carina nebula and the remarkable colliding wind binary Eta Carinae (Eta Car) during the period 2007 July to 2009 January. We detect a gamma-ray source (1AGL J1043-5931) consistent with the position of Eta Car. If 1AGL J1043-5931 is associated with the Eta Car system our data provide the long sought first detection above 100 MeV of a colliding wind binary. The average gamma-ray flux above 100 MeV and integrated over the pre-periastron period 2007 July to 2008 October is F = (37 +/- 5) x 10-8 ph cm-2 s-1 corresponding to an average gamma-ray luminosity of L = 3.4 x 10^34 erg s-1 for a distance of 2.3 kpc. We also report a 2-day gamma-ray flaring episode of 1AGL J1043-5931 on 2008 Oct. 11-13 possibly related to a transient acceleration and radiation episode of the strongly variable shock in the system.

Using gamma-ray data collected by the Astrorivelatore Gamma ad Immagini LEggero (AGILE) satellite over a period of almost one year (from 2007 July to 2008 June), we searched for pulsed signals from 35 potentially interesting radio pulsars, ordered according to $F_{\gamma}\propto \sqrt{\dot{E}} d^{-2}$ and for which contemporary or recent radio data were available. AGILE detected three new top-ranking nearby and Vela-like pulsars with good confidence both through timing and spatial analysis. Among the newcomers we find pulsars with very high rotational energy losses, such as the remarkable PSR B1509-58 with a magnetic field in excess of 10^13 Gauss, and PSR J2229+6114 providing a reliable identification for the previously unidentified EGRET source 3EG 2227+6122. Moreover, the powerful millisecond pulsar B1821-24, in the globular cluster M28, is detected during a fraction of the observations. Four other promising gamma-ray pulsar candidates, among which is the notable J2043+2740 with an age in excess of 1 million years, show a possible detection in the timing analysis only and deserve confirmation.

We report the results of a 3-weeks multi-wavelength campaign on the flat spectrum radio quasar 3C 273 carried out with the AGILE gamma-ray mission, covering the 30 MeV -50 GeV and 18-60 keV, the REM observatory (covering the near-IR and optical), Swift (near-UV/Optical, 0.2-10 keV and 15-50 keV), INTEGRAL (3 - 200 keV) and Rossi XTE (2-12 keV). This is the first observational campaign including gamma-ray data, after the last EGRET observations, more than 8 years ago. Our study was carried out using simultaneous light curves of the source flux from all the involved instruments, in the different energy ranges, in search for correlated variability. Then a time-resolved SED was used for a detailed physical modelling of the emission mechanisms. The source was detected in gamma-rays only in the second week of our campaign. We found indication of a possible anti-correlation between the emission at gamma-rays and at soft and hard X-rays, supported by the complete set of instruments. Instead, optical data do not show short term variability. Only in two EGRET observations (in 1993 and 1997) 3C 273 showed intra-observation variability in gamma-rays. The SED with almost simultaneous data, partially covers the regions of the synchrotron emission, the big blue bump, and the inverse-Compton. The soft X-ray emission is consistent with combined SSC and EC mechanisms, while hard X and gamma-ray emissions are compatible with EC from thermal photons of the disk. The time evolution of the SED is well interpreted and modelled in terms of an acceleration episode of the electrons population, leading to a shift in the IC peak towards higher energies.

Context. We report the detection by the AGILE satellite of an intense gamma-ray flare from the gamma-ray source 3EG J1255-0549, associated to the Flat Spectrum Radio Quasar 3C 279, during the AGILE pointings towards the Virgo Region on 2007 July 9-13. Aims. The simultaneous optical, X-ray and gamma-ray covering allows us to study the spectral energy distribution (SED) and the theoretical models relative to the flaring episode of mid-July. Methods. AGILE observed the source during its Science Performance Verification Phase with its two co-aligned imagers: the Gamma- Ray Imaging Detector (GRID) and the hard X-ray imager (Super-AGILE) sensitive in the 30 MeV - 50 GeV and 18 - 60 keV respectively. During the AGILE observation the source was monitored simultaneously in optical band by the REM telescope and in the X-ray band by the Swift satellite through 4 ToO observations. Results. During 2007 July 9-13 July 2007, AGILE-GRID detected gamma-ray emission from 3C 279, with the source at ~2 deg from the center of the Field of View, with an average flux of (210+-38) 10^-8 ph cm^-2 s^-1 for energy above 100 MeV. No emission was detected by Super-AGILE, with a 3-sigma upper limit of 10 mCrab. During the observation lasted about 4 days no significative gamma-ray flux variation was observed. Conclusions. The Spectral Energy Distribution is modelled with a homogeneous one-zone Synchrotron Self Compton emission plus the contributions by external Compton scattering of direct disk radiation and, to a lesser extent, by external Compton scattering of photons from the Broad Line Region.

AGILE is a small gamma-ray astronomy satellite mission of the Italian Space Agency dedicated to high-energy astrophysics launched in 2007 April. Its 1 microsecond absolute time tagging capability coupled with a good sensitivity in the 30 MeV-30 GeV range, with simultaneous X-ray monitoring in the 18-60 keV band, makes it perfectly suited for the study of gamma-ray pulsars following up on the CGRO/EGRET heritage. In this paper we present the first AGILE timing results on the known gamma-ray pulsars Vela, Crab, Geminga and B 1706-44. The data were collected from 2007 July to 2008 April, exploiting the mission Science Verification Phase, the Instrument Timing Calibration and the early Observing Pointing Program. Thanks to its large field of view, AGILE collected a large number of gamma-ray photons from these pulsars (about 10,000 pulsed counts for Vela) in only few months of observations. The coupling of AGILE timing capabilities, simultaneous radio/X-ray monitoring and new tools aimed at precise photon phasing, exploiting also timing noise correction, unveiled new interesting features at sub-millisecond level in the pulsars' high-energy light-curves.

We report the detection by the AGILE (Astro-rivelatore Gamma a Immagini LEggero) satellite of an intense gamma-ray flare from the source AGL J1511-0909, associated with the powerful quasar PKS 1510-089, during ten days of observations from 23 August to 1 September 2007. During the observation period, the source was in optical decrease following a flaring event monitored by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). The simultaneous gamma-ray, optical, and radio coverage allows us to study the spectral energy distribution and the theoretical models based on the synchrotron and inverse Compton (IC) emission mechanisms. AGILE observed the source with its two co-aligned imagers, the Gamma-Ray Imaging Detector and the hard X-ray imager Super-AGILE sensitive in the 30 MeV - 50 GeV and 18 - 60 keV bands, respectively. Between 23 and 27 August 2007, AGILE detected gamma-ray emission from PKS 1510-089 when this source was located about 50 degrees off-axis, with an average flux of (270 +/- 65) x 10^{-8} photons cm^{-2} s^{-1} for photon energy above 100 MeV. In the following period, 28 August - 1 September, after a satellite re-pointing, AGILE detected the source at about 35 degrees off-axis, with an average flux (E > 100 MeV) of (195 +/- 30) x 10^{-8} photons cm^{-2} s^{-1}. No emission was detected by Super-AGILE, with a 3-sigma upper limit of 45 mCrab in 200 ksec. The spectral energy distribution is modelled with a homogeneous one-zone synchrotron self Compton (SSC) emission plus contributions by external photons: the SSC emission contributes primarily to the X-ray band, whereas the contribution of the IC from the external disc and the broad line region match the hard gamma-ray spectrum observed.

The Mini-Calorimeter (MCAL) instrument on-board the AGILE satellite is a non-imaging gamma-ray scintillation detector sensitive in the 300keV-100MeV energy range with a total on-axis geometrical area of 1400cm^2. Gamma-Ray Bursts (GRBs) are one of the main scientific targets of the AGILE mission and the MCAL design as an independent self-triggering detector makes it a valuable all-sky monitor for GRBs. Furthermore MCAL is one of the very few operative instruments with microsecond timing capabilities in the MeV range. In this paper the results of GRB detections with MCAL after one year of operation in space are presented and discussed. A flexible trigger logic implemented in the AGILE payload data-handling unit allows the on-board detection of GRBs. For triggered events, energy and timing information are sent to telemetry on a photon-by-photon basis, so that energy and time binning are limited by counting statistics only. When the trigger logic is not active, GRBs can be detected offline in ratemeter data, although with worse energy and time resolution. Between the end of June 2007 and June 2008 MCAL detected 51 GRBs, with a detection rate of about 1 GRB/week, plus several other events at a few milliseconds timescales. Since February 2008 the on-board trigger logic has been fully active. Comparison of MCAL detected events and data provided by other space instruments confirms the sensitivity and effective area estimations. MCAL also joined the 3rd Inter-Planetary Network, to contribute to GRB localization by means of triangulation.

[Abridged] We report on a multiwavelength observation of the blazar 3C 454.3 (which we dubbed "crazy diamond") carried out on November 2007 by means of the astrophysical satellites AGILE, INTEGRAL, Swift, the WEBT Consortium, and the optical-NIR telescope REM. 3C 454.3 is detected at a $\sim 19-\sigma$ level during the 3-week observing period, with an average flux above 100 MeV of $F_{\rm E>100MeV} = (170 \pm 13) \times 10^{-8}$ \phcmsec. The gamma-ray spectrum can be fit with a single power-law with photon index $\Gamma_{\rm GRID} = 1.73 \pm 0.16$ between 100 MeV and 1 GeV. We detect significant day-by-day variability of the gamma-ray emission during our observations, and we can exclude that the fluxes are constant at the 99.6% ($\sim 2.9 \sigma$) level. The source was detected typically around 40 degrees off-axis, and it was substantially off--axis in the field of view of the AGILE hard X-ray imager. However, a 5-day long ToO observation by INTEGRAL detected 3C 454.3 at an average flux of about $F_{\rm 20-200 keV} = 1.49 \times 10^{-3}$ \phcmsec with an average photon index of $\Gamma_{\rm IBIS} = 1.75 \pm 0.24$ between 20--200 keV. Swift also detected 3C 454.3 with a flux in the 0.3--10 keV energy band in the range $(1.23-1.40) \times 10^{-2}$ \phcmsec{} and a photon index in the range $\Gamma_{\rm XRT} = 1.56-1.73$. In the optical band, both WEBT and REM show an extremely variable behavior in the $R$ band. A correlation analysis based on the entire data set is consistent with no time-lags between the gamma-ray and the optical flux variations. Our simultaneous multifrequency observations strongly indicate that the dominant emission mechanism between 30 MeV and 30 GeV is dominated by inverse Compton scattering of relativistic electrons in the jet on the external photons from the broad line region.

GRB 080514B is the first gamma ray burst (GRB), since the time of EGRET, for which individual photons of energy above several tens of MeV have been detected with a pair-conversion tracker telescope. This burst was discovered with the Italian AGILE gamma-ray satellite. The GRB was localized with a cooperation by AGILE and the interplanetary network (IPN). The gamma-ray imager (GRID) estimate of the position, obtained before the SuperAGILE-IPN localization, is found to be consistent with the burst position. The hard X-ray emission observed by SuperAGILE lasted about 7 s, while there is evidence that the emission above 30 MeV extends for a longer duration (at least ~13 s). Similar behavior was seen in the past from a few other GRBs observed with EGRET. However, the latter measurements were affected, during the brightest phases, by instrumental dead time effects, resulting in only lower limits to the burst intensity. Thanks to the small dead time of the AGILE/GRID we could assess that in the case of GRB 080514B the gamma-ray to X-ray flux ratio changes significantly between the prompt and extended emission phase.

The concept of "age" as a parameter for the description of the state of development of high energy showers in the atmosphere has been in use in cosmic ray studies for several decades. In this work we briefly discuss how this concept, originally introduced to describe the average behavior of electromagnetic cascades, can be fruitfully applied to describe individual showers generated by primary particles of different nature, including protons, nuclei and neutrinos. Showers with the same age share three different important properties: (i) their electron size has the same fractional rate of change with increasing depth, (ii) the bulk of the electrons and photons in the shower (excluding high energy particles) have energy spectra with shapes and relative normalization uniquely determined by the age parameter, (iii) the electrons and photons in the shower have also the same angular and lateral distributions sufficiently far from the shower axis. In this work we discuss how the properties associated with the shower age can be understood with simple arguments, and how the shapes of the electron and photon spectra and the relative normalization that correspond to a certain age can be calculated analytically.

We report the gamma-ray activity from the intermediate BL Lac S5 0716+714 during observations acquired by the AGILE satellite in September and October 2007. These detections of activity were contemporaneous with a period of intense optical activity, which was monitored by GASP-WEBT. This simultaneous optical and gamma-ray coverage allows us to study in detail the light curves, time lags, gamma-ray photon spectrum, and Spectral Energy Distributions (SEDs) during different states of activity. AGILE observed the source with its two co-aligned imagers, the Gamma-Ray Imaging Detector (GRID) and the hard X-ray imager (Super-AGILE), which are sensitive to the 30 MeV-50 GeV and 18-60 keV energy ranges, respectively. Observations were completed in two different periods, the first between 2007 September 4-23, and the second between 2007 October 24-November 1. Over the period 2007 September 7-12, AGILE detected gamma-ray emission from the source at a significance level of 9.6-sigma with an average flux (E>100 MeV) of (97 +/- 15) x 10^{-8} photons cm^{-2} s^{-1}, which increased by a factor of at least four within three days. No emission was detected by Super-AGILE for the energy range 18-60 keV to a 3-sigma upper limit of 10 mCrab in 335 ksec. In October 2007, AGILE repointed toward S5 0716+714 following an intense optical flare, measuring an average flux of (47 +/- 11) x 10^{-8} photons cm^{-2} s^{-1} at a significance level of 6.0-sigma. The gamma-ray flux of S5 0716+714 detected by AGILE is the highest ever detected for this blazar and one of the most intense gamma-ray fluxes detected from a BL Lac object. The SED of mid-September appears to be consistent with the synchrotron self-Compton (SSC) emission model, but only by including two SSC components of different variabilities.

We present the AGILE gamma-ray observations of the field containing the puzzling gamma-ray source 3EG J1835+5918. This source is one of the most remarkable unidentified EGRET sources. An unprecedentedly long AGILE monitoring of this source yields important information on the positional error box, flux evolution, and spectrum. 3EG J1835+5918 has been in the AGILE field of view several times in 2007 and 2008 for a total observing time of 138 days from 2007 Sept 04 to 2008 June 30 encompassing several weeks of continuous coverage. With an exposure time approximately twice that of EGRET, AGILE confirms the existence of a prominent gamma-ray source (AGL J1836+5926) at a position consistent with that of EGRET, although with a remarkably lower average flux value for photon energies greater than 100 MeV. A 5-day bin temporal analysis of the whole data set of AGL J1836+5926 shows some evidence for variability of the gamma-ray flux. The source spectrum between 100 MeV and 1 GeV can be fitted with a power law with photon index in the range 1.6-1.7, fully consistent with the EGRET value. The faint X-ray source RX J1836.2+5925 that has been proposed as a possible counterpart of 3EG J1835+5918 is well within the AGILE error box. Future continuous monitoring (both by AGILE and GLAST) is needed to confirm the gamma-ray flux variability and to unveil the source origin, a subject that is currently being pursued through a multiwavelength search for counterparts.

This contribution discusses some of the main problems in high energy astrophysics, and the perspectives to solve them using different types of "messengers": cosmic rays, photons and neutrinos

The study of extragalactic sources of high energy radiation via the direct measurement of the proton and neutrino fluxes that they are likely to emit is one of the main goals for the future observations of the recently developed air showers detectors and neutrino telescopes. In this work we discuss the relation between the inclusive proton and neutrino signals from the ensemble of all sources in the universe, and the resolved signals from the closest and brightest objects. We also compare the sensitivities of proton and neutrino telescopes and comment on the relation between these two new astronomies.

We report the first blazar detection by the AGILE satellite. AGILE detected 3C 454.3 during a period of strongly enhanced optical emission in July 2007. AGILE observed the source with a dedicated repointing during the period 2007 July 24-30 with its two co-aligned imagers, the Gamma-Ray Imaging Detector and the hard X-ray imager Super-AGILE sensitive in the 30 MeV-50 GeV and 18-60 keV, respectively. Over the entire period, AGILE detected gamma-ray emission from 3C 454.3 at a significance level of 13.8-$\sigma$ with an average flux (E$>$100 MeV) of $(280 \pm 40) \times 10^{-8}$ photons cm$^{-2}$ s$^{-1}$. The gamma-ray flux appears to be variable towards the end of the observation. No emission was detected by Super-AGILE in the energy range 20-60 keV, with a 3-$\sigma$ upper limit of $2.3 \times 10^{-3}$ photons cm$^{-2}$ s$^{-1}$. The gamma-ray flux level of 3C 454.3 detected by AGILE is the highest ever detected for this quasar and among the most intense gamma-ray fluxes ever detected from Flat Spectrum Radio Quasars.

GRB 070724B is the first Gamma Ray Burst localized by SuperAGILE, the hard X-ray monitor aboard the AGILE satellite. The coordinates of the event were published $\sim 19$ hours after the trigger. The Swift X-Ray Telescope pointed at the SuperAGILE location and detected the X-ray afterglow inside the SuperAGILE error circle. The AGILE gamma-ray Tracker and Minicalorimeter did not detect any significant gamma ray emission associated with GRB 070724B in the MeV and GeV range, neither prompt nor delayed. Searches of the optical afterglow were performed by the Swift UVOT and the Palomar automated 60-inch telescopes without any significant detection. Similarly the Very Large Array did not detect a radio afterglow. This is the first GRB event with a firm upper limit in the 100 MeV -- 30 GeV energy range, associated with an X-ray afterglow.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

The measurement of the flavor composition of the neutrino fluxes from astrophysical sources has been proposed as a method to study not only the nature of their emission mechanisms, but also the neutrino fundamental properties. It is however problematic to reconcile these two goals, since a sufficiently accurate understanding of the neutrino fluxes at the source is needed to extract information about the physics of neutrino propagation. In this work we discuss critically the expectations for the flavor composition and energy spectrum from different types of astrophysical sources, and comment on the theoretical uncertainties connected to our limited knowledge of their structure.

This work discusses the perspectives to observe fluxes of high energy astrophysical neutrinos with the planned km3 telescopes. On the basis of the observations of GeV and TeV gamma-rays, and of ultra high energy cosmic rays, it is possible to construct well motivated predictions that indicate that the discovery of such fluxes is probable. However the range of these predictions is broad, and the very important opening of the ``neutrino window'' on the high energy universe is not guaranteed with the current design of the detectors. The problem of enlarging the detector acceptance using the same (water/ice Cherenkov) or alternative (acoustic/radio) techniques is therefore of central importance.

A Monte-Carlo calculation of the atmospheric neutrino fluxes [1,2] has been extended to take account of the three-dimensional (3D) nature of the problem, including the bending of secondary particles in the geomagnetic field. Emphasis has been placed on minimizing the approximations when introducing the 3D considerations. In this paper, we describe the techniques used and quantify the effects of the small approximations which remain. We compare 3D and 1D calculations using the same physics input in order to evaluate the conditions under which the 3D calculation is required and when the considerably simpler 1D calculation is adequate. We find that the 1D and 3D results are essentially identical for neutrino energy greater than 5 GeV except for small effects in the azimuthal distributions due to bending of the secondary muon by the geomagnetic field during their propagation in the atmosphere.

Research in cosmic rays is now nearly a century old, but most of the fundamental questions in this field remain unanswered, on the other hand the perspectives of future studies in the next decade are very bright. New detectors will provide higher quality data in the entire energy range from 10^8 to 10^{20} eV (or more if particles of higher energy have non negligible fluxes), moreover cosmic ray astrophysics must now be considered, together with gamma, neutrino and gravitational wave astronomy, as one of the subfields of high energy astrophysics, and using information from these four "messengers" there is the potential of a detailed understanding of the origin of the high energy radiation in the universe. High energy cosmic rays are measured indirectly observing the showers they generate in the atmosphere, and a correct and detailed interpretation of these measurements will require an improved understanding of the properties of hadronic interactions. The new collider experiments, and in particular the LHC project at CERN offer the unique possibility to perform measurements of great value for cosmic ray astrophysics. It is of great importance for cosmic research that this possibility is fully exploited with the appropriate instrumentation and analysis.

New measurements of the cosmic ray fluxes (p, e+, e- and Helium) performed by the Alpha Magnetic Spectrometer (AMS) during a ten days flight of space shuttle have revealed the existence of significant fluxes of particles below the geomagnetic cutoff. These fluxes exhibit a number of remarkable properties, such as a He-3/He-4 ratio of order ~10, an e+/e- ratio of order ~4 and production from well defined regions of the Earth that are distinct for positively and negatively charged particles. In this work we show that the natural hypothesis, that these subcutoff particles are generated as secondary products of primary cosmic ray interactions in the atmosphere can reproduce all the observed properties. We also discuss the implications of the subcutoff fluxes for the estimate of the atmospheric neutrino fluxes, and find that they represent a negligibly small correction. On the other hand the AMS results give important confirmations about the assumption of isotropy for the interplanetary cosmic ray fluxes also on large angular scales, and on the validity of the geomagnetic effects that are important elements for the prediction of the atmospheric neutrino fluxes.

The data on atmospheric neutrinos can be explained assuming the existence of oscillations between muon neutrinos and light sterile neutrinos with mixing close to maximal, and Delta_m^2 approximately equal to 3 E-3 (eV^2). This interpretation of the data is in potential conflict with the successes of big bang nucleosynthesis (BBN), since oscillations can result in a too large contribution of the sterile state to the energy density of the universe at the epoch of nucleosynthesis. The possibility to evade these cosmological constraints has been recently the object of some controversy. In this work we rediscuss this problem and find that the inclusion of a small mixing of the sterile state with tau neutrino can result in the generation of a large lepton asymmetry that strongly suppress the muon-sterile neutrino oscillations eliminating the possible conflict with BBN bounds. In this scheme the mass of the tau neutrino must be larger than few eV's and is compatible with cosmological bounds. Our calculation is performed using a Pauli-Boltzmann method. In this approach it is also possible to develop analytic calculations that allow physical insight in the processes considered and give support to the numerical results.

We reconsider neutrino decay as an explanation for atmospheric neutrino observations. We show that if the mass-difference relevant to the two mixed states \nu_\mu and \nu_\tau is very small (< 10^{-4} eV^2), then a very good fit to the observations can be obtained with decay of a component of \nu_\mu to a sterile neutrino and a Majoron. We discuss how the K2K and MINOS long-baseline experiments can distinguish the decay and oscillation scenarios.

Geomagnetic effects distort the zenith angle distribution of sub--GeV and few--GeV atmospheric neutrinos, breaking the up--down symmetry that would be present in the absence of neutrino oscillations and without a geomagnetic field. The geomagnetic effects also produce a characteristic azimuthal dependence of the $\nu$--fluxes, related to the well known east--west effect, that should be detectable in neutrino experiments of sufficiently large mass. We discuss these effects quantitatively. Because the azimuthal dependence is in first order independent of any oscillation effect, it is a useful diagnostic tool for studying possible systematic effects in the search for neutrino oscillations.