1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like “orphan” TeV flares and broad emission in the high-energy regime that are difficult to interpret ...using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the
Fermi
Large Area Telescope (LAT) and
Swift
instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE,
E
> 100 GeV)
γ
-ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching ∼3 times the Crab Nebula flux was measured by MAGIC on the 13 and 14 of June, and 1 July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the
γ
-ray (both
Fermi
-LAT and VHE bands) and the X-ray bands are quite hard. On 13 June and 1 July 2016, the source showed rapid variations in the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to large values of the Doppler factors (
δ
≥ 30−60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power
L
jet
∼ 10
46
erg s
−1
and under high values of the magnetic field strength (∼100 G) and maximum proton energy (∼few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016.
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.
ABSTRACT
We report a characterization of the multiband flux variability and correlations of the nearby (z = 0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, ...FACT, and other collaborations and instruments from 2014 November till 2016 June. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F>1 TeV < 1.7 × 10−12 ph cm−2 s−1) and in the 2–10 keV (X-ray) band (F$_{2-10\, \mathrm{keV}}\lt 3.6\times 10^{-11}$ erg cm−2 s−1), during which the Swift-BAT data suggest an additional spectral component beyond the regular synchrotron emission. The highest flux variability occurs in X-rays and very high-energy (E > 0.1 TeV) γ-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HRkeV and HRTeV show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 d centred at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 d centred at a time lag of $43^{+9}_{-6}$ d. This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a lognormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process.
Abstract
Extreme high-frequency-peaked BL Lac objects (EHBLs) are blazars that exhibit extremely energetic synchrotron emission. They also feature nonthermal gamma-ray emission whose peak lies in the ...very high-energy (VHE,
E
> 100 GeV) range, and in some sources exceeds 1 TeV: this is the case for hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, 10 targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 265 hr of good-quality data. The data were complemented by coordinated
Swift
observations. The X-ray data analysis confirms that all but two sources are EHBLs. The sources show only a modest variability and a harder-when-brighter behavior, typical for this class of objects. At VHE gamma-rays, three new sources were detected and a hint of a signal was found for another new source. In each case, the intrinsic spectrum is compatible with the hypothesis of a hard-TeV nature of these EHBLs. The broadband spectral energy distributions (SEDs) of all sources are built and modeled in the framework of a single-zone, purely leptonic model. The VHE gamma-ray-detected sources were also interpreted with a spine–layer model and a proton synchrotron model. The three models provide a good description of the SEDs. However, the resulting parameters differ substantially in the three scenarios, in particular the magnetization parameter. This work presents the first mini catalog of VHE gamma-ray and multiwavelength observations of EHBLs.
We report the detection of pulsed gamma-ray emission from the Geminga pulsar (PSR J0633+1746) between 15 GeV and 75 GeV. This is the first time a middle-aged pulsar has been detected up to these ...energies. Observations were carried out with the MAGIC telescopes between 2017 and 2019 using the low-energy threshold Sum-Trigger-II system. After quality selection cuts, ∼80 h of observational data were used for this analysis. To compare with the emission at lower energies below the sensitivity range of MAGIC, 11 years of
Fermi
-LAT data above 100 MeV were also analysed. From the two pulses per rotation seen by
Fermi
-LAT, only the second one,
P
2, is detected in the MAGIC energy range, with a significance of 6.3
σ
. The spectrum measured by MAGIC is well-represented by a simple power law of spectral index Γ = 5.62 ± 0.54, which smoothly extends the
Fermi
-LAT spectrum. A joint fit to MAGIC and
Fermi
-LAT data rules out the existence of a sub-exponential cut-off in the combined energy range at the 3.6
σ
significance level. The power-law tail emission detected by MAGIC is interpreted as the transition from curvature radiation to Inverse Compton Scattering of particles accelerated in the northern outer gap.
We report on observations of the pulsar/Be star binary system PSR J2032+4127/MT91 213 in the energy range between and with the Very Energetic Radiation Imaging Telescope Array and Major Atmospheric ...Gamma Imaging Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new point-like gamma-ray source is detected, coincident with the location of PSR J2032+4127/MT91 213. The gamma-ray light curve and spectrum are well characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar/Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around . This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope on board the Neil Gehrels Swift Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission that we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.
ABSTRACT
Extreme high-energy-peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump-structured spectral energy distribution (SED) peaks at higher energies with respect ...to conventional blazars. Multiwavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073326.7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within an MWL context, including an analysis of the Fermi-LAT data over 10 yr of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EHBL. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broad-band SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-band SED.
Abstract
The results of gamma-ray observations of the binary system HESS J0632 + 057 collected during 450 hr over 15 yr, between 2004 and 2019, are presented. Data taken with the atmospheric ...Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the H
α
emission line. A significant detection of the modulation of the very high-energy gamma-ray fluxes with a period of 316.7 ± 4.4 days is reported, consistent with the period of 317.3 ± 0.7 days obtained with a refined analysis of X-ray data. The analysis of data from four orbital cycles with dense observational coverage reveals short-timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over a timescale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but cannot find any correlation of optical H
α
parameters with fluxes at X-ray or gamma-ray energies in simultaneous observations. The key finding is that the emission of HESS J0632 + 057 in the X-ray and gamma-ray energy bands is highly variable on different timescales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.
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.