On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

We present the observations of Lyman continuum (LyC) emission in the afterglow spectra of GRB 191004B at $z=3.5055$, together with those of the other two previously known LyC-emitting long gamma-ray bursts (LGRB) (GRB 050908 at $z=3.3467$, and GRB 060607A at $z=3.0749$), to determine their LyC escape fraction and compare their properties. From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of $\log(N_{\rm HI}/cm^{-2})= 17.2 \pm 0.15$, and negligible extinction ($A_{\rm V}=0.03 \pm 0.02$ mag). The only metal absorption lines detected are CIV and SiIV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Ly$\alpha$ emission. From its Ly$\alpha$ emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to $1 \leq$ SFR $\leq 4.7$ M$_{\odot}\ yr^{-1}$. We fit the Ly$\alpha$ emission with a shell model and find parameters values consistent with the observed ones. The absolute LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of $0.35^{+0.10}_{-0.11}$, $0.08^{+0.05}_{-0.04}$ and $0.20^{+0.05}_{-0.05}$, respectively. We compare the LyC escape fraction of LGRBs to the values of other LyC emitters found from the literature, showing that LGRB afterglows can be powerful tools to study LyC escape for faint high-redshift star-forming galaxies. Indeed we could push LyC leakage studies to much higher absolute magnitudes. The host galaxies of the three LGRB presented here have all $M_{\rm 1600} > -19.5$ mag, with the GRB 060607A host at $M_{\rm 1600} > -16$ mag. LGRB hosts may therefore be particularly suitable for exploring the ionizing escape fraction in galaxies that are too faint or distant for conventional techniques. Furthermore the time investment is very small compared to galaxy studies. [Abridged]

We report on the results of X-ray observations of 4XMM J111816.0-324910, a transient ultra-luminous X-ray source located in the galaxy NGC 3621. This system is characterised by a transient nature and marked variability with characteristic time-scale of ~3500 s, differently from other ULXs, which in the vast majority show limited intra-observation variability. Such a behaviour is very reminiscent of the so-called heartbeats sometimes observed in the Galactic black hole binary GRS 1915+105, where the variability time-scale is ~10-1000 s. We study the spectral and timing properties of this object and find that overall, once the differences in the variability time-scales are taken into account, they match quite closely those of both GRS 1915+105, and of a number of objects showing heartbeats in their light-curves, including a confirmed neutron star and a super-massive black hole powering an active galactic nucleus. We investigate the nature of the compact object in 4XMM J111816.0-324910 by searching for typical neutron star signatures and by attempting a mass estimate based on different methods and assumptions. Based on the current available data, we are not able to unambiguously determine the nature of the accreting compact object responsible for the observed phenomenology.

During a search for X-ray transients in the XMM-Newton archive within the EXTraS project, we discovered a new X-ray source that is detected only during a ~5 min interval of a ~21 h-long observation performed on 2011 June 21 (EXMM 023135.0-603743, probability of a random Poissonian fluctuation: ~$1.4\times10^{-27}$). With dedicated follow-up observations, we found that its position is consistent with a star-forming galaxy (SFR = 1-2 $M_\odot$ yr$^{-1}$) at redshift $z=0.092\pm0.003$ ($d=435\pm15$ Mpc). At this redshift, the energy released during the transient event was $2.8\times10^{46}$ erg in the 0.3-10 keV energy band (in the source rest frame). The luminosity of the transient, together with its spectral and timing properties, make EXMM 023135.0-603743 a gripping analog to the X-ray transient associated to SN 2008D, which was discovered during a Swift/XRT observation of the nearby ($d=27$ Mpc) supernova-rich galaxy NGC 2770. We interpret the XMM-Newton event as a supernova shock break-out or an early cocoon, and show that our serendipitous discovery is compatible with the rate of core-collapse supernovae derived from optical observations and much higher than that of tidal disruption events.

We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$).

We present the first detection of an X-ray flare from an ultracool dwarf of spectral class L. The event was identified in the EXTraS database of XMM-Newton variable sources, and its optical counterpart, J0331-27, was found through a cross-match with the Dark Energy Survey Year 3 release. Next to an earlier four-photon detection of Kelu-1, J0331-27 is only the second L dwarf detected in X-rays, and much more distant than other ultracool dwarfs with X-ray detections (photometric distance of 240 pc). From an optical spectrum with the VIMOS instrument at the VLT, we determine the spectral type of J0331-27 to be L1. The X-ray flare has an energy of E_X,F ~ 2x10^33 erg, placing it in the regime of superflares. No quiescent emission is detected, and from 2.5 Msec of XMM data we derive an upper limit of L_X,qui < 10^27 erg/s. The flare peak luminosity L_X,peak = 6.3x10^29 erg/s, flare duration tau_decay ~ 2400 s, and plasma temperature (~16 MK) are similar to values observed in X-ray flares of M dwarfs. This shows that strong magnetic reconnection events and the ensuing plasma heating are still present even in objects with photospheres as cool as ~2100 K. However, the absence of any other flares above the detection threshold of E_X,F ~2.5x10^32 erg in a total of ~2.5 Ms of X-ray data yields a flare energy number distribution inconsistent with the canonical power law dN/dE ~ E^-2, suggesting that magnetic energy release in J0331-27 -- and possibly in all L dwarfs -- takes place predominantly in the form of giant flares.

We report on the search for the optical counterpart of the gravitational event GW170814, which was carried out with the VLT Survey Telescope (VST) by the GRAvitational Wave Inaf TeAm (GRAWITA). Observations started 17.5 hours after the LIGO/Virgo alert and we covered an area of 99 deg$^2$ that encloses $\sim 77\%$ and $\sim 59\%$ of the initial and refined localization probability regions, respectively. A total of six epochs were secured over nearly two months. The survey reached an average limiting magnitude of 22 AB mag in the $r-$band. After assuming the model described in Perna et al. 2019 that derives as possible optical counterpart of a BBH event a transient source declining in about one day, we have computed a survey efficiency of about $5\%$. This paper describes the VST observational strategy and the results obtained by our analysis pipelines developed to search for optical transients in multi-epoch images. We report the catalogue of the candidates with possible identifications based on light-curve fitting. We have identified two dozens of SNe, nine AGNs, one QSO. Nineteen transients characterized by a single detection were not classified. We have restricted our analysis only to the candidates that fall into the refined localization map. None out of 39 left candidates could be positively associated with GW170814. This result implies that the possible emission of optical radiation from a BBH merger had to be fainter than r $\sim$ 22 ($L_{optical}$ $\sim$ $1.4 \times 10^{42}$ erg/s) on a time interval ranging from a few hours up to two months after the GW event.

Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and strongly affected by high background. In 2018, we obtained a new, deeper (~90 ks) XMM-Newton observation that allowed us to perform a detailed characterization of the ULXs hosted in the galaxy. ULX-1 and ULX-2, the two brightest objects (Lx~(6-10)e39 erg/s), have spectra that can be described by a two-thermal component model as often found in ULXs. ULX-1 shows also one order of magnitude in flux variability on short-term timescales (hundreds to thousand ks). The other sources (ULX-3 and ULX-4) show flux changes of at least an order of magnitude, and these objects may be candidate transient ULXs although longer X-ray monitoring or further studies are required to ascribe them to the ULX population. In addition, we found a previously undetected source that might be a new candidate ULX (labelled as ULX-5) with a luminosity of ~1e39 erg/s and hard power-law spectral shape, whose nature is still unclear and for which a background Active Galactic Nucleus cannot be excluded. We discuss the properties of all the ULXs in NGC 7456 within the framework of super-Eddington accretion onto stellar mass compact objects. Although no pulsations were detected, we cannot exclude that the sources host neutron stars.

Multi-messenger astronomy received a great boost following the discovery of kilonova AT2017gfo, the optical counterpart of the gravitational wave source GW170817 associated with the short gamma-ray burst GRB 170817A. AT2017gfo was the first kilonova that could be extensively monitored in time both photometrically and spectroscopically. Previously, only few candidates have been observed against the glare of short GRB afterglows. In this work, we aim to search the fingerprints of AT2017gfo-like kilonova emissions in the optical/NIR light curves of 39 short GRBs with known redshift. For the first time, our results allow us to study separately the range of luminosity of the blue and red components of AT2017gfo-like kilonovae in short GRBs. In particular, the red component is similar in luminosity to AT2017gfo, while the blue kilonova can be more than 10 times brighter. Finally, we find further evidence to support all the claimed kilonova detections and we exclude an AT2017gfo-like kilonova in GRBs 050509B and 061201.

GRB 190114C is the first GRB for which the detection of very-high energy emission up to the TeV range has been reported. It is still unclear whether environmental properties might have contributed to the production of these very high-energy photons, or if it is solely related to the released GRB emission. The relatively low redshift of the GRB (z=0.425) allows us to study the host galaxy of this event in detail, and to potentially identify idiosyncrasies that could point to progenitor characteristics or environmental properties responsible for such a unique event. We use ultraviolet, optical, infrared and submillimetre imaging and spectroscopy obtained with HST, VLT and ALMA to obtain an extensive dataset on which the analysis of the host galaxy is based. The host system is composed of a close pair of interacting galaxies (Delta v = 50 km s^-1), both of which are well-detected by ALMA in CO(3-2). The GRB occurred within the nuclear region (~170 pc from the centre) of the less massive but more star-forming galaxy of the pair. The host is more massive (log(M/M_odot)=9.3) than average GRB hosts at that redshift and the location of the GRB is rather unique. The enhanced star-formation rate was probably triggered by tidal interactions between the two galaxies. Our ALMA observations indicate that both host galaxy and companion have a high molecular gas fraction, as has been observed before in interacting galaxy pairs. The location of the GRB within the core of an interacting galaxy with an extinguished line-of-sight is indicative of a denser environment than typically observed for GRBs and could have been crucial for the generation of the very-high-energy photons that were observed.

Modern soft X-ray observatories can yield unique insights into time domain astrophysics, and a huge amount of information is stored - and largely unexploited - in data archives. Like a treasure-hunt, the EXTraS project harvested the hitherto unexplored temporal domain information buried in the serendipitous data collected by the European Photon Imaging Camera instrument onboard the ESA XMM-Newton, in 16 years of observations. All results have been released to the scientific community, together with new software analysis tools. This paper presents the architecture of the EXTraS science gateway, that has the goal to provide the software to the scientific community through a Web based portal using the EGI Federated Cloud infrastructure. The main focus is on the light software architecture of the portal and on the technological insights for an effective use of the EGI ecosystem.

Modern soft X-ray observatories can yield unique insights into time domain astrophysics, and a huge amount of information is stored - and largely unexploited - in data archives. Like a treasure-hunt, the EXTraS project harvested the hitherto unexplored temporal domain information buried in the serendipitous data collected by the European Photon Imaging Camera instrument onboard the XMM- Newton satellite in 20 years of observations. The result is a vast catalogue, describing the temporal behaviour of hundreds of thousands of X-ray sources. But the catalogue is just a starting point because it has to be, in its turn, further analysed. During the project an education activity has been defined and run in several workshops for high school students in Italy, Germany and UK. The final goal is to engage the students, and in perspective citizen scientists, to go through the whole validation process: they look into the data and try to discover new sources, or to characterize already known sources. This paper describes how the EXTraS science gateway is used to accomplish these tasks and highlights the first discovery, a flaring X-ray source in the globular cluster NGC 6540.

The number of supernovae known to be connected with long-duration gamma-ray bursts is increasing and the link between these events is no longer exclusively found at low redshift ($z \lesssim 0.3$) but is well established also at larger distances. We present a new case of such a liaison at $z = 0.33$ between GRB\,171010A and SN\,2017htp. It is the second closest GRB with an associated supernova of only three events detected by Fermi-LAT. The supernova is one of the few higher redshift cases where spectroscopic observations were possible and shows spectral similarities with the well-studied SN\,1998bw, having produced a similar Ni mass ($M_{\rm Ni}=0.33\pm0.02 ~\rm{M_{\odot}}$) with slightly lower ejected mass ($M_{\rm ej}=4.1\pm0.7~\rm{M_{\odot}}$) and kinetic energy ($E_{\rm K} = 8.1\pm2.5 \times 10^{51} ~\rm{erg}$). The host-galaxy is bigger in size than typical GRB host galaxies, but the analysis of the region hosting the GRB revealed spectral properties typically observed in GRB hosts and showed that the progenitor of this event was located in a very bright HII region of its face-on host galaxy, at a projected distance of $\sim$ 10 kpc from its galactic centre. The star-formation rate (SFR$_{GRB} \sim$ 0.2 M$_{\odot}$~yr$^{-1}$) and metallicity (12 + log(O/H) $\sim 8.15 \pm 0.10$) of the GRB star-forming region are consistent with those of the host galaxies of previously studied GRB-SN systems.

Ultraluminous X-ray sources (ULXs) are extragalactic X-ray emitters located off-center of their host galaxy and with a luminosity in excess of a few ${10^{39}\text{ erg s}^{-1}}$, if emitted isotropically. The discovery of periodic modulation revealed that in some ULXs the accreting compact object is a neutron star, indicating luminosities substantially above their Eddington limit. The most extreme object in this respect is ${NGC 5907~ULX-1}$ (ULX1), with a peak luminosity that is 500 times its Eddington limit. During a Chandra observation to probe a low state of ULX1, we detected diffuse X-ray emission at the position of ULX1. Its diameter is $2.7 \pm 1.0$ arcsec and contains 25 photons, none below 0.8 keV. We interpret this extended structure as an expanding nebula powered by the wind of ULX1. Its diameter of about ${200\text{ pc}}$, characteristic energy of ${\sim 1.9\text{ keV}}$, and luminosity of ${\sim 2\times10^{38}\text{ erg s}^{-1}}$ imply a mechanical power of ${1.3\times10^{41}\text{ erg s}^{-1}}$ and an age ${\sim 7 \times 10^{4}\text{ yr}}$. This interpretation suggests that a genuinely super-Eddington regime can be sustained for time scales much longer than the spin-up time of the neutron star powering the system. As the mechanical power from a single ULX nebula can rival the injection rate of cosmic rays of an entire galaxy, ULX nebulae could be important cosmic ray accelerators.

We report the detections of molecular hydrogen (H$_2$), vibrationally-excited H$_2$ (H$^*_2$), and neutral atomic carbon (CI), in two new afterglow spectra of GRBs\,181020A ($z=2.938$) and 190114A ($z=3.376$), observed with X-shooter at the Very Large Telescope (VLT). Both host-galaxy absorption systems are characterized by strong damped Lyman-$\alpha$ absorbers (DLAs) and substantial amounts of molecular hydrogen with $\log N$(HI, H$_2$) = $22.20\pm 0.05,~20.40\pm 0.04$ (GRB\,181020A) and $\log N$(HI, H$_2$) = $22.15\pm 0.05,~19.44\pm 0.04$ (GRB\,190114A). The DLA metallicites, depletion levels and dust extinctions are [Zn/H] = $-1.57\pm 0.06$, [Zn/Fe] = $0.67\pm 0.03$, and $A_V = 0.27\pm 0.02$\,mag (GRB\,181020A) and [Zn/H] = $-1.23\pm 0.07$, [Zn/Fe] = $1.06\pm 0.08$, and $A_V = 0.36\pm 0.02$\,mag (GRB\,190114A). We then examine the molecular gas content of all known H$_2$-bearing GRB-DLAs and explore the physical conditions and characteristics of these systems. We confirm that H$_2$ is detected in all CI- and H$^*_2$-bearing GRB absorption systems, but that these rarer features are not necessarily detected in all GRB H$_2$ absorbers. We find that a large molecular fraction of $f_{\rm H_2} \gtrsim 10^{-3}$ is required for CI to be detected. The defining characteristic for H$^*_2$ to be present is less clear, though a large H$_2$ column density is an essential factor. We then derive the H$_2$ excitation temperatures of the molecular gas and find that they are relatively low with $T_{\rm ex} \approx 100 - 300$\,K, however, there could be evidence of warmer components populating the high-$J$ H$_2$ levels in GRBs\,181020A and 190114A. Finally, we demonstrate that the otherwise successful X-shooter GRB afterglow campaign is hampered by a significant dust bias excluding the most dust-obscured H$_2$ absorbers from identification [Abridged].

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.

The low luminosity, X-ray flaring activity, of the sub-class of high mass X-ray binaries called Supergiant Fast X-ray Transients, has been investigated using XMM-Newton public observations, taking advantage of the products made publicly available by the EXTraS project. One of the goals of EXTraS was to extract from the XMM-Newton public archive information on the aperiodic variability of all sources observed in the soft X-ray range with EPIC (0.2-12 keV). Adopting a Bayesian block decomposition of the X-ray light curves of a sample of SFXTs, we picked out 144 X-ray flares, covering a large range of soft X-ray luminosities (1e32-1e36 erg/s). We measured temporal quantities, like the rise time to and the decay time from the peak of the flares, their duration and the time interval between adjacent flares. We also estimated the peak luminosity, average accretion rate and energy release in the flares. The observed soft X-ray properties of low-luminosity flaring activity from SFXTs is in qualitative agreement with what is expected by the application of the Rayleigh-Taylor instability model in accreting plasma near the neutron star magnetosphere. In the case of rapidly rotating neutron stars, sporadic accretion from temporary discs cannot be excluded.

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.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

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.

(Abridged) Long gamma-ray bursts (LGRB) have been suggested as promising tracers of star formation owing to their association with the core-collapse of massive stars. The goal of this work is to characterise the population of host galaxies of LGRBs at 1 < z < 2, investigate the conditions in which LGRBs form at these redshifts and assess their use as tracers of star formation. We perform a spectro-photometric analysis to determine the stellar mass, star formation rate, specific star formation rate and metallicity of the complete, unbiased host galaxy sample of the Swift/BAT6 LGRB sample at 1 < z < 2. We compare the distribution of these properties to the ones of typical star-forming galaxies from the MOSDEF and COSMOS2015 Ultra Deep surveys, within the same redshift range. We find that, similarly to z < 1, LGRBs do not directly trace star formation at 1 < z < 2, and they tend to avoid high-mass, high-metallicity host galaxies. We also find evidence for an enhanced fraction of starbursts among the LGRB host sample with respect to the star-forming population of galaxies. Nonetheless we demonstrate that the driving factor ruling the LGRB efficiency is metallicity. The LGRB host distributions can be reconciled with the ones expected from galaxy surveys by imposing a metallicity upper limit of 12+logOH ~ 8.55. Metallicity rules the LGRB production efficiency, which is stifled at Z > 0.7 Zsun. Under this hypothesis we can expect LGRBs to trace star formation at z > 3, once the bulk of the star forming galaxy population are characterised by metallicities below this limit. The moderately high metallicity threshold found is in agreement with the conditions necessary to rapidly produce a fast-rotating Wolf-Rayet star a in close binary system, and could be accommodated by single star models under chemically homogeneous mixing with very rapid rotation and weak magnetic coupling.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

We present the first X-ray and UV/optical observations of a very bright and fast nova in the disc of M31, M31N 2013-01b. The nova reached a peak magnitude $R\sim$15 mag and decayed by 2 magnitudes in only 3 days, making it one of the brightest and fastest novae ever detected in Andromeda. From archival multi-band data we have been able to trace its fast evolution down to $U>21$ mag in less than two weeks and to uncover for the first time the Super-Soft X-ray phase, whose onset occurred 10-30 days from the optical maximum. The X-ray spectrum is consistent with a blackbody with a temperature of $\sim$50 eV and emitting radius of $\sim$4$\times 10^{9}$ cm, larger than a white dwarf radius, indicating an expanded region. Its peak X-ray luminosity, 3.5$\times 10^{37}$ erg s$^{-1}$, locates M31N 2013-01b among the most luminous novae in M31. We also unambiguously detect a short 1.28$\pm$0.02 h X-ray periodicity that we ascribe to the binary orbital period, possibly due to partial eclipses. This makes M31N 2013-01b the first nova in M31 with an orbital period determined. The short period also makes this nova one of the few known below the 2-3 h orbital period gap. All the observed characteristics strongly indicate that M31N 2013-01b harbours a massive white dwarf and a very low-mass companion, consistent with being a nova belonging to the disc population of the Andromeda Galaxy.

Aim. The unusual transient ATLAS17aeu was serendipitously detected within the sky localisation of the gravitational wave trigger GW170104. The importance of a possible association with gravitational waves coming from a binary black hole merger led to an extensive follow-up campaign, with the aim of assessing a possible connection with GW170104. Methods. With several telescopes, we carried out both photometric and spectroscopic observations of ATLAS17aeu, for several epochs, between $\sim 3$ and $\sim 230$ days after the first detection. Results. We studied in detail the temporal and spectroscopic properties of ATLAS17aeu and its host galaxy. We detected spectral features similar to those of a broad lined supernova superposed to an otherwise typical long-GRB afterglow. Based on analysis of the optical light curve, spectrum and host galaxy SED, we conclude that the redshift of the source is probably $z \simeq 0.5 \pm 0.2$. Conclusions. While the redshift range we have determined is marginally compatible with that of the gravitational wave event, the presence of a supernova component and the consistency of this transient with the E$_{\rm p}$-E$_{\rm iso}$ correlation support the conclusion that ATLAS17aeu was associated with the long gamma-ray burst GRB170105A. This rules out the association of the GRB170105A/ATLAS17aeu transient with the gravitational wave event GW170104, which was due to a binary black hole merger.

The gamma-ray bursts hosts (GRBHs) are excellent targets to study the extinction properties of dust and its effects on the global emission of distant galaxies. The dust extinction curve is measured along the GRB afterglow line of sight and the analysis of the spectral energy distribution (SED) of the host galaxy gives access to the global dust attenuation of the stellar light. We selected a sample of 30 GRBs for which the extinction curve along the GRB afterglow line-of-sight (l.o.s) is measured in the rest-frame ultraviolet (UV) up to optical and we analysed the properties of the extinction curve as a function of the host galaxy properties. From these 30 GRBs, we selected seven GRBHs with a good rest-frame UV to near infra-red spectral coverage for the host. The attenuation curve was derived by fitting the SEDs of the GRBH sample with the CIGALE SED fitting code. Different star formation histories were studied to recover the star formation rates (SFR) derived using H$_{\alpha}$ luminosities. The most extinguished GRBs are preferentially found in the more massive hosts and the UV bump is preferentially found in the most extinguished GRB l.o.s. Five out of seven hosts are best fitted with a recent burst of star formation, leading to lower stellar mass estimates than previously found. The average attenuation in the host galaxies is about 70% of the amount of extinction along the GRB l.o.s. We find a great variety in the derived attenuation curves of GRBHs, the UV slope can be similar, flatter or even steeper than the extinction curve slope. We find that the flatter (steeper) attenuation curves are found in galaxies with the highest (lowest) SFR and stellar masses. The comparison of our results with radiative transfer simulations leads to a uniform distribution of dust and stars in a very clumpy ISM for half the GRBHs and various dust-stars geometries for the second half of the sample.

Broad-line Ic supernovae (SNe Ic-BL) are a very rare class of core-collapse supernovae exhibiting high ejecta velocities and high kinetic energies. They are the only type of SNe that accompany long gamma-ray burst (GRB) explosions. Systematic differences found in the spectra of SNe Ic-BL with and without GRBs (GRB-SNe and SNe Ic-BL, respectively) suggest that either the progenitor or/and explosion mechanism of SNe Ic-BL with and without a GRB differ, or the difference could be only due to the viewing angle of the observer with respect to the orientation of the collimated explosion. We present the systematic comparison of the host galaxies of broad-lined SNe Ic with and without a detected GRB, the latter being detected in untargeted surveys, with the aim to find out whether there are any systematic differences between the environments in which these two classes of SNe preferentially explode. We study photometric properties of the host galaxies of a sample of 8 GRB-SNe and a sample of 28 SNe Ic-BL at z < 0.2. The two galaxy samples have indistinguishable luminosity and proper size distribution. We find indications that GRB-SNe on average occur closer to the centres of their host galaxies, i.e. the samples have a different distribution of projected offsets, normalized by the galaxy sizes. In addition we compare gas-phase metallicities of the GRB-SNe and SNe Ic-BL host samples and find that a larger fraction of super-solar metallicity hosts are found among the SNe Ic-BL without a GRB. Our results are indicative of a genuine difference between the two types of explosions and suggest that the viewing angle is not the main source of difference in the spectra of the two classes. We discuss the implications our results have on our understanding of progenitors of SNe Ic-BL with and without a GRB.

The UV extinction feature at 2175 \AA\ is ubiquitously observed in the Galaxy but is rarely detected at high redshifts. Here we report the spectroscopic detection of the 2175 \AA\ bump on the sightline to the \gamma-ray burst (GRB) afterglow GRB 180325A at z=2.2486, the only unambiguous detection over the past ten years of GRB follow-up, at four different epochs with the Nordic Optical Telescope (NOT) and the Very Large Telescope (VLT)/X-shooter. Additional photometric observations of the afterglow are obtained with the Gamma-Ray burst Optical and Near-Infrared Detector (GROND). We construct the near-infrared to X-ray spectral energy distributions (SEDs) at four spectroscopic epochs. The SEDs are well-described by a single power-law and an extinction law with R_V~4.4, A_V~1.5, and the 2175 \AA\ extinction feature. The bump strength and extinction curve are shallower than the average Galactic extinction curve. We determine a metallicity of [Zn/H]>-0.98 from the VLT/X-shooter spectrum. We detect strong neutral carbon associated with the GRB with an equivalent width of Wr(\lambda 1656) = 0.85+/-0.05. We also detect optical emission lines from the host galaxy. Based on the H\alpha emission line flux, the derived dust-corrected star-formation rate is ~46+/-4 M_sun/yr and the predicted stellar mass is log M*/M_sun~9.3+/-0.4, suggesting the host galaxy is amongst the main-sequence star-forming galaxies.

We report on the serendipitous discovery of a new transient in NGC 5907, at a peak luminosity of 6.4x10^{39} erg/s. The source was undetected in previous 2012 Chandra observations with a 3 sigma upper limit on the luminosity of 1.5x10^{38} erg/s, implying a flux increase of a factor of >35. We analyzed three recent 60ks/50ks Chandra and 50ks XMM-Newton observations, as well as all the available Swift observations performed between August 2017/March 2018. Until the first half of October 2017, Swift observations do not show any emission from the source. The transient entered the ULX regime in less than two weeks and its outburst was still on-going at the end of February 2018. The 0.3-10 keV spectrum is consistent with a single multicolour blackbody disc (kT~1.5 keV). The source might be a ~30 solar mass black hole accreting at the Eddington limit. However, although we did not find evidence of pulsations, we cannot rule-out the possibility that this ULX hosts an accreting neutron star.

We report the discovery of a flaring X-ray source in the globular cluster NGC 6540, obtained during the EXTraS project devoted to a systematic search for variability in archival data of the XMM-Newton satellite. The source had a quiescent X-ray luminosity of the order of ~10^32 erg/s in the 0.5-10 keV range (for a distance of NGC 6540 of 4 kpc) and showed a flare lasting about 300 s. During the flare, the X-ray luminosity increased by more than a factor 40, with a total emitted energy of ~10^36 erg. These properties, as well as Hubble Space Telescope photometry of the possible optical counterparts, suggest the identification with a chromospherically active binary. However, the flare luminosity is significantly higher than what commonly observed in stellar flares of such a short duration, leaving open the possibility of other interpretations.

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.

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.

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

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).

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}}$.

We performed a search for eclipsing and dipping sources in the archive of the EXTraS project - a systematic characterization of the temporal behaviour of XMM-Newton point sources. We discovered dips in the X-ray light curve of 3XMM J004232.1+411314, which has been recently associated with the hard X-ray source dominating the emission of M31. A systematic analysis of XMM-Newton observations revealed 13 dips in 40 observations (total exposure time $\sim$0.8 Ms). Among them, four observations show two dips, separated by $\sim$4.01 hr. Dip depths and durations are variable. The dips occur only during low-luminosity states (L$_{0.2-12}<1\times10^{38}$ erg s$^{-1}$), while the source reaches L$_{0.2-12}\sim2.8\times10^{38}$ erg s$^{-1}$. We propose this system to be a new dipping Low-Mass X-ray Binary in M31 seen at high inclination (60$^{\circ}$-80$^{\circ}$), the observed dipping periodicity is the orbital period of the system. A blue HST source within the Chandra error circle is the most likely optical counterpart of the accretion disk. The high luminosity of the system makes it the most luminous dipper known to date.

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.

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.

The merger of two dense stellar remnants including at least one neutron star (NS) is predicted to produce gravitational waves (GWs) and short duration gamma ray bursts (GRBs). In the process, neutron-rich material is ejected from the system and heavy elements are synthesized by r-process nucleosynthesis. The radioactive decay of these heavy elements produces additional transient radiation termed "kilonova" or "macronova". We report the detection of linear optical polarization P = (0.50 +/- 0.07)% at 1.46 days after detection of the GWs from GW170817, a double neutron star merger associated with an optical macronova counterpart and a short GRB. The optical emission from a macronova is expected to be characterized by a blue, rapidly decaying, component and a red, more slowly evolving, component due to material rich of heavy elements, the lanthanides. The polarization measurement was made when the macronova was still in its blue phase, during which there is an important contribution from a lanthanide-free outflow. The low degree of polarization is consistent with intrinsically unpolarized emission scattered by Galactic dust, suggesting a symmetric geometry of the emitting region and low inclination of the merger system. Stringent upper limits to the polarization degree from 2.45 - 9.48 days post-burst are consistent with the lanthanides-rich macronova interpretation.

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.

We present a comprehensive temporal and spectral analysis of the long Swift GRB 120327A afterglow data to investigate the possible causes of the observed early time colour variations. We collected data from various instruments/telescopes in different bands (X-rays, ultra- violet, optical and near-infrared) and determined the shapes of the afterglow early-time light curves. We studied the overall temporal behaviour and the spectral energy distributions from early to late times. The ultra-violet, optical, and near-infrared light curves can be modelled with a single power-law component between 200 and 2e4 s after the burst event. The X-ray light curve shows a canonical steep-shallow-steep behaviour, typical of long gamma-ray bursts. At early times a colour variation is observed in the ultra-violet/optical bands, while at very late times a hint of a re-brightening is visible. The observed early time colour change can be explained as a variation in the intrinsic optical spectral index, rather than an evolution of the optical extinction.

We report on the X-ray dust-scattering features observed around the afterglow of the gamma ray burst GRB 160623A. With an XMM-Newton observation carried out ~2 days after the burst, we found evidence of at least six rings, with angular size expanding between ~2 and 9 arcmin, as expected for X-ray scattering of the prompt GRB emission by dust clouds in our Galaxy. From the expansion rate of the rings, we measured the distances of the dust layers with extraordinary precision: 528.1 +\- 1.2 pc, 679.2 +\- 1.9 pc, 789.0 +\- 2.8 pc, 952 +\- 5 pc, 1539 +\- 20 pc and 5079 +\- 64 pc. A spectral analysis of the ring spectra, based on an appropriate dust-scattering model (BARE-GR-B from Zubko et al. 2004}) and the estimated burst fluence, allowed us to derive the column density of the individual dust layers, which are in the range 7x10^20-1.5x10^22 cm^-2. The farthest dust-layer (i.e. the one responsible for the smallest ring) is also the one with the lowest column density and it is possibly very extended, indicating a diffuse dust region. The properties derived for the six dust-layers (distance, thickness, and optical depth) are generally in good agreement with independent information on the reddening along this line of sight and on the distribution of molecular and atomic gas.

We use numerical N-body hydrodynamical simulations with varying PopIII stellar models to investigate the possibility of detecting first star signatures with observations of high-redshift damped Ly$\alpha$ absorbers (DLAs). The simulations include atomic and molecular cooling, star formation, energy feedback and metal spreading due to the evolution of stars with a range of masses and metallicities. Different initial mass functions (IMFs) and corresponding metal-dependent yields and lifetimes are adopted to model primordial stellar populations. The DLAs in the simulations are selected according to either the local gas temperature (temperature selected) or the host mass (mass selected). We find that 3\% (40\%) of mass (temperature) selected high-$z$ ($z\ge5.5$) DLAs retain signatures of pollution from PopIII stars, independently from the first star model. Such DLAs have low halo mass ($<10^{9.6}\,\rm M_{\odot}$), metallicity ($<10^{-3}\,\rm Z_{\odot}$) and star formation rate ($<10^{-1.5}\,\rm M_{\odot}\,yr^{-1}$). { Metal abundance ratios of DLAs imprinted in the spectra of QSO} can be useful tools to infer the properties of the polluting stellar generation and to constrain the first star mass ranges. Comparing the abundance ratios derived from our simulations to those observed in DLAs at $z\ge5$, we find that most of these DLAs are consistent within errors with PopII stars dominated enrichment and strongly disfavor the pollution pattern of very massive first stars (i.e. 100~$\rm M_{\odot}$-500~$\rm M_{\odot}$). However, some of them could still result from the pollution of first stars in the mass range [0.1, 100]~$\rm M_{\odot}$. In particular, we find that the abundance ratios from SDSS J1202+3235 are consistent with those expected from PopIII enrichment dominated by massive (but not extreme) first stars.

Ultraluminous x-ray sources (ULXs) in nearby galaxies shine brighter than any X-ray source in our Galaxy. ULXs are usually modeled as stellar-mass black holes (BHs) accreting at very high rates or intermediate-mass BHs. We present observations showing that NGC5907 ULX is instead an x-ray accreting neutron star (NS) with a spin period evolving from 1.43~s in 2003 to 1.13~s in 2014. It has an isotropic peak luminosity of about 1000 times the Eddington limit for a NS at 17.1~Mpc. Standard accretion models fail to explain its luminosity, even assuming beamed emission, but a strong multipolar magnetic field can describe its properties. These findings suggest that other extreme ULXs (x-ray luminosity > 10^{41} erg/s) might harbor NSs.

We investigate the existence of a metallicity threshold for the production of long gamma-ray bursts (LGRBs). We used the host galaxies of the Swift/BAT6 sample of LGRBs. We considered the stellar mass, star formation rate (SFR), and metallicity determined from the host galaxy photometry and spectroscopy up to z = 2 and used them to compare the distribution of host galaxies to that of field galaxies in the mass-metallicity and fundamental metallicity relation plane. We find that although LGRBs also form in galaxies with relatively large stellar masses, the large majority of host galaxies have metallicities below log(O=H)~8.6. The extension to z = 2 results in a good sampling of stellar masses also above Log(Mstar/Msun)~9.5 and provides evidence that LGRB host galaxies do not follow the fundamental metallicity relation. As shown by the comparison with dedicated numerical simulations of LGRB host galaxy population, these results are naturally explained by the existence of a mild (~0.7 Zsun) threshold for the LGRB formation. The present statistics does not allow us to discriminate between different shapes of the metallicity cutoff, but the relatively high metallicity threshold found in this work is somewhat in disagreement to most of the standard single-star models for LGRB progenitors.

NGC 7793 P13 is a variable (luminosity range ~100) ultraluminous X-ray source (ULX) proposed to host a stellar-mass black hole of less than 15 M$_{\odot}$ in a binary system with orbital period of 64 d and a 18-23 M$_{\odot}$ B9Ia companion. Within the EXTraS project we discovered pulsations at a period of ~0.42 s in two XMM-Newton observations of NGC 7793 P13, during which the source was detected at $L_{\mathrm{X}}\sim2.1\times10^{39}$ and $5\times10^{39}$ erg s$^{-1}$ (0.3-10 keV band). These findings unambiguously demonstrate that the compact object in NGC 7793 P13 is a neutron star accreting at super-Eddington rates. While standard accretion models face difficulties accounting for the pulsar X-ray luminosity, the presence of a multipolar magnetic field with $B$ ~ few $\times$ 10$^{13}$ G close to the base of the accretion column appears to be in agreement with the properties of the system.

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.

The EXTraS project to explore the X-ray Transient and variable Sky searches for coherent signals in the X-ray archival data of XMM-Newton. XMM-Newton performed more than 400 pointed observations in the region of the Large Magellanic Cloud (LMC). We inspected the results of the EXTraS period search to systematically look for new X-ray pulsators in our neighbour galaxy. We analysed the XMM-Newton observations of two sources from the 3XMM catalogue which show significant signals for coherent pulsations. 3XMM J051259.8-682640 was detected as source with hard X-ray spectrum in two XMM-Newton observations, revealing a periodic modulation of the X-ray flux with 956~s. As optical counterpart we identify an early-type star with Halpha emission. The OGLE I-band light curve exhibits a regular pattern with three brightness dips which mark a period of ~1350 d. The X-ray spectrum of 3XMM J051034.6-670353 is dominated by a super-soft blackbody-like emission component (kT ~ 70 eV) which is modulated by nearly 100% with a period of ~1418 s. From GROND observations we suggest a star with r' = 20.9 mag as possible counterpart of the X-ray source. 3XMM J051259.8-682640 is confirmed as a new Be/X-ray binary pulsar in the Large Magellanic Cloud. We discuss the long-term optical period as likely orbital period which would be the longest known from a high-mass X-ray binary. The spectral and temporal properties of the super-soft source 3XMM J051034.6-670353 are very similar to those of RX J0806.3+1527 and RX J1914.4+2456 suggesting that it belongs to the class of double-degenerate polars and is located in our Galaxy rather than in the LMC.

We study the possibility to detect and distinguish signatures of enrichment from PopIII stars in observations of PopII GRBs (GRBIIs) at high redshift by using numerical N-body/hydrodynamical simulations including atomic and molecular cooling, star formation and metal spreading from stellar populations with different initial mass functions (IMFs), yields and lifetimes. PopIII and PopII star formation regimes are followed simultaneously and both a top-heavy and a Salpeter-like IMF for pristine PopIII star formation are adopted. We find that the fraction of GRBIIs hosted in a medium previously enriched by PopIII stars (PopIII-dominated) is model independent. Typical abundance ratios, such as [Si/O] vs [C/O] and [Fe/C] vs [Si/C], can help to disentangle enrichment from massive and intermediate PopIII stars, while low-mass first stars are degenerate with regular PopII generations. The properties of galaxies hosting PopIII-dominated GRBIIs are not very sensitive to the particular assumption on the mass of the first stars.

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.

Long gamma-ray bursts (LGRBs) and superluminous supernovae (SLSNe) are both explosive transients with very massive progenitor stars. Clues about the nature of the progenitors can be found by investigating environments in which such transients occur. While studies of LGRB host galaxies have a long history, dedicated observational campaigns have only recently resulted in a high enough number of photometrically and spectroscopically observed SLSN hosts to allow statistically significant analysis of their properties. In this paper we make a comparison of the host galaxies of hydrogen-poor (H-poor) SLSNe and the Swift/BAT6 sample of LGRBs. In contrast to previous studies we use a complete sample of LGRBs and we address a special attention to the comparison methodology and the selection of SLSN sample whose data have been compiled from the available literature. At intermediate redshifts (0.3 < z < 0.7) the two classes of transients select galaxies whose properties (stellar mass, luminosity, star-formation rate, specific star-formation rate and metallicity) do not differ on average significantly. Moreover, the host galaxies of both classes of objects follow the fundamental metallicity relation and the fundamental plane of metallicity. In contrast to previous studies we show that at intermediate redshifts the emission line equivalent widths of the two populations are essentially the same and that the previous claims regarding the higher fraction of SLSN hosts among the extreme emission line galaxies with respect to LGRBs are mostly due to a larger fraction of strong-line emitters among SLSN hosts at z < 0.3, where samples of LGRB hosts are small and poorly defined.

We present the results of a spectroscopic search for narrow emission and absorption features in the X-ray spectra of long gamma-ray burst (GRB) afterglows. Using XMM-Newton data, both EPIC and RGS spectra, of six bright (fluence >10^{-7} erg cm^{-2}) and relatively nearby (z=0.54-1.41) GRBs, we performed a blind search for emission or absorption lines that could be related to a high cloud density or metal-rich gas in the environ close to the GRBs. We detected five emission features in four of the six GRBs with an overall statistical significance, assessed through Monte Carlo simulations, of <3.0 sigma. Most of the lines are detected around the observed energy of the oxygen edge at ~0.5 keV, suggesting that they are not related to the GRB environment but are most likely of Galactic origin. No significant absorption features were detected. A spectral fitting with a free Galactic column density (N_H) testing different models for the Galactic absorption confirms this origin because we found an indication of an excess of Galactic N_H in these four GRBs with respect to the tabulated values.

Current models suggest gamma-ray bursts could be used as a way of probing Population III stars - the first stars in the early Universe. In this paper we use numerical simulations to demonstrate that late time radio observations of gamma-ray burst afterglows could provide a means of identifying bursts that originate from Population III stars, if these were highly massive, independently from their redshift. We then present the results from a pilot study using the Australia Telescope Compact Array at 17 GHz, designed to test the hypothesis that there may be Population III gamma-ray bursts amongst the current sample of known events. We observed three candidates plus a control gamma-ray burst, and make no detections with upper limits of 20-40 uJy at 500-1300 days post explosion.

We study the X-ray absorption of a complete sample of 99 bright Swift gamma-ray bursts. Over the last few years, a strong correlation between the intrinsic X-ray absorbing column density (N_H(z)) and the redshift was found. This absorption excess in high-z GRBs is now thought to be due to the overlooked contribution of the absorption along the intergalactic medium, by means of both intervening objects and the diffuse warm-hot intergalactic medium along the line of sight. In this work we neglect the absorption along the IGM, because our purpose is to study the eventual effect of a radical change in the Galactic absorption model on the N_H(z) distribution. Therefore, we derive the intrinsic absorbing column densities using two different Galactic absorption models, the Leiden Argentine Bonn HI survey and the more recent model including molecular hydrogen. We find that, if on the one hand the new Galactic model considerably affects the single column density values, on the other hand there is no drastic change in the distribution as a whole. It becomes clear that the contribution of Galactic column densities alone, no matter how improved, is not sufficient to change the observed general trend and it has to be considered as a second-order correction. The cosmological increase of N_H(z) as a function of redshift persists and, in order to explain the observed distribution, it is necessary to include the contribution of both the diffuse intergalactic medium and the intervening systems along the line of sight of the GRBs.

Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and could thus be a potentially powerful tool to trace cosmic star formation. However, especially at low redshifts (z < 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the impact of the environment on GRB formation. We study host galaxy spectra of the Swift/BAT6 complete sample of 14 z < 1 bright LGRBs. We use the detected nebular emission lines to measure the dust extinction, star formation rate (SFR) and nebular metallicity (Z) of the hosts and supplement the data set with previously measured stellar masses M$_{\star}$. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M$_{\star}$ relations) are compared to samples of field star-forming galaxies.We find that LGRB hosts at z < 1 have on average lower SFRs than if they were direct star-formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to Z $\sim 8.4-8.5$, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are not inconsistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of increased incidence of starburst galaxies in the Swift/BAT6 z < 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevent us from reaching more quantitative conclusions.

The detection of the first gravitational wave (GW) transient GW150914 prompted an extensive campaign of follow-up observations at all wavelengths. Although no dedicated XMM-Newton observations have been performed, the satellite passed through the GW150914 error box during normal operations. Here we report the analysis of the data taken during these satellite slews performed two hours and two weeks after the GW event. Our data cover 1.1 square degrees and 4.8 square degrees of the final GW localization region. No credible X-ray counterpart to GW150914 is found down to a sensitivity of 6E-13 erg/cm2/s in the 0.2-2 keV band. Nevertheless, these observations show the great potential of XMM-Newton slew observations for the search of the electromagnetic counterparts of GW events. A series of adjacent slews performed in response to a GW trigger would take <1.5 days to cover most of the typical GW credible region. We discuss this scenario and its prospects for detecting the X-ray counterpart of future GW detections.

The power spectrum of Near InfraRed Background (NIRB) fluctuations measured at 3.6 $\mu$m by {\tt Spitzer} shows a clustering excess over the known galaxies signal that has been interpreted in terms of early ($z\simgt 13$), accreting (direct collapse) black holes (DCBH) or low-$z$ intrahalo light (IHL). In addition, these fluctuations correlate with the cosmic X-ray background (CXB) measured at (0.5-2) keV, supporting the black hole explanation. This scenario has been questioned by the recent detection of a correlation between the two {\tt CIBER} 1.1/1.6 $\mu$m bands with the 3.6 $\mu$m {\tt Spitzer} one. This correlation is hardly explained by early DCBHs that, due to intergalactic absorption, cannot contribute to the shortest wavelength bands. Here we show that the new correlation is caused instead by a Diffuse Galactic Light (DGL) component arising from Galactic stellar light scattered by dust. The black hole interpretation of the excess remains perfectly valid and, actually, the inclusion of DGL allows less demanding (by up to about 30\%) requirements on the DCBH abundance/mass.

X-ray emission from Young Stellar Objects (YSOs) is crucial to understand star formation. A very limited amount of X-ray results is available for the protostellar (ClassI) phase. A systematic search of transient X-ray phenomena combined with a careful evaluation of the evolutionary stage offer a widely unexplored window to our understanding of YSOs X-ray properties. Within the EXTraS project, a search for transients and variability in the whole XMM-Newton archive, we discover transient X-ray emission consistent with ISO-Oph 85, a strongly embedded YSO in the rho Ophiuchi region, not detected in previous time-averaged X-ray studies. We extract an X-ray light curve for the flare and determine its spectral parameters from XMM-Newton/EPIC (European Photon Imaging Camera) data using quantile analysis. The X-ray flare ($2500\,s$), the only one detected in the XMM-Newton archive for ISO-Oph 85, has a luminosity of $LogL_X[erg/s]=31.1$ and a spectrum consistent with a highly-absorbed one-component thermal model ($N_H=1.0^{+1.2}_{-0.5}10^{23}\,cm^{-2}$, $kT=1.15^{+2.35}_{-0.65}\,keV)$. We set an upper limit of $LogL_X[erg/s]<29.5$ to the quiescent X-ray luminosity. We build a SED with IR to mm photometry drawn from literature and mid-IR Spitzer and sub-mm Herschel photometry analysed by us, and compare it with pre-computed models. The sub-mm emission peak in the Herschel data suggests that the object is a ClassI protostar. However, the Herschel/IR position offset is larger than for other YSOs in the region, leaving some doubt on the association. This is the first X-ray flare from a YSO recognised as a candidate ClassI YSO via the analysis of its complete SED. This work shows how the analysis of the whole SED is fundamental for the classification of YSOs, and how the X-ray source detection techniques we developed can open a new era in time-resolved analysis of the X-ray emission from stars.

During a search for coherent signals in the X-ray archival data of XMM-Newton, we discovered a modulation at 1.2 s in 3XMM J004301.4+413017 (3X J0043), a source lying in the direction of an external arm of M 31. This short period indicates a neutron star (NS). Between 2000 and 2013, the position of 3X J0043 was imaged by public XMM-Newton observations 35 times. The analysis of these data allowed us to detect an orbital modulation at 1.27 d and study the long-term properties of the source. The emission of the pulsar was rather hard (most spectra are described by a power law with $\Gamma < 1$) and, assuming the distance to M 31, the 0.3-10 keV luminosity was variable, from $\sim$$3\times10^{37}$ to $2\times10^{38}$ erg s$^{-1}$. The analysis of optical data shows that, while 3X J0043 is likely associated to a globular cluster in M 31, a counterpart with $V\gtrsim22$ outside the cluster cannot be excluded. Considering our findings, there are two main viable scenarios for 3X J0043: a peculiar low-mass X-ray binary, similar to 4U 1822-37 or 4U 1626-67, or an intermediate-mass X-ray binary resembling Her X-1. Regardless of the exact nature of the system, 3X J0043 is the first accreting NS in M 31 in which the spin period has been detected.

The EXTraS project (Exploring the X-ray Transient and variable Sky) will characterise the temporal behaviour of the largest ever sample of objects in the soft X-ray range (0.1-12 keV) with a complex, systematic and consistent analysis of all data collected by the European Photon Imaging Camera (EPIC) instrument onboard the ESA XMM-Newton X-ray observatory since its launch. We will search for, and characterize variability (both periodic and aperiodic) in hundreds of thousands of sources spanning more than nine orders of magnitude in time scale and six orders of magnitude in flux. We will also search for fast transients, missed by standard image analysis. Our analysis will be completed by multiwavelength characterization of new discoveries and phenomenological classification of variable sources. All results and products will be made available to the community in a public archive, serving as a reference for a broad range of astrophysical investigations.

LGRBs are associated with massive stars and are therefore linked to star formation. The conditions necessary to produce LGRBs can affect the relation between the LGRB rate and star formation. By using the power of a complete LGRB sample, our aim is to understand whether such a bias exists and, if it does, what is its origin. In this first paper, we build the SED of the z<1 host galaxies of the BAT6 LGRB sample, and determine their stellar masses from SED fitting. We compare the resulting stellar mass distribution (i) with star-forming galaxies observed in deep surveys (UltraVISTA); (ii) with semi-analitical models of the z<1 star forming galaxy population and (iii) with numerical simulations of LGRB hosts having different metallicity thresholds for the progenitor star environment. We find that at z<1 LGRBs tend to avoid massive galaxies and are powerful in selecting faint low-mass star-forming galaxies. The stellar mass distribution of the hosts is not consistent with that of the UltraVISTA star-forming galaxies weighted for their SFR. This implies that, at least at z<1, LGRBs are not unbiased tracers of star formation. To make the two distributions consistent, a much steeper faint-end of the mass function would be required, or a very shallow SFR-Mass relation for the low mass galaxy population. GRB host galaxy simulations indicates that, to reproduce the stellar mass distribution, a metallicity threshold of the order of Z_th=0.3-0.5Z_sun is necessary. Models without a metallicity threshold or with an extreme threshold of Z_th = 0.1Z_sun are excluded at z<1. The use of the BAT6 complete sample makes this result not affected by possible biases which could have influenced past results based on incomplete samples. The preference for low metallicities (Z<~0.5Z_sun) can be a consequence of the particular conditions needed for the progenitor star to produce a GRB. (Abridged)

We derive, adopting a direct method, the luminosity function and the formation rate of long Gamma Ray Bursts through a complete, flux-limited, sample of Swift bursts which has a high level of completeness in redshift z (~82%). We parametrise the redshift evolution of the GRB luminosity as L = L_0(1+ z)^k and we derive k = 2.5, consistently with recent estimates. The de-evolved luminosity function of GRBs can be represented by a broken power law with slopes a = -1.32 +- 0.21 and b = -1.84 +- 0.24 below and above, respectively, a characteristic break luminosity L_0,b = 10^51.45+-0.15 erg/s. Under the hypothesis of luminosity evolution we find that the GRB formation rate increases with redshift up to z~2, where it peaks, and then decreases in agreement with the shape of the cosmic star formation rate. We test the direct method through numerical simulations and we show that if it is applied to incomplete (both in redshift and/or flux) GRB samples it can misleadingly result in an excess of the GRB formation rate at low redshifts.

High-redshift gamma-ray bursts have several advantages for the study of the distant universe, providing unique information about the structure and properties of the galaxies in which they exploded. Spectroscopic identification with large ground-based telescopes has improved our knowledge of the class of such distant events. We present the multi-wavelength analysis of the high-$z$ Swift gamma-ray burst GRB140515A ($z = 6.327$). The best estimate of the neutral hydrogen fraction of the intergalactic medium (IGM) towards the burst is $x_{HI} \leq 0.002$. The spectral absorption lines detected for this event are the weakest lines ever observed in gamma-ray burst afterglows, suggesting that GRB140515A exploded in a very low density environment. Its circum-burst medium is characterised by an average extinction (A$_{\rm V} \sim 0.1$) that seems to be typical of $z \ge 6$ events. The observed multi-band light curves are explained either with a very flat injected spectrum ($p = 1.7$) or with a multi-component emission ($p = 2.1$). In the second case a long-lasting central engine activity is needed in order to explain the late time X-ray emission. The possible origin of GRB140515A from a Pop III (or from a Pop II stars with local environment enriched by Pop III) massive star is unlikely.

Gamma Ray Bursts are detectable in the gamma-ray band if their jets are oriented towards the observer. However, for each GRB with a typical theta_jet, there should be ~2/theta_jet^2 bursts whose emission cone is oriented elsewhere in space. These off-axis bursts can be eventually detected when, due to the deceleration of their relativistic jets, the beaming angle becomes comparable to the viewing angle. Orphan Afterglows (OA) should outnumber the current population of bursts detected in the gamma-ray band even if they have not been conclusively observed so far at any frequency. We compute the expected flux of the population of orphan afterglows in the mm, optical and X-ray bands through a population synthesis code of GRBs and the standard afterglow emission model. We estimate the detection rate of OA by on-going and forthcoming surveys. The average duration of OA as transients above a given limiting flux is derived and described with analytical expressions: in general OA should appear as daily transients in optical surveys and as monthly/yearly transients in the mm/radio band. We find that ~ 2 OA yr^-1 could already be detected by Gaia and up to 20 OA yr^-1 could be observed by the ZTF survey. A larger number of 50 OA yr^-1 should be detected by LSST in the optical band. For the X-ray band, ~ 26 OA yr^-1 could be detected by the eROSITA. For the large population of OA detectable by LSST, the X-ray and optical follow up of the light curve (for the brightest cases) and/or the extensive follow up of their emission in the mm and radio band could be the key to disentangle their GRB nature from other extragalactic transients of comparable flux density.

Ten years of operations of the Swift satellite have allow us to collect a small sample of long Gamma-Ray Bursts (GRBs) at redshift larger than six. I will review here the present status of this research field and discuss the possible use of GRBs as a fundamental new tool to explore the early Universe, complementary to quasar and galaxy surveys.

The EXTraS project (Exploring the X-ray Transient and variable Sky) will harvest the hitherto unexplored temporal domain information buried in the serendipitous data collected by the European Photon Imaging Camera (EPIC) instrument onboard the ESA XMM-Newton X-ray observatory since its launch. This will include a search for fast transients, as well as a search and characterization of variability (both periodic and aperiodic) in hundreds of thousands of sources spanning more than nine orders of magnitude in time scale and six orders of magnitude in flux. X-ray results will be complemented by multiwavelength characterization of new discoveries. Phenomenological classification of variable sources will also be performed. All our results will be made available to the community. A didactic program in selected High Schools in Italy, Germany and the UK will also be implemented. The EXTraS project (2014-2016), funded within the EU/FP7 framework, is carried out by a collaboration including INAF (Italy), IUSS (Italy), CNR/IMATI (Italy), University of Leicester (UK), MPE (Germany) and ECAP (Germany).

We investigate signatures of population III (PopIII) stars in the metal-enriched environment of GRBs originating from population II-I (PopII/I) stars by using abundance ratios derived from numerical simulations that follow stellar evolution and chemical enrichment. We find that at $z>10$ more than $10%$ of PopII/I GRBs explode in a medium previously enriched by PopIII stars (we refer to them as GRBII$\rightarrow$III). Although the formation of GRBII$\rightarrow$III is more frequent than that of pristine PopIII GRBs (GRBIIIs), we find that the expected GRBII$\rightarrow$III observed rate is comparable to that of GRBIIIs, due to the usually larger luminosities of these latter. GRBII$\rightarrow$III events take place preferentially in small proto-galaxies with stellar masses $\rm M_\star \sim 10^{4.5} - 10^7\,\rm M_\odot$, star formation rates $\rm SFR \sim 10^{-3}-10^{-1}\,\rm M_\odot/yr$ and metallicities $Z \sim 10^{-4}-10^{-2}\,\rm Z_\odot$. On the other hand, galaxies with $Z < 10^{-2.8}\,\rm Z_\odot$ are dominated by metal enrichment from PopIII stars and should preferentially host GRBII$\rightarrow$III. Hence, measured GRB metal content below this limit could represent a strong evidence of enrichment by pristine stellar populations. We discuss how to discriminate PopIII metal enrichment on the basis of various abundance ratios observable in the spectra of GRBs' afterglows. By employing such analysis, we conclude that the currently known candidates at redshift $z\simeq 6$ -- i.e. GRB 050904 \cite[][]{2006Natur.440..184K} and GRB 130606A \cite[][]{2013arXiv1312.5631C} -- are likely not originated in environments pre-enriched by PopIII stars.

The structure of Gamma-Ray Burst (GRB) jets impacts on their prompt and afterglow emission properties. Insights into the still unknown structure of GRBs can be achieved by studying how different structures impact on the luminosity function (LF): i) we show that low ($10^{46} < L_{\rm iso} < 10^{48}$ erg/s) and high (i.e. with $L_{\rm iso} > 10^{50}$ erg/s) luminosity GRBs can be described by a unique LF; ii) we find that a uniform jet (seen on- and off-axis) as well as a very steep structured jet (i.e. $\epsilon(\theta) \propto \theta^{-s}$ with $s > 4$) can reproduce the current LF data; iii) taking into account the emission from the whole jet (i.e. including contributions from mildly relativistic, off-axis jet elements) we find that $E_{\rm iso}(\theta_{\rm v})$ (we dub this quantity "apparent structure") can be very different from the intrinsic structure $\epsilon(\theta)$: in particular, a jet with a Gaussian intrinsic structure has an apparent structure which is more similar to a power law. This opens a new viewpoint on the quasi-universal structured jet hypothesis.

The rapidly declining population of bright quasars at z~3 appears to make an increasingly small contribution to the ionising background at the HI Lyman limit. It is then generally though that massive stars in (pre-)galactic systems may provide the additional ionising flux needed to complete HI reionisation by z>6. A galaxy dominated background, however, may require that the escape fraction of Lyman continuum radiation from high redshift galaxies is as high as 10%, a value somewhat at odds with (admittedly scarce) observational constraints. High escape fractions from dwarf galaxies have been advocated, or, alternatively, a so-far undetected (or barely detected) population of unobscured, high-redshift faint AGNs. Here we question the latter hypothesis, and show that such sources, to be consistent with the measured level of the unresolved X-ray background at z=0, can provide a fraction of the HII filling factor not larger than 13% by z=6. The fraction rises to <27% in the somewhat extreme case of a constant comoving redshift evolution of the AGN emissivity. This still calls for a mean escape fraction of ionising photons from high-z galaxies >10%.

Gamma Ray Bursts (GRBs) are a powerful probe of the high redshift Universe. We present a tool to estimate the detection rate of high-z GRBs by a generic detector with defined energy band and sensitivity. We base this on a population model that reproduces the observed properties of GRBs detected by Swift, Fermi and CGRO in the hard X-ray and gamma-ray bands. We provide the expected cumulative distributions of the flux and fluence of simulated GRBs in different energy bands. We show that scintillator detectors, operating at relatively high energies (e.g. tens of keV to the MeV), can detect only the most luminous GRBs at high redshifts due to the link between the peak spectral energy and the luminosity (Ep-Liso) of GRBs. We show that the best strategy for catching the largest number of high-z bursts is to go softer (e.g. in the soft X-ray band) but with a very high sensitivity. For instance, an imaging soft X-ray detector operating in the 0.2-5 keV energy band reaching a sensitivity, corresponding to a fluence of ~10^-8 erg cm^-2, is expected to detect ~40 GRBs yr^-1 sr^-1 at z>5 (~3 GRBs yr^-1 sr^-1 at z>10). Once high-z GRBs are detected the principal issue is to secure their redshift. To this aim we estimate their NIR afterglow flux at relatively early times and evaluate the effectiveness of following them up and construct usable samples of events with any forthcoming GRB mission dedicated to explore the high-z Universe.

We compare the prompt intrinsic spectral properties of a sample of short Gamma--ray Burst (GRB) with the first 0.3 seconds (rest frame) of long GRBs observed by Fermi/GBM. We find that short GRBs and the first part of long GRBs lie on the same E_p--E_iso correlation, that is parallel to the relation for the time averaged spectra of long GRBs. Moreover, they are indistinguishable in the E_p--L_iso plane. This suggests that the emission mechanism is the same for short and for the beginning of long events, and both short and long GRBs are very similar phenomena, occurring on different timescales. If the central engine of a long GRB would stop after ~0.3 * (1+z) seconds the resulting event would be spectrally indistinguishable from a short GRB.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of gamma-ray bursts. For a summary, we refer to the paper.

The census of heavy elements (metals) produced by all stars through cosmic times up to present-day is limited to ~50%; of these only half are still found within their parent galaxy. The majority of metals is expelled from galaxies into the circumgalactic (or even more distant, intergalactic) space by powerful galactic winds, leaving unpleasant uncertainty on the amount, thermal properties and distribution of these key chemical species. These dispersed metals unavoidably absorb soft X-ray photons from distant sources. We show that their integrated contribution can be detected in the form of increasing X-ray absorption with distance, for all kinds of high-energy cosmic sources. Based on extensive cosmological simulations, we assess that $\sim$ 10\% of all cosmic metals reside in the intergalactic medium. Most of the X-ray absorption arises instead from a few discrete structures along the line of sight. These extended structures, possibly pin-pointing galaxy groups, contain million degree, metal-enriched gas, 100-1,000 times denser than the cosmic mean. An additional ~10% of cosmic metals could reside in this phase.

We investigated the rest frame spectral lags of two complete samples of bright long (50) and short (6) gamma-ray bursts (GRB) detected by Swift. We analysed the Swift/BAT data through a discrete cross-correlation function (CCF) fitted with an asymmetric Gaussian function to estimate the lag and the associated uncertainty. We find that half of the long GRBs have a positive lag and half a lag consistent with zero. All short GRBs have lags consistent with zero. The distributions of the spectral lags for short and long GRBs have different average values. Limited by the small number of short GRBs, we cannot exclude at more than 2 sigma significance level that the two distributions of lags are drawn from the same parent population. If we consider the entire sample of long GRBs, we do not find evidence for a lag-luminosity correlation, rather the lag-luminosity plane appears filled on the left hand side, thus suggesting that the lag-luminosity correlation could be a boundary. Short GRBs are consistent with the long ones in the lag-luminosity plane.

Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic outflows (jets) from accreting black holes at cosmological distances. Bright afterglows are produced when the outflow collides with the ambient medium. Afterglow polarization directly probes the magnetic properties of the jet, when measured minutes after the burst, and the geometric properties of the jet and the ambient medium when measured hours to days after the burst. High values of optical polarization detected minutes after burst in GRB 120308A indicate the presence of large-scale ordered magnetic fields originating from the central engine (the power source of the GRB). Theoretical models predict low degrees of linear polarization and negligable circular polarization at late times, when the energy in the original ejecta is quickly transferred to the ambient medium and propagates farther into the medium as a blastwave. Here we report the detection of circularly polarized optical light in the afterglow of GRB 121024A, measured 0.15 days after the burst. We show that the circular polarization is intrinsic to the afterglow and unlikely to be produced by dust scattering or plasma propagation effects. A possible explanation is to invoke anisotropic (rather than the commonly assumed isotropic) electron pitch angle distributions, and we suggest that new models are required to produce the complex microphysics of realistic shocks in relativistic jets.

The structure of Gamma Ray Burst (GRB) jets impacts on their prompt and afterglow emission properties. The jet of GRBs could be uniform, with constant energy per unit solid angle within the jet aperture, or it could instead be structured, namely with energy and velocity that depend on the angular distance from the axis of the jet. We try to get some insight about the still unknown structure of GRBs by studying their luminosity function. We show that low (1e46-1e48 erg/s) and high (i.e. with L > 1e50 erg/s) luminosity GRBs can be described by a unique luminosity function, which is also consistent with current lower limits in the intermediate luminosity range (1e48-1e50} erg/s). We derive analytical expressions for the luminosity function of GRBs in uniform and structured jet models and compare them with the data. Uniform jets can reproduce the entire luminosity function with reasonable values of the free parameters. A structured jet can also fit adequately the current data, provided that the energy within the jet is relatively strongly structured, i.e. E propto theta^{-k} with k > 4. The classical E propto theta^{-2} structured jet model is excluded by the current data.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We describe the GRB and All-sky Monitor Experiment (GAME) mission submitted by a large international collaboration (Italy, Germany, Czech Repubblic, Slovenia, Brazil) in response to the 2012 ESA call for a small mission opportunity for a launch in 2017 and presently under further investigation for subsequent opportunities. The general scientific objective is to perform measurements of key importance for GRB science and to provide the wide astrophysical community of an advanced X-ray all-sky monitoring system. The proposed payload was based on silicon drift detectors (~1-50 keV), CdZnTe (CZT) detectors (~15-200 keV) and crystal scintillators in phoswich (NaI/CsI) configuration (~20 keV-20 MeV), three well established technologies, for a total weight of ~250 kg and a required power of ~240 W. Such instrumentation allows a unique, unprecedented and very powerful combination of large field of view (3-4 sr), a broad energy energy band extending from ~1 keV up to ~20 MeV, an energy resolution as good as ~300 eV in the 1-30 keV energy range, a source location accuracy of ~1 arcmin. The mission profile included a launch (e.g., by Vega) into a low Earth orbit, a baseline sky scanning mode plus pointed observations of regions of particular interest, data transmission to ground via X-band (4.8 Gb/orbit, Alcantara and Malindi ground stations), and prompt transmission of GRB / transient triggers.

We present a carefully selected sample of short gamma-ray bursts (SGRBs) observed by the Swift satellite up to June 2013. Inspired by the criteria we used to build a similar sample of bright long GRBs (the BAT6 sample), we selected SGRBs with favorable observing conditions for the redshift determination on ground, ending up with a sample of 36 events, almost half of which with a redshift measure. The redshift completeness increases up to about 70% (with an average redshift value of z = 0.85) by restricting to those events that are bright in the 15-150 keV Swift Burst Alert Telescope energy band. Such flux-limited sample minimizes any redshift-related selection effects, and can provide a robust base for the study of the energetics, redshift distribution and environment of the Swift bright population of SGRBs. For all the events of the sample we derived the prompt and afterglow emission in both the observer and (when possible) rest frame and tested the consistency with the correlations valid for long GRBs. The redshift and intrinsic X-ray absorbing column density distributions we obtain are consistent with the scenario of SGRBs originated by the coalescence of compact objects in primordial binaries, with a possible minor contribution (~10%-25%) of binaries formed by dynamical capture (or experiencing large natal kicks). This sample is expected to significantly increase with further years of Swift activity.

We investigate the scaling relation between the observed amount of absorption in the X-ray spectra of Gamma Ray Burst (GRB) afterglows and the absorber redshift. Through dedicated numerical simulations of an ideal instrument, we establish that this dependence has a power law shape with index 2.4. However, for real instruments, this value depends on their low energy cut-off, spectral resolution and on the detector spectral response in general. We thus provide appropriate scaling laws for specific instruments. Finally, we discuss the possibility to measure the absorber redshift from X-ray data alone. We find that 10^5-10^6 counts in the 0.3-10 keV band are needed to constrain the redshift with 10% accuracy. As a test case we discuss the XMM-Newton observation of GRB 090618 at z=0.54. We are able to recover the correct redshift of this burst with the expected accuracy.

By combining cosmological simulations with Frontier Fields project lens models we find that, in the most optimistic case, galaxies as faint as $m \approx 33 - 34$ (AB magnitude at $1.6 \rm \mu m$) can be detected in the Frontier Fields. Such faint galaxies are hosted by dark matter halos of mass $\sim10^9 M_\odot$ and dominate the ionizing photon budget over currently observed bright galaxies, thus allowing for the first time the investigation of the dominant reionization sources. In addition, the observed number of these galaxies can be used to constrain the role of feedback in suppressing star formation in small halos: for example, if galaxy formation is suppressed in halos with circular velocity $v_c < 50$ km s$^{-1}$, galaxies fainter than $m=31$ should not be detected in the FFs.

Long-duration gamma-ray bursts (GRBs) have been found to be associated with broad-lined type-Ic supernovae (SNe), but only a handful of cases have been studied in detail. Prompted by the discovery of the exceptionally bright, nearby GRB130427A (redshift z=0.3399), we aim at characterising the properties of its associated SN2013cq. This is the first opportunity to test directly the progenitors of high-luminosity GRBs. We monitored the field of the Swift long duration GRB130427A using the 3.6-m TNG and the 8.2-m VLT during the time interval between 3.6 and 51.6 days after the burst. Photometric and spectroscopic observations revealed the presence of the type Ic SN2013cq. Spectroscopic analysis suggests that SN2013cq resembles two previous GRB-SNe, SN1998bw and SN2010bh associated with GRB980425 and XRF100316D, respectively. The bolometric light curve of SN2013cq, which is significantly affected by the host galaxy contribution, is systematically more luminous than that of SN2010bh ($\sim$ 2 mag at peak), but is consistent with SN1998bw. The comparison with the light curve model of another GRB-connected SN2003dh, indicates that SN2013cq is consistent with the model when brightened by 20%. This suggests a synthesised radioactive $^{56}$Ni mass of $\sim 0.4 M_\odot$. GRB130427A/SN2013cq is the first case of low-z GRB-SN connection where the GRB energetics are extreme ($E_{\rm \gamma, iso} \sim 10^{54}$ erg). We show that the maximum luminosities attained by SNe associated with GRBs span a very narrow range, but those associated with XRFs are significantly less luminous. On the other hand the isotropic energies of the accompanying GRBs span 6 orders of magnitude (10$^{48}$ erg $< E_{\rm \gamma, iso} <$ 10$^{54}$ erg), although this range is reduced when corrected for jet collimation. The GRB total radiated energy is in fact a small fraction of the SN energy budget.

We present the rest-frame light curves in the optical and X-ray bands of an unbiased and complete sample of Swift long Gamma-Ray Bursts (GRBs), namely the BAT6 sample. The unbiased BAT6 sample (consisting of 58 events) has the highest level of completeness in redshift ($\sim$ 95%), allowing us to compute the rest-frame X-ray and optical light curves for 55 and 47 objects, respectively. We compute the X-ray and optical luminosities accounting for any possible source of absorption (Galactic and intrinsic) that could affect the observed fluxes in these two bands. We compare the behaviour observed in the X-ray and in the optical bands to assess the relative contribution of the emission during the prompt and afterglow phases. We unarguably demonstrate that the GRBs rest-frame optical luminosity distribution is not bimodal, being rather clustered around the mean value Log(L$_{R}$) = 29.9 $\pm$ 0.8 when estimated at a rest frame time of 12 hr. This is in contrast with what found in previous works and confirms that the GRB population has an intrinsic unimodal luminosity distribution. For more than 70% of the events the rest-frame light curves in the X-ray and optical bands have a different evolution, indicating distinct emitting regions and/or mechanisms. The X-ray light curves normalised to the GRB isotropic energy (E$_{\rm iso}$), provide evidence for X-ray emission still powered by the prompt emission until late times ($\sim$ hours after the burst event). On the other hand, the same test performed for the E$_{\rm iso}$-normalised optical light curves shows that the optical emission is a better proxy of the afterglow emission from early to late times.

Orphan Afterglows (OA) are slow transients produced by Gamma Ray Bursts seen off-axis that become visible on timescales of days/years at optical/NIR and radio frequencies, when the prompt emission at high energies (X and gamma rays) has already ceased. Given the typically estimated jet opening angle of GRBs theta_jet ~ 3 deg, for each burst pointing to the Earth there should be a factor ~ 700 more GRBs pointing in other directions. Despite this, no secure OAs have been detected so far. Through a population synthesis code we study the emission properties of the population of OA at radio frequencies. OAs reach their emission peak on year-timescales and they last for a comparable amount of time. The typical peak fluxes (which depend on the observing frequency) are of few micro Jy in the radio band with only a few OA reaching the mJy level. These values are consistent with the upper limits on the radio flux of SN Ib/c observed at late times. We find that the OA radio number count distribution has a typical slope -1.7 at high fluxes and a flatter (-0.4) slope at low fluxes with a break at a frequency-dependent flux. Our predictions of the OA rates are consistent with the (upper) limits of recent radio surveys and archive searches for radio transients. Future radio surveys like VAST/ASKAP at 1.4 GHz should detect ~ 3x10^-3 OA deg^-2 yr-1, MeerKAT and EVLA at 8.4 GHz should see ~ 3x10^-1 OA deg-2 yr-1. The SKA, reaching the micro Jy flux limit, could see up to ~ 0.2-1.5 OA deg^-2 yr^-1. These rates also depend on the duration of the OA above a certain flux limit and we discuss this effect with respect to the survey cadence.

It has been proposed that the first, intermediate-mass ($\approx 10^{5-6}~M_\odot$) black holes might form through direct collapse of unpolluted gas in atomic-cooling halos exposed to a strong Lyman-Werner (LW) or near-infrared (NIR) radiation. As these systems are expected to be Compton-thick, photons above 13.6 eV are largely absorbed and re-processed into lower energy bands. It follows that direct collapse black holes (DCBHs) are very bright in the LW/NIR bands, typically outshining small high-redshift galaxies by more than 10 times. Once the first DCBHs form, they then trigger a runaway process of further DCBH formation, producing a sudden rise in their cosmic mass density. The universe enters the "DCBH era" at $z \approx 20$ when a large fraction of atomic-cooling halos are experiencing DCBH formation. By combining the clustering properties of the radiation sources with Monte Carlo simulations we show that in this scenario the DCBH mass density rises from $\sim 5$~$M_\odot$ Mpc$^{-3}$ at $z\sim 30$ to the peak value $\sim5\times10^5 M_\odot$ Mpc$^{-3}$ at $z \sim 14$ in our fiducial model. However, the abundance of \textit{active} (accreting) DCBHs drops after $z \sim 14$, as gas in the potential formation sites (unpolluted halos with virial temperature slightly above $10^4$~K) is photoevaporated. This effect almost completely suppresses DCBH formation after $z\sim 13$. The DCBH formation era lasts only $\approx 150$ Myr, but it might crucially provide the seeds of the supermassive black holes (SMBHs) powering $z\sim6$ quasars.

Long-duration Gamma-Ray Bursts (GRBs) are an extremely rare outcome of the collapse of massive stars, and are typically found in the distant Universe. Because of its intrinsic luminosity ($L\sim 3 \times 10^{53}$ erg s$^{-1}$) and its relative proximity ($z=0.34$), GRB 130427A was a unique event that reached the highest fluence observed in the gamma-ray band. Here we present a comprehensive multiwavelength view of GRB 130427A with Swift, the 2-m Liverpool and Faulkes telescopes and by other ground-based facilities, highlighting the evolution of the burst emission from the prompt to the afterglow phase. The properties of GRB 130427A are similar to those of the most luminous, high-redshift GRBs, suggesting that a common central engine is responsible for producing GRBs in both the contemporary and the early Universe and over the full range of GRB isotropic energies.

We aim at probing the emission mechanism of the accreting super massive black holes in the high redshift Universe. We study the X-ray spectrum of ULAS1120+0641, the highest redshift quasar detected so far at z=7.085, which has been deeply observed (340 ks) by XMM-Newton. Despite the long integration time the spectral analysis is limited by the poor statistics, with only 150 source counts being detected. We measured the spectrum in the 2-80 keV rest-frame (0.3-10 keV observed) energy band. Assuming a simple power law model we find a photon index of 2.0+/-0.3 and a luminosity of 6.7+/-0.3 10^44 erg/s in the 2-10 keV band, while the intrinsic absorbing column can be only loosely constrained (NH< 1E23 cm^-2). Combining our data with published data we calculate that the X-ray-to-optical spectral index alpha_OX is1.8+/-0.1, in agreement with the alpha_OX-UV luminosity correlation valid for lower redshift quasars. We expanded to high energies the coverage of the spectral energy distribution of ULAS1120+0641. This is the second time that a z >6 quasar has been investigated through a deep X-ray observation. In agreement with previous studies of z~6 AGN samples, we do not find any hint of evolution in the broadband energy distribution. Indeed from our dataset ULAS 1120+0641 is indistinguishable from the population of optically bright quasar at lower redshift.

One of the most elusive features of Gamma Ray Bursts (GRBs) is the sporadic emission prior to the main prompt event observed in at least $\sim 15\%$ of cases. These precursors have spectral and temporal properties similar to the main prompt emission, and smaller, but comparable, energetics. They are separated from the main event by a quiescent time that may be extremely long and, in some cases, more than one precursor has been observed in the same burst. Precursors are still a puzzle: despite many attempts none of the proposed models can account for all the observed features. Based on the complete sample of bright long GRBs observed by Swift (BAT6), we propose a new scenario for which precursors are explained by assuming that the central GRB engine is a newly born magnetar. In this model the precursor and the prompt emission arise from accretion of matter onto the surface of the magnetar. The accretion process can be halted by the centrifugal drag exerted by the rotating magnetosphere onto the in-falling matter, allowing for multiple precursors and very long quiescent times.

We estimate the initial bulk Lorentz factors Gamma_0 for GRBs that show the onset of the afterglow in their optical light curves. We find that Gamma_0 is strongly correlated with both the isotropic equivalent luminosity L_iso and energy E_iso and, with a larger scatter, also with the rest frame peak energy E_peak. These new correlations allow us to interpret the spectral energy correlations E_peak-L_iso (-E_iso) as a sequence of Gamma_0 factors. By accounting for the beaming effects, we find that the comoving frame properties of GRBs result clustered around typical values (e.g. L'_iso~5x10^48 erg/s). Moreover, it is theoretically predicted that there should be a link between the jet dynamics (Gamma_0) and its geometry (theta_jet). Through a population synthesis code we reconstruct the Gamma_0 and theta_jet distributions and search for a possible link between them. We find that Gamma_0 and theta_jet in GRBs should have log-normal distributions and they should be anti correlated (i.e. theta_jet^2*Gamma_0=const).

Starting from the Swift sample we define a complete sub-sample of 58 bright long Gamma Ray Bursts (GRB), 55 of them (95%) with a redshift determination, in order to characterize their properties. Our sample (BAT6) allows us to study the properties of the long GRB population and their evolution with cosmic time. We focus in particular on the GRB luminosity function, on the spectral-energy correlations of their prompt emission, on the nature of dark bursts, on possible correlations between the prompt and the X-ray afterglow properties, and on the dust extinction.

We discuss the results of the analysis of multi-wavelength data for the afterglows of GRB 081007 and GRB 090424, two bursts detected by Swift. One of them, GRB 081007, also shows a spectroscopically confirmed supernova, SN 2008hw, which resembles SN 1998bw in its absorption features, while the maximum luminosity is only about half as large as that of SN 1998bw. Bright optical flashes have been detected in both events, which allows us to derive solid constraints on the circumburst-matter density profile. This is particularly interesting in the case of GRB 081007, whose afterglow is found to be propagating into a constant-density medium, yielding yet another example of a GRB clearly associated with a massive star progenitor which did not sculpt the surroundings with its stellar wind. There is no supernova component detected in the afterglow of GRB 090424, likely due to the brightness of the host galaxy, comparable to the Milky Way. We show that the afterglow data are consistent with the presence of both forward- and reverse-shock emission powered by relativistic outflows expanding into the interstellar medium. The absence of optical peaks due to the forward shock strongly suggests that the reverse shock regions should be mildly magnetized. The initial Lorentz factor of outflow of GRB 081007 is estimated to be \Gamma ~ 200, while for GRB 090424 a lower limit of \Gamma > 170 is derived. We also discuss the prompt emission of GRB 081007, which consists of just a single pulse. We argue that neither the external forward-shock model nor the shock-breakout model can account for the prompt emission data and suggest that the single-pulse-like prompt emission may be due to magnetic energy dissipation of a Poynting-flux dominated outflow or to a dissipative photosphere.

Radio observations of Gamma Ray Bursts afterglows are fundamental in providing insights into their physics and environment, and in constraining the true energetics of these sources. Nonetheless, radio observations of GRB afterglows are presently sparse in the time/frequency domain. Starting from a complete sample of 58 bright Swift long bursts (BAT6), we constructed a homogeneous sub-sample of 38 radio detections/upper limits which preserves all the properties of the parent sample. One half of the bursts have detections between 1 and 5 days after the explosion with typical fluxes F>100 muJy at 8.4 GHz. Through a Population SYnthesis Code coupled with the standard afterglow Hydrodynamical Emission model (PSYCHE) we reproduce the radio flux distribution of the radio sub-sample. Based on these results we study the detectability in the time/frequency domain of the entire long GRB population by present and future radio facilities. We find that the GRBs that typically trigger Swift can be detected at 8.4 GHz by JVLA within few days with modest exposures even at high redshifts. The final SKA can potentially observe the whole GRB population provided that there will be a dedicated GRB gamma-ray detector more sensitive than Swift. For a sizable fraction (50%) of these bursts, SKA will allow us to perform radio-calorimetry, after the trans-relativistic transition (occurring ~100 d), providing an estimate of the true (collimation corrected) energetics of GRBs.

Long Gamma Ray Bursts (LGRBs) are related to the final stages of evolution of massive stars. As such, they should follow the star formation rate (SFR) of galaxies. We can use them to probe for star-forming galaxies (SFGs) in the distant universe following this assumption. The relation between the rate of LGRBs in a given galaxy and its SFR (that we call the LGRB "bias") may be complex, as we have indications that the LGRB hosts are not perfect analogues to the general population of SFGs. In this work, we try to quantify the dependence of the LGRB bias on physical parameters of their host galaxy such as the SFR or the stellar mass. We propose an empirical method based on the comparison of stellar mass functions (and SFR distributions) of LGRB hosts and of SFGs in order to find how the bias depends on the stellar mass or the SFR. By applying this method to a sample of LGRB hosts at redshifts lower than 1.1, where the properties of SFGs are well established, and where the properties of LGRB host galaxies can be deduced from observations (for stellar masses larger than 10**9.25 Msun and SFR larger than 1.8 Msun / yr), we find that the LGRB bias depends on both the stellar mass and SFR. We find that the bias decreases with the SFR, i.e. we see no preference for highly SFGs, once taken into account the larger number of massive stars in galaxies with larger SFR. We do not find any trend with the specific star formation rate (SSFR) but the dynamical range in SSFR in our study is narrow. Although through an indirect method, we relate these trends to a possible decrease of the LGRBs rate / SFR ratio with the metallicity. The method we propose suggests trends that may be useful to constrain models of LGRB progenitors, showing a clear decrease of the LGRB bias with the metallicity. This is promising for the future as the number of LGRB hosts studied will increase.

It has been recently suggested that supermassive black holes at z = 5-6 might form from super-fast (\dot M > 10^4 Msun/yr) accretion occurring in unstable, massive nuclear gas disks produced by mergers of Milky-Way size galaxies. Interestingly, such mechanism is claimed to work also for gas enriched to solar metallicity. These results are based on an idealized polytropic equation of state assumption, essentially preventing the gas from cooling. We show that under more realistic conditions, the disk rapidly (< 1 yr) cools, the accretion rate drops, and the central core can grow only to \approx 100 Msun. In addition, most of the disk becomes gravitationally unstable in about 100 yr, further quenching the accretion. We conclude that this scenario encounters a number of difficulties that possibly make it untenable.

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

A crucial challenge in astrophysics over the coming decades will be to understand the origins of supermassive black holes (SMBHs) that lie at the centres of most, if not all, galaxies. The processes responsible for the initial formation of these SMBHs and their early growth via accretion - when they are seen as Active Galactic Nuclei (AGN) - remain unknown. To address this challenge, we must identify low luminosity and obscured z>6 AGNs, which represent the bulk of early SMBH growth. Sensitive X-ray observations are a unique signpost of accretion activity, uncontaminated by star formation processes, which prevent reliable AGN identification at other wavelengths (e.g. optical, infrared). The Athena+ Wide Field Imager will enable X-ray surveys to be carried out two orders of magnitude faster than with Chandra or XMM-Newton, opening a new discovery space and identifying over 400 z>6 AGN, including obscured sources. Athena+ will also play a fundamental role to enhance the scientific return of future multiwavelength facilities that will probe the physical conditions within the host galaxies of early SMBHs, which is vital for understanding how SMBHs form, what fuels their subsequent growth, and to assess their impact on the early Universe. Follow-up of samples of z>6 galaxies with the Athena+ X-ray Integral Field Unit could also reveal the presence of highly obscured AGNs, thanks to the detection of strong iron lines. Thus, Athena+ will enable the first quantitative measurements of the extent and distribution of SMBH accretion in the early Universe.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Angular fluctuations of the Near InfraRed Background (NIRB) intensity are observed up to scales $\simlt 1^{\ensuremath{^{\circ}}}$. Their interpretation is challenging as even after removing the contribution from detected sources, the residual signal is $>10$ times higher than expected from distant galaxies below the detection limit and first stars. We propose here a novel interpretation in which early, intermediate mass, accreting direct collapse black holes (DCBH), which are too faint to be detected individually in current surveys, could explain the observed fluctuations. We find that a population of highly obscured ($N_{\rm H}\simgt 10^{25} \rm cm^{-2}$) DCBHs formed in metal-free halos with virial temperature $10^4$ K at $z\simgt 12$, can explain the observed level $\approx 10^{-3}$ (nW m$^{-2}$ sr$^{-1})^2$ of the 3.6 and 4.5 $\mu$m fluctuations on scales $>100''$. The signal on smaller scales is instead produced by undetected galaxies at low and intermediate redshifts. Albeit Compton-thick, at scales $\theta> 100''$ DCBHs produce a CXB (0.5-2 keV)-NIRB ($4.5 \rm \mu m$) cross-correlation signal of $\simeq 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ nW m$^{-2}$ sr$^{-1}$ slightly dependent on the specific value of the absorbing gas column ($N_{\rm H} \approx 10^{25} \rm cm^{-2}$) adopted and in agreement with the recent measurements by \cite{2012arXiv1210.5302C}. At smaller scales the cross-correlation is dominated by the emission of high-mass X-ray binaries (HMXB) hosted by the same low-$z$, undetected galaxies accounting for small scale NIRB fluctuations. These results outline the great potential of the NIRB as a tool to investigate the nature of the first galaxies and black holes.

In this paper we compute rest-frame extinctions for the afterglows of a sample of gamma-ray bursts complete in redshift. The selection criteria of the sample are based on observational high-energy parameters of the prompt emission and therefore our sample should not be biased against dusty sight-lines. It is therefore expected that our inferences hold for the general population of gamma-ray bursts. Our main result is that the optical/near-infrared extinction of gamma-ray burst afterglows in our sample does not follow a single distribution. 87% of the events are absorbed by less than 2 mag, and 50% suffer from less than 0.3-0.4 mag extinction. The remaining 13% of the afterglows are highly absorbed. The true percentage of gamma-ray burst afterglows showing high absorption could be even higher since a fair fraction of the events without reliable redshift measurement are probably part of this class. These events may be due to highly dusty molecular clouds/star forming regions associated with the gamma-ray burst progenitor or along the afterglow line of sight, and/or to massive dusty host galaxies. No clear evolution in the dust extinction properties is evident within the redshift range of our sample, although the largest extinctions are at z~1.5-2, close to the expected peak of the star formation rate. Those events classified as dark are characterized, on average, by a higher extinction than typical events in the sample. A correlation between optical/near-infrared extinction and hydrogen-equivalent column density based on X-ray studies is shown although the observed NH appears to be well in excess compared to those observed in the Local Group. Dust extinction does not seem to correlate with GRB energetics or luminosity.

LOFT is a satellite mission currently in Assessment Phase for the ESA M3 selection. The payload is composed of the Large Area Detector (LAD), with 2-50 keV energy band, a peak effective area of about 10 m2 and an energy resolution better than 260 eV, and the Wide Field Monitor (WFM), a coded mask imager with a FOV of several steradians, an energy resolution of about 300 eV and a point source location accuracy of 1 arcmin in the 2-30 keV energy range. Based on preliminary considerations and estimates, we show how the scientific performances of the WFM are particularly suited to investigate some of the most relevant open issues in the study of GRBs: the physics of the prompt emission, the spectral absorption features by circum-burst material (and hence the nature of the progenitors), the population and properties of XRFs, and the detection and rate of high-z GRBs. Measurements of the early afterglow emission with the Large Area Detector (LAD) may also be possible depending on the mission slewing capabilities and TOO observations policy. More refined analysis and simulations are in progress.

Several independent measurements have confirmed the existence of fluctuations ($\delta F_{\rm obs}\approx 0.1 \rm nW/m^{2}/sr$ at $3.6 \rm \mu m$) up to degree angular scales in the source-subtracted Near InfraRed Background (NIRB) whose origin is unknown. By combining high resolution cosmological N-body/hydrodynamical simulations with an analytical model, and by matching galaxy Luminosity Functions (LFs) and the constraints on reionization simultaneously, we predict the NIRB absolute flux and fluctuation amplitude produced by high-$z$ ($z > 5$) galaxies (some of which harboring Pop III stars, shown to provide a negligible contribution). This strategy also allows us to make an empirical determination of the evolution of ionizing photon escape fraction: we find $f_{\rm esc} = 1$ at $z \ge 11$, decreasing to $\approx 0.05$ at $z = 5$. In the wavelength range $1.0-4.5 \rm \mu m$, the predicted cumulative flux is $F =0.2-0.04 \rm nW/m^2/sr$. However, we find that the radiation from high-$z$ galaxies (including those undetected by current surveys) is insufficient to explain the amplitude of the observed fluctuations: at $l=2000$, the fluctuation level due to $z > 5$ galaxies is $\delta F = 0.01-0.002 \rm nW/m^2/sr$, with a relative wavelength-independent amplitude $\delta F/F = 4$%. The source of the missing power remains unknown. This might indicate that an unknown component/foreground, with a clustering signal very similar to that of high-$z$ galaxies, dominates the source-subtracted NIRB fluctuation signal.

We investigate the nature of high-z host galaxies of long Gamma-Ray Bursts (LGRBs) by means of state-of-the-art numerical simulations of cosmic structure formation and evolution of galaxies. We combine results from different runs with various box sizes and resolutions. By assigning to each simulated galaxy the probability to host a LGRB, assumed to be proportional to the mass of young stars, we provide a full description of the physical properties of high-z LGRB host galaxy population. We find that LGRBs at z>6 are hosted in galaxies with typical star formation rates SFR \sim 0.03-0.3 Msun yr^{-1}, stellar masses M \sim 10^{6-8} Msun, and metallicities Z \sim 0.01-0.1 Zsun. Furthermore, the ratio between their doubling time and the corresponding cosmic time seems to be universally equal to ~0.1-0.3, independently from the redshift. The distribution of their UV luminosity places LGRB hosts in the faint-end of the galaxy luminosity function, well below the current capabilities of space- or ground-based optical facilities. This is in line with recent reports of non-detection of LGRB hosts using extremely deep HST and VLT observations. In conclusion, high-z LGRBs are found to trace the position of those faint galaxies that are thought to be the major actors in the re-ionization of the Universe.

We place firm upper limits on the global accretion history of massive black holes at z>5 from the recently measured unresolved fraction of the cosmic X-ray background. The maximum allowed unresolved intensity observed at 1.5 keV implies a maximum accreted-mass density onto massive black holes of rho_acc < 1.4E4 M_sun Mpc^{-3} for z>5. Considering the contribution of lower-z AGNs, the value reduces to rho_acc < 0.66E4 M_sun Mpc^{-3}. The tension between the need for the efficient and rapid accretion required by the observation of massive black holes already in place at z>7 and the strict upper limit on the accreted mass derived from the X-ray background may indicate that black holes are rare in high redshift galaxies, or that accretion is efficient only for black holes hosted by rare galaxies.

We recently found that Gamma Ray Burst energies and luminosities, in their comoving frame, are remarkably similar. This, coupled with the clustering of energetics once corrected for the collimation factor, suggests the possibility that all bursts, in their comoving frame, have the same peak energy E'peak (of the order of a few keV) and the same energetics of the prompt emission E'gamma (of the order of 2e48 erg). The large diversity of bursts energies is then due to the different bulk Lorentz factor Gamma and jet aperture angle theta_jet. We investigated, through a population synthesis code, what are the distributions of Gamma and theta_jet compatible with the observations. Both quantities must have preferred values, with log-normal best fitting distributions and <Gamma0> ~ 275 and <\theta_jet> ~ 8.7 degree. Moreover, the peak values of the Gamma and theta_jet distributions must be related - theta_jet^2.5 Gamma =const: the narrower the jet angle, the larger the bulk Lorentz factor. We predict that ~6% of the bursts that point to us should not show any jet break in their afterglow light curve since they have sin(theta_jet)<1/Gamma. Finally, we estimate that the local rate of GRBs is ~0.3% of all local SNIb/c and ~2.5% of local hypernovae, i.e. SNIb/c with broad absorption lines.

We study the spectral properties of the unresolved cosmic X-ray background (CXRB) in the 1.5-7.0 keV energy band with the aim of providing an observational constraint on the statistical properties of those sources that are too faint to be individually probed. We made use of the Swift X-ray observation of the Chandra Deep Field South complemented by the Chandra data. Exploiting the lowest instrument background (Swift) together with the deepest observation ever performed (Chandra) we measured the unresolved emission at the deepest level and with the best accuracy available today. We find that the unresolved CXRB emission can be modeled by a single power law with a very hard photon index Gamma=0.1+/-0.7 and a flux of 5(+/-3)E-12 cgs in the 2.0-10 keV energy band (1 sigma error). Thanks to the low instrument background of the Swift-XRT, we significantly improved the accuracy with respect to previous measurements. These results point towards a novel ingredient in AGN population synthesis models, namely a positive evolution of the Compton-thick AGN population from local Universe to high redshift.

The association of Type Ic SNe with long-duration GRBs is well established. We endeavor, through accurate ground-based observational campaigns, to characterize these SNe at increasingly high redshifts. We obtained a series of optical photometric and spectroscopic observations of the Type Ic SN2012bz associated with the Swift long-duration GRB120422A (z=0.283) using the 3.6-m TNG and the 8.2-m VLT telescopes. The peak times of the light curves of SN2012bz in various optical filters differ, with the B-band and i'-band light curves reaching maximum at ~9 and ~23 rest-frame days, respectively. The bolometric light curve has been derived from individual bands photometric measurements, but no correction for the unknown contribution in the near-infrared (probably around 10-15%) has been applied. Therefore, the present light curve should be considered as a lower limit to the actual UV-optical-IR bolometric light curve. This pseudo-bolometric curve reaches its maximum (Mbol = -18.56 +/- 0.06) at 13 +/- 1 rest-frame days; it is similar in shape and luminosity to the bolometric light curves of the SNe associated with z<0.2 GRBs and more luminous than those of SNe associated with XRFs. A comparison with the model generated for the bolometric light curve of SN2003dh suggests that SN2012bz produced only about 15% less 56Ni than SN2003dh, about 0.35 Msol. Similarly the VLT spectra of SN2012bz, after correction for Galactic extinction and for the contribution of the host galaxy, suggest comparable explosion parameters with those observed in SN2003dh (EK~3.5 x 10^52 erg, Mej~7 Msol) and a similar progenitor mass (~25-40 Msol). GRB120422A is consistent with the Epeak-Eiso and the EX,iso-Egamma,iso-E_peak relations. GRB120422A/SN2012bz shows the GRB-SN connection at the highest redshift so far accurately monitored both photometrically and spectroscopically.

The jet opening angle theta_jet and the bulk Lorentz factor Gamma_0 are crucial parameters for the computation of the energetics of Gamma Ray Bursts (GRBs). From the ~30 GRBs with measured theta_jet or Gamma_0 it is known that: (i) the real energetic E_gamma, obtained by correcting the isotropic equivalent energy E_iso for the collimation factor ~theta_jet^2, is clustered around 10^50-10^51 erg and it is correlated with the peak energy E_p of the prompt emission and (ii) the comoving frame E'_p and E'_gamma are clustered around typical values. Current estimates of Gamma_0 and theta_jet are based on incomplete data samples and their observed distributions could be subject to biases. Through a population synthesis code we investigate whether different assumed intrinsic distributions of Gamma_0 and theta_jet can reproduce a set of observational constraints. Assuming that all bursts have the same E'_p and E'_gamma in the comoving frame, we find that Gamma_0 and theta_jet cannot be distributed as single power-laws. The best agreement between our simulation and the available data is obtained assuming (a) log-normal distributions for theta_jet and Gamma_0 and (b) an intrinsic relation between the peak values of their distributions, i.e theta_jet^2.5*Gamma_0=const. On average, larger values of Gamma_0 (i.e. the "faster" bursts) correspond to smaller values of theta_jet (i.e. the "narrower"). We predict that ~6% of the bursts that point to us should not show any jet break in their afterglow light curve since they have sin(theta_jet)<1/Gamma_0. Finally, we estimate that the local rate of GRBs is ~0.3% of all local SNIb/c and ~4.3% of local hypernovae, i.e. SNIb/c with broad-lines.

We propose to use high-redshift long $\gamma$-ray bursts (GRBs) as cosmological tools to constrain the amount of primordial non-Gaussianity in the density field. By using numerical, N-body, hydrodynamic, chemistry simulations of different cosmological volumes with various Gaussian and non-Gaussian models, we self-consistently relate the cosmic star formation rate density to the corresponding GRB rate. Assuming that GRBs are fair tracers of cosmic star formation, we find that positive local non-Gaussianities, described in terms of the non-linear parameter, \fnl, might boost significantly the GRB rate at high redshift, $z \gg 6$. Deviations with respect to the Gaussian case account for a few orders of magnitude if \fnl$\sim 1000$, one order of magnitude for \fnl$\sim 100$, and a factor of $\sim 2$ for \fnl$\sim 50$. These differences are found only at large redshift, while at later times the rates tend to converge. Furthermore, a comparison between our predictions and the observed GRB data at $z > 6$ allows to exclude large negative \fnl, consistently with previous works. Future detections of any long GRB at extremely high redshift ($z\sim 15-20$) could favor non-Gaussian scenarios with positive \fnl. More stringent constraints require much larger high-$z$ GRB complete samples, currently not available in literature. By distinguishing the contributions to the GRB rate from the metal-poor population III regime, and the metal-enriched population II-I regime, we conclude that the latter is a more solid tracer of the underlying matter distribution, while the former is strongly dominated by feedback mechanisms from the first, massive, short-lived stars, rather than by possible non-Gaussian fluctuations. This holds quite independently of the assumed population III initial mass function.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

We investigate wheter there is any correlation between the X-ray afterglow luminosity and the prompt emission properties of a carefully selected sub-sample of bright Swift long Gamma-Ray Bursts (GRBs) nearly complete in redshift (~90%). Being free of selection effects (except flux limit), this sample provides the possibility to compare the rest frame physical properties of GRB prompt and afterglow emission in an unbiased way. The afterglow X-ray luminosities are computed at four different rest frame times (5 min, 1 hr, 11 hr and 24 hr after trigger) and compared with the prompt emission isotropic energy E_iso, the isotropic peak luminosity L_iso and the rest frame peak energy E_peak. We find that the rest frame afterglow X-ray luminosity do correlate with these prompt emission quantities, but the significance of each correlation decreases over time. This result is in agreement with the idea that the GRB X-ray light curve can be described as the result of a combination of different components whose relative contribution and weight change with time, with the prompt and afterglow emission dominating at early and late time, respectively. In particular, we found evidence that the plateau and the shallow decay phase often observed in GRB X-ray light curves are powered by activity from the central engine. The existence of the L_X-E_iso correlation at late times (t_rf > 11 hr) suggests a similar radiative efficiency among different bursts with on average about 6% of the total kinetic energy powering the prompt emission.

In this paper we present a sample of 10 short gamma-ray bursts (GRBs) with a robust redshift determination, discovered by the Swift satellite up to January 2011. We measure their X-ray absorbing column densities and collect data on the host galaxy offsets. We find evidence for intrinsic absorption and no correlation between the intrinsic absorbing column density and the projected offset of the GRB from its host galaxy center. We find that the properties in the gamma regime (T90, fluence and 1-s peak photon flux) of short GRBs with "bright" and "faint" X-ray afterglow likely disfavour different prompt emission mechanisms. Host galaxy offset and GRB duration (T90) do not correlate. Instead, there is a hint of anti-correlation between the effective radius normalised host galaxy offset and T90. Finally, we examine the properties of short GRBs with short-lived and long-lived X-ray afterglows, finding that some short GRBs with short-lived X-ray afterglows have their optical afterglow detected. In light of this, the X-ray afterglow duration does not seem to be an unique indicator of a specific progenitor and/or environment for short GRBs.

We present the spectroscopic and photometric evolution of the nearby (z = 0.059) spectroscopically confirmed type Ic supernova, SN 2010bh, associated with the soft, long-duration gamma-ray burst (X-ray flash) GRB 100316D. Intensive follow-up observations of SN 2010bh were performed at the ESO Very Large Telescope (VLT) using the X-shooter and FORS2 instruments. Owing to the detailed temporal coverage and the extended wavelength range (3000--24800 A), we obtained an unprecedentedly rich spectral sequence among the hypernovae, making SN 2010bh one of the best studied representatives of this SN class. We find that SN 2010bh has a more rapid rise to maximum brightness (8.0 +/- 1.0 rest-frame days) and a fainter absolute peak luminosity (L_bol~3e42 erg/s) than previously observed SN events associated with GRBs. Our estimate of the ejected (56)Ni mass is 0.12 +/- 0.02 Msun. From the broad spectral features we measure expansion velocities up to 47,000 km/s, higher than those of SNe 1998bw (GRB 980425) and 2006aj (GRB 060218). Helium absorption lines He I lambda5876 and He I 1.083 microm, blueshifted by ~20,000--30,000 km/s and ~28,000--38,000 km/s, respectively, may be present in the optical spectra. However, the lack of coverage of the He I 2.058 microm line prevents us from confirming such identifications. The nebular spectrum, taken at ~186 days after the explosion, shows a broad but faint [O I] emission at 6340 A. The light-curve shape and photospheric expansion velocities of SN 2010bh suggest that we witnessed a highly energetic explosion with a small ejected mass (E_k ~ 1e52 erg and M_ej ~ 3 Msun). The observed properties of SN 2010bh further extend the heterogeneity of the class of GRB supernovae.

Supernovae Ib/c are likely to be associated to long GRBs, therefore it is important to compare the SN rate in galaxies with the GRB rate. To do that we computed Type Ib/c SN rates in galaxies of different morphological type by assuming different histories of star formation and different supernova Ib/c progenitors. We included some recent suggestions about the dependence of the minimum mass of single Wolf-Rayet (WR) stars upon the stellar metallicity and therefore upon galactic chemical evolution. We adopted several cosmic star formation rates as functions of cosmic time, either observationally or theoretically derived, including the one computed with our galaxy models. Then we computed the cosmic Type Ib/c SN rates. We derived the following conclusions: i) the ratio cosmic GRB - Type Ib/c rate varies in the range 10^{-2}-10^{-4} in the whole redshift range, thus suggesting that only a small fraction of all the Type Ib/c SNe gives rise to GRBs. ii) The metallicity dependence of Type Ib/c SN progenitors produces lower cosmic SN Ib/c rates at early times, for any chosen cosmic star formation rate. iii) Different theoretical cosmic star formation rates, computed under different scenarios of galaxy formation, produce SN Ib/c cosmic rates which differ mainly at very high redshift. However, it is difficult to draw firm conclusions on the high redshift trend because of the large uncertainties in the data. iv) GRBs can be important tracers of star formation at high redshift if their luminosity function does not vary with redshift and they can help in discriminating among different galaxy formation models.

The relative importance of metals and dust grains in the formation of the first low-mass stars has been a subject of debate. The recently discovered Galactic halo star SDSS J102915+172927 (Caffau et al. 2011) has a mass less than 0.8 Msun and a metallicity of Z = 4.5 10^{-5} Zsun. We investigate the origin and properties of this star by reconstructing the physical conditions in its birth cloud. We show that the observed elemental abundance trend of SDSS J102915+172927 can be well fitted by the yields of core-collapse supernovae with metal-free progenitors of 20 Msun and 35 Msun. Using these selected supernova explosion models, we compute the corresponding dust yields and the resulting dust depletion factor taking into account the partial destruction by the supernova reverse shock. We then follow the collapse and fragmentation of a star forming cloud enriched by the products of these SN explosions at the observed metallicity of SDSS J102915+172927. We find that [0.05 - 0.1] Msun mass fragments, which then lead to the formation of low-mass stars, can occur provided that the mass fraction of dust grains in the birth cloud exceeds 0.01 of the total mass of metals and dust. This, in turn, requires that at least 0.4 Msun of dust condense in the first supernovae, allowing for moderate destruction by the reverse shock. If dust formation in the first supernovae is less efficient or strong dust destruction does occur, the thermal evolution of the SDSS J102915+172927 birth cloud is dominated by molecular cooling, and only > 8 Msun fragments can form. We conclude that the observed properties of SDSS J102915+172927 support the suggestion that dust must have condensed in the ejecta of the first supernovae and played a fundamental role in the formation of the first low-mass stars.

We study the possible effects of selection biases on the Ep-Liso correlation caused by the unavoidable presence of flux-limits in the existing samples of Gamma Ray Bursts (GRBs). We consider a well defined complete sample of Swift GRBs and perform Monte Carlo simulations of the GRB population under different assumptions for their luminosity functions. If we assume that there is no correlation between the peak energy Ep and the isotropic luminosity Liso, we are unable to reproduce it as due to the flux limit threshold of the Swift complete sample. We can reject the null hypothesis that there is no intrinsic correlation between Ep and Liso at more than 2.7 sigma level of confidence. This result is robust against the assumptions of our simulations and it is confirmed if we consider, instead of Swift, the trigger threshold of the Batse instrument. Therefore, there must be a physical relation between these two quantities. Our simulations seem to exclude, at a lower confidence level of 1.6 sigma, the possibility that the observed Ep-Liso correlation among different bursts is caused by a boundary, i.e. such that for any given Ep, we see only the largest Liso, which has a flux above the threshold of the current instruments.

A complete sample of bright Swift Gamma-ray Bursts (GRBs) has been recently selected by Salvaterra et al. (2011). The sample has a high level of completeness in redshift (91%). We derive here the intrinsic absorbing X-ray column densities of these GRBs making use of the Swift X-ray Telescope data. This distribution has a mean value of log(NH/cm-2)=21.7+-0.5. This value is consistent with the distribution of the column densities derived from the total sample of GRBs with redshift. We find a mild increase of the intrinsic column density with redshift. This can be interpreted as due to the contribution of intervening systems along the line of sight. Making use of the spectral index connecting optical and X-ray fluxes at 11 hr (beta_OX), we investigate the relation of the intrinsic column density and the GRB `darkness'. We find that there is a very tight correlation between dark GRBs and high X-ray column densities. This clearly indicates that the dark GRBs are formed in a metal-rich environment where dust must be present.

We use a nearly complete sample of Gamma Ray Bursts (GRBs) detected by the Swift satellite to study the correlations between the spectral peak energy Ep of the prompt emission, the isotropic energetics Eiso and the isotropic luminosity Liso. This GRB sample is characterized by a high level of completeness in redshift (90%). This allows us to probe in an unbiased way the issue related to the physical origin of these correlations against selection effects. We find that one burst, GRB 061021, is an outlier to the Ep-Eiso correlation. Despite this case, we find strong Ep-Eiso and Ep-Liso correlations for the bursts of the complete sample. Their slopes, normalisations and dispersions are consistent with those found with the whole sample of bursts with measured redshift and Ep. This means that the biases present in the total sample commonly used to study these correlations do not affect their properties. Finally, we also find no evolution with redshift of the Ep-Eiso and Ep-Liso correlations.

We study a model in which supermassive black holes (SMBHs) can grow by the combined action of gas accretion on heavy seeds and mergers of both heavy (m_s^h=10^5 Msol) and light (m_s^l = 10^2 Msol) seeds. The former result from the direct collapse of gas in T_s^h >1.5x10^4K, H_2-free halos; the latter are the endproduct of a standard H_2-based star formation process. The H_2-free condition is attained by exposing halos to a strong (J_21 > 10^3) Lyman-Werner UV background produced by both accreting BHs and stars, thus establishing a self-regulated growth regime. We find that this condition is met already at z close to 18 in the highly biased regions in which quasars are born. The key parameter allowing the formation of SMBHs by z=6-7 is the fraction of halos that can form heavy seeds: the minimum requirement is that f_heavy>0.001; SMBH as large as 2x10^10 Msol can be obtained when f_heavy approaches unity. Independently of f_heavy, the model produces a high-z stellar bulge-black hole mass relation which is steeper than the local one, implying that SMBHs formed before their bulge was in place. The formation of heavy seeds, allowed by the Lyman-Werner radiative feedback in the quasar-forming environment, is crucial to achieve a fast growth of the SMBH by merger events in the early phases of its evolution, i.e. z>7. The UV photon production is largely dominated by stars in galaxies, i.e. black hole accretion radiation is sub-dominant. Interestingly, we find that the final mass of light BHs and of the SMBH in the quasar is roughly equal by z=6; by the same time only 19% of the initial baryon content has been converted into stars. The SMBH growth is dominated at all epochs z > 7.2 by mergers (exceeding accretion by a factor 2-50); at later times accretion becomes by far the most important growth channel. We finally discuss possible shortcomings of the model.

We present a carefully selected sub-sample of Swift Long Gamma-ray Bursts (GRBs), that is complete in redshift. The sample is constructed by considering only bursts with favorable observing conditions for ground-based follow-up searches, that are bright in the 15-150 keV Swift/BAT band, i.e. with 1-s peak photon fluxes in excess to 2.6 ph s^-1 cm^-2. The sample is composed by 58 bursts, 52 of them with redshift for a completeness level of 90%, while another two have a redshift constraint, reaching a completeness level of 95%. For only three bursts we have no constraint on the redshift. The high level of redshift completeness allows us for the first time to constrain the GRB luminosity function and its evolution with cosmic times in a unbiased way. We find that strong evolution in luminosity (d_l=2.3\pm 0.6) or in density (d_d=1.7\pm 0.5) is required in order to account for the observations. The derived redshift distribution in the two scenarios are consistent with each other, in spite of their different intrinsic redshift distribution. This calls for other indicators to distinguish among different evolution models. Complete samples are at the base of any population studies. In future works we will use this unique sample of Swift bright GRBs to study the properties of the population of long GRBs.

We study the properties of the population of optically dark events present in a carefully selected complete sample of bright Swift long gamma-ray bursts. The high level of completeness in redshift of our sample (52 objects out of 58) allow us to establish the existence of a genuine dark population and we are able to estimate the maximum fraction of dark burst events (~30%) expected for the whole class of long gamma-ray burst. The redshift distribution of this population of dark bursts is similar to the one of the whole sample. Interestingly, the rest-frame X-ray luminosity (and the de-absorbed X-ray flux) of the sub-class of dark bursts is slightly higher than the average luminosity of the non-dark events. At the same time the prompt properties do not differ and the optical flux of dark events is at the lower tail of the optical flux distribution, corrected for Galactic absorption. All these properties suggest that dark bursts events generate in much denser environments with respect to normal bright events. We can therefore exclude the high-z and the low-density scenarios and conclude that the major cause of the origin of optically dark events is the dust extinction.

GRBs generate an afterglow emission that can be detected from radio to X-rays during days, or even weeks after the initial explosion. The peak of this emission crosses the mm/submm range during the first hours to days, making their study in this range crucial for constraining the models. Observations have been limited until now due to the low sensitivity of the observatories in this range. We present observations of 10 GRB afterglows obtained from APEX and SMA, as well as the first detection of a GRB with ALMA, and put them into context with all the observations that have been published until now in the spectral range that will be covered by ALMA. The catalogue of mm/submm observations collected here is the largest to date and is composed of 102 GRBs, of which 88 had afterglow observations, whereas the rest are host galaxy searches. With our programmes, we contributed with data of 11 GRBs and the discovery of 2 submm counterparts. In total, the full sample, including data from the literature, has 22 afterglow detections with redshift ranging from 0.168 to 8.2. GRBs have been detected in mm/submm wavelengths with peak luminosities spanning 2.5 orders of magnitude, the most luminous reaching 10^33erg s^-1 Hz^-1. We observe a correlation between the X-ray brightness at 0.5 days and the mm/submm peak brightness. Finally we give a rough estimate of the distribution of peak flux densities of GRB afterglows, based on the current mm/submm sample. Observations in the mm/submm bands have been shown to be crucial for our understanding of the physics of GRBs, but have until now been limited by the sensitivity of the observatories. With the start of the operations at ALMA, the sensitivity will be increased by more than an order of magnitude. Our estimates predict that, once completed, ALMA will detect up to 98% of the afterglows if observed during the passage of the peak synchrotron emission.

The tidal disruption of a solar mass star around a supermassive black hole has been extensively studied analytically and numerically. In these events the star develops into an elongated banana-shaped structure. After completing an eccentric orbit, the bound debris fall onto the black hole, forming an accretion disk and emitting radiation. The same process may occur on planetary scales, if a minor body orbits too close to its star. In the Solar System, comets fall directly onto our Sun or onto planets. If the star is a compact object, the minor body can become tidally disrupted. Indeed, one of the first mechanisms invoked to produce strong gamma-ray emission involved accretion of comets onto neutron stars in our Galaxy. Here we report that the peculiarities of the `Christmas' gamma-ray burst (GRB 101225A) can be explained by a tidal disruption event of a minor body around a Galactic isolated neutron star. This result would indicate either that minor bodies can be captured by compact stellar remnants more frequently than it occurs in our Solar System or that minor body formation is relatively easy around millisecond radio pulsars. A peculiar Supernova associated to a GRB may provide an alternative explanation.

We compute the contribution of blazars (Flat Spectrum Radio Quasars and BL-Lacs) to the Fermi-LAT Extragalactic gamma-ray Background, assuming a gamma-ray luminosity function proportional to the luminosity function of radio galaxies and the Spectral Energy Distribution of blazars given by the so-called blazar sequence. We find that blazars account for the 45% of the total Fermi-LAT EGB. However, at low (<10 GeV) and high (>50 GeV) energies the blazar contribution falls short to account for the observed background light. At E<10 GeV, the Extragalactic gamma-ray Background can be fully explained by considering the additional contribution of star-forming galaxies, whereas at very high energies, where the absorption due to Extragalactic Background Light dominates, the observed background calls for an additional unknown component.

We perform a detailed study of the gamma-ray burst GRB091127/SN2009nz host galaxy at z=0.490 using the VLT/X-shooter spectrograph in slit and integral-field unit (IFU). From the analysis of the optical and X-ray afterglow data obtained from ground-based telescopes and Swift-XRT we confirm the presence of a bump associated with SN2009nz and find evidence of a possible jet break in the afterglow lightcurve. The X-shooter afterglow spectra reveal several emission lines from the underlying host, from which we derive its integrated properties. These are in agreement with those of previously studied GRB-SN hosts and, more generally, with those of the long GRB host population. We use the Hubble Space Telescope and ground based images of the host to determine its stellar mass (M_star). Our results extend to lower M_star values the M-Z plot derived for the sample of long GRB hosts at 0.3<z<1.0 adding new information to probe the faint end of the M-Z relation and the shift of the LGRB host M-Z relation from that found from emission line galaxy surveys. Thanks to the IFU spectroscopy we can build the 2D velocity, velocity dispersion and star formation rate (SFR) maps. They show that the host galaxy has a perturbed rotation kinematics with evidence of a SFR enhancement consistent with the afterglow position.

Because massive, low-metallicity population III (PopIII) stars may produce very powerful long gamma-ray bursts (LGRBs), high-redshift GRB observations could probe the properties of the first stars. We analyze the correlation between early PopIII stars and LGRBs by using cosmological N-body/hydrodynamical simulations, which include detailed chemical evolution, cooling, star formation, feedback effects and the transition between PopIII and more standard population I/II (PopII/I) stars. From the Swift observed rate of LGRBs, we estimate the fraction of black holes that will produce a GRB from PopII/I stars to be in the range 0.028<f_{GRB}<0.140, depending on the assumed upper metallicity of the progenitor. Assuming that as of today no GRB event has been associated to a PopIII star, we estimate the upper limit for the fraction of LGRBs produced by PopIII stars to be in the range 0.006<f_{GRB}<0.022. When we apply a detection threshold compatible with the BAT instrument, we find that the expected fraction of PopIII GRBs (GRB3) is ~10% of the full LGRB population at z>6, becoming as high has 40% at z>10. Finally, we study the properties of the galaxies hosting our sample of GRB3. We find that the average metallicity of the galaxies hosting a GRB3 is typically higher than the critical metallicity used to select the PopIII stars, due to the efficiency in polluting the gas above such low values. We also find that the highest probability of finding a GRB3 is within galaxies with a stellar mass <10^7 Msun, independently from the redshift.

We present a two-parameter model of the extragalactic gamma-ray background (EGB) in the 0.1-100 GeV range as measured by the Large Area Telescope (LAT) onboard the Fermi satellite. The EGB can be fully explained as the sum of three distinct components, namely blazars, non-beamed AGNs (Seyfert galaxies and QSOs), and cosmic rays from star-forming galaxies. The contribution to the background from beamed sources is obtained by fitting the Fermi-LAT blazar differential number counts assuming that the gamma-ray luminosity function is directly proportional to the radio luminosity function of FRI and FRII galaxies. The high energy emission from non-beamed AGNs is instead determined by popular synthesis models of the observed X-ray background. Finally, the EGB is fit by adding a third component arising from pion decay in cosmic rays, assuming that such component is closely linked to the cosmic star formation history. We find that blazars dominate at energies \geq 10 GeV, for E \leq 0.2 GeV the main contribution is from non-beamed AGNs, while cosmic rays are required in between. Because of absorption due to interaction of gamma-rays with the extra-galactic background light, our model falls short at the highest energies probed by LAT, (\geq 70 GeV), leaving room to a possible contribution from dark matter particle annihilation. As an example, a particle of mass ~ 0.5 TeV and cross section <{\sigma}v>~ 5*10^{-26} cm^3 s^-1 can accomodate the data.

We study the implication of the collapsar model for Long Gamma-Ray Bursts (LGRBs) on the metallicity properties of the host galaxies, by combining high-resolution N-body simulations with semi-analytic models of galaxy formation. The cosmological model that we use reproduces the Fundamental Metallicity Relation recently discovered for the SDSS galaxies, whereby the metallicity decreases with increasing Star Formation Rate for galaxies of a given stellar mass. We select host galaxies housing pockets of gas-particles, young and with different thresholds in metallicities, that can be sites of LRGB events, according to the collapsar model. The simulated samples are compared with 18 observed LGRB hosts in the aim at discriminating whether the metallicity is a primary parameter. We find that a threshold in metallicity for the LGRB progenitors, within the model galaxies, is not necessary in order to reproduce the observed distribution of host metallicities. The low metallicities of observed LGRB hosts is a consequence of the high star formation environment. The star formation rate appears to be the primary parameter to generate a burst event. Finally, we show that only a few LGRBs are observed in massive, highly extincted galaxies, while these galaxies are expected to produce many such events. We identify these missing events with the fraction of dark LGRBs.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

We investigate the metallicity properties of host galaxies of long Gamma-ray Bursts (GRBs) in the light of the Fundamental Metallicity Relation (FMR), the tight dependence of metallicity on mass and SFR recently discovered for SDSS galaxies with stellar masses above 10^9.2 Msun. As most of the GRB hosts have masses below this limit, the FMR can only be used after an extension towards lower masses. At this aim, we study the FMR for galaxies with masses down to ~10^8.3 Msun, finding that the FMR does extend smoothly at lower masses, albeit with a much larger scatter. We then compare the resulting FMR with the metallicity properties of 18 host galaxies of long GRBs. While the GRB hosts show a systematic offset with respect to the mass-metallicity relation, they are fully consistent with the FMR. This shows that the difference with the mass-metallicity relation is due to higher than average SFRs, and that GRBs with optical afterglows do not preferentially select low-metallicity hosts among the star-forming galaxies. The apparent low metallicity is therefore a consequence of the occurrence of long GRB in low mass, actively star-forming galaxies, known to dominate the current cosmic SFR.

Recent observations have gathered a considerable sample of high redshift galaxy candidates and determined the evolution of their luminosity function (LF). To interpret these findings, we use cosmological SPH simulations including, in addition to standard physical processes, a detailed treatment of the Pop III-Pop II transition in early objects. The simulated high-z galaxies match remarkably well the amplitude and slope of the observed LF in the redshift range 5<z<10. The LF shifts towards fainter luminosities with increasing redshift, while its faint-end slope keeps an almost constant value, \alpha ~-2. The stellar populations of high-z galaxies have ages of 100-300 (40-130) Myr at z=5 (z=7-8), implying an early (z>9.4) start of their star formation activity; the specific star formation rate is almost independent of galactic stellar mass. These objects are enriched rapidly with metals and galaxies identified by HST/WFC3 (M_UV < -18) show metallicities ~0.1 Zsun even at z=7-8. Most of the simulated galaxies at z~7 (noticeably the smallest ones) are virtually dust-free, and none of them has an extinction larger than E(B-V) = 0.01. The bulk (50%) of the ionizing photons is produced by objects populating the faint-end of the LF (M_UV < -16), which JWST will resolve up to z=7.3. PopIII stars continue to form essentially at all redshifts; however, at z=6 (z=10) the contribution of Pop III stars to the total galactic luminosity is always less than 5% for M_UV < -17 (M_UV < -16). The typical high-z galaxies closely resemble the GRB host galaxy population observed at lower redshifts, strongly encouraging the use of GRBs to detect the first galaxies.

Gamma-Ray Bursts (GRBs) are the most violent explosions in the Universe. Long duration GRBs are associated with the collapse of massive stars, rivaling their host galaxies in luminosity. The discovery of the most distant spectroscopically confirmed object in the Universe, GRB090423, opened a new window on the high redshift Universe, making it possible to study the cosmic 're-ionization' epoch and the preceding dark ages, as well as the generation of the first stars (Population III) using GRBs. Obviously this enables a wealth of new studies using the near infrared (nIR) characteristics of GRB afterglows. Here we explore a different path, focusing on the next generation of X-ray missions with large area focusing telescopes and fast repointing capabilities. We found that X-ray data can complement nIR observations and for the brightest GRBs can provide an accurate and independent redshift determination. Metallicity studies can also be carried out profitably once the redshift is known. Finally we discuss observational signatures of GRBs arising from Population III stars in the X-ray band.

Gamma-ray bursts have been proved to be detectable up to distances much larger than any other astrophysical object, providing the most effective way, complementary to ordinary surveys, to study the high redshift universe. To this end, we present here the results of an observational campaign devoted to the study of the high-z GRB 090205. We carried out optical/NIR spectroscopy and imaging of GRB 090205 with the ESO-VLT starting from hours after the event up to several days later to detect the host galaxy. We compared the results obtained from our optical/NIR observations with the available Swift high-energy data of this burst. Our observational campaign led to the detection of the optical afterglow and host galaxy of GRB 090205 and to the first measure of its redshift, z=4.65. Similar to other, recent high-z GRBs, GRB 090205 has a short duration in the rest-frame with T_{90,rf}=1.6 s, which suggests the possibility that it might belong to the short GRBs class. The X-ray afterglow of GRB 090205 shows a complex and interesting behaviour with a possible rebrightening at 500-1000s from the trigger time and late flaring activity. Photometric observations of the GRB 090205 host galaxy argue in favor of a starburst galaxy with a stellar population younger than ~ 150 Myr. Moreover, the metallicity of Z > 0.27 Z_Sun derived from the GRB afterglow spectrum is among the highest derived from GRB afterglow measurement at high-z, suggesting that the burst occurred in a rather enriched envirorment. Finally, a detailed analysis of the afterglow spectrum shows the existence of a line corresponding to Lyman-alpha emission at the redshift of the burst. GRB 090205 is thus hosted in a typical Lyman-alpha emitter (LAE) at z=4.65. This makes the GRB 090205 host the farthest GRB host galaxy, spectroscopically confirmed, detected to date.

The approach of Observational Astronomy is mainly aimed at the construction of larger aperture telescopes, more sensitive detectors and broader wavelength coverage. Certainly fruitful, this approach turns out to be not completely fulfilling the needs when phenomena related to the formation of black holes (BH), neutron stars (NS) and relativistic stars in general are concerned. Recently, mainly through the Vela, Beppo-SAX and Swift satellites, we reached a reasonable knowledge of the most violent events in the Universe and of some of the processes we believe are leading to the formation of black holes (BH). We plan to open a new window of opportunity to study the variegated physics of very fast astronomical transients, particularly the one related to extreme compact objects. The innovative approach is based on three cornerstones: 1) the design (the conceptual design has been already completed) of a 3m robotic telescope and related focal plane instrumentation characterized by the unique features: "No telescope points faster"; 2) simultaneous multi-wavelengths observations (photometry, spectroscopy o\& polarimetry); 3) high time resolution observations. The conceptual design of the telescope and related instrumentation is optimized to address the following topics: High frequency a-periodic variability, Polarization, High z GRBs, Short GRBs, GRB-Supernovae association, Multi-wavelengths simultaneous photometry and rapid low dispersion spectroscopy. This experiment will turn the "exception" (like the optical observations of GRB 080319B) to "routine".

Short Gamma-Ray Bursts (SGRBs) are expected to form from the coalescence of compact binaries, either of primordial origin or from dynamical interactions in globular clusters. In this paper, we investigate the possibility that the offset and afterglow brightness of a SGRB can help revealing the origin of its progenitor binary. We find that a SGRB is likely to result from the primordial channel if it is observed within 10 kpc from the center of a massive galaxy and shows a detectable afterglow. The same conclusion holds if it is 100 kpc away from a small, isolated galaxy and shows a weak afterglow. On the other hand, a dynamical origin is suggested for those SGRBs with observable afterglow either at a large separation from a massive, isolated galaxy or with an offset of 10-100 kpc from a small, isolated galaxy. We discuss the possibility that SGRBs from the dynamical channel are hosted in intra-cluster globular clusters and find that GRB 061201 may fall within this scenario.

We discuss expectations for the absorption of high-energy gamma-rays by gamma-gamma pair production with intergalactic radiation fields (IRFs) at very high redshifts (z~5-20), and the prospects thereof for probing the cosmic reionization era. For the evolving IRF, a semi-analytical model incorporating both Population II and Population III stars is employed, which is consistent with a wide variety of existing high-z observations including QSO spectral measurements, WMAP Thomson depth constraints, near-IR source count limits, etc. We find that the UV IRF below the Lyman edge energy with intensities in the range of a few times 10^{-19} erg cm^{-2} s^{-1} Hz^{-1} sr^{-1} can cause appreciable attenuation above ~12 GeV at z~5, down to ~6-8 GeV at z>~8-10. This may be observable in the spectra of blazars or gamma-ray bursts by the Fermi Gamma-ray Space Telescope or next generation facilities such as the Cherenkov Telescope Array, Advanced Gamma-ray Imaging System or 5@5, providing invaluable insight into early star formation and cosmic reionization.

We investigate the cosmic reionization history by comparing semi-analytical models of the Lya forest with observations of high-z quasars and gamma ray bursts absorption spectra. In order to constrain the reionization epoch z_rei, we consider two physically motivated scenarios in which reionization ends either early (ERM, z_rei>= 7) or late (LRM, z_rei~6). We analyze the transmitted flux in a sample of 17 quasars spectra at 5.7<z<6.4 and in the spectrum of the gamma ray burst 050904 at z=6.3, studying the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with our models, we conclude that current observational data do not require any sudden change in the ionization state of the IGM at z~6, favouring indeed a highly ionized Universe at these epochs, as predicted by the ERM. Moreover, we test the predictions of this model through Lya emitters observations, finding that the ERM provide a good fit to the evolution of the luminosity function of Lya emitting galaxies in the redshift range z=5.7-6.5. The overall result points towards an extended reionization process which starts at z>=11 and completes at z_rei>=7, in agreement with the recent WMAP5 data.

Gamma-ray bursts (GRBSs) are produced by rare types of massive stellar explosions. Their rapidly fading afterglows are often bright enough at optical wavelengths, that they are detectable up to cosmological distances. Hirtheto, the highest known redshift for a GRB was z=6.7, for GRB 080913, and for a galaxy was z=6.96. Here we report observations of GRB 090423 and the near-infrared spectroscopic measurement of its redshift z=8.1^{+0.1}_{-0.3}. This burst happened when the Universe was only ~4% of its current age. Its properties are similar to those of GRBs observed at low/intermediate redshifts, suggesting that the mechanisms and progenitors that gave rise to this burst about 600 million years after the Big Bang are not markedly different from those producing GRBs ~10 billion years later.

Combining cosmological SPH simulations with a previously developed Lyman Alpha production/transmission model and the Early Reionization Model (ERM, reionization ends at redshift z~7), we obtain Lyman Alpha and UV Luminosity Functions (LFs) for Lyman Alpha Emitters (LAEs) for redshifts between 5.7 and 7.6. Matching model results to observations at z~5.7 requires escape fractions of Lyman Alpha, f_alpha=0.3, and UV (non-ionizing) continuum photons, f_c=0.22, corresponding to a color excess, E(B-V)=0.15. We find that (i) f_c increases towards higher redshifts, due the decreasing mean dust content of galaxies, (ii) the evolution of f_alpha/f_c hints at the dust content of the ISM becoming progressively inhomogeneous/clumped with decreasing redshift. The clustering photoionization boost is important during the initial reionization phases but has little effect on the Lyman Alpha LF for a highly ionized IGM. Halo (stellar) masses are in the range 10.0 < \log M_h < 11.8 (8.1 < \log M_* < 10.4) with M_h \propto M_*^{0.64}. The star formation rates are between 3-120 solar masses per year, mass-weighted mean ages are greater than 20 Myr at all redshifts, while the mean stellar metallicity increases from Z=0.12 to 0.22 solar metallicity from z~7.6 to z~5.7; both age and metallicity positively correlate with stellar mass. The brightest LAEs are all characterized by large star formation rates and intermediate ages (~200 Myr), while objects in the faint end of the Lyman Alpha LF show large age and star formation rate spreads. With no more free parameters, the Spectral Energy Distributions of three LAE at z~5.7 observed by Lai et al. (2007) are well reproduced by an intermediate age (182-220 Myr) stellar population and the above E(B-V) value.

Aims. With an observed and rest-frame duration of < 2s and < 0.5s, respectively, GRB090426 could be classified as a short GRB. The prompt detection, both from space and ground-based telescopes, of a bright optical counterpart to this GRB offered a unique opportunity to complete a detailed study. Methods. Based on an extensive ground-based observational campaign, we obtained the spectrum of the optical afterglow of GRB090426, measuring its redshift and obtaining information about the medium in which the event took place. We completed follow-up observation of the afterglow optical light curve down to the brightness level of the host galaxy that we firmly identified and studied. We also retrieved and analyzed all the available high-energy data of this event, and compared the results with our findings in the optical. This represents one of the most detailed studies of a short-duration event presented so far. Results. The time properties qualify GRB090426 as a short burst. In this case, its redshift of z = 2.61 would be the highest yet found for a GRB of this class. On the other hand, the spectral and energy properties are more similar to those of long bursts. LBT late-time deep imaging identifies a star-forming galaxy at a redshift consistent with that of the GRB. The afterglow lies within the light of its host and shows evidence of local absorption.

In this paper we examine the effect of X-ray and Lyalpha photons on the intergalactic medium temperature. We calculate the photon production from a population of stars and micro-quasars in a set of cosmological hydrodynamic simulations which self-consistently follow the dark matter dynamics, radiative processes as well as star formation, black hole growth and associated feedback processes. We find that, (i) IGM heating is always dominated by X-rays unless the Lyalpha photon contribution from stars in objects with mass M<10^8 Msun becomes significantly enhanced with respect to the X-ray contribution from BHs in the same halo (which we do not directly model). (ii) Without overproducing the unresolved X-ray background, the gas temperature becomes larger than the CMB temperature, and thus an associated 21 cm signal should be expected in emission, at z<11.5. We discuss how in such a scenario the transition redshift between a 21 cm signal in absorption and in emission could be used to constraint BHs accretion and associated feedback processes.

Gamma-ray bursts are cosmological sources emitting radiation from the gamma-rays to the radio band. Substantial observational efforts have been devoted to the study of GRBs during the prompt phase, i.e. the initial burst of high-energy radiation, and during the longer-lasting afterglows. In spite of many successes in interpreting these phenomena there are still several open key questions about the fundamental emission processes, their energetics and the environment. Moreover, independently of their modeling, GRB spectra are remarkably simple, being satisfactorily fitted with power-laws, and therefore offer a very valuable tool to probe the extragalactic background light distribution affecting all high-energy observations of cosmological sources. Observations carried out with Cherenkov telescopes, as MAGIC, can be fundamental for all these scientific topics. GRB080430, being at a rather moderate redshift, z~0.76, and well-studied in the optical, although observed only a few hours after the high-energy event, is a good test case to evaluate the perspective for late-afterglow observations with ground based GeV/TeV observatories.

Short GRBs are commonly thought to originate from the merging of double compact object binaries but direct evidence for this scenario is still missing. Optical observations of short GRBs allow us to measure redshifts, firmly identify host galaxies, characterize their properties, and accurately localize GRBs within them. Multiwavelength observations of GRB afterglows provide useful information on the emission mechanisms at work. These are all key issues that allow one to discriminate among different models of these elusive events. We carried out photometric observations of the short/hard GRB 051227, GRB 061006, and GRB 071227 with the ESO-VLT starting from several hours after the explosion down to the host galaxy level several days later. For GRB 061006 and GRB 071227 we also obtained spectroscopic observations of the host galaxy. We compared the results obtained from our optical observations with the available X-ray data of these bursts. For all the three above bursts, we discovered optical afterglows and firmly identified their host galaxies. About half a day after the burst, the optical afterglows of GRB 051227 and GRB 061006 present a decay significatly steeper than in the X-rays. In the case of GRB 051227, the optical decay is so steep that it likely indicates different emission mechanisms in the two wavelengths ranges. The three hosts are blue, star forming galaxies at moderate redshifts and with metallicities comparable to the Solar one. The projected offsets of the optical afterglows from their host galaxies centers span a wide range, but all afterglows lie within the light of their hosts and present evidence for local absorption in their X-ray spectra. We discuss our findings in light of the current models of short GRB progenitors.

The first structures in the Universe formed at z>7, at higher redshift than all currently known galaxies. Since GRBs are brighter than other cosmological sources at high redshift and exhibit simple power-law afterglow spectra that is ideal for absorption studies, they serve as powerful tools for studying the early universe. New facilities planned for the coming decade will be able to obtain a large sample of high-redshift GRBs. Such a sample would constrain the nature of the first stars, galaxies, and the reionization history of the Universe.

We compute the luminosity function (LF) and the formation rate of long gamma ray bursts (GRBs) by fitting the observed differential peak flux distribution obtained by the BATSE experiment in two different scenarios: i) the GRB luminosity evolves with redshift and ii) GRBs form preferentially in low-metallicity environments. In both cases, model predictions are consistent with the Swift number counts and with the number of detections at z>2.5 and z>3.5. To discriminate between the two evolutionary scenarios, we compare the model results with the number of luminous bursts (i.e. with isotropic peak luminosity in excess of 10^53 erg s^-1) detected by Swift in its first three years of mission. Our sample conservatively contains only bursts with good redshift determination and measured peak energy. We find that pure luminosity evolution models can account for the number of sure identifications. In the case of a pure density evolution scenario, models with Z_th>0.3 Zsun are ruled out with high confidence. For lower metallicity thresholds, the model results are still statistically consistent with available lower limits. However, many factors can increase the discrepancy between model results and data, indicating that some luminosity evolution in the GRB LF may be needed also for such low values of Z_th. Finally, using these new constraints, we derive robust upper limits on the bright-end of the GRB LF, showing that this cannot be steeper than ~2.6.

We provide measurements of the neutral hydrogen fraction xHI at z~6, by comparing semi-analytical models of the Lyalpha forest with observations of high-z quasars and Gamma Ray Bursts absorption spectra. We analyze the transmitted flux in a sample of 17 QSOs spectra at 5.74<zem<6.42 studying separately the narrow transmission windows (peaks) and the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with our models, we conclude that xHI evolves smoothly from 10^{-4.4} at z=5.3 to 10^{-4.2} at z=5.6, with a robust upper limit xHI<0.36 at z=6.3. We show the results of the first-ever detected transverse proximity effect in the HI Lyalpha forest, produced by the HII region of the faint quasar RD J1148+5253 at z=5.70 intervening along the LOS of SDSS J1148+5251 at z=6.42. Moreover, we propose a novel method to study cosmic reionization using absorption line spectra of high-redshift GRBs afterglows. We show that the time evolution and the statistics of gaps in the observed spectra represent exquisite tools to discriminate among different reionization models. By applying our methods to GRB 050904 detected at z=6.29, we show that the observation of this burst provides strong indications of a highly ionized intergalactic medium at z~6, with an estimated mean neutral hydrogen fraction xHI=6.4\pm 0.3\times 10^{-5} along that line of sight.

We compute the luminosity function (LF) and the formation rate of long gamma ray bursts (GRBs) in three different scenarios: i) GRBs follow the cosmic star formation and their LF is constant in time; ii) GRBs follow the cosmic star formation but the LF varies with redshift; iii) GRBs form preferentially in low-metallicity environments. We then test model predictions against the Swift 3-year data, showing that scenario i) is robustly ruled out. Moreover, we show that the number of bright GRBs detected by Swift suggests that GRBs should have experienced some sort of luminosity evolution with redshift, being more luminous in the past. Finally we propose to use the observations of the afterglow spectrum of GRBs at z>5.5 to constrain the reionization history and we applied our method to the case of GRB 050904.

We explore the effect of cosmic radiative feedback from the sources of reionization on the thermal evolution of the intergalactic medium. We find that different prescriptions for this feedback predict quite different thermal and reionization histories. In spite of this, current data can not discriminate among different reionization scenarios. We find that future observations both from 21-cm and CMB experiments can be used to break the degeneracy among model parameters provided that we will be able to remove the foreground signal at the percent (or better) level.

We propose to study cosmic reionization using absorption line spectra of high-redshift Gamma Ray Burst (GRB) afterglows. We show that the statistics of the dark portions (gaps) in GRB absorption spectra represent exquisite tools to discriminate among different reionization models. We then compute the probability to find the largest gap in a given width range [Wmax, Wmax + dW] at a flux threshold Fth for burst afterglows at redshifts 6.3 < z < 6.7. We show that different reionization scenarios populate the (Wmax, Fth) plane in a very different way, allowing to distinguish among different reionization histories. We provide here useful plots that allow a very simple and direct comparison between observations and model results. Finally, we apply our methods to GRB 050904 detected at z = 6.29. We show that the observation of this burst strongly favors reionization models which predict a highly ionized intergalactic medium at z~6, with an estimated mean neutral hydrogen fraction xHI = 6.4 \pm 0.3 \times 10^-5 along the line of sight towards GRB 050904.

We describe the scientific motivations, the mission concept and the instrumentation of SPACE, a class-M mission proposed for concept study at the first call of the ESA Cosmic-Vision 2015-2025 planning cycle. SPACE aims to produce the largest three-dimensional evolutionary map of the Universe over the past 10 billion years by taking near-IR spectra and measuring redshifts for more than half a billion galaxies at 0<z<2 down to AB~23 over 3\pi sr of the sky. In addition, SPACE will also target a smaller sky field, performing a deep spectroscopic survey of millions of galaxies to AB~26 and at 2<z<10+. These goals are unreachable with ground-based observations due to the ~500 times higher sky background. To achieve the main science objectives, SPACE will use a 1.5m diameter Ritchey-Chretien telescope equipped with a set of arrays of Digital Micro-mirror Devices (DMDs) covering a total field of view of 0.4 deg2, and will perform large-multiplexing multi-object spectroscopy (e.g. ~6000 targets per pointing) at a spectral resolution of R~400 as well as diffraction-limited imaging with continuous coverage from 0.8mum to 1.8mum.

Short-hard Gamma Ray Bursts (SGRBs) are currently thought to arise from gravitational wave driven coalescences of double neutron star systems forming either in the field or dynamically in globular clusters. For both channels we fit the peak flux distribution of BATSE SGRBs to derive the local burst formation rate and luminosity function. We then compare the resulting redshift distribution with Swift 2-year data, showing that both formation channels are needed in order to reproduce the observations. Double neutron stars forming in globular clusters are found to dominate the distribution at z<0.3, whereas the field population from primordial binaries can account for the high-z SGRBs. This result is not in contradiction with the observed host galaxy type of SGRBs.

We compute the imprints left on the CMB by two cosmic reionization models consistent with current observations but characterized by alternative radiative feedback prescriptions (suppression and filtering) resulting in a different suppression of star formation in low-mass halos. The models imply different ionization and thermal histories and 21 cm background signals. The derived Comptonization, u, and free-free distortion, y_B, parameters are below current observational limits for both models. However, the value of u = 1.69 * 10^-7 (9.65 * 10^-8) for the suppression (filtering) model is in the detectability range of the next generation of CMB spectrum experiments. Through the dedicated Boltzmann code CMBFAST, modified to include the above ionization histories, we compute the CMB angular power spectrum (APS) of the TT, TE, and EE modes. For the EE mode the differences between these models are significantly larger than the cosmic and sampling variance over the multipole range l = 5-15, leaving a good chance of discriminating between these feedback mechanisms with forthcoming/future CMB polarization experiments. The main limitations come from foreground contamination: it should be subtracted at per cent level in terms of APS, a result potentially achievable by novel component separation techniques and mapping of Galactic foreground.

We use a semi-analytical model to study the impact of reionization, and the associated radiative feedback, on galaxy formation. Two feedback models have been considered: (i) a standard prescription, according to which star formation is totally suppressed in galaxies with circular velocity below a critical threshold (model CF06) and (ii) a characterization based on the filtering scale (model G00), allowing for a gradual reduction of the gas available for star formation in low-mass galaxies. In model CF06 reionization starts at z ~ 15-20, is 85% complete by z ~ 10; at the same z, the ionized fraction is 16% in model G00. The models match SDSS constraints on the evolution of the neutral hydrogen fraction at z < 7, but predict different Thomson optical depths, tau_e = 0.1017 (CF06), and 0.0631 (G00); such values are within 1 sigma of the WMAP 3-yr determination. Both models are in remarkable good agreement with additional existing data (evolution of Lyman-limit systems, cosmic star formation history, high-z galaxy counts, IGM thermal history), which therefore cannot be used to discriminate among different feedback models. Deviations among radiative feedback prescriptions emerge when considering the expected HI 21 cm background signal, where a ~ 15 mK absorption feature in the range 75-100 MHz is present in model G00 and a global shift of the emission feature preceding reionization towards larger frequencies occurs in the same model. Single dish observations with existing or forthcoming low-frequency radio telescopes can achieve mK sensitivity, allowing the identification of these features provided that foregrounds can be accurately subtracted.

Long Gamma Ray Bursts (GRBs) constitute an important tool to study the Universe near and beyond the epoch of reionization. We delineate here the characteristics of an 'ideal' instrument for the search of GRBs at z>6-10. We find that the detection of these objects requires soft band detectors with a high sensitivity and moderately large FOV. In the light of these results, we compare available and planned GRB missions, deriving conservative predictions on the number of high-z GRBs detectable by these instruments along with the maximum accessible redshift. We show that the Swift satellite will be able to detect various GRBs at z>6, and likely at z>10 if the trigger threshold is decreased by a factor of ~2. Furthermore, we find that INTEGRAL and GLAST are not the best tool to detect bursts at z>6: the former being limited by the small FOV, and the latter by its hard energy band and relatively low sensitivity. Finally, future missions (SVOM, EDGE, but in particular EXIST) will provide a good sample of GRBs at z>6 in a few years of operation.

In this paper we have calculated the effect of Lyalpha photons emitted by the first stars on the evolution of the IGM temperature. We have considered both a standard Salpeter IMF and a delta-function IMF for very massive stars with mass 300 M_sun. We find that the Lyalpha photons produced by the stellar populations considered here are able to heat the IGM at z<25, although never above ~100 K. Stars with a Salpeter IMF are more effective as, due to the contribution from small-mass long-living stars, they produce a higher Lyalpha background. Lyalpha heating can affect the subsequent formation of small mass objects by producing an entropy floor that may limit the amount of gas able to collapse and reduce the gas clumping.We find that the gas fraction in halos of mass below ~ 5 x 10^6 M_sun is less than 50% (for the smallest masses this fraction drops to 1% or less) compared to a case without Lyalpha heating. Finally, Lyalpha photons heat the IGM temperature above the CMB temperature and render the 21cm line from neutral hydrogen visible in emission at z<15.

We compute the probability to detect long Gamma Ray Bursts (GRBs) at z>5 with Swift, assuming that GRBs form preferentially in low-metallicity environments. The model fits well both the observed BATSE and Swift GRB differential peak flux distribution and is consistent with the number of z>2.5 detections in the 2-year Swift data. We find that the probability to observe a burst at z>5 becomes larger than 10% for photon fluxes P<1 ph s^{-1} cm^{-2}, consistent with the number of confirmed detections. The corresponding fraction of z>5 bursts in the Swift catalog is ~10%-30% depending on the adopted metallicity threshold for GRB formation. We propose to use the computed probability as a tool to identify high redshift GRBs. By jointly considering promptly-available information provided by Swift and model results, we can select reliable z>5 candidates in a few hours from the BAT detection. We test the procedure against last year Swift data: only three bursts match all our requirements, two being confirmed at z>5. Other three possible candidates are picked up by slightly relaxing the adopted criteria. No low-z interloper is found among the six candidates.

We compute the luminosity function (LF) and the formation rate of long gamma ray bursts (GRBs) by fitting the observed differential peak flux distribution obtained by the BATSE satellite in three different scenarios: i) GRBs follow the cosmic star formation and their LF is constant in time; ii) GRBs follow the cosmic star formation but the LF varies with redshift; iii) GRBs form preferentially in low-metallicity environments. We find that the differential peak flux number counts obtained by BATSE and by Swift can be reproduced using the same LF and GRB formation rate, indicating that the two satellites are observing the same GRB population. We then check the resulting redshift distributions in the light of Swift 2-year data, focusing in particular on the relatively large sample of GRBs detected at z>2.5. We show that models in which GRBs trace the cosmic star formation and are described by a constant LF are ruled out by the number of high-z Swift detections. This conclusion does not depend on the redshift distribution of bursts that lack of optical identification, nor on the existence of a decline in star formation rate at z>2, nor on the adopted faint-end of the GRB LF. Swift observations can be explained by assuming that the LF varies with redshift and/or that GRB formation is limited to low-metallicity environments.

We study, by means of dedicated simulations of massive black hole build-up, the possibility to constraint the existence and nature of the AGN population at z>6 with available and planned X-ray and near infrared space telescopes. We find that X-ray deep-field observations can set important constraints to the faint-end of the AGN luminosity function at very high redshift. Planned X-ray telescopes should be able to detect AGN hosting black holes with masses down to >10^5 Msun (i.e., X-ray luminosities in excess of 10^42 erg s^-1), and can constrain the evolution of the population of massive black hole at early times (6<z<10). We find that this population of AGN should contribute substantially (~25%) to the unresolved fraction of the cosmic X-ray background in the 0.5-10 keV range, and that a significant fraction (~3-4%) of the total background intensity would remain unaccounted even after future X-ray observations. As byproduct, we compute the expected UV background from AGN at z>6 and we discuss the possible role of AGN in the reionization of the Universe at these early epochs, showing that AGN alone can provide enough ionizing photons only in the (improbable) case of an almost completely homogeneous inter-galactic medium. Finally, we show that super-Eddington accretion, suggested by the observed QSOs at z>6, must be a very rare event, confined to black holes living in the highest density peaks.

We investigate how hierarchical models for the co-evolution of the massive black hole (MBH) and AGN population can reproduce the observed faint X-ray counts. We find that the main variable influencing the theoretical predictions is the Eddington ratio of accreting sources. We compare three different models proposed for the evolution of AGN Eddington ratio, f_Edd: constant f_Edd=1, f_Edd decreasing with redshift, and f_Edd depending on the AGN luminosity, as suggested by simulations of galactic mergers including BHs and AGN feedback. We follow the full assembly of MBHs and host halos from early times to the present in a LambdaCDM cosmology. AGN activity is triggered by halo major mergers and MBHs accrete mass until they satisfy the observed correlation with velocity dispersion. We find that all three models can reproduce fairly well the total faint X-ray counts. The redshift distribution is however poorly matched in the first two models. The Eddington ratios suggested by merger simulations predicts no turn-off of the faint end of the AGN optical luminosity function at redshifts z>=1 down to very low luminosity.

We explore the nature of possible electromagnetic counterparts of coalescences of massive black hole binaries at cosmological distances detectable by the Laser Interferometer Space Antenna (LISA). An electromagnetic precursor, during the last year of gravitational wave (GW)-driven inspiral, or an afterglow within few years after coalescence, may highlight the position in the sky of galaxies hosting LISA sources. We show that observations of precursors and afterglows are mutually exclusive, depending on the mass of the primary black hole. Precursors are expected to occur in binaries where the primary (more massive) black hole is heavier than ~10^7 \Msun. They may correspond to on-off states of accretion, i.e., to a bright X-ray source decaying into quiescence before black hole coalescence, and are likely associated to disturbed galaxies showing signs of ongoing starbursts. Coalescences of lighter binaries, with masses <5x10^6 \Msun, lack of any precursor, as gas is expected to be consumed long before the GW-driven orbital decay. Such events would not be hosted by (massive) galaxies with an associated starburst, given the slow binary inspiral time compared to the typical time scale of starbursts. By contrast, coalescence, for such light binaries, is followed by an electromagnetic afterglow, i.e., an off-on accretion state rising in <20 yrs. Using a cosmological merger tree algorithm, we show that future X-ray missions such as XEUS will be able to identify, in 20 yrs operation, almost all the massive BH binary detectable by LISA, and, in only 5 yrs, all the LISA sources at z>6.

Kashlinsky et al. (2005) find a significant cosmic infrared background fluctuation excess on angular scales >50 arcsec that cannot be explained by instrumental noise or local foregrounds. The excess has been tentatively attributed to emission from primordial very massive (PopIII) stars formed <200 Myr after the Big Bang. Using an evolutionary model motivated by independent observations and including various feedback processes, we find that PopIII stars can contribute <40% of the total background intensity (\nu J_\nu ~ 1-2 nW m^-2 sr^-1 in the 0.8-8 \mum range) produced by all galaxies (hosting both PopIII and PopII stars) at z>5. The infrared fluctuation excess is instead very precisely accounted by the clustering signal of galaxies at z>5, predominantly hosting PopII stars with masses and properties similar to the present ones.

The observed near infrared background excess over light from known galaxies is commonly ascribed to redshifted radiation from early, very massive, PopIII stars. We show here that this interpretation must be discarded as it largely overpredicts the number of J-dropouts and Ly\alpha emitters in ultra deep field searches. Independently of the detailed physics of Ly\alpha line emission, J-dropouts limit the background excess fraction due to PopIII sources to be (at best) < 1/24. As alternative explanations can either be rejected (e.g. miniquasars, decaying neutrinos) or appear unlikely (zodiacal light), whereas the reality of the excess is supported by the interpretation of the angular fluctuations, the origin of this component remains very puzzling. We briefly discuss possible hints to solve the problem.

Motivated by the idea that the recently detected near-infrared (1.2-4 micron) excess over the contribution of known galaxies is due to redshifted light from the first cosmic stars (Salvaterra & Ferrara 2003), we have used the effect caused by photon-photon absorption on gamma-ray spectra of blazars to put constraints on extragalactic background light (EBL) from the optical to the far-IR bands. Our analysis is mainly based on the blazar H 1426+428, for which we assume a power-law unabsorbed spectrum. We find that an EBL model with no excess over known galaxies in the near-infrared background (NIRB) is in agreement with all the considered blazars; however, it implies a very peculiar intrinsic spectrum for H 1426+428. Additional data on the blazars 1ES1101-232, H 2356-309 and PKS 2155-304 exclude the existence of a strong NIRB excess consistent with Kelsall's model of zodiacal light subtraction (ZL); the COBE/DIRBE measurements, after Wright's model ZL subtraction, represent a firm NIRB upper limit. The constraints on the optical EBL are weaker, due to the fact that predictions from different optical EBL models are often comparable to the experimental errors. In the mid-infrared the SPITZER measurement at 24 micron gives a good fit for all the considered blazars.

Motivated by theoretical predictions that first stars were predominantly very massive, we investigate the physics of the transition from an early epoch dominated by massive Pop III stars to a later epoch dominated by familiar low-mass Pop II/I stars by means of a numerically-generated catalogue of dark matter halos coupled with a self-consistent treatment of chemical and radiative feedback. Depending on the strength of the chemical feedback, Pop III stars can contribute a substantial fraction (several percent) of the cosmic star formation activity even at moderate redshifts, z = 5. We find that the three z = 10 sources tentatively detected in NICMOS UDFs should be powered by Pop III stars, if these are massive; however, this scenario fails to reproduce the derived WMAP electron scattering optical depth. Instead, both the UDFs and WMAP constraints can be fulfilled if stars at any time form with a more standard, slightly top-heavy, Larson IMF in the range 1 Msun < M < 100 Msun.

A large population of Intermediate Mass Black Holes (IMBHs) might be produced at early cosmic times as a left over of the evolution of the very massive first stars. Accretion onto IMBHs provides a source of (re)ionizing radiation. We show that the baryon mass fraction locked into IMBHs and their growth is strongly constrained by the observed residual Soft X-ray Background (SXRB) intensity. Thus, unless they are extremely X-ray quiet, miniquasars must be quite rare and/or have a short shining phase. As a byproduct, we find that miniquasars can not be the only source of reionization and that their alleged contribution to the near infrared bands is completely negligible.

We study CMB secondary anisotropies produced by inhomogeneous reionization by means of cosmological simulations coupled with the radiative transfer code CRASH. The reionization history is consistent with the WMAP Thomson optical depth determination. We find that the signal arising from this process dominates over the primary CMB component for l > 4000 and reaches a maximum amplitude of l(l+1)C_l/2\pi ~ 1.6 x 10^{-13} on arcmin scale, i.e. l as large as several thousands. We then cross-correlate secondary CMB anisotropy maps with neutral hydrogen 21cm line emission fluctuations obtained from the same simulations. The two signals are highly anti-correlated on angular scales corresponding to the typical size of HII regions (including overlapping) at the 21cm map redshift. We show how the CMB/21cm cross-correlation can be used to: (a) study the nature of the reionization sources, (b) reconstruct the cosmic reionization history, (c) infer the mean cosmic ionization level at any redshift. We discuss the feasibility of the proposed experiment with forthcoming facilities.

The WMAP first year detection of a high redshift reionization through its imprints on CMB anisotropy T and TE mode angular power spectra calls for a better comprehension of the universe ionization and thermal history after the standard recombination. Different reionization mechanisms predict different signatures in the CMB, both in temperature and polarization anisotropies and in spectral distortions. The Planck capability to distinguish among different scenarios through its sensitivity to T, TE, and E mode angular power spectra is discussed. Perspectives open by future high sensitivity experiments on the CMB polarization anisotropy and spectrum are also presented.

The study of the Extragalactic Background Light (EBL) is crucial to understand many astrophysical problems (as the formation of first stars, the evolution of galaxies and the role of dust emission). At present, one of the most powerful ways to put constraints on EBL is represented by the study of the photon-photon absorption on gamma-ray spectra of TeV blazars. Adopting this method, we found that, if the only contribution to the optical and Near Infrared (NIR) background is given by galaxies, the spectrum of the blazar H1426+428 cannot be fitted. To reproduce the observational data of H1426+428 a Near Infrared excess with respect to galaxy counts is required, with amplitude consistent with both the Matsumoto et al. (2000) data with Kelsall's model of zodiacal light (ZL) subtraction and the DIRBE data with Wright's model of ZL subtraction. The derived constraints on the optical EBL are weaker, because the experimental errors on blazar data are still bigger than the differences among various optical EBL models. In the mid-infrared the SPITZER measurement at 24 micron provides the best fit of the blazar spectrum.

The consequences of energy injection associated with stellar activity in galaxies have been often referred to as ``feedback effects''. Assessing the detailed physics and the role of such effects are central issues in modern theories of galaxy formation and evolution. In these Lectures we present a brief review of various aspects related to feedback with particular emphasis to those occurring when the universe was only 100 Myr old, an epoch now commonly known as the the Cosmic Dawn, when the first galaxies and the stars within them brought the universe back into the visible ages.

We show that the explosion of the first supernovae can trigger low-mass star formation via gravitational fragmentation of the supernova-driven gas shell. If the shell mass does not exceed the host galaxy gas mass, all explosions with energies E_SN > 10^{51} erg can lead to shell fragmentation. However, the minimum ambient density required to induce such fragmentation is much larger, n_0 > 300 cm^{-3}, for Type II supernovae than for pair-instability ones, which can induce star formation even at lower ambient densities. The typical mass of the unstable fragments is 10^{4-7} Msun; their density is in the range 110-6x10^7 cm{-3}. Fragments have a metallicity strictly lower than 10^{-2.6} Zsun and large values of the gravitational-to-pressure force ratio f ~ 8. Based on these findings, we conclude that the second generation of stars produced by such self-propagating star formation is predominantly constituted by low-mass, long-living, extremely metal-poor (or even metal-free, if mixing is suppressed) stars. We discuss the implications of such results for Pop III star formation scenarios and for the most iron-poor halo star HE0107-5240.

At redshift z>5 Type II supernovae (SNII) are the only known dust sources with evolutionary timescales shorter than the Hubble time. We extend the model of dust formation in the ejecta of SNII by Todini & Ferrara (2001) to investigate the same process in pair-instability supernovae (PISN), which are though to arise from the explosion of the first, metal free, very massive (140-260 Msun) cosmic stars. We find that 15%-30% of the PISN progenitor mass is converted into dust, a value >10 times higher than for SNII; PISN dust depletion factors (fraction of produced metals locked into dust grains) range between 0.3 and 0.7. These conclusions depend very weakly on the mass of the PISN stellar progenitor, which instead affects considerably the composition and size distribution. For the assumed temperature evolution, grain condensation starts 150-200 days after the explosion; the dominant compounds for all progenitor masses are SiO2 and Mg2SiO4 while the contribution of amorphous carbon and magnetite grains grows with progenitor mass; typical grain sizes range between 0.001 and a few 0.1 micron and are always smaller than 1 micron. We give a brief discussion of the implications of dust formation for the IMF evolution of the first stars, cosmic reionization and the intergalactic medium.

We have studied the implications of the new generation of CMB spectrum space experiments for our knowledge of the thermal history of the Universe. The combination of two experiments with the sensitivity and the frequency coverage jointly forseen for DIMES and FIRAS II will be able to greatly change our vision of the capability of the CMB spectrum information to constrain physical processes at different cosmic ages. The limits on the energy dissipations at the all cosmic times accessible to CMB spectrum investigations (z < z_therm) could be improved by about two order of magnitudes and even dissipation processes with DE/E ~ 10^-6 could be detected and possibly accurately studied. These results are possible because such levels of accuracy on a so wide frequency range allow to remove the approximate degeneracy both between free-free and Bose-Einstein (BE) like distortions and between Comptonization and BE-like distortions. Finally, we discussed the different signatures imprinted on the CMB spectrum by some late astrophysical and particle decay models recently proposed in the literature and possibly related to the reionization of the Universe indicated by WMAP, and compared them with the sensitivity of such classes of CMB space spectrum experiments.

We show that the emission from the first, metal-free stars inside Population III objects (PopIIIs) are needed to explain the level of fluctuations in the Near Infrared Background (NIRB) recently discovered by Kashlinsky et al. (2002), at least at the shortest wavelengths. Clustering of (unresolved) Pop IIIs can in fact account for the entire signal at almost all the ~1-30 arcsec scales probed by observations in the J band. Their contribution fades away at shorter frequencies and becomes negligible in the K band. ``Normal'', highly clustered, <z>~3 galaxies undergoing intense star-formation such as those found in the Hubble Deep Fields can 'fill in' this gap and provide for the missing signal. It is in fact found that their contribution to the intensity fluctuations is the dominant one at lambda=2.17 mum, while it gradually looses importance in the H and J bands. The joint contribution from these two populations of cosmic objects is able, within the errors, to reproduce the observed power spectrum in the whole Near Infrared range on small angular scales (theta < 200" for Pop III protogalaxies). Signals on larger scales detected by other experiments instead require the presence of more local sources.

The earliest stars to form in the Universe were the first sources of light, heat and metals after the Big Bang. The products of their evolution will have had a profound impact on subsequent generations of stars. Recent studies of primordial star formation have shown that, in the absence of metals (elements heavier than helium), the formation of stars with masses 100 times that of the Sun would have been strongly favoured, and that low-mass stars could not have formed before a minimum level of metal enrichment had been reached. The value of this minimum level is very uncertain, but is likely to be between 10^{-6} and 10^{-4} that of the Sun. Here we show that the recent discovery of the most iron-poor star known indicates the presence of dust in extremely low-metallicity gas, and that this dust is crucial for the formation of lower-mass second-generation stars that could survive until today. The dust provides a pathway for cooling the gas that leads to fragmentation of the precursor molecular cloud into smaller clumps, which become the lower-mass stars.

Population III stars are now believed to contribute to the observed Near Infrared Background (NIRB) and heavy element pollution of the intergalactic medium. Here we show that a Pop III contribution to the primordial He abundance consistent with NIRB constraints, might mask a lower value for the Big Bang Nucleosynthesis (BBN) He abundance.

We analyse the implications of future observations of the CMB absolute temperature at centimeter and decimeter wavelengths, necessary to complement the accurate COBE/FIRAS data. Our analysis shows that forthcoming ground and balloon measures will allow a better understanding of free-free distortions but will not be able to significantly improve the constraints already provided by the FIRAS data on the possible energy exchanges in the primeval plasma. The same holds even improving the sensitivity up to ~10 times. Thus, we have studied the impact of very high quality data, such those in principle achievable with a space experiment like DIMES planned to measure the CMB absolute temperature at 0.5 - 15 cm with a sensitivity of ~0.1 mK, close to that of FIRAS. Such high quality data would improve by a factor ~50 the FIRAS results on the fractional energy exchanges associated to dissipation processes possibly occurred in a wide range of cosmic epochs, at intermediate and high redshifts (y_h > 1). The energy dissipation epoch could be also significantly constrained, while the limits on the energy possibly dissipated at late epochs can be improved by a factor ~2. In addition, such measures can provide an independent and very accurate cross-check of FIRAS calibration. Finally, a DIMES-like experiment will be able to provide indicative independent estimates of the baryon density.

The redshifted light of the first (Pop III) stars might substantially contribute to the near infrared background (NIRB). By fitting recent data with models including up-to-date Pop III stellar spectra, we find that such stars can indeed account for the whole NIRB residual (i.e. after `normal' galaxy contribution subtraction) if the high redshift star formation efficiency is f = 10%-50%, depending on the IMF (the top-heaviest requiring lowest efficiency) and on the unknown galaxy contribution in the L band (our models, however, suggest it to be negligible). Such epoch of Pop III star formation ends in all models by z_end ~ 8.8, with a hard limit z_end < 9 set by J band observations. To prevent an associated IGM over-enrichment with heavy elements compared to observed levels in the IGM, pair-instability supernovae must be the dominant heavy element sources. Alternative explanations must break the light-metal production link by advocating very massive stars (M > 260 M_sun) locking their nucleosynthetic products in the compact remnant or by postulating an extremely inhomogeneous metal enrichment of the Ly alpha forest. We discuss these possibilities in detail along with the uncertainties related to the adopted zodiacal light model.

We have compared the absolute temperature data of the CMB spectrum with models for CMB spectra distorted by a single or two heating processes at different cosmic times. The constraints on the fractional energy injected in the radiation field, DE/E, are mainly provided by the FIRAS instrument aboard the COBE satellite. Under the hypothesis that two heating processes have occurred at different epochs, the limits on DE/E are relaxed by a factor 2 both for the earlier and the later process with respect to the case in which a single energy injection in the thermal history of the universe is considered. In general, the constraints on DE/E are weaker for early processes than for relatively late processes, because of the wavelength coverage of FIRAS data. We considered also the FIRAS calibration as revised by Battistelli et al. 2000, that, in the case of the favourite calibrator emissivity law proposed by the authors, implies significant deviations from a planckian spectrum. An astrophysical explanation of this, although intriguing, seems difficult, both in terms of CMB spectral distortions and in terms of a relevant millimetric foreground. Future precise measurements at longer wavelengths as well as current and future CMB anisotropy space missions will provide independent, direct or indirect, cross checks.

We present the complete set of the measures of the CMB absolute temperature on the basis on the works published in the literature, recognize the main causes of error with the aim of estimating separately the magnitude of the systematic and statistical errors, by focussing in particular on the most recent observations. The main purpose of this work is to create a complete and reasoned database of CMB absolute temperatures. The simple database format permits a reading of the data through any text editor or through programs for data handling. This is the first step for a statistical analysis of CMB spectrum data and their comparison with the theoretical predictions for the distorted spectra in order to derive constraints on physical processes at very high redshifts.