The Giants Arcs as Modeled by the Superbubbles

Show more

1. Introduction

The giant arcs in the cluster of galaxies start to be observed as narrow-like shape by [1] [2] [3] and the first theoretical explanation was the gravitational lensing, see [4] [5] [6] [7] . The determination of the statistical parameters of the giant arcs has been analyzed in order to derive the cosmological parameters, see [8] , in connection with the ΛCDM cosmology, see [9] , in the framework of the triaxiality and substructure of CDM halo, see [10] , in connection with the Wilkinson Microwave Anisotropy Probe (WMAP) data, see [11] , including the effects of baryon cooling in dark matter N-body simulations, see [12] and in order to derive the photometric properties of 105 giant arcs in the Second Red-Sequence Cluster Survey (RCS-2), see [13] . The gravitational lensing is the most common theoretical explanation; we select some approaches among others: [14] evaluated the mass distribution inside distant clusters; [15] studied the statistics of giant arcs in flat cosmologies with and without a cosmological constant; [16] analyzed how the gravitational lensing influences the surface brightness of giant luminous arcs and [17] used the warm dark matter (WDM) cosmologies to explain the lensing in galaxy clusters. Another theoretical line of research explains the giant arcs as shells originated by the galaxies in the cluster: [18] [19] analyzed the limb-brightened shell model, the gravitational lens model and the echo model; [20] suggested that the Gamma-ray burst (GRB) explosions are the sources of the shells with sizes of many kpc, and [21] suggested a connection between the Einstein ring associated to SDP.81 and the evolution of a superbubble (SB) in the intracluster medium.

This paper analyzes in Section 2 three cosmologies in order to calibrate the transversal distance which allows to convert the arcsec in pc. Section 3 is devoted to the evolution of a SB in the intracluster medium. Section 4 reports the observations of the giant arcs and the first phase of a SB. Section 5 reports the various steps which allow to reproduce the shape of the giant arc A2267 and the multiple arcs visible in the cluster of galaxies. Section 6 is dedicated to theory of the image: analytical formulae explain the hole in the central part of the SBs and numerical results reproduce the details of a giant arc.

2. Adopted Cosmologies

In the following we review three cosmological theories.

2.1. ΛCDM Cosmology

The basic parameters of ΛCDM cosmology are: the Hubble constant,
${H}_{0}$ , expressed in km∙s^{−1}∙Mpc^{−1}, the velocity of light, c, expressed in kms^{−1}, and the three numbers
${\Omega}_{\text{M}}$ ,
${\Omega}_{K}$ , and
${\Omega}_{\Lambda}$ , see [22] for more details. In the case of the Union 2.1 compilation, see [23] , the parameters are
${H}_{0}=69.81\text{\hspace{0.17em}}\text{km}\text{\hspace{0.05em}}\cdot {\text{s}}^{-\text{1}}\cdot {\text{Mpc}}^{-\text{1}}$ ,
${\Omega}_{\text{M}}=0.239$ and
${\Omega}_{\Lambda}=0.651$ . To have the luminosity distance,
${D}_{L}\left(z\mathrm{;}{H}_{0}\mathrm{,}c\mathrm{,}{\Omega}_{\text{M}}\mathrm{,}{\Omega}_{\Lambda}\right)$ , as a function of the redshift only, we apply the minimax rational approximation, which is characterized by two parameters, p and q. The luminosity distance,
${D}_{L\mathrm{,3,2}}$ , when
$p=3$ and
$q=2$

${D}_{L,3,2}=\frac{-7.761-1788.53z-3203.06{z}^{2}-65.8463{z}^{3}}{-0.438025-0.334872z+0.0203996{z}^{2}}\text{Mpc}\text{\hspace{1em}}\text{for}\text{\hspace{0.17em}}\text{\hspace{0.05em}}0.001<z<4.$ (1)

The transversal distance in ΛCDM cosmology, ${D}_{T\mathrm{,3,2}}$ , which corresponds to the angle $\delta $ expressed in arcsec is

${D}_{T\mathrm{,3,2}}=\frac{4.84813\delta \left(2.328+502.067z+113.03{z}^{2}\right)}{0.124085+0.149501z+0.0932928{z}^{2}}\text{pc}$ (2)

2.2. Flat Cosmology

The two parameters of the flat cosmology are
${H}_{0}$ , the Hubble constant expressed in km∙s^{−1}∙Mpc^{−1}, and
${\Omega}_{\text{M}}$ which is

${\Omega}_{\text{M}}=\frac{8\text{\pi}\text{\hspace{0.05em}}G{\rho}_{0}}{3{H}_{0}^{2}}\mathrm{,}$ (3)

where G is the Newtonian gravitational constant and ${\rho}_{0}$ is the mass density at the present time. In the case of m = 2 and n = 2 the minimax rational expression for the luminosity distance, ${d}_{L\mathrm{,}m\mathrm{,2,2}}$ , when ${H}_{0}=70\text{\hspace{0.17em}}\text{km}\cdot {\text{s}}^{-\text{1}}\cdot {\text{Mpc}}^{-\text{1}}$ and ${\Omega}_{\text{M}}=0.277$ , is

${d}_{L\mathrm{,}m\mathrm{,2,2}}=\frac{0.0889+748.555z+5.58311{z}^{2}}{0.175804+0.206041z+0.068685{z}^{2}}\text{Mpc}$ (4)

The transversal distance in flat cosmology, ${D}_{Tf\mathrm{,3,2}}$ , which corresponds to the angle $\delta $ expressed in arcsec is

${D}_{Tf\mathrm{,3,2}}=\frac{4.84813\delta \left(0.0889+748.555z+5.58311{z}^{2}\right)}{0.175804+0.206041z+0.0686854{z}^{2}}\text{pc}$ (5)

2.3. Modified Tired Light

In an Euclidean static framework the modified tired light (MTL) has been introduced in Section 2.2 in [24] . The distance in MTL is

$d=\frac{c}{{H}_{0}}\mathrm{ln}\left(1+z\right)\mathrm{.}$ (6)

The distance modulus in the modified tired light (MTL) is

$m-M=\frac{5}{2}\text{\hspace{0.05em}}\frac{\beta \mathrm{ln}\left(z+1\right)}{\mathrm{ln}\left(10\right)}+5\frac{1}{\mathrm{ln}\left(10\right)}\mathrm{ln}\left(\frac{\mathrm{ln}\left(z+1\right)c}{{H}_{0}}\right)+25$ (7)

Here $\beta $ is a parameter comprised between 1 and 3 which allows to match theory with observations. The number of free parameters in MTL is two: ${H}_{0}$ and $\beta $ . The fit of the distance modulus with the data of Union 2.1 compilation gives $\beta =2.37$ , ${H}_{0}=69.32\pm 0.34$ , see [22] , which means the following distance

$d=4324.761\mathrm{ln}\left(1+z\right){\left(1+z\right)}^{1.185}\text{Mpc}.$ (8)

The transversal distance in MTL, ${d}_{T}$ , which corresponds to the angle $\delta $ expressed in arcsec is

${d}_{T}=20967\delta \mathrm{ln}\left(1+z\right){\left(1+z\right)}^{1.185}\text{pc}$ (9)

We report the angular distance for a fixed delta as function of redshift for the three cosmologies, see Figure 1. The angular distance in flat and ΛCDM cosmology does not increase with z, see [25] , in contrast with the modified tired light.

3. The Motion of a SB

We now summarize the adopted profile of density and the equation of motion for a SB.

3.1. The Profile

The density is assumed to have the following hyperbolic dependence on Z which is the third Cartesian coordinate,

Figure 1. Angular distances in kpc for the three cosmologies here considered when $\delta =38.913456$ .

$\rho \left(Z\mathrm{;}{Z}_{0}\mathrm{,}{\rho}_{0}\right)=\{\begin{array}{ll}{\rho}_{0}\hfill & \text{if}\text{\hspace{0.17em}}z\le {Z}_{0}\text{\hspace{0.05em}}\hfill \\ {\rho}_{0}\frac{{Z}_{0}}{z}\hfill & \text{if}\text{\hspace{0.17em}}z>{Z}_{0}\text{\hspace{0.05em}}\hfill \end{array}$ (10)

where the parameter ${Z}_{0}$ fixes the scale and ${\rho}_{0}$ is the density at $Z={Z}_{0}$ . In spherical coordinates the dependence on the polar angle is

$\rho \left(r\mathrm{;}\theta \mathrm{,}{Z}_{0}\mathrm{,}{\rho}_{0}\right)=\{\begin{array}{ll}{\rho}_{0}\hfill & \text{if}\text{\hspace{0.17em}}\mathrm{cos}\left(\theta \right)\le {Z}_{0}\hfill \\ {\rho}_{0}\frac{{Z}_{0}}{r\mathrm{cos}\left(\theta \right)}\hfill & \text{if}\text{\hspace{0.17em}}r\text{\hspace{0.17em}}\mathrm{cos}\left(\theta \right)>{Z}_{0}\hfill \end{array}$ (11)

Given a solid angle $\Delta \Omega $ the mass ${M}_{0}$ swept in the interval $\left[\mathrm{0,}{r}_{0}\right]$ is

${M}_{0}=\frac{1}{3}\text{\hspace{0.05em}}{\rho}_{0}\text{\hspace{0.05em}}{r}_{0}^{3}\Delta \Omega \text{\hspace{0.17em}}.$ (12)

The total mass swept, $M\left(r\mathrm{;}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}\alpha \mathrm{,}\theta \mathrm{,}{\rho}_{0}\right)$ , in the interval $\left[\mathrm{0,}r\right]$ is

$M\left(r\mathrm{;}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}\alpha \mathrm{,}\theta \mathrm{,}{\rho}_{0}\right)=\left(\frac{1}{3}{\rho}_{0}{r}_{0}^{3}+\frac{1}{2}\frac{{\rho}_{0}{Z}_{0}\left({r}^{2}-{r}_{0}^{2}\right)}{\mathrm{cos}\left(\theta \right)}\right)\Delta \Omega \text{\hspace{0.17em}}\mathrm{.}$ (13)

and its approximate value at high values of r is

$M\left(r\mathrm{;}{Z}_{0}\mathrm{,}\alpha \mathrm{,}\theta \mathrm{,}{\rho}_{0}\right)\approx \frac{1}{2}\frac{{r}^{2}{\rho}_{0}{Z}_{0}}{\mathrm{cos}\left(\theta \right)}\Delta \Omega \text{\hspace{0.17em}}\mathrm{.}$ (14)

The density
${\rho}_{0}$ can be obtained by introducing the number density,
${n}_{0}$ , expressed in particles cm^{−3}, the mass of hydrogen,
${m}_{H}$ , and a multiplicative factor f, which is chosen to be 1.4, see [26] ,

${\rho}_{0}=f{m}_{H}{n}_{0}.$ (15)

The astrophysical version of the total approximate swept mass as given by Equation (14), expressed in solar mass units, ${M}_{\odot}$ , is

$M\left({r}_{pc}\mathrm{;}{Z}_{\mathrm{0,}pc}\mathrm{,}{n}_{0}\mathrm{,}\theta \right)\approx \frac{0.0172\text{\hspace{0.05em}}{n}_{0}{z}_{\mathrm{0,}pc}{r}_{pc}^{2}}{\mathrm{cos}\left(\theta \right)}\text{\hspace{0.05em}}{M}_{\odot}\Delta \Omega \mathrm{,}$ (16)

where ${Z}_{\mathrm{0,}pc}$ , and ${r}_{\mathrm{0,}pc}$ are ${Z}_{0}$ , and r expressed in pc.

3.2. The Equation of Motion

The conservation of the classical momentum in spherical coordinates along the solid angle $\Delta \Omega $ in the framework of the thin layer approximation states that

${M}_{0}\left({r}_{0}\right){v}_{0}=M\left(r\right)v\mathrm{,}$ (17)

where ${M}_{0}\left({r}_{0}\right)$ and $M\left(r\right)$ are the swept masses at ${r}_{0}$ and r, and ${v}_{0}$ and v are the velocities of the thin layer at ${r}_{0}$ and r. This conservation law can be expressed as a differential equation of the first order by inserting $v=\frac{\text{d}r}{\text{d}t}$ :

$M\left(r\right)\frac{\text{d}r}{\text{d}t}-{M}_{0}\text{\hspace{0.05em}}{v}_{0}=0.$ (18)

The velocity as a function of the radius r is

$v\left(r\mathrm{;}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)=2\frac{{r}_{0}^{3}{v}_{0}\mathrm{cos}\left(\theta \right)}{2{r}_{0}^{3}\mathrm{cos}\left(\theta \right)-3{r}_{0}^{2}{Z}_{0}+3{r}^{2}{Z}_{0}}\mathrm{.}$ (19)

The differential equation which models the momentum conservation in the case of a hyperbolic profile is

$\left(\frac{1}{3}{r}_{0}^{3}+\frac{1}{2}\frac{{Z}_{0}\left(-{r}_{0}^{2}+{\left(r\left(t\right)\right)}^{2}\right)}{\mathrm{cos}\left(\theta \right)}\right)\frac{\text{d}}{\text{d}t}r\left(t\right)-\frac{1}{3}{r}_{0}^{3}{v}_{0}=\mathrm{0,}$ (20)

where the initial conditions are $r={r}_{0}$ and $v={v}_{0}$ when $t={t}_{0}$ .

The variables can be separated and the radius as a function of the time is

$r\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)=\frac{HN}{HD}\mathrm{,}$ (21)

where

$\begin{array}{l}HN=-\sqrt[3]{3}(2\mathrm{cos}\left(\theta \right)\sqrt[3]{3}{r}_{0}-3\sqrt[3]{3}{Z}_{0}-(-9{Z}_{0}^{3/2}+\left(\left(9t-9{t}_{0}\right){v}_{0}+9{r}_{0}\right)\mathrm{cos}\left(\theta \right)\sqrt{{Z}_{0}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+{\sqrt{3}\sqrt{27}\sqrt{AHN})}^{2/3}){r}_{0}\mathrm{,}\end{array}$ (22)

with

$\begin{array}{l}AHN=(\frac{8{\left(\mathrm{cos}\left(\theta \right)\right)}^{2}{r}_{0}^{3}}{27}+{Z}_{0}\left({\left(t-{t}_{0}\right)}^{2}{v}_{0}^{2}+2{r}_{0}\left(t-{t}_{0}\right){v}_{0}-\frac{1}{3}{r}_{0}^{2}\right)\mathrm{cos}\left(\theta \right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{}{\overset{{}^{}}{}}-2{v}_{0}{Z}_{0}^{2}\left(t-{t}_{0}\right))\mathrm{cos}\left(\theta \right)\end{array}$ (23)

and

$HD=3\sqrt{{Z}_{0}}\times \sqrt[3]{-9{Z}_{0}^{3/2}+\left(\left(9t-9{t}_{0}\right){v}_{0}+9\text{\hspace{0.05em}}{r}_{0}\right)\mathrm{cos}\left(\theta \right)\sqrt{{Z}_{0}}+9\sqrt{BHD}}\mathrm{,}$ (24)

with

$\begin{array}{l}BHD=(\frac{8{\left(\mathrm{cos}\left(\theta \right)\right)}^{2}{r}_{0}^{3}}{27}+{Z}_{0}\left({\left(t-{t}_{0}\right)}^{2}{v}_{0}^{2}+2{r}_{0}\left(t-{t}_{0}\right){v}_{0}-\frac{1}{3}{r}_{0}^{2}\right)\mathrm{cos}\left(\theta \right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{}{\overset{{}^{}}{}}-2\text{\hspace{0.05em}}{v}_{0}{Z}_{0}^{2}\left(t-{t}_{0}\right))\mathrm{cos}\left(\theta \right)\mathrm{.}\end{array}$ (25)

As a consequence the velocity as function of the time is

$v\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)=\frac{\text{d}r\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)}{\text{d}t}\mathrm{.}$ (26)

More details as well the exploration of other profiles of density can be found in [27] . We now continue evaluating the flux of kinetic energy, ${F}_{ek}$ , in the thin emitting layer which is supposed to have density ${\rho}_{l}$

${F}_{ek}\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)=\frac{1}{2}{\rho}_{l}4\text{\pi}r{\left(t\right)}^{2}v{\left(t\right)}^{3}\mathrm{.}$ (27)

The volume of the thin emitting layer, ${V}_{l}$ , is approximated by

${V}_{l}=4\Delta \text{\pi}{r}^{2}\mathrm{,}$ (28)

where $\Delta $ is thickness of the layer; as an example [28] quotes $\Delta =\frac{r}{12}$ . The two approximations for mass, Equation (14), and volume, Equation (28), allows to derive an approximate value for the density in the thin layer

${\rho}_{l}=\frac{1}{8}\frac{{\rho}_{0}{Z}_{0}f}{\mathrm{cos}\left(\theta \right)r\text{\pi}}\mathrm{.}$ (29)

Inserted in Equation (27) the radius, velocity and density as given by Equation (21), Equation (26) and Equation (29), we obtain

${F}_{ek}\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)=\frac{FN}{FD}\mathrm{,}$ (30)

where

$\begin{array}{l}FN=-\sqrt{27}{\left(-3\sqrt{3}{Z}_{0}^{3/2}+\sqrt{27}\sqrt{{F}_{1}\mathrm{cos}\left(\theta \right)}+3{F}_{5}\right)}^{3}f\\ \text{\hspace{0.05em}}\times {\left(2\sqrt[3]{3}{F}_{3}+{\left(-9{Z}_{0}^{3/2}+{F}_{2}+\sqrt{27}\sqrt{{F}_{1}\mathrm{cos}\left(\theta \right)}\sqrt{3}\right)}^{2/3}\right)}^{3}{r}_{0}^{4}\\ \text{\hspace{0.05em}}\times \mathrm{cos}\left(\theta \right){v}_{0}^{3}\sqrt{{Z}_{0}}\left(2\sqrt[3]{3}{F}_{3}-{\left(-9{Z}_{0}^{3/2}+{F}_{2}+\sqrt{27}\sqrt{{F}_{1}\mathrm{cos}\left(\theta \right)}\sqrt{3}\right)}^{2/3}\right)\sqrt[3]{3}{\rho}_{0}\end{array}$ (31)

and

$FD=108\text{\hspace{0.05em}}\sqrt{{F}_{1}\mathrm{cos}\left(\theta \right)}\times {\left(-9{Z}_{0}^{3/2}+{F}_{2}+\sqrt{27}\sqrt{{F}_{1}\mathrm{cos}\left(\theta \right)}\sqrt{3}\right)}^{13/3}{F}_{4}$ (32)

being

$\begin{array}{l}{F}_{1}=\frac{8{\left(\mathrm{cos}\left(\theta \right)\right)}^{2}{r}_{0}^{3}}{27}+\left({\left(t-{t}_{0}\right)}^{2}{v}_{0}^{2}+2{r}_{0}\left(t-{t}_{0}\right){v}_{0}-\frac{1}{3}{r}_{0}^{2}\right){Z}_{0}\mathrm{cos}\left(\theta \right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-2{v}_{0}{Z}_{0}^{2}\left(t-{t}_{0}\right),\end{array}$ (33)

${F}_{2}=\left(\left(9t-9{t}_{0}\right){v}_{0}+9{r}_{0}\right)\mathrm{cos}\left(\theta \right)\sqrt{{Z}_{0}}\mathrm{,}$ (34)

${F}_{3}=\mathrm{cos}\left(\theta \right){r}_{0}-3/2{Z}_{0}$ (35)

$\begin{array}{l}{F}_{4}=8{\left(\mathrm{cos}\left(\theta \right)\right)}^{2}{r}_{0}^{3}+27\left({\left(t-{t}_{0}\right)}^{2}{v}_{0}^{2}+2\text{\hspace{0.05em}}{r}_{0}\left(t-{t}_{0}\right){v}_{0}-\frac{1}{3}{r}_{0}^{2}\right){Z}_{0}\mathrm{cos}\left(\theta \right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-54\text{\hspace{0.05em}}{v}_{0}{Z}_{0}^{2}\left(t-{t}_{0}\right)\mathrm{,}\end{array}$ (36)

${F}_{5}=\left({v}_{0}\left(t-{t}_{0}\right)+{r}_{0}\right)\mathrm{cos}\left(\theta \right)\sqrt{3}\sqrt{{Z}_{0}}\mathrm{.}$ (37)

We now assumes that the amount of luminosity, ${L}_{theo}$ , reversed in the shocked emission is proportional to the flux of kinetic energy as given by Equation (30)

${L}_{theo}\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)\propto {F}_{ek}\left(t\mathrm{;}{t}_{0}\mathrm{,}{r}_{0}\mathrm{,}{Z}_{0}\mathrm{,}{v}_{0}\mathrm{,}\theta \right)\mathrm{.}$ (38)

The theoretical luminosity is not equal along all the SB but is function of the polar angle $\theta $ . In this framework is useful to introduce the ratio, $\kappa $ , between theoretical luminosity at $\theta $ and that one at $\theta =0$ ,

$\kappa =\frac{{L}_{theo}\left(t;{t}_{0},{r}_{0},{Z}_{0},{v}_{0},\theta \right)}{{L}_{theo}\left(t;{t}_{0},{r}_{0},{Z}_{0},{v}_{0},\theta =0\right)}.$ (39)

The above model for the theoretical luminosity is independent from the image theory, see Section 6, and does not explains the hole of luminosity visible in the shells.

4. Astrophysical Environment

We now analyze the catalogue for the giant arcs, the two giant arcs SDP.81 and A2267 and the initial astrophysical conditions for the SBs.

4.1. The Catalogue

Some parameters of the giant arcs as detected as images by cluster lensing and the supernova survey with Hubble (CLASH) which is available as a catalogue at http://vizier.u-strasbg.fr/viz-bin/VizieR, see [29] . We are interested in the arc length which is given in arcsec, the arc length to width ratio, the photometric redshift, and in the radial distance from the arc center to the cluster center in arcsec. Table 1 reports the statistical parameters of the radial distance from the arc center to the cluster center in kpc and Figures 2-4 the histogram of the frequencies in the framework of flat, ΛCDM and MTL cosmology respectively.

4.2. Single Giant Arcs

The ring associated with the galaxy SDP.81, see [30] , is characterized by a foreground galaxy at $z=0.2999$ and a background galaxy at $z=0.3042$ . This ring has been studied with the Atacama Large Millimeter/sub-millimeter Array (ALMA) by [31] - [36] and has the observed parameters as in Table 2.

Another giant arc is that in A2667 which is made by three pieces: A, B and C, see Figure 1 in [37] . The radius can be found from the equation of the circle

Table 1. Statistical parameters of the radial distance from the arc center to the cluster center in kpc.

Table 2. Observed parameters of the giants arcs.

Figure 2. Histogram of the radial distance from the arc center to the cluster center in Mpc in flat cosmology with conversion formula (5).

Figure 3. Histogram of the radial distance from the arc center to the cluster center in Mpc in ΛCDM cosmology with conversion formula (2).

Figure 4. Histogram of the radial distance from the arc center to the cluster center in Mpc in MTL cosmology with conversion formula (9).

given the three points A, B and C, see Table 3. The three pieces can be digitalized for a further comparison with a simulation, see empty red stars in Figure 8.

4.3. The Initial Conditions

We review the starting equations for the evolution of the SB [38] [39] [40] which can be derived from the momentum conservation applied to a pyramidal section. The parameters of the thermal model are
${N}^{\mathrm{*}}$ , the number of SN explosions in
$5.0\times {10}^{7}$ yr,
${Z}_{\text{OB}}$ , the distance of the OB associations from the galactic plane,
${E}_{51}$ , the energy in 10^{51} erg usually chosen equal to one,
${v}_{0}$ , the initial velocity which is fixed by the bursting phase,
${t}_{0}$ , the initial time in yr which is equal to the bursting time, and t the proper time of the SB. With the above definitions the radius of the SB is

$R=111.56{\left(\frac{{E}_{51}{t}_{7}^{3}{N}^{\mathrm{*}}}{{n}_{0}}\right)}^{\frac{1}{5}}\text{pc}\mathrm{,}$ (40)

and its velocity

$V=6.567\frac{1}{{t}_{7}^{2/5}}\sqrt[5]{\frac{{E}_{51}{N}^{\mathrm{*}}}{{n}_{0}}}\text{\hspace{0.17em}}\text{km}\cdot {\text{s}}^{-1}\mathrm{.}$ (41)

In the following, we will assume that the bursting phase ends at
$t={t}_{7,0}$ (the bursting time is expressed in units of 10^{7} yr) when
${N}_{SN}$ SN are exploded

${N}_{SN}={N}^{*}\frac{{t}_{7,0}\times {10}^{7}}{5\times {10}^{7}}.$ (42)

The two following inverted formula allows to derive the parameters of the initial conditions for the SB with ours
${r}_{0}$ expressed in pc and
${v}_{0}$ expressed in km∙s^{−1} are

${t}_{\mathrm{7,0}}=0.05878095238\frac{{r}_{0}}{{v}_{0}}\mathrm{,}$ (43)

Table 3. Radius of the giant arcs in kpc.

and

${N}^{\mathrm{*}}=2.8289\times {10}^{-7}\frac{{r}_{0}^{2}{n}_{0}{v}_{0}^{3}}{{E}_{51}}\mathrm{.}$ (44)

5. Astrophysical Simulation

We simulate a single giant arc, A2267, and then we simulate the statistics of many giant arcs.

5.1. Simulation of A2667

The final stage of the SB connected with A2267 is simulated with the parameters reported in Table 4; in particular Figure 5 displays the 3D shape and Figure 6 reports the 2D section.

Figure 7 reports the 2D section of the SB as well the three pieces of the giant arc connected with A2267. The similarity between the observed radius of curvature of the giant arc as well the theoretical one is reported in a zoom, see Figure 8.

We can understand the reason for which the giant arc A2267 has a limited angular extension of ≈31˚ by plotting the ratio $\kappa $ , equation (37), between the theoretical luminosity as function of $\theta $ and the theoretical luminosity at $\theta =0$ with parameters as in Table 4, see Figure 9. As a practical example at ≈31˚/2, where the factor two arises from the symmetry of the framework, the theoretical luminosity is decreased of a factor $\kappa =0.987$ in respect to the value at $\theta =0$ . We now introduce the threshold luminosity, ${L}_{tr}$ , which is an observational parameter. The theoretical luminosity will scale as function of the polar angle as ${L}_{theo}\left(\theta \right)\propto {L}_{0}\ast r$ where ${L}_{0}$ is the theoretical luminosity at $\theta =0$ and $\kappa $ has been defined in Equation (39). When the inequality ${L}_{theo}<{L}_{tr}$ is verified the giant arc is impossible to detect and only the zone characterized by low values of the polar angle will be detected.

In our model the velocity with parameters as in Table 4 is function of the polar angle, see Figure 10, and has range $37\text{\hspace{0.17em}}\text{km}/\text{s}<v\left(\theta \right)<142\text{\hspace{0.17em}}\text{km}/\text{s}$ . As a comparison a velocity $50\text{\hspace{0.17em}}\text{km}/\text{s}<v<75\text{\hspace{0.17em}}\text{km}/\text{s}$ is measured in A2267, see Figure 5 in [37] .

5.2. Simulation of Many Giants Arcs

The presence of multiple giants arcs in the CLASH cluster, see as an example Figure 11 in [29] , can be simulated adopting the following steps:

Table 4. Theoretical parameters of the SB connected with A2267.

Figure 5. 3D surface of the SB connected with A2267, parameters as in Table 4 and axes in pc.

• A given number of SBs, as an example 15, are generated with variable lifetime, t, see Figure 11.

• For each SB we select a section around polar angle equal to zero characterized by a fixed angle of ≈31˚ and we randomly rotate it around the origin, see Figure 12.

Figure 6. 2D section in the z = 0 plane of the SB connected with A2267, parameters as in Table 4 and axes in kpc.

Figure 7. 2D section in the z = 0 plane of the SB connected with A2267, parameters as in Table 4 (full green points) and the three pieces of the giant arc in A2667 (empty red stars); axes in arcsec.

Figure 8. Enlarged view of the three pieces of the giant arc in A2667 (empty red stars) and the theoretical radius (full green points); axes in arcsec.

Figure 9. Normalized luminosity as function of the polar angle in deg, parameters as in Table 4.

Figure 10. Velocity in km/s as function of the polar angle in deg.

Figure 11. Multiple sections of the SB with time, t comprised in [10^{6} yr, 10^{8} yr] and other parameters as in Table 4.

Figure 12. Multiple sections of the SB as in Figure 11 with angular extension of the polar angle, θ, of ≈31˚ and progressive rotation of the selected piece of section.

• The centers of the SBs are randomly placed in a squared box with side of 300 kpc, see Figure 13.

Table 5 reports the theoretical statistical parameters of the above simulation for the radial distance from the arc center to the cluster center in kpc. A comparison should be done with the astronomical parameters for the CLASH clusters of Table 1.

6. Theory of the Image

We now review the theory of the image for the case of optically thin medium both from an analytical and an analytical point of view.

6.1. The Elliptical Shell

A real ellipsoid represents a first approximation of the asymmetric giants arcs and has equation

$\frac{{z}^{2}}{{a}^{2}}+\frac{{x}^{2}}{{b}^{2}}+\frac{{y}^{2}}{{d}^{2}}=1,$ (45)

in which the polar axis is the z-axis.

We are interested in the section of the ellipsoid $y=0$ which is defined by the following external ellipse

$\frac{{z}^{2}}{{a}^{2}}+\frac{{x}^{2}}{{b}^{2}}=1.$ (46)

We assume that the emission takes place in a thin layer comprised between the external ellipse and the internal ellipse defined by

$\frac{{z}^{2}}{{\left(a-c\right)}^{2}}+\frac{{x}^{2}}{{\left(b-c\right)}^{2}}=1,$ (47)

Table 5. Statistical parameters of the radial distance from the theoretical arc center to the cluster center in kpc.

Figure 13. Multiple sections of SB as in Figure 12 with random shift of the origin of the selected SB (green empty stars). The random shift denotes the galaxies (red crosses).

see Figure 14. We therefore assume that the number density C is constant and in particular rises from 0 at (0, a) to a maximum value ${C}_{m}$ , remains constant up to (0, a-c) and then falls again to 0. The length of sight, when the observer is situated at the infinity of the x-axis, is the locus parallel to the x-axis which crosses the position z in a Cartesian x-z plane and terminates at the external ellipse. The locus length is

${l}_{I}=2\frac{\sqrt{{a}^{2}-{z}^{2}}b}{a}$ (48)

when $\left(a-c\right)\le z<a$

${l}_{II}=2\frac{\sqrt{{a}^{2}-{z}^{2}}b}{a}-2\frac{\sqrt{{a}^{2}-2ac+{c}^{2}-{z}^{2}}\left(b-c\right)}{a-c}$ (49)

when $0\le z<\left(a-c\right)$

In the case of optically thin medium, the intensity is split in two cases

${I}_{I}\left(z;a,b\right)={I}_{m}\times 2\frac{\sqrt{{a}^{2}-{z}^{2}}b}{a}$ (50)

when $\left(a-c\right)\le z<a$

${I}_{II}\left(z;a,c\right)={I}_{m}\times \left(2\frac{\sqrt{{a}^{2}-{z}^{2}}b}{a}-2\frac{\sqrt{{a}^{2}-2ac+{c}^{2}-{z}^{2}}\left(b-c\right)}{a-c}\right)$ (51)

when $0\le z<\left(a-c\right)$

Figure 14. Internal and external ellipses when $a=347\text{\hspace{0.17em}}\text{kpc}$ , $b=237\text{\hspace{0.17em}}\text{kpc}$ and $c=\frac{a}{12}\text{\hspace{0.17em}}\text{kpc}$ .

where ${I}_{m}$ is a constant which allows to compare the theoretical intensity with the observed one. A typical profile in intensity along the z-axis is reported in Figure 15. The ratio, $\kappa $ , between the theoretical intensity at the maximum, ( $z=a-c$ ), and at the minimum, ( $z=0$ ), is given by

$\frac{{I}_{I}\left(z=a-c\right)}{{I}_{II}\left(z=0\right)}=\kappa =\frac{\sqrt{2\text{\hspace{0.05em}}a-c}b}{\sqrt{c}a}.$ (52)

As an example the values $a=6\text{\hspace{0.17em}}\text{kpc}$ , $b=4\text{\hspace{0.17em}}\text{kpc}$ , $c=\frac{a}{12}\text{\hspace{0.17em}}\text{kpc}$ gives $\kappa =3.19$ . The knowledge of the above ratio from the observations allows to deduce c once a and b are given by the observed morphology

$c=2\frac{a{b}^{2}}{{a}^{2}{r}^{2}+{b}^{2}}.$ (53)

The above analytical model explains the hole in luminosity visible in the astrophysical shells such as supernovae and SBs. More details can be found in [41] .

Figure 15. The intensity profile along the z-axis when $a=347\text{\hspace{0.17em}}\text{kpc}$ , $b=237\text{\hspace{0.17em}}\text{kpc}$ $c=\frac{a}{12}\text{kpc}$ and ${I}_{m}=1$ .

6.2. The Numerical Shell

The source of luminosity is assumed here to be the flux of kinetic energy, ${L}_{m}$ ,

${L}_{m}=\frac{1}{2}\rho A{V}^{3},$ (54)

where A is the considered area, V is the velocity and $\rho $ is the density. In our case $A={r}^{2}\Delta \Omega $ , where $\Delta \Omega $ is the considered solid angle and $r\left(\theta \right)$ the temporary radius along the chosen direction. The observed luminosity along a given direction can be expressed as

$L=\u03f5{L}_{m}\mathrm{,}$ (55)

where $\u03f5$ is a constant of conversion from the mechanical luminosity to the observed luminosity.

We review the algorithm that allows to build the image, see [42] :

• An empty memory grid
$\mathcal{M}\left(i\mathrm{,}j\mathrm{,}k\right)$ which contains NDIM^{3} pixels is considered.

• We first generate an internal 3D surface of revolution by rotating the ideal image of 360˚ around the polar direction and a second external surface of revolution at a fixed distance $\Delta R$ from the first surface. As an example, we fixed $\Delta R=R/12$ , where R is the momentary radius of expansion. The points on the memory grid which lie between the internal and external surfaces are memorized on $\mathcal{M}\left(i\mathrm{,}j\mathrm{,}k\right)$ by a variable integer number according to formula (52) and density $\rho $ proportional to the swept mass.

• Each point of $\mathcal{M}\left(i\mathrm{,}j\mathrm{,}k\right)$ has spatial coordinates $x\mathrm{,}y\mathrm{,}z$ which can be represented by the following $1\times 3$ matrix, A,

$A=\left[\begin{array}{c}x\\ y\\ z\end{array}\right].$ (56)

The orientation of the object is characterized by the Euler angles $\left(\Phi \mathrm{,}\Theta \mathrm{,}\Psi \right)$ and therefore by a total $3\times 3$ rotation matrix, E. The matrix point is represented by the following $1\times 3$ matrix, B,

$B=E\cdot A\text{\hspace{0.05em}}\text{\hspace{0.05em}}.$ (57)

• The intensity map is obtained by summing the points of the rotated images along a particular direction.

The image of A2267 built with the above algorithm is shown in Figure 16.

The threshold intensity, ${I}_{tr}$ , is

${I}_{\mathrm{max}}\text{\hspace{0.05em}}\kappa ={I}_{\mathrm{max}},$ (58)

where ${I}_{\mathrm{max}}$ , is the maximum value of intensity characterizing the ring and $\kappa $ is a parameter which allows matching theory with observations and was previously defined in Equation (39). A typical image with a hole is visible in Figure 17. The opening angle of the visible arc can be parametrized as function of the ratio $\kappa $ , see Figure 18. An opening of ≈31˚ is reached at $\kappa \approx 0.95$ .

7. Conclusions

7.1. The Equation of Motion

The giants arcs are connected with the visible part of the SBs which advance in the intracluster medium surrounding the host galaxies. The chosen profile of density is hyperbolic, see Equation (10), and the momentum conservation along

Figure 16. Contour map of I for A2267, the x and y axes are in arcsec. The three Euler angles characterizing the orientation are $\Phi ={0}^{\circ}$ , $\Theta ={90}^{\circ}$ and $\Psi ={90}^{\circ}$ , and NDIM = 400.

Figure 17. The same as Figure 16 parameters as in Figure 16 and $\kappa =0.5$ .

Figure 18. The opening angle of the visible simulated arc as function of $\kappa $ ; parameters as in Figure 16.

a given direction allows to derive the equation of motion as function of the polar angle, see Equation (21).

7.2. The Image

According to the theory here presented the giants arcs are the visible part of an advancing SB. An analytical explanation for the limited angular extent of the giant arcs is represented by the theoretical luminosity as function of the polar angle, see Equation (39). An increase in the polar angle produces a decrease of the theoretical luminosity and the arc becomes invisible. Selecting a given numbers of SBs with variable lifetime and randomly inserting them in a cubic box of side ≈ 600 kpc is possible to simulate the giants arcs visible in the clusters of galaxies, see Figure 13 and relative statistical parameters in Table 5.

Acknowledgements

This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France.

References

[1] Lynds, R. and Petrosian, V. (1986) Giant Luminous Arcs in Galaxy Clusters. Bulletin of the American Astronomical Society, 18, 1014.

[2] Paczynski, B. (1987) Giant Luminous Arcs Discovered in Two Clusters of Galaxies. Nature, 325, 572-573.

https://doi.org/10.1038/325572a0

[3] Soucail, G., Fort, B., Mellier, Y. and Picat, J. (1987) A Blue Ring-Like Structure, in the Center of the A 370 Cluster of Galaxies. Astronomy and Astrophysics, 172, L14-L16.

[4] Kovner, I. (1987) Giant Luminous Arcs from Gravitational Lensing. Nature, 327, 193-194.

https://doi.org/10.1038/327193c0

[5] Waldrop, M.M. (1987) The Giant Arcs Are Gravitational Mirages. Science, 238, 1351-1352.

https://doi.org/10.1126/science.238.4832.1351

[6] Soucail, G., Mellier, Y., Fort, B., Mathez, G. and Cailloux, M. (1988) The Giant Arc in a 370—Spectroscopic Evidence for Gravitational Lensing from a Source at Z = 0.724. Astronomy and Astrophysics, 191, L19-L21.

[7] Narasimha, D. and Chitre, S.M. (1988) Giant Luminous Arcs in Galaxy Clusters. The Astrophysical Journal, 332, 75-80.

https://doi.org/10.1086/166630

[8] Kaufmann, R. and Straumann, N. (2000) Giant Arc Statistics and Cosmological Parameters. Annalen der Physik, 9, 384-394.

https://doi.org/10.1002/(SICI)1521-3889(200005)9:3/5<384::AID-ANDP384>3.0.CO;2-Y

[9] Wambsganss, J., Bode, P. and Ostriker, J.P. (2004) Giant Arc Statistics in Concord with a Concordance Lambda Cold Dark Matter Universe. The Astrophysical Journal, 606, L93-L96.

[10] Dalal, N., Holder, G. and Hennawi, J.F. (2004) Statistics of Giant Arcs in Galaxy Clusters. The Astrophysical Journal, 609, 50-60.

https://doi.org/10.1086/420960

[11] Li, G.L., Mao, S., Jing, Y.P., Mo, H.J., Gao, L. and Lin, W.P. (2006) The Giant Arc Statistics in the Three-Year Wilkinson Microwave Anisotropy Probe Cosmological Model. Monthly Notices of the Royal Astronomical Society, 372, L73-L77.

https://doi.org/10.1111/j.1745-3933.2006.00230.x

[12] Wambsganss, J., Ostriker, J.P. and Bode, P. (2008) The Effect of Baryon Cooling on the Statistics of Giant Arcs and Multiple Quasars. The Astrophysical Journal, 676, 753-760.

https://doi.org/10.1086/527529

[13] Bayliss, M.B. (2012) Broadband Photometry of 105 Giant Arcs: Redshift Constraints and Implications for Giant Arc Statistics. The Astrophysical Journal, 744 156-172.

https://doi.org/10.1088/0004-637X/744/2/156

[14] Hammer, F. and Rigaut, F. (1989) Giant Luminous Arcs from Lensing—Determ- ination of the Mass Distribution inside Distant Cluster Cores. Astronomy and Astrophysics, 226, 45-46.

[15] Wu, X.P. and Mao, S. (1996) The Cosmological Constant and Statistical Lensing of Giant Arcs. The Astrophysical Journal, 463, 404-408.

https://doi.org/10.1086/177256

[16] Lewis, G.F. (2001) Gravitational Microlensing of Giant Luminous Arcs: A Test for Compact Dark Matter in Clusters of Galaxies. Publications of the Astronomical Society of Australia, 18, 182-185.

https://doi.org/10.1071/AS01025

[17] Mahdi, H.S., van Beek, M., Elahi, P.J., Lewis, G.F., Power, C. and Killedar, M. (2014) Gravitational Lensing in WDM Cosmologies: The Cross-Section for Giant Arcs. Monthly Notices of the Royal Astronomical Society, 441, 1954-1963.

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

[18] Dekel, A. and Braun, E. (1988) Giant Arcs-Spherical Shells? In: Audouze, J., Pelletan, M.C., Szalay, A., Zel’dovich, Y.B. and Peebles, P.J.E., Eds., Large Scale Structures of the Universe, Vol. 130, IAU Symposium, 598.

[19] Braun, E. and Dekel, A. (1988) On the Giant Arcs in Clusters of Galaxies—Can They Be Shells? Comments on Astrophysics, 12, 233-244.

[20] Efremov, Y.N., Elmegreen, B.G. and Hodge, P.W. (1998) Giant Shells and Stellar Arcs as Relics of Gamma-Ray Burst Explosions. The Astrophysical Journal, 501, L163-L165.

https://doi.org/10.1086/311468

[21] Zaninetti, L. (2017) The Ring Produced by an Extra-Galactic Superbubble in Flat Cosmology. Journal of High Energy Physics, Gravitation and Cosmology, 3, 339-359.

https://doi.org/10.4236/jhepgc.2017.32029

[22] Zaninetti, L. (2016) Padé Approximant and Minimax Rational Approximation in Standard Cosmology. Galaxies, 4, 4-28.

http://www.mdpi.com/2075-4434/4/1/4

[23] Suzuki, N., Rubin, D., Lidman, C., Aldering, G., Amanullah, R., Barbary, K. and Barrientos, L.F. (2012) The Hubble Space Telescope Cluster Supernova Survey: V. Improving the Dark Energy Constraints above z > 1 and Building an Early-Type-Hosted Supernova Sample. The Astrophysical Journal, 746.

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

[24] Zaninetti, L. (2015) On the Number of Galaxies at High Redshift. Galaxies, 3, 129-155.

https://doi.org/10.3390/galaxies3030129

[25] Peebles, P.J.E. (1993) Principles of Physical Cosmology. Princeton University Press, Princeton, NJ.

[26] McKee, C.F. (1987) Astrophysical Shocks in Diffuse Gas. In: Dalgarno, A. and Layzer, D., Eds., Spectroscopy of Astrophysical Plasmas, Cambridge University Press, Cambridge, 226-252.

https://doi.org/10.1017/CBO9780511564659.010

[27] Zaninetti, L. (2018) The Physics of Asymmetric Supernovae and Supernovae Remnants. Advances in Applied Physics, 6, 25-73.

https://doi.org/10.12988/aap.2018.897

[28] McCray, R.A. (1987) Coronal Interstellar Gas and Supernova Remnants. In: Dalgarno, A. and Layzer, D., Eds., Spectroscopy of Astrophysical Plasmas, Cambridge University Press, Cambridge, 255-278.

https://doi.org/10.1017/CBO9780511564659.011

[29] Xu, B., Postman, M., Meneghetti, M., Seitz, S., Zitrin, A., Merten, J., Maoz, D., Frye, B., Umetsu, K., Zheng, W., Bradley, L., Vega, J. and Koekemoer, A. (2016) The Detection and Statistics of Giant Arcs behind CLASH Clusters. The Astrophysical Journal, 817, 85-110.

https://doi.org/10.3847/0004-637X/817/2/85

[30] Eales, S., Dunne, L., Clements, D. and Cooray, A. (2010) The Herschel ATLAS. Astronomical Society of the Pacific, 122, 499-515.

https://doi.org/10.1086/653086

[31] Tamura, Y., Oguri, M., Iono, D., Hatsukade, B., Matsuda, Y. and Hayashi, M. (2015) High-Resolution ALMA Observations of SDP. 81. I. The Innermost Mass Profile of the Lensing Elliptical Galaxy Probed by 30 Milli-Arcsecond Images. Publications of the Astronomical Society of Japan, 67, 72-78.

https://doi.org/10.1093/pasj/psv040

[32] ALMA Partnership, Vlahakis, C., Hunter, T.R. and Hodge, J.A. (2015) The 2014 ALMA Long Baseline Campaign: Observations of the Strongly Lensed Submillimeter Galaxy HATLAS J090311.6 + 003906 at z = 3.042. The Astrophysical Journal, 808.

[33] Rybak, M., Vegetti, S., McKean, J.P., Andreani, P. and White, S.D.M. (2015) ALMA Imaging of SDP. 81—II. A Pixelated Reconstruction of the CO Emission Lines. Monthly Notices of the Royal Astronomical Society, 453, L26-L30.

[34] Hatsukade, B., Tamura, Y., Iono, D., Matsuda, Y., Hayashi, M. and Oguri, M. (2015) High-Resolution ALMA Observations of SDP. 81. II. Molecular Clump Properties of a Lensed Submillimeter Galaxy at z = 3.042. Publications of the Astronomical Society of Japan, 67, 93.

https://doi.org/10.1093/pasj/psv061

[35] Wong, K.C., Suyu, S.H. and Matsushita, S. (2015) The Innermost Mass Distribution of the Gravitational Lens SDP. 81 from ALMA Observations. The Astrophysical Journal, 811, 115-125.

https://doi.org/10.1088/0004-637X/811/2/115

[36] Hezaveh, Y.D., Dalal, N. and Marrone, D.P. (2016) Detection of Lensing Substructure Using ALMA Observations of the Dusty Galaxy SDP. 81. The Astrophysical Journal, 823, 37-56.

https://doi.org/10.3847/0004-637X/823/1/37

[37] Yuan, T.T., Kewley, L.J., Swinbank, A.M. and Richard, J. (2012) The A2667 Giant Arc at z = 1.03: Evidence for Large-Scale Shocks at High Redshift. The Astrophysical Journal, 759, 66-77.

https://doi.org/10.1088/0004-637X/759/1/66

[38] Dyson, J.E. and Williams, D.A. (1997) The Physics of the Interstellar Medium. Institute of Physics Publishing, Bristol.

[39] McCray, R. and Kafatos, M. (1987) Supershells and Propagating Star Formation. The Astrophysical Journal, 317, 190-196.

https://doi.org/10.1086/165267

[40] Zaninetti, L. (2004) On the Shape of Superbubbles Evolving in the Galactic Plane. Publications of the Astronomical Society of Japan, 56, 1067-1082.

https://doi.org/10.1093/pasj/56.6.1067

[41] Zaninetti, L. (2018) The Fermi Bubbles as a Superbubble. International Journal of Astronomy and Astrophysics, 8, 200-217.

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

[42] Zaninetti, L. (2013) Three Dimensional Evolution of SN 1987A in a Self-Gravitating Disk. International Journal of Astronomy and Astrophysics, 3, 93-98.

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