Search

Full list of authors: Makhathini, S.; Mooley, K. P.; Brightman, M.; Hotokezaka, K.; Nayana, A. J.; Intema, H. T.; Dobie, D.; Lenc, E.; Perley, D. A.; Fremling, C.; Moldòn, J.; Lazzati, D.; Kaplan, D. L.; Balasubramanian, A.; Brown, I. S.; Carbone, D.; Chandra, P.; Corsi, A.; Camilo, F.; Deller, A.; Frail, D. A.; Murphy, T.; Murphy, E. J.; Nakar, E.; Smirnov, O.; Beswick, R. J.; Fender, R.; Hallinan, G.; Heywood, I.; Kasliwal, M.; Lee, B.; Lu, W.; Rana, J.; Perkins, S.; White, S. V.; Józsa, G. I. G.; Hugo, B. ; Kamphuis, P., We present the full panchromatic afterglow light-curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups and attempted to mitigate these differences to provide a more uniform data set. Simple power-law fits to the uniform afterglow light curve indicate a t 0.86±0.04 rise, a t -1.92±0.12 decline, and a peak occurring at 155 ± 4 days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index (-0.584 ± 0.002) across all epochs. This gives a precise and updated estimate of the electron power-law index, p = 2.168 ± 0.004. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized interstellar medium density, The MeerKAT telescope is operated by the South African Radio Astronomy Observatory, (SARAO), which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Technology. K.P.M. is a Jansky Fellow of the National Radio Astronomy Observatory. K.P.M. and G.H. acknowledge support from the National Science Foundation Grant AST-1911199. D.D. is supported by an Australian Government Research Training Program Scholarship. T.M. acknowledges the support of the Australian Research Council through grant DP190100561. Parts of this research were conducted by the Australian Research Council Centre of Excellence for Gravitational Wave Discovery (OzGrav), project number CE170100004. We acknowledge support by the GROWTH (Global Relay of Observatories Watching Transients Happen) project funded by the National Science Foundation PIRE (Partnership in International Research and Education) program under grant No. 1545949. D.L. acknowledges support from NASA grants 80NSSC18K1729 (Fermi) and NNX17AK42G (ATP), Chandra grant TM9-20002X, and NSF grant AST-1907955. This research has made use of NASA's Astrophysics Data System Bibliographic Services. C.F. gratefully acknowledges support of his research by the Heising-Simons Foundation. JM acknowledges financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709) and from the grant RTI2018-096228-B-C31 (MICIU/FEDER, EU). D.L.K. was supported by NSF grant AST-1816492. P.K. is partially supported by the BMBF project 05A17PC2for D-MeerKAT. A.B. and A.C. acknowledge support from the National Science Foundation via grant #1907975.

GRB 190114C was a bright burst that occurred in the local Universe (z=0.425). It was the first gamma-ray burst (GRB) ever detected at TeV energies, thanks to MAGIC. We characterize the ambient medium properties of the host galaxy through the study of the absorbing X-ray column density. Joining Swift, XMM-Newton, and NuSTAR observations, we find that the GRB X-ray spectrum is characterized by a high column density that is well in excess of the expected Milky Way value and decreases, by a factor of ~2, around ~$10^5$ s. Such a variability is not common in GRBs. The most straightforward interpretation of the variability in terms of photoionization of the ambient medium is not able to account for the decrease at such late times, when the source flux is less intense. Instead, we interpret the decrease as due to a clumped absorber, denser along the line of sight and surrounded by lower-density gas. After the detection at TeV energies of GRB 190114C, two other GRBs were promptly detected. They share a high value of the intrinsic column density and there are hints for a decrease of the column density, too. We speculate that a high local column density might be a common ingredient for TeV-detected GRBs., 6 pages, 4 figures. Accepted for publication in A&A

We present new radio observations of the binary neutron star merger GW170817 carried out with the Karl G. Jansky Very large Array (VLA) more than 3\,yrs after the merger. Our combined dataset is derived by co-adding more than $\approx32$\,hours of VLA time on-source, and as such provides the deepest combined observation (rms sensitivity $\approx 0.99\,\mu$Jy) of the GW170817 field obtained to date at 3\,GHz. We find no evidence for a late-time radio re-brightening at a mean epoch of $t\approx 1200$\,d since merger, in contrast to a $\approx 2.1\,\sigma$ excess observed at X-ray wavelengths at the same mean epoch. Our measurements agree with expectations from the post-peak decay of the radio afterglow of the GW170817 structured jet. Using these results, we constrain the parameter space of models that predict a late-time radio re-brightening possibly arising from the high-velocity tail of the GW170817 kilonova ejecta, which would dominate the radio and X-ray emission years after the merger (once the structured jet afterglow fades below detection level). Our results point to a steep energy-speed distribution of the kilonova ejecta (with energy-velocity power law index $\alpha \gtrsim 5$). We suggest possible implications of our radio analysis, when combined with the recent tentative evidence for a late-time re-brightening in the X-rays, and highlight the need for continued radio-to-X-ray monitoring to test different scenarios., Comment: 6 pages, 2 figures. Version published in ApJL (includes minor updates to figures, data and refs)

A comprehensive understanding of Gamma Ray Bursts (GRBs) has been elusive due to the variety of questions surrounding the radiation mechanism at play in these events. Polarization measurements of GRBs can heavily constrain the relevant radiation mechanisms and the structure of the GRB jet; however, there is a limited number of theoretical predictions that observed GRB polarizations can be compared against. Here, we conduct radiative transfer calculations of a set of two dimensional relativistic hydrodynamic long GRB (LGRB) jet simulations, of a constant and a variable jet, using the Monte Carlo Radiation Transport (MCRaT) code. MCRaT has been enhanced by the inclusion of polarization; it has been first verified by reproducing a variety of results in the literature and then used to obtain the time integrated and resolved polarization degrees and angles of the synthetic LGRBs. While the obtained time-integrated polarization degrees {($\lesssim 1$\%)} are consistent with the constraints from the POLAR experiment, they are lower than other theoretical studies due to the lack of strong gradients in the model jet profiles that we use. The time resolved results suggests that GRBs with wide jets observed on axis will have small polarization degrees ($\lesssim 2\%$) and {constant polarization angles}, during the brightest portion of the light curve. GRBs observed off axis will have larger polarization degrees and polarization angles that change with the temporal structure of radiating shells in the outflow. We then place our results in the context of GRB prompt emission models and future LEAP and POLAR-2 GRB polarimetry detections., 14 pages, 7 figures, final version after it has been accepted to ApJ

The association of GRB170817A with GW170817 has confirmed the long-standing hypothesis that binary neutron star (BNS) mergers are the progenitors of at least some short gamma-ray bursts (SGRBs). This connection has ushered in an era in which broadband observations of SGRBs, together with measurements of the time delay between the gravitational waves and the electromagnetic radiation, allow to probe the properties of the emitting outflow and its engine to an unprecedented detail. Since the structure of the radiating outflow is molded by the interaction of a relativistic jet with the binary ejecta, it is of paramount importance to study the system in a realistic setting. Here we present a three-dimensional hydrodynamic simulation of a relativistic jet propagating in the ejecta of a BNS merger, which were computed with a general relativistic magnetohydrodynamic simulation. We find that the jet's centroid oscillates around the axis of the system, due to inhomogeneities encountered in the propagation. These oscillations allow the jet to find the path of least resistance and travel faster than an identical jet in smooth ejecta. In our setup the breakout time is ~0.6 sec, comparable to the expected central engine duration in SGRBs and possibly a non-negligible fraction of the total delay between the gravitational and gamma-ray signals. Our simulation also shows that energy is carried in roughly equal amounts by the jet and by the cocoon, and that about 20 per cent of the injected energy is transferred to the ejecta via mechanical work., Under review for publication in ApJL

It has been proposed that blazar jets are structured, with a fast spine surrounded by a slower sheath or layer. This structured jet model explains some properties of their emission and morphology. Because of their relative motion, the radiation produced by one component is seen amplified by the other, thus enhancing the inverse Compton emission of both. Radiation is emitted anisotropically in the comoving frames, and causes the emitting plasma to recoil. As seen in the observer frame, this corresponds to a deceleration of the fastest component (the spine) and an acceleration of the slower one (the layer). While the deceleration of the spine has already been investigated, here we study for the first time the acceleration of the sheath and find self-consistent velocity profile solutions for both the spine and the sheath while accounting for radiative cooling. We find that the sheath can be accelerated to the velocities required by the observations if its leptons remain energetic in the acceleration region, assumed to be of the order of 100 Schwarzschild radii, demanding continuous injection of energetic particles in that region., Comment: 14 pages, 16 figures, submitted to MNRAS. Comments are welcome!

Classical novae commonly show evidence of rapid dust formation within months of the outburst. However, it is unclear how molecules and grains are able to condense within the ejecta, given the potentially harsh environment created by ionizing radiation from the white dwarf. Motivated by the evidence for powerful radiative shocks within nova outflows, we propose that dust formation occurs within the cool, dense shell behind these shocks. We incorporate a simple molecular chemistry network and classical nucleation theory with a model for the thermodynamic evolution of the post-shock gas, in order to demonstrate the formation of both carbon and forsterite ($\rm Mg_2SiO_4$) grains. The high densities due to radiative shock compression ($n \sim 10^{14}$ cm$^{-3}$) result in CO saturation and rapid dust nucleation. Grains grow efficiently to large sizes $\gtrsim 0.1\mu$m, in agreement with IR observations of dust-producing novae, and with total dust masses sufficient to explain massive extinction events such as V705 Cas. As in dense stellar winds, dust formation is CO-regulated, with carbon-rich flows producing carbon-rich grains and oxygen-rich flows primarily forming silicates. CO is destroyed by non-thermal particles accelerated at the shock, allowing additional grain formation at late times, but the efficiency of this process appears to be low. Given observations showing that individual novae produce both carbonaceous and silicate grains, we concur with previous works attributing this bimodality to chemical heterogeneity of the ejecta. Nova outflows are diverse and inhomogeneous, and the observed variety of dust formation events can be reconciled by different abundances, the range of shock properties, and the observer viewing angle. The latter may govern the magnitude of extinction, with the deepest extinction events occurring for observers within the binary equatorial plane.

The detection of GW170817, it's extensive multi-wavelength follow-up campaign, and the large amount of theoretical development and interpretation that followed, have resulted in a significant step forward in the understanding of the binary neutron star merger phenomenon as a whole. One of its aspects is seeing the merger as a progenitor of short gamma-ray bursts (SGRB), which will be the subject of this review. On the one hand, GW170817 observations have confirmed some theoretical expectations, exemplified by the confirmation that binary neutron star mergers are the progenitors of SGRBs. In addition, the multimessenger nature of GW170817 has allowed for gathering of unprecedented data, such as the trigger time of the merger, the delay with which the gamma-ray photons were detected, and the brightening afterglow of an off-axis event. All together, the incomparable richness of the data from GW170817 has allowed us to paint a fairly detailed picture of at least one SGRB. I will detail what we learned, what new questions have arisen, and the perspectives for answering them when a sample of GW170817-comparable events have been studied., Comment: Mini-review for Frontiers. Final published version

We study the ensemble of linear polarization measurement in the optical afterglows of long-duration gamma-ray bursts. We assume a non sideways-expanding top-hat jet geometry and use the relatively large number of measurements under the assumption that they represent a statistically unbiased sample. This allows us to constrain the ratio between the maximum predicted polarization and the measured one, which is an indicator of the geometry of the magnetic field in the downstream region of the external shock. We find that the measured polarization is substantially suppressed with respect to the maximum possible for either a completely ordered magnetic field parallel to the shock normal or to a field that is entirely contained in the shock plane. The measured polarization is limited, on average, to between 25 and 30% of the maximum theoretically possible value. This reduction requires the perpendicular component of the magnetic field to be dominant in energy with respect to the component parallel to the shock front, as expected for a shock generated and/or shock compressed field. We find, however, that the data only marginally support the assumption of a simple top-hat jet, pointing towards a more complex geometry for the outflow., Comment: 6 pages, submitted to AAS journals

We review the current observational status and theoretical interpretations for the class of broad lines type Ic supernovae. They are characterized by fast photospheric expansion and lack of H and He absorption. They have a larger than normal energy budget, suggesting that they are powered or, at least, augmented by a central engine, like a magnetar or an accreting black hole. There appears therefore to be a link between these supernovae and long-duration gamma-ray bursts. However, its nature has not been satisfactorily demystified.

The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multi-wavelength electromagnetic radiation. The almost simultaneous observation of a pulse of gamma-rays proved that BNS mergers are associated with at least some short gamma-ray bursts (GRBs). However, the gamma-ray pulse was faint, casting doubts on the association of BNS mergers with the luminous, highly relativistic outflows of canonical short GRBs. Here we show that structured jets with a relativistic, energetic core surrounded by slower and less energetic wings produce afterglow emission that brightens characteristically with time, as recently seen in the afterglow of GW170817. Initially, we only see the relatively slow material moving towards us. As time passes, larger and larger sections of the outflow become visible, increasing the luminosity of the afterglow. The late appearance and increasing brightness of the multi-wavelength afterglow of GW170817 allow us to constrain the geometry of its ejecta and thus reveal the presence of an off-axis jet pointing about 30 degrees away from Earth. Our results confirm a single origin for BNS mergers and short GRBs: GW170817 produced a structured outflow with a highly relativistic core and a canonical short GRB. We did not see the bright burst because it was beamed away from Earth. However, approximately one in 20 mergers detected in gravitational waves will be accompanied by a bright, canonical short GRB., Models updated with new data and added references. Accepted for publication in PRL, 8 pages, 7 figures and 1 table. A grid of models, jet properties, and python interpolating routine is available at http://www.science.oregonstate.edu/~lazzatid/cocoon.html

The classical theory of grain nucleation suffers from both theoretical and predictive deficiencies. We strive to alleviate these deficiencies in our understanding of dust formation and growth by utilizing an atomistic model of nucleation. Carbon cluster geometries are determined with a set of global minimization algorithms. Using density functional theory, the binding energies of carbon clusters from $n=2$ to $n=99$ are then calculated. These energies are used to calculate the critical size and nucleation rate of carbon clusters. We find that the critical cluster size is largely determined by the changes in geometry of the clusters. Clusters with size $n=27$ and $n=8$, roughly corresponding to the transition from ring-to-fullerene geometry and chain-to-ring geometry respectively, are the critical sizes across the range of temperature and saturation where nucleation is significant. In contrast to the classical theory, nucleation is enhanced at low-temperatures, and suppressed at high temperatures. These results will be applied to a modified chemical evolution code using results from supernova simulations., 13 pages, 7 figures. Proceedings of Cosmic Dust VIII conference in Tokyo, 2015

The disks of active galactic nuclei (AGNs), traditionally studied as feeders of the supermassive black holes (SMBHs) at their centers, are also hosts to massive stars and hence their neutron star (NS) and black hole (BH) remnants. Migration traps and gas torques in these disks favor binary formation and enhance the rate of compact object mergers. In these environments both long gamma-ray bursts (GRBs) from the death of massive stars and short GRBs from NS–NS to NS–BH mergers are expected. However, their properties in the environment of AGN disks have never been studied. Here we show that GRBs in AGN disks can display unique features, owing to the unusual relative position of the shocks that characterize the burst evolution and the Thomson photosphere of the AGN disk. In dense environments, for example, a relativistic reverse shock develops early, likely powering the prompt emission instead of internal shocks. The transient’s time evolution is also compressed, yielding afterglow emission that is brighter and may peak earlier than for GRBs in the interstellar medium. Additionally, in regions of the disk that are sufficiently dense and extended, the light curves are dominated by diffusion, since the fireball remains inside the disk photosphere throughout the entire evolution. These diffusion-dominated transients emerge on timescales of days in disks around SMBHs of ∼ 106 M ⊙ to years for SMBHs of ∼ 108 M ⊙. Finally, a large fraction of events, especially in AGNs with SMBHs ≲ 107 M ⊙, display time-variable absorption in the X-ray band.

GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet-wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant., Comment: Accepted for publication in ApJ

The discovery of GW170817, the merger of a binary neutron star (NS) triggered by a gravitational wave detection by LIGO and Virgo, has opened a new window of exploration in the physics of NSs and their cosmological role. Among the important quantities to measure are the mass and velocity of the ejecta produced by the tidally disrupted NSs and the delay - if any - between the merger and the launching of a relativistic jet. These encode information on the equation of state of the NS, the nature of the merger remnant, and the jet launching mechanism, as well as yielding an estimate of the mass available for r-process nucleosynthesis. Here we derive analytic estimates for the structure of jets expanding in environments with different density, velocity, and radial extent. We compute the jet-cocoon structure and the properties of the broadband afterglow emission as a function of the ejecta mass, velocity, and time delay between merger and launch of the jet. We show that modeling of the afterglow light curve can constrain the ejecta properties and, in turn, the physics of neutron density matter. Our results increase the interpretative power of electromagnetic observations by allowing for a direct connection with the merger physics., Comment: ApJ in press. Updated references and discussion of BNS ejecta

Gamma-ray bursts are a complex, non-linear system that evolves very rapidly through stages of vastly different conditions. They evolve from scales of few hundred kilometers where they are very dense and hot to cold and tenuous on scales of parsecs. As such, our understanding of such a phenomenon can truly increase by combining theoretical and numerical studies adopting different numerical techniques to face different problems and deal with diverse conditions. In this review, we will describe the tremendous advancement in our comprehension of the bursts phenomenology through numerical modeling. Though we will discuss studies mainly based on jet dynamics across the progenitor star and the interstellar medium, we will also touch upon other problems such as the jet launching, its acceleration, and the radiation mechanisms. Finally, we will describe how combining numerical results with observations from Swift and other instruments resulted in true understanding of the bursts phenomenon and the challenges still lying ahead., Comment: 14 Pages, Journal of High Energy Astrophysics for the dedicated issue: "Swift: Ten Years of Discovery", in press

The mechanism behind the prompt emission of GRBs is a longstanding mystery yet to be resolved. To tackle this issue, we present Monte Carlo simulations of a scattering dominated photon-lepton plasma. Initially in equilibrium, the photons and leptons are driven out of equilibrium by energy injection which mimics energy dissipation through shocks/magnetic reconnection. We discuss non-thermal spectral evolution (in agreement with the Band function) and how spectral correlations occur among the Band parameters, which future observations can verify. We then unveil an improved Monte Carlo code that dynamically evolves and tracks the positions and four-momenta of photons and leptons (which interact via Compton scattering) in an expanding relativistic jet. We discuss preliminary results and potential applications of this novel code.

We present the results of numerical simulations of the prompt emission of short-duration gamma-ray bursts. We consider emission from the relativistic jet, the mildly relativistic cocoon, and the non-relativistic shocked ambient material. We find that the cocoon material is confined between off-axis angles 15, Comment: 6 pages, 8 figures, accepted for publication in ApJL

Short Gamma-Ray Bursts (GRBs) are explosions of cosmic origin believed to be associated with the merger of two compact objects, either two neutron stars, or a neutron star and a black hole. The presence of at least one neutron star has long been thought to be an essential element of the model: its tidal disruption provides the needed baryonic material whose rapid accretion onto the post-merger black hole powers the burst. The recent tentative detection by the Fermi satellite of a short GRB in association with the gravitational wave signal GW150914 produced by the merger of two black holes has challenged this standard paradigm. Here we show that the evolution of two high-mass, low-metallicity stars with main sequence rotational speeds a few tens of percent of the critical speed eventually undergoing a weak supernova explosion {\em can} produce a short gamma-ray burst. The outer layers of the envelope of the last exploding star remain bound and circularize at large radii. With time, the disk cools and becomes neutral, suppressing the magneto-rotational instability, and hence the viscosity. The disk remains 'long-lived dead' until tidal torques and shocks during the pre-merger phase heat it up and re-ignite accretion, rapidly consuming the disk and powering the short gamma-ray burst., Matches version accepted to ApJL

XIPE, the X-ray Imaging Polarimetry Explorer, is a mission dedicated to X-ray Astronomy. At the time of writing XIPE is in a competitive phase A as fourth medium size mission of ESA (M4). It promises to reopen the polarimetry window in high energy Astrophysics after more than 4 decades thanks to a detector that efficiently exploits the photoelectric effect and to X-ray optics with large effective area. XIPE uniqueness is time-spectrally-spatially- resolved X-ray polarimetry as a breakthrough in high energy astrophysics and fundamental physics. Indeed the payload consists of three Gas Pixel Detectors at the focus of three X-ray optics with a total effective area larger than one XMM mirror but with a low weight. The payload is compatible with the fairing of the Vega launcher. XIPE is designed as an observatory for X-ray astronomers with 75 % of the time dedicated to a Guest Observer competitive program and it is organized as a consortium across Europe with main contributions from Italy, Germany, Spain, United Kingdom, Poland, Sweden.

In this paper, we report on the available X-ray data collected by INTEGRAL, Swift, and XMM-Newton during the first outburst of the INTEGRAL transient IGR J17451-3022, discovered in 2014 August. The monitoring observations provided by the JEM-X instruments on-board INTEGRAL and the Swift/XRT showed that the event lasted for about 9 months and that the emission of the source remained soft for the entire period. The source emission is dominated by a thermal component (kT~1.2 keV), most likely produced by an accretion disk. The XMM-Newton observation carried out during the outburst revealed the presence of multiple absorption features in the soft X-ray emission that could be associated to the presence of an ionized absorber lying above the accretion disk, as observed in many high-inclination low mass X-ray binaries. The XMM-Newton data also revealed the presence of partial and rectangular X-ray eclipses (lasting about 820 s), together with dips. The latter can be associated with increases in the overall absorption column density in the direction of the source. The detection of two consecutive X-ray eclipses in the XMM-Newton data allowed us to estimate the source orbital period at 22620.5(-1.8,+2.0) s (1{\sigma} c.l.)., Comment: Accepted for publication on A&A

Unlike mergers of two compact objects containing a neutron star (NS), binary black hole (BBH) mergers are not accompanied by the production of tidally disrupted material, and hence lack the most direct source of accretion to power a jet and generate electromagnetic (EM) radiation. However, following a tentative detection by the Fermi GBM of a $\gamma$-ray counterpart to GW150914, several alternative ideas were proposed for driving a jet and producing EM radiation. {\em If} such jets were in fact produced, they would however lack the bright cocoon emission that makes jets from binary NSs bright also at large viewing angles. Here, via Monte Carlo simulations of a population of BBH mergers with properties consistent with those inferred from the existing LIGO/Virgo observations, and the angular emission characteristic of jets propagating into the interstellar medium, we derive limits on the allowed energetics and Lorentz factors of such jets from EM follow ups to GW-detected BBH merger events to date, and we make predictions which will help tighten these limits with broadband EM follow ups to events in future LIGO/Virgo runs. The condition that $\lesssim 1$ event out of 10 GW-detected BBH mergers be above the Fermi/GBM threshold imposes that any currently allowed emission model has to satisfy the condition $E_{\rm iso}/(10^{49} {\rm erg})(\theta_{jet}/20^\circ)\lesssim 1$., Comment: Matches version in press in ApJ - extended library of lightcurves publicly available here: http://www.astro.sunysb.edu/rosalba/EMmod/models.html

New information on short/hard gamma-ray bursts (GRBs) is being gathered thanks to the discovery of their optical and X-ray afterglows. However, some key aspects are still poorly understood, including the collimation level of the outflow, the duration of the central engine activity, and the properties of the progenitor systems. We want to constrain the physical properties of the short GRB 050724 and of its host galaxy, and make some inferences on the global short GRB population. We present optical observations of the afterglow of GRB 050724 and of its host galaxy, significantly expanding the existing dataset for this event. We compare our results with models, complementing them with available measurements from the literature. We study the afterglow light curve and spectrum including X-ray data. We also present observations of the host galaxy. The observed optical emission was likely related to the large flare observed in the X-ray light curve. The apparent steep decay was therefore not due to the jet effect. Available data are indeed consistent with low collimation, in turn implying a large energy release, comparable to that of long GRBs. The flare properties also constrain the internal shock mechanism, requiring a large Lorentz factor contrast between the colliding shells. This implies that the central engine was active at late times, rather than ejecting all shells simultaneously. The host galaxy has red colors and no ongoing star formation, consistent with previous findings on this GRB. However, it is not a pure elliptical, and has some faint spiral structure. GRB 050724 provides the most compelling case for association between a short burst and a galaxy with old stellar population. It thus plays a pivotal role in constraining progenitors models, which should allow for long delays between birth and explosion., 8 pages, 5 figures, 4 tables, accepted for publication in A&A, typo fixed

The duration of the prompt emission of long gamma-ray bursts is generally considered to be fairly similar to the duration of the activity of the engine in the center of the progenitor star. Here, we investigate the relation between the duration of the engine activity and that of the observed light curve, using inputs from both numerical simulations and observations. We find that the observed burst duration is a good proxy for the engine duration after the time necessary for the jet to break out the star's surface is subtracted. However, the observed duration is a function of the viewing angle and can be significantly shorter than the duration of the engine activity. We also show that the observed, redshift-corrected burst duration evolves only moderately with redshift for both observations and synthetic light curves. We conclude that the broad distribution of the observed duration of long BATSE gamma-ray bursts is mostly accounted for by an engine lasting ~20 s, the dispersion being due to viewing and redshift effects. Our results do not rule out the existence of engines with very long duration. However, we find that they are constrained to be a small minority of the BATSE detected bursts., Comment: 6 pages, 6 figures, accepted for publication in MNRAS. One typo corrected in the abstract

A peculiar infrared ring-like structure was discovered by {\em Spitzer} around the strongly magnetised neutron star SGR 1900$+$14. This infrared structure was suggested to be due to a dust-free cavity, produced by the SGR Giant Flare occurred in 1998, and kept illuminated by surrounding stars. Using a 3D dust radiative transfer code, we aimed at reproducing the emission morphology and the integrated emission flux of this structure assuming different spatial distributions and densities for the dust, and different positions for the illuminating stars. We found that a dust-free ellipsoidal cavity can reproduce the shape, flux, and spectrum of the ring-like infrared emission, provided that the illuminating stars are inside the cavity and that the interstellar medium has high gas density ($n_H\sim$1000 cm$^{-3}$). We further constrain the emitting region to have a sharp inner boundary and to be significantly extended in the radial direction, possibly even just a cavity in a smooth molecular cloud. We discuss possible scenarios for the formation of the dustless cavity and the particular geometry that allows it to be IR-bright., 13 pages, 8 figures, accepted for publication in The Astrophysical Journal

We present calculations of the wide angle emission of short-duration gamma-ray bursts from compact binary merger progenitors. Such events are expected to be localized by their gravitational wave emission, fairly irrespective of the orientation of the angular momentum vector of the system, along which the gamma-ray burst outflow is expected to propagate. We show that both the prompt and afterglow emission are dim and challenging to detect for observers lying outside of the cone within which the relativistic outflow is propagating. If the jet initially propagates through a baryon contaminated region surrounding the merger site, however, a hot cocoon forms around it. The cocoon subsequently expands quasi-isotropically producing its own prompt emission and external shock powered afterglow. We show that the cocoon prompt emission is detectable by Swift BAT and Fermi GBM, We also show that the cocoon afterglow peaks a few hours to a few days after the burst and is detectable for up to a few weeks at all wavelengths. The timing and brightness of the transient are however uncertain due to their dependence on unknown quantities such as the density of the ambient medium surrounding the merger site, the cocoon energy, and the cocoon Lorentz factor. For a significant fraction of the gravitationally-detected neutron-star-binary mergers, the cocoon afterglow could possibly be the only identifiable electromagnetic counterpart, at least at radio and X-ray frequencies., 10 pages, 6 figures, accepted for publication in MNRAS

Long Gamma-ray Bursts (GRBs) signal the death of massive stars. The afterglow emission can be used to probe the progenitor ambient through a detailed study of the absorption pattern imprinted by the circumburst material as well as the host galaxy interstellar medium on the continuum spectrum. This has been done at optical wavelengths with impressive results. Similar studies can in principle be carried out in the X-ray band, allowing us to shed light on the material metallicity, composition and distance of the absorber. We start exploiting this route through high resolution spectroscopy XMM-Newton observations of GRB 090618. We find a high metallicity absorbing medium (Z> 0.2 Zsun) with possible enhancements of S and Ne with respect to the other elements (improving the fit at a level of >3.4 sigma). Including the metallicity effects on the X-ray column density determination, the X-ray and optical evaluations of the absorption are in agreement for a Small Magellanic Cloud extinction curve.

We study the role of the unknown microphysical properties of carbonaceous dust particles in determining the amount and size distribution of carbonaceous dust condensed in type II supernova explosions. We parametrize the microphysical properties in terms of the shape factor of the grain and the sticking coefficient of gas-phase carbon atoms onto the grain surfaces. We find that the amount of dust formed is fairly independent of these properties, within the parameter range considered, though limited by the available amount of carbon atoms not locked in CO molecules. However, we find that the condensation times and size distributions of dust grains depend sensitively on the microphysical parameters, with the mass distributions being weighted toward larger effective radii for conditions considering grains with higher sticking coefficients and/or more aspherical shapes. We discuss that this leads to important consequences on the predicted extinction law of SN dust and on the survival rate of the formed grains as they pass through the reverse shock of the SN. We conclude that a more detailed understanding of the dust formation process and of the microphysical properties of each dust species needs to be achieved before robust prediction on the SN dust yields can be performed.

Long duration gamma-ray bursts are commonly associated with the deaths of massive stars. Spectroscopic studies using the afterglow as a light source provide a unique opportunity to unveil the medium surrounding it, probing the densest region of their galaxies. This material is usually in a low ionisation state and at large distances from the burst site, hence representing the normal interstellar medium in the galaxy. Here we present the case of GRB 090426 at z=2.609, whose optical spectrum indicates an almost fully ionised medium together with a low column density of neutral hydrogen. For the first time, we also observe variations in the Ly alpha absorption line. Photoionisation modeling shows that we are probing material from the vicinity of the burst (~80 pc). The host galaxy is a complex of two luminous interacting galaxies, which might suggest that this burst could have occurred in an isolated star-forming region outside its host galaxy created in the interaction of the two galaxies.

By coupling radiation transfer calculations to hydrodynamic simulations, there have been major advancements in understanding the long Gamma Ray Burst (LGRB) prompt emission. Building upon these achievements, we present an analysis of photospheric emission acquired by using the Monte Carlo Radiation Transfer (MCRaT) code on hydrodynamic simulations with variable jet profiles. MCRaT propagates and Compton scatters individual photons that have been injected into the collimated outflow in order to produce synthetic light curves and spectra. These light curves and spectra allow us to compare our results to LGRB observational data. We find excellent agreement between our fitted time resolved $\beta$ parameters and those that are observed. Additionally, our simulations show that photospheric emission, under certain conditions, is able to create the observationally expected Band $\alpha$ parameter. Finally, we show that the simulations are consistent with the Golenetskii correlation but exhibit some strain with the Amati and Yonetoku correlations., Comment: 10 pages, 8 figures, this version is the final accepted manuscript to appear in ApJ

We present a numerical simulation of a gamma-ray burst jet from a long-lasting engine in the core of a 16 solar mass Wolf-Rayet star. The engine is kept active for 6000 s with a luminosity that decays in time as a power-law with index -5/3. Even though there is no short time-scale variability in the injected engine luminosity, we find that the jet's kinetic luminosity outside the progenitor star is characterized by fluctuations with relatively short time scale. We analyze the temporal characteristics of those fluctuations and we find that they are consistent with the properties of observed flares in X-ray afterglows. The peak to continuum flux ratio of the flares in the simulation is consistent with some, but not all, the observed flares. We propose that propagation instabilities, rather than variability in the engine luminosity, are responsible for the X-ray flares with moderate contrast. Strong flares such as the one detected in GRB 050502B, instead, cannot be reproduced by this model and require strong variability in the engine activity.

We investigate the role of thermal fluctuations and of the finite number of monomers in small clusters of carbon atoms on the nucleation rate of carbonaceous grains. Thermal fluctuations are due to the quantized nature of the energy exchanges between the clusters, the gas, and the radiation field. Nanoscale effects modify the spontaneous detachment of monomers due to the finite amount of internal energy contained in small clusters. We find that both corrections have a big impact on the stability of the clusters and on the rate of nucleation. We implement our model within a Monte Carlo code to derive the new stability conditions for clusters as well as nucleation rates. Due to computing limitations, we can explore the consequences of this approach only at high temperatures, at which particle interactions are not much less frequent than photon interactions. We found that the combined effect of the detachment correction and the temperature fluctuations produces faster nucleation. We also found that the nucleation rate depends on the composition of the gas and not only on the partial pressure of the compound that condensates into grains. This is a unique result of this model that can be used to prove or disprove it.

We present late-time radio observations of GW170817, the first binary neutron star merger discovered through gravitational waves by the advanced LIGO and Virgo detectors. Our observations, carried out with the Karl G. Jansky Very Large Array, were optimized to detect polarized radio emission, and thus to constrain the linear polarization fraction of GW170817. At an epoch of ~244 days after the merger, we rule out linearly polarized emission above a fraction of ~12% at a frequency of 2.8 GHz (99% confidence). Within the structured jet scenario (a.k.a. successful jet plus cocoon system) for GW170817, the derived upper-limit on the radio continuum linear polarization fraction strongly constrains the magnetic field configuration in the shocked ejecta. We show that our results for GW170817 are compatible with the low level of linear polarization found in afterglows of cosmological long gamma-ray bursts. Finally, we discuss our findings in the context of future expectations for the study of radio counterparts of binary neutron star mergers identified by ground-based gravitational-wave detectors., 5 pages, 2 figures, 1 table

Ultra-long gamma ray bursts (ULGRBs) are a distinct class of GRBs characterized by durations of several thousands of seconds, about two orders of magnitude longer than those of standard long GRBs (LGRBs). The driving engine of these events has not been uncovered yet, and ideas range from magnetars, to tidal disruption events, to extended massive stars, such as blue super giants (BSG). BSGs, a possible endpoint of stellar evolution, are attractive for the relatively long free-fall times of their envelopes, allowing accretion to power a long-lasting central engine. At the same time, their large radial extension poses a challenge to the emergence of a jet. Here we perform an end-to-end simulation aimed at assessing the viability of BSGs as ULGRB progenitors. The evolution to core collapse of a BSG star model is calculated with the MESA code. We then compute the accretion rate for the fraction of envelope material with enough angular momentum to circularize and form an accretion disk, and input the corresponding power into a jet which we evolve through the star envelope with the FLASH code. Our simulation shows that the jet can emerge, and the resulting light curves resemble those observed in ULGRBs, with durations $T_{90}$ ranging from $\approx$4000~s to $\approx10^4$ s depending on the viewing angle., Comment: Matches version in press in ApJ

The mass of the stellar precursor of supernova (SN) 1987A and the burst of neutrinos observed at the moment of the explosion are consistent with the core-collapse formation of a neutron star. However, no compelling evidence for the presence in SN 1987A of a compact object of any kind has been found yet in any band of the electromagnetic spectrum, prompting questions on whether the neutron star survived and, if it did, on its properties. Starting from the analysis of recent Chandra observations, here we appraise the current observational situation. We derived limits on the X-ray luminosity of a compact object with a nonthermal, Crab-pulsar-like spectrum of the order of $\approx$$(1$-$5)\times10^{35}$ erg s$^{-1}$, corresponding to limits on the rotational energy loss of a possible X-ray pulsar in SN 1987A of $\approx$$(0.5$-$1.5)\times10^{38}$ erg s$^{-1}$. However, a much brighter X-ray source cannot be excluded if, as is likely, it is enshrouded in a cloud of absorbing matter with metallicity similar to that expected in the outer layers of a massive star towards the end of its life. We found that other limits obtained from various arguments and observations in other energy ranges either are unbinding or allow a similar maximum luminosity of the order of $\approx$$10^{35}$ erg s$^{-1}$. We conclude that while a pulsar alike the one in the Crab Nebula in both luminosity and spectrum is hardly compatible with the observations, there is ample space for an `ordinary' X-ray-emitting young neutron star, born with normal initial spin period, temperature and magnetic field, to be hiding inside the evolving remnant of SN 1987A., Accepted for publication in the Astrophysical Journal; 10 pages, 3 color figures, 2 tables

The properties of a massive star prior to its final explosion are imprinted in the circumstellar medium (CSM) created by its wind and termination shock. We perform a detailed, comprehensive calculation of the time-variable and angle-dependent transmission spectra of an average-luminosity Gamma-Ray Burst (GRB) which explodes in the CSM structure produced by the collapse of a 20 Msun, rapidly rotating, Z=0.001 progenitor star. We study both the case in which metals are initially in the gaseous phase, as well as the situation in which they are heavily depleted into dust. We find that high-velocity lines from low-ionization states of silicon, carbon, and iron are initially present in the spectrum only if the metals are heavily depleted into dust prior to the GRB explosion. However, such lines disappear on timescales of a fraction of a second for a burst observed on-axis, and of a few seconds for a burst seen at high-latitude, making their observation virtually impossible. Rest-frame lines produced in the termination shock are instead clearly visible in all conditions. We conclude that time-resolved, early-time spectroscopy is not a promising way in which the properties of the GRB progenitor wind can be routinely studied. Previous detections of high velocity features in GRB UV spectra must have been due either due to a superposition of a physically unrelated absorber or to a progenitor star with very unusual properties., Comment: Published in MNRAS; higher resolution figures in published version.

We consider radiation emitted by the jitter mechanism in a Blandford-McKee self-similar blastwave. We assume the magnetic field configuration throughout the whole blastwave meets the condition for the emission of jitter radiation and we compute the ensuing images, light curves and spectra. The calculations are performed for both a uniform and a wind environment. We compare our jitter results to synchrotron results. We show that jitter radiation produces slightly different spectra than synchrotron, in particular between the self-absorption and the peak frequency, where the jitter spectrum is flat, while the synchrotron spectrum grows as \nu^{1/3}. The spectral difference is reflected in the early decay slope of the light curves. We conclude that jitter and synchrotron afterglows can be distinguished from each other with good quality observations. However, it is unlikely that the difference can explain the peculiar behavior of several recent observations, such as flat X-ray slopes and uncorrelated optical and X-ray behavior., Comment: 11 pages, 7 postscript figures. Accepted for publication in MNRAS

We compute the average luminosity of X-ray flares as a function of time, for a sample of 10 long-duration gamma-ray burst afterglows. The mean luminosity, averaged over a timescale longer than the duration of the individual flares, declines as a power-law in time with index ~-1.5. We elaborate on the properties of the central engine that can produce such a decline. Assuming that the engine is an accreting compact object, and for a standard conversion factor between accretion rate and jet luminosity, the switch between a neutrino-cooled thin disk and a non-cooled thick disk takes place at the transition from the prompt to the flaring phase. We discuss the implications of this coincidence under different scenarios for the powering of the GRB outflow. We also show that the interaction of the outflow with the envelope of the progenitor star cannot produce flares out of a continuous relativistic flow, and conclude that it is the dynamics of the disk or the jet-launching mechanism that generates an intrinsically unsteady outflow on timescales much longer than the dynamical timescale of the system. This is consistent with the fact that X-ray flares are observed in short-duration GRBs as well as in long-duration ones.

We use the kinetic theory of nucleation to explore the properties of dust nucleation in subsaturated vapours. Due to radiation losses, the subcritical clusters have a smaller temperature compared to their vapour. This alters the dynamical balance between the attachment and detachment of monomers, allowing for stable nucleation of grains in vapours that are subsaturated for their temperature. We find this effect particularly important at low densities and in the absence of a strong background radiation field. We find new conditions for stable nucleation in the n - T phase diagram. The nucleation in the non-local thermodynamic equilibrium (non-LTE) regions is likely to be at much slower rate than in the supersaturated vapours. We evaluate the nucleation rate, warning the reader that it does depend on poorly substantiated properties of the macro-molecules assumed in the computation. On the other hand, the conditions for nucleation depend only on the properties of the large stable grains and are more robust. We finally point out that this mechanism may be relevant in the early Universe as an initial dust pollution mechanism, since once the interstellar medium is polluted with dust, mantle growth is likely to be dominant over non-LTE nucleation in the diffuse medium.

We present MCRaT, a Monte Carlo Radiation Transfer code for self-consistently computing the light curves and spectra of the photospheric emission from relativistic, unmagnetized jets. We apply MCRaT to a relativistic hydrodynamic simulation of a long duration gamma-ray burst jet, and present the resulting light-curves and time-dependent spectra for observers at various angles from the jet axis. We compare our results to observational results and find that photospheric emission is a viable model to explain the prompt phase of long-duration gamma-ray bursts at the peak frequency and above, but faces challenges in reproducing the flat spectrum below the peak frequency. We finally discuss possible limitations of these results both in terms of the hydrodynamics and the radiation transfer and how these limitations could affect the conclusions that we present., Comment: 10 pages, accepted for publication in ApJ

We study the long-term evolution of relativistic jets in collapsars and examine the effects of viewing angle on the subsequent gamma ray bursts. We carry out a series of high-resolution simulations of a jet propagating through a stellar envelope in 2D cylindrical coordinates using the FLASH relativistic hydrodynamics module. For the first time, simulations are carried out using an adaptive mesh that allows for a large dynamic range inside the star while still being efficient enough to follow the evolution of the jet long after it breaks out from the star. Our simulations allow us to single out three phases in the jet evolution: a precursor phase in which relativistic material turbulently shed from the head of the jet first emerges from the star, a shocked jet phase where a fully shocked jet of material is emerging, and an unshocked jet phase where the jet consists of a free-streaming, unshocked core surrounded by a thin boundary layer of shocked jet material. The appearance of these phases will be different to observers at different angles. The precursor has a wide opening angle and would be visible far off axis. The shocked phase has a relatively narrow opening angle that is constant in time. During the unshocked jet phase the opening angle increases logarithmically with time. As a consequence, some observers see prolonged dead times of emission even for constant properties of the jet injected in the stellar core. We also present an analytic model that is able to reproduce the overall properties of the jet and its evolution. We finally discuss the observational implications of our results, emphasizing the possible ways to test progenitor models through the effects of jet propagation in the star. In an appendix, we present 1D and 2D tests of the FLASH relativistic hydrodynamics module.

We present the results of a set of numerical simulations of long-duration gamma-ray burst jets aimed at studying the effect of a variable engine on the peak frequency of the photospheric emission. Our simulations follow the propagation of the jet inside the progenitor star, its break-out, and the subsequent expansion in the environment out to the photospheric radius. A constant and two step-function models are considered for the engine luminosity. We show that our synthetic light-curves follow a luminosity-peak frequency correlation analogous to the Golenetskii correlation found in long-duration gamma-ray burst observations. Within the parameter space explored, it appears that the central engine luminosity profile does not have a significant effect on the location of a gamma-ray burst in the Luminosity-peak frequency plane, bursts from different central engines being indistinguishable from each other.

High resolution spectroscopy of GRB 021004 revealed a wealth of absorption lines from several intermediate ionization species. The velocity structure of the absorber is complex and material with velocity up to >3000 km/s is observed. Since only the blueshifted component is observed, the absorber is very likely to be material closely surrounding the gamma-ray burst. We use a time-dependent photoionization code to track the abundance of the ions over time. Thanks to the presence of absorption from intermediate ionization states at long times, we can estimate the location and mass of the components of the absorber. We interpret those constraints within the hypernova scenario showing that the mass loss rate of the progenitor must have been ~10^{-4} solar masses per year, suggestive of a very massive star. In addition, the wind termination shock must lie at a distance of at least 100 pc, implying a low density environment. The velocity structure of the absorber also requires clumping of the wind at those large distances., Accepted for publication in MNRAS

Short duration (, 6 pages, 2 figures. Accepted for publication in A&A

Gamma-ray burst (GRB) progenitors are observationally linked to the death of massive stars. X-ray studies of the GRB afterglows can deepen our knowledge of the ionization status and metal abundances of the matter in the GRB environment. Moreover, the presence of local matter can be inferred through its fingerprints in the X-ray spectrum, i.e. the presence of absorption higher than the Galactic value. A few studies based on BeppoSAX and XMM-Newton found evidence of higher than Galactic values for the column density in a number of GRB afterglows. Here we report on a systematic analysis of 17 GRBs observed by Swift up to April 15, 2005. We observed a large number of GRBs with an excess of column density. Our sample, together with previous determinations of the intrinsic column densities for GRBs with known redshift, provides evidence for a distribution of absorption consistent with that predicted for randomly occurring GRB within molecular clouds., Comment: Accepted for publication on A&A (5 pages and 2 figures)

In recent years, detailed observations and accurate numerical simulations have provided support to the idea that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms can extract energy from this system and power the SGRBs. Here we connect observations and numerical simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general-relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By assuming a constant efficiency in converting torus mass into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses >~1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some of the SGRBs, the BH would need to have an initial spin of ~0.9, or higher., 6 pages, 3 figures. Extended the discussion on jet efficiency. Accepted for publication on The Astrophysical Journal Letters

The giant flare observed on Dec. 27th 2004 from SGR 1806-20 has revived the idea that a fraction of short ( 1 sec) and a more complex light curve than the Dec. 27th 2004 event. This implies a stringent limit on the rate of extragalactic SGR giant flares with spectral properties analogous to the Dec. 27th flare. We conclude that up to 4 per cent of the short GRBs could be associated to giant flares (2 sigma confidence). This implies that either the distance to SGR 1806-20 is smaller than 15 kpc or the rate of Galactic giant flares is lower than the estimated 0.033 per year.