The VERITAS array of Cherenkov telescopes has carried out a deep observational program on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of good quality selected data, taken between January 2010 and May 2011. No significant $\gamma$-ray emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal galaxy currently known. We derive stringent bounds on various annihilating and decaying dark matter particle models. The upper limits on the velocity-weighted annihilation cross-section are $\mathrm{<\sigma v >^{95% CL} \lesssim 10^{-23} cm^{3} s^{-1}}$, improving our limits from previous observations of dwarf spheroidal galaxies by at least a factor of two for dark matter particle masses $\mathrm{m_{\chi}\gtrsim 300 GeV}$. The lower limits on the decay lifetime are at the level of $\mathrm{\tau^{95% CL} \gtrsim 10^{24} s}$. Finally, we address the interpretation of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario.

The variability of the blazar Markarian 421 in TeV gamma rays over a 14-year time period has been explored with the Whipple 10 m telescope. It is shown that the dynamic range of its flux variations is large and similar to that in X-rays. A correlation between the X-ray and TeV energy bands is observed during some bright flares and when the complete data sets are binned on long timescales. The main database consists of 878.4 hours of observation with the Whipple telescope, spread over 783 nights. The peak energy response of the telescope was 400 GeV with 20% uncertainty. This is the largest database of any TeV-emitting active galactic nucleus (AGN) and hence was used to explore the variability profile of Markarian 421. The time-averaged flux from Markarian 421 over this period was 0.446$\pm$0.008 Crab flux units. The flux exceeded 10 Crab flux units on three separate occasions. For the 2000-2001 season the average flux reached 1.86 Crab units, while in the 1996-1997 season the average flux was only 0.23 Crab units.

We report the discovery of an unidentified, extended source of very-high-energy (VHE) gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hours of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0.23^{\circ} \pm 0.03^{\circ} (stat)+0.04^{\circ}_{-0.02}^{\circ}(sys) and its spectrum is well-characterized by a differential power law (dN/dE = N_0 \times (E/TeV)^{-\Gamma}) with a photon index of {\Gamma} = 2.37 \pm 0.14 (stat) \pm 0.20 (sys) and a flux normalization of N0 = 1.5 \pm 0.2 (stat) \pm 0.4(sys) \times 10^-12 ph TeV^{-1} cm^{-2} s^{-1}. This yields an integral flux of 5.2 \pm 0.8 (stat) \pm 1.4 (sys) \times 10^-12 ph cm^{-2} s^{-1} above 320 GeV, corresponding to 3.7% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the supernova remnant shock.

We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles, due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 \times 10^(-12) erg cm^(-2) s^(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.

We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high-energy (VHE; E > 100 GeV) gamma-ray emission from the high-frequency-peaked BL Lac object RBS 0413. VERITAS, a ground-based gamma-ray observatory, detected VHE gamma rays from RBS 0413 with a statistical significance of 5.5 standard deviations (sigma) and a gamma-ray flux of (1.5 \pm 0.6stat \pm 0.7syst) \times 10^(-8) photons m^(-2) s^(-1) (\sim 1% of the Crab Nebula flux) above 250 GeV. The observed spectrum can be described by a power law with a photon index of 3.18 \pm 0.68stat \pm 0.30syst. Contemporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope detected HE gamma rays from RBS 0413 with a statistical significance of more than 9 sigma, a power-law photon index of 1.57 \pm 0.12stat +0.11sys -0.12sys and a gamma-ray flux between 300 MeV and 300 GeV of (1.64 \pm 0.43stat +0.31sys -0.22sys) \times 10^(-5) photons m^(-2) s^(-1). We present the results from Fermi-LAT and VERITAS, including a spectral energy distribution modeling of the gamma-ray, quasi-simultaneous X-ray (Swift-XRT), ultraviolet (Swift-UVOT) and R-band optical (MDM) data. We find that, if conditions close to equipartition are required, both the combined synchrotron self-Compton/external-Compton and the lepto-hadronic models are preferred over a pure synchrotron self-Compton model.

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 present the results of sixteen Swift-triggered GRB follow-up observations taken with the VERITAS telescope array from January, 2007 to June, 2009. The median energy threshold and response time of these observations was 260 GeV and 320 s, respectively. Observations had an average duration of 90 minutes. Each burst is analyzed independently in two modes: over the whole duration of the observations and again over a shorter time scale determined by the maximum VERITAS sensitivity to a burst with a t^-1.5 time profile. This temporal model is characteristic of GRB afterglows with high-energy, long-lived emission that have been detected by the Large Area Telescope (LAT) on-board the Fermi satellite. No significant VHE gamma-ray emission was detected and upper limits above the VERITAS threshold energy are calculated. The VERITAS upper limits are corrected for gamma-ray extinction by the extragalactic background light (EBL) and interpreted in the context of the keV emission detected by Swift. For some bursts the VHE emission must have less power than the keV emission, placing constraints on inverse Compton models of VHE emission.

We report on the VERITAS discovery of very-high-energy (VHE) gamma- ray emission above 200 GeV from the high-frequency-peaked BL Lac object RXJ0648.7+1516 (GBJ0648+1516), associated with 1FGLJ0648.8+1516. The photon spectrum above 200 GeV is fit by a power law dN/dE = F0(E/E0)-{\Gamma} with a photon index {\Gamma} of 4.4 {\pm} 0.8stat {\pm}0.3syst and a flux normalization F0 of (2.3 {\pm}0.5stat {\pm} 1.2sys) {\times} 10-11 TeV-1cm-2s-1 with E0 = 300 GeV. No VHE vari- ability is detected during VERITAS observations of RXJ0648.7+1516 between 2010 March 4 and April 15. Following the VHE discovery, the optical identifica- tion and spectroscopic redshift were obtained using the Shane 3-m Telescope at the Lick Observatory, showing the unidentified object to be a BL Lac type with a redshift of z = 0.179. Broadband multiwavelength observations contemporaneous with the VERITAS exposure period can be used to sub-classify the blazar as a high-frequency-peaked BL Lac (HBL) object, including data from the MDM ob- servatory, Swift-UVOT and XRT, and continuous monitoring at photon energies above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the absence of undetected, high-energy rapid variability, the one-zone synchrotron self-Compton model (SSC) overproduces the high-energy gamma-ray emission measured by the Fermi-LAT over 2.3 years. The SED can be parameterized sat- isfactorily with an external-Compton or lepto-hadronic model, which have two and six additional free parameters, respectively, compared to the one-zone SSC model.

We report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 Gigaelectronvolts (GeV) with the VERITAS array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 Megaelectronvolts (MeV) and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.

Giant X-ray outbursts, with luminosities of about $ 10^{37}$ erg s$^{-1}$, are observed roughly every 5 years from the nearby Be/pulsar binary 1A 0535+262. In this article, we present observations of the source with VERITAS at very-high energies (VHE; E$>$100 GeV) triggered by the X-ray outburst in December 2009. The observations started shortly after the onset of the outburst, and they provided comprehensive coverage of the episode, as well as the 111-day binary orbit. No VHE emission is evident at any time. We also examined data from the contemporaneous observations of 1A 0535+262 with the Fermi/LAT at high energy photons (HE; E$>$0.1 GeV) and failed to detect the source at GeV energies. The X-ray continua measured with the Swift/XRT and the RXTE/PCA can be well described by the combination of blackbody and Comptonized emission from thermal electrons. Therefore, the gamma-ray and X-ray observations suggest the absence of a significant population of non-thermal particles in the system. This distinguishes 1A~0535+262 from those Be X-ray binaries (such as PSR B1259--63 and LS I +61$^{\circ}$303) that have been detected at GeV--TeV energies. We discuss the implications of the results on theoretical models.

We report on very-high-energy ($>$100 GeV) gamma-ray observations of Swift J164449.3+573451, an unusual transient object first detected by the {\it Swift} Observatory and later detected by multiple radio, optical and X-ray observatories. A total exposure of 28 hours was obtained on Swift J164449.3+573451 with VERITAS during 2011 March 28 -- April 15. We do not detect the source and place a differential upper limit on the emission at 500 GeV during these observations of $1.4 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ (99% confidence level). We also present time-resolved upper limits and use a flux limit averaged over the X-ray flaring period to constrain various emission scenarios that can accommodate both the radio-through-X-ray emission detected from the source and the lack of detection by VERITAS.

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

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

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

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.

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

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

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

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

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

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.

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

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.

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.

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.

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

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

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

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

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

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

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

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

An eleven channel digital asynchronous transceiver (DAT) employing parallel optical link technology has been developed for trigger signal distribution across the Very Energetic Radiation Imaging Telescope Array System (VERITAS). Combinatorial logic functions are implemented in Xilinx Spartan 3 FPGAs, providing a versatile solution adaptable for use in future atmospheric Cerenkov detectors and other high-energy astroparticle experiments. The device is dead-time free and introduces a minimal skew of 1.6 ns between channels. The jitter on each DAT channel is less than 0.8 ns 95% of the time, allowing communication between telescopes and a central trigger system separated by hundreds of meters, without limiting array performance.

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.

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.

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.

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.

Mrk 421 is the first TeV blazar found to exhibit significant spectral variability during strong flaring activity, showing hardening of the TeV spectrum in high emission states. Mrk 421 is also known to exhibit flux variability on time scales as short as 15 minutes. In this paper we present studies of hourly spectral variability of Mrk 421 in 2001 using data from the Whipple Observatory 10m gamma-ray telescope.

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.

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.

The hunt for cosmic TeV particle accelerators is prospering through Imaging Atmospheric Cerenkov Telescopes. We face challenges such as low light levels and MHz trigger rates, and the need to distinguish between particle air showers stemming from primary gamma rays and those due to the hadronic cosmic ray background. Our test beam is provided by the Crab Nebula, a steady accelerator of particles to energies beyond 20 TeV. Highly variable gamma-ray emission, coincident with flares at longer wavelengths, is revealing the particle acceleration mechanisms at work in the relativistic jets of Active Galaxies. These 200 GeV to 20 TeV photons propagating over cosmological distances allow us to place a limit on the infra-red background linked to galaxy formation and, some speculate, to the decay of massive relic neutrinos. Gamma rays produced in neutralino annihilation or the evaporation of primordial black holes may also be detectable. These phenomena and a zoo of astrophysical objects will be the targets of the next generation multi-national telescope facilities.

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 present the results of a search for pulsed TeV emission from the Crab pulsar using the Whipple Observatory's 10 m gamma-ray telescope. The direction of the Crab pulsar was observed for a total of 73.4 hours between 1994 November and 1997 March. During this period the Whipple 10 m telescope was operated at its lowest energy threshold to date. Spectral analysis techniques were applied to search for the presence of a gamma-ray signal from the Crab pulsar over the energy band 250 GeV to 4 TeV. We do not see any evidence of the 33 ms pulsations present in other energy bands from the Crab pulsar. The 99.9% confidence level upper limit for pulsed emission above 250 GeV is derived to be 4.8x10^-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 60 GeV or lower to accommodate these upper limits.

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.

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

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.

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.

The BL Lacertae object Markarian 501 was identified as a source of gamma-ray emission at the Whipple Observatory in March 1995. Here we present a flux variability analysis on several times-scales of the 233 hour data set accumulated over 213 nights (from March 1995 to July 1998) with the Whipple Observatory 10 m atmospheric Cherenkov imaging telescope. In 1995, with the exception of a single night, the flux from Markarian 501 was constant on daily and monthly time-scales and had an average flux of only 10% that of the Crab Nebula, making it the weakest VHE source detected to date. In 1996, the average flux was approximately twice the 1995 flux and showed significant month-to-month variability. No significant day-scale variations were detected. The average gamma-ray flux above ~350 GeV in the 1997 observing season rose to 1.4 times that of the Crab Nebula -- 14 times the 1995 discovery level -- allowing a search for variability on time-scales shorter than one day. Significant hour-scale variability was present in the 1997 data, with the shortest, observed on MJD 50607, having a doubling time of ~2 hours. In 1998 the average emission level decreased considerably from that of 1997 (to ~20% of the Crab Nebula flux) but two significant flaring events were observed. Thus, the emission from Markarian 501 shows large amplitude and rapid flux variability at very high energies as does Markarian 421. It also shows large mean flux level variations on year-to-year time-scales, behaviour which has not been seen from Markarian 421 so far.

The close relation between ground-based TeV observations and satellite borne $\gamma$-ray measurements has been important for the understanding of blazars. The observations which involve the TeV component in blazar studies are reviewed.

We present the discovery of >350 GeV gamma-ray emission from the BL Lacertae object 1ES 2344+514 with the Whipple Observatory 10m gamma-ray telescope. This is the third BL Lac object detected at gamma-ray energies above 300 Gev, the other two being Markarian 421 (Mrk 421) and Mrk501. These three active galactic nuclei are all X-ray selected and have the lowest known redshifts of any BL Lac objects currently identified. The evidence for emission derives primarily from an apparent flare on December 20, 1995 when a 6 sigma excess was detected with a flux approximately 63% of the very high energy gamma-ray emission from the Crab Nebula, the standard candle for TeV gamma-ray sources. Excluding the flare, observations between October 1995 and January 1996 yield a 4 sigma detection corresponding to 11% of the VHE Crab Nebula flux. Observations spanning September 1996 to January 1997 failed to yield a significant detection of a steady flux or any flaring. For this period, the 99.9% confidence level upper limit is <8% of the Crab Nebula. The low baseline emission level and variations in nightly and yearly flux of 1ES 22344+514 are the same as the VHE emission characteristics of Mrk 421 and Mrk 501

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

We present multiwavelength observations of the BL Lacertae object Markarian 501 (Mrk 501) in 1997 between April 8 and April 19. Evidence of correlated variability is seen in very high energy (VHE, E > 350 GeV) gamma-ray observations taken with the Whipple Observatory gamma-ray telescope, data from the Oriented Scintillation Spectrometer Experiment of the Compton Gamma-Ray Observatory, and quicklook results from the All-Sky Monitor of the Rossi X-ray Timing Explorer while the Energetic Gamma-Ray Experiment Telescope did not detect Mrk 501. Short term optical correlations are not conclusive but the U-band flux observed with the 1.2m telescope of the Whipple Observatory was 10% higher than in March. The average energy output of Mrk 501 appears to peak in the 2 keV to 100 keV range suggesting an extension of the synchrotron emission to at least 100 keV, the highest observed in a blazar and ~100 times higher than that seen in the other TeV-emitting BL Lac object, Mrk 421. The VHE gamma-ray flux observed during this period is the highest ever detected from this object. The VHE gamma-ray energy output is somewhat lower than the 2-100 keV range but the variability amplitude is larger. The correlations seen here do not require relativistic beaming of the emission unless the VHE spectrum extends to >5 TeV.

The ~70 unidentified sources of the EGRET sky survey may be one of its most important legacies. The identification of these sources at other wavelengths is critical to understanding their nature. Many have flat spectra out to 10 GeV which, if extrapolated to TeV energies, would be easily detectable relative to the steeply falling diffuse background. The Whipple Observatory gamma-ray telescope has been used to observe a number of these which were selected based on their position, intensity and spectrum and in some cases based on a possible association with a supernova remnant or pulsar. No significant emission has been detected from these sources, and upper limits are given.

Markarian 501 is only the second extragalactic source to be detected with high statistical certainty at TeV energies; it is similar in many ways to Markarian 421. The Whipple Observatory gamma-ray telescope has been used to observe the AGN Markarian 501 in 1996 and 1997, the years subsequent to its initial detection. The apparent variability on the one-day time-scale observed in TeV gamma rays in 1995 is confirmed and compared with the variability in Markarian 421. Observations at X-ray and optical wavelengths from 1997 are also presented.

If supernova remnants (SNRs) are the sites of cosmic-ray acceleration, the associated nuclear interactions should result in observable fluxes of TeV gamma-rays from the nearest SNRs. Measurements of the gamma-ray flux from six nearby, radio-bright, SNRs have been made with the Whipple Observatory gamma-ray telescope. No significant emission has been detected and upper limits on the $>$300 GeV flux are reported. Three of these SNRs (IC443, gamma-Cygni and W44) are spatially coincident with low latitude unidentified sources detected with EGRET. These upper limits weaken the case for the simplest models of shock acceleration and energy dependent propagation.

Markarian 421 was the first extragalactic source to be detected with high statistical certainty at TeV energies. The Whipple Observatory gamma-ray telescope has been used to observe the Active Galactic Nucleus, Markarian 421 in 1996 and 1997. The rapid variability observed in TeV gamma rays in previous years is confirmed. Doubling times as short as 15 minutes are reported with flux levels reaching 15 photons per minute. The TeV energy spectrum is derived using two independent methods. The implications for the intergalactic infra-red medium of an observed unbroken power law spectrum up to energies of 5 TeV is discussed.

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

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