Free Access
Issue
A&A
Volume 552, April 2013
Article Number A12
Number of page(s) 24
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201220345
Published online 13 March 2013

© ESO, 2013

1. Introduction

The simultaneous study of the chemical composition of planetary nebulae (PNe) and H ii regions in galaxies brings several interesting benefits. While H ii regions have since long been used to probe the present-day chemical composition of the interstellar medium (Peimbert 1975), the use of PNe, in principle, offers the possibility to study its fossil record, since PNe are descendants of stars that were formed up to a few Gyrs ago (Peimbert 1990; Maciel & Köppen 1994). Another point is that, in the case of a high-metallicity medium, abundances derived in H ii regions can be biased due to significant temperature gradients in their interiors, while those derived in PNe are not (Stasińska 2005). On the other hand, the chemical composition of PNe is not necessarily identical to that of the material out of which the progenitors were formed. It has long been known that, for example, the nitrogen observed in PNe has been partly produced at various stages of the progenitor evolution (Renzini & Voli 1981). Other elements may be affected as well, in smaller proportion but in a way that is far less well understood. Therefore, the comparison of abundances in PNe and H ii regions can also provide important clues on the nucleosynthesis in low- and intermediate-mass stars of different metallicities.

In the case of spiral galaxies, H ii regions have revealed the presence of radial abundance gradients (Searle 1971; Shields 1974; McCall et al. 1985), which give clues on the chemical evolution of those galaxies. One of the important questions to sort out is how these gradients evolve with time. Do they steepen? Do they flatten? Previous studies based on PNe in the Milky Way have led to contradictory results (Maciel & Köppen 1994; Maciel et al. 2003; Stanghellini & Haywood 2010). One difficulty in the study of abundance gradients from PNe located in the Milky Way is that their distances are in general very poorly known. This is why turning to other galaxies may be rewarding, despite the increased distance.

The face-on spiral galaxy M 33, at a distance of 840 kpc, has been thoroughly studied from the point of view of its ISM chemical abundance, both in terms of its H ii regions (Vilchez et al. 1988; Rosolowsky & Simon 2008; Magrini et al. 2007, 2010; Bresolin 2011) and, more recently, its PNe (Magrini et al. 2009 for the entire range of galatic radii; and Bresolin et al. 2010 for the inner zones). According to Magrini et al. (2010), PNe and H ii regions appear to show a similar metallicity gradient.

NGC 300 is the next nearest, nearly face-on spiral galaxy, at a distance of 1.88 Mpc (Gieren et al. 2005). It is an Scd galaxy, quite similar to M 33 in many respects, located in the Sculptor group. Its main properties are listed in Table 1. Its giant H ii regions have been recently studied by Bresolin et al. (2009), using the FORS 2 spectrograph at the Very Large Telescope. These authors were able to determine for the first time in this galaxy a radial gas-phase oxygen abundance gradient using only temperature-based abundance determinations,  − 0.077 dex/kpc, with a central abundance 12 + log O/H  ≃  8.57. This agrees very well with the trend in metallicity obtained for A and B supergiant stars by Urbaneja et al. (2005) and Kudritzki et al. (2008). In this paper, we use the same instrument to investigate the chemical composition of the brightest PNe and compact H ii regions in this galaxy. In a separate paper, the analysis of the pre-imaging (on-band and off-band [O iii] λ5007 observations) obtained with the same telescope is presented (Peña et al. 2012), and a list of candidate PNe is given.

This paper is organized as follows. Section 2 presents the observational material and the data reduction. Section 3 presents the determination of the physical conditions, and the ionic and total abundances of our targets. Before getting to the heart of this paper, Sect. 4 revisits the classification of the candidate PNe and H ii regions. Section 5 presents classical abundance ratio diagrams for our NGC 300 objects. Section 6 presents the radial behavior of the abundance ratios derived from PNe and H ii regions. Section 7 provides our interpretation of the results, and Sect. 8 summarizes the main findings of our work.

Table 1

Main properties of NGC 300.

2. Observations and data reduction

2.1. Direct imaging

In 2006 we started an observational program aimed at finding a large number of PNe in NGC 300 and analyzing their chemical properties. The Very Large Telescope (VLT) of the European Southern Observatory (ESO) at Cerro Paranal, equipped with the FORS 2 spectrograph, was used for this purpose. [O iii] λ5007 on- and off-band pre-imaging data, was obtained using the FILT-500-5 and OIII/6000 filters, respectively, to search for emission line objects. We covered two regions with the 6.8 × 6.8 arcmin2 field of view of FORS 2: one, centered at RA = 00h54m49.0 s, Dec =  −37°41′02.0′′, includes the center of the galaxy, and the other, centered at RA = 00h55m22.0s, Dec =  −37°43′00.0′′, samples the galaxy outskirts (see Fig. 1 of Peña et al. 2012).

After subtracting the on-band and off-band [O iii] λ5007 images, we searched for emission line objects in the two fields. Hundreds of such objects were found. We considered as PN candidates those objects that appear stellar with no stellar residual in the off-band image. With such a criterion, we found more than 100 PN candidates: about 68 in the central field and about 34 in the external one. We rediscovered most of the 34 PN candidates reported by Soffner et al. (1996). The analysis of the results obtained from our pre-imaging data is presented in Peña et al. (2012). We selected about 50 objects (both PN candidates and other emission line objects) for follow-up spectroscopy with FORS 2 in multiobject (MXU) mode.

Table 2

Summary of spectroscopic observationsa.

2.2. Spectroscopy

Multiobject spectra were obtained on August 19–22, 2006. Masks with 1 arcsec-wide slit were designed for the FORS2 MXU mode to observe as many targets as possible without affecting the desirable wavelength range, in the two regions of the galaxy covered by our narrow-band images. Of course the brightest PN candidates were chosen first, but a certain number of faint objects were included too. The range of log F(5007) values for the observed PNe covers from  − 14.693 to  − 16.485 which correspond to [O iii] λ5007 magnitudes m(5007) from –23 to –27.5. Three grisms were used in order to cover the widest wavelength range possible: 600B (approximate range 3600–5100 Å, 4.5 Å FWHM spectral resolution), 600RI (approximate range 5000–7500 Å, 5 Å  spectral resolution), and 300I (approximate range 6500–9500 Å, 10 Å spectral resolution). Multiple exposures were carried out for each setup. Our exposures are listed in Table 2. Total exposure times for the central field were 3.5 h for grism 600B, 3.3 h for grism 600RI and 2 h for grism 300I. For the external field the exposure times were 3.33 h for grism 600B, 1.5 h for grism 600RI and 0.5 h for grism 300I. Three spectroscopic standard stars (EG 274, LDS 749B and BMP 16274) were observed each night, through a 5′′-wide slit, for flux calibration. Seeing conditions were in general very good (0.7′′–0.9′′), except for the first night (August 19–20), when the seeing was larger than 1′′. Frames containing the same spectral data were combined and then reduced with standard procedures by using the EsoRex pipeline provided by ESO: bias subtraction, flat fielding, and wavelength calibration were applied at this stage. Afterwards we used standard IRAF1 tasks for spectral extraction and flux calibration.

Table 3

Line intensities for H ii regions and PNe, relative to Hβ = 1.00123.

In Fig. 1 the calibrated spectra (showing the full wavelength range) for all the observed objects are shown.

2.2.1. Measuring and dereddening the fluxes

Line intensities were measured in the flux-calibrated spectra with the task splot in IRAF by integrating between two positions over a continuum level estimated by eye (actually there is no detectable continuum in most of the objects, in particular PNe). Then we proceeded to calculate the logarithmic reddening correction, c(Hβ) from the Balmer decrement. Due to the fact that we have three independent spectral ranges (blue – 600B, visible – 600RI and red – 300I) that can be affected by differences in flux calibration, we determined the reddening from the blue spectra, where Hβ, Hγ, Hδ and the upper lines of the Balmer series are present. The logarithmic reddening correction, c(Hβ), was computed from these lines by assuming case B recombination theory (Storey & Hummer 1995). The lines measured in the blue spectral range were dereddened with this c(Hβ) value and by adopting the Seaton (1979) reddening law. In the 600RI and 300I spectra only Hα is available among the H Balmer series lines, so we adopted the c(Hβ) value obtained in the blue to deredden the lines in these spectral ranges, relative to Hα (again Seaton’s reddening law was used). Afterwards we matched the blue, visible and red portions of the spectra by assuming a dereddened Hα/Hβ flux ratio as given by recombination theory [I(Hα)/I(Hβ)  ~  2.86]. The line flux uncertainties (as well as upper limits) were determined as 2 ×  the rms noise measured on each side of the line.

In Table 3 we present the dereddened line fluxes and associated uncertainties for all the lines detected in each object. For lines that are important abundance diagnostics upper limits are given. At the end of the table we also list the c(Hβ) values and the logarithm of the total observed Hβ flux (since the full table is available at the CDS, here we show just an example, for illustration purposes). The deprojected galactocentric distance of each object, given in terms of the 25th magnitude B-band isophotal radius R25 = 9.75′, and calculated by adopting the parameters given in Table 1, is also included in Table 3.

It is interesting to notice that c(Hβ) is very low in all our targets, in particular in PNe where almost no reddening was found. Only a few compact H ii regions show significant reddening. This is consistent with the low foreground E(B − V) values determined by other authors, ranging between 0.013 (Schlegel et al. 1998) and 0.096 (Gieren et al. 2005).

3. Determination of physical conditions, ionic and total abundances

3.1. Atomic data

Since the atomic data adopted in the computation of abundances translate immediately into the results, we present them first. A number of new atomic calculations has recently become available, attempting to improve on previous results. One emblematic case is that of O+, which is used as a density diagnostic, and for which the collision strenghts computed by McLaughlin & Bell (1998) gave results inconsistent with those of other ions, as shown by Copetti & Writzl (2002) and Wang et al. (2004). This prompted new atomic computations and a thorough discussion to address this problem (Kisielius et al. 2009). The latter authors also recommend the use of the transition probabilities from Zeippen (1982).

New collision strenghts and transition probabilities have been computed for S+ by Tayal & Zatsarinny (2010). We used the [S ii] λ6731/6717 and [O ii] λ3726/3729 line ratios observed by Wang et al. (2004) in Galactic PNe to check the densities obtained using the Tayal & Zatsarinny (2010) atomic data. In Fig. 2 we show the data with their error bars in the [S ii] λ6731/6717 vs. [O ii] λ3726/3729 plane, together with the line of equal electron density at a temperature of 104 K for densities ranging from ne = 20 cm-3 (upper right) to ne = 105 cm-3 (lower left), shown in red. It is clear that, at the high-density end – the most drastic regime in our study, as will be seen later – the observational points lie well below this line, meaning that the [S ii] λ6731/6717 ratio indicates a subsantially higher density than the [O ii] λ3726/3729 one (by factors as large as 10), and many of the observed values of [S ii] λ6731/6717 are outside the high density limit. Even if one does not expect the mean densities in the O+ and in the S+ region to be identical, due to the fact that a density gradient is likely to exist and that the O+ and S+ zones are not exactly coextensive, such a large difference looks unrealistic. We have therefore taken the [S ii] transition probabilities from Mendoza & Zeippen (1982), which lead to the line of equal density shown in blue. This line goes beautifully through the observational points, and therefore we adopted this reference for the transition probabilities of [S ii].

New atomic data are also available for N+ and Ar++, which we use in the present work. McLaughlin et al. (2011) announced new collision strenghts for Ne++ but the data are not available yet. We list in Table 4 the references for the atomic data used in the present work for ions of O, N, Ne, S and Ar. Note that for O+ we preferred the transition probabilities from Zeippen (1982) over more recent estimates, as recommended by Kisielius et al. (2009).

thumbnail Fig. 2

Values of the density indicator [S ii] λ6731/6717 as a function of [O ii] λ3726/3729. Points with error bars: observed values in the sample of PNe by Wang et al. (2004). Continuous lines: lines of equal electron density at a temperature of 104 K from ne = 20 cm-3 (upper right) to ne = 105 cm-3 (lower left). Red: using the [S ii] transition probabilities by Tayal & Zatsarinny (2010). Blue: using the [S ii] transition probabilities from Mendoza & Zeippen (1982).

thumbnail Fig. 3

Logarithms of the emissivities of the [O ii] λ3727, [N ii] λ6584, [S ii] λ6716, λ6731, and [Ar iii] λ7135 lines for densities ne = 102, 103, 104 and 105 cm-3. Red: with the atomic data used in the present work. Blue: with the atomic data used in our previous work.

The emissivities of the different lines were obtained with a 5-level atom in each case. In Fig. 3 we compare the present emissivities of the [O ii] λ3727, [N ii] λ6584, [S ii] λ6716, λ6731, and [Ar iii] λ7135 lines in nebular conditions (red curves), with respect to those obtained with atomic data used in our previous work (Peña et al. 2007; Bresolin et al. 2009, 2010) (blue curves), for densities ne = 102,103,104 and 105 cm-3 (the lower curves corresponding to lower emissivities). Overall, the differences are not very large, the largest ones reaching 0.2 dex for [Ar iii]. Incidentally, this plot also shows how much the [O ii] λ3727 emissivity depends on density as soon as it exceeds 103 cm-3, meaning that it is essential to have a good determination of the density for a proper determination of O/H and N/O.

For H+ and He++ we use the case B emissivities from Hummer & Storey (1987), and for He+ we use the emissivities from Porter et al. (2007), including collisional excitation. In the case of optically thin PNe, case A would be more appropriate, implying i) a possible error in the helium abundance; and ii) a possible underestimate of heavy element abundances with respect to hydrogen by up to 30%. However we do not expect such objects to be present in our sample, which is composed of luminous PNe, which are likely to be rather optically thick (see Stasińska et al. 1998). Many of our objects show a detectable [O i] line giving observational support to such a claim (although the possibility of density-boundedness in some directions cannot be discarded).

3.2. Plasma diagnostics and ionic abundances

In Table 5 we list the values of the electron temperatures derived from the [O iii] λ4363/5007 and [N ii] λ5755/6584 line ratios obtained from Table 3, as well as the densities derived from [S ii] λ6731/6717 and [Ar iv] λ4740/4711. We give in brackets the 15 and 85 percentiles, obtained by a Monte Carlo procedure as explained in Sect. 3.4. This means that there is a 70% probability that the true value lies between those percentiles (to be compared to the 68.3% probability to be within one standard deviation for a Gaussian distribution). Similarly to Table 3, Table 5 is ordered by object number, irrespective of the nature of the object (PN or compact H ii region). The first column indicates the object type as assigned in Sect. 4: PN for planetary nebulae and H ii for compact H ii regions.

The ionic abundances were derived from the following lines (when observed): He iλ5876, He iiλ4686, [O ii] λ3727, [O iii] λ5007, [N ii] λ6584, [Ne iii] λ3869, [S ii] λ6716, λ6731, [S iii] λ9069, and [Ar iii] λ7135 with respect to Hβ. In the few cases where the [O ii] λ3727 line was not observed but the [O ii] λ7320,7330 was, we used the latter to derive the O+ abundance. When only an upper limit was available for a line intensity, the intensity was put to zero for the calculations. We checked that the errors on the abundances resulting form this procedure are smaller than the uncertainties presented in Sect. 3.4.

For PNe, we assumed that the temperature of all the ions is equal to T[O iii], the temperature given by [O iii] λ4363/50072. This is, of course, not entirely correct, as temperature gradients are expected to occur inside the nebulae. One could be tempted to use the empirical relation found to exist between T[O iii]/T[N ii] and He ii 4686/Hβ by Kingsburgh & Barlow (1994) or Wang & Liu (2007). However, this relation shows an important spread and its trend is opposite to that seen in the sample of Peña et al. (2001) and García-Rojas et al. (2012). In the few cases where the temperature from [N ii] λ5755/6584 was available, we decided not to use it for abundance determination, the error bars being too large to present any advantage over the temperature from [O iii] λ4363/5007. For the compact H ii regions in our sample, we adopted the same procedure to relate the temperatures of the different ions as described in Bresolin et al. (2009), since there are many theoretical and observational arguments to justify it, as can be seen in Fig. 2 of Bresolin et al. (2009)3.

In many cases, the only information available on the electron density comes from the [S ii] λ6731/6717 line ratio. We adopt this value as representative for the entire nebula. A good evaluation of the density is crucial for O+ and S+ but not so much for the other ions, since the critical density for collisional deexcitation of [O iii] λ5007, [N ii] λ6584, [Ne iii] λ3869, [S iii] λ9069, and [Ar iii] λ7135 is rather high. There are a few cases where the density could be obtained only from [Ar iv] λ4740/4711, and this was adopted for all the ions. In the cases were the densities could be obtained from both [S ii] λ6731/6717 and [Ar iv] λ4740/4711  we chose the density from [S ii] λ6731/6717 to represent the nebula, since, as just mentioned, it is the density in the low ionization zone which is more critical for abundance determinations. However, there are some objects for which the error bar on the density from [Ar iv] λ4740/4711 is significantly smaller than that from [S ii] λ6731/6717 or the density derived from [Ar iv] λ4740/4711 is more reasonable. For such cases (i.e. for objects 25, 35, 91, 103) we used the density obtained from [Ar iv] λ4740/4711. When the spectra did not give any information at all on density (6 cases), we adopted a value of 3000 cm-3 for PNe, and 100 cm-3 for H ii regions. We note that, for objects 91 and 103, the only ones where the density could be derived from both [S ii] λ6731/6717 and [Ar iv] λ4740/4711 with a certain degree of reliability, the density from [Ar iv] is smaller than that of [S ii] by a factor 5 and 3, respectively. It is not clear whether this reflects a true difference in density between the high and low excitation zones. Given that, in the one hundred Galactic PNe observed by Wang et al. (2004), the differences between the two diagnostics are not that large, we are rather inclined to consider that our error bars might have been underevaluated. Note that in both cases, the [S ii] lines are extremely weak.

Table 6 lists the ionic abundances determined for all our objects.

3.3. Elemental abundances

To obtain the elemental abundances from ionic abundances, one must correct for unseen ions. The present situation with ionization correction factors (ICFs) is far from satisfactory. In view of the lack of better options, for PNe we use the ICFs from Kingsburgh & Barlow (1994), which were obtained from detailed models for 10 PNe. However, these models were never published, and it is not clear whether the modeled objects are similar to the ones studied in the extragalactic case, where it is the brightest PNe that are chosen for spectroscopic studies. In addition, our spectra encompass the entire spatial extent of the PNe, while those of Kingsburgh & Barlow (1994) correspond to only small zones of the objects. For H ii regions, we follow the same policy as Bresolin et al. (2009). It uses the ICFs from Izotov et al. (2006), which are based on sequences of photoionization models aimed at reproducing the observed properties of giant H ii regions. The only exception is represented by neon, for which we adopt Ne/O    =    Ne++/O++, since the ICF proposed by Izotov et al. (2006) leads to a trend in Ne/O as a function of excitation. Lastly, in the case of helium we do not attempt to use any ICF, since there is no reliable way to correct for neutral helium in our objects.

Table 7 lists the nominal abundances of O, N, Ne, S, Ar with respect to H, as well as the sum of the abundances of He+ and He++ with respect to H. We also indicate the 15 and 85 percentiles, obtained from our Monte Carlo analysis of errors, as explained in the next section.

Since we wish to compare the abundances of PNe with those of H ii regions, we recomputed the abundances of the giant H ii regions from Bresolin et al. (2009) with the same atomic parameters and procedures as for the compact H ii regions presented in this paper. Our abundances are in general very similar to the ones presented in Bresolin et al. (2009), but differences can amount to 0.1–0.2 dex in the objects with the highest abundances, due to the different atomic data employed.

To test the adopted ICFs, we plot in Fig. 4 the values of Ne/O, S/O and Ar/O as a function of O++/(O++O++). Throughout the paper, PNe are represented by circles, compact H ii regions by stars, and the giant H ii regions from Bresolin et al. (2009) by green squares. No significant trend with excitation is seen for Ar/O and S/O, although for the latter we find a tendency for PNe to give a lower value than H ii regions. We are here witnessing the sulfur problem discussed at length by Henry et al. (2012). For Ne/O, there is no trend with excitation in the case of PNe, but a slight one for H ii regions – both giant and compact ones – (which had already been noted by Bresolin et al. 2009 for giant H ii regions and probably implies that the ICFs for Ne in H ii regions should be revised).

In the top left panel of Fig. 4 we plot the values of (He++He++)/H as a function of O++/(O++O++). Some of these values are very low, clearly indicating the presence of neutral helium. There is no reliable way to correct for that, since the ionization of helium is mainly governed by the temperature of the ionizing star, while that of the heavy elements also depends on the ionization parameter.

thumbnail Fig. 4

Various abundance ratios as a function of O++/(O+ + O++) for PNe and H ii regions in NGC 300. PNe are represented by circles, compact H ii regions by stars and giant H ii regions by squares. The big red star is object no. 39, the big black circle is object no. 74 (see Sect. 4), and the big black star is object no. 57 (see Sect. 3.3). The objects represented by open circles are objects 12, 88 and 109 (see Sect. 3.4).

3.4. Uncertainties

thumbnail Fig. 5

Histogram of log  O/H + 12 values for object PN 48 from the Monte Carlo simulation. The nominal value is indicated by the continuous line, the 15 and 85 percentiles by the dashed lines.

To estimate the uncertainties in the derived abundances we start from the uncertainties in the line intensities listed in Table 3. Those were estimated as explained in Sect. 2.2, and include the uncertainties involved in the dereddening process. We use a Monte Carlo procedure to follow the error propagation in the determination of the plasma parameters and the abundances. We assume for each line flux a Gaussian distribution centered on the flux effectively measured and having a dispersion equal to the estimated flux uncertainty, and given in Table 34. For each line, we considered 200 independent realizations, and estimated the error bars by taking the 15 and 85 percentiles of the final results5. The reason for considering percentiles and not one sigma deviations is that the final distributions on abundances can be highly non-Gaussian, as shown in Fig. 5 for the case of PN 48, and we are in fact interested in the range of values in which there is a 70% chance for the real value to be found (nb: in the case of a Gaussian distribution, an interval of  ± 1σ corresponds to a probability of 68.3%). Note that the nominal values of the abundances listed in Table 7 do not necessarily correspond to the median or the most probable values of the Monte Carlo simulation. Our estimation for the final uncertainties on the abundances does not include uncertainties on ICFs, for which we have no handle. This is unfortunately the situation for all abundance determinations using ICFs that have been published so far.

For PN 97 and 109 no lines from O+ are measured, therefore the abundance of that ion cannot be estimated. The total oxygen abundance should however not be drastically affected since in that case the contribution of O+ to the O/H abundance ratio is expected to be small. In PN 12, 88 and 109 the He iλ5876 line is not observed and O/H is obtained simply by adding O++ and O+ (when measured), so that the computed O/H is a lower limit to the true oxygen abundance, since it does not account for the presence of higher ionization species of oxygen. In all the plots in this paper these three objects are represented by open symbols, and are not considered in the computation of the radial abundance gradients in Sect. 6. On the other hand, PN 97 is represented with the same symbols as the rest of PNe and is included in the computation of the gradients in Sect. 6.

Concerning H ii regions, Stasińska (2005) made the point that abundances derived using temperature-based methods in metal-rich nebulae are prone to strong errors, due to important temperature gradients with the nebulae. From the behavior of the radial variation of O/H in the giant H ii regions observed by Bresolin et al. (2009), it does not seem that these objects belong to the abundance regime where such effects appear (except perhaps the two innermost H ii regions): the radial variation of O/H is rather smooth, and the dispersion can be explained by observational errors only.

4. How to distinguish PNe from H ii regions in NGC 300

thumbnail Fig. 6

Diagrams illustrating observed properties that allow to differentiate between PNe and H ii regions in NGC 300 (see Sect. 4). The symbols have the same meaning as in Fig. 4. The few PNe for which the electron density cannot be measured have been placed at log Mion = 3.

In Peña et al. (2012), the candidate PNe were defined as point-like objects seen in [O iii] emission but not in the continuum. If the objects were seen in the continuum, they were considered as candidate H ii regions. The spectroscopic data presented in this paper allow us to revise our classification on a more solid basis.

In Fig. 6 we show a few diagrams that are useful for this purpose. In this figure, as well as in all the similar figures in this paper, object number 74 is represented by a big black circle, and object number 39 is represented by a big red star. From our photometric analysis, initially object 74 was considered as a compact H ii region and object 39 as a PN but, as we will see, spectroscopic data argue for the opposite. In all the diagrams of Fig. 6, we can see a clear segregation between candidate H ii regions and candidates PNe. The Hβ luminosities of the latter tend to be lower that those of candidate H ii regions. As expected, many candidate PNe show He iiλ4686 emission while none of the H ii regions do. Many candidate PNe have higher N/O ratios than H ii regions, similarly to what is found for well-known PNe in the Milky Way (Henry et al. 2000). PNe tend to have a higher [O iii] temperature at a given O/H since, on average, the effective temperatures of their exciting stars are much higher that those of massive stars that ionize H ii regions. For the same reason, PNe tend to have a higher excitation, as measured by the O++/(O++ O++) ratio. They also tend to have higher densities. Indeed, PNe with densities similar to those of H ii regions would be too faint to be detected at the distance of NGC 300, since the Hβ luminosity is proportional to the gas density times the nebular mass which, in the case of PNe is at most of the order of a few solar masses. For a temperature of 104 K, relevant for our sample, the mass of ionized gas is given by Mion=37.5L(Hβ)/n\begin{equation} \label{eq:mass} M_{\rm ion} = 37.5 \, L(\Hb) / n \end{equation}(1)where L(Hβ) and Mion are in solar units, and n is the mass-weighted average electron density in cm-3. In practise, we use the [S ii] or [Ar iv] ratios to estimated this density.

The L(Hβ) versus Mion relation (Fig. 6 bottom left) clearly shows the sequence: PNe – compact H ii regions – giant H ii regions. Of course, the estimates of Mion are coarse, the densities being uncertain, as seen in Table 5, but even considering those uncertainties, the trend is clearly present. In all the diagrams of Fig. 6 objects number 39 and number 74 appear at the junction between PNe and compact H ii regions. So why do we consider object 74 as a PN and object number 39 as an H ii region? There are strong arguments for that. If object 74 were a compact H ii region, from its Hβ luminosity one would deduce that the effective temperature of the ionizing star would be about 30 000 K, using the relation between the Lyman continuum output and the effective temperature of main sequence stars from Martins et al. (2005). But the observed value of O++/(O++ O++) for this object is 0.79, clearly incompatible with such a low effective temperature. On the other hand, for object 39, using the equation above, we find an ionized mass of around 100 M, clearly incompatible with a PN.

Such sanity checks on the classification of point-like extragalactic nebulae are very important. In the present case, only 2 objects out of 35 would have been misclassified without the use of spectroscopic data. But the proportion could be much larger in other studies, especially if the preselection is done with a telescope too small to detect the continuum emission of an O-type star.

5. Abundance patterns of PNe and H ii regions in NGC 300

thumbnail Fig. 7

Abundance patterns of PNe and H ii regions in NGC 300. The symbols have the same meaning as in Fig. 4.

If Fig. 7 we present some classical diagrams comparing the abundance patterns of PNe and H ii regions in NGC 300.

The most evident feature is the systematically larger N/O ratio in PNe relative to H ii regions. Such a behavior is observed in other galaxies like, for example, the Magellanic Clouds (Leisy & Dennefeld 2006), NGC 6822 (Hernández-Martínez et al. 2009) or M 33 (Bresolin et al. 2010). As seen in panel e of Fig. 7, N/O appears to be correlated with the helium abundance as estimated from (He+ + He++)/H (discarding the PNe PN 12 and 88 for which no He i line is measured). This is again similar to what was found in previous studies of other galaxies and is in agreement with the idea that the observed enhancement of nitrogen is due to the second dredge-up, which brings CNO-cycle processed material to the stellar surface.

The Ne/O ratio is very constant among both PNe and H ii regions. There is a small offset (by about 0.2 dex) seen in Panels b and f between the two types of objects. Such an offset has already been shown to occur in other galaxies (e.g., Bresolin et al. 2010). It is difficult to say whether it tells us something about the physics of these objects rather than simply being due to inadequate ICFs (see however Sect. 7). Note that the compact H ii region H ii 57 has surprizingly high value of Ne/O (0.29) for a reason that we do not understand.

The Ar/O values also do not show any trend with O/H, both for H ii regions and PNe. For the latter objects, the dispersion is larger than for Ne/O, probably both because the [Ar iii] λ7135 lines are weaker than the [Ne iii] λ3869 lines and because the ICF used for argon is extremely coarse.

For sulfur, we find S/O ratios systematically lower for PNe than for H ii regions, similar to what was found by other authors in many samples of PNe (see Henry et al. 2012 for a recent discussion).

6. Radial abundance gradients in NGC 300

thumbnail Fig. 8

Abundances with respect to hydrogen as a function of the normalized galactocentric radius R/R25. The symbols have the same meaning as in Fig. 4. The continuous blue lines and the dashed red lines represent the straight line fits to the observations of H ii regions and of PNe, respectively.

thumbnail Fig. 9

Abundances with respect to oxygen as a function of the normalized galactocentric radius R/R25. The layout of the figure is the same as in Fig. 4.

We now focus on the radial behavior of the properties of PNe and H ii regions (considering compact H ii regions and giant H ii regions together).

Figure 8 shows the radial distribution of abundance ratios with respect to hydrogen, as a function of R/R25, where R25 is the radius out to the blue surface brightness level of 25 mag arcsec-1, corrected for inclination and galactic extinction. This normalization of the radius is the most widely used since Vila-Costas & Edmunds (1992), and allows one to soundly compare the radial behavior of abundance ratios in different galaxies. The symbols have the same meanings as in all the figures starting from Fig. 4. The estimated error in the line ratios is indicated for each object. One can see that, in the central zones, error bars for giant H ii regions are larger than those for PNe, despite the fact that the former are more luminous by about two orders of magnitude. This is due to the fact that metal-rich H ii regions have low electron temperatures, and an error in their estimated values strongly affects the resulting abundance. Note that the error bars reported in Fig. 8 take into account only error propagation from the observed line fluxes, but consider that the scheme adopted to evaluate the temperatures of the different zones is perfectly valid, so the real uncertainties in the abundances with respect to hydrogen are likely higher, especially in H ii regions.

Figure 9 is analogous to Fig. 8 but shows abundance ratios with respect to oxygen.

The radial trend of the abundances in compact H ii regions is similar to that in the giant H ii regions, although one may note a tendency for compact H ii regions to have systematically lower oxygen and neon abundances than giant H ii regions at similar galactocentric radii. We do not see any convincing explanation for this fact. Larger statistics would be needed to first ascertain whether the effect is real or not. On the other hand, PNe and H ii regions appear to have convincingly different radial behaviors. This is especially clear for abundance ratios involving O, N and Ne (and also S, but for this element, as noted in Sect. 5 the abundances in PNe are uncertain).

Table 8 presents the results of straight-line fitting of various abundance ratios as a function of fractional galactocentric distance, R/R25. The fitting was done taking into account the uncertainties in abundance ratios, but discarding errors in radial distances. For each abundance ratio, we performed one fitting for PNe, and one for H ii regions (merging compact and giant H ii regions). For each fitting, the table gives the values of A (the intercept) and B (the slope) in the equation y=A+B×R/R25,\begin{equation} \label{ eq:gradient} y = A + B \times R/R_{25}, \end{equation}(2)as well as the number of objects considered.

Table 8 clearly shows that the radial abundance gradients of PNe are significantly shallower than for H ii regions in the case of O/H, Ne/H, Ar/H and S/H. As a matter of fact, the abundance gradients for PNe are almost flat. We should note, however, that the abundance dispersion at any galactocentric radius is larger for PNe than for H ii regions. The abundance gradient that we find for O/H in our sample of H ii regions is slightly flatter than found by Bresolin et al. (2009):  − 0.36 ± 0.05 as compared to  − 0.41 ± 0.03. As can be deduced from Fig. 8, this is due to the inclusion of the compact H ii regions.

We also note that the values of A are smaller for PNe than for H ii regions, in a very significant way in the case of O/H.

Concerning the radial behavior of the abundance ratios with respect to oxygen we observe that Ne/O shows no gradient at all, with the Ne/O value being systematically higher by almost 0.2 dex for PNe than for H ii regions at any galactocentric distance. There is no gradient in Ar/O for H ii regions and a mild one for PNe (whose significance is however not obvious, given that the ICFs in the determination of the argon abundance are very crude, as mentioned earlier). Concerning S/O, the regression lines for PNe and H ii regions run almost parallel, with the PNe values being systematically lower than the H ii region ones by 0.2 dex. From those three ratios, given the gentle trends and small differences between PNe and H ii regions, and considering the uncertainties in the abundance derivation, it would be premature to draw any conclusion of astrophysical value.

For N/O, on the contrary, the difference in the observed trends is spectacular. While for H ii regions the gradient in N/O is similar to the one in O/H (as already found in Bresolin et al. 2009), there is almost no gradient in the N/O values determined in PNe – just a large dispersion. Besides, at any galactocentric distance, almost all the observed values of N/O in our PNe are strongly above the H ii region trend! This implies that in most of the PNe in our sample the N/O ratio does not reflect that of the progenitor’s star at birth, but has been heavily modified by dredge-up of nucleosynthesis products.

7. Discussion

We now try to understand the main features of the behavior of abundances in PNe and H ii regions in NGC 300 within a more general framework. The most outstanding findings requiring an explanation are:

  • 1.

    The N/O ratios in most PNe are larger by0.4–1 dex than inH ii regions at similargalactocentric distances.

  • 2.

    The O/H ratio for PNe presents a larger intrinsic scatter than for H ii regions6.

  • 3.

    The O/H radial gradient obtained from PNe is flatter than that obtained from H ii regions.

In order to aid our discussion, we list in Table 9 a few relevant properties of PN progenitors, i.e. of stars with initial masses Mi between 1 M and 7 M. These are: the main sequence lifetimes, tMS, the “final” masses, i.e. the masses of the central stars of the PNe they produce, and some characteristics of the computed chemical yields: the values of Δlog  N/O and Δlog  O/H, which are the predicted changes of the PN chemical composition with respect to the initial composition of the progenitor atmosphere. Those values have been obtained from recent stellar evolutionary models (Marigo 2001; Karakas 2010; Lagrange et al. 2012). They are listed for three different metallicities: solar, Z = 0.008 and Z = 0.004. In Fig. 10 we show the evolution of the PNe of different central star masses, using a simple toy model where a spherical nebula of homogeneous density and constant mass (0.5 M) expands with a velocity of 20 km s-1 around the central star. For the stellar evolution, we use the set of hydrogen-burning post-AGB evolutionary tracks of Blöcker (1995, and references therein), with the same interpolation procedure as in Stasińska et al. (1997). The stars are assumed to radiate like black bodies, and L(Hβ) has been computed analytically assuming an electron temperature of 104 K. While obviously the evolution of real PNe is far more complex than that (for example the covering factor is not necessarily unity, the density is not homogeneous and the expansion velocity is not uniform), this simple toy model serves to fix order of magnitude properties. In the left panel, the black rectangle delimits the region covered by the PNe we observed in NGC 300. The right panel shows the time evolution of L(Hβ) for the same models, starting from the ejection of the nebula.

thumbnail Fig. 10

Evolution of the PN toy model described in Sect. 7 for PNe having progenitor masses (from right to left) of 2 M (red), 2.5 M (green), 3 M (blue), 4 M (cyan), 5 M (magenta). For 1 M, the Hβ luminosity would be too small to be detected. The black box in the left panel indicates the observed range for our PNe. The right panel indicates the time evolution of L(Hβ) since the formation of the PN model. It illustrates that the chance to observe a PN having a progenitor mass of 2 M is much larger than for 2.5 M.

thumbnail Fig. 11

Observed and predicted radial trends of O/H and N/O in the PNe of NGC 300. The thin lines show the predicted trends assuming that the progenitor initial abundances follow the regression line for H ii regions (represented with the thick blue line). The different colors correspond to the different stellar evolution models listed in Table 9: green for Marigo (2001), grey for Karakas (2010), blue for the standard models of Lagarde et al. (2012), purple for the models with rotation of Lagarde et al. (2012). The thick dashed lines shows is regression line for the observed PN abundances.

By considering Fig. 10 and Table 9, one can see that the bulk of our PNe likely arise from stars of initial masses around 2–2.5 M, with masses closer to 2 M being probably more frequent since their PNe are sufficiently bright for a longer time. This implies that the progenitors were born in a rather narrow time interval: roughly 0.5–1 Gyr ago, according to Table 9.

It follows that the oxygen abundance scatter observed among PNe at any galactocentric distance is likely not due to sampling of different epochs of the galaxy chemical evolution. Another way to produce the observed abundance scatter would be the presence of radial stellar motions. These would imply that a PN seen now at a certain galactocentric distance may have been formed in a radially distinct zone. The importance of radial stellar motions in NGC 300 has been assessed by Gogarten et al. (2010) by means of color–magnitude diagrams of stellar populations in different radial bins. With the help of Monte Carlo simulations, these authors found that stars in their model for NGC 300 do not undergo important radial migration. This leaves us with only one alternative to explain the O/H dispersion: the modification of the original abundances by nucleosynthesis and mixing during the progenitor evolution. Indeed, the abundance of oxygen can be increased by dredge up of He-burning products or decreased by hot bottom-burning. For example, the models of Marigo (2001) listed in Table 9 show the effect of third dredge-up on oxygen in a stellar mass range appropriate to NGC 300. However, the models by Karakas (2010) as well as those by Lagarde et al. (2012), which include more complete physics than Marigo (2001), do not predict an important oxygen abundance variation. This is illustrated in Fig. 11, which shows the radial variation of the oxygen abundance in the PNe of NGC 300, assuming that the initial abundance is equal to that of the H ii regions (represented by the thick blue line) and taking the predictions from the different models shown in Table 9 for an inital mass of 2.5 M. Clearly, only the models of Marigo (2001) produce a noticeable trend. However, this trend, which would reverse the gradient with respect to that of the H ii regions, is not in agreement with the observations. We should note, however, that the available grids of models were done using specific parametrization of processes such as mass-loss or convection and do not explore the entire parameter space. One could think that, in different conditions, models could predict a different behavior of the oxygen abundance. For example, Decressin et al. (2009) find, for their solar metallicity models with initial rotation velocity of 300 km s-1 at the surface (i.e. larger than the Lagarde et al. models) that the oxygen abundance has decreased by 0.65 dex at the end of the early AGB phase.

We note that also in M 33 the scatter of the O/H values observed for PNe is much larger than that observed for H ii regions at similar galactocentric distances (Bresolin et al. 2010). Many O/H values observed in PNe in that galaxy are actually significantly higher than any of the O/H values observed in H ii regions. This suggests that oxygen production in the PN progenitors in that galaxy is quite common, too.

Finally, we have found that in NGC 300 the O/H radial gradient derived from PNe is flatter than that obtained from H ii regions. Taken at face value, and ignoring the considerations above, this might be considered as evidence that the metallicity gradient in this galaxy is steepening with time, since PNe are an older population compared to H ii regions. However, we have just argued that the oxygen abundances in PNe are affected by nucleosynthesis in the progenitors and this can happen in both ways. Thus it is difficult to determine the metallicity gradient at birth from the observed O/H ratios. We note, however, that argon – which is not suspected of being modified by progenitor nucleosynthesis – also shows a flatter galactic gradient in PNe than in H ii regions (see Fig. 8 and Table 8), although not so clearly as in the case of O/H, and the dispersion is even larger. However, as mentioned several times in this paper, the ICFs for argon are very crude, so the argon abundances are not fully reliable. On the other hand, there is a priori no reason why uncertainties in the argon abundances would conspire to produce an artificial flattening of the gradient. The available data on PNe would therefore indicate that there is a steepening of the metallicity gradient in NGC 300 over the last Gyr, however the evidence is not really strong. One worry is the fact that at R/R25 = 0.8, the argon abundances seen in PNe are larger than those in the interstellar medium as traced by the H ii regions by about 0.1–0.2 dex. It is not clear whether this is significant. Sulfur is another element that is not suspected of being modified by progenitor nucleosynthesis. However, as mentioned in Sect. 3.4, there is a well documented problem with sulfur abundances in PNe, so until this problem is not solved, it is not recommended to use PN sulfur abundances to test chemical evolution scenarios. Since we have argued that the oxygen abundances of the PNe observed in NGC 300 have been affected by nucleosynthesis in their progenitors, there is one more topic to be discussed. This is the small dispersion in the Ne/O ratio in PNe and the fact that this ratio is higher than the one in H ii regions by about 0.2 dex. It is known that, while 20Ne, the main isotope of neon, does not change its abundance during the evolution of intermediate mass stars, 22Ne can be produced through combined H and He burning during the thermally pulsing evolution (Karakas & Lattanzio 2003). The range of stellar masses in which this production can significantly affect the total Ne abundance in the models of Karakas (2010) is Mi = 2.5–3.5 M at a metallicity Z = 0.008 and Mi > 2 M at a metallicity Z = 0.004. So the PNe of our sample could well be affected by this process, which would explain the fact that their Ne/O ratio is larger than in H ii regions. The question is: why are the O and Ne abundances so strongly correlated, since they are not produced by the same process, as discussed in Karakas & Lattanzio (2003)? (Note that we have argued in Sect. 5 that the ICFs for neon might be inadequate, but from Fig. 4 the Ne/O ratios calculated for our PNe do not depend on excitation.)

As regards the behavior of the N/O ratio in the PNe, Table 9 shows that for such progenitor masses as relevant for our NGC 300 sample, PNe are expected to show a N/O enhancement of about 0.3 to 0.6 dex (depending on model and metalicity). Our observations, however, show that, in the central parts of the Galaxy, N/O ratios in PNe range from –1.2 (i.e. close to the value found in H ii regions) to 0.35. So the observed dispersion is much larger than predicted. In the outer part of the galaxy, they are all above –0.4. The thin lines representing the radial behavior for N/O for a progenitor mass of 2.5 M do not follow the observed behavior very well, for any of the considered models. As in the case of the O/H ratio, the observed behavior of the N/O ratio in PNe argues for a large variety of conditions that affect the mixing processes. Note that the behavior of N/O in NGC 300 is actually similar that in M 33 (Bresolin et al. 2010). This means that the requirements for stellar evolution models are similar for both galaxies, which is reassuring.

8. Summary

We have obtained high signal-to-noise spectroscopy with the VLT of 26 PNe and 9 compact H ii regions in the nearby spiral galaxy NGC 300. These objects were detected in the [O iii] on- and off-band imaging survey of Peña et al. (2012) and supplement the data obtained by Bresolin et al. (2009) for a sample of giant H ii regions. After a discussion of the available atomic data, we have determined the abundances of He, N, O, Ne, S and Ar in all these objects in a consistent way.

We have discussed the nature of the candidate PNe and compact H ii regions in the light of our new spectroscopic information. It turns out that one object that was identified as a candidate PN from the photometric data is in fact a compact H ii region, while one candidate compact H ii region is in fact a PN.

Concerning the general abundance patterns shown in Fig. 7, we find, not surprizingly, that compact H ii regions have chemical abundances that are globally similar to those of giant H ii regions. PNe have systematically larger N/O ratios compared to H ii regions, while they have similar Ne/O and Ar/O ratios. The sulfur abundances in PNe show an unphysical trend with excitation, probably indicating an inadequate ionization correction factor for this element. All theses results are in rough agreement with recent findings in other galaxies such as the Small and Large Magellanic Clouds (Leisy & Dennefeld 2006), NGC 3109 (Peña et al. 2007), NGC 6822 (Hernández-Martínez et al. 2009), and M 33 (Bresolin et al. 2011). The large N/O ratios in PNe as compared to H ii regions is an indication that N has been altered by nucleosynthesis in the progenitors. While similar results obtained in other galaxies are generally interpreted as being the result of second dredge-up, this cannot be the case here, since we have shown that most of the PNe in our sample likely arise from stars of initial masses around 2–2.5 M. An extra mixing process is required to bring the freshly produced N to the stellar surface before the ejection of the PN, likely driven by rotation.

With regard to the radial behavior of the chemical abundances across the galaxy, we find that compact H ii regions follow the giant H ii regions rather well (with however a tendency for compact H ii to have slightly lower abundances with respect to hydrogen than giant H ii rgions). On the other hand, PNe behave very differently. In the central parts of the galaxy PNe have an average O/H abundance ratio 0.15 dex smaller than H ii regions. The abundance dispersion at any galactocentric radius is significantly larger for PNe than for H ii regions. This suggests that the oxygen abundance in PNe is also affected by nucleosynthesis in their progenitors, in a way which depends on several parameters. The formal O/H, Ne/H and Ar/He abundance gradients for PNe are shallower than for H ii regions. In fact, they are almost flat. We have argued that this indicates a steepening of the metallicity gradient in NGC 300 over the last Gyr, rather than the effect of radial stellar motions. However the large observed abundance dispersion makes any conclusion on a possible steepening of the gradients with time only tentative.

We end this study by recalling several problems that are general to the study of the chemical abundances of PNe, which we have not considered here and should be kept in mind. One is the unsatisfactory situation with the ionization correction factors used to compute the abundances: this is known to be important for the sulfur abundance, but all the other elements, including oxygen, might be affected to a certain degree. An additional problem is the amount of depletion into dust grains, which has been brought up by Rodríguez & Delgado-Inglada (2011) for PNe and H ii regions in the solar vicinity. Finally, there is the still unsolved problem of the discrepancy between abundances derived from collisionally excited and from recombination lines (see e.g. the book by Stasińska et al. 2012 for a summary), which casts some doubt on the relevance of the computed abundances. The reliability of the conclusions to which we arrived at in this study – as in any other study of this kind – depends somehow on the solutions that will be found to those problems.

Online material

thumbnail Fig. 1

Calibrated spectra of our observed PNe and compact H ii regions.

thumbnail Fig. 1

continued.

thumbnail Fig. 1

continued.

thumbnail Fig. 1

continued.

thumbnail Fig. 1

continued. Note that for PN 74, the red spectrum has been obtained with grism 300I, which has a much lower resolution than 600B and 600RI.

thumbnail Fig. 1

continued.

thumbnail Fig. 1

continued.

thumbnail Fig. 1

continued.

Table 4

Atomic data used in our abundance determinations.

Table 5

Electron densities and temperatures for the PNe and compact H ii regions in NGC 300.

Table 6

Ionic abundances for the PNe and compact H ii regions in NGC 300.

Table 7

Computed abundances for the PNe and H ii regions in NGC 300a.

Table 8

Straight-line fitting of the radial variations of abundance ratios for the PNe and H ii regions in NGC 300.

Table 9

Typical properties of PNe and their progenitors.


1

IRAF is distributed by the NOAO, which is operated by the AURA, Inc., under cooperative agreement with the NSF.

2

For one of our objects, PN 20, no direct measurement of the electron temperature is possible, so no abundances have been calculated.

3

This way of proceeding assumes that we already know which of our objects are PNe and which are compact H ii regions, and we had to iterate after making these decisions as explained in Sect. 4.

4

This is perhaps not the most realistic distribution, but it is convenient. Wesson et al. (2012) have used more sophisticated distributions, but the strong biases they claim for line intensities with signal-to-noise ratios  ≤ 4 are not supported by observations of line ratios such as [O iii] λ4959/5007 or [N ii] λ6548/6484 which, on average, are consistent with the values predicted by atomic physics (see, e.g., Fig. 10 of Bresolin et al. 2005).

5

Whenever this procedure lead to a negative value for a line intensity, which may happen principally for low signal-to-noise cases, it was replaced by one tenth of the measured value. Since we use percentiles and not standard deviations to estimate the uncertainties, such a way of dealing with the problem is acceptable.

6

Part of the scatter observed in H ii regions at a given galactocentric distance is likely due to abundance uncertainties, as discussed in Sect. 6.

Acknowledgments

M.P. is grateful to DAS, Universidad de Chile, for hospitality during a sabbatical stay when part of this work was done. This work received financial support from grant IN-105511(PAPIIT DGAPA-UNAM). F.B. gratefully acknowledges partial support from the National Science Foundation grants AST-0707911 and AST-1008798. Y.G.T. acknowledges the award of a Marie Curie intra-European Fellowship within the 7th European Community Framework Programme (grant agreement PIEF-GA-2009-236486).

References

  1. Aggarwal, K. M., & Keenan, F. P. 1999, ApJS, 123, 311 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  2. Blöcker, T. 1995, A&A, 299, 755 [NASA ADS] [Google Scholar]
  3. Bresolin, F. 2011, ApJ, 730, 129 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bresolin, F., Schaerer, D., González Delgado, R. M., & Stasińska, G. 2005, A&A, 441, 981 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bresolin, F., Gieren, W., Kudritzki, R.-P., et al. 2009, ApJ, 700, 309 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bresolin, F., Stasińska, G., Vílchez, J. M., Simon, J. D., & Rosolowsky, E. 2010, MNRAS, 404, 1679 [NASA ADS] [Google Scholar]
  7. Copetti, M. V. F., & Writzl, B. C. 2002, A&A, 382, 282 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Decressin, T., Charbonnel, C., Siess, L., et al. 2009, A&A, 505, 727 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Froese Fischer, C., Tachiev, G., & Irimia, A. 2006, ADNDT, 92, 607 [CrossRef] [Google Scholar]
  11. Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1997, A&AS, 123, 159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. García-Rojas, J., Peña, M., Morisset, C., Mesa-Delgado, A., & Ruiz, M. T. 2012, A&A, 538, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Gieren, W., Pietrzyński, G., Soszynski, I., et al. 2005, ApJ, 628, 695 [NASA ADS] [CrossRef] [Google Scholar]
  14. Gogarten, S. M., Dalcanton, J. J., Williams, B. F., et al. 2010, ApJ, 712, 858 [NASA ADS] [CrossRef] [Google Scholar]
  15. Henry, R. B. C., Kwitter, K. B., & Bates, J. A. 2000, ApJ, 531, 928 [NASA ADS] [CrossRef] [Google Scholar]
  16. Henry, R. B. C., Speck, A., Karakas, A. I., Ferland, G. J., & Maguire, M. 2012, ApJ, 749, 61 [NASA ADS] [CrossRef] [Google Scholar]
  17. Hernández-Martínez, L., Peña, M., Carigi, L., & García-Rojas, J. 2009, A&A, 505, 1027 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Hummer, D. G., & Storey, P. J. 1987, MNRAS, 224, 801 [NASA ADS] [CrossRef] [Google Scholar]
  19. Izotov, Y. I., Stasińska, G., Meynet, G., Guseva, N. G., & Thuan, T. X. 2006, A&A, 448, 955 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Karakas, A. I. 2010, MNRAS, 403, 1413 [NASA ADS] [CrossRef] [Google Scholar]
  21. Karakas, A. I., & Lattanzio, J. C. 2003, PASA, 20, 393 [Google Scholar]
  22. Karakas, A. I., van Raai, M. A., Lugaro, M., Sterling, N. C., & Dinerstein, H. L. 2009, ApJ, 690, 1130 [NASA ADS] [CrossRef] [Google Scholar]
  23. Kingsburgh, R. L., & Barlow, M. J. 1994, MNRAS, 271, 257 [NASA ADS] [CrossRef] [Google Scholar]
  24. Kisielius, R., Storey, P. J., Ferland, G. J., & Keenan, F. P. 2009, MNRAS, 397, 903 [Google Scholar]
  25. Kudritzki, R.-P., Urbaneja, M. A., Bresolin, F., et al. 2008, ApJ, 681, 269 [NASA ADS] [CrossRef] [Google Scholar]
  26. Lagarde, N., Decressin, T., Charbonnel, C., et al. 2012, A&A, 543, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Leisy, P., & Dennefeld, M. 2006, A&A, 456, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Maciel, W. J., & Koppen, J. 1994, A&A, 282, 436 [NASA ADS] [Google Scholar]
  29. Maciel, W. J., Costa, R. D. D., & Uchida, M. M. M. 2003, A&A, 397, 667 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Magrini, L., Vílchez, J. M., Mampaso, A., Corradi, R. L. M., & Leisy, P. 2007, A&A, 470, 865 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Magrini, L., Stanghellini, L., & Villaver, E. 2009, ApJ, 696, 729 [NASA ADS] [CrossRef] [Google Scholar]
  32. Magrini, L., Stanghellini, L., Corbelli, E., Galli, D., & Villaver, E. 2010, A&A, 512, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Marigo, P. 2001, A&A, 370, 194 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. McCall, M. L., Rybski, P. M., & Shields, G. A. 1985, ApJS, 57, 1 [NASA ADS] [CrossRef] [Google Scholar]
  36. McLaughlin, B. M., & Bell, K. L. 1998, JPhB, 31, 4317 [Google Scholar]
  37. McLaughlin, B. M., & Bell, K. L. 2000, JPhB, 33, 597 [Google Scholar]
  38. McLaughlin, B. M., Lee, T.-G., Ludlow, J. A., et al. 2011, JPhB, 44, 175206 [NASA ADS] [Google Scholar]
  39. Mendoza, C., & Zeippen, C. J. 1982, MNRAS, 198, 127 [NASA ADS] [Google Scholar]
  40. Munoz Burgos, J. M., Loch, S. D., Ballance, C. P., & Boivin, R. F. 2009, A&A, 500, 1253 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Paturel, G., Petit, C., Prugniel, P., et al. 2003, A&A, 412, 45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Peimbert, M. 1975, ARA&A, 13, 113 [NASA ADS] [CrossRef] [Google Scholar]
  43. Peimbert, M. 1990, RPPh, 53, 1559 [Google Scholar]
  44. Peña, M., Stasińska, G., & Medina, S. 2001, A&A, 367, 983 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Peña, M., Stasińska, G., & Richer, M. G. 2007, A&A, 476, 745 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Peña, M., Reyes-Pérez, J., Hernández-Martínez, L., & Pérez-Guillén, M. 2012, A&A, 547, 78 [Google Scholar]
  47. Porter, R. L., Ferland, G. J., & MacAdam, K. B. 2007, ApJ, 657, 327 [NASA ADS] [CrossRef] [Google Scholar]
  48. Ramsbottom, C. A., Bell, K. L., & Keenan, F. P. 1997, MNRAS, 284, 754 [NASA ADS] [Google Scholar]
  49. Renzini, A., & Voli, M. 1981, A&A, 94, 175 [NASA ADS] [Google Scholar]
  50. Rodríguez, M., & Delgado-Inglada, G. 2011, ApJ, 733, L50 [NASA ADS] [CrossRef] [Google Scholar]
  51. Rosolowsky, E., & Simon, J. D. 2008, ApJ, 675, 1213 [NASA ADS] [CrossRef] [Google Scholar]
  52. Searle, L. 1971, ApJ, 168, 327 [NASA ADS] [CrossRef] [Google Scholar]
  53. Seaton, M. J. 1979, MNRAS, 187, 73 [Google Scholar]
  54. Shields, G. A. 1974, ApJ, 193, 335 [NASA ADS] [CrossRef] [Google Scholar]
  55. Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef] [Google Scholar]
  56. Soffner, T., Mendez, R. H., Jacoby, G. H., et al. 1996, A&A, 306, 9 [NASA ADS] [Google Scholar]
  57. Stanghellini, L., & Haywood, M. 2010, ApJ, 714, 1096 [NASA ADS] [CrossRef] [Google Scholar]
  58. Stasińska, G. 2005, A&A, 434, 507 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Stasińska, G., Gorny, S. K., & Tylenda, R. 1997, A&A, 327, 736 [NASA ADS] [Google Scholar]
  60. Stasińska, G., Richer, M. G., & McCall, M. L. 1998, A&A, 336, 667 [NASA ADS] [Google Scholar]
  61. Stasińska, G., et al. 2012, Oxygen in the Universe, EAS Pub. Ser., 54 (EDP Sciences) [Google Scholar]
  62. Storey, P. J., & Hummer, D. G. 1995, MNRAS, 272, 41 [NASA ADS] [CrossRef] [Google Scholar]
  63. Storey, P. J., & Zeippen, C. J. 2000, MNRAS, 312, 813 [Google Scholar]
  64. Tayal, S. S. 2011, ApJS, 195, 12 [NASA ADS] [CrossRef] [Google Scholar]
  65. Tayal, S. S., & Gupta, G. P. 1999, ApJ, 526, 544 [NASA ADS] [CrossRef] [Google Scholar]
  66. Tayal, S. S., & Zatsarinny, O. 2010, ApJS, 188, 32 [NASA ADS] [CrossRef] [Google Scholar]
  67. Urbaneja, M. A., Herrero, A., Bresolin, F., et al. 2005, ApJ, 622, 862 [NASA ADS] [CrossRef] [Google Scholar]
  68. Vila-Costas, M. B., & Edmunds, M. G. 1992, MNRAS, 259, 121 [NASA ADS] [CrossRef] [Google Scholar]
  69. Vilchez, J. M., Pagel, B. E. J., Diaz, A. I., Terlevich, E., & Edmunds, M. G. 1988, MNRAS, 235, 633 [NASA ADS] [CrossRef] [Google Scholar]
  70. Wang, W., & Liu, X.-W. 2007, MNRAS, 381, 669 [NASA ADS] [CrossRef] [Google Scholar]
  71. Wang, W., Liu, X.-W., Zhang, Y., & Barlow, M. J. 2004, A&A, 427, 873 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Weidemann, V. 2000, A&A, 363, 647 [NASA ADS] [Google Scholar]
  73. Wesson, R., Stock, D. J., & Scicluna, P. 2012, MNRAS, 422, 3516 [NASA ADS] [CrossRef] [Google Scholar]
  74. Zeippen, C. J. 1982, MNRAS, 198, 111 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Main properties of NGC 300.

Table 2

Summary of spectroscopic observationsa.

Table 3

Line intensities for H ii regions and PNe, relative to Hβ = 1.00123.

Table 4

Atomic data used in our abundance determinations.

Table 5

Electron densities and temperatures for the PNe and compact H ii regions in NGC 300.

Table 6

Ionic abundances for the PNe and compact H ii regions in NGC 300.

Table 7

Computed abundances for the PNe and H ii regions in NGC 300a.

Table 8

Straight-line fitting of the radial variations of abundance ratios for the PNe and H ii regions in NGC 300.

Table 9

Typical properties of PNe and their progenitors.

All Figures

thumbnail Fig. 2

Values of the density indicator [S ii] λ6731/6717 as a function of [O ii] λ3726/3729. Points with error bars: observed values in the sample of PNe by Wang et al. (2004). Continuous lines: lines of equal electron density at a temperature of 104 K from ne = 20 cm-3 (upper right) to ne = 105 cm-3 (lower left). Red: using the [S ii] transition probabilities by Tayal & Zatsarinny (2010). Blue: using the [S ii] transition probabilities from Mendoza & Zeippen (1982).

In the text
thumbnail Fig. 3

Logarithms of the emissivities of the [O ii] λ3727, [N ii] λ6584, [S ii] λ6716, λ6731, and [Ar iii] λ7135 lines for densities ne = 102, 103, 104 and 105 cm-3. Red: with the atomic data used in the present work. Blue: with the atomic data used in our previous work.

In the text
thumbnail Fig. 4

Various abundance ratios as a function of O++/(O+ + O++) for PNe and H ii regions in NGC 300. PNe are represented by circles, compact H ii regions by stars and giant H ii regions by squares. The big red star is object no. 39, the big black circle is object no. 74 (see Sect. 4), and the big black star is object no. 57 (see Sect. 3.3). The objects represented by open circles are objects 12, 88 and 109 (see Sect. 3.4).

In the text
thumbnail Fig. 5

Histogram of log  O/H + 12 values for object PN 48 from the Monte Carlo simulation. The nominal value is indicated by the continuous line, the 15 and 85 percentiles by the dashed lines.

In the text
thumbnail Fig. 6

Diagrams illustrating observed properties that allow to differentiate between PNe and H ii regions in NGC 300 (see Sect. 4). The symbols have the same meaning as in Fig. 4. The few PNe for which the electron density cannot be measured have been placed at log Mion = 3.

In the text
thumbnail Fig. 7

Abundance patterns of PNe and H ii regions in NGC 300. The symbols have the same meaning as in Fig. 4.

In the text
thumbnail Fig. 8

Abundances with respect to hydrogen as a function of the normalized galactocentric radius R/R25. The symbols have the same meaning as in Fig. 4. The continuous blue lines and the dashed red lines represent the straight line fits to the observations of H ii regions and of PNe, respectively.

In the text
thumbnail Fig. 9

Abundances with respect to oxygen as a function of the normalized galactocentric radius R/R25. The layout of the figure is the same as in Fig. 4.

In the text
thumbnail Fig. 10

Evolution of the PN toy model described in Sect. 7 for PNe having progenitor masses (from right to left) of 2 M (red), 2.5 M (green), 3 M (blue), 4 M (cyan), 5 M (magenta). For 1 M, the Hβ luminosity would be too small to be detected. The black box in the left panel indicates the observed range for our PNe. The right panel indicates the time evolution of L(Hβ) since the formation of the PN model. It illustrates that the chance to observe a PN having a progenitor mass of 2 M is much larger than for 2.5 M.

In the text
thumbnail Fig. 11

Observed and predicted radial trends of O/H and N/O in the PNe of NGC 300. The thin lines show the predicted trends assuming that the progenitor initial abundances follow the regression line for H ii regions (represented with the thick blue line). The different colors correspond to the different stellar evolution models listed in Table 9: green for Marigo (2001), grey for Karakas (2010), blue for the standard models of Lagarde et al. (2012), purple for the models with rotation of Lagarde et al. (2012). The thick dashed lines shows is regression line for the observed PN abundances.

In the text
thumbnail Fig. 1

Calibrated spectra of our observed PNe and compact H ii regions.

In the text
thumbnail Fig. 1

continued.

In the text
thumbnail Fig. 1

continued.

In the text
thumbnail Fig. 1

continued.

In the text
thumbnail Fig. 1

continued. Note that for PN 74, the red spectrum has been obtained with grism 300I, which has a much lower resolution than 600B and 600RI.

In the text
thumbnail Fig. 1

continued.

In the text
thumbnail Fig. 1

continued.

In the text
thumbnail Fig. 1

continued.

In the text

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.