Imaging Atmospheric Cherenkov Telescopes have long been viewed as potential light collectors to be used for long baseline optical intensity interferometry observations. Intensity interferometry, as implemented with Cherenkov telescopes, is well suited for studying the spatial structure of stars of O/B/A stellar types at short optical wavelengths. Such observations complement those with the current generation of optical amplitude interferometers, which are typically restricted to longer wavelengths. Dedicated intensity interferometry instrumentation has been developed for the VERITAS observatory with engineering tests and observations occurring since October 2018. Here, the first results using two of the VERITAS telescopes are reported. The system has already been extended to the two additional telescopes, enabling SII observations with all four VERITAS telescopes.

In 2018-2019 the VERITAS VHE gamma-ray observatory was augmented with highspeed optical instrumentation and continuous data recording electronics to create a sensitive Stellar Intensity Interferometry (SII) observatory, VERITAS-SII. The primary science goal of VERITAS-SII is to perform stellar diameter measurements and image analysis in the visible wavebands on a selection of bright (m< 6), hot (O/B/A) stars. The VERITAS Collaboration has agreed to the deployment and operation of VERITAS-SII during several days each month around the full moon period when VERITAS does not perform VHE gamma-ray observations. The VERITAS-SII augmentation employs custom high-speed/low-noise focal plane instrumentation using high quantum efficiency photomultiplier tubes, and a battery-powered, fiber-optic controlled High Voltage supply. To reduce engineering time, VERITAS-SII uses commercially available high-speed (250 MS/sec), continuously streaming electronics to record the time dependence of the intensity fluctuations at each VERITAS telescope. VERITAS-SII also uses fast ( < 100 psec) data acquisition clock synchronization over inter-telescope distances (greater than 100 m) using a commercially available White Rabbit based timing solution. VERITAS-SII is now in full operation at the VERITAS observatory, F.L.Whipple Observatory, Amado, AZ USA. This paper describes the design of the instrumentation hardware used for VERITAS-SII augmentation of the VERITAS observatory, the status of initial VERITAS-SII observations, and plans for future improvements to VERITAS-SII.

Recent advances in telescope design, photodetector efficiency, and high-speed electronic data recording and synchronization have created the observational capability to achieve unprecedented angular resolution for several thousand bright (m< 6) and hot (O/B/A) stars by means of a modern implementation of Stellar Intensity Interferometry (SII). This technology, when deployed on future arrays of large diameter optical telescopes, has the ability to image astrophysical objects with an angular resolution better than 40 {\mu} arc-sec. This paper describes validation tests of the SII technique in the laboratory using various optical sensors and correlators, and SII measurements on nearby stars that have recently been completed as a technology demonstrator. The paper describes ongoing and future developments that will advance the impact and instrumental resolution of SII during the upcoming decade.

We present measurements of the second order spatial coherence function of thermal light sources using Hanbury-Brown and Twiss interferometry with a digital correlator. We demonstrate that intensity fluctuations between orthogonal polarizations, or at detector separations greater than the spatial coherence length of the source, are uncorrelated but can be used to reduce systematic noise. The work performed here can readily be applied to existing and future Imaging Air-Cherenkov Telescopes used as star light collectors for Stellar Intensity Interferometry (SII) to measure spatial properties of astronomical objects.

Stellar intensity interferometers will achieve stellar imaging with a tenth of a milli- arcsecond resolution in the optical band by taking advantage of the large light collect- ing area and broad range of inter-telescope distances offered by future gamma-ray Air Cherenkov Telescope (ACT) arrays. Up to now, studies characterizing the capabilities of intensity interferometers using ACTs have not accounted for realistic effects such as telescope mirror extension, detailed photodetector time response, excess noise, and night sky contamination. In this paper, we present the semi-classical quantum optics Monte-Carlo simulation we developed in order to investigate these experimental limi- tations. In order to validate the simulation algorithm, we compare our first results to models for sensitivity and signal degradation resulting from mirror extension, pulse shape, detector excess noise, and night sky contamination.

We present the results of a joint observational campaign between the Green Bank radio telescope and the VERITAS gamma-ray telescope, which searched for a correlation between the emission of very-high-energy (VHE) gamma rays ($E_{\gamma} >$ 150 GeV) and Giant Radio Pulses (GRPs) from the Crab pulsar at 8.9 GHz. A total of 15366 GRPs were recorded during 11.6 hours of simultaneous observations, which were made across four nights in December 2008 and in November and December 2009. We searched for an enhancement of the pulsed gamma-ray emission within time windows placed around the arrival time of the GRP events. In total, 8 different time windows with durations ranging from 0.033 ms to 72 s were positioned at three different locations relative to the GRP to search for enhanced gamma-ray emission which lagged, led, or was concurrent with, the GRP event. Further, we performed separate searches on main pulse GRPs and interpulse GRPs and on the most energetic GRPs in our data sample. No significant enhancement of pulsed VHE emission was found in any of the preformed searches. We set upper limits of 5-10 times the average VHE flux of the Crab pulsar on the flux simultaneous with interpulse GRPs on single-rotation-period time scales. On $\sim$8-second time scales around interpulse GRPs, we set an upper limit of 2-3 times the average VHE flux. Within the framework of recent models for pulsed VHE emission from the Crab pulsar, the expected VHE-GRP emission correlations are below the derived limits.

Using kilometric arrays of air Cherenkov telescopes, intensity interferometry may increase the spatial resolution in optical astronomy by an order of magnitude, enabling images of rapidly rotating stars with structures in their circumstellar disks and winds, or mapping out patterns of nonradial pulsations across stellar surfaces. Intensity interferometry (pioneered by Hanbury Brown and Twiss) connects telescopes only electronically, and is practically insensitive to atmospheric turbulence and optical imperfections, permitting observations over long baselines and through large airmasses, also at short optical wavelengths. The required large telescopes with very fast detectors are becoming available as arrays of air Cherenkov telescopes, distributed over a few square km. Digital signal handling enables very many baselines to be synthesized, while stars are tracked with electronic time delays, thus synthesizing an optical interferometer in software. Simulated observations indicate limiting magnitudes around m(v)=8, reaching resolutions ~30 microarcsec in the violet. The signal-to-noise ratio favors high-temperature sources and emission-line structures, and is independent of the optical passband, be it a single spectral line or the broad spectral continuum. Intensity interferometry provides the modulus (but not phase) of any spatial frequency component of the source image; for this reason image reconstruction requires phase retrieval techniques, feasible if sufficient coverage of the interferometric (u,v)-plane is available. Experiments are in progress; test telescopes have been erected, and trials in connecting large Cherenkov telescopes have been carried out. This paper reviews this interferometric method in view of the new possibilities offered by arrays of air Cherenkov telescopes, and outlines observational programs that should become realistic already in the rather near future.

Recent proposals have been advanced to apply imaging air Cherenkov telescope arrays to stellar intensity interferometry (SII). Of particular interest is the possibility of model-independent image recovery afforded by the good (u, v)-plane coverage of these arrays, as well as recent developments in phase retrieval techniques. The capabilities of these instruments used as SII receivers have already been explored for simple stellar objects, and here the focus is on reconstructing stellar images with non-uniform radiance distributions. We find that hot stars (T > 6000 K) containing hot and/or cool localized regions (T \sim 500 K) as small as \sim 0.1 mas can be imaged at short wavelengths ({\lambda} = 400 nm).

With its unprecedented light-collecting area for night-sky observations, the Cherenkov Telescope Array (CTA) holds great potential for also optical stellar astronomy, in particular as a multi-element intensity interferometer for realizing imaging with sub-milliarcsecond angular resolution. Such an order-of-magnitude increase of the spatial resolution achieved in optical astronomy will reveal the surfaces of rotationally flattened stars with structures in their circumstellar disks and winds, or the gas flows between close binaries. Image reconstruction is feasible from the second-order coherence of light, measured as the temporal correlations of arrival times between photons recorded in different telescopes. This technique (once pioneered by Hanbury Brown and Twiss) connects telescopes only with electronic signals and is practically insensitive to atmospheric turbulence and to imperfections in telescope optics. Detector and telescope requirements are very similar to those for imaging air Cherenkov observatories, the main difference being the signal processing (calculating cross correlations between single camera pixels in pairs of telescopes). Observations of brighter stars are not limited by sky brightness, permitting efficient CTA use during also bright-Moon periods. While other concepts have been proposed to realize kilometer-scale optical interferometers of conventional amplitude (phase-) type, both in space and on the ground, their complexity places them much further into the future than CTA, which thus could become the first kilometer-scale optical imager in astronomy.

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.

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

We describe a Chandra X-ray ToO project designed to isolate the site of TeV flaring in the radio galaxy M87. To date, we have triggered the Chandra observations only once (2010 April) and by the time of the first of our 9 observations, the TeV flare had ended. However, we found that the X-ray intensity of the unresolved nucleus was at an elevated level for our first observation. Of the more than 60 Chandra observations we have made of the M87 jet covering 9 years, the nucleus was measured at a comparably high level only 3 times. Two of these occasions can be associated with TeV flaring, and at the time of the third event, there were no TeV monitoring activities. From the rapidity of the intensity drop of the nucleus, we infer that the size of the emitting region is of order a few light days times the unknown beaming factor; comparable to the same sort of estimate for the TeV emitting region. We also find evidence of spectral evolution in the X-ray band which seems consistent with radiative losses affecting the non-thermal population of the emitting electrons within the unresolved nucleus.

Optical stellar intensity interferometry with air Cherenkov telescope arrays, composed of nearly 100 telescopes, will provide means to measure fundamental stellar parameters and also open the possibility of model-independent imaging. In addition to sensitivity issues, a main limitation of image recovery in intensity interferometry is the loss of phase of the complex degree of coherence during the measurement process. Nevertheless, several model-independent phase reconstruction techniques have been developed. Here we implement a Cauchy-Riemann based algorithm to recover images from simulated data. For bright stars (m_v~6) and exposure times of a few hours, we find that scale features such as diameters, oblateness and overall shapes are reconstructed with uncertainties of a few percent. More complex images are also well reconstructed with high degrees of correlation with the pristine image. Results are further improved by using a forward algorithm.

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.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

In the field of gamma-ray astronomy, irregular and noisy datasets make difficult the characterization of light-curve features in terms of statistical significance while properly accounting for trial factors associated with the search for variability at different times and over different timescales. In order to address these difficulties, we propose a method based on the Haar wavelet decomposition of the data. It allows statistical characterization of possible variability, embedded in a white noise background, in terms of a confidence level. The method is applied to artificially generated data for characterization as well as to the the very high energy M87 light curve recorded with VERITAS in 2008 which serves here as a realistic application example.

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

We report the discovery of TeV gamma-ray emission from the Type Ia supernova remnant (SNR) G120.1+1.4, known as Tycho's supernova remnant. Observations performed in the period 2008-2010 with the VERITAS ground-based gamma-ray observatory reveal weak emission coming from the direction of the remnant, compatible with a point source located at $00^{\rm h} \ 25^{\rm m} \ 27.0^{\rm s},\ +64^{\circ} \ 10^{\prime} \ 50^{\prime\prime}$ (J2000). The TeV photon spectrum measured by VERITAS can be described with a power-law $dN/dE = C(E/3.42\;\textrm{TeV})^{-\Gamma}$ with $\Gamma = 1.95 \pm 0.51_{stat} \pm 0.30_{sys}$ and $C = (1.55 \pm 0.43_{stat} \pm 0.47_{sys}) \times 10^{-14}$ cm$^{-2}$s$^{-1}$TeV$^{-1}$. The integral flux above 1 TeV corresponds to $\sim 0.9%$ percent of the steady Crab Nebula emission above the same energy, making it one of the weakest sources yet detected in TeV gamma rays. We present both leptonic and hadronic models which can describe the data. The lowest magnetic field allowed in these models is $\sim 80 \mu$G, which may be interpreted as evidence for magnetic field amplification.

The X-ray binary LSI +61 303, consisting of a main sequence Be star and a compact object has been detected in the TeV range with MAGIC and VERITAS, and showed a clear intensity modulation as a function of the orbital phase. We describe a gamma-ray attenuation model and apply it to this system. Our first result is that interaction of high energy photons with the background radiation produced by the main sequence star alone does not account for the observed modulation. We then include interactions between very high energy radiation and matter, and are able to constrain fundamental parameters of the system such as the mass of the compact object and the density of circumstellar matter around the Be star. In our analysis of the TeV data, we find that the compact object has mass $M_2>2.5M_{\odot}$ at the 99% confidence level, implying it is most likely a black hole. However, we find a column density which conflicts with results from X-ray observations, suggesting that attenuation may not play an important role in the modulation.

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.

Kilometric-scale optical imagers seem feasible to realize by intensity interferometry, using telescopes primarily erected for measuring Cherenkov light induced by gamma rays. Planned arrays envision 50--100 telescopes, distributed over some 1--4 km$^2$. Although array layouts and telescope sizes will primarily be chosen for gamma-ray observations, also their interferometric performance may be optimized. Observations of stellar objects were numerically simulated for different array geometries, yielding signal-to-noise ratios for different Fourier components of the source images in the interferometric $(u,v)$-plane. Simulations were made for layouts actually proposed for future Cherenkov telescope arrays, and for subsets with only a fraction of the telescopes. All large arrays provide dense sampling of the $(u,v)$-plane due to the sheer number of telescopes, irrespective of their geographic orientation or stellar coordinates. However, for improved coverage of the $(u,v)$-plane and a wider variety of baselines (enabling better image reconstruction), an exact east-west grid should be avoided for the numerous smaller telescopes, and repetitive geometric patterns avoided for the few large ones. Sparse arrays become severely limited by a lack of short baselines, and to cover astrophysically relevant dimensions between 0.1--3 milliarcseconds in visible wavelengths, baselines between pairs of telescopes should cover the whole interval 30--2000 m.

Intensity interferometry permits very long optical baselines and the observation of sub-milliarcsecond structures. Using planned kilometric arrays of air Cherenkov telescopes at short wavelengths, intensity interferometry may increase the spatial resolution achieved in optical astronomy by an order of magnitude, inviting detailed studies of the shapes of rapidly rotating hot stars with structures in their circumstellar disks and winds, or mapping out patterns of nonradial pulsations across stellar surfaces. Signal-to-noise in intensity interferometry favors high-temperature sources and emission-line structures, and is independent of the optical passband, be it a single spectral line or the broad spectral continuum. Prime candidate sources have been identified among classes of bright and hot stars. Observations are simulated for telescope configurations envisioned for large Cherenkov facilities, synthesizing numerous optical baselines in software, confirming that resolutions of tens of microarcseconds are feasible for numerous astrophysical targets.

Experiments are in progress to prepare for intensity interferometry with arrays of air Cherenkov telescopes. At the Bonneville Seabase site, near Salt Lake City, a testbed observatory has been set up with two 3-m air Cherenkov telescopes on a 23-m baseline. Cameras are being constructed, with control electronics for either off- or online analysis of the data. At the Lund Observatory (Sweden), in Technion (Israel) and at the University of Utah (USA), laboratory intensity interferometers simulating stellar observations have been set up and experiments are in progress, using various analog and digital correlators, reaching 1.4 ns time resolution, to analyze signals from pairs of laboratory telescopes.

Sub milli-arcsecond imaging in the visible band will provide a new perspective in stellar astrophysics. Even though stellar intensity interferometry was abandoned more than 40 years ago, it is capable of imaging and thus accomplishing more than the measurement of stellar diameters as was previously thought. Various phase retrieval techniques can be used to reconstruct actual images provided a sufficient coverage of the interferometric plane is available. Planned large arrays of Air Cherenkov telescopes will provide thousands of simultaneously available baselines ranging from a few tens of meters to over a kilometer, thus making imaging possible with unprecedented angular resolution. Here we investigate the imaging capabilities of arrays such as CTA or AGIS used as Stellar Intensity Interferometry receivers. The study makes use of simulated data as could realistically be obtained from these arrays. A Cauchy-Riemann based phase recovery allows the reconstruction of images which can be compared to the pristine image for which the data were simulated. This is first done for uniform disk stars with different radii and corresponding to various exposure times, and we find that the uncertainty in reconstructing radii is a few percent after a few hours of exposure time. Finally, more complex images are considered, showing that imaging at the sub-milli-arc-second scale is possible.

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.

We present results from an intensive VERITAS monitoring campaign of the high-frequency peaked BL Lac object 1ES 1218+304 in 2008/2009. Although 1ES 1218+304 was detected previously by MAGIC and VERITAS at a persistent level of ~6% of the Crab Nebula flux, the new VERITAS data reveal a prominent flare reaching ~20% of the Crab. While very high energy (VHE) flares are quite common in many nearby blazars, the case of 1ES 1218+304 (redshift z = 0.182) is particularly interesting since it belongs to a group of blazars that exhibit unusually hard VHE spectra considering their redshifts. When correcting the measured spectra for absorption by the extragalactic background light, 1ES 1218+304 and a number of other blazars are found to have differential photon indices less than 1.5. The difficulty in modeling these hard spectral energy distributions in blazar jets has led to a range of theoretical gamma-ray emission scenarios, one of which is strongly constrained by these new VERITAS observations. We consider the implications of the observed light curve of 1ES 1218+304, which shows day scale flux variations, for shock acceleration scenarios in relativistic jets, and in particular for the viability of kiloparsec-scale jet emission scenarios.

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.

At the time large air Cherenkov arrays are being discussed for future gamma-ray observatories, we review the relationship between the targeted capabilities and the main design parameters taking into account construction costs. As an example application, we describe a telescope array optimized for observations between 1 TeV and a few 100 TeV and use detailed simulations to estimate its performances in comparison to science objectives.

The imminent availability of large arrays of large light collectors deployed to exploit atmospheric Cherenkov radiation for gamma-ray astronomy at more than 100GeV, motivates the growing interest in application of intensity interferometry in astronomy. Indeed, planned arrays numbering up to one hundred telescopes will offer close to 5,000 baselines, ranging from less than 50m to more than 1000m. Recent and continuing signal processing technology developments reinforce this interest. Revisiting Stellar Intensity Interferometry for imaging is well motivated scientifically. It will fill the short wavelength (B/V bands) and high angular resolution (< 0.1mas) gap left open by amplitude interferometers. It would also constitute a first and important step toward exploiting quantum optics for astronomical observations, thus leading the way for future observatories. In this paper we outline science cases, technical approaches and schedule for an intensity interferometer to be constructed and operated in the visible using gamma-ray astronomy Air Cherenkov Telescopes as receivers.

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.

The Short GAmma Ray Front Air Cherenkov Experiment (SGARFACE) uses the Whipple 10 m telescope to search for bursts of $\gamma$ rays. SGARFACE is sensitive to bursts with duration from a few ns to $\sim$20 $\mu$s and with $\gamma$-ray energy above 100 MeV. SGARFACE began operating in March 2003 and has collected 2.2 million events during an exposure time of 2267 hours. A search for bursts of $\gamma$ rays from explosions of primordial black holes (PBH) was carried out. A Hagedorn-type PBH explosion is predicted to be visible within 60 pc of Earth. Background events were caused by cosmic rays and by atmospheric phenomena and their rejection was accomplished to a large extent using the time-resolved images. No unambiguous detection of bursts of $\gamma$ rays could be made as the remaining background events mimic the expected shape and time development of bursts. Upper limits on the PBH explosion rate were derived from the SGARFACE data and are compared to previous and future experiments. We note that a future array of large wide-field air-Cherenkov telescopes equipped with a SGARFACE-like trigger would be able to operate background-free with a 20 to 30 times higher sensitivity for PBH explosions.

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 dark matter science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses the prospects for dark matter detection with future gamma-ray experiments, and the complementarity of gamma-ray measurements with other indirect, direct or accelerator-based searches. We conclude that any comprehensive search for dark matter should include gamma-ray observations, both to identify the dark matter particle (through the charac- teristics of the gamma-ray spectrum) and to measure the distribution of dark matter in galactic halos.

This is a report on the findings of the extragalactic science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses extragalactic science topics including active galactic nuclei, cosmic ray acceleration in galaxies, galaxy clusters and large scale structure formation shocks, and the study of the extragalactic infrared and optical background radiation. The scientific potential of ground based gamma-ray observations of Gamma-Ray Bursts and dark matter annihilation radiation is covered in other sections of the white paper.

In this poster contribution we highlight the equivalence between an Imaging Air Cherenkov Telescope (IACT) array and an Intensity Interferometer for a range of technical requirements. We touch on the differences between a Michelson and an Intensity Interferometer and give a brief overview of the current IACT arrays, their upgrades and next generation concepts (CTA, AGIS, completion 2015). The latter are foreseen to include 30-90 telescopes that will provide 400-4000 different baselines that range in length between 50m and a kilometre. Intensity interferometry with such arrays of telescopes attains 50 micro-arcseconds resolution for a limiting V magnitude of ~8.5. 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 discuss possibilities for using IACT arrays as optical Intensity Interferometers.

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 SNR/cosmic-ray working group for the white paper on the status and future of ground-based 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 the status of past and current attempts to observe shell-type supernova remnants and diffuse emission from cosmic rays at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to study the acceleration of relativistic charged particles which is one of the main unsolved, yet fundamental, problems in modern astrophysics. The acceleration of particles relies on interactions between energetic particles and magnetic turbulence. In the case of SNRs we can perform spatially resolved studies in systems with known geometry, and the plasma physics deduced from these observations will help us to understand other systems where rapid particle acceleration is believed to occur and where observations as detailed as those of SNRs are not possible.

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.

(Abridged) We present a detailed analysis of week-long simultaneous observations of the blazar Mrk421 at 2-60 keV X-rays (RXTE) and TeV gamma-rays (Whipple and HEGRA) in 2001. The unprecedented quality of this dataset enables us to establish firmly the existence of the correlation between the TeV and X-ray luminosities, and to start unveiling some of its more detailed characteristics, in particular its energy dependence, and time variability. The source shows strong, highly correlated variations in X-ray and gamma-ray. No evidence of X-ray/gamma-ray interband lag is found on the full week dataset (<3 ks). However, a detailed analysis of the March 19 flare reveals that data are not consistent with the peak of the outburst in the 2-4 keV X-ray and TeV band being simultaneous. We estimate a 2.1+/-0.7 ks TeV lag. The amplitudes of the X-ray and gamma-ray variations are also highly correlated, and the TeV luminosity increases more than linearly w.r.t. the X-ray one. The strong correlation supports the standard model in which a unique electrons population produces the X-rays by synchrotron radiation and the gamma-ray component by inverse Compton scattering. However, for the individual best observed flares the gamma-ray flux scales approximately quadratically w.r.t. the X-ray flux, posing a serious challenge to emission models for TeV blazars. Rather special conditions and/or fine tuning of the temporal evolution of the physical parameters of the emission region are required in order to reproduce the quadratic correlation.

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.

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

Gamma-ray burst (GRB) observations at very high energies (VHE, E > 100 GeV) can impose tight constraints on some GRB emission models. Many GRB afterglow models predict a VHE component similar to that seen in blazars and plerions, in which the GRB spectral energy distribution has a double-peaked shape extending into the VHE regime. VHE emission coincident with delayed X-ray flare emission has also been predicted. GRB follow-up observations have had high priority in the observing program at the Whipple 10m Gamma-ray Telescope and GRBs will continue to be high priority targets as the next generation observatory, VERITAS, comes on-line. Upper limits on the VHE emission, at late times (>~4 hours), from seven GRBs observed with the Whipple Telescope are reported here.

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.

The high energy end of gamma-ray source spectra might provide important clues regarding the nature of the processes involved in gamma-ray emission. Several galactic sources with hard emission spectra extending up to more than 30TeV have already been reported. Measurements around 100TeV and above should be an important goal for the next generation of high energy gamma-ray astronomy experiments. Here we present several techniques providing the required exposure (100 km^2.h). We focus our study on three Imaging Atmospheric Cherenkov Technique (IACT) based approaches: low elevation observations, large field of view telescopes, and large telescope arrays. We comment on the advantages and disadvantages of each approach and report simulation based estimates of their energy ranges and sensitivities.

The trajectory of a primary gamma-ray detected with an array of at least four atmospheric Cherenkov imaging telescopes can be reconstructed from the shower image centroid positions and geometrical considerations independent of the primary energy. Using only the image centroid positions some cosmic-ray discrimination is also possible. This minimal approach opens the possibility of pushing the analysis threshold to lower values, close to the hardware threshold.

The first atmospheric Cherenkov telescope of VERITAS (the Very Energetic Radiation Imaging Telescope Array System) has been in operation since February 2005. We present here a technical description of the instrument and a summary of its performance. The calibration methods are described, along with the results of Monte Carlo simulations of the telescope and comparisons between real and simulated data. The analysis of TeV $\gamma$-ray observations of the Crab Nebula, including the reconstructed energy spectrum, is shown to give results consistent with earlier measurements. The telescope is operating as expected and has met or exceeded all design specifications.

Galaxy clusters might be sources of TeV gamma rays emitted by high-energy protons and electrons accelerated by large scale structure formation shocks, galactic winds, or active galactic nuclei. Furthermore, gamma rays may be produced in dark matter particle annihilation processes at the cluster cores. We report on observations of the galaxy clusters Perseus and Abell 2029 using the 10 m Whipple Cherenkov telescope during the 2003-2004 and 2004-2005 observing seasons. We apply a two-dimensional analysis technique to scrutinize the clusters for TeV emission. In this paper we first determine flux upper limits on TeV gamma-ray emission from point sources within the clusters. Second, we derive upper limits on the extended cluster emission. We subsequently compare the flux upper limits with EGRET upper limits at 100 MeV and theoretical models. Assuming that the gamma-ray surface brightness profile mimics that of the thermal X-ray emission and that the spectrum of cluster cosmic rays extends all the way from thermal energies to multi-TeV energies with a differential spectral index of -2.1, our results imply that the cosmic ray proton energy density is less than 7.9% of the thermal energy density for the Perseus cluster.

A calibrated laser pulse propagating through the atmosphere produces a flash of Rayleigh scattered light with an intensity that can be calculated very accurately when atmospheric conditions are good. This is used in a technique developed for the absolute calibration of ultra high energy cosmic ray fluorescence telescopes, and it can also be applied to imaging atmospheric Cherenkov telescopes (IACTs). In this paper we present the absolute calibration system being constructed and tested for the VERITAS project.

The idea of searching for optical signals from extraterrestrial civilizations has become increasingly popular over the last five years, with dedicated projects at a number of observatories. The method relies on the detection of a brief (few ns), intense light pulse with fast photon detectors. Ground-based gamma-ray telescopes such as the Whipple 10m, providing a large mirror area and equipped with an array of photomultiplier tubes (PMTs), are ideal instruments for this kind of observation if the background of cosmic-ray events can be rejected. We report here on a method for searching for optical SETI pulses, using background discrimination techniques based on the image shape.

The Narrabri intensity interferometer was successfully used until 1974 to observe 32 stars, all brighter than B=+2.5, among which some were found to have an angular diameter as small as 0.41+/-0.03 milli-arc-seconds (mas). The technique was then abandoned in favor of Michelson interferometry. Here we consider the technical feasibility and scientific potential of implementing intensity interferometry on Imaging Air Cherenkov Telescope arrays. The scientific motivations are varied, including stellar diameter measurements and investigations of the circumstellar environment. Long baselines and short wavelengths are easily accesible to this technique, making it uniquely suited for some applications.

The Short GAmma Ray Front Air Cherenkov Experiment (SGARFACE) is operated at the Whipple Observatory utilizing the Whipple 10m gamma-ray telescope. SGARFACE is sensitive to gamma-ray bursts of more than 100MeV with durations from 100ns to 35us and provides a fluence sensitivity as low as 0.8 gamma-rays per m^2 above 200MeV (0.05 gamma-rays per m^2 above 2GeV) and allows to record the burst time structure.

The Whipple 10-m telescope was used to observe M87 since 2000. No significant gamma-ray signal was found and upper limits compared to the HEGRA detection suggest the source may be variable. We found weak evidence for a correlation with the X-ray activity in 2000-2001 but this tendency did not persist in 2002-2003.

The Whipple 10m Very High Energy gamma-ray telescope has been used to search for indications of WIMP annihilation in the direction of the globular cluster M15. The upper limits derived constrain the amount of super-symmetric dark matter that may reside in globular clusters.

A decade after the discovery of TeV gamma-rays from the blazar Mrk 421 (Punch et al. 1992), the list of TeV blazars has increased to five BL Lac objects: Mrk 421 (Punch et al. 1992; Petry et al. 1996; Piron et al. 2001), Mrk 501 (Quinn et al. 1996; Aharonian et al. 1999; Djannati-Atai et al. 1999), 1ES2344+514 (Catanese et al. 1998), H1426+428 (Horan et al. 2000, 2002; Aharonian et al. 2002; Djannati-Atai et al. 2002) and 1ES1959+650 (Nishiyama et al. 1999; Konopelko et al. 2002; Holder et al. 2002). In this paper we report results from recent observations of Mrk 421, H1426+428 and 1ES1959+650 using the Whipple Observatory 10 m telescope.

We present the status and motivation of the Short GAmma Ray Front Air Cherenkov Experiment (SGARFACE) which will be operated parallel to standard Very High Energy gamma-ray observations using the Whipple 10m telescope. SGARFACE is sensitive to 100MeV-10GeV gamma-ray bursts with durations ranging from 100ns to 100us providing a fluence sensitivity as low as few $\rm 1E-9 erg/cm^2. Preliminary data taking started in November 2002.

The atmosphere is an intrinsic part of any ground based Cherenkov gamma-ray telescope, and the telescope response is therefore sensitive to unpredictable changes in the atmospheric transparency which are difficult to measure and interpret in the absence of a calibrated beam of high energy gamma-rays. In this paper, we use the detector response to Cherenkov emission from cosmic ray initiated air showers to obtain a relative calibration for data obtained under different instrumental and atmospheric conditions as well as over a range of source angles to the Zenith. We show that such a relative calibration is useful and efficient for data selection, for correcting the measured gamma-ray rate and for inter-calibration between the elements of an array of Cherenkov telescopes.

The BL Lac object H1426+428 was recently detected as a high energy gamma-ray source by the VERITAS collaboration (Horan et al. 2002). We have reanalyzed the 2001 portion of the data used in the detection in order to examine the spectrum of H1426+428 above 250 GeV. We find that the time-averaged spectrum agrees with a power law of the shape dF/dE = 10^(-7.31 +- 0.15(stat) +- 0.16(syst)) x E^(-3.50 +- 0.35(stat) +- 0.05(syst)) m^(-2)s^(-1)TeV^(-1) The statistical evidence from our data for emission above 2.5 TeV is 2.6 sigma. With 95% c.l., the integral flux of H1426+428 above 2.5 TeV is larger than 3% of the corresponding flux from the Crab Nebula. The spectrum is consistent with the (non-contemporaneous) measurement by Aharonian et al. (2002) both in shape and in normalization. Below 800 GeV, the data clearly favours a spectrum steeper than that of any other TeV Blazar observed so far indicating a difference in the processes involved either at the source or in the intervening space.

The detection of spectral variability of the gamma-ray blazar Mrk 421 at TeV energies is reported. Observations with the Whipple Observatory 10m gamma-ray telescope taken in 2000/2001 revealed exceptionally strong and long-lasting flaring activity. Flaring levels of 0.4 to 13 times that of the Crab Nebula flux provided sufficient statistics for a detailed study of the energy spectrum between 380 GeV and 8.2 TeV as a function of flux level. These spectra are well described by a power law with an exponential cutoff. There is no evidence for variation in the cutoff energy with flux, and all spectra are consistent with an average value for the cutoff energy of 4.3 TeV. The spectral index varies between 1.89 in a high flux state and 2.72 in a low state. The correlation between spectral index and flux is tight when averaging over the total 2000/2001 data set. Spectral measurements of Mrk~421 from previous years (1995/96 and 1999) by the Whipple collaboration are consistent with this flux-spectral index correlation, which suggest that this may be a constant or a long-term property of the source. If a similar flux-spectral index correlation were found for other gamma-ray blazars, this universal property could help disentangle the intrinsic emission mechanism from external absorption effects.

Exceptionally strong and long lasting flaring activity of the blazar Markarian 421 (Mrk 421) occurred between January and March 2001. Based on the excellent signal-to-noise ratio of the data we derive the energy spectrum between 260 GeV - 17 TeV with unprecedented statistical precision. The spectrum is not well described by a simple power law even with a curvature term. Instead the data can be described by a power law with exponential cutoff: $\rm {{dN}\over{dE}} \propto \rm E^{-2.14 \pm 0.03_{stat}} \times e^{-E/E_{0}} m^{-2} s^{-1} TeV^{-1}$ with $\rm E_{0} = 4.3 \pm 0.3_{stat} TeV$. Mrk 421 is the second $\gamma$-ray blazar that unambiguously exhibits an absorption-like feature in its spectral energy distribution at 3-6 TeV suggesting that this may be a universal phenomenon, possibly due to the extragalactic infra-red background radiation.

A new imaging atmospheric Cherenkov telescope (CANGAROO-II) with a light-weight reflector has been constructed. Light, robust, and durable mirror facets of containing CFRP (Carbon Fiber Reinforced Plastic) laminates were developed for the telescope. The attitude of each facet can be adjusted by stepping motors. In this paper, we describe the design, manufacturing, alignment procedure, and the performance of the CANGAROO-II optical reflector system.

In 1998 and 1999 the Whipple Observatory 10 m telescope was used to search for diffuse gamma ray emission from the Galactic Plane. No signifiant evidence of emission was found. Assuming the TeV emission profile matches EGRET observations above 1 GeV with a differential spectral index of 2.4, we derive an upper limit of {$\rm {3.0\cdot10^{-8}cm^{-2}s^{-1}sr^{-1}}$} for the average diffuse emission above {{$\rm500 GeV$}} in the galactic latitude range from {{$\rm-2^o$ to $\rm +2^o$}} at galactic longitude {{$\rm 40^o$}}. Comparisons with EGRET observations provide a lower limit of 2.31 for the differential spectral index of the diffuse emission, assuming there is no break in the spectrum between 30 GeV and 500 GeV. This constrains models for diffuse emission with a significant inverse Compton contribution.

We give an overview of the current status and scientific goals of VERITAS, a proposed hexagonal array of seven 10 m aperture imaging Cherenkov telescopes. The selected site is Montosa Canyon (1390 m a.s.l.) at the Whipple Observatory, Arizona. Each telescope, of 12 m focal length, will initially be equipped with a 499 element photomultiplier camera covering a 3.5 degree field of view. A central station will initiate the readout of 500 MHz FADCs upon receipt of multiple telescope triggers. The minimum detectable flux sensitivity will be 0.5% of the Crab Nebula flux at 200 GeV. Detailed simulations of the array's performance are presented elsewhere at this meeting. VERITAS will operate primarily as a gamma-ray observatory in the 50 GeV to 50 TeV range for the study of active galaxies, supernova remnants, pulsars and gamma-ray bursts.

We present the results of a search for pulsed TeV emission from the Crab pulsar using the Whipple Observatory's 10m gamma-ray telescope. The direction of the Crab pulsar was observed for a total of 73.4 hours between 1994 November and 1997 March. Spectral analysis techniques were applied to search for the presence of a gamma-ray signal from the Crab pulsar over the energy band 250GeV to 4TeV. At these energies we do not see any evidence of the 33ms pulsations present at lower energies from the Crab pulsar. The 99.9% confidence level upper limit for pulsed emission above 250GeV is derived to be 4.8 10^-12 cm^-2 s^-1 or <3% of the steady flux from the Crab Nebula. These results imply a sharp cut-off of the power-law spectrum seen by the EGRET instrument on the Compton Gamma-Ray Observatory. If the cut-off is exponential, it must begin at 60GeV or lower to accommodate these upper limits.

Markarian 501, a nearby (z=0.033) X-ray selected BL Lacertae object, is a well established source of Very High Energy (VHE, E>=300 GeV) gamma rays. Dramatic variability in its gamma-ray emission on time-scales from years to as short as two hours has been detected. Multiwavelength observations have also revealed evidence that the VHE gamma-ray and hard X-ray fluxes may be correlated. Here we present results of observations made with the Whipple Collaboration's 10 m Atmospheric Cerenkov Imaging Telescope during 1999 and discuss them in the context of observations made on Markarian 501 during the period from 1996-1998.