Showing 1–50 of 147 results for author: Shore, S N

Shear, writhe and filaments: turbulence in the high latitude molecular cloud MBM 40

Authors: Marco Monaci, Loris Magnani, Steven N. Shore, Henrik Olofsson, Mackenzie R. Joy

Abstract: Context. It is almost banal to say that the interstellar medium (ISM) is structurally and thermodynamically complex. But the variety of the governing processes, including stellar feedback, renders the investigation challenging. High latitude molecular clouds (HLMCs) with no evidence of internal star formation, such as MBM 40, are excellent sites for studying the chemistry and dynamic evolution of… ▽ More Context. It is almost banal to say that the interstellar medium (ISM) is structurally and thermodynamically complex. But the variety of the governing processes, including stellar feedback, renders the investigation challenging. High latitude molecular clouds (HLMCs) with no evidence of internal star formation, such as MBM 40, are excellent sites for studying the chemistry and dynamic evolution of the cold neutral ISM. Aims. We used this high latitude cloud as an exemplar for the dynamical and chemical processes in the diffuse interstellar medium. Methods. We analyzed new and archival $^{12}$CO, $^{13}$CO, CH, HCO$^+$, CS, H$_2$CO, HCN data from Five College Radio Observatory (FCRAO), Onsala Space Observatory (OSO), Arizona Radio Observatory (ARO) and W. Gordon telescope (Arecibo) combined with the Galactic Arecibo L-band Feed Array HI (GALFA-HI) HI 21 cm data set, to study the chemistry, thermal state, and dynamics of MBM 40. A new dynamical analytical approach was adopted by considering each line profile as a line of sight Probability Distribution Function (PDF) of the turbulence weighted by gas emissivity. Results. The atomic and molecular gas are smoothly distributed in space and velocity. No steep transition is seen between circumcloud atomic and cloud molecular gas in either radial velocity or structure. We proposed a topology of the cloud from the molecular tracers, a contorted filamentary structure that is shaped by a broad embedding shear flow in the neutral atomic gas. Comparative examination of different molecular tracers shows that $^{13}$CO, H$_2$CO and CS arise from only denser molecular cores, where $^{12}$CO, CH and HCO$^+$ traces diffuse gas with broader range of dynamics. △ Less

Submitted 11 June, 2023; originally announced June 2023.

Comments: 18 pages, 22 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 676, A138 (2023)

Spectroscopic Signature of a Re-established Accretion Disk in Symbiotic-like Recurrent Nova RS Ophiuchi

Authors: Alessandra Azzollini, Steven Neil Shore, Paul Kuin

Abstract: A novel method is presented which can pin down the time the accretion disk re-established itself in the RS Oph system after it experienced a nova disruption. The method is based on the re-ionisation of the ejecta by photoionisation from the radiation released in the boundary layer from accretion. A novel method is presented which can pin down the time the accretion disk re-established itself in the RS Oph system after it experienced a nova disruption. The method is based on the re-ionisation of the ejecta by photoionisation from the radiation released in the boundary layer from accretion. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 3 pages, 1 figure

Journal ref: Res. Notes AAS 6 92 (2022)

Multiwavelength spectroscopic study of shock driven phenomena in explosive outbursts in symbiotic-like recurrent novae with emphasis on RS Ophiuchi

Authors: Alessandra Azzollini, Steven Neil Shore, Paul Kuin, Kim Page

Abstract: To detail the development of RS Ophiuchi and the other Galactic Symbiotic-like Recurrent Novae throughout their outburst and quiescence, with a particular emphasis on the propagation of the shock wave during the outburst of the binaries. The spectral analysis has been performed using archival data according to the features of the individual datasets. Swift grism spectra were reduced and extracted… ▽ More To detail the development of RS Ophiuchi and the other Galactic Symbiotic-like Recurrent Novae throughout their outburst and quiescence, with a particular emphasis on the propagation of the shock wave during the outburst of the binaries. The spectral analysis has been performed using archival data according to the features of the individual datasets. Swift grism spectra were reduced and extracted using a combination of the pre-existing UVOTPY Python routine and newly written pipelines in Matlab. Other datasets were directly available in reduced form, already corrected for instrumental or background contamination, calibrated in wavelength and flux or intensity. The work on these was done through pipelines suited for reading the data and elaborating them to extract quantities of interest for the analysis. We find striking similarities in different outbursts of the same object and for different novae. For example, RS Oph 2021 was almost identical to the 2006 outburst, despite having occurred at a different orbital phase with the observations made from a different line of sight through the red giant wind. Despite the intrinsically different properties of the binaries, striking similarities are found for different systems of the same class, for instance, that the trend of the electron density over time during outburst appears to follow a general temporal development. △ Less

Submitted 19 June, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: 20 pages, 40 figures

Journal ref: A&A 674, A139 (2023)

The mixing of dust and gas in the high latitude translucent cloud MBM 40

Authors: Marco Monaci, Loris Magnani, Steven N. Shore

Abstract: Context. High latitude molecular clouds (hereafter HLMCs) permit the study of interstellar gas dynamics and astrochemistry with good accuracy due to their proximity, generally clear lines of sight, and lack of internal star-forming activity which can heavily modify the physical context. MBM 40, one of the nearest HLMCs, has been extensively studied, making it a superb target to infer and study the… ▽ More Context. High latitude molecular clouds (hereafter HLMCs) permit the study of interstellar gas dynamics and astrochemistry with good accuracy due to their proximity, generally clear lines of sight, and lack of internal star-forming activity which can heavily modify the physical context. MBM 40, one of the nearest HLMCs, has been extensively studied, making it a superb target to infer and study the dust-to-gas mixing ratio (DGMR). Aims. The mixing of dust and gas in the interstellar medium remains a fundamental issue to keep track of astrochemistry evolution and molecular abundances. Accounting for both molecular and atomic gas is difficult because $H_2$ is not directly observable and HI spectra always show different dynamical profiles blended together which are not directly correlated with the cloud. We used two independent strategies to infer the molecular and atomic gas column densities and compute the dust-to-gas mixing ratio. Methods. We combined $HI$ 21 cm and $^{12}CO$ line observations with the IRAS 100 $μ$m image to infer the dust-to-gas mixing ratio within the cloud. The cloud 21 cm profile was extracted using a hybrid Gaussian decomposition where $^{12}CO$ was used to deduce the total molecular hydrogen column density. Infrared images were used to calculate the dust emission. Results. The dust-to-gas mixing ratio is nearly uniform within the cloud as outlined by the hairpin structure. The total hydrogen column density and 100 $μ$m emissivity are linearly correlated over a range in $N(H_{tot})$ of one order of magnitude. △ Less

Submitted 12 January, 2023; originally announced January 2023.

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

Journal ref: A&A, 668 (2022) L9

A semiempirical approach to low-energy cosmic ray propagation in the diffuse interstellar medium

Authors: Riccardo Franceschi, Steven Neil Shore

Abstract: We investigate the ionization of the diffuse interstellar medium by cosmic rays by modeling their propagation along the wandering magnetic fields using a Monte Carlo method. We study how low-energy cosmic rays propagate in turbulent, translucent molecular clouds, and how they regulate the ionization and both lose and gain energy from the medium. As a test case, we use high spatial resolution (0.03… ▽ More We investigate the ionization of the diffuse interstellar medium by cosmic rays by modeling their propagation along the wandering magnetic fields using a Monte Carlo method. We study how low-energy cosmic rays propagate in turbulent, translucent molecular clouds, and how they regulate the ionization and both lose and gain energy from the medium. As a test case, we use high spatial resolution (0.03 pc) CO maps of a well-studied high latitude translucent cloud, MBM 3, to model turbulence. The propagation problem is solved with a modified Monte Carlo procedure that includes trapping, energization, and ionization losses. In a homogeneous medium, trapping and re-energization do not produce a significant effect. In a nonuniform medium, particles can be trapped for a long time inside the cloud. This modifies the cosmic ray distribution due to stochastic acceleration at the highest energies (about 100 MeV). At lower energies, the re-energization is too weak to produce an appreciable effect. The change in the energy distribution does not significantly affect the ionization losses, so ionization changes are due to trapping effects. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 8 pages, 13 figures

Journal ref: A&A, 667, 8 (2022)

Fermi-LAT Gamma-ray Detection of the Recurrent Nova RS Ophiuchi during its 2021 Outburst

Authors: C. C. Cheung, T. J. Johnson, P. Jean, M. Kerr, K. L. Page, J. P. Osborne, A. P. Beardmore, K. V. Sokolovsky, F. Teyssier, S. Ciprini, G. Marti-Devesa, I. Mereu, S. Razzaque, K. S. Wood, S. N. Shore, S. Korotkiy, A. Levina, A. Blumenzweig

Abstract: We report the Fermi-LAT gamma-ray detection of the 2021 outburst of the symbiotic recurrent nova RS Ophiuchi. In this system, unlike classical novae from cataclysmic binaries, the ejecta from the white dwarf form shocks when interacting with the dense circumstellar wind environment of the red giant companion. We find the LAT spectra from 50 MeV to ~20-23 GeV, the highest-energy photons detected in… ▽ More We report the Fermi-LAT gamma-ray detection of the 2021 outburst of the symbiotic recurrent nova RS Ophiuchi. In this system, unlike classical novae from cataclysmic binaries, the ejecta from the white dwarf form shocks when interacting with the dense circumstellar wind environment of the red giant companion. We find the LAT spectra from 50 MeV to ~20-23 GeV, the highest-energy photons detected in some sub-intervals, are consistent with $π^{\rm 0}$-decay emission from shocks in the ejecta as proposed by Tatischeff & Hernanz (2007) for its previous 2006 outburst. The LAT light-curve displayed a fast rise to its peak >0.1 GeV flux of $\simeq$6x10^-6 ph cm^-2 s^-1 beginning on day 0.745 after its optically-constrained eruption epoch of 2021 August 8.50. The peak lasted for ~1 day, and exhibited a power-law decline up to the final LAT detection on day 45. We analyze the data on shorter timescales at early times and found evidence of an approximate doubling of emission over ~200 minutes at day 2.2, possibly indicating a localized shock-acceleration event. Comparing the data collected by the AAVSO, we measured a constant ratio of ~2.8x10^-3 between the gamma-ray and optical luminosities except for a ~5x smaller ratio within the first day of the eruption likely indicating attenuation of gamma rays by ejecta material and lower high-energy proton fluxes at the earliest stages of the shock development. The hard X-ray emission due to bremsstrahlung from shock-heated gas traced by the Swift-XRT 2-10 keV light-curve peaked at day ~6, later than at GeV and optical energies. Using X-ray derived temperatures to constrain the velocity profile, we find the hadronic model reproduces the observed >0.1 GeV light-curve. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: ApJ, accepted. 21 pages, 10 figures, 4 tables

Environmental Effects on the Dynamical Evolution of Star Clusters in Turbulent Molecular Clouds

Authors: Paolo Suin, Steven N. Shore, Václav Pavlík

Abstract: Context: Star clusters form within giant molecular clouds that are strongly altered by the feedback action of the massive stars, but the cluster still remains embedded in a dense, highly turbulent medium and interactions with ambient structures may modify its dynamical evolution from that expected if it were isolated. Aims: We aim to study coupling mechanisms between the dynamical evolution of t… ▽ More Context: Star clusters form within giant molecular clouds that are strongly altered by the feedback action of the massive stars, but the cluster still remains embedded in a dense, highly turbulent medium and interactions with ambient structures may modify its dynamical evolution from that expected if it were isolated. Aims: We aim to study coupling mechanisms between the dynamical evolution of the cluster, accelerated by the mass segregation process, with harassment effects caused by the gaseous environment. Methods: We simulated the cluster dynamical evolution combining $N$-body and hydrodynamic codes within the Astronomical Multipurpose Software Environment (AMUSE). Conclusions: Tidal harassment produces a sparser configuration more rapidly than the isolated reference simulations. The evolution of the asymptotic power-law density distribution exponent also shows substantially different behaviour in the two cases. The background is more effective on clusters in advanced stages of dynamical development. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 9 pages, 2 appendices, 4 tables, 6 figures, accepted for publication to A&A

Revisiting V1309 Sco 2008 outburst spectra. Observational evidence for theoretical modeling of stellar mergers

Authors: Elena Mason, Steven N. Shore

Abstract: CONTEXT: V1309 Sco is the only certain noncompact stellar merger, due to its indisputable preoutburst light curve matching that of a contact binary of almost equal mass stars. Therefore, anything that can be deduced from the existing observations serves as benchmark constraints for models. AIMS: We present some observational evidences to guide future hydrodynamical simulations and common envelope… ▽ More CONTEXT: V1309 Sco is the only certain noncompact stellar merger, due to its indisputable preoutburst light curve matching that of a contact binary of almost equal mass stars. Therefore, anything that can be deduced from the existing observations serves as benchmark constraints for models. AIMS: We present some observational evidences to guide future hydrodynamical simulations and common envelope studies. METHODS: Using archive spectra taken at high and mid spectral resolution during the V1309 Sco outburst and late decline, together with the inferential methods we developed to study nova ejecta through panchromatic high resolution spectroscopic follow ups, we constrain the physical state, structure, dynamics and geometry of the transient originated in the stellar merger. RESULTS: We found that the emitted spectra arise from two distinct contributions: matter expelled during the 2008 outburst and circumbinary gas produced during historic mass loss episodes. These two components likely have orthogonal geometry with the 2008 mass loss displaying a dust-laden bipolar ejecta produced by a time limited rapidly accelerating wind and the circumbinary gas having a donut-like shape. A central source powers them both, having produced a fluorescent light pulse, but we cannot precisely determine the time it started or its spectral energy distribution. We can, however, place its upper energy cutoff at about 54 eV and the bulk of its emission at $<$20 eV. We also know that the central source turned off within months from the outburst and before the ejecta turned optically thin. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 11 pages, 9 figures (png format). Accepted for publication in A&A

Journal ref: A&A 664, A12 (2022)

The 2021 outburst of the recurrent nova RS Ophiuchi observed in X-rays by the Neil Gehrels Swift Observatory: a comparative study

Authors: K. L. Page, A. P. Beardmore, J. P. Osborne, U. Munari, J. -U. Ness, P. A. Evans, M. F. Bode, M. J. Darnley, J. J. Drake, N. P. M. Kuin, T. J. O'Brien, M. Orio, S. N. Shore, S. Starrfield, C. E. Woodward

Abstract: On 2021 August 8, the recurrent nova RS Ophiuchi erupted again, after an interval of 15.5 yr. Regular monitoring by the Neil Gehrels Swift Observatory began promptly, on August 9.9 (0.37 day after the optical peak), and continued until the source passed behind the Sun at the start of November, 86 days later. Observations then restarted on day 197, once RS Oph emerged from the Sun constraint. This… ▽ More On 2021 August 8, the recurrent nova RS Ophiuchi erupted again, after an interval of 15.5 yr. Regular monitoring by the Neil Gehrels Swift Observatory began promptly, on August 9.9 (0.37 day after the optical peak), and continued until the source passed behind the Sun at the start of November, 86 days later. Observations then restarted on day 197, once RS Oph emerged from the Sun constraint. This makes RS Oph the first Galactic recurrent nova to have been monitored by Swift throughout two eruptions. Here we investigate the extensive X-ray datasets from 2006 and 2021, as well as the more limited data collected by EXOSAT in 1985. The hard X-rays arising from shock interactions between the nova ejecta and red giant wind are similar following the last two eruptions. In contrast, the early super-soft source (SSS) in 2021 was both less variable and significantly fainter than in 2006. However, 0.3-1 keV light-curves from 2021 reveal a 35 s quasi-periodic oscillation consistent in frequency with the 2006 data. The Swift X-ray spectra from 2021 are featureless, with the soft emission typically being well parametrized by a simple blackbody, while the 2006 spectra showed much stronger evidence for superimposed ionized absorption edges. Considering the data after day 60 following each eruption, during the supersoft phase the 2021 spectra are hotter, with smaller effective radii and lower wind absorption, leading to an apparently reduced bolometric luminosity. We explore possible explanations for the gross differences in observed SSS behaviour between the 2006 and 2021 outbursts. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: 21 pages, 16 colour figures, accepted for publication in MNRAS. Table 1 in full is included as an ancillary PDF (will be supplementary online material when published by MNRAS)

Close encounters with the Death Star: Interactions between collapsed bodies and the Solar System

Authors: Václav Pavlík, Steven N. Shore

Abstract: Aims: We aim to investigate the consequences of a fast massive stellar remnant - a black hole (BH) or a neutron star (NS) - encountering a planetary system. Methods: We modelled a close encounter between the actual Solar System (SS) and a $2\,M_\odot$ NS and a $10\,M_\odot$ BH, using a few-body symplectic integrator. We used a range of impact parameters, orbital phases at the start of the simulati… ▽ More Aims: We aim to investigate the consequences of a fast massive stellar remnant - a black hole (BH) or a neutron star (NS) - encountering a planetary system. Methods: We modelled a close encounter between the actual Solar System (SS) and a $2\,M_\odot$ NS and a $10\,M_\odot$ BH, using a few-body symplectic integrator. We used a range of impact parameters, orbital phases at the start of the simulation derived from the current SS orbital parameters, encounter velocities, and incidence angles relative to the plane of the SS. Results: We give the distribution of possible outcomes, such as when the SS remains bound, when it suffers a partial or complete disruption, and in which cases the intruder is able to capture one or more planets, yielding planetary systems around a BH or a NS. We also show examples of the long-term stability of the captured planetary systems. △ Less

Submitted 7 April, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 25 pages, 22 figures 1 table, published in A&A Letters (Free access)

Journal ref: A&A 648, L2 (2021)

The aftermath of nova Cen 2013 (V1369 Cen)

Authors: Elena Mason, Steven N. Shore, Jeremy Drake, Steve B. Howell, Paul Kuin, Enza Magaudda

Abstract: Context: Classical nova progenitors are cataclysmic variables and very old novae are observed to match high mass transfer rate and (relatively) long orbital period systems. However, the aftermath of a classical nova has never been studied in detail. Aims: To probe the aftermath of a classical nova explosion in cataclysmic variables and observe as the binary system relaxes to quiescence. Methods: W… ▽ More Context: Classical nova progenitors are cataclysmic variables and very old novae are observed to match high mass transfer rate and (relatively) long orbital period systems. However, the aftermath of a classical nova has never been studied in detail. Aims: To probe the aftermath of a classical nova explosion in cataclysmic variables and observe as the binary system relaxes to quiescence. Methods: We used multi-wavelength time resolved optical and near-infrared spectroscopy for a bright, well studied classical nova five years after outburst. We were able to disentangle the contribution of the ejecta at this late epoch using its previous characterization, separating the ejecta emission from that of the binary system. Results: We determined the binary orbital period (P=3.76 hr), the system separation and mass ratio (q>=0.17 for an assumed white dwarf mass of 1.2 solar masses). We find evidence of an irradiated secondary star and no unambiguous signature of an accretion disk, although we identify a second emission line source tied to the white dwarf with an impact point. The data are consistent with a bloated white dwarf envelope and the presence of unsettled gas within the white dwarf Roche lobe. Conclusions: At more than 5 years after eruption, it appears that this classical nova has not yet relaxed. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: Accepted for publication in A&A. Main body of manuscript: 10 pages and 10 figures; 16 pages including appendix tables

Journal ref: A&A 649, A28 (2021)

An optical spectroscopic and polarimetric study of the microquasar binary system SS 433

Authors: P. Picchi, S. N. Shore, E. J. Harvey, A. Berdyugin

Abstract: We present a study of the mass transfer and wind outflows of SS433, focusing on the so-called stationary lines based on archival high and low resolution optical spectra, and new optical multifilter polarimetry and low resolution optical spectra spanning an interval of a decade and a broad range of precessional and orbital phases. We derive $\text{E(B-V)}=0.86\pm0.10$ and revised UV and U band pola… ▽ More We present a study of the mass transfer and wind outflows of SS433, focusing on the so-called stationary lines based on archival high and low resolution optical spectra, and new optical multifilter polarimetry and low resolution optical spectra spanning an interval of a decade and a broad range of precessional and orbital phases. We derive $\text{E(B-V)}=0.86\pm0.10$ and revised UV and U band polarizations and polarization angles that yield the same position angle as the optical. The polarization wavelength dependence is consistent with optical-dominating electron scattering with a Rayleigh component in U and the UV filters; no polarization changes were observed during a flare event. Using profile orbital and precessional modulation of multiple lines we derive properties for the accretion disk, present evidence for a strong disk wind, determine its velocity structure, and demonstrate its variability on timescales unrelated to the orbit. We derive a mass ratio $q=0.37\pm0.04$, and masses $\text{M}_X=4.2\pm0.4\ \text{M}_\odot$, $\text{M}_A=11.3\pm 0.6\ \text{M}_\odot$, and show that the A star fills its Roche surface. The O I 7772 Å and 8446 Å lines show different but related orbital modulation and no evidence for a circumbinary disk component. Instead, the spectral line profile variability can be understood with an ionization stratified outflow predicted by thermal wind modeling, which also accounts for an extended equatorial structure detected at long wavelength. △ Less

Submitted 6 August, 2020; v1 submitted 19 July, 2020; originally announced July 2020.

Comments: Typos corrected. Hyperref added. Accepted for publication in A&A (2020 May 27). Acknowledgements modified

Journal ref: A&A 640, A96 (2020)

The ambiguous transient ASASSN-17hx. A possible nova-impostor

Authors: Elena Mason, Steven N. Shore, Paul Kuin, Terry Bohlsen

Abstract: Some transients, although classified as novae based on their maximum and early decline optical spectra, cast doubts on their true nature and whether nova impostors might exist. We monitored a candidate nova which displayed a distinctly unusual light curve at maximum and early decline through optical spectroscopy (3000-10000 Å, 500<R<100000) complemented with Swift UV and AAVSO optical photometry.… ▽ More Some transients, although classified as novae based on their maximum and early decline optical spectra, cast doubts on their true nature and whether nova impostors might exist. We monitored a candidate nova which displayed a distinctly unusual light curve at maximum and early decline through optical spectroscopy (3000-10000 Å, 500<R<100000) complemented with Swift UV and AAVSO optical photometry. We use the spectral line series to characterize the ejecta dynamics, structure, and mass. We found that the ejecta are in free ballistic expansion and structured as typical of classical novae. However, their derived mass is at least an order of magnitude larger than the typical ejecta masses obtained for classical novae. Specifically, we found M$_{ej}\simeq$9$\times$10$^{-3}$ M$_\odot$ independent of the distance for a filling factor $\varepsilon$=1. By constraining the distance we derived $\varepsilon$ in the range 0.08-0.10, giving a mass 7$\times$10$^{-4}\lesssim$ M$_{ej}\lesssim$9$\times$10$^{-4}$ M$_\odot$. The nebular spectrum, characterized by unusually strong coronal emission lines, confines the ionizing source energy to the range 20-250 eV, possibly peaking in the range 75-100 or 75-150 eV. We link this source to other slow novae which showed similar behavior and suggest that they might form a distinct physical sub-group. They may result from a classical nova explosion occurring on a very low mass white dwarf or be impostors for an entirely different type of transient. △ Less

Submitted 31 January, 2020; originally announced January 2020.

Comments: accepted for publication on A&A

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

123-321 Models of Classical Novae

Authors: Jordi Jose, Steven N. Shore, Jordi Casanova

Abstract: High-resolution spectroscopy has revealed large concentrations of CNO and sometimes other intermediate-mass elements in the shells ejected during nova outbursts, suggesting that the solar composition material transferred from the secondary mixes with the outermost layers of the underlying white dwarf during the thermonuclear runaway. Multidimensional simulations have shown that Kelvin-Helmholtz in… ▽ More High-resolution spectroscopy has revealed large concentrations of CNO and sometimes other intermediate-mass elements in the shells ejected during nova outbursts, suggesting that the solar composition material transferred from the secondary mixes with the outermost layers of the underlying white dwarf during the thermonuclear runaway. Multidimensional simulations have shown that Kelvin-Helmholtz instabilities provide self-enrichment of the accreted envelope with material from the outermost layers of the white dwarf, at levels that agree with observations. However, the Eulerian and time-explicit nature of most multidimensional codes used to date and the overwhelming computational load have limited their applicability, and no multidimensional simulation has been conducted for a full nova cycle. This paper explores a new methodology that combines 1-D and 3-D simulations. The early stages of the explosion (i.e., mass-accretion and initiation of the runaway) have been computed with the 1-D hydrodynamic code SHIVA. When convection extends throughout the entire envelope, the structures for each model were mapped into 3-D Cartesian grids and were subsequently followed with the multidimensional code FLASH. Two key physical quantities were extracted from the 3-D simulations and subsequently implemented into SHIVA, which was used to complete the simulation through the late expansion and ejection stages: the time-dependent amount of mass dredged-up from the outer white dwarf layers, and the time-dependent convective velocity profile throughout the envelope. More massive envelopes than those reported from previous models with pre-enrichment have been found. This results in more violent outbursts, characterized by higher peak temperatures and greater ejected masses, with metallicity enhancements in agreement with observations. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: 8 pages, accepted for publication in Astronomy & Astrophysics

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

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

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

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

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

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

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

The January 2016 eruption of recurrent nova LMC 1968

Authors: N. P. M. Kuin, K. L. Page, P. Mróz, M. J. Darnley, S. N. Shore, J. P. Osborne, F. Walter, F. Di Mille, N. Morrell, U. Munari, T. Bohlsen, A. Evans, R. D. Gehrz, S. Starrfield, M. Henze, S. C. Williams, G. J. Schwarz, A. Udalski, M. K. Szymański, R. Poleski, I. Soszyński, V. A. R. M. Ribeiro, R. Angeloni, A. A. Breeveld, A. P. Beardmore , et al. (1 additional authors not shown)

Abstract: We present a comprehensive review of all observations of the eclipsing recurrent Nova LMC 1968 in the Large Magellanic Cloud which was previously observed in eruption in 1968, 1990, 2002, 2010, and most recently in 2016. We derive a probable recurrence time of $6.2 \pm 1.2$ years and provide the ephemerides of the eclipse. In the ultraviolet-optical-IR photometry the light curve shows high variabi… ▽ More We present a comprehensive review of all observations of the eclipsing recurrent Nova LMC 1968 in the Large Magellanic Cloud which was previously observed in eruption in 1968, 1990, 2002, 2010, and most recently in 2016. We derive a probable recurrence time of $6.2 \pm 1.2$ years and provide the ephemerides of the eclipse. In the ultraviolet-optical-IR photometry the light curve shows high variability right from the first observation around two days after eruption. Therefore no colour changes can be substantiated. Outburst spectra from 2016 and 1990 are very similar and are dominated by H and He lines longward of 2000 Angstrom. Interstellar reddening is found to be E(B-V) = $0.07\pm0.01$. The super soft X-ray luminosity is lower than the Eddington luminosity and the X-ray spectra suggest the mass of the WD is larger than 1.3 M$_\odot$. Eclipses in the light curve suggest that the system is at high orbital inclination. On day four after the eruption a recombination wave was observed in Fe II ultraviolet absorption lines. Narrow line components are seen after day 6 and explained as being due to reionisation of ejecta from a previous eruption. The UV spectrum varies with orbital phase, in particular a component of the He II 1640 Angstrom emission line, which leads us to propose that early-on the inner WD Roche lobe might be filled with a bound opaque medium prior to the re-formation of an accretion disk. Both this medium and the ejecta can cause the delay in the appearance of the soft X-ray source. △ Less

Submitted 20 October, 2019; v1 submitted 7 September, 2019; originally announced September 2019.

Comments: MNRAS accepted October 18, 2019

Structure and Expansion Law of HII Regions in structured Molecular Clouds

Authors: Manuel Zamora-Avilés, Enrique Vázquez-Semadeni, Ricardo F. González, José Franco, Steven N. Shore, Lee W. Hartmann, Javier Ballesteros-Paredes, Robi Banerjee, Bastian Körtgen

Abstract: We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventu… ▽ More We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^α$, with $α$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models. △ Less

Submitted 17 May, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS

Observatory science with eXTP

Authors: Jean J. M. in 't Zand, Enrico Bozzo, Jinlu Qu, Xiang-Dong Li, Lorenzo Amati, Yang Chen, Immacolata Donnarumma, Victor Doroshenko, Stephen A. Drake, Margarita Hernanz, Peter A. Jenke, Thomas J. Maccarone, Simin Mahmoodifar, Domitilla de Martino, Alessandra De Rosa, Elena M. Rossi, Antonia Rowlinson, Gloria Sala, Giulia Stratta, Thomas M. Tauris, Joern Wilms, Xuefeng Wu, Ping Zhou, Iván Agudo, Diego Altamirano , et al. (159 additional authors not shown)

Abstract: In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to stu… ▽ More In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

The extreme HBL behaviour of Markarian 501 during 2012

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

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

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

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

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

Two-dimensional simulations of mixing in classical novae: the effect of the white dwarf composition and mass

Authors: Jordi Casanova, Jordi Jose, Steven N. Shore

Abstract: Context. Classical novae are explosive phenomena that take place in stellar binary systems. They are powered by mass transfer from a low-mass main sequence star onto a white dwarf (either CO or ONe). The material accumulates for 10+4 - 10+5 yr until ignition under degenerate conditions, resulting in a thermonuclear runaway. The nuclear energy released produces peak temperatures of about 0.1 - 0.4… ▽ More Context. Classical novae are explosive phenomena that take place in stellar binary systems. They are powered by mass transfer from a low-mass main sequence star onto a white dwarf (either CO or ONe). The material accumulates for 10+4 - 10+5 yr until ignition under degenerate conditions, resulting in a thermonuclear runaway. The nuclear energy released produces peak temperatures of about 0.1 - 0.4 GK. During these events, 10-7 - 10-3 Msun enriched in intermediate-mass elements (with respect to solar abundances) are ejected into the interstellar medium. However, the origin of the large metallicity enhancements and the inhomogeneous distribution of chemical species observed in high-resolution spectra of ejected nova shells is not fully understood. Aims. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities that operate at the core-envelope interface can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, such multidimensional simulations have been performed for a small number of cases, and much of the parameter space remains unexplored. Methods. Here we investigate the dredge-up, driven by Kelvin-Helmholtz instabilities, for white dwarf masses in the range 0.8-1.25 Msun and different core compositions (CO-rich and ONe-rich substrates). We present a set of five numerical simulations performed in two dimensions aimed at analyzing the possible impact of the white dwarf mass (and composition) on the metallicity enhancement and on the explosion characteristics. Results. We observe greater mixing (about 30 per cent higher when measured in the same conditions), at the time we stop the simulations, and more energetic outbursts for ONe-rich substrates than for CO-rich substrates and for more massive white dwarfs. △ Less

Submitted 27 July, 2018; originally announced July 2018.

Comments: Accepted for Publication in Astronomy and Astrophysics

V1369 Cen high resolution panchromatic late nebular spectra in the context of a unified picture for nova ejecta

Authors: Elena Mason, Steven N. Shore, Ivan De Gennaro Aquino, Luca Izzo, Kim Page, Greg J. Schwarz

Abstract: Nova Cen 2013 (V1369 Cen) is the fourth bright nova observed panchromatically through high resolution UV+optical multi epoch spectroscopy. It is also the nova with the richest set of spectra (both in terms of data quality and number of epochs) thanks to its exceptional brightness. Here, we use the late nebular spectra taken between day ~250 and day ~837 after outburst to derive the physical, geome… ▽ More Nova Cen 2013 (V1369 Cen) is the fourth bright nova observed panchromatically through high resolution UV+optical multi epoch spectroscopy. It is also the nova with the richest set of spectra (both in terms of data quality and number of epochs) thanks to its exceptional brightness. Here, we use the late nebular spectra taken between day ~250 and day ~837 after outburst to derive the physical, geometrical and kinematical properties of the nova. We compare the results with those determined for the other panchromatic studies in this series: T Pyx, V339 Del (nova Del 2013), and V959 Mon (nova Mon 2012). From this we conclude that in all these novae the ejecta geometry and phenomenology can be consistently explained by clumpy gas expelled during a single, brief ejection episode and in ballistic expansion, and not by a wind. For V1369 Cen the ejecta mass (about 1E-4 solar masses) and filling factor (0.1<=f<=0.2) are consistent with those of classical novae but larger (by at least an order of magnitude) than those of T Pyx and the recurrent novae. V1369 Cen has an anomalously high relative to solar N/C ratio that is beyond the range currently predicted for a CO nova, and the Ne emission line strengths are dissimilar to those of typical ONe or CO white dwarfs. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Journal ref: Published, 2018, ApJ, 853, 27

Spectroscopic diagnostics of dust formation and evolution in classical nova ejecta

Authors: Steven N. Shore, N. Paul Kuin, Elena Mason, Ivan De Gennaro Aquino

Abstract: A fraction of classical novae form dust during the early stages of their outbursts. The classical CO nova V5668 Sgr (Nova Sgr. 2015b) underwent a deep photometric minimum about 100 days after outburst that was covered across the spectrum. A similar event was observed for an earlier CO nova, V705 Cas (Nova Cas 1993) and a less optically significant event for the more recent CO nova V339 Del (Nova D… ▽ More A fraction of classical novae form dust during the early stages of their outbursts. The classical CO nova V5668 Sgr (Nova Sgr. 2015b) underwent a deep photometric minimum about 100 days after outburst that was covered across the spectrum. A similar event was observed for an earlier CO nova, V705 Cas (Nova Cas 1993) and a less optically significant event for the more recent CO nova V339 Del (Nova Del 2013). This study provides a "compare and contrast" of these events to better understand the very dynamical event of dust formation. We show the effect of dust formation on multiwavelength high resolution line profiles in the interval 1200Å - 9200Å using a biconical ballistic structure that has been applied in our previous studies of the ejecta. We find that both V5668 Sgr and V339 Del can be modeled using a grey opacity for the dust, indicating fairly large grains (at least 0.1 micron) and that the persistent asymmetries of the line profiles in late time spectra, up to 650 days after the event for V5668 Sgr and 866 days for V339 Del, point to the survival of the dust well into the transparent, nebular stage of the ejecta evolution. This is a general method for assessing the properties of dust forming novae well after the infrared is completely transparent in the ejecta. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: 15 pages 14 figures, accepted for publication in A&A, 2018 June 28

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

The Temporal Development of Dust Formation and Destruction in Nova Sagittarii 2015#2 (V5668 Sgr): A Panchromatic Study

Authors: R. D. Gehrz, A. Evans, C. E. Woodward, L. A. Helton, D. P. K. Banerjee, M. K. Srivastava, N. M. Ashok, V. Joshi, S. P. S. Eyres, Joachim Krautter, N. P. M. Kuin, K. L. Page, J. P. Osborne, G. J. Schwarz, D. P. Shenoy, S. N. Shore, S. G. Starrfield, R. M. Wagner

Abstract: We present 5-28 micron SOFIA FORECAST spectroscopy complemented by panchromatic X-ray through infrared observations of the CO nova V5668 Sgr documenting the formation and destruction of dust during 500 days following outburst. Dust condensation commenced by 82 days after outburst at a temperature of 1090 K. The condensation temperature indicates that the condensate was amorphous carbon. There was… ▽ More We present 5-28 micron SOFIA FORECAST spectroscopy complemented by panchromatic X-ray through infrared observations of the CO nova V5668 Sgr documenting the formation and destruction of dust during 500 days following outburst. Dust condensation commenced by 82 days after outburst at a temperature of 1090 K. The condensation temperature indicates that the condensate was amorphous carbon. There was a gradual decrease of the grain size and dust mass during the recovery phase. Absolute parameter values given here are for an assumed distance of 1.2 kpc. We conclude that the maximum mass of dust produced was 1.2 x 10-7 solar masses if the dust was amorphous carbon. The average grain radius grew to a maximum of 2.9 microns at a temperature of 720 K around day 113 when the shell visual optical depth was Tau = 5.4. Maximum grain growth was followed by followed by a period of grain destruction. X-rays were detected with Swift from day 95 to beyond day 500. The Swift X-ray count rate due to the hot white dwarf peaked around day 220, when its spectrum was that of a kT = 35 eV blackbody. The temperature, together with the super-soft X-ray turn-on and turn-off times, suggests a WD mass of 1.1 solar masses. We show that the X-ray fluence was sufficient to destroy the dust. Our data show that the post-dust event X-ray brightening is not due to dust destruction, which certainly occurred, as the dust is optically thin to X-rays. △ Less

Submitted 24 April, 2018; v1 submitted 2 April, 2018; originally announced April 2018.

Comments: Accepted for publication in the Astrophysical Journal. This is the final version after correction of the proofs

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

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

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

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

Comments: HTML page with links

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

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

Abstract: The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capabili… ▽ More The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses and precession of the circumstellar disk periodically covering the central binary system is expected to be at its minimum. The eastern and western SS433/W50 interaction regions are also examined. We aim to constrain some theoretical models previously developed for this system. We made use of dedicated observations from MAGIC and H.E.S.S. from 2006 to 2011 which were combined for the first time and accounted for a total effective observation time of 16.5 h. Gamma-ray attenuation does not affect the jet/medium interaction regions. The analysis of a larger data set amounting to 40-80 h, depending on the region, was employed. No evidence of VHE gamma-ray emission was found. Upper limits were computed for the combined data set. We place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Our findings suggest that the fraction of the jet kinetic power transferred to relativistic protons must be relatively small to explain the lack of TeV and neutrino emission from the central system. At the SS433/W50 interface, the presence of magnetic fields greater 10$μ$G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with energies up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions. △ Less

Submitted 12 July, 2017; originally announced July 2017.

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

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

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

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

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

Submitted 5 July, 2017; originally announced July 2017.

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

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

Periodicity makes galactic shocks unstable - I. Linear analysis

Authors: Mattia C. Sormani, Emanuele Sobacchi, Steven N. Shore, Robin G. Tress, Ralf S. Klessen

Abstract: We study the dynamical stability of stationary galactic spiral shocks. The steady-state equilibrium flow contains a shock of the type derived by Roberts in the tightly wound approximation. We find that boundary conditions are critical in determining whether the solutions are stable or not. Shocks are unstable if periodic boundary conditions are imposed. For intermediate strengths of the spiral pot… ▽ More We study the dynamical stability of stationary galactic spiral shocks. The steady-state equilibrium flow contains a shock of the type derived by Roberts in the tightly wound approximation. We find that boundary conditions are critical in determining whether the solutions are stable or not. Shocks are unstable if periodic boundary conditions are imposed. For intermediate strengths of the spiral potential, the instability disappears if boundary conditions are imposed such that the upstream flow is left unperturbed as in the classic analysis of D'yakov and Kontorovich. This reconciles apparently contradictory findings of previous authors regarding the stability of spiral shocks. This also shows that the instability is distinct from the Kelvin-Helmholtz instability, confirming the findings of Kim et al. We suggest that instability is a general characteristics of periodic shocks, regardless of the presence of shear, and provide a physical picture as to why this is the case. For strong spiral potentials, high post-shock shear makes the system unstable also to parasitic Kelvin-Helmholtz instability regardless of the boundary conditions. Our analysis is performed in the context of a simplified problem that, while preserving all the important characteristics of the original problem, strips it from unnecessary complications, and assumes that the gas is isothermal, non self-gravitating, non-magnetised. △ Less

Submitted 2 July, 2017; originally announced July 2017.

Comments: Accepted for publication in MNRAS

A Systematic Survey of the Effects of Wind Mass Loss Algorithms on the Evolution of Single Massive Stars

Authors: Mathieu Renzo, Christian D. Ott, Steven N. Shore, Selma E. de Mink

Abstract: Mass loss is a key uncertainty in the evolution of massive stars. Stellar evolution calculations must employ parametric algorithms for mass loss, and usually only include stellar winds. We carry out a parameter study of the effects of wind mass loss on massive star evolution using the open-source stellar evolution code MESA. We provide a systematic comparison of wind mass loss algorithms for solar… ▽ More Mass loss is a key uncertainty in the evolution of massive stars. Stellar evolution calculations must employ parametric algorithms for mass loss, and usually only include stellar winds. We carry out a parameter study of the effects of wind mass loss on massive star evolution using the open-source stellar evolution code MESA. We provide a systematic comparison of wind mass loss algorithms for solar-metallicity, nonrotating, single stars in the initial mass range of $15-35\,M_\odot$. We consider combinations drawn from two hot phase algorithms, three cool phase algorithms, and two Wolf-Rayet algorithms. We consider linear wind efficiency scale factors of $1$, $0.33$, and $0.1$ to account for reductions in mass loss rates due to wind inhomogeneities. We find that the initial to final mass mapping for each zero-age main-sequence (ZAMS) mass has a $\sim 50\%$ uncertainty if all algorithm combinations and wind efficiencies are considered. The ad-hoc efficiency scale factor dominates this uncertainty. While the final total mass and internal structure of our models vary tremendously with mass loss treatment, final observable parameters are much less sensitive for ZAMS mass $\lesssim 30\,M_\odot$. This indicates that uncertainty in wind mass loss does not negatively affect estimates of the ZAMS mass of most single-star supernova progenitors from pre-explosion observations. Furthermore, we show that the internal structure of presupernova stars is sensitive to variations in both main sequence and post main-sequence mass loss. We find that the compactness parameter $ξ\propto M/R(M)$ varies by as much as $30\%$ for a given ZAMS mass evolved with different wind efficiencies and mass loss algorithm combinations. [abridged] △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: Accepted for publication on A&A, 22 pages + 2 appendixes, 12 figures, online input parameters available at https://stellarcollapse.org/renzo2017 and data at https://zenodo.org/record/292924#.WK0q2tWi6Wg

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

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

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

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

Submitted 22 March, 2017; originally announced March 2017.

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

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

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

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

Abstract: We present the first detection of the nearby (z=0.084) low-luminosity BL Lac object 1ES 1741+196 in the very high energy (VHE: E$>$100 GeV) band. This object lies in a triplet of interacting galaxies. Early predictions had suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac sources, within the reach of MAGIC detectability. Its detection by MAGIC, later confirmed by VERITAS… ▽ More We present the first detection of the nearby (z=0.084) low-luminosity BL Lac object 1ES 1741+196 in the very high energy (VHE: E$>$100 GeV) band. This object lies in a triplet of interacting galaxies. Early predictions had suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac sources, within the reach of MAGIC detectability. Its detection by MAGIC, later confirmed by VERITAS, helps to expand the small population of known TeV BL Lacs. The source was observed with the MAGIC telescopes between 2010 April and 2011 May, collecting 46 h of good quality data. These observations led to the detection of the source at 6.0 $σ$ confidence level, with a steady flux $\mathrm{F}(> 100 {\rm GeV}) = (6.4 \pm 1.7_{\mathrm{stat}}\pm 2.6_{\mathrm{syst}}) \cdot 10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and a differential spectral photon index $Γ= 2.4 \pm 0.2_{\mathrm{stat}} \pm 0.2_{\mathrm{syst}}$ in the range of $\sim$80 GeV - 3 TeV. To study the broad-band spectral energy distribution (SED) simultaneous with MAGIC observations, we use KVA, Swift/UVOT and XRT, and Fermi/LAT data. One-zone synchrotron-self-Compton (SSC) modeling of the SED of 1ES 1741+196 suggests values for the SSC parameters that are quite common among known TeV BL Lacs except for a relatively low Doppler factor and slope of electron energy distribution. A thermal feature seen in the SED is well matched by a giant elliptical's template. This appears to be the signature of thermal emission from the host galaxy, which is clearly resolved in optical observations. △ Less

Submitted 22 February, 2017; originally announced February 2017.

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

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

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

Abstract: We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux… ▽ More We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties. △ Less

Submitted 30 December, 2016; originally announced December 2016.

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

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

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

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

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

Submitted 21 November, 2016; originally announced November 2016.

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

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

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

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

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

Submitted 14 November, 2016; originally announced November 2016.

Comments: 45 pages, 15 figures

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

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

Abstract: Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In… ▽ More Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models. △ Less

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

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

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

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

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

Abstract: Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components, spatially indistinguishable by gamma-ray instruments, but separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energ… ▽ More Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components, spatially indistinguishable by gamma-ray instruments, but separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energy distribution of QSO B0218+357 can give information on the energetics of z ~ 1 very high energy gamma- ray sources. Moreover the gamma-ray emission can also be used as a probe of the extragalactic background light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during the expected arrival time of the delayed component of the emission. The MAGIC and Fermi-LAT observations were accompanied by quasi-simultaneous optical data from the KVA telescope and X-ray observations by Swift-XRT. We construct a multiwavelength spectral energy distribution of QSO B0218+357 and use it to model the source. The GeV and sub-TeV data, obtained by Fermi-LAT and MAGIC, are used to set constraints on the extragalactic background light. Results. Very high energy gamma-ray emission was detected from the direction of QSO B0218+357 by the MAGIC telescopes during the expected time of arrival of the trailing component of the flare, making it the farthest very high energy gamma-ray sources detected to date. The observed emission spans the energy range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy distribution of QSO B0218+357 is consistent with current extragalactic background light models. The broad band emission can be modeled in the framework of a two zone external Compton scenario, where the GeV emission comes from an emission region in the jet, located outside the broad line region. △ Less

Submitted 5 September, 2016; originally announced September 2016.

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

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

Three-dimensional simulations of turbulent convective mixing in ONe and CO classical nova explosions

Authors: Jordi Casanova, Jordi Jose, Enrique Garcia-Berro, Steven N. Shore

Abstract: Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures about 100 - 400 MK. During these events, about 10-3 - 10-7 Msun, enric… ▽ More Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures about 100 - 400 MK. During these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the large concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. Three-dimensional simulations of mixing at the core-envelope interface during nova outbursts have been performed with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We show that the presence of an ONe-rich substrate, as in "neon novae", yields larger metallicity enhancements in the ejecta, compared to CO,rich substrates (i.e., non-neon novae). A number of requirements and constraints for such 3-D simulations (e.g., minimum resolution, size of the computational domain) are also outlined. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: Accepted for publication in Astronomy & Astrophysics

The panchromatic spectroscopic evolution of the classical CO nova V339 Del (Nova Del 2013) until X-ray turnoff

Authors: S. N. Shore, E. Mason, G. J. Schwarz, F. M. Teyssier, C. Buil, I. De Gennaro Aquino, K. L. Page, J. P. Osborne, S. Scaringi, S. Starrfield, H. van Winckel, R. E. Williams, C. E. Woodward

Abstract: Classical novae are the product of thermonuclear runaway-initiated explosions occurring on accreting white dwarfs. V339 Del (Nova Delphinus 2013) was one of the brightest classical novae of the last hundred years. Spectroscopy and photometry are available from $γ$-rays through infrared at stages that have frequently not been well observed. The complete data set is intended to provide a benchmark f… ▽ More Classical novae are the product of thermonuclear runaway-initiated explosions occurring on accreting white dwarfs. V339 Del (Nova Delphinus 2013) was one of the brightest classical novae of the last hundred years. Spectroscopy and photometry are available from $γ$-rays through infrared at stages that have frequently not been well observed. The complete data set is intended to provide a benchmark for comparison with modeling and for understanding more sparsely monitored historical classical and recurrent novae. This paper is the first in the series of reports on the development of the nova. We report here on the early stages of the outburst, through the X-ray active stage. A time sequence of optical, flux calibrated high resolution spectra was obtained with the Nordic Optical Telescope (NOT) using FIES simultaneously, or contemporaneously, with the Space Telescope Imaging Spectrograph (STIS) aboard the Hubble Space Telescope during the early stages of the outburst. These were supplemented with MERCATOR/HERMES optical spectra. High resolution IUE ultraviolet spectra of OS And 1986, taken during the Fe curtain phase, served as a template for the distance determination. We used standard plasma diagnostics (e.g., [O III] and [N II] line ratios, and the H$β$ line flux) to constrain electron densities and temperatures of the ejecta. Using Monte Carlo modeling of the ejecta, we derived the structure, filling factor, and mass from comparisons of the optical and ultraviolet line profiles. △ Less

Submitted 6 June, 2016; originally announced June 2016.

Comments: 14 pages, 18 figures

Journal ref: 2016 A&A, 590, 123

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

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

Abstract: We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux v… ▽ More We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Metsähovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Metsähovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades. △ Less

Submitted 29 May, 2016; originally announced May 2016.

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

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

Fermi LAT Gamma-ray Detections of Classical Novae V1369 Centauri 2013 and V5668 Sagittarii 2015

Authors: C. C. Cheung, P. Jean, S. N. Shore, L. Stawarz, R. H. D. Corbet, J. Knodlseder, S. Starrfield, D. L. Wood, R. Desiante, F. Longo, G. Pivato, K. S. Wood

Abstract: We report the Fermi Large Area Telescope (LAT) detections of high-energy (>100 MeV) gamma-ray emission from two recent optically bright classical novae, V1369 Centauri 2013 and V5668 Sagittarii 2015. At early times, Fermi target-of-opportunity observations prompted by their optical discoveries provided enhanced LAT exposure that enabled the detections of gamma-ray onsets beginning ~2 days after th… ▽ More We report the Fermi Large Area Telescope (LAT) detections of high-energy (>100 MeV) gamma-ray emission from two recent optically bright classical novae, V1369 Centauri 2013 and V5668 Sagittarii 2015. At early times, Fermi target-of-opportunity observations prompted by their optical discoveries provided enhanced LAT exposure that enabled the detections of gamma-ray onsets beginning ~2 days after their first optical peaks. Significant gamma-ray emission was found extending to 39-55 days after their initial LAT detections, with systematically fainter and longer duration emission compared to previous gamma-ray detected classical novae. These novae were distinguished by multiple bright optical peaks that encompassed the timespans of the observed gamma rays. The gamma-ray light curves and spectra of the two novae are presented along with representative hadronic and leptonic models, and comparisons to other novae detected by the LAT are discussed. △ Less

Submitted 13 May, 2016; originally announced May 2016.

Comments: 13 pages, 6 figures, 4 tables, ApJ accepted

Fermi Reveals New Light on Novae in Gamma rays

Authors: C. C. Cheung, P. Jean, S. N. Shore, J. E. Grove, M. Leising

Abstract: Novae are now firmly established as a high-energy (>100 MeV) gamma-ray source class by the Fermi Large Area Telescope (LAT). In symbiotic binary systems such as V407 Cyg 2010, there is a firm theoretical framework for the production of shock-accelerated particles in the nova ejecta from interactions with the dense wind of the red giant companion. Yet, the high-energy gamma-ray emission detected in… ▽ More Novae are now firmly established as a high-energy (>100 MeV) gamma-ray source class by the Fermi Large Area Telescope (LAT). In symbiotic binary systems such as V407 Cyg 2010, there is a firm theoretical framework for the production of shock-accelerated particles in the nova ejecta from interactions with the dense wind of the red giant companion. Yet, the high-energy gamma-ray emission detected in classical novae involving less evolved stellar companions cannot be explained in the same way and could instead be produced in internal shocks in the ejecta. We summarize the Fermi-LAT gamma-ray observations of novae, highlighting the main properties that will guide further studies. Additionally, we report on the soft gamma-ray (~0.1 MeV) continuum detection of the oxygen-neon type classical nova V382 Vel 1999 with the OSSE detector aboard the Compton Gamma Ray Observatory in light of its Fermi-era analog, V959 Mon 2012. △ Less

Submitted 4 May, 2016; originally announced May 2016.

Comments: 9 pages, 5 figures. The 34th International Cosmic Ray Conference, 30 July - 6 August 2015 (The Hague, The Netherlands), PoS(ICRC2015) 880

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

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

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

Submitted 23 March, 2016; originally announced March 2016.

Comments: 13 pages, accepted for publication in MNRAS

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

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

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

Abstract: The gamma-ray binary LS I +61$^{\circ}$303 is a well established source from centimeter radio up to very high energy (VHE; E$>$100 GeV). Its broadband emission shows a periodicity of $\sim$26.5 days, coincident with the orbital period. A longer (super-orbital) period of 1667 $\pm$ 8 days was discovered in radio and confirmed in optical and high energy (HE; E>100 MeV) gamma-ray observations. We pre… ▽ More The gamma-ray binary LS I +61$^{\circ}$303 is a well established source from centimeter radio up to very high energy (VHE; E$>$100 GeV). Its broadband emission shows a periodicity of $\sim$26.5 days, coincident with the orbital period. A longer (super-orbital) period of 1667 $\pm$ 8 days was discovered in radio and confirmed in optical and high energy (HE; E>100 MeV) gamma-ray observations. We present a four-year campaign performed by MAGIC together with archival data concentrating on a search for a long timescale signature in the VHE emission. We focus on the search for super-orbital modulation of the VHE peak and on the search for correlations between TeV emission and optical determination of the extension of the circumstellar disk. A four-year campaign has been carried out by MAGIC. The source was observed during the orbital phases when the periodic VHE outbursts have occurred ($φ$=0.55-0.75). Additionally, we included archival MAGIC observations and data published by the VERITAS collaboration in these studies. For the correlation studies, LS I +61$^{\circ}$303 has also been observed during the orbital phases where sporadic VHE emission had been detected in the past ($φ$=0.75-1.0). These MAGIC observations were simultaneous with optical spectroscopy from the LIVERPOOL telescope. The TeV flux of the periodical outburst in orbital phases $φ$=0.5--0.75 was found to show yearly variability consistent with the $\sim$4.5 years long-term modulation found in the radio band. This modulation of the TeV flux can be well described by a sine function with the best fit period of $1610\pm 58$ days. The complete dataset span two super-orbital periods. There is no evidence for a correlation between the TeV emission and the mass-loss rate of the Be star but this may be affected by the strong, short timescale (as short as intra-day) variation displayed by the H$α$ fluxes. △ Less

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

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

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

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

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

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

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

Comments: 14 pages, 7 figures

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

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

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

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

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

Comments: 8 pages, 6 figures

MAGIC observations of the February 2014 flare of 1ES 1011+496 and ensuing constraint of the EBL density

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

Abstract: In February-March 2014, the MAGIC telescopes observed the high-frequency peaked BL Lac 1ES 1011+496 (z=0.212) in flaring state at very-high energy (VHE, E>100GeV). The flux reached a level more than 10 times higher than any previously recorded flaring state of the source. We present the description of the characteristics of the flare presenting the light curve and the spectral parameters of the ni… ▽ More In February-March 2014, the MAGIC telescopes observed the high-frequency peaked BL Lac 1ES 1011+496 (z=0.212) in flaring state at very-high energy (VHE, E>100GeV). The flux reached a level more than 10 times higher than any previously recorded flaring state of the source. We present the description of the characteristics of the flare presenting the light curve and the spectral parameters of the night-wise spectra and the average spectrum of the whole period. From these data we aim at detecting the imprint of the Extragalactic Background Light (EBL) in the VHE spectrum of the source, in order to constrain its intensity in the optical band. For this we implement the method developed by the H.E.S.S. collaboration in which the intrinsic energy spectrum of the source is modeled with a simple function, and the EBL-induced optical depth is calculated using a template EBL model. The likelihood of the observed spectrum is then maximized, including a normalization factor for the EBL opacity among the free parameters. From the data collected differential energy spectra was produced for all nights of the observed period. Evaluating the changes in the fit parameters we conclude that the spectral shape for most of the nights were compatible, regardless of the flux level, which enabled us to produce an average spectrum from which the EBL imprint could be constrained. The likelihood ratio test shows that the model with an EBL density 1.07(-0.20,+0.24)_{stat+sys}, relative to the one in the tested EBL template (Dominguez et al.2011), is preferred at the 4.6 sigma level to the no-EBL hypothesis, with the assumption that the intrinsic source spectrum can be modeled as a log-parabola. This would translate into a constraint of the EBL density in the wavelength range [0.24 um,4.25 um], with a peak value at 1.4 um of F=12.27_{-2.29}^{+2.75} nW m^{-2} sr^{-1}, including systematics. △ Less

Submitted 18 February, 2016; v1 submitted 16 February, 2016; originally announced February 2016.

Comments: 8 pages, 8 figures, updated references in section 1

Journal ref: A&A 590, A24 (2016)

Limits to dark matter annihilation cross-section from a combined analysis of MAGIC and Fermi-LAT observations of dwarf satellite galaxies

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

Abstract: We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi-LAT. We obtain limits on the annihilation cross-secti… ▽ More We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi-LAT. We obtain limits on the annihilation cross-section for dark matter particle masses between 10 GeV and 100 TeV - the widest mass range ever explored by a single gamma-ray analysis. These limits improve on previously published Fermi-LAT and MAGIC results by up to a factor of two at certain masses. Our new inclusive analysis approach is completely generic and can be used to perform a global, sensitivity-optimized dark matter search by combining data from present and future gamma-ray and neutrino detectors. △ Less

Submitted 16 February, 2016; v1 submitted 25 January, 2016; originally announced January 2016.

Comments: 19 pages, 3 figures. V2: Few typos corrected and references added. Matches published version JCAP 02 (2016) 039

Journal ref: JCAP 02 (2016) 039

Very-high-energy gamma-rays from the Universe's middle age: detection of the z=0.940 blazar PKS 1441+25 with MAGIC

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

Abstract: The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies… ▽ More The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy. △ Less

Submitted 12 January, 2018; v1 submitted 14 December, 2015; originally announced December 2015.

Comments: Corresponding Authors: J. Becerra (josefa.becerra@nasa.gov), M. Nievas Rosillo (mnievas@ucm.es), M. Manganaro (manganaro@iac.es), F. Tavecchio (fabrizio.tavecchio@brera.inaf.it) Published in The Astrophysical Journal Letters (ApJL), 2015ApJ...815L..23A, DOI: 10.1088/2041-8205/815/2/L23

Journal ref: ApJL, 2015, 815, L23

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

Authors: M. Baloković, D. Paneque, G. Madejski, A. Furniss, J. Chiang, the NuSTAR team, :, M. Ajello, D. M. Alexander, D. Barret, R. Blandford, S. E. Boggs, F. E. Christensen, W. W. Craig, K. Forster, P. Giommi, B. W. Grefenstette, C. J. Hailey, F. A. Harrison, A. Hornstrup, T. Kitaguchi, J. E. Koglin, K. K. Madsen, P. H. Mao, H. Miyasaka , et al. (286 additional authors not shown)

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

Submitted 7 December, 2015; originally announced December 2015.

Comments: 32 pages, 14 figures; accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 819, 156 (2016)

Teraelectronvolt pulsed emission from the Crab pulsar detected by MAGIC

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

Abstract: Aims: To investigate the extension of the very-high-energy spectral tail of the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours of good quality data of Crab with the MAGIC telescope, obtained from February 2007 until April 2014. Results: We report the most energetic pulsed emission ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile shows two na… ▽ More Aims: To investigate the extension of the very-high-energy spectral tail of the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours of good quality data of Crab with the MAGIC telescope, obtained from February 2007 until April 2014. Results: We report the most energetic pulsed emission ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with the ones measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC telescopes the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least $5\times 10^6$. These results strongly suggest IC scattering off low energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder. △ Less

Submitted 5 December, 2018; v1 submitted 23 October, 2015; originally announced October 2015.

First NuSTAR Observations of Mrk 501 within a Radio to TeV Multi-Instrument Campaign

Authors: A. Furniss, K. Noda, S. Boggs, J. Chiang, F. Christensen, W. Craig, P . Giommi, C. Hailey, F. Harisson, G. Madejski, K. Nalewajko, M. Perri, D. Stern, M. Urry, F. Verrecchia, W. Zhang, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio , et al. (251 additional authors not shown)

Abstract: We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the F… ▽ More We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsähovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission. △ Less

Submitted 24 September, 2015; v1 submitted 16 September, 2015; originally announced September 2015.

Comments: Accepted for publication in ApJ