Showing 1–50 of 206 results for author: Colin, P

Study of the variable broadband emission of Markarian 501 during the most extreme Swift X-ray activity

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, G. Bonnoli, G. Busetto, R. Carosi, G. Ceribella, S. Cikota, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras , et al. (163 additional authors not shown)

Abstract: Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios… ▽ More Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$σ$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$σ$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$σ$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade. △ Less

Submitted 16 March, 2021; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: Accepted for publication in A&A. Corresponding authors: J. Becerra González (jbecerra@iac.es), D. Paneque (dpaneque@mppmu.mpg.de), C. Wendel (cwendel@astro.uni-wuerzburg.de)

Journal ref: A&A 637, A86 (2020)

Statistics of VHE gamma-Rays in Temporal Association with Radio Giant Pulses from the Crab Pulsar

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, B. Biasuzzi, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin, E. Colombo , et al. (124 additional authors not shown)

Abstract: The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E>100 GeV) gamma-rays from the Crab pulsar. 16 hours of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV with the Major Atmospheric Gamma-r… ▽ More The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E>100 GeV) gamma-rays from the Crab pulsar. 16 hours of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes were performed. We searched for a statistical correlation between the radio and VHE gamma-ray emission with search windows of different lengths and different time lags to the arrival times of a radio GP. A dedicated search for an enhancement in the number of VHE gamma-rays correlated with the occurrence of radio GPs was carried out separately for the P1 and P2 phase ranges respectively. 99444 radio GPs have been detected in the radio data sample. We find no significant correlation between the GPs and VHE photons in any of the search windows. Depending on phase cuts and the chosen search windows we find upper limits at 95% confidence level on an increase in VHE gamma-ray events correlated with radio GPs between 7% and 61% of the average Crab pulsar VHE flux for the P1 and P2 phase ranges respectively. This puts upper limits on the flux increase during a radio GP of 12% to 2900% (depending on search window duration and phase cuts) of the pulsed VHE flux. This is the most stringent upper limit on a correlation between gamma-ray emission and radio GPs reported so far.} △ Less

Submitted 31 January, 2020; v1 submitted 1 November, 2019; originally announced November 2019.

Comments: 15 pages, 7 figures, accepted for publication by Astronomy and Astrophysics

Journal ref: A&A 634, A25 (2020)

Galactic Cosmic Ray Sun Shadow during the declining phase of cycle 24 observed by HAWC

Authors: Alejandro Lara, Paulina Colin, K. P. Arunbabu, James Ryan

Abstract: The High Altitude Water Cherenkov (HAWC) array is sensitive to high energy Cosmic Rays (CR) in the $\sim 10$ to $\sim 200$ TeV energy range, making it possible to construct maps of the so called "Sun Shadow" ($SS$), i. e. of the deficit of CR coming from the direction of the Sun. In this work, we present the variation of the Relative Intensity of the deficit ($SS_{RI}$) for three years of HAWC obs… ▽ More The High Altitude Water Cherenkov (HAWC) array is sensitive to high energy Cosmic Rays (CR) in the $\sim 10$ to $\sim 200$ TeV energy range, making it possible to construct maps of the so called "Sun Shadow" ($SS$), i. e. of the deficit of CR coming from the direction of the Sun. In this work, we present the variation of the Relative Intensity of the deficit ($SS_{RI}$) for three years of HAWC observations form 2016 to 2018 in which we found a clear decreasing trend of the ($SS_{RI}$) over the studied period, corresponding to the declining phase of the solar cycle 24. By comparing the $SS_{RI}$ with the photospheric magnetic field evolution, we show that there is a linear relationship between the $SS_{RI}$ and the median photospheric magnetic field of the Active Region belt (-40$^\circ \le$ lat $\le$ 40$^\circ$) and a inverse linear relationship with the polar photospheric magnetic field (lat $\ge \pm$ 60$^\circ$). The former relationship is due to the magnetic field causing a deviation of the CR, whereas the latter reflects the change of the heliospheric field topology from multipolar to dipolar configurations. These relationships are valid only when the median magnetic field is lower than 8 G, during the declining and minimum phases of the solar cycle 24. Finally, we show that relativistic charged particles, in the 10 to 200 TeV energy range, are deflected a few degrees. △ Less

Submitted 20 August, 2019; originally announced August 2019.

Journal ref: PoS(ICRC2019)1104

Deep observations of the globular cluster M15 with the MAGIC telescopes

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, G. Busetto, R. Carosi, G. Ceribella, S. Cikota, S. M. Colak , et al. (131 additional authors not shown)

Abstract: A population of globular clusters (GCs) has been recently established by the Fermi-LAT telescope as a new class of GeV $γ$-ray sources. Leptons accelerated to TeV energies, in the inner magnetospheres of MSPs or in their wind regions, should produce $γ$-rays through the inverse Compton scattering in the dense radiation field from the huge population of stars. We have conducted deep observations of… ▽ More A population of globular clusters (GCs) has been recently established by the Fermi-LAT telescope as a new class of GeV $γ$-ray sources. Leptons accelerated to TeV energies, in the inner magnetospheres of MSPs or in their wind regions, should produce $γ$-rays through the inverse Compton scattering in the dense radiation field from the huge population of stars. We have conducted deep observations of the globular cluster M15 with the MAGIC telescopes and used 165 hrs in order to search for $γ$-ray emission. A strong upper limit on the TeV $γ$-ray flux $<3.2\times 10^{-13}\mathrm{cm^{-2}s^{-1}}$ above 300 GeV ($<0.26\%$ of the Crab nebula flux) has been obtained. We interpret this limit as a constraint on the efficiency of the acceleration of leptons in the magnetospheres of the MSPs. We constrain the injection rate of relativistic leptons, $η_{\rm e}$, from the MSPs magnetospheres and their surrounding. We conclude that $η_{\rm e}$ must be lower than expected from the modelling of high energy processes in MSP inner magnetospheres. For leptons accelerated with the power law spectrum in the MSP wind regions, $η_{\rm e}$ is constrained to be much lower than derived for the wind regions around classical pulsars. These constraints are valid for the expected range of magnetic field strengths within the GC and for the range of likely energies of leptons injected from the inner magnetospheres, provided that the leptons are not removed from the globular cluster very efficiently due to advection process. We discuss consequences of these constraints for the models of radiation processes around millisecond pulsars. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 484, Issue 2, 1 April 2019, Pages 2876-2885

MAGIC and Fermi-LAT gamma-ray results on unassociated HAWC sources

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak, P. Colin , et al. (318 additional authors not shown)

Abstract: The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the fir… ▽ More The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail. △ Less

Submitted 13 January, 2019; originally announced January 2019.

Comments: 12 pages, 3 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 485, Issue 1, May 2019, Pages 356-366

A Fast Very High Energy $γ$-ray Flare from BL Lacertae during a Period of Multiwavelength activity in June 2015

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, S. Cikota, S. M. Colak , et al. (141 additional authors not shown)

Abstract: The mechanisms producing fast variability of the $γ$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $γ$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a… ▽ More The mechanisms producing fast variability of the $γ$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $γ$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $γ$-ray bands. In this paper we present the MAGIC VHE $γ$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $γ$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $γ$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $γ$-rays is very similar to two other observed VHE $γ$-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43\,GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution, based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the $γ$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $γ$-ray flare, but all have some tension with the multiwavelength observations. △ Less

Submitted 7 January, 2019; originally announced January 2019.

Comments: Accepted to Astronomy & Astrophysics

Journal ref: A&A 623, A175 (2019)

Discovery of TeV $γ$-ray emission from the neighbourhood of the supernova remnant G24.7+0.6 by MAGIC

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee , et al. (129 additional authors not shown)

Abstract: SNR G24.7+0.6 is a 9.5 kyrs radio and $γ$-ray supernova remnant evolving in a dense medium. In the GeV regime, SNR G24.7+0.6 (3FHL\,J1834.1--0706e/FGES\,J1834.1--0706) shows a hard spectral index ($Γ$$\sim$2) up to $200$\,GeV, which makes it a good candidate to be observed with Cherenkov telescopes such as MAGIC. We observed the field of view of \snr\ with the MAGIC telescopes for a total of 31 ho… ▽ More SNR G24.7+0.6 is a 9.5 kyrs radio and $γ$-ray supernova remnant evolving in a dense medium. In the GeV regime, SNR G24.7+0.6 (3FHL\,J1834.1--0706e/FGES\,J1834.1--0706) shows a hard spectral index ($Γ$$\sim$2) up to $200$\,GeV, which makes it a good candidate to be observed with Cherenkov telescopes such as MAGIC. We observed the field of view of \snr\ with the MAGIC telescopes for a total of 31 hours. We detect very high energy $γ$-ray emission from an extended source located 0.34\degr\ away from the center of the radio SNR. The new source, named \mgc\ is detected up to 5\,TeV, and its spectrum is well-represented by a power-law function with spectral index of $2.74 \pm 0.08$. The complexity of the region makes the identification of the origin of the very-high energy emission difficult, however the spectral agreement with the LAT source and overlapping position at less than 1.5$σ$ point to a common origin. We analysed 8 years of \fermi-LAT data to extend the spectrum of the source down to 60\,MeV. \fermi-LAT and MAGIC spectra overlap within errors and the global broad band spectrum is described by a power-law with exponential cutoff at $1.9\pm0.5$\,TeV. The detected $γ$-ray emission can be interpreted as the results of proton-proton interaction between the supernova and the CO-rich surrounding. △ Less

Submitted 12 December, 2018; originally announced December 2018.

Comments: accepted for publication by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 483, Issue 4, March 2019, Pages 4578-4585

Periastron Observations of TeV Gamma-Ray Emission from a Binary System with a 50-year Period

Authors: The VERITAS Collaboration, A. U. Abeysekara, W. Benbow, R. Bird, A. Brill, R. Brose, J. H. Buckley, A. J. Chromey, M. K. Daniel, A. Falcone, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, P. Kaaret, P. Kar, M. Kertzman, D. Kieda , et al. (191 additional authors not shown)

Abstract: We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric 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 mont… ▽ More We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric 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 0.5 TeV. 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 (XRT) on board the Neil Gehrels \textit{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 we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Comments: 12 pages, 3 figures, 1 table. Accepted for publication in ApJL

Journal ref: The Astrophysical Journal Letters, 867:L19 (8pp), 2018 November 1

Constraining very-high-energy and optical emission from FRB 121102 with the MAGIC telescopes

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee , et al. (133 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous o… ▽ More Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$σ$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $σ$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: 9 pages, 3 figures. Accepted for publication in MNRAS the 2nd September 2018

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 481, Issue 2, 1 December 2018, Pages 2479-2486

The extreme HBL behaviour of Markarian 501 during 2012

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras, J. Cortina, S. Covino , et al. (254 additional authors not shown)

Abstract: A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0… ▽ More A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays. △ Less

Submitted 14 August, 2018; v1 submitted 13 August, 2018; originally announced August 2018.

Comments: 25 pages, 9 figures, accepted for publication in A&A. Corresponding authors: Gareth Hughes (gareth.hughes@cfa.harvard.edu), David Paneque (dpaneque@mppmu.mpg.de), Amit Shukla (amit.shukla@astro.uni-wuerzburg.de)

Journal ref: A&A 620, A181 (2018)

The blazar TXS 0506+056 associated with a high-energy neutrino: insights into extragalactic jets and cosmic ray acceleration

Authors: S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak, P. Colin , et al. (124 additional authors not shown)

Abstract: A neutrino with energy of $\sim$290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS~0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at $\sim 3σ$ level. We monitored the object in the very-high-energy (VHE) band with the MAGIC telescopes for $\sim$41 hours from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability i… ▽ More A neutrino with energy of $\sim$290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS~0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at $\sim 3σ$ level. We monitored the object in the very-high-energy (VHE) band with the MAGIC telescopes for $\sim$41 hours from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of $\sim 10^{45} - 4 \times 10^{46} {\rm erg \ s^{-1}}$. The steep spectrum observed by MAGIC is concordant with internal $γγ$ absorption above a few tens of GeV entailed by photohadronic production of a $\sim$290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton up-scattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet co-moving frame can be in the range $\sim 10^{14}$ to $10^{18}$ eV. △ Less

Submitted 30 July, 2018; v1 submitted 11 July, 2018; originally announced July 2018.

Comments: 11 pages, 3 figues, 3 tables, Accepted for publication in ApJL. Major changes in the revised Version: corrected author list

Journal ref: The Astrophysical Journal Letters, 863 (2018) L10

The broad-band properties of the intermediate synchrotron peaked BL Lac S2 0109+22 from radio to VHE gamma rays

Authors: MAGIC Collaboration, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak , et al. (141 additional authors not shown)

Abstract: The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classifi… ▽ More The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classification. Finally, we compare the source long-term behaviour to other VHE gamma-ray emitting (TeV) blazars. We performed a temporal and spectral analysis of the data centred around the MAGIC interval of observation (MJD 57225-57231). Long-term radio and optical data have also been investigated using the discrete correlation function. The redshift of the source is estimated through optical host-galaxy imaging and also using the amount of VHE gamma-ray absorption. The quasi-simultaneous multi-frequency spectral energy distribution (SED) is modelled with the conventional one-zone synchrotron self-Compton (SSC) model. MAGIC observations resulted in the detection of the source at a significance level of $5.3\,σ$. The VHE gamma-ray emission of S2 0109+22 is variable on a daily time scale. VHE gamma-ray luminosity of the source is lower than the average of TeV BL Lacs. The optical polarization, and long-term optical/radio behaviour of the source are different from the general population of TeV blazars. All these findings agree with the classification of the source as an ISP BL Lac object. We estimate the source redshift as $z = 0.36 \pm 0.07$. The SSC parameters describing the SED are rather typical for blazars. △ Less

Submitted 5 July, 2018; originally announced July 2018.

Comments: 15 pages, 10 figures, Corresponding authors: Fallah Ramazani, V. (vafara@utu.fi), Hovatta, T. (talvikki.hovatta@utu.fi), Lindfors, E. (elilin@utu.fi) and Nilsson, K. (kani@utu.fi)

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 480, Issue 1, 11 October 2018, Pages 879-892

Multi-wavelength characterization of the blazar S5~0716+714 during an unprecedented outburst phase

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak , et al. (165 additional authors not shown)

Abstract: 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. The comprehensive dataset collected is investigated in order to shed light… ▽ More 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. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. 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 VHE. In the constructed SED the \textit{Swift}-XRT+\textit{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 \textit{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 EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out 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 modelling. △ Less

Submitted 1 July, 2018; originally announced July 2018.

Journal ref: Astronomy & Astrophysics 619, A45 (2018)

Constraining Dark Matter lifetime with a deep gamma-ray survey of the Perseus Galaxy Cluster with MAGIC

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, S. Cikota, S. M. Colak , et al. (126 additional authors not shown)

Abstract: Clusters of galaxies are the largest known gravitationally bound structures in the Universe, with masses around $10^{15}\ M_\odot$, most of it in the form of dark matter. The ground-based Imaging Atmospheric Cherenkov Telescope MAGIC made a deep survey of the Perseus cluster of galaxies using almost 400 h of data recorded between 2009 and 2017. This is the deepest observational campaign so far on… ▽ More Clusters of galaxies are the largest known gravitationally bound structures in the Universe, with masses around $10^{15}\ M_\odot$, most of it in the form of dark matter. The ground-based Imaging Atmospheric Cherenkov Telescope MAGIC made a deep survey of the Perseus cluster of galaxies using almost 400 h of data recorded between 2009 and 2017. This is the deepest observational campaign so far on a cluster of galaxies in the very high energy range. We search for gamma-ray signals from dark matter particles in the mass range between 200 GeV and 200 TeV decaying into standard model pairs. We apply an analysis optimized for the spectral and morphological features expected from dark matter decays and find no evidence of decaying dark matter. From this, we conclude that dark matter particles have a decay lifetime longer than $\sim10^{26}$~s in all considered channels. Our results improve previous lower limits found by MAGIC and represent the strongest limits on decaying dark matter particles from ground-based gamma-ray instruments. △ Less

Submitted 28 August, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: 5 figures, corresponding authors: Joaquim Palacio (jpalacio@ifae.es), Michele Doro (michele.doro@unipd.it), Mónica Vazquez-Acosta (monicava@iac.es) and Javier Rico (jrico@ifae.es)

Journal ref: Physics of the Dark Universe Volume 22, December 2018, Pages 38-47

Detection of persistent VHE gamma-ray emission from PKS 1510-089 by the MAGIC telescopes during low states between 2012 and 2017

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee , et al. (143 additional authors not shown)

Abstract: PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize… ▽ More PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modelled within an External Compton scenario with a stationary emission region through which plasma and magnetic field are flowing. We perform also the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hrs of data during times when the Fermi-LAT flux measured above 1 GeV was below 3x10-8cm-2s-1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5σ) VHE gamma-ray emission at the level of (4.27+-0.61stat)x10-12cm-2s-1 above 150GeV, a factor 80 smaller than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broad-band emission is compatible with the EC scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% C.L. on the location of the emission region allows us to place it at the distance >74% of the outer radius of the BLR. △ Less

Submitted 3 September, 2018; v1 submitted 14 June, 2018; originally announced June 2018.

Comments: 12 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 619, A159 (2018)

Gamma-ray flaring activity of NGC 1275 in 2016-2017 measured by MAGIC

Authors: MAGIC collaboration, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak, P. Colin , et al. (130 additional authors not shown)

Abstract: We report on the detection of flaring activity from the Fanaroff-Riley I radio galaxy NGC 1275 in very-high-energy (VHE, E $>$ 100 GeV) gamma rays with the MAGIC telescopes. Observations were performed between 2016 September and 2017 February as part of a monitoring program. The brightest outburst with $\sim1.5$ times the Crab Nebula flux above 100 GeV (C.U.) was observed during the night between… ▽ More We report on the detection of flaring activity from the Fanaroff-Riley I radio galaxy NGC 1275 in very-high-energy (VHE, E $>$ 100 GeV) gamma rays with the MAGIC telescopes. Observations were performed between 2016 September and 2017 February as part of a monitoring program. The brightest outburst with $\sim1.5$ times the Crab Nebula flux above 100 GeV (C.U.) was observed during the night between 2016 December 31 and 2017 January 1 (fifty times higher than the mean previously measured in two observational campaigns between 2009 and 2011). Significant variability of the day-by-day light curve was measured, the shortest flux-doubling time-scales was found to be of $(611\pm101)$ min. The combined spectrum of the MAGIC data during the strongest flare state and simultaneous data from the Fermi-LAT around 2017 January 1 follows a power-law with an exponential cutoff at the energy $(492\pm35)$ GeV. Simultaneous optical flux density measurements in the R-band obtained with the KVA telescope are also presented and the correlation between the optical and gamma-ray emission is investigated. Due to possible internal pair-production, the fast flux variability constrains the Doppler factor to values which are inconsistent with a large viewing angle as observed in the radio band. We investigate different scenarios for the explanation of fast gamma-ray variability, namely emission from: magnetospheric gaps, relativistic blobs propagating in the jet (mini-jets) or external cloud (or star) entering the jet. We find that the only plausible model to account for the luminosities here observed would be the production of gamma rays in a magnetospheric gap around the central black hole only in the eventuality of an enhancement of the magnetic field threading the hole from its equipartition value with the gas pressure in the accretion flow. △ Less

Submitted 5 June, 2018; originally announced June 2018.

Comments: accepted for publications in Astronomy and Astrophysics

Journal ref: A&A 617, A91 (2018)

Limits on the flux of tau neutrinos from 1 PeV to 3 EeV with the MAGIC telescopes

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak , et al. (120 additional authors not shown)

Abstract: A search for tau neutrino induced showers with the MAGIC telescopes is presented. The MAGIC telescopes located at an altitude of 2200 m a.s.l. in the Canary Island of La Palma, can point towards the horizon or a few degrees below across an azimuthal range of about 80 degrees. This provides a possibility to search for air showers induced by tau leptons arising from interactions of tau neutrinos in… ▽ More A search for tau neutrino induced showers with the MAGIC telescopes is presented. The MAGIC telescopes located at an altitude of 2200 m a.s.l. in the Canary Island of La Palma, can point towards the horizon or a few degrees below across an azimuthal range of about 80 degrees. This provides a possibility to search for air showers induced by tau leptons arising from interactions of tau neutrinos in the Earth crust or the surrounding ocean. In this paper we show how such air showers can be discriminated from the background of very inclined hadronic showers by using Monte Carlo simulations. Taking into account the orography of the site, the point source acceptance and the event rates expected have been calculated for a sample of generic neutrino fluxes from photo-hadronic interactions in AGNs. The analysis of about 30 hours of data taken towards the sealeads to a 90\% C.L. point source limit for tau neutrinos in the energy range from $1.0 \times 10^{15}$ eV to $3.0 \times 10^{18}$ eV of about $E_{ν_τ}^{2}\times φ(E_{ν_τ}) < 2.0 \times 10^{-4}$ GeV cm$^{-2}$ s$^{-1}$ for an assumed power-law neutrino spectrum with spectral index $γ$=-2. However, with 300 hours and in case of an optimistic neutrino flare model, limits of the level down to $E_{ν_τ}^{2}\times φ(E_{ν_τ}) < 8.4 \times 10^{-6}$ GeV cm$^{-2}$ s$^{-1}$ can be expected. △ Less

Submitted 7 May, 2018; originally announced May 2018.

Comments: accepted for publication in Astroparticle Physics

Journal ref: Astroparticle Physics Volume 102, November 2018, Pages 77-88

The detection of the blazar S4 0954+65 at very-high-energy with the MAGIC telescopes during an exceptionally high optical state

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak, P. Colin , et al. (147 additional authors not shown)

Abstract: The very-high-energy (VHE, $\gtrsim 100$ GeV) $γ$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physi… ▽ More The very-high-energy (VHE, $\gtrsim 100$ GeV) $γ$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $γ$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $γ$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL). △ Less

Submitted 12 January, 2018; originally announced January 2018.

Comments: 15 pages, 7 figures, submitted to A&A

Journal ref: A&A 617, A30 (2018)

The global and radial stellar mass assembly of Milky Way-sized galaxies

Authors: V. Avila-Reese, A. González-Samaniego, P. Colín, H. Ibarra-Medel, A. Rodríguez-Puebla

Abstract: We study the global and radial stellar mass assembly of eight zoomed-in MW-sized galaxies produced in Hydrodynamics cosmological simulations. The disk-dominated galaxies (4) show a fast initial stellar mass growth in the innermost parts, driven mostly by in-situ SF, but since $z\sim2-1$ the SF enters in a long-term quenching phase. The outer regions follow this trend but more gentle as more extern… ▽ More We study the global and radial stellar mass assembly of eight zoomed-in MW-sized galaxies produced in Hydrodynamics cosmological simulations. The disk-dominated galaxies (4) show a fast initial stellar mass growth in the innermost parts, driven mostly by in-situ SF, but since $z\sim2-1$ the SF enters in a long-term quenching phase. The outer regions follow this trend but more gentle as more external they are. As the result, the radial stellar mass growth is highly inside-out due to both the inside-out structural growth and inside-out SF quenching. The half-mass radius evolves fast; for instance, $R_{0.5}$($z=1$)$<0.5$$R_{0.5}$($z=0$). Two other runs resemble lenticular galaxies. One shows also a pronounced inside-out growth and the other one presents a nearly uniform radial mass assembly. The other two galaxies suffered late major mergers. Their normalized radial mass growth histories (MGHs) are nearly close among them but with periods of outside-in assembly during or after the mergers. For all the simulations, the archaeological radial MGHs calculated from the $z = 0$ stellar-particles age distribution are similar to the current MGHs, which evidences that the mass assembly by ex-situ stars and the radial mass transport do not change significantly their radial mass distributions. Our results agree qualitatively with observational inferences from the fossil record method applied to a survey of local galaxies and from look-back observations of progenitors of MW-sized galaxies. However, the inside-out growth mode is more pronounced and the $R_{0.5}$ growth is faster in simulations than in observational inferences. △ Less

Submitted 31 January, 2018; v1 submitted 5 January, 2018; originally announced January 2018.

Comments: 23 pages, 15 figures. Accepted for publication in ApJ. Minor corrections and more references added. The discussion was extended in section 6.1

Indirect dark matter searches in the dwarf satellite galaxy Ursa Major II with the MAGIC Telescopes

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak , et al. (120 additional authors not shown)

Abstract: The dwarf spheroidal galaxy Ursa Major II (UMaII) is believed to be one of the most dark-matter dominated systems among the Milky Way satellites and represents a suitable target for indirect dark matter (DM) searches. The MAGIC telescopes carried out a deep observation campaign on UMaII between 2014 and 2016, collecting almost one hundred hours of good-quality data. This campaign enlarges the pool… ▽ More The dwarf spheroidal galaxy Ursa Major II (UMaII) is believed to be one of the most dark-matter dominated systems among the Milky Way satellites and represents a suitable target for indirect dark matter (DM) searches. The MAGIC telescopes carried out a deep observation campaign on UMaII between 2014 and 2016, collecting almost one hundred hours of good-quality data. This campaign enlarges the pool of DM targets observed at very high energy (E$\gtrsim$50GeV) in search for signatures of dark matter annihilation in the wide mass range between $\sim$100 GeV and $\sim$100 TeV. To this end, the data are analyzed with the full likelihood analysis, a method based on the exploitation of the spectral information of the recorded events for an optimal sensitivity to the explored dark matter models. We obtain constraints on the annihilation cross-section for different channels that are among the most robust and stringent achieved so far at the TeV mass scale from observations of dwarf satellite galaxies. △ Less

Submitted 9 March, 2018; v1 submitted 7 December, 2017; originally announced December 2017.

Comments: 18 pages, 8 figures

Journal ref: JCAP03(2018)009

Science with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, B. S. Acharya, I. Agudo, I. Al Samarai, R. Alfaro, J. Alfaro, C. Alispach, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, E. Antolini, L. A. Antonelli, C. Aramo, M. Araya, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, M. Ashley, M. Backes, C. Balazs, M. Balbo, O. Ballester , et al. (558 additional authors not shown)

Abstract: The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black ho… ▽ More The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments. The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources. The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document. △ Less

Submitted 21 January, 2018; v1 submitted 22 September, 2017; originally announced September 2017.

Comments: 213 pages, including references and glossary. Version 2: credits and references updated, some figures updated, and author list updated

Cherenkov Telescope Array Contributions to the 35th International Cosmic Ray Conference (ICRC2017)

Authors: F. Acero, B. S. Acharya, V. Acín Portella, C. Adams, I. Agudo, F. Aharonian, I. Al Samarai, A. Alberdi, M. Alcubierre, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, J. -P. Amans, E. Amato, L. Ambrogi, G. Ambrosi, M. Ambrosio, J. Anderson, M. Anduze, E. O. Angüner, E. Antolini, L. A. Antonelli, V. Antonuccio , et al. (1117 additional authors not shown)

Abstract: List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea. List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea. △ Less

Submitted 24 October, 2017; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: Index of Cherenkov Telescope Array conference proceedings at the ICRC2017, Busan, Korea

Constraining Lorentz invariance violation using the Crab Pulsar emission observed up to TeV energies by MAGIC

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin , et al. (120 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Journal ref: The Astrophysical Journal Supplement Series, 232:9, September 2017

Performance of the MAGIC telescopes under moonlight

Authors: D. Guberman, P. Colin

Abstract: MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we present a dedicated Moon-adapted analysis and characterize the performanc… ▽ More MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we present a dedicated Moon-adapted analysis and characterize the performance of MAGIC under moonlight. We evaluate energy threshold, angular resolution and sensitivity of MAGIC under different background light levels, based on Crab Nebula observations and tuned Monte Carlo simulations. This study includes observations taken under non-standard hardware configurations, such as reducing the camera photomultiplier tubes gain by a factor $\sim$1.7 (reduced HV settings) with respect to standard settings (nominal HV) or using UV-pass filters to strongly reduce the amount of moonlight reaching the telescopes cameras. The Crab Nebula spectrum is correctly reconstructed in all the studied illumination levels, that reach up to 30 times brighter than under dark conditions. The main effect of moonlight is an increase in the analysis energy threshold and in the systematic uncertainties on the flux normalization. The sensitivity degradation is constrained to be below 10\%, within 15-30\% and between 60 and 80\% for nominal HV, reduced HV and UV-pass filter observations, respectively. No worsening of the angular resolution was found. Thanks to observations during moonlight, the duty cycle can be doubled, suppressing the need to stop observations around full Moon. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: Proceedings of the 35th International Cosmic Ray Conference (ICRC 2017), Bexco, Busan, Korea (arXiv:1708.05153)

Report number: MAGIC-ICRC/2017/18

MAGIC Contributions to the 35th International Cosmic Ray Conference (ICRC2017)

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras , et al. (121 additional authors not shown)

Abstract: MAGIC (Major Atmospheric Gamma Imaging Cherenkov) is a system of two 17 m diameter, F/1.03 Imaging Atmospheric Cherenkov Telescopes (IACT). They are dedicated to the observation of gamma rays from galactic and extragalactic sources in the very high energy range (VHE, 30 GeV to 100 TeV). This submission contains links to the proceedings for the 35th International Cosmic Ray Conference (ICRC2017), h… ▽ More MAGIC (Major Atmospheric Gamma Imaging Cherenkov) is a system of two 17 m diameter, F/1.03 Imaging Atmospheric Cherenkov Telescopes (IACT). They are dedicated to the observation of gamma rays from galactic and extragalactic sources in the very high energy range (VHE, 30 GeV to 100 TeV). This submission contains links to the proceedings for the 35th International Cosmic Ray Conference (ICRC2017), held in Bexco, Busan, Korea from the 12th to the 17th of July, 2017. △ Less

Submitted 14 September, 2017; v1 submitted 17 August, 2017; originally announced August 2017.

Comments: HTML page with links

Search for very-high-energy gamma-ray emission from the microquasar Cygnus X-1 with the MAGIC telescopes

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, P. Colin, E. Colombo , et al. (122 additional authors not shown)

Abstract: 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 t… ▽ More 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 $Γ=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. △ Less

Submitted 11 August, 2017; originally announced August 2017.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 472, Issue 3, 11 December 2017, Pages 3474-3485

Constraints on particle acceleration in SS433/W50 from MAGIC and H.E.S.S. observations

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, R. Carosi, A. Carosi, A. Chatterjee, P. Colin, E. Colombo , et al. (386 additional authors not shown)

Abstract: 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 capabili… ▽ More 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$μ$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. △ Less

Submitted 12 July, 2017; originally announced July 2017.

Comments: Accepted for publication in A&A, 9 pages, 2 figures

Journal ref: A&A 612, A14 (2018)

A cut-off in the TeV gamma-ray spectrum of the SNR Cassiopeia A

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin , et al. (119 additional authors not shown)

Abstract: It is widely believed that the bulk of the Galactic cosmic rays are accelerated in supernova remnants (SNRs). However, no observational evidence of the presence of particles of PeV energies in SNRs has yet been found. The young historical SNR Cassiopeia A (Cas A) appears as one of the best candidates to study acceleration processes. Between December 2014 and October 2016 we observed Cas A with the… ▽ More It is widely believed that the bulk of the Galactic cosmic rays are accelerated in supernova remnants (SNRs). However, no observational evidence of the presence of particles of PeV energies in SNRs has yet been found. The young historical SNR Cassiopeia A (Cas A) appears as one of the best candidates to study acceleration processes. Between December 2014 and October 2016 we observed Cas A with the MAGIC telescopes, accumulating 158 hours of good-quality data. We derived the spectrum of the source from 100 GeV to 10 TeV. We also analysed $\sim$8 years of $Fermi$-LAT to obtain the spectral shape between 60 MeV and 500 GeV. The spectra measured by the LAT and MAGIC telescopes are compatible within the errors and show a clear turn off (4.6 $σ$) at the highest energies, which can be described with an exponential cut-off at $E_c = 3.5\left(^{+1.6}_{-1.0}\right)_{\textit{stat}} \left(^{+0.8}_{-0.9}\right)_{\textit{sys}}$ TeV. The gamma-ray emission from 60 MeV to 10 TeV can be attributed to a population of high-energy protons with spectral index $\sim$2.2 and energy cut-off at $\sim$10 TeV. This result indicates that Cas A is not contributing to the high energy ($\sim$PeV) cosmic-ray sea in a significant manner at the present moment. A one-zone leptonic model fails to reproduce by itself the multi-wavelength spectral energy distribution. Besides, if a non-negligible fraction of the flux seen by MAGIC is produced by leptons, the radiation should be emitted in a region with a low magnetic field (B$\lessapprox$100$μ$G) like in the reverse shock.) △ Less

Submitted 5 July, 2017; originally announced July 2017.

Comments: 8 pages, 3 figures, submitted to MNRAS. Corresponding authors: Daniel Guberman (dguberman@ifae.es), Emma de Oña Wilhelmi (wilhelmi@ice.csic.es) and Daniel Galindo (dgalindo@am.ub.es)

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 472, Issue 3, 11 December 2017, Pages 2956-2962

MAGIC observations of the microquasar V404 Cygni during the 2015 outburst

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, P. Colin, E. Colombo, J. L. Contreras, J. Cortina , et al. (120 additional authors not shown)

Abstract: The microquasar V404 Cygni underwent a series of outbursts in 2015, June 15-31, during which its flux in hard X-rays (20-40 keV) reached about 40 times the Crab Nebula flux. Because of the exceptional interest of the flaring activity from this source, observations at several wavelengths were conducted. The MAGIC telescopes, triggered by the INTEGRAL alerts, followed-up the flaring source for sever… ▽ More The microquasar V404 Cygni underwent a series of outbursts in 2015, June 15-31, during which its flux in hard X-rays (20-40 keV) reached about 40 times the Crab Nebula flux. Because of the exceptional interest of the flaring activity from this source, observations at several wavelengths were conducted. The MAGIC telescopes, triggered by the INTEGRAL alerts, followed-up the flaring source for several nights during the period June 18-27, for more than 10 hours. One hour of observation was conducted simultaneously to a giant 22 GHz radio flare and a hint of signal at GeV energies seen by Fermi-LAT. The MAGIC observations did not show significant emission in any of the analysed time intervals. The derived flux upper limit, in the energy range 200--1250 GeV, is 4.8$\times 10^{-12}$ ph cm$^{-2}$ s$^{-1}$. We estimate the gamma-ray opacity during the flaring period, which along with our non-detection, points to an inefficient acceleration in the V404\,Cyg jets if VHE emitter is located further than $1\times 10^{10}$ cm from the compact object. △ Less

Submitted 4 July, 2017; originally announced July 2017.

Comments: 7 pages, accepted by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 471, Issue 2, 21 October 2017, Pages 1688-1693

Observation of the Black Widow B1957+20 millisecond pulsar binary system with the MAGIC telescopes

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, R. Carosi, A. Carosi, A. Chatterjee, P. Colin , et al. (123 additional authors not shown)

Abstract: B1957+20 is a millisecond pulsar located in a black widow type compact binary system with a low mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high energy γ rays in the inverse Compton process. We performe… ▽ More B1957+20 is a millisecond pulsar located in a black widow type compact binary system with a low mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high energy γ rays in the inverse Compton process. We performed extensive observations with the MAGIC telescopes of B1957+20. We interpret results in the framework of a few different models, namely emission from the vicinity of the millisecond pulsar, the interaction of the pulsar and stellar companion wind region, or bow shock nebula. No significant steady very high energy γ-ray emission was found. We derived a 95% confidence level upper limit of 3.0 x 10 -12 cm -2 s -1 on the average γ-ray emission from the binary system above 200 GeV. The upper limits obtained with MAGIC constrain, for the first time, different models of the high-energy emission in B1957+20. In particular, in the inner mixed wind nebula model with mono-energetic injection of electrons, the acceleration efficiency of electrons is constrained to be below ~(2-10)% of the pulsar spin down power. For the pulsar emission, the obtained upper limits for each emission peak are well above the exponential cut-off fits to the Fermi-LAT data, extrapolated to energies above 50 GeV. The MAGIC upper limits can rule out a simple power-law tail extension through the sub-TeV energy range for the main peak seen at radio frequencies. △ Less

Submitted 5 June, 2017; originally announced June 2017.

Comments: 12 pages, 10 figures, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 470, Issue 4, October 2017, Pages 4608-4617

Performance of the MAGIC telescopes under moonlight

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, P. Colin, E. Colombo , et al. (121 additional authors not shown)

Abstract: MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we develop a dedicated Moon-adapted analysis to characterize the performance… ▽ More MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we develop a dedicated Moon-adapted analysis to characterize the performance of MAGIC under moonlight. We evaluate energy threshold, angular resolution and sensitivity of MAGIC under different background light levels, based on Crab Nebula observations and tuned Monte Carlo simulations. This study includes observations taken under non-standard hardware configurations, such as reducing the camera photomultiplier tubes gain by a factor ~1.7 (Reduced HV settings) with respect to standard settings (Nominal HV) or using UV-pass filters to strongly reduce the amount of moonlight reaching the cameras of the telescopes. The Crab Nebula spectrum is correctly reconstructed in all the studied illumination levels, that reach up to 30 times brighter than under dark conditions. The main effect of moonlight is an increase in the analysis energy threshold and in the systematic uncertainties on the flux normalization. The sensitivity degradation is constrained to be below 10%, within 15-30% and between 60 and 80% for Nominal HV, Reduced HV and UV-pass filter observations, respectively. No worsening of the angular resolution was found. Thanks to observations during moonlight, the maximal duty cycle of MAGIC can be increased from ~18%, under dark nights only, to up to ~40% in total with only moderate performance degradation. △ Less

Submitted 2 August, 2017; v1 submitted 4 April, 2017; originally announced April 2017.

Comments: 14 pages, 13 figures, accepted in Astroparticle Physics. Corresponding authors: Daniel Guberman (dguberman@ifae.es) and Pierre Colin (colin@mppmu.mpg.de)

Journal ref: Astropart.Phys. 94 (2017) 29-41

First Multi-wavelength Campaign on the Gamma-ray-loud Active Galaxy IC 310

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, R. Carosi, A. Carosi, A. Chatterjee, P. Colin, E. Colombo , et al. (138 additional authors not shown)

Abstract: The extragalactic VHE gamma-ray sky is rich in blazars. These are jetted active galactic nuclei viewed at a small angle to the line-of-sight. Only a handful of objects viewed at a larger angle are known so far to emit above 100 GeV. Multi-wavelength studies of such objects up to the highest energies provide new insights into the particle and radiation processes of active galactic nuclei. We report… ▽ More The extragalactic VHE gamma-ray sky is rich in blazars. These are jetted active galactic nuclei viewed at a small angle to the line-of-sight. Only a handful of objects viewed at a larger angle are known so far to emit above 100 GeV. Multi-wavelength studies of such objects up to the highest energies provide new insights into the particle and radiation processes of active galactic nuclei. We report the results from the first multi-wavelength campaign observing the TeV detected nucleus of the active galaxy IC 310, whose jet is observed at a moderate viewing angle of 10 deg - 20 deg. The multi-instrument campaign was conducted between 2012 Nov. and 2013 Jan., and involved observations with MAGIC, Fermi, INTEGRAL, Swift, OVRO, MOJAVE and EVN. These observations were complemented with archival data from the AllWISE and 2MASS catalogs. A one-zone synchrotron self-Compton model was applied to describe the broad-band spectral energy distribution. IC 310 showed an extraordinary TeV flare at the beginning of the campaign, followed by a low, but still detectable TeV flux. Compared to previous measurements, the spectral shape was found to be steeper during the low emission state. Simultaneous observations in the soft X-ray band showed an enhanced energy flux state and a harder-when-brighter spectral shape behaviour. No strong correlated flux variability was found in other frequency regimes. The broad-band spectral energy distribution obtained from these observations supports the hypothesis of a double-hump structure. The harder-when-brighter trend in the X-ray and VHE emission is consistent with the behaviour expected from a synchrotron self-Compton scenario. The contemporaneous broad-band spectral energy distribution is well described with a one-zone synchrotron self-Compton model using parameters that are comparable to those found for other gamma-ray-emitting misaligned blazars. △ Less

Submitted 22 March, 2017; originally announced March 2017.

Comments: 15 pages, 10 figures. Accepted for publication in A&A

Journal ref: A&A 603, A25 (2017)

Very-high-energy gamma-ray observations of the Type Ia Supernova SN 2014J with the MAGIC telescopes

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, C. Arcaro, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, R. Carosi, A. Carosi, A. Chatterjee , et al. (127 additional authors not shown)

Abstract: 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… ▽ More 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^{-Γ}$, with a spectral index of $Γ=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. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Comments: 5 pages, 2 figures. Accepted for publication in A&A

Journal ref: A&A 602, A98 (2017)

MAGIC detection of very high energy gamma-ray emission from the low-luminosity blazar 1ES 1741+196

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, C. Arcaro, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero , et al. (137 additional authors not shown)

Abstract: 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… ▽ More 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 $σ$ 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 $Γ= 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. △ Less

Submitted 22 February, 2017; originally announced February 2017.

Comments: Accepted for publication in MNRAS, corresponding authors: Nijil Mankuzhiyil, Massimo Persic, Saverio Lombardi, Josefa Becerra

Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (268 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 30 December, 2016; originally announced December 2016.

Comments: 33 pages, 15 figures, accepted for publication in A&A. Corresponding authors: Marlene Doert (marlene.doert@tu-dortmund.de) and David Paneque (dpaneque@mppmu.mpg.de)

Journal ref: A&A 603, A31 (2017)

Observations of Sagittarius A* during the pericenter passage of the G2 object with MAGIC

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, C. Arcaro, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero , et al. (131 additional authors not shown)

Abstract: Context. We present the results of a multi-year monitoring campaign of the Galactic Center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a… ▽ More Context. We present the results of a multi-year monitoring campaign of the Galactic Center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a unique chance to study the effect of in-falling matter on the broad-band emission of a SMBH. Aims. We search for potential flaring emission of very-high-energy (VHE; $\geq$100 GeV) gamma rays from the direction of the SMBH at the GC due to the passage of the G2 object. Using these data we also study the morphology of this complex region. Methods. We observed the GC region with the MAGIC Imaging Atmospheric Cherenkov Telescopes during the period 2012-2015, collecting 67 hours of good-quality data. In addition to a search for variability in the flux and spectral shape of the GC gamma-ray source, we use a point-source subtraction technique to remove the known gamma-ray emitters located around the GC in order to reveal the TeV morphology of the extended emission inside that region. Results. No effect of the G2 object on the VHE gamma-ray emission from the GC was detected during the 4 year observation campaign. We confirm previous measurements of the VHE spectrum of Sagittarius A*, and do not detect any significant variability of the emission from the source. Furthermore, the known VHE gamma-ray emitter at the location of the supernova remnant G0.9+0.1 was detected, as well as the recently discovered VHE source close to the GG radio Arc. △ Less

Submitted 21 November, 2016; originally announced November 2016.

Comments: 11 pages, 8 figures, accepted to A&A. Corresponding authors: Christian Fruck (fruck@mpp.mpg.de), Ievgen Vovk (ievgen.vovk@mpp.mpg.de) and John Ennis Ward (jward@ifae.es)

Journal ref: A&A 601, A33 (2017)

A search for spectral hysteresis and energy-dependent time lags from X-ray and TeV gamma-ray observations of Mrk 421

Authors: A. U. Abeysekara, S. Archambault, A. Archer, W. Benbow, R. Bird, M. Buchovecky, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, X. Chen, L. Ciupik, M. P. Connolly, W. Cui, J. D. Eisch, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, A. Flinders, L. Fortson, A. Furniss, S. Griffin, M. Hütten, N. Håkansson , et al. (198 additional authors not shown)

Abstract: 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" (ToO) observations of Mrk 421,… ▽ More 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" (ToO) 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 X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) 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 $\gtrsim 4\times 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. △ Less

Submitted 14 November, 2016; originally announced November 2016.

Comments: 45 pages, 15 figures

Hierarchical Cluster Assembly in Globally Collapsing Clouds

Authors: Enrique Vazquez-Semadeni, Alejandro Gonzalez-Samaniego, Pedro Colin

Abstract: We discuss the mechanism of cluster formation in a numerical simulation of a molecular cloud (MC) undergoing global hierarchical collapse (GHC). The global nature of the collapse implies that the SFR increases over time. The hierarchical nature of the collapse consists of small-scale collapses within larger-scale ones. The large-scale collapses culminate a few Myr later than the small-scale ones a… ▽ More We discuss the mechanism of cluster formation in a numerical simulation of a molecular cloud (MC) undergoing global hierarchical collapse (GHC). The global nature of the collapse implies that the SFR increases over time. The hierarchical nature of the collapse consists of small-scale collapses within larger-scale ones. The large-scale collapses culminate a few Myr later than the small-scale ones and consist of filamentary flows that accrete onto massive central clumps. The small-scale collapses form clumps that are embedded in the filaments and falling onto the large-scale collapse centers. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps. Thus, the filaments feed both gaseous and stellar material to the massive central clump. This leads to the presence of a few older stars in a region where new protostars are forming, and also to a self-similar structure, in which each unit is composed of smaller-scale sub-units that approach each other and may merge. Because the older stars formed in the filaments share the infall motion of the gas onto the central clump, they tend to have larger velocities and to be distributed over larger areas than the younger stars formed in the central clump. Finally, interpreting the IMF at face-value as a probability distribution implies that massive stars only form once the {\it local} SFR is large enough to sample the IMF up to high masses. In combination with the increase of the SFR, this implies that massive stars tend to appear late in the evolution of the MC, and only in the central massive clumps. We discuss the correspondence of these features with observed properties of young stellar clusters, finding very good qualitative agreement, thus providing support to the scenario of global, hierarchical collapse of MCs, while explaining the origin of the observed cluster structure. △ Less

Submitted 31 October, 2016; originally announced November 2016.

Comments: 17 pages, 12 figures. Submitted to MNRAS. Comments welcome

Multiwavelength observations of a VHE gamma-ray flare from PKS 1510-089 in 2015

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, R. Carosi, A. Carosi, A. Chatterjee, P. Colin , et al. (151 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 28 February, 2017; v1 submitted 28 October, 2016; originally announced October 2016.

Comments: 12 pages, 10 figures, accepted for a publication in A&A

Journal ref: A&A 603, A29 (2017)

Contributions of the Cherenkov Telescope Array (CTA) to the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016)

Authors: The CTA Consortium, :, A. Abchiche, U. Abeysekara, Ó. Abril, F. Acero, B. S. Acharya, C. Adams, G. Agnetta, F. Aharonian, A. Akhperjanian, A. Albert, M. Alcubierre, J. Alfaro, R. Alfaro, A. J. Allafort, R. Aloisio, J. -P. Amans, E. Amato, L. Ambrogi, G. Ambrosi, M. Ambrosio, J. Anderson, M. Anduze, E. O. Angüner , et al. (1387 additional authors not shown)

Abstract: List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany. List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: Index of CTA conference proceedings for the Gamma 2016, Heidelberg, Germany

Very High-Energy Gamma-Ray Follow-Up Program Using Neutrino Triggers from IceCube

Authors: IceCube Collaboration, M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, K. Andeen, T. Anderson, I. Ansseau, G. Anton, M. Archinger, C. Arguelles, J. Auffenberg, S. Axani, X. Bai, S. W. Barwick, V. Baum, R. Bay, J. J. Beatty, J. Becker-Tjus, K. -H. Becker, S. BenZvi , et al. (519 additional authors not shown)

Abstract: We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-f… ▽ More We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015. △ Less

Submitted 12 November, 2016; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: accepted for publication in JINST

Journal ref: 2016 JINST 11 P11009

Detection of very high energy gamma-ray emission from the gravitationally-lensed blazar QSO B0218+357 with the MAGIC telescopes

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, C. Arcaro, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee , et al. (129 additional authors not shown)

Abstract: 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 energ… ▽ More 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. △ Less

Submitted 5 September, 2016; originally announced September 2016.

Comments: 11 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 595, A98 (2016)

Cosmological simulations of Milky Way-sized galaxies

Authors: Pedro Colin, Vladimir Avila-Reese, Santi Roca-Fabrega, Octavio Valenzuela

Abstract: We introduce a new set of eight Milky Way-sized cosmological simulations performed using the AMR code ART + Hydrodynamics in a LCDM cosmology. The set of zoom-in simulations covers present-day virial masses in the 0.83-1.56 x 10^12 msun range and is carried out with our simple but effective deterministic star formation (SF) and ``explosive' stellar feedback prescriptions. The work is focused on sh… ▽ More We introduce a new set of eight Milky Way-sized cosmological simulations performed using the AMR code ART + Hydrodynamics in a LCDM cosmology. The set of zoom-in simulations covers present-day virial masses in the 0.83-1.56 x 10^12 msun range and is carried out with our simple but effective deterministic star formation (SF) and ``explosive' stellar feedback prescriptions. The work is focused on showing the goodness of the simulated set of ``field' Milky Way-sized galaxies. Our results are as follows. (a) The circular velocity curves of our simulated galaxies are nearly flat. (b) Runs ending with a significant disk component, for their stellar masses, have V_max, radius, SF rate, gas fraction, and specific angular momentum values consistent with observations of late-type galaxies. (C) The two most spheroid-dominated galaxies formed in halos with late active merger histories, but other run that ends also as spheroid-dominated, never had major mergers. (d) Our simulations are consistent with the empirical stellar-to-halo mass correlation, and those that end as disk-dominated, evolve mostly along the low-mass branch of this correlation. (e) Moreover, since the last 6.5-10 Gyr, the baryonic/stellar and halo mass growth histories are proportional. (f) Within Rvir ~ 25-50% of the baryons are missed. (g) The z ~ 0 gas velocity dispersion profiles, sigma_z(r), are nearly flat and can be mostly explained by the kinetic energy injected by stars. (h) The average values of sigma_z increase at higher redshifts, following roughly the shape of the SF history. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 17 pages, 14 figures, ApJ accepted

Long-term multi-wavelength variability and correlation study of Markarian 421 from 2007 to 2009

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (154 additional authors not shown)

Abstract: 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 v… ▽ More 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ähovi 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ähovi 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. △ Less

Submitted 29 May, 2016; originally announced May 2016.

Comments: Corresponding authors: Ann-Kristin Overkemping (ann-kristin.overkemping@tu-dortmund.de), Marina Manganaro (manganaro@iac.es), Diego Tescaro (diego.tescaro@gmail.com), To be published in Astronomy&Astrophysics (A&A), 12 pages, 9 figures

Journal ref: A&A 593, A91 (2016)

Multi-Wavelength Observations of the Blazar 1ES 1011+496 in Spring 2008

Authors: The MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra Gonzalez, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (136 additional authors not shown)

Abstract: The BL Lac object 1ES 1011+496 was discovered at Very High Energy gamma-rays by MAGIC in spring 2007. Before that the source was little studied in different wavelengths. Therefore a multi-wavelength (MWL) campaign was organized in spring 2008. Along MAGIC, the MWL campaign included the Metsahovi radio observatory, Bell and KVA optical telescopes and the Swift and AGILE satellites. MAGIC observatio… ▽ More The BL Lac object 1ES 1011+496 was discovered at Very High Energy gamma-rays by MAGIC in spring 2007. Before that the source was little studied in different wavelengths. Therefore a multi-wavelength (MWL) campaign was organized in spring 2008. Along MAGIC, the MWL campaign included the Metsahovi radio observatory, Bell and KVA optical telescopes and the Swift and AGILE satellites. MAGIC observations span from March to May, 2008 for a total of 27.9 hours, of which 19.4 hours remained after quality cuts. The light curve showed no significant variability. The differential VHE spectrum could be described with a power-law function. Both results were similar to those obtained during the discovery. Swift XRT observations revealed an X-ray flare, characterized by a harder when brighter trend, as is typical for high synchrotron peak BL Lac objects (HBL). Strong optical variability was found during the campaign, but no conclusion on the connection between the optical and VHE gamma-ray bands could be drawn. The contemporaneous SED shows a synchrotron dominated source, unlike concluded in previous work based on nonsimultaneous data, and is well described by a standard one zone synchrotron self Compton model. We also performed a study on the source classification. While the optical and X-ray data taken during our campaign show typical characteristics of an HBL, we suggest, based on archival data, that 1ES 1011+496 is actually a borderline case between intermediate and high synchrotron peak frequency BL Lac objects. △ Less

Submitted 23 March, 2016; originally announced March 2016.

Comments: 13 pages, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 459, Issue 3, p.2286-2298 (2016)

Super-orbital variability of LS I +61°303 at TeV energies

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (128 additional authors not shown)

Abstract: 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 pre… ▽ More 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 ($φ$=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 ($φ$=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 $φ$=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$α$ fluxes. △ Less

Submitted 21 April, 2016; v1 submitted 22 March, 2016; originally announced March 2016.

Comments: 7 pages, 4 figures, to be published in A&A

Journal ref: A&A 591, A76 (2016)

Insights into the emission of the blazar 1ES 1011+496 through unprecedented broadband observations during 2011 and 2012

Authors: J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, C. Arcaro, A. Babic, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, P. Colin, E. Colombo, J. L. Contreras, J. Cortina , et al. (139 additional authors not shown)

Abstract: 1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $γ$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to… ▽ More 1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $γ$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to explain the observed broadband spectrum. We analyzed VHE to radio data from 2011 and 2012 collected by MAGIC, $Fermi$-LAT, $Swift$, KVA, OVRO, and Metsähovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. The VHE spectrum was fit with a simple power law with a photon index of $3.69\pm0.22$ and a flux above 150 GeV of $(1.46\pm0.16)\times10^{-11}$ ph cm$^{-2}$ s$^{-1}$. 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of $1.8\pm0.4\,c$, which is the highest speed statistically significantly measured so far in a high-frequency-peaked BL Lac. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field to the optical emission. △ Less

Submitted 6 June, 2016; v1 submitted 22 March, 2016; originally announced March 2016.

Comments: 14 pages, 7 figures

Journal ref: A&A 591, A10 (2016)

Investigating the peculiar emission from the new VHE gamma-ray source H1722+119

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (137 additional authors not shown)

Abstract: 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… ▽ More 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) $γ$-ray band with a statistical significance of 5.9 $σ$. 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) $γ$-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. △ Less

Submitted 30 September, 2016; v1 submitted 21 March, 2016; originally announced March 2016.

Comments: 12 pages, 5 figures, 2 tables

Journal ref: MNRAS, 2016, 459, 3271-3281

Search for VHE gamma-ray emission from Geminga pulsar and nebula with the MAGIC telescopes

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra Gonzalez, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (130 additional authors not shown)

Abstract: The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wi… ▽ More The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 x 10^{-13} TeV cm^{-2} s^{-1} for the Geminga pulsar and 3.5 x 10^{-12} TeV cm^{-2} s^{-1} for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult. △ Less

Submitted 5 March, 2016; v1 submitted 2 March, 2016; originally announced March 2016.

Comments: 8 pages, 6 figures

Cluster Assembly in Hierarchically Collapsing Clouds

Authors: Enrique Vazquez-Semadeni, Alejandro Gonzalez-Samaniego, Manuel Zamora-Aviles, Pedro Colin

Abstract: We discuss the mechanism of cluster formation in hierarchically collapsing molecular clouds. Recent evidence, both observational and numerical, suggests that molecular clouds (MCs) may be undergoing global, hierarchical gravitational collapse. The "hierarchical" regime consists of small-scale collapses within larger-scale ones. The latter implies that the star formation rate increases systematical… ▽ More We discuss the mechanism of cluster formation in hierarchically collapsing molecular clouds. Recent evidence, both observational and numerical, suggests that molecular clouds (MCs) may be undergoing global, hierarchical gravitational collapse. The "hierarchical" regime consists of small-scale collapses within larger-scale ones. The latter implies that the star formation rate increases systematically during the early stages of evolution, and occurs via filamentary flows onto "hubs" of higher density, mass, and velocity dispersion, and culminates a few Myr after than the small-scale collapses have started to form stars. In turn, the small-scale collapses occur in clumps embedded in the filaments, and are themselves falling into the larger potential well of the still-ongoing large-scale collapse. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps towards the larger potential trough, so that the filaments feed both gaseous and stellar material to the hubs. This leads to the presence of older stars in a region where new protostars are still forming, to a scale-free or fractal structure of the clusters, in which each unit is composed of smaller-scale ones, and to the eventual merging of the subunits, explaining the observed structural features of open clusters. △ Less

Submitted 23 February, 2016; originally announced February 2016.

Comments: 6 pages, 7 figures, to appear in the proceedings of IAU Symposium 316: Formation, Evolution, and Survival of Massive Star Clusters