Source author record

Davide Lazzati

Davide Lazzati appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2022arXiv

Dynamics of a relativistic jet through magnetized media

The merger of two neutron stars (NSs) produces the emission of gravitational waves, the formation of a compact object surrounded by a dense and magnetized environment, and the launching of a collimated and relativistic jet, which will eventually produce a short gamma-ray burst (SGRB). The interaction of the jet with the environment has been shown to play a major role in shaping the structure of the outflow that eventually powers the gamma-ray emission. In this paper, we present a set of 2.5 dimensional RMHD simulations that follow the evolution of a relativistic non-magnetized jet through a medium with different magnetization levels, as produced after the merger of two NSs. We find that the predominant consequence of a magnetized ambient medium is that of suppressing instabilities within the jet, and preventing the formation of a series of collimation shocks. One implication of this is that internal shocks lose efficiency, causing bursts with low-luminosity prompt emission. On the other hand, the jet-head velocity and the induced magnetization within the jet are fairly independent from the magnetization of the ambient medium. Future numerical studies with a larger domain are necessary to obtain light curves and spectra in order to better understand the role of magnetized media.

preprint2022arXiv

The emergence of diffused Gamma-Ray Burst afterglows from the disks of Active Galactic Nuclei

The disks of Active Galactic Nuclei (AGNs) have emerged as rich environments for the production and capture of stars and the compact objects that they leave behind. These stars produce long Gamma-Ray Bursts (LGRBs) at their deaths, while frequent interactions among compact objects form binary neutron stars and neutron star-black hole binaries, leading to short GRBs (SGRBs) upon their merger. Predicting the properties of these transients as they emerge from the dense environments of AGN disks is key to their proper identification and to better constrain the star and compact object population in AGN disks. Some of these transients would appear unusual because they take place in much higher densities than the interstellar medium. Others, which are the subject of this paper, would additionally be modified by radiation diffusion since they are generated within optically thick regions of the accretion disks. Here we compute the GRB afterglow light curves for diffused GRB sources for a representative variety of central black-hole masses and disk locations. We find that the radiation from radio to UV and soft X-rays can be strongly suppressed by synchrotron self-absorption in the dense medium of the AGN disk. In addition, photon diffusion can significantly delay the emergence of the emission peak, turning a beamed, fast transient into a slow, isotropic, and dimmer one. These would appear as broadband-correlated AGN variability with dominance at the higher frequencies. Their properties can constrain both the stellar populations within AGN disks as well as the disk structure.

preprint2021arXiv

Host galaxies and electromagnetic counterparts to binary neutron star mergers across the cosmic time: Detectability of GW170817-like events

The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to maximize observational prospects. Here we perform population synthesis calculations of BNS evolution with the code MOBSE, and seed the binaries in galaxies at three representative redshifts (z=0.01,0.1,1) of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the astrophysical parameters of GRB170817A as a prototype, we numerically compute the broadband lightcurves of jets from BNS mergers, with the afterglow brightness depending on the local medium density at the merger sites. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find that ~70-80% of BNSs from z=0.01 should be detectable in gamma-rays. The afterglow detection probabilities of GW-triggered BNS mergers vary between ~0.3-0.7%, with higher values for jets aligned with the BNS angular momentum, and are comparable across the high and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses.

preprint2021arXiv

Photospheric Prompt Emission From Long Gamma Ray Burst Simulations -- II. Spectropolarimetry

Although Gamma Ray Bursts (GRBs) have been detected for many decades, the lack of knowledge regarding the radiation mechanism that produces the energetic flash of radiation, or prompt emission, from these events has prevented the full use of GRBs as probes of high energy astrophysical processes. While there are multiple models that attempt to describe the prompt emission, each model can be tuned to account for observed GRB characteristics in the gamma and X-ray energy bands. One energy range that has not been fully explored for the purpose of prompt emission model comparison is that of the optical band, especially with regards to polarization. Here, we use an improved MCRaT code to calculate the expected photospheric optical and gamma-ray polarization signatures ($Π_\mathrm{opt}$ and $Π_γ$, respectively) from a set of two relativistic hydrodynamic long GRB simulations, which emulate a constant and variable jet. We find that time resolved $Π_\mathrm{opt}$ can be large ($\sim 75\%$) while time-integrated $Π_\mathrm{opt}$ can be smaller due to integration over the asymmetries in the GRB jet where optical photons originate; $Π_γ$ follows a similar evolution as $Π_\mathrm{opt}$ with smaller polarization degrees. We also show that $Π_\mathrm{opt}$ and $Π_γ$ agree well with observations in each energy range. Additionally, we make predictions for the expected polarization of GRBs based on their location within the Yonetoku relationship. While improvements can be made to our analyses and predictions, they exhibit the insight that global radiative transfer simulations of GRB jets can provide with respect to current and future observations.

preprint2020arXiv

Electromagnetic Signatures of Relativistic Explosions in AGN Disks

The disks of Active Galactic Nuclei (AGNs), traditionally studied as the feeders of the supermassive black holes (SMBHs) at their centers, have recently triggered a lot of interest also as 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 and NS-BH mergers are expected. However, their properties in the unique environments of AGN disks have never been studied. Here we show that GRBs in AGNs 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, the external shock develops before the internal shocks, leading to prompt emission powered by a relativistic reverse shock. The transient's time evolution is also compressed, yielding afterglow emission that is much brighter and peaks much 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 is trapped inside the disc photosphere. These diffusion-dominated transients emerge on timescales of days in disks around SMBHs of $\sim 10^6 M_\odot$ to years for SMBHs of $\sim 10^8 M_\odot$. Finally, a large fraction of events, especially in AGNs with SMBHs $\lesssim 10^7 M_\odot$, display time-variable absorption in the X-ray band.

preprint2020arXiv

Intrinsic properties of the engine and jet that powered the short gamma-ray burst associated with GW170817

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.

preprint2020arXiv

Photospheric Polarization Signatures From Long Gamma Ray Burst Simulations

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.

preprint2020arXiv

Polarization constraints on the geometry of the magnetic field in the external shock of gamma-ray bursts

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.

preprint2016arXiv

Formation and properties of astrophysical carbonaceous dust

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.

preprint2016arXiv

On radiative acceleration in spine-sheath structured blazar jets

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.

preprint2016arXiv

Short Gamma-Ray Bursts from the Merger of Two Black Holes

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.

preprint2016arXiv

Three-dimensional simulations of long duration gamma-ray burst jets: time scales from variable engines

Gamma-ray burst light curves are characterized by marked variability, each showing unique properties. The origin of such variability, at least for a fraction of long GRBs, may be the result of an unsteady central engine. It is thus important to study the effects that an episodic central engine has on the jet propagation and, eventually, on the prompt emission within the collapsar scenario. Thus, in this study we follow the interaction of pulsed outflows with their progenitor stars with hydrodynamic numerical simulations in both two and three dimensions. We show that the propagation of unsteady jets is affected by the interaction with the progenitor material well after the break-out time, especially for jets with long quiescent times, comparable to or larger than a second. We also show that this interaction can lead to an asymmetric behavior in which pulse durations and quiescent periods are systematically different. After the pulsed jets drill through the progenitor and the interstellar medium, we find that, on average, the quiescent epochs last longer than the pulses (even in simulations with symmetrical active and quiescent engine times). This could explain the asymmetry detected in the light curves of long-quiescent times duration gamma-ray bursts.

preprint2015arXiv

Gamma-ray burst spectra and spectral correlations from sub-photospheric Comptonization

One of the most important unresolved issues in gamma-ray burst physics is the origin of the prompt gamma-ray spectrum. Its general non-thermal character and the softness in the X-ray band remain unexplained. We tackle these issues by performing Monte Carlo simulations of radiation-matter interactions in a scattering dominated photon-lepton plasma. The plasma -- initially in equilibrium -- is driven to non-equilibrium conditions by a sudden energy injection in the lepton population, mimicking the effect of a shock wave or the dissipation of magnetic energy. Equilibrium restoration occurs due to energy exchange between the photons and leptons. While the initial and final equilibrium spectra are thermal, the transitional photon spectra are characterized by non-thermal features such as power-law tails, high energy bumps, and multiple components. Such non-thermal features are observed at infinity if the dissipation occurs at small to moderate optical depths, and the spectrum is released before thermalization is complete. We model the synthetic spectra with a Band function and show that the resulting spectral parameters are similar to observations for a frequency range of 2-3 orders of magnitude around the peak. In addition, our model predicts correlations between the low-frequency photon index and the peak frequency as well as between the low- and high-frequency indices. We explore baryon and pair dominated fireballs and reach the conclusion that baryonic fireballs are a better model for explaining the observed features of gamma-ray burst spectra.

preprint2015arXiv

Numerical Simulations of Gamma-Ray Burst Explosions

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.

preprint2015arXiv

The interplay between chemistry and nucleation in the formation of carbonaceous dust in supernova ejecta

Core-collapse supernovae are considered to be important contributors to the primitive dust enrichment of the interstellar medium in the high-redshift universe. Theoretical models of dust formation in stellar explosions have so far provided controversial results and a generally poor fit to the observations of dust formation in local supernovae. We present a new methodology for the calculation of carbonaceous dust formation in young supernova remnants. Our new technique uses both the nucleation theory and a chemical reaction network to allow us to compute the dust growth beyond the molecular level as well as to consider chemical erosion of the forming grains. We find that carbonaceous dust forms efficiently in the core of the ejecta, but takes several years to condensate, longer than previously estimated. It forms unevenly and remains concentrated in the inner part of the remnant. These results support the role of core-collapse supernovae as dust factories and provide new insight on the observations of SN 1987A, in which large amounts of dust have been detected to form on a timescale of years after core collapse.

preprint2014arXiv

Formation and properties of astrophysical carbonaceous dust. I: ab-initio calculations of the configuration and binding energies of small carbon clusters

The binding energies of n < 100 carbon clusters are calculated using the ab-initio density functional theory code Quantum Espresso. Carbon cluster geometries are determined using several levels of classical techniques and further refined using density functional theory. The resulting energies are used to compute the work of cluster formation and the nucleation rate in a saturated, hydrogen-poor carbon gas. Compared to classical calculations that adopt the capillary approximation, we find that nucleation of carbon clusters is enhanced at low temperatures and depressed at high temperatures. This difference is ascribed to the different behavior of the critical cluster size. We find that the critical cluster size is at n = 27 or n = 8 for a broad range of temperatures and saturations, instead of being a smooth function of such parameters. The results of our calculations can be used to follow carbonaceous cluster/grain formation, stability, and growth in hydrogen poor environments, such as the inner layers of core-collapse supernova and supernova remnants.

preprint2014arXiv

Photospheric emission from long duration gamma-ray bursts powered by variable engines

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.

preprint2013arXiv

On the observed duration distribution of gamma-ray bursts from collapsars

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.

preprint2012arXiv

Compact Binary Progenitors of Short Gamma-Ray Bursts

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.

preprint2012arXiv

Unifying the Zoo of Jet-Driven Stellar Explosions

We present a set of numerical simulations of stellar explosions induced by relativistic jets emanating from a central engine sitting at the center of compact, dying stars. We explore a wide range of durations of the central engine activity, two candidate stellar progenitors, and two possible values of the total energy release. We find that even if the jets are narrowly collimated, their interaction with the star unbinds the stellar material, producing a stellar explosion. We also find that the outcome of the explosion can be very different depending on the duration of the engine activity. Only the longest-lasting engines result in successful gamma-ray bursts. Engines that power jets only for a short time result in relativistic supernova explosions, akin to observed engine-driven SNe such as SN2009bb. Engines with intermediate durations produce weak gamma-ray bursts, with properties similar to nearby bursts such as GRB 980425. Finally, we find that the engines with the shortest durations, if they exist in nature, produce stellar explosions that lack sizable amounts of relativistic ejecta and are therefore dynamically indistinguishable from ordinary core-collapse supernovae.

preprint2011arXiv

High-efficiency photospheric emission of long-duration gamma-ray burst jets: the effect of the viewing angle

We present the results of a numerical investigation of the spectra and light curves of the emission from the photospheres of long-duration gamma-ray burst jets. We confirm that the photospheric emission has high efficiency and we show that the efficiency increases slightly with the off-axis angle. We show that the peak frequency of the observed spectrum is proportional to the square root of the photosphere's luminosity, in agreement with the Amati relation. However, a quantitative comparison reveals that the thermal peak frequency is too small for the corresponding total luminosity. As a consequence, the radiation must be out of thermal equilibrium with the baryons in order to reproduce the observations. Finally, we show that the spectrum integrated over the emitting surface is virtually indistinguishable from a Planck law, and therefore an additional mechanism has to be identified to explain the non-thermal behavior of the observed spectra at both high and low frequencies.

preprint2011arXiv

On the effects of microphysical grain properties on the yields of carbonaceous dust from type II SNe

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.

preprint2010arXiv

Non-Thermal emission from the photospheres of Gamma-Ray Burst outflows. I: High frequency tails

We study the spectrum of high frequency radiation emerging from mildly dissipative photospheres of long-duration gamma-ray burst outflows. Building on the results of recent numerical investigations, we assume that electrons are heated impulsively to mildly relativistic energies by either shocks or magnetic dissipation at Thomson optical depths of several and subsequently cool by inverse Compton, scattering off the thermal photons of the photosphere. We show that even in the absence of magnetic field and non-thermal leptons, inverse Compton scattering produces power-law tails that extend from the peak of the thermal radiation, at several hundred keV, to several tens of MeV, and possibly up to GeV energies. The slope of the high-frequency power-law is predicted to vary substantially during a single burst, and the model can easily account for the diversity of high-frequency spectra observed by BATSE. Our model works in baryonic as well as in magnetically dominated outflows, as long as the magnetic field component is not overwhelmingly dominant.

preprint2010arXiv

Short-duration gamma-ray bursts from off-axis collapsars

We present 2D high-resolution hydrodynamic simulations of the relativistic outflows of long-duration gamma-ray burst progenitors. We analyze the properties of the outflows at wide off-axis angles, produced by the expansion of the hot cocoon that surrounds the jet inside the progenitor star. We find that the cocoon emission at wide angles may have properties similar to those of the subclass of short-duration gamma-ray bursts with persistent X-ray emission. We compute the predicted duration distribution, redshift distribution, and afterglow brightness and we find that they are all in agreement with the observed properties of short GRBs with persistent emission. We suggest that a SN component, the properties of the host galaxies, and late afterglow observations can be used as a crucial test to verify this model.

preprint2010arXiv

The origin and propagation of variability in the outflows of long duration gamma-ray bursts

We present the results of hydrodynamical simulations of gamma-ray burst jets propagating through their stellar progenitor material and subsequently through the surrounding circumstellar medium. We consider both jets that are injected with constant properties in the center of the star and jets injected with a variable luminosity. We show that the variability properties of the jet outside the star are a combination of the variability injected by the engine and the variability caused by the jet propagation through the star. Comparing power spectra for the two cases shows that the variability injected by the engine is preserved even if the jet is heavily shocked inside the star. Such shocking produces additional variability at long time scales, of order several seconds. Our findings suggest that the broad pulses of several seconds duration typically observed in gamma-ray bursts are due to the interaction of the jet with the progenitor, while the short-timescale variability, characterized by fluctuations on time scales of milliseconds, has to be injected at the base of the jet. Studying the properties of the fast variability in GRBs may therefore provide clues to the nature of the inner engine and the mechanisms of energy extraction from it.

preprint2010arXiv

The rest-frame ultraviolet spectra of GRBs from massive rapidly-rotating stellar progenitors

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.

preprint2010arXiv

Thermal fluctuations and nanoscale effects in the nucleation of carbonaceous dust grains

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.

preprint2010arXiv

X-ray flares from propagation instabilities in long Gamma-Ray Burst jets

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.

preprint2009arXiv

Polarization signature of gamma-ray bursts from fragmented fireballs

We study the polarization properties of the prompt emission of gamma-ray bursts produced by fragmented fireballs. Such fireballs, known in the literature under various names, are made by the superposition of many individual blobs, each of which produces a spike in the light curve. The differences between pulses are due both to the intrinsic diversity in the blobs' properties and to their orientation with respect to the line of sight. We show that the peak flux and the polarization of each pulse are connected through the orientation of the blob that produces the pulse, while the position angle fluctuates randomly from one pulse to the next and is constant within an individual pulse. The most polarized pulses are those with approximately one tenth the peak flux of the brightest pulse. These conclusions do not depend on the assumed radiation mechanism nor on the energy and Lorentz factor of the blobs. We compare the prediction of this model to a simulated set of observations, showing that a limited sample of GRBs with time-resolved polarization measurements would provide a crucial test for this model. We finally show that a hint of the predicted correlation may have been already observed.

preprint2009arXiv

X-ray Polarization of Gamma-Ray Bursts

The degree and the temporal evolution of linear polarization in the prompt and afterglow emission of gamma-ray bursts is a very robust diagnostic of some key features of gamma-ray bursts jets and their micro and macro physics. In this contribution, I review the current status of the theory of polarized emission from GRB jets during the prompt, optical flash, and afterglow emission. I compare the theoretical predictions to the available observations and discuss the future prospect from both the theoretical and observational standpoints.

preprint2007arXiv

Dust Echoes from the Ambient Medium of Gamma-Ray Bursts

Long gamma-ray bursts (GRBs) are likely associated with the collapse of massive stars, which produce dust and are born in dusty environments. Absorption and scattering of ultraviolet/X-ray photons from the prompt, optical flash and afterglow emission of the GRB produce dust echoes. We perform time-dependent calculations of these echoes, accounting for the evolution of the dust grain distribution due to selective grain destruction by the GRB radiation, and for off-axis beaming. We explore cloud configurations of differing density and size -- the echo light curve and spectrum depend on the cloud radius, with larger clouds peaking at longer wavelengths. For a region ~ 3 pc in size with hydrogen density ~ 1000 per cubic centimeter, the echo spectrum peaks at ~ 3.6 microns and ~ 8.8 eV for thermal and scattered components, respectively. Dust echoes should be detectable with the Very Large Telescope up to z ~ 0.1, IRAC onboard the Spitzer Space Telescope up to z ~ 0.2, and NICMOS onboard the Hubble Space Telescope up to z ~ 0.3. Furthermore, the shape of the echo light curve allows one to infer: the jet opening angle; the inclination of the jet axis with respect to the line of sight; the size of the dust-emitting region. For sources with symmetric, bipolar jets, dust echoes exhibit two bumps in the light curve, making them easily distinguishable from the rebrightening due to an underlying supernova.

preprint2007arXiv

Jitter radiation as a possible mechanism for Gamma-Ray Burst afterglows. Spectra and lightcurves

The standard model of GRB afterglows assumes that the radiation observed as a delayed emission is of synchrotron origin, which requires the shock magnetic field to be relatively homogeneous on small scales. An alternative mechanism -- jitter radiation, which traditionally has been applied to the prompt emission -- substitutes synchrotron when the magnetic field is tangled on a microscopic scale. Such fields are produced at relativistic shocks by the Weibel instability. Here we explore the possibility that small-scale fields populate afterglow shocks. We derive the spectrum of jitter radiation under the afterglow conditions. We also derive the afterglow lightcurves for the ISM and Wind profiles of the ambient density. Jitter self-absorption is calculated here for the first time. We find that jitter radiation can produce afterglows similar to synchrotron-generated ones, but with some important differences. We compare the predictions of the two emission mechanisms. By fitting observational data to the synchrotron and jitter afterglow lightcurves, it can be possible to discriminate between the small-scale vs large-scale magnetic field models in afterglow shocks.

preprint2007arXiv

Non-LTE dust nucleation in sub-saturated vapors

We use the kinetic theory of nucleation to explore the properties of dust nucleation in sub-saturated vapors. Due to radiation losses, the sub-critical clusters have a smaller temperature compared to their vapor. This alters the dynamical balance between attachment and detachment of monomers, allowing for stable nucleation of grains in vapors that are sub-saturated 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-LTE regions is likely to be at much slower rate than in the super-saturated vapors. 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.

preprint2007arXiv

Temporal and Angular Properties of GRB Jets Emerging from Massive Stars

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.

preprint2006arXiv

The physics of GRB jets and their interaction with the progenitor star

It is now generally accepted that long gamma-ray bursts are associated with the final evolutionary stages of massive stars. As a consequence, their jets must propagate through the stellar progenitor and break out on their surface, before they can reach the photospheric radius and produce the gamma-ray photons. We investigate the role of the progenitor star in shaping the jet properties. We show that even a jet powered by a steady engine can develop a rich phenomenology at the stellar surface. We present special-relativistic simulations and compare the results to analytic considerations. We show that the jet is complex in the time as well as in the angular domain, so that observers located along different lines of sight detect significantly different bursts.

preprint2005arXiv

GRB Progenitors and Environment

The study and knowledge of the environment of Gamma-Ray Bursts is of great interest from many points of view. For high redshift (z>0.5) events, the structure of the ambient medium is one of the best indicators of the nature and properties of the progenitor. It also tells us about the last stages of the pre-explosion evolution of the progenitor. In addition, it is interesting in its own as a sample of the interstellar medium in a high redshift galaxy. Measures of the density and structure of the GRB environment are however sparse, and different methods yield different, often incompatible, results. I will review the methods and results with particular emphasis on the case of GRB 021004, a puzzling but highly informative event. I will finally underline the advancements that will be possible in the Swift era.

preprint2004arXiv

Gamma-Ray Burst Progenitors Confront Observations

The discovery of a supernova emerging at late times in the afterglow of GRB 030329 has apparently settled the issue on the nature of the progenitor of gamma-ray bursts. We now know that at least a fraction of cosmological GRBs are associated with the death of massive stars, and that the two explosions are most likely simultaneous. Even though the association was already suggested for GRB 980425, the peculiarity of that burst did not allow to extend the association to all GRBs. The issue is now to understand whether GRB 030329 is a "standard burst" or not. I will discuss some peculiarities of GRB 030329 and its afterglow lightcurve showing how, rather than a classical cosmological GRB, it looks more like a transition object linking weak events like GRB 980425 to the classical long duration GRBs. I will also discuss the problems faced by the Hypernova scenario to account for the X-ray features detected in several GRBs and their afterglows.

preprint1998arXiv

On the location of X-ray peaks and dominant galaxies in clusters

We analyze the structure of the X-ray emission of a sample of 22 Abell clusters of galaxies with a cD in their centre, observed with the ROSAT PSPC. Utilizing the multi-scale power of the Wavelet Transform we detect significant (~ 50 h^{-1} kpc) offsets between the large scale centroid and the peak of X-ray emission. Despite the uncertainties on the satellite pointing, the X-ray to optical correlation indicates a likely association between the X-ray peak and the dominant galaxy. We develop a model in which the offset is produced by small amplitude oscillations of the cD galaxy around the bottom of the cluster potential well and successfully compare it to the observed distribution of offsets. Within this scenario the offsets are not due to dynamic instabilities and the number of structured clusters is greatly reduced.