Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed.

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $\gamma$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $\gamma$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $\gamma$-ray (VERITAS, MAGIC), high-energy (HE) $\gamma$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Mets\"ahovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($\delta \gtrsim 33$) and the size of the emission region ($ \delta^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

HESS J0632+057 is a gamma-ray binary composed of a compact object orbiting a Be star with a period of about $315$ days. Extensive X-ray and TeV gamma-ray observations have revealed a peculiar light curve containing two peaks, separated by a dip. We present the results of simultaneous observations in hard X-rays with NuSTAR and in TeV gamma-rays with VERITAS, performed in November and December 2017. These observations correspond to the orbital phases $\phi\approx0.22$ and $0.3$, where the fluxes are rising towards the first light-curve peak. A significant variation of the spectral index from 1.77$\pm$0.05 to 1.56$\pm$0.05 is observed in the X-ray data. The multi-wavelength spectral energy distributions (SED) derived from the observations are interpreted in terms of a leptonic model, in which the compact object is assumed to be a pulsar and non-thermal radiation is emitted by high-energy electrons accelerated at the shock formed by the collision between the stellar and pulsar wind. The results of the SED fitting show that our data can be consistently described within this scenario, and allow us to estimate the magnetization of the pulsar wind at the location of the shock formation. The constraints on the pulsar-wind magnetization provided by our results are shown to be consistent with those obtained from other systems.

Compilation of papers presented by the VERITAS Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to observe cosmic neutrinos (CNs) above 20 PeV and ultra-high energy cosmic rays (UHECRs) above 20 EeV over the full sky. The POEMMA mission calls for two identical satellites flying in loose formation, each comprised of a 4-meter wide field-of-view (45 degrees) Schmidt photometer. The hybrid focal surface includes a fast (1 ${\mu}$s) ultraviolet camera for fluorescence observations and an ultrafast (10 ns) optical camera for Cherenkov observations. POEMMA will provide new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these otherwise inaccessible energies.

Efforts to detect ultrahigh energy neutrinos are driven by several objectives: What is the origin of astrophysical neutrinos detected with IceCube? What are the sources of ultrahigh energy cosmic rays? Do the ANITA detected events point to new physics? Shedding light on these questions requires instruments that can detect neutrinos above $10^7$ GeV with sufficient sensitivity - a daunting task. While most ultrahigh energy neutrino experiments are based on the detection of a radio signature from shower particles following a neutrino interaction, we believe that the detection of Cherenkov and fluorescence light from shower particles is an attractive alternative. Imaging air showers with Cherenkov and fluorescence light is a technique that is successfully used in several ultrahigh energy cosmic ray and very-high energy gamma-ray experiments. We performed a case study of an air-shower imaging system for the detection of earth-skimming tau neutrinos. The detector configuration we consider consists of an imaging system that is located on top of a mountain and is pointed at the horizon. From the results of this study we conclude that a sensitivity of $3\cdot10^{-9}$ GeV cm$^{-2}$s$^{-1}$sr$^{-1}$ can be achieved at $2\cdot10^8$ GeV with a relatively small and modular system after three years of observation. In this presentation we discuss key findings of our study and how they translate into design requirements for an imaging system we dub Trinity.

The detection of astrophysical neutrinos by IceCube and the potential to constrain source models of ultra-high energy cosmic rays provide the motivation to develop instruments for the observation of neutrinos above $10^7$ GeV. Among the different techniques to detect ultra-high energy neutrinos is the Earth-skimming technique. It makes use of the fact that the tau produced in a tau neutrino interaction inside the Earth can emerge from the ground and initiate an upward-going particle shower in the atmosphere. The particle shower and thus the neutrino can be reconstructed by measuring the Cherenkov and radio emission from the shower particles. In this presentation, we discuss our ongoing development of a Cherenkov telescope for the detection of tau neutrinos, which is to be deployed on the Extreme Universe Space Observatory Super Pressure Balloon 2 (EUSO-SPB2) and is a precursor experiment for the proposed Probe of Extreme Multi-Messenger Astrophysics (POEMMA) mission. POEMMA aims at the detection of ultrahigh energy cosmic rays and ultrahigh energy neutrinos from low earth orbit. The 1 m$^2$ Cherenkov telescope for EUSO-SPB2 will use silicon photomultipliers coupled to a 100 MS/s readout based on the ASIC for General Electronics for TPC`s (AGET) switch capacitor ring sampler. We present the optics, results from our studies to qualify the readout concept and the design of the mechanical integration of the photon detectors and the readout into the telescope.

Trinity is a proposed air-shower imaging system optimized for the detection of earth-skimming ultrahigh energy tau neutrinos with energies between $10^7$ GeV and $10^{10}$ GeV. Trinity will pursue three major scientific objectives. 1) It will narrow in on possible source classes responsible for the astrophysical neutrino flux measured by IceCube. 2) It will help find the sources of ultrahigh-energy cosmic rays (UHECR) and understand the composition of UHECR. 3) It will test fundamental neutrino physics at the highest energies. Trinity uses the imaging technique, which is well established and successfully used by the very high-energy gamma-ray community (CTA, H.E.S.S., MAGIC, and VERITAS) and the UHECR community (Telescope Array, Pierre Auger)

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is a NASA Astrophysics probe-class mission designed to observe ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos from space. Astro2020 APC white paper: Medium-class Space Particle Astrophysics Project.

We discuss the acceptance and sensitivity of a small air-shower imaging system to detect earth-skimming ultrahigh-energy tau neutrinos. The instrument we study is located on top of a mountain and has an azimuthal field of view of $360^\circ$. We find that the acceptance and sensitivity of such a system is close to maximal if it is located about 2 km above ground, has a vertical field of view of $5^\circ$, allows the reconstruction of an at least $0.3^\circ$ long air-shower image, and features an effective light-collection area of $10$ m$^2$ in any direction. After three years of operation, an imaging system with these features achieves an all-flavor neutrino flux sensitivity of $5\times10^{-9}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ at $2\times10^8$ GeV.

We conduct a search for periodic emission in the very high-energy gamma-ray band (VHE; E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hours. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in non-detections of pulsed VHE gamma rays from each pulsar. Upper limits on a potential VHE gamma-ray flux are derived at the 95% confidence level above three energy thresholds using two methods. These are the first such searches for pulsed VHE emission from each of the pulsars, and the obtained limits constrain a possible flux component manifesting at VHEs as is seen for the Crab pulsar.

We present characterization results of two silicon photomultipliers; the Hamamatsu LVR-6050-CN and the Ketek PM3325 WB. With our measurements of the bias dependence of the breakdown probability we are able to draw conclusions about the location and spatial extension of the avalanche region. For the KETEK SiPM we find that the avalanche region is located close to the surface. In the Hamamatsu SiPM the high-field region is located $0.5\,\mu$m below the surface, while the volume above is depleted almost until the surface. Furthermore, for the Hamamatsu SiPM we find that charge carriers produced by optical-crosstalk photons enter a cell below the avalanche region as opposed to an earlier device where most of the photons enter a cell from above. In the here tested Hamamatsu device the crosstalk photons probably absorb in the bulk and the generated holes diffuse into the active volume of the cell within 2\,ns and initiate a breakdown. The present paper is an attempt to spur further interest in the use of the bias dependence of the breakdown probability and establish it as a standard tool not only to determine the location of the high-field region but also to determine the origin of charge carriers relative to the high-field region. With the knowledge of where the charges come from it should be possible to further improve the optical crosstalk, dark count, and afterpulsing characteristics of SiPM.

We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new, point-like, gamma-ray source is detected, coincident with the location of PSR J2032+4127 / MT91 213. The gamma-ray light curve and spectrum are well-characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar / Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope (XRT) on board the Neil Gehrels \textit{Swift} Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering processes while propagating within the Galaxy and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance there are only a few known astrophysical objects capable of accelerating electrons to such high energies. It is also possible that the CREs are the products of the annihilation or decay of heavy dark matter (DM) particles. VERITAS, an array of imaging air Cherenkov telescopes in southern Arizona, USA, is primarily utilized for gamma-ray astronomy, but also simultaneously collects CREs during all observations. We describe our methods of identifying CREs in VERITAS data and present an energy spectrum, extending from 300 GeV to 5 TeV, obtained from approximately 300 hours of observations. A single power-law fit is ruled out in VERITAS data. We find that the spectrum of CREs is consistent with a broken power law, with a break energy at 710 $\pm$ 40$_{stat}$ $\pm$ 140$_{syst}$ GeV.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 to 500 TeV. The measurement makes use of a template-based analysis method, which, for the first time, is applied to the energy reconstruction of iron-induced air showers recorded by the VERITAS array of imaging atmospheric Cherenkov telescopes. The event selection makes use of the direct Cherenkov light which is emitted by charged particles before the first interaction, as well as other parameters related to the shape of the recorded air shower images. The measured spectrum is well described by a power law $\frac{\mathrm{d} F}{\mathrm{d} E}=f_0\cdot \left(\frac{E}{E_0}\right)^{-\gamma}$ over the full energy range, with $\gamma = 2.82 \pm 0.30 \mathrm{(stat.)} ^{+0.24}_{-0.27} \mathrm{(syst.)}$ and $f_0 = \left( 4.82 \pm 0.98 \mathrm{(stat.)}^{+2.12}_{-2.70} \mathrm{(syst.)} \right)\cdot 10^{-7}$m$^{-2}$s$^{-1}$sr$^{-1}$TeV$^{-1}$ at $E_0=50$TeV, with no indication of a cutoff or spectral break. The measured differential flux is compatible with previous results, with improved statistical uncertainty at the highest energies.

On 2017 September 22, the IceCube Neutrino Observatory reported the detection of the high-energy neutrino event \icnu, of potential astrophysical origin. It was soon determined that the neutrino direction was consistent with the location of the gamma-ray blazar \txs~(3FGL J0509.4+0541), which was in an elevated gamma-ray emission state as measured by the \emph{Fermi} satellite. VERITAS observations of the neutrino/blazar region started on 2017 September 23 in response to the neutrino alert and continued through 2018 February 6. While no significant very-high-energy (VHE; E $>$ 100 GeV) emission was observed from the blazar by VERITAS in the two-week period immediately following the IceCube alert, TXS 0506+056 was detected by VERITAS with a significance of 5.8 standard deviations ($\sigma$) in the full 35-hour data set. The average photon flux of the source during this period was $(8.9 \pm 1.6) \times 10^{-12} \; \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$, or 1.6\% of the Crab Nebula flux, above an energy threshold of 110 GeV, with a soft spectral index of $4.8 \pm 1.3$.

HESS J1943+213 is a very-high-energy (VHE; $>$100 GeV) $\gamma$-ray source in the direction of the Galactic Plane. Studies exploring the classification of the source are converging towards its identification as an extreme synchrotron BL Lac object. Here we present 38 hours of VERITAS observations of HESS J1943+213 taken over two years. The source is detected with $\sim$20 standard deviations significance, showing a remarkably stable flux and spectrum in VHE $\gamma$-rays. Multi-frequency very-long-baseline array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with European VLBI Network and detect this component in the 4.6 GHz and the 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, $Fermi$-LAT, $Swift$-XRT, FLWO 48$''$ telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron-self-Compton model. The well-measured $\gamma$-ray peak of the SED with VERITAS and $Fermi$-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary $\gamma$-rays from ultra-high-energy cosmic ray-initiated electromagnetic cascades to the $\gamma$-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and $\gamma$-ray bands. No statistically significant flux or spectral variability is detected.

We present results from deep observations towards the Cygnus region using 300 hours of very-high-energy (VHE) $\gamma$-ray data taken with the VERITAS Cherenkov telescope array and over seven years of high-energy $\gamma$-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse $\gamma$-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130), and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL 2021.0+4031e (the Gamma-Cygni supernova remnant) was examined and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE $\gamma$-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula ux. We interpret these observations in a multiwavelength context and present the most detailed $\gamma$-ray view of the region to date.

Combined with very-long-baseline interferometry measurements, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL~Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL~Lacertae observed by VERITAS, with a rise time of $\sim$2.3~hr and a decay time of $\sim$36~min. The peak flux above 200 GeV is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photon} \;\text{m}^{-2}\; \text{s}^{-1}$ measured with a 4-minute-binned light curve, corresponding to $\sim$180\% of the flux which is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in VLBA observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models which invoke relativistic plasma passing stationary shocks.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

Imaging Atmospheric Cherenkov Telescopes study the highest energy (up to tens of TeV) photon emission coming from nearby and distant astrophysical sources, thus providing valuable results from searches for Lorentz Invariance Violation (LIV) effects. Highly variable, energetic and distant sources such as Pulsars and AGNs are the best targets for the Time-of-Flight LIV studies. However, the limited number of observations of AGN flares or of high-energy pulsed emission greatly restricts the potential of such studies, especially any potential LIV effects as a function of redshift. To address these issues, an inter-experiment working group has been established by the three major collaborations taking data with Imaging Atmospheric Cherenkov Telescopes (H.E.S.S., MAGIC and VERITAS) with the aim to increase sensitivity to any effects of LIV, together with an improved control of systematic uncertainties, by sharing data samples and developing joint analysis methods. This will allow an increase in the number of available sources and to perform a sensitive search for redshift dependencies. This presentation reviews the first combined maximum likelihood method analyses using simu- lations of published source observations done in the past with H.E.S.S., MAGIC and VERITAS. The results from analyses based on combined maximum likelihood methods, the strategies to deal with data from different types of sources and instruments, as well as future plans will be presented.

Compilation of papers presented by the VERITAS Collaboration at the 35th International Cosmic Ray Conference (ICRC), held July 12 through July 20, 2017 in Busan, South Korea.

Very-high-energy (VHE; $>$ 100 GeV) gamma-ray emission from the blazar RGB J2243+203 was discovered with the VERITAS Cherenkov telescope array, during the period between 21 and 24 December 2014. The VERITAS energy spectrum from this source can be fit by a power law with a photon index of $4.6 \pm 0.5$, and a flux normalization at 0.15 TeV of $(6.3 \pm 1.1) \times 10^{-10} ~ \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{TeV}^{-1}$. The integrated \textit{Fermi}-LAT flux from 1 GeV to 100 GeV during the VERITAS detection is $(4.1 \pm 0.8) \times 10^{\textrm{-8}} ~\textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$, which is an order of magnitude larger than the four-year-averaged flux in the same energy range reported in the 3FGL catalog, ($4.0 \pm 0.1 \times 10^{\textrm{-9}} ~ \textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$). The detection with VERITAS triggered observations in the X-ray band with the \textit{Swift}-XRT. However, due to scheduling constraints \textit{Swift}-XRT observations were performed 67 hours after the VERITAS detection, not simultaneous with the VERITAS observations. The observed X-ray energy spectrum between 2 keV and 10 keV can be fitted with a power-law with a spectral index of $2.7 \pm 0.2$, and the integrated photon flux in the same energy band is $(3.6 \pm 0.6) \times 10^{-13} ~\textrm{cm}^{-2} \textrm{s}^{-1}$. EBL model-dependent upper limits of the blazar redshift have been derived. Depending on the EBL model used, the upper limit varies in the range from z $<~0.9$ to z $<~1.1$.

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) mission is being designed to establish charged-particle astronomy with ultra-high energy cosmic rays (UHECRs) and to observe cosmogenic tau neutrinos (CTNs). The study of UHECRs and CTNs from space will yield orders-of-magnitude increase in statistics of observed UHECRs at the highest energies, and the observation of the cosmogenic flux of neutrinos for a range of UHECR models. These observations should solve the long-standing puzzle of the origin of the highest energy particles ever observed, providing a new window onto the most energetic environments and events in the Universe, while studying particle interactions well beyond accelerator energies. The discovery of CTNs will help solve the puzzle of the origin of UHECRs and begin a new field of Astroparticle Physics with the study of neutrino properties at ultra-high energies.

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.

We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($\tau^{+}\tau^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.

Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive photomultiplier tube (PMT) cameras. Exposure to high levels of background illumination degrades the efficiency of and potentially destroys these photo-detectors over time, so IACTs cannot be operated in the same configuration in the presence of bright moonlight as under dark skies. Since September 2012, observations have been carried out with the VERITAS IACTs under bright moonlight (defined as about three times the night-sky-background (NSB) of a dark extragalactic field, typically occurring when Moon illumination > 35%) in two observing modes, firstly by reducing the voltage applied to the PMTs and, secondly, with the addition of ultra-violet (UV) bandpass filters to the cameras. This has allowed observations at up to about 30 times previous NSB levels (around 80% Moon illumination), resulting in 30% more observing time between the two modes over the course of a year. These additional observations have already allowed for the detection of a flare from the 1ES 1727+502 and for an observing program targeting a measurement of the cosmic-ray positron fraction. We provide details of these new observing modes and their performance relative to the standard VERITAS observations.

We report about the optical and electrical characterization of three high efficiency and blue sensitive Silicon photomultipliers from FBK, Hamamatsu, and SensL. Key features of the tested devices when operated at 90% breakdown probability are peak photon detection efficiencies between 40% and 55%, temperature dependencies of gain and PDE that are less than 1%/$^{\circ}$C, dark rates of $\sim$50\,kHz/mm$^{2}$ at room temperature, afterpulsing of about 2%, and direct optical crosstalk between 6% and 20%. The characteristics of all three devices impressively demonstrate how the Silicon-photomultiplier technology has improved over the past ten years. It is further demonstrated how the voltage and temperature characteristics of a number of quantities can be parameterized on the basis of physical models. The models provide a deeper understanding of the device characteristics over a wide bias and temperature range. They also serve as examples how producers could provide the characteristics of their SiPMs to users. A standardized parameterization of SiPMs would enable users to find the optimal SiPM for their application and the operating point of SiPMs without having to perform measurements thus significantly reducing design and development cycles.

High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it is a young, type Ia SNR that is well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho's SNR with 147 hours of VERITAS and 84 months of Fermi-LAT observations, which represents about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is $2.92 \pm 0.42_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$. It is also softer than the spectral index in the GeV energy range, $2.14 \pm 0.09_{\mathrm{stat}} \pm 0.02_{\mathrm{sys}}$, measured by this study using Fermi--LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi--LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously.

B2 1215+30 is a BL Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes, and subsequently confirmed by the VERITAS observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2 1215+30 during routine monitoring observations of the blazar 1ES 1218+304, located in the same field of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability timescale of < 3.6 h. Multiwavelength observations with Fermi-LAT, Swift, and the Tuorla observatory revealed a correlated high GeV flux state and no significant optical counterpart to the flare, with a spectral energy distribution where the gamma-ray luminosity exceeds the synchrotron luminosity. When interpreted in the framework of a one-zone leptonic model, the observed emission implies a high degree of beaming, with Doppler factor > 10, and an electron population with spectral index < 2.3.

We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGN), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma rays from AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background (CMB) photons via inverse-Compton scattering to produce gamma rays. Due to the deflection of the electron-positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218+304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around $10^{-14}$G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.

We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015.

The Cherenkov Telescope Array (CTA) is planned to be the next-generation very-high-energy (VHE; E > 100 GeV) gamma-ray observatory. It is anticipated that CTA will improve upon the sensitivity of the current generation of VHE experiments, such as VERITAS, HESS and MAGIC, by an order of magnitude. CTA is planned to consist of two graded arrays of Cherenkov telescopes with three primary-mirror sizes. A proof-of-concept telescope, based on the dual-mirror Schwarzchild-Couder design, is being constructed on the VERITAS site at the F.L. Whipple Observatory in southern Arizona, USA, and is a candidate design for the medium-sized telescopes. The construction of the telescope will be completed in early 2017, and the status of this project is presented here.

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.

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.

The F-type star KIC 8462852 has recently been identified as an exceptional target for SETI (search for extraterrestrial intelligence) observations. We describe an analysis methodology for optical SETI, which we have used to analyse nine hours of serendipitous archival observations of KIC 8462852 made with the VERITAS gamma-ray observatory between 2009 and 2015. No evidence of pulsed optical beacons, above a pulse intensity at the Earth of approximately 1 photon per m^2, is found. We also discuss the potential use of imaging atmospheric Cherenkov telescope arrays in searching for extremely short duration optical transients in general.

The TeV binary system LS I +61$^\circ$ 303 is known for its regular, non-thermal emission pattern which traces the orbital period of the compact object in its 26.5 day orbit around its B0 Ve star companion. The system typically presents elevated TeV emission around apastron passage with flux levels between 5% and 15% of the steady flux from the Crab Nebula (> 300 GeV). In this article, VERITAS observations of LS I +61$^\circ$ 303 taken in late 2014 are presented, during which bright TeV flares around apastron at flux levels peaking above 30% of the Crab Nebula flux were detected. This is the brightest such activity from this source ever seen in the TeV regime. The strong outbursts have rise and fall times of less than a day. The short timescale of the flares, in conjunction with the observation of 10 TeV photons from LS I +61$^\circ$ 303 during the flares, provides constraints on the properties of the accelerator in the source.

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.

We present the development of a novel 11328 pixel silicon photomultiplier (SiPM) camera for use with a ground-based Cherenkov telescope with Schwarzschild-Couder optics as a possible medium-sized telescope for the Cherenkov Telescope Array (CTA). The finely pixelated camera samples air-shower images with more than twice the optical resolution of cameras that are used in current Cherenkov telescopes. Advantages of the higher resolution will be a better event reconstruction yielding improved background suppression and angular resolution of the reconstructed gamma-ray events, which is crucial in morphology studies of, for example, Galactic particle accelerators and the search for gamma-ray halos around extragalactic sources. Packing such a large number of pixels into an area of only half a square meter and having a fast readout directly attached to the back of the sensors is a challenging task. For the prototype camera development, SiPMs from Hamamatsu with through silicon via (TSV) technology are used. We give a status report of the camera design and highlight a number of technological advancements that made this development possible.

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.