Source author record

Fabio De Colle

Fabio De Colle 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

22works
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

22 published item(s)

preprint2026arXiv

First Results from the PanRadio Gamma-Ray Burst Collaboration: The 400-day Afterglow of GRB 230815A

We introduce the PanRadio Gamma-ray Burst (GRB) program carried out on the Australia Telescope Compact Array: a systematic, multi-year, radio survey of all southern $\textit{Swift}$ GRB events, comprehensively following the multi-frequency evolution of their afterglows from within an hour to years post-burst. We present the results of the 400-day observing campaign following the afterglow of the long-duration (collapsar) GRB 230815A, the first one detected through this program. Typically, GRB 230815A would not otherwise receive traditional radio follow-up, given it has no known redshift and lacks comprehensive multi-wavelength follow-up due to its high line-of-sight extinction with $A_V=2.3$. We found its early X-ray jet break at ${\sim}0.1$ days post-burst to be at odds with the evolution of the multi-frequency radio light curves that were traced over an unusually long duration of $400$ days. The radio light curves approximately evolved (with minor deviations) according to the standard self-similar expansion for a relativistic blast wave in a homogeneous environment prior to the jet break, showing no evidence for evolutions of the microphysical parameters describing the electron acceleration processes. We reconcile these features by proposing a two-component jet: the early X-ray break originates from a narrow component with a half-opening angle ${\sim}2.1^{\circ}$, while the evolution of the radio afterglow stems from a wider component with a half-opening angle of $\gtrapprox 35^{\circ}$. The PanRadio GRB program will establish a sample of comprehensively followed GRBs, where a rigorous inspection of their microphysical and dynamical parameters can be performed, thereby revealing the diversity of features in their outflows and environments.

preprint2022arXiv

Jets in Common Envelopes: a low mass main sequence star in a red giant

We present small-scale three-dimensional hydrodynamical simulations of the evolution of a 0.3Msun main sequence star which launches two perpendicular jets within the envelope of a 0.88Msun red giant. Based on previous large-scale simulations, we study the dynamics of the jets either when the secondary star is grazing, when it has plunged-in, or when it is well-within the envelope of the red giant (in each stage for ~11 days). The dynamics of the jets through the common envelope (CE) depend on the conditions of the environment as well as on their powering. In the grazing stage and the commencement of the plunge self-regulated jets need higher efficiencies to break out of the envelope of the RG. Deep inside the CE, on the timescales simulated, jets are choked independently of whether they are self-regulated or constantly powered. Jets able to break out of the envelope of the RG in large-scale simulations, are choked in our small-scale simulations. The accreted angular momentum onto the secondary star is not large enough to form a disk. The mass accretion onto the MS star is 1-10% of the Bondi-Hoyle-Littleton rate (~1e-3 - 1e-1 Msun yr/s). High luminosity emission, from X-rays to UV and optical, is expected if the jets break out of the CE. Our simulations illustrate the need for inclusion of more realistic accretion and jet models in the dynamical evolution of the CEs.

preprint2022arXiv

The large landscape of supernova, GRB and cocoon interactions

Long gamma ray bursts (LGRBs) are associated to the collapse of a massive star and the formation of a relativistic jet. As the jet propagates through the star, it forms an extended, hot cocoon. The dynamical evolution of the jet/cocoon system and its interaction with the environment has been studied extensively both analytically and numerically. On the other hand, the role played by the supernova (SN) explosion associated with LGRBs in determining the outcome of the system has been barely considered. In this paper, we discuss the large landscape of outcomes resulting from the interaction of the SN, jet and cocoon. We show that the outcome depends mainly on three timescales: the times for the cocoon and supernova shock wave to break through the surface of the progenitor star, and the time needed for the cocoon to engulf completely the progenitor star. The delay between the launch of the SN shock moving through the progenitor star and the jet can be related to these three timescales. Depending on the ordering of these time scales, the jet-cocoon might propagate inside the SN ejecta or the other way around, and the outcome for the properties of the explosion would be different. We discuss the imprint of the complex interaction between the jet-cocoon and the supernova shock on the emergent thermal and non-thermal radiation.

preprint2021arXiv

Survival of the Fittest: Numerical Modeling of Supernova 2014C

Initially classified as a supernova (SN) type Ib, $\sim$ 100 days after the explosion SN\,2014C made a transition to a SN type II, presenting a gradual increase in the H$α$ emission. This has been interpreted as evidence of interaction between the supernova shock wave and a massive shell previously ejected from the progenitor star. In this paper, we present numerical simulations of the propagation of the SN shock through the progenitor star and its wind, as well as the interaction of the SN ejecta with the massive shell. To determine with high precision the structure and location of the shell, we couple a genetic algorithm to a hydrodynamic and a bremsstrahlung radiation transfer code. We iteratively modify the density stratification and location of the shell by minimizing the variance between X-ray observations and synthetic predictions computed from the numerical model. By assuming spherical symmetry, we found that the shell has a mass of 2.6 M$_\odot$, extends from 1.6 $\times 10^{16}$ cm to $1.87 \times 10^{17}$ cm, implying that it was ejected $\sim 60/(v_w/100 {\rm \; km \; s^{-1}})$ yrs before the SN explosion, and has a density stratification decaying as $\sim r^{-3}$. We found that the product of metallicity by the ionization fraction (due to photo-ionization by the post-shock X-ray emission) %and/or the SN UV radiation is $\sim$ 0.5. Finally, we predict that, if the density stratification follows the same power-law behaviour, the SN will break out from the shell by mid 2022, i.e. 8.5 years after explosion.

preprint2020arXiv

Disk formation and jet inclination effects in Common Envelopes

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 $η$. Depending on the value of $η$ 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 ($η\sim$0.5). For main sequence stars, there is not enough energy reservoir to launch a successful jet.

preprint2020arXiv

Modelling the prompt optical emission of GRB 180325A: the evolution of a spike from the optical to gamma-rays

The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with the TAROT and RATIR ground-based telescopes. These observations show two strong optical flashes with separate peaks at $\sim50\;$s and $\sim120\;$s, followed by a temporally extended optical emission. We also present X-rays and gamma-ray observations of GRB 180325A, detected by the Burst Alert Telescope (BAT) and X-ray Telescope (XRT), on the Neil Gehrels Swift observatory, which both observed a narrow flash at $\sim80\;$s. We show that the prompt gamma- and X-ray early emission shares similar temporal and spectral features consistent with internal dissipation within the relativistic outflow (e.g. by internal shocks or magnetic reconnection), while the early optical flashes are likely generated by the reverse shock that decelerates the ejecta as it sweeps up the external medium.

preprint2020arXiv

The Fate of the Merger Remnant in GW170817 and its Imprint on the Jet Structure

The first neutron star binary merger detected in gravitational waves, GW170817 and the subsequent detection of its emission across the electromagnetic spectrum showed that these systems are viable progenitors of short $γ$-ray bursts (sGRB). The afterglow signal of GW170817 has been found to be consistent with a structured GRB jet seen off-axis, requiring significant amounts of relativistic material at large angles. This trait can be attributed to the interaction of the relativistic jet with the external wind medium. Here we perform numerical simulations of relativistic jets interacting with realistic wind environments in order to explore how the properties of the wind and central engine affect the structure of successful jets. We find that the angular energy distribution of the jet depends primarily on the ratio between the lifetime of the jet and the time it takes the merger remnant to collapse. We make use of these simulations to constrain the time it took for the merger remnant in GW170817 to collapse into a black hole based on the angular structure of the jet as inferred from afterglow observations. We conclude that the lifetime of the merger remnant in GW170817 was $\approx 1-1.7$s, which, after collapse, triggered the formation of the jet.

preprint2019arXiv

GRB 180620A: Evidence for late-time energy injection

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 $α_\mathrm{O,decay}=1.32\pm 0.01$, and a $Δ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.

preprint2019arXiv

Jets from Tidal Disruption Events

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, $γ/$X-ray radiation mechanism, jet composition, and the multi-messenger picture. Finally, open questions and future prospects in this field are summarized.

preprint2016arXiv

Transverse velocity shifts in protostellar jets: rotation or velocity asymmetries?

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.

preprint2015arXiv

A Study of Radio Polarization in Protostellar Jets

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.

preprint2014arXiv

A Stellar Wind Origin for the G2 Cloud: Three-Dimensional Numerical Simulations

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$γ$ 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$γ$ luminosity, remain puzzling, and are explained here assuming that the emission is dominated by the free-wind region. The observed Br$γ$ 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$γ$ 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$γ$ emission is more compact after periastron passage.

preprint2014arXiv

Diversity Of Short Gamma-Ray Burst Afterglows From Compact Binary Mergers Hosting Pulsars

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.

preprint2014arXiv

Necessary Conditions for Short Gamma-Ray Burst Production in Binary Neutron Star Mergers

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.

preprint2013arXiv

Conditions For Successful Helium Detonations In Astrophysical Environments

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 ...

preprint2013arXiv

The Morphology and Dynamics of Jet-Driven Supernova Remnants: the Case of W49B

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.

preprint2013arXiv

Understanding the Unusual X-Ray Emission Properties of the Massive, Close Binary WR 20a: A High Energy Window into the Stellar Wind Initiation Region

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.

preprint2012arXiv

The dynamics, appearance and demographics of relativistic jets triggered by tidal disruption of stars in quiescent supermassive black holes

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]

preprint2011arXiv

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

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.

preprint2011arXiv

Gamma-Ray Burst Dynamics and Afterglow Radiation from Adaptive Mesh Refinement, Special Relativistic Hydrodynamic Simulations

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 $ρ\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.

preprint2011arXiv

Simulations of GRB Jets in a Stratified External Medium: Dynamics, Afterglow Lightcurves, Jet Breaks and Radio Calorimetry

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 $ρ\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.

preprint2009arXiv

Radiation-magnetohydrodynamic simulations of the photoionisation of magnetised globules

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.