Matthews a,d on behalf of the VERITAS Collaboration

Previous detections of individual astrophysical sources of neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017, we detected a high-energy neutrino, IceCube-170922A, with an energy of ~290 tera-electron volts. Its arrival direction was consistent with the location of a known γ-ray blazar, TXS 0506+056, observed to be in a flaring state. An extensive multiwavelength campaign followed, ranging from radio frequencies to γ-rays. These observations characterize the variability and energetics of the blazar and include the detection of TXS 0506+056 in very-high-energy γ-rays. This observation of a neutrino in spatial coincidence with a γ-ray-emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (>  | | b 10), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL Lacs. The most abundant detected BL Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical, and X-ray flux distributions, which hints that even the faintest known blazars could eventually shine in gamma-rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

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.

We present a catalog of sources detected above 50 GeV by the Fermi-Large Area Telescope (LAT) in 80 months of data. The newly delivered Pass8 event-level analysis allows the detection and characterization of sources in the 50 GeV-2 TeV energy range. In this energy band, Fermi-LAT has detected 360 sources, which constitute the second catalog of hard Fermi-LAT sources (2FHL). The improved angular resolution enables the precise localization of point sources (∼1 7 radius at 68% C.L.) and the detection and characterization of spatially extended sources. We find that 86% of the sources can be associated with counterparts at other wavelengths, of which the majority (75%) are active galactic nuclei and the rest (11%) are Galactic sources. Only 25% of the 2FHL sources have been previously detected by Cherenkov telescopes, implying that the 2FHL provides a reservoir of candidates to be followed up at very high energies. This work closes the energy gap between the observations performed at GeV energies by Fermi-LATon orbit and the observations performed at higher energies by Cherenkov telescopes from the ground.

We report on the γ-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average Large Area Telescope (LAT) γ-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 ± 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of two), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 ± 0.14, and the softest one is 2.51 ± 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3 GeV. In this paper, we also present the first results from the 4.5 month long multifrequency campaign (2009 March 15-August 1) on Mrk 501, which included the Very Long Baseline Array (VLBA), Swift, RXTE, MAGIC, and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The extensive radio to TeV data set from this campaign provides us with the most detailed spectral energy distribution yet collected for this source during its relatively low activity. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton (SSC) model. In the framework of this model, we find that the dominant emission region is characterized by a size 0.1 pc (comparable within a factor of few to the size of the partially resolved VLBA core at 15-43 GHz), and that the total jet power (10 44 erg s −1) constitutes only a small fraction (∼10 −3) of the Eddington luminosity. The energy distribution of the freshly accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3 GeV-10 TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20 GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks. We find that the ultrarelativistic electrons and mildly relativistic protons within the blazar zone, if comparable in number, are in approximate energy equipartition, with their energy dominating the jet magnetic field energy by about two orders of magnitude.

We report on TeV γ -ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10 m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV γ -ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hr of VERITAS and 96 hr of Whipple 10 m data were acquired between 2006 January and 2008 June. We present the results of a study of the TeV γ -ray energy spectra as a function of time and for different flux levels. On 2008 May 2 and 3, bright TeV γ -ray flares were detected with fluxes reaching the level of 10 Crab. The TeV γ -ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio wave band. The combination of the Rossi X-ray Timing Explorer and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV 1 The Astrophysical Journal, 738:25 (19pp), 2011 September 1 Acciari et al.γ -rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone synchrotron self-Compton model.

We present the first installment of HI sources extracted from the Arecibo Legacy Fast ALFA (ALFALFA) extragalactic survey, initiated in 2005. Sources have been extracted from 3-D spectral data cubes exploiting a matched filtering technique and then examined interactively to yield global HI parameters. A total of 730 HI detections are catalogued within the solid angle 11 h 44 m < R.A.(J2000) < 14 h 00 m and +12 • < Dec.(J2000) < +16 • , and redshift range −1600 km s −1 < cz < 18000 km s −1 . In comparison, the HI Parkes All-Sky Survey (HIPASS) detected 40 HI signals in the same region. Optical counterparts are assigned via examination of digital optical imaging databases. ALFALFA HI detections are reported for three distinct classes of signals: (a) detections, typically with S/N > 6.5; (b) high velocity clouds in the Milky Way or its periphery; and (c) signals of lower S/N (to ∼ 4.5) which coincide spatially with an optical object of known similar redshift. Although this region of the sky has been heavily surveyed by previous targeted observations based on optical flux-or size-limited samples, 69% of the extracted sources are newly reported HI detections. The resultant positional accuracy of HI sources is dependent on S/N: it averages 24 ′′ (20 ′′ median) for all sources with S/N > 6.5 and is of order ∼17 ′′ (14 ′′ median) for signals with S/N > 12. The median redshift of the sample is ∼7000 km s −1 and its distribution reflects the known local large scale structure including the Virgo cluster and the void behind it, the A1367-Coma supercluster at cz ∼7000 km s −1 and a third more distant overdensity at cz ∼13000 km s −1 . Distance uncertainties in and around the Virgo cluster perturb the derived HI mass distribution. Specifically, an apparent deficiency of the lowest HI mass objects can be attributed, at least in part, to the incorrect assignment of some foreground objects to the cluster distance. Several extended HI features are found in the vicinity of the Virgo cluster. A small percentage (6%) of HI detections have no identifiable optical counterpart, more than half of which are high velocity clouds in the Milky Way vicinity; the remaining 17 objects do not appear connected to or associated with any known galaxy. Based on these initial results, ALFALFA is expected to fulfill, and even exceed, its predicted performance objectives in terms of the number and quality of HI detections.