A radio-map of the colliding winds in the very massive binary system HD 93129A

Radio observations are an effective tool to discover particle acceleration regions in colliding-wind binaries, through detection of synchrotron radiation; these regions are natural laboratories for the study of relativistic particles. Wind-collision region (WCR) models can reproduce the radio continuum spectra of massive binaries that contain both thermal and non-thermal radio emission; however, key constraints for models come from high-resolution imaging. Only five WCRs have been resolved to date at radio frequencies at milliarcsec (mas) angular scales. The source HD 93129A, prototype of the very few known O2 I stars, is a promising target for study: recently, a second massive, early-type star about 50 mas away was discovered, and a non-thermal radio source detected in the region. Preliminary long-baseline array data suggest that a significant fraction of the radio emission from the system comes from a putative WCR. We sought evidence that HD 93129A is a massive binary system with colliding stellar winds that produce non-thermal radiation, through spatially resolved images of the radio emitting regions. We completed observations with the Australian Long Baseline Array (LBA) to resolve the system at mas angular resolutions and reduced archival Australia Telescope Compact Array (ATCA) data to derive the total radio emission. We also compiled optical astrometric data of the system in a homogeneous way. We reduced historical Hubble Space Telescope data and obtained absolute and relative astrometry with milliarcsec accuracy. The astrometric analysis leads us to conclude that the two stars in HD 93129A form a gravitationally bound system. The LBA data reveal an extended arc-shaped non-thermal source between the two stars, indicative of a WCR. The wind momentum-rate ratio of the two stellar winds is estimated. The ATCA data show a point source with a change in flux level ...Comment: Accepted in Astronomy and Astrophysic

Detection of bridge emission above 50 GeV from the Crab pulsar with the MAGIC telescopes

The Crab pulsar is the only astronomical pulsed source detected at very high energy (VHE, E>100GeV) gamma rays. The emission mechanism of VHE pulsation is not yet fully understood, although several theoretical models have been proposed. Aims. In order to test new models, we measured the light curve and the spectra of the Crab pulsar with high precision by means of deep observations. Methods. We analyzed 135 hours of selected MAGIC data taken between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower energies, 5.5 years of Fermi-LAT data were also analyzed. Results. The known two pulses per period were detected with a significance of 8.0 σ and 12.6 σ. In addition, significant emission was found between the two pulses with 6.2 σ. Conclusions. We discovered the bridge emission above 50 GeV between the two main pulses. This emission can not be explained with the existing theories. These data can be used for testing new theoretical models

Very High Energy γ-Rays from the Universe's Middle Age: Detection of the z = 0.940 Blazar PKS 1441+25 with MAGIC

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σ 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 timescale 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

The major upgrade of the MAGIC telescopes, Part I: The hardware improvements and the commissioning of the system

The MAGIC telescopes are two Imaging Atmospheric Cherenkov Telescopes (IACTs) located on the Canary island of La Palma. The telescopes are designed to measure Cherenkov light from air showers initiated by gamma rays in the energy regime from around 50 GeV to more than 50 TeV. The two telescopes were built in 2004 and 2009, respectively, with different cameras, triggers and readout systems. In the years 2011-2012 the MAGIC collaboration undertook a major upgrade to make the stereoscopic system uniform, improving its overall performance and easing its maintenance. In particular, the camera, the receivers and the trigger of the first telescope were replaced and the readout of the two telescopes was upgraded. This paper (Part I) describes the details of the upgrade as well as the basic performance parameters of MAGIC such as raw data treatment, linearity in the electronic chain and sources of noise. In Part II, we describe the physics performance of the upgraded system

Teraelectronvolt pulsed emission from the Crab Pulsar detected by MAGIC

Aims: We investigate the extension of the very high-energy spectral tail of the Crab Pulsar at energies above 400 GeV. Methods: We analyzed ~320 h of good-quality Crab data obtained with the MAGIC telescope from February 2007 to 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 those 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 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 ± 0.1. Conclusions: Using data from the MAGIC telescopes we measured 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 × 106. 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

MAGIC observations of MWC 656, the only known Be/BH system

Context. MWC 656 has recently been established as the first observationally detected high-mass X-ray binary system containing a Be star and a black hole (BH). The system has been associated with a gamma-ray flaring event detected by the AGILE satellite in July 2010. Aims: Our aim is to evaluate whether the MWC 656 gamma-ray emission extends to very high energy (VHE > 100 GeV) gamma rays. Methods: We observed MWC 656 with the MAGIC telescopes for ~23 h during two observation periods: between May and June 2012, and in June 2013. During the last period, observations were performed contemporaneously with X-ray (XMM-Newton) and optical (STELLA) instruments. Results: We did not detect the MWC 656 binary system at TeV energies with the MAGIC telescopes in either of the two campaigns. Upper limits (ULs) to the integral flux above 300 GeV have been set, as well as differential ULs at a level of ~5% of the Crab nebula flux. The results obtained from the MAGIC observations do not support persistent emission of VHE gamma rays from this system at a level of 2.4% the Crab flux

Multiwavelength Study of Quiescent States of Mrk 421 with Unprecedented Hard X-Ray Coverage Provided by NuSTAR in 2013

We present coordinated multiwavelength observations of the bright, nearby BL Lacertae object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very high energy (VHE) γ-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep ${\rm{\Gamma }}\approx 3$ power law, with no evidence for an exponential cutoff or additional hard components up to ~80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure that relates to the two bumps of the broadband SED. In each bump, the variability increases with energy, which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in situ electron acceleration and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states

A Search for Spectral Hysteresis and Energy-dependent Time Lags from X-Ray and TeV Gamma-Ray Observations of Mrk 421

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three 'target-of-opportunity' observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering the X-ray and optical/ultraviolet bands) and VERITAS (covering the TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-Large Area Telescope) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g., the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at gsim4 × 10−4 Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies

Constraints on the persistent radio source associated with FRB 20190520B using the European VLBI Network

We present very-long-baseline interferometry (VLBI) observations of a continuum radio source potentially associated with the fast radio burst source FRB 20190520B. Using the European VLBI network (EVN), we find the source to be compact on VLBI scales with an angular size of $<2.3$ mas ($3\sigma$). This corresponds to a transverse physical size of $<9$ pc (at the $z=0.241$ redshift of the host galaxy), confirming it to be an FRB persistent radio source (PRS) like that associated with the first-known repeater FRB 20121102A. The PRS has a flux density of $201 \pm 34 \rm{\mu Jy}$ at 1.7 GHz and a spectral radio luminosity of $L_{1.7 \rm GHz} = (3.0 \pm 0.5) \times 10^{29}\,\mathrm{erg s^{-1} Hz^{-1}}$ (also similar to the FRB 20121102A PRS). Comparing to previous lower-resolution observations, we find that no flux is resolved out on milliarcsecond scales. We have refined the PRS position, improving its precision by an order of magnitude compared to previous results. We also report the detection of a FRB 20190520B burst at 1.4 GHz and find the burst position to be consistent with the PRS position, at $\lesssim20$ mas. This strongly supports their direct physical association and the hypothesis that a single central engine powers both the bursts and the PRS. We discuss the model of a magnetar in a wind nebula and present an allowed parameter space for its age and the radius of the putative nebula powering the observed PRS emission. Alternatively, we find that an accretion-powered 'hypernebula' model also fits our observational constraints.Comment: 13 pages, 6 figures, Submitted to ApJ

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

Context. Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy distribution (SED). Aims: The goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail. Methods: A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E> 100 GeV) γ-rays with MAGIC, VERITAS, Whipple, Fermi-LAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population. Results: The activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2-10 keV band and the VHE γ-ray flux above 200 GeV, consistent with the γ-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE γ-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 × 104<γ< 6 × 105, where γ is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands. Conclusions: Both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (δ, B and R), in addition to the parameters related to the particle population. Appendices are available in electronic form at http://www.aanda.orgMulti-wavelength light curves (data in Fig. 1) and broadband spectral energy distributions (the data in Figs. 7, 8a-9f, 12a-13f) are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/578/A2