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.}

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 \AA, 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.

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.

The investigation of the emergence of life is a major endeavour of science. Astronomy is contributing to it in three fundamental manners: (1) by measuring the chemical enrichment of the Universe, (2) by investigating planet formation and searching for exoplanets with signatures of life and, (3) by determining the abundance of aminoacids and the chemical routes to aminoacid and protein growth in astronomical bodies. This proposal deals with the first two. In the Voyage to 2050, the world-wide scientific community is getting equipped with large facilities for the investigation of the emergence of life in the Universe (i.e. VLT, JWST, ELT, GMT, TMT, ALMA, FAST, VLA, ATHENA, SKA) including the ESA's CHEOPS, PLATO and ARIEL missions. This white paper is a community effort to call for the development of a large ultraviolet optical observatory to gather fundamental data for this investigation that will not be accessible through other ranges of the electromagnetic spectrum. A versatile space observatory with UV sensitivity a factor of 50-100 greater than existing facilities will revolutionize our understanding of the pathway to life in the Universe.

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.

Current burning issues in stellar physics, for both hot and cool stars, concern their magnetism. In hot stars, stable magnetic fields of fossil origin impact their stellar structure and circumstellar environment, with a likely major role in stellar evolution. However, this role is complex and thus poorly understood as of today. It needs to be quantified with high-resolution UV spectropolarimetric measurements. In cool stars, UV spectropolarimetry would provide access to the structure and magnetic field of the very dynamic upper stellar atmosphere, providing key data for new progress to be made on the role of magnetic fields in heating the upper atmospheres, launching stellar winds, and more generally in the interaction of cool stars with their environment (circumstellar disk, planets) along their whole evolution. UV spectropolarimetry is proposed on missions of various sizes and scopes, from POLLUX on the 15-m telescope LUVOIR to the Arago M-size mission dedicated to UV spectropolarimetry.

The Hubble Space Telescope has produced astonishing science over the past thirty years. Hubble's productivity can continue to soar for years to come provided some worn out components get upgraded. While powerful new ground-based and space telescopes are expected to come online over the next decade, none of them will have the UV capabilities that make Hubble a unique observatory. Without Hubble, progress in UV and blue optical astrophysics will be halted. Observations at these wavelengths are key for a range of unresolved astrophysics questions, ranging from the characterization of solar system planets to understanding interaction of galaxies with the intergalactic medium and the formation history of the universe. Hubble will remain our only source of high-angular resolution UV imaging and high-sensitivity UV spectroscopy for the next two decades, offering the ability for continued unique science and maximizing the science return from complementary observatories. Therefore, we recommend that NASA, ESA, and the private sector study the scientific merit, technical feasibility, and risk of a new servicing mission to Hubble to boost its orbit, fix aging components, and expand its instrumentation. Doing so would: 1) keep Hubble on its path to reach its unmet full potential, 2) extend the mission's lifetime past the next decade, which will maximize the synergy of Hubble with other upcoming facilities, and 3) enable and enhance the continuation of scientific discoveries in UV and optical astrophysics.

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^\alpha$, with $\alpha$ 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.

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.

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.

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.

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.

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.

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\AA\ - 9200\AA\ 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.

POLLUX is a high-resolution, UV spectropolarimeter proposed for the 15-meter primary mirror option of LUVOIR. The instrument Phase 0 study is supported by the French Space Agency (CNES) and performed by a consortium of European scientists. POLLUX has been designed to deliver high-resolution spectroscopy (R> 120,000) over a broad spectral range (90-390 nm). Its unique spectropolarimetric capabilities will open-up a vast new parameter space, in particular in the unexplored UV domain and in a regime where high-resolution observations with current facilities in the visible domain are severely photon starved. POLLUX will address a range of questions at the core of the LUVOIR Science portfolio. The combination of high resolution and broad coverage of the UV bandpass will resolve narrow UV emission and absorption lines originating in diffuse media, thus permitting the study of the baryon cycle over cosmic time: from galaxies forming stars out of interstellar gas and grains, and stars forming planets, to the various forms of feedback into the interstellar and intergalactic medium (ISM and IGM), and active galactic nuclei (AGN). UV circular and linear polarimetry will reveal the magnetic fields for a wide variety of objects for the first time, from AGN outflows to a diverse range of stars, stellar explosions (both supernovae and their remnants), the ISM and IGM. It will enable detection of polarized light reflected from exoplanets (or their circumplanetary material and moons), characterization of the magnetospheres of stars and planets (and their interactions), and measurements of the influence of magnetic fields at the (inter)galactic scale. In this paper, we outline the key science cases of POLLUX, together with its high-level technical requirements. The instrument design, its estimated performances, and the required technology development are presented in a separated paper.

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.

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

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.

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

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 $\sigma$) 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$\mu$G) like in the reverse shock.)

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.

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 $\xi\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]

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.

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.

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

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

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.

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.

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.

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.

1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $\gamma$-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\"ahovi 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.

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 $\gamma$-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$\beta$ 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.

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

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.

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.

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

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.

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.

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.

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.

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.

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\"ahovi 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.

V745 Sco is a recurrent nova, with the most recent eruption occurring in February 2014. V745 Sco was first observed by Swift a mere 3.7 hr after the announcement of the optical discovery, with the super-soft X-ray emission being detected around four days later and lasting for only ~two days, making it both the fastest follow-up of a nova by Swift and the earliest switch-on of super-soft emission yet detected. Such an early switch-on time suggests a combination of a very high velocity outflow and low ejected mass and, together with the high effective temperature reached by the super-soft emission, a high mass white dwarf (>1.3 M_sun). The X-ray spectral evolution was followed from an early epoch where shocked emission was evident, through the entirety of the super-soft phase, showing evolving column density, emission lines, absorption edges and thermal continuum temperature. UV grism data were also obtained throughout the super-soft interval, with the spectra showing mainly emission lines from lower ionization transitions and the Balmer continuum in emission. V745 Sco is compared with both V2491 Cyg (another nova with a very short super-soft phase) and M31N 2008-12a (the most rapidly recurring nova yet discovered). The longer recurrence time compared to M31N 2008-12a could be due to a lower mass accretion rate, although inclination of the system may also play a part. Nova V745 Sco (2014) revealed the fastest evolving super-soft source phase yet discovered, providing a detailed and informative dataset for study.

The MAGIC stereoscopic system collected 69 hours of Crab Nebula data between October 2009 and April 2011. Analysis of this data sample using the latest improvements in the MAGIC stereoscopic software provided an unprecedented precision of spectral and night-by-night light curve determination at gamma rays. We derived a differential spectrum with a single instrument from 50 GeV up to almost 30 TeV with 5 bins per energy decade. At low energies, MAGIC results, combined with Fermi-LAT data, show a flat and broad Inverse Compton peak. The overall fit to the data between 1 GeV and 30 TeV is not well described by a log-parabola function. We find that a modified log-parabola function with an exponent of 2.5 instead of 2 provides a good description of the data ($\chi^2=35/26$). Using systematic uncertainties of red the MAGIC and Fermi-LAT measurements we determine the position of the Inverse Compton peak to be at (53 $\pm$ 3stat + 31syst -13syst) GeV, which is the most precise estimation up to date and is dominated by the systematic effects. There is no hint of the integral flux variability on daily scales at energies above 300 GeV when systematic uncertainties are included in the flux measurement. We consider three state- of-the-art theoretical models to describe the overall spectral energy distribution of the Crab Nebula. The constant B-field model cannot satisfactorily reproduce the VHE spectral measurements presented in this work, having particular difficulty reproducing the broadness of the observed IC peak. Most probably this implies that the assumption of the homogeneity of the magnetic field inside the nebula is incorrect. On the other hand, the time-dependent 1D spectral model provides a good fit of the new VHE results when considering a 80 {\mu}G magnetic field. However, it fails to match the data when including the morphology of the nebula at lower wavelengths.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

The high-frequency-peaked BL Lac (HBL) 1ES 0806+524 (z = 0.138) was discovered in VHE $\gamma$ rays in 2008. Until now, the broad-band spectrum of 1ES 0806+524 has been only poorly characterized, in particular at high energies. We analysed multiwavelength observations from $\gamma$ rays to radio performed from 2011 January to March, which were triggered by the high activity detected at optical frequencies. These observations constitute the most precise determination of the broad-band emission of 1ES 0806+524 to date. The stereoscopic MAGIC observations yielded a $\gamma$-ray signal above 250 GeV of $(3.7 \pm 0.7)$ per cent of the Crab Nebula flux with a statistical significance of 9.9 $\sigma$. The multiwavelength observations showed significant variability in essentially all energy bands, including a VHE $\gamma$-ray flare that lasted less than one night, which provided unprecedented evidence for short-term variability in 1ES 0806+524. The spectrum of this flare is well described by a power law with a photon index of $2.97 \pm 0.29$ between $\sim$150 GeV and 1 TeV and an integral flux of $(9.3 \pm 1.9)$ per cent of the Crab Nebula flux above 250 GeV. The spectrum during the non-flaring VHE activity is compatible with the only available VHE observation performed in 2008 with VERITAS when the source was in a low optical state. The broad-band spectral energy distribution can be described with a one-zone Synchrotron Self Compton model with parameters typical for HBLs, indicating that 1ES 0806+524 is not substantially different from the HBLs previously detected.

A distinguishing trait of the three known Galactic recurrent novae with the shortest orbital periods, T Pyx, IM Nor, and CI Aql, is that their optical decline time-scales are significantly longer than those of the other recurrent systems. On the other hand, some estimates of the mass of the ejecta, the velocity of the ejecta, and the duration of the soft X-rays emission of these systems are of the order of those of the other recurrent systems and the fast classical novae. We put forth a tentative explanation of this phenomenon. We propose that in these systems part of the material transferred from the companion during the first few days of the eruption remains within the Roche lobe of the white dwarf, preventing the radiation from ionizing the ejecta of the system and increasing the optical decline time-scale. We explain why this phenomenon is more likely in systems with a high mass transfer rate and a short orbital period. Finally, we present a schematic model that shows that the material transferred from the companion is sufficient to absorb the radiation from the white dwarf in these systems, ultimately supporting this scenario as quantitatively realistic.

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ~50GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220GeV is (0.66+/-0.03)% of Crab Nebula flux in 50 h of observations. The angular resolution, defined as the sigma of a 2-dimensional Gaussian distribution, at those energies is < 0.07 degree, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11-18% in flux normalization and +/-0.15 for the energy spectrum power-law slope.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

We present the results of an intensive multiwavelength campaign on nova LMC 2012. This nova evolved very rapidly in all observed wavelengths. The time to fall two magnitudes in the V band was only 2 days. In X-rays the super soft phase began 13$\pm$5 days after discovery and ended around day 50 after discovery. During the super soft phase, the \Swift/XRT and \Chandra\ spectra were consistent with the underlying white dwarf being very hot, $\sim$ 1 MK, and luminous, $\sim$ 10$^{38}$ erg s$^{-1}$. The UV, optical, and near-IR photometry showed a periodic variation after the initial and rapid fading had ended. Timing analysis revealed a consistent 19.24$\pm$0.03 hr period in all UV, optical, and near-IR bands with amplitudes of $\sim$ 0.3 magnitudes which we associate with the orbital period of the central binary. No periods were detected in the corresponding X-ray data sets. A moderately high inclination system, $i$ = 60$\pm$10$^{\arcdeg}$, was inferred from the early optical emission lines. The {\it HST}/STIS UV spectra were highly unusual with only the \ion{N}{5} (1240\AA) line present and superposed on a blue continuum. The lack of emission lines and the observed UV and optical continua from four epochs can be fit with a low mass ejection event, $\sim$ 10$^{-6}$ M$_{\odot}$, from a hot and massive white dwarf near the Chandrasekhar limit. The white dwarf, in turn, significantly illuminated its subgiant companion which provided the bulk of the observed UV/optical continuum emission at the later dates. The inferred extreme white dwarf characteristics and low mass ejection event favor nova LMC 2012 being a recurrent nova of the U Sco subclass.

Supermassive black holes with masses of millions to billions of solar masses are commonly found in the centers of galaxies. Astronomers seek to image jet formation using radio interferometry, but still suffer from insufficient angular resolution. An alternative method to resolve small structures is to measure the time variability of their emission. Here, we report on gamma-ray observations of the radio galaxy IC 310 obtained with the MAGIC telescopes revealing variability with doubling time scales faster than 4.8 min. Causality constrains the size of the emission region to be smaller than 20\% of the gravitational radius of its central black hole. We suggest that the emission is associated with pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base of the radio jet.

The number of known very high energy (VHE) blazars is $\sim\,50$, which is very small in comparison to the number of blazars detected in other frequencies. This situation is a handicap for population studies of blazars, which emit about half of their luminosity in the $\gamma$-ray domain. Moreover, VHE blazars, if distant, allow for the study of the environment that the high-energy $\gamma$-rays traverse in their path towards the Earth, like the extragalactic background light (EBL) and the intergalactic magnetic field (IGMF), and hence they have a special interest for the astrophysics community. We present the first VHE detection of 1ES\,0033+595 with a statistical significance of 5.5\,$\sigma$. The VHE emission of this object is constant throughout the MAGIC observations (2009 August and October), and can be parameterized with a power law with an integral flux above 150 GeV of $(7.1\pm1.3)\times 10^{-12} {\mathrm{ph\,cm^{-2}\,s^{-1}}}$ and a photon index of ($3.8\pm0.7$). We model its spectral energy distribution (SED) as the result of inverse Compton scattering of synchrotron photons. For the study of the SED we used simultaneous optical R-band data from the KVA telescope, archival X-ray data by \textit{Swift} as well as \textit{INTEGRAL}, and simultaneous high energy (HE, $300$\,MeV~--~$10$\,GeV) $\gamma$-ray data from the \textit{Fermi} LAT observatory. Using the empirical approach of Prandini et al. (2010) and the \textit{Fermi}-LAT and MAGIC spectra for this object, we estimate the redshift of this source to be $0.34\pm0.08\pm0.05$. This is a relevant result because this source is possibly one of the ten most distant VHE blazars known to date, and with further (simultaneous) observations could play an important role in blazar population studies, as well as future constraints on the EBL and IGMF.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

We aim to characterize the broadband emission from 2FGL J2001.1+4352, which has been associated with the unknown-redshift blazar MG4 J200112+4352. Based on its gamma-ray spectral properties, it was identified as a potential very high energy (VHE; E > 100 GeV) gamma-ray emitter. The source was observed with MAGIC first in 2009 and later in 2010 within a multi-instrument observation campaign. The MAGIC observations yielded 14.8 hours of good quality stereoscopic data. The object was monitored at radio, optical and gamma-ray energies during the years 2010 and 2011. The source, named MAGIC J2001+439, is detected for the first time at VHE with MAGIC at a statistical significance of 6.3 {\sigma} (E > 70 GeV) during a 1.3-hour long observation on 2010 July 16. The multi-instrument observations show variability in all energy bands with the highest amplitude of variability in the X-ray and VHE bands. We also organized deep imaging optical observations with the Nordic Optical Telescope in 2013 to determine the source redshift. We determine for the first time the redshift of this BL Lac object through the measurement of its host galaxy during low blazar activity. Using the observational evidence that the luminosities of BL Lac host galaxies are confined to a relatively narrow range, we obtain z = 0.18 +/- 0.04. Additionally, we use the Fermi-LAT and MAGIC gamma-ray spectra to provide an independent redshift estimation, z = 0.17 +/- 0.10. Using the former (more accurate) redshift value, we adequately describe the broadband emission with a one-zone SSC model for different activity states and interpret the few-day timescale variability as produced by changes in the high-energy component of the electron energy distribution.

We critically examine the recent claimed detection of Raman scattered O VI at around 6830\AA\ in the iron curtain stage spectra of the classical CO nova V339 Del. The observed line variations are compatible in profile and timing of emission line strength with an excited state transition of neutral carbon. Line formation in classical nova ejecta is physically very different from that in symbiotic binaries, in which the O VI emission line is formed within the wind of the companion red giant at low differential velocity. The ejecta velocity and density structure prevent the scattering from producing analogous features. There might , however, be a broadband spectropolarimetric signature of the Raman process and also Rayleigh scattering at some stage in the expansion. We show that the neutral carbon spectrum, hitherto under-exploited for novae, is especially useful as a probe of the structure of the ejecta during the early, optically thick stages of the expansion