The space-time correlations of streams of photons can provide fundamentally new channels of information about the Universe. Today's astronomical observations essentially measure certain amplitude coherence functions produced by a source. The spatial correlations of wave fields has traditionally been exploited in Michelson-style amplitude interferometry. However the technology of the past was largely incapable of fine timing resolution and recording multiple beams. When time and space correlations are combined it is possible to achieve spectacular measurements that are impossible by any other means. Stellar intensity interferometry is ripe for development and is one of the few unexploited mechanisms to obtain potentially revolutionary new information in astronomy. As we discuss below, the modern use of stellar intensity interferometry can yield unprecedented measures of stellar diameters, binary stars, distance measures including Cepheids, rapidly rotating stars, pulsating stars, and short-time scale fluctuations that have never been measured before.

Intensity interferometry exploits a quantum optical effect in order to measure objects with extremely small angular scales. The first experiment to use this technique was the Narrabri intensity interferometer, which was successfully used in the 1970s to measure 32 stellar diameters at optical wavelengths; some as small as 0.4 milli-arcseconds. The advantage of this technique, in comparison with Michelson interferometers, is that it requires only relatively crude, but large, light collectors equipped with fast (nanosecond) photon detectors. Ground-based gamma-ray telescope arrays have similar specifications, and a number of these observatories are now operating worldwide, with more extensive installations planned for the future. These future instruments (CTA, AGIS, completion 2015) with 30-90 telescopes will provide 400-4000 different baselines that range in length between 50m and a kilometre. Intensity interferometry with such arrays of telescopes attains $50 \mu$-arcsecond resolution for a limiting visual magnitude ~8.5. Phase information can be extracted from the interferometric measurement with phase closure, allowing image reconstruction. This technique opens the possibility of a wide range of studies amongst others, probing the stellar surface activity and the dynamic AU scale circumstellar environment of stars in various crucial evolutionary stages. Here we focuse on the astrophysical potential of an intensity interferometer utilising planned new gamma-ray instrumentation.

The present generation of ground-based Very High Energy (VHE) gamma-ray observatories consist of arrays of up to four large (> 12m diameter) light collectors quite similar to those used by R. Hanbury Brown to measure stellar diameters by Intensity Interferometry in the late 60's. VHE gamma-ray observatories to be constructed over the coming decade will involve several tens of telescopes of similar or greater sizes. Used as intensity interferometers, they will provide hundreds of independent baselines. Now is the right time to re-assess the potential of intensity interferometry so that it can be taken into consideration in the design of these large facilities.

In the 1970s, the Narrabri intensity interferometer was used to measure 32 stellar diameters; some as small as 0.4 milli-arc-seconds (mas). The interferometer consisted of a pair of 6.5m telescopes with relatively crude optics, similar to those currently in use as Atmospheric Cherenkov Telescopes (ACT). We explore the possibility of implementing a modern intensity interferometer on an ACT array. Developments in fast digital signal processing technology now make such a system relatively easy to implement, and provide improved sensitivity. Allowing measurements at short wavelength (<400nm), with long baselines (> 100m), which are still challenging for Michelson interferometers, present ACT arrays could be used to probe angular structures as small as ~0.2mas, and smaller with large array projects already being discussed. This would provide measurements of stellar diameters, binary systems, circumstellar environments and, possibly, stellar surface features. ACT arrays could be used as intensity interferometers during bright moon periods, providing valuable scientific output for little expense and no impact on the gamma-ray observing schedule.

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive gamma-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 100,000 square meter, combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic Ray initiated air showers that resemble gamma-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of sigma[63%]= 0.1 degree (0.2 degree) and sigma[68%]= 15% (22%), respectively. The main emphasis of the paper is the discussion of gamma-hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a gamma-hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.

In April-May 2003, the blazar 1ES 1959+650 showed an increased level of X-ray activity. This prompted a multiwavelength observation campaign with the Whipple 10 m gamma-ray telescope, the Rossi X-ray Timing Explorer, the Bordeaux Optical Observatory, and the University of Michigan Radio Astrophysical Observatory. We present the multiwavelength data taken from May 2, 2003 to June 7, 2003 and compare the source characteristics with those measured during observations taken during the years 2000 and 2002. The X-ray observations gave a data set with high signal-to-noise light curves and energy spectra; however, the gamma-ray observations did not reveal a major TeV gamma-ray flare. Furthermore, we find that the radio and optical fluxes do not show statistically significant deviations from those measured during the 2002 flaring periods. While the X-ray flux and X-ray photon index appear correlated during subsequent observations, the apparent correlation evolved significantly between the years 2000, 2002, and 2003. We discuss the implications of this finding for the mechanism that causes the flaring activity.

We present a method of atmospheric Cherenkov imaging which reconstructs the unique arrival direction of TeV gamma rays using a single telescope. The method is derived empirically and utilizes several features of gamma-ray induced air showers which determine, to a precision of 0.12 degrees, the arrival direction of photons, on an event-by-event basis. Data from the Whipple Observatory's 10 m gamma-ray telescope is utilized to test selection methods based on source location. The results compare these selection methods with traditional techniques and three different camera fields of view. The method will be discussed in the context of a search for a gamma-ray signal from a point source located anywhere within the field of view and from regions of extended emission.

The BL Lac object Mrk 501 went into a very high state of activity during 1997, both in VHE gamma-rays and X-rays. We present here results from observations at energies above 250 GeV carried out between March and October 1997 with the CAT Cerenkov imaging Telescope. The average differential spectrum between 30 GeV and 13 TeV shows significant curvature and is well represented by phi_0 * E_TeV^{-(alpha + beta*log10(E_TeV))}, with: phi_0 = 5.19 +/- 0.13 {stat} +/- 0.12 {sys-MC} +1.66/-1.04 {sys-atm} * 10^-11 /cm^2/s/TeV alpha = 2.24 +/- 0.04 {stat} +/- 0.05 {sys} beta = 0.50 +/- 0.07 {stat} (negligible systematics). The TeV spectral energy distribution of Mrk 501 clearly peaks in the range 500 GeV-1 TeV. Investigation of spectral variations shows a significant hardness-intensity correlation with no measurable effect on the curvature. This can be described as an increase of the peak TeV emission energy with intensity. Simultaneous and quasi-simultaneous CAT VHE gamma-ray and BeppoSAX hard X-ray detections for the highest recorded flare on 16th April and for lower-activity states of the same period show correlated variability with a higher luminosity in X-rays than in gamma-rays. The observed spectral energy distribution and the correlated variability between X-rays and gamma-rays, both in amplitude and in hardening of spectra, favour a two-component emission scheme where the low and high energy components are attributed to synchrotron and inverse Compton (IC) radiation, respectively.

Gamma-ray observations above 200 MeV are conventionally made by satellite-based detectors. The EGRET detector on the Compton Gamma Ray Observatory (CGRO) has provided good sensitivity for the detection of bursts lasting for more than 200 ms. Theoretical predictions of high-energy gamma-ray bursts produced by quantum-mechanical decay of primordial black holes (Hawking 1971) suggest the emission of bursts on shorter time scales. The final stage of a primordial black hole results in a burst of gamma-rays, peaking around 250 MeV and lasting for a tenth of a microsecond or longer depending on particle physics. In this work we show that there is an observational window using ground-based imaging Cherenkov detectors to measure gamma-ray burst emission at energies E greater than 200 MeV. This technique, with a sensitivity for bursts lasting nanoseconds to several microseconds, is based on the detection of multi-photon-initiated air showers.

Only BL Lac objects have been detected as extragalactic sources of very high energy (E > 300 GeV) gamma rays. Using the Whipple Observatory Gamma-ray Telescope, we have attempted to detect more BL Lacs using three approaches. First, we have conducted surveys of nearby BL Lacs, which led to the detections of Mrk 501 and 1ES 2344+514. Second, we have observed X-ray bright BL Lacs when the RXTE All-Sky Monitor identifies high state X-ray emission in an object, in order to efficiently detect extended high emission states. Third, we have conducted rapid observations of several BL Lacs and QSOs located close together in the sky to search for very high flux, short time-scale flare states such as have been seen from Mrk 421. We will present the results of a survey using the third observational technique.

Markarian 421 was observed for about four days with BeppoSAX and the Whipple Observatory gamma-ray telescope in April 1998. A pronounced, well-defined, flare with hour-scale variability was observed simultaneously in X-rays and very high energy gamma-rays. These data provide the first evidence that the X-ray and TeV intensities are well correlated on time-scales of hours. While the rise of the flare occurred on a similar time-scale in the two wavebands, the decay of the flare was much more rapid in gamma rays, providing the first clear indication that the X-ray and gamma-ray emission may not be completely correlated in Markarian 421.

If cosmic rays with energies <100 TeV originate in the galaxy and are accelerated in shock waves in shell-type supernova remnants (SNRs), gamma-rays will be produced as the result of proton and electron interactions with the local interstellar medium, and by inverse Compton emission from electrons scattering soft photon fields. We report on observations of two supernova remnants with the Whipple Observatory's 10 m gamma-ray telescope. No significant detections have been made and upper limits on the >500 GeV flux are reported. Non-thermal X-ray emission detected from one of these remnants (Cassiopeia A) has been interpreted as synchrotron emission from electrons in the ambient magnetic fields. Gamma-ray emission detected from the Monoceros/Rosette Nebula region has been interpreted as evidence of cosmic-ray acceleration. We interpret our results in the context of these observations.

A new method of shower-image analysis is presented which appears very powerful as applied to those Cherenkov Imaging Telescopes with very high definition imaging capability. It provides hadron rejection on the basis of a single cut on the image shape, and simultaneously determines the energy of the electromagnetic shower and the position of the shower axis with respect to the detector. The source location is also reconstructed for each individual gamma-ray shower, even with one single telescope, so for a point source the hadron rejection can be further improved. As an example, this new method is applied to data from the CAT (Cherenkov Array at Themis) imaging telescope, which has been operational since Autumn, 1996.

The CAT (Cherenkov Array at Themis) imaging telescope, equipped with a very-high-definition camera (546 fast phototubes with 0.12 degrees spacing surrounded by 54 larger tubes in two guard rings) started operation in Autumn 1996 on the site of the former solar plant Themis (France). Using the atmospheric Cherenkov technique, it detects and identifies very high energy gamma-rays in the range 250 GeV to a few tens of TeV. The instrument, which has detected three sources (Crab nebula, Mrk 421 and Mrk 501), is described in detail.

The CAT Imaging Telescope has observed the BL Lac object Markarian 501 between March and August 1997. We report here on the variability over this time including several large flares. We present also preliminary spectra for all these data, for the low emission state, and for the largest flare.

The CAT imaging telesope, which uses the Atmospheric Cherenkov Technique to observe TeV gamma-ray sources, has been operating since September 1996. Located in southern France, it features a fine grained camera consisting of 546 PMT with a pixel size of 0.12 degrees. Strong gamma-ray signals from Mkn 501 above 200 GeV were detected during observations made in April 1997. The high signal/noise ratio for these data enables direct characterization of the gamma-ray images as recorded in the camera. These results are used to validate the Monte-Carlo simulations of the response of the telescope to gamma-ray showers and to improve its performances. Observational results on Mkn 180 and the Crab nebula spectrum are presented.

The VHE gamma-ray imaging telescope CAT started taking data in October 1996. Located at the Themis solar site in southern France, it features a 17.7 m^2 Davies-Cotton mirror equipped with 600 PMT camera at the focal plane. The mechanics and optics, the PMTs and the electronics are presented. The performance, based on the first 7 months of operation, is discribed.

The CAT imaging telescope, recently built on the site of the former solar plant Themis (French Pyrenees), observed gamma-rays from the Crab nebula from October 1996 to March 1997. This steady source, often considered as the standard candle of very-high-energy gamma-ray astronomy, is used as a test-beam to probe the performances of the new telescope, particularly its energy threshold (220 GeV at 20 degrees zenith angle) and the stability of its response. Due to the fine-grain camera, an accurate analysis of the longitudinal profiles of shower images is performed, yielding the source position in two dimensions for each individual shower.

The CAT imaging telescope on the site on the former solar plant Themis has been observing gamma-rays from Mrk501 above 220 GeV in March and April 1997. This source is shown to be highly variable and the light curve is presented. The detected gamma-ray rate for the most intense flare is in excess of 10 per minute.