Context.
Certain types of supernova remnants (SNRs) in our Galaxy are assumed to be PeVatrons, capable of accelerating cosmic rays (CRs) to ~ PeV energies. However, conclusive observational evidence ...for this has not yet been found. The SNR G106.3+2.7, detected at 1–100 TeV energies by different γ-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape, which can be divided into a head and a tail region with different physical conditions. However, in which region the 100 TeV emission is produced has not yet been identified because of the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear as to whether the origin of the γ-ray emission is leptonic or hadronic.
Aims.
With the better angular resolution provided by new MAGIC data compared to earlier
γ
-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1° resolution at TeV energies.
Methods.
We observed the SNR G106.3+2.7 using the MAGIC telescopes for 121.7 h in total – after quality cuts – between May 2017 and August 2019. The analysis energy threshold is ~0.2 TeV, and the angular resolution is 0.07−0.1°. We examined the
γ
-ray spectra of different parts of the emission, whilst benefitting from the unprecedented statistics and angular resolution at these energies provided by our new data. We also used measurements at other wavelengths such as radio, X-rays, GeV
γ
-rays, and 10 TeV
γ
-rays to model the emission mechanism precisely.
Results.
We detect extended γ-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detect a significant
γ
-ray emission with energies above 6.0 TeV from only the tail region suggests that the emissions above 10 TeV detected with air shower experiments (Milagro, HAWC, Tibet AS
γ
and LHAASO) are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of ~1 PeV for that region. Such high-energy emission in this middle-aged SNR (4−10 kyr) can be explained by considering a scenario where protons escaping from the SNR in the past interact with surrounding dense gases at present.
Conclusions.
The
γ
-ray emission region detected with the MAGIC telescopes in the SNR G106.3+2.7 is extended and spatially coincident with the radio continuum morphology. The multi-wavelength spectrum of the emission from the tail region suggests proton acceleration up to ~PeV, while the emission mechanism of the head region could either be hadronic or leptonic.
Context. PKS 1510−089 is one of only a few flat spectrum radio quasars detected in the very-high-energy (VHE, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties ...of the PKS 1510−089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510−089 with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes during a long, high gamma-ray state in May 2015. In order to perform broadband modeling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray, and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510−089 was detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, owing to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to those obtained during previous measurements of the source. The observed flux variability sets constraints for the first time on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.
Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon ...in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400 GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of ( ) for a linear, and ( ) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, which are found to worsen the statistical limits by about 36%-42%. Our constraints would have been much more stringent if the intrinsic pulse shape of the pulsar between 200 GeV and 400 GeV was understood in sufficient detail and allowed inclusion of events well below 400 GeV.
Context. The BL Lac object S5 0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations ...were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. Aims.The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Methods. Multi-wavelength light curves have been studied together with the broadband spectral energy distributions (SEDs). The sample includes data from Effelsberg, OVRO, Metsähovi, VLBI, CARMA, IRAM, SMA, Swift-UVOT, KVA, Tuorla, Steward, RINGO3, KANATA, AZT-8+ST7, Perkins, LX-200, Swift-XRT, NuSTAR, Fermi-LAT and MAGIC. Results. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the very-high-energy (VHE, E > 100 GeV). In the constructed SED, the Swift-XRT+NuSTAR data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1 GeV to 600 GeV by Fermi+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to γ-ray data. Instead we use a two-zone model. The electric vector position angle (EVPA) shows an unprecedented fast rotation. An estimation of the redshift of the source by combined high-energy (HE, 0.1 GeV < E < 100 GeV) and VHE data provides a value of z = 0.31 ± 0.02stats ± 0.05sys, confirming the literature value. Conclusions. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out of a recollimation shock in the inner jet. A shock–shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modeling.
Millisecond pulsars, old neutron stars spun up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" ...rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311—3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.
Context.
The BL Lac object 1ES 0647+250 is one of the few distant
γ
-ray emitting blazars detected at very high energies (VHEs; ≳100 GeV) during a non-flaring state. It was detected with the MAGIC ...telescopes during a period of low activity in the years 2009−2011 as well as during three flaring activities in the years 2014, 2019, and 2020, with the highest VHE flux in the last epoch. An extensive multi-instrument data set was collected as part of several coordinated observing campaigns over these years.
Aims.
We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions.
Methods.
We evaluated the variability and correlation of the emission in the different energy bands with the fractional variability and the Z-transformed discrete correlation function, as well as its spectral evolution in X-rays and
γ
rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimated its distance in an indirect manner through a comparison of the GeV and TeV spectra from simultaneous observations with
Fermi
-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data.
Results.
We find significant long-term variability, especially in X-rays and VHE
γ
rays. Furthermore, significant (3−4
σ
) correlations were found between the radio, optical, and high-energy (HE)
γ
-ray fluxes, with the radio emission delayed by about ∼400 days with respect to the optical and
γ
-ray bands. The spectral analysis reveals a harder-when-brighter trend during the non-flaring state in the X-ray domain. However, no clear patterns were observed for either the enhanced states or the HE (30 MeV <
E
< 100 GeV) and VHE
γ
-ray emission of the source. The indirect estimation of the redshift yielded a value of
z
= 0.45 ± 0.05, which is compatible with some of the values reported in the literature. The SEDs related to the low-activity state and the three flaring states of 1ES 0647+250 can be described reasonably well with the both one-component and two-component leptonic scenarios. However, the long-term correlations indicate the need for an additional radio-producing region located about 3.6 pc downstream from the gamma-ray producing region.
Aims. We present a study of the very high-energy (VHE; E> 100 GeV) γ-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength ...spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z ≥ 0.6 and makes it a promising candidate to be the most distant VHE source. Methods. The source has been observed with the MAGIC telescopes in VHE γ rays for a total observation time of ~33.6 h from 2009 to 2011. A detailed analysis of its γ-ray spectrum and time evolution has been carried out. Moreover, we have collected and analyzed simultaneous and quasi-simultaneous multiwavelength data. Results. The source was marginally detected in VHE γ rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2σ. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 ± 1.2 to 5.0 ± 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009–2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z ~ 0.6.