Source author record

A. Neronov

A. Neronov 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

89works
10topics
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

89 published item(s)

preprint2026arXiv

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Identifying Galactic PeVatrons (PeV particle accelerators) from the ultra-high-energy (UHE, >100 TeV) $γ$-ray sources plays a crucial role in revealing the origin of Galactic cosmic rays. The UHE source 1LHAASO J1857+0203u is suggested to be associated with HESS J1858+020, which may be attributed to the possible PeVatron candidate supernova remnant (SNR) G35.6$-$0.4 or HII region G35.6$-$0.5. We perform detailed analysis on the very-high-energy and UHE $γ$-ray emissions towards this region with data from the Large High Altitude Air Shower Observatory (LHAASO). 1LHAASO J1857+0203u is detected with a significance of 11.6$σ$ above 100 TeV, indicating the presence of a PeVatron. It has an extension of $\sim 0.18^\circ$ with a power-law (PL) spectral index of $\sim$2.5 in 1-25 TeV and a point-like emission with a PL spectral index of $\sim$3.2 above 25 TeV. Using the archival CO and HI data, we identify some molecular and atomic clouds that may be associated with the TeV $γ$-ray emissions. Our modelling indicates that the TeV $γ$-ray emissions are unlikely to arise from the clouds illuminated by the protons that escaped from SNR G35.6$-$0.4. In the scenario that HII region G35.6$-$0.5 could accelerate particles to the UHE band, the observed GeV-TeV $γ$-ray emission could be well explained by a hadronic model with a PL spectral index of $\sim$2.0 and cutoff energy of $\sim$450 TeV. However, an evolved pulsar wind nebula origin cannot be ruled out.

preprint2026arXiv

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

preprint2026arXiv

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $γ$-rays detectable by arrays such as LHAASO. These $γ$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $pγ$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $γ$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $γ$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $γ$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

preprint2026arXiv

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude $b\approx 10.5^\circ$. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21$σ$ and 17$σ$ at 8$-$100 TeV and $>$100 TeV, respectively. The corresponding extensions are determined to be 0.23$^{\circ}\pm$0.03$^{\circ}$ and 0.17$^{\circ}\pm$0.03$^{\circ}$. The emission is proposed to originate from the relativistic electrons and positrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function $dN/dE = (42.4\pm4.1)(\frac{E}{20\rm\ TeV})^{-2.31\pm0.11}\exp(-\frac{E}{110\pm25\rm\ TeV})$ $\rm\ TeV^{-1}\ cm^{-2}\ s^{-1}$in the energy range from 8 TeV to 300 TeV, implying a steady-state parent electron spectrum $dN_e/dE_e\propto (\frac{E_e}{100\rm\ TeV})^{-3.13\pm0.16}\exp[(\frac{-E_e}{373\pm70\rm\ TeV})^2]$ at energies above $\approx 50 \rm\ TeV$. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to $\approx 4.5\rm\ μG$. To satisfy the multi-wavelength spectrum and the $γ$-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

preprint2026arXiv

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range.

preprint2026arXiv

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories.

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

preprint2026arXiv

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of $8.2σ$ for $E>25$~TeV assuming a Gaussian template. The best-fit location is (R.A., Dec.)$=(95^{\circ}\!.47\pm0^{\circ}\!.11,\,37^{\circ}\!.92 \pm0^{\circ}\!.09)$, and the extension is $0^{\circ}\!.40\pm0^{\circ}\!.07$. The energy spectrum can be described by a power-law spectrum with an index of ${-2.92 \pm 0.17_{\rm stat} \pm 0.02_{\rm sys} }$. No clear extended multi-wavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of $\sim160$~TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

preprint2026arXiv

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments.

preprint2026arXiv

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

preprint2026arXiv

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

The diffuse Galactic gamma-ray emission is a very important tool used to study the propagation and interaction of cosmic rays in the Milky Way. In this work, we report the measurements of the diffuse emission from the Galactic plane, covering Galactic longitudes from $15^{\circ}$ to $235^{\circ}$ and latitudes from $-5^{\circ}$ to $+5^{\circ}$, in an energy range of 1 TeV to 25 TeV, with the Water Cherenkov Detector Array (WCDA) of the Large High Altitude Air Shower Observatory (LHAASO). After masking the sky regions of known sources, the diffuse emission is detected with $24.6σ$ and $9.1σ$ significance in the inner Galactic plane and outer Galactic plane, respectively. The WCDA spectra in both regions can be well described by a power-law function, with spectral indices of $-2.67\pm0.05_{\rm stat}$ in the inner region and $-2.83\pm0.19_{\rm stat}$ in the outer region, respectively. Combined with the Square Kilometer Array (KM2A) measurements at higher energies, a clear softening of the spectrum is found in the inner region, with change of spectral indices by $\sim0.5$ at a break energy around $30$ TeV. The fluxes of the diffuse emission are higher by a factor of $1.5-2.7$ than the model prediction assuming local CR spectra and the gas column density, which are consistent with those measured by the KM2A. Along Galactic longitude, the spatial distribution of the diffuse emission shows deviation from that of the gas column density. The spectral shape of the diffuse emission are possibly variation in different longitude region. The WCDA measurements bridge the gap between the low-energy measurements by space detectors and the ultra-high-energy observations by LHAASO-KM2A and other experiments. These results suggest that improved modeling of the wide-band diffuse emission is required.

preprint2026arXiv

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties.

preprint2025arXiv

Observation of the $γ$-ray Emission from W43 with LHAASO

In this paper, we report the detection of the very-high-energy (VHE, $ 100{\rm\ GeV} < E < 100{\rm\ TeV} $) and ultra-high-energy (UHE, $E > 100\rm\ TeV$) $γ$-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended $γ$-ray source was detected with a significance of ${\sim}16\,σ$ by KM2A and ${\sim}17\,σ$ by WCDA, respectively. The angular extension of this $γ$-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the $γ$-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

preprint2025arXiv

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4.

preprint2022arXiv

Signatures of anisotropic diffusion around PeVatrons in 100 TeV gamma-ray data

The Tibet AS$γ$ collaboration has reported a diffuse $γ$-ray emission signal from the Galactic Plane. We consider the possibility that the diffuse emission from the outer Galactic Plane at the highest energies is produced by cosmic rays spreading from a single supernova-type source either in the Local or Perseus arm of the Milky Way. We show that anisotropic diffusion of multi-PeV cosmic rays along the Galactic magnetic field can produce an extended source spanning ten(s) of degrees on the sky, with a flux-per-unit-solid-angle consistent with Tibet AS$γ$ measurements. Observations of this new type of very extended sources, and measurements of their morphologies, can be used to characterize the anisotropic diffusion of PeV cosmic rays in the Galactic magnetic field, and to constrain the locations and properties of past PeVatrons.

preprint2022arXiv

Status of the K-EUSO Orbital Detector of Ultra-high Energy Cosmic Rays

K-EUSO (KLYPVE-EUSO) is a planned orbital mission aimed at studying ultra-high energy cosmic rays (UHECRs) by detecting fluorescence and Cherenkov light emitted by extensive air showers in the nocturnal atmosphere of Earth in the ultraviolet (UV) range. The observatory is being developed within the JEM-EUSO collaboration and is planned to be deployed on the International Space Station after 2025 and operated for at least two years. The telescope, consisting of $\sim10^{5}$ independent pixels, will allow a spatial resolution of $\sim0.6$ km on the ground, and, from a 400 km altitude, it will achieve a large and full sky exposure to sample the highest energy range of the UHECR spectrum. We provide a comprehensive review of the current status of the development of the K-EUSO experiment, paying special attention to its hardware parts and expected performance. We demonstrate how results of the K-EUSO mission can complement the achievements of the existing ground-based experiments and push forward the intriguing studies of ultra-high energy cosmic rays, as well as bring new knowledge about other phenomena manifesting themselves in the atmosphere in the UV range.

preprint2021arXiv

Possibility of measurement of cosmic ray electron spectrum up to 100 TeV with two-layer water Cherenkov detector array

Measurements of cosmic ray electron+positron spectrum above 10 TeV with ground-based experiments is challenging because of the difficulty of rejection of hadronic extensive air shower background. We study the efficiency of rejection of the hadronic background with water Cherenkov detector array supplemented by muon detection layer. We show that addition of a continuous muon detection layer to the experimental setup allows to achieve a ~ 1e-5 rejection factor for hadronic background at 10 TeV and enables measurement of electron spectrum in 10-100 TeV energy range. We show that measurements of electron spectrum in this energy range do not require a high-altitude experiment and can be done with a sea-level detector.

preprint2021arXiv

Sensitivity of top-of-the-mountain fluorescence telescope system for astrophysical neutrino flux above 10 PeV

Tau neutrinos with energies in the PeV-EeV range produce up-going extensive air showers (UEAS) if they interact underground close enough to the surface of the Earth. This work studies detectability of the UEAS with a system of fluorescence telescopes overlooking dark, low reflectivity, area on the ground up to the distances 20-30 km from mountain top(s). Such system could provide sensitivity sufficient for accumulation of the astrophysical neutrino signal statistics at the rate ten(s) events per year in the energy range beyond 10 PeV, thus allowing to extend the energy frontier of neutrino astronomy into 10-100 PeV range. Comparison of sensitivities of the top-of-the-mountain telescope and IceCube Generation II shows that the two approaches for neutrino detection are complementary, providing comparable performance in adjacent energy bands below and above 10 PeV. Sensitivity of the top-of-the-mountain fluorescence telescope system is also sufficient for the discovery of theoretically predicted cosmogenic neutrino signal.

preprint2021arXiv

The POEMMA (Probe of Extreme Multi-Messenger Astrophysics) Observatory

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to accurately observe ultra-high-energy cosmic rays (UHECRs) and cosmic neutrinos from space with sensitivity over the full celestial sky. POEMMA will observe the extensive air showers (EASs) from UHECRs and UHE neutrinos above 20 EeV via air fluorescence. Additionally, POEMMA will observe the Cherenkov signal from upward-moving EASs induced by Earth-interacting tau neutrinos above 20 PeV. The POEMMA spacecraft are designed to quickly re-orientate to follow up transient neutrino sources and obtain unparalleled neutrino flux sensitivity. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two identical satellites flying in loose formation in 525 km altitude orbits. Each POEMMA instrument incorporates a wide field-of-view (45$^\circ$) Schmidt telescope with over 6 m$^2$ of collecting area. The hybrid focal surface of each telescope includes a fast (1~$μ$s) near-ultraviolet camera for EAS fluorescence observations and an ultrafast (10~ns) optical camera for Cherenkov EAS observations. In a 5-year mission, POEMMA will provide measurements that open new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these otherwise inaccessible energies.

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

LHAASO sensitivity for diffuse gamma-ray signals from the Galaxy

We estimate the sensitivity of LHAASO telescope for the large angular scale diffuse gamma-ray flux in multi-TeV - multi-PeV energy range. We discuss possible sources of the signal in this energy range including the guaranteed flux from cosmic ray interactions in the interstellar medium and possible flux from decaying dark matter. We show that LHAASO will be able to detect the diffuse cosmic ray induced gamma-ray flux up to high Galactic latitude regions thus providing firm identification of the Galactic cosmic ray component of the astrophysical neutrino signal detected by IceCube and clarification of the nature of the knee feature of the cosmic ray spectrum. Comparing the diffuse flux sensitivity with the diffuse gamma-ray flux expected from the dark matter decays, we show LHAASO will be able to detect the gamma-ray signal from dark matter particles of PeV-EeV mass decaying on the time scale up to 3e29 s.

preprint2020arXiv

Mapping large-scale diffuse gamma-ray emission in 10-100 TeV band with Cherenkov telescopes

Measurement of diffuse gamma-ray emission from the Milky Way with Imaging Atmospheric Cherenkov Telescopes (IACT) is difficult because of high level of charged cosmic ray background and small field-of-view. We show that such a measurement is nevertheless possible in the energy band 10-100 TeV. The minimal charged particle background for IACTs is achieved in event selections used for the analyses of the cosmic ray electrons. Tight cuts on the event quality in these event selections allow to achieve the background level sufficiently low to allow the measurements of the diffuse Galactic gamma-ray flux above 10 TeV. We calculate the sensitivities of different types of IACT arrays for the Galactic diffuse emission measurements and compare them with the diffuse gamma-ray flux from different parts of the sky measured by the Fermi Large Area Telescope below 3 TeV and with the astrophysical neutrino signal measured by IceCube telescope. We show that deep exposure of existing IACT systems is sufficient for detection of the diffuse flux from all the Galactic Plane up to Galactic latitude |b|~ degrees. Medium Size Telescope array of CTA will be able to detect the diffuse flux up 30 degree Galactic latitude. Its sensitivity will be sufficient for detection of the gamma-ray counterpart of the Galactic component of the IceCube astrophysical neutrino signal above 10 TeV. We also propose that a dedicated IACT system composed of small but wide field-of-view telescopes could be used for mapping the 10-100 TeV the diffuse gamma-ray emission from across all the sky.

preprint2020arXiv

Snowmass 2021 Letter of Interest: The Probe Of Multi-Messenger Astrophysics (POEMMA)

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to identify the sources of Ultra-High-Energy Cosmic Rays (UHECRs) and to observe cosmic neutrinos, both with full-sky coverage. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two spacecraft flying in a loose formation at 525 km altitude, 28.5 deg inclination orbits. Each spacecraft hosts a Schmidt telescope with a large collecting area and wide field of view. A novel focal plane is optimized to observe both the UV fluorescence signal from extensive air showers (EASs) and the beamed optical Cherenkov signals from EASs. In POEMMA-stereo fluorescence mode, POEMMA will measure the spectrum, composition, and full-sky distribution of the UHECRs above 20 EeV with high statistics along with remarkable sensitivity to UHE neutrinos. The spacecraft are designed to quickly re-orient to a POEMMA-limb mode to observe neutrino emission from Target-of-Opportunity (ToO) transient astrophysical sources viewed just below the Earth's limb. In this mode, POEMMA will have unique sensitivity to cosmic neutrino tau events above 20 PeV by measuring the upward-moving EASs induced by the decay of the emerging tau leptons following the interactions of neutrino tau inside the Earth.

preprint2020arXiv

The INTEGRAL view of the pulsating hard X-ray sky: from accreting and transitional millisecond pulsars to rotation-powered pulsars and magnetars

In the last 25 years, a new generation of X-ray satellites imparted a significant leap forward in our knowledge of X-ray pulsars. The discovery of accreting and transitional millisecond pulsars proved that disk accretion can spin up a neutron star to a very high rotation speed. The detection of MeV-GeV pulsed emission from a few hundreds of rotation-powered pulsars probed particle acceleration in the outer magnetosphere, or even beyond. Also, a population of two dozens of magnetars has emerged. INTEGRAL played a central role to achieve these results by providing instruments with high temporal resolution up to the hard X-ray/soft gamma-ray band and a large field of view imager with good angular resolution to spot hard X-ray transients. In this article, we review the main contributions by INTEGRAL to our understanding of the pulsating hard X-ray sky, such as the discovery and characterization of several accreting and transitional millisecond pulsars, the generation of the first catalog of hard X-ray/soft gamma-ray rotation-powered pulsars, the detection of polarization in the hard X-ray emission from the Crab pulsar, and the discovery of persistent hard X-ray emission from several magnetars.

preprint2019arXiv

A Mini-Imaging Air Cherenkov Telescope

In this paper we describe the different software and hardware elements of a mini-telescope for the detection of cosmic rays and gamma-rays using the Cherenkov light emitted by their induced particle showers in the atmosphere. We estimate the physics reach of the standalone mini-telescope and present some results of the measurements done at the Sauverny Observatory of the University of Geneva and at the Saint-Luc Observatory, which demonstrate the ability of the telescope to observe cosmic rays with energy above about 100 TeV. Such a mini-telescope can constitute a cost-effective out-trigger array that can surround other gamma-ray telescopes or extended air showers detector arrays. Its development was born out of the desire to illustrate to students and amateurs the cosmic ray and gamma-ray detection from ground, as an example of what is done in experiments using larger telescopes. As a matter of fact, a mini-telescope can be used in outreach night events. While outreach is becoming more and more important in the scientific community to raise interest in the general public, the realisation of the mini-telescope is also a powerful way to train students on instrumentation such as photosensors, their associated electronics, acquisition software and data taking. In particular, this mini-telescope uses silicon photomultipliers (SiPM) and the dedicated ASIC, CITIROC.

preprint2019arXiv

Galactic diffuse gamma-ray emission at TeV energy

Measuring the diffuse Galactic gamma-ray flux in the TeV range is difficult for ground-based gamma-ray telescopes because of the residual cosmic-ray background, which is higher than the gamma-ray flux by several orders of magnitude. Its detection is also challenging for space-based telescopes because of low signal statistics. We characterize the diffuse TeV flux from the Galaxy using decade-long exposures of the Fermi Large Area Telescope. Considering that the level of diffuse Galactic emission in the TeV band approaches the level of residual cosmic-ray background, we estimated the level of residual cosmic-ray background in the SOURCEVETO event selection and verified that the TeV diffuse Galactic emission flux is well above the residual cosmic-ray background up to high Galactic latitude regions. We study spectral and imaging properties of the diffuse TeV signal from the Galactic plane. We find much stronger emission from the inner Galactic plane than in previous HESS telescope estimates (lower bound). We also find a significant difference in the measurement of the Galactic longitude and latitude profiles of the signal measured by Fermi and HESS. These discrepancies are presumably explained by the fact that regions of background estimate in HESS have non-negligible gamma-ray flux. Comparing Fermi measurements with those of ARGO-YBJ, we find better agreement, with the notable exception of the Cygnus region, where we find much higher flux (by a factor 1.5). We also measure the TeV diffuse emission spectrum up to high Galactic latitude and show that the spectra of different regions of the sky have spectral slopes consistent with Gamma=2.34+/- 0.04. We discuss the possible origin of the hard slope of the TeV diffuse emission. Fermi/LAT provides reliable measurements of the diffuse Galactic emission spectrum in the TeV range.

preprint2019arXiv

Search for decaying eV-mass axion-like particles using gamma-ray signal from blazars

Decaying axion-like particles (ALP) with masses in the eV range which might occupy dark matter halos of the Milky Way and other galaxies produce a characteristic "bump" feature in the spectrum of extragalactic background light (EBL). This feature leaves an imprint on the gamma-ray spectra of distant extragalactic sources. We derive constraints on the ALP coupling to photons based on analysis of spectra of very-high-energy gamma-ray loud blazars. We combine gamma-ray spectral measurements by Fermi/LAT and Cherenkov telescopes and fit a model in which the intrinsic source spectrum is modified by pair production on photons produced by ALP decays. We constrain the amplitude of gamma-ray flux suppression by this effect. We find that the combined Fermi/LAT and VERITAS data set for the source 1ES 1218+304 currently provides the tightest constraint on ALP-two-photon coupling which is complementary to the constraints imposed by non-observation of excess energy loss in Horizontal Branch stars, by the high-resolution spectroscopic observations of galaxy clusters with optical telescopes and by the searches of ALP signal with CERN Solar Axion Telescope. Our analysis favours existence of a bump in the EBL spectrum which could be produced by ALPs in the mass range 2-3 eV and axion-photon coupling $\sim 10^{-10}$ GeV$^{-1}$. We discuss possibilities for verification of this hint with deeper Cherenkov telescope observations of large number of blazars with current generation instruments and with the Cherenkov Telescope Array.

preprint2019arXiv

Search for spectral features in extragalactic background light with gamma-ray telescopes

Cumulative optical and infrared emission from galaxies accumulated over cosmological time scales, the extragalactic background light (EBL), could be probed by complementary techniques of direct observations and source counting in the visible and infrared as well as via its imprint on the signal of distant active galactic nuclei in gamma rays. We compare the visible and infrared measurements with the gamma-ray constraints and study if the discrepancies of the measurements with different methods could be due to the presence of features in the EBL spectrum that are localised in the micron wavelength range. We combined data on time-averaged spectra of selected blazars that were obtained by Fermi and ground-based gamma-ray telescopes. We also modelled the effect of absorption on EBL while allowing for the existence of a previously unaccounted spectral feature. We show that a previously reported 'excess' in EBL flux in the $\sim 1$~micron wavelength range is consistent with gamma-ray measurements, that is, if the excess has the form of a narrow feature of the width $δλ< λ$ and an overall flux of up to 15 nW/(m$^2$ sr) above the 'minimal' EBL, which is estimated from the visible and infrared source counts. Such 'bump-like' spectral features could originate, for example, from decaying dark-matter particles, or either axions or peculiar astrophysical processes in the course of star-formation history. We discuss the possibilities for the search of spectral features in the EBL with the Cherenkov Telescope Array (CTA).

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.

preprint2016arXiv

Abundances of 59Co and 59Ni in the cosmic ray flux

Two main hypotheses for the origin of Galactic cosmic rays are the "supernova" and "superbubble" origin hypotheses. We analyse the evidence for the superbubble hypothesis provided by the measurements of the relativive abundances of isotopes of cobalt and nickel in the cosmic ray flux. We compare the measured upper limit on the abundance of 59Ni in the cosmic ray flux with the 59Ni abundance predictions of the up-to-date stellar evolution models. Non-detection of 59Ni in the cosmic ray flux has previously been attributed to a large time delay of the order of 1e5 yr between the moment of supernova explosion and the onset of particle acceleration process. This large time delay was considered as an argument in favour of the "superbubble" scenario. We show that the recent calculation of the 59Ni yield of massive stars, which takes into account the initial mass range up to 120 solar masses and includes stellar rotation, results in prediction of low 59Ni abundance relative to its decay product 59Co. The predicted abundance is consistent with the upper bound on 59Ni abundance in the cosmic ray flux for the supernova parameters assumed. This result removes the necessity of decay of 59Ni in the time interval between the supernova explosion and the onset of acceleration process and restores the consistency of measurements of 59Ni / 59Co abundances with the "supernova" hypothesis of the CR origin.

preprint2016arXiv

Bokeh Mirror Alignment for Cherenkov Telescopes

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures and high image intensities to map the faint Cherenkov light emitted from cosmic ray air showers onto their image sensors. Segmented reflectors fulfill these needs, and composed from mass production mirror facets they are inexpensive and lightweight. However, as the overall image is a superposition of the individual facet images, alignment remains a challenge. Here we present a simple, yet extendable method, to align a segmented reflector using its Bokeh. Bokeh alignment does not need a star or good weather nights but can be done even during daytime. Bokeh alignment optimizes the facet orientations by comparing the segmented reflectors Bokeh to a predefined template. The optimal Bokeh template is highly constricted by the reflector's aperture and is easy accessible. The Bokeh is observed using the out of focus image of a near by point like light source in a distance of about 10 focal lengths. We introduce Bokeh alignment on segmented reflectors and demonstrate it on the First Geiger-mode Avalanche Cherenkov Telescope (FACT) on La Palma, Spain.

preprint2016arXiv

Galactic and extragalactic contributions to the astrophysical muon neutrino signal

Spectral and anisotropy properties of IceCube astrophysical neutrino signal reveal an evidence for a significant Galactic contribution to the neutrino flux in Southern hemisphere. We check if the Galactic contribution is detectable in the astrophysical muon neutrino flux observed from a low positive declinations region of the Northern sky. Estimating the Galactic neutrino flux in this part of the sky from gamma-ray and Southern sky neutrino data, we find that the Northern sky astrophysical muon neutrino signal shows an excess over the Galactic flux. This points to the presence of an additional hard spectrum (extragalactic or large scale Galactic halo) component of astrophysical neutrino flux. We show that the Galactic flux component should still be detectable in the muon neutrino data in a decade long IceCube exposure.

preprint2016arXiv

Hard X-ray view of microlensing events in RX J1131-1231

RX J1131-1231 is a gravitationally lensed system which includes four images of a quasar lensed by an elliptical galaxy. The flux in the individual images is known to be affected by microlensing effect in the visible and X-ray bands. We study the multi-wavelength properties of RX J1131-1231 over a broad energy range, from optical to hard X-ray, during the periods of the microlensing caustic crossings. We aim to constrain the spatial extent of the X-ray emission region at different energies. We combine the data of the source monitoring in the visible band with the X-ray data of the Burst Alert Telescope (BAT) on board of SWIFT satellite and Chandra X-ray observatory. Inspecting the broad band spectrum and lightcurves of the source we identify several microlensing caustic crossing events, and study the details of variability of the source during these events. The caustic crossings of image A on MJD 55150 and 55500 produce strong variations of the overall X-ray flux from the source. In the soft X-ray band, the caustic crossing events are characterised by the exponential growth / decay of the flux ratio between the source images with the exponentiating time scale T_d~60 d, which corresponds to the caustic transit time of the source of the size ~10^{14} cm. Variability induced by the microlensing dominates the overall variability of the sum of the fluxes of four images of the source which should also be detectable in the hard X-ray lightcurve, if the hard X-ray source is comparably compact. Non-detection of the caustic crossing events in the SWIFT/BAT lightcurve of the source indicates that geometry of the region of emission of X-rays with energies above 15 keV is different from that of the lower energy X-rays, with the hard X-ray emission an order of magnitude larger in size.

preprint2016arXiv

Normalized and Asynchronous Mirror Alignment for Cherenkov Telescopes

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures and high image intensities to map the faint Cherenkov light emitted from cosmic ray air showers onto their image sensors. Segmented reflectors fulfill these needs, and as they are composed from mass production mirror facets they are inexpensive and lightweight. However, as the overall image is a superposition of the individual facet images, alignment is a challenge. Here we present a computer vision based star tracking alignment method, which also works for limited or changing star light visibility. Our method normalizes the mirror facet reflection intensities to become independent of the reference star's intensity or the cloud coverage. Using two CCD cameras, our method records the mirror facet orientations asynchronously of the telescope drive system, and thus makes the method easy to integrate into existing telescopes. It can be combined with remote facet actuation, but does not require one to work. Furthermore, it can reconstruct all individual mirror facet point spread functions without moving any mirror. We present alignment results on the 4 meter First Geiger-mode Avalanche Cherenkov Telescope (FACT).

preprint2016arXiv

Particle acceleration in the vacuum gaps in black hole magnetospheres

We consider particle acceleration in vacuum gaps in magnetospheres of black holes powered through Blandford-Znajek mechanism and embedded into radiatively-inefficient accretion flow (RIAF) environment. In such situation the gap height is limited by the onset of gamma-gamma pair production on the infrared photons originating from the RIAF. We numerically calculate acceleration and propagation of charged particles taking into account the detailed structure of electric and magnetic field in the gap and in the entire black hole magnetosphere, radiative energy losses and interactions of gamma rays produced by the propagated charged particles with the background radiation field of RIAF. We show that the presence of the vacuum gap has clear observational signatures. The spectra of emission from gaps embedded into a relatively high luminosity RIAF are dominated by the inverse Compton emission with a sharp, super-exponential cut-off in the very-high-energy gamma-ray band. The cut-off energy is determined by the properties of the RIAF and is largely independent of the structure of magnetosphere and geometry of the gap.The spectra of the gap residing in low-luminosity RIAFs are dominated by synchrotron / curvature emission with the spectra extending into 1-100 GeV energy range. We also consider the effect of possible acceleration of protons in the gap and find that proton energies could reach the ultra-high-energy cosmic ray (UHECR) range only in extremely low luminosity RIAFs.

preprint2016arXiv

Toward the full test of the nuMSM sterile neutrino dark matter model with Athena

We discuss the potential of Athena X-ray telescope, in particular of its X-ray Integral Field Unit (X-IFU), for detection of the signal from the light-weight decaying dark matter with mass in the keV range. We show that high energy resolution and large collection area of X-IFU will provide an improvement of sensitivity which will be sufficient for the full test of the neutrino Minimal extension of the Standard Model (nuMSM). Search for the narrow spectral line produced by the decay of the dark matter sterile neutrino in the spectra of dwarf spheroidal galaxies with X-IFU will explore the whole allowed range masses and mixing angles of the nuMSM lightest sterile neutrino and in this way either to find the dark matter signal or rule out the nuMSM model.

preprint2015arXiv

Constraints on 3.55 keV line emission from stacked observations of dwarf spheroidal galaxies

Several recent works have reported the detection of an unidentified X-ray line at 3.55~keV, which could possibly be attributed to the decay of dark matter (DM) particles in the halos of galaxy clusters and in the M31 galaxy. We analyze all publicly-available XMM-Newton data of dwarf spheroidal galaxies to test the possible DM origin of the line. Dwarf spheroidal galaxies have high mass-to-light ratios and their interstellar medium is not a source of diffuse X-ray emission, thus they are expected to provide the cleanest DM decay line signal. Our analysis shows no evidence for the presence of the line in the stacked spectra of the dwarf galaxies. It excludes the sterile neutrino DM decay origin of the 3.5 keV line reported by Bulbul et al. (2014) at the level of 4.1sigma under standard assumptions about the Galactic DM column density in the direction of selected dwarf galaxies and at the level of 3.2sigma assuming minimal Galactic DM column density. Our analysis is still consistent with the estimate of sterile neutrino DM parameters by Boyarsky et al. (2014), because of its larger uncertainty. However, the central value of their estimate of the mixing angle is inconsistent with our dwarf spheroidals data at 3.4sigma (2.5sigma) level assuming the mean (minimal) Galactic DM column density. As a by-product of our analysis, we provide updated upper limits to the mixing angle of sterile neutrino DM in the mass range between 2 and 20 keV.

preprint2015arXiv

Data compression for the First G-APD Cherenkov Telescope

The First Geiger-mode Avalanche photodiode (G-APD) Cherenkov Telescope (FACT) has been operating on the Canary island of La Palma since October 2011. Operations were automated so that the system can be operated remotely. Manual interaction is required only when the observation schedule is modified due to weather conditions or in case of unexpected events such as a mechanical failure. Automatic operations enabled high data taking efficiency, which resulted in up to two terabytes of FITS files being recorded nightly and transferred from La Palma to the FACT archive at ISDC in Switzerland. Since long term storage of hundreds of terabytes of observations data is costly, data compression is mandatory. This paper discusses the design choices that were made to increase the compression ratio and speed of writing of the data with respect to existing compression algorithms. Following a more detailed motivation, the FACT compression algorithm along with the associated I/O layer is discussed. Eventually, the performances of the algorithm is compared to other approaches.

preprint2015arXiv

DigiCam - Fully Digital Compact Read-out and Trigger Electronics for the SST-1M Telescope proposed for the Cherenkov Telescope Array

The SST-1M is one of three prototype small-sized telescope designs proposed for the Cherenkov Telescope Array, and is built by a consortium of Polish and Swiss institutions. The SST-1M will operate with DigiCam - an innovative, compact camera with fully digital read-out and trigger electronics. A high level of integration will be achieved by massively deploying state-of-the-art multi-gigabit transmission channels, beginning from the ADC flash converters, through the internal data and trigger signals transmission over backplanes and cables, to the camera's server link. Such an approach makes it possible to design the camera to fit the size and weight requirements of the SST-1M exactly, and provide low power consumption, high reliability and long lifetime. The structure of the digital electronics will be presented, along with main physical building blocks and the internal architecture of FPGA functional subsystems.

preprint2015arXiv

FACT - Monitoring Blazars at Very High Energies

The First G-APD Cherenkov Telescope (FACT) was built on the Canary Island of La Palma in October 2011 as a proof of principle for silicon based photosensors in Cherenkov Astronomy. The scientific goal of the project is to study the variability of active galatic nuclei (AGN) at TeV energies. Observing a small sample of TeV blazars whenever possible, an unbiased data sample is collected. This allows to study the variability of the selected objects on timescales from hours to years. Results from the first three years of monitoring will be presented. To provide quick flare alerts to the community and trigger multi-wavelength observations, a quick look analysis has been installed on-site providing results publicly online within the same night. In summer 2014, several flare alerts were issued. Results of the quick look analysis are summarized.

preprint2015arXiv

Hard spectrum of cosmic rays in the Disks of Milky Way and Large Magellanic Cloud

The slope of the locally measured spectrum of cosmic rays varies from 2.8 for protons with energies below 200 GeV down to 2.5 for heavy nuclei with energies in the TeV-PeV range. It is not clear if the locally measured slope values are representative for those of the overall population of Galactic cosmic rays and if the slope of the cosmic ray spectrum varies across the Galaxy. We use the data of Fermi Space gamma-ray Telescope to derive a measurement of the slope of the cosmic ray spectrum across the Galactic Disk and to compare it with that of the Large Magellanic Cloud cosmic rays. A special choice of the background estimation regions allows us to single out the neutral pion decay component of the gamma-ray flux in the energy range above 10 GeV and to separate it from (a) emission from the local interstellar medium around the Solar system and (b) from the inverse Compton emission produced by cosmic ray electrons. Results. The spectrum of the pion decay gamma-ray emission from the Galactic disk in the energy band 10 GeV - 1 TeV has the slope 2.4. There is no evidence for the variation of the slope with Galactic longitude / distance from the Galactic Centre. The slope of the spectrum of cosmic rays derived from the gamma-ray data, 2.45, is harder than the slope of the locally observed cosmic ray proton spectrum. Pion decay emission from a powerlaw distribution of cosmic rays with the same hard slope also provides a fit to the gamma-ray spectrum of the Large Magellanic Cloud. Identical and hard slopes of the spectra of cosmic rays in the Milky Way and in the Large Magellanic Cloud are consistent with a straightforward theoretical model in which cosmic rays are injected by shock acceleration with the spectrum with the slope 2...2.1 which is subsequently modified by 1/3...1/2 by the energy-dependent escape of cosmic rays through the turbulent Galactic magnetic field.

preprint2015arXiv

Leptonic origin of the 100 MeV gamma-ray emission from the Galactic Centre

The Galactic centre is a bright gamma-ray source with the GeV-TeV band spectrum composed of two distinct components in the 1-10 GeV and 1-10 TeV energy ranges. The nature of these two components is not clearly understood. We investigate the gamma-ray properties of the Galactic centre to clarify the origin of the observed emission. We report imaging, spectral, and timing analysis of data from 74 months of observations of the Galactic centre by FERMI/LAT gamma-ray telescope complemented by sub-MeV data from approximately ten years of INTEGRAL/PICsIT observations. We find that the Galactic centre is spatially consistent with the point source in the GeV band. The tightest 3 sigma upper limit on its radius is 0.13 degree in the 10-300 GeV energy band. The spectrum of the source in the 100 MeV energy range does not have a characteristic turnover that would point to the pion decay origin of the signal. Instead, the source spectrum is consistent with a model of inverse Compton scattering by high-energy electrons. In this a model, the GeV bump in the spectrum originates from an episode of injection of high-energy particles, which happened ~300 years ago. This injection episode coincides with the known activity episode of the Galactic centre region, previously identified using X-ray observations. The hadronic model of source activity could be still compatible with the data if bremsstrahlung emission from high-energy electrons was present in addition to pion decay emission.

preprint2015arXiv

Multi-wavelength Observations of the Binary System PSR B1259-63/LS~2883 around the 2014 Periastron Passage

We report on the results of the extensive multi-wavelength campaign from optical to GeV gamma-rays of the 2014 periastron passage of PSR B1259-63, which is a unique high-mass gamma-ray emitting binary system with a young pulsar companion. Observations demonstrate the stable nature of the post-periastron GeV flare and prove the coincidence of the flare with the start of rapid decay of the H$α$ equivalent width, usually interpreted as a disruption of the Be stellar disk. Intensive X-ray observations reveal changes in the X-ray spectral behaviour happening at the moment of the GeV flare. We demonstrate that these changes can be naturally explained as a result of synchrotron cooling of monoenergetic relativistic electrons injected into the system during the GeV flare.

preprint2015arXiv

Signatures of a two million year old supernova in the spectra of cosmic ray protons, antiprotons and positrons

The locally observed cosmic ray spectrum has several puzzling features, such as the excess of positrons and antiprotons above $\sim 20$ GeV and the discrepancy in the slopes of the spectra of cosmic ray protons and heavier nuclei in the TeV-PeV energy range. We show that these features are consistently explained by a nearby source which was active $\sim 2$ Myr ago and has injected $(1-2)\times 10^{50}$ erg in cosmic rays. The transient nature of the source and its overall energy budget point to the supernova origin of this local cosmic ray source. The age of the supernova suggests that the local cosmic ray injection was produced by the same supernova that has deposited $^{60}$Fe isotopes in the deep ocean crust.

preprint2015arXiv

Test of relativistic gravity using microlensing of relativistically broadened lines in gravitationally lensed quasars

We show that observation of the time-dependent effect of microlensing of relativistically broadened emission lines (such as e.g. the Fe Kalpha line in X-rays) in strongly lensed quasars could provide data on celestial mechanics of circular orbits in the direct vicinity of the horizon of supermassive black holes. This information can be extracted from the observation of evolution of red / blue edge of the magnified line just before and just after the period of crossing of the innermost stable circular orbit by the microlensing caustic. The functional form of this evolution is insensitive to numerous astrophysical parameters of the accreting black hole and of the microlensing caustics network system (as opposed to the evolution the full line spectrum). Measurement of the temporal evolution of the red / blue edge could provide a precision measurement of the radial dependence of the gravitational redshift and of velocity of the circular orbits, down to the innermost stable circular orbit. These measurements could be used to discriminate between the General Relativity and alternative models of the relativistic gravity in which the dynamics of photons and massive bodies orbiting the gravitating centre is different from that of the geodesics in the Schwarzschild or Kerr space-times.

preprint2015arXiv

The Atmospheric Monitoring System of the JEM-EUSO Space Mission

An Atmospheric Monitoring System (AMS) is a mandatory and key device of a space-based mission which aims to detect Ultra-High Energy Cosmic Rays (UHECR) and Extremely-High Energy Cosmic Rays (EHECR) from Space. JEM-EUSO has a dedicated atmospheric monitoring system that plays a fundamental role in our understanding of the atmospheric conditions in the Field of View (FoV) of the telescope. Our AMS consists of a very challenging space infrared camera and a LIDAR device, that are being fully designed with space qualification to fulfil the scientific requirements of this space mission. The AMS will provide information of the cloud cover in the FoV of JEM-EUSO, as well as measurements of the cloud top altitudes with an accuracy of 500 m and the optical depth profile of the atmosphere transmittance in the direction of each air shower with an accuracy of 0.15 degree and a resolution of 500 m. This will ensure that the energy of the primary UHECR and the depth of maximum development of the EAS ( Extensive Air Shower) are measured with an accuracy better than 30\% primary energy and 120 $g/cm^2$ depth of maximum development for EAS occurring either in clear sky or with the EAS depth of maximum development above optically thick cloud layers. Moreover a very novel radiometric retrieval technique considering the LIDAR shots as calibration points, that seems to be the most promising retrieval algorithm is under development to infer the Cloud Top Height (CTH) of all kind of clouds, thick and thin clouds in the FoV of the JEM-EUSO space telescope.

preprint2015arXiv

The SST-1M camera for the Cherenkov Telescope Array

The prototype camera of the single-mirror Small Size Telescopes (SST-1M) proposed for the Cherenkov Telescope Array (CTA) project has been designed to be very compact and to deliver high performance over thirty years of operation. The camera is composed of an hexagonal photo-detection plane made of custom designed large area hexagonal silicon photomultipliers and a high throughput, highly configurable, fully digital readout and trigger system (DigiCam). The camera will be installed on the telescope structure at the H. Niewodnicza{ń}ski institute of Nuclear Physics in Krakow in fall 2015. In this contribution, we review the steps that led to the development of the innovative photo-detection plane and readout electronics, and we describe the test and calibration strategy adopted.

preprint2015arXiv

Unified model for cosmic rays above $10^{17}$ eV and the diffuse gamma-ray and neutrino backgrounds

We investigate how the extragalactic proton component derived within the "escape model" can be explained by astrophysical sources. We consider as possible cosmic ray (CR) sources normal/starburst galaxies and radio-loud active galactic nuclei (AGN). We find that the contribution to the total extragalactic proton flux from normal and starburst galaxies is only subdominant and does not fit the spectral shape deduced in the escape model. In the case of radio-loud AGN, we show that the complete extragalactic proton spectrum can be explained by a single source population, BL Lac/FR I, for any of the potential acceleration sites in these sources. We calculate the diffuse neutrino and $γ$-ray fluxes produced by these CR protons interacting with gas inside their sources. For a spectral slope of CRs close to $α=2.1-2.2$ as suggested by shock acceleration, we find that these UHECR sources contribute the dominant fraction of both the isotropic $γ$-ray background and of the extragalactic part of the astrophysical neutrino signal observed by IceCube.

preprint2014arXiv

AtmoHEAD 2013 workshop / Atmospheric Monitoring for High-Energy Astroparticle Detectors

A 3-day international workshop on atmospheric monitoring and calibration for high-energy astroparticle detectors, with a view towards next-generation facilities. The atmosphere is an integral component of many high-energy astroparticle detectors. Imaging atmospheric Cherenkov telescopes and cosmic-ray extensive air shower detectors are the two instruments driving the rapidly evolving fields of very-high- and ultra-high-energy astrophysics. In these instruments, the atmosphere is used as a giant calorimeter where cosmic rays and gamma rays deposit their energy and initiate EASs; it is also the medium through which the resulting Cherenkov light propagates. Uncertainties in real-time atmospheric conditions and in the fixed atmospheric models typically dominate all other systematic errors. With the improved sensitivity of upgraded IACTs such as H.E.S.S.-II and MAGIC-II and future facilities like the Cherenkov Telescope Array (CTA) and JEM-EUSO, statistical uncertainties are expected to be significantly reduced, leaving the atmosphere as the limiting factor in the determination of astroparticle spectra. Varying weather conditions necessitate the development of suitable atmospheric monitoring to be integrated in the overall instrument calibration, including Monte Carlo simulations. With expertise distributed across multiple collaborations and scientific domains, an interdisciplinary workshop is being convened to advance progress on this critical and timely topic.

preprint2014arXiv

Multi-wavelength observations of the binary system PSR B1259-63/LS 2883 around the 2010-2011 periastron passage

We report on broad multi-wavelength observations of the 2010-2011 periastron passage of the gamma-ray loud binary system PSR B1259-63. High resolution interferometric radio observations establish extended radio emission trailing the position of the pulsar. Observations with the Fermi Gamma-ray Space Telescope reveal GeV gamma-ray flaring activity of the system, reaching the spin-down luminosity of the pulsar, around 30 days after periastron. There are no clear signatures of variability at radio, X-ray and TeV energies at the time of the GeV flare. Variability around periastron in the H$α$ emission line, can be interpreted as the gravitational interaction between the pulsar and the circumstellar disk. The equivalent width of the H$α$ grows from a few days before periastron until a few days later, and decreases again between 18 and 46 days after periastron. In near infrared we observe the similar decrease of the equivalent width of Br$γ$ line between the 40th and 117th day after the periastron. For the idealized disk, the variability of the H$α$ line represents the variability of the mass and size of the disk. We discuss possible physical relations between the state of the disk and GeV emission under assumption that GeV flare is directly related to the decrease of the disk size.

preprint2014arXiv

PeV neutrinos from interactions of cosmic rays with the interstellar medium in the Galaxy

We present a self-consistent interpretation of the of very-high-energy neutrino signal from the direction of the inner Galaxy, which is a part of the astronomical neutrino signal reported by IceCube. We demonstrate that an estimate of the neutrino flux in the E>100 TeV energy range lies at the high-energy power-law extrapolation of the spectrum of diffuse gamma-ray emission from the Galactic Ridge, as observed by Fermi telescope. This suggests that IceCube neutrino and Fermi/LAT gamma-ray fluxes are both produced in interactions of cosmic rays with the interstellar medium in the Norma arm and/or in the Galactic Bar. Cosmic rays responsible for the gamma-ray and neutrino flux are characterised by hard spectrum with the slope harder than -2.4 and cut-off energy higher than 10 PeV.

preprint2014arXiv

Remote sensing of clouds and aerosols with cosmic rays

Remote sensing of atmosphere is conventionally done via a study of extinction / scattering of light from natural (Sun, Moon) or artificial (laser) sources. Cherenkov emission from extensive air showers generated by cosmic rays provides one more natural light source distributed throughout the atmosphere. We show that Cherenkov light carries information on three-dimensional distribution of clouds and aerosols in the atmosphere and on the size distribution and scattering phase function of cloud/aerosol particles. Therefore, it could be used for the atmospheric sounding. The new atmospheric sounding method could be implemented via an adjustment of technique of imaging Cherenkov telescopes. The atmospheric sounding data collected in this way could be used both for atmospheric science and for the improvement of the quality of astronomical gamma-ray observations.

preprint2013arXiv

An evaluation of the exposure in nadir observation of the JEM-EUSO mission

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.

preprint2013arXiv

An exploration of hadronic interactions in blazars using IceCube

Context: Hadronic models, involving matter (proton or nuclei) acceleration in blazar jets, imply high energy photon and neutrino emissions due to interactions of high-energy protons with matter and/or radiation in the source environment. Aims: This paper shows that the sensitivity of the IceCube neutrino telescope in its 40-string configuration (IC-40) is already at the level of constraining the parameter space of purely hadronic scenarios of activity of blazars. Methods: Assuming that the entire source power originates from hadronic interactions, and assuming that the models describe the data, we estimate the expected neutrino flux from blazars based on the observed gamma-ray flux by Fermi, simultaneously with IC-40 observations. We consider two cases separately to keep the number of constrainable parameters at an acceptable level: proton-proton or proton-gamma interactions are dominant. Comparing the IC-40 sensitivity to the neutrino flux expected from some of the brightest blazars, we constrain model parameters characterizing the parent high-energy proton spectrum. Results: We find that when pp interactions dominate, some constraints on the primary proton spectrum can be imposed. For instance, for the tightest constrained source 3C 454.3, the very high energy part of the spectra of blazars is constrained to be harder than E^-2 with cut-off energies in the range of Ecut >10^18 eV. When interactions of high-energy protons on soft photon fields dominate, we can find similarly tight constraints on the proton spectrum parameters. [abridged]

preprint2013arXiv

Cosmological Magnetic Fields: Their Generation, Evolution and Observation

We review the possible mechanisms for the generation of cosmological magnetic fields, discuss their evolution in an expanding Universe filled with the cosmic plasma and provide a critical review of the literature on the subject. We put special emphasis on the prospects for observational tests of the proposed cosmological magnetogenesis scenarios using radio and gamma-ray astronomy and ultra high energy cosmic rays. We argue that primordial magnetic fields are observationally testable. They lead to magnetic fields in the intergalactic medium with magnetic field strength and correlation length in a well defined range. We also state the unsolved questions in this fascinating open problem of cosmology and propose future observations to address them.

preprint2013arXiv

Magnetic Fields in the Large Scale Structure from Faraday Rotation measurements

We search for magnetic fields outside galaxies and galaxy clusters by investigating redshift evolution of Faraday rotation measures (RM) of extragalactic radio sources. Our analysis reveals a strong evidence for the redshift dependence of the mean of the absolute value of the RM. The evidence is further strengthened if the Galactic contribution to the RM is subtracted. The hypothesis of the absence of the redshift evolution of residual RM is ruled out at 5 sigma level. The observed redshift dependence of RM is consistent with the possibility of the presence of nano-Gauss strength magnetic fields with correlation length shorter than 0.1 Mpc in the weakly overdense elements of the Large Scale Structure traced by the Ly-alpha clouds.

preprint2013arXiv

Measuring the correlation length of intergalactic magnetic fields from observations of gamma-ray induced cascades

Context. The imaging and timing properties of γ-ray emission from electromagnetic cascades initiated by very-high-energy (VHE) γ-rays in the intergalactic medium depend on the strength B and correlation length λB of intergalactic magnetic fields (IGMF). Aims. We study the possibility of measuring both B and λB via observations of the cascade emission with γ-ray telescopes. Methods. For each measurement method, we find two characteristics of the cascade signal, which are sensitive to the IGMF B and λB values in different combinations. For the case of IGMF measurement using the observation of extended emission around extragalactic VHE γ-ray sources, the two characteristics are the slope of the surface brightness profile and the overall size of the cascade source. For the case of IGMF measurement from the time delayed emission, these two characteristics are the initial slope of the cascade emission light curve and the overall duration of the cascade signal. Results. We show that measurement of the slope of the cascade induced extended emission and/or light curve can both potentially provide measure of the IGMF correlation length, provided it lies within the range 10 kpc< λB <1 Mpc. For correlation lengths outside this range, gamma-ray observations can provide upper or lower bound on λB. The latter of the two methods holds great promise in the near future for providing a measurement/constraint using measurements from present/next-generation γ-ray-telescopes. Conclusions. Measurement of the IGMF correlation length will provide an important constraint on its origin. In particular, it will enable to distinguish between an IGMF of galactic wind origin from an IGMF of cosmological origin.

preprint2013arXiv

Neutrino signal from extended Galactic sources in IceCube

We explore the detectability of the neutrino flux from the entire Galactic Plane or from a part of it with IceCube. We calculate the normalization and the spectral index of the neutrino power law spectrum from different regions of the Galactic plane, based on the observed spectral characteristics of the pion decay gamma-ray diffuse emission observed by the Fermi/LAT telescope in the energy band above 100 GeV. We compare the neutrino flux calculated in this way with the sensitivity of IceCube for the detection of extended sources. Assuming a binned extended source analysis method, we find that the only possible evidence for neutrino emission for sources located in the Northern hemisphere is from the Cygnus region after 20 years of exposure. For other parts of the Galactic Plane even a 20 years exposure with IceCube is not sufficient for the detection. Taking into account marginal significance of the detectable source in the Cygnus region, we find a precise position and size of the source region which optimizes the signal-to-noise ratio for neutrinos. We also calculate the low-energy threshold above which the neutrino signal could be detected with the highest signal-to-noise ratio. This calculation of precise source position, size and energy range, based on the gamma-ray data, could be used to remove the 'trial factor' in the analysis of the real neutrino data of IceCube. We notice that the diffuse neutrino emission from the inner Galactic Plane in the Southern Hemisphere is much brighter. A neutrino detector with characteristics equivalent to IceCube, but placed at the Northern Hemisphere (such as KM3NeT), would detect several isolated neutrino sources in the Galactic Plane within just 5 years exposure at 5σ level. These isolated sources of ~TeV neutrinos would unambiguously localize sources of cosmic rays which operated over the last 10 thousand years in the Galaxy.[abridged]

preprint2013arXiv

Potential of LOFT telescope for the search of dark matter

Large Observatory For X-ray Timing (LOFT) is a next generation X-ray telescope selected by European Space Agency as one of the space mission concepts within the ``Cosmic Vision'' programme. The Large Area Detector on board of LOFT will be a collimator-type telescope with an unprecedentedly large collecting area of about 10 square meters in the energy band between 2 and 100 keV. We demonstrate that LOFT will be a powerful dark matter detector, suitable for the search of the X-ray line emission expected from decays of light dark matter particles in galactic halos. We show that LOFT will have sensitivity for dark matter line search more than an order of magnitude higher than that of all existing X-ray telescopes. In this way, LOFT will be able to provide a new insight into the fundamental problem of the nature of dark matter.

preprint2012arXiv

Degree-scale GeV "jets" from active and dead TeV blazars

We show that images of TeV blazars in the GeV energy band should contain, along with point-like sources, degree-scale jet-like extensions. These GeV extensions are the result of electromagnetic cascades initiated by TeV gamma-rays interacting with extragalactic background light and the deflection of the cascade electrons/positrons in extragalactic magnetic fields (EGMF). Using Monte-Carlo simulations, we study the spectral and timing properties of the degree-scale extensions in simulated GeV band images of TeV blazars. We show that the brightness profile of such degree-scale extensions can be used to infer the lightcurve of the primary TeV gamma-ray source over the past 1e7 yr, i.e. over a time scale comparable to the life-time of the parent active galactic nucleus. This implies that the degree-scale jet-like GeV emission could be detected not only near known active TeV blazars, but also from "TeV blazar remnants", whose central engines were switched off up to ten million years ago. Since the brightness profile of the GeV "jets" depends on the strength and the structure of the EGMF, their observation provides additionally information about the EGMF.

preprint2012arXiv

Extragalactic Very-High-Energy gamma-ray background

We study the origin of the extragalactic diffuse gamma-ray background using the data from the Fermi telescope. To estimate the background level, we count photons at high Galactic latitudes |b|>60 degrees. Subtracting photons associated to known sources and the residual cosmic ray and Galactic diffuse backgrounds, we estimate the Extragalactic Gamma-ray Background (EGB) flux. We find that the spectrum of EGB in the very-high-energy (VHE) band above 30 GeV follows the stacked spectrum of BL Lacs. LAT data reveal the positive (1+z)^k, 1<k<4 cosmological evolution of the BL Lac source population consistent with that of their parent population, FR I radio galaxies. We show that EGB at E>30 GeV could be completely explained by emission from unresolved BL Lacs if k~3.

preprint2012arXiv

High Galactic latitude Fermi sources of gamma-rays with energies above 100 GeV

We present a catalog of sources of very high energy (E>100 GeV) gamma-rays detected by the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope at Galactic latitudes |b|>=10 degrees. We cross-correlate the directions of individual photons with energies above 100 GeV detected by Fermi/LAT with the first year Fermi catalog of sources detected at lower energies. We find a significant correlation between the arrival directions of the highest energy photons and positions of Fermi sources, with the possibility of chance coincidences at the level of 1e-45. We present a list of Fermi sources contributing to the correlation signal. A similar analysis is performed to cross-correlate the catalog of BL Lac objects with the highest energy photons detected by Fermi. We produce a catalog of high Galactic latitude Fermi sources visible at 100-300 GeV energies. The catalog is divided into two parts. The first part contains a list of 50 higher significance sources among which there can be three possible false detections. The second part of the catalog contains a list of 25 lower significance sources, among which six are possibly false detections. Finally we identify eight additional sources from the cross-correlation analysis with the BL Lac catalog. The reported sources of E>100 GeV gamma-rays span a broad range of redshifts, up to z~1. Most of the sources are BL Lac type objects. Only 17 out of 83 objects in our list were previously reported as VHE gamma-ray sources.

preprint2012arXiv

Origin of TeV Galactic Cosmic Rays

We consider a possibility of identification of sources of cosmic rays (CR) of the energy above 1 TeV via observation of degree-scale extended gamma-ray emission which traces the locations of recent sources in the Galaxy. Such emission in the energy band above 100 GeV is produced by CR nuclei and electrons released by the sources and spreading into the interstellar medium. We use the data from the Fermi gamma-ray telescope to locate the degree-scale 100 GeV gamma-ray sources. We find that the number of such sources and their overall power match to those expected when CRs injection events happen every ~100 yr in portions of ~1e50 erg. We find that most of the sources are associated to pulsars with spin down age less than ~30 kyr and hence to the recent supernova explosions. This supports the hypothesis of supernova origin of Galactic CRs. We notice that the degree-scale extended emission does not surround shell-like supernova remnants without pulsars. Based on this observation, we argue that the presence of the pulsar is essential for the CR acceleration process. We expect that a significant fraction of the degree-scale sources should be detectable as extended sources with km3-scale neutrino detectors.

preprint2012arXiv

Search for variable gamma-ray emission from the Galactic plane in the Fermi data

High-energy gamma-ray emission from the Galactic plane above ~100 MeV is composed of three main contributions: diffuse emission from cosmic ray interactions in the interstellar medium, emission from extended sources, such as supernova remnants and pulsar wind nebulae, and emission from isolated compact source populations. The diffuse emission and emission from the extended sources provide the dominant contribution to the flux almost everywhere in the inner Galaxy, preventing the detection of isolated compact sources. In spite of this difficulty, compact sources in the Galactic plane can be singled out based on the variability properties of their gamma-ray emission. Our aim is to find sources in the Fermi data that show long-term variability. We performed a systematic study of the emission variability from the Galactic plane, by constructing the variability maps. We find that emission from several directions along the Galactic plane is significantly variable on a time scale of months. These directions include, in addition to known variable Galactic sources and background blazars, the Galactic ridge region at positive Galactic longitudes and several regions containing young pulsars. We argue that variability on the time scale of months may be common to pulsars, originating from the inner parts of pulsar wind nebulae, similarly to what is observed in the Crab pulsar.

preprint2012arXiv

Superorbital modulation of X-ray emission from gamma-ray binary LSI +61 303

We report the discovery of a systematic constant time lag between the X-ray and radio flares of the gamma-ray binary LSI +61 303, persistent over long, multi-year, time scale. Using the data of monitoring of the system by RXTE we show that the orbital phase of X-ray flares from the source varies from $ϕ_X\simeq 0.35$ to $ϕ_X\simeq 0.75$ on the superorbital 4.6 yr time scale. Simultaneous radio observations show that periodic radio flares always lag the X-ray flare by $Δϕ_{X-R}\simeq 0.2$. We propose that the constant phase lag corresponds to the time of flight of the high-energy particle filled plasma blobs from inside the binary to the radio emission region at the distance ~10 times the binary separation distance. We put forward a hypothesis that the X-ray bursts correspond to the moments of formation of plasma blobs inside the binary system.

preprint2012arXiv

Ultra-High Energy Cosmic Ray production in the polar cap regions of black hole magnetospheres

We develop a model of ultra-high energy cosmic ray (UHECR) production via acceleration in a rotation-induced electric field in vacuum gaps in the magnetospheres of supermassive black holes (BH). We show that if the poloidal magnetic field near the BH horizon is misaligned with the BH rotation axis, charged particles, which initially spiral into the BH hole along the equatorial plane, penetrate into the regions above the BH "polar caps" and are ejected with high energies to infinity. We show that in such a model acceleration of protons near a BH of typical mass 3e8 solar masses is possible only if the magnetic field is almost aligned with the BH rotation axis. We find that the power of anisotropic electromagnetic emission from an UHECR source near a supermassive BH should be at least 10-100 times larger then UHECR power of the source. This implies that if the number of UHECR sources within the 100 Mpc sphere is ~100, the power of electromagnetic emission which accompanies proton acceleration in each source, $10^{42-43}$ erg/s, is comparable to the typical luminosities of active galactic nuclei (AGN) in the local Universe. We also explore the acceleration of heavy nuclei, for which the constraints on the electromagnetic luminosity and on the alignment of magnetic field in the gap are relaxed.

preprint2012arXiv

Very hard gamma-ray emission from a flare of Mrk 501

We investigate the peculiar properties of a large TeV gamma-ray flare from Mrk~501 detected during the 2009 multiwavelength campaign. We identify the counterpart of the flare in the Fermi/LAT telescope data and study its spectral and timing characteristics. A strong order-of-magnitude increase of the very-high-energy gamma-ray flux during the flare was not accompanied by an increase in the X-ray flux, so that the flare was one of the "orphan"-type TeV flares observed in BL Lacs. The flare lasted about 1 month at energies above 10 GeV. The flaring source spectrum in the 10-200 GeV range was very hard, with a photon index 1.1+/-0.2, harder than that observed in any other blazar in the gamma-ray band. No simultaneous flaring activity was detected below 10 GeV. Different variability properties of the emission below and above 10 GeV indicate the existence of two separate components in the spectrum. We investigate possible explanations of the very hard flaring component. We consider, among others, the possibility that the flare is produced by an electromagnetic cascade initiated by very-high-energy gamma-rays in the intergalactic medium. Within such an interpretation, peculiar spectral and temporal characteristics of the flare could be explained if the magnetic field in the intergalactic medium is of the order of 1e-17 - 1e-16 G.

preprint2012arXiv

Very-high-energy gamma-ray emission from high-redshift blazars

We study the possible detection of and properties of very high-energy (VHE) gamma-ray emission (in the energy band above 100 GeV) from high redshift sources. We report on the detection of VHE gamma-ray flux from blazars with redshifts z>0.5. We use the data of Fermi telescope in the energy band above 100 GeV and identify significant sources via cross-correlation of arrival directions of individual VHE gamma-rays with the positions of known Fermi sources. There are thirteen high-redshift sources detected in the VHE band by Fermi/LAT telescope. The present statistics of the Fermi signal from these sources is too low for a sensible study of the effects of suppression of the VHE flux by pair production through interactions with Extragalactic Background Light photons. We find that the detection of these sources with ground-based gamma-ray telescopes would be challenging. However, several sources including BL Lacs PKS 0426-380 at z=1.11, KUV 00311-1938 at z=0.61, B3 1307+433 at z=0.69, PG 1246+586 at z=0.84, Ton 116 at z=1.065 as well as a flat-spectrum radio quasar 4C +55.17 at z=0.89 should be detectable by HESS-II, MAGIC-II and CTA. A high-statistics study of a much larger number of VHE gamma-ray sources at cosmological distances would be possible with the proposed high-altitude Cherenkov telescope 5@5.

preprint2011arXiv

Discovery of High-Energy Gamma-Ray Emission from the Binary System PSR B1259-63/LS 2883 Around Periastron with Fermi

We report on the discovery of \geq 100 MeV γ rays from the binary system PSR B1259-63/LS 2883 using the Large Area Telescope (LAT) on board Fermi. The system comprises a radio pulsar in orbit around a Be star. We report on LAT observations from near apastron to ~ 60 days after the time of periastron, tp, on 2010 December 15. No γ-ray emission was detected from this source when it was far from periastron. Faint γ-ray emission appeared as the pulsar approached periastron. At ~ tp + 30d, the \geq 100 MeV γ-ray flux increased over a period of a few days to a peak flux 20-30 times that seen during the pre-periastron period, but with a softer spectrum. For the following month, it was seen to be variable on daily time scales, but remained at ~ 1 - 4 \times 10^-6 cm^-2 s^-1 before starting to fade at ~ tp + 57d. The total γ-ray luminosity observed during this period is comparable to the spin-down power of the pulsar. Simultaneous radio and X-ray observations of the source showed no corresponding dramatic changes in radio and X-ray flux between the pre-periastron and post-periastron flares. We discuss possible explanations for the observed γ-ray-only flaring of the source.

preprint2011arXiv

EGMF Constraints from Simultaneous GeV-TeV Observations of Blazars

Attenuation of the TeV gamma-ray flux from distant blazars through pair production with extragalactic background light leads to the development of electromagnetic cascades and subsequent, lower energy, GeV secondary gamma-ray emission. Due to the deflection of VHE cascade electrons by extragalactic magnetic fields (EGMF), the spectral shape of this arriving cascade gamma-ray emission is dependent on the strength of the EGMF. Thus, the spectral shape of the GeV-TeV emission from blazars has the potential to probe the EGMF strength along the line of sight to the object. We investigate constraints on the EGMF derived from observations of blazars for which TeV observations simultaneous with those by the Fermi telescope were reported. We study the dependence of the EGMF bound on the hidden assumptions it rests upon. We select blazar objects for which simultaneous Fermi/LAT GeV and Veritas, MAGIC or HESS TeV emission have been published. We model the development of electromagnetic cascades along the gamma-ray beams from these sources using Monte Carlo simulations, including the calculation of the temporal delay incurred by cascade photons, relative to the light propagation time of direct gamma-rays from the source. Constraints on EGMF could be derived from the simultaneous GeV-TeV data on the blazars RGB J0710+591, 1ES 0229+200, and 1ES 1218+304. The measured source flux level in the GeV band is lower than the expected cascade component calculated under the assumption of zero EGMF. Assuming that the reason for the suppression of the cascade component is the extended nature of the cascade emission, we find that B>10^{-15} G (assuming EGMF correlation length of ~1 Mpc) is consistent with the data. Alternatively, the assumption that the suppression of the cascade emission is caused by the time delay of the cascade photons the data are consistent with B>10^{-17} G for the same correlation length.

preprint2011arXiv

Low-energy break in the spectrum of Galactic cosmic rays

Measurements of the low energy spectrum of Galactic Cosmic Rays (GCR) by detectors at or near the Earth are affected by Solar modulation. To overcome this difficulty, we consider nearby molecular clouds as GCR detectors outside the Solar system. Using gamma-ray observations of the clouds by the Fermi telescope we derive the spectrum of GCRs in the clouds from the observed gamma-ray emission spectrum. We find that the GCR spectrum has a low energy break with the spectral slope hardening by 1.1+/-0.3 at an energy of 9+/-3 GeV. Detection of a low-energy break enables a measurement of GCR energy density in the interstellar space 0.9+/-0.3 eV/cm3.

preprint2010arXiv

Detection of very high energy gamma-ray emission from the Perseus cluster head-tail galaxy IC 310 by the MAGIC telescopes

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.

preprint2010arXiv

Galactic sources of E>100 GeV gamma-rays seen by Fermi telescope

We perform a search for sources of gamma-rays with energies E>100 GeV at low Galactic latitudes |b|<10 deg using the data of Fermi telescope. To separate compact gamma-ray sources from the diffuse emission from the Galaxy, we use the Minimal Spanning Tree method with threshold of 5 events in inner Galaxy (Galactic longitude |l|<60 deg) and of 3 events in outer Galaxy. Using this method, we identify 22 clusters of very-high-energy (VHE) gamma-rays, which we consider as "source candidates". 3 out of 22 event clusters are expected to be produced in result of random coincidences of arrival directions of diffuse background photons. To distinguish clusters of VHE events produced by real sources from the background we perform likelihood analysis on each source candidate. We present a list of 19 higher significance sources for which the likelihood analysis in the energy band E>100 GeV gives Test Statistics (TS) values above 25. Only 10 out of the 19 high-significance sources can be readily identified with previously known VHE gamma-ray sources. 4 sources could be parts of extended emission from known VHE gamma-ray sources. Five sources are new detections in the VHE band. Among these new detections we tentatively identify one source as a possible extragalactic source PMN J1603-4904 (a blazar candidate), one as a pulsar wind nebula around PSR J1828-1007. High significance cluster of VHE events is also found at the position of a source coincident with the Eta Carinae nebula. In the Galactic Center region, strong VHE gamma-ray signal is detected from Sgr C molecular cloud, but not from the Galactic Center itself.

preprint2010arXiv

Gamma ray emission from magnetized relativistic GRB outflows

We argue that small pitch angle synchrotron emission provides an important dissipation mechanism which has to be taken into account in the models of formation of relativistic magnetized gamma-ray burst (GRB) outflows from the newborn black holes and/or magnetars. We show that if the GRB outflow is proton loaded, spectral energy distribution of this emission is expected to sharply peak in 0.1-1~MeV energy band. If the small pitch angle synchrotron emission efficiently cools relativistic particles of the outflow, its spectrum below the peak energy is a powerlaw with spectral index alpha ~ -1, close to the typical spectral index of the time-resolved GRB spectra. Otherwise, the low energy spectral index can be as hard as alpha ~ 0, as observed at the beginning of the GRB pulses. We make a conjecture that small pitch angle synchrotron emission from proton-loaded magnetized GRB outflow could significantly contribute to the Band component of the prompt emission of GRBs while electromagnetic cascade initiated by the protons could be responsible for the GeV component.

preprint2010arXiv

INTEGRAL hard X-ray spectra of the cosmic X-ray background and Galactic ridge emission

We derive the spectra of the cosmic X-ray background (CXB) and of the Galactic ridge X-ray emission (GRXE) in the ~20-200 keV range from the data of the IBIS instrument aboard the INTEGRAL satellite obtained during the four dedicated Earth-occultation observations of early 2006. We analyse the modulation of the IBIS/ISGRI detector counts induced by the passage of the Earth through the field of view of the instrument. Unlike previous studies, we do not fix the spectral shape of the various contributions, but model instead their spatial distribution and derive for each of them the expected modulation of the detector counts. The spectra of the diffuse emission components are obtained by fitting the normalizations of the model lightcurves to the observed modulation in different energy bins. The obtained CXB spectrum is consistent with the historic HEAO-1 results and falls slightly below the spectrum derived with Swift/BAT. A 10% higher normalization of the CXB cannot be completely excluded, but it would imply an unrealistically high albedo of the Earth. The derived spectrum of the GRXE confirms the presence of a minimum around 80 keV with improved statistics and yields an estimate of ~0.6 M_Sun for the average mass of white dwarfs in the Galaxy. The analysis also provides updated normalizations for the spectra of the Earth's albedo and the cosmic-ray induced atmospheric emission.

preprint2010arXiv

No evidence for gamma-ray halos around active galactic nuclei resulting from intergalactic magnetic fields

We analyze the gamma-ray halo around stacked AGNs reported in Ap.J.Lett., 2010, 722, L39. First, we show that the angular distribution of gamma-rays around the stacked AGNs is consistent with the angular distribution of the gamma-rays around the Crab pulsar, which is a point source for Fermi/LAT. This makes it unlikely that the halo is caused by an electromagnetic cascade of TeV photons in the intergalactic space. We then compare the angular distribution of gamma-rays around the stacked AGNs with the point-spread function (PSF) of Fermi/LAT and confirm the existence of an excess above the PSF. However, we demonstrate that the magnitude and the angular size of this effect is different for photons converted in the front and back parts of the Fermi/LAT instrument, and thus is an instrumental effect.

preprint2010arXiv

Search for an extended VHE gamma-ray emission from Mrk 421 and Mrk 501 with the MAGIC Telescope

Context: Part of the very high energy $γ$-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.

preprint2010arXiv

The catalog of variable sources detected by INTEGRAL I: Catalog and Techniques

Context: In 6 years of operation, INTEGRAL/ISGRI revealed more than 500 sources. Many of these sources are variable. Taking into account that nearly half of INTEGRAL/ISGRI sources are new and many of them are still unidentified, the variability properties of the sources can serve as additional parameters that may help to classify and identify the unknown sources. Aims: In order to study the variability properties of the sources detected by INTEGRAL/ISGRI we develop a method to quantify the variability of a source. We describe here our techniques and compile a catalog of the sources that fit our criteria of variability. Methods: We use the natural time binning of INTEGRAL observations called Science Window ($\approx 2000$ seconds) and test the hypothesis that the detected sources are constant using a $χ^2$ all-sky map in three energy bands (20-40, 40-100, 100-200 keV). We calculate an intrinsic variance of the flux in individual pixels and use it to define the fractional variability of a source. The method is sensitive to the source variability on time scales of one Science Window and higher. We concentrate only on the sources which were already reported to be detected by INTEGRAL. Results: We present a catalog of 202 sources which are found to be significantly variable. For the catalog sources we give the measure of variability and fluxes with corresponding errors in 20-40, 40-100, 100-200 keV energy bands, and we present some statistics about the population of variable sources. The description of the physical properties of the variable sources will be given in a forthcoming paper.

preprint2010arXiv

Very high-energy gamma-ray emission from IC 310

We construct a systematic survey of extragalactic γ-ray sky at the energies above 100 GeV using the data of Fermi telescope. Such survey has not been previously done by the ground-based Cherenkov gamma-ray telescopes which have, contrary to Fermi, narrow field of view. We study a map of arrival directions of the highest energy photons detected by Fermi at Galactic latitudes |b| > 10 degrees and search for significant point source like excesses above the diffuse Galactic and extragalactic γ-ray backgrounds. We identify eight significant point source like excesses in this map. Seven of the eight sources are known TeV blazars. The previously unknown source is identified with a head-tail radio galaxy IC 310, situated in Perseus cluster of galaxies. The source is detected with significance 6 sigma above 30 GeV. We identify two possible scenaria for gamma-ray emission from this source. One possibility is that emission originates from the base of relativistic outflow from the active nucleus, as in the BL Lacs and FR I type radio galaxies. Otherwise gamma-ray photons could be produced at the bow shock formed in result of fast motion of the galaxy through the intracluster medium. The two models could be distinguished via the study of variability of the γ-ray signal.

preprint2009arXiv

A compact pulsar wind nebula model of the gamma-ray loud binary LS I +61 303

We study a model of of LS I +61 303 in which its radio to TeV emission is due to interaction of a relativistic wind from a young pulsar with the wind from its companion Be star. We assume the fast polar wind is clumpy, which is typical for radiatively-driven winds. The clumpiness cause the two winds to mix. The relativistic electrons from the pulsar wind are retained in the moving clumps by inhomogeneities of the magnetic field, which explains the X-ray variability observed on time scales much shorter than the orbital period. We calculate detailed inhomogeneous spectral models reproducing the average broad-band spectrum from radio to TeVs. Given the uncertainties the form of the distribution of relativistic electrons, the X-ray spectrum could be dominated by either Compton or synchrotron emission. The recent Fermi observations constrain the high-energy cut-off in the electron distribution to be at the Lorentz factor of 2 10^4 or 10^8 in the former and latter model, respectively. We provide formulae comparing the losses of the relativistic electrons due to Compton, synchrotron and Coulomb processes vs. the distance from the Be star. We calculate the optical depth of the wind to free-free absorption, showing that it will suppress most of the radio emission within the orbit. We point out the importance of Compton and Coulomb heating of the stellar wind within and around the gamma-ray emitting region. Then, we find the most likely mechanism explaining the orbital modulation at TeV energies is anisotropy of emission.

preprint2009arXiv

Exceptional flaring activity of the anomalous X-ray pulsar 1E 1547.0-5408

(Abridged) We studied an exceptional period of activity of the anomalous X-ray pulsar 1E 1547.0-5408 in January 2009, during which about 200 bursts were detected by INTEGRAL. The major activity episode happened when the source was outside the field of view of all the INTEGRAL instruments. But we were still able to study the properties of 84 bursts detected simultaneously by the anti-coincidence shield of the spectrometer SPI and by the detector of the imager ISGRI. We find that the luminosity of the 22 January 2009 bursts of 1E 1547.0-5408 was > 1e42 erg/s. This luminosity is comparable to that of the bursts of soft gamma repeaters (SGR) and is at least two orders of magnitude larger than the luminosity of the previously reported bursts from AXPs. Similarly to the SGR bursts, the brightest bursts of 1E 1547.0-5408 consist of a short spike of ~100 ms duration with a hard spectrum, followed by a softer extended tail of 1-10 s duration, which occasionally exhibits pulsations with the source spin period of ~2 s. The observation of AXP bursts with luminosities comparable to the one of SGR bursts strengthens the conjecture that AXPs and SGRs are different representatives of one and the same source type.

preprint2009arXiv

Where are Swift Gamma-ray bursts beyond the "synchrotron deathline"?

We study time-resolved spectra of the prompt emission of Swift Gamma-ray bursts (GRB). Our goal is to see if previous BATSE claims of the existence of a large amount of spectra with the low energy photon indices harder than 2/3 are consistent with Swift data. We perform a systematic search of the episodes of the spectral hardening down to the photon indices below 2/3 in the prompt emission spectra of Swift GRBs. We show that the data of the BAT instrument on board of Swift are consistent with BATSE data, if one takes into account differences between the two instruments. Much lower statistics of the very hard spectra in Swift GRBs is explained by the smaller field of view and narrower energy band of the BAT telescope.

preprint2009arXiv

X-ray observations of PSR B1259-63 near the 2007 periastron passage

PSR B1259-63 is a 48 ms radio pulsar in a highly eccentric 3.4 year orbit with a Be star SS 2883. Unpulsed gamma-ray, X-ray and radio emission components are observed from the binary system. It is likely that the collision of the pulsar wind with the anisotropic wind of the Be star plays a crucial role in the generation of the observed non-thermal emission. The 2007 periastron passage was observed in unprecedented details with Suzaku, Swift, XMM-Newton and Chandra missions. We present here the results of this campaign and compare them with previous observations. With these data we are able, for the first time, to study the details of the spectral evolution of the source over a 2 months period of the passage of the pulsar close to the Be star. New data confirm the pre-periastron spectral hardening, with the photon index reaching a value smaller than 1.5, observed during a local flux minimum. If the observed X-ray emission is due to the inverse Compton (IC) losses of the 10 MeV electrons, then such a hard spectrum can be a result of Coulomb losses, or can be related to the existence of the low-energy cut-off in the electron spectrum. Alternatively, if the X-ray emission is a synchrotron emission of very high energy electrons, the observed hard spectrum can be explained if the high energy electrons are cooled by IC emission in Klein-Nishina regime. Unfortunately the lack of simultaneous data in the TeV energy band prevents us from making a definite conclusion on the nature of the observed spectral hardening and, therefore, on the origin of the X-ray emission.

preprint2005arXiv

TeV signatures of compact UHECR accelerators

We study numerically particle acceleration by the electric field induced near the horizon of a rotating supermassive (M \sim 10^9-10^{10} M_{\odot}) black hole embedded in the magnetic field B. We find that acceleration of protons to energy E > 10^{20} eV is possible only at extreme values of M and B. We also find that the acceleration is very inefficient and is accompanied by a broad band MeV-TeV radiation whose total power exceeds at least by a factor of 1000 the total power emitted in ultra-high energy cosmic rays (UHECR). This implies that if O(10) nearby quasar remnants were sources of proton events with energy E > 10^{20} eV, then each quasar remnant would overshine e.g. the Crab nebula by more than two orders of magnitude in the TeV energy band. Recent TeV observations exclude this possibility. A model in which O(100) sources are situated at 100-1000 Mpc is not ruled out and can be experimentally tested by present TeV gamma-ray telescopes. Such a model can explain the observed UHECR flux at moderate energies E \sim (4-5) 10^{19} eV.

preprint2002arXiv

Large scale extragalactic jets powered by very-high energy gamma rays

The radiative cooling of electrons responsible for the nonthermal synchrotron emission of large scale jets of radiogalaxies and quasars requires quasi continuous (in time and space) production of relativistic electrons throughout the jets over the scales exceeding 100 kpc. While in the standard paradigm of large scale jets this implies in situ acceleration of electrons, in this letter we propose a principally different ``non-acceleration'' origin of these electrons, assuming that they are implemented all over the length of the jet through effective development of electromagnetic cascades initiated by extremely high energy gamma-rays injected into the jet from the central object. This scenario provides a natural and very economic way to power the jets up to distances of 100 kpc and beyond.