J. L. Payne1 - M. D. Filipovic2,3 - T. G. Pannuti4 - P. A. Jones3 - N. Duric5 - G. L. White1 - S. Carpano6
1 - Centre for Astronomy, James Cook University, Townsville,
QLD, 4811, Australia
2 - University of Western Sydney, Locked Bag 1797, Penrith South,
DC, NSW, 1797, Australia
3 - Australia Telescope National Facility, CSIRO, PO Box 76,
Epping, NSW, 1710, Australia
4 - Spitzer Science Center, California Institute of Technology,
Mailstop 220-6, Pasadena, CA 91125, USA
5 - Department of Physics and Astronomy, University of New
Mexico, 800 Yale Bd. N.E., Albuquerque, NM 87131, USA
6 - Institut für Astronomie und Astrophysik, Universität
Tübingen, Sand 1, 72076 Tübingen, Germany
Received 8 April 2004 / Accepted 4 June 2004
Abstract
We present a multi-frequency study of supernova remnants
(SNRs) and H II regions in the nearby Sculptor Group Sd galaxy
NGC 300, based on new ATCA observations at the wavelengths of 13
and 20 cm, XMM-Newton observations, newly-processed ROSAT (PSPC/HRI; Read & Pietsch 2001) and VLA (20/6 cm) images of
this galaxy. We have investigated the physical properties at the
X-ray and radio wavelengths of the 28 optical SNRs found by
Blair & Long (1997) and have expanded on the multi-wavelength work by
Pannuti et al. (2000) on this same galaxy. From a total of 54 radio
sources and 11 X-ray sources, we report 18 SNRs and three (3) SNR candidates (classified by spectral index alone) in
NGC 300. Five of these 18 SNRs are associated with reported
optical SNRs and three have X-ray counterparts. An
additional 12 radio SNRs are seen in the Blair & Long (1997) [S
II] images. We also investigate the luminosity function of our
SNRs. Three background radio sources are confirmed and 12 other
sources could represent additional background objects. Twenty two
radio correlations with OB associations within NGC 300 correspond
to either H II regions or SNRs making them a good tracer of SNRs
near star-forming regions. Additionally, two radio sources
coincide with potential globular clusters of NGC 300 reported by
Kim et al. (2002).
Key words: galaxies: individual: NGC 300 - ISM: supernova remnants - radio continuum: galaxies - X-rays: galaxies - ISM: H II regions
To address these difficulties and to garner more insights into the nature of SNRs, observers have searched for SNRs in nearby galaxies such as the Large Magellanic Cloud (LMC) (Filipovic et al., in preparation; Williams et al. 1999) , the Small Magellanic Cloud (SMC) (Filipovic et al., in preparation), M 31 (Braun & Walterbos 1993), M 33 (Gordon et al. 1998), and the Sculptor Group galaxies NGC 300 and NGC 7793 (Blair & Long 1997, hereafter referred to as BL97; and Pannuti et al. 2002). Most of these surveys have concentrated on optical methods to detect SNRs and have met with considerable success in revealing a large number of new sources. However, such optical observations need to be complemented by observations at other wavelengths associated with emission from SNRs (such as X-ray and radio) to both detect a maximum number of SNRs in a galaxy of interest and to locate SNRs at all stages of their evolution.
As pointed out in Lacey & Duric (2001), there may be a selection effect inherent in optical surveys, which are more sensitive to detecting SNRs away from H II regions. Such SNRs may be weak emitters in the X-ray and radio because of their location within low density environments which are not conducive to high X-ray and radio luminosities. In contrast, SNRs that are powerful sources of X-ray and radio emission are often located within H II regions. It is precisely because of their location in H II regions that they may be missed by optical surveys.
Although there are certainly exceptions, it can be thought that SNRs located in regions of low density were parented by low-mass white dwarf progenitor stars as Type Ia supernovae, while in contrast SNRs that are found in regions of high density were created by high-mass progenitor stars such as Type Ib/Ic/II supernovae. Thus, multi-wavelength studies also give a more complete picture of a galaxy's star formation history (Pannuti et al. 2002).
In this paper, we present results of a multi-wavelength search for SNRs in NGC 300 using new radio and X-ray data. Table 1 lists some of the more important properties of this galaxy. In Sect. 2, we describe our observations and data analysis. We discuss and explain the method used to identify SNRs, H II regions and background sources in Sect. 3. Position and flux density analysis is presented in Sect. 4, and a description of individual sources is given in Sect. 5. Finally, Sect. 6 is a summary and concluding discussion.
Table 1: Important properties of NGC 300. (NED = NASA/IPAC Extragalactic Database).
![]() |
Figure 1:
ATCA 1374 MHz contours overlaying an DSS2-Red (R) image
(grey scale) of NGC 300. An ellipse (21.9 |
| Open with DEXTER | |
![]() |
Figure 2:
The ATCA 2496 MHz image of NGC 300 bounded by our ellipse with its contours.
Flux from two nearby contaminating sources,
NVSS J005403-374636 and NVSS J005353-374020, can
be seen in the bottom right-hand corner of the image.
Contours are: 0.2, 0.3, 0.4,
0.5, 0.7, 1, 1.5, 2, 3, 4 and 5 mJy
|
| Open with DEXTER | |
Radio-continuum observations of NGC 300 were made on 28 February
2000 with the Australia Telescope Compact Array (ATCA). The array
was in the 6C antenna configuration, with baselines ranging
between 153 and 6000 m. All observations were made simultaneously
at 1374 and 2496 MHz (
and 13 cm). Using this baseline
configuration, we achieved a resolution of 6
for each
frequency (Table 2). This was considered as the best
angular resolution with which to study NGC 300 in detail, ensuring
at the same time sufficient sensitivity to any extended emission
component. Primary flux density calibration was achieved using
PKS 1934-638 and for phase calibration we used the
secondary calibrator J0048-427. Data reductions were
performed using the Multi-channel Image Reconstruction, Image
Analysis and Display ( MIRIAD) software package (Sault & Killeen 2003).
Radio-continuum images of these observations are shown in
Figs. 1 and 2.
Table 2: Summary of radio-continuum observations of the region of NGC 300.
NGC 300 was also observed by Pannuti et al. (2000) (hereafter referred to as PDL) using the Very Large Array (VLA) of the National Radio Astronomy Observatory (NRAO) at a wavelength of 6 cm in the hybrid CnB configuration (northern arm in the B array) on 22 May 1993 and at a wavelength of 20 cm in the BnA configuration on 13 June 1998 (Table 2). To minimize bandwidth smearing, which limits the field of view when observing in normal continuum mode, the VLA observed NGC 300 in multi-channel line mode with seven channels per IF and a channel width of 3.125 MHz, thereby synthesizing a total band of 37.5 MHz after dropping the first of the seven channels. The effective observing frequencies for the images made from the remaining six channels were 4860 MHz (6 cm) and 1448 MHz (20 cm). Data reduction were performed using the Astronomical Image Processing System ( AIPS) software package (Greisen 2003).
The VLA images were corrected for primary beam attenuation and then converted to MIRIAD format before undertaking flux density measurements using the MIRIAD package. A summary of all radio observations is given in Table 2.
XMM-Newton observations on NGC 300 occurred December of 2000
and January of 2001 during revolutions 192 (
37 ks) and 195
(
47 ks), respectively. Data was collected using the EPIC
MOS 1, MOS 2 and PN detectors utilizing a medium filter
(Ehle et al. 2003). XMM-Newton's 5
sensitivity limit is
estimated at
ergs s-1 (1 erg s
-1 = 10-7 Watts).
After processing each dataset (using Science Analysis Software,
SAS; Loiseau 2003), an X-ray image of NGC 300
was created (see Fig. 3). To do this, filtered MOS 1
event files from revolution 192 and 195 were first merged. In
similar fashion, merged MOS 2 and PN event files from both
observations were created. MOS 1 and MOS 2 merged event files were
next combined and finally a single event file containing all
events from both observations created the image. Editing provided
"good time interval'' (gti) files by removing times associated
with contamination from flaring particles.
![]() |
Figure 3:
XMM-Newton (energy range 0.3-6.0 keV) image (contours)
of NGC 300 bounded by our ellipse overlaid on DSS2-Red (R) image (gray scale). Radio
identifications are marked and X-ray contours are:
1.9, 5.8, 17.3, 52.2 and 104.3 in units of
10-2 cnts
|
| Open with DEXTER | |
Source detection ( SAS's edetect chain) was accomplished
using the combined (both observations) filtered event lists from
MOS 1, MOS 2 and PN separately. During source detection, lists
were generated using a minimum likelihood threshold of 15 (which
corresponds to a Gaussian significance of
). These
were then combined to produce a source list. This allowed
determination of the coordinates of preliminary X-ray
sources. We found a total of 533 X-ray sources but further
analysis is in preparation by one of us (S. Carpano). Ten (10)
coinciding sources were chosen visually by comparing overlapping
contours from the XMM-Newton and all four radio images using the
KARMA software package.
For each of the XMM-Newton point-like sources that coincided with radio sources, count rates were determined from each individual detector by first extracting its filtered spectrum and background. Backscale calibration was automatically applied to these files. Response and auxiliary matrix files were created to calibrate the data to the instrumental response of the telescope. It is noted that occasional sources fell very near to the detectors grid; this may have affected a source's hardness ratio somewhat, but not as greatly as the count rate itself.
The total counts and net count rates for each source were found with the High Energy Astrophysics Science Archive Research Center's (HEAsoft's) XSPEC package in the full (0.3-6.0 keV), soft (0.3-1.0 keV), medium (1.0-2.0 keV) and hard (2.0-6.0 keV) energy bands.
The ROSAT archive contains several pointed PSPC and HRI observations centered close to the field of NGC 300. More details on the ROSAT mission can be found in Trümper (1983). The two nearest ROSAT PSPC (energy range 0.1-2.4 keV) observations which have the best angular resolution were used in this study. For more details of the PSPC observations see Table 3.
In Fig. 4 we show the ROSAT PSPC image
("hard'' energy range 0.52-2.01 keV) of the NGC 300
region. It was obtained from combining the two PSPC observations and
binned to 5
pixels.
Several ROSAT HRI observations were performed within the
field of NGC 300 but lower exposure times (total of
40 ks)
did not allow a deeper quantitative study than was obtained from
the PSPC observations (total of
46 ks). For more details on
the PSPC and HRI observations, see Read & Pietsch (2001).
We conducted the multi-frequency analysis of ATCA (1374 and 2496 MHz) and VLA (1448 and 4860 MHz) sources by identifying source positions with data from previous catalogues and our XMM-Newton data. These catalogues include optical observations by BL97, Soffner et al. (1996), Kim et al. (2002), Pietrzynski et al. (2001) and Deharveng et al. (1988); radio analysis by PDL and X-ray data from Read & Pietsch (2001).
Initially, we assume all radio images have a positional accuracy
of less than 1
.
Then, utilizing the viewer in the Karma
software package (Gooch 1996), sources were selected using
multi-frequency contour levels and catalogue annotation files
overlaying a DSS2-Red (R)
image of the galaxy. A elliptical visual boundary to
NGC 300 of
,
centered on RA(J2000) =
Dec(J2000) = -37
,
is provided by the
NASA/IPAC Extragalactic Database (Jet Propulsion Laboratory,
California Institute of Technology and National Aeronautics and
Space Administration). A grid annotation file allowed systematic visualisation
of each source field such that no region would be missed.
Table 3: Summary of X-ray observations of the region of NGC 300.
After converting VLA images to MIRIAD format, task IMSTAT was
used to determine the rms sensitivity of all datasets individually
by averaging several large quiet regions. We found sensitivity
values of 0.058, 0.062, 0.066 and 0.037 mJy
for frequencies of
1374, 2496, 1465 and 4885 MHz, respectively. Assuming a distance
of 2.02 Mpc (Freedman et al. 2001), this corresponds to luminosity
limits of
,
,
and
W Hz-1.
In our visual selection
of radio sources in this galaxy, we were also required to account for poorly cleaned sidelobes
from two nearby radio sources, NVSS J005403-374636 and
NVSS J005353-374020.
The position and flux density for each source at a given radio frequency were determined using the two-dimensional point-fitting algorithms in the MIRIAD software package with primary beam corrected images. This involved creating a region using task CGCURS, containing only the selected source, and analyzing this region with task IMFIT to find the source's flux density (with errors) and position. Since 20% of our sources are extended beyond image resolution, a Gaussian fitting process must include integration over the entire area of the object.
Estimates of the spectral index (
)
of
each radio feature/source were made based on flux densities
obtained from corresponding radio-frequencies. (The spectral index
is defined by the relation
,
where
is the integrated flux density
and
is frequency.) This process was accomplished using a
standard line of best fit algorithm. Errors (
)
were deduced given the scatter in flux density and errors in individual
flux density estimates.
We list in Table 4, for each of 54 detected radio
objects; source number, position (RA and Dec), integrated flux
densities (1374, 1448, 2496 and 4860 MHz) and corresponding errors, spectral index and
error, source type,
comparisons to BL97 images
and identifications with sources at other wavelengths.
![]() |
Figure 4:
ROSAT PSPC (energy range 0.52-2.01 keV, 46 ks)
image (contours) of NGC 300 bounded by our ellipse overlaid on DSS2-Red (R) image
(gray scale). Contours are:
1.4, 1.6, 2.0, 2.7, 3.8, 6.9, 12.8, 24.9, 48.9, 97.0
and 193.2 in units of 10-3 cnts
|
| Open with DEXTER | |
Table 4:
Measured radio-continuum properties of objects within NGC 300.
Listed positions (RA and Dec) are from 1374 MHz observations where
possible.
In Col. 10, visual comparisons are made to BL97 [S II]
versus H
(with some contamination from [N II])
images; 0 = No flux, VF = very faint flux, F = faint flux, M = medium flux and S = strong flux.
Prefixes to source numbers used in Col. 11: H = HRI (Read & Pietsch 2001),
P = PSPC (Read & Pietsch 2001), D = Deharveng et al. (1988), S = Soffner et al. (1996),
DSS2-Red(f/m) = STScI Digitized Sky Survey (1993, 1994),
BL = Blair & Long (1997), RPS97 = Read et al. (1997), PSNR = Pannuti et al. (2000),
GC = Kim et al. (2002), AS_ = Pietrzynski et al. (2001) and XMM = this Paper. Source type abbreviations are:
SNR = supernova remnant, HII = H II region, BKG = background object and xrb = X-ray binary
(capital letters denote higher confidence).
- denotes snr candidates.
Using the methods described above we found ten XMM-Newton
X-ray sources in the area of NGC 300 having radio-continuum
counterparts within an estimated XMM-Newton positional
error of
4
.
We present a contour image of XMM-Newton observations of NGC 300 in Fig. 3 and
Table 6 lists the X-ray sources that match our
radio sources.
Table 5:
XMM-Newton photoelectric absorption
model fits to radio counterparts in NGC 300 (1
confidence interval using MOS 1, MOS 2
and PN data). Source XMM5 was also found to closely fit a power law spectrum. The corresponding
values for this
model are shown in parentheses.
For point like sources with more than 200 counts per observation, we were able to fit simple spectral models to the combined data for each observation using XSPEC after grouping the data to a minimum of 25 counts. Using this technique it was possible to fit a spectrum to two X-ray sources, XMMU J005450.0-373853 and XMMU J005510.7-374835 (our XMM5 and XMM8, respectively). As Table 5 shows, we find the best spectral fit for XMM5 to be the photoelectric absorption with bremsstrahlung model although the photoelectric absorption with power law model was a nearly as good. For XMM8, the best fit was clearly the photoelectric absorption with power law model.
Spectral fitting was more difficult for sources having few counts (<200). For these sources, we used hardness ratios to give some indication of their X-ray spectra. Using our energy bands: soft (0.3-1.0 keV), medium (1.0-2.0 keV) and hard (2.0-6.0 keV); HR1 was defined as the ratio of the difference of medium and soft counts to their sum ((med-soft)/(med+soft)). In similar fashion, HR2 was defined as the ratio of the differences of hard and soft counts to their sum ((hard-soft)/(hard+soft)). Table 6 shows the results of these calculations for data obtained from the PN detector for each observation.
Figure 5 shows a colour-colour diagram of HR1 vs. HR2 for each source and observation. Numerous authors including Haberl et al. (2000) have created similar diagrams to establish a relation between various types of sources including supersoft sources, foreground stars, X-ray binaries and SNRs, utilizing ROSAT observations of the Magellanic Clouds. The definitions of HR1 and HR2 vary with different X-ray telescopes and so it is not possible to use ROSAT source type criteria with XMM-Newton data.
Kong & Di Stefano (2003, and reference therein) have created XMM-Newton colour-colour diagrams in their analysis of IC 342 with energy bands and hardness ratios similar to ours. Superimposed models in their Fig. 4 colour-colour diagram shows that power law spectra tend to occupy the top right section while soft thermal models occupy the lower left. For each of their models, the column density increases from left to right. We used XSPEC's "fakeit'' command to create our own simulated simple model spectra for our colour-colour diagram shown in Fig. 5.
Our hardness ratios primarily add insight to previously selected SNRs as the colours of the X-ray sources are not reliable enough to identify source type alone. It must be remembered that this exercise tells us only about the spectrum of these objects. Recent observations of Galactic SNRs such as RX J1713.7-3946 and SN 1006 (Koyama et al. 1997) have shown that they can have non-thermal as well as thermal properties in X-ray.
Table 7 lists the number of sources in common with our data and other (optical and X-ray) catalogues. This information was compiled from both visual inspection of our data and annotation files (created from the literature) using KARMA'S kview.
Positional identifications with published sources in the optical
and X-ray domains were considered high if they fell within
the positional limits of these surveys (
). An
exception was the ROSAT catalogue created by Read & Pietsch (2001) who
cite positional errors for each of their X-ray sources
individually. Annotation circles were created for this catalogue
with a radius corresponding to these errors in order to find each
corresponding radio source. Next, each individual source
identification was carefully reexamined using contour levels from
different wavelengths before being added to our list.
Table 6: XMM-Newton sources found to match our radio sources showing PN detector hardness ratios and count rates. A complete catalogue of XMM-Newton sources can be found in Carpano (in preparation). Column 4 abbreviation Pe represents the estimated positional error.
Using this process, an interesting systematic error is suspected
with respect to data presented by Soffner et al. (1996) and
Deharveng et al. (1988). When a correction of
Dec = -3.9
was applied to each of these
datasets, the respective positions of their sources matched the
DSS2-Red (R) plate almost exactly. It seems reasonable that errors
of this magnitude occur in the literature presumably due to
registration of their optical image on the sky astrometric frame.
As a matter of background, we now discuss catalogues of NGC 300 found in the literature in more detail.
BL97 detected 28 SNRs candidates through CCD imagery with interference filters
to isolate emission from H
+ [N II], [S II] and a continuum near
610 nm, to identify isolated nebulae bright in [S II] relative to H
.
Optical
long-slit CCD spectra of these candidates confirmed that
these sources met the usual criteria of a flux ratio
([S II]/H
)
of 0.4 or greater.
(All five BL97 SNRs which match our radio sources have
ratios exceeding 0.57.) As noted in Gordon et al. (1999), the physical
basis for this is that in photoionized regions, most sulfur is in
the doubly ionized state. By contrast, in shock-heated nebulae,
the cooling time is short so that many ionization states occur,
including S+. This causes us to observe higher [S
II]/H
ratios in SNRs than in H II regions.
In this study we use BL97 1500s images transferred from tape to
help characterize our sources.
BL97 also determined the optical diameter for their sources,
noting that the majority of the objects were reasonably well
resolved. They based their measurements on the earlier distance
value of 2.1 Mpc, noting that the smallest SNRs in NGC 300 are at
the seeing disk value of 1.5
(1
= 10.2 pc). These
apparent diameters should be considered upper limits. The majority
of BL97 SNRs have apparent diameters in the 25-75 pc range, which
they felt reflected a well-evolved ISM-dominated population for
which their search technique worked best. However, BL97 reported
four (4) large SNRs in NGC 300 with diameters in the 75-200 pc
range for which there are no Galactic counterparts; a trend also
seen in other galaxies (e.g. M 33). These objects may represent
superbubbles; remnants formed from more than one supernova
explosion along with stellar winds found in a group of massive
stars. Dunne et al. (2001) have studied X-rays from superbubbles of
about the same size located within the Large Magellanic Clouds and
found that the X-ray emission is brighter than that
theoretically expected for a wind-blown bubble alone, suggesting that
the X-ray emission has been enhanced by interactions
with several SNRs of similar age.
![]() |
Figure 5: Colour-colour diagram for all correlated XMM-Newton sources using PN hardness ratios. Each source is numbered while filled triangles indicate sources from revolution 192 and open triangles represent those from revolution 195. The estimated hardness ratios from different models are shown as lines connected to "fake'' data points (stars). Power law models tend to have hardness ratios in the upper right of the diagram while thermal models have ratios in the lower left of the diagram. |
| Open with DEXTER | |
H
images were used to determine H II regions
in NGC 300 by Deharveng et al. (1988) and Soffner et al. (1996).
Deharveng et al. (1988) cataloged 176 H II regions within NGC 300 while
Soffner et al. (1996) report finding 90. In addition, Soffner et al. (1996)
show 31 counterpart regions reported in Deharveng et al. (1988).
Related to these H II regions, OB associations composed of groups of young, massive, gravitationally unbound stars formed from molecular clouds within NGC 300 were studied by Pietrzynski et al. (2001). Using the Path Linkage Criterion (PLC), they found 117 such associations. These OB associations were reported to overlap with at least 60% of the H II regions found by Deharveng et al. (1988). Pietrzynski et al. (2001) note that these associations are well suited as tracers of recent or ongoing star formation. In our analysis, we find (Table 4) 22 OB associations that correspond to either H II regions or SNRs.
We are mindful to distinguish these H II regions and OB associations from other objects near NGC 300. Kim et al. (2002) found 17 globular cluster candidates in NGC 300 using U BV I CCD photometry. Their search used photometric information, morphological parameters and visual inspection, and they note that the globular clusters were barely resolved in their images.
Results from ROSAT PSPC and HRI X-ray observations of NGC 300 have been presented by Read & Pietsch (2001). A summary of X-ray sources and their properties from the ROSAT PSPC image (Fig. 4) is presented in their Table 2.
Table 7: Results of multi-frequency source comparison. Numbers represent sources in common to two surveys. XMM-Newton sources are in "()'' if not confirmed. Abbreviations used are: ROSAT = Read & Pietsch (2001); D = Deharveng et al. (1988); S = Soffner et al. (1996), DSS2-Red(f/m) = STScI Digitized Sky Survey (1993, 1994), BL = Blair & Long (1997), XMM-Newton = this paper, K = Kim et al. (2002) and P = Pietrzynski et al. (2001).
Table 8: Selection criteria for radio source type. Note use of logical AND and OR for this table.
Multi-frequency analysis enhances our understanding of the nature of discrete sources in NGC 300. We assume that this galaxy is essentially transparent to radio emission, and, radio and X-ray catalogues contain not only sources located within NGC 300, but also background sources lying behind it. We use a similar classification to that used by Filipovic et al. (1998) to divide the radio sources into two major groups:
Therefore, we use our XMM-Newton data, identifications from six other catalogues, DSS2-Red (R) images, original BL97 images and radio spectral index to create our classification. Because of the distance to NGC 300 (2.02 Mpc), we are unable to use source extension and morphology as part of our criteria.
Specifically, we used the following classification scheme (Table 8) to classify each source within the area of NGC 300 defined by our ellipse:
The accuracy of the positions of all of our radio sources found by methods
described in Sect. 3.1 is
limited by the pointing accuracy and beam size of the ATCA and VLA
telescopes. Since our criteria for radio source
cross-identifications are based initially on a radius of less than
1
and radio source cross-identifications with optical
catalogues on less than 4
,
we expect a small amount of
scatter in position differences if there is true positional
agreement.
We have compared the positions of 16 H II regions cited in
Deharveng et al. (1988) (systematically corrected as discussed Sect. 3.1) to
their ATCA counterparts at 1374 MHz and have found no significant
positional bias between our radio positions and the published
optical positions. At the 95 percent confidence level (using the
student-t with n-1 degrees of freedom), the positional differences
in right ascension and declination (
and
,
respectively) are
,
.
Comparisons of 1374 MHz positions with positions at the other
radio frequencies and XMM-Newton
show the radio catalogues to be self-consistent
(Table 9).
Table 9: Comparisons of positions including 95 percent confidence intervals (CI) and standard deviations. D = Deharveng et al. (1988).
In Fig. 6 we compare ATCA (1374 MHz) to VLA (1448 MHz)
flux densities. We find the line of best fit for this data to be
.
There is no evidence for systematic differences in the flux densities from the ATCA and the VLA.
![]() |
Figure 6:
Comparisons of flux density between ATCA (1374 MHz) and VLA (1448 MHz)
sources with individual flux density error bars plotted on a linear scale.
A weighted line of best fit is
(
|
| Open with DEXTER | |
Because of the relative certainty of the SNRs distance, the
luminosity function is directly proportional to our observed flux densities. The
luminosity of each radio source at a specific frequency,
,
is given in W Hz-1 and
is defined by the relation:
![]() |
(1) |
where D is the distance to NGC 300 of 2.02 Mpc
(Freedman et al. 2001) and
is the flux density at a given radio
frequency. In our case, the source flux density at 1374 MHz was
used to estimate the luminosity of each SNR source. The 5
completeness level at this frequency is 0.29 mJy with a
corresponding level in luminosity of
W Hz-1.
In Fig. 7 we show a histogram of 17 SNRs
(ATCA J005450.2-374030
does not have a flux density at 1374 MHz). The
mean luminosity of the SNRs is
W Hz-1with standard deviation
W Hz-1. Here, we
are seeing only selected high luminosity sources just above the
completeness level. Deeper studies are required for a complete
sampling of SNRs in NGC 300.
Of the five BL97 SNRs identified with our sources, two (BL-S10, Fig. 8d and BL-S11, Fig. 8h) have radio spectral indices consistent with SNRs. We cannot determine a spectral index for the remaining three sources (BL-S6, Fig. 8a, BL-S26, Fig. 8v and BL-S28, Fig. 8u). Additionally, three of these BL97 sources (BL-S6, BL-S10 and BL-26) have corresponding X-ray emission. As noted by Chen & Chu (1998), only a small number of optically identified extragalactic SNRs can be confirmed at radio and X-ray wavelengths with these being among the most luminous remnants.
ATCA J005431.9-373825 corresponds to BL-S6 with a [S
II]/H
ratio of 0.60. This SNR coincides with an OB association and
is seen both in ROSAT and XMM-Newton observations. Our
colour-colour diagram is not helpful in determining the
X-ray spectrum of this source and it contains too few
counts for model fitting.
BL97's S10 identifies with ATCA J005440.6-374049 and is
associated with an H II region and an OB association. Read & Pietsch (2001)
suggest that it is a soft source and our XMM-Newton data implies
that it has a thermal spectrum. Our calculated radio spectral
index of -0.5 is also consistent with an SNR.
![]() |
Figure 7: Luminosity function of SNRs in NGC 300. |
| Open with DEXTER | |
ATCA J005515.4-374439 (BL-S26) also has a thermal spectrum by inspection of its XMM-Newton X-ray colour diagram (Fig. 5, XMM9) and is also labeled P49 in ROSAT observations (see Figs. 8v, 8w and 8x). Its radio spectral index cannot be reliably found since we only have flux densities for two very near frequencies (1374 and 1448 MHz).
In Table 7, we note a total of six (6) XMM-Newton
sources with an existence likelihood of 15 which appear to be valid
sources on visual inspection that match BL97 SNRs within
a 4
positional radius. Thus, possibly as many as three
BL97 optical SNRs have X-ray emission but no radio
detection. We also find four (4) and six (6) X-ray
cross-matches within H II regions from Deharveng et al. (1988) and
Soffner et al. (1996), respectively. Four BL97 SNRs with
X-ray emission are associated with H II regions and only
two of these have radio emission.
![]() |
Figure 10:
Spectrum fits for our XMM-Newton sources 5 and 8. In
each diagram, grouped spectral data show model fits from PN data ( top)
with MOS 1 and MOS 2 data below (x-axis shows energy range, 0.3 keV to |
| Open with DEXTER | |
BL97 notes that while the canonical value of 0.4 in the
ratio of [S II] to H
has worked well in the Milky
Way and Local Group galaxies to separate shock-heated nebulae from
photoionized gas, the exact value of this cutoff has not been
important since there is a gap in the distribution. Photoionized
regions tend to show ratios of
0.2 or below. In Col. 10
of Table 4 we show the results of a visual inspection
of original 1500 s plates from BL97. [S II] plates
were also used to create optical backgrounds for our radio and
X-ray contours shown in Fig. 8. Although not
quantitative, we see [S II] emission from 12 sources other
than the SNRs reported by them. H II regions without an SNR
identification tend to show little or no [S II] emission.
Some of the more interesting objects in Fig. 8 not reported in BL97 include
ATCA J005438.1-374144 (Fig. 8g), in which the typical spherical shape of a SNR
can be seen. Radio emission (with a resolution of about 6
compared to an optical resolution of 1.5
)
with a spectral index of -0.8
is seen from the most
intense area at "4 o'clock'', highlighting that non-thermal radio
emission from SNRs is formed at the shock. In fact, the
expanding shell of a SNR may be undergoing different processes at
the same time causing some areas to emit non-thermal radio
radiation while other regions do not (also see the same process in
BL-S28 shown in Fig. 8u). This may also be the case
with sources ATCA J005450.2-374030 and ATCA J005450.7-374022. BL97
noted that there appeared to be 4 large SNRs in their data that
measured up to 200 pc (20
). With these two sources, we may
be detecting radio emission from one large SNR or superbubble (see
Figs. 8k and 8l).
Radio emission from what appears to be a shell shock front is also seen in J005451.1-373826 (Fig. 8n) and J005500.5-374037 (Fig. 8p), with the latter having a more complete shell appearance. ATCA J005443.1-374311 (Fig. 8i), J005445.3-373847 (Fig. 8j), J005450.3-373822 (Fig. 8m), J005451.7-373939 (Fig. 8o), J005503.5-374246 (Fig. 8r) and J005512.7-374140 (Fig. 8t) have less defined shells. ATCA J005503.6-374320 (Fig. 8s) may show only a half shell in [S II], but this would put the diameter of the remnant at about 100 pc and its shape is not well defined.
We cannot determine [S II]: H
ratios of the above sources from
existing BL97 observations.
We feel that optical long-slit CCD spectra from these SNRs initially found in our radio data and now suspected from
visual inspection of BL97's [S II]
images, would be helpful in their final determination as SNRs. This information might
also fine tune the [S II]: H
ratio cutoff between H II regions and SNRs.
We found thirteen additional SNRs than reported by BL97 based on radio spectral index and identification with H II regions or OB associations (Table 4). Although most of them have some optical emission as discussed above, there is a general lack of X-ray emission from this group of objects. This emphasizes that radio wavelengths are especially important for the detection these objects.
SNR ATCA J005451.7-373939 has a positional identification with both a globular cluster (GC6) and an OB association. It seems likely in this case, that the radio emission is associated with the OB association and the globular cluster is simply superimposed, since cluster SNRs most likely dissipated into the interstellar medium long ago.
We find three background sources based on positional identifications and spectral index. Although most of these sources represent distant galaxies, one of these, ATCA J005453.3-374311, is coincident with a globular cluster as reported by Kim et al. (2002). Its X-ray colour (XMM7 in Fig. 5) and radio spectral index is consistent with a power law spectrum.
ATCA J005425.2-374441 also has a very steep spectral index of -1.8, but its X-ray colour (XMM1) is a bit softer. Still, it is certainly possible that this object also has a X-ray power law spectrum as evidenced by its X-ray colour proximity to ATCA J005510.8-374835 (Fig. 9d) which clearly corresponds to an AGN (XMM8) as reported by Read & Pietsch (2001). Using XMM-Newton data, we were able to fit a photoelectric absorption power law model to the latter as shown in Figs. 10c and 10d (also see Table 5).
Based on spectral index alone, there are at least two other probable background sources. It is not possible to tell if the remaining 10 "bkg'' objects are really background objects as they could represent a different source type.
ATCA J005450.3-373850 (Fig. 9a) is possibly
associated with an X-ray binary
(XMM5) as reported by Read & Pietsch (2001) which they base on its
distance of 30
from the nearest SNR and its HR2 value of
0.75. Using HRI data, they also found a 41% probability that the
X-ray source is variable. We have labeled this source
"xrb'' since the spectrum we found is also consistent with their
conclusion. Our X-ray model fit shows the object to have a
photoelectric absorption Bremsstrahlung spectrum (Figs. 10a and 10b).
However, the photoelectric absorption power law
fit is also very good (see XMM5 in Table 5). Perhaps
this object has both thermal and power law properties. Its radio
emission may be the result of synchrotron emission from jets,
making this a rare example of a microquasar. Only one other
suspected extragalactic microquasar has been reported by
Dubus (2003) (in the nucleus of M 33). We expect this type of
source to be rarely found in extragalactic studies since detection
of the relatively faint radio emission is unlikely at far
distances.
![]() |
Figure 11: Venn diagram showing the intersection of selected sets of SNRs for NGC 300. Optical refers to the objects listed in BL97. Parenthesis show the total number of SNRs in each category if our additional [S II] optical identifications are added to those reported by BL97. The "?'' denotes that at this time, we do not know the total number of X-ray SNRs within this galaxy. That information is being prepared by S. Carpano for later independent publication. |
| Open with DEXTER | |
PDL could place only limits on the radio spectral
indices from eight of their 14 radio sources because they had only
one flux density (
= 20 cm) for each of them. One of
the goals of the present study is to expand on those observations
using new ATCA and XMM-Newton data. Our findings do support the
six SNRs identifications that they made for radio sources having
VLA flux densities at both 6 and 20 cm wavelengths and we pick up a 6 cm
flux density for their R1 allowing us to agree with a total of seven of their sources.
We also name their R2 as a possible "snr'' only, since its association with Read & Pietsch (2001) P44
has been called
into doubt by those authors (see page 980 of Pannuti et al. 2002).
While we disagree with six of PDL's SNRs, with the addition of ATCA flux densities in
the 13 and 20 cm wavelength range, we find an additional six
radio SNRs. Updating Fig. 9 in Pannuti et al. (2002) we present a Venn
diagram depicting the intersection of sets of SNRs within NGC 300
in Fig. 11. The diagrams are very similar except that we
do not include "possible'' X-ray sources and although the
number of radio sources is almost identical, they do not
represent the same sources between the two studies.
The optical and radio numbers in parentheses represent the number of
SNRs in these categories if our additional BL97 identifications with radio sources
are included.
In summary, out of a total of 54 radio sources, we have found and described 18 SNRs in NGC 300; 13 in addition to those found in the optical by BL97. Three of the five SNRs identified by BL97 that have radio counterparts also have X-ray counterparts. We find evidence for [S II] emission from an additional 12 of our radio sources, implying that a total of 17 SNRs have some optical component.
There are three additional SNR candidates in NGC 300 which will await further studies for verification. These sources do not appear to have any optical component, although ATCA J005528.2-374903 is not within the boundary of the BL97 images.
All identifications with ROSAT sources are included in our XMM-Newton X-ray identifications with radio sources except source number 2 (coincident with PDL's R2 discussed above) reported in Read et al. (1997). Although this source could be a supernova remnant also, we do not have a radio spectral index to prove it.
We have investigated the luminosity functions of our SNRs and find that conclusions about them as a group cannot be made since we are selecting for only the most luminous SNRs.
While it has been understood that many SNRs lie embedded within H II regions, we also note a strong correlation of our radio sources with OB associations within this galaxy. This is another good place to look for SNRs from massive progenitor stars.
We also find the expected number of background sources within the region covered by NGC 300. We note that other objects, for example, globular clusters near NGC 300, may lie in the same direction, thus adding some confusion with identification of sources. An additional seven XMM-Newton X-ray sources match our radio sources other than the SNRs mentioned above. Three are identified with background objects and one is associated with a likely X-ray binary as discussed above. Also, two X-ray sources could be background objects while one identifies with another X-ray supernova candidate (P36) discussed in Read & Pietsch (2001). The latter's counterpart, XMM10, shows hardness ratios consistent with a thermal spectrum but we have no additional information to prove it as a SNR.
Acknowledgements
We thank Wolfgang Pietsch and Andy Reid for allowing us to use their ROSAT image of NGC 300 for our analysis. T.G.P. would like to thank the staff of the ATNF for their hospitality during the observing run. We also used the Karma software package developed by the ATNF and the EXSAS/MIDAS software package developed by the MPE. E. M. Berkhuijsen deserves thanks for her advice and kind patience.
We also thank William Blair who graciously provided us with the images of NGC 300 used in BL97. These Los Campanas Observatory observations were transferred from tape at our request.
DSS2-Red (R) plates were based on photographic data obtained using the UK Schmidt Telescope. The UK Schmidt Telescope was operated by the Royal Observatory Edinburgh, with funding from the UK Science and Engineering Research Council, until 1988 June, and thereafter by the Anglo-Australian Observatory. Original plate material is copyright (c) of the Royal Observatory Edinburgh and the Anglo-Australian Observatory. The plates were processed into the present compressed digital form with their permission. The Digitized Sky Survey was produced at the Space Telescope Science Institute (STScI) under US Government grant NAG W-2166. TGP gratefully acknowledges support from a Chandra General Observer Award issued by the Smithsonian Astrophysical Observatory.
![]() |
Figure 8:
SNRs with radio and X-ray contours
in NGC 300 overlaying BL97 1500s [S II]
images.
Radio contours are |
| Open with DEXTER | |
![]() |
Figure 8: continued. |
| Open with DEXTER | |
![]() |
Figure 8: continued. |
| Open with DEXTER | |
![]() |
Figure 8: continued. |
| Open with DEXTER | |
![]() |
Figure 9:
Other interesting objects in NGC 300 with radio and X-ray contours
overlaying the DSS2-Red (R) image. 1374 MHz contours are 0.174 and 0.291 mJy
|
| Open with DEXTER | |