Can all breaks in gamma-ray burst afterglows be explained by jet effects?

Whether gamma-ray bursts are highly beamed or not is a very important question, as it has been pointed out that the beaming will lead to a sharp break in the afterglow light curves during the ultrarelativistic phase, with the breaking point determined by  Γ∼1/ θ ₀  , where Γ is the bulk Lorentz factor and θ ₀ is the initial half opening angle of the ejecta, and such a break is claimed to be present in the light curves of some GRBs. In this paper we will examine whether all the observed breaks in GRB afterglow light curves can be explained by jet effects. Here we present a detailed calculation of the jet evolution and emission, and have obtained a simple formula of bulk Lorentz factor evolution. We show that the light curves are very smoothly steepened by jet effect, and the shape of the light curve is determined by only one parameter –     , where E and n are the fireball energy and surrounding medium density, respectively. We find that for GRB 990123 and GRB 991216, the jet model can approximately fit their light curves, and the values of     are about 0.17 and 0.22, respectively. On the other hand, the light curves of GRB 990510, GRB 000301c, GRB 000926 and GRB 010222 cannot be fitted by the jet model, which suggests that the breaks may be caused by some other reasons, and the jet effect should be not the unique reason.     

The nature and structure of the emission line nebula K 3-35: a very young planetary nebula with precessing bipolar jet-like outflows?

We present Hα, [N  II ]6583 and 6-cm continuum images of the emission line nebula K 3-35. The optical images reveal an extended nebula (size ≃ 11 × 9 arcsec² in [N  II ]) in which most of the emission originates in a very narrow (width 0.7–1.3 arcsec) S-shaped region which extends almost all along the nebula (≃ 7 arcsec). The 6-cm continuum emission also arises in this narrow region, which is characterized by an exceedingly high point-symmetry and systematic and continuous changes of the orientation with respect to the nebular centre. The properties of the narrow region suggest that it represents a system of precessing bipolar jet-like components. Two low-excitation, compact bipolar knots near the tips of the jet-like components are observed in the deduced [N  II ]/Hα image ratio. These knots may be generated by the interaction of the collimated outflows with surrounding material. A comparison of the optical and radio images shows the existence of differential extinction within the nebula. Maximum extinction is observed in a disc-like region which traces the equator of the elliptical shell previously observed at 20-cm continuum. All available data strongly suggest that K 3-35 is a very young planetary nebula in which we could be observing the first stages of the formation of collimated outflows and point-symmetric structures typically observed in planetary nebulae. The properties of the jet-like components in K 3-35 are in good agreement with models of binary central stars in which highly collimated outflows originate either from a precessing accretion disc or via magnetic collimation in a precessing star.     

Numerical simulations of astrophysical jets from Keplerian discs – III. The effects of mass loading

We present 2.5D time-dependent simulations of the non-linear evolution of non-relativistic outflows from the surface of Keplerian accretion discs. The gas is accelerated from the surface of the disc (which is a fixed platform in these simulations) into a cold corona in stable hydrostatic equilibrium. We explore the dependence of the resulting jet characteristics upon the mass loading of the winds. Two initial configurations of the threading disc magnetic field are studied: a potential field and a uniform vertical field configuration.
We show that the nature of the resulting highly collimated, jet-like outflows (steady or episodic) is determined by the mass load of the disc wind. The mass load controls the interplay between the collimating effects of the toroidal field and the kinetic energy density in the outflow. In this regard, we demonstrate that the onset of episodic behaviour of jets appears to be determined by the quantity     which compares the speed for a toroidal Alfvén wave to cross the diameter of the jet, with the flow speed v _p along the jet. This quantity decreases with increasing load. For sufficiently large N (small mass loads), disturbances appear to grow leading to instabilities and shocks. Knots are then generated and the outflow becomes episodic. These effects are qualitatively independent of the initial magnetic configuration that we employed and are probably generic to a wide variety of magnetized accretion disc models.     

Chimneys in the starburst galaxy M82

We present MERLIN and VLA observations at 1.4 and 5 GHz of the diffuse radio emission in the centre of M82. We detect a large expanding shell of ionized gas surrounding the brightest supernova remnant 41.95+57.5 with a diameter of ∼100 pc and an expansion velocity of ∼100 km s⁻¹. We observe a 50-pc-scale 'blow-out' from this region, in the form of a 'cone' of missing 5-GHz continuum emission, which appears to be an excellent example of a galactic chimney.
On larger radio scales, we observe four chimney structures extending out to the north ∼100–200 pc along the minor axis. One of these features is almost certainly related to the 50-pc-scale blow-out from 41.95+57.5, although this is not the most prominent feature. The other features have also been traced to expulsion of material from the very centre by using an 'unsharp masked' image from 5-GHz VLA B-array observations, with the supernova remnant removed.
We interpret these chimney features as the base of the superwind, which implies that the ejection of material into the halo does not occur smoothly over the starburst region. Instead, very localized venting of hot gas is clearly in evidence.     

Metallicities of high redshift GRB hosts: The case of GRB 100219A

GRB 100219A at z = 4.667 has been the highest redshift gamma‐ray burst observed with the X‐shooter spectrograph up to now. The spectrum covering the range from 5000 to 24000 Å and a large number of absorption lines allows to make a detailed study of the interstellar medium in a high redshift galaxy. The ISM in the low ionisation state and the kinematics of the absorption line components reveal a complex velocity field. The metallicity measured from different absorption lines is around 0.1 solar. Other GRB hosts at redshift beyond ∼3 have similar metallicities albeit with a large scatter in the metallicity distribution. X‐shooter will allow us to determine metallicities of a larger number of GRB hosts beyond redshift 5, to probe the early chemical enrichment of the Universe and to study its evolution from redshift 2 to beyond 10 (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Galactic winds and circulation of the interstellar medium in dwarf galaxies

We study, through 2D hydrodynamical simulations, the feedback of a starburst on the ISM of typical gas-rich dwarf galaxies. The main goal is to address the circulation of the ISM and metals following the starburst. We assume a single-phase rotating ISM in equilibrium in the galactic potential generated by a stellar disc and a spherical dark halo. The starburst is assumed to occur in a small volume in the centre of the galaxy, and it generates a mechanical power of 3.8×10³⁹ or 3.8×10⁴⁰ erg s⁻¹ for 30 Myr. We find, in accordance with previous investigations, that the galactic wind is not very effective in removing the ISM. The metal-rich stellar ejecta, however, can be efficiently expelled from the galaxy and dispersed in the intergalactic medium.
Moreover, we find that the central region of the galaxy is always replenished with cold and dense gas a few 100 million years after the starburst, achieving the requisite for a new star formation event in ≈0.5–1 Gyr. The hydrodynamical evolution of galactic winds is thus consistent with the episodic star formation regime suggested by many chemical evolution studies.
We also discuss the X-ray emission of these galaxies and find that the observable (emission-averaged) abundance of the hot gas underestimates the real one if thermal conduction is effective. This could explain the very low hot-gas metallicities estimated in starburst galaxies.     

A structure and energy dissipation efficiency of relativistic reconfinement shocks

We present a semi-analytical hydrodynamical model for the structure of reconfinement shocks formed in astrophysical relativistic jets interacting with external medium. We take into account exact conservation laws, both across the shock front and in the zone of the shocked matter, and exact angular relations. Our results confirm a good accuracy of the approximate formulae derived by Komissarov & Falle. However, including the transverse pressure gradient in the shocked jet, we predict an absolute size of the shock to be about twice larger. We calculate the efficiency of the kinetic energy dissipation in the shock and show a strong dependence on both the bulk Lorentz factor and opening angle of the jet.     

The 'Carina Flare' supershell: probing the atomic and molecular ISM in a Galactic chimney

The 'Carina Flare' supershell, GSH 287+04−17, is a molecular supershell originally discovered in  ¹²CO( J = 1–0)  with the NANTEN 4 m telescope. We present the first study of the shell's atomic ISM, using H  i 21-cm line data from the Parkes 64-m telescope Southern Galactic Plane Survey. The data reveal a gently expanding,  ∼230 × 360  pc H  i supershell that shows strong evidence of Galactic Plane blowout, with a break in its main body at   z ∼ 280  pc and a capped high-latitude extension reaching   z ∼ 450  pc. The molecular clouds form comoving parts of the atomic shell, and the morphology of the two phases reflects the supershell's influence on the structure of the ISM. We also report the first discovery of an ionized component of the supershell, in the form of delicate, streamer-like filaments aligned with the proposed direction of blowout. The distance estimate to the shell is re-examined, and we find strong evidence to support the original suggestion that it is located in the Carina Arm at a distance of  2.6 ± 0.4 kpc  . Associated H  i and H₂ masses are estimated as   M _H I≈ 7 ± 3 × 10⁵ M_⊙  and     , and the kinetic energy of the expanding shell as   E _K ∼ 1 × 10⁵¹  erg. We examine the results of analytical and numerical models to estimate a required formation energy of several 10⁵¹ to  ∼10⁵²  erg, and an age of  ∼10⁷ yr  . This age is compatible with molecular cloud formation time-scales, and we briefly consider the viability of a supershell-triggered origin for the molecular component.     

Failed gamma-ray bursts and orphan afterglows

It is believed that orphan afterglow searches can help to measure the beaming angle in gamma-ray bursts (GRBs). Great expectations have been put on this method. We point out that the method is in fact not as simple as we originally expected. As a result of the baryon-rich environment that is common to almost all popular progenitor models, there should be many failed gamma-ray bursts, i.e. fireballs with Lorentz factor much less than  100–1000  , but still much larger than unity. In fact, the number of failed gamma-ray bursts may even be much larger than that of successful bursts. Owing to the existence of these failed gamma-ray bursts, there should be many orphan afterglows even if GRBs are due to isotropic fireballs, then the simple discovery of orphan afterglows never means that GRBs are collimated. Unfortunately, to distinguish between a failed-GRB orphan and a jetted but off-axis GRB orphan is not an easy task. The major problem is that the trigger time is unknown. Some possible solutions to the problem are suggested.     

Gemini GMOS/integral field unit spectroscopy of NGC 1569 – II. Mapping the roots of the galactic outflow

We present a set of four Gemini-North Multi-Object Spectrograph/integral field unit (IFU) observations of the central disturbed regions of the dwarf irregular starburst galaxy NGC 1569, surrounding the well-known superstar clusters A and B. This continues on directly from a companion paper, in which we describe the data reduction and analysis techniques employed and present the analysis of one of the IFU pointings. By decomposing the emission-line profiles across the IFU fields, we map out the properties of each individual component identified and identify a number of relationships and correlations that allow us to investigate in detail the state of the ionized interstellar medium (ISM). Our observations support and expand on the main findings from the analysis of the first IFU position, where we conclude that a broad (≲400 km s⁻¹) component underlying the bright nebular emission lines is produced in a turbulent mixing layer on the surface of cool gas knots, set up by the impact of the fast-flowing cluster winds. We discuss the kinematic, electron-density and excitation maps of each region in detail and compare our results to previous studies. Our analysis reveals a very complex environment with many overlapping and superimposed components, including dissolving gas knots, rapidly expanding shocked shells and embedded ionizing sources, but no evidence for organized bulk motions. We conclude that the four IFU positions presented here lie well within the starburst region where energy is injected, and, from the lack of substantial ordered gas flows, within the quasi-hydrostatic zone of the wind interior to the sonic point. The net outflow occurs at radii beyond 100–200 pc, but our data imply that mass-loading of the hot ISM is active even at the roots of the wind.     

The visibility of gamma-ray burst afterglows in dusty star-forming regions

Recent observations of the environments of gamma-ray bursts (GRBs) favour massive stars as their progenitors, which are likely to be surrounded by gas and dust. The visibility of the optical and UV emission of a GRB is expected to depend on the characteristics of both the dust and the GRB emission itself. A reasonable distribution of surrounding dust is capable of absorbing all the optical and UV emission of the optical flash and afterglow of a GRB, unless the optical flash has a peak isotropic luminosity L _peak≳10⁴⁹ erg s⁻¹ . This means that dark bursts should exist and these bursts will have to be studied at infrared rather than optical wavelengths. In this paper details will be given about the infrared GRB dust emission. The reprocessed dust emission peaks at a rest-frame wavelength of about 8 μm. Forthcoming space telescopes, in particular the IRAC camera on board the Space Infrared Telescope Facility , could detect this emission out to a redshift of about two. However, an accurate position of the GRB afterglow must be provided for this emission to be identified, because the light curve of the reprocessed dust emission does not vary on time-scales less than several years.     

Three-dimensional hydrodynamical simulations of the large-scale structure of W50−SS433

We present 3D hydrodynamical simulations of a precessing jet propagating inside a supernova remnant (SNR) shell, particularly applied to the W50−SS433 system in a search for the origin of its peculiar elongated morphology. Several runs were carried out with different values for the mass-loss rate of the jet, the initial radius of the SNR, and the opening angle of the precession cone. We found that our models successfully reproduce the scale and morphology of W50 when the opening angle of the jets is set to 10° or if this angle linearly varies with time. For these models, more realistic runs were made considering that the remnant is expanding into an interstellar medium with an exponential density profile (as H  i observations suggest). Taking into account all these ingredients, the large-scale morphology of the W50−SS433 system, including the asymmetry between the lobes (formed by the jet–SNR interaction), is well reproduced.     

Shattering and coagulation of dust grains in interstellar turbulence

We investigate shattering and coagulation of dust grains in turbulent interstellar medium (ISM). The typical velocity of dust grain as a function of grain size has been calculated for various ISM phases based on a theory of grain dynamics in compressible magnetohydrodynamic turbulence. In this paper, we develop a scheme of grain shattering and coagulation and apply it to turbulent ISM by using the grain velocities predicted by the above turbulence theory. Since large grains tend to acquire large velocity dispersions as shown by earlier studies, large grains tend to be shattered. Large shattering effects are indeed seen in warm ionized medium within a few Myr for grains with radius   a ≳ 10⁻⁶  cm. We also show that shattering in warm neutral medium can limit the largest grain size in ISM  ( a ∼ 2 × 10⁻⁵ cm)  . On the other hand, coagulation tends to modify small grains since it only occurs when the grain velocity is small enough. Coagulation significantly modifies the grain size distribution in dense clouds (DC), where a large fraction of the grains with   a < 10⁻⁶ cm  coagulate in 10 Myr. In fact, the correlation among   R_V   , the carbon bump strength and the ultraviolet slope in the observed Milky Way extinction curves can be explained by the coagulation in DC. It is possible that the grain size distribution in the Milky Way is determined by a combination of all the above effects of shattering and coagulation. Considering that shattering and coagulation in turbulence are effective if dust-to-gas ratio is typically more than ∼1/10 of the Galactic value, the regulation mechanism of grain size distribution should be different between metal-poor and metal-rich environments.     

G309.2–00.6 and jets in supernova remnants

We present Australia Telescope Compact Array observations of the supernova remnant (SNR) G309.2–00.6. In a 1.3-GHz continuum image the remnant appears as a near-circular shell, but with two brightened and distorted arcs of emission on opposite sides. H  i absorption against the SNR yields a distance in the range 5.4 to 14.1 kpc, corresponding to an age (1−20) × 10³ yr. On the basis of the morphology of the SNR we argue that it is a younger analogue of the W 50/SS 433 system, and that its unusual appearance is a result of opposed jets or outflows from a central source. A jet-like feature and breaks in the shell can both be seen along the axis of proposed outflow, providing further support for this interpretation; the central source itself is not detected. The SNR may be interacting with the adjacent H  ii region RCW 80 through an extension of the proposed outflow beyond its shell. This would put the SNR at the lower limit of its distance range and would imply an age 4000 yr. We consider other SNRs similar to G309.2–00.6, and propose remnants whose shells are affected by jets as one of several classes of SNR from which the presence of a central source can be inferred.     

The influence of gravitational acceleration on the supernova-driven Parker instability

Within a framework of 2D magnetohydrodynamic (MHD) simulations, we explore the dynamical regimes initiated by a supernova explosion in a magnetized stratified interstellar medium (ISM). We concentrate on the formation of large-scale magnetic structures and outflows connected with the Parker instability. For the sake of simplicity we only show models with a fixed explosion energy corresponding to a single supernova (SN) occurring in host galaxies with different fixed values of the gravitational acceleration g and different ratios of specific heats. We show that in general, depending on these two parameters, three different regimes are possible: a slowly growing Parker instability on time-scales much longer than the galactic rotation period for small g; the Parker instability growing at roughly the rotation period, which for ratios of specific heats larger than one is accompanied by an outflow resulting from the explosion for intermediate g; and a rapidly growing instability and a strong blow-out flow for large g . By means of numerical simulations and analytical estimates we show that the explosion energy and gravitational acceleration which separate the three regimes scale as Eg ²∼constant in the 2D case. We expect that in the 3D case this scaling law is Eg ³∼constant . Our simulations demonstrate furthermore that a single SN explosion can lead to the growth of multiple Parker loops in the disc and large-scale magnetic field loops in the halo, extending over 2–3 kpc horizontally and up to 3 kpc vertically above the mid-plane of the disc.     

Ly+ and ultraviolet emission from high-redshift gamma-ray burst hosts: to what extent do gamma-ray bursts trace star formation?

We report the result of a search for Lyα emission from the host galaxies of the gamma-ray bursts  (GRBs) 030226 ( z = 1.986), 021004 ( z = 2.335)  and  020124 ( z = 3.198)  . We find that the host galaxy of GRB 021004 is an extended (around 8 kpc) strong Lyα emitter with a rest-frame equivalent width (EW) of 68⁺¹²₋₁₁Å, and a star formation rate of  10.6 ± 2.0 M_⊙ yr⁻¹  . We do not detect the hosts of GRB 030226 and GRB 020124, but the upper limits on their Lyα fluxes do not rule out large rest-frame EWs. In the fields of GRB 021004 and GRB 030226 we find seven and five other galaxies, respectively, with excess emission in the narrow-band filter. These galaxies are candidate Lyα-emitting galaxies in the environment of the host galaxies. We have also compiled a list of all   z ≳ 2  GRB hosts, and demonstrate that a scenario where they trace star formation in an unbiased way is compatible with current observational constraints. Fitting the   z = 3  luminosity function (LF) under this assumption results in a characteristic luminosity of   R *= 24.6  and a faint-end slope of  α=−1.55  , consistent with the LF measured for Lyman-break galaxies.     

The energetics and mass loss of the dwarf starburst Markarian 33

We present ROSAT HRI X-ray data and optical imaging of the important dwarf starburst Markarian 33. We find an extended, complex, shell-like morphology in the X-ray emission, with an extent of ∼     , coincident with the bright star-forming regions at the centre of the galaxy. The physical extent of this X-ray emission from Mrk 33 is very similar to the observed H α emission, and suggests that the bulk of the X-ray emission is coming from an expanding superbubble.
We estimate the age and mass of Mrk 33's starburst to be 5.8 Myr and     respectively, with the energy injection rate in the central regions of the galaxy being ∼10⁴¹ erg s⁻¹, while the associated mass-loss rate from the star-forming regions is estimated to be ∼0.2 M_⊙ yr⁻¹. We suggest that the X-ray emission is predominantly powered by starburst-type activity, and argue that a blow-out in the form of a galactic wind is the most likely fate for Mrk 33, resulting in the loss of most of the galaxy's metal-enriched material and a small fraction (<1 per cent) of the ISM.     

MERLIN polarization observations of compact steep-spectrum sources at 5 GHz

We present polarization observations of 28 compact steep-spectrum sources made with the upgraded MERLIN at 5 GHz. With an angular resolution of 60 milliarcsec and rms noise levels of about 0.1 mJy beam⁻¹, the total intensity images reveal new details in many of these sources. A few of the more extended lobes and jets are quite highly polarized, but more than half of the components are completely unpolarized. Comparison with published data implies that this is due to Faraday depolarization, probably occurring in the surrounding medium with nB  ∼ 1 cm^−3μG. The high resolution of the present observations implies that the variations in Faraday rotation, probably due to magnetic field tangling, occur on scales of less than about 100 pc.     

Cuspy dark matter haloes and the Galaxy

The microlensing optical depth to Baade's Window constrains the minimum total mass in baryonic matter within the Solar circle to be greater than ∼     , assuming the inner Galaxy is barred with viewing angle ∼20°. From the kinematics of solar neighbourhood stars, the local surface density of dark matter is ∼     . We construct cuspy haloes normalized to the local dark matter density and calculate the circular-speed curve of the halo in the inner Galaxy. This is added in quadrature to the rotation curve provided by the stellar and ISM discs, together with a bar sufficiently massive so that the baryonic matter in the inner Galaxy reproduces the microlensing optical depth. Such models violate the observational constraint provided by the tangent-velocity data in the inner Galaxy (typically at radii     . The high baryonic contribution required by the microlensing is consistent with implications from hydrodynamical modelling and the pattern speed of the Galactic bar. We conclude that the cuspy haloes favoured by the cold dark matter cosmology (and its variants) are inconsistent with the observational data on the Galaxy.     

Dynamical and chemical evolution of gas-rich dwarf galaxies

We study the effect of a single, instantaneous starburst on the dynamical and chemical evolution of a gas-rich dwarf galaxy, the potential well of which is dominated by a dark matter halo. We follow the dynamical and chemical evolution of the interstellar medium (ISM) by means of an improved two-dimensional hydrodynamical code coupled with detailed chemical yields originating from type II SNe, type Ia SNe and single low- and intermediate-mass stars (IMS). In particular we follow the evolution of the abundances of H, He, C, N, O, Mg, Si and Fe. We find that for a galaxy resembling IZw18, a galactic wind develops as a consequence of the starburst and it carries out of the galaxy mostly the metal-enriched gas. In addition, we find that different metals are lost differentially in the sense that the elements produced by type Ia SNe are lost more efficiently than others. As a consequence of that, we predict larger [ α /Fe] ratios for the gas inside the galaxy than for the gas leaving the galaxy. A comparison of our predicted abundances of C, N, O and Si in the case of a burst occurring in a primordial gas shows a very good agreement with the observed abundances in IZw18 as long as the burst has an age of ∼31 Myr and IMS produce some primary nitrogen. However, we cannot exclude that a previous burst of star formation had occurred in IZw18, especially if the pre-enrichment produced by the older burst was lower than Z =0.01 Z_⊙. Finally, at variance with previous studies, we find that most of the metals reside in the cold gas phase already after a few Myr. This result is mainly caused by the assumed low SN II heating efficiency, and justifies the generally adopted homogeneous and instantaneous mixing of gas in chemical evolution models.