Recent results obtained with gamma-ray satellites have established supernova remnants as accelerators of GeV hadronic cosmic rays. In such processes, CRs accelerated in SNR shocks interact with particles from gas clouds in their surrounding. In particular, the rich medium in which core-collapse SNRs explode provides a large target density to boost hadronic gamma-rays. SNR G39.2-0.3 is one of the brightest SNR in infrared wavelengths, and its broad multi-wavelength coverage allows detailed modeling of its radiation from radio to high energies. We reanalyzed the Fermi-LAT data on this region and compare it with new radio observations from the MWISP survey. The modeling of the spectral energy distribution from radio to GeV energies favors an hadronic origin of the gamma-ray emission and constrains the SNR magnetic field to be at least ~100 uG. Despite the large magnetic field, the present acceleration of protons seems to be limited to ~10 GeV, which points to a drastic slow down of the shock velocity due to the dense wall traced by the CO observations, surrounding the remnant. Further investigation of the gamma-ray spectral shape points to a dynamically old remnant subjected to severe escape of CRs and a decrease of acceleration efficiency. The low-energy peak of the gamma-ray spectrum also suggests that that the composition of accelerated particles might be enriched by heavy nuclei which is certainly expected for a core-collapse SNR. Alternatively, the contribution of the compressed pre-existing Galactic cosmic rays is discussed, which is, however, found to not likely be the dominant process for gamma-ray production.

Depending on their type, supernova remnants may have ejecta material with high abundance of heavy elements such as carbon or oxygen. In addition, core-collapse supernovae explode in the wind material of their progenitor star that may also have a high abundance of heavy elements. Hadronic collisions in these enriched media spawn the production of gamma rays, neutrinos, and secondary electrons with spectra that cannot be scaled from those calculated for pp collisions, potentially leading to erroneous results. We used Monte-Carlo event generators to calculate the differential production rate of particles such as gamma rays, neutrinos, and secondary electrons for H, He, C, and O nuclei as projectiles and as target material. The cross sections and the multiplicity spectra are separately computed for each of the 16 combinations of projectile and target. We describe characteristic effects of heavy nuclei in the shape and normalization of the spectra of various particles produced.

We present a study of $\gamma$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from \textit{Fermi}-LAT and VERITAS deviates significantly from a simple power-law, and is best described by a power-law with spectral index of $2.17\pm 0.02_{stat}$ with a cut-off energy of $2.3 \pm 0.5_{stat}$ TeV. These results, along with radio, X-ray and $\gamma$-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed $\gamma$-ray spectrum. From modeling of the entire multi-wavelength spectrum, a minimum magnetic field inside the remnant of $B_{\mathrm{min}}\approx150\,\mathrm{\mu G}$ is deduced.

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.

With 8 hours of observations, VERITAS confirms the detection of two very high energy gamma-ray sources. The gamma-ray binary LS 5039 is detected with a statistical significance of $8.8\sigma$. The measured flux above 1 TeV is $(2.5 \pm 0.4) \times 10^{-12} \rm \, cm^{-2} \, s^{-1}$ near inferior conjunction and $(7.8 \pm 2.8) \times 10^{-13} \rm \, cm^{-2} \, s^{-1}$ near superior conjunction. The pulsar wind nebula HESS J1825-137 is detected with a statistical significance of $6.7\sigma$ and a measured flux above 1 TeV of $(3.9 \pm 0.8) \times 10^{-12} \rm \, cm^{-2} \, s^{-1}$.

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.

We present 1 to 10GHz radio continuum flux density, spectral index, polarisation and Rotation Measure (RM) images of the youngest known Galactic Supernova Remnant (SNR) G1.9+0.3, using observations from the Australia Telescope Compact Array (ATCA). We have conducted an expansion study spanning 8 epochs between 1984 and 2017, yielding results consistent with previous expansion studies of G1.9+0.3. We find a mean radio continuum expansion rate of ($0.78 \pm 0.09$) per cent year$^{-1}$ (or $\sim8900$ km s$^{-1}$ at an assumed distance of 8.5 kpc), although the expansion rate varies across the SNR perimeter. In the case of the most recent epoch between 2016 and 2017, we observe faster-than-expected expansion of the northern region. We find a global spectral index for G1.9+0.3 of $-0.81\pm0.02$ (76 MHz$-$10 GHz). Towards the northern region, however, the radio spectrum is observed to steepen significantly ($\sim -$1). Towards the two so called (east & west) "ears" of G1.9+0.3, we find very different RM values of 400-600 rad m$^{2}$ and 100-200 rad m$^{2}$ respectively. The fractional polarisation of the radio continuum emission reaches (19 $\pm$ 2)~per~cent, consistent with other, slightly older, SNRs such as Cas~A.

Supernova remnants are known to accelerate cosmic rays on account of their non-thermal emission of radio waves, X-rays, and gamma rays. Although there are many models for the acceleration of cosmic rays in Supernova remnants, the escape of cosmic rays from these sources is yet understudied. We use our time-dependent acceleration code RATPaC to study the acceleration of cosmic rays and their escape in post-adiabatic Supernova remnants and calculate the subsequent gamma-ray emission from inverse-Compton scattering and Pion decay. We performed spherically symmetric 1-D simulations in which we simultaneously solve the transport equations for CRs, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in a volume large enough to keep all CRs in the simulation. The transport equations for cosmic-rays and magnetic turbulence are coupled via the cosmic-ray gradient and the spatial diffusion coefficient of the cosmic rays, while the cosmic-ray feedback onto the shock structure can be ignored. Our simulations span 100kyrs, thus covering the remnants evolution until the beginning of the post-adiabatic phase. At later stages of the evolution cosmic rays over a wide range of energy can reside outside of the remnant, creating spectra that are softer than predicted by standard DSA and feature breaks in the 10-100 GeV-range. The total spectrum of cosmic rays released into the interstellar medium has a spectral index of s~2.4 above roughly 10 GeV which is close to that required by Galactic propagation models. We further find the gamma-ray luminosity to peak around an age of 4,000 years for inverse-Compton-dominated high-energy emission. Remnants expanding in low-density media emit generally more inverse-Compton radiation matching the fact that the brightest known supernova remnants - RCW86, Vela Jr, HESSJ1731-347 and RXJ1713.7-3946 - are all expanding in low density environments.

The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the Universe's history. Spectral measurements of blazars at very high energies ($>$100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the gamma-ray photons with the EBL. The spectra of 14 VERITAS-detected blazars are included in a new measurement of the EBL SED that is independent of EBL SED models. The resulting SED covers an EBL wavelength range of 0.56--56 $\mu$m, and is in good agreement with lower limits obtained by assuming that the EBL is entirely due to radiation from cataloged galaxies.

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.

Context. The youngest Galactic supernova remnant G1.9+0.3 is an interesting target for next generation gamma-ray observatories. So far, the remnant is only detected in the radio and the X-ray bands, but its young age of ~100 yrs and inferred shock speed of ~14,000 km/s could make it an efficient particle accelerator. Aims. We aim to model the observed radio and X-ray spectra together with the morphology of the remnant. At the same time, we aim to estimate the gamma-ray flux from the source and evaluated the prospects of its detection with future gamma-ray experiments. Methods. We performed spherical symmetric 1-D simulations with the RATPaC code, in which we simultaneously solve the transport equation for cosmic rays, the transport equation for magnetic turbulence, and the hydro-dynamical equations for the gas flow. Separately computed distributions of the particles accelerated at the forward and the reverse shock are then used to calculate the spectra of synchrotron, inverse Compton, and pion-decay radiation from the source. Results. The emission from G1.9+0.3 can be self-consistently explained within the test-particle limit. We find that the X-ray flux is dominated by emission from the forward shock while most of the radio emission originates near the reverse shock, which makes G1.9+0.3 the first remnant with non-thermal radiation detected from the reverse shock. The flux of very-high-energy gamma-ray emission from G1.9+0.3 is expected to be close to the sensitivity threshold of the Cherenkov Telescope Array, CTA. The limited time available to grow large-scale turbulence limits the maximum energy of particles to values below 100 TeV, hence G1.9+0.3 is not a PeVatron.

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.

The angular size of a star is a critical factor in determining its basic properties. Direct measurement of stellar angular diameters is difficult: at interstellar distances stars are generally too small to resolve by any individual imaging telescope. This fundamental limitation can be overcome by studying the diffraction pattern in the shadow cast when an asteroid occults a star, but only when the photometric uncertainty is smaller than the noise added by atmospheric scintillation. Atmospheric Cherenkov telescopes used for particle astrophysics observations have not generally been exploited for optical astronomy due to the modest optical quality of the mirror surface. However, their large mirror area makes them well suited for such high-time-resolution precision photometry measurements. Here we report two occultations of stars observed by the VERITAS Cherenkov telescopes with millisecond sampling, from which we are able to provide a direct measurement of the occulted stars' angular diameter at the $\leq0.1$ milliarcsecond scale. This is a resolution never achieved before with optical measurements and represents an order of magnitude improvement over the equivalent lunar occultation method. We compare the resulting stellar radius with empirically derived estimates from temperature and brightness measurements, confirming the latter can be biased for stars with ambiguous stellar classifications.

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.

Vela Jr. (RX J0852.0$-$4622) is one of just a few known supernova remnants (SNRs) with a resolved shell across the whole electromagnetic spectrum from radio to very-high-energy ($>100$ GeV; VHE) gamma-rays. Its proximity and large size allow for detailed spatially resolved observations of the source making Vela Jr. one of the primary sources used for the study of particle acceleration and emission mechanisms in SNRs. High-resolution X-ray observations reveal a steepening of the spectrum toward the interior of the remnant. In this study we aim for a self-consistent radiation model of Vela Jr. which at the same time would explain the broadband emission from the source and its intensity distribution. We solve the full particle transport equation combined with the high-resolution 1D hydrodynamic simulations (using Pluto code) and subsequently calculate the radiation from the remnant. The equations are solved in the test particle regime. We test two models for the magnetic field profile downstream of the shock: damped magnetic field which accounts for the damping of strong magnetic turbulence downstream, and transported magnetic field. Neither of these scenarios can fully explain the observed radial dependence of the X-ray spectrum under spherical symmetry. We show, however, that the softening of the spectrum and the X-ray intensity profile can be explained under the assumption that the emission is enhanced within a cone.

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

We present very-high-resolution 1D MHD simulations of the late-stage supernova remnants (SNR). In the post-adiabatic stage the magnetic field has an important and significant dynamical effect on the shock dynamics, the flow structure, and hence the acceleration and emission of cosmic rays. We find that the tangential component of the magnetic field provides pressure support that to a fair degree prevents the collapse of the radiative shell and thus limits the total compression ratio of the partially or fully radiative forward shock. A consequence is that the spectra of cosmic rays would not be as hard as in hydrodynamic simulations. We also investigated the effect on the flow profiles of the magnetic-field inclination and a large-scale gradient in the gas density and/or the magnetic field. A positive density gradient shortens the evolutionary stages whereas a shock obliquity lowers the shock compression. The compression of the tangential component of the magnetic field leads to its dominance in the downstream region of post-adiabatic shocks for a wide range of orientation of the upstream field, which may explain why one preferentially observes tangential radio polarization in old SNRs. As most cosmic rays are produced at late stages of SNR evolution, the post-adiabatic phase and the influence of the magnetic field during it are most important for modeling the cosmic-ray acceleration at old SNRs and the gamma-ray emission from late-stage SNRs interacting with clouds.

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.

On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant since the very-high-energy (VHE) observation started only $\sim$2 minutes after the prompt emission peaked, and $Fermi$-LAT observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB~150323A ($z=0.593$) limits the attenuation by the extragalactic background light to $\sim 50$ % at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below $\sim100$ GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be $A\gtrsim 3\times 10^{11}$ g cm$^{-1}$, consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the ISM, which therefore cannot be ruled out as the environment of GRB 150323A.

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.

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

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.

Context. Supernova remnants are known as sources of galactic cosmic rays for their non-thermal emission of radio waves, X-rays, and gamma-rays. However, the observed soft broken power-law spectra are hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvenic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a co-moving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces non-linearity in addition to that imposed by cosmic-ray feedback.