-
Aims: To investigate the extension of the very-high-energy spectral tail of
the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours
of good quality data of Crab with the MAGIC telescope, obtained from February
2007 until April 2014. Results: We report the most energetic pulsed emission
ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile
shows two narrow peaks synchronized with the ones measured in the GeV energy
range. The spectra of the two peaks follow two different power-law functions
from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above
10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When
making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices
of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC
telescopes the most energetic pulsed photons from a pulsar to date. Such TeV
pulsed photons require a parent population of electrons with a Lorentz factor
of at least $5\times 10^6$. These results strongly suggest IC scattering off
low energy photons as the emission mechanism and a gamma-ray production region
in the vicinity of the light cylinder.
-
The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is
detected between 40 and 250 GeV with a significance of 25.5 {\sigma} using the
MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357
(z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar
detected to date. The observations were triggered by an outburst in 2015 April
seen at GeV energies with the Large Area Telescope on board Fermi.
Multi-wavelength observations suggest a subdivision of the high state into two
distinct flux states. In the band covered by MAGIC, the variability time scale
is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy
distribution with an external Compton model, the location of the emitting
region is understood as originating in the jet outside the broad line region
(BLR) during the period of high activity, while being partially within the BLR
during the period of low (typical) activity. The observed VHE spectrum during
the highest activity is used to probe the extragalactic background light at an
unprecedented distance scale for ground-based gamma-ray astronomy.
-
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 MAGIC collaboration. A profile likelihood analysis of pulsar
events reconstructed for energies above 400GeV 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 $E_{\mathrm{QG}_1} > 5.5\cdot 10^{17}$GeV ($4.5\cdot 10^{17}$GeV) for
a linear, and $E_{\mathrm{QG}_2} > 5.9\cdot 10^{10}$GeV ($5.3\cdot 10^{10}$GeV)
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, found to worsen the
statistical limits by about 36--42\%. Our constraints would have resulted much
more competitive if the intrinsic pulse shape of the pulsar between 200GeV and
400GeV was understood in sufficient detail and allowed inclusion of events well
below 400GeV.
-
The microquasar Cygnus X-1 displays the two typical soft and hard X-ray
states of a black-hole transient. During the latter, Cygnus X-1 shows a
one-sided relativistic radio-jet. Recent detection of the system in the high
energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT
associates this emission with the outflow. Former MAGIC observations revealed a
hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV)
regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken
with the MAGIC telescopes between July 2007 and October 2014. To shed light on
the correlation between hard X-ray and VHE gamma rays as previously suggested,
we study each main X-ray state separately. We perform an orbital phase-folded
analysis to look for variability in the VHE band. Additionally, to place this
variability behavior in a multiwavelength context, we compare our results with
\textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI},
\textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the
VHE regime. We establish upper limits for each X-ray state, assuming a
power-law distribution with photon index $\Gamma=3.2$. For steady emission in
the hard and soft X-ray states, we set integral upper limits at 95\% confidence
level for energies above 200 GeV at $2.6\times10^{-12}$~photons
cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$,
respectively. We rule out steady VHE gamma-ray emission above this energy
range, at the level of the MAGIC sensitivity, originating in the interaction
between the relativistic jet and the surrounding medium, while the emission
above this flux level produced inside the binary still remains a valid
possibility.
-
The large jet kinetic power and non-thermal processes occurring in the
microquasar SS 433 make this source a good candidate for a very high-energy
(VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the
central binary and the interaction regions between jets and surrounding nebula.
Also, non-thermal emission at lower energies has been previously reported. We
explore the capability of SS 433 to emit VHE gamma rays during periods in which
the expected flux attenuation due to periodic eclipses and precession of the
circumstellar disk periodically covering the central binary system is expected
to be at its minimum. The eastern and western SS433/W50 interaction regions are
also examined. We aim to constrain some theoretical models previously developed
for this system. We made use of dedicated observations from MAGIC and H.E.S.S.
from 2006 to 2011 which were combined for the first time and accounted for a
total effective observation time of 16.5 h. Gamma-ray attenuation does not
affect the jet/medium interaction regions. The analysis of a larger data set
amounting to 40-80 h, depending on the region, was employed. No evidence of VHE
gamma-ray emission was found. Upper limits were computed for the combined data
set. We place constraints on the particle acceleration fraction at the inner
jet regions and on the physics of the jet/medium interactions. Our findings
suggest that the fraction of the jet kinetic power transferred to relativistic
protons must be relatively small to explain the lack of TeV and neutrino
emission from the central system. At the SS433/W50 interface, the presence of
magnetic fields greater 10$\mu$G is derived assuming a synchrotron origin for
the observed X-ray emission. This also implies the presence of high-energy
electrons with energies up to 50 TeV, preventing an efficient production of
gamma-ray fluxes in these interaction regions.
-
Context. PKS 1510-089 is one of only a few flat spectrum radio quasars
detected in the VHE (very-high-energy, > 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 MAGIC telescopes during a long high
gamma-ray state in May 2015. In order to perform broadband modelling 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 has been
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, due 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 the ones obtained
during previous measurements of the source. The observed flux variability sets
for the first time constraints 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.
-
In this work we present data from observations with the MAGIC telescopes of
SN 2014J detected in January 21 2014, the closest Type Ia supernova since
Imaging Air Cherenkov Telescopes started to operate. We probe the possibility
of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early
stages of Type Ia supernova explosions. We performed follow-up observations
after this supernova explosion for 5 days, between January 27 and February 2 in
2014. We search for gamma-ray signal in the energy range between 100 GeV and
several TeV from the location of SN 2014J using data from a total of $\sim5.5$
hours of observations. Prospects for observing gamma-rays of hadronic origin
from SN 2014J in the near future are also being addressed. No significant
excess was detected from the direction of SN 2014J. Upper limits at 95$\%$
confidence level on the integral flux, assuming a power-law spectrum,
d$F/$d$E\propto E^{-\Gamma}$, with a spectral index of $\Gamma=2.6$, for
energies higher than 300 GeV and 700 GeV, are established at
$1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$,
respectively. For the first time, upper limits on the VHE emission of a Type Ia
supernova are established. The energy fraction isotropically emitted into TeV
gamma rays during the first $\sim10$ days after the supernova explosion for
energies greater than 300 GeV is limited to $10^{-6}$ of the total available
energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario,
the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be
detectable in the future by any current or planned generation of Imaging
Atmospheric Cherenkov Telescopes.
-
We present the first detection of the nearby (z=0.084) low-luminosity BL Lac
object 1ES 1741+196 in the very high energy (VHE: E$>$100 GeV) band. This
object lies in a triplet of interacting galaxies. Early predictions had
suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac
sources, within the reach of MAGIC detectability. Its detection by MAGIC, later
confirmed by VERITAS, helps to expand the small population of known TeV BL
Lacs. The source was observed with the MAGIC telescopes between 2010 April and
2011 May, collecting 46 h of good quality data. These observations led to the
detection of the source at 6.0 $\sigma$ confidence level, with a steady flux
$\mathrm{F}(> 100 {\rm GeV}) = (6.4 \pm 1.7_{\mathrm{stat}}\pm
2.6_{\mathrm{syst}}) \cdot 10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and a differential
spectral photon index $\Gamma = 2.4 \pm 0.2_{\mathrm{stat}} \pm
0.2_{\mathrm{syst}}$ in the range of $\sim$80 GeV - 3 TeV. To study the
broad-band spectral energy distribution (SED) simultaneous with MAGIC
observations, we use KVA, Swift/UVOT and XRT, and Fermi/LAT data. One-zone
synchrotron-self-Compton (SSC) modeling of the SED of 1ES 1741+196 suggests
values for the SSC parameters that are quite common among known TeV BL Lacs
except for a relatively low Doppler factor and slope of electron energy
distribution. A thermal feature seen in the SED is well matched by a giant
elliptical's template. This appears to be the signature of thermal emission
from the host galaxy, which is clearly resolved in optical observations.
-
We present an extensive study of the BL Lac object Mrk 501 based on a data
set collected during the multi-instrument campaign spanning from 2009 March 15
to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS,
Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase
in the fractional variability with energy, while no significant interband
correlations of flux changes are found in the acquired data set. The higher
variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the
lack of correlation with the X-ray emission indicate that the highest-energy
electrons that are responsible for the VHE gamma-rays do not make a dominant
contribution to the ~1 keV emission. Alternatively, there could be a very
variable component contributing to the VHE gamma-ray emission in addition to
that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC
model parameters is probed following a dedicated grid-scan strategy, allowing
for a wide range of models to be tested and offering a study of the degeneracy
of model-to-data agreement in the individual model parameters. We find that
there is some degeneracy in both the one-zone and the two-zone SSC scenarios
that were probed, with several combinations of model parameters yielding a
similar model-to-data agreement, and some parameters better constrained than
others. The SSC model grid-scan shows that the flaring activity around 2009 May
22 cannot be modeled adequately with a one-zone SSC scenario, while it can be
suitably described within a two-independent-zone SSC scenario. The observation
of an electric vector polarization angle rotation coincident with the gamma-ray
flare from 2009 May 1 resembles those reported previously for low frequency
peaked blazars, hence suggesting that there are many similarities in the
flaring mechanisms of blazars with different jet properties.
-
The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed
the BL Lac object H1722+119 (redshift unknown) for six consecutive nights
between 2013 May 17 and 22, for a total of 12.5 h. The observations were
triggered by high activity in the optical band measured by the KVA (Kungliga
Vetenskapsakademien) telescope. The source was for the first time detected in
the very high energy (VHE, $E > 100$ GeV) $\gamma$-ray band with a statistical
significance of 5.9 $\sigma$. The integral flux above 150 GeV is estimated to
be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous
high energy (HE, 100 MeV $ < E < 100$ GeV) $\gamma$-ray observations from
Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within
the framework of the current extragalactic background light models, we estimate
the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous
X-ray to radio data collected by the instruments on board the Swift satellite,
the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study
multifrequency characteristics of the source. We found no significant temporal
variability of the flux in the HE and VHE bands. The flux in the optical and
radio wavebands, on the other hand, did vary with different patterns. The
spectral energy distribution (SED) of H1722+119 shows surprising behaviour in
the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled
using an inhomogeneous helical jet synchrotron self-Compton model.
-
Context. QSO B0218+357 is a gravitationally lensed blazar located at a
redshift of 0.944. The gravitational lensing splits the emitted radiation into
two components, spatially indistinguishable by gamma-ray instruments, but
separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a
violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes.
Aims. The spectral energy distribution of QSO B0218+357 can give information on
the energetics of z ~ 1 very high energy gamma- ray sources. Moreover the
gamma-ray emission can also be used as a probe of the extragalactic background
light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during
the expected arrival time of the delayed component of the emission. The MAGIC
and Fermi-LAT observations were accompanied by quasi-simultaneous optical data
from the KVA telescope and X-ray observations by Swift-XRT. We construct a
multiwavelength spectral energy distribution of QSO B0218+357 and use it to
model the source. The GeV and sub-TeV data, obtained by Fermi-LAT and MAGIC,
are used to set constraints on the extragalactic background light. Results.
Very high energy gamma-ray emission was detected from the direction of QSO
B0218+357 by the MAGIC telescopes during the expected time of arrival of the
trailing component of the flare, making it the farthest very high energy
gamma-ray sources detected to date. The observed emission spans the energy
range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy
distribution of QSO B0218+357 is consistent with current extragalactic
background light models. The broad band emission can be modeled in the
framework of a two zone external Compton scenario, where the GeV emission comes
from an emission region in the jet, located outside the broad line region.
-
We study the multi-band variability and correlations of the TeV blazar Mrk
421 on year time scales, which can bring additional insight on the processes
responsible for its broadband emission. We observed Mrk 421 in the very high
energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March
2007 to June 2009 for a total of 96 hours of effective time after quality cuts.
The VHE flux variability is quantified with several methods, including the
Bayesian Block algorithm, which is applied to data from Cherenkov telescopes
for the first time. The 2.3 year long MAGIC light curve is complemented with
data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO,
and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an
excellent characterisation of the multi-band variability and correlations over
year time scales. Mrk 421 was found in different gamma-ray emission states
during the 2.3 year long observation period. Flares and different levels of
variability in the gamma-ray light curve could be identified with the Bayesian
Block algorithm. The same behaviour of a quiet and active emission was found in
the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct
correlation in time. The behaviour of the optical light curve of GASP-WEBT and
the radio light curves by OVRO and Mets\"ahovi are different as they show no
coincident features with the higher energetic light curves and a less variable
emission. The fractional variability is overall increasing with energy. The
comparable variability in the X-ray and VHE bands and their direct correlation
during both high- and low-activity periods spanning many months show that the
electron populations radiating the X-ray and gamma-ray photons are either the
same, as expected in the Synchrotron-Self-Compton mechanism, or at least
strongly correlated, as expected in electromagnetic cascades.
-
The gamma-ray binary LS I +61$^{\circ}$303 is a well established source from
centimeter radio up to very high energy (VHE; E$>$100 GeV). Its broadband
emission shows a periodicity of $\sim$26.5 days, coincident with the orbital
period. A longer (super-orbital) period of 1667 $\pm$ 8 days was discovered in
radio and confirmed in optical and high energy (HE; E>100 MeV) gamma-ray
observations. We present a four-year campaign performed by MAGIC together with
archival data concentrating on a search for a long timescale signature in the
VHE emission. We focus on the search for super-orbital modulation of the VHE
peak and on the search for correlations between TeV emission and optical
determination of the extension of the circumstellar disk. A four-year campaign
has been carried out by MAGIC. The source was observed during the orbital
phases when the periodic VHE outbursts have occurred ($\phi$=0.55-0.75).
Additionally, we included archival MAGIC observations and data published by the
VERITAS collaboration in these studies. For the correlation studies, LS I
+61$^{\circ}$303 has also been observed during the orbital phases where
sporadic VHE emission had been detected in the past ($\phi$=0.75-1.0). These
MAGIC observations were simultaneous with optical spectroscopy from the
LIVERPOOL telescope. The TeV flux of the periodical outburst in orbital phases
$\phi$=0.5--0.75 was found to show yearly variability consistent with the
$\sim$4.5 years long-term modulation found in the radio band. This modulation
of the TeV flux can be well described by a sine function with the best fit
period of $1610\pm 58$ days. The complete dataset span two super-orbital
periods. There is no evidence for a correlation between the TeV emission and
the mass-loss rate of the Be star but this may be affected by the strong, short
timescale (as short as intra-day) variation displayed by the H$\alpha$ fluxes.
-
Context. In the last five years the Fermi Large Area Telescope (LAT)
instrument detected GeV {\gamma}-ray emission from five novae. The GeV emission
can be interpreted in terms of an inverse Compton process of electrons
accelerated in a shock. In this case it is expected that protons in the same
conditions can be accelerated to much higher energies. Consequently they may
produce a second component in the {\gamma}-ray spectrum at TeV energies. Aims.
We aim to explore the very-high-energy domain to search for {\gamma}-ray
emission above 50 GeV and to shed light on the acceleration process of leptons
and hadrons in nova explosions. Methods. We have performed observations with
the MAGIC telescopes of the classical nova V339 Del shortly after the 2013
outburst, triggered by optical and subsequent GeV {\gamma}-ray detec- tions. We
also briefly report on VHE observations of the symbiotic nova YY Her and the
dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the
analysis of con- temporaneous Fermi-LAT data of the sources. The TeV and GeV
observations are compared in order to evaluate the acceleration parameters for
leptons and hadrons. Results. No significant TeV emission was found from the
studied sources. We computed upper limits on the spectrum and night-by-night
flux. The combined GeV and TeV observations of V339 Del limit the ratio of
proton to electron luminosities to Lp<~0.15 Le.
-
The high-frequency-peaked BL Lac (HBL) 1ES 0806+524 (z = 0.138) was
discovered in VHE $\gamma$ rays in 2008. Until now, the broad-band spectrum of
1ES 0806+524 has been only poorly characterized, in particular at high
energies. We analysed multiwavelength observations from $\gamma$ rays to radio
performed from 2011 January to March, which were triggered by the high activity
detected at optical frequencies. These observations constitute the most precise
determination of the broad-band emission of 1ES 0806+524 to date. The
stereoscopic MAGIC observations yielded a $\gamma$-ray signal above 250 GeV of
$(3.7 \pm 0.7)$ per cent of the Crab Nebula flux with a statistical
significance of 9.9 $\sigma$. The multiwavelength observations showed
significant variability in essentially all energy bands, including a VHE
$\gamma$-ray flare that lasted less than one night, which provided
unprecedented evidence for short-term variability in 1ES 0806+524. The spectrum
of this flare is well described by a power law with a photon index of $2.97 \pm
0.29$ between $\sim$150 GeV and 1 TeV and an integral flux of $(9.3 \pm 1.9)$
per cent of the Crab Nebula flux above 250 GeV. The spectrum during the
non-flaring VHE activity is compatible with the only available VHE observation
performed in 2008 with VERITAS when the source was in a low optical state. The
broad-band spectral energy distribution can be described with a one-zone
Synchrotron Self Compton model with parameters typical for HBLs, indicating
that 1ES 0806+524 is not substantially different from the HBLs previously
detected.