Energy Conservation in the thin layer approximation: V. The surface brightness in supernova remnants

Two new equations of motion for a supernova remnant (SNR) are derived in the framework of energy conservation for the thin-layer approximation. The first one is based on an inverse square law for the surrounding density and the second one on a non-cubic dependence of the swept mass. Under the assumption that the observed radio-flux scales as the flux of kinetic energy, two scaling laws are derived for the temporal evolution of the surface brightness of SNRs. The astrophysical applications cover two galactic samples of surface brightness and an extragalactic one.


Introduction
The surface brightness versus diameter, (Σ − D ), for supernova remnants (SNRs) were 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 catalogue 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.

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 where M 0 (r 0 ) and M (r) 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.

The inverse square law
The medium around the SN is assumed to scale as an inverse square law where ρ 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 The trajectory, i.e. the radius as a function of time, is and the velocity as a function of time is More details can be found in [15]. The rate of transfer of mechanical energy, L m , is where ρ(t),r(t) and v(t) are the instantaneous density, radius and velocity of the SN. We now assume that the density in front of the advancing expansion scales as where r 0 is the radius at t 0 and d is a parameter which allows matching the observations. The mechanical luminosity is now In the case here analysed of the inverse square profile for density we have where DLM = ρ 0 6 d r 0 We now assume that the observed luminosity, L ν , in a given band denoted by the frequency ν is proportional to the mechanical luminosity where L ν 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 which is

Non cubic dependence
The swept mass is assumed to scale as where M 0 is the swept mass at r = r 0 , and δ 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 (r) in equation (1) dr which has as the solution where The velocity is Figure 1 reports the analytical solution for the NCD case as given by equation (16) for SN 1993J . The mechanical luminosity is assumed to scale as in equation (8) and therefore in the NCD case is where and and The surface brightness is derived according to equation (12) and in the NCD case is where and A comparison of the two models here implemented for the trajectory is reported in Figure 2.

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 − ρ s v 2 4 πr 2 , where is a constant and rho s is the density in the thin advancing layer. The volume, V , of the advancing layer is with ∆r = r/12, therefore the above density for the advancing layer is Inserting in the above equation the velocity to the first order as given by equation (18) the radiative losses, Q(r; r 0 , v 0 , δ, ), are The sum of the radiative losses between r 0 and r is given by the following integral, L, The conservation of energy in the presence of the back reaction due to the radiative losses is An analytical solution for the velocity to second order, v c (r; r 0 , c 0 , δ, ), is v c (r; r 0 , v 0 , δ, ) = r − δ 2 r 0 δ 2 24 ln (r 0 ) − 24 ln (r) + 1v 0 .
The inclusion of the back reaction allows the evaluation of the SRS's maximum length, r back (r 0 , δ, ) , which can be derived by setting the above velocity equal to zero.
(36) Figure 3 reports the finite radius of the advancing SNR as a function of .

Galactic Application
In the following we will process a sample of data x i , y i with i varying between 1 and N by a power law fit of the the type where C and α 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, D, the maximum diameter, D max , the minimum Σ − D , Σ − D min , the average Σ − D , Σ − 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, θ, are reported for 295 SNRs. The surface-brightness in SI is   Table 1.
which will also be adopted here.
The probability density function (PDF), p(z), to have an SNR as a function of the galactic height z is characterized by an exponential PDF with b = 83 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 the galactic height where z 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 km 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 (D(pc), Σ − D, α) 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.   Table 1. Table 2. Statistics of the observed [16] and simulated Σ in the framework of the inverse square model.   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. 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 r0 = 1 pc, r0,min = 2 r0, r0,max = 1.4 r0,min, v0 = 4000 km/s, δ = 1.3, d = 1.1, t0 = 10 yr, tmin = 2 t0 and tmax = 2.9 10 4 yr.  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,min , and a maximum value, r 0,max , see Table 3 and Figure 10.

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: 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. 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.