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.

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi Large Area Telescope (LAT) catalog of γ-ray sources. Based on the first 12 years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss the degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

Axions are well motivated particles that could make up most or all of the dark matter if they have masses below 100 µeV. Microwave cavity techniques comprised of closed resonant structures immersed in solenoid magnets are sensitive to dark matter axions with masses of a few µeV, but face difficulties scaling to higher masses. We present the a novel detector architecture consisting of an open, Fabry-Pérot resonator and a series of current-carrying wire planes, and demonstrate this technique with a search for dark matter axion-like particles called Orpheus. This search excludes dark matter axion-like particles with masses between 68.2 and 76.5 µeV and axion-photon couplings greater than 4 × 10 −7 GeV −1 . We project that the fundamental sensitivity of this technique could be extended to be sensitive to couplings below 1 × 10 −15 GeV −1 , consistent with the DFSZ model of QCD axions.

A deep survey of the Large Magellanic Cloud at ∼ 0.1−100 TeV photon energies with the Cherenkov Telescope Array is planned. We assess the detection prospects based on a model for the emission of the galaxy, comprising the four known TeV emitters, mock populations of sources, and interstellar emission on galactic scales. We also assess the detectability of 30 Doradus and SN 1987A, and the constraints that can be derived on the nature of dark matter. The survey will allow for fine spectral studies of N 157B, N 132D, LMC P3, and 30 Doradus C, and half a dozen other sources should be revealed, mainly pulsar-powered objects. The remnant from SN 1987A could be detected if it produces cosmic-ray nuclei with a flat power-law spectrum at high energies, or with a steeper index 2.3 − 2.4 pending a flux increase by a factor > 3 − 4 over ∼ 2015 − 2035. Large-scale interstellar emission remains mostly out of reach of the survey if its > 10 GeV spectrum has a soft photon index ∼ 2.7, but degree-scale 0.1 − 10 TeV pion-decay emission could be detected if the cosmic-ray spectrum hardens above >100 GeV. The 30 Doradus star-forming region is detectable if acceleration efficiency is on the order of 1 − 10% of the mechanical luminosity and diffusion is suppressed by two orders of magnitude within < 100 pc. Finally, the survey could probe the canonical velocity-averaged cross section for self-annihilation of weakly interacting massive particles for cuspy Navarro-Frenk-White profiles.

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.

The Cherenkov Telescope Array (CTA) is the major next-generation observa-7 tory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

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 Very Energetic Radiation Imaging Telescope Array and Major Atmospheric Gamma Imaging 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 on board the Neil Gehrels 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 that we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

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.