In this paper we present the updated empirical radio surface-brightness-to-diameter ($\Sigma$--$D$) relation for Galactic supernova remnants (SNRs) calibrated using $110$ SNRs with reliable distances. We apply orthogonal fitting procedure and kernel density smoothing in $\Sigma-D$ plane and compare the results with the latest theoretical $\Sigma-D$ relations derived from simulations of radio evolution of SNRs. We argue that the best agreement between the empirical and simulated $\Sigma-D$ relations is achieved if the mixed-morphology SNRs and SNRs of both, low brightness and small diameter, are filtered out from the calibration sample. The distances to $5$ newly discovered remnants and $27$ new candidates for shell SNRs are estimated from our full and filtered calibration samples.

We present 1 to 10GHz radio continuum flux density, spectral index, polarisation and Rotation Measure (RM) images of the youngest known Galactic Supernova Remnant (SNR) G1.9+0.3, using observations from the Australia Telescope Compact Array (ATCA). We have conducted an expansion study spanning 8 epochs between 1984 and 2017, yielding results consistent with previous expansion studies of G1.9+0.3. We find a mean radio continuum expansion rate of ($0.78 \pm 0.09$) per cent year$^{-1}$ (or $\sim8900$ km s$^{-1}$ at an assumed distance of 8.5 kpc), although the expansion rate varies across the SNR perimeter. In the case of the most recent epoch between 2016 and 2017, we observe faster-than-expected expansion of the northern region. We find a global spectral index for G1.9+0.3 of $-0.81\pm0.02$ (76 MHz$-$10 GHz). Towards the northern region, however, the radio spectrum is observed to steepen significantly ($\sim -$1). Towards the two so called (east & west) "ears" of G1.9+0.3, we find very different RM values of 400-600 rad m$^{2}$ and 100-200 rad m$^{2}$ respectively. The fractional polarisation of the radio continuum emission reaches (19 $\pm$ 2)~per~cent, consistent with other, slightly older, SNRs such as Cas~A.

This review presents the fundamentals of the particle acceleration processes active in interstellar medium (ISM), which are essentially based on the so-called Fermi mechanism theory. More specifically, the review presents here in more details the first order Fermi acceleration process -- also known as diffusive shock acceleration (DSA) mechanism. In this case, acceleration is induced by the interstellar (IS) shock waves. These IS shocks are mainly associated with emission nebulae (HII regions, planetary nebulae and supernova remnants). Among all types of emission nebulae, the strongest shocks are associated with supernova remnants (SNRs). Due to this fact they also provide the most efficient manner to accelerate ISM particles to become {high energy particles}, i.e.~cosmic-rays (CRs). The review therefore focuses on the particle acceleration at the strong shock waves of supernova remnants.

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies.

In this paper we present new optical observations of the galaxy NGC 185 intended to reveal the status of supernova remnants (SNRs) in this dwarf companion of the Andromeda galaxy. Previously, it was reported that this galaxy hosts one SNR. Our deep photometric study with the 2m telescope at Rozhen National Astronomical Observatory using narrow-band H$\alpha$ and [SII] filters revealed complex structure of the interstellar medium in the center of the galaxy. To confirm the classification and to study the kinematics of the detected nebulae, we carried out spectroscopic observations using the SCORPIO multi-mode spectrograph at the 6m telescope at the Special Astrophysical Observatory of the Russian Academy of Science, both in low- and high-resolution modes. We also searched the archival X-ray and radio data for counterparts of the candidate SNRs identified by our optical observations. Our observations imply the presence of one more SNR, one possible HII region previously cataloged as part of an SNR, and the presence of an additional source of shock ionization in one low-brightness PN. We detected enhanced [SII]/H_alpha and [NII]/H_alpha line ratios, as well as relatively high (up to 90 km s$^{-1}$) expansion velocities of the two observed nebulae, motivating their classification as SNRs (with diameters of 45 pc and 50 pc), confirmed by both photometric and spectral observations. The estimated electron density of emission nebulae is 30 - 200 cm$^{-3}$. Archival XMM-Newton observations indicate the presence of an extended, low-brightness, soft source in projection of one of the optical SNRs, whereas the archival VLA radio image shows weak, unresolved emission in the center of NGC 185.

In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$.

A widely accepted paradigm is that equipartition (eqp) between the energy density of cosmic rays (CRs) and the energy density of the magnetic field cannot be sustained in supernova remnants (SNRs). However, our 3D hydrodynamic supercomputer simulations, coupled with a non-linear diffusive shock acceleration model (NLDSA), provide evidence that eqp may be established at the end of the Sedov phase of evolution in which most SNRs spend the longest portions of their lives. We introduce term "constant partition" for any constant ratio between the CR energy density and the energy density of magnetic field in an SNR, while term "equipartition" should be reserved for the case of approximately the same values of the energy density (also, it is {constant partition} in the order of magnitude) of ultrarelativistic electrons only (or CRs in total) and the energy density of magnetic field. Our simulations suggest that this approximate constant partition does exist in all but the youngest SNRs. We speculate that since evolved SNRs at the end of the Sedov phase of evolution can reach eqp between CRs and magnetic fields, they may be responsible for initializing this type of eqp in the interstellar medium. Additionally, we show that eqp between the electron component of CRs and magnetic field may be used for calculating the magnetic field strength directly from observations of synchrotron emission from SNRs. The values of magnetic field strengths in SNRs given here are approximately 2.5 times lower than values calculated by Arbutina et al. (2012, 2013).

We present a model for the radio evolution of supernova remnants (SNRs) obtained by using three-dimensional (3D) hydrodynamic simulations, coupled with nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. We model the radio evolution of SNRs on a global level, by performing simulations for wide range of the relevant physical parameters, such as the ambient density, the supernova (SN) explosion energy, the acceleration efficiency and the magnetic field amplification (MFA) efficiency. We attribute the observed spread of radio surface brightnesses for corresponding SNR diameters to the spread of these parameters. In addition to our simulations of type Ia SNRs, we also considered SNR radio evolution in denser, nonuniform circumstellar environment, modified by the progenitor star wind. These simulations start with the mass of the ejecta substantially higher than in the case of a type Ia SN and presumably lower shock speed. The magnetic field is understandably seen as very important for the radio evolution of SNRs. In terms of MFA, we include both resonant and non-resonant modes in our large scale simulations, by implementing models obtained from first-principles, particle-in-cell (PIC) simulations and non-linear magnetohydrodynamical (MHD) simulations. We test the quality and reliability of our models on a sample consisting of Galactic and extragalactic SNRs. Our simulations give $\Sigma-D$ slopes between -4 and -6 for the full Sedov regime. Recent empirical slopes obtained for the Galactic samples are around -5, while for the extragalactic samples are around -4.

We analyze the influence of fractal structure of the interstellar matter (ISM) density on the parameter values for the radio surface brightness to diameter ($\Sigma-D$) relation for supernovae remnants (SNRs). We model a dense ISM as a molecular cloud with fractal density structure. SNRs are modelled as spheres of different radius scattered in the modelled ISM. The surface brightness of the SNRs is calculated from the simple relation $\Sigma \propto \rho^{0.5}D^{-3.5}$ and also from the parametrized more general form $\Sigma \propto \rho^{\eta}D^{-\beta_0}$. Our results demonstrate that empirical $\Sigma-D$ slopes that are steeper than the ones derived from theory, might be partly explained with the fractal structure of the ambient medium into which SNRs expand. The slope of the $\Sigma-D$ relation steepens if the density of the regions where SNRs are formed is higher. The simple geometrical effects combined with the fractal structure of the ISM can contribute to a steeper empirical $\Sigma-D$ slopes, especially for older remnants, and this is more pronounced if $\Sigma$ has a stronger dependence on ambient density.

Integrated radio-spectrum of Cas A in continuum was analyzed with special emphasis on possible high frequency spectral curvature. We conclude that the most probable scenario is that Planck's new data reveal the imprint of non-linear particle acceleration in the case of this young Galactic supernova remnant (SNR).

We present the analysis of supernova remnants (SNRs) in the Large Magellanic Cloud (LMC) and their influence on the environment at far-infrared (FIR) and submillimeter wavelengths. We use new observations obtained with the {\it Herschel} Space Observatory and archival data obtained with the {\it Spitzer} Space Telescope, to make the first FIR atlas of these objects. The SNRs are not clearly discernible at FIR wavelengths, however their influence becomes apparent in maps of dust mass and dust temperature, which we constructed by fitting a modified black-body to the observed spectral energy distribution in each sightline. Most of the dust that is seen is pre-existing interstellar dust in which SNRs leave imprints. The temperature maps clearly reveal SNRs heating surrounding dust, while the mass maps indicate the removal of 3.7$^{+7.5}_{-2.5}$ M$_{\odot}$ of dust per SNR. This agrees with the calculations by others that significant amounts of dust are sputtered by SNRs. Under the assumption that dust is sputtered and not merely pushed away, we estimate a dust destruction rate in the LMC of $0.037^{+0.075}_{-0.025}$ M$_\odot$ yr$^{-1}$ due to SNRs, yielding an average lifetime for interstellar dust of $2^{+4.0}_{-1.3}\times10^7$ yr. We conclude that sputtering of dust by SNRs may be an important ingredient in models of galactic evolution, that supernovae may destroy more dust than they produce, and that they therefore may not be net producers of long lived dust in galaxies.

We present updated empirical radio surface-brightness-to-diameter ($\Sigma-D$) relation for supernova remnants (SNRs) in our Galaxy. Our original calibration sample of Galactic SNRs with independently determined distances (Pavlovic et al. 2013, hereafter Paper I) is reconsidered and updated with data which became available in the past two years. The orthogonal fitting procedure and probability-density-function-based (PDF) method are applied to the calibration sample in the $\log \Sigma - \log D$ plane. Orthogonal regression keeps $\Sigma-D$ and $D-\Sigma$ relations invariant. Our previous Monte Carlo simulations verified that the slopes of the empirical $\Sigma-D$ relation should be determined by using orthogonal regression. Updated calibration sample contains 65 shell SNRs. 6 new Galactic SNRs are added to the sample from Paper I, one is omitted and distances are changed for 10 SNRs. The slope derived here is slightly steeper ($\beta \approx 5.2$) than $\Sigma-D$ slope in Paper I ($\beta \approx 4.8$). The PDF method relies on data points density maps which can provide more reliable calibrations that preserve more information contained in the calibration sample. We estimate distances to five new faint Galactic SNRs discovered for the first time by Canadian Galactic Plane Survey and obtained distances of 2.3, 4.0, 1.3, 2.9 and 4.7 kiloparsecs for G108.5+11.0, G128.5+2.6, G149.5+3.2, G150.8+3.8 and G160.1$-$1.1, respectively. The updated empirical relation is used to estimate distances of 160 shell Galactic SNRs and new results change their distance scales up to 15 per cent, compared to results from Paper I. The PDF calculation can provide even few times higher or lower values in comparison with orthogonal fit, as it uses totally different approach. However, in average, this difference is 32, 24 and 18 per cent for mode, median and mean distances.

The theoretical fundamentals of formation of the supernova remnant (SNR) continuum radio spectra are presented in this review. Mainly based on the Fermi 1 theory - also known as diffuse shock acceleration (DSA) - the different shapes (linear or curved in log-log scale) of SNR radio spectra are predicted for both young and evolved SNRs. On the other hand, some particular forms of spectra of older SNRs can be predicted by including the additional processes such as Fermi 2 acceleration mechanism or thermal bremsstrahlung radiation. Also, all of these theoretically predicted forms of radio spectra are compared with real spectra obtained from observations. Finally this review can represent some kind of "atlas" with initial patterns for the different kinds of SNR radio spectra - it should be helpful for radio astronomers in their interpretation of the observed radio spectra.

We present a detailed radio, X-ray and optical study of a newly discovered Large Magellanic Cloud (LMC) supernova remnant (SNR) which we denote MCSNR J0508-6902. Observations from the Australian Telescope Compact Array (ATCA) and the $\textit{XMM-Newton}$ X-ray observatory are complemented by deep H$\alpha$ images and Anglo Australian Telescope AAOmega spectroscopic data to study the SNR shell and its shock-ionisation. Archival data at other wavelengths are also examined. The remnant follows a filled-in shell type morphology in the radio-continuum and has a size of $\sim$74 pc $\times$ 57 pc at the LMC distance. The X-ray emission exhibits a faint soft shell morphology with Fe-rich gas in its interior $-$ indicative of a Type Ia origin. The remnant appears to be mostly dissipated at higher radio-continuum frequencies leaving only the south-eastern limb fully detectable while in the optical it is the western side of the SNR shell that is clearly detected. The best-fit temperature to the shell X-ray emission ($kT = 0.41^{+0.05}_{-0.06}$ keV) is consistent with other large LMC SNRs. We determined an O/Fe ratio of $<21$ and an Fe mass of 0.5-1.8$~M_{\odot}$ in the interior of the remnant, both of which are consistent with the Type Ia scenario. We find an equipartition magnetic field for the remnant of $\sim$28 $\mu$G, a value typical of older SNRs and consistent with other analyses which also infer an older remnant.

We present a detailed study of Australia Telescope Compact Array (ATCA) observations of a newly discovered Large Magellanic Cloud (LMC) supernova remnant (SNR), SNR J0533-7202. This object follows a horseshoe morphology, with a size 37 pc x 28 pc (1-pc uncertainty in each direction). It exhibits a radio spectrum with the intrinsic synchrotron spectral index of alpha= -0.47+-0.06 between 73 and 6 cm. We report detections of regions showing moderately high fractional polarisation at 6 cm, with a peak value of 36+-6% and a mean fractional polarisation of 12+-7%. We also estimate an average rotation measure across the remnant of -591 rad m^-2. The current lack of deep X-ray observation precludes any conclusion about high-energy emission from the remnant. The association with an old stellar population favours a thermonuclear supernova origin of the remnant.

We report on new Australian Telescope Compact Array (ATCA) observations of SNR J0519-6902. The Supernova Remnant (SNR) is small in size (~8 pc) and exhibits a typical SNR spectrum of alpha = -0.53 +- 0.07, with steeper spectral indices found towards the northern limb of the remnant. SNR J0519-6902 contains a low level of radially orientated polarisation at wavelengths of 3 & 6 cm, which is characteristic of younger SNRs. A fairly strong magnetic field was estimated of ~171 microG. The remnant appears to be the result of a typical Type Ia supernovae, sharing many properties as another small and young Type Ia LMC SNR, J0509-6731.

We use observations of the continuum radio emission at 1420, 820, 408, 34.5 and 22 MHz to estimate the mean brightness temperatures of the HB 21 supernova remnant (SNR) at five frequencies. We also presented mean spectral index of HB 21. The spectra of HB 21 are estimated for mean temperatures versus frequency for 1420, 820, 408, 34.5 and 22 MHz. We also presented $T-T$ plots of three frequency pairs: between 1420--34.5, 1420--22, 34.5--22 MHz. We noticed flatter spectral indices at frequencies below 408 MHz. Probably this is due to the absorption by thermal plasma at low frequencies.

The average brightness temperatures and surface brightnesses at 22 MHz are derived for the four main Galactic radio-continuum loops (Loops I, II, III and IV, hereafter radio-loops). Also the radio-continuum spectra for the radio-loops are presented. Adding the average brightness temperatures at 22 MHz derived here with the average brightness temperatures derived from spectra published previously at 408, 820 and 1420 MHz we obtained clearly non-thermal spectral indices for all radio-loops. Our derived spectral indices are slightly shallower than those measured by previous works.

The estimated brightness of the Cygnus loop supernova remnant (SNR) at 2720, 1420, 820, 408 and 34.5 MHz are presented. The observations of the continuum radio emission are used to calculate the mean brightness temperatures and surface brightnesses of this loop at the five frequencies in wide spectral range, using the method we have previously developed for large radio loops. The spectrum for mean temperatures versus frequency between the five frequencies is estimated and the spectral index of Cygnus loop is also obtained. Also, from our results can be concluded that Cygnus loop evolves in the low density environment and the initial energy of supernova explosion was relatively low. The obtained results confirm non-thermal origin of the Cygnus radio loop and show that our method is applicable to almost all remnants.

In this paper we estimated the temperatures and brightnesses of the Monoceros radio loop at 1420, 820 and 408 MHz. Linear spectrum is estimated for mean temperatures versus frequency between 1420, 820 and 408 MHz. The spectral index of Monoceros loop is also obtained. The brightness temperatures and surface brightnesses of the loop are computed using data taken from radio-continuum surveys at the three frequencies. The spectral index of the loop is also obtained from $T-T$ plots between 1420 - 820, 1420 - 408 and 820 - 408 MHz. The obtained results confirm non-thermal origin of the Monoceros radio loop.

We present an extended analysis of the relation between radio surface brightness and diameter -- the so-called $\Sigma-D$ relation for planetary nebulae (PNe). We revise our previous derivation of the theoretical $\Sigma-D$ relation for the evolution of bremsstrahlung surface brightness in order to include the influence of the fast wind from the central star. Different theoretical forms are derived: $\Sigma \propto D^{-1}$ for the first and second phases of evolution and $\Sigma\propto D^{-3}$ for the final stage of evolution. Also, we analyzed several different Galactic PN samples. All samples are influenced by severe selection effects, but Malmquist bias seems to be less influential here than in the supernova remnant (SNR) samples. We derived empirical $\Sigma-D$ relations for 27 sample sets using 6 updated PN papers from which an additional 21 new sets were extracted. Twenty four of these have a trivial form of $\beta \approx 2$. However, we obtain one empirical $\Sigma-D$ relation that may be useful for determining distances to PNe. This relation is obtained by extracting a recent nearby (< 1 kpc) Galactic PN sample.

It has recently been reported that some of the flux density values for an evolved supernova remnant (SNR) HB 3 (G132.7 + 1.3) are not accurate. In this work we revised an analysis of the possible thermal emission at radio frequencies from the SNR HB 3 using the recently published, corrected, flux density values. A model including a sum of non-thermal (purely synchrotron) and thermal (bremsstrahlung) component is applied for fitting integrated radio spectrum of the SNR. The contribution of thermal component in total volume emissivity at $1 \mathrm{GHz}$ was estimated to be $\approx37 %$. The ambient density was also estimated to be $n\approx 9 \mathrm{cm}^{-3}$ for the $\mathrm{T}=10^{4} \mathrm{K}$. Again, we obtained the relatively significant presence of thermal emission at radio frequencies from the SNR so we could support interaction between SNR HB 3 and adjacent molecular cloud associated with the H {\sc ii} region W3. Our model estimates for thermal component contribution to total volume emissivity at $1 \mathrm{GHz}$ and ambient density are similar to those obtained earlier ($\approx40 %$, $\approx10 \mathrm{cm^{-3}}$). It is clear that the corrected flux density values do not change the basic conclusions.

In this article we use 1420 MHz data to demonstrate the likely reality of Galactic radio Loops V and VI. We further estimate distances and spectral indices for both these and the four main radio loops. In the cases of Loops I - IV, radio spectral indices are calculated from the mean brightnesses at 1420 and 820/404 MHz. The spectral indices of Loops V and VI are obtained from $T - T$ plots between 1420 and 408 MHz. Using the supernova remnant (SNR) hypothesis for the origin of radio loops, distances are calculated from the surface brightnesses and the angular diameters at 1420 MHz. We also study how results for brightnesses and distances of radio loops agree with current theories of SNR evolution. For this purpose, the ambient density and initial explosion energy of the loops are discussed. We also discuss applications of different $\Sigma - D$ relations. The results obtained confirm a non-thermal origin and nearby locations for the Galactic radio loops. Therefore, we have indications that they are very old SNRs that evolve in low ambient densities, with high initial explosion energies.