Observations by the Fermi Large Area Telescope of gamma-ray millisecond pulsar light curves imply copious pair production in their magnetospheres, and not exclusively in those of younger pulsars. Such pair cascades may be a primary source of Galactic electrons and positrons, contributing to the observed enhancement in positron flux above ~10 GeV. Fermi has also uncovered many new millisecond pulsars, impacting Galactic stellar population models. We investigate the contribution of Galactic millisecond pulsars to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day millisecond pulsars. We simulate their pair spectra invoking an offset-dipole magnetic field. We also consider positrons and electrons that have been further accelerated to energies of several TeV by strong intrabinary shocks in black widow and redback systems. Since millisecond pulsars are not surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs freely escape and undergo losses only in the intergalactic medium. We compute the transported pair spectra at Earth, following their diffusion and energy loss through the Galaxy. The predicted particle flux increases for non-zero offsets of the magnetic polar caps. Pair cascades from the magnetospheres of millisecond pulsars are only modest contributors around a few tens of GeV to the lepton fluxes measured by AMS-02, PAMELA, and Fermi, after which this component cuts off. The contribution by black widows and redbacks may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.

Pulsars are believed to be sources of relativistic electrons and positrons. The abundance of detections of gamma-ray millisecond pulsars by Fermi Large Area Telescope coupled with their light curve characteristics that imply copious pair production in their magnetospheres, motivated us to investigate this old pulsar population as a source of Galactic electrons and positrons and their contribution to the enhancement in cosmic-ray positron flux at GeV energies. We use a population synthesis code to predict the source properties (number, position, and power) of the present-day Galactic millisecond pulsars, taking into account the latest Fermi and radio observations to calibrate the model output. Next, we simulate pair cascade spectra from these pulsars using a model that invokes an offset-dipole magnetic field. We assume free escape of the pairs from the pulsar environment. We then compute the cumulative spectrum of transported electrons and positrons at Earth, following their diffusion and energy losses as they propagate through the Galaxy. Our results indicate that the predicted particle flux increases for non-zero offsets of the magnetic polar caps. Comparing our predicted spectrum and positron fraction to measurements by AMS-02, PAMELA, and Fermi, we find that millisecond pulsars are only modest contributors at a few tens of GeV, after which this leptonic spectral component cuts off. The positron fraction is therefore only slightly enhanced above 10 GeV relative to a background flux model. This implies that alternative sources such as young, nearby pulsars and supernova remnants should contribute additional primary positrons within the astrophysical scenario.

Recent observations by the Fermi Large Area Telescope (LAT) and the High Energy Stereoscopic System (H.E.S.S.) have revealed globular clusters (GC) to be sources of high-energy (HE) and very-high-energy (VHE) gamma-rays. It has been suggested that the presence of large numbers of millisecond pulsars (MSPs) within these clusters may be either directly responsible for these gamma-ray fluxes through emission of pulsed curvature radiation, or indirectly through the injection of relativistic leptons into the cluster. These relativistic particles are plausibly re-accelerated in shocks, created by the collision of stellar winds, before interacting with the soft-photon radiation field set up by the stellar population of the host cluster. Inverse Compton (IC) scattering then produces gamma-radiation in the TeV band. In order to calculate the IC spectrum, an accurate profile for the energy density of the soft-photon field is required. We construct such a profile by deriving a radially-dependent expression for the stellar energy density, and then solve it numerically. As a next step, the average energy density values for three different regions of the cluster (demarcated by its core, half-mass, and tidal radii) are determined, which we consequently import into an existing radiation code to predict the TeV gamma-ray spectrum. As an application, we consider the case of Terzan 5, boasting a population of 34 radio MSPs, and compare our predicted spectrum with that measured by H.E.S.S.

Given that the bulk of the galactic cosmic rays (CRs) originates from transient point-like sources, such as supernova remnants, the flux of the CR primary component measured at Earth depends on the local source history. Whereas the secondary omponent shows little to no variation due to nearby sources, and instead, should depend on the global distribution of CRs. Steady-state, rotational symmetric models (2D) of CR propagation assume smeared-out CR sources in the Galaxy and cannot properly take into account the influence of nearby point sources. In this paper, we infer evidence of nearby sources by illustrating that a 2D propagation model will not describe CR primaries and secondaries equally well. We adapted the 2D version of the GALPROP code to a compute-cluster environment and perform parameter studies comparing CR spectra with primary and secondary CR data separately. Doing a parameter study, one may expect different best fit values looking at the primary and secondary CR components separately, as it is unlikely that the source history mimicked by the 2D models coincides with the real local source history. We find that the primaries and secondaries are fitted with differing best fit parameters, with the results indicating that the 2D model being more suited to modelling secondary CRs. The insufficient modelling of primary CRs can be contributed to the possible presence of local CR sources.

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

Galactic transport models for cosmic rays involve the diffusive motion of these particles in the interstellar medium. Due to the large-scale structured galactic magnetic field this diffusion is anisotropic with respect to the local field direction. We included this transport effect along with continuous loss processes in a quantitative model of galactic propagation for cosmic ray protons which is based on stochastic differential equations. We calculated energy spectra at different positions along the Sun's galactic orbit and compared them to the isotropic diffusion case. The results show that a larger amplitude of variation as well as different spectral shapes are obtained in the introduced anisotropic diffusion scenario and emphasize the need for accurate galactic magnetic field models.

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.

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.

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

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.

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.

The Fermi Large Area Telescope (LAT) has recently detected a population of globular clusters (GCs) in high-energy (HE) gamma-rays. Their spectral properties and energetics are consistent with cumulative emission from a population of millisecond pulsars (MSPs) hosted by these clusters. For example, the HE spectra exhibit fairly hard power-law indices and cutoffs around a few GeV, typical of pulsed spectra measured for the gamma-ray pulsar population. The energetics may be used to constrain the number of visible MSPs in the cluster (N_vis), assuming canonical values for the average gamma-ray efficiency and spin-down power. This interpretation is indeed strengthened by the fact that the first gamma-ray MSP has now been identified in the GC NGC 6624, and this MSP is responsible for almost all of the HE emission from this cluster. On the other hand, it has been argued that the MSPs are also sources of relativistic leptons which may be reaccelerated in shocks originating in collisions of stellar winds in the cluster core, and may upscatter bright starlight and cosmic microwave background photons to very high energies. Therefore, this unpulsed component may give an independent constraint on the total number of MSPs (N_tot) hosted in the GC, for a given cluster magnetic field B and diffusion coefficient k_0. Lastly, the transport properties of the energetic leptons may be further constrained using multiwavelength data, e.g., to infer the radial dependence of k_0 and B. We present results on our modeling of the pulsed and unpulsed gamma-ray fluxes from the GC Terzan 5.

Recent data reported by the PAMELA and ATIC experiments, as well as cosmic ray lepton results from FERMI and H.E.S.S. collaborations sparked a series of papers explaining these results either by contributions of electron positron pairs to the local interstellar cosmic ray (CR) spectrum by dark matter (DM) or pulsars. Focusing here on pulsars, we argue that at the present, our knowledge about particle acceleration at pulsars as well as of the local Galactic CR propagation is still limited, i.e. the recent results for CR electrons and positrons constrain pulsar and propagation models. We will thus not present another attempt to explain the data by contributions of pulsars to the local CR lepton flux but rather to highlight the caveats in doing so.

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)

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.

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.

Aims: To provide a significantly improved probability distribution for the H-test for periodicity in X-ray and $\gamma$-ray arrival times, which is already extensively used by the $\gamma$-ray pulsar community. Also, to obtain an analytical probability distribution for stacked test statistics in the case of a search for pulsed emission from an ensemble of pulsars where the significance per pulsar is relatively low, making individual detections insignificant on their own. This information is timely given the recent rapid discovery of new pulsars with the Fermi-LAT t $\gamma$-ray telescope. Methods: Approximately $10^{14}$ realisations of the H-statistic ($H$) for random (white) noise is calculated from a random number generator for which the repitition cycle is $\gg 10^{14}$. From these numbers the probability distribution $P(>H)$ is calculated. Results: The distribution of $H$ is is found to be exponential with parameter $\lambda=0.4$ so that the cumulative probability distribution $P(>H)=\exp{(-\lambda H)}$. If we stack independent values for $H$, the sum of $K$ such values would follow the Erlang-K distribution with parameter $\lambda$ for which the cumulative probability distribution is also a simple analytical expression. Conclusion: Searches for weak pulsars with unknown pulse profile shapes in the Fermi-LAT, Agile or other X-ray data bases should benefit from the {\it H-test} since it is known to be powerful against a broad range of pulse profiles, which introduces only a single statistical trial if only the {\it H-test} is used. The new probability distribution presented here favours the detection of weaker pulsars in terms of an improved sensitivity relative to the previously known distribution.

We present an X-ray analysis and a model of the nonthermal emission of the pulsar wind nebula (PWN) MSH15-52. We analyzed XMM-Newton data to obtain the spatially resolved spectral parameters around the pulsar PSRB1509-58. A steepening of the fitted power-law spectra and decrease in the surface brightness is observed with increasing distance from the pulsar. In the second part of this paper, we introduce a model for the nonthermal emission, based on assuming the ideal magnetohydrodynamic limit. This model is used to constrain the parameters of the termination shock and the bulk velocity of the leptons in the PWN. Our model is able to reproduce the spatial variation of the X-ray spectra. The parameter ranges that we found agree well with the parameter estimates found by other authors with different approaches. In the last part of this paper, we calculate the inverse Compton emission from our model and compare it to the emission detected with the H.E.S.S. telescope system. Our model is able to reproduce the flux level observed with H.E.S.S., but not the spectral shape of the observed TeV {\gamma}-ray emission.

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 recent discovery of diffuse, VHE gamma radiation from the Galactic center ridge by the H.E.S.S. telescope allow for the first time the direct determination of parameters of galactic cosmic ray propagation models. In this paper we show that the diffuse gamma radiation near the Galactic center may be explained by the interaction of VHE cosmic ray (CR) protons with the interstellar gas located in several giant molecular clouds leading to a measurement of the cosmic ray diffusion coefficient for the galactic center region of k = 1.3 kpc^2 Myr^-1 for a mean proton energy of about 3 TeV, if we assume that the CR protons originated from a supernova event (Sgr A East), which took off about 10 kyr ago. This value of k is about 5 to 10 times smaller than the locally measured value.

Recent EGRET observations of the diffuse Galactic gamma-ray emission reveal a spectrum which is incompatible with the assumption that the cosmic ray spectra measured locally hold throughout the Galaxy: the spectrum above 1 GeV, where the emission is supposedly dominated by pi^0-decay, is harder than that derived from the local cosmic ray proton spectrum. We demonstrate that in case of a SNR origin of cosmic ray nucleons part of this gamma-ray excess may be attributed to the dispersion of the spectral indices in these objects. In global averages, as are gamma-ray line-of-sight integrals, this dispersion leads to a positive curvature in the composite spectrum, and hence to modified pi^0-decay gamma-ray spectra.