The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259--63/LS 2883, making it an ideal candidate for the study of high-energy non-thermal emission. It has been the subject of multi-wavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high-energy gamma-ray observations carried out by VERITAS before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that very high-energy gamma rays are produced via an inverse-Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than $\sim$2 G before the disappearance of the radio pulsar and greater than $\sim$10 G afterwards.

The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85 GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes. The observations, taken between 2008 May and 2010 June and totalling 53 hours, resulted in the discovery of $\gamma$-ray emission from the blazar, which has a redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of 5.5 standard deviations above the background with an integral flux of (2.8$\pm0.7_{\mathrm{stat}}\pm0.8_{\mathrm{sys}}$) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ (1.2\% of the Crab Nebula's flux) above 200 GeV. The measured spectrum is described well by a power law from 0.2 TeV to 1.3 TeV with a photon index of 3.1 $\pm$ 0.4$_{\mathrm{stat}}$ $\pm$ 0.2$_{\mathrm{sys}}$. Quasi-simultaneous multi-wavelength data from the Fermi Large Area Telescope (0.3--300 GeV) and the Swift X-ray Telescope (0.2--10 keV) are additionally used to model the properties of the emission region. A synchrotron self-Compton model produces a good representation of the multi-wavelength data. Adding an external-Compton or a hadronic component also adequately describes the data.

The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $\gamma$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation with an increase of a factor $\sim$4 in 25 minutes, and a falling time of $\sim$50 minutes. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.

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.

During moonlit nights, observations with ground-based Cherenkov telescopes at very high energies (VHE, $E>100$ GeV) are constrained since the photomultiplier tubes (PMTs) in the telescope camera are extremely sensitive to the background moonlight. Observations with the VERITAS telescopes in the standard configuration are performed only with a moon illumination less than 35$\%$ of full moon. Since 2012, the VERITAS collaboration has implemented a new observing mode under bright moonlight, by either reducing the voltage applied to the PMTs (reduced-high-voltage configuration, RHV), or by utilizing UV-transparent filters. While these operating modes result in lower sensitivity and increased energy thresholds, the extension of the available observing time is useful for monitoring variable sources such as blazars and sources requiring spectral measurements at the highest energies. In this paper we report the detection of $\gamma$-ray flaring activity from the BL Lac object 1ES 1727+502 during RHV observations. This detection represents the first evidence of VHE variability from this blazar. The integral flux is $(1.1\pm0.2)\times10^{-11}\mathrm{cm^{-2}s^{-1}}$ above 250 GeV, which is about five times higher than the low-flux state. The detection triggered additional \veritas\ observations during standard dark-time. Multiwavelength observations with the FLWO 48" telescope, and the Swift and Fermi satellites are presented and used to produce the first spectral energy distribution (SED) of this object during $\gamma$-ray flaring activity. The SED is then fitted with a standard synchrotron-self-Compton model, placing constraints on the properties of the emitting region and of the acceleration mechanism at the origin of the relativistic particle population in the jet.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

We present the results of 71.6 hours of observations of the Geminga pulsar (PSR J0633+1746) with the VERITAS very-high-energy gamma-ray telescope array. Data taken with VERITAS between November 2007 and February 2013 were phase-folded using a Geminga pulsar timing solution derived from data recorded by the XMM-\emph{Newton} and \emph{Fermi}-LAT space telescopes. No significant pulsed emission above 100 GeV is observed, and we report upper limits at the 95% confidence level on the integral flux above 135 GeV (spectral analysis threshold) of 4.0$\times10^{-13}$ s$^{-1}$ cm$^{-2}$ and 1.7$\times10^{-13}$ s$^{-1}$ cm$^{-2}$ for the two principal peaks in the emission profile. These upper limits, placed in context with phase-resolved spectral energy distributions determined from five years of data from the \emph{Fermi}-LAT, constrain possible hardening of the Geminga pulsar emission spectra above $\sim$50 GeV.

We summarize broadband observations of the TeV-emitting blazar 1ES 1959+650, including optical R-band observations by the robotic telescopes Super-LOTIS and iTelescope, UV observations by Swift UVOT, X-ray observations by the Swift X-ray Telescope (XRT), high-energy gamma-ray observations with the Fermi Large Area Telescope (LAT) and very-high-energy (VHE) gamma-ray observations by VERITAS above 315 GeV, all taken between 17 April 2012 and 1 June 2012 (MJD 56034 and 56079). The contemporaneous variability of the broadband spectral energy distribution is explored in the context of a simple synchrotron self Compton (SSC) model. In the SSC emission scenario, we find that the parameters required to represent the high state are significantly different than those in the low state. Motivated by possible evidence of gas in the vicinity of the blazar, we also investigate a reflected-emission model to describe the observed variability pattern. This model assumes that the non-thermal emission from the jet is reflected by a nearby cloud of gas, allowing the reflected emission to re-enter the blob and produce an elevated gamma-ray state with no simultaneous elevated synchrotron flux. The model applied here, although not required to explain the observed variability pattern, represents one possible scenario which can describe the observations. As applied to an elevated VHE state of 66% of the Crab Nebula flux, observed on a single night during the observation period, the reflected-emission scenario does not support a purely leptonic non-thermal emission mechanism. The reflected emission model does, however, predict a reflected photon field with sufficient energy to enable elevated gamma-ray emission via pion production with protons of energies between 10 and 100 TeV.

We present results from VERITAS observations of the BL Lac object PG 1553+113 spanning the years 2010, 2011, and 2012. The time-averaged spectrum, measured between 160 and 560\,GeV, is well described by a power law with a spectral index of $4.33 \pm 0.09$. The time-averaged integral flux above $200\,$GeV measured for this period was $(1.69 \pm 0.06) \times 10^{-11} \, \mathrm{ph} \, \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$, corresponding to 6.9\% of the Crab Nebula flux. We also present the combined $\gamma$-ray spectrum from the Fermi Large Area Telescope and VERITAS covering an energy range from 100~MeV to 560~GeV. The data are well fit by a power law with an exponential cutoff at $\rm {101.9 \pm 3.2 \, \mathrm{GeV}} $. The origin of the cutoff could be intrinsic to PG~1553+113 or be due to the $\gamma$-ray opacity of our universe through pair production off the extragalactic background light (EBL). Given lower limits to the redshift of $\rm z \negthinspace > \negthinspace 0.395$ based on optical/UV observations of PG~1553+113, the cutoff would be dominated by EBL absorption. Conversely, the small statistical uncertainties of the VERITAS energy spectrum have allowed us to provide a robust upper limit on the redshift of PG 1553+113 of $z \negthinspace \leq \negthinspace 0.62$. A strongly-elevated mean flux of $(2.50 \pm 0.14) \times 10^{-11} \, \mathrm{ph} \, \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$ (10.3\% of the Crab Nebula flux) was observed during 2012, with the daily flux reaching as high as $(4.44 \pm 0.71) \times 10^{-11} \, \mathrm{ph} \, \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$ (18.3\% of the Crab Nebula flux) on MJD 56048. The light curve measured during the 2012 observing season is marginally inconsistent with a steady flux, giving a $\chi^2$ probability for a steady flux of 0.03\%.

We report on multifrequency observations performed during 2012 December-2013 August of the first narrow-line Seyfert 1 galaxy detected in gamma rays, PMN J0948+0022 ($z$ = 0.5846). A gamma-ray flare was observed by the Large Area Telescope on board Fermi during 2012 December-2013 January, reaching a daily peak flux in the 0.1-100 GeV energy range of (155 $\pm$ 31) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 2013 January 1, corresponding to an apparent isotropic luminosity of about 1.5$\times$10$^{48}$ erg s$^{-1}$. The gamma-ray flaring period triggered Swift and VERITAS observations in addition to radio and optical monitoring by OVRO, MOJAVE, and CRTS. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the gamma-ray flare, reaching a record flux for this source from optical to gamma rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of F (> 0.2 TeV) < 4.0$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$. We compared the spectral energy distribution (SED) of the flaring state in 2013 January with that of an intermediate state observed in 2011. The two SEDs, modelled as synchrotron emission and an external Compton scattering of seed photons from a dust torus, can be modelled by changing both the electron distribution parameters and the magnetic field.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

Prompt emission from the very fluent and nearby (z=0.34) gamma-ray burst GRB 130427A was detected by several orbiting telescopes and by ground-based, wide-field-of-view optical transient monitors. Apart from the intensity and proximity of this GRB, it is exceptional due to the extremely long-lived high-energy (100 MeV to 100 GeV) gamma-ray emission, which was detected by the Large Area Telescope on the Fermi Gamma-ray Space Telescope for ~70 ks after the initial burst. The persistent, hard-spectrum, high-energy emission suggests that the highest-energy gamma rays may have been produced via synchrotron self-Compton processes though there is also evidence that the high-energy emission may instead be an extension of the synchrotron spectrum. VERITAS, a ground-based imaging atmospheric Cherenkov telescope array, began follow-up observations of GRB 130427A ~71 ks (~20 hr) after the onset of the burst. The GRB was not detected with VERITAS; however, the high elevation of the observations, coupled with the low redshift of the GRB, make VERITAS a very sensitive probe of the emission from GRB 130427A for E > 100 GeV. The non-detection and consequent upper limit derived place constraints on the synchrotron self-Compton model of high-energy gamma-ray emission from this burst.

We report on deep observations of the extended TeV gamma-ray source MGRO J1908+06 made with the VERITAS very high energy (VHE) gamma-ray observatory. Previously, the TeV emission has been attributed to the pulsar wind nebula (PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a significance level of 14 standard deviations (14 sigma) and measure a photon index of 2.20 +/- 0.10_stat +/- 0.20_sys. The TeV emission is extended, covering the region near PSR J1907+0602 and also extending towards SNR G40.5--0.5. When fitted with a 2-dimensional Gaussian, the intrinsic extension has a standard deviation of sigma_src = 0.44 +/- 0.02 degrees. In contrast to other TeV PWNe of similar age in which the TeV spectrum softens with distance from the pulsar, the TeV spectrum measured near the pulsar location is consistent with that measured at a position near the rim of G40.5--0.5, 0.33 degrees away.

We present deep VERITAS observations of the blazar PKS 1424+240, along with contemporaneous Fermi Large Area Telescope, Swift X-ray Telescope and Swift UV Optical Telescope data between 2009 February 19 and 2013 June 8. This blazar resides at a redshift of $z\ge0.6035$, displaying a significantly attenuated gamma-ray flux above 100 GeV due to photon absorption via pair-production with the extragalactic background light. We present more than 100 hours of VERITAS observations from three years, a multiwavelength light curve and the contemporaneous spectral energy distributions. The source shows a higher flux of (2.1$\pm0.3$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV in 2009 and 2011 as compared to the flux measured in 2013, corresponding to (1.02$\pm0.08$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV. The measured differential very high energy (VHE; $E\ge100$ GeV) spectral indices are $\Gamma=$3.8$\pm$0.3, 4.3$\pm$0.6 and 4.5$\pm$0.2 in 2009, 2011 and 2013, respectively. No significant spectral change across the observation epochs is detected. We find no evidence for variability at gamma-ray opacities of greater than $\tau=2$, where it is postulated that any variability would be small and occur on longer than year timescales if hadronic cosmic-ray interactions with extragalactic photon fields provide a secondary VHE photon flux. The data cannot rule out such variability due to low statistics.

TeV J2032+4130 was the first unidentified source discovered at very high energies (VHE; E $>$ 100 GeV), with no obvious counterpart in any other wavelength. It is also the first extended source to be observed in VHE gamma rays. Following its discovery, intensive observational campaigns have been carried out in all wavelengths in order to understand the nature of the object, which have met with limited success. We report here on a deep observation of TeV J2032+4130, based on 48.2 hours of data taken from 2009 to 2012 by the VERITAS (Very Energetic Radiation Imaging Telescope Array System) experiment. The source is detected at 8.7 standard deviations ($\sigma$) and is found to be extended and asymmetric with a width of 9.5$^{\prime}$$\pm$1.2$^{\prime}$ along the major axis and 4.0$^{\prime}$$\pm$0.5$^{\prime}$ along the minor axis. The spectrum is well described by a differential power law with an index of 2.10 $\pm$ 0.14$_{stat}$ $\pm$ 0.21$_{sys}$ and a normalization of (9.5 $\pm$ 1.6$_{stat}$ $\pm$ 2.2$_{sys}$) $\times$ 10$^{-13}$TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. We interpret these results in the context of multiwavelength scenarios which particularly favor the pulsar wind nebula (PWN) interpretation.

In March 2013, a flaring episode from the Crab Nebula lasting ~2 weeks was detected by the Fermi-LAT (Large Area Telescope on board the Fermi Gamma-ray Space Telescope). VERITAS provides simultaneous observations throughout this period. During the flare, the Fermi-LAT detected a 20-fold increase in flux above the average synchrotron flux >100 MeV seen from the Crab Nebula. Simultaneous measurements with VERITAS are consistent with the non-variable long-term average Crab Nebula flux at TeV energies. Assuming a linear correlation between the very-high-energy flux change >1 TeV and the flux change seen in the Fermi-LAT band >100 MeV during the period of simultaneous observations, the linear correlation factor can be constrained to be at most 8.6 * 10^-3 with 95% confidence.

The high-frequency-peaked BL Lacertae object 1ES 0229+200 is a relatively distant (z = 0.1396), hard-spectrum (Gamma ~ 2.5), very-high-energy-emitting (E > 100 GeV) gamma-ray blazar. Very-high-energy measurements of this active galactic nucleus have been used to place constraints on the intensity of the extragalactic background light and the intergalactic magnetic field. A multi-wavelength study of this object centered around very-high-energy observations by VERITAS is presented. This study obtained, over a period of three years, an 11.7 standard deviation detection and an average integral flux F(E>300 GeV) = (23.3 +- 2.8_stat +- 5.8_sys) x 10^-9 photons m^-2 s^-1, or 1.7% of the Crab Nebula's flux (assuming the Crab Nebula spectrum measured by H.E.S.S). Supporting observations from Swift and RXTE are analyzed. The Swift observations are combined with previously published Fermi observations and the very-high-energy measurements to produce an overall spectral energy distribution which is then modeled assuming one-zone synchrotron-self-Compton emission. The chi^2 probability of the TeV flux being constant is 1.6%. This, when considered in combination with measured variability in the X-ray band, and the demonstrated variability of many TeV blazars, suggests that the use of blazars such as 1ES 0229+200 for intergalactic magnetic field studies may not be straightforward and challenges models that attribute hard TeV spectra to secondary gamma-ray production along the line of sight.

HESS J0632+057 is the only gamma-ray binary known so far whose position in the sky allows observations with ground-based observatories both in the northern and southern hemispheres. Here we report on long-term observations of HESS J0632+057 conducted with the VERITAS and H.E.S.S. Cherenkov Telescopes and the X-ray Satellite Swift, spanning a time range from 2004 to 2012 and covering most of the system's orbit. The VHE emission is found to be variable, and is correlated with that at X-ray energies. An orbital period of $315 ^{+6}_{-4}$ days is derived from the X-ray data set, which is compatible with previous results, $P = (321 \pm 5$) days. The VHE light curve shows a distinct maximum at orbital phases close to 0.3, or about 100 days after periastron passage, which coincides with the periodic enhancement of the X-ray emission. Furthermore, the analysis of the TeV data shows for the first time a statistically significant ($> 6.5 \sigma$) detection at orbital phases 0.6--0.9. The obtained gamma-ray and X-ray light curves and the correlation of the source emission at these two energy bands are discussed in the context of the recent ephemeris obtained for the system. Our results are compared to those reported for other gamma-ray binaries.

We report results from TeV gamma-ray observations of the microquasar Cygnus X-3. The observations were made with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) over a time period from 2007 June 11 to 2011 November 28. VERITAS is most sensitive to gamma rays at energies between 85 GeV to 30 TeV. The effective exposure time amounts to a total of about 44 hours, with the observations covering six distinct radio/X-ray states of the object. No significant TeV gamma-ray emission was detected in any of the states, nor with all observations combined. The lack of a positive signal, especially in the states where GeV gamma rays were detected, places constraints on TeV gamma-ray production in Cygnus X-3. We discuss the implications of the results.

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 present the results of a multiwavelength observational campaign on the TeV binary system LS I +61 303 with the VERITAS telescope array (>200 GeV), Fermi-LAT (0.3-300 GeV), and Swift-XRT (2-10 keV). The data were taken from December 2011 through January 2012 and show a strong detection in all three wavebands. During this period VERITAS obtained 24.9 hours of quality selected livetime data in which LS I +61 303 was detected at a statistical sig- nificance of 11.9 sigma. These TeV observations show evidence for nightly variability in the TeV regime at a post-trial significance of 3.6 sigma. The combination of the simultaneously obtained TeV and X-ray fluxes do not demonstrate any evidence for a correlation between emission in the two bands. For the first time since the launch of the Fermi satellite in 2008, this TeV detection allows the construction of a detailed MeV-TeV spectral energy distribution from LS I +61 303. This spectrum shows a distinct cutoff in emission near 4 GeV, with emission seen by the VERITAS observations following a simple power-law above 200 GeV. This feature in the spectrum of LS I +61 303, obtained from overlapping observations with Fermi-LAT and VERITAS, may indicate that there are two distinct populations of accelerated particles producing the GeV and TeV emission.

We report on VERITAS observations of the BL Lac object B2 1215+30 between 2008 and 2012. During this period, the source was detected at very high energies (VHE; E > 100 GeV) by VERITAS with a significance of $8.9\sigma$ and showed clear variability on time scales larger than months. In 2011, the source was found to be in a relatively bright state and a power-law fit to the differential photon spectrum yields a spectral index of $3.6 \pm 0.4_{\mathrm{stat}} \pm 0.3_{\mathrm{syst}}$ with an integral flux above 200 GeV of $(8.0 \pm 0.9_{\mathrm{stat}} \pm 3.2_{\mathrm{syst}}) \times 10^{-12}\, \mathrm{cm}^{-2} \mathrm{s}^{-1}$. No short term variability could be detected during the bright state in 2011. Multi-wavelength data were obtained contemporaneous with the VERITAS observations in 2011 and cover optical (Super-LOTIS, MDM, Swift-UVOT), X-ray (Swift-XRT), and gamma-ray (Fermi-LAT) frequencies. These were used to construct the spectral energy distribution (SED) of B2 1215+30. A one-zone leptonic model is used to model the blazar emission and the results are compared to those of MAGIC from early 2011 and other VERITAS-detected blazars. The SED can be well reproduced with model parameters typical for VHE-detected BL Lacs.

Compilation of papers contributed by the VERITAS Collaboration to the 33rd International Cosmic Ray Conference, held 2-9 July, 2013, in Rio de Janeiro, Brazil.

We report the detection of a new TeV gamma-ray source, VER J0521+211, based on observations made with the VERITAS imaging atmospheric Cherenkov telescope array. These observations were motivated by the discovery of a cluster of >30GeV photons in the first year of Fermi-LAT observations. VER J0521+211 is relatively bright at TeV energies, with a mean photon flux of 1.93 +/- 0.13_stat +/- 0.78_sys 10^-11 cm-2 s-1 above 0.2 TeV during the period of the VERITAS observations. The source is strongly variable on a daily timescale across all wavebands, from optical to TeV, with a peak flux corresponding to ~0.3 times the steady Crab Nebula flux at TeV energies. Follow-up observations in the optical and X-ray bands classify the newly-discovered TeV source as a BL Lac-type blazar with uncertain redshift, although recent measurements suggest z=0.108. VER J0521+211 exhibits all the defining properties of blazars in radio, optical, X-ray, and gamma-ray wavelengths.

We report on the VERITAS observations of the high-frequency peaked BL Lac object 1ES 1959+650 in the period 2007-2011. This source is detected at TeV energies by VERITAS at 16.4 standard deviation (sigma) significance in 7.6 hours of observation in a low flux state. A multiwavelength spectral energy distribution (SED) is constructed from contemporaneous data from VERITAS, Fermi-LAT, RXTE PCA, and Swift UVOT. Swift XRT data is not included in the SED due to a lack of simultaneous observations with VERITAS. In contrast to the orphan gamma-ray flare exhibited by this source in 2002, the X-ray flux of the source is found to vary by an order of magnitude, while other energy regimes exhibit less variable emission. A quasi-equilibrium synchrotron self-Compton model with an additional external radiation field is used to describe three SEDs corresponding to the lowest, highest, and average X-ray states. The variation in the X-ray spectrum is modeled by changing the electron injection spectral index, with minor adjustments of the kinetic luminosity in electrons. This scenario produces small-scale flux variability of order >~2 in the HE (E>1 MeV) and VHE (E>100 GeV) gamma-ray regimes, which is corroborated by the Fermi-LAT, VERITAS, and Whipple 10m telescope light curves.

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 the discovery of TeV gamma-ray emission coincident with the shell-type radio supernova remnant (SNR) CTA 1 using the VERITAS gamma-ray observatory. The source, VER J0006+729, was detected as a 6.5 standard deviation excess over background and shows an extended morphology, approximated by a two-dimensional Gaussian of semi-major (semi-minor) axis 0.30 degree (0.24 degree) and a centroid 5' from the Fermi gamma-ray pulsar PSR J0007+7303 and its X-ray pulsar wind nebula (PWN). The photon spectrum is well described by a power-law dN/dE = N_0 (E/3 TeV)^(-\Gamma), with a differential spectral index of \Gamma = 2.2 +- 0.2_stat +- 0.3_sys, and normalization N_0 = (9.1 +- 1.3_stat +- 1.7_sys) x 10^(-14) cm^(-2) s^(-1) TeV^(-1). The integral flux, F_\gamma = 4.0 x 10^(-12) erg cm^(-2) s^(-1) above 1 TeV, corresponds to 0.2% of the pulsar spin-down power at 1.4 kpc. The energetics, co-location with the SNR, and the relatively small extent of the TeV emission strongly argue for the PWN origin of the TeV photons. We consider the origin of the TeV emission in CTA 1.

We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.

We report on VERITAS very-high-energy (VHE; E>100 GeV) observations of six blazars selected from the Fermi Large Area Telescope First Source Catalog (1FGL). The gamma-ray emission from 1FGL sources was extrapolated up to the VHE band, taking gamma-ray absorption by the extragalactic background light into account. This allowed the selection of six bright, hard-spectrum blazars that were good candidate TeV emitters. Spectroscopic redshift measurements were attempted with the Keck Telescope for the targets without Sloan Digital Sky Survey (SDSS) spectroscopic data. No VHE emission is detected during the observations of the six sources described here. Corresponding TeV upper limits are presented, along with contemporaneous Fermi observations and non-concurrent Swift UVOT and XRT data. The blazar broadband spectral energy distributions (SEDs) are assembled and modeled with a single-zone synchrotron self-Compton model. The SED built for each of the six blazars show a synchrotron peak bordering between the intermediate- and high-spectrum-peak classifications, with four of the six resulting in particle-dominated emission regions.

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.

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated, suggesting that clusters may also be sources of very high-energy (VHE; E>100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cherenkov telescopes, with complementing Fermi-LAT observations at GeV energies. No significant gamma-ray emission from the Coma cluster was detected. Integral flux upper limits at the 99% confidence level were measured to be on the order of (2-5)*10^-8\ ph. m^-2 s^-1 (VERITAS, >220 GeV} and ~2*10^-6 ph. m^-2 s^-1 (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain CRs and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be < 16% from VERITAS data and < 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be <50%. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of (2 - 5.5) muG, depending on the radial magnetic-field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this {renders} the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark-matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally-averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, <\sigma v>. (abr.)

We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between January 2008 and February 2011, resulting in 56.2 hours of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as 5.2 +/- 1.1_stat +/- 2.6_sys * 10^-12 photons cm^-2 s^-1 above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an photon index of Gamma 3.4 +/- 0.5_stat +/- 0.3_sys and a flux normalization of 1.6 +/- 0.3_stat +/- 0.8_sys * 10^-11 photons cm^-2 s^-1 at 300 GeV. We also present multiwavelength results taken in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.

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.

Compilation of papers contributed by the VERITAS Collaboration to the 32nd International Cosmic Ray Conference, held 11-18 August 2011 in Beijing, China.

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.

The VERITAS telescope array has been operating smoothly since 2007, and has detected gamma-ray emission above 100 GeV from 40 astrophysical sources. These include blazars, pulsar wind nebulae, supernova remnants, gamma-ray binary systems, a starburst galaxy, a radio galaxy, the Crab pulsar, and gamma-ray sources whose origin remains unidentified. In 2009, the array was reconfigured, greatly improving the sensitivity. We summarize the current status of the observatory, describe some of the scientific highlights since 2009, and outline plans for the future.

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.

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

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

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.

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.

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

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

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

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

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

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

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

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

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

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

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