Issue |
A&A
Volume 503, Number 3, September I 2009
|
|
---|---|---|
Page(s) | 855 - 858 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200911807 | |
Published online | 15 July 2009 |
The
analysis of recently detected radio planetary nebulae in the Magellanic Clouds
(Research Note)
B. Vukotic1 - D. Urosevic2,4 - M. D. Filipovic3 - J. L. Payne3
1 - Astronomical Observatory Belgrade,
Volgina 7, 11160 Belgrade-74, Serbia
2 -
University of Belgrade, Faculty of Mathematics, Studentski trg 16, 11000 Belgrade, Serbia
3 -
University of Western Sydney, Locked Bag 1797, Penrith South, DC, NSW 1797, Australia
4 - Isaac Newton Institute of Chile, Yugoslavia Branch, Yugoslavia
Received 6 February 2009 / Accepted 26 June 2009
Abstract
Aims. We investigate and analyze the radio surface brightness to diameter ()
relation for recently detected, bright radio-continuum planetary nebulae (PNe) in the Magellanic Clouds (MC).
Methods. We applied a Monte Carlo analysis in order to account for sensitivity selection effects on measured
relation slopes for bright radio PNe in the MCs.
Results. In the
plane, these radio MCs PNe are positioned among the brightest of the nearby Galactic PNe and are close to the D-2 sensitivity line of the MCs radio maps. The fitted Large Magellanic Cloud (LMC) data slope appears to be influenced by survey sensitivity. This suggests the MCs radio PN sample represents just the ``tip of the iceberg'' of the actual luminosity function. Specifically, our results imply that sensitivity selection tends to flatten the slope of the
relation. Although MCs PNe appear to share the evolution properties of Galactic PNe, the small number of data points prevented us from further constraining their evolution properties.
Key words: methods: statistical - ISM: planetary nebulae: general - galaxies: Magellanic Clouds - radio continuum: ISM
1 Introduction
Filipovic et al. (2009) has recently reported 15 radio-continuum PNe detected at the
level and above from various radio Magellanic Clouds (MCs) surveys. Their detections are mainly based on positional coincidences with optically detected planetary nebulae (PNe) and optical spectroscopy (Payne et al. 2008a,b). Included are Large Magellanic Cloud (LMC) mosaics (Dickel et al. 2005) with sensitivities of
0.5 mJy beam-1 at both 4.8 and 8.64 GHz with resolutions of 33 and 20
,
respectively. The Small Magellanic Cloud (SMC) mosaics have 0.5 mJy beam-1 sensitivities at 4.8 and 8.64 GHz with resolutions of 30 and 15
(Dickel 2009). To compare these sources with Galactic PNe, we use radio surface brightness (
)
since this quantity is distance independent:
![]() |
(1) |
where S is the flux density and



All of the detected radio PNe are unresolved, despite five of them being observed in ``snap-shot'' mode with the resolution of up to 1
and sensitivity of
(for details see Filipovic et al. 2009). The optical diameters are available for some 10 out of 15 (see Table 1 in Filipovic et al. 2009). The MCs PNe data at 4.8 GHz are plotted in Fig. 1 with respective LMC map
sensitivity line at
(hereafter referred to as sensitivity line). Five LMC PNe are above the sensitivity line and one is below (``snap-shot'' mode). Five PNe with estimated upper flux density limit of
are plotted on the sensitivity line and the remaining one without optical diameter is not plotted. Four SMC PNe are plotted on the SMC map resolution limit of 30
.
We use the PNe sample of Galactic PNe with reliable individual distances from Stanghellini et al. 2008; hereafter SSV, and divide them into two subsamples (nearby PNe with distances smaller than 1 kpc and those with distances greater than 1 kpc). Urosevic et al. (2009), in an extensive analysis of the empirical
relation for PNe, argued that the SSV sample of nearby PNe is least influenced by selection effects. For comparison, we plot two radio PNe from the Sagittarius dwarf galaxy (Dudziak et al. 2000), also representative of extragalactic radio PNe. An arrow in Fig. 1 represents the D-2 direction in which unresolved objects should move if their diameters were known.
![]() |
Figure 1:
A plot of 4.8 GHz PNe from the LMC, SMC and Sagittarius Dwarf Galaxy.
The solid line represents the LMC map sensitivity ( |
Open with DEXTER |
Inspection of Fig. 1 reveal that the detected radio-continuum MCs PNe are very close to their sensitivity lines and they are as bright as the brightest Galactic PNe. When compared with Galactic PNe it appears that current data set is only the very peak of the actual MCs PNe luminosity function. In this paper, using
analysis and Monte Carlo generated data samples, we investigate the sensitivity related selection effects on the
properties of these objects and compare them with the properties of Galactic PNe.
As a convenient platform for distance determination, the
relation
has a
50 year long history of exploitation for supernova remnants (SNRs). The study of the relation for these very luminous radio synchrotron sources began with the work of Shklovskii (1960), while the first empirical
relation was presented by Poveda & Woltjer (1968). Decades of analysis resulted in a good understanding of the
relation for SNRs and various data sample selection effects that influence that relation.
The previous analysis of PNe radio data samples (e.g. Phillips 2002, and references therein) were mainly focused on radius vs radio brightness temperature (
)
empirical relations. These were sparse in the consideration of data selection effects. Recent works by Urosevic et al. (2007) and Urosevic et al. (2009) use a SNR formalism in the analysis of PNe data along with a discussion of these selection effects. Although they note that all PN samples suffer from selection effects caused by limitations in survey sensitivity and resolution, the nearby SSV sample, having a
slope of
,
is least influenced by these effects.
Table 1: Fit parameters of the LMC data sample.
2 Analysis
Samples for Monte Carlo simulations are made using sources from Table 1 in Filipovic et al. (2009) that have a known optical diameter (LMC sources only) and are above the map sensitivity line for a given frequency. This includes 5 PNe at both frequencies (8.64 and 4.8 GHz) though not necessarily the same ones. Using adopted distances to the LMC of 50 kpc and to the SMC of 60 kpc (Alves 2004; Hilditch et al. 2005), the smallest sources in the selected samples have a mean geometrical diameter of 0.09 pc at both frequencies. The largest source has a mean geometrical diameter of 0.12 pc at 4.8 GHz and 0.15 pc at 8.64 GHz. Surface brightness ranges from
to
Wm-2 Hz-1 sr-1 at 8.64 GHz and from
to
Wm-2 Hz-1 sr-1 at 4.8 GHz. For the LMC maps, the 4.8 GHz sensitivity line is at
Wm-2 Hz-1 sr-1 with a D-2 break at D=8.00 pc and a 8.64 GHz sensitivity line at
and D=4.85 pc, respectively.
The relative flux density errors at both frequencies of the LMC maps are <10% (Filipovic et al. 2009) and they are approximately the size of the symbols shown in Fig. 1. If we assume that the optical diameters have no significant errors and that the absolute flux density error equals the flux density standard deviation
,
according to error propagation theory it follows from Eq. (1). that the standard deviation for
(
)
is
.
This gives
.
The resulting parameters of the
fits are given in Table 1 for each frequency. The parameters a and
,
their standard uncertainties
and
,
the linear correlation coefficient r, the fit quality r2, the probability Q of obtaining larger WSSR (weighted sum of square residuals) and a ratio of WSSR to the number of degrees of freedom (ndof) are given in this table. A value of
is expected for a statistically acceptable fit (the case when errors of the dependent variable have a non-Gaussian distribution), if the data is were well approximated by the model. This should apply for fits presented in this paper since we have transformed flux density errors to
errors. When introduced in the least squares fitting procedure, our values of Q implies a statistically acceptable fit.
Although the fits are statistically believable the resulting
at 4.8 GHz substantially differs from
at 8.64 GHz, likely in debt to a small number of data points. Sampling only the peak of the luminosity function and with the majority of the sample likely hidden below the sensitivity line, the resulting values for
should not be taken seriously. Comparing the LMC radio PN sample and the sample of a nearby Galactic PNe, we made an attempt to access a more meaningful
slope using Monte Carlo simulations as described below.
We performed Monte Carlo simulations similar to those described by Urosevic et al. (2005). The results of their simulations (with 21 data points above the sensitivity line) showed that slopes steeper than -2 are under influence of sensitivity related selection effects. First, we determined the empirical
standard deviation from the best fit line, assuming
as the independent variable. We then selected an interval in
five times as long as that of the real data. This interval is then sprinkled with random points of the same
density as that of the real data.
The simulated points, that lie on the log D axis, are then projected onto a series of lines at different slopes (in steps of 0.1 from -3.5 to -1.5). Each of these lines passes through the extreme upper left hand end of the best fit line to the real data. We also added Gaussian noise
in
,
which is related to the scatter of the real data by a parameter called ``scatter''. A scatter of 1 corresponds to the same standard deviation as that of the real data.
An appropriate sensitivity cutoff is applied to the simulated data points, selecting points above the sensitivity line. This is done 100 times for each simulated slope and a least squares best fit line is generated for artificial samples that have five or more selected data points. The number of such samples is given in the last column of Tables 2 and 3.
In Tables 2 and 3, the first column lists the value of the simulated slope, while the mean and standard deviation of the best fit slopes for the generated samples are given in second and third column, respectively. The fourth and fifth column gives the mean and standard deviation of the best fit slopes for sensitivity selected generated samples, respectively. In Fig. 2 we present one of our Monte Carlo generated samples at 8.64 GHz for a scatter of 1 and the simulated slope -2.7.
Table 2: Monte Carlo simulation results at 4.8 GHz.
Table 3: Monte Carlo simulation results at 8.64 GHz.
![]() |
Figure 2: The Monte Carlo generated sample at 8.64 GHz for a scatter of 2 and simulated slope of -2.7. The LMC data points (asterisks) are plotted along with the sensitivity (thick solid) line; artificially generated points are plotted above (filled circles) and below (open circles) this line. The selection of points above the sensitivity line give a fit with a flatter slope (thin solid line) than the slope of a fit using all simulated points (dotted line). |
Open with DEXTER |
3 Discussion and conclusions
From Tables 2 and 3, it is evident that the sensitivity cutoff tends to flatten the intrinsic (simulated) slopes less than -2 (the most likely case for the real data). At 4.8 GHz, a slope of
(shallower than the sensitivity line) is not influenced by sensitivity related selection effects. From the other side, this slope does not have any physical interpretation and is only the consequence of scatter in the
plane. The 8.64 GHz slope of
appears to correspond to somewhat steeper selection free slopes but due to a large error could have any value smaller than -2.3. This is within the range of the SSV slope at 4.8 GHz. With current LMC data samples it is not possible to constrain the lower limit, due to a smaller number of selected points for steeper slopes. The value of a mean slope after selection at 8.64 GHz starts to oscillate for the slopes
-2.4 (the standard deviation of slope after selection exceeds the interval between two successive simulated slopes). Table 2 is given here rather for completeness, but nevertheless shows that for the samples with small number of data points and a slope greater than -2 it is not possible to extract any meaningful information with this kind of simulations.
With these recently confirmed MCs radio PNe being amongst brightest Galactic PNe in the
plane and given their proximity to the sensitivity line, it is likely they represent only the peak of the actual MCs PNe luminosity function, which is probably the ``bright-end'' extension of the Galactic PNe luminosity function (for a more detailed discussion of the nature of these MCs PNe see Filipovic et al. 2009). The results of our LMC Monte Carlo simulations suggest that current sparse data samples cannot give meaningful
slope. However, they confirmed that the corresponding apparent
slope flattens toward
because of the sensitivity related selection effects..
This should be kept in mind as a serious measure of caution in forthcoming studies of extragalactic PN samples. Future high sensitivity images of the MCs will surely provide even better and more complete radio samples of PN population and enable more robust constraining of PN evolution parameters.
Acknowledgements
The authors would like to thank the referee Prof. John Dickel for useful comments that have improved the manuscript. This research has been supported by the Ministry of Science and Technological Development of the Republic of Serbia through the Project #146012.
References
- Alves, D. R. 2004, New Astron. Rev. 48, 659
- Dickel, J. R., McIntyre, V. J., Gruendl, R. A., & Milne, D. K. 2005, AJ, 129, 790 [NASA ADS] [CrossRef] (In the text)
- Dickel, J. R., et al. 2009, in The Magellanic System: Stars, Gas, and Galaxies, ed. J. van Loon, & J. M. Oliveira (Cambridge University Press), in press (In the text)
- Dudziak, G., Péquignot, D., Zijlstra, A. A., & Walsh, J. R. 2000, A&A, 363, 717 [NASA ADS] (In the text)
- Filipovic, M. D., Cohen, M., Reid, W. A., et al. 2009, MNRAS, in press, [arXiv:0906.4588v1] (In the text)
- Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 [NASA ADS]
- Payne, J. L., Filipovic, M. D., Crawford, E. J., et al. 2008a, Serb. Astron. J., 176, 65 [NASA ADS]
- Payne, J. L., Filipovic, M. D., Millar, W. C., et al. 2008b, Serb. Astron. J., 177, 53 [NASA ADS]
- Phillips, J. P. 2002, ApJS, 139, 199 [NASA ADS] [CrossRef] (In the text)
- Poveda, A., & Woltjer, L. 1968, AJ, 73, 65 [NASA ADS] [CrossRef] (In the text)
- Shklovskii, I. S. 1960, AZh, 37, 256 [NASA ADS] (In the text)
- Stanghellini, L., Shaw, R. A., & Villaver, E. 2008, ApJ, 689, 194 [NASA ADS] [CrossRef] (In the text)
- Urosevic, D., Pannuti, T. G., Duric, N., & Theodorou, A. 2005, A&A, 435, 437 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Urosevic, D., Vukotic, B., Arbutina, B., & Ilic, D. 2007, Serb. Astron. J., 174, 73 [NASA ADS]
- Urosevic, D., Vukotic, B., Arbutina, B., Ilic, D., et al. 2009, A&A, 495, 537 [NASA ADS] [CrossRef] [EDP Sciences]
Footnotes
- ... relation
- In this paper, we use the form
, where
represents slope.
All Tables
Table 1: Fit parameters of the LMC data sample.
Table 2: Monte Carlo simulation results at 4.8 GHz.
Table 3: Monte Carlo simulation results at 8.64 GHz.
All Figures
![]() |
Figure 1:
A plot of 4.8 GHz PNe from the LMC, SMC and Sagittarius Dwarf Galaxy.
The solid line represents the LMC map sensitivity ( |
Open with DEXTER | |
In the text |
![]() |
Figure 2: The Monte Carlo generated sample at 8.64 GHz for a scatter of 2 and simulated slope of -2.7. The LMC data points (asterisks) are plotted along with the sensitivity (thick solid) line; artificially generated points are plotted above (filled circles) and below (open circles) this line. The selection of points above the sensitivity line give a fit with a flatter slope (thin solid line) than the slope of a fit using all simulated points (dotted line). |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.