Effective absorbing column density in the gamma-ray burst afterglow X-ray spectra

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.Comment: MNRAS accepted. 6 figures. 3 table

A complete sample of bright Swift short Gamma-Ray Bursts

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.Comment: 17 pages, 8 figures, 7 tables. Accepted for publication in MNRA

GRB orphan afterglows in present and future radio transient surveys

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.Comment: (10 pages, 5 figures, 1 table) Accepted for publication by PAS

Unveiling the population of orphan Gamma Ray Bursts

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.Comment: 9 pages, 4 figures, 2 tables. Accepted for publication by Astronomy and Astrophysic

Optical and X-ray Rest-frame Light Curves of the BAT6 sample

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.Comment: Accepted for publication in A&A: 10 pages, 5 figures, 2 table

Short GRBs at the dawn of the gravitational wave era

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

Co-Occurring Psychiatric and Substance Use Disorders: Clinical Survey Among a Rural Cohort of Italian Patients

Purpose: Dual diagnosis (DD) is the co-occurrence of both a mental illness and a substance use disorder (SUD). Lots of studies have analysed the integrated clinical approach, which involves both psychiatry and toxicology medical experts. The purpose of this study is to analyse the socio-demographic characteristics and treatment strategies of patients with DD in a rural area of Italy. Patients and Methods: Clinical data of 750 patients were collected in 2016 through the analysis of health plan records. Results: The rate of co-occurring disorders is highly variable among people with SUD. In the considered area, patients with DD are 24%, of these only 46.1% have been treated with an integrated clinical program. Moreover, this percentage is further reduced (35.8%) if only patients with heroin use disorder are considered. Conclusion: A comprehensive revision of DD treatment is needed, especially for people suffering from heroin use disorder and living in remote areas. Meticulous data analysis from other addiction health services of rural areas could be necessary to identify a science-based clinical intervention

Dust extinction for an unbiased sample of GRB afterglows

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.Comment: 18 pages, 7 figures, 10 tables, MNRAS, in pres