Showing 1–37 of 37 results for author: Somero, A

SN2017egm: A Helium-rich Superluminous Supernova with Multiple Bumps in the Light Curves

Authors: Jiazheng Zhu, Ning Jiang, Subo Dong, Alexei V. Filippenko, Richard J. Rudy, A. Pastorello, Christopher Ashall, Subhash Bose, R. S. Post, D. Bersier, Stefano Benetti, Thomas G. Brink, Ping Chen, Liming Dou, N. Elias-Rosa, Peter Lundqvist, Seppo Mattila, Ray W. Russell, Michael L. Sitko, Auni Somero, M. D. Stritzinger, Tinggui Wang, Peter J. Brown, E. Cappellaro, Morgan Fraser , et al. (6 additional authors not shown)

Abstract: When discovered, SN~2017egm was the closest (redshift $z=0.03$) hydrogen-poor superluminous supernova (SLSN-I) and a rare case that exploded in a massive and metal-rich galaxy. Thus, it has since been extensively observed and studied. We report spectroscopic data showing strong emission at around He~I $λ$10,830 and four He~I absorption lines in the optical. Consequently, we classify SN~2017egm as… ▽ More When discovered, SN~2017egm was the closest (redshift $z=0.03$) hydrogen-poor superluminous supernova (SLSN-I) and a rare case that exploded in a massive and metal-rich galaxy. Thus, it has since been extensively observed and studied. We report spectroscopic data showing strong emission at around He~I $λ$10,830 and four He~I absorption lines in the optical. Consequently, we classify SN~2017egm as a member of an emerging population of helium-rich SLSNe-I (i.e., SLSNe-Ib). We also present our late-time photometric observations. By combining them with archival data, we analyze high-cadence ultra-violet, optical, and near-infrared light curves spanning from early pre-peak ($\sim -20\,d$) to late phases ($\sim +300\,d$). We obtain its most complete bolometric light curve, in which multiple bumps are identified. None of the previously proposed models can satisfactorily explain all main light-curve features, while multiple interactions between the ejecta and circumstellar material (CSM) may explain the undulating features. The prominent infrared excess with a blackbody luminosity of $10^7$--$10^8\,L_{sun}$ detected in SN~2017egm could originate from the emission of either an echo of a pre-existing dust shell, or newly-formed dust, offering an additional piece of evidence supporting the ejecta-CSM interaction model. Moreover, our analysis of deep $Chandra$ observations yields the tightest-ever constraint on the X-ray emission of an SLSN-I, amounting to an X-ray-to-optical luminosity ratio $\lesssim 10^{-3}$ at late phases ($\sim100-200\,d$), which could help explore its close environment and central engine. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 25 pages, 14 Figures, 4 Tables; accepted for publication in ApJ (Mar. 2023)

V838 Mon: A slow waking up of Sleeping Beauty?

Authors: T. Liimets, I. Kolka, M. Kraus, T. Eenmäe, T. Tuvikene, T. Augusteijn, L. Antunes Amaral, A. A. Djupvik, J. H. Telting, B. Deshev, E. Kankare, J. Kankare, J. E. Lindberg, T. M. Amby, T. Pursimo, A. Somero, A. Thygesen, P. A. Strøm

Abstract: Context. V838 Monocerotis is a peculiar binary that underwent an immense stellar explosion in 2002, leaving behind an expanding cool supergiant and a hot B3V companion. Five years after the outburst, the B3V companion disappeared from view, and so far did not recover. Aims. We investigate the changes in the light curve and spectral features Methods. A monitoring campaign has been performed during… ▽ More Context. V838 Monocerotis is a peculiar binary that underwent an immense stellar explosion in 2002, leaving behind an expanding cool supergiant and a hot B3V companion. Five years after the outburst, the B3V companion disappeared from view, and so far did not recover. Aims. We investigate the changes in the light curve and spectral features Methods. A monitoring campaign has been performed during the past 13 years with the Nordic Optical Telescope to obtain optical photometric and spectroscopic data. The data sets are used to analyse the temporal evolution of the spectral features and the spectral energy distribution, and to characterize the object. Results. Our photometric data show a steady brightening in all bands during the past 13 years, which is particularly prominent in the blue. This rise is also reflected in the spectra, showing a gradual relative increase in the continuum flux at shorter wavelengths. In addition, a slow brightening of the Ha emission line starting in 2015 was detected. These changes might imply that the B3V companion is slowly reappearing. During the same time interval, our analysis reveals a considerable change in the observed colours of the object along with a steady decrease in the strength and width of molecular absorption bands in our low-resolution spectra. These changes suggest a rising temperature of the cool supergiant along with a weakening of its wind, most likely combined with a slow recovery of the secondary due to the evaporation of the dust and accretion of the material from the shell in which the hot companion is embedded. From our medium-resolution spectra, we find that the heliocentric radial velocity of the atomic absorption line of TiI 6556.06 A has been stable for more than a decade. We propose that TiI lines are tracing the velocity of the red supergiant in V838 Mon, and not representing the infalling matter as previously stated. △ Less

Submitted 30 November, 2022; v1 submitted 12 November, 2022; originally announced November 2022.

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

Journal ref: A&A 670, A13 (2023)

The nuclear transient AT 2017gge: a tidal disruption event in a dusty and gas-rich environment and the awakening of a dormant SMBH

Authors: F. Onori, G. Cannizzaro, P. G. Jonker, M. Kim, M. Nicholl, S. Mattila, T. M. Reynolds, M. Fraser, T. Wevers, E. Brocato, J. P. Anderson, R. Carini, P. Charalampopoulos, P. Clark, M. Gromadzki, C. P. Gutiérrez, N. Ihanec, C. Inserra, A. Lawrence, G. Leloudas, P. Lundqvist, T. E. Müller-Bravo, S. Piranomonte, M. Pursiainen, K. A. Rybicki , et al. (6 additional authors not shown)

Abstract: We present the results from a dense multi-wavelength (optical/UV, near-infrared (IR), and X-ray) follow-up campaign of the nuclear transient AT2017gge, covering a total of 1698 days from the transient's discovery. The bolometric lightcurve, the black body temperature and radius, the broad H and He I $λ$5876 emission lines and their evolution with time, are all consistent with a tidal disruption ev… ▽ More We present the results from a dense multi-wavelength (optical/UV, near-infrared (IR), and X-ray) follow-up campaign of the nuclear transient AT2017gge, covering a total of 1698 days from the transient's discovery. The bolometric lightcurve, the black body temperature and radius, the broad H and He I $λ$5876 emission lines and their evolution with time, are all consistent with a tidal disruption event (TDE) nature. A soft X-ray flare is detected with a delay of $\sim$200 days with respect to the optical/UV peak and it is rapidly followed by the emergence of a broad He II $λ$4686 and by a number of long-lasting high ionization coronal emission lines. This indicate a clear connection between a TDE flare and the appearance of extreme coronal line emission (ECLEs). An IR echo, resulting from dust re-radiation of the optical/UV TDE light is observed after the X-ray flare and the associated near-IR spectra show a transient broad feature in correspondence of the He I $λ$10830 and, for the first time in a TDE, a transient high-ionization coronal NIR line (the [Fe XIII] $λ$10798) is also detected. The data are well explained by a scenario in which a TDE occurs in a gas and dust rich environment and its optical/UV, soft X-ray, and IR emission have different origins and locations. The optical emission may be produced by stellar debris stream collisions prior to the accretion disk formation, which is instead responsible for the soft X-ray flare, emitted after the end of the circularization process. △ Less

Submitted 9 September, 2022; v1 submitted 31 May, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS

Intermediate-luminosity red transients: Spectrophotometric properties and connection to electron-capture supernova explosions

Authors: Y. -Z. Cai, A. Pastorello, M. Fraser, M. T. Botticella, N. Elias-Rosa, L. -Z. Wang, R. Kotak, S. Benetti, E. Cappellaro, M. Turatto, A. Reguitti, S. Mattila, S. J. Smartt, C. Ashall, S. Benitez, T. -W. Chen, A. Harutyunyan, E. Kankare, P. Lundqvist, P. A. Mazzali, A. Morales-Garoffolo, P. Ochner, G. Pignata, S. J. Prentice, T. M. Reynolds , et al. (34 additional authors not shown)

Abstract: We present the spectroscopic and photometric study of five intermediate-luminosity red transients (ILRTs), namely AT 2010dn, AT 2012jc, AT 2013la, AT 2013lb, and AT 2018aes. They share common observational properties and belong to a family of objects similar to the prototypical ILRT SN~2008S. These events have a rise time that is less than 15 days and absolute peak magnitudes of between $-11.5$ an… ▽ More We present the spectroscopic and photometric study of five intermediate-luminosity red transients (ILRTs), namely AT 2010dn, AT 2012jc, AT 2013la, AT 2013lb, and AT 2018aes. They share common observational properties and belong to a family of objects similar to the prototypical ILRT SN~2008S. These events have a rise time that is less than 15 days and absolute peak magnitudes of between $-11.5$ and $-14.5$ mag. Their pseudo-bolometric light curves peak in the range $0.5$ - $9.0 \times10^{40}~\mathrm{erg~s}^{-1}$ and their total radiated energies are on the order of (0.3 - 3) $\times$~10$^{47}$~erg. After maximum brightness, the light curves show a monotonic decline or a plateau, resembling those of faint supernovae IIL or IIP, respectively. At late phases, the light curves flatten, roughly following the slope of the $^{56}$Co decay. If the late-time power source is indeed radioactive decay, these transients produce $^{56}$Ni masses on the order of $10^{-4}$ to $10^{-3}$~\msun. The spectral energy distribution of our ILRT sample, extending from the optical to the mid-infrared (MIR) domain, reveals a clear IR excess soon after explosion and non-negligible MIR emission at very late phases. The spectra show prominent H lines in emission with a typical velocity of a few hundred km~s$^{-1}$, along with Ca~II features. In particular, the [Ca~II] $λ$7291,7324 doublet is visible at all times, which is a characteristic feature for this family of transients. The identified progenitor of SN~2008S, which is luminous in archival Spitzer MIR images, suggests an intermediate-mass precursor star embedded in a dusty cocoon. We propose the explosion of a super-asymptotic giant branch star forming an electron-capture supernova as a plausible explanation for these events. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: 31 pages, 17 figures; accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 654, A157 (2021)

ASASSN-18am/SN 2018gk : An overluminous Type IIb supernova from a massive progenitor

Authors: Subhash Bose, Subo Dong, C. S. Kochanek, M. D. Stritzinger, Chris Ashall, Stefano Benetti, E. Falco, Alexei V. Filippenko, Andrea Pastorello, Jose L. Prieto, Auni Somero, Tuguldur Sukhbold, Junbo Zhang, Katie Auchettl, Thomas G. Brink, J. S. Brown, Ping Chen, A. Fiore, Dirk Grupe, T. W. -S. Holoien, Peter Lundqvist, Seppo Mattila, Robert Mutel, David Pooley, R. S. Post , et al. (7 additional authors not shown)

Abstract: ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of $M_V \approx -20$ mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am d… ▽ More ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of $M_V \approx -20$ mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am declines extremely rapidly for a Type II SN, with a photospheric-phase decline rate of $\sim6.0~\rm mag~(100 d)^{-1}$. Owing to the weakening of HI and the appearance of HeI in its later phases, ASASSN-18am is spectroscopically a Type IIb SN with a partially stripped envelope. However, its photometric and spectroscopic evolution show significant differences from typical SNe IIb. Using a radiative diffusion model, we find that the light curve requires a high synthesised $\rm ^{56}Ni$ mass $M_{\rm Ni} \sim0.4~M_\odot$ and ejecta with high kinetic energy $E_{\rm kin} = (7-10) \times10^{51} $ erg. Introducing a magnetar central engine still requires $M_{\rm Ni} \sim0.3~M_\odot$ and $E_{\rm kin}= 3\times10^{51} $ erg. The high $\rm ^{56}Ni$ mass is consistent with strong iron-group nebular lines in its spectra, which are also similar to several SNe Ic-BL with high $\rm ^{56}Ni$ yields. The earliest spectrum shows "flash ionisation" features, from which we estimate a mass-loss rate of $ \dot{M}\approx 2\times10^{-4}~\rm M_\odot~yr^{-1} $. This wind density is too low to power the luminous light curve by ejecta-CSM interaction. We measure expansion velocities as high as $ 17,000 $ km/s for $H_α$, which is remarkably high compared to other SNe II. We estimate an oxygen core mass of $1.8-3.4$ $M_\odot$ using the [OI] luminosity measured from a nebular-phase spectrum, implying a progenitor with a zero-age main sequence mass of $19-26$ $M_\odot$. △ Less

Submitted 24 February, 2021; v1 submitted 30 June, 2020; originally announced July 2020.

Observation of inverse Compton emission from a long $γ$-ray burst

Authors: 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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella, Y. Chai , et al. (279 additional authors not shown)

Abstract: Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the ex… ▽ More Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Journal ref: Nature 575 (2019) 459-463

Teraelectronvolt emission from the $γ$-ray burst GRB 190114C

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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, A. Carosi, R. Carosi , et al. (155 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to… ▽ More Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to GeV gamma rays, due to synchrotron radiation by energetic electrons in accompanying shock waves. Although emission of gamma rays at even higher, TeV energies by other radiation mechanisms had been theoretically predicted, it had never been detected previously. Here we report the clear detection of GRB 190114C in the TeV band, achieved after many years of dedicated searches for TeV emission from GRBs. Gamma rays in the energy range 0.2--1 TeV are observed from about 1 minute after the burst (at more than 50 standard deviations in the first 20 minutes). This unambiguously reveals a new emission component in the afterglow of a GRB, whose power is comparable to that of the synchrotron component. The observed similarity in the radiated power and temporal behaviour of the TeV and X-ray bands points to processes such as inverse Compton radiation as the mechanism of the TeV emission, while processes such as synchrotron emission by ultrahigh-energy protons are disfavoured due to their low radiative efficiency. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Journal ref: Nature 575 (2019) 455-458

MAGIC observations of the diffuse $γ$-ray emission in the vicinity of the Galactic Centre

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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (151 additional authors not shown)

Abstract: Aims: $γ$ rays can be used as a tracer in the search of sources of Galactic cosmic rays (CRs). We present deep observations of the Galactic Centre (GC) region with the MAGIC telescopes, which we use for inferring the underlying CR distribution. Methods: We observed the GC region for ${\approx}100$ hours with the MAGIC telescopes from 2012 to 2017, at high zenith angles (58-70~deg). This implies a… ▽ More Aims: $γ$ rays can be used as a tracer in the search of sources of Galactic cosmic rays (CRs). We present deep observations of the Galactic Centre (GC) region with the MAGIC telescopes, which we use for inferring the underlying CR distribution. Methods: We observed the GC region for ${\approx}100$ hours with the MAGIC telescopes from 2012 to 2017, at high zenith angles (58-70~deg). This implies a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using new software, we derive instrument response and background models, enabling us to study the diffuse emission in the region. We use pre-existing data of the gas distribution in the GC region to derive the underlying distribution of CRs. Results: We obtain a significant detection for all four model components used to fit our data (Sgr~A*, ``Arc'', G0.9+0.1, and an extended component for the Galactic Ridge). We find that the diffuse component is best described as a power-law with index 2 and an exponential cut-off at around 20~TeV with the significance of the cut-off being only 2~$σ$. The derived cosmic-ray profile hints to a peak at the GC position, with a measured profile index of $1.2 \pm 0.3$, supporting the hypothesis of a CR accelerator at the GC. We argue that the measurements of this profile are presently limited by our knowledge of the gas distribution in the GC vicinity. △ Less

Submitted 31 May, 2020; originally announced June 2020.

Comments: 10 pages, 10 figures. Accepted for publication by A&A. Corresponding authors: Christian Fruck (fruck@mpp.mpg.de), Ievgen Vovk (vovk@icrr.u-tokyo.ac.jp), Yuki Iwamura (iwamura@icrr.u-tokyo.ac.jp) and Marcel Strzys (strzys@icrr.u-tokyo.ac.jp)

Journal ref: A&A 642, A190 (2020)

A search for dark matter in Triangulum II 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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (150 additional authors not shown)

Abstract: We present the first results from very-high-energy observations of the dwarf spheroidal satellite candidate Triangulum II with the MAGIC telescopes from 62.4 hours of good-quality data taken between August 2016 and August 2017. We find no gamma-ray excess in the direction of Triangulum II, and upper limits on both the differential and integral gamma-ray flux are presented. Currently, the kinematic… ▽ More We present the first results from very-high-energy observations of the dwarf spheroidal satellite candidate Triangulum II with the MAGIC telescopes from 62.4 hours of good-quality data taken between August 2016 and August 2017. We find no gamma-ray excess in the direction of Triangulum II, and upper limits on both the differential and integral gamma-ray flux are presented. Currently, the kinematics of Triangulum II are affected by large uncertainties leading to a bias in the determination of the properties of its dark matter halo. Using a scaling relation between the annihilation J-factor and heliocentric distance of well-known dwarf spheroidal galaxies, we estimate an annihilation J-factor for Triangulum II for WIMP dark matter of $\log[J_{\text{ann}}({0.5^{\circ}})/$ GeV$^{2}$ cm$^{-5}] = 19.35 \pm 0.37$. We also derive a dark matter density profile for the object relying on results from resolved simulations of Milky Way sized dark matter halos. We obtain 95% confidence-level limits on the thermally averaged annihilation cross section for WIMP annihilation into various Standard Model channels. The most stringent limits are obtained in the $τ^{+}τ^{-}$ final state, where a cross section for annihilation down to $\langle σ_{\text{ann}} v \rangle = 3.05 \times 10^{-24}$ cm$^{3}$ s$^{-1}$ is excluded. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: 14 pages, 8 figures, accepted for publication in Physics of the Dark Universe

Journal ref: Physics of the Dark Universe Volume 28, May 2020, 100529

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

Authors: A. U. Abeysekara, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. K. Daniel, J. Dumm, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, N. Galante, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (234 additional authors not shown)

Abstract: We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a co… ▽ More We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $γ$-ray (VERITAS, MAGIC), high-energy (HE) $γ$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Metsähovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($δ\gtrsim 33$) and the size of the emission region ($ δ^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: 41 pages including 3 appendices, 13 figures; version accepted to Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 890, Number 2 (2020)

Broadband characterisation of the very intense TeV flares of the blazar 1ES 1959+650 in 2016

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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Z. Bosnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (142 additional authors not shown)

Abstract: 1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a broad emission in the high-energy regime, that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavel… ▽ More 1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a broad emission in the high-energy regime, that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the Fermi Large Area Telescope (LAT) and Swift instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE, E >100 GeV) gamma-ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching ~ 3 times the Crab Nebula flux was measured by MAGIC on the 13th, 14th of June and 1st July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the gamma-ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13th June and 1st July 2016, the source showed rapid variations of the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to high values of the Doppler factors (>=30-60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power L_{jet}~10^{46} erg/s and under high values of the magnetic field strength (~100 G) and maximum proton energy (~few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016. △ Less

Submitted 31 January, 2020; originally announced February 2020.

Comments: 18 pages, 9 figures, accepted to A&A

Journal ref: A&A 638, A14 (2020)

MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV

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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (147 additional authors not shown)

Abstract: Aims: We aim to measure the Crab Nebula gamma-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies. Methods: We use the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also develop an auxiliary procedure of monitoring atmospheric… ▽ More Aims: We aim to measure the Crab Nebula gamma-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies. Methods: We use the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also develop an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording of optical images of the stellar field next to the source position, which provides a better than 10% accuracy for the transmission measurements. Results: We demonstrate that MAGIC very large zenith angle observations yield a collection area larger than a square kilometer. In only ~56 hr of observations, we detect the gamma-ray emission from the Crab Nebula up to 100 TeV, thus providing the highest energy measurement of this source to date with Imaging Atmospheric Cherenkov Telescopes. Comparing accumulated and archival MAGIC and Fermi/LAT data with some of the existing emission models, we find that none of them provides an accurate description of the 1 GeV to 100 TeV gamma-ray signal. △ Less

Submitted 26 January, 2020; originally announced January 2020.

Comments: 10 pages, 10 figures, accepted to A&A

Journal ref: A&A 635, A158 (2020)

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)

Monitoring of the radio galaxy M87 during a low emission state from 2012 to 2015 with 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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto , et al. (166 additional authors not shown)

Abstract: M87 is one of the closest (z=0.00436) extragalactic sources emitting at very-high-energies (VHE, E > 100 GeV). The aim of this work is to locate the region of the VHE gamma-ray emission and to describe the observed broadband spectral energy distribution (SED) during the low VHE gamma-ray state. The data from M87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed a… ▽ More M87 is one of the closest (z=0.00436) extragalactic sources emitting at very-high-energies (VHE, E > 100 GeV). The aim of this work is to locate the region of the VHE gamma-ray emission and to describe the observed broadband spectral energy distribution (SED) during the low VHE gamma-ray state. The data from M87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed and combined with multi-wavelength data from Fermi-LAT, Chandra, HST, EVN, VLBA and the Liverpool Telescope. The averaged VHE gamma-ray spectrum can be fitted from 100GeV to 10TeV with a simple power law with a photon index of (-2.41 $\pm$ 0.07), while the integral flux above 300GeV is $(1.44 \pm 0.13) \times 10^{-12} cm^{-2} s^{-1}$. During the campaign between 2012 and 2015, M87 is generally found in a low emission state at all observed wavelengths. The VHE gamma-ray flux from the present 2012-2015 M87 campaign is consistent with a constant flux with some hint of variability ($\sim3σ$) on a daily timescale in 2013. The low-state gamma-ray emission likely originates from the same region as the flare-state emission. Given the broadband SED, both a leptonic synchrotron self Compton and a hybrid photo-hadronic model reproduce the available data well, even if the latter is preferred. We note, however, that the energy stored in the magnetic field in the leptonic scenario is very low suggesting a matter dominated emission region. △ Less

Submitted 1 August, 2020; v1 submitted 6 January, 2020; originally announced January 2020.

Comments: 12 pages, 4 figures, Published on MNRAS. Corresponding authors (alphabetical ordered): C. Arcaro, P. Bangale, M. Manganaro, D. Mazin

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 492, Issue 4 (2020), p.5354-5365

New hard-TeV extreme blazars detected with the MAGIC telescopes

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

Abstract: Extreme high-frequency peaked BL Lac objects (EHBLs) are blazars which exhibit extremely energetic synchrotron emission. They also feature non-thermal gamma-ray emission whose peak lies in the very high-energy (VHE, E > 100 GeV) range, and in some sources exceeds 1TeV: this is the case of hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, ten targets were observed… ▽ More Extreme high-frequency peaked BL Lac objects (EHBLs) are blazars which exhibit extremely energetic synchrotron emission. They also feature non-thermal gamma-ray emission whose peak lies in the very high-energy (VHE, E > 100 GeV) range, and in some sources exceeds 1TeV: this is the case of hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, ten targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 262 h of good quality data. The data were complemented by coordinated Swift observations. The X-ray data analysis confirms that all the sources but two are EHBLs. The sources show only a modest variability and a harder-when-brighter behavior, typical for this class of objects. At VHE gamma rays, three new sources were detected and a hint of signal was found for another new source. In each case the intrinsic spectrum is compatible with the hypothesis of a hard-TeV nature of these EHBLs. The broadband spectral energy distributions (SEDs) of all sources are built and modeled in the framework of a single-zone purely leptonic model. The VHE gamma-ray detected sources were also interpreted with a spine-layer and a proton synchrotron models. The three models provide a good description of the SEDs. However, the resulting parameters differ substantially in the three scenarios, in particular the magnetization parameter. This work presents a first mini-catalog of VHE gamma-ray and multi-wavelength observations of EHBLs. △ Less

Submitted 20 August, 2020; v1 submitted 15 November, 2019; originally announced November 2019.

Comments: 30 pages, 9 figures. Corresponding author: E. Prandini. Main authors: C. Arcaro, V. Fallah Ramazani

Journal ref: The Astrophysical Journal Supplement Series, Volume 247, Issue 1, id.16, 24 pp. (2020)

SN 2016gsd: An unusually luminous and linear type II supernova with high velocities

Authors: T. M. Reynolds, M. Fraser, S. Mattila, M. Ergon, P. Lundqvist, L. Dessart, Subo Dong, N. Elias-Rosa, L. Galbany, C. P. Gutiérrez, T. Kangas, E. Kankare, R. Kotak, H. Kuncarayakti, A. Pastorello, O. Rodriguez, S. J. Smartt, M. Stritzinger, L. Tomasella, Ping Chen, J. Harmanen, G. Hozzeinzadeh, D. Andrew Howell, C. Inserra, M. Nicholl , et al. (5 additional authors not shown)

Abstract: We present observations of the unusually luminous Type II supernova (SN) 2016gsd. With a peak absolute magnitude of V = $-$19.95 $\pm$ 0.08, this object is one of the brightest Type II SNe, and lies in the gap of magnitudes between the majority of Type II SNe and the superluminous SNe. Its light curve shows little evidence of the expected drop from the optically thick phase to the radioactively po… ▽ More We present observations of the unusually luminous Type II supernova (SN) 2016gsd. With a peak absolute magnitude of V = $-$19.95 $\pm$ 0.08, this object is one of the brightest Type II SNe, and lies in the gap of magnitudes between the majority of Type II SNe and the superluminous SNe. Its light curve shows little evidence of the expected drop from the optically thick phase to the radioactively powered tail. The velocities derived from the absorption in H$α$ are also unusually high with the blue edge tracing the fastest moving gas initially at 20000 km s$^{-1}$, and then declining approximately linearly to 15000 km s$^{-1}$ over $\sim$100 d. The dwarf host galaxy of the SN indicates a low-metallicity progenitor which may also contribute to the weakness of the metal lines in its spectra. We examine SN 2016gsd with reference to similarly luminous, linear Type II SNe such as SNe 1979C and 1998S, and discuss the interpretation of its observational characteristics. We compare the observations with a model produced by the JEKYLL code and find that a massive star with a depleted and inflated hydrogen envelope struggles to reproduce the high luminosity and extreme linearity of SN 2016gsd. Instead, we suggest that the influence of interaction between the SN ejecta and circumstellar material can explain the majority of the observed properties of the SN. The high velocities and strong H$α$ absorption present throughout the evolution of the SN may imply a circumstellar medium configured in an asymmetric geometry. △ Less

Submitted 3 June, 2020; v1 submitted 30 September, 2019; originally announced September 2019.

Comments: 25 pages, 17 figures, Published by MNRAS. Re-uploaded to reflect published version. Changes made to Figures 1,8,11 to give more information, new figure 15 added. Changes to text to improve clarity and add additional information and analysis requested by the referee

Testing emission models on the extreme blazar 2WHSP J073326.7+515354 detected at very high energies 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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (151 additional authors not shown)

Abstract: Extreme high-energy peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multi-wavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up t… ▽ More Extreme high-energy peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multi-wavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073326.7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within a MWL context, including an analysis of the Fermi-LAT data over ten years of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EHBL. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broad-band SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-band SED. △ Less

Submitted 19 November, 2019; v1 submitted 25 September, 2019; originally announced September 2019.

Comments: MNRAS published

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 490, Issue 2, December 2019, Pages 2284-2299

The long-lived Type IIn SN 2015da: Infrared echoes and strong interaction within an extended massive shell

Authors: L. Tartaglia, A. Pastorello, J. Sollerman, C. Fransson, S. Mattila, M. Fraser, F. Taddia, L. Tomasella, M. Turatto, A. Morales-Garoffolo, N. Elias-Rosa, P. Lundqvist, J. Harmanen, T. Reynolds, E. Cappellaro, C. Barbarino, A. Nyholm, E. Kool, E. Ofek, X. Gao, Z. Jin, H. Tan, D. J. Sand, F. Ciabattari, X. Wang , et al. (19 additional authors not shown)

Abstract: In this paper we report the results of the first $\sim$four years of spectroscopic and photometric monitoring of the Type IIn supernova SN 2015da (also known as PSN J13522411+3941286, or iPTF16tu). The supernova exploded in the nearby spiral galaxy NGC 5337 in a relatively highly extinguished environment. The transient showed prominent narrow Balmer lines in emission at all times and a slow rise t… ▽ More In this paper we report the results of the first $\sim$four years of spectroscopic and photometric monitoring of the Type IIn supernova SN 2015da (also known as PSN J13522411+3941286, or iPTF16tu). The supernova exploded in the nearby spiral galaxy NGC 5337 in a relatively highly extinguished environment. The transient showed prominent narrow Balmer lines in emission at all times and a slow rise to maximum in all bands. In addition, early observations performed by amateur astronomers give a very well-constrained explosion epoch. The observables are consistent with continuous interaction between the supernova ejecta and a dense and extended H-rich circumstellar medium. The presence of such an extended and dense medium is difficult to reconcile with standard stellar evolution models, since the metallicity at the position of SN 2015da seems to be slightly subsolar. Interaction is likely the mechanism powering the light curve, as confirmed by the analysis of the pseudo bolometric light curve, which gives a total radiated energy $\gtrsim10^{51}\,\rm{erg}$. Modeling the light curve in the context of a supernova shock breakout through a dense circumstellar medium allowed us to infer the mass of the prexisting gas to be $\simeq8\,\rm{M_{\odot}}$, with an extreme mass-loss rate for the progenitor star $\simeq0.6\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, suggesting that most of the circumstellar gas was produced during multiple eruptive events. Near- and mid-infrared observations reveal a flux excess in these domains, similar to those observed in SN 2010jl and other interacting transients, likely due to preexisting radiatively heated dust surrounding the supernova. By modeling the infrared excess, we infer a mass $\gtrsim0.4\times10^{-3}\,\rm{M_{\odot}}$ for the dust. △ Less

Submitted 21 January, 2020; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 23 pages (19+appendix), 15 figures, 11 tables (Tables A.1-A.4 at CDS only). Accepted for publication on Astronomy and Astrophysics

Journal ref: A&A 635, A39 (2020)

Constraints on gamma-ray and neutrino emission from NGC 1068 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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (150 additional authors not shown)

Abstract: Starburst galaxies and star-forming active galactic nuclei (AGN) are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the do… ▽ More Starburst galaxies and star-forming active galactic nuclei (AGN) are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hours of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very high energy band. We did not detect significant gamma-ray emission, and set upper limits at 95\% confidence level to the gamma-ray flux above 200 GeV f<5.1x10^{-13} cm^{-2} s ^{-1} . This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr^{-1}. △ Less

Submitted 12 August, 2019; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 11 pages, 4 figures, accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal 883 (2019) 135

Measurement of the Extragalactic Background Light using MAGIC and Fermi-LAT gamma-ray observations of blazars up to z = 1

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, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, G. Busetto, R. Carosi, G. Ceribella, Y. Chai , et al. (135 additional authors not shown)

Abstract: We present a measurement of the extragalactic background light (EBL) based on a joint likelihood analysis of 32 gamma-ray spectra for 12 blazars in the redshift range z = 0.03 to 0.944, obtained by the MAGIC telescopes and Fermi-LAT. The EBL is the part of the diffuse extragalactic radiation spanning the ultraviolet, visible and infrared bands. Major contributors to the EBL are the light emitted b… ▽ More We present a measurement of the extragalactic background light (EBL) based on a joint likelihood analysis of 32 gamma-ray spectra for 12 blazars in the redshift range z = 0.03 to 0.944, obtained by the MAGIC telescopes and Fermi-LAT. The EBL is the part of the diffuse extragalactic radiation spanning the ultraviolet, visible and infrared bands. Major contributors to the EBL are the light emitted by stars through the history of the universe, and the fraction of it which was absorbed by dust in galaxies and re-emitted at longer wavelengths. The EBL can be studied indirectly through its effect on very-high energy photons that are emitted by cosmic sources and absorbed via photon-photon interactions during their propagation across cosmological distances. We obtain estimates of the EBL density in good agreement with state-of-the-art models of the EBL production and evolution. The 1-sigma upper bounds, including systematic uncertainties, are between 13% and 23% above the nominal EBL density in the models. No anomaly in the expected transparency of the universe to gamma rays is observed in any range of optical depth.We also perform a wavelength-resolved EBL determination, which results in a hint of an excess of EBL in the 0.18 - 0.62 $μ$m range relative to the studied models, yet compatible with them within systematics. △ Less

Submitted 29 March, 2019; originally announced April 2019.

Comments: Accepted by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 486, Issue 3, July 2019, Pages 4233-4251

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

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

Strongly Bipolar Inner Ejecta of the Normal Type IIP Supernova ASASSN-16at

Authors: Subhash Bose, Subo Dong, N. Elias-Rosa, B. J. Shappee, David Bersier, Stefano Benetti, M. D. Stritzinger, D. Grupe, C. S. Kochanek, J. L. Prieto, Ping Chen, H. Kuncarayakti, Seppo Mattila, Antonia Morales-Garoffolo, Nidia Morrell, F. Onori, Thomas M Reynolds, A. Siviero, Auni Somero, K. Z. Stanek, Giacomo Terreran, Todd A. Thompson, L. Tomasella, C. Ashall, Christa Gall , et al. (2 additional authors not shown)

Abstract: We report distinctly double-peakedH-alpha and H-beta emission lines in the late-time, nebular-phase spectra (>~200 d) of the otherwise normal at early phases (<~ 100 d) Type IIP supernova ASASSN-16at (SN 2016X). Such distinctly double-peaked nebular Balmer lines have never been observed for a Type II SN. The nebular-phase Balmer emission is driven by the radioactive Co56 decay, so the observed lin… ▽ More We report distinctly double-peakedH-alpha and H-beta emission lines in the late-time, nebular-phase spectra (>~200 d) of the otherwise normal at early phases (<~ 100 d) Type IIP supernova ASASSN-16at (SN 2016X). Such distinctly double-peaked nebular Balmer lines have never been observed for a Type II SN. The nebular-phase Balmer emission is driven by the radioactive Co56 decay, so the observed line-profile bifurcation suggests a strong bipolarity in the Ni56 distribution or in the line-forming region of the inner ejecta. The strongly bifurcated blue- and red-shifted peaks are separated by ~3x10^3 km/s and are roughly symmetrically positioned with respect to the host-galaxy rest frame, implying that the inner ejecta are composed of two almost detached blobs. The red peak progressively weakens relative to the blue peak, and disappears in the 740 d spectrum. One possible reason for the line-ratio evolution is increasing differential extinction from continuous formation of dust within the envelope, which is also supported by the near-infrared flux excess that develops after ~100 d. △ Less

Submitted 3 March, 2019; v1 submitted 29 October, 2018; originally announced October 2018.

Comments: Published in ApJL

Journal ref: Bose et al 2019 ApJL 873 L3

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

AT 2017be - a new member of the class of Intermediate-Luminosity Red Transients

Authors: Y-Z. Cai, A. Pastorello, M. Fraser, M. T. Botticella, C. Gall, I. Arcavi, S. Benetti, E. Cappellaro, N. Elias-Rosa, J. Harmanen, G. Hosseinzadeh, D. A. Howell, J. Isern, T. Kangas, E. Kankare, H. Kuncarayakti, P. Lundqvist, S. Mattila, C. McCully, T. M. Reynolds, A. Somero, M. D. Stritzinger, G. Terreran

Abstract: We report the results of our spectrophotometric monitoring campaign for AT~2017be in NGC~2537. Its lightcurve reveals a fast rise to an optical maximum, followed by a plateau lasting about 30 days, and finally a fast decline. Its absolute peak magnitude ($M_{r}$ $\simeq$ $-$12 $\rm{mag}$) is fainter than that of core-collapse supernovae, and is consistent with those of supernova impostors and othe… ▽ More We report the results of our spectrophotometric monitoring campaign for AT~2017be in NGC~2537. Its lightcurve reveals a fast rise to an optical maximum, followed by a plateau lasting about 30 days, and finally a fast decline. Its absolute peak magnitude ($M_{r}$ $\simeq$ $-$12 $\rm{mag}$) is fainter than that of core-collapse supernovae, and is consistent with those of supernova impostors and other Intermediate-Luminosity Optical Transients. The quasi-bolometric lightcurve peaks at $\sim$ 2 $\times$ 10$^{40}$ erg s$^{-1}$, and the late-time photometry allows us to constrain an ejected $^{56}$Ni mass of $\sim$ 8 $\times$ 10$^{-4}$\msun. The spectra of AT~2017be show minor evolution over the observational period, a relatively blue continuum showing at early phases, which becomes redder with time. A prominent H$α$ emission line always dominates over other Balmer lines. Weak Fe {\sc ii} features, Ca~{\sc ii} H$\&$K and the Ca {\sc ii} NIR triplet are also visible, while P-Cygni absorption troughs are found in a high resolution spectrum. In addition, the [Ca~{\sc ii}] $λ$7291,7324 doublet is visible in all spectra. This feature is typical of Intermediate-Luminosity Red Transients (ILRTs), similar to SN~2008S. The relatively shallow archival Spitzer data are not particularly constraining. On the other hand, a non-detection in deeper near-infrared HST images disfavours a massive Luminous Blue Variable eruption as the origin for AT~2017be. As has been suggested for other ILRTs, we propose that AT~2017be is a candidate for a weak electron-capture supernova explosion of a super-asymptotic giant branch star, still embedded in a thick dusty envelope. △ Less

Submitted 31 July, 2018; originally announced July 2018.

Comments: 21 pages, 15 figures, accepted by MNRAS

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)

SN 2017dio: a type-Ic supernova exploding in a hydrogen-rich circumstellar medium

Authors: Hanindyo Kuncarayakti, Keiichi Maeda, Christopher J. Ashall, Simon J. Prentice, Seppo Mattila, Erkki Kankare, Claes Fransson, Peter Lundqvist, Andrea Pastorello, Giorgos Leloudas, Joseph P. Anderson, Stefano Benetti, Melina C. Bersten, Enrico Cappellaro, Regis Cartier, Larry Denneau, Massimo Della Valle, Nancy Elias-Rosa, Gaston Folatelli, Morgan Fraser, Lluis Galbany, Christa Gall, Avishay Gal-Yam, Claudia P. Gutierrez, Aleksandra Hamanowicz , et al. (20 additional authors not shown)

Abstract: SN 2017dio shows both spectral characteristics of a type-Ic supernova (SN) and signs of a hydrogen-rich circumstellar medium (CSM). Prominent, narrow emission lines of H and He are superposed on the continuum. Subsequent evolution revealed that the SN ejecta are interacting with the CSM. The initial SN Ic identification was confirmed by removing the CSM interaction component from the spectrum and… ▽ More SN 2017dio shows both spectral characteristics of a type-Ic supernova (SN) and signs of a hydrogen-rich circumstellar medium (CSM). Prominent, narrow emission lines of H and He are superposed on the continuum. Subsequent evolution revealed that the SN ejecta are interacting with the CSM. The initial SN Ic identification was confirmed by removing the CSM interaction component from the spectrum and comparing with known SNe Ic, and reversely, adding a CSM interaction component to the spectra of known SNe Ic and comparing them to SN 2017dio. Excellent agreement was obtained with both procedures, reinforcing the SN Ic classification. The light curve constrains the pre-interaction SN Ic peak absolute magnitude to be around $M_g = -17.6$ mag. No evidence of significant extinction is found, ruling out a brighter luminosity required by a SN Ia classification. These pieces of evidence support the view that SN 2017dio is a SN Ic, and therefore the first firm case of a SN Ic with signatures of hydrogen-rich CSM in the early spectrum. The CSM is unlikely to have been shaped by steady-state stellar winds. The mass loss of the progenitor star must have been intense, $\dot{M} \sim 0.02$ $(ε_{Hα}/0.01)^{-1}$ $(v_\textrm{wind}/500$ km s$^{-1}$) $(v_\textrm{shock}/10 000$ km s$^{-1})^{-3}$ $M_\odot$~yr$^{-1}$, peaking at a few decades before the SN. Such a high mass loss rate might have been experienced by the progenitor through eruptions or binary stripping. △ Less

Submitted 6 February, 2018; v1 submitted 30 November, 2017; originally announced December 2017.

Comments: ApJL, 11 pages, 5 figures (accepted 22 Jan 2018)

Gaia17biu/SN 2017egm in NGC 3191: The closest hydrogen-poor superluminous supernova to date is in a "normal", massive, metal-rich spiral galaxy

Authors: Subhash Bose, Subo Dong, A. Pastorello, Alexei V. Filippenko, C. S. Kochanek, Jon Mauerhan, C. Romero-Canizales, Thomas Brink, Ping Chen, J. L. Prieto, R. Post, Christopher Ashall, Dirk Grupe, L. Tomasella, Stefano Benetti, B. J. Shappee, K. Z. Stanek, Zheng Cai, E. Falco, Peter Lundqvist, Seppo Mattila, Robert Mutel, Paolo Ochner, David Pooley, M. D. Stritzinger , et al. (33 additional authors not shown)

Abstract: Hydrogen-poor superluminous supernovae (SLSNe-I) have been predominantly found in low-metallicity, star-forming dwarf galaxies. Here we identify Gaia17biu/SN 2017egm as an SLSN-I occurring in a "normal" spiral galaxy (NGC 3191) in terms of stellar mass (several times 10^10 M_sun) and metallicity (roughly Solar). At redshift z=0.031, Gaia17biu is also the lowest redshift SLSN-I to date, and the abs… ▽ More Hydrogen-poor superluminous supernovae (SLSNe-I) have been predominantly found in low-metallicity, star-forming dwarf galaxies. Here we identify Gaia17biu/SN 2017egm as an SLSN-I occurring in a "normal" spiral galaxy (NGC 3191) in terms of stellar mass (several times 10^10 M_sun) and metallicity (roughly Solar). At redshift z=0.031, Gaia17biu is also the lowest redshift SLSN-I to date, and the absence of a larger population of SLSNe-I in dwarf galaxies of similar redshift suggests that metallicity is likely less important to the production of SLSNe-I than previously believed. With the smallest distance and highest apparent brightness for an SLSN-I, we are able to study Gaia17biu in unprecedented detail. Its pre-peak near-ultraviolet to optical color is similar to that of Gaia16apd and among the bluest observed for an SLSN-I while its peak luminosity (M_g = -21 mag) is substantially lower than Gaia16apd. Thanks to the high signal-to-noise ratios of our spectra, we identify several new spectroscopic features that may help to probe the properties of these enigmatic explosions. We detect polarization at the ~0.5% level that is not strongly dependent on wavelength, suggesting a modest, global departure from spherical symmetry. In addition, we put the tightest upper limit yet on the radio luminosity of an SLSN-I with <5.4x10^26 erg/s/Hz (at 10 GHz), which is almost a factor of 40 better than previous upper limits and one of the few measured at an early stage in the evolution of an SLSN-I. This limit largely rules out an association of this SLSNe-I with known populations of gamma-ray burst (GRB) like central engines. △ Less

Submitted 25 December, 2017; v1 submitted 2 August, 2017; originally announced August 2017.

Comments: Accepted for publication in ApJ. Ancillary ASCII tables added: TRL.txt -- blackbody temperature, radius and luminosity; uvw2uvm2uvw1uvu.txt -- UV photometry; BgVri.txt -- optical photometry; zJHK.txt -- NIR photometry

The 67P/Churyumov-Gerasimenko observation campaign in support of the Rosetta mission

Authors: C. Snodgrass, M. F. A'Hearn, F. Aceituno, V. Afanasiev, S. Bagnulo, J. Bauer, G. Bergond, S. Besse, N. Biver, D. Bodewits, H. Boehnhardt, B. P. Bonev, G. Borisov, B. Carry, V. Casanova, A. Cochran, B. C. Conn, B. Davidsson, J. K. Davies, J. de León, E. de Mooij, M. de Val-Borro, M. Delacruz, M. A. DiSanti, J. E. Drew , et al. (90 additional authors not shown)

Abstract: We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond t… ▽ More We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: Author prepared version; final published version available at journal. 22 pages

Journal ref: Phil. Trans. R. Soc. A 375, 20160249 (2017)

Gaia16apd -- a link between fast-and slowly-declining type I superluminous supernovae

Authors: T. Kangas, N. Blagorodnova, S. Mattila, P. Lundqvist, M. Fraser, U. Burgaz, E. Cappellaro, J. M. Carrasco Martínez, N. Elias-Rosa, L. K. Hardy, J. Harmanen, E. Y. Hsiao, J. Isern, E. Kankare, Z. Kołaczkowski, M. B. Nielsen, T. M. Reynolds, L. Rhodes, A. Somero, M. D. Stritzinger, Ł. Wyrzykowski

Abstract: We present ultraviolet, optical and infrared photometry and optical spectroscopy of the type Ic superluminous supernova (SLSN) Gaia16apd (= SN 2016eay), covering its evolution from 26 d before the $g$-band peak to 234.1 d after the peak. Gaia16apd was followed as a part of the NOT Unbiased Transient Survey (NUTS). It is one of the closest SLSNe known ($z = 0.102\pm0.001$), with detailed optical an… ▽ More We present ultraviolet, optical and infrared photometry and optical spectroscopy of the type Ic superluminous supernova (SLSN) Gaia16apd (= SN 2016eay), covering its evolution from 26 d before the $g$-band peak to 234.1 d after the peak. Gaia16apd was followed as a part of the NOT Unbiased Transient Survey (NUTS). It is one of the closest SLSNe known ($z = 0.102\pm0.001$), with detailed optical and ultraviolet (UV) observations covering the peak. Gaia16apd is a spectroscopically typical type Ic SLSN, exhibiting the characteristic blue early spectra with O II absorption, and reaches a peak $M_{g} = -21.8 \pm 0.1$ mag. However, photometrically it exhibits an evolution intermediate between the fast- and slowly-declining type Ic SLSNe, with an early evolution closer to the fast-declining events. Together with LSQ12dlf, another SLSN with similar properties, it demonstrates a possible continuum between fast- and slowly-declining events. It is unusually UV-bright even for a SLSN, reaching a non-$K$-corrected $M_{uvm2} \simeq -23.3$ mag, the only other type Ic SLSN with similar UV brightness being SN 2010gx. Assuming that Gaia16apd was powered by magnetar spin-down, we derive a period of $P = 1.9\pm0.2$ ms and a magnetic field of $B = 1.9\pm0.2 \times 10^{14}$ G for the magnetar. The estimated ejecta mass is between 8 and 16 $\mathrm{M}_{\odot}$ and the kinetic energy between 1.3 and $2.5 \times 10^{52}$ erg, depending on opacity and assuming that the entire ejecta is swept up into a thin shell. Despite the early photometric differences, the spectra at late times are similar to slowly-declining type Ic SLSNe, implying that the two subclasses originate from similar progenitors. △ Less

Submitted 5 June, 2017; v1 submitted 30 November, 2016; originally announced November 2016.

Comments: 14 pages; 8 figures. Updated on June 5 2017; now published in MNRAS after revision

Journal ref: 2017MNRAS.469.1246K

High-resolution optical spectroscopy of RS Ophiuchi during 2008 -- 2009

Authors: A. Somero, P. Hakala, G. A. Wynn

Abstract: RS Ophiuchi is a symbiotic variable and a recurrent nova. We have monitored it with the Nordic Optical Telescope and obtained 30 high resolution (R=46 000) optical spectra over one orbital cycle during quiescence. To our knowledge this is the best-sampled high resolution spectroscopic dataset of RS Oph over one orbital period. We do not detect any direct signatures of an accretion disc such as dou… ▽ More RS Ophiuchi is a symbiotic variable and a recurrent nova. We have monitored it with the Nordic Optical Telescope and obtained 30 high resolution (R=46 000) optical spectra over one orbital cycle during quiescence. To our knowledge this is the best-sampled high resolution spectroscopic dataset of RS Oph over one orbital period. We do not detect any direct signatures of an accretion disc such as double peaked emission lines, but many line profiles are complex consisting of superimposed emission and absorption components. We measure the spin of the red giant and conclude that it is tidally locked to the binary orbit. We observe Na I absorption features, probably arising from the circumbinary medium, that has been shaped by previous recurrent nova outbursts. We do not detect any intrinsic polarisation in the optical wavelengths. △ Less

Submitted 4 October, 2016; originally announced October 2016.

Comments: 9 pages, 7 figures + Appendix with trailed spectrograms. Accepted for publication in MNRAS

Phase-resolved optical and X-ray spectroscopy of low-mass X-ray binary X1822-371

Authors: A. Somero, P. Hakala, P. Muhli, P. Charles, O. Vilhu

Abstract: (Abridged) X1822-371 is the prototypical accretion disc corona X-ray source, a low-mass X-ray binary viewed at very high inclination, thereby allowing the disc structure and extended disc coronal regions to be visible. We study the structure of the accretion disc in X1822-371 by modelling the phase-resolved spectra both in optical and X-ray regime. We analyse high time resolution optical ESO/VLT s… ▽ More (Abridged) X1822-371 is the prototypical accretion disc corona X-ray source, a low-mass X-ray binary viewed at very high inclination, thereby allowing the disc structure and extended disc coronal regions to be visible. We study the structure of the accretion disc in X1822-371 by modelling the phase-resolved spectra both in optical and X-ray regime. We analyse high time resolution optical ESO/VLT spectra of X1822-371 to study the variability in the emission line profiles. In addition, we use data from XMM-Newton space observatory to study phase-resolved as well as high resolution X-ray spectra. We apply the Doppler tomography technique to reconstruct a map of the optical emission distribution in the system. We fit multi-component models to the X-ray spectra. We find that our results from both the optical and X-ray analysis can be explained with a model where the accretion disc has a thick rim in the region where the accretion stream impacts the disc. The behaviour of the H_beta line complex implies that some of the accreting matter creates an outburst around the accretion stream impact location and that the resulting outflow of matter moves both away from the accretion disc and towards the centre of the disc. Such behaviour can be explained by an almost isotropic outflow of matter from the accretion stream impact region. The optical emission lines of HeII 4686 and 5411 show double peaked profiles, typical for an accretion disc at high inclination. However, their velocities are slower than expected for an accretion disc in a system like X1822-371. This, combined with the fact that the HeII emission lines do not get eclipsed during the partial eclipse in the continuum, suggests that the line emission does not originate in the orbital plane and is more likely to come from above the accretion disc, for example the accretion disc wind. △ Less

Submitted 17 January, 2012; originally announced January 2012.

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

First Kepler results on compact pulsators VI. Targets in the final half of the survey phase

Authors: R. H. Østensen, R. Silvotti, S. Charpinet, R. Oreiro, S. Bloemen, A. S. Baran, M. D. Reed, S. D. Kawaler, J. H. Telting, E. M. Green, S. J. O'Toole, C. Aerts, B. T. Gänsicke, T. R. Marsh, E. Breedt, U. Heber, D. Koester, A. C. Quint, D. W. Kurtz, C. Rodríguez-López, M. Vučković, T. A. Ottosen, S. Frimann, A. Somero, P. A. Wilson , et al. (7 additional authors not shown)

Abstract: We present results from the final six months of a survey to search for pulsations in white dwarfs and hot subdwarf stars with the Kepler spacecraft. Spectroscopic observations are used to separate the objects into accurate classes, and we explore the physical parameters of the subdwarf B (sdB) stars and white dwarfs in the sample. From the Kepler photometry and our spectroscopic data, we find that… ▽ More We present results from the final six months of a survey to search for pulsations in white dwarfs and hot subdwarf stars with the Kepler spacecraft. Spectroscopic observations are used to separate the objects into accurate classes, and we explore the physical parameters of the subdwarf B (sdB) stars and white dwarfs in the sample. From the Kepler photometry and our spectroscopic data, we find that the sample contains 5 new pulsators of the V1093 Her type, one AM CVn type cataclysmic variable, and a number of other binary systems. This completes the survey for compact pulsators with Kepler. No V361 Hya type of short-period pulsating sdB stars were found in this half, leaving us with a total of one single multiperiodic V361 Hya and 13 V1093 Her pulsators for the full survey. Except for the sdB pulsators, no other clearly pulsating hot subdwarfs or white dwarfs were found, although a few low-amplitude candidates still remain. The most interesting targets discovered in this survey will be observed throughout the remainder of the Kepler Mission, providing the most long-term photometric datasets ever made on such compact, evolved stars. Asteroseismic investigations of these datasets will be invaluable in revealing the interior structure of these stars, and will boost our understanding of their evolutionary history. △ Less

Submitted 21 January, 2011; originally announced January 2011.

Comments: 11 pages, accepted for publication in MNRAS

A planetary nebula around nova V458 Vul undergoing flash ionization

Authors: R. Wesson, M. J. Barlow, R. L. M. Corradi, J. E. Drew, P. J. Groot, C. Knigge, D. Steeghs, B. T. Gaensicke, R. Napiwotzki, P. Rodriguez-Gil, A. A. Zijlstra, M. F. Bode, J. J. Drake, D. J. Frew, E. A. Gonzalez-Solares, R. Greimel, M. J. Irwin, L. Morales-Rueda, G. Nelemans, Q. A. Parker, S. E. Sale, J. L. Sokoloski, A. Somero, H. Uthas, N. A. Walton , et al. (3 additional authors not shown)

Abstract: Nova V458 Vul erupted on 2007 August 8th and reached a visual magnitude of 8.1 a few days later. H$α$ images obtained six weeks before the outburst as part of the IPHAS galactic plane survey reveal an 18th magnitude progenitor surrounded by an extended nebula. Subsequent images and spectroscopy of the nebula reveal an inner nebular knot increasing rapidly in brightness due to flash ionization by… ▽ More Nova V458 Vul erupted on 2007 August 8th and reached a visual magnitude of 8.1 a few days later. H$α$ images obtained six weeks before the outburst as part of the IPHAS galactic plane survey reveal an 18th magnitude progenitor surrounded by an extended nebula. Subsequent images and spectroscopy of the nebula reveal an inner nebular knot increasing rapidly in brightness due to flash ionization by the nova event. We derive a distance of 13 kpc based on light travel time considerations, which is supported by two other distance estimation methods. The nebula has an ionized mass of 0.2 M$_{\odot}$ and a low expansion velocity: this rules it out as ejecta from a previous nova eruption, and is consistent with it being a ~14,000 year old planetary nebula, probably the product of a prior common envelope (CE) phase of evolution of the binary system. The large derived distance means that the mass of the erupting WD component of the binary is high. We identify two possible evolutionary scenarios, in at least one of which the system is massive enough to produce a Type Ia supernova on merging. △ Less

Submitted 2 October, 2008; originally announced October 2008.

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