Source author record

H. Kubo

H. Kubo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE hep-ex astro-ph.IM astro-ph.CO physics.ins-det nucl-ex astro-ph cond-mat.mtrl-sci

Catalog footprint

What is connected

121works

8topics

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

Constraint on the Matter-Antimatter Symmetry-Violating Phase in Neutrino Oscillations

The current laws of physics do not explain the observed imbalance of matter and antimatter in the universe. Sakharov proposed that an explanation would require the violation of CP symmetry between matter and antimatter. The only CP violation observed so far is in the weak interactions of quarks, and it is too small to explain the matter-antimatter imbalance of the universe. It has been shown that CP violation in the lepton sector could generate the matter-antimatter disparity through the process called leptogenesis. The quantum mixing of neutrinos, the neutral leptons in the Standard Model, provides a potential source of CP violation through a complex phase dCP, which may have consequences for theoretical models of leptogenesis. This CP violation can be measured in muon neutrino to electron neutrino oscillations and the corresponding antineutrino oscillations, which are experimentally accessible with accelerator-produced beams as established by the T2K experiment. Until now, the value of dCP has not been significantly constrained by neutrino oscillation experiments. Here the T2K collaboration reports a measurement that favors large enhancement of the neutrino oscillation probability, excluding values of dCP which result in a large enhancement of the observed anti-neutrino oscillation probability at three standard deviations (3 sigma). The 3 sigma confidence level interval for dCP, which is cyclic and repeats every 2pi, is [-3.41,-0.03] for the so-called normal mass ordering, and [-2.54,-0.32] for the inverted mass ordering. Our results show an indication of CP violation in the lepton sector. Herein we establish methods for sensitive searches for matter-antimatter asymmetry in neutrino oscillations using accelerator-produced neutrino beams. Future measurements with larger data samples will determine whether the leptonic CP violation is larger than the quark sector CP violation.

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.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]"

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

First combined measurement of the muon neutrino and antineutrino charged-current cross section without pions in the final state at T2K

This paper presents the first combined measurement of the double-differential muon neutrino and antineutrino charged-current cross sections with no pions in the final state on hydrocarbon at the off-axis near detector of the T2K experiment. The data analyzed in this work comprise 5.8$\times$10$^{20}$ and 6.3$\times$10$^{20}$ protons on target in neutrino and antineutrino mode respectively, at a beam energy peak of 0.6 GeV. Using the two measured cross sections, the sum, difference and asymmetry were calculated with the aim of better understanding the nuclear effects involved in such interactions. The extracted measurements have been compared with the prediction from different Monte Carlo generators and theoretical models showing that the difference between the two cross sections have interesting sensitivity to nuclear effects.

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

Search for heavy neutrinos with the T2K near detector ND280

This paper reports on the search for heavy neutrinos with masses in the range $140 < M_N < 493$ MeV/c$^2$ using the off-axis near detector ND280 of the T2K experiment. These particles can be produced from kaon decays in the standard neutrino beam and then subsequently decay in ND280. The decay modes under consideration are $N \to \ell^{\pm}_α π^{\mp}$ and $N \to \ell^+_α \ell^-_β ν(\barν)$ ($α,β=e,μ$). A search for such events has been made using the Time Projection Chambers of ND280, where the background has been reduced to less than two events in the current dataset in all channels. No excess has been observed in the signal region. A combined Bayesian statistical approach has been applied to extract upper limits on the mixing elements of heavy neutrinos to electron-, muon- and tau- flavoured currents ($U_e^2$, $U_μ^2$, $U_τ^2$) as a function of the heavy neutrino mass, e.g. $U_e^2 < 10^{-9}$ at $90\%$ C.L. for a mass of $390$ MeV/c$^2$. These constraints are competitive with previous experiments.

preprint2020arXiv

Simultaneous measurement of the muon neutrino charged-current cross section on oxygen and carbon without pions in the final state at T2K

This paper reports the first simultaneous measurement of the double differential muon neutrino charged-current cross section on oxygen and carbon without pions in the final state as a function of the outgoing muon kinematics, made at the ND280 off-axis near detector of the T2K experiment. The ratio of the oxygen and carbon cross sections is also provided to help validate various models' ability to extrapolate between carbon and oxygen nuclear targets, as is required in T2K oscillation analyses. The data are taken using a neutrino beam with an energy spectrum peaked at 0.6 GeV. The extracted measurement is compared with the prediction from different Monte Carlo neutrino-nucleus interaction event generators, showing particular model separation for very forward-going muons. Overall, of the models tested, the result is best described using Local Fermi Gas descriptions of the nuclear ground state with RPA suppression.

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

preprint2019arXiv

First Measurement of the Charged Current $\overlineν_μ$ Double Differential Cross Section on a Water Target without Pions in the final state

This paper reports the first differential measurement of the charged-current $\overlineν_μ$ interaction cross section on water with no pions in the final state. The unfolded flux-averaged measurement using the T2K off-axis near detector is given in double differential bins of $μ^+$ momentum and angle. The integrated cross section in a restricted phase space is $σ=\left(1.11\pm0.18\right)\times10^{-38}$ cm$^{2}$ per water molecule. Comparisons with several nuclear models are also presented.

preprint2019arXiv

Measurement of neutrino and antineutrino neutral-current quasielastic-like interactions on oxygen by detecting nuclear de-excitation $γ$-rays

Neutrino- and antineutrino-oxygen neutral-current quasielastic-like interactions are measured at Super-Kamiokande using nuclear de-excitation $γ$-rays to identify signal-like interactions in data from a $14.94 \ (16.35)\times 10^{20}$ protons-on-target exposure of the T2K neutrino (antineutrino) beam. The measured flux-averaged cross sections on oxygen nuclei are $\langle σ_{ν{\rm -NCQE}} \rangle = 1.70 \pm 0.17 ({\rm stat.}) ^{+ {\rm 0.51}}_{- {\rm 0.38}} ({\rm syst.}) \times 10^{-38} \ {\rm cm^2/oxygen}$ with a flux-averaged energy of 0.82 GeV and $\langle σ_{\barν {\rm -NCQE}} \rangle = 0.98 \pm 0.16 ({\rm stat.}) ^{+ {\rm 0.26}}_{- {\rm 0.19}} ({\rm syst.}) \times 10^{-38} \ {\rm cm^2/oxygen}$ with a flux-averaged energy of 0.68 GeV, for neutrinos and antineutrinos, respectively. These results are the most precise to date, and the antineutrino result is the first cross section measurement of this channel. They are compared with various theoretical predictions. The impact on evaluation of backgrounds to searches for supernova relic neutrinos at present and future water Cherenkov detectors is also discussed.

preprint2019arXiv

Measurement of the muon neutrino charged-current single $π^+$ production on hydrocarbon using the T2K off-axis near detector ND280

We report the measurements of single and double differential cross section of muon neutrino charged-current interactions on carbon with a single positively charged pion in the final state at the T2K off-axis near detector using $5.56\times10^{20}$ protons on target. The analysis uses data control samples for the background subtraction and the cross section signal, defined as a single negatively charged muon and a single positively charged pion exiting from the target nucleus, is extracted using an unfolding method. The model dependent cross section, integrated over the T2K off-axis neutrino beam spectrum peaking at $0.6$~GeV, is measured to be $σ= (11.76 \pm 0.44 \text{(stat)} \pm 2.39 \text{(syst)}) \times 10^{-40} \text{cm}^2$~$\text{nucleon}^{-1}$. Various differential cross sections are measured, including the first measurement of the Adler angles for single charged pion production in neutrino interactions with heavy nuclei target.

preprint2019arXiv

Search for Electron Antineutrino Appearance in a Long-baseline Muon Antineutrino Beam

Electron antineutrino appearance is measured by the T2K experiment in an accelerator-produced antineutrino beam, using additional neutrino beam operation to constrain parameters of the PMNS mixing matrix. T2K observes 15 candidate electron antineutrino events with a background expectation of 9.3 events. Including information from the kinematic distribution of observed events, the hypothesis of no electron antineutrino appearance is disfavored with a significance of 2.40σ and no discrepancy between data and PMNS predictions is found. A complementary analysis that introduces an additional free parameter which allows non-PMNS values of electron neutrino and antineutrino appearance also finds no discrepancy between data and PMNS predictions.

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

Prospects for CTA observations of the young SNR RX J1713.7-3946

We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We present a series of simulated images of RX J1713.7$-$3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the non-thermal X-ray emission observed by {\it XMM-Newton}, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H{\sc i} observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic vs leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.

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

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

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

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

Fermi LAT Discovery of Extended Gamma-Ray Emissions in the Vicinity of the HB3 Supernova Remnant

We report the discovery of extended gamma-ray emission measured by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope in the region of the supernova remnant (SNR) HB3 (G132.7+1.3) and the W3 HII complex adjacent to the southeast of the remnant. W3 is spatially associated with bright 12CO (J=1-0) emission. The gamma-ray emission is spatially correlated with this gas and the SNR. We discuss the possibility that gamma rays originate in interactions between particles accelerated in the SNR and interstellar gas or radiation fields. The decay of neutral pions produced in nucleon-nucleon interactions between accelerated hadrons and interstellar gas provides a reasonable explanation for the gamma-ray emission. The emission from W3 is consistent with irradiation of the CO clouds by the cosmic rays accelerated in HB3.

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

Measurements of $π^{\pm}$ differential yields from the surface of the T2K replica target for incoming 31 GeV/c protons with the NA61/SHINE spectrometer at the CERN SPS

Measurements of particle emission from a replica of the T2K 90 cm-long carbon target were performed in the NA61/SHINE experiment at CERN SPS, using data collected during a high-statistics run in 2009. An efficient use of the long-target measurements for neutrino flux predictions in T2K requires dedicated reconstruction and analysis techniques. Fully-corrected differential yields of $π^\pm$-mesons from the surface of the T2K replica target for incoming 31 GeV/c protons are presented. A possible strategy to implement these results into the T2K neutrino beam predictions is discussed and the propagation of the uncertainties of these results to the final neutrino flux is performed.

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

Readout technologies for directional WIMP Dark Matter detection

The measurement of the direction of WIMP-induced nuclear recoils is a compelling but technologically challenging strategy to provide an unambiguous signature of the detection of Galactic dark matter. Most directional detectors aim to reconstruct the dark-matter-induced nuclear recoil tracks, either in gas or solid targets. The main challenge with directional detection is the need for high spatial resolution over large volumes, which puts strong requirements on the readout technologies. In this paper we review the various detector readout technologies used by directional detectors. In particular, we summarize the challenges, advantages and drawbacks of each approach, and discuss future prospects for these technologies.

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

An Electron-Tracking Compton Telescope for a Survey of the Deep Universe by MeV gamma-rays

Photon imaging for MeV gammas has serious difficulties due to huge backgrounds and unclearness in images, which are originated from incompleteness in determining the physical parameters of Compton scattering in detection, e.g., lack of the directional information of the recoil electrons. The recent major mission/instrument in the MeV band, Compton Gamma Ray Observatory/COMPTEL, which was Compton Camera (CC), detected mere $\sim30$ persistent sources. It is in stark contrast with $\sim$2000 sources in the GeV band. Here we report the performance of an Electron-Tracking Compton Camera (ETCC), and prove that it has a good potential to break through this stagnation in MeV gamma-ray astronomy. The ETCC provides all the parameters of Compton-scattering by measuring 3-D recoil electron tracks; then the Scatter Plane Deviation (SPD) lost in CCs is recovered. The energy loss rate (dE/dx), which CCs cannot measure, is also obtained, and is found to be indeed helpful to reduce the background under conditions similar to space. Accordingly the significance in gamma detection is improved severalfold. On the other hand, SPD is essential to determine the point-spread function (PSF) quantitatively. The SPD resolution is improved close to the theoretical limit for multiple scattering of recoil electrons. With such a well-determined PSF, we demonstrate for the first time that it is possible to provide reliable sensitivity in Compton imaging without utilizing an optimization algorithm. As such, this study highlights the fundamental weak-points of CCs. In contrast we demonstrate the possibility of ETCC reaching the sensitivity below $1\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$ at 1 MeV.

preprint2015arXiv

CTA Contributions to the 34th International Cosmic Ray Conference (ICRC2015)

List of contributions from the CTA Consortium presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

preprint2015arXiv

Development and operational experience of magnetic horn system for T2K experiment

A magnetic horn system to be operated at a pulsed current of 320 kA and to survive high-power proton beam operation at 750 kW was developed for the T2K experiment. The first set of T2K magnetic horns was operated for over 12 million pulses during the four years of operation from 2010 to 2013, under a maximum beam power of 230 kW, and $6.63\times10^{20}$ protons were exposed to the production target. No significant damage was observed throughout this period. This successful operation of the T2K magnetic horns led to the discovery of the $ν_μ\rightarrowν_e$ oscillation phenomenon in 2013 by the T2K experiment. In this paper, details of the design, construction, and operation experience of the T2K magnetic horns are described.

preprint2015arXiv

Development of Slow Control Boards for the Large Size Telescopes of the Cherenkov Telescope Array

The camera of the Large Size Telescopes (LSTs) of the Cherenkov Telescope Array (CTA) consists of 265 photosensor modules, each of them containing 7 photomultiplier tubes (PMTs), a slow control board (SCB), a readout board, and a trigger logic. We have developed the SCB, which is installed between the 7 PMTs and the readout board. The main task for SCBs is the controlling of the high voltages for the PMTs and the monitoring of their anode currents. In addition, the SCB provides the functionality to create test pulses that can be injected at the input of the PMT preamplifier in order to emulate a PMT signal without the need of setting a high voltage, or even without the PMT itself. The test pulses have a very similar width as the PMT pulses (less than 3 ns FWHM) and their amplitude can be adjusted in a wide dynamic range. These features allow us not only to test the functionality of the camera modules but also to fully characterize these. We report on the design and the functions of the SCB together with the results of test measurements.

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

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

Measurement of the $ν_μ$ CCQE cross section on carbon with the ND280 detector at T2K

The Charged-Current Quasi-Elastic (CCQE) interaction, $ν_{l} + n \rightarrow l^{-} + p$, is the dominant CC process at $E_ν\sim 1$ GeV and contributes to the signal in accelerator-based long-baseline neutrino oscillation experiments operating at intermediate neutrino energies. This paper reports a measurement by the T2K experiment of the $ν_μ$ CCQE cross section on a carbon target with the off-axis detector based on the observed distribution of muon momentum ($p_μ$) and angle with respect to the incident neutrino beam ($θ_μ$). The flux-integrated CCQE cross section was measured to be $(0.83 \pm 0.12) \times 10^{-38}\textrm{ cm}^{2}$ in good agreement with NEUT MC value of ${0.88 \times 10^{-38}} \textrm{ cm}^{2}$. The energy dependence of the CCQE cross section is also reported. The axial mass, $M_A^{QE}$, of the dipole axial form factor was extracted assuming the Smith-Moniz CCQE model with a relativistic Fermi gas nuclear model. Using the absolute (shape-only) $p_μcosθ_μ$ distribution, the effective $M_A^{QE}$ parameter was measured to be ${1.26^{+0.21}_{-0.18} \textrm{ GeV}/c^{2}}$ (${1.43^{+0.28}_{-0.22} \textrm{ GeV}/c^{2}}$).

preprint2015arXiv

Measurement of the $ν_μ$ charged current quasi-elastic cross-section on carbon with the T2K on-axis neutrino beam

We report a measurement of the $ν_μ$ charged current quasi-elastic cross-sections on carbon in the T2K on-axis neutrino beam. The measured charged current quasi-elastic cross-sections on carbon at mean neutrino energies of 1.94 GeV and 0.93 GeV are $(11.95\pm 0.19(stat.)_{-1.47}^{+1.82} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$ and $(10.64\pm 0.37(stat.)_{-1.65}^{+2.03} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$, respectively. These results agree well with the predictions of neutrino interaction models. In addition, we investigated the effects of the nuclear model and the multi-nucleon interaction.

preprint2015arXiv

Measurement of the Electron Neutrino Charged-current Interaction Rate on Water with the T2K ND280 pi-zero Detector

This paper presents a measurement of the charged current interaction rate of the electron neutrino beam component of the beam above $1.5$~GeV using the large fiducial mass of the T2K $π^0$ detector. The predominant poriton of the $ν_e$ flux ($\sim$85 %) at these energies comes from kaon decays. The measured ratio of the observed beam interaction rate to the predicted rate in the detector with water targets filled is 0.89 $\pm$ 0.08 (stat.) $\pm$ 0.11 (sys.), and with the water targets emptied is 0.90 $\pm$ 0.09 (stat.) $\pm$ 0.13 (sys.). The ratio obtained for the interactions on water only from an event subtraction method is 0.87 $\pm$ 0.33 (stat.) $\pm$ 0.21 (sys.). This is the first measurement of the interaction rate of electron neutrinos on water, which is particularly of interest to experiments with water Cherenkov detectors.

preprint2015arXiv

Measurements of neutrino oscillation in appearance and disappearance channels by the T2K experiment with 6.6E20 protons on target

We report on measurements of neutrino oscillation using data from the T2K long-baseline neutrino experiment collected between 2010 and 2013. In an analysis of muon neutrino disappearance alone, we find the following estimates and 68% confidence intervals for the two possible mass hierarchies: Normal Hierarchy: $\sin^2θ_{23}=0.514^{+0.055}_{-0.056}$ and $Δm^2_{32}=(2.51\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ Inverted Hierarchy: $\sin^2θ_{23}=0.511\pm0.055$ and $Δm^2_{13}=(2.48\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ The analysis accounts for multi-nucleon mechanisms in neutrino interactions which were found to introduce negligible bias. We describe our first analyses that combine measurements of muon neutrino disappearance and electron neutrino appearance to estimate four oscillation parameters and the mass hierarchy. Frequentist and Bayesian intervals are presented for combinations of these parameters, with and without including recent reactor measurements. At 90% confidence level and including reactor measurements, we exclude the region: $δ_{CP}=[0.15,0.83]π$ for normal hierarchy and $δ_{CP}=[-0.08,1.09]π$ for inverted hierarchy. The T2K and reactor data weakly favor the normal hierarchy with a Bayes Factor of 2.2. The most probable values and 68% 1D credible intervals for the other oscillation parameters, when reactor data are included, are: $\sin^2θ_{23}=0.528^{+0.055}_{-0.038}$ and $|Δm^2_{32}|=(2.51\pm0.11)\times 10^{-3}$ eV$^2$/c$^4$.

preprint2015arXiv

Neutrino Oscillation Physics Potential of the T2K Experiment

The observation of the recent electron neutrino appearance in a muon neutrino beam and the high-precision measurement of the mixing angle $θ_{13}$ have led to a re-evaluation of the physics potential of the T2K long-baseline neutrino oscillation experiment. Sensitivities are explored for CP violation in neutrinos, non-maximal $\sin^22θ_{23}$, the octant of $θ_{23}$, and the mass hierarchy, in addition to the measurements of $δ_{CP}$, $\sin^2θ_{23}$, and $Δm^2_{32}$, for various combinations of $ν$-mode and $\barν$-mode data-taking. With an exposure of $7.8\times10^{21}$~protons-on-target, T2K can achieve 1-$σ$ resolution of 0.050(0.054) on $\sin^2θ_{23}$ and $0.040(0.045)\times10^{-3}~\rm{eV}^2$ on $Δm^2_{32}$ for 100\%(50\%) neutrino beam mode running assuming $\sin^2θ_{23}=0.5$ and $Δm^2_{32} = 2.4\times10^{-3}$ eV$^2$. T2K will have sensitivity to the CP-violating phase $δ_{\rm{CP}}$ at 90\% C.L. or better over a significant range. For example, if $\sin^22θ_{23}$ is maximal (i.e $θ_{23}$=$45^\circ$) the range is $-115^\circ<δ_{\rm{CP}}<-60^\circ$ for normal hierarchy and $+50^\circ<δ_{\rm{CP}}<+130^\circ$ for inverted hierarchy. When T2K data is combined with data from the NO$ν$A experiment, the region of oscillation parameter space where there is sensitivity to observe a non-zero $δ_{CP}$ is substantially increased compared to if each experiment is analyzed alone.

preprint2015arXiv

Performance of a new electron-tracking Compton camera under intense radiations from a water target irradiated with a proton beam

We have developed an electron-tracking Compton camera (ETCC) for use in next-generation MeV gamma ray telescopes. An ETCC consists of a gaseous time projection chamber (TPC) and pixel scintillator arrays (PSAs). Since the TPC measures the three dimensional tracks of Compton-recoil electrons, the ETCC can completely reconstruct the incident gamma rays. Moreover, the ETCC demonstrates efficient background rejection power in Compton-kinematics tests, identifies particle from the energy deposit rate (dE/dX) registered in the TPC, and provides high quality imaging by completely reconstructing the Compton scattering process. We are planning the "Sub-MeV gamma ray Imaging Loaded-on-balloon Experiment" (SMILE) for our proposed all-sky survey satellite. Performance tests of a mid-sized 30 cm-cubic ETCC, constructed for observing the Crab nebula, are ongoing. However, observations at balloon altitudes or satellite orbits are obstructed by radiation background from the atmosphere and the detector itself. The background rejection power was checked using proton accelerator experiments conducted at the Research Center for Nuclear Physics, Osaka University. To create the intense radiation fields encountered in space, which comprise gamma rays, neutrons, protons, and other energetic entities, we irradiated a water target with a 140 MeV proton beam and placed a SMILE-II ETCC near the target. In this situation, the counting rate was five times than that expected at the balloon altitude. Nonetheless, the ETCC stably operated and identified particles sufficiently to obtain a clear gamma ray image of the checking source. Here, we report the performance of our detector and demonstrate its effective background rejection based in electron tracking experiments.

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.

preprint2014arXiv

Black hole lightning due to particle acceleration at subhorizon scales

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

preprint2014arXiv

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

The Crab pulsar is the only astronomical pulsed source detected at very high energy (VHE, E>100GeV) gamma-rays. The emission mechanism of VHE pulsation is not yet fully understood, although several theoretical models have been proposed. In order to test the new models, we measured the light curve and the spectra of the Crab pulsar with high precision by means of deep observations. We analyzed 135 hours of selected MAGIC data taken between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower energies, 4.6 years of {\it Fermi}-LAT data were also analyzed. The known two pulses per period were detected with a significance of $8.0 σ$ and $12.6 σ$. In addition, significant emission was found between the two pulses with $6.2 σ$. We discovered the bridge emission above 50 GeV between the two main pulses. This emission can not be explained with the existing theories. These data can be used for testing new theoretical models.

preprint2014arXiv

Development of a 30 cm-cube Electron-Tracking Compton Camera for the SMILE-II Experiment

To explore the sub-MeV/MeV gamma-ray window for astronomy, we have developed the Electron-Tracking Compton Camera (ETCC), and carried out the first performance test at room condition using several gamma-ray sources in the sub-MeV energy band. Using a simple track analysis for a quick first test of the performance, the gamma-ray imaging capability was demonstrated by clear images and 5.3 degrees of angular resolution measure (ARM) measured at 662 keV. As the greatest impact of this work, a gamma-ray detection efficiency on the order of $10^{-4}$ was achieved at the sub-MeV gamma-ray band, which is one order of magnitude higher than our previous experiment. This angular resolution and detection efficiency enables us to detect the Crab Nebula at the 5 sigma level with several hours observation at balloon altitude in middle latitude. Furthermore, good consistency of efficiencies between this performance test and simulation including only physical processes has a large importance; it means we achieve nearly 100% detection of Compton recoil electrons. Thus, our estimation of enhancements by upgrades of the detector is more dependable. We are planning to confirm the imaging capability of the ETCC by observation of celestial objects in the SMILE-II (Sub-MeV gamma ray Imaging Loaded-on-balloon Experiment II). The SMILE-II and following SMILE-III project will be an important key of sub-MeV/MeV gamma-ray astronomy.

preprint2014arXiv

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

The pulsar wind nebula (PWN) 3C 58 is one of the historical very-high-energy (VHE; E>100 GeV) gamma-ray source candidates. It is energized by one of the highest spin-down power pulsars known (5% of Crab pulsar) and it has been compared to the Crab Nebula due to their morphological similarities. This object was previously observed by imaging atmospheric Cherenkov telescopes (Whipple, VERITAS and MAGIC), although not detected, with an upper limit of 2.4% Crab Unit (C.U.) at VHE. It was detected by Fermi-LAT with a spectrum extending beyond 100 GeV. We analyzed 81 hours of 3C 58 data taken with the MAGIC telescopes and we detected VHE gamma-ray emission with a significance of 5.7 sigma and an integral flux of 0.65% C.U. above 1 TeV. The differential energy spectrum between 400 GeV and 10 TeV is well described by a power-law function dϕ/dE=f_0(E/1TeV)^{-Gamma} with f_0=(2.0\pm0.4_{stat}\pm0.6_{sys})\times10^{-13}cm^{-2}s^{-1}TeV^{-1} and Gamma=2.4\pm0.2_{stat}\pm0.2_{sys}. The skymap is compatible with an unresolved source. We report the first significant detection of PWN 3C 58 at TeV energies. According to our results 3C 58 is the least luminous VHE gamma-ray PWN ever detected at VHE and the one with the lowest flux at VHE to date. We compare our results with the expectations of time-dependent models in which electrons up-scatter photon fields. The best representation favors a distance to the PWN of 2 kpc and Far Infrared (FIR) comparable to CMB photon fields. If we consider an unexpectedly high FIR density, the data can also be reproduced by models assuming a 3.2 kpc distance. A low magnetic field, far from equipartition, is required to explain the VHE data. Hadronic contribution from the hosting supernova remnant (SNR) requires unrealistic energy budget given the density of the medium, disfavoring cosmic ray acceleration in the SNR as origin of the VHE gamma-ray emission.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

preprint2014arXiv

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

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

preprint2014arXiv

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

HESS J1857+026 is an extended TeV gamma-ray source that was discovered by H.E.S.S. as part of its Galactic plane survey. Given its broadband spectral energy distribution and its spatial coincidence with the young energetic pulsar PSR J1856+0245, the source has been put forward as a pulsar wind nebula (PWN) candidate. MAGIC has performed follow-up observations aimed at mapping the source down to energies approaching 100 GeV in order to better understand its complex morphology. HESS J1857+026 was observed by MAGIC in 2010, yielding 29 hours of good quality stereoscopic data that allowed us to map the source region in two separate ranges of energy. We detected very-high-energy gamma-ray emission from HESS J1857+026 with a significance of $12 σ$ above $150$ GeV. The differential energy spectrum between $100$ GeV and $13$ TeV is well described by a power law function $dN/dE = N_0(E/1\textrm{TeV})^{-Γ}$ with $N_0 = (5.37 \pm0.44_{stat} \pm1.5_{sys}) \times 10^{-12} (\textrm{TeV}^{-1} \textrm{cm}^{-2}$ $\textrm{ s}^{-1})$ and $Γ= 2.16\pm0.07_{stat} \pm0.15_{sys}$, which bridges the gap between the GeV emission measured by Fermi-LAT and the multi-TeV emission measured by H.E.S.S.. In addition, we present a detailed analysis of the energy-dependent morphology of this region. We couple these results with archival multi-wavelength data and outline evidence in favor of a two-source scenario, whereby one source is associated with a PWN, while the other could be linked with a molecular cloud complex containing an HII region and a possible gas cavity.

preprint2014arXiv

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

It has been claimed that the nova-like cataclysmic variable (CV) AE Aquarii (AE Aqr) is a very-high-energy (VHE, $E>$100 GeV) source both on observational and theoretical grounds. We aim to search for VHE gamma-ray emission from AE Aqr during different states of the source at several wavelengths to confirm or rule out previous claims of detection of gamma-ray emission from this object. We report on observations of AE Aqr performed by MAGIC. The source was observed during 12 hours as part of a multiwavelength campaign carried out between May and June 2012 covering the optical, X-ray, and gamma-ray ranges. Besides MAGIC, the other facilities involved were the KVA, Skinakas, and Vidojevica telescopes in the optical and Swift in X-rays. We calculated integral upper limits coincident with different states of the source in the optical. We computed upper limits to the pulsed emission limiting the signal region to 30% of the phaseogram and we also searched for pulsed emission at different frequencies applying the Rayleigh test. AE Aqr was not detected at VHE energies during the multiwavelength campaign. We establish integral upper limits at the 95\% confidence level for the steady emission assuming the differential flux proportional to a power-law function dϕ/dE \propto E^{-Gamma}, with a Crab-like photon spectral index of Gamma=2.6. The upper limit above 200 GeV is 6.4\times10^{-12} cm^{-2}s^{-1} and above 1 TeV is 7.4\times10^{-13} cm^{-2}s^{-1}. We obtained an upper limit for the pulsed emission of 2.6\times10^{-12} cm^{-2}s^{-1} for energies above 200 GeV. Applying the Rayleigh test for pulsed emission at different frequencies we did not find any significant signal. Our results indicate that AE Aqr is not a VHE gamma-ray emitter at the level of emission previously claimed. We have established the most constraining upper limits for the VHE gamma-ray emission of AE Aqr.

preprint2014arXiv

Measurement of the inclusive $ν_μ$ charged current cross section on iron and hydrocarbon in the T2K on-axis neutrino beam

We report a measurement of the $ν_μ$ inclusive charged current cross sections on iron and hydrocarbon in the T2K on-axis neutrino beam. The measured inclusive charged current cross sections on iron and hydrocarbon averaged over the T2K on-axis flux with a mean neutrino energy of 1.51 GeV are $(1.444\pm0.002(stat.)_{-0.157}^{+0.189}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, and $(1.379\pm0.009(stat.)_{-0.147}^{+0.178}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, respectively, and their cross section ratio is $1.047\pm0.007(stat.)\pm0.035(syst.)$. These results agree well with the predictions of the neutrino interaction model, and thus we checked the correct treatment of the nuclear effect for iron and hydrocarbon targets in the model within the measurement precisions.

preprint2014arXiv

Measurement of the Inclusive Electron Neutrino Charged Current Cross Section on Carbon with the T2K Near Detector

The T2K off-axis near detector, ND280, is used to make the first differential cross-section measurements of electron neutrino charged current interactions at energies ~1 GeV as a function of electron momentum, electron scattering angle and four-momentum transfer of the interaction. The total flux-averaged $ν_e$ charged current cross-section on carbon is measured to be $1.11\pm0.09~(stat)\pm0.18~(syst)\times10^{-38} cm^2/nucleon$. The differential and total cross-section measurements agree with the predictions of two leading neutrino interaction generators, NEUT and GENIE. The NEUT prediction is $1.23\times10^{-38} cm^2/nucleon$ and the GENIE prediction is $1.08\times10^{-38} cm^2/nucleon$. The total $ν_e$ charged current cross-section result is also in agreement with data from the Gargamelle experiment.

preprint2014arXiv

Measurement of the intrinsic electron neutrino component in the T2K neutrino beam with the ND280 detector

The T2K experiment has reported the first observation of the appearance of electron neutrinos in a muon neutrino beam. The main and irreducible background to the appearance signal comes from the presence in the neutrino beam of a small intrinsic component of electron neutrinos originating from muon and kaon decays. In T2K, this component is expected to represent 1.2% of the total neutrino flux. A measurement of this component using the near detector (ND280), located 280 m from the target, is presented. The charged current interactions of electron neutrinos are selected by combining the particle identification capabilities of both the time projection chambers and electromagnetic calorimeters of ND280. The measured ratio between the observed electron neutrino beam component and the prediction is 1.01+-0.10 providing a direct confirmation of the neutrino fluxes and neutrino cross section modeling used for T2K neutrino oscillation analyses. Electron neutrinos coming from muons and kaons decay are also separately measured, resulting in a ratio with respect to the prediction of 0.68+-0.30 and 1.10+-0.14, respectively.

preprint2014arXiv

Measurement of the muon beam direction and muon flux for the T2K neutrino experiment

The Tokai-to-Kamioka (T2K) neutrino experiment measures neutrino oscillations by using an almost pure muon neutrino beam produced at the J-PARC accelerator facility. The T2K muon monitor was installed to measure the direction and stability of the muon beam which is produced together with the muon neutrino beam. The systematic error in the muon beam direction measurement was estimated, using data and MC simulation, to be 0.28 mrad. During beam operation, the proton beam has been controlled using measurements from the muon monitor and the direction of the neutrino beam has been tuned to within 0.3 mrad with respect to the designed beam-axis. In order to understand the muon beam properties,measurement of the absolute muon yield at the muon monitor was conducted with an emulsion detector. The number of muon tracks was measured to be $(4.06\pm0.05)\times10^4$ cm$^{-2}$ normalized with $4\times10^{11}$ protons on target with 250 kA horn operation. The result is in agreement with the prediction which is corrected based on hadron production data.

preprint2014arXiv

Measurement of the neutrino-oxygen neutral-current interaction cross section by observing nuclear deexcitation $γ$ rays

We report the first measurement of the neutrino-oxygen neutral-current quasielastic (NCQE) cross section. It is obtained by observing nuclear deexcitation $γ$-rays which follow neutrino-oxygen interactions at the Super-Kamiokande water Cherenkov detector. We use T2K data corresponding to $3.01 \times 10^{20}$ protons on target. By selecting only events during the T2K beam window and with well-reconstructed vertices in the fiducial volume, the large background rate from natural radioactivity is dramatically reduced. We observe 43 events in the $4-30$ MeV reconstructed energy window, compared with an expectation of 51.0, which includes an estimated 16.2 background events. The background is primarily nonquasielastic neutral-current interactions and has only 1.2 events from natural radioactivity. The flux-averaged NCQE cross section we measure is $1.55 \times 10^{-38}$ cm$^2$ with a 68\% confidence interval of $(1.22, 2.20) \times 10^{-38}$ cm$^2$ at a median neutrino energy of 630 MeV, compared with the theoretical prediction of $2.01 \times 10^{-38}$ cm$^2$.

preprint2014arXiv

Multifrequency Studies of the Peculiar Quasar 4C +21.35 During the 2010 Flaring Activity

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

preprint2014arXiv

Optimized dark matter searches in deep observations of Segue 1 with MAGIC

We present the results of stereoscopic observations of the satellite galaxy Segue 1 with the MAGIC Telescopes, carried out between 2011 and 2013. With almost 160 hours of good-quality data, this is the deepest observational campaign on any dwarf galaxy performed so far in the very high energy range of the electromagnetic spectrum. We search this large data sample for signals of dark matter particles in the mass range between 100 GeV and 20 TeV. For this we use the full likelihood analysis method, which provides optimal sensitivity to characteristic gamma-ray spectral features, like those expected from dark matter annihilation or decay. In particular, we focus our search on gamma-rays produced from different final state Standard Model particles, annihilation with internal bremsstrahlung, monochromatic lines and box-shaped signals. Our results represent the most stringent constraints to the annihilation cross-section or decay lifetime obtained from observations of satellite galaxies, for masses above few hundred GeV. In particular, our strongest limit (95% confidence level) corresponds to a ~500 GeV dark matter particle annihilating into tau+tau-, and is of order <sigma v> ~ 1.2x10^{-24} cm^3 s^{-1} - a factor ~40 above the <sigma v> thermal value.

preprint2014arXiv

Precise Measurement of the Neutrino Mixing Parameter θ_{23} from Muon Neutrino Disappearance in an Off-axis Beam

New data from the T2K neutrino oscillation experiment produce the most precise measurement of the neutrino mixing parameter theta_{23}. Using an off-axis neutrino beam with a peak energy of 0.6 GeV and a data set corresponding to 6.57 x 10^{20} protons on target, T2K has fit the energy-dependent nu_mu oscillation probability to determine oscillation parameters. Marginalizing over the values of other oscillation parameters yields sin^2 (theta_{23}) = 0.514 +0.055/-0.056 (0.511 +- 0.055), assuming normal (inverted) mass hierarchy. The best-fit mass-squared splitting for normal hierarchy is Delta m^2_{32} = (2.51 +- 0.10) x 10^{-3} eV^2/c^4 (inverted hierarchy: Delta m^2_{13} = (2.48 +- 0.10) x 10^{-3} eV^2/c^4). Adding a model of multinucleon interactions that affect neutrino energy reconstruction is found to produce only small biases in neutrino oscillation parameter extraction at current levels of statistical uncertainty.

preprint2014arXiv

Search for Proton Decay via $p \rightarrow νK^{+}$ using 260 kiloton$\cdot$year data of Super-Kamiokande

We have searched for proton decay via $p \rightarrow νK^{+}$ using Super-Kamiokande data from April 1996 to February 2013, 260 kiloton$\cdot$year exposure in total. No evidence for this proton decay mode is found. A lower limit of the proton lifetime is set to $5.9 \times 10^{33}$ years at 90% confidence level.

preprint2014arXiv

Search for short baseline $ν_e$ disappearance with the T2K near detector

The T2K experiment has performed a search for $ν_e$ disappearance due to sterile neutrinos using $5.9 \times 10^{20}$ protons on target for a baseline of $280 m$ in a neutrino beam peaked at about $500 MeV$. A sample of ν_e CC interactions in the off-axis near detector has been selected with a purity of 63\% and an efficiency of 26\%. The p-value for the null hypothesis is 0.085 and the excluded region at 95\% CL is approximately $sin^2 2 θ_{ee} > 0.3$ for $Δm^2_{eff} > 7 eV^2 / c^4$.

preprint2014arXiv

Search for Very-High-Energy Gamma Rays from the z = 0.896 Quasar 4C +55.17 with the MAGIC telescopes

The bright gamma-ray quasar 4C +55.17 is a distant source ($z = 0.896$) with a hard spectrum at GeV energies as observed by the Large Area Telescope (LAT) on board the {\it Fermi} satellite. This source is identified as a good source candidate for very-high-energy (VHE; $> 30$ GeV) gamma rays. In general VHE gamma rays from distant sources provide an unique opportunity to study the extragalactic background light (EBL) and underlying astrophysics. The flux intensity of this source in the VHE range is investigated. Then, constraints on the EBL are derived from the attenuation of gamma-ray photons coming from the distant blazar. We searched for a gamma-ray signal from this object using the 35-hour observations taken by the MAGIC telescopes between November 2010 and January 2011. No significant VHE gamma-ray signal was detected. We computed the upper limits of the integrated gamma-ray flux at $95\%$ confidence level of $9.4 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ and $2.5 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ above $100$ GeV and $200$ GeV, respectively. The differential upper limits in four energy bins in the range from $80$ GeV to $500$ GeV are also derived. The upper limits are consistent with the attenuation predicted by low-flux EBL models on the assumption of a simple power-law spectrum extrapolated from LAT data.

preprint2013arXiv

Calibration of the Super-Kamiokande Detector

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

preprint2013arXiv

CTA contributions to the 33rd International Cosmic Ray Conference (ICRC2013)

Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil

preprint2013arXiv

Measurement of Neutrino Oscillation Parameters from Muon Neutrino Disappearance with an Off-axis Beam

The T2K collaboration reports a precision measurement of muon neutrino disappearance with an off-axis neutrino beam with a peak energy of 0.6 GeV. Near detector measurements are used to constrain the neutrino flux and cross section parameters. The Super-Kamiokande far detector, which is 295 km downstream of the neutrino production target, collected data corresponding to $3.01 \times 10^{20}$ protons on target. In the absence of neutrino oscillations, $205 \pm 17$ (syst.) events are expected to be detected and only 58 muon neutrino event candidates are observed. A fit to the neutrino rate and energy spectrum assuming three neutrino flavors, normal mass hierarchy and $θ_{23}\leq π/4$ yields a best-fit mixing angle $\sin^2(2θ_{23})=1.000$ and mass splitting $|Δm^2_{32}| =2.44 \times 10^{-3}$ eV$^2$/c$^4$. If $θ_{23}\geq π/4$ is assumed, the best-fit mixing angle changes to $\sin^2(2θ_{23})=0.999$ and the mass splitting remains unchanged.

preprint2013arXiv

Measurement of the Inclusive NuMu Charged Current Cross Section on Carbon in the Near Detector of the T2K Experiment

T2K has performed the first measurement of νμ inclusive charged current interactions on carbon at neutrino energies of ~1 GeV where the measurement is reported as a flux-averaged double differential cross section in muon momentum and angle. The flux is predicted by the beam Monte Carlo and external data, including the results from the NA61/SHINE experiment. The data used for this measurement were taken in 2010 and 2011, with a total of 10.8 x 10^{19} protons-on-target. The analysis is performed on 4485 inclusive charged current interaction candidates selected in the most upstream fine-grained scintillator detector of the near detector. The flux-averaged total cross section is <σ_CC>_ϕ=(6.91 +/- 0.13 (stat) +/- 0.84 (syst)) x10^{-39} cm^2/nucleon for a mean neutrino energy of 0.85 GeV.

preprint2013arXiv

The T2K Neutrino Flux Prediction

The Tokai-to-Kamioka (T2K) experiment studies neutrino oscillations using an off-axis muon neutrino beam with a peak energy of about 0.6 GeV that originates at the J-PARC accelerator facility. Interactions of the neutrinos are observed at near detectors placed at 280 m from the production target and at the far detector -- Super-Kamiokande (SK) -- located 295 km away. The flux prediction is an essential part of the successful prediction of neutrino interaction rates at the T2K detectors and is an important input to T2K neutrino oscillation and cross section measurements. A FLUKA and GEANT3 based simulation models the physical processes involved in the neutrino production, from the interaction of primary beam protons in the T2K target, to the decay of hadrons and muons that produce neutrinos. The simulation uses proton beam monitor measurements as inputs. The modeling of hadronic interactions is re-weighted using thin target hadron production data, including recent charged pion and kaon measurements from the NA61/SHINE experiment. For the first T2K analyses the uncertainties on the flux prediction are evaluated to be below 15% near the flux peak. The uncertainty on the ratio of the flux predictions at the far and near detectors is less than 2% near the flux peak.

preprint2012arXiv

Dual baseline search for muon antineutrino disappearance at 0.1 eV^2 < Δm^2 < 100 eV^2

The MiniBooNE and SciBooNE collaborations report the results of a joint search for short baseline disappearance of \bar{ν_μ} at Fermilab's Booster Neutrino Beamline. The MiniBooNE Cherenkov detector and the SciBooNE tracking detector observe antineutrinos from the same beam, therefore the combined analysis of their datasets serves to partially constrain some of the flux and cross section uncertainties. Uncertainties in the ν_μ background were constrained by neutrino flux and cross section measurements performed in both detectors. A likelihood ratio method was used to set a 90% confidence level upper limit on \bar{ν_μ} disappearance that dramatically improves upon prior limits in the Δm^2=0.1-100 eV^2 region.

preprint2012arXiv

First Muon-Neutrino Disappearance Study with an Off-Axis Beam

We report a measurement of muon-neutrino disappearance in the T2K experiment. The 295-km muon-neutrino beam from Tokai to Kamioka is the first implementation of the off-axis technique in a long-baseline neutrino oscillation experiment. With data corresponding to 1.43 10**20 protons on target, we observe 31 fully-contained single muon-like ring events in Super-Kamiokande, compared with an expectation of 104 +- 14 (syst) events without neutrino oscillations. The best-fit point for two-flavor nu_mu -> nu_tau oscillations is sin**2(2 theta_23) = 0.98 and |Δm**2_32| = 2.65 10**-3 eV**2. The boundary of the 90 % confidence region includes the points (sin**2(2 theta_23),|Δm**2_32|) = (1.0, 3.1 10**-3 eV**2), (0.84, 2.65 10**-3 eV**2) and (1.0, 2.2 10**-3 eV**2).

preprint2012arXiv

Neutron imaging detector based on the muPIC micro-pixel chamber

We have developed a prototype time-resolved neutron imaging detector employing the micro-pixel chamber (muPIC), a micro-pattern gaseous detector, coupled with a field programmable gate array-based data acquisition system for applications in neutron radiography at high-intensity neutron sources. The prototype system, with an active area of 10cm x 10cm and operated at a gas pressure of 2 atm, measures both the energy deposition (via time-over-threshold) and 3-dimensional track of each neutron-induced event, allowing the reconstruction of the neutron interaction point with improved accuracy. Using a simple position reconstruction algorithm, a spatial resolution of 349 +/- 36 microns was achieved, with further improvement expected. The detailed tracking allows strong rejection of background gamma-rays, resulting in an effective gamma sensitivity of 10^-12 or less, coupled with stable, robust neutron identification. The detector also features a time resolution of 0.6 microseconds.

preprint2012arXiv

Pion emission from the T2K replica target: method, results and application

The T2K long-baseline neutrino oscillation experiment in Japan needs precise predictions of the initial neutrino flux. The highest precision can be reached based on detailed measurements of hadron emission from the same target as used by T2K exposed to a proton beam of the same kinetic energy of 30 GeV. The corresponding data were recorded in 2007-2010 by the NA61/SHINE experiment at the CERN SPS using a replica of the T2K graphite target. In this paper details of the experiment, data taking, data analysis method and results from the 2007 pilot run are presented. Furthermore, the application of the NA61/SHINE measurements to the predictions of the T2K initial neutrino flux is described and discussed.

preprint2012arXiv

Search for Proton Decay into Muon plus Neutral Kaon in Super-Kamiokande I, II, and III

We have searched for proton into muon plus neutral kaon using data from a 91.7 kiloton-year exposure of Super-Kamiokande-I, a 49.2 kiloton-year exposure of Super-Kamiokande-II, and a 31.9 kiloton-year exposure of Super-Kamiokande-III. The number of candidate events in the data was consistent with the atmospheric neutrino background expectation and no evidence for proton decay in this mode was found. We set a partial lifetime lower limit of 1.6x10^33 years at the 90% confidence level.

preprint2011arXiv

CANGAROO-III observation of TeV gamma rays from the unidentified gamma-ray source HESS J1614-518

We report the detection, with the CANGAROO-III imaging atmospheric Cherenkov telescope array, of a very high energy gamma-ray signal from the unidentified gamma-ray source HESS J1614-518, which was discovered in the H.E.S.S. Galactic plane survey. Diffuse gamma-ray emission was detected above 760 GeV at the 8.9 sigma level during an effective exposure of 54 hr from 2008 May to August. The spectrum can be represented by a power-law: 8.2+-2.2_{stat}+-2.5_{sys}x10^{-12}x (E/1TeV)^{-Gamma} cm^{-2} s^{-1} TeV^{-1} with a photon index Gamma of 2.4+-0.3_{stat}+-0.2_{sys}, which is compatible with that of the H.E.S.S. observations. By combining our result with multi-wavelength data, we discuss the possible counterparts for HESS J1614-518 and consider radiation mechanisms based on hadronic and leptonic processes for a supernova remnant, stellar winds from massive stars, and a pulsar wind nebula. Although a leptonic origin from a pulsar wind nebula driven by an unknown pulsar remains possible, hadronic-origin emission from an unknown supernova remnant is preferred.

preprint2011arXiv

Dual baseline search for muon neutrino disappearance at 0.5 eV^2 < Δm^2 < 40 eV^2

The SciBooNE and MiniBooNE collaborations report the results of a ν_μdisappearance search in the Δm^2 region of 0.5-40 eV^2. The neutrino rate as measured by the SciBooNE tracking detectors is used to constrain the rate at the MiniBooNE Cherenkov detector in the first joint analysis of data from both collaborations. Two separate analyses of the combined data samples set 90% confidence level (CL) limits on ν_μdisappearance in the 0.5-40 eV^2 Δm^2 region, with an improvement over previous experimental constraints between 10 and 30 eV^2.

preprint2011arXiv

Indication of Electron Neutrino Appearance from an Accelerator-produced Off-axis Muon Neutrino Beam

The T2K experiment observes indications of $ν_μ\rightarrow ν_e$ appearance in data accumulated with $1.43\times10^{20}$ protons on target. Six events pass all selection criteria at the far detector. In a three-flavor neutrino oscillation scenario with $|Δm_{23}^2|=2.4\times10^{-3}$ eV$^2$, $\sin^2 2θ_{23}=1$ and $\sin^2 2θ_{13}=0$, the expected number of such events is 1.5$\pm$0.3(syst.). Under this hypothesis, the probability to observe six or more candidate events is 7$\times10^{-3}$, equivalent to 2.5$σ$ significance. At 90% C.L., the data are consistent with 0.03(0.04)$<\sin^2 2θ_{13}<$ 0.28(0.34) for $δ_{\rm CP}=0$ and a normal (inverted) hierarchy.

preprint2011arXiv

Measurement of inclusive charged current interactions on carbon in a few-GeV neutrino beam

The SciBooNE Collaboration reports a measurement of inclusive charged current interactions of muon neutrinos on carbon with an average energy of 0.8 GeV using the Fermilab Booster Neutrino Beam. We compare our measurement with two neutrino interaction simulations: NEUT and NUANCE. The charged current interaction rates (product of flux and cross section) are extracted by fitting the muon kinematics, with a precision of 6-15% for the energy dependent and 3% for the energy integrated analyses. We also extract CC inclusive interaction cross sections from the observed rates, with a precision of 10-30% for the energy dependent and 8% for the energy integrated analyses. This is the first measurement of the CC inclusive cross section on carbon around 1 GeV. These results can be used to convert previous SciBooNE cross section ratio measurements to absolute cross section values.

preprint2011arXiv

Measurement of K+ production cross section by 8 GeV protons using high energy neutrino interactions in the SciBooNE detector

The SciBooNE Collaboration reports K+ production cross section and rate measurements using high energy daughter muon neutrino scattering data off the SciBar polystyrene (C8H8) target in the SciBooNE detector. The K+ mesons are produced by 8 GeV protons striking a beryllium target in Fermilab Booster Neutrino Beam line (BNB). Using observed neutrino and antineutrino events in SciBooNE, we measure d2σ/dpdΩ = (5.34 \times 0.76) mb/(GeV/c \times sr) for p + Be -> K+ + X at mean K+ energy of 3.9 GeV and angle (with respect to the proton beam direction) of 3.7 degrees, corresponding to the selected K+ sample. Compared to Monte Carlo predictions using previous higher energy K+ production measurements, this measurement, which uses the NUANCE neutrino interaction generator, is consistent with a normalization factor of 0.85\times0.12. This agreement is evidence that the extrapolation of the higher energy K+ measurements to an 8 GeV beam energy using Feynman scaling is valid. This measurement reduces the error on the K+ production cross section from 40% to 14%.

preprint2011arXiv

Measurements of the T2K neutrino beam properties using the INGRID on-axis near detector

Precise measurement of neutrino beam direction and intensity was achieved based on a new concept with modularized neutrino detectors. INGRID (Interactive Neutrino GRID) is an on-axis near detector for the T2K long baseline neutrino oscillation experiment. INGRID consists of 16 identical modules arranged in horizontal and vertical arrays around the beam center. The module has a sandwich structure of iron target plates and scintillator trackers. INGRID directly monitors the muon neutrino beam profile center and intensity using the number of observed neutrino events in each module. The neutrino beam direction is measured with accuracy better than 0.4 mrad from the measured profile center. The normalized event rate is measured with 4% precision.

preprint2010arXiv

Design and performance of the muon monitor for the T2K neutrino oscillation experiment

This article describes the design and performance of the muon monitor for the T2K (Tokaito-Kamioka) long baseline neutrino oscillation experiment. The muon monitor consists of two types of detector arrays: ionization chambers and silicon PIN photodiodes. It measures the intensity and profile of muons produced, along with neutrinos, in the decay of pions. The measurement is sensitive to the intensity and direction of the neutrino beam. The linearity and stability of the detectors were measured in beam tests to be within 2.4% and 1.5%, respectively. Based on the test results, the precision of the beam direction measured by the muon monitor is expected to be 0.25 mrad.

preprint2010arXiv

Measurement of Inclusive Neutral Current Neutral Pion Production on Carbon in a Few-GeV Neutrino Beam

The SciBooNE Collaboration reports inclusive neutral current neutral pion production by a muon neutrino beam on a polystyrene target (C8H8). We obtain (7.7 \pm 0.5(stat.) \pm 0.5 (sys.)) x 10^(-2) as the ratio of the neutral current neutral pion production to total charged current cross section; the mean energy of neutrinos producing detected neutral pions is 1.1 GeV. The result agrees with the Rein-Sehgal model implemented in our neutrino interaction simulation program with nuclear effects. The spectrum shape of the neutral pion momentum and angle agree with the model. We also measure the ratio of the neutral current coherent pion production to total charged current cross section to be (0.7 \pm 0.4) x 10^(-2).

preprint2010arXiv

Measurement of neutral current coherent neutral pion production on carbon in a few-GeV neutrino beam

The SciBooNE Collaboration reports a measurement of neutral current coherent neutral pion production on carbon by a muon neutrino beam with average energy 0.8 GeV. The separation of coherent from inclusive neutral pion production has been improved by detecting recoil protons from resonant neutral pion production. We measure the ratio of the neutral current coherent neutral pion production to total charged current cross sections to be (1.16 +/- 0.24) x 10-2. The ratio of charged current coherent pion to neutral current coherent pion production is calculated to be 0.14+0.30 -0.28, using our published charged current coherent pion measurement.

preprint2009arXiv

CANGAROO-III Observation of TeV Gamma Rays from the vicinity of PSR B1 706-44

Observation by the CANGAROO-III stereoscopic system of the Imaging Cherenkov Telescope has detected extended emission of TeV gamma rays in the vicinity of the pulsar PSR B1706$-$44. The strength of the signal observed as gamma-ray-like events varies when we apply different ways of emulating background events. The reason for such uncertainties is argued in relevance to gamma-rays embedded in the "off-source data", that is, unknown sources and diffuse emission in the Galactic plane, namely, the existence of a complex structure of TeV gamma-ray emission around PSR B1706$-$44.

preprint2009arXiv

CANGAROO-III search for TeV Gamma-rays from two clusters of galaxies

Because accretion and merger shocks in clusters of galaxies may accelerate particles to high energies, clusters are candidate sites for the origin of ultra-high-energy (UHE) cosmic-rays. A prediction was presented for gamma-ray emission from a cluster of galaxies at a detectable level with the current generation of imaging atmospheric Cherenkov telescopes. The gamma-ray emission was produced via inverse Compton upscattering of cosmic microwave background (CMB) photons by electron-positron pairs generated by collisions of UHE cosmic rays in the cluster. We observed two clusters of galaxies, Abell 3667 and Abell 4038, searching for very-high-energy gamma-ray emission with the CANGAROO-III atmospheric Cherenkov telescope system in 2006. The analysis showed no significant excess around these clusters, yielding upper limits on the gamma-ray emission. From a comparison of the upper limit for the north-west radio relic region of Abell 3667 with a model prediction, we derive a lower limit for the magnetic field of the region of ~0.1 micro G. This shows the potential of gamma-ray observations in studies of the cluster environment. We also discuss the flux upper limit from cluster center regions using a model of gamma-ray emission from neutral pions produced in hadronic collisions of cosmic-ray protons with the intra-cluster medium (ICM). The derived upper limit of the cosmic-ray energy density within this framework is an order of magnitude higher than that of our Galaxy.

preprint2009arXiv

Search for VHE gamma rays from SS433/W50 with the CANGAROO-II telescope

SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate sub-TeV gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 system have been performed to detect sub-TeV gamma-rays using the 10m CANGAROO-II telescope in August and September, 2001, and July and September, 2002. The total observation times are 85.2 hours for ON source, and 80.8 hours for OFF source data. No significant excess of sub-TeV gamma rays has been found at 3 regions of the western X-ray lobe of SS433/W50 system. We have derived 99% confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to TeV. The derived lower limits are 4.3 microgauss for the center of the brightest X-ray emission region and 6.3 microgauss for the far end from SS433 in the western X-ray lobe. In addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe.

preprint2009arXiv

The case for a directional dark matter detector and the status of current experimental efforts

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

preprint2006arXiv

Monte-Carlo Simulator and Ancillary Response Generator of Suzaku XRT/XIS System for Spatially Extended Source Analysis

We have developed a framework for the Monte-Carlo simulation of the X-Ray Telescopes (XRT) and the X-ray Imaging Spectrometers (XIS) onboard Suzaku, mainly for the scientific analysis of spatially and spectroscopically complex celestial sources. A photon-by-photon instrumental simulator is built on the ANL platform, which has been successfully used in ASCA data analysis. The simulator has a modular structure, in which the XRT simulation is based on a ray-tracing library, while the XIS simulation utilizes a spectral "Redistribution Matrix File" (RMF), generated separately by other tools. Instrumental characteristics and calibration results, e.g., XRT geometry, reflectivity, mutual alignments, thermal shield transmission, build-up of the contamination on the XIS optical blocking filters (OBF), are incorporated as completely as possible. Most of this information is available in the form of the FITS (Flexible Image Transport System) files in the standard calibration database (CALDB). This simulator can also be utilized to generate an "Ancillary Response File" (ARF), which describes the XRT response and the amount of OBF contamination. The ARF is dependent on the spatial distribution of the celestial target and the photon accumulation region on the detector, as well as observing conditions such as the observation date and satellite attitude. We describe principles of the simulator and the ARF generator, and demonstrate their performance in comparison with in-flight data.

preprint2003arXiv

Development of a time projection chamber with micro pixel electrodes

A time projection chamber (TPC) based on a gaseous chamber with micro pixel electrodes (micro-PIC) has been developed for measuring three-dimensional tracks of charged particles. The micro-PIC with a detection area of 10 cm square consists of a double-sided printing circuit board. Anode pixels are formed with 0.4 mm pitch on strips aligned perpendicular to the cathode strips in order to obtain a two-dimensional position. In the TPC with drift length of 8 cm, 4 mm wide field cage electrodes are aligned at 1mm spaces and a uniform electric field of about 0.4 kV/cm is produced. For encoding of the three-dimensional position a synchronous readout system has been developed using Field Programmable Gate Arrays with 40 MHz clock. This system enables us to reconstruct the three-dimensional track of the particle at successive points like a cloud chamber even at high event rate. The drift velocity of electrons in the TPC was measured with the tracks of cosmic muons for three days, during which the TPC worked stably with the gas gain of 3000. With a radioisotope of gamma-ray source the three-dimensional track of a Compton scattered electron was taken successfully.

preprint2003arXiv

Performance of large area Micro Pixel Chamber

A novel gaseous two-dimensional imaging detector "Micro Pixel Chamber (micro-PIC)" has been developed. This detector is based on double sided printed circuit board (PCB). We have developed large area (10cm x 10cm) micro-PICs with 65536 pixel anodes of 400um pitch on a 100um thick insulating substrate. Achieved energy resolution was 30% (FWHM) at 5.9keV, and a gas gain of 7000 was obtained with argon ethane (8:2) gas mixture. This gain is high enough to detect minimum ionizing particles with such a small electrode pitch. Although several discharges occurred during 65 hours continuous operation, the detectors have kept stable operation with high gain. The micro-PIC is a useful detector for many applications e.g. X-ray, gamma ray, and charged particle imaging. The micro electrode structure allows us to measure directions of primary electrons due to incident X-rays or gamma rays, which provide a strong method for X-ray polarimetry and gamma-ray imaging.

preprint1998arXiv

ASCA Observations of Blazars and Multiband Analysis

We present data for 18 blazars observed with ASCA, half of which were also observed contemporaneously with EGRET as parts of multi-wavelength campaigns. The observations show a clear difference in the spectra between three subclasses of blazars, namely the High-energy peaked BL Lac objects (HBLs), Low-energy peaked BL Lac objects (LBLs), and quasar-hosted blazars (QHBs). We find that the radiation process responsible for the HE peak for HBLs can be explained solely by Synchrotron-Self-Compton (SSC) emission, with the Doppler factor consistent with the VLBI and/or gamma-ray variability data. For many QHBs, on the other hand, the gamma-rays cannot be solely due to the SSC mechanism unless the Doppler factor is significantly in excess of that inferred from VLBI data. We consider an alternative scenario consistent with the measured values of the Doppler factor, where the SSC component is still present in QHBs and it dominates in the X-ray band, but it is below the observed gamma-ray spectrum. With an assumption that the peak of the SSC emission is on the extrapolation of the X-ray spectrum, and adopting the Doppler factor of 10, we infer the magnetic field to be 0.1 - 1 Gauss, and Lorentz factors of electrons radiating at the peak of the nu F(nu) spectrum of 10^3 for QHBs; this is much lower than 10^5 for HBLs. This difference is most likely due to the large photon density expected in QHBs.

H. Kubo

What is connected

Connect this record

See the researcher in context

Building this map preview

121 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

Constraint on the Matter-Antimatter Symmetry-Violating Phase in Neutrino Oscillations

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

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

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

First combined measurement of the muon neutrino and antineutrino charged-current cross section without pions in the final state at T2K

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

Search for heavy neutrinos with the T2K near detector ND280

Simultaneous measurement of the muon neutrino charged-current cross section on oxygen and carbon without pions in the final state at T2K

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

First Measurement of the Charged Current $\overlineν_μ$ Double Differential Cross Section on a Water Target without Pions in the final state

Measurement of neutrino and antineutrino neutral-current quasielastic-like interactions on oxygen by detecting nuclear de-excitation $γ$-rays

Measurement of the muon neutrino charged-current single $π^+$ production on hydrocarbon using the T2K off-axis near detector ND280

Search for Electron Antineutrino Appearance in a Long-baseline Muon Antineutrino Beam

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

Prospects for CTA observations of the young SNR RX J1713.7-3946

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

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

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

Fermi LAT Discovery of Extended Gamma-Ray Emissions in the Vicinity of the HB3 Supernova Remnant

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

Measurements of $π^{\pm}$ differential yields from the surface of the T2K replica target for incoming 31 GeV/c protons with the NA61/SHINE spectrometer at the CERN SPS

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

Readout technologies for directional WIMP Dark Matter detection

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

An Electron-Tracking Compton Telescope for a Survey of the Deep Universe by MeV gamma-rays

CTA Contributions to the 34th International Cosmic Ray Conference (ICRC2015)

Development and operational experience of magnetic horn system for T2K experiment

Development of Slow Control Boards for the Large Size Telescopes of the Cherenkov Telescope Array

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

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

Measurement of the $ν_μ$ CCQE cross section on carbon with the ND280 detector at T2K

Measurement of the $ν_μ$ charged current quasi-elastic cross-section on carbon with the T2K on-axis neutrino beam

Measurement of the Electron Neutrino Charged-current Interaction Rate on Water with the T2K ND280 pi-zero Detector

Measurements of neutrino oscillation in appearance and disappearance channels by the T2K experiment with 6.6E20 protons on target

Neutrino Oscillation Physics Potential of the T2K Experiment

Performance of a new electron-tracking Compton camera under intense radiations from a water target irradiated with a proton beam

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