Showing 1–50 of 66 results for author: Archer, A

Investigation of the Microquasar SS 433 with VERITAS

Authors: The VERITAS Collaboration, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, N. R. Bond, Y. Chen, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, M. E. Godoy, J. E. Pedrosa, S. Feldman, Q. Feng, S. Filbert, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, T. B. Humensky, M. Iskakova, W. Jin, M. N. Johnson , et al. (35 additional authors not shown)

Abstract: Microquasars such as SS 433 are considered potential contributors to cosmic rays up to the knee of the cosmic ray energy spectrum ($\sim4\,\mathrm{PeV}$), where a transition in the dominant acceleration processes is expected. The SS 433 system, located within the W50 supernova remnant, is a Galactic microquasar with relativistic jets interacting with the surrounding medium over parsec scales, prov… ▽ More Microquasars such as SS 433 are considered potential contributors to cosmic rays up to the knee of the cosmic ray energy spectrum ($\sim4\,\mathrm{PeV}$), where a transition in the dominant acceleration processes is expected. The SS 433 system, located within the W50 supernova remnant, is a Galactic microquasar with relativistic jets interacting with the surrounding medium over parsec scales, providing an example for studying jet-driven particle acceleration. A deep morphological and spectral study of SS 433 is performed using more than 150 hours of observations with VERITAS, sensitive to $γ$-ray energies $>100\,\mathrm{GeV}$. With an angular resolution better than $0.1^°$, extended TeV $γ$-ray emission is resolved from both the eastern and western jet lobes, located tens of parsecs from the central binary. The emission appears elongated along the jet axis and coincides with regions where the jets interact with the surrounding supernova remnant. No TeV emission is detected from the central binary, nor is significant emission observed between the central binary and the jet lobes. Phase-resolved analyses show no evidence for variability with orbital or precessional phase, supporting a steady emission scenario. The observed morphology and spectra are consistent with scenarios where particles are accelerated in the lobes of the jets, possibly through shocks or alternative processes such as magnetic reconnection. The extended TeV emission from the jet lobes of SS 433 favors a leptonic origin in the VERITAS energy range, suggesting any hadronic acceleration is subdominant. The results offer valuable constraints on how microquasar jets may contribute to the Galactic cosmic-ray population toward the knee. △ Less

Submitted 23 March, 2026; originally announced March 2026.

Comments: 18 pages, 12 figures, Accepted for publication in the Astrophysical Journal (ApJ)

VERITAS Observations Contemporaneous with the LHAASO Detection of NGC 4278

Authors: The VERITAS Collaboration, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, J. H. Buckley, Y. Chen, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, M. Escobar Godoy, S. Feldman, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, Z. Hughes, T. B. Humensky, M. Iskakova, W. Jin , et al. (44 additional authors not shown)

Abstract: Significant gamma-ray emission between 1 TeV and 20 TeV from a point source, 1LHAASO J1219+2915, consistent with the location of the LINER/LLAGN galaxy NGC 4278 was recently reported by the LHAASO collaboration. These data were later split into active and quasi-quiet states, with most of the LHAASO significance coming from the active state (MJD 59449-59589). Subsequent analysis of Fermi-LAT and Sw… ▽ More Significant gamma-ray emission between 1 TeV and 20 TeV from a point source, 1LHAASO J1219+2915, consistent with the location of the LINER/LLAGN galaxy NGC 4278 was recently reported by the LHAASO collaboration. These data were later split into active and quasi-quiet states, with most of the LHAASO significance coming from the active state (MJD 59449-59589). Subsequent analysis of Fermi-LAT and Swift-XRT observations have been used to explore the double-peaked broad-band emission. Models of the spectral energy distribution (SED) are currently unconstrained due to the lack of contemporaneous multi-wavelength data at either peak. Here we report serendipitous observations of NGC 4278 with VERITAS, made possible by the contemporaneous observations of the nearby blazars 1ES 1218+304, 1ES 1215+303, and W Comae, each of which are located within $2^\circ$ of NGC 4278. VERITAS did not detect any gamma-ray emission and a flux upper limit was calculated. The flux upper limits constrain the photon spectrum of the quasi-quiet period, and together with Fermi-LAT, indicate a peak in the SED between 100 GeV and 2 TeV. We present an interpretation of the broadband SED that is based on acceleration of protons in the corona of the AGN, followed by p-$γ$ interactions and optically thin $γ$-ray emission. Within this framework, the implied neutrino signal is slightly below the current sensitivity of IceCube. △ Less

Submitted 18 March, 2026; originally announced March 2026.

Comments: 11 pages, 2 figures, Accepted for publication in the Astrophysical Journal (ApJ)

VERITAS contributions to the 39th International Cosmic Ray Conference

Authors: A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, Y. Chen, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Escobar Godoy, J. Escudero Pedrosa, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, T. B. Humensky, M. Iskakova, W. Jin, M. N. Johnson, E. Joshi, M. Kertzman, M. Kherlakian , et al. (33 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 39th International Cosmic Ray Conference (ICRC), held July 14 through July 24, 2025 in Geneva, Switzerland. Compilation of papers presented by the VERITAS Collaboration at the 39th International Cosmic Ray Conference (ICRC), held July 14 through July 24, 2025 in Geneva, Switzerland. △ Less

Submitted 6 January, 2026; v1 submitted 22 December, 2025; originally announced December 2025.

Prompt Searches for Very-High-Energy γ-Ray Counterparts to IceCube Astrophysical Neutrino Alerts

Authors: J. Abhir, A. Biland, K. Brand, T. Bretz, D. Dorner, L. Eisenberger, D. Elsaesser, P. Günther, S. Hasan, D. Hildebrand, K. Mannheim, M. Linhoff, F. Pfeifle, W. Rhode, B. Schleicher, V. Sliusar, M. Vorbrugg, R. Walter, F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin , et al. (809 additional authors not shown)

Abstract: The search for sources of high-energy astrophysical neutrinos can be significantly advanced through a multi-messenger approach, which seeks to detect the gamma rays that accompany neutrinos as they are produced at their sources. Multi-messenger observations have so far provided the first evidence for a neutrino source, illustrated by the joint detection of the flaring blazar TXS 0506+056 in highen… ▽ More The search for sources of high-energy astrophysical neutrinos can be significantly advanced through a multi-messenger approach, which seeks to detect the gamma rays that accompany neutrinos as they are produced at their sources. Multi-messenger observations have so far provided the first evidence for a neutrino source, illustrated by the joint detection of the flaring blazar TXS 0506+056 in highenergy (HE, E > 1 GeV) and very-high-energy (VHE, E > 100 GeV) gamma rays in coincidence with the high-energy neutrino IceCube-170922A, identified by IceCube. Imaging atmospheric Cherenkov telescopes (IACTs), namely FACT, H.E.S.S., MAGIC, and VERITAS, continue to conduct extensive neutrino target-of-opportunity follow-up programs. These programs have two components: followup observations of single astrophysical neutrino candidate events (such as IceCube-170922A), and observation of known gamma-ray sources after the identification of a cluster of neutrino events by IceCube. Here we present a comprehensive analysis of follow-up observations of high-energy neutrino events observed by the four IACTs between September 2017 (after the IceCube-170922A event) and January 2021. Our study found no associations between gamma-ray sources and the observed neutrino events. We provide a detailed overview of each neutrino event and its potential counterparts. Furthermore, a joint analysis of all IACT data is included, yielding combined upper limits on the VHE gamma-ray flux. △ Less

Submitted 18 December, 2025; originally announced December 2025.

Comments: accepted for publication in ApJ

Probing a cosmogenic origin of astrophysical neutrinos and cosmic rays using gamma-ray observations of TXS 0506+056

Authors: A. Acharyya, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, J. H. Buckley, Y. Chen, J. L. Christiansen, A. Duerr, M. Errando, M. Escobar Godoy, A. Falcone, S. Feldman, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, Z. Hughes, M. Iskakova, W. Jin, P. Kaaret , et al. (36 additional authors not shown)

Abstract: In September 2017, a high-energy neutrino event detected by the IceCube Neutrino Observatory (IceCube-170922A) was associated, at the $3σ$ level, with a gamma-ray flare from the blazar TXS 0506+056. Cosmic rays that are accelerated in astrophysical sources can escape from their jets and interact with background radiation fields. Interactions with the extragalactic background light can produce pion… ▽ More In September 2017, a high-energy neutrino event detected by the IceCube Neutrino Observatory (IceCube-170922A) was associated, at the $3σ$ level, with a gamma-ray flare from the blazar TXS 0506+056. Cosmic rays that are accelerated in astrophysical sources can escape from their jets and interact with background radiation fields. Interactions with the extragalactic background light can produce pions and hence neutrinos, while interactions with the cosmic microwave background predominantly drive inverse Compton scattering, contributing to electromagnetic cascades in intergalactic space. The resulting secondary gamma-ray emission can be detected with high-energy gamma-ray telescopes. Here, we report on a new search for such cosmogenic cascade emission from the blazar TXS 0506+056, using a combined data set from the Fermi-Large Area Telescope and VERITAS. We compare the gamma-ray spectrum and neutrino observations with the predictions of cosmic-ray induced cascades in intergalactic space. The observed gamma-ray spectrum is modeled as a combination of the primary spectrum and the cascade spectrum. We apply a Monte Carlo simulation with a $Δχ^2$-based likelihood analysis to jointly determine the best-fit parameters of a proton emission spectrum describing the data and derive constraints on the proton escape luminosity. Assuming a log-parabola primary photon spectrum, we find consistency with a proton injection spectral index of $α_{p} \simeq 2.0$ and a cutoff energy of $E_{p,\text{max}} \simeq 1.3 \times 10^{16}$ eV, and constrain the isotropic proton escape luminosity to $1 \times 10^{44}$ erg s$^{-1}$ $\lesssim L_{p, esc} \lesssim 3 \times 10^{45}$ erg s$^{-1}$ at the 90 % confidence level. △ Less

Submitted 8 November, 2025; originally announced November 2025.

Comments: 14 pages, 5 figures, Accepted for publication in the Astrophysical Journal (ApJ)

Constraints on Axion-Like Particles from VERITAS Observations of a Flaring Radio Galaxy in the Perseus Cluster

Authors: C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, Y. Chen, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, M. Escobar Godoy, J. Escudero Pedrosa, S. Feldman, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, Z. Hughes, T. B. Humensky, M. Iskakova, W. Jin , et al. (40 additional authors not shown)

Abstract: Background: Axion-like particles (ALPs) are hypothetical particles that emerge in numerous theoretical extensions to the Standard Model. Their coupling to electromagnetic field implies that ALPs would mix with photons in the presence of external magnetic fields. As ALP phenomenology is governed by the mass and strength of its coupling, there is a subset of this parameter space in which this mixing… ▽ More Background: Axion-like particles (ALPs) are hypothetical particles that emerge in numerous theoretical extensions to the Standard Model. Their coupling to electromagnetic field implies that ALPs would mix with photons in the presence of external magnetic fields. As ALP phenomenology is governed by the mass and strength of its coupling, there is a subset of this parameter space in which this mixing would be expected to leave an imprint on the spectra of TeV gamma-ray sources. Data: In 2017, the VERITAS gamma-ray observatory recorded the second day of a dramatic flare of the radio galaxy NGC 1275, embedded at the center of the Perseus galaxy cluster. This serendipitous locale provides a spatially-extended magnetic field of strength O(10$μ$G) through which escaping photons traverse, making it an excellent target to study ALPs. Methods: We analyze the VERITAS data of NGC 1275's 2017 flare with the gammapy analysis package. Extensive fitting and modeling are performed to ultimately conduct a likelihood analysis used to search for any evidence of a preference for ALPs and to explore the confidence with which constraints can be set. We adopt the CLs method for this study for its conservative approach to setting limits in regimes where the search has limited sensitivity. Results: No evidence for the existence of ALPs is found, and no combination of mass and coupling strength can be excluded at or above 95% confidence level. We provide a map showing the strength of our exclusions in the mass and coupling parameter space. The strongest exclusions are found in the mass range $2 \times 10^{-7}$eV $\lesssim m_a \lesssim 4 \times 10^{-7}$eV and at the coupling strength of $g_{aγ} \gtrsim 3 \times 10^{-11}$ GeV$^{-1}$ up to 80% confidence level, which are consistent with previous studies. Conclusions: We find the CLs method to be a trustworthy approach, and advocate for its... △ Less

Submitted 21 October, 2025; originally announced October 2025.

Comments: 18 pages, 6 figures, accepted for publication in the Physical Review D (PRD)

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

Authors: Fermi-LAT Collaboration, :, S. Abdollahi, L. Baldini, R. Bellazzini, B. Berenji, E. Bissaldi, R. Bonino, P. Bruel, S. Buson, E. Charles, A. W. Chen, S. Ciprini, M. Crnogorcevic, A. Cuoco, F. D'Ammando, A. de Angelis, M. Di Mauro, N. Di Lalla, L. Di Venere, A. Domínguez, S. J. Fegan, A. Fiori, P. Fusco, V. Gammaldi , et al. (582 additional authors not shown)

Abstract: Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect dark matter (DM) searches using gamma-ray telescopes because they are thought to have high DM content and a low astrophysical background. The sensitivity of these searches is improved by combining the observations of dSphs made by different gamma-ray telescopes. We present the results of a combined search by the most sensitive cu… ▽ More Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect dark matter (DM) searches using gamma-ray telescopes because they are thought to have high DM content and a low astrophysical background. The sensitivity of these searches is improved by combining the observations of dSphs made by different gamma-ray telescopes. We present the results of a combined search by the most sensitive currently operating gamma-ray telescopes, namely: the satellite-borne Fermi-LAT telescope; the ground-based imaging atmospheric Cherenkov telescope arrays H.E.S.S., MAGIC, and VERITAS; and the HAWC water Cherenkov detector. Individual datasets were analyzed using a common statistical approach. Results were subsequently combined via a global joint likelihood analysis. We obtain constraints on the velocity-weighted cross section $\langle σ\mathit{v} \rangle$ for DM self-annihilation as a function of the DM particle mass. This five-instrument combination allows the derivation of up to 2-3 times more constraining upper limits on $\langle σ\mathit{v} \rangle$ than the individual results over a wide mass range spanning from 5 GeV to 100 TeV. Depending on the DM content modeling, the 95% confidence level observed limits reach $1.5\times$10$^{-24}$ cm$^3$s$^{-1}$ and $3.2\times$10$^{-25}$ cm$^3$s$^{-1}$, respectively, in the $τ^+τ^-$ annihilation channel for a DM mass of 2 TeV. △ Less

Submitted 27 August, 2025; originally announced August 2025.

Journal ref: JCAP 03 (2026) 035

Multi-wavelength Study of HESS J0632+057: New Insights into Pulsar-Disk Interaction

Authors: Jaegeun Park, Hongjun An, Chanho Kim, Natalie Matchett, Kaya Mori, Brian van Soelen, VERITAS Collaboration, :, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, J. H. Buckley, Y. Chen, A. J. Chromey, A. Duerr, M. Errando, M. Escobar Godoy, A. Falcone, S. Feldman, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon , et al. (38 additional authors not shown)

Abstract: We present an analysis of new multi-wavelength observations of the TeV gamma-ray binary HESS J0632+057, conducted using SALT, Swift, NuSTAR, and VERITAS in 2023--2024. By combining these new data with archival observations, we confirm previous suggestions of orbital variability in the source's X-ray spectrum, including increased X-ray absorption at the orbital phase interval of… ▽ More We present an analysis of new multi-wavelength observations of the TeV gamma-ray binary HESS J0632+057, conducted using SALT, Swift, NuSTAR, and VERITAS in 2023--2024. By combining these new data with archival observations, we confirm previous suggestions of orbital variability in the source's X-ray spectrum, including increased X-ray absorption at the orbital phase interval of $φ\approx0.3\textrm{--}0.4$. The source's X-ray flux within this phase interval seems to have exhibited a significant change on an orbital timescale. Additionally, occasional short-term variations in the X-ray band on a timescale of less than 3 days have been observed. The measured duration of the increased absorbing column density and the flux variability timescales can provide clues about the interaction between the putative pulsar and the Be companion's disk if, as previously suggested, the pulsar crosses the disk at this phase interval. Moreover, the new contemporaneous X-ray and TeV observations around the pulsar-crossing phases revealed independent variability in the X-ray and TeV fluxes, contrary to a previous observation of concurrent flux increases. While these observations alone cannot provide definitive conclusions, we discuss our results in the context of pulsar-disk interaction and intrabinary shock emission scenarios. △ Less

Submitted 31 July, 2025; originally announced July 2025.

Comments: 15 pages, 7 figures. Accepted for publication in ApJ

Measurement of the photosphere oblateness of $γ$ Cassiopeiae via Stellar Intensity Interferometry with the VERITAS Observatory

Authors: A. Archer, J. P. Aufdenberg, P. Bangale, J. T. Bartkoske, W. Benbow, J. H. Buckley, Y. Chen, N. B. Y. Chin, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Escobar Godoy, S. Feldman, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, Z. Hughes, T. B. Humensky, W. Jin, M. N. Johnson , et al. (36 additional authors not shown)

Abstract: We use the stellar intensity interferometry system implemented with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) at Fred Lawrence Whipple Observatory (FLWO) as a light collector to obtain measurements of the rapid rotator star $γ$ Cassiopeiae, at a wavelength of 416 nm. Using data from baselines sampling different position angles, we extract the size, oblateness, and proje… ▽ More We use the stellar intensity interferometry system implemented with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) at Fred Lawrence Whipple Observatory (FLWO) as a light collector to obtain measurements of the rapid rotator star $γ$ Cassiopeiae, at a wavelength of 416 nm. Using data from baselines sampling different position angles, we extract the size, oblateness, and projected orientation of the photosphere. Fitting the data with a uniform ellipse model yields a minor-axis angular diameter of $0.43\pm0.02$ mas, a major-to-minor-radius ratio of $1.28\pm0.04$, and a position angle of $116^\circ\pm5^\circ$ for the axis of rotation. A rapidly-rotating stellar atmosphere model that includes limb and gravity darkening describes the data well with a fitted angular diameter of $0.604^{+0.041}_{-0.034}$ mas corresponding to an equatorial radius of 10.9$^{+0.8}_{-0.6}~R_\odot$, a rotational velocity with a $1~σ$ lower limit at $97.7\%$ that of breakup velocity, and a position angle of $114.7^{+6.4}_{-5.7}$ degrees. These parameters are consistent with H$α$ line spectroscopy and infrared-wavelength Michelson interferometric measurements of the star's decretion disk. This is the first measurement of an oblate photosphere using intensity interferometry. △ Less

Submitted 17 June, 2025; originally announced June 2025.

Multiwavelength observation of a candidate pulsar halo LHAASO J0621+3755 and the first X-ray detection of PSR J0622+3749

Authors: C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, J. H. Buckley, Y. Chen, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, M. Escobar Godoy, A. Falcone, S. Feldman, Q. Feng, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, T. B. Humensky, W. Jin, M. N. Johnson, P. Kaaret , et al. (49 additional authors not shown)

Abstract: Pulsar halos are regions around middle-aged pulsars extending out to tens of parsecs. The large extent of the halos and well-defined central cosmic-ray accelerators make this new class of Galactic sources an ideal laboratory for studying cosmic-ray transport. LHAASO J0621+3755 is a candidate pulsar halo associated with the middle-aged gamma-ray pulsar PSR J0622+3749. We observed LHAASO J0621+3755… ▽ More Pulsar halos are regions around middle-aged pulsars extending out to tens of parsecs. The large extent of the halos and well-defined central cosmic-ray accelerators make this new class of Galactic sources an ideal laboratory for studying cosmic-ray transport. LHAASO J0621+3755 is a candidate pulsar halo associated with the middle-aged gamma-ray pulsar PSR J0622+3749. We observed LHAASO J0621+3755 with VERITAS and XMM-Newton in the TeV and X-ray bands, respectively. For this work, we developed a novel background estimation technique for imaging atmospheric Cherenkov telescope observations of such extended sources. No halo emission was detected with VERITAS (0.3--10 TeV) or XMM-Newton (2--7 keV) within 1 degree and 10 arcmin around PSR J0622+3749, respectively. Combined with the LHAASO-KM2A and Fermi-LAT data, VERITAS flux upper limits establish a spectral break at ~1--10 TeV, a unique feature compared with Geminga, the most studied pulsar halo. We model the gamma-ray spectrum and LHAASO-KM2A surface brightness as inverse Compton emission and find suppressed diffusion around the pulsar, similar to Geminga. A smaller diffusion suppression zone and harder electron injection spectrum than Geminga are necessary to reproduce the spectral cutoff. A magnetic field <= 1 uG is required by our XMM-Newton observation and synchrotron spectral modeling, consistent with Geminga. Our findings support slower diffusion and lower magnetic field around pulsar halos than the Galactic averages, hinting at magnetohydrodynamic turbulence around pulsars. Additionally, we report the detection of an X-ray point source spatially coincident with PSR J0622+3749, whose periodicity is consistent with the gamma-ray spin period of 333.2 ms. The soft spectrum of this source suggests a thermal origin. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: Accepted for publication in ApJ. 28 pages, 19 figures

Constraints on the X-ray and Very High Energy $γ$-ray Flux from Supernova Remnant W44

Authors: A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, J. H. Buckley, Y. Chen, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, M. Escobar Godoy, S. Feldman, Q. Feng, J. Foote, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, T. B. Humensky, W. Jin, M. N. Johnson, P. Kaaret, M. Kertzman , et al. (33 additional authors not shown)

Abstract: Observations of GeV gamma-ray emission from the well-studied mixed-morphology supernova remnant (SNR) W44 by Fermi-LAT and AGILE imply that it is a site of significant cosmic ray acceleration. The spectral energy distribution (SED) derived from the GeV data suggest that the gamma-ray emission likely originates from the decay of neutral pions generated by cosmic-ray interactions. It is essential to… ▽ More Observations of GeV gamma-ray emission from the well-studied mixed-morphology supernova remnant (SNR) W44 by Fermi-LAT and AGILE imply that it is a site of significant cosmic ray acceleration. The spectral energy distribution (SED) derived from the GeV data suggest that the gamma-ray emission likely originates from the decay of neutral pions generated by cosmic-ray interactions. It is essential to measure the SED of W44 in the X-ray and very high energy (VHE) gamma-ray bands to verify the hadronic origin of the emission and to gauge the potential contributions from leptonic emission. We report an upper-limit of the nonthermal X-ray flux from W44 of 5 $\times$ 10$^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 0.5 - 8.0 keV band based on $\sim$ 300 ks of XMM-Newton observations. The X-ray upper limit is consistent with previously estimated hadronic models, but in tension with the leptonic models. We estimate the VHE flux upper limit of $\sim$ 1.2 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$ in the 0.5 - 5.0 TeV range from W44 using data from the Very Energetic Radiation Imaging Telescope Array System (VERITAS). Our non-detection of W44 at VHE wavlengths is in agreemnent with observations from other imaging atmospheric Cherenkov telescopes (IACTs) and is perhaps consistent with the evolutionary stage of the SNR. △ Less

Submitted 12 March, 2025; originally announced March 2025.

Comments: 9 pages, 3 figures, 1 table, accepted for publication in ApJ

A multi-wavelength study to decipher the 2017 flare of the blazar OJ 287

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. Brill, J. P. Caldwell, M. Carini, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, Q. Feng, J. P. Finley, J. Foote, L. Fortson, A. Furniss, G. Gallagher, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Hoang , et al. (49 additional authors not shown)

Abstract: In February 2017, the blazar OJ~287 underwent a period of intense multiwavelength activity. It reached a new historic peak in the soft X-ray (0.3-10 keV) band, as measured by Swift-XRT. This event coincides with a very-high-energy (VHE) $γ$-ray outburst that led VERITAS to detect emission above 100 GeV, with a detection significance of $10σ$ (from 2016 December 9 to 2017 March 31). The time-averag… ▽ More In February 2017, the blazar OJ~287 underwent a period of intense multiwavelength activity. It reached a new historic peak in the soft X-ray (0.3-10 keV) band, as measured by Swift-XRT. This event coincides with a very-high-energy (VHE) $γ$-ray outburst that led VERITAS to detect emission above 100 GeV, with a detection significance of $10σ$ (from 2016 December 9 to 2017 March 31). The time-averaged VHE $γ$-ray spectrum was consistent with a soft power law ($Γ= -3.81 \pm 0.26$) and an integral flux corresponding to $\sim2.4\%$ that of the Crab Nebula above the same energy. Contemporaneous data from multiple instruments across the electromagnetic spectrum reveal complex flaring behavior, primarily in the soft X-ray and VHE bands. To investigate the possible origin of such an event, our study focuses on three distinct activity states: before, during, and after the February 2017 peak. The spectral energy distributions during these periods suggest the presence of at least two non-thermal emission zones, with the more compact one responsible for the observed flare. Broadband modeling results and observations of a new radio knot in the jet of OJ~287 in 2017 are consistent with a flare originating from a strong recollimation shock outside the radio core. △ Less

Submitted 26 August, 2024; v1 submitted 16 July, 2024; originally announced July 2024.

Comments: 21 pages, 7 figures, accepted for publication in The Astrophysical Journal

VERITAS contributions to the 38th International Cosmic Ray Conference

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, Z. Hughes, T. B. Humensky, W. Jin, M. N. Johnson, M. Kertzman , et al. (39 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: html page. ICRC 2023, Nagoya, Japan

Multiwavelength Observations of the Blazar PKS 0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. Brill, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, T. B. Humensky , et al. (185 additional authors not shown)

Abstract: We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ra… ▽ More We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution (SED), and the gamma-ray data from Fermi -LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both X-ray and gamma-ray measurements provide strong constraints on the leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed SED. The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models and allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. We show a numerical lepto-hadronic model with external target photons that reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 21 pages, 3 figures, accepted by ApJ

A VERITAS/Breakthrough Listen Search for Optical Technosignatures

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Pedro Batista, Wystan Benbow, Aryeh Brill, M Capasso, Manel Errando, Abraham Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Sean Griffin, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, T. Humensky, Weidong Jin, Philip Kaaret , et al. (43 additional authors not shown)

Abstract: The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectabl… ▽ More The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectable over interstellar distances. We report here on the analysis and results of dedicated VERITAS observations of Breakthrough Listen targets conducted in 2019 and 2020 and of archival VERITAS data collected since 2012. Thirty hours of dedicated observations of 136 targets and 249 archival observations of 140 targets were analyzed and did not reveal any signals consistent with a technosignature. The results are used to place limits on the fraction of stars hosting transmitting civilizations. We also discuss the minimum-pulse sensitivity of our observations and present VERITAS observations of CALIOP: a space-based pulsed laser onboard the CALIPSO satellite. The detection of these pulses with VERITAS, using the analysis techniques developed for our technosignature search, allows a test of our analysis efficiency and serves as an important proof-of-principle. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 15 pages, 7 figures

VERITAS discovery of very high energy gamma-ray emission from S3 1227+25 and multiwavelength observations

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Wystan Benbow, Aryeh Brill, Jodi Christiansen, Alisha Chromey, Manel Errando, Abe Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Greg Gallagher, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, Weidong Jin, Madalyn Johnson, Philip Kaaret , et al. (46 additional authors not shown)

Abstract: We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th… ▽ More We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th resulted in a strong 13$σ$ detection with a differential photon spectral index, $Γ$ = 3.8 $\pm$ 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target of opportunity observations with Swift, optical photometry, polarimetry and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of $τ_{obs}$ = 6.2 $\pm$ 0.9 hours, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cut-off. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 18 pages, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

Search for Ultraheavy Dark Matter from Observations of Dwarf Spheroidal Galaxies with VERITAS

Authors: A. Acharyya, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, M. Baumgart, W. Benbow, J. H. Buckley, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. F. Hanlon, O. Hervet, J. Hoang, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (29 additional authors not shown)

Abstract: Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has bee… ▽ More Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has been suggested as an under-explored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS gamma-ray observatory. With 216 hours of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a $γ$-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle. △ Less

Submitted 24 February, 2023; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 10 pages, 7 figures

VTSCat: The VERITAS Catalog of Gamma-Ray Observations

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. H. Buckley, A. Brill, M. Capasso, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, K. A Farrell, Q. Feng, J. P. Finley, G. M Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, C. Giuri, O. Gueta , et al. (64 additional authors not shown)

Abstract: The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four… ▽ More The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hours 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. △ Less

Submitted 13 January, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: Version with corrected author list

Journal ref: Research Notes of the AAS, Volume 7, Number 1, 2023

VERITAS contributions to the 37th International Cosmic Ray Conference

Authors: C. B. Adams, A. Archer, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, K. A. Farrell, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, O. Hervet, J. Holder, B. Hona, T. B. Humensky, W. Jin , et al. (36 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 37th International Cosmic Ray Conference (ICRC), held July 12 through July 23, 2021 (online) in Berlin, Germany. Compilation of papers presented by the VERITAS Collaboration at the 37th International Cosmic Ray Conference (ICRC), held July 12 through July 23, 2021 (online) in Berlin, Germany. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: html page. 2021 ICRC, Berlin, Germany

VERITAS Discovery of VHE Emission from the Radio Galaxy 3C 264: A Multi-Wavelength Study

Authors: A. Archer, W. Benbow, R. Bird, A. Brill, M. Buchovecky, J. H. Buckley, M. T. Carini, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, Q. Feng, P. Fortin, L. Fortson, A. Furniss, A. Gent, M. Georganopoulos, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder , et al. (45 additional authors not shown)

Abstract: The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $γ$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a stat… ▽ More The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $γ$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a statistical significance of 7.8 standard deviations above background. The observed VHE flux is variable on monthly time scales, with an elevated flux seen in 2018 observations. The VHE emission during this elevated state is well-characterized by a power-law spectrum with a photon index $Γ= 2.20 \pm 0.27$ and flux F($>315$ GeV) = ($7.6\pm 1.2_{\mathrm stat} \pm 2.3_{\mathrm syst})\times 10^{-13}$ cm$^{-2}$ s$^{-1}$, or approximately 0.7% of the Crab Nebula flux above the same threshold. 3C 264 ($z = 0.0217$) is the most distant radio galaxy detected at VHE, and the elevated state is thought to be similar to that of the famously outbursting jet in M 87. Consequently, extensive contemporaneous multi-wavelength data were acquired in 2018 at the time of the VHE high state. An analysis of these data, including VLBA, VLA, HST, Chandra and Swift observations in addition to the VERITAS data, is presented, along with a discussion of the resulting spectral energy distribution. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: 19 pages, 11 figures, Accepted for publication in Astrophysical Journal

Evidence for proton acceleration up to TeV energies based on VERITAS and Fermi-LAT observations of the Cas A SNR

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, A. J. Chromey, W. Cui, M. K. Daniel, S. Das, V. V. Dwarkadas, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (38 additional authors not shown)

Abstract: We present a study of $γ$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the… ▽ More We present a study of $γ$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from \textit{Fermi}-LAT and VERITAS deviates significantly from a simple power-law, and is best described by a power-law with spectral index of $2.17\pm 0.02_{stat}$ with a cut-off energy of $2.3 \pm 0.5_{stat}$ TeV. These results, along with radio, X-ray and $γ$-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed $γ$-ray spectrum. From modeling of the entire multi-wavelength spectrum, a minimum magnetic field inside the remnant of $B_{\mathrm{min}}\approx150\,\mathrm{μG}$ is deduced. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 33 pages, 9 Figures, 6 Tables

A decade of multi-wavelength observations of the TeV blazar 1ES 1215+303: Extreme shift of the synchrotron peak frequency and long-term optical-gamma-ray flux increase

Authors: Janeth Valverde, Deirdre Horan, Denis Bernard, Stephen Fegan, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. H. Buckley, J. L. Christiansen, W. Cui, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan , et al. (64 additional authors not shown)

Abstract: Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radi… ▽ More Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed. △ Less

Submitted 12 February, 2020; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: ApJ 891 (2020) 170

Probing the Properties of the Pulsar Wind in the Gamma-Ray Binary HESS J0632+057 with NuSTAR and VERITAS Observations

Authors: A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. L. Christiansen, A. J. Chromey, W. Cui, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, P. Kaaret , et al. (38 additional authors not shown)

Abstract: HESS J0632+057 is a gamma-ray binary composed of a compact object orbiting a Be star with a period of about $315$ days. Extensive X-ray and TeV gamma-ray observations have revealed a peculiar light curve containing two peaks, separated by a dip. We present the results of simultaneous observations in hard X-rays with NuSTAR and in TeV gamma-rays with VERITAS, performed in November and December 2017… ▽ More HESS J0632+057 is a gamma-ray binary composed of a compact object orbiting a Be star with a period of about $315$ days. Extensive X-ray and TeV gamma-ray observations have revealed a peculiar light curve containing two peaks, separated by a dip. We present the results of simultaneous observations in hard X-rays with NuSTAR and in TeV gamma-rays with VERITAS, performed in November and December 2017. These observations correspond to the orbital phases $φ\approx0.22$ and $0.3$, where the fluxes are rising towards the first light-curve peak. A significant variation of the spectral index from 1.77$\pm$0.05 to 1.56$\pm$0.05 is observed in the X-ray data. The multi-wavelength spectral energy distributions (SED) derived from the observations are interpreted in terms of a leptonic model, in which the compact object is assumed to be a pulsar and non-thermal radiation is emitted by high-energy electrons accelerated at the shock formed by the collision between the stellar and pulsar wind. The results of the SED fitting show that our data can be consistently described within this scenario, and allow us to estimate the magnetization of the pulsar wind at the location of the shock formation. The constraints on the pulsar-wind magnetization provided by our results are shown to be consistent with those obtained from other systems. △ Less

Submitted 12 December, 2019; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: Accepted for publication in The Astrophysical Journal

Measurement of the extragalactic background light spectral energy distribution with VERITAS

Authors: VERITAS collaboration, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. L. Christiansen, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, M. Fernandez-Alonso, J. P. Finley, L. Fortson, A. Furniss, A. Gent, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (37 additional authors not shown)

Abstract: The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the Universe's history. Spectral measurements of blazars at very high energies ($>$100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the… ▽ More The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the Universe's history. Spectral measurements of blazars at very high energies ($>$100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the gamma-ray photons with the EBL. The spectra of 14 VERITAS-detected blazars are included in a new measurement of the EBL SED that is independent of EBL SED models. The resulting SED covers an EBL wavelength range of 0.56--56 $μ$m, and is in good agreement with lower limits obtained by assuming that the EBL is entirely due to radiation from cataloged galaxies. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: Accepted for publication in The Astrophysical Journal

VERITAS contributions to the 36th International Cosmic Ray Conference

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. K. Daniel, S. Das, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (43 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin. Compilation of papers presented by the VERITAS Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin. △ Less

Submitted 23 September, 2019; v1 submitted 17 September, 2019; originally announced September 2019.

Comments: html page. 2019 ICRC, Madison, Wisconsin

A Search for Pulsed Very High-Energy Gamma Rays from Thirteen Young Pulsars in Archival VERITAS Data

Authors: A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, A. J. Chromey, W. Cui, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret, P. Kar, N. Kelley-Hoskins , et al. (36 additional authors not shown)

Abstract: We conduct a search for periodic emission in the very high-energy gamma-ray band (VHE; E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hours. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in non-detections of pulsed VHE gamma rays from each pulsar. Up… ▽ More We conduct a search for periodic emission in the very high-energy gamma-ray band (VHE; E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hours. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in non-detections of pulsed VHE gamma rays from each pulsar. Upper limits on a potential VHE gamma-ray flux are derived at the 95% confidence level above three energy thresholds using two methods. These are the first such searches for pulsed VHE emission from each of the pulsars, and the obtained limits constrain a possible flux component manifesting at VHEs as is seen for the Crab pulsar. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Comments: 32 pages, 20 figures, accepted for publication in ApJ

VERITAS and Fermi-LAT observations of new HAWC sources

Authors: VERITAS Collaboration, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, V. Bugaev, A. J. Chromey, M. P. Connolly, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, M. Hutten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson , et al. (259 additional authors not shown)

Abstract: The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detect… ▽ More The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: Accepted for publication in the ApJ, Corresponding author: Nahee Park (VERITAS Collaboration), John W. Hewitt (Fermi-LAT Collaboration), Ignacio Taboada (HAWC Collaboration), 30 pages, 12 figures

Journal ref: ApJ 866 (2018) no.1, 24

Measurement of Cosmic-ray Electrons at TeV Energies by VERITAS

Authors: VERITAS Collaboration, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, V. Bugaev, M. P. Connolly, W. Cui, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, G. Gillanders, M. Hütten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret, P. Kar, N. Kelley-Hoskins , et al. (36 additional authors not shown)

Abstract: Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering processes while propagating within the Galaxy and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance there a… ▽ More Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering processes while propagating within the Galaxy and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance there are only a few known astrophysical objects capable of accelerating electrons to such high energies. It is also possible that the CREs are the products of the annihilation or decay of heavy dark matter (DM) particles. VERITAS, an array of imaging air Cherenkov telescopes in southern Arizona, USA, is primarily utilized for gamma-ray astronomy, but also simultaneously collects CREs during all observations. We describe our methods of identifying CREs in VERITAS data and present an energy spectrum, extending from 300 GeV to 5 TeV, obtained from approximately 300 hours of observations. A single power-law fit is ruled out in VERITAS data. We find that the spectrum of CREs is consistent with a broken power law, with a break energy at 710 $\pm$ 40$_{stat}$ $\pm$ 140$_{syst}$ GeV. △ Less

Submitted 29 August, 2018; originally announced August 2018.

Comments: 17 pages, 2 figures, accepted for publication in PRD

The extreme HBL behaviour of Markarian 501 during 2012

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras, J. Cortina, S. Covino , et al. (254 additional authors not shown)

Abstract: A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0… ▽ More A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays. △ Less

Submitted 14 August, 2018; v1 submitted 13 August, 2018; originally announced August 2018.

Comments: 25 pages, 9 figures, accepted for publication in A&A. Corresponding authors: Gareth Hughes (gareth.hughes@cfa.harvard.edu), David Paneque (dpaneque@mppmu.mpg.de), Amit Shukla (amit.shukla@astro.uni-wuerzburg.de)

Journal ref: A&A 620, A181 (2018)

Measurement of the Iron Spectrum in Cosmic Rays by VERITAS

Authors: The VERITAS collaboration, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, V. Bugaev, M. P. Connolly, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, M. Hütten, C. A. Johnson, P. Kaaret, N. Kelley-Hoskins , et al. (33 additional authors not shown)

Abstract: We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 to 500 TeV. The measurement makes use of a template-based analysis method, which, for the first time, is applied to the energy reconstruction of iron-induced air showers recorded by the VERITAS array of imaging atmospheric Cherenkov telescopes. The event selection makes use of the direct Cherenkov light which… ▽ More We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 to 500 TeV. The measurement makes use of a template-based analysis method, which, for the first time, is applied to the energy reconstruction of iron-induced air showers recorded by the VERITAS array of imaging atmospheric Cherenkov telescopes. The event selection makes use of the direct Cherenkov light which is emitted by charged particles before the first interaction, as well as other parameters related to the shape of the recorded air shower images. The measured spectrum is well described by a power law $\frac{\mathrm{d} F}{\mathrm{d} E}=f_0\cdot \left(\frac{E}{E_0}\right)^{-γ}$ over the full energy range, with $γ= 2.82 \pm 0.30 \mathrm{(stat.)} ^{+0.24}_{-0.27} \mathrm{(syst.)}$ and $f_0 = \left( 4.82 \pm 0.98 \mathrm{(stat.)}^{+2.12}_{-2.70} \mathrm{(syst.)} \right)\cdot 10^{-7}$m$^{-2}$s$^{-1}$sr$^{-1}$TeV$^{-1}$ at $E_0=50$TeV, with no indication of a cutoff or spectral break. The measured differential flux is compatible with previous results, with improved statistical uncertainty at the highest energies. △ Less

Submitted 20 July, 2018; originally announced July 2018.

Comments: 16 pages, 5 figures, submitted to Phys. Rev. D

Journal ref: Phys. Rev. D 98, 022009 (2018)

VERITAS observations of the BL Lac object TXS 0506+056

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, G. H. Gillanders, O. Gueta, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret, P. Kar , et al. (35 additional authors not shown)

Abstract: On 2017 September 22, the IceCube Neutrino Observatory reported the detection of the high-energy neutrino event \icnu, of potential astrophysical origin. It was soon determined that the neutrino direction was consistent with the location of the gamma-ray blazar \txs~(3FGL J0509.4+0541), which was in an elevated gamma-ray emission state as measured by the \emph{Fermi} satellite. VERITAS observation… ▽ More On 2017 September 22, the IceCube Neutrino Observatory reported the detection of the high-energy neutrino event \icnu, of potential astrophysical origin. It was soon determined that the neutrino direction was consistent with the location of the gamma-ray blazar \txs~(3FGL J0509.4+0541), which was in an elevated gamma-ray emission state as measured by the \emph{Fermi} satellite. VERITAS observations of the neutrino/blazar region started on 2017 September 23 in response to the neutrino alert and continued through 2018 February 6. While no significant very-high-energy (VHE; E $>$ 100 GeV) emission was observed from the blazar by VERITAS in the two-week period immediately following the IceCube alert, TXS 0506+056 was detected by VERITAS with a significance of 5.8 standard deviations ($σ$) in the full 35-hour data set. The average photon flux of the source during this period was $(8.9 \pm 1.6) \times 10^{-12} \; \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$, or 1.6\% of the Crab Nebula flux, above an energy threshold of 110 GeV, with a soft spectral index of $4.8 \pm 1.3$. △ Less

Submitted 12 July, 2018; originally announced July 2018.

Journal ref: ApJL 2018

HESS J1943+213: An Extreme Blazar Shining Through The Galactic Plane

Authors: The VERITAS Collaboration, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, V. Bugaev, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, A. Flinders, L. Fortson, A. Furniss, G. H. Gillanders, M. Hütten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret, P. Kar , et al. (38 additional authors not shown)

Abstract: HESS J1943+213 is a very-high-energy (VHE; $>$100 GeV) $γ$-ray source in the direction of the Galactic Plane. Studies exploring the classification of the source are converging towards its identification as an extreme synchrotron BL Lac object. Here we present 38 hours of VERITAS observations of HESS J1943+213 taken over two years. The source is detected with $\sim$20 standard deviations significan… ▽ More HESS J1943+213 is a very-high-energy (VHE; $>$100 GeV) $γ$-ray source in the direction of the Galactic Plane. Studies exploring the classification of the source are converging towards its identification as an extreme synchrotron BL Lac object. Here we present 38 hours of VERITAS observations of HESS J1943+213 taken over two years. The source is detected with $\sim$20 standard deviations significance, showing a remarkably stable flux and spectrum in VHE $γ$-rays. Multi-frequency very-long-baseline array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with European VLBI Network and detect this component in the 4.6 GHz and the 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, $Fermi$-LAT, $Swift$-XRT, FLWO 48$''$ telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron-self-Compton model. The well-measured $γ$-ray peak of the SED with VERITAS and $Fermi$-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary $γ$-rays from ultra-high-energy cosmic ray-initiated electromagnetic cascades to the $γ$-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and $γ$-ray bands. No statistically significant flux or spectral variability is detected. △ Less

Submitted 11 June, 2018; originally announced June 2018.

Comments: Accepted for publication in The Astrophysical Journal

A Very High Energy $γ$-Ray Survey towards the Cygnus Region of the Galaxy

Authors: The VERITAS Collaboration, A. U. Abeysekara, A. Archer, T. Aune, W. Benbow, R. Bird, R. Brose, M. Buchovecky, V. Bugaev, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, A. Flinders, L. Fortson, A. Furniss, E. V. Gotthelf, J. Grube, D. Hanna, O. Hervet, J. Holder, K. Huang, G. Hughes , et al. (46 additional authors not shown)

Abstract: We present results from deep observations towards the Cygnus region using 300 hours of very-high-energy (VHE) $γ$-ray data taken with the VERITAS Cherenkov telescope array and over seven years of high-energy $γ$-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse $γ$-ray emission in the northern sky, the Cygnus region provides a promising area to… ▽ More We present results from deep observations towards the Cygnus region using 300 hours of very-high-energy (VHE) $γ$-ray data taken with the VERITAS Cherenkov telescope array and over seven years of high-energy $γ$-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse $γ$-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130), and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL 2021.0+4031e (the Gamma-Cygni supernova remnant) was examined and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE $γ$-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula ux. We interpret these observations in a multiwavelength context and present the most detailed $γ$-ray view of the region to date. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: 49 Pages, 22 Figures, 13 Tables, Accepted to ApJ

A strong limit on the very-high-energy emission from GRB 150323A

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, V. Bugaev, M. P. Connolly, W. Cui, M. Errando, A. Falcone, Q. Feng, J. P. Finley, A. Flinders, L. Fortson, A. Furniss, G. H. Gillanders, M. Hütten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret , et al. (41 additional authors not shown)

Abstract: On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit i… ▽ More On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant since the very-high-energy (VHE) observation started only $\sim$2 minutes after the prompt emission peaked, and $Fermi$-LAT observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB~150323A ($z=0.593$) limits the attenuation by the extragalactic background light to $\sim 50$ % at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below $\sim100$ GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be $A\gtrsim 3\times 10^{11}$ g cm$^{-1}$, consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the ISM, which therefore cannot be ruled out as the environment of GRB 150323A. △ Less

Submitted 3 March, 2018; originally announced March 2018.

Comments: 6 pages, 1 figure. Accepted for publication in The Astrophysical Journal

VERITAS contributions to the 35th International Cosmic Ray Conference

Authors: A. U. Abeysekara, S. Archambault, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. L. Christiansen, M. P. Connolly, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, M. Fernandez-Alonso, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, G. H. Gillanders, M. Hütten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky , et al. (41 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 35th International Cosmic Ray Conference (ICRC), held July 12 through July 20, 2017 in Busan, South Korea. Compilation of papers presented by the VERITAS Collaboration at the 35th International Cosmic Ray Conference (ICRC), held July 12 through July 20, 2017 in Busan, South Korea. △ Less

Submitted 22 September, 2017; originally announced September 2017.

Comments: html page. 2017 ICRC, Busan, South Korea. A full list of author affiliations can be found at this link: https://veritas.sao.arizona.edu/conferences/ICRC2017_AuthList.html

Discovery of very-high-energy emission from RGB J2243+203 and derivation of its redshift upper limit

Authors: A. U. Abeysekara, S. Archambault, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, V. Bugaev, M. Cerruti, M. P. Connolly, W. Cui, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, G. H. Gillanders, S. Griffin, J. Grube, M. Hutten, D. Hanna, O. Hervet, J. Holder , et al. (45 additional authors not shown)

Abstract: Very-high-energy (VHE; $>$ 100 GeV) gamma-ray emission from the blazar RGB J2243+203 was discovered with the VERITAS Cherenkov telescope array, during the period between 21 and 24 December 2014. The VERITAS energy spectrum from this source can be fit by a power law with a photon index of $4.6 \pm 0.5$, and a flux normalization at 0.15 TeV of… ▽ More Very-high-energy (VHE; $>$ 100 GeV) gamma-ray emission from the blazar RGB J2243+203 was discovered with the VERITAS Cherenkov telescope array, during the period between 21 and 24 December 2014. The VERITAS energy spectrum from this source can be fit by a power law with a photon index of $4.6 \pm 0.5$, and a flux normalization at 0.15 TeV of $(6.3 \pm 1.1) \times 10^{-10} ~ \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{TeV}^{-1}$. The integrated \textit{Fermi}-LAT flux from 1 GeV to 100 GeV during the VERITAS detection is $(4.1 \pm 0.8) \times 10^{\textrm{-8}} ~\textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$, which is an order of magnitude larger than the four-year-averaged flux in the same energy range reported in the 3FGL catalog, ($4.0 \pm 0.1 \times 10^{\textrm{-9}} ~ \textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$). The detection with VERITAS triggered observations in the X-ray band with the \textit{Swift}-XRT. However, due to scheduling constraints \textit{Swift}-XRT observations were performed 67 hours after the VERITAS detection, not simultaneous with the VERITAS observations. The observed X-ray energy spectrum between 2 keV and 10 keV can be fitted with a power-law with a spectral index of $2.7 \pm 0.2$, and the integrated photon flux in the same energy band is $(3.6 \pm 0.6) \times 10^{-13} ~\textrm{cm}^{-2} \textrm{s}^{-1}$. EBL model-dependent upper limits of the blazar redshift have been derived. Depending on the EBL model used, the upper limit varies in the range from z $<~0.9$ to z $<~1.1$. △ Less

Submitted 15 September, 2017; originally announced September 2017.

Very-High-Energy $γ$-Ray Observations of the Blazar 1ES 2344+514 with VERITAS

Authors: C. Allen, S. Archambault, A. Archer, W. Benbow, R. Bird, E. Bourbeau, R. Brose, M. Buchovecky, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, X. Chen, J. L. Christiansen, M. P. Connolly, W. Cui, M. K. Daniel, J. D. Eisch, A. Falcone, Q. Feng, M. Fernandez-Alonso, J. P. Finley, H. Fleischhack, A. Flinders, L. Fortson , et al. (57 additional authors not shown)

Abstract: We present very-high-energy $γ$-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above background of $20.8σ$ in $47.2$ hours (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Usi… ▽ More We present very-high-energy $γ$-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above background of $20.8σ$ in $47.2$ hours (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Using these observations the temporal properties of 1ES 2344+514 are studied on short and long times scales. We fit a constant flux model to nightly- and seasonally-binned light curves and apply a fractional variability test, to determine the stability of the source on different timescales. We reject the constant-flux model for the 2007-2008 and 2014-2015 nightly-binned light curves and for the long-term seasonally-binned light curve at the $> 3σ$ level. The spectra of the time-averaged emission before and after correction for attenuation by the extragalactic background light are obtained. The observed time-averaged spectrum above 200 GeV is satisfactorily fitted (${χ^2/NDF = 7.89/6}$) by a power-law function with index $Γ= 2.46 \pm 0.06_{stat} \pm 0.20_{sys} $ and extends to at least 8 TeV. The extragalactic-background-light-deabsorbed spectrum is adequately fit (${χ^2/NDF = 6.73/6}$) by a power-law function with index $Γ= 2.15 \pm 0.06_{stat} \pm 0.20_{sys} $ while an F-test indicates that the power-law with exponential cutoff function provides a marginally-better fit ($χ^2/NDF $ = $2.56 / 5 $) at the 2.1$σ$ level. The source location is found to be consistent with the published radio location and its spatial extent is consistent with a point source. △ Less

Submitted 9 August, 2017; originally announced August 2017.

Comments: 7 pages, 2 figures. Published in Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 471, Issue 2, 21 October 2017, Pages 2117-2123

Dark Matter Constraints from a Joint Analysis of Dwarf Spheroidal Galaxy Observations with VERITAS

Authors: VERITAS Collaboration, S. Archambault, A. Archer, W. Benbow, R. Bird, E. Bourbeau, T. Brantseg, M. Buchovecky, J. H. Buckley, V. Bugaev, K. Byrum, M. Cerruti, J. L. Christiansen, M. P. Connolly, W. Cui, M. K. Daniel, Q. Feng, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, A. Geringer-Sameth, S. Griffin, J. Grube, M. Hütten , et al. (47 additional authors not shown)

Abstract: We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We… ▽ More We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($τ^{+}τ^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($γγ$) final state. △ Less

Submitted 8 May, 2017; v1 submitted 15 March, 2017; originally announced March 2017.

Comments: 14 pages, 9 figures, published in PRD, Ascii tables containing annihilation cross sections limits are available for download as ancillary files with readme.txt file description of limits

Journal ref: Physical Review D 95, 082001 (2017)

Gamma-ray Observations Under Bright Moonlight with VERITAS

Authors: S. Archambault, A. Archer, W. Benbow, R. Bird, E. Bourbeau, A. Bouvier, M. Buchovecky, V. Bugaev, J. V Cardenzana, M. Cerruti, L. Ciupik, M. P. Connolly, W. Cui, M. K. Daniel, M. Errando, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, G. H. Gillanders, S. Griffin, D. Hanna, O. Hervet , et al. (40 additional authors not shown)

Abstract: Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive photomultiplier tube (PMT) cameras. Exposure to high levels of background illumination degrades the efficiency of and potentially destroys these photo-detectors over time, so IACTs cannot be operated in the same configuration in the presence of bright moonlight as under dark skies. Since September 2012, observations have… ▽ More Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive photomultiplier tube (PMT) cameras. Exposure to high levels of background illumination degrades the efficiency of and potentially destroys these photo-detectors over time, so IACTs cannot be operated in the same configuration in the presence of bright moonlight as under dark skies. Since September 2012, observations have been carried out with the VERITAS IACTs under bright moonlight (defined as about three times the night-sky-background (NSB) of a dark extragalactic field, typically occurring when Moon illumination > 35%) in two observing modes, firstly by reducing the voltage applied to the PMTs and, secondly, with the addition of ultra-violet (UV) bandpass filters to the cameras. This has allowed observations at up to about 30 times previous NSB levels (around 80% Moon illumination), resulting in 30% more observing time between the two modes over the course of a year. These additional observations have already allowed for the detection of a flare from the 1ES 1727+502 and for an observing program targeting a measurement of the cosmic-ray positron fraction. We provide details of these new observing modes and their performance relative to the standard VERITAS observations. △ Less

Submitted 3 March, 2017; originally announced March 2017.

Multiwavelength follow-up of a rare IceCube neutrino multiplet

Authors: M. G. Aartsen, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, I. Al Samarai, D. Altmann, K. Andeen, T. Anderson, I. Ansseau, G. Anton, M. Archinger, C. Argüelles, J. Auffenberg, S. Axani, X. Bai, S. W. Barwick, V. Baum, R. Bay, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, D. Berley , et al. (479 additional authors not shown)

Abstract: On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at… ▽ More On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns. △ Less

Submitted 28 November, 2017; v1 submitted 20 February, 2017; originally announced February 2017.

Comments: 23 pages, 13 figures, accepted by A&A on July 30 2017

Journal ref: A&A 607, A115 (2017)

Gamma-ray observations of Tycho's SNR with VERITAS and Fermi

Authors: VERITAS Collaboration, S. Archambault, A. Archer, W. Benbow, R. Bird, E. Bourbeau, M. Buchovecky, J. H. Buckley, V. Bugaev, M. Cerruti, M. P. Connolly, W. Cui, V. V. Dwarkadas, M. Errando, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, S. Griffin, M. Hutten, D. Hanna, J. Holder, C. A. Johnson , et al. (42 additional authors not shown)

Abstract: High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it is a young, type Ia SNR that is well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-… ▽ More High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it is a young, type Ia SNR that is well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho's SNR with 147 hours of VERITAS and 84 months of Fermi-LAT observations, which represents about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is $2.92 \pm 0.42_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$. It is also softer than the spectral index in the GeV energy range, $2.14 \pm 0.09_{\mathrm{stat}} \pm 0.02_{\mathrm{sys}}$, measured by this study using Fermi--LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi--LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: Accepted for publication in The Astrophysical Journal

A luminous and isolated gamma-ray flare from the blazar B2 1215+30

Authors: VERITAS Collaboration, A. U. Abeysekara, S. Archambault, A. Archer, W. Benbow, R. Bird, M. Buchovecky, J. H. Buckley, V. Bugaev, K. Byrum, M. Cerruti, X. Chen, L. Ciupik, W. Cui, H. J. Dickinson, J. D. Eisch, M. Errando, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, G. H. Gillanders, S. Griffin , et al. (62 additional authors not shown)

Abstract: B2 1215+30 is a BL Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes, and subsequently confirmed by the VERITAS observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2 1215+30 during routine monitoring observations of the blazar 1ES 1218+304, located in the same field of view. T… ▽ More B2 1215+30 is a BL Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes, and subsequently confirmed by the VERITAS observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2 1215+30 during routine monitoring observations of the blazar 1ES 1218+304, located in the same field of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability timescale of < 3.6 h. Multiwavelength observations with Fermi-LAT, Swift, and the Tuorla observatory revealed a correlated high GeV flux state and no significant optical counterpart to the flare, with a spectral energy distribution where the gamma-ray luminosity exceeds the synchrotron luminosity. When interpreted in the framework of a one-zone leptonic model, the observed emission implies a high degree of beaming, with Doppler factor > 10, and an electron population with spectral index < 2.3. △ Less

Submitted 4 January, 2017; originally announced January 2017.

Comments: 8 pages, 3 figures, accepted for publication in The Astrophysical Journal

Journal ref: ApJ, 836, 205 (2017)

Search for Magnetically Broadened Cascade Emission From Blazars with VERITAS

Authors: VERITAS Collaboration, S. Archambault, A. Archer, W. Benbow, M. Buchovecky, V. Bugaev, M. Cerruti, M. P. Connolly, W. Cui, A. Falcone, M. Fernández Alonso, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, S. Griffin, M. Hütten, O. Hervet, J. Holder, T. B. Humensky, C. A. Johnson, P. Kaaret, P. Kar, D. Kieda, M. Krause , et al. (31 additional authors not shown)

Abstract: We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGN), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma rays from AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background (CMB) photons via inverse-Comp… ▽ More We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGN), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma rays from AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background (CMB) photons via inverse-Compton scattering to produce gamma rays. Due to the deflection of the electron-positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218+304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around $10^{-14}$G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission. △ Less

Submitted 2 January, 2017; originally announced January 2017.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 835:288, 2017 February 1

Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (268 additional authors not shown)

Abstract: We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux… ▽ More We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties. △ Less

Submitted 30 December, 2016; originally announced December 2016.

Comments: 33 pages, 15 figures, accepted for publication in A&A. Corresponding authors: Marlene Doert (marlene.doert@tu-dortmund.de) and David Paneque (dpaneque@mppmu.mpg.de)

Journal ref: A&A 603, A31 (2017)

A search for spectral hysteresis and energy-dependent time lags from X-ray and TeV gamma-ray observations of Mrk 421

Authors: A. U. Abeysekara, S. Archambault, A. Archer, W. Benbow, R. Bird, M. Buchovecky, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, X. Chen, L. Ciupik, M. P. Connolly, W. Cui, J. D. Eisch, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, A. Flinders, L. Fortson, A. Furniss, S. Griffin, M. Hütten, N. Håkansson , et al. (198 additional authors not shown)

Abstract: Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421,… ▽ More Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies. △ Less

Submitted 14 November, 2016; originally announced November 2016.

Comments: 45 pages, 15 figures

Very High-Energy Gamma-Ray Follow-Up Program Using Neutrino Triggers from IceCube

Authors: IceCube Collaboration, M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, K. Andeen, T. Anderson, I. Ansseau, G. Anton, M. Archinger, C. Arguelles, J. Auffenberg, S. Axani, X. Bai, S. W. Barwick, V. Baum, R. Bay, J. J. Beatty, J. Becker-Tjus, K. -H. Becker, S. BenZvi , et al. (519 additional authors not shown)

Abstract: We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-f… ▽ More We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015. △ Less

Submitted 12 November, 2016; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: accepted for publication in JINST

Journal ref: 2016 JINST 11 P11009

A Search for Very High-Energy Gamma Rays from the Missing Link Binary Pulsar J1023+0038 with VERITAS

Authors: E. Aliu, S. Archambault, A. Archer, W. Benbow, R. Bird, J. Biteau, M. Buchovecky, J. H. Buckley, V. Bugaev, K. Byrum, J. V Cardenzana, M. Cerruti, X. Chen, L. Ciupik, M. P. Connolly, W. Cui, H. J. Dickinson, J. D. Eisch, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, A. Flinders, P. Fortin, L. Fortson , et al. (60 additional authors not shown)

Abstract: The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259--63/LS 2883, making it an ideal candidate for the study of high-energy non-thermal emission. It has been the subject of multi-wavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk a… ▽ More The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259--63/LS 2883, making it an ideal candidate for the study of high-energy non-thermal emission. It has been the subject of multi-wavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high-energy gamma-ray observations carried out by VERITAS before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that very high-energy gamma rays are produced via an inverse-Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than $\sim$2 G before the disappearance of the radio pulsar and greater than $\sim$10 G afterwards. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: 7 pages, 3 figures, accepted for publication in ApJ

Very-high-energy observations of the binaries V 404 Cyg and 4U 0115+634 during giant X-ray outbursts

Authors: A. Archer, W. Benbow, R. Bird, E. Bourbeau, M. Buchovecky, J. H. Buckley, V. Bugaev, K. Byrum, M. Cerruti, M. P. Connolly, W. Cui, M. Errando, A. Falcone, Q. Feng, M. Fernandez-Alonso, J. P. Finley, H. Fleischhack, A. Flinders, L. Fortson, A. Furniss, S. Griffin, J. Grube, M. Hütten, D. Hanna, O. Hervet , et al. (40 additional authors not shown)

Abstract: Transient X-ray binaries produce major outbursts in which the X-ray flux can increase over the quiescent level by factors as large as $10^7$. The low-mass X-ray binary V 404 Cyg and the high-mass system 4U 0115+634 underwent such major outbursts in June and October 2015, respectively. We present here observations at energies above hundreds of GeV with the VERITAS observatory taken during some of t… ▽ More Transient X-ray binaries produce major outbursts in which the X-ray flux can increase over the quiescent level by factors as large as $10^7$. The low-mass X-ray binary V 404 Cyg and the high-mass system 4U 0115+634 underwent such major outbursts in June and October 2015, respectively. We present here observations at energies above hundreds of GeV with the VERITAS observatory taken during some of the brightest X-ray activity ever observed from these systems. No gamma-ray emission has been detected by VERITAS in 2.5 hours of observations of the microquasar V 404 Cyg from 2015, June 20-21. The upper flux limits derived from these observations on the gamma-ray flux above 200 GeV of F $< 4.4\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ correspond to a tiny fraction (about $10^{-6}$) of the Eddington luminosity of the system, in stark contrast to that seen in the X-ray band. No gamma rays have been detected during observations of 4U 0115+634 in the period of major X-ray activity in October 2015. The flux upper limit derived from our observations is F $< 2.1\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ for gamma rays above 300 GeV, setting an upper limit on the ratio of gamma-ray to X-ray luminosity of less than 4%. △ Less

Submitted 23 August, 2016; originally announced August 2016.

Comments: Accepted for publication in the Astrophysical Journal

Discovery of Very High Energy Gamma Rays from 1ES 1440+122

Authors: VERITAS Collaboration, S. Archambault, A. Archer, A. Barnacka, B. Behera, M. Beilicke, W. Benbow, K. Berger, R. Bird, M. Boettcher, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, X. Chen, J. L. Christiansen, L. Ciupik, E. Collins-Hughes, M. P. Connolly, W. Cui, H. J. Dickinson, J. Dumm, J. D. Eisch, M. Errando, A. Falcone , et al. (65 additional authors not shown)

Abstract: The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85 GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes. The observations, taken between 2008 May and 2010 June and totalling 53 hours, resulted in the discovery of $γ$-ray emission from the blazar, which has a redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of 5.5 standar… ▽ More The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85 GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes. The observations, taken between 2008 May and 2010 June and totalling 53 hours, resulted in the discovery of $γ$-ray emission from the blazar, which has a redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of 5.5 standard deviations above the background with an integral flux of (2.8$\pm0.7_{\mathrm{stat}}\pm0.8_{\mathrm{sys}}$) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ (1.2\% of the Crab Nebula's flux) above 200 GeV. The measured spectrum is described well by a power law from 0.2 TeV to 1.3 TeV with a photon index of 3.1 $\pm$ 0.4$_{\mathrm{stat}}$ $\pm$ 0.2$_{\mathrm{sys}}$. Quasi-simultaneous multi-wavelength data from the Fermi Large Area Telescope (0.3--300 GeV) and the Swift X-ray Telescope (0.2--10 keV) are additionally used to model the properties of the emission region. A synchrotron self-Compton model produces a good representation of the multi-wavelength data. Adding an external-Compton or a hadronic component also adequately describes the data. △ Less

Submitted 9 August, 2016; originally announced August 2016.

Comments: 8 pages, 4 figures. Accepted for publication in MNRAS

Very High Energy outburst of Markarian 501 in May 2009

Authors: E. Aliu, S. Archambault, A. Archer, T. Arlen, T. Aune, A. Barnacka, B. Behera, M. Beilicke, W. Benbow, K. Berger, R. Bird, A. Bouvier, M. Böttcher, M. Buchovecky, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, A. Cesarini, X. Chen, L. Ciupik, E. Collins-Hughes, M. P. Connolly, W. Cui, J. Dumm , et al. (86 additional authors not shown)

Abstract: The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $γ$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux (… ▽ More The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $γ$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation with an increase of a factor $\sim$4 in 25 minutes, and a falling time of $\sim$50 minutes. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data. △ Less

Submitted 4 August, 2016; originally announced August 2016.

Comments: 13 pages, 10 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 594, A76 (2016)