The VERITAS array of Cherenkov telescopes has carried out a deep observational program on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of good quality selected data, taken between January 2010 and May 2011. No significant $\gamma$-ray emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal galaxy currently known. We derive stringent bounds on various annihilating and decaying dark matter particle models. The upper limits on the velocity-weighted annihilation cross-section are $\mathrm{<\sigma v >^{95% CL} \lesssim 10^{-23} cm^{3} s^{-1}}$, improving our limits from previous observations of dwarf spheroidal galaxies by at least a factor of two for dark matter particle masses $\mathrm{m_{\chi}\gtrsim 300 GeV}$. The lower limits on the decay lifetime are at the level of $\mathrm{\tau^{95% CL} \gtrsim 10^{24} s}$. Finally, we address the interpretation of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario.

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.

We report the discovery of an unidentified, extended source of very-high-energy (VHE) gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hours of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0.23^{\circ} \pm 0.03^{\circ} (stat)+0.04^{\circ}_{-0.02}^{\circ}(sys) and its spectrum is well-characterized by a differential power law (dN/dE = N_0 \times (E/TeV)^{-\Gamma}) with a photon index of {\Gamma} = 2.37 \pm 0.14 (stat) \pm 0.20 (sys) and a flux normalization of N0 = 1.5 \pm 0.2 (stat) \pm 0.4(sys) \times 10^-12 ph TeV^{-1} cm^{-2} s^{-1}. This yields an integral flux of 5.2 \pm 0.8 (stat) \pm 1.4 (sys) \times 10^-12 ph cm^{-2} s^{-1} above 320 GeV, corresponding to 3.7% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the supernova remnant shock.

Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated, suggesting that clusters may also be sources of very high-energy (VHE; E>100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cherenkov telescopes, with complementing Fermi-LAT observations at GeV energies. No significant gamma-ray emission from the Coma cluster was detected. Integral flux upper limits at the 99% confidence level were measured to be on the order of (2-5)*10^-8\ ph. m^-2 s^-1 (VERITAS, >220 GeV} and ~2*10^-6 ph. m^-2 s^-1 (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain CRs and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be < 16% from VERITAS data and < 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be <50%. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of (2 - 5.5) muG, depending on the radial magnetic-field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this {renders} the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark-matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally-averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, <\sigma v>. (abr.)

We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between January 2008 and February 2011, resulting in 56.2 hours of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as 5.2 +/- 1.1_stat +/- 2.6_sys * 10^-12 photons cm^-2 s^-1 above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an photon index of Gamma 3.4 +/- 0.5_stat +/- 0.3_sys and a flux normalization of 1.6 +/- 0.3_stat +/- 0.8_sys * 10^-11 photons cm^-2 s^-1 at 300 GeV. We also present multiwavelength results taken in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.

We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles, due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 \times 10^(-12) erg cm^(-2) s^(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.

We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high-energy (VHE; E > 100 GeV) gamma-ray emission from the high-frequency-peaked BL Lac object RBS 0413. VERITAS, a ground-based gamma-ray observatory, detected VHE gamma rays from RBS 0413 with a statistical significance of 5.5 standard deviations (sigma) and a gamma-ray flux of (1.5 \pm 0.6stat \pm 0.7syst) \times 10^(-8) photons m^(-2) s^(-1) (\sim 1% of the Crab Nebula flux) above 250 GeV. The observed spectrum can be described by a power law with a photon index of 3.18 \pm 0.68stat \pm 0.30syst. Contemporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope detected HE gamma rays from RBS 0413 with a statistical significance of more than 9 sigma, a power-law photon index of 1.57 \pm 0.12stat +0.11sys -0.12sys and a gamma-ray flux between 300 MeV and 300 GeV of (1.64 \pm 0.43stat +0.31sys -0.22sys) \times 10^(-5) photons m^(-2) s^(-1). We present the results from Fermi-LAT and VERITAS, including a spectral energy distribution modeling of the gamma-ray, quasi-simultaneous X-ray (Swift-XRT), ultraviolet (Swift-UVOT) and R-band optical (MDM) data. We find that, if conditions close to equipartition are required, both the combined synchrotron self-Compton/external-Compton and the lepto-hadronic models are preferred over a pure synchrotron self-Compton model.

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.

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.

Multiwavelength (MWL) observations of the blazar PKS 2155-304 during two weeks in July and August 2006, the period when two exceptional flares at very high energies (VHE, E>= 100 GeV) occurred, provide a detailed picture of the evolution of its emission. The complete data set from this campaign is presented, including observations in VHE gamma-rays (H.E.S.S.), X-rays (RXTE, CHANDRA, SWIFT XRT), optical (SWIFT UVOT, Bronberg, Watcher, ROTSE), and in the radio band (NRT, HartRAO, ATCA). Optical and radio light curves from 2004 to 2008 are compared to the available VHE data from this period, to put the 2006 campaign into the context of the long-term evolution of the source. The X-ray and VHE gamma-ray emission are correlated during the observed high state of the source, but show no direct connection with longer wavelengths. The long-term flux evolution in the optical and radio bands is found to be correlated and shows that the source reaches a high state at long wavelengths after the occurrence of the VHE flares. Spectral hardening is seen in the SWIFT XRT data. The nightly averaged high-energy spectra of the non-flaring nights can be reproduced by a stationary one-zone SSC model, with only small variations in the parameters. The spectral and flux evolution in the high-energy band during the night of the second VHE flare is modelled with multi-zone SSC models, which can provide relatively simple interpretations for the hour time-scale evolution of the high-energy emission, even for such a complex data set. For the first time in this type of source, a clear indication is found for a relation between high activity at high energies and a long-term increase in the low frequency fluxes.

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. ...

VERITAS has been monitoring the very-high-energy (VHE; >100GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +- 0.16) X 10^{-11} cm^{-2} s^{-1} at energies above 350GeV. In 2010 April, VERITAS detected a flare from M87 with peak flux of (2.71 +- 0.68) X 10^{-11} cm^{-2} s^{-1} for E>350GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90^{+0.22}_{-0.15} days. The shortest detected exponential rise time is three times as long, at 2.87^{+1.65}_{-0.99} days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M87 and the constraints they provide on the physical size of the emitting region.

We present the results of sixteen Swift-triggered GRB follow-up observations taken with the VERITAS telescope array from January, 2007 to June, 2009. The median energy threshold and response time of these observations was 260 GeV and 320 s, respectively. Observations had an average duration of 90 minutes. Each burst is analyzed independently in two modes: over the whole duration of the observations and again over a shorter time scale determined by the maximum VERITAS sensitivity to a burst with a t^-1.5 time profile. This temporal model is characteristic of GRB afterglows with high-energy, long-lived emission that have been detected by the Large Area Telescope (LAT) on-board the Fermi satellite. No significant VHE gamma-ray emission was detected and upper limits above the VERITAS threshold energy are calculated. The VERITAS upper limits are corrected for gamma-ray extinction by the extragalactic background light (EBL) and interpreted in the context of the keV emission detected by Swift. For some bursts the VHE emission must have less power than the keV emission, placing constraints on inverse Compton models of VHE emission.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

Compilation of papers contributed by the VERITAS Collaboration to the 32nd International Cosmic Ray Conference, held 11-18 August 2011 in Beijing, China.

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 VERITAS discovery of very-high-energy (VHE) gamma- ray emission above 200 GeV from the high-frequency-peaked BL Lac object RXJ0648.7+1516 (GBJ0648+1516), associated with 1FGLJ0648.8+1516. The photon spectrum above 200 GeV is fit by a power law dN/dE = F0(E/E0)-{\Gamma} with a photon index {\Gamma} of 4.4 {\pm} 0.8stat {\pm}0.3syst and a flux normalization F0 of (2.3 {\pm}0.5stat {\pm} 1.2sys) {\times} 10-11 TeV-1cm-2s-1 with E0 = 300 GeV. No VHE vari- ability is detected during VERITAS observations of RXJ0648.7+1516 between 2010 March 4 and April 15. Following the VHE discovery, the optical identifica- tion and spectroscopic redshift were obtained using the Shane 3-m Telescope at the Lick Observatory, showing the unidentified object to be a BL Lac type with a redshift of z = 0.179. Broadband multiwavelength observations contemporaneous with the VERITAS exposure period can be used to sub-classify the blazar as a high-frequency-peaked BL Lac (HBL) object, including data from the MDM ob- servatory, Swift-UVOT and XRT, and continuous monitoring at photon energies above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the absence of undetected, high-energy rapid variability, the one-zone synchrotron self-Compton model (SSC) overproduces the high-energy gamma-ray emission measured by the Fermi-LAT over 2.3 years. The SED can be parameterized sat- isfactorily with an external-Compton or lepto-hadronic model, which have two and six additional free parameters, respectively, compared to the one-zone SSC model.

The VERITAS telescope array has been operating smoothly since 2007, and has detected gamma-ray emission above 100 GeV from 40 astrophysical sources. These include blazars, pulsar wind nebulae, supernova remnants, gamma-ray binary systems, a starburst galaxy, a radio galaxy, the Crab pulsar, and gamma-ray sources whose origin remains unidentified. In 2009, the array was reconfigured, greatly improving the sensitivity. We summarize the current status of the observatory, describe some of the scientific highlights since 2009, and outline plans for the future.

We report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 Gigaelectronvolts (GeV) with the VERITAS array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 Megaelectronvolts (MeV) and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.

Several newly discovered very-high-energy (VHE; E > 100 GeV) gamma-ray sources in the Galaxy are thought to be associated with energetic pulsars. Among them, middle-aged (> 1E+4 yr) systems exhibit large centre-filled VHE nebulae, offset from the pulsar position, which result from the complex relationship between the pulsar wind and the surrounding medium, and reflect the past evolution of the pulsar. Imaging Atmospheric Cherenkov Telescopes (IACTs) have been successful in revealing extended emission from these sources in the VHE regime. Together with radio and X-ray observations, this observational window allows one to probe the energetics and magnetic field inside these large-scale nebulae. H.E.S.S., with its large field of view, angular resolution of < 0.1deg and unprecedented sensitivity, has been used to discover a large population of such VHE sources. In this paper, the H.E.S.S. data from the continuation of the Galactic Plane Survey (-80deg < l < 60deg, |b| < 3deg), together with the existing multi-wavelength observations, are used. A new VHE gamma-ray source was discovered at R.A. (J2000) = 13h56m00s, Dec. (J2000) = -64d30m00s with a 2' statistical error in each coordinate, namely HESS J1356-645. The source is extended, with an intrinsic Gaussian width of (0.20 +/- 0.02)deg. Its integrated energy flux between 1 and 10 TeV of 8E-12 erg cm-2 s-1 represents ~ 11% of the Crab Nebula flux in the same energy band. The energy spectrum between 1 and 20 TeV is well described by a power law dN/dE ~ E-Gamma with photon index Gamma = 2.2 +/- 0.2stat +/- 0.2sys. The inspection of archival radio images at three frequencies and the analysis of X-ray data from ROSAT/PSPC and XMM-Newton/MOS reveal the presence of faint non-thermal diffuse emission coincident with HESS J1356-645. HESS J1356-645 is most likely associated with the young and energetic pulsar PSR J1357-6429 (Abridged)

We report on very-high-energy ($>$100 GeV) gamma-ray observations of Swift J164449.3+573451, an unusual transient object first detected by the {\it Swift} Observatory and later detected by multiple radio, optical and X-ray observatories. A total exposure of 28 hours was obtained on Swift J164449.3+573451 with VERITAS during 2011 March 28 -- April 15. We do not detect the source and place a differential upper limit on the emission at 500 GeV during these observations of $1.4 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ (99% confidence level). We also present time-resolved upper limits and use a flux limit averaged over the X-ray flaring period to constrain various emission scenarios that can accommodate both the radio-through-X-ray emission detected from the source and the lack of detection by VERITAS.

The H.E.S.S. very-high-energy (VHE, E > 0.1 TeV) gamma-ray telescope system has discovered a new source, HESS J1747-248. The measured integral flux is (1.2 +/- 0.3) \times 10^-12 cm-2 s-1 above 440 GeV for a power-law photon spectral index of 2.5 +/- 0.3 stat +/- 0.2 sys. The VHE gamma-ray source is located in the close vicinity of the Galactic globular cluster Terzan 5 and extends beyond the H.E.S.S. point spread function (0.07 degree). The probability of a chance coincidence with Terzan 5 and an unrelated VHE source is quite low (~ 10^-4). With the largest population of identified millisecond pulsars (msPSRs), a very high core stellar density and the brightest GeV range flux as measured by Fermi-LAT, Terzan 5 stands out among Galactic globular clusters. The properties of the VHE source are briefly discussed in the context of potential emission mechanisms, notably in relation to msPSRs. Interpretation of the available data accommodates several possible origins for this VHE gamma-ray source, although none of them offers a satisfying explanation of its peculiar morphology.

The recent discovery of the radio shell-type supernova remnant (SNR), G353.6-0.7, in spatial coincidence with the unidentified TeV source HESS J1731-347 has motivated further observations of the source with the High Energy Stereoscopic System (H.E.S.S.) Cherenkov telescope array to test a possible association of the gamma-ray emission with the SNR. With a total of 59 hours of observation, representing about four times the initial exposure available in the discovery paper of HESS J1731-347, the gamma-ray morphology is investigated and compared with the radio morphology. An estimate of the distance is derived by comparing the interstellar absorption derived from X-rays and the one obtained from 12CO and HI observations. The deeper gamma-ray observation of the source has revealed a large shell-type structure with similar position and extension (r~0.25{\deg}) as the radio SNR, thus confirming their association. By accounting for the H.E.S.S. angular resolution and projection effects within a simple shell model, the radial profile is compatible with a thin, spatially unresolved, rim. Together with RX J1713.7-3946, RX J0852.0-4622 and SN 1006, HESS J1731-347 is now the fourth SNR with a significant shell morphology at TeV energies. The derived lower limit on the distance of the SNR of 3.2 kpc is used together with radio and X-ray data to discuss the possible origin of the gamma-ray emission, either via inverse Compton scattering of electrons or the decay of neutral pions resulting from proton-proton interaction.

We present the results of observations of the TeV binary LS I +61 303 with the VERITAS telescope array between 2008 and 2010, at energies above 300 GeV. In the past, both ground-based gamma-ray telescopes VERITAS and MAGIC have reported detections of TeV emission near the apastron phases of the binary orbit. The observations presented here show no strong evidence for TeV emission during these orbital phases; however, during observations taken in late 2010, significant emission was detected from the source close to the phase of superior conjunction (much closer to periastron passage) at a 5.6 standard deviation (5.6 sigma) post-trials significance. In total, between October 2008 and December 2010 a total exposure of 64.5 hours was accumulated with VERITAS on LS I +61 303, resulting in an excess at the 3.3 sigma significance level for constant emission over the entire integrated dataset. The flux upper limits derived for emission during the previously reliably active TeV phases (i.e. close to apastron) are less than 5% of the Crab Nebula flux in the same energy range. This result stands in apparent contrast to previous observations by both MAGIC and VERITAS which detected the source during these phases at >10% of the Crab Nebula flux. During the two year span of observations, a large amount of X-ray data were also accrued on LS I +61 303 by the Swift X-ray Telescope (XRT) and the Rossi X-ray Timing Explorer Timing (RXTE) Proportional Counter Array (PCA). We find no evidence for a correlation between emission in the X-ray and TeV regimes during 20 directly overlapping observations. We also comment on data obtained contemporaneously by the Fermi Large Area Telescope (LAT).

Observations of the globular clusters NGC 6388 and M 15 were carried out by the H.E.S.S. array of Cherenkov telescopes for a live time of 27.2 and 15.2 hours respectively. No gamma-ray signal is found at the nominal target position of NGC 6388 and M 15. In the primordial formation scenario, globular clusters are formed in a dark matter halo and dark matter could still be present in the baryon-dominated environment of globular clusters. This opens the possibility of observing a dark matter self-annihilation signal. The dark matter content of the globular clusters NGC 6388 and M 15 is modelled taking into account the astrophysical processes that can be expected to influence the dark matter distribution during the evolution of the globular cluster: the adiabatic contraction of dark matter by baryons, the adiabatic growth of a black hole in the dark matter halo and the kinetic heating of dark matter by stars. 95% confidence level exclusion limits on the dark matter particle velocity-weighted annihilation cross section are derived for these dark matter haloes. In the TeV range, the limits on the velocity-weighted annihilation cross section are derived at the 10-25 cm3 s-1 level and a few 10-24 cm3 s-1 for NGC 6388 and M 15 respectively.

A search for a very-high-energy (VHE; >= 100 GeV) gamma-ray signal from self-annihilating particle Dark Matter (DM) is performed towards a region of projected distance r ~ 45-150 pc from the Galactic Center. The background-subtracted gamma-ray spectrum measured with the High Energy Stereoscopic System (H.E.S.S.) gamma-ray instrument in the energy range between 300 GeV and 30 TeV shows no hint of a residual gamma-ray flux. Assuming conventional Navarro-Frenk-White (NFW) and Einasto density profiles, limits are derived on the velocity-weighted annihilation cross section < \sigma v> as a function of the DM particle mass. These are among the best reported so far for this energy range. In particular, for the DM particle mass of ~1 TeV, values for <\sigma v> above 3 * 10^(-25) cm^3 s^(-1) are excluded for the Einasto density profile. The limits derived here differ much less for the chosen density profile parametrizations, as opposed to limits from gamma-ray observations of dwarf galaxies or the very center of the Milky Way, where the discrepancy is significantly larger.

We report on a newly detected point-like source, HESS J1943+213 located in the Galactic plane. This source coincides with an unidentified hard X-ray source IGR J19443+2117, which was proposed to have radio and infrared counterparts. HESS J1943+213 is detected at the significance level of 7.9 \sigma (post-trials) at RA(J2000)=19h 43m 55s +- 1s (stat) +- 1s (sys), DEC(J2000) = +21deg 18' 8" +- 17" (stat) +- 20" (sys). The source has a soft spectrum with photon index Gamma = 3.1 +- 0.3 (stat) +- 0.2 (sys) and a flux above 470 GeV of 1.3 +- 0.2 (stat) +- 0.3 (sys) x 10^{-12} cm^{-2} s^{-1}. There is no Fermi/LAT counterpart down to a flux limit of 6 x 10^{-9} cm^{-2} s^{-1} in the 0.1-100 GeV energy range (95% confidence upper limit calculated for an assumed power-law model with a photon index Gamma=2.0). The data from radio to VHE gamma-rays do not show any significant variability. We combine new H.E.S.S., Fermi/LAT and Nancay Radio Telescope observations with pre-existing non-simultaneous multi-wavelength observations of IGR J19443+2117 and discuss the likely source associations as well as the interpretation as an active galactic nucleus, a gamma-ray binary or a pulsar wind nebula. The lack of a massive stellar counterpart disfavors the binary hypothesis, while the soft VHE spectrum would be very unusual in case of a pulsar wind nebula. In addition, the distance estimates for Galactic counterparts places them outside of the Milky Way. All available observations favor an interpretation as an extreme, high-frequency peaked BL Lac object with a redshift z>0.14. This would be the first time a blazar is detected serendipitously from ground-based VHE observations, and the first VHE AGN detected in the Galactic Plane.

The gamma-ray pulsar PSR B1706-44 and the adjacent supernova remnant (SNR) candidate G343.1-2.3 were observed by H.E.S.S. during a dedicated observation campaign in 2007. As a result of this observation campaign, a new source of very-high-energy (VHE; E > 100 GeV) gamma-ray emission, HESS J1708-443, was detected with a statistical significance of 7 sigma, although no significant point-like emission was detected at the position of the energetic pulsar itself. In this paper, the morphological and spectral analyses of the newly-discovered TeV source are presented. The centroid of HESS J1708-443 is considerably offset from the pulsar and located near the apparent center of the SNR, at RA(J2000) = 17h08m11s +/- 17s and Dec(J2000) = -44d20' +/- 4'. The source is found to be significantly more extended than the H.E.S.S. point spread function (~0.1 deg), with an intrinsic Gaussian width of 0.29 deg +/- 0.04 deg. Its integral flux between 1 and 10 TeV is ~ 3.8 x 10^-12 ph cm^-2 s^-1, equivalent to 17% of the Crab Nebula flux in the same energy range. The measured energy spectrum is well-fit by a power law with a relatively hard photon index Gamma = 2.0 +/- 0.1 (stat) +/- 0.2 (sys). Additional multi-wavelength data, including 330 MHz VLA observations, were used to investigate the VHE gamma-ray source's possible associations with the pulsar wind nebula of PSR B1706-44 and/or with the complex radio structure of the partial shell-type SNR G343.1-2.3.

Several models of Quantum Gravity predict Lorentz Symmetry breaking at energy scales approaching the Planck scale (10^{19} GeV). With present photon data from the observations of distant astrophysical sources, it is possible to constrain the Lorentz Symmetry breaking linear term in the standard photon dispersion relations. Gamma-ray Bursts (GRB) and flaring Active Galactic Nuclei (AGN) are complementary to each other for this purpose, since they are observed at different distances in different energy ranges and with different levels of variability. Following a previous publication of the High Energy Stereoscopic System (H.E.S.S.) collaboration, a more sensitive event-by-event method consisting of a likelihood fit is applied to PKS 2155-304 flare data of MJD 53944 (July 28, 2006) as used in the previous publication. The previous limit on the linear term is improved by a factor of ~3 up to M^{l}_{QG} > 2.1x10^{18} GeV and is currently the best result obtained with blazars. The sensitivity to the quadratic term is lower and provides a limit of M^{q}_{QG} > 6.4x10^10 GeV, which is the best value obtained so far with an AGN and similar to the best limits obtained with GRB.

The Sculptor and Carina Dwarf spheroidal galaxies were observed with the H.E.S.S. Cherenkov telescope array between January 2008 and December 2009. The data sets consist of a total of 11.8 and 14.8 hours of high quality data, respectively. No gamma-ray signal was detected at the nominal positions of these galaxies above 220 GeV and 320 GeV, respectively. Upper limits on the gamma-ray fluxes at 95% C.L. assuming two forms for the spectral energy distribution (a power law shape and one derived from dark matter annihilation) are obtained at the level of 10^-13 to 10^-12 cm^-2s^-1 in the TeV range. Constraints on the velocity weighted dark matter particle annihilation cross section for both Sculptor and Carina dwarf galaxies range from <\sigma v> ~ 10^-21 cm^3s^-1 down to <\sigma v> ~ 10^-22 cm^3s^-1 depending on the dark matter halo model used. Possible enhancements of the gamma-ray flux are studied: the Sommerfeld effect, which is found to exclude some dark matter particle masses, the internal Bremsstrahlung and clumps in the dark-matter halo distributions.

Context: The detection of gamma-rays in the very-high-energy (VHE) range (100 GeV-100 TeV) offers the possibility of studying the parent population of ultrarelativistic particles found in astrophysical sources, so it is useful for understanding the underlying astrophysical processes in nonthermal sources. Aim: The discovery of the VHE gamma-ray source HESS J1507-622 is reported and possibilities regarding its nature are investigated. Methods: The H.E.S.S. array of imaging atmospheric Cherenkov telescopes (IACTs) has a high sensitivity compared with previous instruments (~1% of the Crab flux in 25 hours observation time for a 5 sigma point-source detection) and has a large field of view (~5 deg in diameter). HESS J1507-622 was discovered within the ongoing H.E.S.S. survey of the inner Galaxy, and the source was also studied by means of dedicated multiwavelength observations. Results: A Galactic gamma-ray source, HESS J1507-622, located ~3.5 deg from the Galactic plane was detected with a statistical significance > 9 sigma. Its energy spectrum is well fitted by a power law with spectral index \Gamma = 2.24 +/- 0.16_{stat} +/- 0.20_{sys} and a flux above 1 TeV of (1.5 +/- 0.4_{stat} +/- 0.3_{sys}) X 10^{-12} cm^{-2} s^{-1}. Possible interpretations (considering both hadronic and leptonic models) of the VHE gamma-ray emission are discussed in the absence of an obvious counterpart.

Aims. Previous observations with the H.E.S.S. telescope array revealed the existence of extended very-high-energy (VHE; E>100 GeV) {\gamma}-ray emission, HESS J1023-575, coincident with the young stellar cluster Westerlund 2. At the time of discovery, the origin of the observed emission was not unambiguously identified, and follow-up observations have been performed to further investigate the nature of this {\gamma}-ray source. Methods. The Carina region towards the open cluster Westerlund 2 has been re-observed, increasing the total exposure to 45.9 h. The combined dataset includes 33 h of new data and now permits a search for energy-dependent morphology and detailed spectroscopy. Results. A new, hard spectrum VHE {\gamma}-ray source, HESSJ1026-582, was discovered with a statistical significance of 7{\sigma}. It is positionally coincident with the Fermi LAT pulsar PSR J1028-5819. The positional coincidence and radio/{\gamma}-ray characteristics of the LAT pulsar favors a scenario where the TeV emission originates from a pulsar wind nebula. The nature of HESS J1023-575 is discussed in light of the deep H.E.S.S. observations and recent multi-wavelength discoveries, including the Fermi LAT pulsar PSRJ1022-5746 and giant molecular clouds in the region. Despite the improved VHE dataset, a clear identification of the object responsible for the VHE emission from HESS J1023-575 is not yet possible, and contribution from the nearby high-energy pulsar and/or the open cluster remains a possibility.

PKS 0548-322 (z=0.069) is a ``high-frequency-peaked'' BL Lac object and a candidate very high energy (VHE, E>100 GeV) gamma-ray emitter, due to its high X-ray and radio flux. Observations at the VHE band provide insights into the origin of very energetic particles present in this source and the radiation processes at work. We report observations made between October 2004 and January 2008 with the H.E.S.S. array, a four imaging atmospheric-Cherenkov telescopes. Contemporaneous UV and X-ray observations with the Swift satellite in November 2006 are also reported. PKS 0548-322 is detected for the first time in the VHE band with H.E.S.S. We measure an excess of 216 gamma-rays corresponding to a significance of 5.6 standard deviations. The photon spectrum of the source is described by a power-law, with a photon index of Gamma=2.86 +/- 0.34 (stat) +/- 0.10 (sys). The integral flux above 200 GeV is 1.3 % of the flux of the Crab Nebula, and is consistent with being constant in time. Contemporaneous Swift/XRT observations reveal an X-ray flux between 2 and 10 keV of F_{2-10 keV}=2.3 +/- 0.2 x 10^{-11} erg.cm^{-2}. s^{-1}, an intermediate intensity state with respect to previous observations. The spectral energy distribution can be reproduced using a simple one-zone synchrotron self Compton model, with parameters similar those observed for other sources of this type.

Observations of very high energy gamma-rays from blazars provide information about acceleration mechanisms occurring in their innermost regions. Studies of variability in these objects allow a better understanding of the mechanisms at play. To investigate the spectral and temporal variability of VHE (>100 GeV) gamma-rays of the well-known high-frequency-peaked BL Lac object PKS 2155-304 with the H.E.S.S. imaging atmospheric Cherenkov telescopes over a wide range of flux states. Data collected from 2005 to 2007 are analyzed. Spectra are derived on time scales ranging from 3 years to 4 minutes. Light curve variability is studied through doubling timescales and structure functions, and is compared with red noise process simulations. The source is found to be in a low state from 2005 to 2007, except for a set of exceptional flares which occurred in July 2006. The quiescent state of the source is characterized by an associated mean flux level of 4.32 +/-0.09 x 10^-11 cm^-2 s^-1 above 200 GeV, or approximately 15% of the Crab Nebula, and a power law photon index of 3.53 +/-0.06. During the flares of July 2006, doubling timescales of ~2 min are found. The spectral index variation is examined over two orders of magnitude in flux, yielding different behaviour at low and high fluxes,which is a new phenomenon in VHE gamma-ray emitting blazars. The variability amplitude characterized by the fractional r.m.s. is strongly energy-dependent and is proportional to E^(0.19 +/- 0.01). The light curve r.m.s. correlates with the flux. This is the signature of a multiplicative process which can be accounted for as a red noise with a Fourier index of ~2. This unique data set shows evidence for a low level gamma-ray emission state from PKS 2155-304, which possibly has a different origin than the outbursts. The discovery of the light curve lognormal behaviour might be an indicator ..

AIMS: The properties of the broad-band emission from the high-frequency peaked BL Lac H 2356-309 (z=0.165) are investigated. METHODS: Very High Energy (VHE; E > 100 GeV) observations of H 2356-309 were performed with the High Energy Stereoscopic System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM optical monitoring observations were also made in 2007. RESULTS: A strong VHE signal, ~13 sigma total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The integral flux above 240 GeV is I(>240 GeV) = (3.06 +- 0.26 {stat} +- 0.61 {syst}) x 10^{-12} cm^{-2} s^{-1}, corresponding to ~1.6% of the flux observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a power law (photon index of Gamma = 3.06 +- 0.15 {stat} +- 0.10 {syst}). Significant small-amplitude variations in the VHE flux from H 2356-309 are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by a hard, broken-power-law spectrum on both nights. CONCLUSIONS: The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed X-ray flux is at a historically low level.

Recent theoretical predictions of the lowest very high energy (VHE) luminosity of SN 1006 are only a factor 5 below the previously published H.E.S.S. upper limit, thus motivating further in-depth observations of this source. Deep observations at VHE energies (above 100 GeV) were carried out with the High Energy Stereoscopic System (H.E.S.S.) of Cherenkov Telescopes from 2003 to 2008. More than 100 hours of data have been collected and subjected to an improved analysis procedure. Observations resulted in the detection of VHE gamma-rays from SN 1006. The measured gamma-ray spectrum is compatible with a power-law, the flux is of the order of 1% of that detected from the Crab Nebula, and is thus consistent with the previously established H.E.S.S. upper limit. The source exhibits a bipolar morphology, which is strongly correlated with non-thermal X-rays. Because the thickness of the VHE-shell is compatible with emission from a thin rim, particle acceleration in shock waves is likely to be the origin of the gamma-ray signal. The measured flux level can be accounted for by inverse Compton emission, but a mixed scenario that includes leptonic and hadronic components and takes into account the ambient matter density inferred from observations also leads to a satisfactory description of the multi-wavelength spectrum.

The full Klein-Nishina cross section for the inverse Compton scattering interactions of electrons implies a significant reduction of the electron energy loss rate compared to the Thomson limit when the electron energy exceeds the critical Klein-Nishina energy E_K = gamma_K * m_e c^2 = 0.27 m_e^2 c^2/(k_BT), where T denotes the temperature of the photon graybody distribution. As a consequence the total radiative energy loss rate of single electrons exhibits sudden drops in the overall energy loss rate ~ gamma ^2 - dependence when the electron energy reaches the critical Klein-Nishina energy. The strength of the drop is proportional to the energy density of the photon radiation field. The diffuse galactic optical photon fields from stars of spectral type B and G-K lead to critical Klein-Nishina energies of 40 and 161 GeV, respectively. Associated with the drop in the loss rate are sudden increases (Klein-Nishina steps) in the equilibrium spectrum of cosmic ray electrons. Because the radiative loss rate of electrons is the main ingredient in any transport model of high-energy cosmic ray electrons, Klein-Nishina steps will modify any calculated electron equilibrium spectrum irrespective of the electron sources and spatial transport mode. To delineate most clearly the consequences of the Klein-Nishina drops in the radiative loss rate, we chose as illustrative example the simplest realistic model for cosmic ray electron dynamics in the Galaxy, consisting of the competition of radiative losses and secondary production by inelastic hadron-hadron collisions. We demonstrate that the spectral structure in the FERMI and H.E.S.S. data is well described and even the excess measured by ATIC might be explained by Klein-Nishina steps.