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.

The evolution and physics of the common envelope (CE) phase are still not well understood. Jets launched from a compact object during this stage may define the evolutionary outcome of the binary system. We focus on the case in which jets are launched from a neutron star (NS) engulfed in the outer layers of a red giant (RG). We run a set of three-dimensional hydrodynamical simulations of jets with different luminosities and inclinations. The luminosity of the jet is self-regulated by the mass accretion rate and an efficiency $\eta$. Depending on the value of $\eta$ the jet can break out of the BHL bulge ("successful jet") and aligns against the incoming wind, in turn, it will realign in favour of the direction of the wind. The jet varies in size and orientation and may present quiescent and active epochs. The inclination of the jet and the Coriolis and centrifugal forces, only slightly affect the global evolution. As the accretion is hypercritical, and the specific angular momentum is above the critical value for the formation of a disk, we infer the formation of a disk and launching of jets. The disks' mass and size would be $\sim$10$^{-2}$~M$_\odot$ and $\gtrsim 10^{10}$ cm, and it may have rings with different rotation directions. In order to have a successful jet from a white dwarf, the ejection process needs to be very efficient ($\eta\sim$0.5). For main sequence stars, there is not enough energy reservoir to launch a successful jet.

The production of both gravitational waves and short gamma-ray bursts (sGRBs) is widely associated with the merger of compact objects. Several studies have modelled the evolution of the electromagnetic emission using the synchrotron emission produced by the deceleration of both a relativistic top-hat jet seen off-axis, and a wide-angle quasi-spherical outflow (both using numerical studies). In this study we present an analytical model of the synchrotron and synchrotron self-Compton (SSC) emission for an off-axis top-hat jet and a quasi-spherical outflow. We calculate the light curves obtained from an analytic model in which the synchrotron and SSC emission (in the fast- or slow-cooling regime) of an off-axis top-hat jet and a quasi-spherical outflow are decelerated in either a homogeneous or a wind-like circumburst medium. We show that the synchrotron emission of the quasi-spherical outflow is stronger than that of the off-axis jet during the first $\sim$ 10 - 20 days, and weaker during the next $\gtrsim$ 80 days. Moreover, we show that if the off-axis jet is decelerated in a wind-like medium, then the SSC emission is very likely to be detected. Applying a MCMC code to our model (for synchrotron emission only), we find the best-fit values for the radio, optical and X-ray emission of GRB 170817A which are in accordance with observations. For GRB 170817A, we find using our model that the synchrotron emission generated by the quasi-spherical outflow and off-axis top-hat jet increase as $F_\nu\propto t^{\alpha}$ with $\alpha\lesssim 0.8$ and $\alpha>3$, respectively. Finally, we obtain the correspondent SSC light curves which are in accordance with the very-high-energy gamma-ray upper limits derived with the GeV - TeV observatories.

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.

Jets launched from a compact object (CO) during a common envelope (CE) may play a key role in the evolution of the system, and may also be an efficient removal channel for its material. In this work we study, through a large set of three-dimensional hydrodynamic simulations, the effects that jets launched from either a black hole (BH) or a neutron star (NS) have during a CE phase. The jet is self-consistently powered by a fraction ($\eta$) of the mass accretion rate that reaches the CO. For low mass accretion efficiencies ($\eta$<0.1\%), the jet is not able to drill through the material accreting onto the CO and forms a bulge around it. For higher efficiencies, the jet is able to drill through the bulge and in intermediate efficiencies ($\eta \sim$1-5\%) may produce an oscillating mass accretion rate behavior. We find that the jet may deposit enough energy to unbind the outer layers of the CE. The self-regulated jets present variability in size and orientation while their cocoons expand smoothly over the CE. The mass accretion rate initially decreases due to the ram pressure of the cocoon. If the launched jet is not able to drill through the bulge, the accretion is such that it may convert a NS into a BH within a decade or double the mass of a BH in a few years. If a jet is present, it deposits enough energy to unbind the outer layers of the CE and the system evolves into a grazing envelope, configuration.

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.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

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.

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.

We report on the development of Mezcal-SRHD, a new adaptive mesh refinement, special relativistic hydrodynamics (SRHD) code, developed with the aim of studying the highly relativistic flows in Gamma-Ray Burst sources. The SRHD equations are solved using finite volume conservative solvers. The correct implementation of the algorithms is verified by one-dimensional (1D) shock tube and multidimensional tests. The code is then applied to study the propagation of 1D spherical impulsive blast waves expanding in a stratified medium with $\rho \propto r^{-k}$, bridging between the relativistic and Newtonian phases, as well as to a two-dimensional (2D) cylindrically symmetric impulsive jet propagating in a constant density medium. It is shown that the deceleration to non-relativistic speeds in one-dimension occurs on scales significantly larger than the Sedov length. This transition is further delayed with respect to the Sedov length as the degree of stratification of the ambient medium is increased. This result, together with the scaling of position, Lorentz factor and the shock velocity as a function of time and shock radius, is explained here using a simple analytical model based on energy conservation. The method used for calculating the afterglow radiation by post-processing the results of the simulations is described in detail. The light curves computed using the results of 1D numerical simulations during the relativistic stage correctly reproduce those calculated assuming the self-similar Blandford-McKee solution for the evolution of the flow. The jet dynamics from our 2D simulations and the resulting afterglow lightcurves, including the jet break, are in good agreement with those presented in previous works. Finally, we show how the details of the dynamics critically depend on properly resolving the structure of the relativistic flow.

The dynamics of GRB jets during the afterglow phase is most reliably and accurately modelled using hydrodynamic simulations. All published simulations, 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 $\rho \propto r^{-k}$ for k = 0, 1, 2. The simulations are performed in 2D using the special relativistic version of the Mezcal code. The dynamics for stratified external media are broadly similar to those derived for expansion into a uniform external medium. The jet half-opening angle start increasing logarithmically with time once the Lorentz factor drops below 1/theta_0. For larger k values the lateral expansion is faster at early times 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 lightcurve for k = 2 although the shape of the break is affected more by the viewing angle than by the slope of the external density profile. Steeper density profiles are found to produce more gradual jet breaks while larger viewing angles cause smoother and later appearing jet breaks. The counter-jet 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 is likely to consistently over-estimate the true energy by up to a factor of a few for k=2, and either over-predict or under-predict it by a smaller factor for k = 0,1.

We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and 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 Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.

Numerical simulations of HH jets never show side-entrainment of environmental material into the jet beam. This is because the bow shock associated with the jet head pushes the surrounding environment into a dense shell, which is never in direct contact with the sides of the jet beam. We present 3D simulations in which a side-streaming motion (representing the motion of the outflow source through the surrounding medium) pushes the post-bow shock shell into direct contact with the jet beam. This is a possible mechanism for modelling well collimated "molecular jets" as an atomic/ionic flow which entrains molecules initially present only in the surrounding environment.

The rotation rate in pre-supernova cores is an important ingredient which 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.

Even thought a large amount of long gamma ray bursts (LGRBs) present quiescent periods, their origin remains unclear. In this talk, it is shown how different angular momentum distributions, as a function of the stellar radius), can lead to neutrino luminosity variability and the possibility of quiescent epochs in LGRBs.

We consider the long term evolution of debris following the tidal disruption of compact stars in the context of short gamma ray bursts (SGRBs). 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 alpha 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 time scale of the new accretion event can reach tens of seconds to minutes, depending on the details of the system. The associated energies and time scales are consistent with those occurring in X-ray flares.

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.