We present very-high-energy $\gamma$-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above background of $20.8\sigma$ in $47.2$ hours (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Using these observations the temporal properties of 1ES 2344+514 are studied on short and long times scales. We fit a constant flux model to nightly- and seasonally-binned light curves and apply a fractional variability test, to determine the stability of the source on different timescales. We reject the constant-flux model for the 2007-2008 and 2014-2015 nightly-binned light curves and for the long-term seasonally-binned light curve at the $> 3\sigma$ level. The spectra of the time-averaged emission before and after correction for attenuation by the extragalactic background light are obtained. The observed time-averaged spectrum above 200 GeV is satisfactorily fitted (${\chi^2/NDF = 7.89/6}$) by a power-law function with index $\Gamma = 2.46 \pm 0.06_{stat} \pm 0.20_{sys} $ and extends to at least 8 TeV. The extragalactic-background-light-deabsorbed spectrum is adequately fit (${\chi^2/NDF = 6.73/6}$) by a power-law function with index $\Gamma = 2.15 \pm 0.06_{stat} \pm 0.20_{sys} $ while an F-test indicates that the power-law with exponential cutoff function provides a marginally-better fit ($\chi^2/NDF $ = $2.56 / 5 $) at the 2.1$\sigma$ level. The source location is found to be consistent with the published radio location and its spatial extent is consistent with a point source.

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

We present results from multiwavelength observations of the BL Lacertae object 1ES 1741+196, including results in the very-high-energy $\gamma$-ray regime using the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The VERITAS time-averaged spectrum, measured above 180 GeV, is well-modelled by a power law with a spectral index of $2.7\pm0.7_{\mathrm{stat}}\pm0.2_{\mathrm{syst}}$. The integral flux above 180 GeV is $(3.9\pm0.8_{\mathrm{stat}}\pm1.0_{\mathrm{syst}})\times 10^{-8}$ m$^{-2}$ s$^{-1}$, corresponding to 1.6% of the Crab Nebula flux on average. The multiwavelength spectral energy distribution of the source suggests that 1ES 1741+196 is an extreme-high-frequency-peaked BL Lacertae object. The observations analysed in this paper extend over a period of six years, during which time no strong flares were observed in any band. This analysis is therefore one of the few characterizations of a blazar in a non-flaring state.

Between the beginning of its full-scale scientific operations in 2007 and 2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of these, 26 were detected as very-high-energy (VHE; E>100 GeV) {\gamma}-ray sources. In this work, we present the analysis results of a sample of 114 undetected objects. The observations constitute a total live-time of ~570 hours. The sample includes several unidentified Fermi-Large Area Telescope (LAT) sources (located at high Galactic latitude) as well as all the sources from the second Fermi-LAT catalog which are contained within the field of view of the VERITAS observations. We have also performed optical spectroscopy measurements in order to estimate the redshift of some of these blazars that do not have a spectroscopic distance estimate. We present new optical spectra from the Kast instrument on the Shane telescope at the Lick observatory for 18 blazars included in this work, which allowed for the successful measurement or constraint on the redshift of four of them. For each of the blazars included in our sample we provide the flux upper limit in the VERITAS energy band. We also study the properties of the significance distributions and we present the result of a stacked analysis of the data-set, which shows a 4 {\sigma} excess.

Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EBL), or the radiation near the supermassive black hole lying at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS 1441+25 (z=0.939) during April 2015, a period of high activity across all wavelengths. This observation of PKS 1441+25 suggests that the emission region is located thousands of Schwarzschild radii away from the black hole. The gamma-ray detection also sets a stringent upper limit on the near-ultraviolet to near-infrared EBL intensity, suggesting that galaxy surveys have resolved most, if not all, of the sources of the EBL at these wavelengths.

Compilation of papers presented by the VERITAS Collaboration at the 34th International Cosmic Ray Conference (ICRC), held July 30 through August 6, 2015 in The Hague, The Netherlands.

The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is a ground-based array located at the Fred Lawrence Whipple Observatory in southern Arizona and is one of the world's most sensitive gamma-ray instruments at energies of 85 GeV to $>$30 TeV. VERITAS has a wide scientific reach that includes the study of extragalactic and Galactic objects as well as the search for astrophysical signatures of dark matter and the measurement of cosmic rays. In this paper, we will summarize the current status of the VERITAS observatory and present some of the scientific highlights from the last two years, focusing in particular on those results shown at the 2015 ICRC in The Hague, Netherlands.

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.

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

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.

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

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.

The Crab Nebula has long been considered a standard candle in high energy astrophysics, but in recent years this assumption has been strongly contradicted in keV-GeV wavebands. In light of these developments, a search for variability is being performed on the nebula at Very High Energies (VHE; E>300 GeV), the preliminary results of which are presented here. This initial study is based on 10 years (2001-2011) of archival data from the Whipple 10 m telescope. The data set was searched for evidence of variability on the timescales of 1, 7, and 14 days. To date, no significant flaring activity has been found, but simulations are in progress to determine the level of variability that would be detected.

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