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.

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.

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.

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.

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

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.