Demonstration of stellar intensity interferometry with the four VERITAS telescopes

High angular resolution observations at optical wavelengths provide valuable insights in stellar astrophysics, directly measuring fundamental stellar parameters, and probing stellar atmospheres, circumstellar disks, elongation of rapidly rotating stars, and pulsations of Cepheid variable stars. 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 $\beta$ Canis Majoris and $\epsilon$ 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.Comment: Accepted for publication in Nature Astronomy (2020

An Archival Search for Neutron-Star Mergers in Gravitational Waves and Very-High-Energy Gamma Rays

The recent discovery of electromagnetic signals in coincidence with neutron-star mergers has solidified the importance of multimessenger campaigns in studying the most energetic astrophysical events. Pioneering multimessenger observatories, such as LIGO/Virgo and IceCube, record many candidate signals below the detection significance threshold. These sub-threshold event candidates are promising targets for multimessenger studies, as the information provided by them may, when combined with contemporaneous gamma-ray observations, lead to significant detections. Here we describe a new method that uses such candidates to search for transient events using archival very-high-energy gamma-ray data from imaging atmospheric Cherenkov telescopes (IACTs). We demonstrate the application of this method to sub-threshold binary neutron star (BNS) merger candidates identified in Advanced LIGO's first observing run. We identify eight hours of archival VERITAS observations coincident with seven BNS merger candidates and search them for TeV emission. No gamma-ray emission is detected; we calculate upper limits on the integral flux and compare them to a short gamma-ray burst model. We anticipate this search method to serve as a starting point for IACT searches with future LIGO/Virgo data releases as well as in other sub-threshold studies for multimessenger transients, such as IceCube neutrinos. Furthermore, it can be deployed immediately with other current-generation IACTs, and has the potential for real-time use that places minimal burden on experimental operations. Lastly, this method may serve as a pilot for studies with the Cherenkov Telescope Array, which has the potential to observe even larger fields of view in its divergent pointing mode

Deep Learning Transient Detection with VERITAS

Ground-based $\gamma$-ray observatories, such as the VERITAS array of imaging atmospheric Cherenkov telescopes, provide insight into very-high-energy (VHE, $\mathrm{E}>100\,\mathrm{GeV}$) astrophysical transient events. Examples include the evaporation of primordial black holes, gamma-ray bursts and flaring blazars. Identifying such events with a serendipitous location and time of occurrence is difficult. Thus, employing a robust search method becomes crucial. An implementation of a transient detection method based on deep-learning techniques for VERITAS will be presented. This data-driven approach significantly reduces the dependency on the characterization of the instrument response and the modelling of the expected transient signal. The response of the instrument is affected by various factors, such as the elevation of the source and the night sky background. The study of these effects allows enhancing the deep learning method with additional parameters to infer their influences on the data. This improves the performance and stability for a wide range of observational conditions. We illustrate our method for an historic flare of the blazar BL Lac that was detected by VERITAS in October 2016. We find a promising performance for the detection of such a flare in timescales of minutes that compares well with the VERITAS standard analysis

A Search for TeV Gamma-Ray Emission from Pulsar Tails by VERITAS

We report on the search for very-high-energy gamma-ray emission from the regions around three nearby supersonic pulsars (PSR B0355+54, PSR J0357+3205, and PSR J1740+1000) that exhibit long X-ray tails. To date there is no clear detection of TeV emission from any pulsar tail that is prominent in X-ray or radio. We provide upper limits on the TeV flux, and luminosity, and also compare these limits with other pulsar wind nebulae detected in X-rays and the tail emission model predictions. We find that at least one of the three tails is likely to be detected in observations that are a factor of 2–3 more sensitive. The analysis presented here also has implications for deriving the properties of pulsar tails, for those pulsars whose tails could be detected in TeV

Variability and Spectral Characteristics of Three Flaring Gamma-Ray Quasars Observed by VERITAS and Fermi-LAT

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

Variability and Spectral Characteristics of Three Flaring Gamma-Ray Quasars Observed by VERITAS and Fermi-LAT

International audienceFlat-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

VTSCat: The VERITAS Catalog of Gamma-Ray Observations

The ground-based gamma-ray observatory Very Energetic Radiation Imaging Telescope Array System (VERITAS, https://veritas.sao.arizona.edu/) is sensitive to photons of astrophysical origin with energies in the range between ≈85 GeV and ≈30 TeV. The instrument consists of four 12 m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hr of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star-forming galaxies, dwarf-spheroidal galaxies, and highly variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications

Gamma-ray flaring activity of NGC 1275 in 2016-2017 measured by MAGIC

We report on the detection of flaring activity from the Fanaroff-Riley I radio galaxy NGC 1275 in very-high-energy (VHE, E $>$ 100 GeV) gamma rays with the MAGIC telescopes. Observations were performed between 2016 September and 2017 February as part of a monitoring program. The brightest outburst with $\sim1.5$ times the Crab Nebula flux above 100 GeV (C.U.) was observed during the night between 2016 December 31 and 2017 January 1 (fifty times higher than the mean previously measured in two observational campaigns between 2009 and 2011). Significant variability of the day-by-day light curve was measured, the shortest flux-doubling time-scales was found to be of $(611\pm101)$ min. The combined spectrum of the MAGIC data during the strongest flare state and simultaneous data from the Fermi-LAT around 2017 January 1 follows a power-law with an exponential cutoff at the energy $(492\pm35)$ GeV. Simultaneous optical flux density measurements in the R-band obtained with the KVA telescope are also presented and the correlation between the optical and gamma-ray emission is investigated. Due to possible internal pair-production, the fast flux variability constrains the Doppler factor to values which are inconsistent with a large viewing angle as observed in the radio band. We investigate different scenarios for the explanation of fast gamma-ray variability, namely emission from: magnetospheric gaps, relativistic blobs propagating in the jet (mini-jets) or external cloud (or star) entering the jet. We find that the only plausible model to account for the luminosities here observed would be the production of gamma rays in a magnetospheric gap around the central black hole only in the eventuality of an enhancement of the magnetic field threading the hole from its equipartition value with the gas pressure in the accretion flow

Multiwavelength Observations of the Blazar VER J0521+211 during an Elevated TeV Gamma-Ray State

We report on a long-lasting, elevated gamma-ray flux state from VER J0521+211 observed by VERITAS, MAGIC, and Fermi-LAT in 2013 and 2014. The peak integral flux above 200 GeV measured with the nightly binned light curve is (8.8 ± 0.4) × 10$^{−7}$ photons m$^{−2}$ s$^{−1}$, or ∼37% of the Crab Nebula flux. Multiwavelength observations from X-ray, UV, and optical instruments are also presented. A moderate correlation between the X-ray and TeV gamma-ray fluxes was observed, and the X-ray spectrum appeared harder when the flux was higher. Using the gamma-ray spectrum and four models of the extragalactic background light (EBL), a conservative 95% confidence upper limit on the redshift of the source was found to be z ≤ 0.31. Unlike the gamma-ray and X-ray bands, the optical flux did not increase significantly during the studied period compared to the archival low-state flux. The spectral variability from optical to X-ray bands suggests that the synchrotron peak of the spectral energy distribution (SED) may become broader during flaring states, which can be adequately described with a one-zone synchrotron self-Compton model varying the high-energy end of the underlying particle spectrum. The synchrotron peak frequency of the SED and the radio morphology of the jet from the MOJAVE program are consistent with the source being an intermediate-frequency-peaked BL Lac object

Combined dark matter searches towards dwarf spheroidal galaxies with Fermi-LAT, HAWC, H.E.S.S., MAGIC, and VERITAS

Cosmological and astrophysical observations suggest that 85\% of the total matter of the Universe is made of Dark Matter (DM). However, its nature remains one of the most challenging and fundamental open questions of particle physics. Assuming particle DM, this exotic form of matter cannot consist of Standard Model (SM) particles. Many models have been developed to attempt unraveling the nature of DM such as Weakly Interacting Massive Particles (WIMPs), the most favored particle candidates. WIMP annihilations and decay could produce SM particles which in turn hadronize and decay to give SM secondaries such as high energy $\gamma$ rays. In the framework of indirect DM search, observations of promising targets are used to search for signatures of DM annihilation. Among these, the dwarf spheroidal galaxies (dSphs) are commonly favored owing to their expected high DM content and negligible astrophysical background. In this work, we present the very first combination of 20 dSph observations, performed by the Fermi-LAT, HAWC, H.E.S.S., MAGIC, and VERITAS collaborations in order to maximize the sensitivity of DM searches and improve the current results. We use a joint maximum likelihood approach combining each experiment's individual analysis to derive more constraining upper limits on the WIMP DM self-annihilation cross-section as a function of DM particle mass. We present new DM constraints over the widest mass range ever reported, extending from 5 GeV to 100 TeV thanks to the combination of these five different $\gamma$-ray instruments