Search

For stars hosting Circumstellar Habitable Zone (CHZ) exoplanets, we investigate the time-evolution of their ultraviolet habitable zone (UHZ), the annular region around a star where an exoplanet could experience a suitable ultraviolet environment for the presence and emergence of life, and the possible intersection of the UHZ with the CHZ. To estimate their UV luminosity evolution, and therefore the evolution of their UHZ, we analyse Swift-UV/Optical telescope observations and adopt the near-UV luminosity evolutionary tracks derived using GALEX observations of young moving groups. We find that an intersection between CHZ and UHZ could exist (or have existed) around all stars of our sample at different epochs, except for the coldest M-dwarfs (temperature < 2800 K, e.g. Trappist-1). For hotter M-dwarfs the formation of RNA precursors through cyanosulfidic chemistry triggered by near-UV radiation could occur during the first 1-2 Gyrs. The radial-extension and time-duration of the CHZ-UHZ intersection increase with the stellar effective temperature and the exoplanet atmospheric transmissivity at near-UV wavelengths. Within our sample, Proxima Centauri represents a golden target for the quest of life outside the Solar system because it experienced a long-lasting and more extended, compared to similar M-dwarfs, CHZ-UHZ intersection., Comment: published in MNRAS Letters, 7 pages, 4 figures

The study of the high-energy (MeV-GeV) part of GRBs spectrum can play a crucial role in investigating the physics of the prompt emission, but it is often hampered by low statistic and the paucity of GeV observations. In this work, we analyze the prompt emission spectra of the 22 brightest GRBs which have been simultaneously observed by Fermi/GBM and Fermi/LAT, spanning 6 orders of magnitude in energy. The high-energy photon spectra can be modelled with a power-law $N(E)\propto E^{-\beta}$ possibly featuring an exponential cutoff. We find that, with the inclusion of the LAT data, the spectral index $\beta$ is softer than what typically inferred from the analysis of Fermi/GBM data alone. Under the assumption that the emission is synchrotron, we derive the index $p\sim2.79$ of the power-law energy distribution of accelerated particles ($N(\gamma)\propto \gamma^{-p}$). In 9 out of 22 GRB spectra, we find a significant presence of an exponential cut-off at high-energy, ranging between 14 and 298 MeV. By interpreting the observed cut-off as a sign of pair-production opacity, we estimate the jet bulk Lorentz factor $\Gamma$, finding values in the range 130-330. These values are consistent with those inferred from the afterglow light curve onset time. Finally, by combining the information from the high-energy prompt emission spectrum with the afterglow lightcurve, we provide a method to derive the distance R from the central engine where the prompt emission occurs. These results highlight the importance of including high-energy data, when available, in the study of prompt spectra and their role in addressing the current challenges of the GRB standard model., Comment: 17 pages, 17 figures. Accepted for publication in A&A

Studies of the most distant AGNs allow us to test our current understanding of the physics present in radio-jetted AGNs across a range of environments. The decrease in apparent luminosity with distance is the primary difficulty to overcome in the study of these distant AGNs, which requires highly sensitive instruments. Our goal is to employ new long wavelength radio data to better parametrise the broad-band SED of GB 1508+5714, a high-redshift (z=4.30) AGN. Its high redshift, high intrinsic luminosity, and classification as a blazar allow us to test emission models that consider the efficient cooling of jet electrons via inverse Compton losses in interactions with the dense CMB photon field at high redshifts. A significant detection of this effect in GB 1508+5714 may partly explain the apparent sparsity of high-redshift radio galaxies in wide-field surveys; detections of this kind are only becoming possible with the current generation of SKA precursors. We used international LOFAR telescope to image the long wavelength radio emission around the high-redshift blazar GB 1508+5714 on arcsecond scales at frequencies between 128 MHz and 160 MHz. This allowed us to compare the spatially resolved structure with higher frequency observations, and to construct spectral index maps. The LOFAR image shows a compact unresolved core and two resolved emission regions around 2 arcsec to the east and to the west of the radio core. We find structure consistent with previous VLA observations, as well as a previously unreported emission region to the east. We interpret the arcsecond-scale radio structure of GB 1508+5714 as a FR II-like radio galaxy at a small viewing angle. Our SED modelling shows that a scenario featuring significant quenching effects caused by interaction with the CMB provides a good description of the data, and notably explains the suppressed radio emission., Comment: 11 pages, 10 figures, 2 tables

We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $\nu L_\nu\propto \nu^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $\nu L_\nu \propto \nu^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $\gamma_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $\gamma_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $\gamma_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $\gamma_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare., Comment: accepted for publication in Astronomy & Astrophysics; preview abstract shortened due to arXiv requirements

Gamma-ray Bursts (GRBs) prompt emission spectra are often fitted with the empirical ''Band" function, namely two power laws smoothly connected. The typical slope of the low energy (sub-MeV) power law is $\alpha_{B}\simeq -1$. In a small fraction of long GRBs this power law splits into two components such that the spectrum presents, in addition to the typical $\sim$ MeV $\nu F_{\nu}$ peak, a break at the order of a few keV or hundreds keV. The typical power law slopes below and above the break are -0.6 and -1.5 respectively. If the break is a common feature, the value of $\alpha_{B}$ could be an ''average'' of the spectral slopes below and above the break in GRBs fitted with Band function. We analyze the spectra of 27 (9) bright long (short) GRBs detected by the Fermi satellite finding a low energy break between 80 keV and 280 keV in 12 long GRBs, but in none of the short events. Through spectral simulations we show that if the break is moved closer (farther) to the peak energy a relatively harder (softer) $\alpha_{B}$ is found by fitting the simulated spectra with the Band function. The hard average slope $\alpha_{B}\simeq-0.38$ found in short GRBs suggests that the break is close to the peak energy. We show that for 15 long GRBs best fitted by the Band function only, the break could be present, but it is not identifiable in the Fermi/GBM spectrum, because either at low energies, close to the detector limit for relatively soft $\alpha_{B}\lesssim-1$, or in the proximity of the energy peak for relatively hard $\alpha_{B}\gtrsim-1$. A spectrum with two breaks could be typical of GRB prompt emission, though hard to identify with current detectors. Instrumental design such that conceived for the THESEUS space mission, extending from 0.3 keV to several MeV and featuring a larger effective area with respect to Fermi/GBM, can reveal a larger fraction of GRBs with a spectral energy break.

THESEUS is a medium size space mission of the European Space Agency, currently under evaluation for a possible launch in 2032. Its main objectives are to investigate the early Universe through the observation of gamma-ray bursts and to study the gravitational waves electromagnetic counterparts and neutrino events. On the other hand, its instruments, which include a wide field of view X-ray (0.3-5 keV) telescope based on lobster-eye focusing optics and a gamma-ray spectrometer with imaging capabilities in the 2-150 keV range, are also ideal for carrying out unprecedented studies in time domain astrophysics. In addition, the presence onboard of a 70 cm near infrared telescope will allow simultaneous multi-wavelegth studies. Here we present the THESEUS capabilities for studying the time variability of different classes of sources in parallel to, and without affecting, the gamma-ray bursts hunt., Comment: Submitted to Experimental Astronomy

Counted among the most powerful cosmic events, supernovae (SNe) and gamma-ray bursts (GRBs) can be highly disruptive for life: Their radiation can be harmful for biota or induce extinction by removing most of the protective atmospheric ozone layer from terrestrial planets (TPs). Nearby high-energy transient astrophysical events have been proposed as possible triggers of mass extinctions on Earth. We assess the habitability of the Milky Way (MW) throughout its cosmic history against potentially disruptive astrophysical transients with the aim of identifying the safest places and epochs within our Galaxy. We also test the hypothesis that one long GRB played a leading role in the late Ordovician mass-extinction event (~445 Myr ago). We characterised the habitability of the MW throughout its cosmic history as a function of galactocentric distance of TPs. We estimated the dangerous effects of transient astrophysical events (long and short GRBs and SNe) with a model that connects their rate to the specific star formation and metallicity evolution within the Galaxy throughout its cosmic history. Our model also accounts for the probability that TPs form around FGK and M stars. Until about six billion years ago, the outskirts of the Galaxy were the safest places to live, despite the relatively low density of TPs. In the last about four billion years, regions between 2 and 8 kpc from the center, which had a higher density of TPs, became the best places for a relatively safer biotic life growth. We confirm the hypothesis that one long GRB played a leading role in the late Ordovician mass-extinction event. In the last 500 Myr, the safest neighborhood in the Galaxy was a region at a distance of 2 to 8 kpc from the Galactic center, whereas the MW outskirts were sterilized by two to five long GRBs., Comment: 11 pages, 10 figures, 4 tables. Accepted for pubblication in Astronomy & Astrophysics

We present the discovery of PSO J030947.49+271757.31, the radio brightest (23.7 mJy at 1.4 GHz) active galactic nucleus (AGN) at z>6.0. It was selected by cross-matching the NRAO VLA Sky Survey and the Panoramic Survey Telescope and Rapid Response System PS1 databases and its high-z nature was confirmed by a dedicated spectroscopic observation at the Large Binocular Telescope. A pointed Neil Gehrels $Swift$ Observatory XRT observation allowed us to measure a flux of $\sim$3.4$\times$10$^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the [0.5-10] keV energy band, which also makes this object the X-ray brightest AGN ever observed at z>6.0. Its flat radio spectrum ($\alpha_{\nu r}$<0.5), very high radio loudness (R>10$^3$), and strong X-ray emission, compared to the optical, support the hypothesis of the blazar nature of this source. %i.e. a radio-loud (RL) AGN with the relativistic jet pointed toward us. Assuming that this is the only blazar at this redshift in the surveyed area of sky, we derive a space density of blazars at z$\sim$6 and with M$_{1450 \mbox{\AA}}$ < -25.1 of 5.5$^{+11.2}_{-4.6}$$\times$10$^{-3}$ Gpc$^{-3}$. From this number, and assuming a reasonable value of the bulk velocity of the jet ($\Gamma$=10), we can also infer a space density of the entire radio-loud AGN population at z$\sim$6 with the same optical/UV absolute magnitude of 1.10$^{+2.53}_{-0.91}$ Gpc$^{-3}$. Larger samples of blazars will be necessary to better constrain these estimates., Comment: 7 pages, 8 figures, 3 tables, Accepted to publication in A&A Letters. Typo corrections

We discuss the new surprising observational results that indicate quite convincingly that the prompt emission of Gamma-Ray Bursts (GRBs) is due to synchrotron radiation produced by a particle distribution that has a low energy cut-off. The evidence of this is provided by the low energy part of the spectrum of the prompt emission, that shows the characteristic F(nu) \propto nu^(1/3) shape followed by F(nu) \propto nu^(-1/2) up to the peak frequency. This implies that although the emitting particles are in fast cooling, they do not cool completely. This poses a severe challenge to the basic ideas about how and where the emission is produced, because the incomplete cooling requires a small value of the magnetic field, to limit synchrotron cooling, and a large emitting region, to limit the self-Compton cooling, even considering Klein-Nishina scattering effects. Some new and fundamental ingredient is required for understanding the GRBs prompt emission. We propose proton-synchrotron as a promising mechanism to solve the incomplete cooling puzzle., Comment: 7 pages, 3 figures, accepted for publication in A&A

GrailQuest (Gamma Ray Astronomy International Laboratory for QUantum Exploration of Space-Time) is a mission concept based on a constellation (hundreds/thousands) of nano/micro/small-satellites in low (or near) Earth orbits. Each satellite hosts a non-collimated array of scintillator crystals coupled with Silicon Drift Detectors with broad energy band coverage (keV-MeV range) and excellent temporal resolution ( below or equal 100 nanoseconds) each with effective area around 100 cm2. This simple and robust design allows for mass-production of the satellites of the fleet. This revolutionary approach implies a huge reduction of costs, flexibility in the segmented launching strategy, and an incremental long-term plan to increase the number of detectors and their performance: a living observatory for next-generation, space-based astronomical facilities. GrailQuest is conceived as an all-sky monitor for fast localisation of high signal-to-noise ratio transients in the X/gamma-ray band, e.g. the elusive electromagnetic counterparts of gravitational wave events. Robust temporal triangulation techniques will allow unprecedented localisation capabilities, in the keV-MeV band, of a few arcseconds or below, depending on the temporal structure of the transient event. The ambitious ultimate goal of this mission is to perform the first experiment, in quantum gravity, to directly probe space-time structure down to the minuscule Planck scale, by constraining or measuring a first order dispersion relation for light in vacuo. This is obtained by detecting delays between photons of different energies in the prompt emission of Gamma-ray Bursts., Comment: White paper submitted in response to the ESA call Voyage 2050 - Submitted to Experimental Astronomy

The Gamma Ray Burst (GRB) 180720B is one of the brightest events detected by the Fermi satellite and the first GRB detected by the H.E.S.S. telescope above 100 GeV. We analyse the Fermi (GBM and LAT) and Swift (XRT and BAT) data and describe the evolution of the burst spectral energy distribution in the 0.5 keV - 10 GeV energy range over the first 500 seconds of emission. We reveal a smooth transition from the prompt phase, dominated by synchrotron emission in a moderately fast cooling regime, to the afterglow phase whose emission has been observed from the radio to the GeV energy range. The LAT (0.1 - 100 GeV) light curve initially rises ($F_{\rm LAT}\propto t^{2.4}$), peaks at $\sim$78 s, and falls steeply ($F_{\rm LAT}\propto t^{-2.2}$) afterwards. The peak, which we interpret as the onset of the fireball deceleration, allows us to estimate the bulk Lorentz factor $\Gamma_{0}\sim 150 \ (300)$ under the assumption of a wind-like (homogeneous) circum-burst medium density. We derive a flux upper limit in the LAT energy range at the time of H.E.S.S. detection, but this does not allow us to unveil the nature of the high energy component observed by H.E.S.S. We fit the prompt spectrum with a physical model of synchrotron emission from a non-thermal population of electrons. The 0 - 35 s spectrum after its $E F(E)$ peak (at 1 - 2 MeV) is a steep power law extending to hundreds of MeV. We derive a steep slope of the injected electron energy distribution $N(\gamma)\propto \gamma^{-5}$. Our fit parameters point towards a very low magnetic field ($B'\sim 1 $ G) in the emission region., Comment: 10 pages, 6 figures, submitted to A&A

We present the discovery and properties of DESJ014132.4-542749.9 (DES0141-54), a new powerful radio-loud active galactic nucleus (AGN) in the early Universe (z=5.0). It was discovered by cross-matching the first data release of the Dark Energy Survey (DES DR1) with the Sidney University Molonglo Survey (SUMSS) radio catalog at 0.843 GHz. This object is the first radio-loud AGN at high redshift discovered in the DES. The radio properties of DES0141-54, namely its very large radio-loudness (R>10$^{4}$), the high radio luminosity (L$_{0.8 GHz}$=1.73$\times$10$^{28}$ W Hz$^{-1}$), and the flatness of the radio spectrum ($\alpha$=0.35) up to very high frequencies (120 GHz in the source's rest frame), classify this object as a blazar, meaning, a radio-loud AGN observed along the relativistic jet axis. However, the X--ray luminosity of DESJ0141-54 is much lower compared to those of the high redshift (z$\geq$4.5) blazars discovered so far. Moreover its X-ray-to-radio luminosity ratio (log($\frac{L_{[0.5-10]keV}}{L_{1.4GHz}}$)=9.96$\pm$0.30 Hz) is small also when compared to lower redshift blazars: only 2\% of the low-z population has a similar ratio. By modeling the spectral energy distribution we found that this peculiar X--ray weakness and the powerful radio emission could be related to a particularly high value of the magnetic field. Finally, the mass of the central black hole is relatively small (M$_{BH}$ = 3-8 $\times$10$^8$ M$_{\odot}$) compared to other confirmed blazars at similar redshift, making DES0141-54 the radio-loud AGN that host the smallest supermassive black hole ever discovered at z$\geq$5., Comment: 12 pages, 10 figures, 6 tables. Accepted for publication in Astronomy and Astrophysics

In the last few years many exoplanets in the habitable zone (HZ) of M-dwarfs have been discovered, but the X-ray/UV activity of cool stars is very different from that of our Sun. The high-energy radiation environment influences the habitability, plays a crucial role for abiogenesis, and impacts planetary atmospheres. LHS 1140b is a super-Earth-size planet orbiting in the HZ of LHS 1140, an M4.5 dwarf at ~15 parsecs. We present the results of a Swift X-ray/UV observing campaign. We characterize for the first time the X-ray/UV radiation environment of LHS 1140b. We measure the variability of the near ultraviolet (NUV) flux and estimate the far ultraviolet (FUV) flux with a correlation between FUV and NUV flux of a sample of low-mass stars in the GALEX archive. We highlight the presence of a dominating X-ray source close to the J2000 coordinates of LHS 1140, characterize its spectrum, and derive an X-ray flux upper limit for LHS 1140. We find that this contaminant source could have influenced the previously estimated spectral energy distribution. No significant variation of the NUV flux of LHS 1140 is found over 3 months, and we do not observe any flare during the 38 ks on the target. LHS 1140 is in the 25th percentile of least variable M4-M5 dwarfs of the GALEX sample. Analyzing the UV flux experienced by the HZ planet LHS 1140b, we find that outside the atmosphere it receives a NUV flux <2% with respect to that of the present-day Earth, while the FUV/NUV ratio is ~100-200 times higher. This represents a lower limit to the true FUV/NUV ratio since the GALEX FUV band does not include Lyman-alpha, which dominates the FUV output of low-mass stars. This is a warning for future searches for biomarkers, which must take into account this high ratio. The relatively low level and stability of UV flux experienced by LHS 1140b should be favorable for its present-day habitability., Comment: 6 pages, 4 figures, 2 tables. Accepted for pubblication in Astronomy & Astrophysics

We observed with the NuSTAR satellite 3 blazars at z>2, detected in gamma-rays by Fermi/LAT and in the soft X-rays, but not yet observed above 10 keV. The flux and slope of their X-ray continuum, together with Fermi/LAT data allows us to estimate their total electromagnetic output and peak frequency. For some of them we can study the source in different states, and investigate the main cause of the observed different spectral energy distribution. We then collected all blazars at redshift greater than 2 observed by NuSTAR, and confirm that these hard and luminous X-ray blazars are among the most powerful persistent sources in the Universe. We confirm the relation between the jet power and the disk luminosity, extending it at the high energy end., Comment: 10 pages, 13 figures, accepted for publication in A&A

The peculiar short gamma-ray burst GRB 170817A has been firmly associated to the gravitational wave event GW170817, which has been unaninmously interpreted as due to the coalescence of a double neutron star binary. The unprecedented behaviour of the non-thermal afterglow led to a debate about its nature, which was eventually settled by high-resolution VLBI observations, which strongly support the off-axis structured jet scenario. Using information on the jet structure derived from multi-wavelength fitting of the afterglow emission and of the apparent VLBI image centroid motion, we compute the appearance of a GRB 170817A-like jet as seen by an on-axis observer and we compare it to the previously observed population of SGRB afterglows and prompt emission events. We find that the intrinsic properties of the GRB 170817A jet are representative of a typical event in the SGRB population, hinting at a quasi-universal jet structure. The diversity in the SGRB afterglow population could therefore be ascribed in large part to extrinsic (redshift, density of the surrounding medium, viewing angle) rather than intrinsic properties. Although more uncertain, the comparison can be extended to the prompt emission properties, leading to similar conclusions., Comment: 9 pages, 3 figures, submitted to A&A. Comments and suggestions are welcome

We review the physics potential of a next generation search for solar axions: the International Axion Observatory (IAXO). Endowed with a sensitivity to discover axion-like particles (ALPs) with a coupling to photons as small as $g_{a\gamma}\sim 10^{-12}$ GeV$^{-1}$, or to electrons $g_{ae}\sim$10$^{-13}$, IAXO has the potential to find the QCD axion in the 1 meV$\sim$1 eV mass range where it solves the strong CP problem, can account for the cold dark matter of the Universe and be responsible for the anomalous cooling observed in a number of stellar systems. At the same time, IAXO will have enough sensitivity to detect lower mass axions invoked to explain: 1) the origin of the anomalous "transparency" of the Universe to gamma-rays, 2) the observed soft X-ray excess from galaxy clusters or 3) some inflationary models. In addition, we review string theory axions with parameters accessible by IAXO and discuss their potential role in cosmology as Dark Matter and Dark Radiation as well as their connections to the above mentioned conundrums., Comment: Prepared for submission to JCAP; new version with very minor corrections after referee review

We update a flux-limited complete sample of Swift-based SGRBs (SBAT4, D'Avanzo et al. 2014), bringing it to 25 events and doubling its previous redshift range. We then evaluate the column densities of the events in the updated sample, in order to compare them with the NH distribution of LGRBs, using the sample BAT6ext (Arcodia et al. 2016). We rely on Monte Carlo simulations of the two populations and compare the computed NH distributions with a two sample Kolmogorov Smirnov (K-S) test. We then study how the K-S probability varies with respect to the redshift range we consider. We find that the K-S probability keeps decreasing as redshift increases until at z$\sim$1.8 the probability that short and long GRBs come from the same parent distribution drops below 1$\%$. This testifies for an observational difference among the two populations. This difference may be due to the presence of highly absorbed LGRBs above z$\sim$1.3, which have not been observed in the SGRB sample yet, although this may be due to our inability to detect them, or to the relatively small sample size., Comment: 7 pages, 3 figures, 1 table. Accepted for publication on A&A

The long-lasting tension between the observed spectra of gamma ray bursts (GRBs) and the predicted synchrotron emission spectrum might be solved if electrons do not completely cool. Evidence for incomplete cooling was recently found in Swift GRBs with prompt observations down to 0.1 keV and in one bright Fermi burst, GRB 160625B. Here we systematically search for evidence of incomplete cooling in the spectra of the ten brightest short and long GRBs observed by Fermi. We find that in 8/10 long GRBs there is compelling evidence of a low energy break (below the peak energy) and good agreement with the photon indices of the synchrotron spectrum (respectively -2/3 and -3/2 below the break and between the break and the peak energy). Interestingly, none of the ten short GRBs analysed shows a break but the low energy spectral slope is consistent with -2/3. In a standard scenario, these results imply a very low magnetic field in the emission region (B' ~ 10 G in the comoving frame), at odd with expectations., Comment: 14 pages, 15 figures, in press, accepted for publication in A&A

GRB 190114C is the first gamma-ray burst detected at Very High Energies (VHE, i.e. >300 GeV) by the MAGIC Cherenkov telescope. The analysis of the emission detected by the Fermi satellite at lower energies, in the 10 keV -- 100 GeV energy range, up to ~ 50 seconds (i.e. before the MAGIC detection) can hold valuable information. We analyze the spectral evolution of the emission of GRB 190114C as detected by the Fermi Gamma-Ray Burst Monitor (GBM) in the 10 keV -- 40 MeV energy range up to ~60 sec. The first 4 s of the burst feature a typical prompt emission spectrum, which can be fit by a smoothly broken power-law function with typical parameters. Starting on ~4 s post-trigger, we find an additional nonthermal component, which can be fit by a power law. This component rises and decays quickly. The 10 keV -- 40 MeV flux of the power-law component peaks at ~ 6 s; it reaches a value of 1.7e-5 erg cm-2 s-1. The time of the peak coincides with the emission peak detected by the Large Area Telescope (LAT) on board Fermi. The power-law spectral slope that we find in the GBM data is remarkably similar to that of the LAT spectrum, and the GBM+LAT spectral energy distribution seems to be consistent with a single component. This suggests that the LAT emission and the power-law component that we find in the GBM data belong to the same emission component, which we interpret as due to the afterglow of the burst. The onset time allows us to estimate the initial jet bulk Lorentz factor Gamma_0 is about 500, depending on the assumed circum-burst density., Comment: 7 pages, 2 figures, in press, accepted for publication in A&A

The binary neutron star merger event GW170817 was detected through both electromagnetic radiation and gravitational waves. Its afterglow emission may have been produced by either a narrow relativistic jet or an isotropic outflow. High spatial resolution measurements of the source size and displacement can discriminate between these scenarios. We present Very Long Baseline Interferometry observations, performed 207.4 days after the merger, using a global network of 32 radio telescopes. The apparent source size is constrained to be smaller than 2.5 milliarcseconds at the 90% confidence level. This excludes the isotropic outflow scenario, which would have produced a larger apparent size, indicating that GW170817 produced a structured relativistic jet. Our rate calculations show that at least 10% of neutron star mergers produce such a jet., Comment: 28 pages, 10 figures, 2 tables. This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science , (2019-02-21), doi: 10.1126/science.aau8815

e-ASTROGAM is a concept for a breakthrough observatory space mission carrying a gamma-ray telescope dedicated to the study of the non-thermal Universe in the photon energy range from 0.15 MeV to 3 GeV. The lower energy limit can be pushed down to energies as low as 30 keV for gamma-ray burst detection with the calorimeter. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with remarkable polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous and current generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will be a major player of the multiwavelength, multimessenger time-domain astronomy of the 2030s, and provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LISA, LIGO, Virgo, KAGRA, the Einstein Telescope and the Cosmic Explorer, IceCube-Gen2 and KM3NeT, SKA, ALMA, JWST, E-ELT, LSST, Athena, and the Cherenkov Telescope Array., Comment: 15 pages, 7 figures. Submitted to the proceedings of the conference SPIE Astronomical Telescopes and Instrumentation 2018: Ultraviolet to Gamma Ray. v2: corrections of authors' affiliations

The role played by the intergalactic medium (IGM) in the X-ray absorption towards high-redshift sources has recently drawn more attention in spectral analysis studies. Here, we study the X-ray absorption towards 15 flat-spectrum radio quasars at $z>2$, relying on high counting statistic ($\gtrsim10\,000$ photons) provided by XMM-Newton, with additional NuSTAR (and simultaneous Swift-XRT) observations when available. Blazars can be confidently considered to have negligible X-ray absorption along the line of sight within the host galaxy, likely swept by the kpc-scale relativistic jet. This makes our sources ideal for testing the absorption component along the IGM. Our new approach is to revisit the origin of the soft X-ray spectral hardening observed in high-z blazars in terms of X-ray absorption occurring along the IGM, with the help of a low-z sample used as comparison. We verify that the presence of absorption in excess of the Galactic value is the preferred explanation to explain the observed hardening, while intrinsic energy breaks, predicted by blazars' emission models, can easily occur out of the observing energy band in most sources. First, we perform an indirect analysis comparing the inferred amount of absorption in excess of the Galactic value with a simulated IGM absorption contribution, that increases with redshift and includes both a minimum component from diffuse IGM metals, and the additional contribution of discrete denser intervening regions. Then, we directly investigate the warm-hot IGM with a spectral model on the best candidates of our sample, obtaining an average IGM density of $n_0=1.01^{+0.53}_{-0.72}\times10^{-7}$ cm$^{-3}$ and temperature of $\log(T/\text{K})=6.45^{+0.51}_{-2.12}$. A more dedicated study is currently beyond reach, but our results can be used as a stepping stone for future more accurate analysis, involving Athena., Comment: 32 pages, 7 figures. Accepted, to be published in A&A

On May 15th, 2017, the \emph{FERMI}/LAT gamma-ray telescope observed a transient source not present in any previous high-energy catalogue: J1544-0649. It was visible for two consecutive weeks, with a flux peak on May 21st. Subsequently observed by a \emph{Swift}/XRT follow-up starting on May 26, the X-ray counterpart position was coincident with the optical transient ASASSN-17gs = AT2017egv, detected on May 25, with a potential host galaxy at $z$=0.171. We conducted a 4-months follow-up in radio (Effelsberg-100m) and optical (San Pedro M\'artir, 2.1m) bands, in order to build the overall Spectral Energy Distribution (SED) of this object. The radio data from 5 to 15 GHz confirmed the flat spectrum of the source, favoring a line of sight close to jet axis, not showing significant variability in the explored post-burst time-window. The Rx ratio, common indicator of radio loudness, gives a value at the border between the radio-loud and radio-quiet AGN populations. The Ca$_{\rm{II}}$ H\&K break value (0.29$\pm$0.05) is compatible with the range expected for the long-sought intermediate population between BL Lacs and FRI radio galaxies. An overall SED fitting from Radio to $\gamma$-ray band shows properties typical of a low-power BL Lac. As a whole, these results suggest that this transient could well be a new example of the recently discovered class of radio-weak BL Lac, showing for the first time a flare in the gamma/X-ray bands., Comment: Accepted for publication in ApJ Letters

The prompt emission of Gamma Ray Bursts extends from the early pulses observed in gamma-rays (>15 keV) to very late flares of X-ray photons (0.3-10 keV). The duration of prompt gamma-ray pulses is rather constant while the width of X-ray flares correlates with their peak time suggesting a possible different origin. However, pulses and flares have similar spectral properties. Considering internal and external shock scenarios, we derive how the energy and duration of pulses scale with their time of occurrence and we compare with observations. The absence of an observed correlation between prompt emission pulse duration and its time of occurrence favours an "internal" origin and confirms the earlier results of Ramirez-Ruiz & Fenimore. We show that also the energetic and temporal properties of X-ray flares are consistent with being produced by internal shocks between slow fireballs with a small contrast between their bulk Lorentz factors. These results relax the requirement of a long lasting central engine to explain the latest X-ray flares., Comment: 8 pages, 5 figures, 1 table, submitted to Astronomy & Astrophysics

We report our observation of the short GRB 170817A, associated to the binary neutron star merger event GW 170817, perfomed in the X-ray band with XMM-Newton 135 d after the event (on the 29th December 2017). We find evidence for a flattening of the X-ray light curve with respect to the previously observed brightening. This is also supported by a nearly simultaneous optical Hubble Space Telescope and successive X-ray Chandra and low-frequency radio observations recently reported in the literature. Since the optical-to-X-ray spectral slope did not change with respect to previous observations, we exclude that the change in the temporal evolution of the light curve is due to the passage of the cooling frequency: its origin must be geometric or dynamical. We interpret all the existing afterglow data with two models: i) a structured jet and ii) a jet-less isotropic fireball with some stratification in its radial velocity structure. Both models fit the data and predict that the radio flux must decrease simultaneously with the optical and the X-ray one, making hard to distinguish between them at the present stage. Polarimetric measures and the rate of short GRB-GW association in future LIGO/Virgo runs will be key to disentangle these two geometrically different scenarios., Comment: Accepted as a Letter in Astronomy & Astrophysics

Knowledge of the bulk Lorentz factor $\Gamma_{0}$ of GRBs allows us to compute their comoving frame properties shedding light on their physics. Upon collisions with the circumburst matter, the fireball of a GRB starts to decelerate, producing a peak or a break (depending on the circumburst density profile) in the light curve of the afterglow. Considering all bursts with known redshift and with an early coverage of their emission, we find 67 GRBs with a peak in their optical or GeV light curves at a time $t_{\rm p}$. For another 106 GRBs we set an upper limit $t_{\rm p}^{\rm UL}$. We show that $t_{\rm p}$ is due to the dynamics of the fireball deceleration and not to the passage of a characteristic frequency of the synchrotron spectrum across the optical band. Considering the $t_{\rm p}$ of 66 long GRBs and the 85 most constraining upper limits, using censored data analysis methods, we reconstruct the most likely distribution of $t_{\rm p}$. All $t_{\rm p}$ are larger than the time $t_{\rm p,g}$ when the prompt emission peaks, and are much larger than the time $t_{\rm ph}$ when the fireball becomes transparent. The reconstructed distribution of $\Gamma_0$ has median value $\sim$300 (150) for a uniform (wind) circumburst density profile. In the comoving frame, long GRBs have typical isotropic energy, luminosity, and peak energy $\langle E_{\rm iso}\rangle=3(8)\times 10^{50}$ erg, $\langle L_{\rm iso}\rangle=3(15) \times 10^{47}$ erg s$^{-1}$ , and $\langle E_{\rm peak}\rangle =1(2)$ keV in the homogeneous (wind) case. We confirm that the significant correlations between $\Gamma$ and the rest frame isotropic energy ($E_{\rm iso}$), luminosity ($L_{\rm iso}$) and peak energy ($E_{\rm peak}$) are not due to selection effects. Assuming a typical opening angle of 5 degrees, we derive the distribution of the jet baryon loading which is centered around a few $10^{-6} {\rm M_{\odot}}$., Comment: 19 pages, 11 figures, 6 tables. Accepted for publication on Astronomy & Astrophysics

Hard-TeV BL Lacs are a new type of blazars characterized by a hard intrinsic TeV spectrum, locating the peak of their gamma-ray emission in the spectral energy distribution (SED) above 2-10 TeV. Such high energies are problematic for the Compton emission, using a standard one-zone leptonic model. We study six examples of this new type of BL Lacs in the hard X-ray band with the NuSTAR satellite. Together with simultaneous observations with the SWIFT satellite, we fully constrain the peak of the synchrotron emission in their SED, and test the leptonic synchrotron self-Compton (SSC) model. We confirm the extreme nature of 5 objects also in the synchrotron emission. We do not find evidence of additional emission components in the hard X-ray band. We find that a one-zone SSC model can in principle reproduce the extreme properties of both peaks in the SED, from X-ray up to TeV energies, but at the cost of i) extreme electron energies with very low radiative efficiency, ii) conditions heavily out of equipartition (by 3 to 5 orders of magnitude), and iii) not accounting for the simultaneous UV data, which then should belong to a different emission component, possibly the same as the far-IR (WISE) data. We find evidence of this separation of the UV and X-ray emission in at least two objects. In any case, the TeV electrons must not "see" the UV or lower-energy photons, even if coming from different zones/populations, or the increased radiative cooling would steepen the VHE spectrum., Comment: 13 pages, 2 figures. Version accepted for publication in MNRAS. Fig. 2 corrected for a small plotting error

The sources responsible for the emission of high-energy ($\gtrsim$ 100 TeV) neutrinos detected by IceCube are still unknown. Among the possible candidates, active galactic nuclei with relativistic jets are often examined, since the outflowing plasma seems to offer the ideal environment to accelerate the required parent high-energy cosmic rays. The non-detection of single point sources or -- almost equivalently -- the absence, in the IceCube events, of multiplets originating from the same sky position, constrains the cosmic density and the neutrino output of these sources, pointing to a numerous population of faint sources. Here we explore the possibility that FR0 radiogalaxies, the population of compact sources recently identified in large radio and optical surveys and representing the bulk of radio-loud AGN population, can represent suitable candidates for neutrino emission. Modeling the spectral energy distribution of a FR0 radiogalaxy recently associated to a $\gamma$-ray source detected by the Large Area Telescope onboard Fermi, we derive the physical parameters of its jet, in particular the power carried by it. We consider the possible mechanisms of neutrino production, concluding that $p\gamma$ reactions in the jet between protons and ambient radiation is too inefficient to sustain the required output. We propose an alternative scenario, in which protons, accelerated in the jet, escape from it and diffuse in the host galaxy, producing neutrinos as a result of $pp$ scattering with the interstellar gas, in strict analogy with the processes taking place in star-forming galaxies., Comment: 5 pages, 3 figures, accepted for publication in MNRAS

We present time resolved spectral analysis of prompt emission from GRB 160625B, one of the brightest bursts ever detected by Fermi in its nine years of operations. Standard empirical functions fail to provide an acceptable fit to the GBM spectral data, which instead require the addition of a low-energy break to the fitting function. We introduce a new fitting function, called 2SBPL, consisting of three smoothly connected power laws. Fitting this model to the data, the goodness of the fits significantly improves and the spectral parameters are well constrained. We also test a spectral model that combines non-thermal and thermal (black body) components, but find that the 2SBPL model is systematically favoured. The spectral evolution shows that the spectral break is located around $E_{\rm break}\sim$ 100 keV, while the usual $\nu F_{\nu}$ peak energy feature $E_{\rm peak}$ evolves in the 0.5-6 MeV energy range. The slopes below and above $E_{\rm break}$ are consistent with the values -0.67 and -1.5, respectively, expected from synchrotron emission produced by a relativistic electron population with a low energy cut-off. If $E_{\rm break}$ is interpreted as the synchrotron cooling frequency, the implied magnetic field in the emitting region is $\sim$ 10 Gauss, i.e. orders of magnitudes smaller than the value expected for a dissipation region located at $\sim 10^{13-14}$ cm from the central engine. The low ratio between $E_{\rm peak}$ and $E_{\rm break}$ implies that the radiative cooling is incomplete, contrary to what is expected in strongly magnetized and compact emitting regions., Comment: 11 pages, 11 figures, 3 tables. Accepted for publication in A&A

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA., Comment: Published on Journal of High Energy Astrophysics (Elsevier)

The delay in the arrival times between high and low energy photons from cosmic sources can be used to test the violation of the Lorentz invariance (LIV), predicted by some quantum gravity theories, and to constrain its characteristic energy scale ${\rm E_{QG}}$ that is of the order of the Planck energy. Gamma-ray bursts (GRBs) and blazars are ideal for this purpose thanks to their broad spectral energy distribution and cosmological distances: at first order approximation, the constraints on ${\rm E_{QG}}$ are proportional to the photon energy separation and the distance of the source. However, the LIV tiny contribution to the total time delay can be dominated by intrinsic delays related to the physics of the sources: long GRBs typically show a delay between high and low energy photons related to their spectral evolution (spectral lag). Short GRBs have null intrinsic spectral lags and are therefore an ideal tool to measure any LIV effect. We considered a sample of $15$ short GRBs with known redshift observed by Swift and we estimate a limit on ${\rm E_{QG}}\gtrsim 1.5\times 10^{16}$ GeV. Our estimate represents an improvement with respect to the limit obtained with a larger (double) sample of long GRBs and is more robust than the estimates on single events because it accounts for the intrinsic delay in a statistical sense., Comment: 6 pages, 2 figures, 1 table, accepted for publication on Astronomy and Astrophysics

We report the results of deep optical follow-up surveys of the first two gravitational-wave sources, GW150914 and GW151226, done by the GRAvitational Wave Inaf TeAm Collaboration (GRAWITA). The VLT Survey Telescope (VST) responded promptly to the gravitational-wave alerts sent by the LIGO and Virgo Collaborations, monitoring a region of $90$ deg$^2$ and $72$ deg$^2$ for GW150914 and GW151226, respectively, and repeated the observations over nearly two months. Both surveys reached an average limiting magnitude of about 21 in the $r-$band. The paper describes the VST observational strategy and two independent procedures developed to search for transient counterpart candidates in multi-epoch VST images. Several transients have been discovered but no candidates are recognized to be related to the gravitational-wave events. Interestingly, among many contaminant supernovae, we find a possible correlation between the supernova VSTJ57.77559-59.13990 and GRB150827A detected by {\it Fermi}-GBM. The detection efficiency of VST observations for different types of electromagnetic counterparts of gravitational-wave events are evaluated for the present and future follow-up surveys., Comment: 26 pages, 11 figures- Submitted to MNRAS

The merger of two neutron stars is predicted to give rise to three major detectable phenomena: a short burst of gamma-rays, a gravitational wave signal, and a transient optical/near-infrared source powered by the synthesis of large amounts of very heavy elements via rapid neutron capture (the r-process). Such transients, named "macronovae" or "kilonovae", are believed to be centres of production of rare elements such as gold and platinum. The most compelling evidence so far for a kilonova was a very faint near-infrared rebrightening in the afterglow of a short gamma-ray burst at z = 0.356, although findings indicating bluer events have been reported. Here we report the spectral identification and describe the physical properties of a bright kilonova associated with the gravitational wave source GW 170817 and gamma-ray burst GRB 170817A associated with a galaxy at a distance of 40 Mpc from Earth. Using a series of spectra from ground-based observatories covering the wavelength range from the ultraviolet to the near-infrared, we find that the kilonova is characterized by rapidly expanding ejecta with spectral features similar to those predicted by current models. The ejecta is optically thick early on, with a velocity of about 0.2 times light speed, and reaches a radius of about 50 astronomical units in only 1.5 days. As the ejecta expands, broad absorption-like lines appear on the spectral continuum indicating atomic species produced by nucleosynthesis that occurs in the post-merger fast-moving dynamical ejecta and in two slower (0.05 times light speed) wind regions. Comparison with spectral models suggests that the merger ejected 0.03-0.05 solar masses of material, including high-opacity lanthanides., Comment: version accepted for publication in Nature. Some minor changes are expected with respect to the journal version

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA)., Comment: Accepted for publication in Advances in Space Research. Partly based on the proposal submitted on October 2016 in response to the ESA Call for next M5 mission, with expanded and updated science sections

Blazars are often characterized by a spectral break at soft X-rays, whose origin is still debated. While most sources show a flattening, some exhibit a blackbody-like soft excess with temperatures of the order of $\sim$0.1 keV, similar to low-luminosity, non-jetted Seyferts. Here we present the analysis of the simultaneous XMM-Newton and NuSTAR observation of the luminous FSRQ 4C+25.05 ($z=2.368$). The observed 0.3-30 keV spectrum is best described by the sum of a hard X-ray power law ($\Gamma = 1.38_{-0.03}^{+0.05}$) and a soft component, approximated by a blackbody with $kT_{\rm BB} = 0.66_{-0.04}^{+0.05}$ keV (rest frame). If the spectrum of 4C+25.05 is interpreted in the context of bulk Comptonization by cold electrons of broad-line region photons emitted in the direction of the jet, such an unusual temperature implies a bulk Lorentz factor of the jet of $\Gamma_{\rm bulk}\sim 11.7$. Bulk Comptonization is expected to be ubiquitous on physical grounds, yet no clear signature of it has been found so far, possibly due to its transient nature and the lack of high-quality, broad-band X-ray spectra., Comment: 6 pages, 4 figures, 2 tables. Accepted for publication in MNRASL

A brief Chandra observation of the ultraluminous quasar, SDSS J010013.02+280225.8 at redshift 6.326, showed it to be a relatively bright, soft X-ray source with a count rate of about 1 ct/ks. In this paper we present results for the quasar from a 65ks XMM-Newton observation, which well constrains its spectral shape. The quasar is clearly detected with a total of $\sim$ 460 net counts in the 0.2-10 keV band. The spectrum is characterised by a simple power-law model with photon index of $\Gamma = 2.30^{+0.10}_{-0.10}$, and the intrinsic 2-10 keV luminosity is $3.14\times10^{45}$ erg $\text{s}^{-1}$. The 1 $\sigma$ upper limit to any intrinsic absorption column density is $N_{H} = 6.07\times 10^{22} {\text{cm}}^{-2}$. No significant iron emission lines were detected. We derive the X-ray-to-optical flux ratio $\alpha_{\text{ox}}$ of $-1.74\pm$0.01, consistent with the values found in other quasars of comparable ultraviolet luminosity. We did not detect significant flux variations either in the XMM-Newton exposure or between XMM-Newton and XMM-Newton observations, which are separated by $\sim$ 8 months. The X-ray observation enables the bolometric luminosity to be calculated after modelling the spectral energy distribution: the accretion rate is found to be sub-Eddington., Comment: 6 pages, 6 figures, accepted by MNRAS

PBC J2333.9-2343 is a giant radio galaxy which shows different characteristics at different wavebands that are difficult to explain within the actual generic schemes of unification of active galactic nuclei (AGN), thus being a good candidate to host different phases of nuclear activity. We aim at disentangling the nature of this AGN by using simultaneous multiwavelength data. We obtained data in 2015 from the Very Long Baseline Array (VLBA), the San Pedro M\'{a}rtir telescope, and the XMM-Newton}observatories. This allows the study of the nuclear parts of the galaxy through its morphology and spectra, as well as the analysis of the spectral energy distribution (SED). We also reanalysed optical data from the San Pedro M\'{a}rtir telescope from 2009 previously presented in Parisi et al. (2012) for a homogeneous comparison. At X-ray frequencies, the source is unabsorbed. The optical spectra are of a type 1.9 AGN, both in 2009 and 2015, although showing a broader component in 2015. The VLBA radio images show an inverted spectrum with self-absorbed, optically thick compact core and steep spectrum, optically thin jet. The SED resembles that of typical blazars and is best represented by an external Compton (EC) model with a viewing angle of $\sim$ 3-6 degrees. The apparent size of the large scale structure of PBC J2333.9-2343 must correspond to an intrinsic deprojected value of $\sim 7$ Mpc for $\theta_v<10^\circ$, and to $> 13$ Mpc for $\theta_v<5^\circ$, a value much larger than the bigger giant radio galaxy known (4.5 Mpc). The above arguments suggest that PBC J2333.9-2343 has undergone through a new episode of nuclear activity and that the direction of the new jet has changed, in the plane of the sky and is now pointing towards us, making this source from being a radio galaxy to become a blazar, a very exceptional case of restarting activity., Comment: 11 pages, 6 figures, 5 tables, accepted for publication in A&A

We revisit the blazar sequence exploiting the complete, flux limited sample of blazars with known redshift detected by the Fermi satellite after 4 years of operations (the 3LAC sample). We divide the sources into gamma-ray luminosity bins, collect all the archival data for all blazars, and construct their spectral energy distribution (SED). We describe the average SED of blazars in the same luminosity bin through a very simple, completely phenomenological function consisting of two broken power laws connecting with a power law of fixed slope describing the radio emission. We do that separately for BL Lacs and for flat spectrum radio quasars (FSRQs) and also for all blazars together. The main results are: i) FSRQs display approximately the same SED as the luminosity increases, except for the fact that the relative importance of the high energy peak increases; ii) as a consequence, X-ray spectra of FSRQs become harder for larger luminosities; iii) BL Lacs form indeed a sequence: they become redder (i.e. the peak frequencies becomes smaller) for increasing luminosities, with a steeper gamma-ray slope and a larger dominance of the high energy peak; iv) for all blazars (BL Lacs+FSRQs) these properties becomes more prominent, as the highest luminosity bin is populated mostly by FSRQs and the lowest luminosity bin mostly by BL Lacs. This agrees with the original blazar sequence, although BL Lacs never have an average gamma-ray slope as hard as found in the original sequence. v) At high luminosities, a large fraction of FSRQs shows signs of thermal emission from the accretion disc, contributing in the optical-UV., Comment: 13 pages, 11 figures, accepted for publication in MNRAS

We investigate the global properties of the radiation emitted by the accretion disk around Kerr black holes. Using the Kerr blackbody (KERRBB) numerical model, we build an analytic approximation of the disk emission features focusing on the pattern of the produced radiation as a function of the black hole spin, mass, accretion rate and viewing angle. The assumption of having a geometrically thin disk limits our analysis to systems emitting below ~0.3 of the Eddington luminosity. We apply this analytical model to four blazars (whose jets are pointing at us) at high redshift that show clear signatures of disk emission. For them, we derive the black hole masses as a function of spin. If these jetted sources are powered by the black hole rotation, they must have high spin values, further constraining their masses., Comment: 13 pages, 14 figures, 6 tables, accepted for publications in A&A

Blazars are the most luminous and variable AGNs, and thus excellent probes of accretion and emission processes close to the central engine. We focus on PKS 1510-089 ($z=0.36$), one of the brightest gamma-ray sources in the Fermi LAT catalog, to study its complex multi-wavelength variability. PKS 1510-089 was observed twice in hard X-rays with the IBIS instrument onboard INTEGRAL during the flares of Jan 2009 and Jan 2010, and simultaneously with Swift and NOT, in addition to the constant Fermi monitoring. The optical polarization was measured in several bands on 18 Jan 2010 at the NOT. Using our and archival data we constructed historical light curves at gamma-to-radio wavelengths covering nearly 20 years and applied variability tests. We assembled SEDs in 2009 and 2010 and compared them with those at two previous epochs and with a model based on synchrotron and inverse Compton (IC) radiation. The SED modeling suggests that the physical quantities that undergo the largest variations are the total power injected into the emitting region and the random Lorentz factor of the electron distribution cooling break, that are higher in the higher gamma-ray states. This suggests a correlation of the injected power with enhanced activity of the acceleration mechanism. The cooling likely takes place at a much smaller distance ($\sim$0.03 pc) than the BLR radius. The emission at a few hundred GeV can be reproduced with IC scattering of highly relativistic electrons off FIR photons at $\sim$0.2 pc, presumably in a dusty torus. DCF analysis between the long-term optical and gamma-ray light curves yields a good correlation with no measurable delay. Our time analysis of the RXTE PCA and Fermi LAT light curves reveals no obvious (quasi-)periodicities, up to the maximum time scale (a few years) probed by the light curves, which are severely affected by red noise., Comment: 23 pages, 7 figures, 6 tables, accepted for publication in Astronomy & Astrophysics

We present QSFit, a new software package to automatically perform the analysis of Active Galactic Nuclei (AGN) optical spectra. The software provides luminosity estimates for the AGN continuum, the Balmer continuum, both optical and UV iron blended complex, host galaxy and emission lines, as well as width, velocity offset and equivalent width of 20 emission lines. Improving on a number of previous studies on AGN spectral analysis, QSFit fits all the components simultaneously, using an AGN continuum model which extends over the entire available spectrum, and is thus a probe of the actual AGN continuum whose estimates are scarcely influenced by localized features (e.g. emission lines) in the spectrum. We used QSFit to analyze 71,251 optical spectra of Type 1 AGN at z < 2 (obtained by the Sloan Digital Sky Survey, SDSS) and to produce a publicly available catalog of AGN spectral properties. Such catalog allowed us (for the first time) to estimate the AGN continuum slope and the Balmer continuum luminosity on a very large sample, and to show that there is no evident correlation between these quantities the redshift. All data in the catalog, the plots with best fitting model and residuals, and the IDL code we used to perform the analysis, are available on a dedicated website. The whole fitting process is customizable for specific needs, and can be extended to analyze spectra from other data sources. The ultimate purpose of QSFit is to allow astronomers to run standardized recipes to analyze the AGN data, in a simple, replicable and shareable way., Comment: 33 pages, 25 figures, 3 appendices. Replaced with accepted version in MNRAS. Reference website is: http://qsfit.inaf.it/

We derive the luminosity function and redshift distribution of short Gamma Ray Bursts (SGRBs) using (i) all the available observer-frame constraints (i.e. peak flux, fluence, peak energy and duration distributions) of the large population of Fermi SGRBs and (ii) the rest-frame properties of a complete sample of Swift SGRBs. We show that a steep $\phi(L)\propto L^{-a}$ with a>2.0 is excluded if the full set of constraints is considered. We implement a Monte Carlo Markov Chain method to derive the $\phi(L)$ and $\psi(z)$ functions assuming intrinsic Ep-Liso and Ep-Eiso correlations or independent distributions of intrinsic peak energy, luminosity and duration. To make our results independent from assumptions on the progenitor (NS-NS binary mergers or other channels) and from uncertainties on the star formation history, we assume a parametric form for the redshift distribution of SGRBs. We find that a relatively flat luminosity function with slope ~0.5 below a characteristic break luminosity ~3$\times10^{52}$ erg/s and a redshift distribution of SGRBs peaking at z~1.5-2 satisfy all our constraints. These results hold also if no Ep-Liso and Ep-Eiso correlations are assumed. We estimate that, within ~200 Mpc (i.e. the design aLIGO range for the detection of GW produced by NS-NS merger events), 0.007-0.03 SGRBs yr$^{-1}$ should be detectable as gamma-ray events. Assuming current estimates of NS-NS merger rates and that all NS-NS mergers lead to a SGRB event, we derive a conservative estimate of the average opening angle of SGRBs: $\theta_{jet}$~3-6 deg. Our luminosity function implies an average luminosity L~1.5$\times 10^{52}$ erg/s, nearly two orders of magnitude higher than previous findings, which greatly enhances the chance of observing SGRB "orphan" afterglows. Efforts should go in the direction of finding and identifying such orphan afterglows as counterparts of GW events., Comment: 13 pages, 5 figures, 2 tables. Accepted for publication in Astronomy & Astrophysics. Figure 5 and angle ranges corrected in revised version

We present a multiwavelength study of GRB 151027A. This is the 999th GRB detected by the Swift satellite and it has a densely sampled emission in the X-ray and optical band and has been observed and detected in the radio up to 140 days after the prompt. The multiwavelength light curve from 500 s to 140 days can be modelled through a standard forward shock afterglow but requires an additional component to reproduce the early X-ray and optical emission. We present TNG and LBT optical observations performed 19.6, 33.9 and 92.3 days after the trigger which show a bump with respect to a standard afterglow flux decay and are possibly interpreted as due to the underlying SN and host galaxy (of 0.4 uJy in the R band). Radio observations, performed with SRT, Medicina, EVN and VLBA between day 4 and 140, suggest that the burst exploded in an environment characterised by a density profile scaling with the distance from the source (wind profile). A remarkable feature of the prompt emission is the presence of a bright flare 100 s after the trigger, lasting 70 seconds in the soft X-ray band, which was simultaneously detected from the optical band up to the MeV energy range. By combining Swift-BAT/XRT and Fermi-GBM data, the broadband (0.3-1000 keV) time resolved spectral analysis of the flare reveals the coexistence of a non-thermal (power law) and thermal blackbody components. The BB component contributes up to 35% of the luminosity in the 0.3-1000 keV band. The gamma-ray emission observed in Swift-BAT and Fermi-GBM anticipates and lasts less than the soft X-ray emission as observed by Swift-XRT, arguing against a Comptonization origin. The BB component could either be produced by an outflow becoming transparent or by the collision of a fast shell with a slow, heavy and optically thick fireball ejected during the quiescent time interval between the initial and later flares of the burst., Comment: 15 pages, 12 figures, 7 tables, accepted for publication A&A

MeV blazars are a sub--population of the blazar family, exhibiting larger--than--average jet powers, accretion luminosities and black hole masses. Because of their extremely hard X--ray continua, these objects are best studied in the X-ray domain. Here, we report on the discovery by the $Fermi$ Large Area Telescope and subsequent follow-up observations with $NuSTAR$, $Swift$ and GROND of a new member of the MeV blazar family: PMN J0641$-$0320. Our optical spectroscopy provides confirmation that this is a flat--spectrum radio quasar located at a redshift of $z=1.196$. Its very hard $NuSTAR$ spectrum (power--law photon index of $\sim$1 up to $\sim$80 keV) indicates that the emission is produced via inverse Compton scattering off photons coming from outside the jet.The overall spectral energy distribution of PMN J0641$-$0320 is typical of powerful blazars and by reproducing it with a simple one-zone leptonic emission model we find the emission region to be located either inside the broad line region or within the dusty torus., Comment: Submitted to ApJ, comments welcome

The haunt of high redshift BL Lacerate objects is day by day more compelling, to firmly understand their intrinsic nature and evolution. SDSS J004054.65-0915268 is, at the moment, one of the most distant BL Lac candidate at z \sim 5 (Plotkin et al 2010). We present a new optical-near IR spectrum obtained with ALFOSC-NOT with a new, custom designed dispersive grating aimed to detect broad emission lines that could disprove this classification. In the obtained spectra we do not detect any emission features and we provide an upper limit to the luminosity of the C IV broad emission line. Therefore, the nature of the object is then discussed, building the overall spectral energy distribution and fitting it with three different models. Our fits, based on the SED modeling with different possible scenarios, cannot rule out the possibility that this source is indeed a BL Lac object although, the absence of optical variability and lack of strong radio flux, they seems to suggest that the observed optical emission originate from a thermalized accretion disk., Comment: (18 pages, 5 figures) - Accepted for publication in Astronomical Journal