ABSTRACT
In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC ...telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC’s standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory’s Northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths.
Context.
It has become evident that one-zone synchrotron self-Compton models are not always adequate for very high-energy (VHE) gamma-ray-emitting blazars. While two-component models perform better, ...they are difficult to constrain due to the large number of free parameters.
Aims.
In this work, we make a first attempt at taking into account the observational constraints from very long baseline interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays), and optical polarisation to limit the parameter space for a two-component model and test whether or not it can still reproduce the observed spectral energy distribution (SED) of the blazars.
Methods.
We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modelling of the optical polarisation variability of the objects. We selected eight epochs for these five objects based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modelled with a two-component model.
Results.
We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behaviour of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempt to use optical polarisation data to shed new light on the behaviour of the two components in the optical band. Our observationally constrained two-component model allows explanation of the entire SED from radio to VHE with two co-located emission regions.
ABSTRACT
PG 1553 + 113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high energy (VHE; >100 GeV) gamma-ray emitter. To better ...understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low- and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required.
Starburst galaxies and star-forming active galactic nuclei are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the ...brightest of the star-forming galaxies found to emit gamma-rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hr of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very-high-energy band. We did not detect significant gamma-ray emission, and set upper limits at the 95% confidence level to the gamma-ray flux above 200 GeV f < 5.1 × 10−13 cm−2 s−1. This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr−1.
Aims. We have performed the first broadband study of Mrk 421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. Methods. The data were collected as part ...of an extensive multiwavelength campaign carried out between May and June 2022 using MAGIC, Fermi -LAT, NuSTAR , XMM-Newton , Swift , and several optical and radio telescopes to complement IXPE data. Results. During the IXPE exposures, the measured 0.2–1 TeV flux was close to the quiescent state and ranged from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the very high-energy (VHE) and X-ray emission are positively correlated at a 4 σ significance level. The IXPE measurements reveal an X-ray polarization degree that is a factor of 2–5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift -XRT monitoring reveals an X-ray flux increase with a clear spectral hardening. This suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counterclockwise), implying important changes in the particle acceleration efficiency on approximately hour timescales.
ABSTRACT
A deep survey of the Large Magellanic Cloud at ∼0.1–100 TeV photon energies with the Cherenkov Telescope Array is planned. We assess the detection prospects based on a model for the emission ...of the galaxy, comprising the four known TeV emitters, mock populations of sources, and interstellar emission on galactic scales. We also assess the detectability of 30 Doradus and SN 1987A, and the constraints that can be derived on the nature of dark matter. The survey will allow for fine spectral studies of N 157B, N 132D, LMC P3, and 30 Doradus C, and half a dozen other sources should be revealed, mainly pulsar-powered objects. The remnant from SN 1987A could be detected if it produces cosmic-ray nuclei with a flat power-law spectrum at high energies, or with a steeper index 2.3–2.4 pending a flux increase by a factor of >3–4 over ∼2015–2035. Large-scale interstellar emission remains mostly out of reach of the survey if its >10 GeV spectrum has a soft photon index ∼2.7, but degree-scale 0.1–10 TeV pion-decay emission could be detected if the cosmic-ray spectrum hardens above >100 GeV. The 30 Doradus star-forming region is detectable if acceleration efficiency is on the order of 1−10 per cent of the mechanical luminosity and diffusion is suppressed by two orders of magnitude within <100 pc. Finally, the survey could probe the canonical velocity-averaged cross-section for self-annihilation of weakly interacting massive particles for cuspy Navarro–Frenk–White profiles.
Aims. The BL Lac 1ES 2344+514 is known for temporary extreme properties characterised by a shift of the synchrotron spectral energy distribution (SED) peak energy ν synch, p above 1 keV. While those ...extreme states have only been observed during high flux levels thus far, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to very high energy (VHE) performed so far, focussing on a systematic characterisation of the intermittent extreme states. Methods. We organised a monitoring campaign covering a 3-year period from 2019 to 2021. More than ten instruments participated in the observations in order to cover the emission from radio to VHE. In particular, sensitive X-ray measurements by XMM-Newton , NuSTAR , and AstroSat took place simultaneously with multi-hour MAGIC observations, providing an unprecedented constraint of the two SED components for this blazar. Results. While our results confirm that 1ES 2344+514 typically exhibits ν synch, p > 1 keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase in the electron acceleration efficiency without a change in the electron injection luminosity. On the other hand, we also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of ν synch, p . During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. By combining Swift -XRT and Swift -UVOT measurements during a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions significantly contribute to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton , NuSTAR , and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3−2 keV band. Using time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model.
Aims.
Large-Sized Telescope 1 (LST-1), the prototype for the Large-Sized Telescope at the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning phase at the Observatorio del ...Roque de los Muchachos on the island of La Palma. The proximity of LST-1 to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes makes it possible to carry out observations of the same gamma-ray events with both systems.
Methods.
We describe the joint LST-1+MAGIC analysis pipeline and used simultaneous Crab Nebula observations and Monte Carlo simulations to assess the performance of the three-telescope system. The addition of the LST-1 telescope allows for the recovery of events in which one of the MAGIC images is too dim to survive analysis quality cuts.
Results.
Thanks to the resulting increase in the collection area and stronger background rejection, we found a significant improvement in sensitivity, allowing for the detection of 30% weaker fluxes in the energy range between 200 GeV and 3 TeV. The spectrum of the Crab Nebula, reconstructed in the energy range between ~60 GeV and ~10 TeV, is in agreement with previous measurements.
ABSTRACT
Extreme high-frequency BL Lacs (EHBL) feature their synchrotron peak of the broad-band spectral energy distribution (SED) at νs ≥ 1017 Hz. The BL Lac object 1ES 2344+514 was included in the ...EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In 2016 August, 1ES 2344+514 was detected with the ground-based γ-ray telescope FACT during a high γ-ray state, triggering multiwavelength (MWL) observations. We studied the MWL light curves of 1ES 2344+514 during the 2016 flaring state, using data from radio to very-high-energy (VHE) γ-rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, Swift-UVOT, Swift-XRT, Fermi-LAT, FACT, and MAGIC. With simultaneous observations of the flare, we built the broad-band SED and studied it in the framework of a leptonic and a hadronic model. The VHE γ-ray observations show a flux level of 55 per cent of the Crab Nebula flux above 300 GeV, similar to the historical maximum of 1995. The combination of MAGIC and Fermi-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The Γ index of the intrinsic spectrum in the VHE γ-ray band is 2.04 ± 0.12stat ± 0.15sys. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to 1018 Hz.
Context. Diffusive shock acceleration (DSA) is the most promising mechanism that accelerates Galactic cosmic rays (CRs) in the shocks of supernova remnants (SNRs). It is based on particles scattering ...caused by turbulence ahead and behind the shock. The turbulence upstream is supposedly generated by the CRs, but this process is not well understood. The dominant mechanism may depend on the evolutionary state of the shock and can be studied via the CRs escaping upstream into the interstellar medium (ISM). Aims. Previous observations of the γ Cygni SNR showed a difference in morphology between GeV and TeV energies. Since this SNR has the right age and is at the evolutionary stage for a significant fraction of CRs to escape, our aim is to understand γ-ray emission in the vicinity of the γ Cygni SNR. Methods. We observed the region of the γ Cygni SNR with the MAGIC Imaging Atmospheric Cherenkov telescopes between 2015 May and 2017 September recording 87 h of good-quality data. Additionally, we analysed Fermi-LAT data to study the energy dependence of the morphology as well as the energy spectrum in the GeV to TeV range. The energy spectra and morphology were compared against theoretical predictions, which include a detailed derivation of the CR escape process and their γ-ray generation. Results. The MAGIC and Fermi-LAT data allowed us to identify three emission regions that can be associated with the SNR and that dominate at different energies. Our hadronic emission model accounts well for the morphology and energy spectrum of all source components. It constrains the time-dependence of the maximum energy of the CRs at the shock, the time-dependence of the level of turbulence, and the diffusion coefficient immediately outside the SNR shock. While in agreement with the standard picture of DSA, the time-dependence of the maximum energy was found to be steeper than predicted, and the level of turbulence was found to change over the lifetime of the SNR.