Multiwavelength observations of the high-frequency-peaked blazar 1ES2344+514 were performed from 2007 October to 2008 January. The campaign represents the first contemporaneous data on the object at very high energy (VHE, E >100 GeV) {\gamma}-ray, X-ray, and UV energies. Observations with VERITAS in VHE {\gamma}-rays yield a strong detection of 20 {\sigma} with 633 excess events in a total exposure of 18.1 hours live-time. A strong VHE {\gamma}-ray flare on 2007 December 7 is measured at F(>300 GeV) = (6.76 \pm 0.62) \times 10-11 ph cm-2 s-1, corresponding to 48% of the Crab Nebula flux. Excluding this flaring episode, nightly variability at lower fluxes is observed with a time-averaged mean of F(>300 GeV) = (1.06 \pm 0.09) \times 10-11 ph cm-2 s-1 (7.6% of the Crab Nebula flux). The differential photon spectrum between 390 GeV and 8.3 TeV for the time-averaged observations excluding 2007 December 7 is well described by a power law with a photon index of {\Gamma} = 2.78 \pm 0.09stat \pm 0.15syst. Over the full period of VERITAS observations contemporaneous X-ray and UV data were taken with Swift and RXTE. The measured 2-10 keV flux ranged by a factor of ~7 during the campaign. On 2007 December 8 the highest ever observed X-ray flux from 1ES 2344+514 was measured by Swift XRT at a flux of F(2-10 keV) = (6.28 \pm 0.31) \times 10-11 erg cm-2 s-1. Evidence for a correlation between the X-ray flux and VHE {\gamma}-ray flux on nightly time-scales is indicated with a Pearson correlation coefficient of r = 0.60 \pm 0.11. Contemporaneous spectral energy distributions (SEDs) of 1ES 2344+514 are presented for two distinct flux states. A one-zone synchrotron self-Compton (SSC) model describes both SEDs using parameters consistent with previous SSC modeling of 1ES 2344+514 from non-contemporaneous observations.

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

The intermediate-frequency peaked BL Lacertae (IBL) object 3C 66A is detected during 2007 - 2008 in VHE (very high energy: E > 100 GeV) gamma-rays with the VERITAS stereoscopic array of imaging atmospheric Cherenkov telescopes. An excess of 1791 events is detected, corresponding to a significance of 21.2 standard deviations (sigma), in these observations (32.8 hours live time). The observed integral flux above 200 GeV is 6% of the Crab Nebula's flux and shows evidence for variability on the time-scale of days. The measured energy spectrum is characterized by a soft power law with photon index Gamma = 4.1 +- 0.4_stat +- 0.6_sys. The radio galaxy 3C 66B is excluded as a possible source of the VHE emission.

Indirect dark matter searches with ground-based gamma-ray observatories provide an alternative for identifying the particle nature of dark matter that is complementary to that of direct search or accelerator production experiments. We present the results of observations of the dwarf spheroidal galaxies Draco, Ursa Minor, Bootes 1, and Willman 1 conducted by VERITAS. These galaxies are nearby dark matter dominated objects located at a typical distance of several tens of kiloparsecs for which there are good measurements of the dark matter density profile from stellar velocity measurements. Since the conventional astrophysical background of very high energy gamma rays from these objects appears to be negligible, they are good targets to search for the secondary gamma-ray photons produced by interacting or decaying dark matter particles. No significant gamma-ray flux above 200 GeV was detected from these four dwarf galaxies for a typical exposure of ~20 hours. The 95% confidence upper limits on the integral gamma-ray flux are in the range 0.4-2.2x10^-12 photons cm^-2s^-1. We interpret this limiting flux in the context of pair annihilation of weakly interacting massive particles and derive constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the WIMPs. The limits are obtained under conservative assumptions regarding the dark matter distribution in dwarf galaxies and are approximately three orders of magnitude above the generic theoretical prediction for WIMPs in the minimal supersymmetric standard model framework. However significant uncertainty exists in the dark matter distribution as well as the neutralino cross sections which under favorable assumptions could further lower the limits.

We report the discovery of very high energy gamma-ray emission from the direction of the SNR G54.1+0.3 using the VERITAS ground-based gamma-ray observatory. The TeV signal has an overall significance of 6.8$\sigma$ and appears point-like given the 5$^{arcminute}$ resolution of the instrument. The integral flux above 1 TeV is 2.5% of the Crab Nebula flux and significant emission is measured between 250 GeV and 4 TeV, well described by a power-law energy spectrum dN/dE $\sim$ E$^{-\Gamma}$ with a photon index $\Gamma= 2.39\pm0.23_{stat}\pm0.30_{sys}$. We find no evidence of time variability among observations spanning almost two years. Based on the location, the morphology, the measured spectrum, the lack of variability and a comparison with similar systems previously detected in the TeV band, the most likely counterpart of this new VHE gamma-ray source is the PWN in the SNR G54.1+0.3. The measured X-ray to VHE gamma-ray luminosity ratio is the lowest among all the nebulae supposedly driven by young rotation-powered pulsars, which could indicate a particle-dominated PWN.

M87 is a nearby radio galaxy that is detected at energies ranging from radio to VHE gamma-rays. Its proximity and its jet, misaligned from our line-of-sight, enable detailed morphological studies and extensive modeling at radio, optical, and X-ray energies. Flaring activity was observed at all energies, and multi-wavelength correlations would help clarify the origin of the VHE emission. In this paper, we describe a detailed temporal and spectral analysis of the VERITAS VHE gamma-ray observations of M87 in 2008 and 2009. In the 2008 observing season, VERITAS detected an excess with a statistical significance of 7.2 sigma from M87 during a joint multi-wavelength monitoring campaign conducted by three major VHE experiments along with the Chandra X-ray Observatory. In February 2008, VERITAS observed a VHE flare from M87 occurring over a 4-day timespan. The peak nightly flux above 250GeV was 7.7% of the Crab Nebula flux. M87 was marginally detected before this 4-day flare period, and was not detected afterwards. Spectral analysis of the VERITAS observations showed no significant change in the photon index between the flare and pre-flare states. Shortly after the VHE flare seen by VERITAS, the Chandra X-ray Observatory detected the flux from the core of M87 at a historical maximum, while the flux from the nearby knot HST-1 remained quiescent. Acciari et al. (2009) presented the 2008 contemporaneous VHE gamma-ray, Chandra X-ray, and VLBA radio observations which suggest the core as the most likely source of VHE emission, in contrast to the 2005 VHE flare that was simultaneous with an X-ray flare in the HST-1 knot. In 2009, VERITAS continued its monitoring of M87 and marginally detected a 4.2 sigma excess corresponding to a flux of ~1% of the Crab Nebula. No VHE flaring activity was observed in 2009.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

We report on observations of very high-energy gamma rays from the shell-type supernova remnant Cassiopeia A with the VERITAS stereoscopic array of four imaging atmospheric Cherenkov telescopes in Arizona. The total exposure time for these observations is 22 hours, accumulated between September and November of 2007. The gamma-ray source associated with the SNR Cassiopeia A was detected above 200 GeV with a statistical significance of 8.3 s.d. The estimated integral flux for this gamma-ray source is about 3% of the Crab-Nebula flux. The photon spectrum is compatible with a power law dN/dE ~ E^(-Gamma) with an index Gamma = 2.61 +/- 0.24(stat) +/- 0.2(sys). The data are consistent with a point-like source. We provide a detailed description of the analysis results, and discuss physical mechanisms that may be responsible for the observed gamma-ray emission.

VERITAS is a state-of-the-art ground-based gamma-ray observatory that operates in the very high-energy (VHE) region of 100 GeV to 50 TeV. The observatory consists of an array of four 12m-diameter imaging atmospheric Cherenkov telescopes located in southern Arizona, USA. The four-telescope array has been fully operational since September 2007, and over the last two years, VERITAS has been operating with high efficiency and with excellent performance. This talk summarizes the recent results from VERITAS, including the discovery of eight new VHE gamma-ray sources.

VERITAS is the high-sensitivity instrument of latest generation. It is often used for the short AGN monitoring exposures evenly distributed over entire observational season of a source of interest. Each of these exposures is long enough to detect the source at the flux level of about 1 Crab. During the 2009 observing season a number of exposures of Mkn 501 with VERITAS revealed variable TeV gamma-ray emission at the flux level eventually exceeding 2 Crab. The spectral and flux variability measurements in TeV gamma rays for the 2009 data sample of Mkn 501 are summarized in this paper.

We report the detection of very-high-energy (VHE) gamma-ray emission from supernova remnant (SNR) G106.3+2.7. Observations performed in 2008 with the VERITAS atmospheric Cherenkov gamma-ray telescope resolve extended emission overlapping the elongated radio SNR. The 7.3 sigma (pre-trials) detection has a full angular extent of roughly 0.6deg by 0.4deg. Most notably, the centroid of the VHE emission is centered near the peak of the coincident 12CO (J = 1-0) emission, 0.4deg away from the pulsar PSR J2229+6114, situated at the northern end of the SNR. Evidently the current-epoch particles from the pulsar wind nebula are not participating in the gamma-ray production. The VHE energy spectrum measured with VERITAS is well characterized by a power law dN/dE = N_0(E/3 TeV)^{-G} with a differential index of G = 2.29 +/- 0.33stat +/- 0.30sys and a flux of N_0 = (1.15 +/- 0.27stat +/- 0.35sys)x 10^{-13} cm^{-2} s^{-1} TeV^{-1}. The integral flux above 1 TeV corresponds to ~5 percent of the steady Crab Nebula emission above the same energy. We describe the observations and analysis of the object and briefly discuss the implications of the detection in a multiwavelength context.

Although Galactic cosmic rays (protons and nuclei) are widely believed to be dominantly accelerated by the winds and supernovae of massive stars, definitive evidence of this origin remains elusive nearly a century after their discovery [1]. The active regions of starburst galaxies have exceptionally high rates of star formation, and their large size, more than 50 times the diameter of similar Galactic regions, uniquely enables reliable calorimetric measurements of their potentially high cosmic-ray density [2]. The cosmic rays produced in the formation, life, and death of their massive stars are expected to eventually produce diffuse gamma-ray emission via their interactions with interstellar gas and radiation. M 82, the prototype small starburst galaxy, is predicted to be the brightest starburst galaxy in gamma rays [3, 4]. Here we report the detection of >700 GeV gamma rays from M 82. From these data we determine a cosmic-ray density of 250 eV cm-3 in the starburst core of M 82, or about 500 times the average Galactic density. This result strongly supports that cosmic-ray acceleration is tied to star formation activity, and that supernovae and massive-star winds are the dominant accelerators.

The recent detection by the Fermi gamma-ray space telescope of high-energy gamma-rays from the radio galaxy NGC 1275 makes the observation of the very high energy (VHE: E > 100 GeV) part of its broadband spectrum particularly interesting, especially for the understanding of active galactic nuclei (AGN) with misaligned multi-structured jets. The radio galaxy NGC 1275 was recently observed by VERITAS at energies above 100 GeV for about 8 hours. No VHE gamma-ray emission was detected by VERITAS from NGC 1275. A 99% confidence level upper limit of 2.1% of the Crab Nebula flux level is obtained at the decorrelation energy of approximately 340 GeV, corresponding to 19% of the power-law extrapolation of the Fermi Large Area Telescope (LAT) result.

We report results from an intensive multiwavelength campaign on the intermediate-frequency-peaked BL Lacertae object W Com (z=0.102) during a strong outburst of very high energy gamma-ray emission in June 2008. The very high energy gamma-ray signal was detected by VERITAS on 2008 June 7-8 with a flux F(>200 GeV) = (5.7+-0.6)x10^-11 cm-2s-1, about three times brighter than during the discovery of gamma-ray emission from W Com by VERITAS in 2008 March. The initial detection of this flare by VERITAS at energies above 200 GeV was followed by observations in high energy gamma-rays (AGILE, E>100 MeV), and X-rays (Swift and XMM-Newton), and at UV, and ground-based optical and radio monitoring through the GASP-WEBT consortium and other observatories. Here we describe the multiwavelength data and derive the spectral energy distribution (SED) of the source from contemporaneous data taken throughout the flare.

We present evidence that the very-high-energy (VHE, E > 100 GeV) gamma-ray emission coincident with the supernova remnant IC 443 is extended. IC 443 contains one of the best-studied sites of supernova remnant/molecular cloud interaction and the pulsar wind nebula CXOU J061705.3+222127, both of which are important targets for VHE observations. VERITAS observed IC 443 for 37.9 hours during 2007 and detected emission above 300 GeV with an excess of 247 events, resulting in a significance of 8.3 standard deviations (sigma) before trials and 7.5 sigma after trials in a point-source search. The emission is centered at 06 16 51 +22 30 11 (J2000) +- 0.03_stat +- 0.08_sys degrees, with an intrinsic extension of 0.16 +- 0.03_stat +- 0.04_sys degrees. The VHE spectrum is well fit by a power law (dN/dE = N_0 * (E/TeV)^-Gamma) with a photon index of 2.99 +- 0.38_stat +- 0.3_sys and an integral flux above 300 GeV of (4.63 +- 0.90_stat +- 0.93_sys) * 10^-12 cm^-2 s^-1. These results are discussed in the context of existing models for gamma-ray production in IC 443.

HESS J0632+057 is one of only two unidentified very-high-energy gamma-ray sources which appear to be point-like within experimental resolution. It is possibly associated with the massive Be star MWC 148 and has been suggested to resemble known TeV binary systems like LS I +61 303 or LS 5039. HESS J0632+057 was observed by VERITAS for 31 hours in 2006, 2008 and 2009. During these observations, no significant signal in gamma rays with energies above 1 TeV was detected from the direction of HESS J0632+057. A flux upper limit corresponding to 1.1% of the flux of the Crab Nebula has been derived from the VERITAS data. The non-detection by VERITAS excludes with a probability of 99.993% that HESS J0632+057 is a steady gamma-ray emitter. Contemporaneous X-ray observations with Swift XRT reveal a factor of 1.8+-0.4 higher flux in the 1-10 keV range than earlier X-ray observations of HESS J0632+057. The variability in the gamma-ray and X-ray fluxes supports interpretation of the ob ject as a gamma-ray emitting binary.

We present results from a long-term monitoring campaign on the TeV binary LSI +61 303 with VERITAS at energies above 500 GeV, and in the 2-10 keV hard X-ray bands with RXTE and Swift, sampling nine 26.5 day orbital cycles between September 2006 and February 2008. The binary was observed by VERITAS to be variable, with all integrated observations resulting in a detection at the 8.8 sigma (2006/2007) and 7.3 sigma (2007/2008) significance level for emission above 500 GeV. The source was detected during active periods with flux values ranging from 5 to 20% of the Crab Nebula, varying over the course of a single orbital cycle. Additionally, the observations conducted in the 2007-2008 observing season show marginal evidence (at the 3.6 sigma significance level) for TeV emission outside of the apastron passage of the compact object around the Be star. Contemporaneous hard X-ray observations with RXTE and Swift show large variability with flux values typically varying between 0.5 and 3.0*10^-11 ergs cm^-2 s^-1 over a single orbital cycle. The contemporaneous X-ray and TeV data are examined and it is shown that the TeV sampling is not dense enough to detect a correlation between the two bands.

The VERITAS collaboration reports the detection of very-high-energy (VHE) gamma-ray emission from the high-frequency-peaked BL Lac object 1ES 1218+304 located at a redshift of z=0.182. A gamma-ray signal was detected with a statistical significance of 10.4 standard deviations (10.4 sigma) for the observations taken during the first three months of 2007, confirming the discovery of this object made by the MAGIC collaboration. The photon spectrum between ~160 GeV and ~1.8 TeV is well described by a power law with an index of Gamma = 3.08 +/- 0.34_stat +/- 0.2_sys. The integral flux is Phi(E > 200 GeV) = (12.2 +/- 2.6) X 10^-12 cm^-2 s^-1, which corresponds to ~6% of that of the Crab Nebula. The light curve does not show any evidence for VHE flux variability. Using lower limits on the density of the extragalactic background light in the near to mid-infrared we are able to limit the range of intrinsic energy spectra for 1ES 1218+304. We show that the intrinsic photon spectrum has an index that is harder than Gamma = 2.32 +/- 0.37_stat. When including constraints from the spectra of 1ES 1101-232 and 1ES 0229+200, the spectrum of 1ES 1218+304 is likely to be harder than Gamma = 1.86 +/- 0.37_stat.

Since September 2005, the Whipple 10m Gamma-ray Telescope has been operated primarily as a blazar monitor. The five Northern Hemisphere blazars that have already been detected at the Whipple Observatory, Markarian 421, H1426+428, Markarian 501, 1ES 1959+650 and 1ES 2344+514, are monitored routinely each night that they are visible. We report on the Markarian 421 observations taken from November 2005 to June 2006 in the gamma-ray, X-ray, optical and radio bands. During this time, Markarian 421 was found to be variable at all wavelengths probed. Both the variability and the correlations among different energy regimes are studied in detail here. A tentative correlation, with large spread, was measured between the X-ray and gamma-ray bands, while no clear correlation was evident among the other energy bands. In addition to this, the well-sampled spectral energy distribution of Markarian 421 (1101+384) is presented for three different activity levels. The observations of the other blazar targets will be reported separately.

The high-frequency-peaked BL-Lacertae object \objectname{1ES 0806+524}, at redshift z=0.138, was observed in the very-high-energy (VHE) gamma-ray regime by VERITAS between November 2006 and April 2008. These data encompass the two-, and three-telescope commissioning phases, as well as observations with the full four-telescope array. \objectname{1ES 0806+524} is detected with a statistical significance of 6.3 standard deviations from 245 excess events. Little or no measurable variability on monthly time scales is found. The photon spectrum for the period November 2007 to April 2008 can be characterized by a power law with photon index $3.6 \pm 1.0_{\mathrm{stat}} \pm 0.3_{\mathrm{sys}}$ between $\sim$300 GeV and $\sim$700 GeV. The integral flux above 300 GeV is $(2.2\pm0.5_{\mathrm{stat}}\pm0.4_{\mathrm{sys}})\times10^{-12}\:\mathrm{cm}^{2}\:\mathrm{s}^{-1}$ which corresponds to 1.8% of the Crab Nebula flux. Non contemporaneous multiwavelength observations are combined with the VHE data to produce a broadband spectral energy distribution that can be reasonably described using a synchrotron-self Compton model.

This is a report on the findings of the technology working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses different technology opportunities and the technical feasibility for substantially improving IACTS and ground based particle detectors to achieve an order of magnitude better sensitivity than the instruments employed today as well as their planned upgrades. A technology roadmap for improving IACTS and ground based particle detectors is presented.

This is a report on the findings of the Galactic compact objects working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper is an APS commissioned document, and the full version has also been released and can be found on astro-ph. This section of the white paper discusses the potential of future ground-based gamma-ray experiments to advance our understanding of the physics of Galactic compact objects including pulsars, pulsar wind nebulae, and X-ray binaries.

VERITAS, an Imaging Atmospheric Cherenkov Telescope (IACT) system for gammma-ray astronomy in the GeV-TeV range, has recently completed its first season of observations with a full array of four telescopes. A number of astrophysical gamma-ray sources have been detected, both galactic and extragalactic, including sources previously unknown at TeV energies. We describe the status of the array and some highlight results, and assess the technical performance, sensitivity and shower reconstruction capabilities.

We report the detection of very high-energy gamma-ray emission from the intermediate-frequency-peaked BL Lacertae object W Comae (z=0.102) by VERITAS. The source was observed between January and April 2008. A strong outburst of gamma-ray emission was measured in the middle of March, lasting for only four days. The energy spectrum measured during the two highest flare nights is fit by a power-law and is found to be very steep, with a differential photon spectral index of Gamma = 3.81 +- 0.35_stat +- 0.34_syst. The integral photon flux above 200GeV during those two nights corresponds to roughly 9% of the flux from the Crab Nebula. Quasi-simultaneous Swift observations at X-ray energies were triggered by the VERITAS observations. The spectral energy distribution of the flare data can be described by synchrotron-self-Compton (SSC) or external-Compton (EC) leptonic jet models, with the latter offering a more natural set of parameters to fit the data.

Many of the recently discovered TeV gamma-ray sources are associated with pulsar wind nebulae (PWNe). In fact, they represent the most populous class of Galactic sources at TeV energies. In addition, HESS has also discovered, in a survey of the Galactic plane, a population of TeV gamma-ray sources that are still without definitive counterparts at longer wavelengths. For a number of these sources, a pulsar is an evident association, which is often located within an extended region of the TeV gamma-ray emission. These particular HESS sources are promising candidates for yet not resolved pulsar wind nebulae. Here we have undertaken a systematic search for X-ray counterparts of the sources, using the archival Chandra data, within the spatial bounds of the unidentified HESS sources. A number of X-ray sources have been detected in the Chandra fields. Two of them, CXOU J161729.3-505512 and CXOU J170252.4-412848, are of a special interest because of their excellent positional coincidence with the pulsars PSR J1617-5055 and PSR J1702-4128, respectively. The first source is extended, with a bright core of 2.6 arcsec (FWHM) in radius but the emission can be seen up to roughly 20 arcsec. The second one is much fainter and detected only with marginal significance (4.6 sigma). It might also be slightly extended, although the situation is quite uncertain due to very limited statistics. The analysis of the archival Chandra data for a middle-aged pulsar (PSR J1913+1011) does not reveal any statistically significant excess at and around the position of the pulsar or the center of gravity of its plausible TeV gamma-ray counterpart (HESS J1912+101). We discuss the implications of the results.

LS I +61 303 is one of only a few high-mass X-ray binaries currently detected at high significance in very high energy gamma-rays. The system was observed over several orbital cycles (between September 2006 and February 2007) with the VERITAS array of imaging air-Cherenkov telescopes. A signal of gamma-rays with energies above 300 GeV is found with a statistical significance of 8.4 standard deviations. The detected flux is measured to be strongly variable; the maximum flux is found during most orbital cycles at apastron. The energy spectrum for the period of maximum emission can be characterized by a power law with a photon index of Gamma=2.40+-0.16_stat+-0.2_sys and a flux above 300 GeV corresponding to 15-20% of the flux from the Crab Nebula.

The multiwavelength observation of the nearby radio galaxy M87 provides a unique opportunity to study in detail processes occurring in Active Galactic Nuclei from radio waves to TeV gamma-rays. Here we report the detection of gamma-ray emission above 250 GeV from M87 in spring 2007 with the VERITAS atmospheric Cherenkov telescope array and discuss its correlation with the X-ray emission. The gamma-ray emission is measured to be point-like with an intrinsic source radius less than 4.5 arcmin. The differential energy spectrum is fitted well by a power-law function: dPhi/dE=(7.4+-1.3_{stat}+-1.5_{sys})(E/TeV)^{-2.31+-0.17_{stat}+-0.2_{sys}} 10^{-9}m^{-2}s^{-1}TeV^{-1}. We show strong evidence for a year-scale correlation between the gamma-ray flux reported by TeV experiments and the X-ray emission measured by the ASM/RXTE observatory, and discuss the possible short-time-scale variability. These results imply that the gamma-ray emission from M87 is more likely associated with the core of the galaxy than with other bright X-ray features in the jet.

The HEGRA gamma-ray source TeV J2032+4130 is considered the prototypical 'dark accelerator', since it was the first TeV source detected with no firm counterparts at lower frequencies. The Whipple collaboration observed this source in 2003-5 and the emission hotspot appears displaced about 9 arcminutes to the northeast of the HEGRA position, though given the large positional uncertainties the HEGRA and Whipple positions are consistent. Here we report on Westerbork Synthesis Radio Telescope (WSRT), Very Large Array (VLA), Chandra and INTEGRAL data covering the locations of the Whipple and HEGRA hotspots. We confirm a dual-lobed radio source (also see Marti et al., 2007) coincident with the Whipple hotspot, as well as a weak, partially non-thermal shell-like object, with a location and morphology very similar to the HEGRA source, in our WSRT and mosaicked VLA datasets, respectively. Due to its extended nature, it is likely that the latter structure is a more plausible counterpart of the reported very high energy (VHE) gamma-ray emissions in this region. If so, TeV J2032+4130 may not be a 'dark accelerator' after all. Further observations with the new generation of imaging Cherenkov telescopes are needed to pin down the precise location and morphology of the TeV emission region and thus clear up the confusion over its possible lower frequency counterparts.

We present observations of the dwarf galaxies Draco and Ursa Minor, the local group galaxies M32 and M33, and the globular cluster M15 conducted with the Whipple 10m gamma-ray telescope to search for the gamma-ray signature of self-annihilating weakly interacting massive particles (WIMPs) which may constitute astrophysical dark matter (DM). We review the motivations for selecting these sources based on their unique astrophysical environments and report the results of the data analysis which produced upper limits on excess rate of gamma rays for each source. We consider models for the DM distribution in each source based on the available observational constraints and discuss possible scenarios for the enhancement of the gamma-ray luminosity. Limits on the thermally averaged product of the total self-annihilation cross section and velocity of the WIMP, <\sigma v>, are derived using conservative estimates for the magnitude of the astrophysical contribution to the gamma-ray flux. Although these limits do not constrain predictions from the currently favored theoretical models of supersymmetry (SUSY), future observations with VERITAS will probe a larger region of the WIMP parameter phase space, <\sigma v> and WIMP particle mass (m_\chi).

Many of the recently discovered galactic very high-energy (VHE) gamma-ray sources are associated with Pulsar Wind Nebulae, which is the most populous Galactic source category at TeV energies. The extended synchrotron nebulae of these objects observed in the X-ray band are a hallmark of the relativistic winds, generated by the young, energetic neutron stars, that interact with the matter ejected by the supernova explosion and the surrounding interstellar gas. Relativistic electrons, or protons, accelerated in the pulsar winds, or at their shock boundaries, interact with the magnetic field and low energy seed photons to produce the observed VHE gamma-ray emission. The VERITAS array of four imaging atmospheric Cherenkov telescopes was designed to study astrophysical sources of gamma rays in the energy domain from about 100 GeV up to several tens of TeV. The sensitivity of the VERITAS array allows detailed studies of the morphology and spectral features of gamma-ray emission from PWNe. Three northern sky PWNe, G75.2+0.1, G106.6+2.9, and 3C58, were observed with VERITAS during 2006. No evidence for TeV gamma-ray emission at the position of the pu lsar associated with these PWNe is demonstrated.

The BL Lac object (blazar) Mrk 421 was observed during its outburst in April 2004 with the Whipple 10 m telescope for a total of about 24.5 hours. The measured gamma-ray rate varied substantially over the range from 4 to 10 gamma's/min and eventually exceeded the steady gamma-ray rate of the Crab Nebula (standard candle) by a factor of 3. The overall significance of the gamma-ray signal exceeded 70 sigma and the total number of excess events was more than 10,000. The signal light curve does not show any particular variability pattern. This unique Mrk 421 outburst enabled the measurement of a high quality spectrum of very high-energy gamma rays in a high state of emission. This spectrum is a power-law and it extends beyond 10 TeV.

Recently, ground-based very high-energy (VHE) gamma-ray astronomy achieved a remarkable advancement in the development of the observational technique for the registration and study of gamma-ray emission above 100 GeV. Construction of telescopes of substantially larger sizes than the currently used 12 m class telescopes can drastically improve the sensitivity of ground-based detectors to gamma rays of energy from 10 GeV to 100 GeV. Based on Monte Carlo simulations we have studied the response of an array of three large area imaging atmospheric Cherenkov telescopes (IACT) as a prototype for a future large-scale low energy ground-based experiment. The sensitivity of a three-telescope array as a function of optical reflector size was investigated here in detail.

We carried out a 5-pointing mosaic observation of TeV J2032+4130 at 1.4 and 4.8 GHz with the VLA in April of 2003. The analysis of the 4.8GHz data indicate weak wispy shell-like radio structure(s) which are at least partially non-thermal. The radio data is compatible with one or more young supernova remnants or perhaps the signature of large scale cluster shocks in this region induced by the violent action of the many massive stars in Cyg OB2.

During the last decade ground-based very high-energy gamma-ray astronomy achieved a remarkable advancement in the development of the observational technique for the registration and study of gamma-ray emission above 100 GeV. It is widely believed that the next step in its future development will be the construction of telescopes of substantially larger size than the currently used 10 m class telescopes. This can drastically improve the sensitivity of the ground-based detectors for gamma rays of energy from 10 to 100 GeV. Based on Monte Carlo simulations of the response of a single stand-alone 30 m imaging atmospheric Cherenkov telescope (IACT) the maximal rejection power against background cosmic ray showers for low energy gamma-rays was investigated in great detail. An advanced Bayesian multivariate analysis has been applied to the simulated Cherenkov light images of the gamma-ray- and proton-induced air showers. The results obtained here quantitatively testify that the separation between the signal and background images degrades substantially at low energies, and consequently the maximum overall quality factor can only be about 3.1 for gamma rays in the 10-30 GeV energy range. Various selection criteria as well as optimal combinations of the standard image parameters utilized for effective image separation have been also evaluated.

We report on observations of the sky region around the unidentified TeV gamma-ray source TeV J2032+4130 carried out with the Whipple Observatory 10 m atmospheric Cherenkov telescope for a total of 65.5 hrs between 2003 and 2005. The standard two-dimensional analysis developed by the Whipple collaboration for a stand-alone telescope reveals an excess in the field of view at a pre-trials significance level of 6.1 standard deviations. The measured position of this excess is alpha(2000) =20 h 32 m 27 s, delta(2000) = 41 deg 39 min 17 s. The estimated integral flux for this gamma-ray source is about 8% of the Crab-Nebula flux. The data are consistent with a point-like source. Here we present a detailed description of the standard two-dimensional analysis technique used for the analysis of data taken with the Whipple Observatory 10 m telescope and the results for the TeV J2032+4130 campaign. We include a short discussion of the physical mechanisms that may be responsible for the observed gamma-ray emission, based on possible association with known astrophysical objects, in particular Cygnus OB2.

We present high-resolution X-ray images taken with the Chandra X-ray Observatory of the field that contains the unidentified TeV gamma-ray source HESS J1804-216. A total of eleven discrete sources were detected with a posteriori significance of > 5 sigma over the entire field of view. Among them, only one, designated as CXOU J180351.4-213707, is significantly extended. The source is about 40" away from the radio pulsar PSR J1803-2137, which was the target of the Chandra observation but was not detected in X-rays. A natural question is whether the two sources are physically related. While it is conceivable that CXOU J180351.4-213707 could be associated with a previously unknown supernova remnant (SNR), in which the pulsar was born, it seems equally plausible that it might be a pulsar wind nebula (PWN) that is powered by a different pulsar whose emission is beamed away from us. In either case, we argue that CXOU J180351.4-213707 is likely the X-ray counterpart of HESS J1804-216, based on the fact that the Galactic TeV gamma-ray sources are predominantly SNRs or PWNe. The X-ray spectrum of the source can be fitted well with a power law, although the model is not well constrained due to large statistical uncertainties. The spectrum seems to be very hard, with the best-fit photon index about 1.2. Under the assumption that CXOU J180351.4-213707 is the X-ray counterpart of HESS J1804-216, we attempted to model the X-ray and TeV emission as synchrotron and inverse Compton scattered radiation from relativistic electrons. We briefly discuss the results.

The H.E.S.S. collaboration recently performed the accurate measurements of the gamma-ray spectra above 200 GeV for two BL Lac-type AGN in the Southern sky: PKS 2155-304 and PKS 2005-489. The TeV spectrum of the BL Lac object 1ES 2344+514 has also been recently measured by the Whipple collaboration. Using the results of phenomenological calculations of the SED of the EBL we have extracted the intrinsic gamma-ray spectra of these blazars. Consequently we have shown that these spectra can be fitted well by means of a synchrotron self Compton model.

Pulsars are commonly regarded as highly magnetized neutron stars, rotating up to several hundred times per second. Over 1,500 radio pulsars have been found so far, about 70 of which are X-ray pulsars, but only a handful have been observed in gamma rays. This high-energy emission is believed to be produced by the electrons accelerated to TeV energies in the pulsar magnetosphere, resulting in cascades of secondary particles. The H.E.S.S. experiment is a system of four imaging atmospheric Cherenkov telescopes in Namibia, which is sensitive to gamma rays above 100 GeV. Three young pulsars, Crab, Vela, and PSR B1706-44 have been observed with H.E.S.S. The results of the search for pulsed emission for these targets, and the constraints on the theories of pulsed very high energy gamma-ray emission, are summarized here.

The construction of the H.E.S.S. (High Energy Stereoscopic System), a superior system of four 12 m imaging atmospheric Cherenkov telescopes, has been completed recently in the Namib desert close to Windhoek, at 1800 m above sea level. This new generation of ground-based gamma-ray detectors has an energy threshold of about 100 GeV in observations at zenith and a sensitivity of about 1% Crab flux for a point-like gamma-ray source. Such high sensitivity was achieved due to 0.1 degree angular resolution, as well as, severe cosmic ray background rejection acquired from multi-fold imaging of individual atmospheric showers. In addition H.E.S.S. has a rather good energy resolution of 15%. H.E.S.S. has been taking routine observations of the gamma-ray sources since December 2003. The outstanding physics results achieved with H.E.S.S. already in a first year of its exploitation strongly encourage further development of imaging atmospheric Cherenkov technique for high-quality gamma-ray observations, which is basically driven by father reduction of the energy threshold of a forthcoming major future gamma-ray detector. Here we are dealing with such a detector. Basic results on performance and sensitivity for a single stand-alone 30 m imaging atmospheric Cherenkov telescope, as well as for a system of two and five telescopes, are discussed here.

At present the ground-based Very High Energy (VHE) gamma-ray astronomy is racing to complete construction of a number of modern gamma-ray detectors, i.e. CANGAROO III, MAGIC, H.E.S.S., and VERITAS. They should be fully operational in a year's time. After much debate, the further development of this gamma-ray astronomy in the foreseeable future must be widely anticipated to proceed with the designing and building of a new instrumentation, which is primarily intended for the further drastic reduction of the energy threshold in gamma-ray observations down to about 10 GeV. On the ground one can hardly reach such low energy thresholds without considerably larger, up to 30 m diameter, optical telescopes, which might be able to collect sufficient amount of Cherenkov light from the atmospheric gamma-ray showers of that low energy. If not taken off the ground entirely (like GLAST), then it seems to be profitable to mount future low energy Cherenkov telescopes at higher altitudes in the atmosphere in order that they will be able to detect substantially more unabsorbed Cherenkov light from a shower. There are a few sites up on the high mountains of roughly 5 km height worldwide, which can be used for such a venture. However, one has to remember that actual time profile, and in particular a two-dimensional distribution (image) of Cherenkov light flash from atmospheric gamma-ray showers, undergo a rapid change after an increase in the observational level. This paper briefly describes the results of a topological analysis of Cherenkov light images calculated at both conventional and desirable altitudes of 2.2 and 5 km above the sea, respectively. A discussion on major upgrades of image topology at high altitude is also given.

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.

Routine observations of the Crab Nebula for a total of about 250 hrs, performed with the HEGRA stereoscopic system of 5 imaging atmospheric Cerenkov telescopes in the standard operational mode, have proven the energy threshold of the system to be 500 GeV for small zenith angles (< 20 deg). A topological trigger applied along with the convergent observational mode allows to reduce noticeably the energy threshold of the system down to 350 GeV. Here we present the relevant Monte Carlo simulations as well as the analysis results of 15 hrs Crab Nebula data taken in such an observational mode. From the Crab Nebula data, the final energy threshold was found to be 350 GeV. The estimated gamma ray flux from the Crab Nebula above 350 GeV is (8.1 +- 0.1(stat) +- 0.2(syst)) 10^(-11) ph cm^(-2) s^(-1), which is consistent with recent measurements reported by the STACEE, CELESTE, CAT, and Whipple groups.

In this paper, we present the results of Monte Carlo simulations of atmospheric showers induced by diffuse gamma rays as detected by the high-energy gamma ray astronomy (HEGRA) system of five imaging atmospheric Cerenkov telescopes (IACTs). We have investigated the sensitivity of observations on extended gamma ray emission over the entire field of view of the instrument. We discuss a technique to search for extended gamma ray sources within the field of view of the instrument. We give estimates for HEGRA sensitivity of observations on extended TeV gamma ray sources.

Since their discovery at radio wavelengths pulsars have been persistent targets for widespread multi-wave observations throughout optics, radio, X-rays, and high-energy gamma-rays. Observations with the EGRET gamma-ray telescope, on board Compton GRO satellite, confirmed the expectation of a pulsed high-energy emission up to a few GeV. Presently, at least seven objects are known as well established high-energy gamma-ray pulsars. A few of those emit gamma-rays well above 1 GeV. Forthcoming ground-based Cherenkov telescopes will enable observations of gamma-rays well below 100 GeV, finally reaching the yet unexplored energy gap at tens of GeV. H.E.S.S. (High Energy Stereoscopic System) is one of such instruments which is planned to be operational in 2004. Here I summarize the basic scientific motivations, the H.E.S.S. sensitivity, and the first targets for future pulsar observations at high energies from the ground.

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.

The HEGRA-IACT-system with a FoV of ~1.5 deg radius has been used for surveying one quater of the Galactic disc in respect to point sources, moderately extended sources and for diffuse emission in the energy range above 1 TeV. In total 140 h of good observation time were accumulated. No new source has been discoverd. Limits on the level of 20% or lower of the Crab flux on about 87 potential sources like SNR, Pulsars and EGRET sources are derived. A limit on the diffuse emission is given on the level of dPhi/dE(E=1 TeV) = 6.1 10^-15 ph cm-2 s-1 sr-1 Mev-1 resulting in a lower limit of 2.5 on the spectral index for the extrapolation of the meassured EGRET flux for the diffuse emission.

At present, the ground-based astronomy of very high energy (VHE) (E>100 GeV) gamma-rays experiences the challenging transition phase caused by the substantial upgrading of its observational instrumentation. Recently the HEGRA collaboration has evidently demonstrated the power of arrays of imaging atmospheric Cherenkov telescopes and initiated a number of complementary projects, such as H.E.S.S., CANGAROO IV, VERITAS. The common philosophy of those projects is based on the stereoscopic observations of air showers with a number of imaging telescopes simultaneously. Such observations allow to reduce substantially the effective energy threshold of detected gamma-ray showers; to improve the angular and energy resolution for individual gamma-rays; to gain in suppression of the cosmic ray background. In my talk I have summarized the general advantages and achieved sensitivity of the stereoscopic observations using the HEGRA data taken with an array of 5 imaging telescopes having the energy threshold of 500 GeV. These results completed with the relevant simulations may be extrapolated on a firm ground towards the future H.E.S.S. (phase I) array of four 12 m imaging telescopes with the estimated energy threshold below 100 GeV. Here I discuss in short the anticipated performance of such an array.

We discuss the concept and the performance of a powerful future ground-based astronomical instrument - a stereoscopic array of several large imaging atmospheric Cherenkov telescopes installed at a very high mountain elevation of about 5 km a.s.l. or more - for the study of the gamma-ray sky at energies from several GeV to 100 GeV.

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

Ground-based atmospheric Cherenkov telescopes are proven to be effective instruments for observations of very high energy (VHE) gamma-radiation from celestial objects. For effective use of such technique one needs detailed Monte Carlo simulations of gamma-ray- and proton/nuclei-induced air showers in Earth atmosphere. Here we discuss in detail the algorithms used in the numerical code ALTAI developed particularly for the simulations of Cherenkov light emission from air showers of energy below 50 TeV. The specific scheme of sampling the charged particle transport in the atmosphere allows the performance of very fast and accurate simulations used for interpretation of the VHE gamma-ray observations.

We discuss an effect of the intergalactic absorption of the TeV gamma-rays in time-averaged spectrum of Mkn 501 measured by the HEGRA Collaboration. Analysis of the spectral behavior, variability time scale and relevant calculations of TeV gamma-ray emission allow to conclude the presence of a noticeable absorption of the TeV gamma-rays in the Mkn 501 energy spectrum.

For stereoscopic systems of imaging atmospheric Cherenkov telescopes (IACTs), a key parameter to optimize the sensitivity for VHE gamma-ray point sources is the intertelescope spacing. Using pairs of telescopes of the HEGRA IACT system, the sensitivity of two-telescope stereo IACT systems is studied as a function of the telescope spacing, ranging from 70 m to 140 m. Data taken during the 1997 outburst of Mrk 501 are used to evaluate both the detection rates before cuts, and the sensitivity for weak signals after cuts to optimize the significance of signals. Detection rates decrease by about 1/3 between the minimum and maximum spacing. The significance of signals is essentially independent of distance, in the range investigated.

We present analysis techniques to improve the energy resolution of stereoscopic systems of imaging atmospheric Cherenkov telescopes, using the HEGRA telescope system as an example. The techniques include (i) the determination of the height of the shower maximum, which is then taken into account in the energy determination, and (ii) the determination of the location of the shower core with the additional constraint that the direction of the gamma rays is known a priori. This constraint can be applied for gamma-ray point sources, and results in a significant improvement in the localization of the shower core, which translates into better energy resolution. Combining both techniques, the HEGRA telescopes reach an energy resolution between 9% and 12%, over the entire energy range from 1 TeV to almost 100 TeV. Options for further improvements of the energy resolution are discussed.

The sensitivity of imaging atmospheric Cherenkov telescopes (IACTs) in TeV gamma-ray observations reachs its maximum at small zenith angles (< 30 degree) which provide the minimum attainable energy threshold of an instrument. However, for a specific telescope site a number of gamma-ray sources, or source candidates, can only be observed at much larger zenith angles (< 60 degree). Moreover the observations at large zenith angles allow to extend the observation time window for any object seen at small zenith angles, as well as to enlarge the dynamic energy range of an instrument towards the highest observable energies of gamma-rays. Based on Monte Carlo simulations we present here the results on the sensitivity of a stereoscopic system of 5 IACTs in observations at large zenith angles. We point out some important parameters of the telescope design which could substantially improve the efficiency of such observations with forthcoming IACT arrays like CANGAROO III, HESS and VERITAS.

For air showers observed simultaneously by more than two imaging atmospheric Cherenkov telescopes, the shower geometry is overconstrained by the images and image information should be combined taking into account the quality of the images. Different algorithms are discussed and tested experimentally using data obtained from observations of Mkn 501 with the HEGRA IACT system. Most of these algorithms provide an estimate of the accuracy of the reconstruction of shower geometry on an event-by-event basis, allowing, e.g., to select higher-quality subsamples for precision measurements.

The recent HEGRA observations of the blazar Mkn 501 show strong curvature in the very high energy gamma-ray spectrum. Applying the gamma-ray opacity derived from an empirically based model of the intergalactic infrared background radiation field (IIRF), to these observations, we find that the intrinsic spectrum of this source is consistent with a power-law: dN/dE~ E^-alpha with alpha=2.00 +/- 0.03 over the range 500 GeV - 20 TeV. Within current synchrotron self-Compton scenarios, the fact that the TeV spectral energy distribution of Mkn 501 does not vary with luminosity, combined with the correlated, spectrally variable emission in X-rays, as observed by the BeppoSAX and RXTE instruments, also independently implies that the intrinsic spectrum must be close to alpha=2. Thus, the observed curvature in the spectrum is most easily understood as resulting from intergalactic absorption.

Arrays of imaging air Cherenkov telescopes (IACTs) like VERITAS, HESS have been recently proposed as the instruments of the next generation for ground based very high energy gamma-ray astronomy invading into 50-100 GeV energy range. Here we present results of design studies for the future IACT arrays which have been performed by means of Monte Carlo simulations. We studied different trigger strategies, abilities of cosmic ray rejection for arrays of 4 and 16 telescopes with 10 m reflectors, equipped with cameras comprising 271 and 721 pixels of 0.25 and 0.15 degree, respectively. The comparative analysis of the performance of such telescope arrays has been done for both camera options, providing almost the same field of view of 4.3 degree. An important issue is the choice of the optimum spacing between the telescopes in such an array. In order to maximize the signal-to-noise ratio in observations at the small zenith angles of 20 degree as well as at large zenith angles of 60 degree, different arrangements of IACT array have been examined. Finally, we present a major recommendations regarding the optimum configuration.

Based on the Monte Carlo simulations we have studied the performance of the HEGRA system of imaging air Cherenkov telescopes (IACTs) in its present configuration of 4 IACTs as well as in its future final configuration of 5 IACTs. Here we present the results on the basic characteristics of the IACT system which are used in the standard data analysis procedure, i.e., the collection areas, the detection rates, the angular resolution, the energy resolution, and the gamma/hadron-separation efficiency. By comparing several key Monte Carlo predictions with experimental results it is possible to check the accuracy of the simulations. The Monte Carlo results concerning hadron-nuclear showers are tested with the recorded cosmic ray events and the results concerning photon-induced showers are tested with a large data sample of gamma-rays observed from BL Lac object Mkn 501 during its high flaring activity in 1997. Summarizing the simulations and current observations we give the basic recommendations of using the instrument and the major values of its sensitivity.

The BL Lac object Mkn 501 has shown very high emission in TeV gamma-rays from March to October, 1997. During this period the source was continuously monitored with the HEGRA stereoscopic system of 4 imaging air Cherenkov telescopes for a total exposure time of 110 hours. The unprecedented statistics of about 38,000 TeV photons, combined with the good energy resolution of ~20 % over the entire energy range and with detailed simulations of the detector response, allowed a determination of the average energy spectrum from 500 GeV to 24 TeV. Although the gamma-ray flux varied strongly with time, the daily energy spectra remained rather constant in their shape. Therefore it is justified to derive a time averaged spectrum and thus to extend, for the first time, the spectral measurement well beyond 10 TeV. This TeV spectrum of Mkn 501 shows a gradual steepening which could be caused by a number of physical processes, such as a limited energy range of the radiating particles, intrinsic gamma-ray absorption inside the source, the Klein-Nishina effect (in the case of an inverse Compton origin of the radiation) and finally, an absorption of the TeV gamma-rays propagating in the intergalactic medium.

The Crab Nebula has been observed by the HEGRA stereoscopic system of 4 imaging air Cherenkov telescopes (IACTs) for about 100 hours from 1997 September to 1998 March. The recent detailed studies on the system performance give an energy threshold and an energy resolution for gamma-rays of 500 GeV and \~20 %, respectively. The Crab spectrum was measured with the IACTs system in a very broad energy range from 500 GeV to 25 TeV, using the observations at zenith angles up to 50 degree. The pure power-law spectrum in TeV gamma-rays constrains the physics parameters of the nebula environment as well as the models of photon emission.

We describe the features of the stereoscopic approach in the imaging atmospheric Cherenkov technique, and discuss the performance of future low threshold telescope arrays.

The HEGRA 5 m^2 air Cherenkov telescope (CT1) was used to search for gamma-ray emission above 1.5 TeV from a series of low redshift AGN (Mkn 421, Mkn 501, MS 0116+319, PKS 2209+236, NGC 315 and W Comae). Here we present results from a total of 383 hours of ON-source observations at zenith angles up to 35$^\circ$ made between February 1996 and April 1997 showing positive detections of Mkn 421 and Mkn 501 and upper limits on the other objects' emission. More recent results from observations of Mkn 421 und Mkn 501 will be added at the conference.

A detection of TeV $\gamma$-rays from Mkn 501 is reported, based on observations made between March and August 1996 with the first HEGRA Cherenkov telescope (CT1). From the image analysis, 351 excess candidate $\gamma$-ray events are obtained from the 147 h dataset. The statistical significance of the excess is 5.2 $\sigma$. The average excess rate is $2.4 \pm 0.5$ h$^{-1}$ above the $\approx$ 1.5 TeV threshold of CT1. Under the assumption that the spectrum of Mkn 501 follows a power law we find a differential spectral index of 2.6$\pm0.5$ and obtain a time-averaged integral flux above 1.5 TeV of $2.3 (\pm0.4)_{Stat} (+1.5-0.6)_{Syst} \times 10^{-12}$ cm$^{-2}$s$^{-1}$. Comparison with our near contemporary observations of the Crab Nebula, used as a standard candle to test CT1 after upgrading to a 127 pixel camera, indicates that Mkn 501 has a spectrum similar to that of the Crab Nebula above 1.5 TeV. The integral flux above 1.5 TeV from Mkn 501 is found to have been between 2.2 and 3.6 times smaller than that from the Crab Nebula. HEGRA is the second experiment to have detected Mkn 501 in the TeV range.

A detection of Gamma-rays from Mkn 421 at energies above 1 TeV is reported, based on observations made in December 1994 - May 1995 with the first two HEGRA Cherenkov telescopes. >From the image analysis, 111 excess gamma candidates are obtained from the 26 h Telescope #1 (CT1) dataset (significance 4.0 sigma) and 218 from the 41 h Telescope #2 (CT2) dataset (significance 4.2 sigma) at zenith angles theta < 25 degrees. The combined significance is approx. 5.8 sigma. This is the second detection of Mkn 421 at TeV energies. The average excess rate is 4.3 +- 1.0 h^-1 for CT1 and 5.4 +- 1.3 h^-1 for CT2. Comparison with our contemporary observations of the Crab Nebula indicates that Mkn 421 has a steeper spectrum than the Crab Nebula above 1 TeV. Under the assumption that the spectrum of Mkn 421 follows a power law, we obtain a differential spectral index of 3.6 +- 1.0 and an integral flux above 1 TeV of 8 (+-2)_(Stat) (+6-3)_(Syst) x 10^-12 cm^-2 s^-1 from a comparison with Monte Carlo data. This flux is smaller than the Crab Nebula flux by a factor of 2.0 +- 0.8.