Near-infrared [Fe ii] emission from supernova remnants and the supernova rate of starburst galaxies

In an effort better to calibrate the supernova rate of starburst galaxies as determined from near-infrared [Fe  ii ] features, we report on a [Fe  ii ] λ 1.644 μm line-imaging survey of a sample of 42 optically selected supernova remnants (SNRs) in M33. A wide range of [Fe  ii ] luminosities are observed within our sample (from less than 6 to 695 L_⊙). Our data suggest that the bright [Fe  ii ] SNRs are entering the radiative phase and that the density of the local interstellar medium (ISM) largely controls the amount of [Fe  ii ] emission. We derive the following relation between the [Fe  ii ] λ 1.644 μm line luminosity of radiative SNRs and the electronic density of the post-shock gas, n _e: L _{[Fe  ii ]}     (cm⁻³). We also find a correlation in our data between L _{[Fe  ii ]} and the metallicity of the shock-heated gas, but the physical interpretation of this result remains inconclusive, as our data also show a correlation between the metallicity and n _e. The dramatically higher level of [Fe  ii ] emission from SNRs in the central regions of starburst galaxies is most likely due to their dense environments, although metallicity effects might also be important. The typical [Fe  ii ]-emitting lifetime of a SNR in the central regions of starburst galaxies is found to be of the order of 10⁴ yr. On the basis of these results, we provide a new empirical relation allowing the determination of the current supernova rate of starburst galaxies from their integrated near-infrared [Fe  ii ] luminosity.     

VLBI observations of supernova remnants in Messier 82

We have used the European VLBI Network (EVN) at 18 cm to study five of the more compact radio sources in the starburst galaxy M82. The angular resolution of the observations is 15 mas, corresponding to 0.2 pc at the distance of M82. The observations reveal shells ranging in diameter from 40 to 90 mas (0.6 to 1.4 pc), although the strongest source (41.95+575) is only marginally resolved by these measurements (∼20×10 mas²).
We have found clear evidence for expansion in one of the shell sources (43.31+592) by re-analysing, in wide-field mode, EVN data taken in 1986. Between 1986 and 1997 this source has increased its diameter by 13.6±2 mas, corresponding to an average expansion velocity of 9850±1500 km s⁻¹. If we assume that the remnant is in free expansion, this is consistent with a supernova event in the early 1960s. Hence this remnant is almost certainly younger than the strongest, most compact source (41.95+575) which was known to be present in the 1960s. 41.95+575 shows no clear evidence for expansion (<4000 km s⁻¹), consistent with a greater age; this is further evidence of its anomalous status. Comparison of the EVN images with earlier MERLIN data is also consistent with expansion in at least two more of the sources. We discuss the flux density variability of the compact sources in M82 and conclude that, with the exception of 41.95+575 and two transient sources, there is little evidence for significant changes in flux density of most of the remnants since the early 1980s.     

Neutral hydrogen absorption against supernova remnants in M82

We present MERLIN neutral hydrogen absorption measurements against supernova remnants in the central starburst region of M82 with an angular resolution of ∼ 0.4 arcsec. We detect H  I absorption or set significant upper limits against 33 supernova remnants from which we have been able to deduce column densities. Hence, using these measurements, we are able to probe the neutral hydrogen distribution and dynamics of the interstellar medium in M82 along 33 lines of sight on linear scales of order 1 pc.   Our results show column densities ranging from <1.6 to >30 × 10²¹ atom cm⁻² with the highest values seen towards the edge of the 250-pc 'ring'. The absorption velocities show a gradient of 7.3 ± 4 km s⁻¹ arcsec⁻¹, consistent with rotation parameters of this 'ring' inferred from other measurements. The absorption velocities against individual remnants show deviations of typically 30 km s⁻¹ from simple solid body rotation, and a number show multiple velocity absorption features. Although some of these deviations may be the result of the remnants being embedded at different depths within the neutral gas, the velocities cannot be explained by a simple rotating ring.     

Radio emission at 1.4 GHz from luminous compact galaxies

We analyse the properties of the 1.4 GHz continuum emission in 52 star‐forming compact luminous galaxies (LCGs). The fluxes of the 1.4 GHz thermal (free‐free) component are derived from the extinction‐ and aperture‐corrected fluxes of the Hα emission line. The fraction of the thermal 1.4 GHz emission is shown tobe in the range of 2–52 % with a median value of about 15–17 %, and its distribution is close to a log‐normal one. Both the thermal and non‐thermal radio continuum luminosities are shown to be proportional to the masses of the young stellar population of the LCGs. We discuss the approximation of the ratio of the 1.4 GHz luminosity‐to‐the mass of the young stellar population as a function of starburst age. The derived approximation is shown to be similar but less clearly defined compared to ones obtained previously for LCGs from their Hα and UV continuum luminosities. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Revised equipartition and minimum energy formula for magnetic field strength estimates from radio synchrotron observations

The commonly used classical equipartition or minimum‐energy estimate of total magnetic fields strengths from radio synchrotron intensities is of limited practical use because it is based on the hardly known ratio K of the total energies of cosmic ray protons and electrons and also has inherent problems. We present a revised formula, using the number density ratio K for which we give estimates. For particle acceleration in strong shocks K is about 40 and increases with decreasing shock strength. Our revised estimate for the field strength gives larger values than the classical estimate for flat radio spectra with spectral indices of about 0.5–0.6, but smaller values for steep spectra and total fields stronger than about 10 µG. In very young supernova remnants, for example, the classical estimate may be too large by up to 10×. On the other hand, if energy losses of cosmic ray electrons are important, K increases with particle energy and the equipartition field may be underestimated significantly. Our revised larger equipartition estimates in galaxy clusters and radio lobes are consistent with independent estimates from Faraday rotation measures, while estimates from the ratio between radio synchrotron and X‐ray inverse Compton intensities generally give much weaker fields. This may be explained e.g. by a concentration of the field in filaments. Our revised field strengths may also lead to major revisions of electron lifetimes in jets and radio lobes estimated from the synchrotron break frequency in the radio spectrum. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulations of mixed-morphology supernova remnants with anisotropic thermal conduction

We explore the role of anisotropic thermal conduction on the evolution of supernova remnants (SNRs) through interstellar media with a range of densities via numerical simulations. We find that a remnant expanding in a dense environment can produce centre-bright hard X-ray emission within 20 kyr, and centre-bright soft X-ray emission within 60 kyr of the supernova event. In a more tenuous environment, the appearance of a centre-bright structure in hard X-rays is delayed until about 60 kyr. The soft X-ray emission from such a remnant may not become centre bright during its observable lifetime. This can explain the observations that show that mixed-morphology SNRs preferentially occur close to denser, molecular environments. Remnants expanding into denser environments tend to be smaller, making it easier for thermal conduction to make large changes in the temperatures of their hot gas bubbles. We show that the lower temperatures make it very favourable to use high-stage ions as diagnostics of the hot gas bubbles in SNRs. In particular, the distribution of O  viii transitions from shell bright at early epochs to centre bright at later epochs in the evolution of an SNR expanding in a dense interstellar medium when the physics of thermal conduction is included.     

Saturated magnetic field amplification at supernova shocks

Cosmic ray streaming instabilities at supernova shocks are discussed in the quasi-linear diffusion formalism which takes into account the feedback effect of wave growth on the cosmic ray streaming motion. In particular, the non-resonant instability that leads to magnetic field amplification in the short wavelength regime is considered. The linear growth rate is calculated using kinetic theory for a streaming distribution. We show that the non-resonant instability is actually driven by a compensating current in the background plasma. The non-resonant instability can develop into a non-linear regime generating turbulence. The saturation of the amplified magnetic fields due to particle diffusion in the turbulence is derived analytically. It is shown that the evolution of parallel and perpendicular cosmic ray pressures is predominantly determined by non-resonant diffusion. However, the saturation is determined by resonant diffusion which tends to reduce the streaming motion through pitch angle scattering. The saturated level can exceed the mean background magnetic field.     

Σ–D relation for supernova remnants and its dependence on the density of the interstellar medium

During the last couple of decades of work on the  Σ– D   (radio surface brightness to diameter) relation for supernova remnants (SNRs), it has been generally accepted that no single  Σ– D   relation can be constructed for all SNRs. However, it may still be possible to construct the relations for some classes of SNRs. In our previous paper we analysed  Σ– D   relation(s) for remnants in the dense environments of molecular clouds. The aim of this paper is to examine, in the same context, a class of oxygen-rich SNRs, and to extend the analysis to remnants evolving in lower-density interstellar media, namely Balmer-dominated SNRs. We have obtained good relations with certain similarities to our previous findings – similarities that emphasize, again, the role of ambient density in the evolution of SNRs.     

Modeling the total and polarized emission in evolving galaxies: “Spotty” magnetic structures

Future radio observations with the Square Kilometre Array (SKA) and its precursors will be sensitive to trace spiral galaxies and their magnetic field configurations up to redshift z ≈ 3. We suggest an evolutionary model for the magnetic configuration in star‐forming disk galaxies and simulate the magnetic field distribution, the total and polarized synchrotron emission, and the Faraday rotation measures for disk galaxies at z ≲ 3. Since details of dynamo action in young galaxies are quite uncertain, we model the dynamo action heuristically relying only on well‐established ideas of the form and evolution of magnetic fields produced by the mean‐field dynamo in a thin disk. We assume a small‐scale seed field which is then amplified by the small‐scale turbulent dynamo up to energy equipartition with kinetic energy of turbulence. The large‐scale galactic dynamo starts from seed fields of 100 pc and an averaged regular field strength of 0.02 μG, which then evolves to a “spotty” magnetic field configuration in about 0.8 Gyr with scales of about one kpc and an averaged regular field strength of 0.6 μG. The evolution of these magnetic spots is simulated under the influence of star formation, dynamo action, stretching by differential rotation of the disk, and turbulent diffusion. The evolution of the regular magnetic field in a disk of a spiral galaxy, as well as the expected total intensity, linear polarization and Faraday rotation are simulated in the rest frame of a galaxy at 5GHz and 150 MHz and in the rest frame of the observer at 150 MHz. We present the corresponding maps for several epochs after disk formation. Dynamo theory predicts the generation of large‐scale coherent field patterns (“modes”). The timescale of this process is comparable to that of the galaxy age. Many galaxies are expected not to host fully coherent fields at the present epoch, especially those which suffered from major mergers or interactions with other galaxies. A comparison of our predictions with existing observations of spiral galaxies is given and discussed (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A search for young Galactic supernova remnants

A sample of eight small-diameter radio sources has been selected from the Molonglo Galactic Plane Survey (MGPS) as candidates for young Galactic supernova remnants. The sources have been identified in the IRAS and Midcourse Space Experiment infrared data bases and imaged in the H107α radio recombination line (RRL) using the Australia Telescope Compact Array (ATCA). Seven of the sources display high ratios of infrared-to-radio-continuum flux density and/or detectable RRLs and are almost certainly H ii regions. One source (G282.8−1.2) is identified as a possible new young Galactic supernova remnant, based on its relatively weak infrared emission, steep radio spectrum, and possible X-ray emission. The adopted method for distinguishing thermal and non-thermal Galactic radio sources seems promising and could be fruitfully applied to more than 100 small-diameter sources listed in the MGPS.     

Synchrotron emission from anisotropic disordered magnetic fields

We derive expressions for the total and linearly polarized synchrotron emissivity of an element of plasma containing relativistic particles and disordered magnetic field that has been sheared or compressed along three independent directions. Our treatment follows that given by Matthews & Scheuer in the special case of a power-law electron energy spectrum. We show that the emissivity integrals depend on a single parameter, making it straightforward to generate one-dimensional look-up tables. We also demonstrate that our formulae give identical results to those in the literature in special cases.     

Influence of accelerated particles on the large-scale magnetic field in young supernova remnants

We offer a possible explanation for the observational data on the magnetic-field structure in young supernova remnants (SN 1006, Tycho, Kepler, Cas A) that have been obtained by analyzing the polarizations of electromagnetic radiation in the radio, infrared, and other wavelength ranges. The authors of observational works interpret these data as evidence that the ordered magnetic-field component is predominantly radial, but it can be much smaller in amplitude than the stochastic field component that accounts for the bulk of the total magnetic energy. We calculate the magnetic field in supernova remnants by taking into account the shock compression of the primary field and the generation of a large-scale magnetic field by the particles accelerated at the shock front. The assumption that the field in the supernova remnant is the explosion-compressed primary field near the star is inconsistent with observational data, because the tangential (relative to the shock front) field component perpendicular to the radius must prevail in this case. However, allowing for the generation of an additional magnetic field by the electric current of the particles accelerated by a strong shock front leads us to conclude that the field components parallel to the front are suppressed by accelerated particles by several orders of magnitude. Only the component perpendicular to the front remains. Such a field configuration for uniform injection does not lead to the generation of an additional magnetic field, and, in this sense, it is stable. This explains the data on the radial direction of the ordered field component. As regards the stochastic field component, we show that it is effectively generated by accelerated particles if their injection into acceleration at the shock front is nonuniform along the front. Injection nonuniformity can be caused by upstream density nonuniformities. A relative density nonuniformity of the order of several percent is enough for an observable magnetic field with scales on the order of the density nonuniformity scales to be generated.     

Non-thermal emission from old supernova remnants

We study the non-thermal emission from old shell-type supernova remnants (SNRs) on the frame of a time-dependent model. In this model, the time-dependent non-thermal spectra of both primary electrons and protons as well as secondary electron/positron (e^±) pairs can be calculated numerically by taking into account the evolution of the secondary e^± pairs produced from proton–proton (p–p) interactions as accelerated protons collide with the ambient matter in an SNR. The multiwavelength photon spectrum for a given SNR can be produced through leptonic processes such as electron/positron synchrotron radiation, bremsstrahlung and inverse Compton scattering as well as hadronic interaction. Our results indicate that the non-thermal emission of the secondary e^± pairs is becoming more and more prominent when the SNR ages in the radiative phase because the source of the primary electrons has been cut off and the electron synchrotron energy loss is significant for a radiative SNR, whereas the secondary e^± pairs can be produced continuously for a long time in the phase due to the large energy-loss time for the p–p interaction. We apply the model to two old SNRs, G8.7−0.1 and G23.3−0.3, and the predicted results can explain the observed multiwavelength photon spectra for the two sources.