The pulsar wind nebula (PWN) associated with the Vela pulsar is a bright source in the radio, X-ray and gamma-ray bands, but not in the optical. This source is very near, lying at a distance of 290 pc, as inferred from the radio and optical parallax measurements of the pulsar. Knowledge of the brightness and structure of the Vela PWN in optical is important in order to constrain the underlying particle spectrum (and possibly the B-field properties and particle losses) associated with this extended source. We use results from the Digital Sky Survey, as well as results obtained using the SAAO 1.0 m telescope equipped with an imaging CCD (STE4) and BV filters, in an attempt to measure optical radiation from Vela X. To enlarge our field of view, we constructed a mosaic consisting of 3 x 3 frames around the pulsar position. We present spectral measurements from the High Energy Stereoscopic System (H.E.S.S.), Fermi Large Area Telescope (LAT), ASCA, Hubble Space Telescope (HST), Very Large Telescope (VLT), New Technology Telescope (NTT), and Wilkinson Microwave Anisotropy Probe (WMAP), in addition to our optical results. Using these data, we investigate whether or not the radio synchrotron component can be smoothly extrapolated to the optical band. This would constrain the electron population to consist of either a single or multiple components, representing a significant advancement in our understanding of this complex multiwavelength source.

Globular clusters (GCs) are astronomical tapestries embroidered with an abundance of exotic stellar-type objects. Their high age promises a rich harvest of evolved stellar products, while the deep potential wells and high mass densities at their centres probably facilitate the formation of multiple-member stellar systems via increased stellar encounter rates. The ubiquity of low-mass X-ray binaries, thought to be the progenitors of millisecond pulsars (MSPs), explain the large number of observed GC radio pulsars and X-ray counterparts. The Fermi Large Area Telescope (LAT) recently unveiled the first gamma-ray GC pulsar (PSR J1823-3021A). The first observations of GCs in the GeV and TeV bands furthermore created much excitement, and in view of the above, it seems natural to explain these high-energy lanterns by investigating an MSP origin. An MSP population is expected to radiate several pulsed spectral components in the radio through gamma-ray wavebands, in addition to being sources of relativistic particles. The latter may interact with background photons in the clusters producing TeV excesses, while they may also radiate synchrotron photons as they traverse the cluster magnetic field. We present multiwavelength modelling results for Terzan 5. We also briefly discuss some alternative interpretations for the observed GC gamma-ray signals.

Composite supernova remnants consist of a pulsar wind nebula located inside a shell-type remnant. The presence of a shell has implications on the evolution of the nebula, although the converse is generally not true. The purpose of this paper is two-fold. The first aim is to determine the effect of the pulsar's initial luminosity and spin-down rate, the supernova ejecta mass, and density of the interstellar medium on the evolution of a spherically-symmetric, composite supernova remnant expanding into a homogeneous medium. The second aim is to investigate the evolution of the magnetic field in the pulsar wind nebula when the the composite remnant expands into a non-uniform interstellar medium. The Euler conservation equations for inviscid flow, together with the magnetohydrodynamic induction law in the kinematic limit, are solved numerically for a number of scenarios where the ratio of magnetic to particle energy is $\sigma < 0.01$. The simulations in the first part of the paper is solved in a one-dimensional configuration. In the second part of the paper, the effect of an inhomogeneous medium on the evolution is studied using a two-dimensional, axis-symmetric configuration.

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.

Diffuse X-ray emission has recently been detected from the globular cluster (GC) Terzan 5, extending out to ~2.5' from the cluster centre. This emission may arise from synchrotron radiation (SR) by energetic leptons being injected into the cluster by the resident millisecond pulsar (MSP) population that interact with the cluster field. These leptons may also be reaccelerated in shocks created by collisions of pulsar winds, and may interact with bright starlight and cosmic microwave background (CMB) photons, yielding gamma rays at very high energies (VHE) through the inverse Compton (IC) process. In the GeV range, Fermi Large Area Telescope (LAT) has detected a population of GCs, very plausibly including Terzan 5, their spectral properties and energetics being consistent with cumulative magnetospheric emission from a population of MSPs. H.E.S.S. has furthermore detected a VHE excess in the direction of Terzan 5. One may derive constraints on the number of MSPs, N_tot, and the radial profiles of the GC B-field, stellar energy density, as well as the diffusion coefficient using the spatially-resolved X-ray, high-energy (HE), and VHE fluxes. If the Fermi LAT flux is due to magnetospheric processes, it will scale with the number of visible gamma-ray MSPs, N_vis. The HE spectrum therefore provides an independent way of constraining the number of MSPs (since N_tot >= N_vis). Consequently, the synthesis of available multiwavelength data presents a unique opportunity to constrain several parameters of the GC Terzan 5.

In a Pulsar Wind Nebula (PWN), the lifetime of inverse Compton (IC) emitting electrons exceeds the lifetime of its progenitor pulsar (as well as its shell-type remnant), but it also exceeds the age of those that emit via synchrotron radiation. Therefore, during its evolution, the PWN can remain bright in IC so that its GeV-TeV gamma-ray flux remains high for timescales much larger (for 10^5 - 10^6 yrs) than the pulsar lifetime and the X-ray PWN lifetime. In this scenario, the magnetic field in the cavity induced by the wind of the progenitor star plays a crucial role. This scenario is in line with the discovery of several unidentified or "dark" sources in the TeV gamma-ray band without X-ray counterparts; and it is also finding confirmation in the recent discoveries at GeV gamma rays. Moreover, these consequences could be also important for reinterpreting the detection of starburst galaxies in the TeV gamma-ray band when considering a leptonic origin of the gamma-ray signal. Both theoretical aspects and their observational proofs will be discussed, as well as the first results of our new modeling code.

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

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.

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

The Fermi Large Area Telescope (LAT) has recently reported the detection of pulsed gamma-rays from 6 young pulsars (J0631+1036, J0659+1414, J0742-2822, J1420-6048, J1509-5850, and J1718-3825), all exhibiting single-peaked pulse profiles (Weltevrede et al., 2010). High-quality radio polarization data are also available for 5 of these pulsars, allowing derivation of constraints on their viewing geometries. We obtain independent constraints on the viewing geometries of these pulsars by using a geometric pulsar emission code to model the Fermi LAT and radio light curves. We find fits for the magnetic inclination and observer angles alpha and zeta with typical errors of ~ 5deg. Our results are generally consistent with those obtained by Weltevrede et al. (2010), although we do find differences in some cases. Our model may lastly provide a framework to constrain the radio emission altitude.

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.

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

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.

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.

Galactic very high energy (VHE, > 100 GeV) gamma ray sources in the inner Galaxy H.E.S.S. survey tend to cluster within 1 degree in latitude around the Galactic plane. HESS J1507-622 instead is unique, since it is located at latitude of ~3.5 degrees. HESS J1507-622 is slightly extended over the PSF of the instrument and hence its Galactic origin is clear. The search for counterparts in other wavelength regimes (radio, infrared and X-rays) failed to show any plausible counterparts; and given its position off the Galactic plane and hence the absorption almost one order of magnitude lower, it is very surprising to not see any counterparts especially at X-rays wavelengths (by ROSAT, XMM Newton and Chandra). Its latitude implies that it is either rather close, within about 1 kpc, or is located well off the Galactic plane. And also the models reflect the uniqueness of this object: a leptonic PWN scenario would place this source due to its quite small extension to multi-kpc distance whereas a hadronic scenario would preferentially locate this object at distances of < 1 kpc where the density of target material is higher.

H.E.S.S. is one of the most sensitive instruments in the very high energy (VHE; > 100 GeV) gamma-ray domain and has revealed many new sources along the Galactic Plane. After the successful first VHE Galactic Plane Survey of 2004, H.E.S.S. has continued and extended that survey in 2005-2008, discovering a number of new sources, many of which are unidentified. Some of the unidentified H.E.S.S. sources have several positional counterparts and hence several different possible scenarios for the origin of the VHE gamma-ray emission; their identification remains unclear. Others have so far no counterparts at any other wavelength. Particularly, the lack of an X-ray counterpart puts serious constraints on emission models. Several newly discovered and still unidentified VHE sources are reported here.

Low-Mass X-ray Binaries (LMXRBs), believed to be the progenitors of recycled millisecond pulsars (MSPs), occur abundantly in globular clusters (GCs). GCs are therefore expected to host large numbers of MSPs. This is also confirmed observationally. The MSPs continuously inject relativistic electrons into the ambient region beyond their light cylinders, and these relativistic particles produce unpulsed radiation via the synchrotron and inverse Compton (IC) processes. It is thus possible, in the context of General Relativistic (GR) frame-dragging MSP models, to predict unpulsed very high energy radiation expected from nearby GCs. We use a period-derivative cleaned sample of MSPs in 47 Tucanae, where the effects of the cluster potential on the individual period derivatives have been removed. Using a Monte Carlo process to include effects of pulsar geometry, we obtain average injection spectra (with relatively small errors) of particles leaving the MSPs. These spectra are next used to predict cumulative synchrotron and IC spectra expected from 47 Tucanae, which is a lower limit, as no reacceleration is assumed. We find that the IC radiation from 47 Tucanae may be visible for H.E.S.S., depending on the nebular field B as well as the number of MSPs N in the GC. Telescopes such as Chandra and Hubble may find it difficult to test the SR component prediction of diffuse radiation if there are many unresolved sources in the field of view. These results may be rescaled for other GCs where less information is available, assuming universal GC MSP characteristics.

The energetic pulsar PSR B1706-44 and the adjacent supernova remnant (SNR) candidate G 343.1-2.3 were observed by H.E.S.S. during a dedicated observational campaign in 2007. A new source of very-high-energy (VHE; E > 100 GeV) gamma-ray emission, HESS J1708-443, was discovered with its centroid at RA(J2000) = 17h08m10s and Dec(J2000) = -44d21', with a statistical error of 3 arcmin on each axis. The VHE gamma-ray source is significantly more extended than the H.E.S.S. point-spread function, with an intrinsic Gaussian width of 0.29 +/- 0.04 deg. Its energy spectrum can be described by a power law with a photon index Gamma = 2.0 +/- 0.1 (stat) +/- 0.2 (sys). The integral flux measured between 1-10 TeV is ~17% of the Crab Nebula flux in the same energy range. The possible associations with PSR B1706-44 and SNR G343.1-2.3 are discussed.

H.E.S.S. is one of the most sensitive instruments in the very high energy (VHE; > 100 GeV) gamma-ray domain and has revealed many new sources along the Galactic Plane. After the successful first VHE Galactic Plane Survey of 2004, H.E.S.S. has continued and extended that survey in 2005-2008, discovering a number of new sources, many of which are unidentified. Some of the unidentified H.E.S.S. sources have several positional counterparts and hence several different possible scenarios for the origin of the VHE gamma-ray emission; their identification remains unclear. Others have so far no counterparts at any other wavelength. Particularly, the lack of an X-ray counterpart puts serious constraints on emission models. Several newly discovered and still unidentified VHE sources are reported here.

In this paper we explore the evolution of a PWN while the pulsar is spinning down. An MHD approach is used to simulate the evolution of a composite remnant. Particular attention is given to the adiabatic loss rate and evolution of the nebular field strength with time. By normalising a two component particle injection spectrum (which can reproduce the radio and X-ray components) at the pulsar wind termination shock to the time dependent spindown power, and keeping track with losses since pulsar/PWN/SNR birth, we show that the average field strength decreases with time as $t^{-1.3}$, so that the synchrotron flux decreases, whereas the IC gamma-ray flux increases, until most of the spindown power has been dumped into the PWN. Eventually adiabatic and IC losses will also terminate the TeV visibility and then eventually the GeV visibility.

The recent Fermi detection of the globular cluster (GC) 47 Tucanae highlighted the importance of modeling collective gamma-ray emission of millisecond pulsars (MSPs) in GCs. Steady flux from such populations is also expected in the very high energy (VHE) domain covered by ground-based Cherenkov telescopes. We present pulsed curvature radiation (CR) as well as unpulsed inverse Compton (IC) calculations for an ensemble of MSPs in the GCs 47 Tucanae and Terzan 5. We demonstrate that the CR from these GCs should be easily detectable for Fermi, while constraints on the total number of MSPs and the nebular B-field may be derived using the IC flux components.

Vela X is a pulsar wind nebula (PWN) associated with the active pulsar B0833-45 and contained within the Vela supernova remnant (SNR). A collimated X-ray filament ("cocoon") extends south-southwest from the pulsar to the center of Vela X. VLA observations uncovered radio emission coincident with the eastern edge of the cocoon and H.E.S.S. has detected TeV $\gamma$-ray emission from this region as well. Using XMM-\textit{Newton} archival data, covering the southern portion of this feature, we analyze the X-ray properties of the cocoon. The X-ray data are best fit by an absorbed nonequilibrium plasma model with a powerlaw component. Our analysis of the thermal emission shows enhanced abundances of O, Ne, and Mg within the cocoon, indicating the presence of ejecta-rich material from the propagation of the SNR reverse shock, consistent with Vela X being a disrupted PWN. We investigate the physical processes that excite the electrons in the PWN to emit in the radio, X-ray and $\gamma$-ray bands. The radio and non-thermal X-ray emission can be explained by synchrotron emission. We model the $\gamma$-ray emission by Inverse Compton scattering of electrons off of cosmic microwave background (CMB) photons. We use a 3-component broken power law to model the synchrotron emission, finding an intrinsic break in the electron spectrum at $\sim5 \times 10^{6}$ keV and a cooling break at $\sim$ 5.5 $\times 10^{10}$ keV. This cooling break along with a magnetic field strength of 5 $\times 10^{-6}$ G indicate that the synchrotron break occurs at $\sim$1 keV.

Morphologically it appears as if the Vela X PWN consists of two emission regions: whereas X-ray (1 keV) and very high energy (VHE) H.E.S.S. gamma-ray observations appear to define a cocoon type shape south of the pulsar, radio observations reveal an extended area of size 2 deg by 3 deg (including the cocoon area), also south of the Vela pulsar. Since no wide field of view (FoV) observations of the synchrotron emission between radio and X-rays are available, we do not know how the lepton (e+/-) spectra of these two components connect and how the morphology changes with energy. Currently we find that two distinct lepton spectra describe the respective radio and X-ray/VHE gamma-ray spectra, with a field strength of 5 muG self-consistently describing a radiation spectral break (or energy maximum) in the multi-TeV domain as observed by H.E.S.S. (if interpreted as IC radiation), while predicting the total hard X-ray flux above 20 keV (measured by the wide FoV INTEGRAL instrument) within a factor of two. If this same field strength is also representative of the radio structure (including filaments), the implied IC component corresponding to the highest radio frequencies should reveal a relatively bright high energy gamma-ray structure and Fermi LAT should be able to resolve it. A higher field strength in the filaments would however imply fewer leptons in Vela X and hence a fainter Fermi LAT signal.

Although only 22 millisecond pulsars (MSPs) are currently known to exist in the globular cluster (GC) 47 Tucanae, this cluster may harbor 30-60 MSPs, or even up to ~200. In this Letter, we model the pulsed curvature radiation (CR) gamma-ray flux expected from a population of MSPs in 47 Tucanae. These MSPs produce gamma-rays in their magnetospheres via accelerated electron primaries which are moving along curved magnetic field lines. A GC like 47 Tucanae containing a large number of MSPs provides the opportunity to study a randomized set of pulsar geometries. Geometry-averaged spectra make the testing of the underlying pulsar model more reliable, since in this case the relative flux uncertainty is reduced by one order of magnitude relative to the variation expected for individual pulsars (if the number of visible pulsars N=100). Our predicted spectra violate the EGRET upper limit at 1 GeV, constraining the product of the number of visible pulsars N and the average integral flux above 1 GeV per pulsar. GLAST/LAT should place even more stringent constraints on this product, and may also limit the maximum average accelerating potential by probing the CR spectral tail. For N=22-200, a GLAST/LAT non-detection will lead to the constraints that the average integral flux per pulsar should be lower by factors 0.03-0.003 than current model predictions.

Geminga and B0656+14 are the closest pulsars with characteristic ages in the ran ge of 100 kyr to 1 Myr. They both have spindown powers of the order 3e34 erg/s at present. The winds of these pulsars had most probably powered pulsar wind nebulae (PWNe) that broke up less than about 100 kyr after the birth of the pulsars. Assuming that leptonic particles accelerated by the pulsars were confined in th e PWNe and were released into the interstellar medium (ISM) on breakup of the PW Ne, we show that, depending on the pulsar parameters, both pulsars make a non-ne gligible contribution to the local cosmic ray (CR) positron spectrum, and they m ay be the main contributors above several GeV. The relatively small angular dist ance between Geminga and B0656+14 thus implies an anisotropy in the local CR po sitron flux at these energies. We calculate the contribution of these pulsars to the locally observed CR electr on and positron spectra depending on the pulsar birth period and the magnitude o f the local CR diffusion coefficient. We further give an estimate of the expecte d anisotropy in the local CR positron flux. Our calculations show that within the framework of our model, the local CR posit ron spectrum imposes constraints on pulsar parameters for Geminga and B0656+14, notably the pulsar period at birth, and also the local interstellar diffusion co efficient for CR leptons.

In this paper we show that the high energy $\gamma$-ray flux in the GeV domain from mature pulsar wind nebulae (PWN) scales as the change in rotational kinetic energy $I(\Omega_0^2-\Omega^2)/2$ since birth, rather than the present day spindown power $I\Omega\dot{\Omega}$. This finding holds as long as the lifetime of inverse Compton emitting electrons exceeds the age of the system. For a typical $\gamma^{-2}$ electron spectrum, the predicted flux depends mostly on the pulsar birth period, conversion efficiency of spindown power to relativistic electrons and distance to the PWN, so that first order estimates of the birth period can be assessed from {\it GLAST/LAT} observations of PWN. For this purpose we derive an analytical expression. The associated (``uncooled'') photon spectral index in the GeV domain is expected to cluster around $\sim 1.5$, which is bounded at low energies by an intrinsic spectral break, and at higher energies by a second spectral break where the photon index steepens to $\sim 2$ due to radiation losses. Mature PWN are expected to have expanded to sizes larger than currently known PWN, resulting in relatively low magnetic energy densities and hence survival of GeV inverse Compton emitting electrons. Whereas such a PWN may be radio and X-ray quiet in synchrotron radiation, it may still be detectable as a {\it GLAST/LAT} source as a result of the relic electrons in the PWN.

In this review paper on pulsar wind nebulae (PWN) we discuss the properties of such nebulae within the context of containment against cross-field diffusion (versus normal advection), the effect of reverse shocks on the evolution of offset ``Vela-like'' PWN, constraints on maximum particle energetics, magnetic field strength estimates based on spectral and spatial properties, and the implication of such field estimates on the composition of the wind. A significant part of the discussion is based on the High Energy Stereoscopic System ({\it H.E.S.S.} or {\it HESS}) detection of the two evolved pulsar wind nebulae Vela X (cocoon) and HESS J1825-137. In the case of Vela X (cocoon) we also review evidence of a hadronic versus a leptonic interpretation, showing that a leptonic interpretation is favored for the {\it HESS} signal. The constraints discussed in this review paper sets a general framework for the interpretation of a number of offset, filled-center nebulae seen by {\it HESS}. These sources are found along the galactic plane with galactic latitudes $|b|\sim 0$, where significant amounts of molecular gas is found. In these regions, we find that the interstellar medium is inhomogeneous, which has an effect on the morphology of supernova shock expansion. One consequence of this effect is the formation of offset pulsar wind nebulae as observed.

The deeper and more extended survey of the central parts of the Galactic Plane by H.E.S.S. during 2005-2007 has revealed a number of new point-like, as well as, extended sources. Two point-like sources can be associated to two remarkable objects around "Crab-like" young and energetic pulsars in our Galaxy: G21.5-0.9 and Kes 75. The characteristics of each of the sources are presented and possible interpretations are briefly discussed.

Combined X-ray synchrotron and inverse-Compton gamma-ray observations of pulsar wind nebulae (PWN) may help to elucidate the processes of acceleration and energy loss in these systems. In particular, such observations provide constraints on the particle injection history and the magnetic field strength in these objects. The newly discovered TeV gamma-ray source HESS J1718-385 has been proposed as the likely PWN of the high spin-down luminosity pulsar PSR J1718-3825. The absence of previous sensitive X-ray measurements of this pulsar, and the unusual energy spectrum of the TeV source, motivated observations of this region with XMM-Newton. The data obtained reveal a hard spectrum X-ray source at the position of PSR 1718-3825 and evidence for diffuse emission in the vicinity of the pulsar. We derive limits on the keV emission from the centroid of HESS J1718-385 and discuss the implications of these findings for the PWN nature of this object.

Since 2002 the VHE (>100 GeV) gamma-ray flux of the high-frequency peaked BL Lac PKS 2155-304 has been monitored with the High Energy Stereoscopic System (HESS). An extreme gamma-ray outburst was detected in the early hours of July 28, 2006 (MJD 53944). The average flux above 200 GeV observed during this outburst is ~7 times the flux observed from the Crab Nebula above the same threshold. Peak fluxes are measured with one-minute time scale resolution at more than twice this average value. Variability is seen up to ~600 s in the Fourier power spectrum, and well-resolved bursts varying on time scales of ~200 seconds are observed. There are no strong indications for spectral variability within the data. Assuming the emission region has a size comparable to the Schwarzschild radius of a ~10^9 solar mass black hole, Doppler factors greater than 100 are required to accommodate the observed variability time scales.

Radio, X-ray, and HESS gamma-ray observations of the Galactic Center (GC) composite supernova remnant SNR G0.9+0.1 are used to constrain a time-dependent injection model of the downstream electron spectrum responsible for the total multiwavelength spectrum. The effect of spindown power evolution aswell as nebular field evolution is employed to reproduce the present-day multiwavelength spectrum. Assuming a nebular magnetic field decay model of typical HESS-type pulsar wind nebulae (PWN), ending with a present-day field strength of 6muG, we obtain an initial spindown power of ~ 10E+38 ergs/s if we assume a birth period and age of 43ms and 6,500 yr respectively to reproduce the properties of the SNR shell. This gives a present-day spindown power of ~ 10E+37 ergs/s, which agrees well with the present-day spindown power derived from X-ray observations.

We present the highly significant detection of a quasi-periodic flux modulation with a period of 22.2 min seen in the X-ray data of the Sgr A* flare of 2004 August 31. This flaring event, which lasted a total of about three hours, was detected simultaneously by EPIC on XMM-Newton and the NICMOS near-infrared camera on the HST. Given the inherent difficulty in, and the lack of readily available methods for quantifying the probability of a periodic signal detected over only several cycles in a data set where red noise can be important, we developed a general method for quantifying the likelihood that such a modulation is indeed intrinsic to the source and does not arise from background fluctuations. We here describe this Monte Carlo based method, and discuss the results obtained by its application to a other XMM-Newton data sets. Under the simplest hypothesis that we witnessed a transient event that evolved, peaked and decayed near the marginally stable orbit of the supermassive black hole, this result implies that for a mass of 3.5 x 10^{6} Msun, the central object must have an angular momentum corresponding to a spin parameter of a=0.22.

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.

The high resolution capabilities of the High Energy Stereoscopic System (HESS) introduced a new era in Gamma-Ray Astronomy, and opens a new window on pulsar wind nebula (PWN) research. A rotationally induced jet (associated with PSR B1509-58) is resolved for the first time in gamma-rays, allowing us to trace the particle transport directly, without having the complicating effect of spatially varying field distributions on the synchrotron emissivity. For PWN older or more extended than Crab (i.e. those with lower field strengths), HESS also reveals the properties of electrons contributing to the EUV/soft X-ray synchrotron bands, whereas EUV/soft X-rays suffer from severe interstellar absorption effects. Finally, HESS morphological studies of evovled PWN also allow us to directly measure the effects of assymetric reverse shock interactions due to SNR forward shock expansion into the inhomogeneous interstellar medium.

The excellent sensitivity and high resolution capability of wide FoV ground-based imaging atmospheric Cerenkov telescopes allow us for the first time to resolve the morphological structures of pulsar wind nebulae (PWN) which are older and more extended than the Crab Nebula. VHE gamma-ray observations of such extended nebulae (with field strengths below ~ 20 micro Gauss) probe the electron component corresponding to the unseen extreme ultraviolet (EUV) synchrotron component, which measures electron injection from earlier evolutionary epochs. VHE observations of PWN therefore introduce a new window on PWN research. This review paper also identifies conditions for maximal VHE visbility of PWN. Regarding pulsar pulsed emission, it is becoming clear that the threshold energies of current telescopes are not sufficient to probe the pulsed gamma-ray component from canonical pulsars. Theoretical estimates of pulsed gamma-ray emission from millisecond pulsars seem to converge and it becomes clear that such detections with current 3rd generation telescopes will not be possible, unless the geometry is favourable.

HESS J1825-137 was detected with a significance of 8.1 $\sigma$ in the Galactic Plane survey conducted with the H.E.S.S. instrument in 2004. Both HESS J1825-137 and the X-ray pulsar wind nebula G18.0--0.7 (associated with the Vela-like pulsar PSR B1823-13) are offset south of the pulsar, which may be the result of the SNR expanding into an inhomogeneous medium. The TeV size ($\sim 35$ pc, for a distance of 4 kpc) is $\sim 6$ times larger than the X-ray size, which may be the result of propagation effects as a result of the longer lifetime of TeV emitting electrons, compared to the relatively short lifetime of keV synchrotron emitting electrons. The TeV photon spectral index of $\sim 2.4$ can also be related to the extended PWN X-ray synchrotron photon index of $\sim 2.3$, if this spectrum is dominated by synchrotron cooling. The anomalously large size of the pulsar wind nebula can be explained if the pulsar was born with a relatively large initial spindown power and braking index $n\sim 2$, provided that the SNR expanded into the hot ISM with relatively low density ($\sim 0.003$ cm$^{-3}$).

PSR J0437-4715 is a millisecond pulsar (MSP) thought to be ``pair formation starved'' (having limited pair cascades due to magnetic photon absorption). Fortunately the general relativistic (GR) electrodynamical model under consideration applicable to this pulsar have few free parameters. We model PSR J0437-4715's visibility, using a 3D model which incorporates the variation of the GR E-field over the polar cap (PC), taking different observer and inclination angles into account. Using this pulsar as a case study, one may generalize to conducting a pulsar population visibility study. We lastly comment on the role of the proposed South African SKA (Square Kilometre Array) prototype, KAT (Karoo Array Telescope), for GLAST gamma-ray pulsar identification.

We simulate the magnetosphere of the nearby millisecond pulsar PSR J0437-4715, which is expected to have an unscreened electric potential due to the lack of magnetic pair production. We incorporate General Relativistic (GR) effects and study curvature radiation (CR) by primary electrons, but neglect inverse Compton (IC) scattering of thermal X-ray photons by these electrons. We find that the CR spectrum cuts off at energies below ~ 17 GeV, well below the threshold of the H.E.S.S. telescope <~ 100 GeV, while other models predict a much higher cutoff of >~ 100 GeV. GR theory also predicts a relatively narrow pulse of ~ 0.2 phase width centered on the magnetic axis. EGRET observations above 100 MeV significantly constrain the application of the Muslimov & Harding (1997) model for gamma-ray production as a result of GR frame dragging, and ultimately its polar cap (PC) current and accelerating potential. Whereas the standard prediction of this pulsar's gamma-ray luminosity due to GR frame dragging is ~10% of the spindown power, a non-detection by forthcoming H.E.S.S. observations will constrain it to <~ 0.3%, enforcing an even more severe revision of the accelerating electric field and PC current.

The poor angular resolution of imaging gamma-ray telescopes is offset by the large reflector areas of next generation telescopes such as MAGIC (17~m diameter), which makes the study of optical emission associated with some gamma-ray sources feasible. Furthermore, the extremely fast time response of photomultipliers (PMs) makes them ideal detectors for fast (subsecond) optical transients and periodic sources like pulsars. The optical pulse of the Crab pulsar was detected with the HEGRA CT1 central pixel using a modified PM, similar to the future MAGIC camera PMs. The purpose of these periodic observations was to determine the light of the night sky (LONS) for the galactic anticenter Crab region.Our results are between 2.5 and 3 times larger than the previously measured LONS (outside the galactic plane), as expected since the Crab pulsar is in the galactic plane, which implies a slightly higher energy threshold for Crab observations, if the higher value of CT1 measured LONS rate for galactic sources is used.

We present a new model for the synchrotron compact nebular emissions for Vela and PSR B1706-44, and derive fundamental pulsar/plerion parameters such as the pair production multiplicity, M, and wind magnetization parameter, sigma. The pair cascade above the pulsar polar cap, combined with the energy from the pulsar wind, injects a pair plasma into the surrounding environment. This wind, consisting of particles and fields, is shocked by the environment, resulting in synchrotron emission from "thermalized" pairs in a compact nebula. The broad band nebular spectrum depends crucially on the spindown power, distance, pair multiplicity, pulsar wind shock radius and sigma at this shock. We construct such a model for the particle spectra in the compact X-ray nebula of Vela and PSR B1706-44. Fits to the multiwavelength spectra of these sources indicate that 300<M<1000, whereas 0.05<sigma<0.5 for Vela. The same sigma interval (as for Vela) was independently derived from the radial gradients of the X-ray compact nebular emission from both Vela and PSR B1706-44, giving us confidence in our results. The M we derive for Vela is too small to be explained by curvature losses in the open magnetoshere of Vela, but the presence of optical pulsed photons could modify the predicted multiplicity.

We discuss the results of modelling the TeV gamma-ray spectra of two AGNs, Mkn 501 and Mkn 421 that have almost the same redshifts: z=0.031 and z=0.034, respectively. The effect of intergalactic gamma-ray absorption is treated as an uncertainty in the measurement of the intrinsic spectrum. Although the objects differ, we obtain satisfactory fits for both of them in a synchrotron self-Compton scenario. Compared to previous models, our fits are characterised by higher values of the Doppler factor (>= 50) and an electron injection spectrum extending to higher energies (Gmax = 1.5x10^5). In the case of Mkn 421, the observed difference in spectral slope in X-rays and TeV gamma-rays between the high and low states can be explained as a variation of a single parameter - the maximum energy Gmax mc^2 at which electrons are injected.

Using the recent models of Malkan & Stecker (2001) for the infrared background radiation and extrapolating them into the optical and UV range using recent galaxy count data, we rederive the optical depth of the Universe to high energy gamma-rays as a function of energy and redshift for energies between 50 GeV and 100 TeV and redshifts between 0.03 and 0.3. We then use these results to derive the intrinsic gamma-ray spectrum of Mkn 501 during its 1997 high state. We find that the time averaged spectral energy distribution of Mkn 501 while flaring had a broad, flat peak in the 5 to 10 TeV range which corresponds to a broad, flat time averaged X-ray peak in the 50 to 100 keV range observed during the flare. The spectral index of our derived intrinsic differential photon spectrum for Mkn 501 at energies below about 2 TeV was found to be apx. 1.6 to 1.7. This corresponds to a time averaged spectral index of 1.76 found in soft X-rays at energies below the X-ray (synchrotron) peak. These results appear to favor a synchrotron-self Compton origin for the TeV emission together wuth jet parameters which are consistent with time variability constraints within the context of a simple SSC model.

The non-detection of pulsed sub-TeV gamma-rays from EGRET pulsars proves that the EGRET pulsed spectra of all gamma-ray pulsars should terminate at energies below a few hundred GeV. The spectrum of a typical integrated pulse profile predicted by the polar cap model resemble typically a hard component, followed by a super exponential cutoff between 1 MeV (PSR B1509-58) and tens of GeV (e.g. Crab, PSR B1951+32 etc). Using a topological trigger for the H.E.S.S. Telescope System in the non-imaging mode, we show that H.E.S.S. should be able to detect pulsed emission from PSR B1706-44 within a few hours if the cutoff energy is above 30 GeV as suggested by EGRET observations. The recently detected radio pulsar PSR J1837-0604 (pulsar period: 96 ms) associated with the unidentified EGRET source GeV J1837-06010 should also be detectable within a few hours if the source is pulsed and if its cutoff is similar to that of PSR B1706-44. H.E.S.S. should even be able to image middle-aged, low-multiplicity pulsars for which the mean photon energy is expected to be well above 10 GeV. Such observations should provide important constraints on the final evolutionary status of gamma-ray pulsars and millisecond pulsars in general.

Seventeen southern sky BL Lacs were observed in UBVRI using the CCD Camera on the 1.0m telescope at the South African Astronomical Observatory (SAAO) in Aug and Nov 1999. The analyses of all the seventeen sources are now complete, and are available via anonymous ftp (ftp pukrs1.puk.ac.za/pub/Blazars). A few examples of our results are however given in this paper. Whereas PKS 2005-489 and PKS 2155-304 apear to have been in a high state, PKS 0048-097 and PKS 0521-365 showed evidence of variability on a time-scale of a few days, with the amplitude of variability increasing towards short wavelengths. This is consistent with observations of gamma-ray BL Lacs, which show similar behaviour in optical and X-rays.

Observations and theoretical studies have demonstrated that the pulsed spectra of all gamma-ray pulsars terminate at energies below a few hundred GeV. In most cases we expect this cutoff energy E_o to be around 10 GeV. Only with next generation ground-based gamma-ray telescopes, which are expected to have non-zero trigger probabilities near 10 GeV, can we expect to detect pulsations. The large gamma-ray fluxes below E_o, together with the associated sharp pulse profiles, compensate for the lack of imaging capability near threshold. For H.E.S.S. we find that the pulsed component of PSR B1706-44 should be detectable near threshold, whereas the unidentified GeV EGRET sources should be detectable if the superexponential cutoff energy is larger that approximately 30 GeV for relatively hard pulsar photon spectra (~E^-1.5).

The detection of unpulsed TeV emission from PSR B1706-44 raised the question if this VHE emission results from inverse Compton scattering of multi-TeV electrons on the CMBR and other soft photon fields. This hints at the existence of an unseen synchrotron nebula, which is larger than the observed 1 arcminute compact X-ray nebula. In a search for synchrotron counterpart of the TeV nebula, we have taken V-band CCD images of the region around the PSR B1706-44 pulsar position. By rejecting all point source contributions down to the diffuse sky limit, we were able to search for such extended plerionic synchrotron emission on scale sizes limited by the angular resolution of the CANGAROO imaging TeV observations. By combining the V-band flux limits with the observed TeV flux, we were able to obtain upper limits for the magnetic field strength as a function of the radius of the assumed TeV nebula. Assuming that the 2 arcminute radio plerion also defines the size of the TeV nebula, we constrain the steady state particle spectral index to the values between 4 and 5.5, with a field strength limit of about 3x10^-6 G. This compares with an index of 3.6+/-0.6 inside the 1 arcminute compact X-ray nebula.

We propose a new hypothesis for the origin of non-solar hadronic cosmic rays (CRs) at all energies: Highly relativistic, narrowly collimated jets from the birth or collapse of neutron stars (NSs) in our Galaxy accelerate ambient disk and halo matter to CR energies and disperse it in ``hot spots'' which they form when they stop in the Galactic halo. Such events - ``Galactic Gamma-Ray Bursters'' (GGRBs) - are proposed to cause cosmological gamma-ray bursts (GRBs) in other galaxies when their beamed radiation happens to point in our direction. Our hypothesis naturally explains some observations which are difficult to understand with the currently popular ideas about CR origin - e.g. the small Galacto-centric gradient of the cosmic-ray density and the absence of the Greisen-Zatsepin-Kuzmin cutoff. Our idea stands or falls with the existence of the ``hot spots'' (``GGRB remnants'') in the Galactic halo. We discuss their expected observational signatures and find that they could appear as EGRET unidentified high-latitude sources.

The BL Lac Object Mkn 501 was in a state of high activity in the TeV range in 1997. During this time Mkn 501 was observed by all Cherenkov-Telescopes of the HEGRA-Collaboration. Part of the data were also taken during moonshine thus providing a nearly continuous coverage for this object in the TeV-range. We have carried out a QPO analysis and found evidence for a 23 day periodicity. We applied the same analysis on the 'data by dwell' x-ray lightcurve from the RXTE/ASM database and found also evidence for the 23 day periodicity. The combined probability was P = 2.8 *10^-4.

We present a new calculation of the intergalactic gamma-ray absorption coefficient as a function of both energy and redshift. In reexamining this problem, we make use of a new, empirically based calculation (as opposed to previous model calculations) of the intergalactic infrared radiation field. We find smaller opacities than those given previously (Stecker & De Jager 1997). We apply our results to the new observations of the flaring gamma-ray spectra of Mrk421 and Mrk501, both at a redshift of apx. 0.03. Our new calculations indicate that there should be no significant curvature in the spectra of these sources for energies below 10 TeV, as indicated by recent observations. However, the intrinsic spectra of these sources should be harder by apx. 0.2 to 0.45 in the spectral index in the 1 to 10 TeV range with an intergalactic absorption cutoff above apx. 20 TeV.

Due to observations made by the Compton Gamma-Ray Observatory over the last six years, it appears that a number of galactic supernova remnants may be candidates for sources of cosmic gamma-rays. These include shell-type remnants such as IC443 and $\gamma$ Cygni, which have no known parent pulsars, but have significant associations with unidentified EGRET sources, and others that appear to be composite, where a pulsar is embedded in a shell (e.g. W44 and Vela), or are purely pulsar-driven, such as the Crab Nebula. This review discusses our present understanding of gamma-ray production in plerionic and non-plerionic supernova remnants, and explores the relationship between such emission and that in other wavebands. Focuses include models of the Crab and Vela nebulae, the composite nature of W44, the relationship of shell-type remnants to cosmic ray production, the relative importance of shock-accelerated protons and electrons, constraints on models placed by TeV, X-ray and radio observations, and the role of electrons injected directly into the remnants by parent pulsars.

We discuss the problem of the absorption of very high-energy gamma-rays by pair production interactions with extragalactic photons which originate from stellar emission in the near IR-UV and reradiation of starlight in the mid- and far-IR. The absorption of gamma-rays above 1 TeV is dominated by interactions with infrared photons. We make a new determination of the optical depth of the universe to multi-TeV photons as a function of energy and redshift and use the results to compare with recent spectral data of Mrk 421 and Mrk 501, sources that have been observed in the flaring state up to apx. 10 TeV energy. For the optical depth calculations, we have made use of a new, empirically based calculation of the intergalactic radiation field by Malkan & Stecker which we consider to be more accurate than that based on previous theoretical modeling. We also discuss the absorption of sub-TeV gamma-rays by starlight photons at high redshifts.

We suggest that low-redshift XBLs (X-ray selected BL Lacertae objects) may be the only extragalactic gamma-ray sources observable at TeV energies. We use simple physical considerations involving synchrotron and Compton component spectra for blazars to suggest why the observed TeV sources are XBLs, whereas mostly RBLs and FSRQs are seen at GeV energies. These considerations indicate that the differences between XBLs and RBLs cannot be explained purely as relativistic jet orientation effects. We note that the only extragalactic TeV sources which have been observed are XBLs and that a nearby RBL with a very hard spectrum in the GeV range has not been seen at TeV energies. We also note that of the 14 BL Lacs observed by EGRET, 12 are RBLs, whereas only 2 are XBLs. We give a list of nearby XBLs which we consider to be good candidate TeV sources and predict estimated TeV fluxes for these objects.

We present a new calculation of the intergalactic $\gamma$-ray pair-production absorption coefficient as a function of both energy and redshift up to the redshift of 3C279, z = 0.54. In reexamining this problem, we make use of new observational data on the intergalactic infrared radiation field (IIRF), together with recent theoretical models of the galactic spectral energy distributions of the IIRF from stars and dust reradiation and estimates of the IIRF from galaxy counts and {\it COBE} results. We present our results for two fairly well defined IIRF spectral energy distributions, one of which is within $1 \sigma $ of our previous estimate of the IIRF at $ \sim 20$ $\mu$m. We then apply our results to the $\gamma$-ray spectrum of Mrk 421, and obtain good agreement with the observational data, including the recent results of the {\it HEGRA} group.

Supernova 1006 is the first shell type supernova remnant to show evidence of particle acceleration to TeV energies. In the present paper we examine this possibility by modeling the observed X-ray non-thermal emission in terms of synchrotron radiation from Fermi accelerated electrons. The predicted synchrotron spectrum fits the radio and non-thermal component of the observed soft X-ray to hard X-ray emission quite well. These particles can produce TeV gamma rays by inverse Compton scattering on the microwave radiation and other ambient fields, and the derived electron distribution is also used to calculate the expected inverse Compton flux. We find that if the remnant is characterised by a magnetic field strength lower than $\sim 7\mu$G, then the TeV flux can be higher than that of the Crab Nebula. About 75\% of the TeV emission from SN 1006 is expected to be concentrated in the synchrotron bright NE and SW rims (the ``hard aegis") of the remnant, which would allow a sensitive search if the Atmospheric Imaging Cherenkov Technique is used.

Following our previously proposed technique, we have used the recent gamma-ray observations of Mrk 421 to place theoretically significant constraints on and possible estimates of the intergalactic infrared radiation field (IIRF) which are consistent with normal galactic IR production by stars and dust and rule out exotic mechanisms proposed to produce a larger IIRF. Using models for the low energy intergalactic photon spectrum from microwave to UV energies, we calculate the opacity of inter- galactic space to gamma-rays as a function of energy and redshift. These calculations indicate that the GeV gamma-ray burst recently observed by the EGRET experiment on CGRO originates at a redshift less than approximately 1.5.