3D numerical simulations of a radiative Herbig–Haro jet in a supersonic side wind

We have performed 3D numerical simulations of an over-pressurized Herbig–Haro-type jet which propagates into a sidestreaming environment. The interaction between the jet and the sidewind results in a perpendicular acceleration of the jet material, and a consequent curvature of the jet as it moves into the anisotropic medium. We find that an approximately steady configuration is achieved both for a sidewind that is perpendicular to the jet and for a sidewind inclined at 45° towards the jet source. The curvature obtained in both these models is consistent with analytic models of the jet/sidewind problem.   We have also calculated Hα maps, which show an emitting sheath around the upwind (with respect to the sidewind) side of the jet beam. This emitting sheath may explain part of the observed emission from curved stellar jets.     

Chemistry in luminous AGN and starburst galaxies

Molecular line emission is a useful tool for probing the highly obscured inner kpc of starburst galaxies and buried AGNs. Molecular line ratios serve as diagnostic tools of the physical conditions of the gas—but also of its chemical properties. Both provide important clues to the type and evolutionary stage of the nuclear activity. While CO emission remains the main tracer for molecular distribution and dynamics, molecules such as HCN, HNC, HCO⁺, CN and HC₃N are useful for probing the properties of the denser (n≳10⁴ cm⁻³), star-forming gas. Here I discuss current views on how line emission from these species can be interpreted in luminous galaxies. HNC, HCO⁺ and CN are all species that can be associated both with photon dominated regions (PDRs) in starbursts—as well as X-ray dominated regions (XDRs) associated with AGN activity. HC₃N line emission may identify galaxies where the starburst is in the early stage of its evolution.     

Dust extinction and X-ray emission from the starburst galaxy NGC 1482

We present the results based on multiwavelength imaging observations of the prominent dust lane starburst galaxy NGC 1482 aimed to investigate the extinction properties of dust existing in the extreme environment. (B-V) colour-index map derived for the starburst galaxy NGC 1482 confirms two prominent dust lanes running along its optical major axis and are found to extend up to ∼11 kpc. In addition to the main lanes, several filamentary structures of dust originating from the central starburst are also evident. Though, the dust is surrounded by exotic environment, the average extinction curve derived for this target galaxy is compatible with the Galactic curve, with R_V = 3.05, and imply that the dust grains responsible for the optical extinction in the target galaxy are not really different than the canonical grains in the Milky Way. Our estimate of total dust content of NGC 1482 assuming screening effect of dust is ∼2.7 × 10⁵ M_⊙, and provide lower limit due to the fact that our method is not sensitive to the intermix component of dust. Comparison of the observed dust in the galaxy with that supplied by the SNe to the ISM, imply that this supply is not sufficient to account for the observed dust and hence point towards the origin of dust in this galaxy through a merger like event.Our multiband imaging analysis reveals a qualitative physical correspondence between the morphologies of the dust and Hα emission lines as well as diffuse X-ray emission in this galaxy. Spatially resolved spectral analysis of the hot gas along outflows exhibit a gradient in the temperature. Similar gradient was also noticed in the measured values of metallicity, indicating that the gas in the halo is not yet enriched. High resolution, 2-8 keV Chandra image reveals a pair of point sources in the nuclear region with their luminosities equal to 2.27 × 10³⁹ erg s⁻¹ and 9.34 × 10³⁹ erg s⁻¹, and are in excess of the Eddington-limit of 1.5 M_⊙ accreting source. Spectral analysis of these sources exhibit an absorbed-power law with the hydrogen column density higher than that derived from the optical measurements.     

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.     

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.     

Jet simulations extending radially self-similar magnetohydrodynamics models

We perform a numerical simulation of magnetohydrodynamics (MHD) radially self-similar jets, whose prototype is the Blandford & Payne analytical example. The final steady state that is reached is valid close to the rotation axis and also at large distances above the disc where the classical analytical model fails to provide physically acceptable solutions. The outflow starts with a subslow magnetosonic speed, which subsequently crosses all relevant MHD critical points and corresponding magnetosonic separatrix surfaces. The characteristics are plotted together with the Mach cones and the superfast magnetosonic outflow satisfies MHD causality. The final solution remains close enough to the analytical one, which is thus shown to be topologically stable and robust for various boundary conditions.     

XMM–Newton Reflection Grating Spectrometer observations of the prototypical starburst galaxy M82

We present results from XMM–Newton Reflection Grating Spectrometer observations of the prototypical starburst galaxy M82. These high-resolution spectra represent the best X-ray spectra to date of a starburst galaxy. A complex array of lines from species over a wide range of temperatures is seen, the most prominent being due to Lyman α emission from abundant low- Z elements such as N, O, Ne, Mg and Si. Emission lines from helium-like charge states of the same elements are also seen in emission, as are strong lines from the entire Fe L series. Further, the O vii line complex is resolved and is seen to be consistent with gas in collisional ionization equilibrium.
Spectral fitting indicates emission from a large mass of gas with a differential emission measure over a range of temperatures (from ∼ 0.2 to ∼ 1.6 keV, peaking at ∼ 0.7 keV), and evidence for super-solar abundances of several elements is indicated. Spatial analysis of the data indicates that low-energy emission is more extended to the south and east of the nucleus than to the north and west. Higher energy emission is far more centrally concentrated.     

Near-infrared integral-field spectroscopy of violent starburst environments

Near-infrared (NIR) integral-field spectroscopy (IFS) of violent starburst environments at high spatial (and spectral) resolution has the potential to revolutionise our ideas regarding the local interactions between the newly formed massive stars and the interstellar medium (ISM) of their host galaxies. To illustrate this point, I present NIR IFS analysis of the central starburst region of NGC 1140, obtained with CIRPASS on Gemini-South. While strong [Feii] emission is found throughout the galaxy, higher-order Brackett emission is predominantly associated with the northern starburst region. Based on the spatial distributions of the [Feii] versus Brackett line emission, I conclude that a galaxy-wide starburst was induced several ×10⁷ yr ago, with more recent starburst activity concentrated around the northern starburst region. I look forward and discuss the exciting prospects that IFS at higher spatial (and spectral) resolution will allow us trace (i) the massive outflows (“superwinds”) expected to originate in the dense, young massive star clusters commonly found in intense starburst environments, and (ii) their impact on the galaxy’s ISM.     

AMRVAC and relativistic hydrodynamic simulations for gamma-ray burst afterglow phases

We apply a novel adaptive mesh refinement (AMR) code, AMRVAC (Adaptive Mesh Refinement version of the Versatile Advection Code), to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold interstellar medium (ISM). We do this for both 1D isotropic and full 2D jet-like fireball models. This is relevant for gamma-ray bursts (GRBs), and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell–ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor  γ= 100  up to the almost Newtonian     phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jet-like GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.     

Hα emission line morphologies in Markarian starburst galaxies

We present broad bandR and narrow band Hα emission line images of a sample of optically selected starburst galaxies from the Markarian lists. The emission line morphology is studied and global properties like luminosities, equivalent widths and star formation rates are derived. The radial distribution of Ha flux and the EW are determined using concentric aperture photometry on the emission line and the continuum images. Ha flux is generally found to peak in the nuclear region and fall off outwards. The EW is found to peak off-center in most of the cases implying that though the intensity of emission is maximum at the nucleus, the star formation activity relative to the underlying continuum often peaks away from the center in Markarian starburst galaxies.     

GRMHD/RMHD simulations & stability of magnetized spine-sheath relativistic jets

A new general relativistic magnetohydrodynamics (GRMHD) code “RAISHIN” used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, c/2, on the KH instability associated with a relativistic, γ=2.5, jet spine-sheath interaction. In the simulations sound speeds up to and Alfvén wave speeds up to ∼0.56c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from c/2 sheath speeds is found.     

Hydrodynamic simulations of rotating molecular jets

Molecular outflows and the jets which may drive them can be expected to display signatures associated with rotation if they are the channels through which angular momentum is extracted from material accreting on to protostars. Here, we determine some basic signatures of rapidly rotating flows through three-dimensional numerical simulations of hydrodynamic jets with molecular cooling and chemistry. We find that these rotating jets generate a broad advancing interface which is unstable and develops into a large swarm of small bow features. In comparison to precessing jets, there is no stagnation point along the axis. The greater the rotation rate, the greater the instability. On the other hand, velocity signatures are only significant close to the jet inlet since jet expansion rapidly reduces the rotation speed. We present predictions for atomic, H₂ and CO submillimetre images and spectroscopy including velocity channel maps and position–velocity diagrams. We also include simulated images corresponding to Spitzer IRAC band images and CO emission, relevant for APEX and eventual ALMA observations. We conclude that protostellar jets often show signs of slow precession but only a few sources display properties which could indicate jet rotation.     

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.     

Time-dependent simulations of steady C-type shocks

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging, and studied the effect of single- and multiple-sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionization structures including appropriate microphysics. At low densities, the grains do not play any significant rôle in the shock dynamics. At high densities, the ionization fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allow a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena.