Search for: All records

Context.Blazars, which include BL Lacs and flat-spectrum radio quasars, represent the brightest persistent extragalactic sources in the high-energy (HE; 10 MeV–100 GeV) and very-high-energy (VHE;E > 100 GeV)γ-ray sky. Due to their almost featureless optical/UV spectra, it is challenging to measure the redshifts of BL Lacs. As a result, about 50% ofγ-ray BL Lacs lack a firm measurement of this property, which is fundamental for population studies, indirect estimates of the extragalactic background light, and fundamental physics probes (e.g., searches for Lorentz-invariance violation or axion-like particles). Aims.This paper is the third in a series of papers aimed at determining the redshift of a sample of blazars selected as prime targets for future observations with the next generation, ground-based VHEγ-ray astronomy observatory, Cherenkov Telescope Array Observatory (CTAO). The accurate determination of the redshift of these objects is an important aid in source selection and planning of future CTAO observations. Methods.Promising targets were selected following a sample selection obtained with Monte Carlo simulations of CTAO observations. The selected targets were expected to be detectable with CTAO in observations of 30 h or less. We performed deep spectroscopic observations of 41 of these blazars using the Keck II, Lick, SALT, GTC, and ESO/VLT telescopes. We carefully searched for spectral lines in the spectra and whenever features of the host galaxy were detected, we attempted to model the properties of the host galaxy. The magnitudes of the targets at the time of the observations were also compared to their long-term light curves. Results.Spectra from 24 objects display spectral features or a high signal-to-noise ratio (S/N). From these, 12 spectroscopic redshifts were determined, ranging from 0.2223 to 0.7018. Furthermore, 1 tentative redshift (0.6622) and 2 redshift lower limits atz > 0.6185 andz > 0.6347 were obtained. The other 9 BL Lacs showed featureless spectra, despite the high S/N (≥100) observations. Our comparisons with long-term optical light curves tentatively suggest that redshift measurements are more straightforward during an optical low state of the active galactic nucleus. Overall, we have determined 37 redshifts and 6 spectroscopic lower limits as part of our programme thus far. 

Aims.We introduce the TELAMON program which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, specifically TeV blazars and candidate neutrino-associated AGN. Here, we present and characterize our main sample of TeV-detected blazars. Methods.We analyzed the data sample from the first ∼2.5 yr of observations between August 2020 and February 2023 in the range from 14 GHz to 45 GHz. During this pilot phase, we observed all 59 TeV-detected blazars in the Northern Hemisphere (i.e., Dec > 0°) known at the time of observation. We discuss the basic data reduction and calibration procedures used for all TELAMON data and introduce a sub-band averaging method used to calculate average light curves for the sources in our sample. Results.The TeV-selected sources in our sample exhibit a median flux density of 0.12 Jy at 20 mm, 0.20 Jy at 14 mm, and 0.60 Jy at 7 mm. The spectrum for most of the sources is consistent with a flat radio spectrum and we found a median spectral index (S(ν)∝ν^α) ofα = −0.11. Our results on flux density and spectral index are consistent with previous studies of TeV-selected blazars. Compared to the GeV-selected F-GAMMA sample, TELAMON sources are significantly fainter in the radio band. This is consistent with the double-humped spectrum of blazars being shifted towards higher frequencies for TeV-emitters (in particular for high-synchrotron peaked BL Lac type objects), which results in a lower radio flux density. The spectral index distribution of our TeV-selected blazar sample is not significantly different from the GeV-selected F-GAMMA sample. Moreover, we present a strategy to track the light curve evolution of sources in our sample for future variability and correlation analysis. 

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. 

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 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 Active galactic nuclei (AGNs) make up about 35 per cent of the more than 250 sources detected in very high-energy (VHE) gamma rays to date with the imaging atmospheric Cherenkov telescopes. Apart from four nearby radio galaxies and two AGNs of unknown type, all known VHE AGNs are blazars. Knowledge of the cosmological redshift of gamma-ray blazars is key to enabling the study of their intrinsic emission properties, as the interaction between gamma rays and the extragalactic background light (EBL) results in a spectral softening. Therefore, the redshift determination exercise is crucial to indirectly placing tight constraints on the EBL density, and to studying blazar population evolution across cosmic time. Due to the powerful relativistic jets in blazars, most of their host galaxies’ spectral features are outshined, and dedicated high signal-to-noise (S/N) spectroscopic observations are required. Deep medium- to high-resolution spectroscopy of 33 gamma-ray blazar optical counterparts was performed with the European Southern Observatory, New Technology Telescope, Keck II telescope, Shane 3-metre telescope, and the Southern African Large Telescope. From the sample, spectra from 25 objects display spectral features or are featureless and have high S/N. The other eight objects have low-quality featureless spectra. We systematically searched for absorption and emission features and estimated, when possible, the fractional host galaxy flux in the measured total flux. Our measurements yielded 14 firm spectroscopic redshifts, ranging from 0.0838 to 0.8125, one tentative redshift, and two lower limits: one at $z > 0.382$ and the other at z > 0.629. 

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. 

Context. Standing and moving shocks in relativistic astrophysical jets are very promising sites for particle acceleration to large Lorentz factors and for the emission from the radio up to the γ -ray band. They are thought to be responsible for at least part of the observed variability in radio-loud active galactic nuclei. Aims. We aim to simulate the interactions of moving shock waves with standing recollimation shocks in structured and magnetized relativistic jets and to characterize the profiles of connected flares in the radio light curve. Methods. Using the relativistic magneto-hydrodynamic code MPI-AMRVAC and a radiative transfer code in post-processing, we explore the influence of the magnetic-field configuration and transverse stratification of an over-pressured jet on its morphology, on the moving shock dynamics, and on the emitted radio light curve. First, we investigate different large-scale magnetic fields with their effects on the standing shocks and on the stratified jet morphology. Secondly, we study the interaction of a moving shock wave with the standing shocks. We calculated the synthetic synchrotron maps and radio light curves and analyze the variability at two frequencies 1 and 15.3 GHz and for several observation angles. Finally, we compare the characteristics of our simulated light curves with radio flares observed from the blazar 3C 273 with the Owens Valley Radio Observatory and Very Long Baseline Array in the MOJAVE survey between 2008 and 2019. Results. We find that in a structured over-pressured relativistic jet, the presence of the large-scale magnetic field structure changes the properties of the standing shock waves and leads to an opening in the jet. The interaction between waves from inner and outer jet components can produce strong standing shocks. When crossing such standing shocks, moving shock waves accompanying overdensities injected in the base of the jet cause very luminous radio flares. The observation of the temporal structure of these flares under different viewing angles probes the jet at different optical depths. At 1 GHz and for small angles, the self-absorption caused by the moving shock wave becomes more important and leads to a drop in the observed flux after it interacts with the brightest standing knot. A weak asymmetry is seen in the shape of the simulated flares, resulting from the remnant emission of the shocked standing shocks. The characteristics of the simulated flares and the correlation of peaks in the light curve with the crossing of moving and standing shocks favor this scenario as an explanation of the observed radio flares of 3C 273. 

Understanding the nature and identity of dark matter is a key goal in the physics community. In the case that TeV-scale dark matter particles decay or annihilate into standard model particles, very-high-energy (VHE) gamma rays (greater than 100 GeV) will be present in the final state. The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an imaging atmospheric Cherenkov telescope array that can indirectly detect VHE gamma rays in an energy range of 100 GeV to > 30 TeV. Dwarf spheroidal galaxies (dSphs) are ideal candidates in the search for dark matter due to their high dark matter content, high mass-to-light ratios, and their low gamma-ray fluxes from astrophysical processes. This study uses a legacy data set of 638 hours collected on 17 dSphs, built over 11 years with an observing strategy optimized according to the dark matter content of the targets. The study addresses a broad dark matter particle mass range, extending from 200 GeV to 30 PeV. In the absence of a detection, we set the upper limits on the dark matter velocity-weighted annihilation cross section 

Extragalactic background light (EBL) plays an important role in cosmology since it traces the history of galaxy formation and evolution. Such diffuse radiation from near-UV to far-infrared wavelengths can interact with γ -rays from distant sources such as active galactic nuclei (AGNs), and is responsible for the high-energy absorption observed in their spectra. However, probing the EBL from γ -ray spectra of AGNs is not trivial due to internal processes that can mimic its effect. Such processes are usually taken into account in terms of curvature of the intrinsic spectrum. Hence, an improper choice of parametrization for the latter can seriously affect EBL reconstruction. In this paper, we propose a statistical approach that avoids a priori assumptions on the intrinsic spectral curvature and that, for each source, selects the best-fit model on a solid statistical basis. By combining the Fermi -LAT observations of 490 blazars, we determine the γ -ray-inferred level of EBL for various state-of-the-art EBL models. We discuss the EBL level obtained from the spectra of both BL Lacs and flat spectrum radio quasars (FSRQ) in order to investigate the impact of internal absorption in different classes of objects. We further scrutinize constraints on the EBL evolution from γ -ray observations by reconstructing the EBL level in four redshift ranges, up to z  ∼ 2.5. The approach implemented in this paper, carefully addressing the question of the modeling of the intrinsic emission at the source, can serve as a solid stepping stone for studies of hundreds of high-quality spectra acquired by next-generation γ -ray instruments.