Source author record

P. G. Prada Moroni

P. G. Prada Moroni appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph.EP nucl-th

Catalog footprint

What is connected

129works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.

preprint2021arXiv

Combined searches for dark matter in dwarf spheroidal galaxies observed with the MAGIC telescopes, including new data from Coma Berenices and Draco

Milky Way dwarf spheroidal galaxies (dSphs) are among the best candidates to search for signals of dark matter annihilation with Imaging Atmospheric Cherenkov Telescopes, given their high mass-to-light ratios and the fact that they are free of astrophysical gamma-ray emitting sources. Since 2011, MAGIC has performed a multi-year observation program in search for Weakly Interacting Massive Particles (WIMPs) in dSphs. Results on the observations of Segue 1 and Ursa Major II dSphs have already been published and include some of the most stringent upper limits (ULs) on the velocity-averaged cross-section $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP annihilation from observations of dSphs. In this work, we report on the analyses of 52.1 h of data of Draco dSph and 49.5 h of Coma Berenices dSph observed with the MAGIC telescopes in 2018 and in 2019 respectively. No hint of a signal has been detected from either of these targets and new constraints on the $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP candidates have been derived. In order to improve the sensitivity of the search and reduce the effect of the systematic uncertainties due to the $J$-factor estimates, we have combined the data of all dSphs observed with the MAGIC telescopes. Using 354.3 h of dSphs good quality data, 95 % CL ULs on $\langle σ_{\mathrm{ann}} v \rangle$ have been obtained for 9 annihilation channels. For most of the channels, these results reach values of the order of $10^{-24} $cm$^3$/s at ${\sim}1$ TeV and are the most stringent limits obtained with the MAGIC telescopes so far.

preprint2021arXiv

Gaia-ESO Survey: Role of magnetic activity and starspots on pre-main sequence lithium evolution

Pre-main sequence models with inflated radii are needed to simultaneously reproduce the colour-magnitude diagram and the lithium depletion pattern in young open clusters. We tested a new set of PMS models including radius inflation due to starspots or magnetic inhibition of convection, using five clusters observed by the Gaia-ESO Survey, spanning the age range ~10-100 Myr where such effects could be important. Gaia-ESO radial velocities were combined with Gaia EDR3 astrometry to obtain clean lists of high-probability members for the five clusters. A Bayesian maximum likelihood method was adopted to derive the best model parameters and the cluster reddening and age. Models were calculated for different values of the mixing length parameter ($α_{ML}=2.0$, 1.5 and 1.0), without spots or with effective spot coverage $β_{spot}=0.2$ and 0.4. To reproduce the CMD and the Li depletion pattern in Gamma Vel A and B and in 25 Ori we need both a reduced convection efficiency $α_{ML}=1.0$ and an effective spot coverage of ~20%. We obtained ages of 18 Myr and 21 Myr for Gamma Vel A and B, respectively, and 19 Myr for 25 Ori. However, a single isochrone is not sufficient to account for the Li dispersion, and an increasing level of spot coverage as mass decreases seems to be required. The older clusters (NGC2451B at 30 Myr, NGC2547 at 35 Myr, and NGC2516 at 138 Myr) are consistent with standard models ($α_{ML}=2.0$ and no spots) except at low masses: a 20% spot coverage seems to reproduce the sequence of M-type stars better and might explain the observed abundance spread. The quality of Gaia-ESO data combined with Gaia allows us to gain important insights on PMS evolution. Models including starspots can provide a consistent explanation of the cluster sequences and Li abundances of young clusters, although a range of starspot coverage is required to fully reproduce the data.

preprint2021arXiv

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet.

preprint2021arXiv

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical Hα parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.

preprint2020arXiv

A search for dark matter in Triangulum II with the MAGIC telescopes

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

preprint2020arXiv

An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $ν_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $γ$-ray telescope FACT during a high $γ$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $γ$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $γ$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $Γ$ index of the intrinsic spectrum in the VHE $γ$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

preprint2020arXiv

Bounds on Lorentz invariance violation from MAGIC observation of GRB 190114C

On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification.

preprint2020arXiv

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

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

preprint2020arXiv

MAGIC observations of the nearby short gamma-ray burst GRB 160821B

The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma rays are still poorly constrained, with only $\sim$20 events detected in the GeV band, and none in the TeV band. GRB~160821B is one of the nearest short GRBs known at $z=0.162$. Recent analyses of the multiwavelength observational data of its afterglow emission revealed an optical-infrared kilonova component, characteristic of heavy-element nucleosynthesis in a BNS merger. Aiming to better clarify the nature of short GRBs, this burst was automatically followed up with the MAGIC telescopes, starting from 24 seconds after the burst trigger. Evidence of a gamma-ray signal is found above $\sim$0.5 TeV at a significance of $\sim3\,σ$ during observations that lasted until 4 hours after the burst. Assuming that the observed excess events correspond to gamma-ray emission from GRB 160821B, in conjunction with data at other wavelengths, we investigate its origin in the framework of GRB afterglow models. The simplest interpretation with one-zone models of synchrotron-self-Compton emission from the external forward shock has difficulty accounting for the putative TeV flux. Alternative scenarios are discussed where the TeV emission can be relatively enhanced. The role of future GeV-TeV observations of short GRBs in advancing our understanding of BNS mergers and related topics is briefly addressed.

preprint2020arXiv

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

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

preprint2020arXiv

Monitoring of the radio galaxy M87 during a low emission state from 2012 to 2015 with MAGIC

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

preprint2020arXiv

New hard-TeV extreme blazars detected with the MAGIC telescopes

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

preprint2020arXiv

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

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

preprint2020arXiv

On the Metamorphosis of the Bailey diagram for RR Lyrae stars

We collected over 6000 high-resolution spectra of four dozen field RR Lyrae (RRL) variables pulsating either in the fundamental (39 RRab) or in the first overtone (9 RRc) mode. We measured radial velocities (RVs) of four strong metallic and four Balmer lines along the entire pulsational cycle and derived RV amplitudes with accuracies better than 1$-$2~\kmsec. The new amplitudes were combined with literature data for 23~RRab and 3~RRc stars (total sample 74 RRLs) which allowed us to investigate the variation of the Bailey diagram (photometric amplitude versus period) when moving from optical to mid-infrared bands and to re-cast the Bailey diagram in terms of RV amplitudes. We found that RV amplitudes for RRab are minimally affected by nonlinear phenomena (shocks) and multi-periodicity (Blazhko effect). The RV slope ($\log P$--A(V$_r$)) when compared with the visual slope ($\log P$--A($V$)) is shallower and the dispersion, at fixed period, decreases by a factor of two. We constructed homogeneous sets of Horizontal Branch evolutionary models and nonlinear, convective pulsation models of RRLs to constrain the impact of evolutionary effects on their pulsation properties. Evolution causes, on the Bailey diagram based on RV amplitudes, a modest variation in pulsation period and a large dispersion in amplitude. The broad dispersion in period of the Bailey diagram is mainly caused by variation in RRL intrinsic parameters (stellar mass, chemical composition). Empirical evidence indicates that RV amplitudes are an optimal diagnostic for tracing the mean effective temperature across the RRab instability strip.

preprint2020arXiv

Protostellar accretion and the cosmological lithium problem

The cosmological lithium problem, i.e. the discrepancy between the lithium abundance predicted by the Big Bang Nucleosynthesis and the one observed for the stars of the "Spite plateau", is one of the long standing problems of modern astrophysics. A possible astrophysical solution involves lithium burning due to protostellar mass accretion on Spite plateau stars. In present work, for the first time, we investigate with accurate evolutionary computations the impact of accretion on the lithium evolution in the metal-poor regime, that relevant for stars in the Spite plateau.

preprint2020arXiv

Protostellar accretion in low mass metal poor stars and the cosmological lithium problem

The cosmological lithium problem, that is, the discrepancy between the lithium abundance predicted by the Big Bang nucleosynthesis and the one observed for the stars of the 'Spite plateau', is one of the long standing problems of modern astrophysics. Recent hints for a possible solution involve lithium burning induced by protostellar mass accretion on Spite plateau stars. The purpose of this paper is to analyze the effect of protostellar accretion on low metallicity low-mass stars with a focus on PMS lithium evolution. We computed the evolution from the protostar to the MS phase of accreting models with final masses of 0.7 and 0.8 M$_\odot$, and three metallicities Z=0.0001, Z=0.0010, and Z=0.0050. The effects of changing the main parameters affecting accreting models, that is the accretion energy (cold versus hot accretion), the initial seed mass $M_{seed}$ and radius $R_{seed}$, and the mass accretion rate $\dot{m}$, have been investigated in detail. As for the main stellar properties and the surface $^7 Li$ abundance, hot accretion models converge to standard non-accreting ones within 1 Myr, regardless of the actual value of $M_{seed}$, $R_{seed}$, and $\dot{m}$. Also, cold accretion models with a relatively large $M_{seed}$ ($\gtrsim 10~M_{jup}$) or $R_{seed}$ ($\gtrsim 1~R_\odot$) converge to standard non-accreting ones in less than about 10-20~Myr. A drastically different evolution occurs whenever a cold protostellar accretion process starts from small values of $M_{seed}$ and $R_{seed}$ ($M_{seed}\sim 1~M_{jup}$, $R_{seed} \lesssim 1~R_\odot$). These models almost entirely skip the standard Hayashi track evolution and deplete Li before the end of the accretion phase. The exact amount of depletion depends on the actual combination of the accretion parameters ($\dot{m}$, $M_{seed}$, and $R_{seed}$), achieving in some cases the complete exhaustion of Li in the whole star.

preprint2020arXiv

Relevance of the small frequency separation for asteroseismic stellar age, mass, and radius. A statistical investigation for main-sequence low-mass stars

We performed a theoretical analysis on the relevance of the small frequency separation $δν$ in determining stellar ages, masses, and radii. We adopted the SCEPtER pipeline for low-mass stars, [0.7, 1.05] Msun. Synthetic stars were generated and reconstructed assuming different relative precision in $δν$ (5% and 2%). The quantification of the systematic errors arising from a mismatch between synthetic stars and the recovery grid was also performed. The biases were negligible. The statistical error in age estimates was strongly dependent on the stellar evolutionary phase. It is at its maximum at ZAMS and decreases to about 11% and 6% ($δν$ at 5% and 2% level) when stars reach the 30% of their evolutionary MS lifetime, then vanishes in the last 20% of the MS. For mass and radius estimates we detected nearly identical behavior for an observational uncertainty of 5%. No benefit was detected for mass and radius determinations from a reduction $δν$ precision to 2%. The impact on the age of the initial helium abundance resulted in negligible results. On the other hand greater bias (2% and 1%) in mass and radius estimates were detected whenever $δν$ is in the observational pool. This result, together with the presence of further unexplored uncertainty sources, suggest that precision in the derived stellar quantities below these thresholds may possibly be overoptimistic. [...] Finally, we compared the estimates by the SCEPtER pipeline for 13 Kepler asteroseismic LEGACY sample stars with those given by six different pipelines. Overall, on the LEGACY sample data, we obtained a multi-pipeline precision of about 4.4%, 1.7%, and 11% on the estimated masses, radii, and ages, respectively.

preprint2020arXiv

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

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

preprint2020arXiv

Studying the nature of the unidentified gamma-ray source HESS J1841-055 with the MAGIC telescopes

We investigate the physical nature and origin of the gamma-ray emission from the extended source HESS J1841-055 observed at TeV and GeV energies. We observed HESS J1841-055 at TeV energies for a total effective time of 43 hours with the MAGIC telescopes, in 2012 and 2013. Additionally, we analysed the GeV counterpart making use of about 10 years of Fermi-LAT data. Using both Fermi-LAT and MAGIC, we study both the spectral and energy-dependent morphology of the source for almost four decades of energy. The origin of the gamma-ray emission from this region is investigated using multi-waveband information on sources present in this region, suggested to be associated with this unidentified gamma-ray source. We find that the extended emission at GeV-TeV energies is best described by more than one source model. We also perform the first energy-dependent analysis of the HESS J1841-055 region at GeV-TeV. We find that the emission at lower energies comes from a diffuse or extended component, while the major contribution of gamma rays above 1 TeV arises from the southern part of the source. Moreover, we find that a significant curvature is present in the combined observed spectrum of MAGIC and Fermi-LAT. The first multi-wavelength spectral energy distribution of this unidentified source shows that the emission at GeV-TeV energies can be well explained with both leptonic and hadronic models. For the leptonic scenario, bremsstrahlung is the dominant emission compared to inverse Compton. On the other hand, for the hadronic model, gamma-ray resulting from the decay of neutral pions ($π^0$) can explain the observed spectrum. The presence of dense molecular clouds overlapping with HESS J1841-055 makes both bremsstrahlung and $π^0$-decay processes the dominant emission mechanisms for the source.

preprint2020arXiv

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

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

preprint2020arXiv

Testing two-component models on very-high-energy gamma-ray emitting BL Lac objects

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

preprint2020arXiv

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

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

preprint2020arXiv

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

preprint2019arXiv

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

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

preprint2019arXiv

Constraints on gamma-ray and neutrino emission from NGC 1068 with the MAGIC telescopes

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

preprint2018arXiv

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

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

preprint2018arXiv

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

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $γ$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

preprint2018arXiv

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

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

preprint2018arXiv

The extreme HBL behaviour of Markarian 501 during 2012

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.

preprint2017arXiv

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

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

preprint2017arXiv

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

Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the MAGIC collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent~CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of $E_{\mathrm{QG}_1} > 5.5\cdot 10^{17}$GeV ($4.5\cdot 10^{17}$GeV) for a linear, and $E_{\mathrm{QG}_2} > 5.9\cdot 10^{10}$GeV ($5.3\cdot 10^{10}$GeV) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, found to worsen the statistical limits by about 36--42\%. Our constraints would have resulted much more competitive if the intrinsic pulse shape of the pulsar between 200GeV and 400GeV was understood in sufficient detail and allowed inclusion of events well below 400GeV.

preprint2017arXiv

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

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.

preprint2017arXiv

MAGIC observations of the microquasar V404 Cygni during the 2015 outburst

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

preprint2017arXiv

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

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.

preprint2017arXiv

Performance of the MAGIC telescopes under moonlight

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

preprint2017arXiv

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

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black-hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV) regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behavior in a multiwavelength context, we compare our results with \textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI}, \textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index $Γ=3.2$. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95\% confidence level for energies above 200 GeV at $2.6\times10^{-12}$~photons cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$, respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of the MAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility.

preprint2017arXiv

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

In this work we present data from observations with the MAGIC telescopes of SN 2014J detected in January 21 2014, the closest Type Ia supernova since Imaging Air Cherenkov Telescopes started to operate. We probe the possibility of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early stages of Type Ia supernova explosions. We performed follow-up observations after this supernova explosion for 5 days, between January 27 and February 2 in 2014. We search for gamma-ray signal in the energy range between 100 GeV and several TeV from the location of SN 2014J using data from a total of $\sim5.5$ hours of observations. Prospects for observing gamma-rays of hadronic origin from SN 2014J in the near future are also being addressed. No significant excess was detected from the direction of SN 2014J. Upper limits at 95$\%$ confidence level on the integral flux, assuming a power-law spectrum, d$F/$d$E\propto E^{-Γ}$, with a spectral index of $Γ=2.6$, for energies higher than 300 GeV and 700 GeV, are established at $1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$, respectively. For the first time, upper limits on the VHE emission of a Type Ia supernova are established. The energy fraction isotropically emitted into TeV gamma rays during the first $\sim10$ days after the supernova explosion for energies greater than 300 GeV is limited to $10^{-6}$ of the total available energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario, the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be detectable in the future by any current or planned generation of Imaging Atmospheric Cherenkov Telescopes.

preprint2016arXiv

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

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (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.

preprint2016arXiv

A statistical test on the reliability of the non-coevality of stars in binary systems

We develop a statistical test on the expected difference in age estimates of two coeval stars in detached double-lined eclipsing binary systems that are only caused by observational uncertainties. We focus on stars in the mass range [0.8; 1.6] Msun, and on stars in the main-sequence phase. The ages were obtained by means of the maximum-likelihood SCEPtER technique. The observational constraints used in the recovery procedure are stellar mass, radius, effective temperature, and metallicity [Fe/H]. We defined the statistic W computed as the ratio of the absolute difference of estimated ages for the two stars over the age of the older one. We determined the critical values of this statistics above which coevality can be rejected. The median expected difference in the reconstructed age between the coeval stars of a binary system -- caused alone by the observational uncertainties -- shows a strong dependence on the evolutionary stage. This ranges from about 20% for an evolved primary star to about 75% for a near ZAMS primary. The median difference also shows an increase with the mass of the primary star from 20% for 0.8 Msun stars to about 50% for 1.6 Msun stars. The reliability of these results was checked by repeating the process with a grid of stellar models computed by a different evolutionary code. We show that the W test is much more sensible to age differences in the binary system components than the alternative approach of comparing the confidence interval of the age of the two stars. We also found that the distribution of W is, for almost all the examined cases, well approximated by beta distributions. The proposed method improves upon the techniques that are commonly adopted for judging the coevality of an observed system. It also provides a result founded on reliable statistics that simultaneously accounts for all the observational uncertainties.

preprint2016arXiv

Calibrating convective-core overshooting with eclipsing binary systems. The case of low-mass main-sequence stars

In a robust statistical way, we quantify the uncertainty that affects the calibration of the overshooting efficiency parameter $β$ that is owing to the uncertainty on the observational data in double-lined eclipsing binary systems. We also quantify the bias that is caused by the lack of constraints on the initial helium content and on the efficiencies of the superadiabatic convection and microscopic diffusion. We adopted a modified grid-based SCEPtER pipeline using as observational constraints the effective temperatures, [Fe/H], masses, and radii of the two stars. In a reference scenario of mild overshooting $β= 0.2$ for the synthetic data, we found both large statistical uncertainties and biases on the estimated $β$. For the first 80% of the MS evolution, $β$ is biased and practically unconstrained in the whole explored range [0.0; 0.4]. In the last 5% of the MS the bias vanishes and the $1 σ$ error is about 0.05. For synthetic data computed with $β= 0.0$, the estimated $β$ is biased by about 0.12 in the first 80% of the MS evolution, and by 0.05 afterwards. Assuming an uncertainty of $\pm 1$ in the helium-to-metal enrichment ratio $ΔY/ΔZ$, we found that in the terminal part of the MS evolution the error on the estimated $β$ values ranges from -0.05 to +0.10, while $β$ is basically unconstrained throughout the explored range at earlier evolutionary stages. A uniform variation of $\pm 0.24$ in the mixing-length parameter around the solar-calibrated value causes in last 5% of the MS an uncertainty from -0.09 to +0.15. A complete neglect of diffusion in the stellar evolution computations produces a $1 σ$ uncertainty of $\pm 0.08$ in the last 5% of the MS, while $β$ is practically unconstrained in the first 80% of the MS. Overall, the calibration appears poorly reliable.

preprint2016arXiv

Deep observation of the NGC 1275 region with MAGIC: search of diffuse gamma-ray emission from cosmic rays in the Perseus cluster

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse gamma-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E>100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6$^\circ$ from the centre, have been detected as point-like VHE gamma-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is well described by a power law with a photon index of $3.6\pm0.2_{stat}\pm0.2_{syst}$ between 90 GeV and 1.2 TeV. We do not detect any diffuse gamma-ray emission from the cluster and set stringent constraints on its CR population. In order to bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes $α$. For $α=2.2$, the CR-to-thermal pressure within the cluster virial radius is constrained to be below 1-2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to <20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, $B_0$, that depend on the CR distribution. For $α=2.2$, $B_0\gtrsim5-8 μ$G, which is below the 25 $μ$G inferred from Faraday rotation measurements, whereas, for $α\lesssim2.1$, the hadronic interpretation of the diffuse radio emission is in contrast with our gamma-ray flux upper limits independently of the magnetic field strength.

preprint2016arXiv

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

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.

preprint2016arXiv

Dust photophoretic transport around a T Tauri star: Implications for comets composition

There is a growing body of evidences for the presence of crystalline material in comets. These crystals are believed to have been annealed in the inner part of the proto-solar nebula, while comets should have been formed in the outer regions. Several transport processes have been proposed to reconcile these two facts; among them a migration driven by photophoresis. The primarily goal of this work is to assess whether disk irradiation by a Pre-Main Sequence star would influence the photophoretic transport. To do so, we have implemented an evolving 1+1D model of an accretion disk, including advanced numerical techniques, undergoing a time-dependent irradiation, consistent with the evolution of the proto-Sun along the Pre-Main Sequence. The photophoresis is described using a formalism introduced in several previous works. Adopting the opacity prescription used in these former studies, we find that the disk irradiation enhances the photophoretic transport: the assumption of a disk central hole of several astronomical units in radius is no longer strictly required, whereas the need for an ad hoc introduction of photoevaporation is reduced. However, we show that a residual trail of small particles could annihilate the photophoretic driven transport via their effect on the opacity. We have also confirmed that the thermal conductivity of transported aggregates is a crucial parameter which could limit or even suppress the photophoretic migration and generate several segregation effects.

preprint2016arXiv

Effect of planet ingestion on low-mass stars evolution: the case of 2MASS J08095427--4721419 star in the Gamma Velorum cluster

We analysed the effects of planet ingestion on the characteristics of a pre-MS star similar to the Gamma Velorum cluster member 2MASS J08095427--4721419 (#52). We discussed the effects of changing the age $t_0$ at which the accretion episode occurs, the mass of the ingested planet and its chemical composition. We showed that the mass of the ingested planet required to explain the current [Fe/H]^#52 increases by decreasing the age $t_0$ and/or by decreasing the Iron content of the accreted matter. We compared the predictions of a simplified accretion method -- where only the variation of the surface chemical composition is considered -- with that of a full accretion model that properly accounts for the modification of the stellar structure. We showed that the two approaches result in different convective envelope extension which can vary up to 10 percent. We discussed the impact of the planet ingestion on a stellar model in the colour-magnitude diagram, showing that a maximum shift of about 0.06 dex in the colour and 0.07 dex in magnitude are expected and that such variations persist even much later the accretion episode. We also analysed the systematic bias in the stellar mass and age inferred by using a grid of standard non accreting models to recover the characteristics of an accreting star. We found that standard non accreting models can safely be adopted for mass estimate, as the bias is <= 6 percent, while much more caution should be used for age estimate where the differences can reach about 60 percent.

preprint2016arXiv

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

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.

preprint2016arXiv

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

The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, $E > 100$ GeV) $γ$-ray band with a statistical significance of 5.9 $σ$. The integral flux above 150 GeV is estimated to be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous high energy (HE, 100 MeV $ < E < 100$ GeV) $γ$-ray observations from Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution (SED) of H1722+119 shows surprising behaviour in the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Accreting pre-main sequence models and abundance anomalies in globular clusters

We investigated the possibility of producing helium enhanced stars in globular clusters by accreting polluted matter during the pre-main sequence phase. We followed the evolution of two different classes of pre-main sequence accreting models, one which neglects and the other that takes into account the protostellar evolution. We analysed the dependence of the final central helium abundance, of the tracks position in the HR diagram and of the surface lithium abundance evolution on the age at which the accretion of polluted material begins and on the main physical parameters that govern the protostellar evolution. The later is the beginning of the late accretion and the lower are both the central helium and the surface lithium abundances at the end of the accretion phase and in ZAMS (Zero Age Main Sequence). In order to produce a relevant increase of the central helium content the accretion of polluted matter should start at ages lower than 1 Myr. The inclusion of the protostellar evolution has a strong impact on the ZAMS models too. The adoption of a very low seed mass (i.e. 0.001 M$_{\odot}$) results in models with the lowest central helium and surface lithium abundances. The higher is the accretion rate and the lower is the final helium content in the core and the residual surface lithium. In the worst case -- i.e. seed mass 0.001 M$_\odot$ and accretion rate $\ge 10^{-5}$ M$_\odot$ yr$^{-1}$ -- the central helium is not increased at all and the surface lithium is fully depleted in the first few million years.

preprint2015arXiv

Astrophysical implications of the proton-proton cross section updates

The p(p,e^+ ν_e)^2H reaction rate is an essential ingredient for theoretical computations of stellar models. In the past several values of the corresponding S-factor have been made available by different authors. Prompted by a recent evaluation of S(E), we analysed the effect of the adoption of different proton-proton reaction rates on stellar models, focusing, in particular, on the age of mid and old stellar clusters (1-12 Gyr) and on standard solar model predictions. By comparing different widely adopted p(p,e^+ ν_e)^2H reaction rates, we found a maximum difference in the temperature regimes typical of main sequence hydrogen-burning stars (5x10^6 - 3x10^7 K) of about 3%. Such a variation translates into a change of cluster age determination lower than 1%. A slightly larger effect is observed in the predicted solar neutrino fluxes with a maximum difference, in the worst case, of about 8%. Finally we also notice that the uncertainty evaluation of the present proton-proton rate is at the level of few \permil, thus the p(p,e^+ ν_e)^2H reaction rate does not constitute anymore a significant uncertainty source in stellar models.

preprint2015arXiv

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

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.

preprint2015arXiv

Grid-based estimates of stellar ages in binary systems. SCEPtER: Stellar CharactEristics Pisa Estimation gRid

We investigate the performance of grid-based techniques in estimating the age of stars in detached eclipsing binary systems. We evaluate the precision of the estimates due to the uncertainty in the observational constraints, and the systematic bias caused by the uncertainty in convective core overshooting, element diffusion, mixing-length value, and initial helium content. We adopted the SCEPtER grid, which includes stars with mass in the range [0.8; 1.6] Msun and evolutionary stages from the ZAMS to the central hydrogen depletion. Age estimates have been obtained by a generalisation of the technique described in our previous work. We showed that the typical 1 sigma random error in age estimates - due to the uncertainty on the observational constraints - is about +- 7%, which is nearly independent of the masses of the two stars. However, such an error strongly depends on the evolutionary phase and becomes larger and asymmetric for stars near the ZAMS where it ranges from about +90% to -25%. The systematic bias due to the including mild and strong convective core overshooting is about 50% and 120% of the error due to observational uncertainties. A variation of +- 1 in the helium-to-metal enrichment ratio accounts for about +- 150% of the random error. The neglect of microscopic diffusion accounts for a bias of about 60% of the random error. We also introduced a statistical test of the expected difference in the recovered age of two coeval stars in a binary system. We find that random fluctuations within the current observational uncertainties can lead genuine coeval binary components to appear to be non-coeval with a difference in age as high as 60%.

preprint2015arXiv

Kepler423b: a half-Jupiter mass planet transiting a very old solar-like star

We report the spectroscopic confirmation of the Kepler object of interest KOI-183.01 (Kepler-423b), a half-Jupiter mass planet transiting an old solar-like star every 2.7 days. Our analysis is the first to combine the full Kepler photometry (quarters 1-17) with high-precision radial velocity measurements taken with the FIES spectrograph at the Nordic Optical Telescope. We simultaneously modelled the photometric and spectroscopic data-sets using Bayesian approach coupled with Markov chain Monte Carlo sampling. We found that the Kepler pre-search data conditioned (PDC) light curve of KOI-183 exhibits quarter-to-quarter systematic variations of the transit depth, with a peak-to-peak amplitude of about 4.3 % and seasonal trends reoccurring every four quarters. We attributed these systematics to an incorrect assessment of the quarterly variation of the crowding metric. The host star KOI-183 is a G4 dwarf with $M_\star=0.85\pm0.04$ M$_\rm{Sun}$, $R_\star=0.95\pm0.04$ R$_\rm{Sun}$, $T_\mathrm{eff}=5560\pm80$ K, $[M/H]=-0.10\pm0.05$ dex, and with an age of $11\pm2$ Gyr. The planet KOI-183b has a mass of $M_\mathrm{p}=0.595\pm0.081$ M$_\mathrm{Jup}$ and a radius of $R_\mathrm{p}=1.192\pm0.052$ R$_\mathrm{Jup}$, yielding a planetary bulk density of $ρ_\mathrm{p}=0.459\pm0.083$ g/cm$^{3}$. The radius of KOI-183b is consistent with both theoretical models for irradiated coreless giant planets and expectations based on empirical laws. The inclination of the stellar spin axis suggests that the system is aligned along the line of sight. We detected a tentative secondary eclipse of the planet at a 2-$σ$ confidence level ($ΔF_{\mathrm{ec}}=14.2\pm6.6$ ppm) and found that the orbit might have a small non-zero eccentricity of $e=0.019^{+0.028}_{-0.014}$. With a Bond albedo of $A_\mathrm{B}=0.037\pm0.019$, KOI-183b is one of the gas-giant planets with the lowest albedo known so far.

preprint2015arXiv

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

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

preprint2015arXiv

On the age of Galactic bulge microlensed dwarf and subgiant stars

Recent results by Bensby and collaborators on the ages of microlensed stars in the Galactic bulge have challenged the picture of an exclusively old stellar population. However, these age estimates have not been independently confirmed. In this paper we verify these results by means of a grid-based method and quantify the systematic biases that might be induced by some assumptions adopted to compute stellar models. We explore the impact of increasing the initial helium abundance, neglecting the element microscopic diffusion, and changing the mixing-length calibration in theoretical stellar track computations. We adopt the SCEPtER pipeline with a novel stellar model grid for metallicities [Fe/H] from -2.00 to 0.55 dex, and masses in the range [0.60; 1.60] Msun from the ZAMS to the helium flash at the red giant branch tip. We show for the considered evolutionary phases that our technique provides unbiased age estimates. Our age results are in good agreement with Bensby and collaborators findings and show 16 stars younger than 5 Gyr and 28 younger than 9 Gyr over a sample of 58. The effect of a helium enhancement as large as Delta Y/Delta Z = 5 is quite modest, resulting in a mean age increase of metal rich stars of 0.6 Gyr. Even simultaneously adopting a high helium content and the upper values of age estimates, there is evidence of 4 stars younger than 5 Gyr and 15 younger than 9 Gyr. For stars younger than 5 Gyr, the use of stellar models computed by neglecting microscopic diffusion or by assuming a super-solar mixing-length value leads to a mean increase in the age estimates of about 0.4 Gyr and 0.5 Gyr respectively. Even considering the upper values for the age estimates, there are four stars estimated younger than 5 Gyr is in both cases. Thus, the assessment of a sizeable fraction of young stars among the microlensed sample in the Galactic bulge appears robust.

preprint2015arXiv

The 2009 multiwavelength campaign on Mrk 421: Variability and correlation studies

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.

preprint2015arXiv

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

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

preprint2015arXiv

The major upgrade of the MAGIC telescopes, Part II: A performance study using observations of the Crab Nebula

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.

preprint2015arXiv

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

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 ($δ$, $B$ and $R$), in addition to the parameters related to the particle population.

preprint2015arXiv

Very-high-energy γ-ray observations of novae and dwarf novae with the MAGIC telescopes

Context. In the last five years the Fermi Large Area Telescope (LAT) instrument detected GeV γ-ray emission from five novae. The GeV emission can be interpreted in terms of an inverse Compton process of electrons accelerated in a shock. In this case it is expected that protons in the same conditions can be accelerated to much higher energies. Consequently they may produce a second component in the γ-ray spectrum at TeV energies. Aims. We aim to explore the very-high-energy domain to search for γ-ray emission above 50 GeV and to shed light on the acceleration process of leptons and hadrons in nova explosions. Methods. We have performed observations with the MAGIC telescopes of the classical nova V339 Del shortly after the 2013 outburst, triggered by optical and subsequent GeV γ-ray detec- tions. We also briefly report on VHE observations of the symbiotic nova YY Her and the dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the analysis of con- temporaneous Fermi-LAT data of the sources. The TeV and GeV observations are compared in order to evaluate the acceleration parameters for leptons and hadrons. Results. No significant TeV emission was found from the studied sources. We computed upper limits on the spectrum and night-by-night flux. The combined GeV and TeV observations of V339 Del limit the ratio of proton to electron luminosities to Lp<~0.15 Le.

P. G. Prada Moroni

What is connected

Connect this record

See the researcher in context

Building this map preview

129 published item(s)

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

Combined searches for dark matter in dwarf spheroidal galaxies observed with the MAGIC telescopes, including new data from Coma Berenices and Draco

Gaia-ESO Survey: Role of magnetic activity and starspots on pre-main sequence lithium evolution

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

A search for dark matter in Triangulum II with the MAGIC telescopes

An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state

Bounds on Lorentz invariance violation from MAGIC observation of GRB 190114C

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

MAGIC observations of the nearby short gamma-ray burst GRB 160821B

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

Monitoring of the radio galaxy M87 during a low emission state from 2012 to 2015 with MAGIC

New hard-TeV extreme blazars detected with the MAGIC telescopes

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

On the Metamorphosis of the Bailey diagram for RR Lyrae stars

Protostellar accretion and the cosmological lithium problem

Protostellar accretion in low mass metal poor stars and the cosmological lithium problem

Relevance of the small frequency separation for asteroseismic stellar age, mass, and radius. A statistical investigation for main-sequence low-mass stars

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

Studying the nature of the unidentified gamma-ray source HESS J1841-055 with the MAGIC telescopes

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

Testing two-component models on very-high-energy gamma-ray emitting BL Lac objects

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

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

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

Constraints on gamma-ray and neutrino emission from NGC 1068 with the MAGIC telescopes

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

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

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

The extreme HBL behaviour of Markarian 501 during 2012

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

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

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

MAGIC observations of the microquasar V404 Cygni during the 2015 outburst

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

Performance of the MAGIC telescopes under moonlight

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

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

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

A statistical test on the reliability of the non-coevality of stars in binary systems

Calibrating convective-core overshooting with eclipsing binary systems. The case of low-mass main-sequence stars

Deep observation of the NGC 1275 region with MAGIC: search of diffuse gamma-ray emission from cosmic rays in the Perseus cluster

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

Dust photophoretic transport around a T Tauri star: Implications for comets composition

Effect of planet ingestion on low-mass stars evolution: the case of 2MASS J08095427--4721419 star in the Gamma Velorum cluster

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

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

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

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

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

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

Accreting pre-main sequence models and abundance anomalies in globular clusters

Astrophysical implications of the proton-proton cross section updates

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

Grid-based estimates of stellar ages in binary systems. SCEPtER: Stellar CharactEristics Pisa Estimation gRid

Kepler423b: a half-Jupiter mass planet transiting a very old solar-like star

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

On the age of Galactic bulge microlensed dwarf and subgiant stars

The 2009 multiwavelength campaign on Mrk 421: Variability and correlation studies

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

The major upgrade of the MAGIC telescopes, Part II: A performance study using observations of the Crab Nebula

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

Very-high-energy γ-ray observations of novae and dwarf novae with the MAGIC telescopes

Black hole lightning due to particle acceleration at subhorizon scales

Contemporaneous observations of the radio galaxy NGC 1275 from radio to very high energy gamma-rays

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

Discovery of TeV gamma-ray emission from the pulsar wind nebula 3C 58 by MAGIC

Discovery of very high energy gamma-ray emission from the blazar 1ES 0033+595 by the MAGIC telescopes

Evolution of the habitable zone of low-mass stars. Detailed stellar models and analytical relationships for different masses and chemical compositions

First broadband characterization and redshift determination of the VHE blazar MAGIC J2001+439

MAGIC gamma-ray and multifrequency observations of flat spectrum radio quasar PKS 1510-089 in early 2012

MAGIC long-term study of the distant TeV blazar PKS 1424+240 in a multiwavelength context

MAGIC reveals a complex morphology within the unidentified gamma-ray source HESS J1857+026

MAGIC search for VHE $γ$-ray emission from AE Aquarii in a multiwavelength context