Search for: All records

Abstract Microquasars such as SS 433 are considered potential contributors to cosmic rays up to the knee of the cosmic-ray energy spectrum (∼4 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 hr of observations with VERITAS, sensitive toγ-ray energies >100 GeV. With an angular resolution better than 0 $\stackrel{\text{°}}{.}$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 in the 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 that 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. 

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 Large High Altitude Air Shower Observatory (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 broadband emission. Models of the spectral energy distribution (SED) are currently unconstrained due to the lack of contemporaneous multiwavelength data at either peak. Here, we report serendipitous observations of NGC 4278 with the Very Energetic Radiation Imaging Telescope Array System (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° 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 active galactic nuclei, followed byp–γinteractions and optically thinγ-ray emission. Within this framework, the implied neutrino signal is slightly below the current sensitivity of IceCube. 

Abstract In 2017 September, 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 Δχ²-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 an LP primary photon spectrum, we find consistency with a proton injection spectral index ofα_p ≃ 2.0 and a cutoff energy ofE_p,max ≃ 1.3 × 10¹⁶eV and constrain the isotropic proton escape luminosity to 1 × 10⁴⁴erg s⁻¹≲L_p,esc ≲ 3 × 10⁴⁵erg s⁻¹at the 90% confidence level. 

Abstract We use the stellar intensity interferometry system implemented with the Very Energetic Radiation Imaging Telescope Array System at Fred Lawrence Whipple Observatory 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 ± 0.02 mas, a major-to-minor-radius ratio of 1.28 ± 0.04, and a position angle of 116° ± 5° 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.034}^{+0.041}$mas corresponding to an equatorial radius of ${10.9}_{-0.6}^{+0.8}{R}_{\odot }$, a rotational velocity with a 1σlower limit at 97.7% of breakup velocity, and a position angle of $114\stackrel{\mathrm{°}}{.}{7}_{-5.7}^{+6.4}$. 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. 

Abstract Observations of GeV gamma-ray emission from the well-studied mixed-morphology supernova remnant (SNR) W44 by Fermi-Large Area Telescope and AGILE imply that it is a site of significant cosmic-ray acceleration. The spectral energy distribution (SED) derived from the GeV data suggests 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  × 10⁻¹³erg cm⁻²s⁻¹in the 0.5–8.0 keV band based on  ∼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  ∼1.2  × 10⁻¹²erg s⁻¹cm⁻²in the 0.5–5.0 TeV range from W44 using data from the Very Energetic Radiation Imaging Telescope Array System. Our nondetection of W44 at VHE wavelengths is in agreement with observations from other imaging atmospheric Cherenkov telescopes and is perhaps consistent with the evolutionary stage of the SNR. 

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. 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 $\mathcal{O}\left(10\mathrm{\mu }\mathrm{G}\right)$through which escaping photons traverse, making it an excellent target to study ALPs. 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 ${\mathrm{CL}}_{\mathrm{s}}$method for this study for its conservative approach to setting limits in regimes where the search has limited sensitivity. 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}\mathrm{eV}\lesssim m_a\lesssim 4\times {10}^{-7}\mathrm{eV}$and at the coupling strength of $g_{a\gamma }\gtrsim 3\times {10}^{-11}{\mathrm{GeV}}^{-1}$up to 80% confidence level, which are consistent with previous studies. We find the ${\mathrm{CL}}_{\mathrm{s}}$to be a trustworthy approach, and advocate for its continued usage in future studies. We note that many of the limitations contributing to the limited sensitivity seen by VERITAS in this study will be improved with next-generation $\gamma$-ray instruments, such as the Cherenkov Telescope Array Observatory. 

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 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^∘and $10^{\prime }$around PSR J0622+3749, respectively. Combined with the LHAASO Kilometer Square Array (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μG 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. 

Abstract We present an analysis of new multiwavelength 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ϕ ≈ 0.3–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. 

Abstract Assuming Galactic cosmic rays originate in supernovae and the winds of massive stars, starburst galaxies should produce very-high-energy (VHE;E > 100 GeV) gamma-ray emission via the interaction of their copious quantities of cosmic rays with the large reservoirs of dense gas within the galaxies. Such VHE emission was detected by VERITAS from the starburst galaxy M82 in 2008–09. An extensive, multiyear campaign followed these initial observations, yielding a total of 254 hr of good-quality VERITAS data on M82. Leveraging modern analysis techniques and the larger exposure, these VERITAS data show a more statistically significant VHE signal (∼6.5 standard deviations,σ). The corresponding photon spectrum is well fit by a power law (Γ = 2.3 ± 0.3_stat ± 0.2_sys), and the observed integral flux isF(>450 GeV) = (3.2 ± 0.6_stat ± 0.6_sys) × 10⁻¹³cm⁻²s⁻¹, or ∼0.4% of the Crab Nebula flux above the same energy threshold. The improved VERITAS measurements, when combined with various multiwavelength data, enable modeling of the underlying emission and transport processes. A purely leptonic scenario is found to be a poor representation of the gamma-ray spectral energy distribution (SED). A lepto-hadronic scenario with cosmic rays following a power-law spectrum in momentum (indexs ≃ 2.25) and with significant bremsstrahlung below 1 GeV provides a good match to the observed SED. The synchrotron emission from the secondary electrons indicates that efficient nonradiative losses of cosmic-ray electrons may be related to advective escape from the starburst core. 

Abstract We use the Very Energetic Radiation Imaging telescope Array System (VERITAS) imaging air Cherenkov telescope array to obtain the first measured angular diameter ofβUMa at visual wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak (βUMa), an A1-type subgiant, by comparing effective temperature and luminosity constraints to model stellar evolution tracks. Previous interferometric limb-darkened angular-diameter measurements ofβUMa in the near-infrared (Center for High Angular Resolution Astronomy (CHARA) Array, 1.149 ± 0.014 mas) and mid-infrared (Keck Nuller, 1.08 ± 0.07 mas), together with the measured parallax and bolometric flux, have constrained the effective temperature. This paper presents current VERITAS-SII observation and analysis procedures to derive squared visibilities from correlation functions. We fit the resulting squared visibilities to find a limb-darkened angular diameter of 1.07 ± 0.04 (stat) ± 0.05 (sys) mas, using synthetic visibilities from a stellar atmosphere model that provides a good match to the spectrum ofβUMa in the optical wave band. The VERITAS-SII limb-darkened angular diameter yields an effective temperature of 9700 ± 200 ± 200 K, consistent with ultraviolet spectrophotometry, and an age of 390 ± 29 ± 32 Myr, using MESA Isochrones and Stellar Tracks. This age is consistent with 408 ± 6 Myr from the CHARA Array angular diameter.