Source author record

K. Nishijima

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

128works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

128 published item(s)

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.

preprint2022arXiv

Pre-Supernova Alert System for Super-Kamiokande

In 2020, the Super-Kamiokande (SK) experiment moved to a new stage (SK-Gd) in which gadolinium (Gd) sulfate octahydrate was added to the water in the detector, enhancing the efficiency to detect thermal neutrons and consequently improving the sensitivity to low energy electron anti-neutrinos from inverse beta decay (IBD) interactions. SK-Gd has the potential to provide early alerts of incipient core-collapse supernovae through detection of electron anti-neutrinos from thermal and nuclear processes responsible for the cooling of massive stars before the gravitational collapse of their cores. These pre-supernova neutrinos emitted during the silicon burning phase can exceed the energy threshold for IBD reactions. We present the sensitivity of SK-Gd to pre-supernova stars and the techniques used for the development of a pre-supernova alarm based on the detection of these neutrinos in SK, as well as prospects for future SK-Gd phases with higher concentrations of Gd. For the current SK-Gd phase, high-confidence alerts for Betelgeuse could be issued up to nine hours in advance of the core-collapse itself.

preprint2022arXiv

Search for neutrinoless quadruple beta decay of $^{136}$Xe in XMASS-I

A search for the neutrinoless quadruple beta decay of $^{136}$Xe was conducted with the liquid-xenon detector XMASS-I using $\rm 327\; kg \times 800.0 \; days$ of the exposure. The pulse shape discrimination based on the scintillation decay time constant which distinguishes $γ$-rays including the signal and $β$-rays was used to enhance the search sensitivity. No significant signal excess was observed from the energy spectrum fitting with precise background evaluation, and we set a lower limit of the half life of 3.7 $\times$ 10$^{24}$ years at 90$\%$ confidence level. This is the first experimental constraint of the neutrinoless quadruple beta decay of $^{136}$Xe.

preprint2022arXiv

Search for proton decay via $p\rightarrow μ^+K^0$ in 0.37 megaton-years exposure of Super-Kamiokande

We searched for proton decay via $p\toμ^+K^0$ in 0.37\,Mton$\cdot$years of data collected between 1996 and 2018 from the Super-Kamiokande water Cherenkov experiment. The selection criteria were defined separately for $K^0_S$ and $K^0_L$ channels. No significant event excess has been observed. As a result of this analysis, which extends the previous search by an additional 0.2\,Mton$\cdot$years of exposure and uses an improved event reconstruction, we set a lower limit of $3.6\times10^{33}$ years on the proton lifetime.

preprint2022arXiv

Search for solar electron anti-neutrinos due to spin-flavor precession in the Sun with Super-Kamiokande-IV

Due to a very low production rate of electron anti-neutrinos ($\barν_e$) via nuclear fusion in the Sun, a flux of solar $\barν_e$ is unexpected. An appearance of $\barν_e$ in solar neutrino flux opens a new window for the new physics beyond the standard model. In particular, a spin-flavor precession process is expected to convert an electron neutrino into an electron anti-neutrino (${ν_e\to\barν_e}$) when neutrino has a finite magnetic moment. In this work, we have searched for solar $\barν_e$ in the Super-Kamiokande experiment, using neutron tagging to identify their inverse beta decay signature. We identified 78 $\barν_e$ candidates for neutrino energies of 9.3 to 17.3 MeV in 2970.1 live days with a fiducial volume of 22.5 kiloton water (183.0 kton$\cdot$year exposure). The energy spectrum has been consistent with background predictions and we thus derived a 90% confidence level upper limit of ${4.7\times10^{-4}}$ on the $ν_e\to\barν_e$ conversion probability in the Sun. We used this result to evaluate the sensitivity of future experiments, notably the Super-Kamiokande Gadolinium (SK-Gd) upgrade.

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

First Gadolinium Loading to Super-Kamiokande

In order to improve Super-Kamiokande's neutron detection efficiency and to thereby increase its sensitivity to the diffuse supernova neutrino background flux, 13 tons of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ (gadolinium sulfate octahydrate) was dissolved into the detector's otherwise ultrapure water from July 14 to August 17, 2020, marking the start of the SK-Gd phase of operations. During the loading, water was continuously recirculated at a rate of 60 m$^3$/h, extracting water from the top of the detector and mixing it with concentrated $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ solution to create a 0.02% solution of the Gd compound before injecting it into the bottom of the detector. A clear boundary between the Gd-loaded and pure water was maintained through the loading, enabling monitoring of the loading itself and the spatial uniformity of the Gd concentration over the 35 days it took to reach the top of the detector. During the subsequent commissioning the recirculation rate was increased to 120 m$^3$/h, resulting in a constant and uniform distribution of Gd throughout the detector and water transparency equivalent to that of previous pure-water operation periods. Using an Am-Be neutron calibration source the mean neutron capture time was measured to be $115\pm1$ $μ$s, which corresponds to a Gd concentration of $111\pm2$ ppm, as expected for this level of Gd loading. This paper describes changes made to the water circulation system for this detector upgrade, the Gd loading procedure, detector commissioning, and the first neutron calibration measurements in SK-Gd.

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.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

Neutron-Antineutron Oscillation Search using a 0.37 Megaton$\cdot$Year Exposure of Super-Kamiokande

As a baryon number violating process with $ΔB=2$, neutron-antineutron oscillation ($n\to\bar n$) provides a unique test of baryon number conservation. We have performed a search for $n\to\bar n$ oscillation with bound neutrons in Super-Kamiokande, with the full data set from its first four run periods, representing an exposure of 0.37~Mton-years. The search used a multivariate analysis trained on simulated $n\to\bar n$ events and atmospheric neutrino backgrounds and resulted in 11 candidate events with an expected background of 9.3 events. In the absence of statistically significant excess, we derived a lower limit on $\bar n$ appearance lifetime in $^{16}$O nuclei of $3.6\times{10}^{32}$ years and on the neutron-antineutron oscillation time of $τ_{n\to\bar n} > 4.7\times10^{8}$~s at 90\% C.L..

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 event bursts in XMASS-I associated with gravitational-wave events

We performed a search for event bursts in the XMASS-I detector associated with 11 gravitational-wave events detected during LIGO/Virgo's O1 and O2 periods. Simple and loose cuts were applied to the data collected in the full 832 kg xenon volume around the detection time of each gravitational-wave event. The data were divided into four energy regions ranging from keV to MeV. Without assuming any particular burst models, we looked for event bursts in sliding windows with various time width from 0.02 to 10 s. The search was conducted in a time window between $-$400 and $+$10,000 s from each gravitational-wave event. For the binary neutron star merger GW170817, no significant event burst was observed in the XMASS-I detector and we set 90% confidence level upper limits on neutrino fluence for the sum of all the neutrino flavors via coherent elastic neutrino-nucleus scattering. The obtained upper limit was (1.3-2.1)$\times 10^{11}$ cm$^{-2}$ under the assumption of a Fermi-Dirac spectrum with average neutrino energy of 20 MeV. The neutrino fluence limits for mono-energetic neutrinos in the energy range between 14 and 100 MeV were also calculated. Among the other 10 gravitational wave events detected as the binary black hole mergers, a burst candidate with a 3.0$σ$ significance was found at 1801.95-1803.95 s in the analysis for GW151012. However, no significant deviation from the background in the reconstructed energy and position distributions was found. Considering the additional look-elsewhere effect of analyzing the 11 GW events, the significance of finding such a burst candidate associated with any of them is 2.1$σ$.

preprint2020arXiv

Search for exotic neutrino-electron interactions using solar neutrinos in XMASS-I

We have searched for exotic neutrino-electron interactions that could be produced by a neutrino millicharge, by a neutrino magnetic moment, or by dark photons using solar neutrinos in the XMASS-I liquid xenon detector. We observed no significant signals in 711 days of data. We obtain an upper limit for neutrino millicharge of 5.4$\times$10$^{-12} e$ at 90\% confidence level assuming all three species of neutrino have common millicharge. We also set flavor dependent limits assuming the respective neutrino flavor is the only one carrying a millicharge, $7.3 \times 10^{-12} e$ for $ν_e$, $1.1 \times 10^{-11} e$ for $ν_μ$, and $1.1 \times 10^{-11} e$ for $ν_τ$. These limits are the most stringent yet obtained from direct measurements. We also obtain an upper limit for the neutrino magnetic moment of 1.8$\times$10$^{-10}$ Bohr magnetons. In addition, we obtain upper limits for the coupling constant of dark photons in the $U(1)_{B-L}$ model of 1.3$\times$10$^{-6}$ if the dark photon mass is 1$\times 10^{-3}$ MeV$/c^{2}$, and 8.8$\times$10$^{-5}$ if it is 10 MeV$/c^{2}$.

preprint2020arXiv

Search for proton decay into three charged leptons in 0.37 megaton-years exposure of the Super-Kamiokande

A search for proton decay into three charged leptons has been performed by using 0.37$\,$Mton$\cdot$years of data collected in Super-Kamiokande. All possible combinations of electrons, muons and their anti-particles consistent with charge conservation were considered as decay modes. No significant excess of events has been found over the background, and lower limits on the proton lifetime divided by the branching ratio have been obtained. The limits range between $9.2\times10^{33}$ to $3.4\times10^{34}$ years at 90$\,$% confidence level, improving by more than an order of magnitude upon limits from previous experiments. A first limit has been set for the $p\rightarrowμ^-e^+e^+$ mode.

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

Search for Astronomical Neutrinos from Blazar TXS0506+056 in Super-Kamiokande

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

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 measurement of the time profile of scintillation induced by low energy gamma-rays in liquid xenon with the XMASS-I detector

We report the measurement of the emission time profile of scintillation from gamma-ray induced events in the XMASS-I 832 kg liquid xenon scintillation detector. Decay time constant was derived from a comparison of scintillation photon timing distributions between the observed data and simulated samples in order to take into account optical processes such as absorption and scattering in liquid xenon. Calibration data of radioactive sources, $^{55}$Fe, $^{241}$Am, and $^{57}$Co were used to obtain the decay time constant. Assuming two decay components, $τ_1$ and $τ_2$, the decay time constant $τ_2$ increased from 27.9 ns to 37.0 ns as the gamma-ray energy increased from 5.9 keV to 122 keV. The accuracy of the measurement was better than 1.5 ns at all energy levels. A fast decay component with $τ_1 \sim 2$ ns was necessary to reproduce data. Energy dependencies of $τ_2$ and the fraction of the fast decay component were studied as a function of the kinetic energy of electrons induced by gamma-rays. The obtained data almost reproduced previously reported results and extended them to the lower energy region relevant to direct dark matter searches.

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

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

Direct dark matter search by annual modulation in XMASS-I

A search for dark matter was conducted by looking for an annual modulation signal due to the Earth's rotation around the Sun using XMASS, a single phase liquid xenon detector. The data used for this analysis was 359.2 live days times 832 kg of exposure accumulated between November 2013 and March 2015. When we assume Weakly Interacting Massive Particle (WIMP) dark matter elastically scattering on the target nuclei, the exclusion upper limit of the WIMP-nucleon cross section 4.3$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ was obtained and we exclude almost all the DAMA/LIBRA allowed region in the 6 to 16 GeV/c$^2$ range at $\sim$10$^{-40}$cm$^2$. The result of a simple modulation analysis, without assuming any specific dark matter model but including electron/$γ$ events, showed a slight negative amplitude. The $p$-values obtained with two independent analyses are 0.014 and 0.068 for null hypothesis, respectively. we obtained 90\% C.L. upper bounds that can be used to test various models. This is the first extensive annual modulation search probing this region with an exposure comparable to DAMA/LIBRA.

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 the atmospheric neutrino flux by Super-Kamiokande: energy spectra, geomagnetic effects, and solar modulation

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric $ν_e+{\barν}_e$ and $ν_μ+{\barν}_μ$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the $ν_e$ and $ν_μ$ samples at 8.0 σ and 6.0 σ significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 σ level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 σ level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 σ level.

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

Real-Time Supernova Neutrino Burst Monitor at Super-Kamiokande

We present a real-time supernova neutrino burst monitor at Super-Kamiokande (SK). Detecting supernova explosions by neutrinos in real time is crucial for giving a clear picture of the explosion mechanism. Since the neutrinos are expected to come earlier than light, a fast broadcasting of the detection may give astronomers a chance to make electromagnetic radiation observations of the explosions right at the onset. The role of the monitor includes a fast announcement of the neutrino burst detection to the world and a determination of the supernova direction. We present the online neutrino burst detection system and studies of the direction determination accuracy based on simulations at SK.

preprint2016arXiv

Search for Neutrinos in Super-Kamiokande associated with Gravitational Wave Events GW150914 and GW151226

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

preprint2016arXiv

Search for two-neutrino double electron capture on $^{124}$Xe with the XMASS-I detector

Double electron capture is a rare nuclear decay process in which two orbital electrons are captured simultaneously in the same nucleus. Measurement of its two-neutrino mode would provide a new reference for the calculation of nuclear matrix elements whereas observation of its neutrinoless mode would demonstrate lepton number violation. A search for two-neutrino double electron capture on $^{124}$Xe is performed using 165.9 days of data collected with the XMASS-I liquid xenon detector. No significant excess above background was observed and we set a lower limit on the half-life as $4.7 \times 10^{21}$ years at 90% confidence level. The obtained limit has ruled out parts of some theoretical expectations. We obtain a lower limit on the $^{126}$Xe two-neutrino double electron capture half-life of $4.3 \times 10^{21}$ years at 90% confidence level as well.

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

Blazar Variability and Evolution in the GeV Regime

One of the most important problem of the blazar astrophysics is to understand the physical origin of the blazar sequence. In this study, we focus on the GeV gamma-ray variability of blazars and evolution perspective we search the relation between the redshift and the variability amplitude of blazars for each blazar subclass. We analyzed the Fermi-LAT data of the TeV blazars and the bright AGNs (flux $\geq$ 4$\times10^{-9}$ cm$^{-2}$s$^{-1}$) selected from the 2LAC (the 2nd LAT AGN catalog) data base. As a result, we found a hint of the correlation between the redshift and the variability amplitude in the FSRQs. Furthermore the BL Lacs which have relatively lower peak frequency of the synchrotron radiation and relatively lower redshift, have a tendency to have a smaller variability amplitude.

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

Limits on sterile neutrino mixing using atmospheric neutrinos in Super-Kamiokande

We present limits on sterile neutrino mixing using 4,438 live-days of atmospheric neutrino data from the Super-Kamiokande experiment. We search for fast oscillations driven by an eV$^2$-scale mass splitting and for oscillations into sterile neutrinos instead of tau neutrinos at the atmospheric mass splitting. When performing both these searches we assume that the sterile mass splitting is large, allowing $\sin^2(Δm^2 L/4E)$ to be approximated as $0.5$, and we assume that there is no mixing between electron neutrinos and sterile neutrinos ($|U_{e4}|^2 = 0$). No evidence of sterile oscillations is seen and we limit $|U_{\mu4}|^2$ to less than 0.041 and $|U_{\tau4}|^2$ to less than 0.18 for $Δm^2 > 0.8$ eV$^2$ at the 90% C.L. in a 3+1 framework. The approximations that can be made with atmospheric neutrinos allow these limits to be easily applied to 3+N models, and we provide our results in a generic format to allow comparisons with other sterile neutrino models.

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

Probing Very Early Stage of Radio Source Evolution in NGC 1275 with VERA

We report on VLBI observations of NGC 1275 with VLBI Exploration of Radio Astrometry (VERA) and Very Long Baseline Array (VLBA) during 2008 and 2013 at 15, 22 and 43 GHz. Our observations provide long-term variations in the radio flux and structure of the inner jet in NGC 1275 on sub-parsec scales, particularly at C3 which is a new radio component accountable for the recent active state in radio regime since early-2000s. We found the apparent velocity of C3 was $β_{\rm{app}} = 0.267 \pm 0.007 $, which was consistent with typical velocities of the hot spot and/or radio lobes in the young radio sources. Furthermore, the radio flux and size of C3 had simultaneously increased during September 2007 and May 2012 with keeping optically-thin spectra, suggesting that the particle acceleration occurred inside C3 as generally seen in the hot spots of the radio lobes. Our new observations suggest that C3 is a new-born hot spot and/or radio lobe at a very early stage of its evolution. Our results imply the radio lobes in the radio galaxies might be already formed on sub-parsec scales from the central SMBH and might play an important role to produce $\rm γ$-ray emission in the young radio sources.

preprint2015arXiv

Search for dinucleon decay into pions at Super-Kamiokande

A search for dinucleon decay into pions with the Super-Kamiokande detector has been performed with an exposure of 282.1 kiloton-years. Dinucleon decay is a process that violates baryon number by two units. We present the first search for dinucleon decay to pions in a large water Cherenkov detector. The modes $^{16}$O$(pp) \rightarrow$ $^{14}$C$π^{+}π^{+}$, $^{16}$O$(pn) \rightarrow$ $^{14}$N$π^{+}π^{0}$, and $^{16}$O$(nn) \rightarrow$ $^{14}$O$π^{0}π^{0}$ are investigated. No significant excess in the Super-Kamiokande data has been found, so a lower limit on the lifetime of the process per oxygen nucleus is determined. These limits are: $τ_{pp\rightarrowπ^{+}π^{+}} > 7.22 \times 10^{31}$ years, $τ_{pn\rightarrowπ^{+}π^{0}} > 1.70 \times 10^{32}$ years, and $τ_{nn\rightarrowπ^{0}π^{0}} > 4.04 \times 10^{32}$ years. The lower limits on each mode are about two orders of magnitude better than previous limits from searches for dinucleon decay in iron.

preprint2015arXiv

Search for n-nbar oscillation in Super-Kamiokande

A search for neutron-antineutron ($n-\bar{n}$) oscillation was undertaken in Super-Kamiokande using the 1489 live-day or $2.45 \times 10^{34}$ neutron-year exposure data. This process violates both baryon and baryon minus lepton numbers by an absolute value of two units and is predicted by a large class of hypothetical models where the seesaw mechanism is incorporated to explain the observed tiny neutrino masses and the matter-antimatter asymmetry in the Universe. No evidence for $n-\bar{n}$ oscillation was found, the lower limit of the lifetime for neutrons bound in ${}^{16}$O, in an analysis that included all of the significant sources of experimental uncertainties, was determined to be $1.9 \times 10^{32}$~years at the 90\% confidence level. The corresponding lower limit for the oscillation time of free neutrons was calculated to be $2.7 \times 10^8$~s using a theoretical value of the nuclear suppression factor of $0.517 \times 10^{23}$~s$^{-1}$ and its uncertainty.

preprint2015arXiv

Search for neutrinos from annihilation of captured low-mass dark matter particles in the Sun by Super-Kamiokande

Super-Kamiokande (SK) can search for weakly interacting massive particles (WIMPs) by detecting neutrinos produced from WIMP annihilations occurring inside the Sun. In this analysis, we include neutrino events with interaction vertices in the detector in addition to upward-going muons produced in the surrounding rock. Compared to the previous result, which used the upward-going muons only, the signal acceptances for light (few-GeV/$c^2$ $\sim$ 200-GeV/$c^2$) WIMPs are significantly increased. We fit 3903 days of SK data to search for the contribution of neutrinos from WIMP annihilation in the Sun. We found no significant excess over expected atmospheric-neutrino background and the result is interpreted in terms of upper limits on WIMP-nucleon elastic scattering cross sections under different assumptions about the annihilation channel. We set the current best limits on the spin-dependent (SD) WIMP-proton cross section for WIMP masses below 200 GeV/$c^2$ (at 10 GeV/$c^2$, 1.49$\times 10^{-39}$ cm$^2$ for $χχ\rightarrow b \bar{b}$ and 1.31$\times 10^{-40}$ cm$^2$ for $χχ\rightarrowτ^+τ^-$ annihilation channels), also ruling out some fraction of WIMP candidates with spin-independent (SI) coupling in the few-GeV/$c^2$ mass range.

preprint2015arXiv

Search for Nucleon and Dinucleon Decays with an Invisible Particle and a Charged Lepton in the Final State at the Super-Kamiokande Experiment

Search results for nucleon decays $p \rightarrow e^+X$, $p \rightarrow μ^+X$, $n \rightarrow νγ$ (where $X$ is an invisible, massless particle) as well as dinucleon decays $np \rightarrow e^+ν$, $np \rightarrow μ^+ν$ and $np \rightarrow τ^+ν$ in the Super-Kamiokande experiment are presented. Using single-ring data from an exposure of 273.4 kton $\cdot$ years, a search for these decays yields a result consistent with no signal. Accordingly, lower limits on the partial lifetimes of $τ_{p \rightarrow e^+X} > 7.9 \times 10^{32}$ years, $τ_{p \rightarrow μ^+X} > 4.1 \times 10^{32}$ years, $τ_{n \rightarrow νγ} > 5.5 \times 10^{32}$ years, $τ_{np \rightarrow e^+ν} > 2.6 \times 10^{32}$ years, $τ_{np \rightarrow μ^+ν} > 2.2 \times 10^{32}$ years and $τ_{np \rightarrow τ^+ν} > 2.9 \times 10^{31}$ years at a $90 \% $ confidence level are obtained. Some of these searches are novel.

preprint2015arXiv

Test of Lorentz invariance with atmospheric neutrinos

A search for neutrino oscillations induced by Lorentz violation has been performed using 4,438 live-days of Super-Kamiokande atmospheric neutrino data. The Lorentz violation is included in addition to standard three-flavor oscillations using the non-perturbative Standard Model Extension (SME), allowing the use of the full range of neutrino path lengths, ranging from 15 to 12,800 km, and energies ranging from 100 MeV to more than 100 TeV in the search. No evidence of Lorentz violation was observed, so limits are set on the renormalizable isotropic SME coefficients in the $eμ$, $μτ$, and $eτ$ sectors, improving the existing limits by up to seven orders of magnitude and setting limits for the first time in the neutrino $μτ$ sector of the SME.

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

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

First Indication of Terrestrial Matter Effects on Solar Neutrino Oscillation

We report an indication that the elastic scattering rate of solar $^8$B neutrinos with electrons in the Super-Kamiokande detector is larger when the neutrinos pass through the Earth during nighttime. We determine the day/night asymmetry, defined as the difference of the average day rate and average night rate divided by the average of those two rates, to be $(-3.2\pm1.1(\text{stat})\pm0.5(\text{syst}))\%$, which deviates from zero by 2.7 $σ$. Since the elastic scattering process is mostly sensitive to electron-flavored solar neutrinos, a non-zero day/night asymmetry implies that the flavor oscillations of solar neutrinos are affected by the presence of matter within the neutrinos' flight path. Super-Kamiokande's day/night asymmetry is consistent with neutrino oscillations for $4\times10^{-5}$eV$^2\leqΔm^2_{21}\leq7\times10^{-5}$eV$^2$ and large mixing values of $θ_{12}$, at the $68\%$ C.L.

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

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

Search for bosonic superweakly interacting massive dark matter particles with the XMASS-I detector

Bosonic superweakly interacting massive particles (super-WIMPs) are a candidate for warm dark matter. With the absorption of such a boson by a xenon atom these dark matter candidates would deposit an energy equivalent to their rest mass in the detector. This is the first direct detection experiment exploring the vector super-WIMPs in the mass range between 40 and 120 keV. Using 165.9 days of data no significant excess above background was observed in the fiducial mass of 41 kg. The present limit for the vector super-WIMPs excludes the possibility that such particles constitute all of dark matter. The absence of a signal also provides the most stringent direct constraint on the coupling constant of pseudoscalar super-WIMPs to electrons. The unprecedented sensitivity was achieved exploiting the low background at a level $10^{-4}$ kg$^{-1}$keV$_{ee}^{-1}$day$^{-1}$ in the detector.

preprint2014arXiv

Search for inelastic WIMP nucleus scattering on $^{129}$Xe in data from the XMASS-I experiment

A search for inelastic scattering of Weakly Interacting Massive Particles (WIMPs) on the isotope $^{129}$Xe was done in data taken with the single phase liquid xenon detector XMASS at the Kamioka Observatory. Using a restricted volume containing 41 kg of LXe at the very center of our detector we observed no significant excess of events in 165.9 days of data. Our background reduction allowed us to derive our limits without explicitly subtracting the remaining events which are compatible with background expectations and derive for e.g. a 50 GeV WIMP an upper limit for its inelastic cross section on $^{129}$Xe nuclei of 3.2 pb at the 90% confidence level.

preprint2014arXiv

Search for Trilepton Nucleon Decay via $p \rightarrow e^+ νν$ and $p \rightarrow μ^+ νν$ in the Super-Kamiokande Experiment

The trilepton nucleon decay modes $p \rightarrow e^+ νν$ and $p \rightarrow μ^+ νν$ violate $|Δ(B - L)|$ by two units. Using data from a 273.4 kiloton year exposure of Super-Kamiokande a search for these decays yields a fit consistent with no signal. Accordingly, lower limits on the partial lifetimes of $τ_{p \rightarrow e^+ νν} > 1.7 \times 10^{32}$ years and $τ_{p \rightarrow μ^+ νν} > 2.2 \times 10^{32}$ years at a $90 \% $ confidence level are obtained. These limits can constrain Grand Unified Theories which allow for such processes.

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

A Search for Nucleon Decay via $n \rightarrow \barν π^{0}$ and $p \rightarrow \barν π^{+}$ in Super-Kamiokande

We present the results of searches for nucleon decay via bound neutron to antineutrino plus pizero and proton to antineutrino plus piplus using data from a combined 172.8 kiloton-years exposure of Super-Kamiokande-I, -II, and -III. We set lower limits on the partial lifetime for each of these modes. For antineutrino pizero, the partial lifetime is >1.1x10^{33} years; for antineutrino piplus, the partial lifetime is >3.9x10^{32} years at 90% confidence level.

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

Light WIMP search in XMASS

A search for light dark matter using low-threshold data from the single phase liquid xenon scintillation detector XMASS, has been conducted. Using the entire 835 kg inner volume as target, the analysis threshold can be lowered to 0.3 keVee (electron-equivalent) to search for light dark matter. With low-threshold data corresponding to a 5591.4 kg$\cdot$day exposure of the detector and without discriminating between nuclear-recoil and electronic events, XMASS excludes part of the parameter space favored by other experiments.

preprint2013arXiv

Search for solar axions in XMASS, a large liquid-xenon detector

XMASS, a low-background, large liquid-xenon detector, was used to search for solar axions that would be produced by bremsstrahlung and Compton effects in the Sun. With an exposure of 5.6ton days of liquid xenon, the model-independent limit on the coupling for mass $\ll$ 1keV is $|g_{aee}|< 5.4\times 10^{-11}$ (90% C.L.), which is a factor of two stronger than the existing experimental limit. The bounds on the axion masses for the DFSZ and KSVZ axion models are 1.9 and 250eV, respectively. In the mass range of 10-40keV, this study produced the most stringent limit, which is better than that previously derived from astrophysical arguments regarding the Sun to date.

preprint2012arXiv

An Indirect Search for WIMPs in the Sun using 3109.6 days of upward-going muons in Super-Kamiokande

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

preprint2012arXiv

Search for GUT Monopoles at Super-Kamiokande

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce ν_{e}, ν_μ and \barν_μ. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(σ_0/1 mb) < 6.3 \times 10^{-24} (β_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where β_M is the monopole velocity in units of the speed of light and σ_0 is the catalysis cross section at β_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for β_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

preprint2012arXiv

Search for Nucleon Decay into Charged Anti-lepton plus Meson in Super-Kamiokande I and II

Searches for a nucleon decay into a charged anti-lepton (e^+ or μ^+) plus a light meson (π^0, π^-, η, ρ^0, ρ^-, ω) were performed using the Super-Kamiokande I and II data. Twelve nucleon decay modes were searched for. The total exposure is 140.9 kiloton \cdot years, which includes a 91.7 kiloton \cdot year exposure (1489.2 live days) of Super-Kamiokande-I and a 49.2 kiloton \cdot year exposure (798.6 live days) of Super-Kamiokande-II. The number of candidate events in the data was consistent with the atmospheric neutrino background expectation. No significant evidence for a nucleon decay was observed in the data. Thus, lower limits on the nucleon partial lifetime at 90% confidence level were obtained. The limits range from 3.6 \times 10^31 to 8.2 \times 10^33 years, depending on the decay modes.

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.

preprint2012arXiv

Searches for very high energy gamma rays from blazars with CANGAROO-III telescope in 2005-2009

We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in the field of view. No significant VHE gamma ray emission was detected from any of the four blazars. The GeV gamma-ray spectra of these objects were obtained by analyzing Fermi/LAT archival data. Non-simultaneous wide range (radio to VHE gamma-ray bands) spectral energy distributions (SEDs) including CANGAROO-III upper limits, GeV gamma-ray spectra, and archival data are discussed using a one-zone synchrotron self-Compton (SSC) model in combination with a external Compton (EC) radiation. The HBLs (H 2356-309 and PKS 2155-304) can be explained by a simple SSC model, and PKS 0537-441 and 3C 279 are well modeled by a combination of SSC and EC model. We find a consistency with the blazar sequence in terms of strength of magnetic field and component size.

preprint2012arXiv

Solar neutrino results in Super-Kamiokande-III

The results of the third phase of the Super-Kamiokande solar neutrino measurement are presented and compared to the first and second phase results. With improved detector calibrations, a full detector simulation, and improved analysis methods, the systematic uncertainty on the total neutrino flux is estimated to be ?2.1%, which is about two thirds of the systematic uncertainty for the first phase of Super-Kamiokande. The observed 8B solar flux in the 5.0 to 20 MeV total electron energy region is 2.32+/-0.04 (stat.)+/-0.05 (sys.) *10^6 cm^-2sec^-1, in agreement with previous measurements. A combined oscillation analysis is carried out using SK-I, II, and III data, and the results are also combined with the results of other solar neutrino experiments. The best-fit oscillation parameters are obtained to be sin^2 θ12 = 0.30+0.02-0.01(tan^2 θ12 = 0.42+0.04 -0.02) and Δm2_21 = 6.2+1.1-1.9 *10^-5eV^2. Combined with KamLAND results, the best-fit oscillation parameters are found to be sin^2 θ12 = 0.31+/-0.01(tan^2 θ12 = 0.44+/-0.03) and Δm2_21 = 7.6?0.2*10^-5eV^2 . The 8B neutrino flux obtained from global solar neutrino experiments is 5.3+/-0.2(stat.+sys.)*10^6cm^-2s^-1, while the 8B flux becomes 5.1+/-0.1(stat.+sys.)*10^6cm^-2s^-1 by adding KamLAND result. In a three-flavor analysis combining all solar neutrino experiments, the upper limit of sin^2 θ13 is 0.060 at 95% C.L.. After combination with KamLAND results, the upper limit of sin^2 θ13 is found to be 0.059 at 95% C.L..

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

Letter of Intent: The Hyper-Kamiokande Experiment --- Detector Design and Physics Potential ---

We propose the Hyper-Kamiokande (Hyper-K) detector as a next generation underground water Cherenkov detector. It will serve as a far detector of a long baseline neutrino oscillation experiment envisioned for the upgraded J-PARC, and as a detector capable of observing -- far beyond the sensitivity of the Super-Kamiokande (Super-K) detector -- proton decays, atmospheric neutrinos, and neutrinos from astronomical origins. The baseline design of Hyper-K is based on the highly successful Super-K, taking full advantage of a well-proven technology. (to be continued)

preprint2011arXiv

Scintillation-only Based Pulse Shape Discrimination for Nuclear and Electron Recoils in Liquid Xenon

In a dedicated test setup at the Kamioka Observatory we studied pulse shape discrimination (PSD) in liquid xenon (LXe) for dark matter searches. PSD in LXe was based on the observation that scintillation light from electron events was emitted over a longer period of time than that of nuclear recoil events, and our method used a simple ratio of early to total scintillation light emission in a single scintillation event. Requiring an efficiency of 50% for nuclear recoil retention we reduced the electron background to 7.7\pm1.1(stat)\pm1.2 0.6(sys)\times10-2 at energies between 4.8 and 7.2 keVee and to 7.7\pm2.8(stat)\pm2.5 2.8(sys)\times10-3 at energies between 9.6 and 12 keVee for a scintillation light yield of 20.9 p.e./keV. Further study was done by masking some of that light to reduce this yield to 4.6 p.e./keV, the same method results in an electron event reduction of 2.4\pm0.2(stat)\pm0.3 0.2(sys)\times10-1 for the lower of the energy regions above. We also observe that in contrast to nuclear recoils the fluctuations in our early to total ratio for electron events are larger than expected from statistical fluctuations.

preprint2011arXiv

Search for Differences in Oscillation Parameters for Atmospheric Neutrinos and Antineutrinos at Super-Kamiokande

We present a search for differences in the oscillations of antineutrinos and neutrinos in the Super-Kamiokande -I, -II, and -III atmospheric neutrino sample. Under a two-flavor disappearance model with separate mixing parameters between neutrinos and antineutrinos, we find no evidence for a difference in oscillation parameters. Best fit antineutrino mixing is found to be at (dm2bar, sin2 2 thetabar) = (2.0x10^-3 eV^2, 1.0) and is consistent with the overall Super-K measurement.

preprint2011arXiv

Study of Non-Standard Neutrino Interactions with Atmospheric Neutrino Data in Super-Kamiokande I and II

In this paper we study non-standard neutrino interactions as an example of physics beyond the standard model using atmospheric neutrino data collected during the Super-Kamiokande I(1996-2001) and II(2003-2005) periods. We focus on flavor-changing-neutral-currents (FCNC), which allow neutrino flavor transitions via neutral current interactions, and effects which violate lepton non-universality (NU) and give rise to different neutral-current interaction-amplitudes for different neutrino flavors. We obtain a limit on the FCNC coupling parameter, varepsilon_{mu tau}, |varepsilon_{mu tau}|<1.1 x 10^{-2} at 90%C.L. and various constraints on other FCNC parameters as a function of the NU coupling, varepsilon_{e e}. We find no evidence of non-standard neutrino interactions in the Super-Kamiokande atmospheric data.

preprint2011arXiv

Supernova Relic Neutrino Search at Super-Kamiokande

A new Super-Kamiokande (SK) search for Supernova Relic Neutrinos (SRNs) was conducted using 2853 live days of data. Sensitivity is now greatly improved compared to the 2003 SK result, which placed a flux limit near many theoretical predictions. This more detailed analysis includes a variety of improvements such as increased efficiency, a lower energy threshold, and an expanded data set. New combined upper limits on SRN flux are between 2.8 and 3.0 nu_e cm^-2 s^-1 > 16 MeV total positron energy (17.3 MeV E_nu).

preprint2011arXiv

Towards an Axiomatic Formulation of Noncommutative Quantum Field Theory

We propose new Wightman functions as vacuum expectation values of products of field operators in the noncommutative space-time. These Wightman functions involve the $\star$-product among the fields, compatible with the twisted Poincaré symmetry of the noncommutative quantum field theory (NC QFT). In the case of only space-space noncommutativity ($θ_{0i}=0$), we prove the CPT theorem using the noncommutative form of the Wightman functions. We also show that the spin-statistics theorem, demonstrated for the simplest case of a scalar field, holds in NC QFT within this formalism.

preprint2010arXiv

Atmospheric neutrino oscillation analysis with sub-leading effects in Super-Kamiokande I, II, and III

We present a search for non-zero theta_{13} and deviations of sin^2 theta_{23} from 0.5 in the oscillations of atmospheric neutrino data from Super-Kamiokande -I, -II, and -III. No distortions of the neutrino flux consistent with non-zero theta_{13} are found and both neutrino mass hierarchy hypotheses are in agreement with the data. The data are best fit at Delta m^2 = 2.1 x 10^-3 eV^2, sin^2 theta_{13} = 0.0, and sin^2 theta_{23} =0.5. In the normal (inverted) hierarchy theta_{13} and Delta m^2 are constrained at the one-dimensional 90% C.L. to sin^2 theta_{13} < 0.04 (0.09) and 1.9 (1.7) x 10^-3 < Delta m^2 < 2.6 (2.7) x 10^-3 eV^2. The atmospheric mixing angle is within 0.407 <= sin^2 theta_{23} <= 0.583 at 90% C.L.

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

Self-shielding effect of a single phase liquid xenon detector for direct dark matter search

Liquid xenon is a suitable material for a dark matter search. For future large scale experiments, single phase detectors are attractive due to their simple configuration and scalability. However, in order to reduce backgrounds, they need to fully rely on liquid xenon's self-shielding property. A prototype detector was developed at Kamioka Observatory to establish vertex and energy reconstruction methods and to demonstrate the self-shielding power against gamma rays from outside of the detector. Sufficient self-shielding power for future experiments was obtained.

preprint2008arXiv

Noncommutative Quantum Field Theory: A Confrontation of Symmetries

The concept of a noncommutative field is formulated based on the interplay between twisted Poincaré symmetry and residual symmetry of the Lorentz group. Various general dynamical results supporting this construction, such as the light-wedge causality condition and the integrability condition for Tomonaga-Schwinger equation, are presented. Based on this analysis, the claim of the identity between commutative QFT and noncommutative QFT with twisted Poincaré symmetry is refuted.

preprint2007arXiv

CANGAROO-III Observations of the 2006 Outburst of PKS2155-304

We have used the CANGAROO-III imaging atmospheric Cherenkov telescopes to observe the high-frequency-peaked BL Lacertae (HBL) object PKS 2155-304 between 2006 July 28 (MJD 53944) and August 2, triggered by the H.E.S.S. report that the source was in a high state of TeV gamma-ray emission. A signal was detected at the 4.8 sigma level in an effective live time of 25.1 hours during the outburst period. The flux of Very High Energy gammarays from the CANGAROO-III observations shows the variability on the time scale of less than a few hours. The averaged integral flux above 660 GeV is (1.6 +/- 0.3_stat +/- 0.5_syst) x 10^-11 cm^-2 sec^-1 which corresponds to ~45% of the flux observed from the Crab nebula. Follow-up observations between August 17 (MJD 53964) and 25 indicate the source activity had decreased.

preprint2004arXiv

Search for Dark Matter WIMPs using Upward Through-going Muons in Super-Kamiokande

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

preprint2001arXiv

Renormalization Constant of the Color Gauge Field as a Probe of Confinement

The mechanism of color confinement as a consequence of an unbroken non-abelian gauge symmetry and asymptotic freedom is elucidated and compared with that of other models based on an analogy with the type II superconductor. It is demonstrated that a sufficient condition for color confinement is given by $Z_3^{-1}=0$ where $Z_3$ denotes the renormalization constant of the color gauge field. It is shown that this condition is actually satisfied in quantum chromodynamics and that some of the characteristic features of other models follow from it.

preprint1999arXiv

TeV gamma-ray observations of three X-ray selected BL Lacs

Despite extensive surveys of extragalactic TeV gamma-ray candidates only 3 sources have so far been detected. All three are northern hemisphere objects and all three are low-redshift X-ray selected BL Lacs (XBLs). In this paper we present the results of observations of the three nearest southern hemisphere XBLs (PKS0548-322, PKS2005-489 and PKS2155-304) with the CANGAROO 3.8m imaging telescope. During the period of observation we estimate that the threshold of the 3.8m telescope was around 1.5TeV. Searches for both steady and short timescale emission have been performed for each source. Additionally, we are able to monitor the X-ray state of each source on a daily basis and we have made contemporaneous measurements of optical activity for PKS0548-322 and PKS2155-304.

preprint1999arXiv

TeV observations of Centaurus A

We have searched for TeV gamma-rays from Centaurus A and surrounding region out to +/- 1.0 deg using the CANGAROO 3.8m telescope. No evidence for TeV gamma-ray emission was observed from the search region, which includes a number of interesting features located away from the tracking centre of our data. The 3 sigma upper limit to the flux of gamma-rays above 1.5 TeV from an extended source of radius 14' centred on Centaurus A is 1.28e-11 photons cm^-2 s^-1.

preprint1998arXiv

Discovery of TeV Gamma Rays from SN1006: Further Evidence for the SNR Origin of Cosmic Rays

This paper reports the first discovery of TeV gamma-ray emission from a supernova remnant made with the CANGAROO 3.8 m Telescope. TeV gamma rays were detected at the sky position and extension coincident with the north-east (NE) rim of shell-type Supernova remnant (SNR) SN1006 (Type Ia). SN1006 has been a most likely candidate for an extended TeV Gamma-ray source, since the clear synchrotron X-ray emission from the rims was recently observed by ASCA (Koyama et al. 1995), which is a strong evidence of the existence of very high energy electrons up to hundreds of TeV in the SNR. The observed TeV gamma-ray flux was $(2.4\pm 0.5(statistical) \pm 0.7(systematic)) \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ ($\ge 3.0\pm 0.9$ TeV) and $ (4.6\pm 0.6 \pm 1.4) \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ ($\ge 1.7\pm 0.5$ TeV) from the 1996 and 1997 observations, respectively. Also we set an upper limit on the TeV gamma-ray emission from the SW rim, estimated to be $ 1.1 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ ($\ge 1.7\pm 0.5$ TeV, 95% CL) in the 1997 data. The TeV gamma rays can be attributed to the 2.7 K cosmic background photons up-scattered by electrons of energies up to about 10$^{14}$ eV by the inverse Compton (IC) process. The observed flux of the TeV gamma rays, together with that of the non-thermal X-rays, gives firm constraints on the acceleration process in the SNR shell; a magnetic field of $6.5\pm2$ $μ$G is inferred from both the synchrotron X-rays and inverse Compton TeV gamma-rays, which gives entirely consistent mechanisms that electrons of energies up to 10$^{14}$ eV are produced via the shock acceleration in SN1006.

preprint1998arXiv

TeV gamma-ray observations of southern BL Lacs with the CANGAROO 3.8m Imaging Telescope

Observational and theoretical results indicate that low-redshift BL Lacertae objects are the most likely extragalactic sources to be detectable at TeV energies. In this paper we present the results of observations of 4 BL Lacertae objects (PKS0521-365, EXO0423.4-0840, PKS2005-489 and PKS2316-423) made between 1993 and 1996 with the CANGAROO 3.8m imaging Cherenkov telescope. During the period of these observations the gamma-ray energy threshold of the 3.8m telescope was ~2TeV. Searches for steady long-term emission have been made, and, inspired by the TeV flares detected from Mkn421 and Mkn501, a search on a night-by-night timescale has also been performed for each source. Comprehensive Monte Carlo simulations are used to estimate upper limits for both steady and short timescale emission.

preprint1997arXiv

Detection of Gamma Rays of Up to 50 TeV From the Crab Nebula

Gamma rays with energies greater than 7 TeV from the Crab pulsar/nebula have been observed at large zenith angles, using the Imaging Atmospheric Technique from Woomera, South Australia. CANGAROO data taken in 1992, 1993 and 1995 indicate that the energy spectrum extends up to at least 50 TeV, without a change of the index of the power law spectrum. The observed differential spectrum is \noindent $(2.01\pm 0.36)\times 10^{-13}(E/{7 TeV})^{-2.53 \pm 0.18} TeV^{-1}cm^{-2}s^{-1}$ between 7 TeV and 50 TeV. There is no apparent cut-off. The spectrum for photon energies above $\sim$10 TeV allows the maximum particle acceleration energy to be inferred, and implies that this unpulsed emission does not originate near the light cylinder of the pulsar, but in the nebula where the magnetic field is not strong enough to allow pair creation from the TeV photons. The hard gamma-ray energy spectrum above 10 TeV also provides information about the varying role of seed photons for the inverse Compton process at these high energies, as well as a possible contribution of $π^{\circ}$-gamma rays from proton collisions.

preprint1997arXiv

Very High Energy Gamma Rays from the Vela Pulsar Direction

We have observed the Vela pulsar region at TeV energies using the 3.8 m imaging Cherenkov telescope near Woomera, South Australia between January 1993 and March 1995. Evidence of an unpulsed gamma-ray signal has been detected at the 5.8 sigma level. The detected gamma-ray flux is (2.9 +/- 0.5 +/- 0.4) x 10^-12 photons cm^-2 sec^-1 above 2.5 +/- 1.0 TeV and the signal is consistent with steady emission over the two years. The gamma-ray emission region is offset from the Vela pulsar position to the southeast by about 0.13 deg. No pulsed emission modulated with the pulsar period has been detected and the 95 % confidence flux upper limit to the pulsed emission from the pulsar is (3.7 +/- 0.7) x 10^-13 photons cm^-2 sec^-1 above 2.5 +/- 1.0 TeV.