Search

Gamma-ray bursts (GRBs) are intense pulses of high-energy emission associated with massive stars' death or compact objects' coalescence. Their multi-wavelength observations help verify the reliability of the standard fireball model. We analyze 14 GRBs observed contemporaneously in gamma-rays by the \textit{Fermi} Large Area Telescope (LAT), in X-rays by the \textit{Swift} Telescope, and in the optical bands by \textit{Swift} and many ground-based telescopes. We study the correlation between the spectral and temporal indices using closure relations according to the synchrotron forward-shock model in the stratified medium ($n \propto r^{-k}$) with $k$ ranging from 0 to 2.5. We find that the model without energy injection is preferred over the one with energy injection in all the investigated wavelengths. In gamma-rays, we only explored the $\nu > $ max\{$\nu_c,\nu_m$\} (SC/FC) cooling condition (where $\nu_c$ and $\nu_m$ are the cooling and characteristic frequencies, namely the frequencies at the spectral break). In the X-ray and optical bands, we explored all the cooling conditions, including $\nu_m < \nu < \nu_c$ (SC), $\nu_c < \nu < \nu_m$ (FC), and SC/FC, and found a clear preference for SC for X-rays and SC/FC for optical. Within these cooling conditions, X-rays exhibit the highest rate of occurrence for the density profile with $k = 0$, while the optical band has the highest occurrence for $k$ = 2.5 when considering no energy injection. Although we can pinpoint a definite environment for some GRBs, we find degeneracies in other GRBs., Comment: 34 pages, 12 figures (two have 2 panels, two have a single panel, one has 14 panels, one has 17 panels, one has 10 panels, four has 15 panels, one has 6 panels), 10 tables (six have single panel, two have 2 panels, two have 3 panels). Accepted for publication in ApJ

$\Lambda$CDM provides a leading framework in the interpretation of modern cosmology. Nevertheless, the scientific community still struggles with many open problems in cosmology. Among the most noticeable ones, the tension in the Hubble constant $H_0$ is particularly intriguing, prompting a wide range of possible solutions. In the present work, the flat scale-free cosmology ($S$CDM) of Maeder (2017) is tested for $H_0$ tension in fits to the Pantheon sample of Supernovae Ia. The Pantheon sample is a collection of 1048 SNe Ia, which formally defines $H_0=H(0)$ by extrapolation to redshift zero of data over positive redshifts $z>0$. Here, we consider $H_{0,k}$ in fits of $S$CDM over $k$ equally sized bins of sub-samples of mean redshift $z_k$. To quantify a trend in $H_{0,k}$ with $z_k$, the results are fit by $f(z)=H'_0/(1+z)^\alpha$ in the two parameters $H'_0$ and $\alpha$. This approach tests for model imperfections or data biases by trends that are inconsistent with zero. Our findings show a decreasing trend inconsistent with zero at $5.3 \sigma$ significance, significantly more so than in $\Lambda$CDM. These results are further confirmed in Pantheon+. It appears that a solution to $H_0$ tension is to be found in models with a deceleration parameter $q_0

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide an unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium. We present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of 1 < p < 2 and 2 $\leq$ p. We show that the SSC scenario plays a relevant role in the radiative parameter $\epsilon$, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies $\geq$ 100 MeV in a sample of nine bursts from the second Fermi/LAT GRB catalog that exhibited temporal and spectral indices with $\geq$ 1.5 and $\approx$ 2, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modeled the multi-wavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor and the mechanism responsible for explaining the energetic GeV photons., Comment: Accepted in MNRAS. 26 pages, 7 figures, 6 tables

We present the largest optical photometry compilation of Gamma-Ray Bursts (GRBs) with redshifts ($z$). We include 64813 observations of 535 events (including upper limits) from 28 February 1997 up to 18 August 2023. We also present a user-friendly web tool \textit{grbLC} which allows users the visualization of photometry, coordinates, redshift, host galaxy extinction, and spectral indices for each event in our database. Furthermore, we have added a Gamma Ray Coordinate Network (GCN) scraper that can be used to collect data by gathering magnitudes from the GCNs. The web tool also includes a package for uniformly investigating colour evolution. We compute the optical spectral indices for 138 GRBs for which we have at least 4 filters at the same epoch in our sample and craft a procedure to distinguish between GRBs with and without colour evolution. By providing a uniform format and repository for the optical catalogue, this web-based archive is the first step towards unifying several community efforts to gather the photometric information for all GRBs with known redshifts. This catalogue will enable population studies by providing light curves (LCs) with better coverage since we have gathered data from different ground-based locations. Consequently, these LCs can be used to train future LC reconstructions for an extended inference of the redshift. The data gathering also allows us to fill some of the orbital gaps from Swift in crucial points of the LCs, e.g., at the end of the plateau emission or where a jet break is identified., Comment: 20 pages, 16 figures, 2 tables. Submitted to MNRAS, this version matches the third revision. The Online Materials and data will be available after the publication

Gamma-ray bursts (GRBs) are one of the most promising transient events for studying multi-wavelength observations in extreme conditions. Observation of GeV photons from bursts would provide crucial information on GRB physics, including the off-axis emission. The Second Gamma-ray Burst Catalog (2FLGC) was announced by the Fermi Large Area Telescope (Fermi-LAT) Collaboration. This catalog includes 29 bursts with photon energy higher than 10 GeV. While the synchrotron forward-shock model has well explained the afterglow data of GRBs, photon energies greater than 10 GeV are very difficult to interpret within this framework. To study the spectral and temporal indices of those bursts described in 2FLGC, Fraija et al. (2022a) proposed the closure relations (CRs) of the synchrotron self-Compton (SSC) emitted from an on-axis jet which decelerates in stellar-wind and the constant-density medium. In this paper, we extend the CRs of the SSC afterglow from an on-axis scenario to an off-axis, including the synchrotron afterglow radiation that seems off-axis. In order to investigate the spectral and temporal index evolution of those bursts reported in 2FLGC, we consider the hydrodynamical evolution with energy injection in the adiabatic and radiative regime for an electron distribution with a spectral index of $1

The LambdaCDM model is the most commonly accepted framework in modern cosmology. However, the local measurements of the Hubble constant, H0, via the Supernovae Type Ia (SNe Ia) calibrated on Cepheids provide a value which is in significant disagreement, from 4 to 6 sigma, with the value of H0 inferred from the Cosmic Microwave Background (CMB) observed by Planck. This disagreement is the so-called Hubble constant tension. To find out the reason for this discrepancy, we analyze the behaviour of the H0 in the Pantheon sample of SNe Ia through a binning approach: we divide the Pantheon into 3 and 4 bins ordered with redshift (z), and for each of them, we estimate the H0. After the H0 estimation, we fit the H0 values with a decreasing function of z, finding out that H0 undergoes a slow decreasing trend compatible with the evolution scenario in 2.0 sigma. [...] Together with SNe Ia, more astrophysical probes such as quasars (QSO) [...] and Gamma-Ray Bursts (GRBs) [...], are needed to tackle the H0 tension. In the realm of GRB-cosmology, one of the most promising correlations is the fundamental plane relation [...]. In the context of applying this relation as a cosmological tool, we also compute how many GRBs must be gathered to reach the same precision as the SNe Ia. Since we are about two decades away from reaching such precision, we also attempt to find additional correlations for the GRBs associated with SNe Ibc that could be exploited to standardize the class of GRB-SNe Ibc in the future. We find a hint of a correlation between the GRBs' end-of-plateau optical luminosity and the SNe's rest-frame peak time, suggesting that the GRBs with the most luminous optical plateau emission are associated with SNe with the most delayed peaks in their light curves. So far, it is the fundamental plane relation to be the most promising candle for exploring the high-z universe., Comment: 8 pages, 3 figures with 2, 1, and 2 panels respectively. Accepted for publication as a proceeding for The 38th International Cosmic Ray Conference (ICRC2023) through Proceedings of Science (PoS)

The longstanding search for the cosmological model that best describes the Universe has been made more intriguing since the recent discovery of the Hubble constant, $H_{0}$, tension observed between the value of $H_{0}$ from the Cosmic Microwave Background and from type Ia supernovae (SNe Ia). Hence, the commonly trusted flat $\Lambda$CDM model is under investigation. In this scenario, cosmography is a very powerful technique to investigate the evolution of the Universe without any cosmological assumption, thus revealing tensions between observational data and predictions from cosmological models in a completely model-independent way. We here employ a robust cosmographic technique based on an orthogonal logarithmic polynomial expansion of the luminosity distance to fit quasars (QSOs) alone and QSOs combined with Gamma-Ray Bursts (GRBs), SNe Ia, and Baryon Acoustic Oscillations. To apply QSOs and GRBs as probes we use, respectively, the Risaliti-Lusso relation between ultraviolet and X-ray luminosities and the ``Dainotti GRB 3D relation" among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity. We also correct QSOs and GRBs for selection biases and redshift evolution and we employ both the traditional Gaussian likelihood and the newly discovered best-fit likelihoods for each probe investigated. This comprehensive analysis reveals a strong tension ($> 4 \, \sigma$) between our data sets and the flat $\Lambda$CDM model proving the power of both the cosmographic approach and high-redshift sources, such as QSOs and GRBs, which can probe the Universe at early epochs., Comment: 13 pages, 4 figures, 4 tables. Comments are welcome. Accepted for publication in MNRAS

The recent $\sim 4 \, \sigma$ Hubble constant, $H_{0}$, tension is observed between the value of $H_{0}$ from the Cosmic Microwave Background (CMB) and Type Ia Supernovae (SNe Ia). It is a decade since this tension is excruciating the modern astrophysical community. To shed light on this problem is key to consider probes at intermediate redshifts between SNe Ia and CMB and reduce the uncertainty on $H_0$. Toward these goals, we fill the redshift gap by employing Gamma-Ray Bursts (GRBs) and Quasars (QSOs), reaching $z=9.4$ and $z=7.6$, respectively, combined with Baryonic Acoustic Oscillations (BAO) and SNe Ia. To this end, we employ the ``Dainotti GRB 3D relation" among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity, and the ``Risaliti-Lusso" QSO relation between ultraviolet and X-ray luminosities. We inquire the commonly adopted Gaussianity assumption on GRBs, QSOs, and BAO. With the joint sample, we fit the flat $\Lambda$ Cold Dark Matter model with both the Gaussian and the newly discovered likelihoods. We also investigate the impact of the calibration assumed for \textit{Pantheon} and \textit{Pantheon +} SNe Ia on this analysis. Remarkably, we show that only GRBs fulfill the Gaussianity assumption. We achieve small uncertainties on the matter density parameter $\Omega_M$ and $H_0$. We find $H_0$ values compatible within 2 $\sigma$ with the one from the Tip of the Red Giant Branch. Finally, we show that the cosmological results are heavily biased against the arbitrary calibration choice for SNe Ia., Comment: 17 pages, 6 figures, 1 table, accepted in MNRAS

While the dominant radiation mechanism gamma-ray bursts (GRBs) remains a question of debate, synchrotron emission is one of the foremost candidates to describe the multi-wavelength afterglow observations. As such, it is expected that GRBs should present some degree of polarization across their evolution - presenting a feasible means of probing these bursts' energetic and angular properties. Although obtaining polarization data is difficult due to the inherent complexities regarding GRB observations, advances are being made, and theoretical modeling of synchrotron polarization is now more relevant than ever. In this manuscript, we present the polarization for a fiduciary model where the synchrotron forward-shock emission evolving in the radiative-adiabatic regime is described by a radially stratified off-axis outflow. This is parameterized with a power-law velocity distribution and decelerated in a constant-density and wind-like external environment. We apply this theoretical polarization model for selected bursts presenting evidence of off-axis afterglow emission, including the nearest orphan GRB candidates observed by the Neil Gehrels Swift Observatory and a few Gravitational Wave (GWs) events that could generate electromagnetic emission. In the case of GRB 170817A, we require the available polarimetric upper limits in radio wavelengths to constrain its magnetic field geometry., Comment: In submission. 18 pages, 7 figures, 3 tables

As we further our studies on Gamma-ray bursts (GRBs), both on theoretical models and observational tools, more and more options begin to open for exploration of its physical properties. As transient events primarily dominated by synchrotron radiation, it is expected that the synchrotron photons emitted by GRBs should present some degree of polarization throughout the evolution of the burst. Whereas observing this polarization can still be challenging due to the constraints on observational tools, especially for short GRBs, it is paramount that the groundwork is laid for the day we have abundant data. In this work, we present a polarization model linked with an off-axis spreading top-hat jet synchrotron scenario in a stratified environment with a density profile $n(r)\propto r^ {-k}$. We present this model's expected temporal polarization evolution for a realistic set of afterglow parameters constrained within the values observed in the GRB literature for four degrees of stratification $k=0,1,1.5 {\rm \, and\,} 2$ and two magnetic field configurations with high extreme anisotropy. We apply this model and predict polarization from a set of GRBs exhibiting off-axis afterglow emission. In particular, for GRB 170817A, we use the available polarimetric upper limits to rule out the possibility of a extremely anisotropic configuration for the magnetic field., Comment: 17 pages, 3 tables, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

The GRB 190829A has been widely studied due to its nature and the high energy emission presented. Due to the detection of a very-high-energy component by the High Energy Stereoscopic System and the event's atypically middling luminosity, it has been categorized in a select, limited group of bursts bordering classic GRBs and nearby sub-energetic events. Given the range of models utilized to adequately characterize the afterglow of this burst, it has proven challenging to identify the most probable explanation. Nevertheless, the detection of polarization data provided by the MASTER collaboration has added a new aspect to GRB 190829A that permits us to attempt to explore this degeneracy. In this paper, we present a polarization model coupled with a synchrotron forward-shock model -- a component in all models used to describe GRB 190829A's afterglow -- in order to fit the polarization's temporal evolution with the existing upper limits ($\Pi < 6\%$). We find that the polarization generated from an on-axis emission is favored for strongly anisotropic magnetic field ratios, while an off-axis scenario cannot be fully ruled out when a more isotropic framework is taken into account., Comment: 16 pages, 7 figures, 1 table

Gamma-ray bursts (GRBs) are extremely high-energy events that can be observed at very high redshift. In addition to gamma rays, they can emit in X-ray, optical, and sometimes radio wavelengths. Here, following the approach in Srinivasaragavan et al. (2020); Dainotti et al. (2021b,c), and Dainotti et al (2022, submitted), we consider 82 GRBs from Dainotti et al. (2022a) that have been observed in optical wavelengths and fitted with a broken power law (BPL). We consider the relations between the spectral and temporal indices (closure relations; CRs) according to the synchrotron forward-shock model evolving in the constant-density interstellar medium (ISM; k = 0) and the stellar Wind environment (k = 2) in both slow- and fast-cooling regimes, where the density profile is defined as n is proportional to r to the power minus k. We find the nu > max nuc, nu m regime is most favored, where nu c and nu m are the cooling and characteristic frequencies, respectively. Finally, we test the 2D Dainotti correlation between the rest-frame end time of the plateau and the luminosity at that time on GRBs that fulfill the most-favored CRs. When we compare the intrinsic scatter sigma int of those 2D correlations to the scatter presented in Dainotti et al. (2020b, 2022a), we see the scatters of our correlations generally agree with the previous values within 1sigma, both before and after correction for selection bias. This new information has helped us to pinpoint subsamples of GRBs with features that could drive the GRB emission mechanism, and eventually allow for GRBs to be used as standard candles, Comment: 18 pages , 6 figures, Accepted in The Astrophysical Journal

Cosmological models and their corresponding parameters are widely debated because of the current discrepancy between the results of the Hubble constant, $H_{0}$, obtained by SNe Ia, and the Planck data from the Cosmic Microwave Background Radiation. Thus, considering high redshift probes like Gamma-Ray Bursts (GRBs) is a necessary step. However, using GRB correlations between their physical features to infer cosmological parameters is difficult because GRB luminosities span several orders of magnitude. In our work, we use a 3-dimensional relation between the peak prompt luminosity, the rest-frame time at the end of the X-ray plateau, and its corresponding luminosity in X-rays: the so-called 3D Dainotti fundamental plane relation. We correct this relation by considering the selection and evolutionary effects with a reliable statistical method, obtaining a lower central value for the intrinsic scatter, $\sigma_{int}=0.18 \pm 0.07$ (47.1 \%) compared to previous results, when we adopt a particular set of GRBs with well-defined morphological features, called the platinum sample. We have used the GRB fundamental plane relation alone with both Gaussian and uniform priors on cosmological parameters and in combination with SNe Ia and BAO measurements to infer cosmological parameters like $H_{0}$, the matter density in the universe ($\Omega_{M}$), and the dark energy parameter $w$ for a $w$CDM model. Our results are consistent with the parameters given by the $\Lambda$CDM model but with the advantage of using cosmological probes detected up to $z=5$, much larger than the one observed for the furthest SNe Ia., Comment: 43 pages, 29 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Gamma-Ray Bursts (GRBs) are very energetic cosmological transients. Long GRBs are usually associated with Type Ib/c Supernovae (SNe), and we refer to them as GRB-SNe. Since the associated SN for a given GRB is observed only at low redshift, a possible selection effect exists when we consider intrinsically faint sources which cannot be observed at high redshift. Thus, it is important to explore the possible relationships between GRB and SN parameters after these have been corrected for astrophysical biases due to the instrumental selection effects and redshift evolution of the variables involved. So far, only GRB prompt emission properties have been checked against the SNe Ib/c properties without considering the afterglow (AG). This work investigates the existence of relationships among GRB's prompt and AG and associated SN properties. We investigate 91 bidimensional correlations among the SN and GRB observables before and after their correction for selection biases and evolutionary effects. As a result of this investigation, we find hints of a new correlation with a Pearson correlation coefficient > 0.50 and a probability of being drawn by chance < 0.05. This correlation is between the luminosity at the end of the GRB optical plateau emission and the rest-frame peak time of the SN. According to this relation, the brightest optical plateaus are accompanied by the largest peak times. This correlation is corrected for selection biases and redshift evolution and may provide new constraints for the astrophysical models associated with the GRB-SNe connection., Comment: 40 pages, 6 figures, 7 tables; accepted for publication in ApJ

Gamma-ray bursts (GRBs) are known to have the most relativistic jets, with initial Lorentz factors in the order of a few hundreds. Many GRBs display an early X-ray light-curve plateau, which was not theoretically expected and therefore puzzled the community for many years. Here, we show that this observed signal is naturally obtained within the classical GRB "fireball" model, provided that the initial Lorentz factor is rather a few tens, and the expansion occurs into a medium-low density "wind". The range of Lorentz factors in GRB jets is thus much wider than previously thought and bridges an observational gap between mildly relativistic jets inferred in active galactic nuclei, to highly relativistic jets deduced in few extreme GRBs. Furthermore, long GRB progenitors are either not Wolf-Rayet stars, or the wind properties during the final stellar evolution phase are different than at earlier times. We discuss several testable predictions of this model., Comment: 61 pages, 24 figures, 9 tables. A final edited version will appear in Nature Communications

The Fermi Large Area Telescope (Fermi-LAT) Collaboration reported the Second Gamma-ray Burst Catalog (2FLGC), which comprises a subset of 29 bursts with photon energies above 10 GeV. Although the standard synchrotron forward-shock model has successfully explained the Gamma-ray burst (GRB) afterglow observations, energetic photons higher than 10 GeV from these transient events can hardly be described in this scenario. We present the closure relations (CRs) of synchrotron self-Compton (SSC) afterglow model in the adiabatic and radiative scenario and when the central engine injects continuous energy into the blastwave to study the evolution of the spectral and temporal indexes of those bursts reported in 2FLGC. We consider the SSC afterglow model evolving in stellar-wind and interstellar medium, and the CRs as a function of the radiative parameter, the energy injection index, and the electron spectral index for $1

Sub-relativistic materials launched during the merger of binary compact objects and the core-collapse of massive stars acquire velocity structures when expanding in a stratified environment. The remnant (either a spinning magnetized neutron star (NS) or a central black hole) from the compact-object or core-collapse could additionally inject energy into the afterglow via spin-down luminosity or/and by accreting fall-back material, producing a refreshed shock, modifying the dynamics, and leading to rich radiation signatures at distinct timescales and energy bands with contrasting intensities. We derive the synchrotron light curves evolving in a stratified environment when a power-law velocity distribution parametrizes the energy of the shock, and the remnant continuously injects energy into the blastwave. As the most relevant case, we describe the latest multi-wavelength afterglow observations ($\gtrsim 900$ days) of the GW170817/GRB 170817A event via a synchrotron afterglow model with energy injection of a sub-relativistic material. The features of the remnant and the synchrotron emission of the sub-relativistic material are consistent with a spinning magnetized NS and the faster "blue" kilonova afterglow, respectively. Using the multi-band observations of some short-bursts with evidence of kilonova, we provide constraints on the expected afterglow emission., Comment: 36 pages, 14 figures, accepted for publication in ApJ

The Neil Gehrels Swift observatory observe Gamma-Ray bursts (GRBs) plateaus in X-rays. We test the reliability of the closure relations through the fireball model when dealing with the GRB plateau emission. We analyze 455 X-ray lightcurves (LCs) collected by \emph{Swift} from 2005 (January) until 2019 (August) for which the redshift is both known and unknown using the phenomenological Willingale 2007 model. Using these fits, we analyze the emission mechanisms and astrophysical environments of these GRBs through the closure relations within the time interval of the plateau emission. Finally, we test the 3D fundamental plane relation (Dainotti relation) which connects the prompt peak luminosity, the time at the end of the plateau (rest-frame), and the luminosity at that time, on the GRBs with redshift, concerning groups determined by the closure relations. This allows us to check if the intrinsic scatter \sigma_{int} of any of these groups is reduced compared to previous literature. The most fulfilled environments for the electron spectral distribution, p>2, are Wind Slow Cooling (SC) and ISM Slow Cooling for cases in which the parameter q, which indicates the flatness of the plateau emission and accounts for the energy injection, is =0 and =0.5, respectively, both in the cases with known and unknown redshifts. We also find that for the sGRBs All ISM Environments with $q=0$ have the smallest \sigma_{int}=0.04 \pm 0.15 in terms of the fundamental plane relation holding a probability of occurring by chance of p=0.005. We have shown that the majority of GRBs presenting the plateau emission fulfil the closure relations, including the energy injection, with a particular preference for the Wind SC environment. The subsample of GRBs that fulfil given relations can be used as possible standard candles and can suggest a way to reduce the intrinsic scatter of these studied relationships., Comment: 44 pages, 23 figures; Accepted to the PASJ to be published soon

Gamma-Ray Bursts (GRBs) are fascinating events due to their panchromatic nature. Their afterglow emission is observed from sub-TeV energies to radio wavelengths. We investigate GRBs that present an optical plateau, leveraging on the resemblance with the X-ray plateau shown in many GRB light curves (LCs). We comprehensively analyze all published GRBs with known redshifts and optical plateau observed mostly by the Neil Gehrels Swift Observatory (Swift). We fit 267 optical LCs and show the existence of the plateau in 102 cases, which is the largest compilation so far of optical plateaus. For 56 Swift GRBs with optical and X-ray plateaus, we compare the rest-frame end time at both wavelengths (T*_opt , T*_X), and conclude that the plateau is achromatic between T*_opt and T*_X. We also confirm the existence of the two-dimensional relations between T*_opt and the optical luminosity at the end of the plateau emission, which resembles the same luminosity-time correlation in X-rays (Dainotti et al. 2013). The existence of this optical correlation has been demonstrated for the largest sample of optical plateaus in the literature to date. The squared scatter in this optical correlation is smallest for the subset of the Gold GRBs with a decrease in the scatter equivalent to 52.4% when compared to the scatter of the entire GRB sample., Comment: 10 pages, 3 figures, 2 tables

The Fermi-LAT collaboration presented the second gamma-ray burst (GRB) catalog covering its first 10 years of operations. A significant fraction of afterglow-phase light curves in this catalog cannot be explained by the closure relations of the standard synchrotron forward-shock model, suggesting that there could be an important contribution from another process. In view of the above, we derive the synchrotron self-Compton (SSC) light curves from the reverse shock in the thick- and thin-shell regime for a uniform-density medium. We show that this emission could explain the GeV flares exhibited in some LAT light curves. Additionally, we demonstrate that the passage of the forward shock synchrotron cooling break through the LAT band from jets expanding in a uniform-density environment may be responsible for the late time ($\approx10^2$ s) steepening of LAT GRB afterglow light curves. As a particular case, we model the LAT light curve of GRB 160509A that exhibited a GeV flare together with a break in the long-lasting emission, and also two very high energy photons with energies of 51.9 and 41.5 GeV observed 76.5 and 242 s after the onset of the burst, respectively. Constraining the microphysical parameters and the circumburst density from the afterglow observations, we show that the GeV flare is consistent with a SSC reverse-shock model, the break in the long-lasting emission with the passage of the synchrotron cooling break through the Fermi-LAT band and the very energetic photons with SSC emission from the forward shock when the outflow carries a significant magnetic field ($R_{\rm B} \simeq 30$) and it decelerates in a uniform-density medium with a very low density ($n=4.554^{+1.128}_{-1.121}\times 10^{-4}\,{\rm cm^{-3}}$)., Comment: 19 pages, 5 figures. Accepted for publication in ApJ. New affiliations to the author list are added

We present the afterglow light curves produced by the deceleration of the non-relativistic ejecta mass in a stratified circumstellar medium with a density profile $n(r)\propto r^{-k}$ with $k=0$, $1$, $1.5$, $2$ and $2.5$. Once the ejecta mass is launched with equivalent kinetic energy parametrized by $E(>\beta)\propto \beta^{-\alpha}$ (where beta is the ejecta velocity) and propagates into the surrounding circumstellar medium, it first moves with constant velocity (the free-coasting phase), and later it decelerates (the Sedov-Taylor expansion). We present the predicted synchrotron and synchrotron-self Compton light curves during the free-coasting phase, and the subsequent Sedov-Taylor expansion. In particular cases, we show the corresponding light curves generated by the deceleration of several ejecta masses with different velocities launched during the coalescence of binary compact objects and the core-collapse of dying massive stars which will contribute at distinct timescales, frequencies, and intensities. Finally, using the multi-wavelength observations and upper limits collected by a large campaign of orbiting satellites and ground telescopes, we constrain the parameter space of both the KN afterglow in GW170817 and the possibly generated KN afterglow in S190814bv. Further observations on timescales of years post-merger are needed to derive tighter constraints., Comment: 30 pages, 11 figures. Accepted for publication in ApJ. The synchrotron break frequencies in the self-absorption regime were added

The X-ray afterglow plateau emission observed in many Gamma-ray Bursts (GRBs) has been interpreted as either being fueled by fallback onto a newly formed black hole, or by the spin-down luminosity of an ultra-magnetized millisecond neutron star. If the latter model is assumed, GRB X-ray afterglow light curves can be analytically reproduced. We fit a sample of GRB X-ray plateaus, interestingly yielding a distribution in the magnetic field versus spin period (B-P) diagram consistent with $B\propto P^{7/6}$. This is expected from the well-established physics of the spin-up line minimum period for Galactic millisecond pulsars. The normalisation of the relation we obtain perfectly matches spin-up line predictions for the expected masses ($\sim 1 M_{\odot}$) and radii ($\sim 10 {\rm ~km}$) of newly born magnetars, and mass accretion rates consistent with GRB expectations of $10^{-4} M_{\odot}/{\rm s} <\dot{M}< 10^{-1} M_{\odot}/{\rm s}$. Short GRBs with extended emission (SEE) appear towards the high period end of the distribution, while the long GRBs (LGRBs) towards the short period end. This result is consistent with spin-up limit expectations where the total accreted mass determines the position of the neutron star in the B-P diagram. The P-B distribution for LGRBs and SEE are statistically different, further supporting the idea that the fundamental plane relation \citep{dainotti16c,Dainotti2017} is a powerful discriminant among those populations. Our conclusions are robust against suppositions regarding the GRB collimation angle and magnetar breaking index, which shifts the resulting magnetar properties parallel to the spin-up line, and strongly support a magnetar origin for GRBs presenting X-ray plateaus., Comment: 14 pages, 7 figures, accepted to ApJ (v2: minor revisions, 4 figures added, fixed typos)

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., Comment: 25 pages, 13 figures, accepted to Advances in Space Research with minor revisions. Details on the THESEUS instrumentation, science case and expected performances can also be found in Amati et al. 2017 (arXiv:1710.04638) and in the presentations of the THESEUS Workshop 2017 (http://www.isdc.unige.ch/theseus/workshop2017-programme.html). v2 few typos corrected

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

Matter mixing is one important topic in the study of core-collapse supernova (CCSN) explosions. In this paper, we perform two-dimensional hydrodynamic simulations to reproduce the high velocity $^{56}$Ni clumps observed in SN 1987A. This is the first time that large density perturbation is proposed in the CCSN progenitor to generate Rayleigh-Taylor (RT) instability and make the effective matter mixing. In the case of a spherical explosion, RT instability is efficient at both C+O/He and He/H interfaces of the SN progenitor. Radial coherent structures shown in perturbation patterns are important for obtaining high velocity $^{56}$Ni clumps. We can also obtain matter mixing features and high velocity $^{56}$Ni clumps in some cases of aspherical explosion. We find that one of the most favorable models in our work has a combination of bipolar and equatorially asymmetric explosions in which at least 25\% of density perturbation is introduced at different composition interfaces of the CCSN progenitor. These simulation results are comparable to the observational findings of SN 1987A., Comment: ApJ accepted

We present an analysis of 123 Gamma-ray bursts (GRBs) with known redshifts possessing an afterglow plateau phase. We reveal that $L_a-T^{*}_a$ correlation between the X-ray luminosity $L_a$ at the end of the plateau phase and the plateau duration, $T^*_a$, in the GRB rest frame has a power law slope different, within more than 2 $\sigma$, from the slope of the prompt $L_{f}-T^{*}_{f}$ correlation between the isotropic pulse peak luminosity, $L_{f}$, and the pulse duration, $T^{*}_{f}$, from the time since the GRB ejection. Analogously, we show differences between the prompt and plateau phases in the energy-duration distributions with the afterglow emitted energy being on average $10\%$ of the prompt emission. Moreover, the distribution of prompt pulse versus afterglow spectral indexes do not show any correlation. In the further analysis we demonstrate that the $L_{peak}-L_a$ distribution, where $L_{peak}$ is the peak luminosity from the start of the burst, is characterized with a considerably higher Spearman correlation coefficient, $\rho=0.79$, than the one involving the averaged prompt luminosity, $L_{prompt}-L_a$, for the same GRB sample, yielding $\rho=0.60$. Since some of this correlation could result from the redshift dependences of the luminosities, namely from their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic. We apply a partial correlation coefficient to the new de-evolved luminosities showing that the intrinsic correlation exists., Comment: 12 pages, 17 pictures accepted in MNRAS, 28 May 2015

Relativistic jets propagating through an ambient medium must produce some observational effects along their side boundaries because of interactions across the large velocity gradient. One possible effect of such an interaction would be a sheared magnetic field structure at the jet boundaries, leading to a characteristic radio polarization pattern. As proposed by Ostrowski, another effect can come from the generation of a high energy cosmic ray component at the boundary, producing dynamic effects on the medium surrounding the jet and forming a cocoon dominated by cosmic rays with a decreased thermal gas emissivity. We selected this process for our first attempt to look for the effects of this type of interaction. We analyzed the Chandra X-ray data for the radio galaxy M87 in order to verify if the expected regions of diminished emissivity may be present near the spectacular X-ray jet in this source. The detailed analysis of the data, merged from 42 separate observations, shows signatures of lower emissivity surrounding the jet. In particular we detect an intensity dip along the part of the jet, which would be approximately 150 pc x 2 kpc in size, if situated along the jet which is inclined toward us. Due to a highly non-uniform X-ray background in the central region we are not able to claim the discovery of a cosmic ray cocoon around the M87 jet: we only have demonstrated that the data show morphological structures which could be accounted for if a cosmic ray cocoon exists., Comment: 8 pages, 8 pictures accepted for publication in MNRAS

Thanks to their enormous energy release, Gamma Rays Bursts (GRBs) have recently attracted a lot of interest to probe the Hubble diagram (HD) deep into the matter dominated era and hence complement Type Ia Supernovae (SNeIa). We consider here three different calibration methods based on the use of a fiducial LCDM model, on cosmographic parameters and on the local regression on SNeIa to calibrate the scaling relations proposed as an equivalent to the Phillips law to standardize GRBs finding any significant dependence. We then investigate the evolution of these parameters with the redshift to obtain any statistical improvement. Under this assumption, we then consider possible systematics effects on the HDs introduced by the calibration method, the averaging procedure and the homogeneity of the sample arguing against any significant bias., Comment: 5 pages, 1 figure to appear in the Proceedings of Frascati Workshop 2011 (Vulcano, Italy, 23th - 28th May)

Gamma Ray Bursts (GRBs) are the most energetic sources in the universe and among the farthest known astrophysical sources. These features make them appealing candidates as standard candles for cosmological applications so that studying the physical mechanisms for the origin of the emission and correlations among their observable properties is an interesting task. We consider here the luminosity L*X - break time Ta* (hereafter LT) correlation and investigate whether there are systematics induced by selection effects or redshift dependent calibra- tion. We perform this analysis both for the full sample of 77 GRBs with known redshift and for the subsample of GRBs having canonical X-ray light curves, hereafter called U0095 sample. We do not find any systematic bias thus con- firming the existence of physical GRB subclasses revealed by tight correlations of their afterglow properties. Furthermore, we study the possibility of applying the LT correlation as a redshift estimator both for the full distribution and for the canonical lightcurves. The large uncertainties and the non negligible intrin- sic scatter make the results not so encouraging, but there are nevertheless some hints motivating a further analysis with an increased U0095 sample., Comment: 23 pages, 9 figures, to appear in ApJ

Gamma Ray Bursts (GRB) observed up to redshifts $z>8$ are fascinating objects to study due to their still unexplained relativistic outburst mechanisms and a possible use to test cosmological models. Our analysis of 77 GRB afterglows with known redshifts revealed a physical subsample of long GRBs with canonical {\it plateau breaking to power-law} light curves with a significant {\it luminosity $L^*_X$ - break time $T^*_a$} correlation in the GRB rest frame. This subsample forms approximately the {\it upper envelope} of the studied distribution. We have also found a similar relation for a small sample of GRB afterglows that belong to the intermediate class (IC) between the short and the long ones. It proves that within the full sample of afterglows there exist physical subclasses revealed here by tight correlations of their afterglow properties. The afterglows with regular (`canonical') light curves obey not only a mentioned tight physical scaling, but -- for a given $T^*_a$ -- the more regular progenitor explosions lead to preferentially brighter afterglows., Comment: 15 pages, 5 figures accepted to ApJL

We present the Hubble diagram (HD) of 66 Gamma Ray Bursts (GRBs) derived using only data from their X - ray afterglow lightcurve. To this end, we use the recently updated L_X - T_a correlation between the break time T_a and the X - ray luminosity L_X measured at T_a calibrated from a sample of Swift GRBs with lightcurves well fitted by the Willingale et al. (2007) model. We then investigate the use of this HD to constrain cosmological parameters when used alone or in combination with other data showing that the use of GRBs leads to constraints in agreement with previous results in literature. We finally argue that a larger sample of high luminosity GRBs can provide a valuable information in the search for the correct cosmological model., Comment: 6 pages, 2 figures, 2 tables, submitted

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way to extend the Hubble diagram to very high redshift. To this aim, the luminosity (or isotropic emitted energy) of a GRB at redshift z must be evaluated from a correlation with a distance independent quantity so that one can then solve for the luminosity distance D_L(z) and hence the distance modulus mu(z). Averaging over five different two parameters correlations and using a fiducial cosmological model to calibrate them, Schaefer (2007) has compiled a sample of 69 GRBs with measured mu(z) which has since then been widely used to constrain cosmological parameters. We update here that sample by many aspects. First, we add a recently found correlation for the X - ray afterglow and use a Bayesian inspired fitting method to calibrate the different GRBs correlations known insofar assuming a fiducial LCDM model in agreement with the recent WMAP5 data. Averaging over six correlations, we end with a new GRBs Hubble diagram comprising 83 objects. We also extensively explore the impact of varying the fiducial cosmological model considering how the estimated mu(z) change as a function of the $(\Omega_M, w_0, w_a)$ parameters of the Chevallier - Polarski - Linder phenomenological dark energy equation of state. In order to avoid the need of assuming an {\it a priori} cosmological model, we present a new calibration procedure based on a model independent local regression estimate of mu(z) using the Union SNeIa sample to calibrate the GRBs correlations. This finally gives us a GRBs Hubble diagram made out of 69 GRBs whose estimated distance modulus mu(z) is almost independent on the underlying cosmological model., Comment: version accepted for publication on MNRAS; 19 pages, 8 figures, 4 tables; contains two new sections on the analysis of systematic errors and cosmological constraints using present day and forecasts on future data

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way to extend the Hubble diagram to very high redshift. However, the large scatter in their intrinsic properties prevents directly using them as distance indicator so that the hunt is open for a relation involving an observable property to standardize GRBs in the same way as the Phillips law makes it possible to use Type Ia Supernovae (SNeIa) as standardizable candles. We use here the data on the X - ray decay curve and spectral index of a sample of GRBs observed with the Swift satellite. These data are used as input to a Bayesian statistical analysis looking for a correlation between the X - ray luminosity L_X(T_a) and the time constant T_a of the afterglow curve. We find a linear relation between \log{[L_X(T_a)]} and \log{[T_a/(1+z)]} with an intrinsic scatter sigma_{int} = 0.33 comparable to previously reported relations. Remarkably, both the slope and the intrinsic scatter are almost independent on the matter density Omega_M and the constant equation of state w of the dark energy component thus suggesting that the circularity problem is alleviated for the $L_X - T_a$ relation., Comment: 5 pages, 2 figures, accepted for publication on MNRAS Letters