The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $\gamma$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation with an increase of a factor $\sim$4 in 25 minutes, and a falling time of $\sim$50 minutes. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.

We summarize broadband observations of the TeV-emitting blazar 1ES 1959+650, including optical R-band observations by the robotic telescopes Super-LOTIS and iTelescope, UV observations by Swift UVOT, X-ray observations by the Swift X-ray Telescope (XRT), high-energy gamma-ray observations with the Fermi Large Area Telescope (LAT) and very-high-energy (VHE) gamma-ray observations by VERITAS above 315 GeV, all taken between 17 April 2012 and 1 June 2012 (MJD 56034 and 56079). The contemporaneous variability of the broadband spectral energy distribution is explored in the context of a simple synchrotron self Compton (SSC) model. In the SSC emission scenario, we find that the parameters required to represent the high state are significantly different than those in the low state. Motivated by possible evidence of gas in the vicinity of the blazar, we also investigate a reflected-emission model to describe the observed variability pattern. This model assumes that the non-thermal emission from the jet is reflected by a nearby cloud of gas, allowing the reflected emission to re-enter the blob and produce an elevated gamma-ray state with no simultaneous elevated synchrotron flux. The model applied here, although not required to explain the observed variability pattern, represents one possible scenario which can describe the observations. As applied to an elevated VHE state of 66% of the Crab Nebula flux, observed on a single night during the observation period, the reflected-emission scenario does not support a purely leptonic non-thermal emission mechanism. The reflected emission model does, however, predict a reflected photon field with sufficient energy to enable elevated gamma-ray emission via pion production with protons of energies between 10 and 100 TeV.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

We report on deep observations of the extended TeV gamma-ray source MGRO J1908+06 made with the VERITAS very high energy (VHE) gamma-ray observatory. Previously, the TeV emission has been attributed to the pulsar wind nebula (PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a significance level of 14 standard deviations (14 sigma) and measure a photon index of 2.20 +/- 0.10_stat +/- 0.20_sys. The TeV emission is extended, covering the region near PSR J1907+0602 and also extending towards SNR G40.5--0.5. When fitted with a 2-dimensional Gaussian, the intrinsic extension has a standard deviation of sigma_src = 0.44 +/- 0.02 degrees. In contrast to other TeV PWNe of similar age in which the TeV spectrum softens with distance from the pulsar, the TeV spectrum measured near the pulsar location is consistent with that measured at a position near the rim of G40.5--0.5, 0.33 degrees away.

In March 2013, a flaring episode from the Crab Nebula lasting ~2 weeks was detected by the Fermi-LAT (Large Area Telescope on board the Fermi Gamma-ray Space Telescope). VERITAS provides simultaneous observations throughout this period. During the flare, the Fermi-LAT detected a 20-fold increase in flux above the average synchrotron flux >100 MeV seen from the Crab Nebula. Simultaneous measurements with VERITAS are consistent with the non-variable long-term average Crab Nebula flux at TeV energies. Assuming a linear correlation between the very-high-energy flux change >1 TeV and the flux change seen in the Fermi-LAT band >100 MeV during the period of simultaneous observations, the linear correlation factor can be constrained to be at most 8.6 * 10^-3 with 95% confidence.

HESS J0632+057 is the only gamma-ray binary known so far whose position in the sky allows observations with ground-based observatories both in the northern and southern hemispheres. Here we report on long-term observations of HESS J0632+057 conducted with the VERITAS and H.E.S.S. Cherenkov Telescopes and the X-ray Satellite Swift, spanning a time range from 2004 to 2012 and covering most of the system's orbit. The VHE emission is found to be variable, and is correlated with that at X-ray energies. An orbital period of $315 ^{+6}_{-4}$ days is derived from the X-ray data set, which is compatible with previous results, $P = (321 \pm 5$) days. The VHE light curve shows a distinct maximum at orbital phases close to 0.3, or about 100 days after periastron passage, which coincides with the periodic enhancement of the X-ray emission. Furthermore, the analysis of the TeV data shows for the first time a statistically significant ($> 6.5 \sigma$) detection at orbital phases 0.6--0.9. The obtained gamma-ray and X-ray light curves and the correlation of the source emission at these two energy bands are discussed in the context of the recent ephemeris obtained for the system. Our results are compared to those reported for other gamma-ray binaries.

We report results from TeV gamma-ray observations of the microquasar Cygnus X-3. The observations were made with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) over a time period from 2007 June 11 to 2011 November 28. VERITAS is most sensitive to gamma rays at energies between 85 GeV to 30 TeV. The effective exposure time amounts to a total of about 44 hours, with the observations covering six distinct radio/X-ray states of the object. No significant TeV gamma-ray emission was detected in any of the states, nor with all observations combined. The lack of a positive signal, especially in the states where GeV gamma rays were detected, places constraints on TeV gamma-ray production in Cygnus X-3. We discuss the implications of the results.

We present the results of a multiwavelength observational campaign on the TeV binary system LS I +61 303 with the VERITAS telescope array (>200 GeV), Fermi-LAT (0.3-300 GeV), and Swift-XRT (2-10 keV). The data were taken from December 2011 through January 2012 and show a strong detection in all three wavebands. During this period VERITAS obtained 24.9 hours of quality selected livetime data in which LS I +61 303 was detected at a statistical sig- nificance of 11.9 sigma. These TeV observations show evidence for nightly variability in the TeV regime at a post-trial significance of 3.6 sigma. The combination of the simultaneously obtained TeV and X-ray fluxes do not demonstrate any evidence for a correlation between emission in the two bands. For the first time since the launch of the Fermi satellite in 2008, this TeV detection allows the construction of a detailed MeV-TeV spectral energy distribution from LS I +61 303. This spectrum shows a distinct cutoff in emission near 4 GeV, with emission seen by the VERITAS observations following a simple power-law above 200 GeV. This feature in the spectrum of LS I +61 303, obtained from overlapping observations with Fermi-LAT and VERITAS, may indicate that there are two distinct populations of accelerated particles producing the GeV and TeV emission.

Compilation of papers contributed by the VERITAS Collaboration to the 33rd International Cosmic Ray Conference, held 2-9 July, 2013, in Rio de Janeiro, Brazil.

The Fornax galaxy cluster was observed with the High Energy Stereoscopic System (H.E.S.S.) for a total live time of 14.5 hours, searching for very-high-energy (VHE, E>100 GeV) gamma-rays from dark matter (DM) annihilation. No significant signal was found in searches for point-like and extended emissions. Using several models of the DM density distribution, upper limits on the DM velocity-weighted annihilation cross-section <sigma v> as a function of the DM particle mass are derived. Constraints are derived for different DM particle models, such as those arising from Kaluza-Klein and supersymmetric models. Various annihilation final states are considered. Possible enhancements of the DM annihilation gamma-ray flux, due to DM substructures of the DM host halo, or from the Sommerfeld effect, are studied. Additional gamma-ray contributions from internal bremsstrahlung and inverse Compton radiation are also discussed. For a DM particle mass of 1 TeV, the exclusion limits at 95% of confidence level reach values of <sigma v> ~ 10^-23cm^3s^-1, depending on the DM particle model and halo properties. Additional contribution from DM substructures can improve the upper limits on <\sigma v> by more than two orders of magnitude. At masses around 4.5 TeV, the enhancement by substructures and the Sommerfeld resonance effect results in a velocity-weighted annihilation cross-section upper limit at the level of <\sigma v> ~ 10^-26cm^3s^-1.

The quasar PKS 1510-089 (z=0.361) was observed with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes during high states in the optical and GeV bands, to search for very high energy (VHE, defined as E >= 0.1 TeV) emission. VHE \gamma-rays were detected with a statistical significance of 9.2 standard deviations in 15.8 hours of H.E.S.S. data taken during March and April 2009. A VHE integral flux of I(0.15 TeV < E < 1.0 TeV) = (1.0 +- 0.2 (stat) +- 0.2 (sys) x 10^{-11} cm^{-2}s^{-1} is measured. The best-fit power law to the VHE data has a photon index of \Gamma=5.4 +- 0.7 (stat) +- 0.3 (sys). The GeV and optical light curves show pronounced variability during the period of H.E.S.S. observations. However, there is insufficient evidence to claim statistically significant variability in the VHE data. Because of its relatively high redshift, the VHE flux from PKS 1510-089 should suffer considerable attenuation in the intergalactic space due to the extragalactic background light (EBL). Hence, the measured \gamma-ray spectrum is used to derive upper limits on the opacity due to EBL, which are found to be comparable with the previously derived limits from relatively-nearby BL Lac objects. Unlike typical VHE-detected blazars where the broadband spectrum is dominated by non-thermal radiation at all wavelengths, the quasar PKS 1510-089 has a bright thermal component in the optical to UV frequency band. Among all VHE detected blazars, PKS 1510-089 has the most luminous broad line region (BLR). The detection of VHE emission from this quasar indicates a low level of \gamma-\gamma absorption on the internal optical to UV photon field.

Very high-energy gamma-ray emission from PKS 0447-439 was detected with the H.E.S.S. Cherenkov telescope array in December 2009. This blazar is one of the brightest extragalactic objects in the Fermi Bright Source List and has a hard spectrum in the MeV to GeV range. In the TeV range, a photon index of 3.89 +- 0.37 (stat) +- 0.22 (sys) and a flux normalisation at 1 TeV, Phi(1 TeV) = (3.5 +- 1.1 (stat) +- 0.9 (sys)) x 10^{-13} cm^{-2} s^{-1} TeV^{-1}, were found. The detection with H.E.S.S. triggered observations in the X-ray band with the Swift and RXTE telescopes. Simultaneous UV and optical data from Swift UVOT and data from the optical telescopes ATOM and ROTSE are also available. The spectrum and light curve measured with H.E.S.S. are presented and compared to the multi-wavelength data at lower energies. A rapid flare is seen in the Swift XRT and RXTE data, together with a flux variation in the UV band, at a time scale of the order of one day. A firm upper limit of z < 0.59 on the redshift of PKS 0447-439 is derived from the combined Fermi-LAT and H.E.S.S. data, given the assumptions that there is no upturn in the intrinsic spectrum above the Fermi-LAT energy range and that absorption on the Extragalactic Background Light (EBL) is not weaker than the lower limit provided by current models. The spectral energy distribution is well described by a simple one-zone Synchrotron Self-Compton (SSC) scenario, if the redshift of the source is less than z <~ 0.4.

Aim. In this paper we present very high energy (VHE; E>100 GeV) data from the \gamma-ray binary system PSR B1259-63/LS 2883 taken around its periastron passage (15th of December 2010) with the High Energy Stereoscopic System (H.E.S.S.) of Cherenkov Telescopes. We aim to search for a possible TeV counterpart of the GeV flare detected by the Fermi LAT. In addition, we aim to study the current periastron passage in the context of previous observations taken at similar orbital phases, testing the repetitive behavior of the source. Methods. Observations at VHE were conducted with H.E.S.S. from 9th to 16th of January 2011. The total dataset amounts to around 6 h of observing time. Results. The source is detected in the 2011 data at a significance level of 11.5\sigma\ revealing an averaged integral flux above 1 TeV of (1.01 \pm 0.18_{stat} \pm 0.20_{sys}) \times 10^{-12} cm^{-2}s^{-1}. The differential energy spectrum follows a power-law shape with a spectral index \Gamma = 2.92 \pm 0.30_{stat} \pm 0.20_{sys} and a flux normalisation at 1 TeV of N_{0} = 1.95 \pm 0.32_{stat} \pm 0.39_{sys}) \times 10^{-12} TeV^{-1} cm^{-2} s^{-1}. The measured lightcurve does not show any evidence for variability of the source on the daily scale. Conclusions. The measured integral flux and the spectral shape of the 2011 data are compatible with the results obtained around previous periastron passages. The absence of variability in the H.E.S.S. data indicates that the GeV flare observed by Fermi LAT in the time period covered also by H.E.S.S. observations originates in a different physical scenario than the TeV emission. Additionaly, new results compared to those obtained in the observations which were performed in 2004 at a similar orbital phase, further support the hypothesis of the repetitive behavior of the source.

Globular clusters (GCs) are established emitters of high-energy (HE, 100 MeV<E<100 GeV) \gamma-ray radiation which could originate from the cumulative emission of the numerous millisecond pulsars (msPSRs) in the clusters' cores or from inverse Compton (IC) scattering of relativistic leptons accelerated in the GC environment. These stellar clusters could also constitute a new class of sources in the very-high-energy (VHE, E>100 GeV) \gamma-ray regime, judging from the recent detection of a signal from the direction of Terzan 5 with the H.E.S.S. telescope array. We searched for point-like and extended VHE \gamma-ray emission from 15 GCs serendipitously covered by H.E.S.S observations and also performed a stacking analysis combining the data from all GCs to investigate the hypothesis of a population of faint emitters. Assuming IC emission as the origin of the VHE \gamma-ray signal from the direction of Terzan 5, we calculated the expected \gamma-ray flux from each of the 15 GCs, based on their number of millisecond pulsars, their optical brightness and the energy density of background photon fields. We did not detect significant VHE \gamma-ray emission from any of the 15 GCs in either of the two analyses. Given the uncertainties related to the parameter determinations, the obtained flux upper limits allow to rule out the simple IC/msPSR scaling model for NGC 6388 and NGC 7078. The upper limits derived from the stacking analyses are factors between 2 and 50 below the flux predicted by the simple leptonic scaling model, depending on the assumed source extent and the dominant target photon fields. Therefore, Terzan 5 still remains exceptional among all GCs, as the VHE \gamma-ray emission either arises from extra-ordinarily efficient leptonic processes, or from a recent catastrophic event, or is even unrelated to the GC itself.

Gamma-ray line signatures can be expected in the very-high-energy (VHE; E_\gamma > 100 GeV) domain due to self-annihilation or decay of dark matter (DM) particles in space. Such a signal would be readily distinguishable from astrophysical \gamma-ray sources that in most cases produce continuous spectra which span over several orders of magnitude in energy. Using data collected with the H.E.S.S. \gamma-ray instrument, upper limits on line-like emission are obtained in the energy range between ~500 GeV and ~25 TeV for the central part of the Milky Way halo and for extragalactic observations, complementing recent limits obtained with the Fermi-LAT instrument at lower energies. No statistically significant signal could be found. For monochromatic \gamma-ray line emission, flux limits of (2x10^-7 - 2x10^-5) m^-2 s^-1 sr^-1 and (1x10^-8 - 2x10^-6) m^-2 s^-1 sr^-1 are obtained for the central part of the Milky Way halo and extragalactic observations, respectively. For a DM particle mass of 1 TeV, limits on the velocity-averaged DM annihilation cross section < \sigma v >(\chi\chi -> \gamma\gamma) reach ~10^-27 cm^3 s^-1, based on the Einasto parametrization of the Galactic DM halo density profile.

The extragalactic background light (EBL) is the diffuse radiation with the second highest energy density in the Universe after the cosmic microwave background. The aim of this study is the measurement of the imprint of the EBL opacity to gamma-rays on the spectra of the brightest extragalactic sources detected with the High Energy Stereoscopic System (H.E.S.S.). The originality of the method lies in the joint fit of the EBL optical depth and of the intrinsic spectra of the sources, assuming intrinsic smoothness. Analysis of a total of ~10^5 gamma-ray events enables the detection of an EBL signature at the 8.8 std dev level and constitutes the first measurement of the EBL optical depth using very-high energy (E>100 GeV) gamma-rays. The EBL flux density is constrained over almost two decades of wavelengths (0.30-17 microns) and the peak value at 1.4 micron is derived as 15 +/- 2 (stat) +/- 3 (sys) nW / m^2 sr.

The previously unidentified very high-energy (VHE; E > 100 GeV) \gamma-ray source HESS J1303-631, discovered in 2004, is re-examined including new data from the H.E.S.S. Cherenkov telescope array. Archival data from the XMM-Newton X-ray satellite and from the PMN radio survey are also examined. Detailed morphological and spectral studies of VHE \gamma-ray emission as well as of the XMM-Newton X-ray data are performed. Significant energy-dependent morphology of the \gamma-ray source is detected with high-energy emission (E > 10 TeV) positionally coincident with the pulsar PSR J1301-6305 and lower energy emission (E <2 TeV) extending \sim 0.4^{\circ} to the South-East of the pulsar. The spectrum of the VHE source can be described with a power-law with an exponential cut-off N_{0} = (5.6 \pm 0.5) X 10^{-12} TeV^-1 cm^-2 s^-1, \Gamma = 1.5 \pm 0.2) and E_{\rm cut} = (7.7 \pm 2.2) TeV. The PWN is also detected in X-rays, extending \sim 2-3' from the pulsar position towards the center of the \gamma-ray emission region. The spectral energy distribution (SED) is well described by a one zone leptonic scenario which, with its associated caveats, predicts a very low average magnetic field for this source. Significant energy-dependent morphology of this source, as well as the identification of an associated X-ray PWN from XMM-Newton observations enable identification of the VHE source as an evolved PWN associated to the pulsar PSR J1303-6305. However, the large discrepancy in emission region sizes and the low level of synchrotron radiation suggest a multi-population leptonic nature. The low implied magnetic field suggests that the PWN has undergone significant expansion. This would explain the low level of synchrotron radiation and the difficulty in detecting counterparts at lower energies, the reason this source was originally classified as a "dark" VHE \gamma-ray source.

Vela X is a region of extended radio emission in the western part of the Vela constellation: one of the nearest pulsar wind nebulae (PWNe), and associated with the energetic Vela pulsar (PSR B0833-45). Extended very-high-energy (VHE) $\gamma$-ray emission (HESS $\mathrm{J0835\mhyphen 455}$) was discovered using the H.E.S.S. experiment in 2004. The VHE $\gamma$-ray emission was found to be coincident with a region of X-ray emission discovered with ${\it ROSAT}$ above 1.5 keV (the so-called \textit{Vela X cocoon}): a filamentary structure extending southwest from the pulsar to the centre of Vela X. A deeper observation of the entire Vela X nebula region, also including larger offsets from the cocoon, has been performed with H.E.S.S. This re-observation was carried out in order to probe the extent of the non-thermal emission from the Vela X region at TeV energies and to investigate its spectral properties. In order to increase the sensitivity to the faint $\gamma$-ray emission from the very extended Vela X region, a multivariate analysis method combining three complementary reconstruction techniques of Cherenkov-shower images is applied for the selection of $\gamma$-ray events. The analysis is performed with the On/Off background method, which estimates the background from separate observations pointing away from Vela X; towards regions free of $\gamma$-ray sources but with comparable observation conditions. The $\gamma$-ray surface brightness over the large Vela X region reveals that the detection of non-thermal VHE $\gamma$-ray emission from the PWN HESS $\mathrm{J0835\mhyphen 455}$ is statistically significant over a region of radius 1.2$^{\circ}$ around the position $\alpha$ = 08$^{\mathrm{h}}$ 35$^{\mathrm{m}}$ 00$^{\mathrm{s}}$, $\delta$ = -45$^{\circ}$ 36$^{\mathrm{\prime}}$ 00$^{\mathrm{\prime}\mathrm{\prime}}$ (J2000).

We present the significant detection of the first extragalactic pulsar wind nebula (PWN) detected in gamma rays, N157B, located in the large Magellanic Cloud (LMC). Pulsars with high spin-down luminosity are found to power energised nebulae that emit gamma rays up to energies of several tens of TeV. N157B is associated with PSRJ0537-6910, which is the pulsar with the highest known spin-down luminosity. The High Energy Stereoscopic System telescope array observed this nebula on a yearly basis from 2004 to 2009 with a dead-time corrected exposure of 46 h. The gamma-ray spectrum between 600 GeV and 12 TeV is well-described by a pure power-law with a photon index of 2.8 \pm 0.2(stat) \pm 0.3(syst) and a normalisation at 1 TeV of (8.2 \pm 0.8(stat) \pm 2.5(syst)) \times 10^-13 cm^-2s^-1TeV^-1. A leptonic multi-wavelength model shows that an energy of about 4 \times 10^49erg is stored in electrons and positrons. The apparent efficiency, which is the ratio of the TeV gamma-ray luminosity to the pulsar's spindown luminosity, 0.08% \pm 0.01%, is comparable to those of PWNe found in the Milky Way. The detection of a PWN at such a large distance is possible due to the pulsar's favourable spin-down luminosity and a bright infrared photon-field serving as an inverse-Compton-scattering target for accelerated leptons. By applying a calorimetric technique to these observations, the pulsar's birth period is estimated to be shorter than 10 ms.

Very-high-energy (VHE; E >100 GeV) and high-energy (HE; 100 MeV < E < 100 GeV) data from \gamma-ray observations performed with the H.E.S.S. telescope array and the Fermi-LAT instrument, respectively, are analysed in order to investigate the non-thermal processes in the starburst galaxy NGC 253. The VHE \gamma-ray data can be described by a power law in energy with differential photon index \Gamma=2.14 \pm 0.18_stat \pm 0.30_sys and differential flux normalisation at 1 TeV of F_0 = (9.6 \pm 1.5_stat (+5.7,-2.9)_sys) x 10^{-14} TeV^{-1} cm^{-2} s^{-1}. A power-law fit to the differential HE \gamma-ray spectrum reveals a photon index of \Gamma=2.24 \pm 0.14_stat \pm 0.03_sys and an integral flux between 200 MeV and 200 GeV of F(0.2-200 GeV) = (4.9 \pm 1.0_stat \pm 0.3_sys) x 10^{-9} cm^{-2} s^{-1}. No evidence for a spectral break or turnover is found over the dynamic range of both the LAT instrument and the H.E.S.S. experiment: a combined fit of a power law to the HE and VHE \gamma-ray data results in a differential photon index \Gamma=2.34 \pm 0.03 with a p-value of 30%. The \gamma-ray observations indicate that at least about 20% of the energy of the cosmic rays capable of producing hadronic interactions is channeled into pion production. The smooth alignment between the spectra in the HE and VHE \gamma-ray domain suggests that the same transport processes dominate in the entire energy range. Advection is most likely responsible for charged particle removal from the starburst nucleus from GeV to multiple TeV energies. In a hadronic scenario for the \gamma-ray production, the single overall power-law spectrum observed would therefore correspond to the mean energy spectrum produced by the ensemble of cosmic-ray sources in the starburst region.

In some galaxy clusters powerful AGN have blown bubbles with cluster scale extent into the ambient medium. The main pressure support of these bubbles is not known to date, but cosmic rays are a viable possibility. For such a scenario copious gamma-ray emission is expected as a tracer of cosmic rays from these systems. Hydra A, the closest galaxy cluster hosting a cluster scale AGN outburst, located at a redshift of 0.0538, is investigated for being a gamma-ray emitter with the High Energy Stereoscopic System (H.E.S.S.) array and the Fermi Large Area Telescope (Fermi-LAT). Data obtained in 20.2 hours of dedicated H.E.S.S. observations and 38 months of Fermi-LAT data, gathered by its usual all-sky scanning mode, have been analyzed to search for a gamma-ray signal. No signal has been found in either data set. Upper limits on the gamma-ray flux are derived and are compared to models. These are the first limits on gamma-ray emission ever presented for galaxy clusters hosting cluster scale AGN outbursts. The non-detection of Hydra A in gamma-rays has important implications on the particle populations and physical conditions inside the bubbles in this system. For the case of bubbles mainly supported by hadronic cosmic rays, the most favorable scenario, that involves full mixing between cosmic rays and embedding medium, can be excluded. However, hadronic cosmic rays still remain a viable pressure support agent to sustain the bubbles against the thermal pressure of the ambient medium. The largest population of highly-energetic electrons which are relevant for inverse-Compton gamma-ray production is found in the youngest inner lobes of Hydra A. The limit on the inverse-Compton gamma-ray flux excludes a magnetic field below half of the equipartition value of 16 muG in the inner lobes.

The massive binary system Eta Carinae and the surrounding HII complex, the Carina Nebula, are potential particle acceleration sites from which very-high-energy (VHE; E > 100 GeV) \gamma-ray emission could be expected. This paper presents data collected during VHE \gamma-ray observations with the H.E.S.S. telescope array from 2004 to 2010, which cover a full orbit of Eta Carinae. In the 33.1-hour data set no hint of significant \gamma-ray emission from Eta Carinae has been found and an upper limit on the \gamma-ray flux of 7.7 x 10-13 ph cm-2 s-1 (99% confidence level) is derived above the energy threshold of 470 GeV. Together with the detection of high-energy (HE; 0.1 GeV > E > 100 GeV) \gamma-ray emission by the Fermi-LAT up to 100 GeV, and assuming a continuation of the average HE spectral index into the VHE domain, these results imply a cut-off in the \gamma-ray spectrum between the HE and VHE \gamma-ray range. This could be caused either by a cut-off in the accelerated particle distribution or by severe \gamma-\gamma\ absorption losses in the wind collision region. Furthermore, the search for extended \gamma-ray emission from the Carina Nebula resulted in an upper limit on the \gamma-ray flux of 4.2 x 10-12 ph cm-2 s-1 (99% confidence level). The derived upper limit of ~23 on the cosmic-ray enhancement factor is compared with results found for the old-age mixed-morphology supernova remnant W 28.

1RXS J101015.9-311909 is a galaxy located at a redshift of z=0.14 hosting an active nucleus belonging to the class of bright BL Lac objects. Observations at high (HE, E > 100 MeV) and very high (VHE, E > 100 GeV) energies provide insights into the origin of very energetic particles present in such sources and the radiation processes at work. We report on results from VHE observations performed between 2006-10 with H.E.S.S. H.E.S.S. data have been analysed with enhanced analysis methods, making the detection of faint sources more significant. VHE emission at a position coincident with 1RXS J101015.9-311909 is detected with H.E.S.S. for the first time. In a total good-quality livetime of about 49 h, we measure 263 excess counts, corresponding to a significance of 7.1\sigma. The photon spectrum above 0.2 TeV can be described by a power-law with a photon index of \Gamma\ = 3.08\pm0.42_{stat}\pm0.20_{sys}. The integral flux above 0.2 TeV is about 0.8% of the flux of the Crab nebula and shows no significant variability over the time reported. In addition, public Fermi/LAT data are analysed to search for high energy emission from the source. The Fermi/LAT HE emission is significant at 8.3\sigma\ in the chosen 25-month dataset. UV and X-ray contemporaneous observations with the Swift satellite in May 2007 are also reported, together with optical observations performed with the ATOM telescope located at the H.E.S.S. site. Swift observations reveal an absorbed X-ray flux of F_{0.3-7 keV} = 1.04^{+0.04}_{-0.05} \times 10^{-11} erg.cm^{-2}.s^{-1} in the 0.3-7 keV range. Finally, all the available data are used to study the source's multi-wavelength properties. The SED can be reproduced using a simple one-zone SSC model with emission from a region with a Doppler factor of 30 and a magnetic field between 0.025 and 0.16 G. These parameters are similar to those obtained for other sources of this type.

The observational coverage with HESS of the Carina region in VHE gamma-rays benefits from deep exposure (40 h) of the neighboring open cluster Westerlund 2. The observations have revealed a new extended region of VHE gamma-ray emission. The new VHE source HESS J1018-589 shows a bright, point-like emission region positionally coincident with SNR G284.3-1.8 and 1FGL J1018.6 - 5856 and a diffuse extension towards the direction of PSR J1016-5857. A soft Gamma=2.7+-0.5 photon index, with a differential flux at 1TeV of N0=(4.2+-1.1)10^-13 TeV^-1 cm^-2 s^-1 is found for the point-like source, whereas the total emission region including the diffuse emission region is well fit by a power-law function with spectral index Gamma=2.9+-0.4 and differential flux at 1TeV of N0=(6.8+-1.6) 10^-13 TeV^-1 cm^-2 s^-1. This H.E.S.S. detection motivated follow-up X-ray observations with the XMM-Newton satellite to investigate the origin of the VHE emission. The analysis of the XMM-Newton data resulted in the discovery of a bright, non-thermal point-like source (XMMU J101855.4-58564) with a photon index of Gamma=1.65+-0.08 in the center of SNRG284.3-1.8, and a thermal, extended emission region coincident with its bright northern filament. The characteristics of this thermal emission are used to estimate the plasma density in the region as n~0.5 cm^-3(2.9kpc/d)^2. The position of XMMUJ101855.4-58564 is compatible with the position reported by the Fermi-LAT collaboration for the binary system 1FGL J1018.6-5856 and the variable Swift XRT source identified with it. The new X-ray data are used alongside archival multi-wavelength data to investigate the relationship between the VHE gamma-ray emission from HESSJ1018-589 and the various potential counterparts in the Carina arm region.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

1ES 0414+009 (z = 0.287) is a distant high-frequency-peaked BL Lac object, and has long been considered a likely emitter of very-high energy (VHE, E>100 GeV) gamma-rays due to its high X-ray and radio flux. Observations in the VHE gamma-ray band and across the electromagnetic spectrum can provide insights into the origin of highly energetic particles present in the source and the radiation processes at work. Because of the distance of the source, the gamma-ray spectrum might provide further limits on the level of the Extragalactic Background Light (EBL). We report observations made between October 2005 and December 2009 with H.E.S.S., an array of four imaging atmospheric Cherenkov telescopes. Observations at high energies (HE, 100 MeV - 100 GeV) with the Fermi-LAT instrument in the first 20 months of its operation are also reported. To complete the multi-wavelength picture, archival UV and X-ray observations with the Swift satellite and optical observations with the ATOM telescope are also used. Based on the observations with H.E.S.S., 1ES 0414+009 is detected for the first time in the VHE band. An excess of 224 events is measured, corresponding to a significance of 7.8 sigma. The photon spectrum of the source is well described by a power law, with photon index of 3.45 \pm 0.25stat \pm 0.20syst. The integral flux above 200 GeV is (1.88 \pm 0.20stat \pm 0.38syst) \times10-12 cm-2 s-1. Observations with the Fermi-LAT in the first 20 months of operation show a flux between 200 MeV and 100 GeV of (2.3 \pm 0.2stat) \times 10-9 erg cm-2 s-1, and a spectrum well described by a power-law function with a photon index 1.85 \pm 0.18. Swift/XRT observations show an X-ray flux between 2 and 10 keV of (0.8 - 1) \times 10-11 erg cm-2 s-1, and a steep spectrum (2.2 - 2.3). Combining X-ray with optical-UV data, a fit with a log-parabolic function locates the synchrotron peak around 0.1 keV. ...

Results obtained in very-high-energy (VHE; E > 100 GeV) \gamma-ray observations performed with the H.E.S.S. telescope array are used to investigate particle acceleration processes in the vicinity of the young massive stellar cluster Westerlund 1 (Wd 1). Imaging of Cherenkov light from \gamma-ray induced particle cascades in the Earth's atmosphere is used to search for VHE \gamma\ rays from the region around Wd 1. Possible catalogued counterparts are searched for and discussed in terms of morphology and energetics of the H.E.S.S. source. The detection of the degree-scale extended VHE \gamma-ray source HESS J1646-458 is reported based on 45 hours of H.E.S.S. observations performed between 2004 and 2008. The VHE \gamma-ray source is centred on the nominal position of Wd 1 and detected with a total statistical significance of ~20\sigma. The emission region clearly extends beyond the H.E.S.S. point-spread function (PSF). The differential energy spectrum follows a power law in energy with an index of \Gamma=2.19 \pm 0.08_{stat} \pm 0.20_{sys} and a flux normalisation at 1 TeV of \Phi_0 = (9.0 \pm 1.4_{stat} \pm 1.8_{sys}) x 10^{-12} TeV^{-1} cm^{-2} s^{-1}. The integral flux above 0.2 TeV amounts to (5.2 \pm 0.9) x 10^{-11} cm^{-2} s^{-1}. Four objects coincident with HESS J1646-458 are discussed in the search of a counterpart, namely the magnetar CXOU J164710.2-455216, the X-ray binary 4U 1642-45, the pulsar PSR J1648-4611 and the massive stellar cluster Wd 1. In a single-source scenario, Wd 1 is favoured as site of VHE particle acceleration. Here, a hadronic parent population would be accelerated within the stellar cluster. Beside this, there is evidence for a multi-source origin, where a scenario involving PSR J1648-4611 could be viable to explain parts of the VHE \gamma-ray emission of HESS J1646-458.

We report on the latest discovery of an extended Very High Energy (VHE) gamma-ray source near the 67 ms pulsar PSR J1831-0952 during the H.E.S.S. Galactic Plane Survey (GPS). The dispersion measure distance of the pulsar (4.3 kpc) would imply that less than 1% of its spin-down energy is required to provide the observed VHE luminosity of the source. No other plausible counterparts have yet been found through preliminary multi- wavelength searches. The most likely scenario is that the VHE emission originates from the -yet unseen at other wavelengths-wind nebula of PSR J1831-0952. If so this would constitute another case of a gamma-ray discovered pulsar wind nebula.

In this study we set our system in order to study the energy spectrum of single, double and triple particles, detected in a scintillation detector. The goal of doing this experiment was to determine the probability of number of particles (single, double or triple) detected, from the energy spectrum in any given energy spectrum. The results of experiment will be used in our extensive air shower array.

A new technique has been devised for the analysis of extensive air shower data in observing the effect of the moon on this data. In this technique the number of EAS events with arrival directions falling in error circles centered about the moving moon is compared to the mean number of events falling in error circles with centers randomly chosen in the sky. For any assumed angular radius of the error circle the deficit in EAS event count in the direction of moon i.e., Nsky-Nmoon which is a moon-related effect is interpreted as the shadow of the moon. A simple theoretical model has been developed to relate Nsky to the angular radius of the error circle and has been applied to the counts from the moon's direction in order to extract the physical parameters of the shadow of the moon. The technique and the theoretical model has been used on 1.7 *10^5 EAS events recorded at Alborz observatory.

The study of the anisotropy of the arrival directions is an essential tool to investigate the origin and propagation of cosmic rays primaries. A simple way of recording many cosmic rays is to record coincidences between a number of detectors. We have monitored multi-TeV cosmic rays by a small array of water cherenkov detectors in Tehran(35 43 N, 51 20 E, 1200m a.s.l). More than 1.1*10^6 extensive air shower events were recorded. In addition to the Compton- Getting effect due to the motion of the earth in the Galaxy, an anisotropy has been observed which is due to a unidirectional anisotropy of cosmic ray flow along the Galactic arms.

Ultra-High-Energy (UHE) ($E>100 $TeV) Extensive Air Showers (EASs) have been monitored for a period of five years (1997-2003), using a small array of scintillation detectors in Tehran, Iran. The data have been analyzed to take in to account of the dependence of source counts on zenith angle. Because of varying thickness of the overlaying atmosphere, the shower count rate is extremely dependent on zenith angle. During a calendar year different sources come in the field of view of the array at varying zenith angles and have different effective observation time equivalent to zenith in a day. High energy gamma-ray sources from the EGRET third catalogue where observed and the data were analyzed using an excess method. Upper limits were obtained for 10 EGRET sources. Then we investigated the EAS event rates for these 10 sources and obtained a flux for each of them using parameters of our experiment results and simulations. Finally we investigated the gamma-ray spectrum in the UHE range using these fluxes with reported fluxes of the EGRET sources.