Mapping the emission line strengths and kinematics of Supernova Remnant S147 with extensive LAMOST spectroscopic observations
We present extensive spectroscopic observations of supernova remnant (SNR) S147 collected with the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST). The spectra were carefully sky-subtracted taking into account the complex filamentary structure of S147. We have utilized all available LAMOST spectra toward S147, including sky and stellar spectra. By measuring the prominent optical emission lines including , N II 6584, and S II 6717, 6731, we present maps of radial velocity and line intensity ratio covering the whole nebula of S147 with unprecedented detail. The maps spatially correlated well with the complex filamentary structure of S147. For the central 2 of S147, the radial velocity varies from 100 to 100 and peaks between 0 and 10 . The intensity ratios of /S II 6717, 6731, S II 6717/ 6731 and /N II 6584 peak at about 0.77, 1.35 and 1.48, respectively, with a scatter of 0.17, 0.19 and 0.37, respectively. The intensity ratios are consistent with the literature values. However, the range of variations of line intensity ratios estimated here and representative of the whole nebula, are larger than previously estimated.
2018 Vol. X No. XX, 000–000
Supernova remnants (SNRs) are interesting objects to study for a number of reasons. They provide insights into the mechanisms of supernova explosions and are possible sources of the Galactic cosmic rays. Moreover, SNRs probe the immediate surroundings of supernovae, shaped by their progenitors. Currently 294 Galactic SNRs (Green 2014) have been discovered. A great majority ( 79%) of the known Galactic SNRs are believed to be relatively old (with ages greater than 10 years and radii larger than 5 pc) and evolved objects in either the adiabatic or the early radiative stages of their evolutionary developement (Woltjer 1972), and generally show characteristic shell structures (Fesen et al. 1985). About 30% of Galactic SNRs are known to have optical emission associated with their nonthermal radio emission. For old remnants, the optical emission arises from the cooling of shocked interstellar cloud material following the passage of the remnant’s blast waves as they expand outward into the ambient medium. Optical spectra of filaments of evolved SNRs show strong emission including that of , O II 3726,3729, O III 4959,5007, N II 6548,6584 and S II 6717,6731. The line intensity ratios have important physical implications, such as /S II 6717,6731 ratio can be used to distinguish shocked nebulae from photoionized nebulae; the S II 6717/S II 6731 ratio is electron density sensitive; the /N II 6584 ratio is a widely used tool for investigating nitrogen-to-hydrogen abundance variations among SNRs. Thus it important to obtain complete sampling of the optical line emissions in SNRs. Although several evolved remnants have been well studied morphologically in the optical, only a limited amount of spectroscopic data are available for those faint optical remnants, especially for those with large angular sizes including SNR S147.
S147 (also named G180.01.7) is an optically faint, highly filamentary shell-type SNR in the direction of Galactic-anticentre. It was first identified as a SNR candidate by Minkowski (1958). S147 is now believed to be one of the most evolved SNRs (with an estimated age 10 yr) in the Galaxy. Except for the spurs associated with a SNR, S147 has a shell-like structure of diameter 200. In the optical bands, the shell structure is dominated by filamentary H emission. The emission is bright at the northern and southern edges but concentrates mainly in the southern parts (Dinçel et al. 2015). Despite of its large age, it conserves well its spherical symmetry except for the blowout regions in the east and west. The optical and radio shell morphologies are well defined (Minkowski 1958; Fuerst & Reich 1986), and coincide with each other. Fig. 1 shows images of S147 in different bands. Although no X-ray emission has been detected by EXOSAT (Sauvageot et al. 1990), Sun et al. (1996) report unambiguous X-ray emission from S147 based on data from the ROSAT All Sky Survey (RASS, Voges et al. 1999). Chen et al. (2017) claim that the X-ray emission and the dust extinction is anti-correlated. A spatially extended gamma-ray source detected in the energy range 0.2 – 10 GeV is found to coincide with S147. The gamma-ray emission exhibits a possible spatial correlation with the prominent filaments of S147 (Katsuta et al. 2012).
A compact object, the radio pulsar J0538+2817 (Anderson et al. 1996), located 40 west of the centre of S147, is believed to be associated with S147. An early type (B0.5V) runaway star, HD 37424, has also been found inside S147, believed to be the pre-supernova binary companion to the progenitor of the pulsar and the SNR (Dinçel et al. 2015). The distance estimated of S147 in the literature spans a wide range, 0.6 – 1.9 kpc, depending on the methods used (Dinçel et al. 2015). The surface brightness and distance (–D) relation typically yields a distance of 1 kpc from several studies (Clark & Caswell 1976; Kundu et al. 1980; Guseinov et al. 2003). The interstellar reddening toward S147 suggests a smaller distance of 0.8 kpc (Fesen et al. 1985). On the other hand, distances derived from the parallax of the pulsar are significantly larger, ranging between 1.3 – 1.5 kpc (Chatterjee et al. 2009; Ng et al. 2007). A recent study yields distance kpc and extinction mag (Dinçel et al. 2015). Finally, from a dust extinction analysis based on data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-centre (XSTPS-GAC; Liu et al. 2014; Zhang et al. 2014), Chen et al. (2017) obtain a new measurement of distance to S147, kpc. The expansion velocity estimated for the S147 shell ranges between 80 – 120 (Kirshner & Arnold 1979; Phillips et al. 1981).
Previous optical spectrophotometric observations of S147 (Kirshner & Arnold 1979; Fesen et al. 1985) obtain only dozens of spectra, covering a small portion of the extremely large and faint S147 SNR. The observations focused on several bright filaments. S147 consists of numerous filaments embedded in diffuse emission. The two components (i.e., bright filaments and diffuse emission) have quite different physical properties. It is thus important to sample fully the whole extent of S147 in order to obtain a comprehensive understanding of an evolved SNR such as S147.
The Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) is a quasi-meridian reflecting Schmidt telescope of 4-m effective aperture, with a field of view of 5 diameter (Cui et al. 2012). It is located in Xinglong station of National Astronomical Observatories of Chinese Academy of Sciences. Being a dedicated survey telescope, LAMOST takes spectra with 4 000 fibers, targeting interested celestial objects including sky background and calibration sources in one single exposure with 16 fiber-fed spectrographs, each accommodates 250 fibers. The spectra cover the entire optical wavelength range ( 3700 – 9000Å) at a resolving power 1800. Since 2012 September, LAMOST has been carrying out a five-year Regular Surveys. Before that there was a one-year Pilot Surveys preceded by a two-year commissioning phase. The LAMOST Regular Surveys consist of two components (Zhao et al. 2012): the LAMOST Extra-Galactic Survey of galaxies (LEGAS) that aims at studying the large scale structure of the universe, and the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) that aims at obtaining millions of stellar spectra in order to study the structure and evolution of the Milky Way galaxy (Deng et al. 2012).
The LEGUE has three sub-components: the spheroid, the anticentre and the disc surveys. The LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC; Liu et al. 2014; Yuan et al. 2015; Xiang et al. 2017) aims to survey a significant volume of the Galactic thin/thick discs and halo for a contiguous sky area of over 3400 deg centred on the Galactic anticentre (, 150 210), down to a limiting magnitude of 17.8 mag (to 18.5 mag for limited fields). S147 is located within the LSS-GAC footprint. A large number of foreground and background stars toward S147 have been targeted by LAMOST accidentally, some by several times. S147 has a large diameter and contains many filaments. It would be extremely time consuming if one attempts to obtain spectra covering its full spatial extent using regular telescopes. The wide field of view, multi-fiber survey telescope LAMOST thus serves as a perfect tool for mapping the details of this extremely extended, evolved SNR. It provides us a great opportunity to study global properties of S147 with unprecedented detail for the first time.
In this paper, we focus on the analysis of the LAMOST optical spectra of S147. In particular, we concentrate on the kinematic properties revealed by strong emission lines including , N II 6548,6584 and S II 6717,6731. In Section 2 we describe the observation and data reduction. In Section 3 we present parameter determinations and result analysis. Section 4 summarize our findings.
2 Observation and Data Reduction
The LAMOST Data Release 2 (DR2) includes data from 1 934 plates collected by June 2014, amongst them 401 plates are from the Pilot Surveys (from October 2011 to June 2012), 808 plates from the first year Regular Surveys (from September 2012 to June 2013), and 725 plates from the second year Regular Surveys (from September 2013 to June 2014). Totally, the DR2 general catalogue contains 4 132 782 spectra, including 3 779 674 stellar, 37 665 of galaxies, 8 633 of QSOs, and 306 810 of unknown objects. The LAMOST DR2 data was released internally on December 2014 and publicly on June 2016 111http://dr2.lamost.org/.
Being a large-size SNR in the direction of Galactic anticentre, many foreground and background stars in the vicinity region of S147 have been targeted, some by several times by LAMOST since 2011. There are 50 LAMOST plates overlapping with (or very close to) S147. Fig. 2 plots the positions of these 50 plates included in the LAMOST DR2, overlaid with an image from the IPHAS survey (Drew et al. 2005).
Amongst the 50 LAMOST plates, 16 have a centre of field of view within 3.55 of the geometric centre of S147 (10 B plates, and 6 M plates 222LAMOST plates are grouped into four categories based on the magnitude ranges of the targets: Very Bright (VB; 14 mag), Bright (B; 14 16.8 mag), Medium-bright (M; 16.8 17.8 mag), and Faint (F; 17.8 mag) (Luo et al. 2015; Liu et al. 2014; Yuan et al. 2015).). Of those 16 plates, 9 are from the Pilot Surveys, 7 from the first year Regular Surveys. The other 34 plates have centres of field of view between 3.55 and 4.54 from the geometric centre of S147 (5 VB plates, 26 B plates, 16 M plates and 3 F plates). Amongst those 34 plates, 28 are from the Pilot Surveys, 18 from the first year Regular Surveys, and 4 from the second year Regular Surveys.
2.2 Data Reduction
LAMOST raw spectra are processed with the LAMOST two-dimensional (2D) pipeline (Luo et al. 2012, 2015), including dark and bias subtractions, cosmic ray removal, one-dimensional (1D) spectral extraction, flat-field correction, wavelength calibration, sky subtraction, merging sub-exposures, and finally, splicing the sub-spectra from the blue- and red-arm of the spectrograph, respectively. Xiang et al. (2015a) develop an iterative algorithm for the flux calibration in order to overcome the effect of often heavy interstellar reddening that hinders the selection of suitable standard stars (F stars). For a given spectrograph, the spectra are first flux-calibrated using the nominal spectral response curve (SRC) and the initial stellar atmospheric parameters are derived with the LAMOST Stellar Parameter Pipeline at Peking University (LSP3; Xiang et al. 2015b). Then based on those initial estimates of stellar parameters, stars of effective temperature between 5750 and 6750 K (mostly of F spectral type) are selected as the flux-calibration standard stars and used to deduce an updated SRC after corrected for the interstellar reddening, derived by comparing the observed and synthetic photometric colours, assuming the Fitzpatrick (1999) extinction law for a total to selective extinction ratio = 3.1. The new SRC is then used to flux-calibrate the spectra and the stellar parameters are updated accordingly from the recalibrated spectra. The above process is repeated until a convergence is achieved. Usually more than four standard stars per spectrograph can be identified and selected for the majority of the LSS-GAC plates. In cases where not enough standard stars can be selected for a given spectrograph (plate), the SRC(s) deduced from other plates, usually collected on the same night, are used to do the calibration. A comparison of spectra obtained at different epochs for duplicate targets indicates that the relative flux calibration has achieved an accuracy of about 10 per cent for spectra of a signal to noise (SN) ratio per pixel higher than 30 at 4650 Å (Xiang et al. 2015a).
2.2.2 Data toward S147
In total, there are 19 736 spectra (including those from the sky fibers) within 2 of the geometric centre of S147, of which 4 667 spectra (4 349 unique stars) have available LSP3 stellar atmospheric parameters. Fig. 3 shows the histogram distributions of spectral SN ratio in the red (7450 Å) regions of the 19 736 spectra within 2 of S147.
2.2.3 Sky subtraction
For each LAMOST plate, 10 – 15 per cent of the fibers are assigned as object-free sky fibers pointing toward areas of blank sky for modeling the night sky background. A method of B-spline fitting (Stoughton et al. 2002; Luo et al. 2015) is applied to make a supersky for each sub-exposure of a spectrograph and then the principal component analysis (PCA) is used to recorrect the supersky. The fiber assignment strategy works well for fields at high Galactic latitudes, yielding results comparable to those of the SDSS survey. However, for sky areas near the Galactic plane, such as of S147 where the sky is plagued by diffuse and highly filamentary nebular emission, the strategy fails to build a supersky of adequate accuracy. In the current work, we develop a new method for the sky subtraction, which is described below.
Overall S147 has a quite regular and nearly spherical shape, with a diameter of 3, smaller than the field of view of LAMOST (5 in diameter). For the individual LAMOST plates taken in the vicinity field of S147, there are always some regions not contaminated by the emission of S147. Thus for those plates, we are able to use the sky fibers positioned at points uncontaminated by S147 emission to rebuild the supersky, and redo the sky subtraction. This procedure is carried out for the blue- and red-arm spectra separately. The uncontaminated sky spectra from the non-S147 region are selected following the criteria: 1) At least 2 away from the S147 geometric centre; 2) Not contaminated by S147 emission, i.e., those with a detectable H flux (larger than 2, for strictly cut) are also excluded. The sky spectra thus selected and satisfying the above criteria are then averaged to construct the so-called “uncontaminated non-S147 supersky”.
LAMOST has 16 spectrographs, each feeds by 250 fibers, where 10 to 20 of them are normally assigned as sky fibers. For the individual spectrographs, the “uncontaminated non-S147 supersky” are scaled to match each sky spectrum. The scaling procedure includes two steps, scaling by the continuum and by the telluric emission lines. The continuum scaling is carried out first, and then the line scaling. In the current work intensities of the sky emission lines 6498, 6533 in red-arm are used for the line scaling. We then obtain a “new initial supersky” for each of the 16 spectrographs (abnormal sky spectra are excluded during this process). We note that continuum and line scaling works much better for a limited wavelength region than for the whole spectrum. In this paper, we focus only on selected wavelength regions of interest here: 6400 – 6800 Å in the red that contains , and the N II 6548,6584 and S II 6717,6731 nebular lines. The sky subtraction is only carried out for the two selected wavelength regions. Once the “new initial supersky” has been obtained, segments of this spectrum are replaced by the “uncontaminated non-S147 supersky” for the following wavelength ranges: 6490 – 6620 Å and 6695 – 6750 Å . The supersky spectra thus obtained are then used for the final sky subtraction.
For an uncontaminated spectrograph, sky fibers fall outside the S147 emission region and are thus free of S147 emission. The final adopted supersky differs little from the supersky yielded by the LAMOST 2D pipeline. For uncontaminated spectrographs, the reconstructed new supersky agrees well with the original LAMOST supersky, indicating that the LAMOST default supersky works well for uncontaminated sky regions. For contaminated spectrographs, the differences between our new adopted supersky and that given by the LAMOST 2D pipeline can be significant, especially for spectral regions affected by prominent nebular emissions such as , N II 6548,6584, and S II 6717,6731. Our sky subtraction procedure has a typical uncertainty of approximately 10 per cent. After the sky subtraction, nebular features, especially , N II 6548,6584, and S II 6717,6731 lines become clearly visible.
|20111215-GAC082N27B1-7-145||b||05:59:44.65||26:17:10.8||31.2||46.7||36.6||12.6||27.1||N II 6584|
|20111215-GAC082N27B1-6-109||b||05:59:43.78||26:23:31.0||13.9||12.5||0.8||12.6||102.7||N II 6584|
|20111215-GAC082N27B1-7-134||b||05:59:30.64||26:07:59.3||12.4||19.9||21.5||11.3||163.9||N II 6584|
|20111126-GAC082N27M1-7-140||a||05:59:25.41||25:59:41.9||2.9||0.9||-17.2||11.4||42.6||N II 6584|
|20111215-GAC082N27B1-7-139||b||05:59:23.28||26:14:29.6||14.6||20.7||11.6||11.7||169.6||N II 6584|
|20111126-GAC082N27M1-7-94||a||05:59:03.98||25:59:41.9||8.1||2.8||-28.3||13.3||33.7||N II 6584|
|20111126-GAC082N27M1-7-135||b||05:58:43.02||26:07:58.6||11.6||29.4||8.1||11.3||99.4||N II 6584|
|20111215-GAC082N27B1-7-135||b||05:58:42.08||26:11:06.4||18.1||35.1||36.2||10.8||185.9||N II 6584|
|20111215-GAC082N27B1-7-181||b||05:58:24.94||26:11:54.3||33.5||71.1||22.4||15.1||104.8||N II 6584|
|20111215-GAC082N27B1-7-186||b||05:58:12.60||26:07:15.9||16.5||28.4||28.5||11.4||127.1||N II 6584|
3 Parameter determinations
To investigate the kinematic properties of S147, we determine by Gaussian fitting the radial velocities and line intensities of the prominent emission lines including , N II 6548, 6584, and S II 6717, 6731. Only nebular lines in the red-arm spectra are analyzed, given that they are much more stronger than those in the blue. Usually, the red-arm spectra have higher SN ratio than the blue-arm ones.
There are typically 300 sky spectra per plate. To utilize as many spectra as possible, all spectra of foreground and background stars in the S147 vicinity field are analyzed as well. A small portion of those stellar spectra have good SN ratio, allowing stellar parameter determinations. For the remaining ones, because of the low SN ratio, no stellar parameters are available from LSP3 (as shown in Fig. 3). Since the stellar features may affect the measurements of nebular emission lines of S147, they need to be subtracted before determining the properties of nebular emission lines.
For high quality stellar spectra with available stellar atmospheric parameters determined with LSP3, the subtraction is done using high quality spectra of “paired stars” that have almost the same stellar atmospheric parameters in regions uncontaminated by nebular emission (Yuan et al. 2013). The “paired stars” are selected from the LSS-GAC Value-added Catalogues DR1 (Yuan et al. 2015) following the criteria: 1) Located at high Galactic latitudes ( 20); 2) Suffer from low extinction (ebvsp 0.05 mag, where ebvsp is the derived from the star pair method, see Yuan et al. (2015) for details); 3) Have high SN ratio spectra (S/N 20 or S/N 20; In a few cases where no matching stars are found, this is relaxed to S/N 10 or S/N 10, where S/N and S/N represent the signal to noise ratio in and band respectively); 4) Have stellar atmospheric parameters of small uncertainties (effective temperature error tefferr 150 K, surface gravity error loggerr 0.25 dex, metallicity error feherr 0.15 dex); 5) Have stellar atmospheric parameters that differ from the targeted stars by small amounts (specifically, the effective temperature difference 150 K, surface gravity difference 0.25 dex, metallicity difference 0.15 dex). For a given target star, the spectrum of a “paired star” that satisfies all above criteria and has the smallest value of is used to subtract the stellar spectrum from the target spectrum. Fig. 4 gives an example of stellar spectrum subtraction.
Most of those low SN ratio stellar spectra with no stellar atmospheric parameters available from LSP3 are of faint stars. The spectra are thus only marginally affected by stellar features. For those spectra, we have ignored the possible contamination of stellar features and measure directly properties of the nebular emission lines without stellar spectrum subtraction. Fig. 5 shows an example of a low SN ratio stellar spectrum. Although the (stellar) continuum is quite noisy and has flux close to zero, the nebular emission lines are quite prominent.
We determine radial velocities and intensities of , N II 6548, 6584, and S II 6717, 6731 nebular emission lines from all available spectra in the vicinity field of S147. Part of the results are given in Table 4. The full Table can be found online. Fig. 6 shows example red-arm spectra for several selected filaments.
4.1 The radial velocity fields
As shown in Fig. 6, the N II 6548 line is in general rather weak and seen only in a few high quality spectra. is the strongest nebular line and visible in almost all the spectra. The N II 6584 line has a typical strength comparable to the S II 6717, 6731 lines. The latter three lines have comparable occurrence rates, and are often less stronger than . For spectra of comparable SN ratio, the stronger the nebular lines are, the higher the detection rates will be, and vice verse. In the current work, only , N II 6584, and S II 6717, 6731 lines are used to measure the kinematic properties of S147.
Fig. 7 shows the radial velocity field of S147 as revealed by the emission lines measured. One can see that radial velocity fields mapped by emission lines exhibit spatial features well correlated. The radial velocities of different filaments of S147 show complex structure and variations. Some of the filaments are moving towards from us (such as the segment marked by blue ellipse in Fig. 7), some are moving away from us (such as the segment marked by red ellipse in Fig. 7). Probably due to the projection effects or a direct consequence of the complex structure of the filaments, many segments of filaments of blue- and red-shifted velocities are tangled together. One of such examples is marked by green ellipse in Fig. 7. The radial velocity fields presented here cover the whole extent of S147 with unprecedented sampling. Fig. 8 shows histogram distributions of all derived radial velocities of uncertainties less than 20 determined from various lines for spectra within 2 of the geometric centre of S147. The distributions peak between 0 – 10 .
4.2 Maps of line intensity ratios
Fig. 9 shows the distributions of line intensity ratio /N II 6584, /S II 6717, 6731 and S II 6717/S II 6731. Similar to the distributions of radial velocities (Fig. 7), the line intensity ratios are well measured along prominent filaments. Fig. 10 shows the histogram distributions of the above line intensity ratios, derived from spectra within 2 of the centre of S147. The ratios /S II 6717, 6731 and S II 6717/S II 6731 peak at 0.77 (dispersion = 0.17) and 1.35 ( = 0.19; i.e. close to the low electron density limit), respectively. Our results are consistent with previous measurements (/S II 6717, 6731 ranges between 0.7 – 1.08; S II 6717/S II 6731 1.4) in the literature (Fesen et al. 1985; Lozinskaia 1976). Since the N II 6548 line is weak and only detected in a very limited sample of spectra, we use only the N II 6584 line to calculate the /N II ratio, i.e. /N II 6584. The peak value of /N II 6584 is 1.48 ( = 0.37), again in good agreement with the literature values of 1.20 – 1.58 (D’Odorico & Sabbadin 1977). Figs.11 and 12 show example spectra of filaments with high and low line intensity ratios and of strong and weak line intensities, respectively.
Generally SNRs differ from H II regions by showing stronger S II 6717, 6731 emission. Quantitatively, /S II 6717, 6731 2.0 for SNRs (Fesen et al. 1985). The middle panels of Fig. 9 and Fig. 10 show that in S147, the /S II 6717, 6731 ratio has values less than 1.4, consistent with S147 being a typical SNR.
It has been suggested that the intensity ratios /N II 6548,6584, /S II 6717, 6731 and S II 6717/S II 6731 vary less from region to region of a given SNR than they do amongst different SNRs (Daltabuit et al. 1976). Based on this, Daltabuit et al. (1976) develop an evolutionary scheme of SNRs in which /N II 6548,6584 and S II 6717/S II 6731 vary systematically as a function of the SNR diameter. In addition, the /N II 6548,6584 ratio combined with electron density indicated by the S II 6717/S II 6731 ratio has been used to infer the elemental abundance gradient in the Milky Way, M31, and M33 (Fesen et al. 1985; Binette et al. 1982; Dopita et al. 1980; Blair et al. 1981, 1982; Blair & Kirshner 1985). This approach is only valid if the aforementioned line ratio does not vary significantly in a given SNR. Only a very limited filaments of Galactic SNRs have been studied spectroscopically. The extensive data presented in the current work provide an opportunity to investigate the intrinsic variations of line intensity ratios in a given SNR and test the legitimacy of the aforementioned approach.
From Fig. 9 and 10, we estimate that the intrinsic dispersions of line ratios amount to be approximately 14% for S II 6717/S II 6731, 22% for /S II 6717, 6731, and 25% for /N II 6584. The dispersions are larger than those previously estimated by Fesen et al. (1985) based on only five observations (the mean values of /N II 6548,6584, /S II 6717, 6731 and S II 6717/S II 6731 are respectively 1.28, 1.00, 1.61, with corresponding dispersions of 11%, 8%, and 12%) and those estimated by Dodorico & Sabbadin (1976) based on only six observations (the mean values of /N II 6548,6584, /S II 6717, 6731 and S II 6717/S II 6731 are respectively 1.44, 0.79, 1.35, with corresponding dispersions of 9%, 10%, and 6%). On the other hand, the results presented here are comparable to the dispersions estimated for 9 typical SNRs as listed in Table 5 of Fesen et al. (1985). The large intrinsic dispersions of line intensity ratios seen in S147 cast doubt on the viability of using those ratios to infer the evolutionary stage of a SNR or to estimate the Galactic elemental abundance gradient using SNRs (see also Fesen et al. 1985).
Being an electron density diagnostic, the S II 6717/S II 6731 ratio has a mean value of 1.35 for S147, i.e. close to the low density limit of 1.4 (see Fig 6. in Blair & Kirshner (1985) for the S II 6717/S II 6731 ratio as a function of electron density, the electron density goes to the lowest when the S II 6717/S II 6731 ratio has a value of 1.4). The result is consistent with typical values found for evolved SNRs (larger than 1.10, thus implies an electron density lower than 300 cm; Fesen et al. 1985). We estimated that S II 6717/S II 6731 has an intrinsic dispersion of 14% in S147, slightly larger than those found previously by Fesen et al. (1985) and Dodorico & Sabbadin (1976) (12% and 6%, respectively). The result suggests a small range of electron densities and cloud pressures in S147. A relation of electron density and the S II 6717 / S II 6731 line ratio has previously been presented by Blair & Kirshner (1985) (see Blair & Kirshner 1985, Fig. 6) assuming a five level atomic model, the collisional strenghths of Pradhan (1978) and the transition probabilities of Mendoza & Zeippen (1982). From the plot of Blair & Kirshner (1985), the mean value of S II 6717/S II 6731 ratio seen in S147, 1.35 along with a standard deviation 0.19, implies an the electron density of only 50 cm, with an upper limit of about 200 cm.
4.3 Comparison with previous spectroscopic observations
Kirshner & Arnold (1979) observed 15 positions of S147 (see Table 1 in their paper) and measured their radial velocities. Their Fig. 1 shows the slit positions used. The exact coordinates of the slit positions are however not given in the paper. So it is hard to make an exact comparison. Referring to their Fig. 1, we select some of their slit positions that overlap approximately with (or close to) some of the fiber positions measured with LAMOST. Table 2 compares the radial velocities determined by Kirshner & Arnold (1979) and by us. Although the differences between our measurements and literature values for position ‘c’ and ‘k’ seem large ( 30 ), we need to remember that systematic uncertainties of radial velocity is 10 (Luo et al. 2015) and the comparison here is not exact for the spectra obtained in same position. So finally we conclude that the agreement is reasonable.
|Position||Radial Velocity||Radial Velocity|
|(Kirshner & Arnold 1979)||(here)|
|Position||Distance||/S II||/S II|
|(arcmin)||(Fesen et al. 1985)||(here)|
Fesen et al. (1985) carried out spectroscopic observations of five positions of S147. Two of them are close to the filaments also observed with LAMOST (P2 and P5 in Fig.13). Table 3 compares the /S II line intensity ratios as determined by Fesen et al. (1985) and by us for P2 and P5. Because the exact positions of their observations are about 4 arcmin away from our observations, the small differences between the two sets of measurements seem to be acceptable.
We have analyzed all spectra available in the LAMOST DR2 collected in the vicinity field of SNR S147. The spectra are carefully sky subtracted. Both the spectra of sky background and of foreground and background stars are used. For high-quality stellar spectra with atmospheric parameters determined, the underlying stellar features are subtracted before measuring the nebular emission line properties. For stellar spectra of low SN ratio and with no stellar atmospheric parameters determined, the spectra are directly used to measure the nebular line properties.
By measuring the prominent emission lines , N II 6584, and S II 6717, 6731, we obtain fundamental kinematic properties of S147, including maps of radial velocity and intensity ratio of exquisite detail. As it is not a dedicated filamentary survey of LAMOST in S147 region, the number of filamentary spectra is very low comparing with the total spectral number in this region. But still this is first such as a large SNR has been fully mapped spectroscopically.
The radial velocity distribution peaks 0 – 10 . The intensity ratios /S II 6717, 6731 and S II 6717/S II 6731 peak respectively at 0.77 ( = 0.17) and 1.35 ( = 0.19), consistent with previous determinations in the literature. The /N II 6584 ratios peak at 1.48 ( = 0.37), and agrees well with the literature value ( 1.20 – 1.58). The intrinsic dispersions of line intensity ratios estimated here are larger than previous estimates in the literature, based on observations of only a few filaments.
Acknowledgements.This work is supported by Joint Funds of the National Natural Science Foundation of China (Grant No. U1531244) and the National Key Basic Research Program of China 2014CB845700. JJR acknowledges support from the Young Researcher Grant of National Astronomical Observatories, Chinese Academy of Sciences. JJR thanks Albert Zijlstra, Nick Wright, and Janet Drew for kindly providing the IPHAS image of S147. The LAMOST FELLOWSHIP is supported by Special Funding for Advanced Users, budgeted and administrated by Centre for Astronomical Mega-Science, Chinese Academy of Sciences (CAMS). This work has made use of data products from the Guoshoujing Telescope (the Large Sky Area Multi-Object Fibre Spectroscopic Telescope, LAMOST). LAMOST is a National Major Scientific Project built by the Chinese Academy of Sciences. Funding for the project has been provided by the National Development and Reform Commission. LAMOST is operated and managed by the National Astronomical Observatories, Chinese Academy of Sciences.
- Anderson et al. (1996) Anderson, S. B., Cadwell, B. J., Jacoby, B. A., et al. 1996, \apjl, 468, L55
- Binette et al. (1982) Binette, L., Dopita, M. A., Dodorico, S., & Benvenuti, P. 1982, \aap, 115, 315
- Blair & Kirshner (1985) Blair, W. P., & Kirshner, R. P. 1985, \apj, 289, 582
- Blair et al. (1981) Blair, W. P., Kirshner, R. P., & Chevalier, R. A. 1981, \apj, 247, 879
- Blair et al. (1982) Blair, W. P., Kirshner, R. P., & Chevalier, R. A. 1982, \apj, 254, 50
- Chatterjee et al. (2009) Chatterjee, S., Brisken, W. F., Vlemmings, W. H. T., et al. 2009, \apj, 698, 250
- Chen et al. (2017) Chen, B.-Q., Liu, X.-W., Ren, J.-J., et al. 2017, \mnras, 472, 3924
- Clark & Caswell (1976) Clark, D. H., & Caswell, J. L. 1976, \mnras, 174, 267
- Cui et al. (2012) Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Research in Astronomy and Astrophysics, 12, 1197
- Daltabuit et al. (1976) Daltabuit, E., Dodorico, S., & Sabbadin, F. 1976, \aap, 52, 93
- Deng et al. (2012) Deng, L.-C., Newberg, H. J., Liu, C., et al. 2012, Research in Astronomy and Astrophysics, 12, 735
- Dinçel et al. (2015) Dinçel, B., Neuhäuser, R., Yerli, S. K., et al. 2015, \mnras, 448, 3196
- Dodorico & Sabbadin (1976) Dodorico, S., & Sabbadin, F. 1976, \aap, 53, 443
- D’Odorico & Sabbadin (1977) D’Odorico, S., & Sabbadin, F. 1977, \aaps, 28, 439
- Dopita et al. (1980) Dopita, M. A., Dodorico, S., & Benvenuti, P. 1980, \apj, 236, 628
- Drew et al. (2005) Drew, J. E., Greimel, R., Irwin, M. J., et al. 2005, \mnras, 362, 753
- Fesen et al. (1985) Fesen, R. A., Blair, W. P., & Kirshner, R. P. 1985, \apj, 292, 29
- Fitzpatrick (1999) Fitzpatrick, E. L. 1999, \pasp, 111, 63
- Fuerst & Reich (1986) Fuerst, E., & Reich, W. 1986, \aap, 163, 185
- Green (2014) Green, D. A. 2014, Bulletin of the Astronomical Society of India, 42, 47
- Guseinov et al. (2003) Guseinov, O. H., Ankay, A., Sezer, A., & Tagieva, S. O. 2003, Astronomical and Astrophysical Transactions, 22, 273
- Katsuta et al. (2012) Katsuta, J., Uchiyama, Y., Tanaka, T., et al. 2012, \apj, 752, 135
- Kirshner & Arnold (1979) Kirshner, R. P., & Arnold, C. N. 1979, \apj, 229, 147
- Kundu et al. (1980) Kundu, M. R., Angerhofer, P. E., Fuerst, E., & Hirth, W. 1980, \aap, 92, 225
- Liu et al. (2014) Liu, X.-W., Yuan, H.-B., Huo, Z.-Y., et al. 2014, in Setting the scene for Gaia and LAMOST, IAU Symposium, vol. 298, edited by S. Feltzing, G. Zhao, N. A. Walton, & P. Whitelock, 310–321
- Lozinskaia (1976) Lozinskaia, T. A. 1976, \azh, 53, 38
- Luo et al. (2012) Luo, A.-L., Zhang, H.-T., Zhao, Y.-H., et al. 2012, Research in Astronomy and Astrophysics, 12, 1243
- Luo et al. (2015) Luo, A.-L., Zhao, Y.-H., Zhao, G., et al. 2015, Research in Astronomy and Astrophysics, 15, 1095
- Mendoza & Zeippen (1982) Mendoza, C., & Zeippen, C. J. 1982, \mnras, 198, 127
- Minkowski (1958) Minkowski, R. 1958, Reviews of Modern Physics, 30, 1048
- Ng et al. (2007) Ng, C.-Y., Romani, R. W., Brisken, W. F., Chatterjee, S., & Kramer, M. 2007, \apj, 654, 487
- Phillips et al. (1981) Phillips, A. P., Gondhalekar, P. M., & Blades, J. C. 1981, \mnras, 195, 485
- Pradhan (1978) Pradhan, A. K. 1978, \mnras, 183, 89P
- Sauvageot et al. (1990) Sauvageot, J. L., Ballet, J., & Rothenflug, R. 1990, \aap, 227, 183
- Stoughton et al. (2002) Stoughton, C., Lupton, R. H., Bernardi, M., et al. 2002, \aj, 123, 485
- Sun et al. (1996) Sun, X., Anderson, S., Aschenbach, B., et al. 1996, in Roentgenstrahlung from the Universe, edited by H. U. Zimmermann, J. Trümper, & H. Yorke, 195–196
- Voges et al. (1999) Voges, W., Aschenbach, B., Boller, T., et al. 1999, \aap, 349, 389
- Woltjer (1972) Woltjer, L. 1972, \araa, 10, 129
- Xiang et al. (2015a) Xiang, M. S., Liu, X. W., Yuan, H. B., et al. 2015a, \mnras, 448, 90
- Xiang et al. (2015b) Xiang, M. S., Liu, X. W., Yuan, H. B., et al. 2015b, \mnras, 448, 822
- Xiang et al. (2017) Xiang, M.-S., Liu, X.-W., Yuan, H.-B., et al. 2017, \mnras, 467, 1890
- Xiao et al. (2008) Xiao, L., Fürst, E., Reich, W., & Han, J. L. 2008, \aap, 482, 783
- Yuan et al. (2015) Yuan, H.-B., Liu, X.-W., Huo, Z.-Y., et al. 2015, \mnras, 448, 855
- Yuan et al. (2013) Yuan, H. B., Liu, X. W., & Xiang, M. S. 2013, \mnras, 430, 2188
- Zhang et al. (2014) Zhang, H.-H., Liu, X.-W., Yuan, H.-B., et al. 2014, Research in Astronomy and Astrophysics, 14, 456-470
- Zhao et al. (2012) Zhao, G., Zhao, Y.-H., Chu, Y.-Q., Jing, Y.-P., & Deng, L.-C. 2012, Research in Astronomy and Astrophysics, 12, 723