Matthews a,d on behalf of the VERITAS Collaboration

LS I +61°303 is one of around ten gamma-ray binaries detected so far which has a spectral energy distribution dominated by MeV-GeV photons. It is located at a distance of 2 kpc and consists of a compact object (black hole or neutron star) in an eccentric orbit around a 10-15Be star, with an orbital period of 26.496 days. The binary orbit modulates the emission ranging from radio to TeV energies. A second, longer, modulation period of 1667 days (the super-orbital period) has also been detected from radio to TeV observations. The VERITAS imaging atmospheric Cherenkov telescope array has been observing LS I +61°303 since 2006, and has accumulated a dataset that fully covers the entire orbit. Increased coverage of the source in the very-high-energy band is currently underway to provide more results on the modulation pattern, super-orbital period, and orbit-to-orbit variability at the highest energies. The spectral measurements at the highest energies will reveal more information about gamma-ray production/absorption mechanisms, the nature of the compact object, and the particle acceleration mechanism. Using >150 hrs of VERITAS data, we present a detailed study of the spectral energy distribution and periodic behavior of this rare gamma-ray source type at very-high energy.

Matthews a,d on behalf of the VERITAS Collaboration

High angular resolution observations at optical wavelengths provide valuable insights in stellar astrophysics [1, 2], directly measuring fundamental stellar parameters [3, 4], and probing stellar atmospheres, circumstellar disks [5], elongation of rapidly rotating stars [6], and pulsations of Cepheid variable stars [7]. The angular size of most stars are of order one milli-arcsecond or less, and to spatially resolve stellar disks and features at this scale requires an optical interferometer using an array of telescopes with baselines on the order of hundreds of meters. We report on the successful implementation of a stellar intensity interferometry system developed for the four VERITAS imaging atmospheric-Cherenkov telescopes. The system was used to measure the angular diameter of the two sub-mas stars β Canis Majoris and Orionis with a precision better than 5%. The system utilizes an off-line approach where starlight intensity fluctuations recorded at each telescope are correlated post-observation. The technique can be readily scaled onto tens to hundreds of telescopes, providing a capability that has proven technically challenging to current generation optical amplitude interferometry observatories. This work demonstrates the feasibility of performing astrophysical measurements with imaging atmospheric-Cherenkov telescope arrays as intensity interferometers and the promise for integrating an intensity interferometry system within future observatories such as the Cherenkov Telescope Array.

Flat-spectrum radio quasars (FSRQs) are the most luminous blazars at GeV energies but only rarely emit detectable fluxes of TeV gamma rays, typically during bright GeV flares. We explore the gamma-ray variability and spectral characteristics of three FSRQs that have been observed at GeV and TeV energies by Fermi-LAT and VERITAS, making use of almost 100 hr of VERITAS observations spread over 10 yr: 3C 279, PKS 1222+216, and Ton 599. We explain the GeV flux distributions of the sources in terms of a model derived from a stochastic differential equation describing fluctuations in the magnetic field in the accretion disk and estimate the timescales of magnetic flux accumulation and stochastic instabilities in their accretion disks. We identify distinct flares using a procedure based on Bayesian blocks and analyze their daily and subdaily variability and gamma-ray energy spectra. Using observations from VERITAS, as well as Fermi, Swift, and the Steward Observatory, we model the broadband spectral energy distributions of PKS 1222+216 and Ton 599 during very high energy (VHE)–detected flares in 2014 and 2017, respectively, strongly constraining the jet Doppler factors and gamma-ray emission region locations during these events. Finally, we place theoretical constraints on the potential production of PeV-scale neutrinos during these VHE flares.