Search

We present the results of a cluster search in the gamma-ray sky images of the Large Magellanic Cloud (LMC) region by means of the Minimum Spanning Tree (MST) and DBSCAN algorithms, at energies higher than 6 and 10 GeV, using 12 years of Fermi-LAT data. Several significant clusters were found, the majority of which associated with previously known gamma-ray sources. We confirm our previous detection of the Supernova Remnants N 49B and N 63A and found new significant clusters associated with the SNRs N 49, N 186D and N 44. These sources are among the brightest X-ray remnants in the LMC and corresponds to core-collapse supernovae interacting with dense HII regions, indicating that an hadronic origin of high energy photons is the most likely process., Comment: 14 pages, 4 figures. Accepted for publication in MNRAS

We present a new elastic LIDAR concept, based on a bi-axially mounted Nd:YAG laser and a telescope with HPD readout, combined with fast FADC signal digitization and offline pulse analysis. The LIDAR return signals have been extensively quality checked and absolutely calibrated. We analyse seven years of quasi-continuous LIDAR data taken during those nights when the MAGIC telescopes were operating. Characterization of the nocturnal ground layer yields zenith and azimuth angle dependent aerosol extinction scaleheights for clear nights. We derive aerosol transmission statistics for light emitted from various altitudes throughout the year and separated by seasons. We find further seasonal dependencies of cloud base and top altitudes, but none for the LIDAR ratios of clouds. Finally, the night sky background light is characterized using the LIDAR photon backgrounds., The financial support of the Spanish grant PID2019-107847RB-C42, funded by the Spanish Ministry of Science and Innovation (MCIN)/AEI/10.13039/501100011033, the German Federal Ministry of Education and ResearchBMBF and the German Max-Planck-Gesellschaft (MPG), and by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project uniri-prirod-18-48 is gratefully acknowledged.

The gamma-ray emission of RX J1713–3946, despite being extensively studied in the GeV and TeV domains, remains poorly understood. This is mostly because, in this range, two competing mechanisms can efficiently produce gamma-rays: the inverse Compton scattering of accelerated electrons, and interactions of accelerated protons with the nuclei of the interstellar medium (ISM). In addition to the acceleration of particles from the thermal pool, the re-acceleration of pre-existing cosmic rays is often overlooked, and has in fact also been taken into account. Especially, because of the distance to the SNR (∼1 kpc), and the low density in which the shock is currently expanding (∼10−2 cm−3), the re-acceleration of cosmic-ray electrons pre-existing in the ISM can account for a significant fraction of the observed gamma-ray emission, and contribute to the shaping of the spectrum in the GeV–TeV range. Remarkably, this emission of leptonic origin is found to be close to the level of the gamma-ray signal in the TeV range, provided that the spectrum of pre-exisiting cosmic-ray electrons is similar to that observed in the local ISM. The overall gamma-ray spectrum of RX J1713–3946 is naturally produced as the sum of leptonic emission from re-accelerated cosmic-ray electrons, and a subdominant hadronic emission from accelerated protons. We also argue that neutrino observations with next-generation detectors might lead to a detection even in the case of a lepto–hadronic origin of the gamma-ray emission.

We report on results of a multi-band monitoring campaign from radio to gamma rays of the high-redshift flat spectrum radio quasar S5 0836+710 during a high activity period detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Two major flares were detected, in 2015 August and November. In both episodes, the apparent isotropic gamma-ray luminosity exceeds 10^50 erg/s, with a doubling time scale of about 3 hours. The high gamma-ray activity may be related to a superluminal knot that emerged from the core in 2015 April at the peak of the radio activity and is moving downstream along the jet. The low variability observed in X-rays may indicate that X-ray emission is produced by the low-energy tail of the same electron population that produces the gamma-ray emission. The analysis of full-polarization pc-scale radio observations suggests the presence of a limb-brightened polarization structure at about 1 mas from the core in which a rotation measure gradient with a sign change is observed transverse to the jet direction. These characteristics are consistent with a scenario in which Faraday rotation is produced by a sheath of thermal electrons with a toroidal magnetic field surrounding the emitting jet., Comment: 26 pages, 13 figures, and 7 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

TeV photons from extragalactic sources are absorbed in the intergalactic medium and initiate electromagnetic cascades. These cascades offer a unique tool to probe the properties of the universe at cosmological scales. We present a new Monte Carlo code dedicated to the physics of such cascades. This code has been tested against both published results and analytical approximations, and is made publicly available. Using this numerical tool, we investigate the main cascade properties (spectrum, halo extension, time delays), and study in detail their dependence on the physical parameters (extra-galactic magnetic field, extra-galactic background light, source redshift, source spectrum and beaming emission). The limitations of analytical solutions are emphasised. In particular, analytical approximations account only for the first generation of photons and higher branches of the cascade tree are neglected., Comment: 18 pages, 21 figures, Accepted by MNRAS; this preprint matches accepted version

We investigate the prompt emission and the afterglow properties of short-duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multiwavelength facilities including the Gran Canarias Telescope 10.4 m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL-SPI-ACS, Swift-BAT, and Fermi-GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A, and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB 130603B, our new afterglow photometric data constrain the pre-jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new and previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for eight other sGRBs we are able to: place tight constraints on the non-detection of the afterglow, host galaxy, or any underlying 'kilonova' emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, AJCT acknowledges support from the Junta de Andalucia (Project P07-TIC-03094) and support from the Spanish Ministry Projects AYA2012-39727-C03-01 and 201571718R. This work has been supported by the Spanish Science Ministry 'Centro de Excelencia SeveroOchoa' Program under grant SEV-2017-0709. FEDER funds are acknowledged. E.S. acknowledges assistance from the Scientific and Technological Research Council of Turkey (TUBITAK) through project 112T224. We thank TUBITAK for a partial support in using T100 telescope with project number 10CT100-95. A.S.P acknowledges partial support grants RFBR 17-02-01388, 17-51-44018, and 1752-80139. E.D.M., A.A.V., and P.Yu.M. are grateful to RSCF grant 18-12-00522 for support. B.-B.Z. acknowledges support from National Thousand Young Talents program of China and National Key Research and Development Program of China (2018YFA0404204). R.Ya.I. is grateful for partial support by the grant RUSTAVELI/FR/379/6-300/14. R.S.R. acknowledges support from ASI (Italian Space Agency) through the Contract No. 2015-046R.0 and from European Union Horizon 2020 Programme under the AHEAD project (grant agreement No. 654215). SJ acknowledges the support of the Korea Basic Science Research Program through NRF-2015R1D1A4A01020961.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi (γ -rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented γ -ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 ± 63) × 10−8 ph cm−2 s−1, corresponding to a luminosity of (2.2 ± 0.1) × 1050 erg s−1, was reached on 2016 December 28. These four γ -ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and γ -ray activity is found. The γ -ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This γ -ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission. © 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, The data collected by the WEBT collaboration are stored in the WEBT archive at the Osservatorio Astrofisico di Torino INAF (http://www.oato.inaf.it/blazars/webt/); for questions regarding their availability, contact the WEBT President Massimo Villata(massimo.villata@inaf.it).We acknowledge financial contribution from the agreementASI-INAFn. 2017-14-H.0 and from the contract PRIN-SKA-CTA-INAF 2016. This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants DN 08-1/2016, DN 18-13/2017, and KP-06-H28/3 (2018). The Skinakas Observatory is a collaborative project of the University of Crete, the Foundation for Research and Technology -Hellas, and the Max-Planck-Institut fur Extraterrestrische Physik. The Abastumani team acknowledges financial support by the Shota Rustaveli National Science Foundation under contract FR/217554/16. The St. Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. GD and OV gratefully acknowledge the observing grant support from the Institute of Astronomy and Rozhen National Astronomical Observatory, Bulgarian Academy of Sciences, via bilateral joint research project 'Study of ICRF radio-sources and fast variable astronomical objects' (head -G.Damljanovic). This work is a part of the Projects No. 176011 ('Dynamics and Kinematics of Celestial Bodies and Systems'), No. 176004 ('Stellar Physics'), and No. 176021 ('Visible and Invisible Matter in Nearby Galaxies: Theory and Observations') supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia. JE is indebted to DGAPA (Universidad Nacional Autonoma de M ' exico) for financial support, PAPIIT project IN114917. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. Data from the Steward Observatory blazar monitoring project were used. This program is supported by NASA/Fermi Guest Investigator grants NNX12AO93G and NNX15AU81G. We acknowledge support by Bulgarian National Science Programme 'Young Scientists and Postdoctoral Students 2019', Bulgarian National Science Fund under grant DN18-10/2017 and National RI Roadmap Projects DO1-157/28.08.2018 and DO1-153/28.08.2018 of the Ministry of Education and Science of the Republic of Bulgaria. This publication makes use of data obtained at Mets ahovi Radio Observatory, operated by Aalto University in Finland. The Astronomical Observatory of the Autonomous Region of theAostaValley (OAVdA) is managed by the Fondazione Clment Fillietroz-ONLUS, which is supported by the Regional Government of the Aosta Valley, the TownMunicipality of Nus and the 'Unit des Communes valdtaines Mont-milius'. The research at the OAVdA was partially funded by two 'Research and Education' grants from Fondazione CRT. RR acknowledges support from CONICYT project Basal AFB-170002. MM and TM acknowledge support through the Russian Government Program of Competitive Growth of Kazan Federal University. The Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat a l'Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucleaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Etudes Spatiales in France. This work performed in part under DOE Contract DE-AC02-76SF00515. The OVRO 40-m monitoring program is supported in part by NASA grants NNX08AW31G, NNX11A043G, and NNX14AQ89G, and NSF grants AST-0808050 and AST-1109911. We thank the Swift team for making these observations possible, the duty scientists, and science planners. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. FD thanks S. Covino for his help with the REM data reduction. This research has made use of data obtained from the high-energy Astrophysics Science Archive Research Center (HEASARC) provided by NASA's Goddard Space Flight Center.

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 $\delta \lambda < \lambda$ 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)., Comment: V2: 7 pages, 4 figures. Considerable changes in the blazar sample, results unchanged; discussions extended - version accepted by A&A

We make a first attempt to find dwarf galaxies in eight \Fermi-LAT extended, unassociated, source fields using \Gaia\ DR2. We probe previously unexplored heliocentric distances of $d, 9 pages, 6 figures, accepted for publication in MNRAS

PKS 0521-36 is an active galactic nucleus (AGN) with uncertain classification. We investigate the properties of this source from radio to γ-rays. The broad emission lines in the optical and ultraviolet bands and steep radio spectrum indicate a possible classification as an intermediate object between broad-line radio galaxies (BLRG) and steep spectrum radio quasars (SSRQ). On pc-scales PKS 0521-36 shows a knotty structure similar to misaligned AGN. The core dominance and the γ-ray properties are similar to those estimated for other SSRQ and BLRG detected in γ-rays, suggesting an intermediate viewing angle with respect to the observer. In this context the flaring activity detected from this source by Fermi-Large Area Telescope between 2010 June and 2012 February is very intriguing. We discuss the γ-ray emission of this source in the framework of the structured jet scenario, comparing the spectral energy distribution (SED) of the flaring state in 2010 June with that of a low state. We present three alternative models corresponding to three different choices of the viewing angles θv = 6°, 15°, and 20°. We obtain a good fit for the first two cases, but the SED obtained with θv = 15° if observed at a small angle does not resemble that of a typical blazar since the synchrotron emission should dominate by a large factor (̃100) the inverse Compton component. This suggests that a viewing angle between 6° and 15° is preferred, with the rapid variability observed during γ-ray flares favouring a smaller angle. However, we cannot rule out that PKS 0521-36 is the misaligned counterpart of a synchrotron-dominated blazar.

The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85 GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes. The observations, taken between 2008 May and 2010 June and totalling 53 hours, resulted in the discovery of $\gamma$-ray emission from the blazar, which has a redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of 5.5 standard deviations above the background with an integral flux of (2.8$\pm0.7_{\mathrm{stat}}\pm0.8_{\mathrm{sys}}$) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ (1.2\% of the Crab Nebula's flux) above 200 GeV. The measured spectrum is described well by a power law from 0.2 TeV to 1.3 TeV with a photon index of 3.1 $\pm$ 0.4$_{\mathrm{stat}}$ $\pm$ 0.2$_{\mathrm{sys}}$. Quasi-simultaneous multi-wavelength data from the Fermi Large Area Telescope (0.3--300 GeV) and the Swift X-ray Telescope (0.2--10 keV) are additionally used to model the properties of the emission region. A synchrotron self-Compton model produces a good representation of the multi-wavelength data. Adding an external-Compton or a hadronic component also adequately describes the data., Comment: 8 pages, 4 figures. Accepted for publication in MNRAS

The source-count distribution as a function of their flux, dN/dS, is one of the main quantities characterizing gamma-ray source populations. We employ statistical properties of the Fermi Large Area Telescope (LAT) photon counts map to measure the composition of the extragalactic gamma-ray sky at high latitudes (|b|>30 deg) between 1 GeV and 10 GeV. We present a new method, generalizing the use of standard pixel-count statistics, to decompose the total observed gamma-ray emission into (a) point-source contributions, (b) the Galactic foreground contribution, and (c) a truly diffuse isotropic background contribution. Using the 6-year Fermi-LAT data set (P7REP), we show that the dN/dS distribution in the regime of so far undetected point sources can be consistently described with a power law of index between 1.9 and 2.0. We measure dN/dS down to an integral flux of ~2x10^{-11} cm^{-2}s^{-1}, improving beyond the 3FGL catalog detection limit by about one order of magnitude. The overall dN/dS distribution is consistent with a broken power law, with a break at 2.1^{+1.0}_{-1.3}x10^{-8} cm^{-2}s^{-1}. The power-law index n_1=3.1^{+0.7}_{-0.5} for bright sources above the break hardens to n_2=1.97+-0.03 for fainter sources below the break. A possible second break of the dN/dS distribution is constrained to be at fluxes below 6.4x10^{-11} cm^{-2}s^{-1} at 95% confidence level. The high-latitude gamma-ray sky between 1 GeV and 10 GeV is shown to be composed of ~25% point sources, ~69.3% diffuse Galactic foreground emission, and ~6% isotropic diffuse background., 28 pages, 20 figures, 6 tables; v2: minor changes and additions, accepted for publication in ApJS

The current velocity of the Smith Cloud indicates that it has undergone at least one passage of the Galactic disc. Using hydrodynamic simulations we examine the present day structure of the Smith Cloud. We find that a dark matter supported cloud is able to reproduce the observed present day neutral hydrogen mass, column density distribution and morphology. In this case the dark matter halo becomes elongated, owing to the tidal interaction with the Galactic disc. Clouds in models neglecting dark matter confinement are destroyed upon disc passage, unless the initial cloud mass is well in excess of what is observed today. We then determine integrated flux upper limits to the gamma-ray emission around such a hypothesised dark matter core in the Smith Cloud. No statistically significant core or extended gamma-ray emission are detected down to a 95% confidence level upper limit of $1.4\times10^{-10}$ ph cm$^{-2}$ s$^{-1}$ in the 1-300 GeV energy range. For the derived distance of 12.4 kpc, the Fermi upper limits set the first tentative constraints on the dark matter cross sections annihilating into $\tau^{+}{\tau}^{-}$ and $b\bar{b}$ for a high-velocity cloud., Comment: 10 pages, 8 figures. Submitted to MNRAS

We report on multifrequency observations of the gamma-ray emitting narrow-line Seyfert 1 galaxy PKS 1502+036 performed from radio to gamma-rays during 2008 August-2012 November by Fermi-LAT, Swift (XRT and UVOT), OVRO, VLBA, and VLA. No significant variability has been observed in gamma-rays, with 0.1-100 GeV flux that ranged between (3-7)x10^-8 ph/cm^2/s using 3-month time bins. The photon index of the LAT spectrum (Gamma=2.60+/-0.06) and the apparent isotropic gamma-ray luminosity, L(0.1-100 GeV)= 7.8x10^45 erg/s, over 51 months are typical of a flat spectrum radio quasar. The radio spectral variability and the one-sided structure, in addition to the observed gamma-ray luminosity, suggest a relativistic jet with a high Doppler factor. In contrast to SBS 0846+513, the VLBA at 15 GHz did not observe superluminal motion for PKS 1502+036. Despite having the optical characteristics typical of a narrow-line Seyfert 1 galaxy, radio and gamma-ray properties of PKS 1502+036 are found to be similar to those of a blazar at the low end of the black hole mass distribution for blazars. This is in agreement with what has been found in the case of the other gamma-ray emitting narrow-line Seyfert 1 SBS 0846+513., Comment: 11 pages, 9 figures, 6 tables. MNRAS in press - Minor editing changes

High redshift blazars are among the most powerful objects in the Universe. Although they represent a significant fraction of the extragalactic hard X-ray sky, they are not commonly detected in γ-rays. High redshift ( z > 2 ) objects represent < 10 per cent of the active galactic nuclei (AGN) population observed by Fermi so far, and γ-ray flaring activity from these sources is even more uncommon. The characterization of the radio-to-γ-ray properties of high redshift blazars represents a powerful tool for the study of the energetics of such extreme objects and the Extragalactic Background Light. This contribution will present results of multi-band campaigns, from radio to γ-rays, on PKS 0836+710, PKS 2149-306, and TXS 0536+145. The latter is the highest redshift detection of a flaring γ-ray blazar so far. At the peaks of their respective flares these sources reached an apparent isotropic gamma-ray luminosity of about 10 50 erg·s - 1 , which is comparable with the luminosity observed from the most powerful blazars. The physical properties derived from the multi-wavelength observations of these sources are then compared with those shown by the high redshift population.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties., 17 pages, 9 figures, 5 tables, Accepted for publication in ApJS