We report on the gamma-ray activity of the high-synchrotron-peaked BL Lacertae object Mrk 421 during the first 1.5 years of Fermi operation, from 2008 August 5 to 2010 March 12. We find that the Large Area Telescope (LAT) gamma-ray spectrum above 0.3 GeV can be well-described by a power-law function with photon index Gamma=1.78 +/- 0.02 and average photon flux F(>0.3 GeV)=(7.23 +/- 0.16) x 10^{-8} ph cm^{-2} s^{-1}. Over this time period, the Fermi-LAT spectrum above 0.3 GeV was evaluated on 7-day-long time intervals, showing significant variations in the photon flux (up to a factor ~3 from the minimum to the maximum flux), but mild spectral variations. The variability amplitude at X-ray frequencies measured by RXTE/ASM and Swift/BAT is substantially larger than that in gamma-rays measured by Fermi-LAT, and these two energy ranges are not significantly correlated. We also present the first results from the 4.5-month-long multifrequency campaign on Mrk 421, which included the VLBA, Swift, RXTE, MAGIC, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign (2009 January 19 to 2009 June 1). During this campaign, Mrk 421 showed a low activity at all wavebands. The extensive multi-instrument (radio to TeV) data set provides an unprecedented, complete look at the quiescent spectral energy distribution (SED) for this source. The broad band SED was reproduced with a leptonic (one-zone Synchrotron Self-Compton) and a hadronic model (Synchrotron Proton Blazar). Both frameworks are able to describe the average SED reasonably well, implying comparable jet powers but very different characteristics for the blazar emission site.

We present the light curves and spectral data of two exceptionally luminous gamma-ray outburts observed by the Large Area Telescope (LAT) experiment on board Fermi Gamma-ray Space Telescope from 3C 273 in September 2009. During these flares, having a duration of a few days, the source reached its highest gamma-ray flux ever measured. This allowed us to study in some details their spectral and temporal structures. The rise and decay are asymmetric on timescales of 6 hours, and the spectral index was significantly harder during the flares than during the preceding 11 months. We also found that short, very intense flares put out the same time-integrated energy as long, less intense flares like that observed in August 2009.

We report on the gamma-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average LAT gamma-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 +/- 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of 2), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 +/- 0.14, and the softest one is 2.51 +/- 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3GeV. In this paper, we also present the first results from the 4.5-month-long multifrequency campaign (2009 March 15 - August 1) on Mrk 501, which included the VLBA, Swift, RXTE, MAGIC and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton model. In the framework of this model, we find that the dominant emission region is characterized by a size <~ 0.1 pc (comparable within a factor of few to the size of the partially-resolved VLBA core at 15-43 GHz), and that the total jet power (~10^{44} erg s^{-1}) constitutes only a small fraction (~10^{-3}) of the Eddington luminosity. The energy distribution of the freshly-accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3GeV-10TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks.

The BL Lacertae object 3C 66A was detected in a flaring state by the Fermi Large Area Telescope (LAT) and VERITAS in 2008 October. In addition to these gamma-ray observations, F-GAMMA, GASP-WEBT, PAIRITEL, MDM, ATOM, Swift, and Chandra provided radio to X-ray coverage. The available light curves show variability and, in particular, correlated flares are observed in the optical and Fermi-LAT gamma-ray band. The resulting spectral energy distribution can be well fit using standard leptonic models with and without an external radiation field for inverse-Compton scattering. It is found, however, that only the model with an external radiation field can accommodate the intra-night variability observed at optical wavelengths.

The sensitivity of a search for sources of TeV neutrinos can be improved by grouping potential sources together into generic classes in a procedure that is known as source stacking. In this paper, we define catalogs of Active Galactic Nuclei (AGN) and use them to perform a source stacking analysis. The grouping of AGN into classes is done in two steps: first, AGN classes are defined, then, sources to be stacked are selected assuming that a potential neutrino flux is linearly correlated with the photon luminosity in a certain energy band (radio, IR, optical, keV, GeV, TeV). Lacking any secure detailed knowledge on neutrino production in AGN, this correlation is motivated by hadronic AGN models, as briefly reviewed in this paper. The source stacking search for neutrinos from generic AGN classes is illustrated using the data collected by the AMANDA-II high energy neutrino detector during the year 2000. No significant excess for any of the suggested groups was found.

The extreme synchrotron BL Lac object H2356-309, located at a redshift of z = 0.165, was observed from June to December 2004 with a total exposure of approx. 40 h live-time with the H.E.S.S. (High Energy Stereoscopic System) array of atmospheric-Cherenkov telescopes (ACTs). Analysis of this data set yields, for the first time, a strong excess of 453 gamma-rays (10 standard deviations above background) from H2356-309, corresponding to an observed integral flux above 200 GeV of I(>200GeV) = (4.1+-0.5) 10^12 cm^-2.s^-1 (statistical error only). The differential energy spectrum of the source between 200 GeV and 1.3 TeV is well-described by a power law with a normalisation (at 1 TeV) of N_0 = (3.00 +- 0.80_stat +- 0.31_sys) 10^-13 cm^-2.s^-1.TeV^-1 and a photon index of Gamma = 3.09 +- 0.24_stat +- 0.10_sys. H2356-309 is one of the most distant BL Lac objects detected at very-high-energy gamma-rays so far. Results from simultaneous observations from ROTSE-III (optical), RXTE (X-rays) and NRT (radio) are also included and used together with the H.E.S.S. data to constrain a single-zone homogeneous synchrotron self-Compton (SSC) model. This model provides an adequate fit to the H.E.S.S. data when using a reasonable set of model parameters.

The Vela supernova remnant (SNR) is a complex region containing a number of sources of non-thermal radiation. The inner section of this SNR, within 2 degrees of the pulsar PSR B0833-45, has been observed by the H.E.S.S. gamma-ray atmospheric Cherenkov detector in 2004 and 2005. A strong signal is seen from an extended region to the south of the pulsar, within an integration region of radius 0.8 deg. around the position (RA = 08h 35m 00s, dec = -45 deg. 36' J2000.0). The excess coincides with a region of hard X-ray emission seen by the ROSAT and ASCA satellites. The observed energy spectrum of the source between 550 GeV and 65 TeV is well fit by a power law function with photon index = 1.45 +/- 0.09(stat) +/- 0.2(sys) and an exponential cutoff at an energy of 13.8 +/- 2.3(stat) +/- 4.1(sys) TeV. The integral flux above 1 TeV is (1.28 +/- 0.17 (stat) +/- 0.38(sys)) x 10^{-11} cm^{-2} s^{-1}. This result is the first clear measurement of a peak in the spectral energy distribution from a VHE gamma-ray source, likely related to inverse Compton emission. A fit of an Inverse Compton model to the H.E.S.S. spectral energy distribution gives a total energy in non-thermal electrons of ~2 x 10^{45} erg between 5 TeV and 100 TeV, assuming a distance of 290 parsec to the pulsar. The best fit electron power law index is 2.0, with a spectral break at 67 TeV.

The high-frequency peaked BL Lac PG 1553+113 was observed in 2005 with the H.E.S.S. stereoscopic array of imaging atmospheric-Cherenkov telescopes in Namibia. Using the H.E.S.S. standard analysis, an excess was measured at the 4.0 sigma level in these observations (7.6 hours live time). Three alternative, lower-threshold analyses yield >5 sigma excesses. The observed integral flux above 200 GeV is (4.8 +/- 1.3_{stat} +/- 1.0_{syst}) x10^{-12} cm^{-2} s^{-1}, and shows no evidence for variability. The measured energy spectrum is characterized by a very soft power law (photon index of Gamma=4.0 +/- 0.6). Although the redshift of PG 1553+113 is unknown, there are strong indications that it is greater than z=0.25 and possibly larger than z=0.78. The observed spectrum is interpreted in the context of VHE gamma-ray absorption by the Extragalactic Background Light, and is used to place an upper limit on the redshift of z<0.74.

A brief review is given on the astrophysics of cosmic rays in the PeV primary energy range, i.e. the region of the knee.

We fit the data for the binary-lens microlensing event MACHO 98-SMC-1 from 5 different microlensing collaborations and find two distinct solutions characterized by binary separation d and mass ratio q: (d,q)=(0.54,0.50) and (d,q)=(3.65,0.36), where d is in units of the Einstein radius. However, the relative proper motion of the lens is very similar in the two solutions, 1.30 km/s/kpc and 1.48 km/s/kpc, thus confirming that the lens is in the Small Magellanic Cloud. The close binary can be either rotating or approximately static but the wide binary must be rotating at close its maximum allowed rate to be consistent with all the data. We measure limb-darkening coefficients for five bands ranging from I to V. As expected, these progressively decrease with rising wavelength. This is the first measurement of limb darkening for a metal-poor A star.