Long-duration gamma-ray bursts (GRBs) are almost unequivocally associated with very energetic, broad-lined supernovae (SNe) of Type Ic-BL. While the gamma-ray emission is emitted in narrow jets, the SN emits radiation isotropically. Therefore, some SN Ic-BL not associated with GRBs have been hypothesized to arise from events with inner engines such as off-axis GRBs or choked jets. Here we present observations of the nearby ($d = 120$ Mpc) SN 2020bvc (ASAS-SN 20bs) which support this scenario. \textit{Swift} UVOT observations reveal an early decline (up to two days after explosion) while optical spectra classify it as a SN Ic-BL with very high expansion velocities ($\approx$ 70,000 km/s), similar to that found for the jet-cocoon emission in SN 2017iuk associated with GRB 171205A. Moreover, \textit{Swift} X-Ray Telescope and \textit{CXO} X-ray Observatory detected X-ray emission only three days after the SN and decaying onwards, which can be ascribed to an afterglow component. Cocoon and X-ray emission are both signatures of jet-powered GRBs. In the case of SN 2020bvc, we find that the jet is off axis (by $\approx$ 23 degrees), as also indicated by the lack of early ($\approx 1$ day) X-ray emission which explains why no coincident GRB was detected promptly or in archival data. These observations suggest that SN 2020bvc is the first orphan GRB detected through its associated SN emission.

The detection of gravitational waves together with their electromagnetic counterpart, in the gamma-ray burst GRB 170817A, marked a new era of multi-messenger astronomy. Several theoretical models have been proposed to explain the atypical behavior of this event. Recently, it was shown that the multi-wavelength afterglow of GRB 170817A was consistent with a synchrotron forward-shock model when the outflow was viewed off-axis, decelerated in a uniform medium and parametrized through a power-law velocity distribution. Motivated by the upper limits on the very-high-energy emission, and the stratified medium in the close vicinity of a binary neutron star merger proposed to explain the gamma-ray flux in the short GRB 150101B, we extend the mechanism proposed to explain GRB 170817A to a more general scenario deriving the synchrotron self-Compton (SSC) and synchrotron forward-shock model when the off-axis outflow is decelerated in a uniform and stratified circumburst density. As particular cases, we show that the delayed and long-lasting afterglow emission observed in GRB 080503, GRB140903A, GRB 150101B, and GRB 160821B could be interpreted by a similar scenario to the one used to describe GRB 170817A. In addition, we show that the proposed scenario agrees with the MAGIC, Fermi-LAT and H.E.S.S upper limits on gamma-ray emission from GRB 160821B and GRB 170817A.

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 discovery of jets from tidal disruption events (TDEs) rejuvenated the old field of relativistic jets powered by accretion onto supermassive black holes. In this Chapter, we first review the extensive multi-wavelength observations of jetted TDEs. Then, we show that these events provide valuable information on many aspects of jet physics from a new prospective, including the on-and-off switch of jet launching, jet propagation through the ambient medium, $\gamma/$X-ray radiation mechanism, jet composition, and the multi-messenger picture. Finally, open questions and future prospects in this field are summarized.

The early optical emission of gamma-ray bursts gives an opportunity to understand the central engine and first stages of these events. About 30\% of GRBs present flares whose origin is still a subject of discussion. We present optical photometry of GRB 180620A with the COATLI telescope and RATIR instrument. COATLI started to observe from the end of prompt emission at $T+39.3$~s and RATIR from $T+121.4$~s. We supplement the optical data with the X-ray light curve from \emph{Swift}/XRT. %The optical and X-ray light curves show very unusual behavior with features clearly beyond the standard fireball model. We observe an optical flare from $T+110$ to $T+550$~s, with a temporal index decay $\alpha_\mathrm{O,decay}=1.32\pm 0.01$, and a $\Delta t/t=1.63$, which we interpret as the signature of a reverse shock component. After the initial normal decay the light curves show a long plateau from $T+500$ to $T+7800$~s both in X-rays and the optical before decaying again after an achromatic jet break at $T+7800$~s. Fluctuations are seen during the plateau phase in the optical. Adding to the complexity of GRB afterglows, the plateau phase (typically associated with the coasting phase of the jet) is seen in this object after the ``normal'' decay phase (emitted during the deceleration phase of the jet) and the jet break phase occurs directly after the plateau. We suggest that this sequence of events can be explained by a rapid deceleration of the jet with $t_d\lesssim 40$ s due to the high density of the environment ($\approx 100$ cm$^{-3}$) followed by reactivation of the central engine which causes the flare and powers the plateau phase.

The afterglow of GRB$\,$170817A/GW$\,$170817 was very unusual, slowly rising as $F_\nu\propto{}t_{\rm{}obs}^{0.8}\nu^{-0.6}$, peaking at $t_{\rm{obs,pk}}\sim\,150\;$days, and sharply decaying as $\sim{}t_{\rm{}obs}^{-2.2}$. VLBI observations revealed an unresolved radio afterglow image whose flux centroid moved superluminally with $v_{\rm{app}}\approx4c$, clearly indicating that the afterglow was dominated by a relativistic jet's compact core. Different jet angular structures explained the afterglow lightcurves: Gaussian and steep power-law profiles with narrow core angles $\theta_c\lesssim5^\circ$ and larger viewing angles $\theta_{\rm{}obs}/\theta_c\sim3-5$. However, a top-hat jet (with sharp edges at $\theta=\theta_0$) was ruled out since it appeared to produce an early flux rise much steeper than observed. Using 2D relativistic hydrodynamic simulations we show that the initial steep flux rise is an artifact caused by the simulation's finite start time, $t_0$, missing its flux contributions from $t<t_0$ and sometimes "compensated" using an analytic top-hat jet. While an initially top-hat jet is not very physical, such simulations are particularly useful at $t_{\rm{}obs}\gtrsim{}t_{\rm{obs,pk}}$ when the afterglow emission is dominated by the jet's core and becomes insensitive to its exact initial angular profile if it drops off sharply outside of the core. We demonstrate that an initially top-hat jet fits GW$\,$170817/GRB$\,$170817A's afterglow lightcurves and flux centroid motion at $t_{\rm{}obs}\gtrsim{}t_{\rm{obs,pk}}$, for $\theta_{\rm{}obs}/\theta_0\approx3$ and may also fit the earlier lightcurves for $\Gamma_0=\Gamma(t_0)\gtrsim10^{2.5}$. We analytically express the degeneracies between the model parameters, and find a minimal jet energy of $E_{\rm{}min}\approx5.3\times10^{48}\;$erg and circum-burst medium density of $n_{\min}\approx5.3\times10^{-6}~{\rm cm}^{-3}$.

The short GRB 170817A, detected by the Fermi Gamma-ray Burst Monitor, orbiting satellites and ground-based telescopes, was the electromagnetic counterpart of a gravitational-wave transient (GW170817) from a binary neutron star merger. After this merger the $\gamma$-ray light curve exhibited a faint peak at $\sim$ 1.7s and the X-ray, optical and radio light curves displayed an extended emission which increased in brightness up to $\sim$ 160 days. In this paper, we show that the X-ray, optical and radio fluxes are consistent with the synchrotron forward-shock model viewed off-axis when the matter in the outflow is parametrized through a power law velocity distribution. We discuss the origin of the $\gamma$-ray peak in terms of internal and external shocks. We show that the $\gamma$-ray flux might be consistent with a synchrotron self-Compton reverse-shock model observed at high latitudes. Comparing the best-fit values obtained after describing the $\gamma$-ray, X-ray, optical and radio fluxes with our model, we find that the afterglow and $\gamma$-ray emission occurred in different regions and also evidence to propose that the progenitor environment was entrained with magnetic fields and therefore, we argue for the presence of the magnetic field amplification in the binary neutron star merger.

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.

Gamma-ray burst (GRB) jets are narrow, and thus typically point away from us. They are initially ultra-relativistic, causing their prompt $\gamma$-ray and early afterglow emission to be beamed away from us. However, as the jet gradually decelerates its beaming cone widens and eventually reaches our line of sight and the afterglow emission may be detected. Such orphan afterglows were not clearly detected so far. Nevertheless, they should be detected in upcoming optical or radio surveys, and it would be challenging to clearly distinguish between them and other types of transients. Therefore, we perform detailed, realistic calculations of the expected afterglow emission from GRB jets viewed at different angles from the jet's symmetry axis. The dynamics are calculated using 2D relativistic hydrodynamics simulations of jets propagating into different power-law external density profiles, $\rho_{\rm ext}\propto{}R^{-k}$ for $k=0,\,1,\,1.5,\,2$, ranging from a uniform ISM-like medium ($k=0$) to a stratified steady stellar-wind like profile ($k=2$). We calculate radio, optical and X-ray lightcurves, and the evolution of the radio afterglow image size, shape and flux centroid. This may help identify misaligned relativistic jets, whether initially ultra-relativistic and producing a GRB for observers within their beam, or (possibly intrinsically more common) moderately relativistic, in either (i) nearby supernovae Ib/c (some of which are associated with long duration GRBs), or (ii) in binary neutron star mergers, which may produce short duration GRBs, and may also be detected in gravitational waves (e.g. GW$\,$170827/GRB$\,$170817A with a weak prompt $\gamma$-ray emission may harbor an off-axis jet).

Relativistic supernovae constitute a sub-class of type Ic supernovae (SNe). Their non-thermal, radio emission differs notably from that of regular type Ic supernovae as they have a fast expansion speed (with velocities $\sim$ 0.6-0.8 c) which can not be explained by a "standard", spherical SN explosion but advocates for a quickly evolving, mildly relativistic ejecta associated with the SN. In this paper, we compute the synchrotron radiation emitted by the cocoon of a long gamma-ray burst jet (GRB). We show that the energy and velocity of the expanding cocoon, and the radio non-thermal light curves and spectra are consistent with those observed in relativistic SNe. Thus, the radio emission from this events is not coming from the SN shock front, but from the mildly relativistic cocoon produced by the passage of a GRB jet through the progenitor star. We also show that the cocoon radio emission dominates the GRB emission at early times for GRBs seen off-axis, and the flux can be larger at late times compared with on-axis GRBs if the cocoon energy is at least comparable with respect to the GRB energy.

We present hydrodynamic simulations of the hot cocoon produced when a relativistic jet passes through the gamma-ray burst (GRB) progenitor star and its environment, and we compute the lightcurve and spectrum of the radiation emitted by the cocoon. The radiation from the cocoon has a nearly thermal spectrum with a peak in the X-ray band, and it lasts for a few minutes in the observer frame; the cocoon radiation starts at roughly the same time as when $\gamma$-rays from a burst trigger detectors aboard GRB satellites. The isotropic cocoon luminosity ($\sim 10^{47}$ erg s$^{-1}$) is of the same order of magnitude as the X-ray luminosity of a typical long-GRB afterglow during the plateau phase. This radiation should be identifiable in the Swift data because of its nearly thermal spectrum which is distinct from the somewhat brighter power-law component. The detection of this thermal component would provide information regarding the size and density stratification of the GRB progenitor star. Photons from the cocoon are also inverse-Compton (IC) scattered by electrons in the relativistic jet. We present the IC lightcurve and spectrum, by post-processing the results of the numerical simulations. The IC spectrum lies in 10 keV--MeV band for typical GRB parameters. The detection of this IC component would provide an independent measurement of GRB jet Lorentz factor and it would also help to determine the jet magnetisation parameter.

The short-duration ($\lesssim2\;$s) GRB 170817A in the nearby ($D=40\;$Mpc) elliptical galaxy NGC 4993 is the first electromagnetic counterpart of the first gravitational wave (GW) detection of a binary neutron-star (NS-NS) merger. It was followed by optical, IR, and UV emission from half a day up to weeks after the event, as well as late time X-ray and radio emission. The early UV, optical, and IR emission showed a quasi-thermal spectrum suggestive of radioactive-decay powered kilonova-like emission. Comparison to kilonova models favors the formation of a short-lived ($\sim1\;$s) hypermassive NS, which is also supported by the $\Delta t\approx1.74\;$s delay between the GW chirp signal and the prompt GRB onset. However, the late onset of the X-ray (8.9$\;$days) and radio (16.4$\;$days) emission, together with the low isotropic equivalent $\gamma$-ray energy output ($E_{\rm\gamma,iso}\approx5\times10^{46}\;$erg), strongly suggest emission from a narrow relativistic jet viewed off-axis. Here we set up a general framework for off-axis GRB jet afterglow emission, comparing analytic and numerical approaches, and showing their general predictions for short-hard GRBs that accompany binary NS mergers. The prompt GRB emission suggests a viewing angle well outside the jet's core, and we compare the afterglow lightcurves expected in such a case to the X-ray to radio emission from GRB 170817A. We fit an afterglow off-axis jet model to the X-ray and radio data and find that the observations are explained by a viewing angle $\theta_{\rm obs}\approx16^\circ-26^\circ$, GRB jet energy $E\sim10^{48.5}-10^{49.5}~{\rm erg}$, and external density $n\sim10^{-5}-10^{-1}~{\rm cm}^{-3}$ for a $\xi_e\sim 0.1$ non-thermal electron acceleration efficiency.

In the present work we have modeled the supernova remnant (SNR) G296.5+10, by means of 3D magnetohydrodynamics (MHD) simulations. This remnant belongs to the bilateral SNR group and has an additional striking feature: the rotation measure (RM) in its eastern and western parts are very different. In order to explain both the morphology observed in radio-continuum and the RM, we consider that the remnant expands into a medium shaped by the superposition of the magnetic field of the progenitor star with a constant Galactic magnetic field. We have also carried out a polarization study from our MHD results, obtaining synthetic maps of the linearly polarized intensity and the Stokes parameters. This study reveals that both the radio morphology and the reported RM for G$296.5+10$ can be explained if the quasi-parallel acceleration mechanism is taking place in the shock front of this remnant.

Common envelope (CE) is an important phase in the evolution of many binary systems. Giant star / compact object interaction in binaries plays an important role in high-energy phenomena as well as in the evolution of their environment. Material accreted onto the compact object may form a disk and power a jet. We study analytically and through numerical simulations the interaction between the jet and the CE. We determine the conditions under which accreting material quenches the jet or allows it to propagate successfully, in which case even the envelope may be ejected. Close to the stellar core of the companion the compact object accretes at a larger rate. A jet launched from this region needs a larger accretion-to-ejection efficiency to successfully propagate through the CE compared to a jet launched far from the stellar core, and is strongly deflected by the orbital motion. The energy deposited by the jet may be larger than the binding energy of the envelope. The jet can, thus, play a fundamental role in the CE evolution. We find that the energy dissipation of the jet into the CE may stop accretion onto the disk. We expect the jet to be intermittent, unless the energy deposited is large enough to lead to the unbinding of the outer layers of the CE. Given that the energy and duration of the jet are similar to those of ultra-long GRBs, we suggest this as a new channel to produce these events.

Three dimensional magnetohydrodynamical simulations were carried out in order to perform a new polarization study of the radio emission of the supernova remnant SN 1006. These simulations consider that the remnant expands into a turbulent interstellar medium (including both magnetic field and density perturbations). Based on the referenced-polar angle technique, a statistical study was done on observational and numerical magnetic field position-angle distributions. Our results show that a turbulent medium with an adiabatic index of 1.3 can reproduce the polarization properties of the SN 1006 remnant. This statistical study reveals itself as a useful tool for obtaining the orientation of the ambient magnetic field, previous to be swept up by the main supernova remnant shock.

Observations of several protostellar jets show systematic differences in radial velocity transverse to the jet propagation direction, which have been interpreted as evidence of rotation in the jets. In this paper we discuss the origin of these velocity shifts, and show that they could be originated by rotation in the flow, or by side to side asymmetries in the shock velocity, which could be due to asymmetries in the jet ejection velocity/density or in the ambient medium. For typical poloidal jet velocities (~ 100-200 km/s), an asymmetry >~ 10% can produce velocity shifts comparable to those observed. We also present three dimensional numerical simulations of rotating, precessing and asymmetric jets, and show that, even though for a given jet there is a clear degeneracy between these effects, a statistical analysis of jets with different inclination angles can help to distinguish between the alternative origins of transverse velocity shifts. Our analysis indicate that side to side velocities asymmetries could represent an important contribution to transverse velocity shifts, being the most important contributor for large jet inclination angles (with respect the the plane of the sky), and can not be neglected when interpreting the observations.

The most popular model for short gamma-ray bursts (sGRBs) involves the coalescence of binary neutron stars. Because the progenitor is actually hidden from view, we must consider under which circumstances such merging systems are capable of producing a successful sGRB. Soon after coalescence, winds are launched from the merger remnant. In this paper, we use realistic wind profiles derived from global merger simulations in order to investigate the interaction of sGRB jets with these winds using numerical simulations. We analyze the conditions for which these axisymmetric winds permit relativistic jets to breakout and produce a sGRB. We find that jets with luminosities comparable to those observed in sGRBs are only successful when their half-opening angles are below ~20{\deg}. This jet collimation mechanism leads to a simple physical interpretation of the luminosities and opening angles inferred for sGRBs. If wide, low luminosity jets are observed, they might be indicative of a different progenitor avenue such as the merger of a neutron star with a black hole. We also use the observed durations of sGRB to place constraints on the lifetime of the wind phase, which is determined by the time it takes the jet to breakout. In all cases we find that the derived limits argue against completely stable remnants for binary neutron star mergers that produce sGRBs.

Synchrotron radiation is commonly observed in connection with shocks of different velocities, ranging from relativistic shocks associated with active galactic nuclei, gamma-ray bursts or microquasars to weakly- or non-relativistic flows as those observed in supernova remnants. Recent observations of synchrotron emission in protostellar jets are important not only because they extend the range over which the acceleration process works, but also because they allow to determine the jet and/or interstellar magnetic field structure, thus giving insights on the jet ejection and collimation mechanisms. In this paper, we compute for the first time polarized (synchrotron) and non polarized (thermal X-ray) synthetic emission maps from axisymmetrical simulations of magnetized protostellar jets. We consider models with different jet velocities and variability, as well as toroidal or helical magnetic field. Our simulations show that variable, low-density jets with velocities $\sim$ 1000 km s$^{-1}$ and $\sim$ 10 times lighter than the environment can produce internal knots with significant synchrotron emission, and thermal X-rays in the shocked region of the leading bow shock moving in a dense medium. While models with a purely toroidal magnetic field show a very large degree of polarization, models with helical magnetic field show lower values and a decrease of the degree of polarization, in agreement with observations of protostellar jets.

A very important issue concerning protostellar jets is the mechanism behind their formation. Obtaining information on the region at the base of a jet can shed light into the subject and some years ago this has been done through a search for a rotational signature at the jet line spectrum. The existence of such signatures, however, remains controversial. In order to contribute to the clarification of this issue, in this paper we show that the Principal Component Analysis (PCA) can potentially help to distinguish between rotation and precession effects in protostellar jet images. We apply the PCA to synthetic spectro-imaging datacubes generated as an output of numerical simulations of protostellar jets. In this way we generate a benchmark to which a PCA diagnostics of real observations can be confronted. Using the computed emission line profiles for [O I]6300A and [S II]6716A, we recover and analyze the effects of rotation and precession in tomograms generated by PCA. We show that different combinations of the eigenvectors can be used to enhance and to identify the rotation features present in the data. Our results indicate that the PCA can be useful for disentangling rotation from precession in jets with an inclination of the jet with respect to the plane of the sky as high as 45 degrees. We have been able to recover the initially imposed rotation jet profile for models at moderate inclination angle (< 15 degrees) and without precession (abridged).

We use three dimensional magnetohydrodynamic (MHD) simulations to model the supernova remnant SN 1006. From our numerical results, we have carried out a polarization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter $Q$, and the polar-referenced angle, which can be compared with observational results. Synthetic maps were computed considering two possible particle acceleration mechanisms: quasi-parallel and quasi-perpendicular. The comparison of synthetic maps of the Stokes parameter $Q$ maps with observations proves to be a valuable tool to discern unambiguously which mechanism is taking place in the remnant of SN 1006, giving strong support to the quasi-parallel model.

Short gamma-ray bursts (sGRBs) are widely believed to result from the mergers of compact binaries. This model predicts an afterglow that bears the characteristic signatures of a constant, low density medium, including a smooth prompt-afterglow transition, and a simple temporal evolution. However, these expectations are in conflict with observations for a non-negligible fraction of sGRB afterglows. In particular, the onset of the afterglow phase for some of these events appears to be delayed and, in addition, a few of them exhibit late- time rapid fading in their lightcurves. We show that these peculiar observations can be explained independently of ongoing central engine activity if some sGRB progenitors are compact binaries hosting at least one pulsar. The Poynting flux emanating from the pulsar companion can excavate a bow-shock cavity surround- ing the binary. If this cavity is larger than the shock deceleration length scale in the undisturbed interstellar medium, then the onset of the afterglow will be delayed. Should the deceleration occur entirely within the swept-up thin shell, a rapid fade in the lightcurve will ensue. We identify two types of pulsar that can achieve the conditions necessary for altering the afterglow: low field, long lived pulsars, and high field pulsars. We find that a sizable fraction (~20-50%) of low field pulsars are likely to reside in neutron star binaries based on observations, while their high field counterparts are not. Hydrodynamical calculations motivated by this model are shown to be in good agreement with observations of sGRB afterglow lightcurves.

We present 3D, adaptive mesh refinement simulations of G2, a cloud of gas moving in a highly eccentric orbit towards the galactic center. We assume that G2 originates from a stellar wind interacting with the environment of the Sgr A* black hole. The stellar wind forms a cometary bubble which becomes increasingly elongated as the star approaches periastron. A few months after periastron passage, streams of material begin to accrete on the central black hole with accretion rates $\dot{M} \sim 10^{-8}$ M$_\odot$ yr$^{-1}$. Predicted Br$\gamma$ emission maps and position-velocity diagrams show an elongated emission resembling recent observations of G2. A large increase in luminosity is predicted by the emission coming from the shocked wind region during periastron passage. The observations, showing a constant Br$\gamma$ luminosity, remain puzzling, and are explained here assuming that the emission is dominated by the free-wind region. The observed Br$\gamma$ luminosity ($\sim 8 \times 10^{30}$ erg s$^{-1}$) is reproduced by a model with a $v_w=50$ km s$^{-1}$ wind velocity and a $10^{-7}$ M$_\odot$ yr$^{-1}$ mass loss rate if the emission comes from the shocked wind. A faster and less dense wind reproduces the Br$\gamma$ luminosity if the emission comes from the inner, free wind region. The extended cometary wind bubble, largely destroyed by the tidal interaction with the black hole, reforms a few years after periastron passage. As a result, the Br$\gamma$ emission is more compact after periastron passage.

The central engine of short gamma-ray bursts (sGRBs) is hidden from direct view, operating at a scale much smaller than that probed by the emitted radiation. Thus we must infer its origin not only with respect to the formation of the {\it trigger} - the actual astrophysical configuration that is capable of powering a sGRB - but also from the consequences that follow from the various evolutionary pathways that may be involved in producing it. Considering binary neutron star mergers we critically evaluate, analytically and through numerical simulations, whether the neutrino-driven wind produced by the newly formed hyper-massive neutron star can allow the collimated relativistic outflow that follows its collapse to actually produce a sGRB or not. Upon comparison with the observed sGRB duration distribution, we find that collapse cannot be significantly delayed (<=100 ms) before the outflow is choked, thus limiting the possibility that long-lived hyper-massive remnants can account for these events. In the case of successful breakthrough of the jet through the neutrino-driven wind, the energy stored in the cocoon could contribute to the precursor and extended emission observed in sGRBs.

We report on a VLA survey for late-time radio emission from 59 supernovae (SNe) of Type I b/c, which have been associated with long-duration gamma-ray bursts (GRBs). An "off-axis" GRB burst (i.e. whose relativistic jet points away from us) is expected to have late-time radio emission even in the absence of significant prompt gamma-ray emission. From our sample, we detected only SN 2003gk with an 8.4-GHz flux density of $2260 \pm 130 \,\mu$Jy. Our subsequent VLBI observations of SN 2003gk, at an age of $\sim$8 yr, allowed us to determine its radius to be $(2.4 \pm 0.4) \times 10^{17}$ cm, or $94 \pm 15$ light days. This radius rules out relativistic expansion as expected for an off-axis GRB jet, and instead suggests an expansion speed of $\sim 10\:000$ km s$^{-1}$ typical for non-relativistic core-collapse supernovae. We attribute the late-onset radio emission to interaction of the ejecta with a dense shell caused by episodic mass-loss from the progenitor. In addition, we present new calculations for the expected radio lightcurves from GRB jets at various angles to the line of sight, and compare these to our observed limits on the flux densities of the remainder of our SN sample. From this comparison we can say that only a fraction of broadlined Type I b/c SNe have a radio-bright jet similar to those seen for GRB afterglows at cosmological distances. However, we also find that for a reasonable range of parameters, as might be representative of the actual population of GRB events rather than the detected bright ones, the radio emission from the GRB jets can be quite faint, and that at present, radio observations do not place strong constraints on off-axis GRB jets.

The problem of explaining the X-ray emission properties of the massive, close binary WR 20a is discussed. Located near the cluster core of Westerlund 2, WR 20a is composed of two nearly identical Wolf- Rayet stars of 82 and 83 solar masses orbiting with a period of only 3.7 days. Although Chandra observations were taken during the secondary optical eclipse, the X-ray light curve shows no signs of a flux decrement. In fact, WR 20a appears slightly more X-ray luminous and softer during the optical eclipse, opposite to what has been observed in other binary systems. To aid in our interpretation of the data, we compare with the results of hydrodynamical simulations using the adaptive mesh refinement code Mezcal that includes radiative cooling and a radiative acceleration force term. It is shown that the X-ray emission can be successfully explained in models where the wind-wind collision interface in this system occurs while the outflowing material is still being accelerated. Consequently, WR 20a serves as a critical test-case for how radiatively-driven stellar winds initiate and interact. Our models not only procure a robust description of current Chandra data, which cover the orbital phases between 0.3 to 0.6, but provide detailed predictions over the entire orbit.

The circumstellar medium (CSM) of a massive star is modified by its winds before a supernova (SN) explosion occurs, and thus the evolution of the resulting supernova remnant (SNR) is influenced by both the geometry of the explosion as well as the complex structure of the CSM. Motivated by recent work suggesting the SNR W49B was a jet-driven SN expanding in a complex CSM, we explore how the dynamics and the metal distributions in a jet-driven explosion are modified by the interaction with the surrounding environment. In particular, we perform hydrodynamical calculations to study the dynamics and explosive nucleosynthesis of a jet-driven SN triggered by the collapse of a 25 solar masses Wolf-Rayet star and its subsequent interaction with the CSM up to several hundred years following the explosion. We find that although the CSM has small-scale effects on the structure of the SNR, the overall morphology and abundance patterns are reflective of the initial asymmetry of the SN explosion. Thus, we predict that jet-driven SNRs, such as W49B, should be identifiable based on morphology and abundance patterns at ages up to several hundred years, even if they expand into a complex CSM environment.

Several models for type Ia-like supernovae events rely on the production of a self-sustained detonation powered by nuclear reactions. In the absence of hydrogen, the fuel that powers these detonations typically consists of either pure helium (He) or a mixture of carbon and oxygen (C/O). Studies that systematically determine the conditions required to initiate detonations in C/O material exist, but until now no analogous investigation of He matter has been conducted. We perform one-dimensional reactive hydrodynamical simulations at a variety of initial density and temperature combinations and find critical length scales for the initiation of He detonations that range between 1 - $10^{10}$ cm. A simple estimate of the length scales over which the total consumption of fuel will occur for steady-state detonations is provided by the Chapman-Jouguet (CJ) formalism. Our initiation lengths are consistently smaller than the corresponding CJ length scales by a factor of $\Sim 100$, providing opportunities for thermonuclear explosions in a wider range of low-mass white dwarfs (WDs) than previously thought possible. We find that virialized WDs with as little mass as 0.24 $M_\odot$ can be detonated, and that even less massive WDs can be detonated if a sizable fraction of their mass is raised to a higher adiabat. That the initiation length is exceeded by the CJ length implies that certain systems may not reach nuclear statistical equilibrium within the time it takes a detonation to traverse the object. In support of this hypothesis, we demonstrate that incomplete burning will occur in the majority of He WD detonations and that $^{40}$Ca, $^{44}$Ti, or $^{48}$Cr, rather than $^{56}$Ni, is the predominant burning product for many of these events. We anticipate that ...

We examine the consequences of a model in which relativistic jets can be triggered in quiescent massive black holes when a geometrically thick and hot accretion disk forms as a result of the tidal disruption of a star. To estimate the power, thrust and lifetime of the jet, we use the mass accretion history onto the black hole as calculated by detailed hydrodynamic simulations of the tidal disruption of stars. We go on to determine the states of the interstellar medium in various types of quiescent galactic nuclei, and describe how this external matter can affect jets propagating through it. We use this information, together with a two-dimensional hydrodynamic model of the structure of the relativistic flow, to study the dynamics of the jet, the propagation of which is regulated by the density stratification of the environment and by its injection history. The breaking of symmetry involved in transitioning from one to two dimensions is crucial and leads to qualitatively new phenomena. Many of the observed properties of the Swift 1644+57/GRB 110328A event can be understood as resulting from accretion onto and jets driven by a $10^6 M_\odot$ central mass black hole following the disruption of sun-like star. With the inclusion of a stochastic contribution to the luminosity due to variations in the feeding rate driven by instabilities near the tidal radius, we find that our model can explain the X-ray light curve without invoking a rarely-occurring deep encounter. In conjunction with the number density of black holes in the local universe, we hypothesize that the conditions required to produce the Swift event are not anomalous, but are in fact representative of the jet-driven flare population arising from tidal disruptions. [abridged]

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.

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

We present the results of radiation-magnetohydrodynamic simulations of the expansion of HII regions and surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs, including the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border. The field in the neutral gas lies along the ionisation front, while the field in the ionised gas tends to be perpendicular to this. The highest pressure compressed neutral/molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition, but with magnetic energy densities 1 to 3 orders of magnitude lower. A high velocity dispersion (approx 8 km/s) is maintained in the ionised gas throughout our simulations, despite the mean expansion velocity being significantly lower. The B field does not significantly brake the HII region expansion on the length and timescales accessible to our simulations, but it does tend to suppress the small-scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the HII region. However, the relative luminosity of ionising and non-ionising radiation has a much larger influence than the presence or absence of the B field. When the radiation field is relatively soft (as in the case of a lower mass cluster, with earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms. Movies available at http://youtube.com/user/divBequals0

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.

The physical parameters of Herbig-Haro jets are usually determined from emission line ratios, obtained from spectroscopy or narrow band imaging, assuming that the emitting region is homogeneous along the line of sight. Under the more general hypothesis of axisymmetry, we apply tomographic reconstruction techniques to the analysis of Herbig-Haro jets. We use data of the HH30 jet taken by Hartigan & Morse (2007) with the Hubble space telescope using the slitless spectroscopy technique. Using a non-parametric Tikhonov regularization technique, we determine the volumetric emission line intensities of the [SII]6716,6731, [OI]6300 and [NII]6583 forbidden emission lines. From our tomographic analysis of the corresponding line ratios, we produce "three-dimensional" images of the physical parameters. The reconstructed density, temperature and ionization fraction present much steeper profiles than those inferred using the assumption of homogeneity. Our technique reveals that the reconstructed jet is much more collimated than the observed one close to the source (a width ~ 5 AU vs. ~ 20 AU at a distance of 10 AU from the star), while they have similar widths at larger distances. In addition, our results show a much more fragmented and irregular jet structure than the classical analysis, suggesting that the the ejection history of the jet from the star-disk system has a shorter timescale component (~ some months) superimposed on a longer, previously observed timescale (of a few years). Finally, we discuss the possible application of the same technique to other stellar jets and planetary nebulae.

In this work we employ an MHD numerical code to reproduce the morphology observed for SN 1006 in radio synchrotron and thermal X-ray emission. We introduce a density discontinuity, in the form of a flat cloud parallel to the Galactic Plane, in order to explain the NW filament observed in optical wavelengths and in thermal X-rays. We compare our models with observations. We also perform a test that contrasts the radio emitting bright limbs of the SNR against the central region, finding additional support to our results. Our main conclusion is that the most probable direction of the ambient magnetic field is on average perpendicular to the Galactic Plane.

We present the first three-dimensional radiation-magnetohydrodynamic simulations of the photoionisation of a dense, magnetised molecular globule by an external source of ultraviolet radiation. We find that, for the case of a strong ionising field, significant deviations from the non-magnetic evolution are seen when the initial magnetic field threading the globule has an associated magnetic pressure that is greater than one hundred times the gas pressure. In such a strong-field case, the photoevaporating globule will adopt a flattened or "curled up" shape, depending on the initial field orientation, and magnetic confinement of the ionised photoevaporation flow can lead to recombination and subsequent fragmentation during advanced stages of the globule evolution. We find suggestive evidence that such magnetic effects may be important in the formation of bright, bar-like emission features in H II regions. We include simple but realistic fits to heating and cooling rates in the neutral and molecular gas in the vicinity of a high-mass star cluster and show that the frequently used isothermal approximation can lead to an overestimate of the importance of gravitational instability in the radiatively imploded globule. For globules within 2 parsecs of a high-mass star cluster, we find that heating by stellar x rays prevents the molecular gas from cooling below 50 K.

The dynamical effects of magnetic fields in models of radiative, Herbig-Haro (HH) jets have been studied in a number of papers. For example, magnetized, radiative jets from variable sources have been studied with axisymmetric and 3D numerical simulations. In this paper, we present an analytic model describing the effect of a toroidal magnetic field on the internal working surfaces that result from a variability in the ejection velocity. We find that for parameters appropriate for HH jets the forces associated with the magnetic field dominate over the gas pressure force within the working surfaces. Depending on the ram pressure radial cross section of the jet, the magnetic field can produce a strong axial pinch, or, alternatively, a broadening of the internal working surfaces. We check the validity of the analytic model with axisymmetric numerical simulations of variable, magnetized jets.

Stellar jets are normally constituted by chains of knots with some periodicity in their spatial distribution, corresponding to a variability of order of several years in the ejection from the protostar/disk system. A widely accepted theory for the presence of knots is related to the generation of internal working surfaces due to variations in the jet ejection velocity. In this paper we study the effect of variations in the inner disk-wind radius on the jet ejection velocity. We show that a small variation in the inner disk-wind radius produce a variation in the jet velocity large enough to generate the observed knots. We also show that the variation in the inner radius may be related to a variation of the stellar magnetic field.

We present a 3D numerical simulation of the recently discovered cometary structure produced as Mira travels through the galactic ISM. In our simulation, we consider that Mira ejects a steady, latitude-dependent wind, which interacts with a homogeneous, streaming environment. The axisymmetry of the problem is broken by the lack of alignment between the direction of the relative motion of the environment and the polar axis of the latitude-dependent wind. With this model, we are able to produce a cometary head with a ``double bow shock'' which agrees well with the structure of the head of Mira's comet. We therefore conclude that a time-dependence in the ejected wind is not required for reproducing the observed double bow shock.

CONTEXT: In the interpretation of stellar jet observations, the physical parameters are usually determined from emission line ratios, obtained from spectroscopic observations or using the information contained in narrow band images. The basic hypothesis in the interpretation of the observations is that the emitting region is homogeneous along the line of sight. Actually, stellar jets are in general not homogeneous, and therefore line of sight convolution effects may lead to the main uncertainty in the determination of the physical parameters. AIMS: This paper is aimed at showing the systematic errors introduced when assuming an homogeneous medium, and studying the effect of an inhomogeneous medium on plasma diagnostics for the case of a stellar jet. In addition, we explore how to reconstruct the volumetric physical parameters of the jet (i. e., with dependence both across and along the line of sight). METHODS: We use standard techniques to determine the physical parameters, i. e., the electron density, temperature and hydrogen ionisation fraction across the jet, and a multi-Gaussian method to invert the Abel transform and determine the reconstructed physical structure. RESULTS: When assuming an homogeneous medium the physical parameters, integrated along the line of sight, do not represent the average of the true values, and do not have a clear physical interpretation. We show that when some information is available on the emissivity profile across the jet, it is then possible to obtain appropriate derivations of the electron density, temperature and ionisation fraction.

Structure and spectral energy distributions (SEDs) of externally irradiated circumstellar disks are often computed on the basis of the two-temperature model of Chiang & Goldreich. We refine these calculations by using a more realistic temperature profile which is continuous at all optical depths and thus goes beyond the two-temperature model. It is based on the approximate solution of the radiation transfer in the disk obtained from the frequency-integrated moment equations in the Eddington approximation. We come up with a simple procedure (``constant g_z approximation'') for treating the vertical structure of the disk in regions where its optical depth to stellar radiation is high. This allows us to obtain expressions for the vertical profiles of density and pressure at every point in the disk and to determine the shape of its surface. Armed with these analytical results we calculate the full radial structure of the disk and demonstrate that it favorably agrees with the results of direct numerical calculations. We also describe a simple and efficient way of the SED calculation based on our adopted temperature profile. Resulting spectra provide very good match (especially at short wavelengths) to the results of more detailed (but also more time-consuming) SED calculations solving the full frequency- and angle-dependent radiation transfer within the disk.

We study the H$\alpha$ emission from jets using two-dimensional axisymmetrical simulations. We compare the emission obtained from hydrodynamic (HD) simulations with that obtained from magnetohydrodynamics (MHD) simulations. The magnetic field is supposed to be present in the jet only, and with a toroidal configuration. The simulations have time-dependent ejection velocities and different intensities for the initial magnetic field. The results show an increase in the H$\alpha$ emission along the jet for the magnetized cases with respect to the HD case. The increase in the emission is due to a better collimation of the jet in the MHD case, and to a small increase in the shock velocity. These results could have important implications for the interpretation of the observations of jets from young stellar objects.

Using the Yguazu-a three-dimensional hydrodynamic code, we have computed a set of numerical simulations of heavy, supersonic, radiatively cooling jets including variabilities in both the ejection direction (precession) and the jet velocity (intermittence). In order to investigate the effects of jet rotation on the shape of the line profiles, we also introduce an initial toroidal rotation velocity profile, in agreement with some recent observational evidence found in jets from T Tauri stars which seems to support the presence of a rotation velocity pattern inside the jet beam, near the jet production region. Since the Yguazu-a code includes an atomic/ionic network, we are able to compute the emission coefficients for several emission lines, and we generate line profiles for the H, [O I]6300, [S II]6716 and [N II]6548 lines. Using initial parameters that are suitable for the DG Tau microjet, we show that the computed radial velocity shift for the medium-velocity component of the line profile as a function of distance from the jet axis is strikingly similar for rotating and non-rotating jet models. These findings lead us to put forward some caveats on the interpretation of the observed radial velocity distribution from a few outflows from young stellar objects, and we claim that these data should not be directly used as a doubtless confirmation of the magnetocentrifugal wind acceleration models.