Diego López-Cámara

Phase Transitions and He-Synthesis-Driven Winds in Neutrino Cooled Accretion Disks: Prospects for Late Flares in Short Gamma-Ray Bursts

We consider the long-term evolution of debris following the tidal disruption of compact stars in the context of short gamma ray bursts. The initial encounter impulsively creates a hot, dense, neutrino-cooled disk capable of powering the prompt emission. After a long delay, we find that powerful winds are launched from the surface of the disk, driven by the recombination of free nucleons into α-particles. The associated energy release depletes the mass supply and eventually shuts off activity of the central engine. As a result, the luminosity and mass accretion rate deviate from the earlier self-similar behavior expected for an isolated ring with efficient cooling. This then enables a secondary episode of delayed activity to become prominent as an observable signature, when material in the tidal tails produced by the initial encounter returns to the vicinity of the central object. The timescale of the new accretion event can reach tens of seconds to minutes, depending on the details of the system. The associated energies and timescales are consistent with those occurring in X-ray flares.

SIMULATIONS OF GAMMA-RAY BURST JETS IN A STRATIFIED EXTERNAL MEDIUM: DYNAMICS, AFTERGLOW LIGHT CURVES, JET BREAKS, AND RADIO CALORIMETRY

The dynamics of gamma-ray burst (GRB) jets during the afterglow phase is most reliably and accurately modeled using hydrodynamic simulations. All published simulations so far, however, have considered only a uniform external medium, while a stratified external medium is expected around long duration GRB progenitors. Here, we present simulations of the dynamics of GRB jets and the resulting afterglow emission for both uniform and stratified external media with ρext ∝ r −k for k = 0, 1, 2. The simulations are performed in two dimensions using the special relativistic version of the Mezcal code. Common to all calculations is the initiation of the GRB jet as a conical wedge of half-opening angle θ0 = 0.2 whose radial profile is taken from the self-similar Blandford‐McKee solution. The dynamics for stratified external media (k = 1, 2) are broadly similar to those derived for expansion into a uniform external medium (k = 0). The jet half-opening angle is observed to start increasing logarithmically with time (or radius) once the Lorentz factor Γ drops below θ −1 0 . For larger k values, however, the lateral expansion is faster at early times (when Γ >θ −1 0 ) and slower at late times with the jet expansion becoming Newtonian and slowly approaching spherical symmetry over progressively longer timescales. We find that, contrary to analytic expectations, there is a reasonably sharp jet break in the light curve for k = 2 (a wind-like external medium), although the shape of the break is affected more by the viewing angle (for θobs θ0) than by the slope of the external density profile (for 0 k 2). Steeper density profiles (i.e., increasing k values) are found to produce more gradual jet breaks while larger viewing angles cause smoother and later appearing jet breaks. The counterjet becomes visible as it becomes sub-relativistic, and for k = 0 this results in a clear bump-like feature in the light curve. However, for larger k values the jet decelerates more gradually, causing only a mild flattening in the radio light curve that might be hard to discern when k = 2. Late-time radio calorimetry, which makes use of a spherical flow approximation near the non-relativistic transition, is likely to consistently overestimate the true energy by up to a factor of a few for k = 2, but likely to either overpredict or underpredict it by a smaller factor for k = 0, 1.

THREE-DIMENSIONAL ADAPTIVE MESH REFINEMENT SIMULATIONS OF LONG-DURATION GAMMA-RAY BURST JETS INSIDE MASSIVE PROGENITOR STARS

We present the results of special relativistic, adaptive mesh refinement, 3D simulations of gamma-ray burst jets expanding inside a realistic stellar progenitor. Our simulations confirm that relativistic jets can propagate and break out of the progenitor star while remaining relativistic. This result is independent of the resolution, even though the amount of turbulence and variability observed in the simulations is greater at higher resolutions. We find that the propagation of the jet head inside the progenitor star is slightly faster in 3D simulations compared to 2D ones at the same resolution. This behavior seems to be due to the fact that the jet head in 3D simulations can wobble around the jet axis, finding the spot of least resistance to proceed. Most of the average jet properties, such as density, pressure, and Lorentz factor, are only marginally affected by the dimensionality of the simulations and therefore results from 2D simulations can be considered reliable.

Late time afterglow observations reveal a collimated relativistic jet in the ejecta of the binary neutron star merger GW170817

The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multiwavelength 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 doubt 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 multiwavelength 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° 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.

GAMMA-RAY BURST PRODUCTION AND SUPERNOVA SIGNATURES IN SLOWLY ROTATING COLLAPSARS

We consider accretion onto newborn black holes following the collapse of rotating massive stellar cores, at the threshold where a centrifugally supported disk gives way to nearly radial inflow for low angular momentum. For realistic initial conditions taken from pre-supernova (pre-SN) evolution calculations, the densities and temperatures involved require the use of a detailed equation of state and neutrino cooling processes, as well as a qualitative consideration of the effects of general relativity. Through two-dimensional dynamical calculations we show how the energy release is affected by the rotation rate and the strength of angular momentum transport, giving rise to qualitatively different solutions in limits of high and low angular momentum, each being capable of powering a gamma-ray burst (GRB). We explore the likelihood of producing Fe-group elements in the two regimes and suggest that while large and massive centrifugally supported disks are capable of driving strong outflows with a possible SN-like signature, quasi-radial flows lack such a feature and may produce a GRB without such an accompanying feature, as seen in GRB060505.

Off-axis Prompt X-Ray Transients from the Cocoon of Short Gamma-Ray Bursts

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<theta<45 degrees and gives origin to X-ray transients with a duration of a few to ~10 seconds, delayed by a few seconds from the time of the merger. We also discuss the distance at which such transients can be detected, finding that it depends sensitively on the assumptions that are made about the radiation spectrum. Purely thermal cocoon transients are detectable only out to a few Mpc, Comptonized transients can instead be detected by the FERMI GBM out to several tens of Mpc.

CRITICAL ANGULAR MOMENTUM DISTRIBUTIONS IN COLLAPSARS: QUIESCENT PERIODS FROM ACCRETION STATE TRANSITIONS IN LONG GAMMA-RAY BURSTS

The rotation rate in pre-supernova cores is an important ingredient that can profoundly affect the post-collapse evolution and associated energy release in supernovae and long gamma ray bursts (LGRBs). Previous work has focused on whether the specific angular momentum is above or below the critical value required for the creation of a centrifugally supported disk around a black hole. Here, we explore the effect of the distribution of angular momentum with radius in the star and show that qualitative transitions between high and low angular momentum flow, corresponding to high and low luminosity accretion states, can effectively be reflected in the energy output, leading to variability and the possibility of quiescent times in LGRBs.

A MODEL OF MIRA'S COMETARY HEAD/TAIL ENTERING THE LOCAL BUBBLE

We model the cometary structure around Mira as the interaction of an asymptotic giant branch stellar wind from Mira A with a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment, Mira entered the Local Bubble (low-density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a >100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Miras comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.

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.

THREE-DIMENSIONAL SIMULATIONS OF LONG DURATION GAMMA-RAY BURST JETS: TIMESCALES 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.