Energy Conservation in the Thin-Layer Approximation: V. The Surface Brightness in Supernova Remnants

Show more

1. Introduction

The surface brightness versus diameter, (Σ-*D*), for supernova remnants (SNRs) was initially analysed from a theoretical point of view in the framework of the initial conditions for the time evolution of SNRs [1] [2] [3] [4]. Some generic information on Σ-*D* can be found in reviews of SNRs [5] [6].

The astrophysical approach to the Σ-*D* has always mixed the observations with the theory and the statistics. We select some items among others: a catalog of 25 SNRs has been compiled by [7] with the conclusion that SNRs in the galactic halo have diameters greater than those in the disk, the evolutionary properties of SNRs have been deduced from observations with the Molonglo and Parkes radio telescopes [8], the Σ-*D* relationship was used to fix the scale for distances of SNRs [9], the diameters, luminosities, surface brightness, galactic heights for 231 SNRs were processed in the framework of the Sedov solution [10], a new high-latitude SNR [11] was analysed, an updated radio Σ-*D* relationship was derived [12] and the phases of some supernova remnant have been determined from the observations [13]. In order to derive the Σ-*D* relationship, some hypotheses on the single equation of motion should be made, e.g., the conservation of the momentum [14]. Here conversely we will apply the energy conservation in the framework of the thin-layer approximation and we will derive an analytical expression for the Σ-*D* relationship. This paper derives two new equations of motion in the framework of energy conservation for the thin-layer approximation, see Section 2, applies the two analytical expressions for the surface brightness to two galactic samples, see Section 3, and to an extragalactic catalog, see Section 4.

2. The Equations of Motion

The conservation of kinetic energy in spherical coordinates in the framework of the thin-layer approximation, here taken to be an assumption, states that

$\frac{1}{2}{M}_{0}\left({r}_{0}\right){v}_{0}^{2}=\frac{1}{2}M\left(r\right){v}^{2}\mathrm{,}$ (1)

where
${M}_{0}\left({r}_{0}\right)$ and
$M\left(r\right)$ are the swept masses at *r*_{0} and *r*, while *v*_{0} and *v* are the velocities of the thin layer at *r*_{0} and *r*. We now present two equations of motion for SNRs and the back-reaction for one of the two.

2.1. The Inverse Square Law

The medium around the SN is assumed to scale as an inverse square law

$\rho \left(r;{r}_{0}\right)=\{\begin{array}{ll}{\rho}_{c}\hfill & \text{\hspace{0.05em}}\text{if}\text{\hspace{0.17em}}r\le {r}_{0}\hfill \\ {\rho}_{c}{\left(\frac{{r}_{0}}{r}\right)}^{2}\hfill & \text{\hspace{0.05em}}\text{if}\text{\hspace{0.17em}}r>{r}_{0}\hfill \end{array}$ (2)

where ${\rho}_{c}$ is the density at $r=0$ and ${r}_{0}$ is the radius after which the density starts to decrease. When the conservation of energy is applied, the velocity as a function of the radius is

$v\left(r\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\right)=-\frac{\sqrt{-\left(2{r}_{0}-3r\right){r}_{0}}{v}_{0}}{2{r}_{0}-3r}\mathrm{.}$ (3)

The trajectory, *i.e.* the radius as a function of time, is

$r\left(t;{t}_{0},{r}_{0},{v}_{0}\right)=\frac{1}{6}\sqrt[3]{2}\sqrt[3]{{r}_{0}}{\left(\left(9t-9{t}_{0}\right){v}_{0}+2\text{\hspace{0.05em}}{r}_{0}\right)}^{2/3}+\frac{2}{3}{r}_{0},$ (4)

and the velocity as a function of time is

$v\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{v}_{0}\right)=\frac{\sqrt[3]{2}\sqrt[3]{{r}_{0}}{v}_{0}}{\sqrt[3]{\left(9t-9{t}_{0}\right){v}_{0}+2{r}_{0}}}\mathrm{.}$ (5)

More details can be found in [15]. The rate of transfer of mechanical energy, ${L}_{m}$, is

${L}_{m}\left(t\right)=\frac{1}{2}\rho \left(t\right)4\pi r{\left(t\right)}^{2}v{\left(t\right)}^{3}\mathrm{,}$ (6)

where $\rho \left(t\right)$, $r\left(t\right)$ and $v\left(t\right)$ are the instantaneous density, radius and velocity of the SN. We now assume that the density in front of the advancing expansion scales as

$\rho \left(t\right)={\rho}_{0}{\left(\frac{{r}_{0}}{r\left(t\right)}\right)}^{d}\mathrm{,}$ (7)

where *r*_{0} is the radius at *t*_{0} and *d* is a parameter which allows matching the observations. The mechanical luminosity is now

${L}_{m}\left(t\right)=\frac{1}{2}{\rho}_{0}{\left(\frac{{r}_{0}}{r\left(t\right)}\right)}^{d}4\pi r{\left(t\right)}^{2}v{\left(t\right)}^{3}\mathrm{.}$ (8)

In the case here analysed of the inverse square profile for density we have

${L}_{m}\left(t;{v}_{0},{t}_{0},{r}_{0}\right)=\frac{DLM}{81{v}_{0}t-81{v}_{0}{t}_{0}+18{r}_{0}},$ (9)

where

$\begin{array}{c}DLM={\rho}_{0}{6}^{d}{\left({r}_{0}{\left(\sqrt[3]{2}\sqrt[3]{{r}_{0}}{\left(\left(9t-9{t}_{0}\right){v}_{0}+2{r}_{0}\right)}^{\frac{2}{3}}+4{r}_{0}\right)}^{-1}\right)}^{d}\\ \text{\hspace{0.17em}}\times \pi {\left(\sqrt[3]{2}\sqrt[3]{{r}_{0}}{\left(\left(9t-9{t}_{0}\right){v}_{0}+2\text{\hspace{0.05em}}{r}_{0}\right)}^{\frac{2}{3}}+4{r}_{0}\right)}^{2}{r}_{0}{v}_{0}^{3}\mathrm{.}\end{array}$ (10)

We now assume that the observed luminosity, ${L}_{\nu}$, in a given band denoted by the frequency $\nu $ is proportional to the mechanical luminosity

${L}_{\nu}\left(t\right)=cost\ast {L}_{m}\left(t\right)\mathrm{,}$ (11)

where ${L}_{\nu}$ is the observed radio luminosity in a given band and $cost$ a constant which resolves the mismatch between theory and observations. The surface brightness is the luminosity divided by the interested area

$\Sigma =\frac{{L}_{\nu}\left(t\right)}{\pi r{\left(t\right)}^{2}}\mathrm{,}$ (12)

which is

$\Sigma \left(t\mathrm{;}{v}_{0}\mathrm{,}{t}_{0}\mathrm{,}{r}_{0}\right)=cost\frac{4{\rho}_{0}{6}^{d}{\left(\frac{{r}_{0}}{\sqrt[3]{2}\sqrt[3]{{r}_{0}}{\left(\left(9t-9{t}_{0}\right){v}_{0}+2{r}_{0}\right)}^{2/3}+4{r}_{0}}\right)}^{d}{r}_{0}{v}_{0}^{3}}{9{v}_{0}t-9{v}_{0}{t}_{0}+2{r}_{0}}\mathrm{.}$ (13)

2.2. Non-Cubic Dependence

The swept mass is assumed to scale as

$M\left(r;{r}_{0},\delta \right)=\{\begin{array}{ll}{M}_{0}\hfill & \text{\hspace{0.05em}}\text{if}\text{\hspace{0.17em}}r\le {r}_{0}\hfill \\ {M}_{0}{\left(\frac{r}{{r}_{0}}\right)}^{\delta}\hfill & \text{\hspace{0.05em}}\text{if}\text{\hspace{0.17em}}r>{r}_{0}\hfill \end{array}$ (14)

where *M*_{0} is the swept mass at
$r={r}_{0}$, and
$\delta $ is a regulating parameter less than 3. The differential equation of the first order which regulates the motion is obtained by inserting the above
$M\left(r\right)$ in Equation (1)

$\frac{\text{d}r\left(t\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \right)}{\text{d}t}=\frac{{v}_{0}}{\sqrt{{r}^{\delta}{r}_{0}^{-\delta}}}\mathrm{,}$ (15)

which has as the solution

$r\left(t\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \right)=\mathrm{exp}\left(\frac{ER}{\delta +2}\right)$ (16)

where

$\begin{array}{c}ER=\mathrm{ln}({r}_{0}^{\delta +1}\delta t{v}_{0}-{r}_{0}^{\delta +1}\delta {t}_{0}{v}_{0}+2{r}_{0}^{\delta +1}t{v}_{0}-2{r}_{0}^{\delta +1}{t}_{0}{v}_{0}+{r}_{0}^{\delta +2}\\ \text{\hspace{0.17em}}+\frac{{r}_{0}^{\delta}{\delta}^{2}{t}^{2}{v}_{0}^{2}}{4}-\frac{{r}_{0}^{\delta}{\delta}^{2}t{t}_{0}{v}_{0}^{2}}{2}+\frac{{r}_{0}^{\delta}{\delta}^{2}{t}_{0}^{2}{v}_{0}^{2}}{4}+{r}_{0}^{\delta}\delta {t}^{2}{v}_{0}^{2}-2{r}_{0}^{\delta}\delta t{t}_{0}{v}_{0}^{2}\\ \text{\hspace{0.17em}}+{r}_{0}^{\delta}\delta {t}_{0}^{2}{v}_{0}^{2}+{r}_{0}^{\delta}{t}^{2}{v}_{0}^{2}-2{r}_{0}^{\delta}t{t}_{0}{v}_{0}^{2}+{r}_{0}^{\delta}{t}_{0}^{2}{v}_{0}^{2})\mathrm{.}\end{array}$ (17)

The velocity is

$v\left(t\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \right)=\frac{NV}{\left(\delta +2\right){\left(\delta t{v}_{0}-\delta {t}_{0}{v}_{0}+2t{v}_{0}-2{t}_{0}{v}_{0}+2{r}_{0}\right)}^{2}}$ (18)

where

$\begin{array}{c}NV={v}_{0}{\left({r}_{0}^{\delta}{\left({v}_{0}\left(\delta +2\right)\left(t-{t}_{0}\right)+2\text{\hspace{0.05em}}{r}_{0}\right)}^{2}\right)}^{{\left(\delta +2\right)}^{-1}}\times ({24}^{-{\left(\delta +2\right)}^{-1}}{\delta}^{2}t{v}_{0}\begin{array}{c}\text{\hspace{0.05em}}\\ \text{\hspace{0.05em}}\end{array}\\ \text{\hspace{0.17em}}-{24}^{-{\left(\delta +2\right)}^{-1}}{\delta}^{2}{t}_{0}{v}_{0}+{84}^{-{\left(\delta +2\right)}^{-1}}\delta \text{\hspace{0.05em}}t{v}_{0}-{84}^{-{\left(\delta +2\right)}^{-1}}\delta {t}_{0}{v}_{0}\\ \text{\hspace{0.17em}}+{4}^{\frac{\delta +1}{\delta +2}}\delta {r}_{0}+{84}^{-{\left(\delta +2\right)}^{-1}}t{v}_{0}-{84}^{-{\left(\delta +2\right)}^{-1}}{t}_{0}{v}_{0}+{84}^{-{\left(\delta +2\right)}^{-1}}{r}_{0}).\end{array}$ (19)

Figure 1 reports the analytical solution for the NCD case as given by Equation (16) for SN 1993J.

Figure 1. Analytical solution for the NCD case with parameters ${r}_{0}=0.006\text{\hspace{0.17em}}\text{pc}$, ${v}_{0}=10000\text{\hspace{0.17em}}\text{km}/\text{s}$, ${t}_{0}=0.026\text{\hspace{0.17em}}\text{yr}$ and $\delta =0.2$, which gives ${\chi}^{2}=3751$. The astronomical data of SN 1993J are represented with vertical error bars.

The mechanical luminosity is assumed to scale as in Equation (8) and therefore in the NCD case is

$\begin{array}{c}{L}_{m}\left(t;{r}_{0},{v}_{0},{t}_{0},\delta \right)=16\frac{{\text{e}}^{EA}{v}_{0}^{3}{\left({r}_{0}{\text{e}}^{EB}\right)}^{d}\pi {\rho}_{0}}{{\left(\delta +2\right)}^{3}{\left({v}_{0}\left(\delta +2\right)\left(t-{t}_{0}\right)+2{r}_{0}\right)}^{6}}{S}_{2}^{3{\left(\delta +2\right)}^{-1}}\\ \text{\hspace{0.17em}}\times {\left({S}_{1}{4}^{\frac{\delta +1}{\delta +2}}+{4}^{-{\left(\delta +2\right)}^{-1}}{v}_{0}{\delta}^{2}\left(t-{t}_{0}\right)\right)}^{3},\end{array}$ (20)

where

$\begin{array}{c}EA=\frac{1}{\delta +2}(-4\mathrm{ln}\left(2\right)+2\text{\hspace{0.05em}}\mathrm{ln}(4{v}_{0}\left(\delta +2\right)\left(t-{t}_{0}\right){r}_{0}^{\delta +1}+4{r}_{0}^{\delta +2}\\ \text{\hspace{0.17em}}+\left(\left({t}^{2}-2t{t}_{0}+{t}_{0}^{2}\right){\delta}^{2}+\left(4{t}^{2}-8t{t}_{0}+4{t}_{0}^{2}\right)\delta +4{t}^{2}-8t{t}_{0}+4{t}_{0}^{2}\right){v}_{0}^{2}{r}_{0}^{\delta})),\end{array}$ (21)

and

$EB=\frac{2\text{\hspace{0.05em}}\mathrm{ln}\left(2\right)-\mathrm{ln}\left({\left({v}_{0}\left(\delta +2\right)\left(t-{t}_{0}\right)+2{r}_{0}\right)}^{2}{r}_{0}^{\delta}\right)}{\delta +2}\mathrm{,}$ (22)

${S}_{1}=\left(\delta +1\right)\left(t-{t}_{0}\right){v}_{0}+\frac{{r}_{0}\left(\delta +2\right)}{2}\mathrm{,}$ (23)

and

${S}_{2}={\left({v}_{0}\left(\delta +2\right)\left(t-{t}_{0}\right)+2{r}_{0}\right)}^{2}{r}_{0}^{\delta}\mathrm{.}$ (24)

The surface brightness is derived according to Equation (12) and in the NCD case is

$\begin{array}{l}\Sigma \left(t\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}{t}_{0}\mathrm{,}\delta \right)\\ =cost\frac{16\text{\hspace{0.05em}}{\rho}_{0}{\left({4}^{-{\left(\delta +2\right)}^{-1}}{v}_{0}{\delta}^{2}t-{4}^{-{\left(\delta +2\right)}^{-1}}{v}_{0}{\delta}^{2}{t}_{0}+SC{4}^{\frac{\delta +1}{\delta +2}}\right)}^{3}{\text{e}}^{SA}{v}_{0}^{3}{r}_{0}^{\frac{\left(d+3\right)\delta +2d}{\delta +2}}S{B}^{\frac{-6-6\delta}{\delta +2}}}{{\left(\delta +2\right)}^{3}}\mathrm{,}\end{array}$ (25)

where

$SA=-\frac{d\left(-2\mathrm{ln}\left(2\right)+2\mathrm{ln}\left({v}_{0}\left(\delta +2\right)\left(t-{t}_{0}\right)+2{r}_{0}\right)+\delta \mathrm{ln}\left({r}_{0}\right)\right)}{\delta +2}\mathrm{,}$ (26)

$SB={v}_{0}\left(\delta +2\right)\left(t-{t}_{0}\right)+2{r}_{0}$ (27)

and

$SC=\left(\delta +1\right)\left(t-{t}_{0}\right){v}_{0}+\frac{{r}_{0}\left(\delta +2\right)}{2}\mathrm{.}$ (28)

A comparison of the two models here implemented for the trajectory is reported in Figure 2.

2.3. Non-Cubic Dependence and Back Reaction

The radiative losses per unit length are assumed to be proportional to the flux of momentum as an assumption

Figure 2. Analytical solution for the inverse square model (dashed line) when ${r}_{0}=1\text{\hspace{0.17em}}\text{pc}$, ${v}_{0}=4000\text{\hspace{0.17em}}\text{km}/\text{s}$, ${t}_{0}=10\text{\hspace{0.17em}}\text{yr}$, and NCD model (full line) with the same parameters and $\delta =1.3$.

$-\epsilon {\rho}_{s}{v}^{2}4\pi {r}^{2}\mathrm{,}$ (29)

where $\epsilon $ is a constant and ${\rho}_{s}$ is the density in the thin advancing layer.

The volume, *V*, of the advancing layer is

$V=4\pi {r}^{2}\Delta r$ (30)

with $\Delta r=r/12$, therefore the above density for the advancing layer is

${\rho}_{s}=\frac{M\left(r\mathrm{;}{r}_{0}\mathrm{,}\delta \right)}{V}\mathrm{.}$ (31)

Inserting in the above equation the velocity to the first order as given by Equation (18) the radiative losses, $Q\left(r\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)$, are

$Q\left(r\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)=-12\frac{\epsilon {M}_{0}{v}_{0}^{2}}{r}\mathrm{.}$ (32)

The sum of the radiative losses between *r*_{0} and *r* is given by the following integral, *L*,

$L\left(r\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)={\displaystyle {\int}_{{r}_{0}}^{r}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}Q\left(r\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)\text{d}r=-12\epsilon {M}_{0}{v}_{0}^{2}\mathrm{ln}\left(r\right)+12\epsilon {M}_{0}{v}_{0}^{2}\mathrm{ln}\left({r}_{0}\right)\mathrm{.}$ (33)

The conservation of energy in the presence of the back reaction due to the radiative losses is

$\frac{2\pi \text{\hspace{0.05em}}{r}_{0}^{3}{v}^{2}}{3}{\left(\frac{r}{{r}_{0}}\right)}^{\delta}+16\epsilon \text{\hspace{0.05em}}\pi \text{\hspace{0.05em}}{r}_{0}^{3}{v}_{0}^{2}\mathrm{ln}\left(r\right)-16\text{\hspace{0.05em}}\epsilon \text{\hspace{0.05em}}\pi \text{\hspace{0.05em}}{r}_{0}^{3}{v}_{0}^{2}\mathrm{ln}\left({r}_{0}\right)=\frac{2\pi \text{\hspace{0.05em}}{r}_{0}^{3}{v}_{0}^{2}}{3}\mathrm{.}$ (34)

An analytical solution for the velocity to second order, ${v}_{c}\left(r\mathrm{;}{r}_{0}\mathrm{,}{c}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)$, is

${v}_{c}\left(r\mathrm{;}{r}_{0}\mathrm{,}{v}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)={r}^{-\frac{\delta}{2}}{r}_{0}^{\frac{\delta}{2}}\sqrt{24\mathrm{ln}\left({r}_{0}\right)\epsilon -24\mathrm{ln}\left(r\right)\epsilon +1}{v}_{0}\mathrm{.}$ (35)

The inclusion of the back reaction allows the evaluation of the SRS’s maximum length, ${r}_{back}\left({r}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)$, which can be derived by setting the above velocity equal to zero.

${r}_{back}\left({r}_{0}\mathrm{,}\delta \mathrm{,}\epsilon \right)={\text{e}}^{\frac{24\mathrm{ln}\left({r}_{0}\right)\epsilon +1}{24\epsilon}}\mathrm{.}$ (36)

Figure 3 reports the finite radius of the advancing SNR as a function of $\epsilon $.

3. Galactic Application

In the following we will process a sample of data
${x}_{i}\mathrm{,}{y}_{i}$ with *i* varying between 1 and *N* by a power law fit of the type

$y\left(x\right)=C{x}^{\alpha}\mathrm{,}$ (37)

where *C* and
$\alpha $ are two constants to be found from the sample.

The *first* source for the observed Σ-*D* relationship for the SNRs of our galaxy can be found in Figure 3 of [10] which is now digitized, see Figure 4.

Table 1 reports the minimum diameter, *D*_{min}, the average diameter,
$\stackrel{\xaf}{D}$, the maximum diameter, *D*_{max}, the minimum Σ-*D*, Σ-*D*_{min}, the average Σ-*D*,
$\stackrel{\xaf}{\Sigma \text{-}D}$ and the maximum Σ-*D*, Σ-*D*_{max} as well as the two parameters of the power law fit.

A *second* source for the Σ-D relationship is the Green’s catalog [16] where the flux density in Jy at 1 GHz, *S*_{1}, and the mayor and minor angular size in arcmin,
$\theta $, are reported for 295 SNRs. The surface-brightness in SI is

Figure 3. Length of the SNR, ${r}_{back}$, when ${r}_{0}=0.6\text{\hspace{0.17em}}\text{pc}$ as function of $\epsilon $.

Figure 4. Observed Σ-*D* relationship (empty stars) and fitted power law (full line) with parameters as in Table 1.

Table 1. Statistics of the observed Σ-*D* galactic relationship and the two parameters of the power law fit. The theoretical parameters of the inverse square solution are
${r}_{0}=1\text{\hspace{0.17em}}\text{pc}$,
${v}_{0}=4000\text{\hspace{0.17em}}\text{km}/\text{s}$,
$d=1.1$,
${t}_{0}=10\text{\hspace{0.17em}}\text{yr}$,
${t}_{\mathrm{min}}=2{t}_{0}$ and
${t}_{\mathrm{max}}=2.05\times {10}^{5}\text{\hspace{0.05em}}\text{yr}$.

$\Sigma =1.181\times {10}^{-19}\frac{{S}_{1}}{\theta}\frac{\text{W}}{{\text{m}}^{\text{2}}\cdot \text{Hz}\cdot \text{sr}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{SI}\mathrm{,}$ (38)

but after [8] the radio astronomers use the following conversion

$\Sigma =1.505\times {10}^{-19}\frac{{S}_{1}}{\theta}\frac{\text{W}}{{\text{m}}^{\text{2}}\cdot \text{Hz}\cdot \text{sr}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{astronomy},$ (39)

which will also be adopted here.

The probability density function (PDF),
$p\left(z\right)$, to have an SNR as a function of the galactic height *z* is characterized by an exponential PDF

$p\left(z\right)=\frac{1}{b}\mathrm{exp}\text{\hspace{0.05em}}\text{\hspace{0.05em}}-\frac{z}{b},$ (40)

with
$b=83\text{\hspace{0.17em}}\text{pc}$ [10]. The linear diameter (*D*) of an SNR increases with the galactic height see Figure 2 in [10] and in the framework of the model with an inverse square law model for density the parameter *r*_{0} is chosen to have the following dependence with the galactic height

${r}_{0}=0.01+0.5\left(\frac{z}{{z}_{\mathrm{max}}}\right),$ (41)

where ${z}_{\mathrm{max}}=687$ is the maximum galactic height which belongs to an SNR. The above assumption coupled with an initial velocity for the inverse square law model for density of ${v}_{0}=4000\text{\hspace{0.17em}}\text{km}/\text{s}$ for all the SNRs will ensure a longer radius for SNRs with higher galactic heights, see Figure 5.

We now simulate the galactic Σ-*D* relationship for a number of theoretical SNRs, *N*, equal to that observed according to the following rules

1) We randomly generate *N* parameters *z* according to the exponential PDF (40).

2) At each random value of *z* we associate an initial parameter *r*_{0} according to the empirical Equation (41).

3) We randomly generate *N* times, *t*, according to the uniform distribution between a minimum value of time, *t*_{min} and a maximum value of time, *t*_{max}.

4) Given the parameters *r*_{0}, *v*_{0}, *t*_{0} and *t* we evaluate the radius according to Equation (4) for an inverse power law profile for density. The diameter is obtained doubling the above result.

5) Σ is now generated according to Equation (13) once the regulating parameter *d* is provided.

The theoretical display of the results is reported in Figure 6, were we matched the three main parameters (
$\stackrel{\xaf}{D}\left(\text{pc}\right)\mathrm{,}\stackrel{\xaf}{\Sigma \text{-}D}\mathrm{,}\alpha $) which for the observations are (37.17, 2.16 × 10^{−}^{19}, −1.32) and for our simulation are (37.4, 2.24 × 10^{−}^{19}, −1.33).

The surface brightness of the Green’s catalog is reported in Figure 7 and our simulation in Figure 8; see also Table 2 for a comparison between the data of the catalog and the simulation.

4. Extragalactic Application

A sample of SNRs in nearby galaxies with diameter in pc and flux at 1.4 GHz in mJy has been collected [17]. The connected catalog is available at http://cdsweb.u-strasbg.fr/. The statistics of the Σ-*D* relationship for this extragalactic sample is reported in Table 3 and displayed in Figure 9.

Figure 5. Theoretical radius as given by the inverse square model, see Equation (2), as function of the galactic height; ${v}_{0}=4000\text{\hspace{0.17em}}\text{km}/\text{s}$, ${t}_{0}=10\text{\hspace{0.17em}}\text{yr}$ and $t=5\times {10}^{4}\text{yr}$.

Figure 6. Theoretical Σ-D relationship (empty stars) and fitted power law (full line) with parameters as in Table 1.

Figure 7. Histogram of the Σ-*D* for 294 SNR as derived from the Green’s catalog [16].

Figure 8. Histogram for the Σ-*D* of our simulation in the framework of the inverse square model.

Figure 9. Observed Σ-*D* relationship for the extragalactic case (empty stars) and fitting power law (full line) with parameters as in Table 3.

The theoretical Σ-*D* is now evaluated in the framework of the NCD model, see Equation (25), with a numerical procedure which is similar to that of the galactic case but with the difference that there is no dependence of *r*_{0} on *z*. The numerical value of *r*_{0} in the extragalactic case is now randomly generated according to the uniform distribution between a minimum value,
${r}_{0,\mathrm{min}}$, and a maximum value,
${r}_{0,\mathrm{max}}$, see Table 3 and Figure 10.

Figure 10. Theoretical Σ-*D* relationship for the extragalactic case in the framework of the NCD model (empty stars) and fitting power law (full line) with parameters as in Table 3.

Table 2. Statistics of the observed [16] and simulated Σ in the framework of the inverse square model.

Table 3. Statistics of the observed Σ-*D* extragalactic relationship and the two parameters of the power law fit. The theoretical parameters for the NCD case are
${r}_{0}=1\text{\hspace{0.17em}}\text{pc}$,
${r}_{0,\mathrm{min}}=2{r}_{0}$,
${r}_{0,\mathrm{max}}=1.4{r}_{0,\mathrm{min}}$,
${v}_{0}=4000\text{\hspace{0.17em}}\text{km}/\text{s}$,
$\delta =1.3$,
$d=1.1$,
${t}_{0}=10\text{\hspace{0.17em}}\text{yr}$,
${t}_{\mathrm{min}}=2{t}_{0}$ and
${t}_{\mathrm{max}}=2.9\times {10}^{4}\text{yr}$.

5. Conclusions

When an analytical law of motion is available, a theoretical Σ-*D* relationship can be derived as a function of time. Here, in the framework of the energy conservation in the thin-layer approximation, we have derived one formula for Σ-*D* when the density of the surrounding ISM scales as an inverse square law, see Equation (13), and another formula for the NCD case, see Equation (25). The two formulae allow simulating:

1) The galactic Σ-*D* relationship as given by the data of [10], see Figure 6 and Table 1.

2) The galactic Σ-*D* relationship as given by a catalog of SNRs [16], see Figure 8 and Table 2.

3) An extragalactic Σ-*D* catalog [17], see Figure 10 and Table 3.

References

[1] Shklovskii, I.S. (1960) Secular Variation of the Flux and Intensity of Radio Emission from Discrete Sources. Soviet Astronomy, 4, 243-249.

[2] van der Laan, H. (1962) Intense Shell Sources of Radio Emission. Monthly Notices of the Royal Astronomical Society, 124, 179-188.

https://doi.org/10.1093/mnras/124.2.179

[3] Gull, S.F. (1973) A Numerical Model of the Structure and Evolution of Young Supernova Remnants. Monthly Notices of the Royal Astronomical Society, 161, 47-69.

https://doi.org/10.1093/mnras/161.1.47

[4] Caswell, J.L. and Lerche, I. (1979) Supernova Remnants—A Generalized Theoretical Approach to the Radio Evolution. Proceedings of the Astronomical Society of Australia, 3, 343-347.

https://doi.org/10.1017/S1323358000026503

[5] Woltjer, L. (1972) Supernova Remnants. Annual Review of Astronomy and Astrophysics, 10, 129-158.

https://doi.org/10.1146/annurev.aa.10.090172.001021

[6] Long, K.S. (2017) Galactic and Extragalactic Samples of Supernova Remnants: How They Are Identified and What They Tell Us. In: Alsabti, A.W. and Murdin, P., Eds., Handbook of Supernovae, Springer, Berlin, 1-36.

[7] Poveda, A. and Woltjer, L. (1968) Supernovae and Supernova Remnants. Astronomical Journal, 73, 65.

https://doi.org/10.1086/110600

[8] Clark, D.H. and Caswell, J.L. (1976) A Study of Galactic Supernova Remnants, Based on Molonglo-Parkes Observational Data. Monthly Notices of the Royal Astronomical Society, 174, 267-305.

https://doi.org/10.1093/mnras/174.2.267

[9] Huang, Y.L. and Thaddeus, P. (1985) The Sigma-D Relation for Shell-Like Supernova Remnants. Astrophysical Journal, 295, L13-L16.

[10] Xu, J.W., Zhang, X.Z. and Han, J.L. (2005) Statistics of Galactic Supernova Remnants. Chinese Astronomy and Astrophysics, 5, 165-174.

https://doi.org/10.1088/1009-9271/5/2/007

[11] Kothes, R., Reich, P., Foster, T.J. and Reich, W. (2017) G181.1+9.5, a New High-Latitude Low-Surface Brightness Supernova Remnant. Astronomy & Astrophysics, 597, Article No. A116. (Preprint)

https://doi.org/10.1051/0004-6361/201629848

[12] Vukotic, B., Ciprijanovic, A., Vucetic, M., Onic, D. and Urosevic, D. (2019) Updated Radio Sigma-D Relation for Galactic Supernova Remnants—II. Serbian Astronomical Journal, 199, 23-37.

https://doi.org/10.2298/SAJ1999023V

[13] Hu, Y.P., Zeng, H.A., Fang, J., Hou, J.P. and Xu, J.W. (2019) Theoretical *Σ-D* Relations for Shell-Type Galactic Supernova Remnants. Journal of Astrophysics and Astronomy, 40, Article No. 7.

https://doi.org/10.1007/s12036-019-9574-5

[14] Zaninetti, L. (2012) Analytical and Monte Carlo Results for the Surface-Brightness-Diameter Relationship in Supernova Remnants. Astrophysical Journal, 746, Article No. 56. (Preprint)

https://doi.org/10.1088/0004-637X/746/1/56

[15] Zaninetti, L. (2020) Energy Conservation in the Thin Layer Approximation: I. The Spherical Classic Case for Supernovae Remnants International Journal of Astronomy and Astrophysics, 10, 71-88. (Preprint)

https://doi.org/10.4236/ijaa.2020.102006

[16] Green, D.A. (2019) A Revised Catalogue of 294 Galactic Supernova Remnants. Journal of Astrophysics and Astronomy, 40, Article No. 36. (Preprint )

https://doi.org/10.1007/s12036-019-9601-6

[17] Urosević, D., Pannuti, T.G., Duric, N. and Theodorou, A. (2005) The *Σ-D* Relation for Supernova Remnants in Nearby Galaxies. Astronomy & Astrophysics, 435, 437-447.

https://doi.org/10.1051/0004-6361:20042535