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.

In this white paper, we present the scientific cases for adding narrowband optical filters to the Large Synoptic Survey Telescope (LSST). LSST is currently planning to observe the southern sky in 6 broadband optical filters. Three of the four LSST science themes would benefit from adding narrowband filter observations. We discuss the technical considerations of using narrowband filters with the LSST and lay out the scientific impact that would result on the study of AGB stars, emission line nebula (e.g., supernova remnants and planetary nebulae), photometric redshifts of galaxies, and the determination of stellar parameters.

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 collisionless shock wave is treated as a streaming plasma instability in the interstellar medium (ISM). We assume that in a steady state, this instability propagates through the ISM as a self-driven plasma instability, whereby the parameters of the instability are determined such that it causes density and velocity jumps as well as isotropization of the particle velocities, which then must be in accordance with MHD theory. The process of resonant interaction and particle scattering off such instability and their backscattering to an upstream region is simulated. We find that some ISM particles bounce off a shock and thus become suprathermal and eligible to enter the process of acceleration to cosmic ray energies by the mechanism of diffusive shock acceleration.

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 construct the most complete sample of supernova remnants (SNRs) in any galaxy - the Large Magellanic Cloud (LMC) SNR sample. We study their various properties such as spectral index ($\alpha$), size and surface-brightness. We suggest an association between the spatial distribution, environment density of LMC SNRs and their tendency to be located around supergiant shells. We find evidence that the 16 known type Ia LMC SNRs are expanding in a lower density environment compared to the Core-Collapse (CC) type. The mean diameter of our entire population (74) is 41 pc, which is comparable to nearby galaxies. We didn't find any correlation between the type of SN explosion, ovality or age. The $N(<D)$ relationship of $a={0.96}$ implies that the randomised diameters are readily mimicking such an exponent. The rate of SNe occurring in the LMC is estimated to be $\sim$1 per 200 yr. The mean $\alpha$ of the entire LMC SNR population is $\alpha=-0.52$, which is typical of most SNRs. However, our estimates show a clear flattening of the synchrotron $\alpha$ as the remnants age. As predicted, our CC SNRs sample are significantly brighter radio emitters than the type Ia remnants. We also estimate the $\Sigma - D$ relation for the LMC to have a slope $\sim3.8$ which is comparable with other nearby galaxies. We also find the residency time of electrons in the galaxy ($4.0-14.3$ Myr), implying that SNRs should be the dominant mechanism for the production and acceleration of CRs.

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.

In this paper we present the most up-to-date list of nearby galaxies with optically detected supernova remnants (SNRs). We discuss the contribution of the H{\alpha} flux from the SNRs to the total H{\alpha} flux and its influence on derived star formation rate (SFR) for 18 galaxies in our sample. We found that the contribution of SNRs' flux to the total H{\alpha} flux is 5 $\pm$ 5 per cent. Due to the observational selection effects, the SNRs contamination of SFRs derived herein represents only a lower limit.

We present a method for distance calibration without using standard fitting procedures. Instead we use random resampling to reconstruct the probability density function (PDF) of calibration data points in the fitting plane. The resulting PDF is then used to estimate distance-related properties. The method is applied to samples of radio surface brightness to diameter (\Sigma-D) data for the Galactic supernova remnants (SNRs) and planetary nebulae (PNe), and period-luminosity (PL) data for the Large Magellanic Cloud (LMC) fundamental mode classical Cepheids. We argue that resulting density maps can provide more accurate and more reliable calibrations than those obtained by standard linear fitting procedures. For the selected sample of the Galactic SNRs, the presented PDF method of distance calibration results in a smaller average distance fractional error of up to $\approx 16$ percentage points. Similarly, the fractional error is smaller for up to $\approx 8$ and $\approx 0.5$ percentage points, for the samples of Galactic PNe and LMC Cepheids, respectively. In addition, we provide a PDF-based calibration data for each of the samples.

We present observations of the portion of the nearby spiral galaxy IC342 using narrow band [SII] and H$\alpha$ filters. These observations were carried out in November 2011 with the 2m RCC telescope at Rozhen National Astronomical Observatory in Bulgaria. In this paper we report coordinates, diameters, H$\alpha$ and [SII] fluxes for 203 HII regions detected in two fields of view in IC342 galaxy. The number of detected HII regions is 5 times higher than previously known in these two parts of the galaxy.

By applying a method of virtual quanta we derive formulae for relativistic non-thermal bremsstrahlung radiation from relativistic electrons as well as from protons and heavier particles with power-law momentum distribution $N(p)dp = k p^{-q} dp.$ We show that emission which originates from an electron scattering on an ion, represents the most significant component of relativistic non-thermal bremsstrahlung. Radiation from an ion scattering on electron, known as inverse bremsstrahlung, is shown to be negligible in overall non-thermal bremsstrahlung emission. These results arise from theory refinement, where we introduce the dependence of relativistic kinetic energy of an incident particle, upon the energy of scattered photon. In part, it is also a consequence of a different mass of particles and relativistic effects.

The equipartition or minimum-energy calculation is a well-known procedure for estimating magnetic field strength and total energy in the magnetic field and cosmic ray particles by using only the radio synchrotron emission. In one of our previous papers we have offered a modified equipartition calculation for supernova remnants (SNRs) with spectral indices 0.5<\alpha <1. Here we extend the analysis to SNRs with \alpha =0.5 and \alpha =1.

In this paper we present new empirical radio surface brightness-to-diameter ({\Sigma} - D) relations for supernova remnants (SNRs) in our Galaxy. We also present new theoretical derivations of the {\Sigma} - D relation based on equipartition or on constant ratio between cosmic rays and magnetic field energy. A new calibration sample of 60 Galactic SNRs with independently determined distances is created. Instead of (standard) vertical regression, used in previous papers, different fitting procedures are applied to the calibration sample in the log {\Sigma} - log D plane. Non-standard regressions are used to satisfy the requirement that values of parameters obtained from the fitting of {\Sigma} - D and D - {\Sigma} relations should be invariant within estimated uncertainties. We impose symmetry between {\Sigma} - D and D - {\Sigma} due to the existence of large scatter in both D and {\Sigma}. Using four fitting methods which treat {\Sigma} and D symmetrically, different {\Sigma} - D slopes {\beta} are obtained for the calibration sample. Monte Carlo simulations verify that the slopes of the empirical {\Sigma} - D relation should be determined by using orthogonal regression, because of its good performance for data sets with severe scatter. The slope derived here ({\beta} = 4.8) is significantly steeper than those derived in previous studies. This new slope is closer to the updated theoretically predicted surface brightness-diameter slope in the radio range for the Sedov phase. We also analyze the empirical {\Sigma} - D relations for SNRs in the dense environment of molecular clouds and for SNRs evolving in lower-density interstellar medium. Applying the new empirical relation to estimate distances of Galactic SNRs results in a dramatically changed distance scale.

In this paper, we investigate the possibility of significant production of thermal bremsstrahlung radiation at radio continuum frequencies that could be linked to some Galactic supernova remnants (SNRs). The main targets for this investigation are SNRs expanding in high density environments. There are several indicators of radio thermal bremsstrahlung radiation from SNRs, such as a flattening at higher frequencies and thermal absorption at lower frequencies intrinsic to an SNR. In this work we discuss the radio continuum properties of 3 SNRs that are the best candidates for testing our hypothesis of significant thermal emission. In the case of SNRs IC443 and 3C391, thermal absorption has been previously detected. For IC443, the contribution of thermal emission at 1 GHz, from our model fit is 3-57%. It is similar to the estimate obtained from the thermal absorption properties (10-40% at 1 GHz). In the case of the 3C391 the conclusions are not so clear. The results from our model fit (thermal emission contribution of 10-25% at 1 GHz) and results obtained from the low frequency absorption (thermal contribution of 0.15-7% at 1 GHz) do not overlap. For the SNR 3C396 we suggest that if previously detected thermal absorption could be intrinsic to the SNR then the thermal emission (<47% at 1 GHz from our model fit) could be significant enough to shape the radio continuum spectrum at high frequencies. Polarization observations for these SNRs can constrain the strength of a thermal component. Reliable observations at low frequencies (<100 MHz) are needed as well as more data at high radio frequencies (>1 GHz), in order to make stronger conclusions about the existence of "radio thermally active" SNRs.

Expressions for surface brightness distribution and for flux density have been theoretically derived in the case of two simple models of a shell supernova remnant. The models are: a homogenous optically thin emitting shell with constant emissivity and a synchrotron shell source with radial magnetic field. Interactive Data Language (IDL) codes for fitting theoretically derived emission profiles assuming these two models to mean profiles of shell supernova remnants obtained from radio observations have been written.

We present observations of the HI tidal arm near dwarf galaxy NGC 3077 (member of the M81 galaxy group) in narrow band [SII] and H_alpha filters. Observations were carried out in March 2011 with the 2m RCC telescope at NAO Rozhen, Bulgaria. Our search for possible supernova remnant candidates (identified as sources with enhanced [SII] emission relative to their H_alpha emission) in this region yielded no sources of this kind. Nevertheless, we found a number of objects with significant H_alpha emission that probably represent uncatalogued, low brightness HII regions.

Determination of the magnetic field strength in the interstellar medium is one of the most complex tasks of contemporary astrophysics. We can only estimate the order of magnitude of the magnetic field strength by using a few very limited methods. Besides Zeeman effect and Faraday rotation, the equipartition or the minimum-energy calculation is a widespread method for estimating magnetic field strength and energy contained in the magnetic field and cosmic ray particles by using only the radio synchrotron emission. Despite of its approximate character, it remains a useful tool, especially when there is no other data about the magnetic field in a source. In this paper we give a modified calculation which we think is more appropriate for estimating magnetic field strengths and energetics in supernova remnants (SNRs). Finally, we present calculated estimates of the magnetic field strengths for all Galactic SNRs for which the necessary observational data are available. The web application for calculation of the magnetic field strength of SNRs is available at http://poincare.matf.bg.ac.rs/~arbo/eqp/.

The radio surface brightness-to-diameter ({\Sigma} - D) relation for supernova remnants (SNRs) in the starburst galaxy M82 is analyzed in a statistically more robust manner than in the previous studies that mainly discussed sample quality and related selection effects. The statistics of data fits in log {\Sigma} - log D plane are analyzed by using vertical (standard) and orthogonal regressions. As the parameter values of D - {\Sigma} and {\Sigma} - D fits are invariant within the estimated uncertainties for orthogonal regressions, slopes of the empirical {\Sigma} - D relations should be determined by using the orthogonal regression fitting procedure. Thus obtained {\Sigma} - D relations for samples which are not under severe influence of the selection effects could be used for estimating SNR distances. Using the orthogonal regression fitting procedure {\Sigma} - D slope {\beta} \approx 3.9 is obtained for the sample of 31 SNRs in M82. The results of implemented Monte Carlo simulations show that the sensitivity selection effect does not significantly influence the slope of M82 relation. This relation could be used for estimation of distances to SNRs that evolve in denser interstellar environment, with number denisty up to 1000 particles per cm3 .

We present observations of the nearby tidal dwarf galaxy Holmberg IX in M81 galaxy group in narrow band [SII] and H$\alpha$ filters, carried out in March and November 2008 with the 2m RCC telescope at NAO Rozhen, Bulgaria. Our search for resident supernova remnants (identified as sources with enhanced [SII] emission relative to their H$\alpha$ emission) in this galaxy yielded no sources of this kind, besides M&H 10-11 or HoIX X-1. Nevertheless, we found a number of objects with significant H$\alpha$ emission that probably represent uncatalogued HII regions.

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.

The problem of extinction is the most important issue to be dealt with in the process of obtaining true absolute magnitudes of core-collapse (including stripped-envelope) supernovae (SNe). The plane-parallel model, widely used in the past, was shown not to describe extinction adequately. We try to apply an alternative model which introduces radial dependance of extinction in parent galaxies of supernovae. For calculating extinction in our Galaxy we use two different methods and compare the results obtained. Our analysis is primarily focused on a chosen sample of stripped-envelope SNe (Ib/c) for which we find intrinsic peak absolute magnitude $\mathrm{M}_{\mathrm{B}}^{0}=-17.80\pm 0.43$.

The problem of extinction is the most important issue to be dealt with in the process of obtaining true absolute magnitudes of core-collapse supernovae (SNe). The plane-parallel model which gives absorption dependent on galaxy inclination, widely used in the past, was shown not to describe extinction adequately. We try to apply an alternative model which introduces radial ependence of extinction. A certain trend of dimmer SNe with decreasing radius from the center of a galaxy was found, for a chosen sample of stripped-envelope SNe.

An analysis of the relation between radio surface brightness and diameter, so-called Sigma-D relation, for planetary nebulae (PNe) is presented: i) the theoretical Sigma-D relation for the evolution of bremsstrahlung surface brightness is derived; ii) contrary to the results obtained earlier for the Galactic supernova remnant (SNR) samples, our results show that the updated sample of Galactic PNe does not severely suffer from volume selection effect - Malmquist bias (same as for the extragalactic SNR samples) and; iii) we conclude that the empirical Sigma-D relation for PNe derived in this paper is not useful for valid determination of distances for all observed PNe with unknown distances.

Generally accepted scheme distinguishes two main classes of supernovae (SNe): Ia resulting from the old stellar population (deflagration of a white dwarf in close binary systems), and SNe of type II and Ib/c whose ancestors are young massive stars (died in a core-collapse explosion). Concerning the latter, there are suggestions that the SNe II are connected to early B stars, and SNe Ib/c to isolated O or Wolf-Rayet (W-R) stars. However, little or no effort was made to further separate SNe Ib from Ic. We have used assumed SN rates for different SN types in spiral galaxies in an attempt to perform this task. If isolated progenitor hypothesis is correct, our analysis indicates that SNe Ib result from stars of main-sequence mass $23 \mathcal{M}_{\odot} \lesssim \mathcal{M} \lesssim 30 \mathcal{M}_{\odot}$, while the progenitors of SNe Ic are more massive stars with $\mathcal{M} \gtrsim 30 \mathcal{M}_{\odot}$. Alternatively, if the majority of SNe Ib/c appear in close binary systems (CBs) then they would result from the same progenitor population as most of the SNe II, i.e. early B stars with initial masses of order $\mathcal{M} \sim 10 \mathcal{M}_{\odot}$. Future observations of SNe at high-redshift ($z$) and their rate will provide us with unique information on SN progenitors and star-formation history of galaxies. At higher-$z$ (deeper in the cosmic past) we expect to see the lack of type Ia events, i.e. the dominance of core-collapse SNe. Better understanding of the stripped-envelope SNe (Ib/c), and their potential use as distance indicators at high-$z$, would therefore be of great practical importance.

In this paper we apply and discuss a method for the determination of the magnetic field ($H$) evolution in supernova remnants (SNRs) from radio luminosity at given frequency $\nu$ to diameter ($L_\mathrm{\nu}-D$) correlation. We assumed that $H$ evolves as $H\propto D^{-\delta}$, where $D$ is the diameter of the remnant. Value $\delta\approx1.2$ is obtained under the equipartition assumption from the equations for revised equipartition calculation (REC) and by using the data sample from the nearby starburst galaxy M82. We try to investigate whether or not SNRs in M82 are in the equipartition state. This is done by comparison of our empirically obtained $\delta$ with the theoretical value expected for the equipartition conditions. Inconsistence between the value obtained for equipartition conditions and adiabatic expansion $(\delta=1.5)$ and the value empirically obtained herein, can be explained mainly by the influence of sensitivity selection effects which tend to flatten the slope of the $L_\nu-D$ relations for extragalactic samples.