Matthews a,d on behalf of the VERITAS Collaboration

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.

The results of gamma-ray observations of the binary system HESS J0632 + 057 collected during 450 hr over 15 yr, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the very high-energy gamma-ray fluxes with a period of 316.7 ± 4.4 days is reported, consistent with the period of 317.3 ± 0.7 days obtained with a refined analysis of X-ray data. The analysis of data from four orbital cycles with dense observational coverage reveals short-timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over a timescale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but cannot find any correlation of optical Hα parameters with fluxes at X-ray or gamma-ray energies in simultaneous observations. The key finding is that the emission of HESS J0632 + 057 in the X-ray and gamma-ray energy bands is highly variable on different timescales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.

In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass ∼6.5 × 10 9 M e . The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87's spectrum. We can exclude that the simultaneous γ-ray emission is produced via inverse Compton emission in the same region producing the

An accurate modeling of neutrino flux attenuation and the distribution of leptons they produce in transit through the Earth is an essential component to determine neutrino flux sensitivities of underground, sub-orbital and space-based detectors. Through neutrino oscillations over cosmic distances, astrophysical neutrino sources are expected to produce nearly equal fluxes of electron, muon and tau neutrinos. Of particular interest are tau neutrinos that interact in the Earth at modest slant depths to produce -leptons. Some -leptons emerge from the Earth and decay in the atmosphere to produce extensive air showers. Future balloon-borne and satellite-based optical Cherenkov neutrino telescopes will be sensitive to upward air showers from tau neutrino induced -lepton decays. We present nuPyProp, a python code that is part of the nuSpaceSim package.nuPyProp generates look-up tables for exit probabilities and energy distributions for → and → propagation in the Earth. This flexible code runs with either stochastic or continuous electromagnetic energy losses for the lepton transit through the Earth. Current neutrino cross section models and energy loss models are included along with templates for user input of other models. Results from nuPyProp are compared with other recent simulation packages for neutrino and charged lepton propagation. Sources of modeling uncertainties are described and quantified.