Compilation of papers presented by the JEM-EUSO Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

Mini-EUSO will observe the Earth in the UV range (300 - 400 nm) offering the opportunity to study a variety of atmospheric events such as Transient Luminous Events (TLEs), meteors and marine bioluminescence. Furthermore it aims to search for Ultra High Energy Cosmic Rays (UHECR) above $10^{21}$ eV and Strange Quark Matter (SQM). The detector is expected to be launched to the International Space Station in August 2019 and look at the Earth in nadir mode from the UV-transparent window of the Zvezda module of the International Space Station. The instrument comprises a compact telescope with a large field of view ($44^{\circ}$), based on an optical system employing two Fresnel lenses for light collection. The light is focused onto an array of 36 multi-anode photomultiplier tubes (MAPMT), for a total of 2304 pixels and the resulting signal is converted into digital, processed and stored via the electronics subsystems on-board. In addition to the main detector, Mini-EUSO contains two ancillary cameras for complementary measurements in the near infrared (1500 - 1600 nm) and visible (400 - 780 nm) range and also a 8x8 SiPM imaging array.

JEM-EUSO is an international program for the development of space based Ultra High Cosmic Ray observatories. The program consists of a series of missions which are either under development or in the data analysis phase. The instruments consist of wide field of view telescopes, which operate in the UV range, designed to detect the fluorescence light emitted by extended air showers in the atmosphere. We describe the simulation software ESAF and show the physical assumptions done in it. We present here the implementation of the POEMMA, K-EUSO, TUS, Mini-EUSO, EUSO-SPB and EUSO-TA detectors in the simulation software ESAF.

The JEM-EUSO program aims at detecting Ultra High Energy Cosmic Rays (UHECRs) by observing the fluorescence light produced by extensive air showers (EAS) in the Earth's atmosphere. Within this program, a new generation of missions is being built, including (i) Mini-EUSO that will be installed on board the ISS in August 2019, (ii) an upgrade to the ground-based telescope EUSO-TA and (iii) the second super pressure balloon flight (EUSO-SPB2). All these detectors will have a dedicated trigger system based on a board equipped with a Xilinx Zynq device that will be able to detect different types of events on three different time-scales: a microsecond time-scale for cosmic ray detection (L1), a hundreds of microsecond time-scale for slower events like transient luminous events (TLEs) (L2), and a tens of millisecond time-scale used to produce a continuous monitoring, for even slower events like meteors or nuclearites. The L1 trigger logic for the upgrade of EUSO-TA and EUSO-SPB2 are being developed taking into account the peculiarity of each detector (optic system, FOV, frame length) starting from the logic already developed for Mini-EUSO. In particular, every pixel will have an independent threshold that will be dynamically adapted to the level of the background; a predetermined condition on the number, the position and the time distribution of pixels above threshold has to be satisfied in order to issue a trigger. This contribution will summarize the L1 trigger logics and the tests currently performed.

Space Debris (SD) consist of non-operational artificial objects orbiting around the Earth, which could possibly damage space vehicles, such as the International Space Station (ISS) or other manned spacecrafts. The vast majority of such objects are cm-sized, not catalogued and usually the tracking data are not precise enough. Here we present the feasibility study of SD detection and tracking with techniques usually employed in cosmic-ray physics. For this purpose, we have evaluated the possibility of using Mini-EUSO, a space-borne fluorescence telescope to be deployed on the ISS, to track SD illuminated by the Sun. By means of ESAF (EUSO Simulation and analysis Framework) simulation and by developing the trigger algorithms, we estimated the minimum size and maximum distances of detectable SD. We then studied the number of possible SD detections using an ESA software called MASTER (Meteoroid and SD Terrestrial Environment Reference). With the Mini-EUSO Engineering Model (Mini-EUSO EM), we performed some measurements to estimate the reflectance of the most common SD materials and to demonstrate the ability of Mini-EUSO to detect SD events. We also performed some tests in open-sky conditions, identifying and tracking fast-moving objects. In particular, the detection of a rocket body allowed us to confirm the simulation outcomes predictions and the expected performance of the detector.

Mini-EUSO is a UV telescope that will look downwards to the Earth's atmosphere onboard the International Space Station. With the design of the ultra-high energy cosmic ray fluorescence detectors belonging to the JEM-EUSO program, it will make the first UV map of the Earth by observing atmospheric phenomena such as transient luminous events, sprites and lightning, as well as meteors and bioluminescence from earth. Diffused light from laser shots from the ground, which mimic the fluorescence light emitted by Nitrogen molecules when extensive air showers pass through the atmosphere, can be used to verify the capability of this kind of detector to observe ultra-high energy cosmic rays. To validate the electronics and the trigger algorithms developed for Mini-EUSO, a scaled down version of the telescope with 1:9 of the original focal surface and a lens of 2.5 cm diameter has been built. Tests of the Mini-EUSO engineering model have been made in laboratory and in open sky condition. In this paper, we report results of observations of the night sky, which include the detection of stars, meteors, a planet and a rocket body reflecting the sunlight. Interesting results of the observation of city lights are also reported.

The TurLab facility is a laboratory, equipped with a 5 m diameter and 1 m depth rotating tank, located in the Physics Department of the University of Turin. Originally, it was mainly built to study systems of different scales where rotation plays a key role in the fluid behavior such as in atmospheric and oceanic flows. In the past few years the TurLab facility has been used to perform experiments related to the observation of Extreme Energy Cosmic Rays (EECRs) from space using the fluorescence technique. For example, in the case of the JEM-EUSO mission, where the diffuse night brightness and artificial light sources can vary significantly in time and space inside the Field of View of the telescope. The Focal Surface of Mini-EUSO Engineering Model (Mini-EUSO EM) with the level 1 (L1) and 2 (L2) trigger logics implemented in the Photo-Detector Module (PDM) has been tested at TurLab. Tests related to the possibility of using an EUSO-like detector for other type of applications such as Space Debris (SD) monitoring and imaging detector have also been pursued. The tests and results obtained within the EUSO@TurLab Project on these different topics are presented.

162173 Ryugu, the target of Hayabusa2, has a round shape with an equatorial ridge, which is known as a spinning top-shape. A strong centrifugal force is a likely contributor to Ryugu's top-shaped features. Observations by Optical Navigation Camera onboard Hayabusa2 show a unique longitudinal variation in geomorphology; the western side of this asteroid, later called the western bulge, has a smooth surface and a sharp equatorial ridge, compared to the other side. Here, we propose a structural deformation process that generated the western bulge. Applying the mission-derived shape model, we employ a finite element model technique to analyze the locations that experience structural failure within the present shape. Assuming that materials are uniformly distributed, our model shows the longitudinal variation in structurally failed regions when the spin period is shorter than ~3.75 h. Ryugu is structurally intact in the subsurface region of the western bulge while other regions are sensitive to structural failure. We infer that this variation is indicative of the deformation process that occurred in the past, and the western bulge is more relaxed structurally than the other region. Our analysis also shows that this deformation process might occur at a spin period between ~3.5 h and ~3.0 h, providing the cohesive strength ranging between ~4 Pa and ~10 Pa.

EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on.

The Mini-EUSO telescope is designed by the JEM-EUSO Collaboration to observe the UV emission of the Earth from the vantage point of the International Space Station (ISS) in low Earth orbit. The main goal of the mission is to map the Earth in the UV, thus increasing the technological readiness level of future EUSO experiments and to lay the groundwork for the detection of Extreme Energy Cosmic Rays (EECRs) from space. Due to its high time resolution of 2.5 us, Mini-EUSO is capable of detecting a wide range of UV phenomena in the Earth's atmosphere. In order to maximise the scientific return of the mission, it is necessary to implement a multi-level trigger logic for data selection over different timescales. This logic is key to the success of the mission and thus must be thoroughly tested and carefully integrated into the data processing system prior to the launch. This article introduces the motivation behind the trigger design and details the integration and testing of the logic.

The TurLab facility is a laboratory, equipped with a 5 m diameter and 1 m depth rotating tank, located in the Physics Department of the University of Turin. The tank has been built mainly to study problems where system rotation plays a key role in the fluid behaviour such as in atmospheric and oceanic flows at different scales. The tank can be filled with different fluids of variable density, which enables studies in layered conditions such as sea waves. The tank can be also used to simulate the terrestrial surface with the optical characteristics of different environments such as snow, grass, ocean, land with soil, stones etc., fogs and clouds. As it is located in an extremely dark place, the light intensity can be controlled artificially. Such capabilities of the TurLab facility are applied to perform experiments related to the observation of Extreme Energy Cosmic Rays (EECRs) from space using the fluorescence technique, as in the case of the JEM-EUSO mission, where the diffuse night brightness and artificial light sources can vary significantly in time and space inside the Field of View (FoV) of the telescope. Here we will report the currently ongoing activity at the TurLab facility in the framework of the JEM-EUSO mission (EUSO@TurLab).

The regolith-covered surfaces of asteroids preserve records of geophysical processes that have occurred both at their surfaces and sometimes also in their interiors. As a result of the unique micro-gravity environment that these bodies posses, a complex and varied geophysics has given birth to fascinating features that we are just now beginning to understand. The processes that formed such features were first hypothesised through detailed spacecraft observations and have been further studied using theoretical, numerical and experimental methods that often combine several scientific disciplines. These multiple approaches are now merging towards a further understanding of the geophysical states of the surfaces of asteroids. In this chapter we provide a concise summary of what the scientific community has learned so far about the surfaces of these small planetary bodies and the processes that have shaped them. We also discuss the state of the art in terms of experimental techniques and numerical simulations that are currently being used to investigate regolith processes occurring on small-body surfaces and that are contributing to the interpretation of observations and the design of future space missions.

We presented optical and near-infrared multi-band linear polarimetry of the highly reddened Type Ia SN~2014J appeared in M82. SN~2014J exhibits large polarization at shorter wavelengths, e.g., $4.8$\% in $B$ band, and the polarization decreases rapidly at longer wavelengths, with the position angle of the polarization remaining at approximately $40^{\circ}$ over the observed wavelength range. These polarimetric properties suggest that the observed polarization is likely to be caused predominantly by the interstellar dust within M82. Further analysis shows that the polarization peaks at a wavelengths much shorter than those obtained for the Galactic dust. The wavelength dependence of the polarization can be better described by an inverse power law rather than by Serkowski law for Galactic interstellar polarization. These suggests that the nature of the dust in M82 may be different from that in our Galaxy, with polarizing dust grains having a mean radius of $<0.1\ \mu$m.

The gamma-ray-detected blazar 3C 454.3 exhibits dramatic flux and polarization variations in the optical and near-infrared bands. In December 2010, the object emitted a very bright outburst. We monitored it for approximately four years (including the 2010 outburst) by optical and near-infrared photopolarimetry. During the 2010 outburst, the object emitted two rapid, redder brightenings, at which the polarization degrees (PDs) in both bands increased significantly and the bands exhibited a frequency-dependent polarization. The observed frequency-dependent polarization leads us to propose that the polarization vector is composed of two vectors. Therefore, we separate the observed polarization vectors into short and long-term components that we attribute to the emissions of the rapid brightenings and the outburst that varied the timescale of days and months, respectively. The estimated PD of the short-term component is greater than the maximum observed PD and is close to the theoretical maximum PD. We constrain the bulk Lorentz factors and inclination angles between the jet axis and the line of sight from the estimated PDs. In this case, the inclination angle of the emitting region of short-term component from the first rapid brightening should be equal to 90$^{\circ}$, because the estimated PD of the short-term component was approximately equal to the theoretical maximum PD. Thus, the Doppler factor at the emitting region of the first rapid brightening should be equal to the bulk Lorentz factor.

We present the results of the search for a correlation between giant radio pulses (GRPs) at 1.4 GHz and hard X-rays at 15-75 keV from the Crab pulsar. We made simultaneous ground and satellite observations of the Crab pulsar over 12 hours in three occasions in April 2010, March and September 2011, and got a sample of 1.3*10^4 main-pulse phase GRPs. From these samples we have found statistically marginal enhancement (21.5%, 2.70 sigma) of hard X-ray flux within +/- 1.5 degree phase angle of the synchronous peak of main-pulse phase GRPs. This enhancement, if confirmed, implicates that GRPs may accompany plasma density increases in the pulsar magnetosphere.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

Context: The blazar Markarian 421 is one of the brightest TeV gamma-ray sources of the northern sky. From December 2007 until June 2008 it was intensively observed in the very high energy (VHE, E > 100 GeV) band by the single-dish Major Atmospheric Gamma-ray Imaging Cherenkov telescope (MAGIC-I). Aims: We aimed to measure the physical parameters of the emitting region of the blazar jet during active states. Methods: We performed a dense monitoring of the source in VHE with MAGIC-I, and also collected complementary data in soft X-rays and optical-UV bands; then, we modeled the spectral energy distributions (SED) derived from simultaneous multi-wavelength data within the synchrotron self--compton (SSC) framework. Results: The source showed intense and prolonged gamma-ray activity during the whole period, with integral fluxes (E > 200 GeV) seldom below the level of the Crab Nebula, and up to 3.6 times this value. Eight datasets of simultaneous optical-UV (KVA, Swift/UVOT), soft X-ray (Swift/XRT) and MAGIC-I VHE data were obtained during different outburst phases. The data constrain the physical parameters of the jet, once the spectral energy distributions obtained are interpreted within the framework of a single-zone SSC leptonic model. Conclusions: The main outcome of the study is that within the homogeneous model high Doppler factors (40 <= delta <= 80) are needed to reproduce the observed SED; but this model cannot explain the observed short time-scale variability, while it can be argued that inhomogeneous models could allow for less extreme Doppler factors, more intense magnetic fields and shorter electron cooling times compatible with hour or sub-hour scale variability.

We present the results of a long M87 monitoring campaign in very high energy $\gamma$-rays with the MAGIC-I Cherenkov telescope. We aim to model the persistent non-thermal jet emission by monitoring and characterizing the very high energy $\gamma$-ray emission of M87 during a low state. A total of 150\,h of data were taken between 2005 and 2007 with the single MAGIC-I telescope, out of which 128.6\,h survived the data quality selection. We also collected data in the X-ray and \textit{Fermi}--LAT bands from the literature (partially contemporaneous). No flaring activity was found during the campaign. The source was found to be in a persistent low-emission state, which was at a confidence level of $7\sigma$. We present the spectrum between 100\,GeV and 2\,TeV, which is consistent with a simple power law with a photon index $\Gamma=2.21\pm0.21$ and a flux normalization at 300\,GeV of $(7.7\pm1.3) \times 10^{-8}$ TeV$^{-1}$ s$^{-1}$ m$^{-2}$. The extrapolation of the MAGIC spectrum into the GeV energy range matches the previously published \textit{Fermi}--LAT spectrum well, covering a combined energy range of four orders of magnitude with the same spectral index. We model the broad band energy spectrum with a spine layer model, which can satisfactorily describe our data.

We report that the optical polarization in the afterglow of GRB 091208B is measured at t = 149 - 706 s after the burst trigger, and the polarization degree is P = 10.4% +/- 2.5%. The optical light curve at this time shows a power-law decay with index -0.75 +/- 0.02, which is interpreted as the forward shock synchrotron emission, and thus this is the first detection of the early-time optical polarization in the forward shock (rather than that in the reverse shock reported by Steele et al. (2009). This detection disfavors the afterglow model in which the magnetic fields in the emission region are random on the plasma skin depth scales, such as amplified by the plasma instabilities, e.g., Weibel instability. We suggest that the fields are amplified by the magnetohydrodynamic instabilities, which would be tested by future observations of the temporal changes of the polarization degrees and angles for other bursts.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

We study the non-thermal jet emission of the BL Lac object B3 2247+381 during a high optical state. The MAGIC telescopes observed the source during 13 nights between September 30th and October 30th 2010, collecting a total of 14.2 hours of good quality very high energy (VHE) $\gamma$-ray data. Simultaneous multiwavelength data was obtained with X-ray observations by the Swift satellite and optical R-band observations at the KVA-telescope. We also use high energy $\gamma$-ray (HE, 0.1 GeV-100 GeV) data from the Fermi satellite. The BL Lac object B3 2247+381 (z=0.119) was detected, for the first time, at VHE $\gamma$-rays at a statistical significance of 5.6 $\sigma$. A soft VHE spectrum with a photon index of -3.2 $\pm$ 0.6 was determined. No significant short term flux variations were found. We model the spectral energy distribution using a one-zone SSC-model, which can successfully describe our data.

We present the results of five years (2005-2009) of MAGIC observations of the BL Lac object PG 1553+113 at very high energies (VHEs, E > 100 GeV). Power law fits of the individual years are compatible with a steady mean photon index \Gamma = 4.27 $\pm$ 0.14. In the last three years of data, the flux level above 150 GeV shows a clear variability (probability of constant flux < 0.001%). The flux variations are modest, lying in the range from 4% to 11% of the Crab Nebula flux. Simultaneous optical data also show only modest variability that seems to be correlated with VHE gamma ray variability. We also performed a temporal analysis of (all available) simultaneous Fermi/LAT data of PG 1553+113 above 1 GeV, which reveals hints of variability in the 2008-2009 sample. Finally, we present a combination of the mean spectrum measured at very high energies with archival data available for other wavelengths. The mean spectral energy distribution can be modeled with a one-zone Synchrotron Self Compton (SSC) model, which gives the main physical parameters governing the VHE emission in the blazar jet.

We present very high energy (VHE, E > 100 GeV) {\gamma}-ray observations of the {\gamma}-ray binary system LS I+61 303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We detect a 6.3{\sigma} {\gamma}-ray signal above 400 GeV in the combined data set. The integral flux above an energy of 300 GeV is F(E>300 GeV)=(1.4 +- 0.3stat +- 0.4syst) * 10^{-12} cm^{-2} s^{-1}, which corresponds to about 1.3% of the Crab Nebula flux in the same energy range. The orbit-averaged flux of LS I +61 303 in the orbital phase interval 0.6--0.7, where a maximum of the TeV flux is expected, is lower by almost an order of magnitude compared to our previous measurements between 2005 September and 2008 January. This provides evidence for a new low emission state in LS I +61 303. We find that the change to the low flux state cannot be solely explained by an increase of photon-photon absorption around the compact star.

We report the results of the observation of the nearby satellite galaxy Segue 1 performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the background. Differential upper limits on the gamma-ray flux are derived assuming various power-law slopes for the possible emission spectrum. Integral upper limits are also calculated for several power-law spectra and for different energy thresholds. The values are of the order of 10^{-11} ph cm^{-2}$ s^{-1} above 100 GeV and 10^{-12} ph cm^{-2} s^{-1} above 200 GeV. Segue 1 is currently considered one of the most interesting targets for indirect dark matter searches. In these terms, the upper limits have been also interpreted in the context of annihilating dark matter particles. For such purpose, we performed a grid scan over a reasonable portion of the parameter space for the minimal SuperGravity model and computed the flux upper limit for each point separately, taking fully into account the peculiar spectral features of each model. We found that in order to match the experimental upper limits with the model predictions, a minimum flux boost of 10^{3} is required, and that the upper limits are quite dependent on the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we compared the upper limits with the predictions of some dark matter models able to explain the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate annihilating into tau+tau-. A complete exclusion however is not possible due to the uncertainties in the Segue 1 astrophysical factor.

The acceleration of particles up to GeV or higher energies in microquasars has been the subject of considerable theoretical and observational efforts in the past few years. Sco X-1 is a microquasar from which evidence of highly energetic particles in the jet has been found when it is in the so-called Horizontal Branch (HB), a state when the radio and hard X-ray fluxes are higher and a powerful relativistic jet is present. Here we present the first very high energy gamma-ray observations of Sco X-1 obtained with the MAGIC telescopes. An analysis of the whole dataset does not yield a significant signal, with 95% CL flux upper limits above 300 GeV at the level of 2.4x10^{-12} ph/cm^2/s. Simultaneous RXTE observations were conducted to search for TeV emission during particular X-ray states of the source. A selection of the gamma-ray data obtained during the HB based on the X-ray colors did not yield a signal either, with an upper limit of 3.4x10^{-12} ph/cm^2/s. These upper limits place a constraint on the maximum TeV luminosity to non-thermal X-ray luminosity of L_{VHE}/L_{ntX}<0.02, that can be related to a maximum TeV luminosity to jet power ratio of L_{VHE}/L_{j}<10^{-3}. Our upper limits indicate that the underlying high-energy emission physics in Sco X-1 must be inherently different from that of the hitherto detected gamma-ray binaries.

We report on the observation of the region around supernova remnant G65.1+0.6 with the stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL J1954.3+2836 and 1FGL J1958.6+2845. They are identified as GeV pulsars and both have a possible counterpart detected at about 35 TeV by the Milagro observatory. MAGIC collected 25.5 hours of good quality data, and found no significant emission in the range around 1 TeV. We therefore report differential flux upper limits, assuming the emission to be point-like (<0.1 deg) or within a radius of 0.3 deg. In the point-like scenario, the flux limits around 1 TeV are at the level of 3 % and 2 % of the Crab Nebula flux, for the two sources respectively. This implies that the Milagro emission is either extended over a much larger area than our point spread function, or it must be peaked at energies beyond 1 TeV, resulting in a photon index harder than 2.2 in the TeV band.

Very High Energy (VHE) gamma-ray emission from the flat spectrum radio quasar (FSRQ) PKS 1222+21 (4C 21.35, z=0.432) was detected with the MAGIC Cherenkov telescopes during a short observation (~0.5 hr) performed on 2010 June 17. The MAGIC detection coincides with high energy MeV/GeV gamma-ray activity measured by the Large Area Telescope (LAT) on board the Fermi satellite. The VHE spectrum measured by MAGIC extends from about 70 GeV up to at least 400 GeV and can be well described by a power law dN/dE \propto E^-Gamma with a photon index Gamma= 3.75+/-0.27stat +/-0.2syst. The averaged integral flux above 100 GeV is (4.56+/-0.46)x10^(-10) cm^-2 s^-1 (~1 Crab Nebula flux). The VHE flux measured by MAGIC varies significantly within the 30 min exposure implying a flux doubling time of about 10 min. The VHE and MeV/GeV spectra, corrected for the absorption by the extragalactic background light (EBL), can be described by a single power law with photon index 2.72+/-0.34 between 3 GeV and 400 GeV, and is consistent with emission belonging to a single component in the jet. The absence of a spectral cutoff constrains the gamma-ray emission region outside the Broad Line Region, which would otherwise absorb the VHE gamma-rays. Together with the detected fast variability, this challenges present emission models from jets in FSRQ. Moreover, the combined Fermi/LAT and MAGIC spectral data yield constraints on the density of the Extragalactic Background Light in the UV-optical to near-infrared range that are compatible with recent models.

Context. 3C 279, the first quasar discovered to emit VHE gamma-rays by the MAGIC telescope in 2006, was reobserved by MAGIC in January 2007 during a major optical flare and from December 2008 to April 2009 following an alert from the Fermi space telescope on an exceptionally high gamma -ray state. Aims. The January 2007 observations resulted in a detection on January 16 with significance 5.2 sigma, corresponding to a F(> 150 GeV) (3.8 \pm 0.8) \cdot 10^-11 ph cm^-2 s^-1 while the overall data sample does not show significant signal. The December 2008 - April 2009 observations did not detect the source. We study the multiwavelength behavior of the source at the epochs of MAGIC observations, collecting quasi-simultaneous data at optical and X-ray frequencies and for 2009 also gamma-ray data from Fermi. Methods. We study the light curves and spectral energy distribution of the source. The spectral energy distributions of three observing epochs (including the February 2006, which has been previously published in Albert et al. 2008a) are modeled with one-zone inverse Compton models and the emission on January 16, 2007 also with two zone model and with a lepto-hadronic model. Results. We find that the VHE gamma-ray emission detected in 2006 and 2007 challenges standard one-zone model, based on relativistic electrons in a jet scattering broad line region photons, while the other studied models fit the observed spectral energy distribution more satisfactorily.

Between 2004 and 2009 a sample of 28 X-ray selected high- and intermediate-frequency peaked blazars with a X-ray flux larger than 2 uJy at 1 keV in the redshift range from 0.018 to 0.361 was observed with the MAGIC telescope at energies above 100 GeV. Seven among them were detected and the results of these observations are discussed elsewhere. Here we concentrate on the remaining 21 blazars which were not detected during this observation campaign and present the 3 sigma (99.7 %) confidence upper limits on their flux. The individual flux upper limits lie between 1.6 % and 13.6 % of the integral flux from the Crab Nebula. Applying a stacking method to the sample of non-detections with a total of 394.1 hours exposure time, we find evidence for an excess with a cumulative significance of 4.9 standard deviations. It is not dominated by individual objects or flares, but increases linearly with the observation time as for a constant source with an integral flux level of ~1.5 % of that observed from the Crab Nebula above 150 GeV.

We report new observations of the intermediate-frequency peaked BL Lacertae object 3C 66A with the MAGIC telescopes. The data sample we use were taken in 2009 December and 2010 January, and comprises 2.3 hr of good quality data in stereoscopic mode. In this period, we find a significant signal from the direction of the blazar 3C 66A. The new MAGIC stereoscopic system is shown to play an essential role for the separation between 3C 66A and the nearby radio galaxy 3C 66B, which is at a distance of only $6^\prime$. The derived integral flux above $100\eh{GeV}$ is 8.3% of Crab Nebula flux and the energy spectrum is reproduced by a power law of photon index $3.64 \pm 0.39_{\rm stat} \pm 0.25_{\rm sys}$. Within errors, this is compatible with the one derived by VERITAS in 2009. From the spectra corrected for absorption by the extragalactic background light, we only find small differences between the four models that we applied, and constrain the redshift of the blazar to $z < 0.68$.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

We report on the detection with the MAGIC telescopes of very high energy gamma-rays from IC 310, a head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC at a high statistical significance of 7.6sigma in 20.6 hr of stereo data. The observed spectral energy distribution is flat with a differential spectral index of -2.00 \pm 0.14. The mean flux above 300 GeV, between October 2009 and February 2010, (3.1 \pm 0.5)x10^{-12} cm^{-2} s^{-1}, corresponds to (2.5 \pm 0.4)% of Crab Nebula units. Only an upper limit, of 1.9% of Crab Nebula units above 300 GeV, was obtained with the 2008 data. This, together with strong hints (>3sigma) of flares in the middle of October and November 2009, implies that the emission is variable. The MAGIC results favour a scenario with the very high energy emission originating from the inner jet close to the central engine. More complicated models than a simple one-zone SSC scenario, e.g. multi-zone SSC, external Compton or hadronic, may be required to explain the very flat spectrum and its extension over more than three orders of magnitude in energy.

Cygnus X-3 is a microquasar consisting of an accreting compact object orbiting around a Wolf-Rayet star. It has been detected at radio frequencies and up to high-energy gamma rays (above 100 MeV). However, many models also predict a very high energy (VHE) emission (above hundreds of GeV) when the source displays relativistic persistent jets or transient ejections. Therefore, detecting such emission would improve the understanding of the jet physics. The imaging atmospheric Cherenkov telescope MAGIC observed Cygnus X-3 for about 70 hours between 2006 March and 2009 August in different X-ray/radio spectral states and also during a period of enhanced gamma-ray emission. MAGIC found no evidence for a VHE signal from the direction of the microquasar. An upper limit to the integral flux for energies higher than 250 GeV has been set to 2.2 x 10-12 photons cm-2 s-1 (95% confidence level). This is the best limit so far to the VHE emission from this source. The non-detection of a VHE signal during the period of activity in the high-energy band sheds light on the location of the possible VHE radiation favoring the emission from the innermost region of the jets, where absorption is significant. The current and future generations of Cherenkov telescopes may detect a signal under precise spectral conditions.

The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope participated in three multiwavelength (MWL) campaigns, observing the blazar Markarian (Mkn) 421 during the nights of 2006 April 28, 29, and 2006 June 14. We analyzed the corresponding MAGIC very-high energy observations during 9 nights from 2006 April 22 to 30 and on 2006 June 14. We inferred light curves with sub-day resolution and night-by-night energy spectra. A strong gamma-ray signal was detected from Mkn 421 on all observation nights. The flux (E > 250 GeV) varied on night-by-night basis between (0.92+-0.11)10^-10 cm^-2 s^-1 (0.57 Crab units) and (3.21+-0.15)10^-10 cm^-2 s^-1 (2.0 Crab units) in 2006 April. There is a clear indication for intra-night variability with a doubling time of 36+-10(stat) minutes on the night of 2006 April 29, establishing once more rapid flux variability for this object. For all individual nights gamma-ray spectra could be inferred, with power-law indices ranging from 1.66 to 2.47. We did not find statistically significant correlations between the spectral index and the flux state for individual nights. During the June 2006 campaign, a flux substantially lower than the one measured by the Whipple 10-m telescope four days later was found. Using a log-parabolic power law fit we deduced for some data sets the location of the spectral peak in the very-high energy regime. Our results confirm the indications of rising peak energy with increasing flux, as expected in leptonic acceleration models.

Context: Part of the very high energy $\gamma$-ray radiation coming from extragalactic sources is absorbed through the pair production process on the extragalactic background light photons. Extragalactic magnetic fields alter the trajectories of these cascade pairs and, in turn, convert cosmic background photons to gamma-ray energies by inverse Compton scattering. These secondary photons can form an extended halo around bright VHE sources. Aims: We searched for an extended emission around the bright blazars Mrk 421 and Mrk 501 using the MAGIC telescope data. Methods: If extended emission is present, the angular distribution of reconstructed gamma-ray arrival directions around the source is broader than for a point-like source. In the analysis of a few tens of hours of observational data taken from Mrk 421 and Mrk 501 we used a newly developed method that provides better angular resolution. This method is based on the usage of multidimensional decision trees. Comparing the measured shapes of angular distributions with those expected from a point-like source one can detect or constrain possible extended emission around the source. We also studied the influence of different types of systematic errors on the shape of the distribution of reconstructed gamma-ray arrival directions for a point source. Results: We present upper limits for an extended emission calculated for both sources for various source extensions and emission profiles. We discuss possible constraints on the extragalactic magnetic fields strength. We obtain upper limits on the extended emission around the Mrk~421 (Mrk~501) on the level of < 5% (< 4%) of the Crab Nebula flux above the energy threshold of 300 GeV. Using these results we discuss possible constraints on the extragalactic magnetic fields strength around a few times $10^{-15}$~G.

The blazar PG 1553+113 is a well known TeV gamma-ray emitter. In this paper, we determine its spectral energy distribution using simultaneous multi-frequency data in order to study its emission processes. An extensive campaign was carried out between March and April 2008, where optical, X-ray, high-energy (HE) gamma-ray, and very-high-energy (VHE) gamma-ray data were obtained with the KVA, Abastumani, REM, RossiXTE/ASM, AGILE and MAGIC telescopes, respectively. This is the first simultaneous broad-band (i.e., HE+VHE) gamma-ray observation, though AGILE did not detect the source. We combine data to derive source's spectral energy distribution and interpret its double peaked shape within the framework of a synchrotron self compton model

Context: Gamma-ray bursts are cosmological sources emitting radiation from the gamma-rays to the radio band. Substantial observational efforts have been devoted to the study of gamma-ray bursts during the prompt phase, i.e. the initial burst of high-energy radiation, and during the long-lasting afterglows. In spite of many successes in interpreting these phenomena, there are still several open key questions about the fundamental emission processes, their energetics and the environment. Aim: Independently of specific gamma-ray burst theoretical recipes, spectra in the GeV/TeV range are predicted to be remarkably simple, being satisfactorily modeled with power-laws, and therefore offer a very valuable tool to probe the extragalactic background light distribution. Furthermore, the simple detection of a component at very-high energies, i.e. at $\sim 100$\,GeV, would solve the ambiguity about the importance of various possible emission processes, which provide barely distinguishable scenarios at lower energies. Methods: We used the results of the MAGIC telescope observation of the moderate resdhift ($z\sim0.76$) \object{GRB\,080430} at energies above about 80\,GeV, to evaluate the perspective for late-afterglow observations with ground based GeV/TeV telescopes. Results: We obtained an upper limit of $F_{\rm 95%\,CL} = 5.5 \times 10^{-11}$\,erg\,cm$^{-2}$\,s$^{-1}$ for the very-high energy emission of \object{GRB\,080430}, which cannot set further constraints on the theoretical scenarios proposed for this object also due to the difficulties in modeling the low-energy afterglow. Nonetheless, our observations show that Cherenkov telescopes have already reached the required sensitivity to detect the GeV/TeV emission of GRBs at moderate redshift ($z \lesssim 0.8$), provided the observations are carried out at early times, close to the onset of their afterglow phase.

IceCube Collaboration Contributions to the 2009 International Cosmic Ray Conference

Over 5,000 PMTs are being deployed at the South Pole to compose the IceCube neutrino observatory. Many are placed deep in the ice to detect Cherenkov light emitted by the products of high-energy neutrino interactions, and others are frozen into tanks on the surface to detect particles from atmospheric cosmic ray showers. IceCube is using the 10-inch diameter R7081-02 made by Hamamatsu Photonics. This paper describes the laboratory characterization and calibration of these PMTs before deployment. PMTs were illuminated with pulses ranging from single photons to saturation level. Parameterizations are given for the single photoelectron charge spectrum and the saturation behavior. Time resolution, late pulses and afterpulses are characterized. Because the PMTs are relatively large, the cathode sensitivity uniformity was measured. The absolute photon detection efficiency was calibrated using Rayleigh-scattered photons from a nitrogen laser. Measured characteristics are discussed in the context of their relevance to IceCube event reconstruction and simulation efforts.

The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies.

We report the observation of a new dwarf nova, OT J055717+683226, during its first-ever recorded superoutburst in December 2006. Our observation shows that this object is an SU UMa-type dwarf nova having a very short superhump period of 76.67+/- 0.03 min (0.05324+/-0.00002 d). The next superoutburst was observed in March 2008. The recurrence time of superoutbursts (supercycle) is, hence, estimated to be ~480 d. The supercycle is much shorter than those of WZ Sge-type dwarf novae having supercycles of >~ 10 yr, which are a major population of dwarf novae in the shortest orbital period regime (<~85 min). Using a hierarchical cluster analysis, we identified seven groups of dwarf novae in the shortest orbital period regime. We identified a small group of objects that have short supercycles, small outburst amplitudes, and large superhump period excesses, compared with those of WZ Sge stars. OT J055717+683226 probably belongs to this group.

The discovery of very high energy (VHE) gamma-ray emitting X-ray binaries has triggered an intense effort to better understand the particle acceleration, absorption, and emission mechanisms in compact binary systems, which provide variable conditions along eccentric orbits. Despite this, the nature of some of these systems, and of the accelerated particles producing the VHE emission, is unclear. To answer some of these open questions, we conducted a multiwavelength campaign of the VHE gamma-ray emitting X-ray binary LS I +61 303 including the MAGIC telescope, XMM-Newton, and Swift during 60% of an orbit in 2007 September. We detect a simultaneous outburst at X-ray and VHE bands, with the peak at phase 0.62 and a similar shape at both wavelengths. A linear fit to the simultaneous X-ray/VHE pairs obtained during the outburst yields a correlation coefficient of r=0.97, while a linear fit to all simultaneous pairs provides r=0.81. Since a variable absorption of the VHE emission towards the observer is not expected for the data reported here, the correlation found indicates a simultaneity in the emission processes. Assuming that they are dominated by a single particle population, either hadronic or leptonic, the X-ray/VHE flux ratio favors leptonic models. This fact, together with the detected photon indices, suggests that in LS I +61 303 the X-rays are the result of synchrotron radiation of the same electrons that produce VHE emission as a result of inverse Compton scattering of stellar photons.

We present early phase observations in optical and near-infrared wavelengths for the extremely luminous Type Ia supernova (SN Ia) 2009dc. The decline rate of the light curve is $\Delta m_{15}(B)=0.65\pm 0.03$, which is one of the slowest among SNe Ia. The peak $V$-band absolute magnitude is $M_{V}=-19.90\pm 0.15$ mag even if the host extinction is $A_{V}=0$ mag. It reaches $M_{V}=-20.19\pm 0.19$ mag for the host extinction of $A_{V}=0.29$ mag as inferred from the observed Na {\sc i} D line absorption in the host. Our $JHK_{s}$-band photometry shows that the SN is one of the most luminous SNe Ia also in near-infrared wavelengths. These results indicate that SN 2009dc belongs to the most luminous class of SNe Ia, like SN 2003fg and SN 2006gz. We estimate the ejected $^{56}$Ni mass of $1.2\pm 0.3$ $\Msun$ for no host extinction case (or 1.6$\pm$ 0.4 M$_{\odot}$ for the host extinction of $A_{V}=0.29$ mag). The C {\sc ii} $\lambda$6580 absorption line keeps visible until a week after maximum, which diminished in SN 2006gz before its maximum brightness. The line velocity of Si {\sc ii} $\lambda$6355 is about 8000 km s$^{-1}$ around the maximum, being considerably slower than that of SN 2006gz, while comparable to that of SN 2003fg. The velocity of the C {\sc ii} line is almost comparable to that of the Si {\sc ii}. The presence of the carbon line suggests that thick unburned C+O layers remain after the explosion. SN 2009dc is a plausible candidate of the super-Chandrasekhar mass SNe Ia.

Index entry listing the contributed papers of the MAGIC collaboration to the 31th International Cosmic Ray Conference (ICRC 2009), July 7-15 2009, \L\'od\'z, Poland. The individual papers are sorted by subject: Overview and Highlight Papers; MAGIC-II Status and Components; Software and Analysis Techniques; Technical Developments; Scientific Results. This HTML document includes clickable links to the papers that exist on the astro-ph arXiv. We hope that this will make it easy to access the MAGIC contributions in a systematic way.

Geigermode avalanche photodiodes (G-APD) are novel photodetectors, which can detect single photons. This type of diodes might become an alternative to photomultipliers (PMT) in next-generation Imaging Air Cherenkov Telescopes. Prospects, limitations and development directions are be discussed. Results from first tests are reported.

We present the results of numerical calculations on propagation of Extremely-high energy (EHE) neutrinos and charged leptons in the earth for trajectories in all phase space of nadir angles. Our comprehensive calculation has shown that not only the secondary produced muons but also taus survive without decaying in energy range of 10PeV-100PeV with intensity approximately three orders of magnitude lower than the neutrino flux regardless of EHE neutrino production models. They form detectable horizontal or downgoing events in a 1km$^3$ underground neutrino telescope such as the IceCube detector. The event rate and the resultant detectability of EHE signals in comparison with the atmospheric muon background are also evaluated. The 90 % C.L. upperlimit of EHE neutrino fluxes by a km$^2$ detection area would be placed at $E^2dF/dE \simeq 3.7\times 10^{-8}$ GeV/cm$^2$ sec sr for $\nu_{\mu}$ and $4.6\times 10^{-8}$ for $\nu_{\tau}$ with energies of $10^9$ GeV in absence of signals with energy-loss in a detection volume of 10PeV or greater.