─юъєьхэЄ тч Є шч ъ¤°р яюшёъютющ ьр°шэ√. └фЁхё юЁшушэры№эюую фюъєьхэЄр : http://jet.sao.ru/hq/lrs/AL179.PDF
─рЄр шчьхэхэш : Wed Feb 24 17:19:48 2010
─рЄр шэфхъёшЁютрэш : Mon Oct 1 20:35:14 2012
╩юфшЁютър:
╧юшёъют√х ёыютр: m 5

Astronomy Letters, Vol. 31, No. 3, 2005, pp. 179н185. Translated from Pis'ma v Astronomicheski Zhurnal, Vol. 31, No. 3, 2005, pp. 198н205. i Original Russian Text Copyright c 2005 by Gosachinskii.

Interaction of the Supernova Remnant 3 with the Ambient Interstellar Gas
I. V. Gosachinskii
*

Special Astrophysical Observatory, St. Petersburg Branch, Russian Academy of Sciences, Pulkovo, St. Petersburg, 196140 Russia
Received August 4, 2004; in final form,

Abstract--The well-known shell supernova remnant (SNR) HB3 is part of a feature-rich star-forming region together with the nebulae W3, W4, and W5. We study the H I structure around this SNR using five RATAN-600 drift curves obtained at a wavelength of 21 cm with an angular resolution of 2 in one coordinate over the radial-velocity range -183 to +60 km s-1 in a wider region of the sky and with a higher sensitivity than in previous works of other authors. The spatialнkinematic distribution of HI features around the SNR clearly shows two concentric expanding shells of gas that surround the SNR and coincide with it in all three coordinates (, , and V ). The outer shell has a radius of 133 pc, a thickness of 24 pc, and an expansion velocity of 48 km s-1 . The mass of the gas in it is 2.3 т 105 M . For the inner shell, these parameters are 78 pc, 36 pc, 24 km s-1 , and 0.9 т 105 M , respectively. The inner shell is immediately adjacent to the SNR. Assuming that the outer shell was produced by the stellar wind and the inner shell arose from the shock wave of the SNR proper, we estimated the age of the outer shell, 1.7 т 106 yr, and the mechanical luminosity of the stellar wind, 1.5 т 1038 erg s-1 . The inner shell has an age of 106 yr and corresponds to a total supernova explosion energy of 1052 erg. c 2005 Pleiades Publishing, Inc. Key words: interstellar medium, gaseous nebulae, neutral hydrogen, supernovae and supernova remnants.

poor sensitivity (3н15 K) and found nothing new in this region compared to Routledge et al. (1991). Note The supernova remnant (SNR) HB3 (G132.6+1.5) also that Routledge et al. (1991) gave a brief review is well known, since it is part of a rich star-forming of the then available observations of HB3 in the radio, region in the Cas OB6 stellar association together X-ray, and optical ranges and in molecular lines. with the H II regions W3, W4, and W5. This SNR is In the period 2001н2004, we obtained five drift commonly assumed to be physically associated with curves through the SNR HB3 in the HI line on the them. At radio wavelengths, it has a significant angusouthern sector of the RATAN-600 radio telescope in lar size ( 80 ) and a shell structure, suggesting that a region that was twice that in Routledge et al. (1991) its evolution time is fairly long. This SNR has been and with a sensitivity that was a factor of 20 higher extensively studied in the continuum, infrared, X rays, than that in Normandeau et al. (1997). We confirmed and molecular lines (see Landecker et al. (1987) and several conclusions of the former paper and found new references therein). features in the distribution of gas around the SNR. Neutral hydrogen in this region has been studNote also that the distribution of neutral hydroied by several authors. The papers by Routledge et gen emission in the Galactic plane is highly nonunial. (1991) and Normandeau et al. (1997), who used form and rich in various features. It thus follows that the DRAO aperture-synthesis radio telescope with an searching for evidence of the physical association of angular resolution of 2 and 1 , are particularly note- HI gas with such objects as SNRs (i.e., evidence worthy. Note, however, that Routledge et al. (1991) that this is not a chance projection) is not an easy restricted themselves to a field of view of one primary task. In our case, however, the problem is simplified, beam (3 at the 10% level) with an angular resolution since the SNR is an active expanding object, so this of 2 and a sensitivity of 1.34 K, so the sensitivity at the must also be reflected in the kinematics of the ambient edges of the field deteriorated sharply during the data gas. Unfortunately, many authors pay little attention processing. Normandeau et al. (1997) had a rather to this circumstance and restrict their analysis to the distribution of gas only in the plane of the sky. * E-mail: gos@fsao.spb.su The goal of our work is to search for expanding HI
1063-7737/05/3103-0179$26.00 c 2005 Pleiades Publishing, Inc.

INTRODUCTION

180

GOSACHINSKII

shells around the SNR, which manifest themselves as ringlike structures in the coordinateнradial velocity plane. INSTRUMENTATION AND TECHNIQUE To study the distribution of interstellar neutral hydrogen in the region of the SNR HB3, we obtained . five drift curves in right ascension at 06 declination . intervals in the range of declinations from +6125 . to +6375 on the southern sector of RATAN-600 . at lower culmination (an elevation of 173). In this elevation range, the RATAN-600 antenna has an angular resolution of 2.6 т 65 and an effective area of about 900 m2 (Esepkina et al. 1979). An uncooled HEMT amplifier (Il'in et al. 1997) was used at the input. The system noise temperature was about 70 K, the 39-channel filter-bank spectral analyzer had a channel width of 30 kHz (6.3 km s-1 ), and the channel spacing was also 30 kHz (Venger et al. 1982). The system control, data acquisition, and primary data processing were performed with an IBM PC computer (Alferova et al. 1986). The drift curve at each declination consisted of two series with three observations each obtained with a receiver tuning shift by half the channel width, so each drift curve had 78 spectral channels following at 3.15 km s-1 intervals. This technique also allowed us to effectively eliminate the interference. On the averaged record, the mean square of the antenna temperature fluctuations in spectral channels was 0.2 . In each series of observations, the antenna and instrumentation parameters were controlled using measurements of a series of reference sources (Venger et al. 1981). For the southern sector, we used 3C 147 and 3C 430. We determined the parameters of features in each channel using a Gaussian fitting program and then attempted to establish a relationship between HI emission features at different radial velocities and at different declinations. This is the only procedure where a significant subjective factor could be introduced. The measured parameters of HI features have the following errors. The radial velocity of an isolated medium-brightness HI feature is measured with an accuracy of no lower than 1 km s-1 . In certain cases, the accuracy deteriorates due to the difficulties in separating an object from the background or from adjacent features. Given the antenna calibration errors, the measurement error of the brightness temperature of the HI line is about 0.3 K, and the error in the . estimated angular sizes in right ascension is 01. In declination, the antenna resolution is much lower, and, therefore, the accuracy of measuring the angular

sizes is lower. The accuracy of estimating the distances depends on the method of their determination and, in each case, must be considered separately. As a result, the accuracy of estimating the HI mass in an isolated cloud is no better than 0.5н1 order of magnitude. RESULTS OF THE OBSERVATIONS Figures 1н5 show the нV diagrams for each drift curve without elimination of the extended HI emission background. The latter was done to avoid the systematic errors that appear in any technique of eliminating an extended background. The declinations (1950.0) of the drift curves are indicated at the top. To save space, the drift curves are given only in a limited range of radial velocities, from about +9 to -100 km s-1 , and no zero lines of each curve are given. The radial velocity with respect to the Local Standard of Rest is indicated near each curve on the right. The spectrometer channel number and the antenna temperature scale are indicated on the left. Compact HI features in the SNR region were identified using linear interpolation, which, of course, can cause their antenna temperature to be slightly overestimated. In each figure, the upper curve represents the drift curve in the continuum at a frequency of 1420 MHz. The SNR is highlighted by dark shading; in the remaining part of the figure, the heavy vertical lines corresponding to the SNR boundaries at minimum brightness in right ascension were drawn. The same dark shading in some of the figures highlights the absorption line from bright thermal continuum sources (W3, W4). The thin closed lines represent the results of our attempts to combine the HI features at adjacent radial velocities into single structures surrounding the SNR. In our opinion, two such structures are observed, outer and inner. The features of the drift curves pertaining to different structures are highlighted by shading of different densities. The presence of such ringlike (нV ) structures generally reflects the expanding shells of gas around the object. The table lists the measured and calculated parameters of neutral hydrogen in these structures under the assumption that they are indeed expanding shells around the SNR. We adopt the distance to the HI shells, 3 ▒ 1 kpc, in accordance with their mean radial velocities and the IAU 1986 model of Galactic rotation (Kerr et al. 1986). The mean radial velocities that we obtained are very close to the optical velocities of the HB3 features inferred from the H observations by Lozinskaya and Sitnik (1980), and the adopted distance corresponds to this work. We see that the HI distribution in the inner shell is highly asymmetric; therefore, the table gives two values of the HI density, separately for the eastern and western parts of the
ASTRONOMY LETTERS Vol. 31 No. 3 2005

INTERACTION OF THE SUPERNOVA REMNANT 3
Ta, K 3 2 1 21 22
23 24 25 26 27 28 н 12.9 29 н 16.0 30 н 19.2 31 н 22.4 32 33 34 35 36 37 н 25.5

181

+63.75
30 25 20 15

Continuum
10 5 55 9.3 6.1 2.9 н0.2 н3.4 н6.5 н9.7

2н 00

Ta, K 3 2 1 0
н1 40 17 18 19 20 21 22 35 30 25

+63.13 Continuum

20

15

10

5

2н 00

55 9.3 6.2 3.0 н 0.2 н 3.3 н 6.5

23 н 9.7 24 25
н 28.7 н 31.9 н 35.0 н 38.2

н 12.8 н 16.0

26 27

н 19.1 н 22.3

28

н 25.5 н н н н н 28.6 31.8 31.8 35.0 38.1

н 41.4 н 44.5 н 47.7

Radial velocity, km sн1

29 30 30 31 32 33 34 35 36 37

Ta, K 20 10 0

38 н 50.9 39 40 41 42 н 54.0

н 41.3 н 44.5 н 47.6 н 50.8 н 54.0 н57.1

н57.2

43 44 45 46 47 48 49 50 51 52 53

н 60.4 н 63.5 н 66.7 н 69.9 н 73.0 н 76.2 н 79.4 н 82.5 н 85.7 н 88.9

38 н 60.3

Ta, K 20 10 0

39 40 41 42 43 44 45 46

н 63.5 н 66.6 н 69.8 н 73.0 н 76.1 н 79.3 н 82.5 н 85.6 н 88.8 н 92.0 н 95.1

н 92.0

н 95.2 54 55 35 30 25 20 15 10 5 2н00 55 н 98.4

47 48 49 50 51

Right ascension (1950.0)
Fig. 1. Drift curves of the HI line emission in the region . of the SNR HB3 at a declination of +6375. The SNR continuum is highlighted by the black color on the upper curve. The radial velocities of each curve are given on the right. The gray shading marks the HI features that may be associated with the SNR.

н 98.3 35 30 25 20 15 10 5 2н00 55

Right ascension (1950.0)
Fig. 2. Same as Fig. 1 at a declination of +6313. .

shell. The knife-edge beam of the radio telescope makes it impossible to estimate the size of the HI features in declination with an acceptable accuracy. However, it can be seen that the observed radialASTRONOMY LETTERS Vol. 31 No. 3 2005

velocity range of the HI ringlike structures decreases northward and southward of the central drift curve. The simplest physical interpretation of the HI shell structure that we detected around HB3 is that the

Radial velocity, km sн1

182

GOSACHINSKII

Ta, K 3 2 1 0
35 30 31 32 33 34 35 36 37 38 39 40 41 30

+62.5 W3 Continuum HB3
25 20 15 10 5 2н 00 55 н 0.1 30 19 н 9.6 н 12.8 н 15.9 н 19.0 н 22.3 н 25.4 н 28.6 н 31.8 н 34.9 н 38.1 н 41.3 н 44.4 24 25 26 27 28 29 30 31 н 12.8 н 16.0 н 19.1 н 22.3 н н н н 25.5 28.6 31.8 35.0 20 21 22 23 н 6.5 н 9.7 н 0.2 н 3.3 25 20 15 10 5 55 3.0

Ta, K 5 4 3 2 1 W3 0
н1

+61.88

Continuum

н 6.4

2н 00

Ta, K 20 10 0

42 43

н 47.6 44 45 46 47 н57.1 48 н 60.3 49 н 63.4 50 51 52 53 н 72.9 54 55 56 57 58 59 60 61 62 63 64 65 35 30 25 20 15 10 5 2н00 55 н 76.1 н 79.3 н 82.4 н 85.6 н 88.8 н 91.9 н 95.1 н 98.3 н 101.4 н 104.6 48 н 107.8 н 110.9 49 50 30 25 20 15 10 5 2н00 55 40 39 н 66.6 н 69.8 н 50.7

Radial velocity, km sн1

н 38.1

н 53.9

Ta, K 20 10 0
33 34 35 36 37 38

н 44.5 н 47.6 н 50.8 н 54.0 н57.1 н 60.3

н 63.5 н 66.6

41 42 43 44 45 46 47 н 69.8 н 73.0 н 76.1 н 79.3 н 82.5 н 85.6 н 88.8 н 92.0 н 95.1

Right ascension (1950.0)
Fig. 3. Same as Fig. 1 at a declination of +625. .

Right ascension (1950.0)
Fig. 4. Same as Fig. 1 at a declination of +6188. .

outer shell may be assumed to be produced by the wind from the star over its lifetime on the main sequence, while the inner shell arose from the SNR

shock. In this case, we can estimate several physical properties of the star and its explosion using theoretical estimates of the effect of shock waves from the
ASTRONOMY LETTERS Vol. 31 No. 3 2005

Radial velocity, km sн1

32

н 41.3

INTERACTION OF THE SUPERNOVA REMNANT 3

183

Ta, K 5 W3 4 3 2 1 0 W4

+61.25

Continuum HB-3
20 15 10

30 23 24 25 26 27 28 29 30 31 32 33 34 35

25

5

2н 00

55 2.9 н 0.2 н 3.4 н 6.5 н 9.7 н 12.9 н 16.0

н 19.2 н 22.4 н 25.5 н 28.7

н 35.0 н 38.2

Ta, K 20

36

н 41.4 н 44.5 н 47.7 н 50.9

37

10 0

38 39

н 54.0 40 41 42 н 60.4 43 44 45 46 47 н 73.0 48 49 50 51 н 76.1 н 79.4 н 82.5 н 85.7 н 63.5 н 66.7 н 69.9 н57.2

30

25

20

15

10

5

2н00

55

Right ascension (1950.0)
Fig. 5. Same as Fig. 1 at a declination of +6125. .

stellar wind and SNRs on the interstellar medium. The initial density of the medium needed for such estimates can be obtained by "spreading" the total amount of gas in the shells over the entire current volume.
ASTRONOMY LETTERS Vol. 31 No. 3 2005

To estimate the mechanical luminosity of the stellar wind from the observed parameters of the outer shell, we can use the model by Weaver et al. (1977). Our estimations yield an age of the shell 1.7 т 106 yr and a mechanical luminosity of the stellar wind 1.5 т

Radial velocity, km sн1

н 31.9

184 Measured parameters of the HI shells Parameters Coordinates (1950.0) Radial velocity, km s Distance, kpc Outer radius, pc Thickness, pc HI mass, M HI density, cm
-3 -1 -1

GOSACHINSKII

Outer shell = 2h 15 -44.4 47.5 3 133 24 2.3 т 10 2.1
-3 5 m

Inner shell = 2h 15 -42.3 24.3 3 78 36 0.9 т 10
5 m

= +625 .

= +625 .

Expansion velocity, km s

(E) 2.46 (W) 1.23 1.85

Initial HI density, cm

0.95

1038 erg s-1 . At a typical stellar-wind particle velocity of 2000 km s-1 , this mechanical luminosity requires a relatively high mass-loss rate, 6.8 т 10-5 M yr-1 , which is characteristic of massive Otype stars or M-type supergiants (de Jager 1984). It is well known that for WolfнRayet stars, the wind velocity can be several times higher; however, this stage is not too long and unlikely to play a significant role in forming the outer shell. The kinematic parameters of the inner shell make it possible to estimate the parameters of the supernova explosion by using the model from Wheeler et al. (1980). The SNR age is 106 yr, and the total explosion energy is 1052 erg. CONCLUSIONS Thus, in the simplest model, our observations of the neutral-hydrogen distribution around the SNR HB3 have shown that, first, we see here the results of the explosion of one of the most massive stars in the Galaxy, and, second, the shells from the stellar wind and from the SNR shock are observed simultaneously. This case seems extremely rare, although traces of the action of the wind and the SNR were found in several other objects, and the model used may not be unique. Note that the presence of a dense inner HI shell around HB3 was also pointed out by Routledge et al. (1991), although these authors could not determine the full range of its velocities. In addition, they also pointed out traces of the outer HI shell, but failed to get its complete picture due to the limited field of view. In general, our data may be considered to be consistent with the model of Routledge et al. (1991).

However, the quantitative estimates of the shell age and energetics should be considered approximate, bearing in mind the semiquantitative nature of the gas-dynamical models of Weaver et al. (1977) and Wheeler et al. (1980) used for this purpose. However, energetics of this scale is now not considered as something extraordinary. It is rather difficult to estimate the effect of the wind from the stars of the Cas OB6 association on the kinematics of the outer HI shell surrounding the SNR HB3, since the SNR, although being close to the association, is nevertheless outside it (see Fig. 1 from Routledge et al 1991). As regards the inner shell, it probably has a much more complex structure than that assumed above in the simple model of two expanding HI shells. Nevertheless, all of the authors consider it to be the result of the action of the SNR shock on the interstellar gas. In our drift curves, especially at the southern. . most declinations +6125 and +6188, we clearly see that the bulk of the HI is observed in the eastern and close (to the Sun) parts of the shell (the radial velocities are negative with respect to the mean velocity of the object). Incidentally, according to Huang and Thaddeus (1986), it is here that a CO cloud is located. In the western and far (from us) parts of the HB3 shell, where the SNR radio brightness is very low, Routledge et al. (1991) placed the highvelocity features of gas detected, in particular, by Lozinskaya and Sitnik (1980) in their model. According to Lozinskaya and Sitnik, the high-velocity features have a very low mass. As regards the bulk of the H-emitting filaments, their expansion velocity obtained by Lozinskaya and Sitnik (1980), 35 ▒ 25 km s-1 , agrees well with our data on the expansion of the HI shell (see the table). The significant decrease in the density of the cold neutral gas in the shell
ASTRONOMY LETTERS Vol. 31 No. 3 2005

INTERACTION OF THE SUPERNOVA REMNANT 3

185

around the SNR in its northwestern and farther (from the Sun) part is also confirmed by our data shown in Figs. 3н5. Probably, the regions of very hot gas responsible for the X-ray emission could be located there as well (Seward 1990). A strong nonuniformity of the distribution of interstellar gas in the region of the SNR HB3 and a significant amount of the molecular component dissociated by the shock could explain the rather strange mean initial densities that we obtained for the outer and inner shells. Finally, note that we presented here the results of RATAN-600 HI observations around one of the 104 SNRs selected from the catalog of Green (1998). The HI survey around them has now been completed, and the data are now processed. Later, drift curves in the нV plane similar to those presented here will be available on the Web page of the Special Astrophysical Observatory, Russian Academy of Sciences, for free use. ACKNOWLEDGMENTS I wish to thank Z.A. Alferova, G.N. Il'in, and T.M. Monastyreva from the Special Astrophysical Observatory, Russian Academy of Sciences, for maintaining the instrumentation in an operational state and help with the observations and data processing. This work was supported by the Russian Foundation for Basic Research (project no. 01-0217154). REFERENCES
1. Z. A. Alferova, I. V. Gosachinskii , S.R.Zhelenkov, et al., Izv. Spets. Astron. Obs. 23, 89 (1986).

2. N. A. Esepkina, N. S. Bakhvalov, B. A. Vasil'ev, et al., Izv. Spets. Astron. Obs. 11, 182 (1979). 3. D. A. Green, Astrophys. Space Sci. 148, 3 (1998). 4. Y.-L. Huang and P. Thaddeus, Astrophys. J. 309, 804 (1986). 5. G. N. Il'in, V. A. Prozorov, A. M. Pilipenko, Problems of Modern Radio astronomy, XXVII Radioastronomical Conference, Theses of Papers (Inst. PA. Ross. Akad. Nauk, St. Petersburg, 1997), p. 128. 6. F. G.Kerr and D. Linden Bell, MNRAS 221, 1023 (1986). 7. T. L. Landecker, J. F. Vaneldik, P. E. Dewdney, and D. Routledge, Astron. J.94, 111 (1987). 8. T. A. Lozinskaya, T. G. Sitnik, Astron. Zh. 57, 997 (1980) [Sov. Astron. 24, 572 (1980)]. 9. M. Normandeau, A. R. Taylor, and P. E. Dewdney, Astrophys. J. 108, 279 (1997). 10. D. Routledge, P. E. Dewdney, T. L. Landecker, et al., Astron. Astrophys. 247, 529 (1991). 11. F. D. Seward, Astrophys. J. Suppl. Ser. 73,781 (1990). 12. A. P. Venger, I.V.Gosachinskii ,V. G. Grachev, et al., Izv. Spets. Astron. Obs. 14, 118 (1981). 13. A. P. Venger, V. G. Grachev, T. M. Egorova, et al., Soobshch. Spets. Astron. Obs. 35, 5 (1982). 14. R. Weaver, R. McCray, J. Castor, et al., Astrophys. J. 218, 377 (1977). 15. J. C. Wheeler, T. J. Masurek, and A. Sivaramakrishnan Astrophys. J. 237, 781 (1980). 16. K. de Yager, The Most Luminous Stars (Mir, Moscow, 1984)[in Russian].

Translated by V. Astakhov

ASTRONOMY LETTERS

Vol. 31 No. 3 2005