Free Access
Issue
A&A
Volume 615, July 2018
Article Number A29
Number of page(s) 19
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201730861
Published online 06 July 2018

© ESO 2018

1 Introduction

As Sun-like stars (of ~0.8–8 M) evolve, they eventually become planetary nebulae (PNe). In this evolutionary phase, the star has lost enough of its convective envelope through stellar winds to expose its inner, hotter regions, causing this ejected material to become ionised. We can determine the ionic and elemental abundances of this PN ejecta by analysing the strong forbidden emission lines of the ions in the stellar spectrum.

Much of the observable PN emission comes from ionisation states of α-process elements, such as neon, sulphur, and argon. The main stages of ionisation for these three elements are observed at mid-infrared wavelengths. While the abundances of elements such as helium and carbon change significantly throughout the course of stellar evolution, those of neon, sulphur, and argon are unchanged (e.g. Marigo et al. 2003), making them useful probes of metallicity at the epoch of stellar formation. Abundance studies carried out at optical wavelengths commonly use oxygen as a metric for metallicity, as the observed emission lines of O+ and O2+ are always strong. However, during the evolution of these low- to intermediate-mass stars, the abundances of oxygen are known to change. This particularly occurs within the asymptotic giant branch (AGB) phase, in which the third dredge-up brings helium, carbon and a small amount of oxygen to the outer envelope of the star. For stars with M ≳ 4 M, some oxygen will also be destroyed by hot bottom burning (e.g. Karakas & Lattanzio 2014; Delgado-Inglada et al. 2015). An empirical study by Delgado-Inglada et al. (2015) has also shown that oxygen enrichment can occur in Galactic PNe with carbonaceous dust. It has been proposed by García-Hernández et al. (2016) that this can be explained by diffusive convective overshooting processes, in which core mixing is extended beyond the Schwarzschild boundary of main sequence stars (Böhm-Vitense 1958; Herwig et al. 1997), producing significant increases in oxygen abundances around sub-solar and solar metallicities (Marigo 2001; Pignatari et al. 2016).

While most abundances from observational studies have been measured using optical spectra, there are some advantages to analysing PNe using infrared spectra. These are described in several studies (e.g. Rubin et al. 1988; Pottasch & Beintema 1999; Bernard-Salas 2003), but can be summarised as follows: extinction corrections are greatly reduced at IR wavelengths compared to those in the optical and ultraviolet regions; many ionic emission lines are observable for Ne, S and Ar within this wavelength range, and hence the need for ionisation correction factors (ICFs) in calculating elemental abundances is reduced; as these IR lines also originate from energy levels close to the ground state, both the uncertainties in the electron temperatures of any ion measured at IR wavelengths and temperature fluctuations within the PN can have little effect on the overall abundances. For this study in particular, the extinction corrections are reduced further as we have analysed the Galactic anti-centre, a region with much less extinction than the bulge (e.g. Pottasch & Bernard-Salas 2015).

The presence of the Galactic metallicity gradient was made clear in a sample of H II regions by Shaver et al. (1983) for nitrogen, oxygen, sulphur, and argon. Since then, it has been further studied not only in H II regions (e.g. Martín-Hernández et al. 2002; Esteban et al. 2017; Fernández-Martín et al. 2017) but also in PNe (e.g. Maciel & Quireza 1999; Pottasch & Bernard-Salas 2006; Maciel et al. 2015), young B-type stars (e.g. Fitzsimmons et al. 1992; Rolleston et al. 2000), Cepheid variables (e.g. Andrievsky et al. 2002a,b,c; Luck et al. 2003; Lemasle et al. 2013; Genovali et al. 2015), open clusters (e.g. Friel 1995) and young stars (e.g. Magrini et al. 2017).

While the presence of the metallicity gradient in the Galactic disk is agreed upon over Galactocentric distances (Rg) in the range4–10 kpc, its continuation towards the anti-centre is debated. Studies of H II regions (Esteban et al. 2017; Fernández-Martín et al. 2017) and B-type stars (Smartt 2000) have found that there is little variation in the gradient far from the Galactic centre, yet a study of Cepheid variables from Andrievsky et al. (2002c) showed the gradient flattening with Rg. Samples of PNe have also been previously analysed in the anti-centre with conflicting results. Costa et al. (2004) showed that the oxygen abundances of a group of PNe, 8–15 kpc away from the Galactic centre, did not directly follow the gradient but instead flattened beyond 10 kpc. This has also been observed in the nearby spiral galaxies M31, M33, M81 and NGC 300 (Magrini et al. 2016). However, the sample of Henry et al. (2010) suggested that the gradient steepened beyond this distance.

Chemical evolution models of the Milky Way have predicted that the radial abundance gradient will flatten over time due to several factors, such as the death of massive stars, which causes the metallicity to increase over time (e.g. Minchev et al. 2013), and radial migration (e.g. Minchev et al. 2012, 2014; Vera-Ciro et al. 2014; Kubryk et al. 2015), in which the angular momentum from stars is redistributed, leading to the movement of stars from the Galactic disk and hence contributing to a flattening radial metallicity gradient within the disk (e.g. Sellwood & Binney 2002).

Investigations into the time evolution of the radial metallicity gradient have given varying results, with several studies of PNe finding an overall steepening with time (e.g. Maciel & Quireza 1999; Chiappini et al. 2001; Stanghellini & Haywood 2010; Kubryk et al. 2015), suggesting that the Galactic disk formed slowly (Chiappini et al. 1997). However, Maciel et al. (2003) showed the gradient flattening over time. Studies of open clusters and field stars have also given varying conclusions on this matter (Anders et al. 2017).

In this paper, we have derived the abundances of neon, sulphur, and argon in a sample of 23 PNe located towards the Galactic anti-centre using IR data in order to study the metallicity gradient beyond 10 kpc, and compared them to other IR spectroscopic samples from the Milky Way that were analysed in the same way.

The layout of this paper is as follows: In Sect. 2 we discuss the source selection and the basic data reduction and extraction methods. In Sect. 3 we explain the methods used to calculate flux, intensity and abundance values as well as Galactocentric distances. The implications of these data on the metallicity gradient in the further regions of the Milky Way are considered in Sect. 4. Finally, in Sect. 5 we present our conclusions and summarise our main results.

2 Data

2.1 Observations

The observations were made with the Infrared Spectrograph on board the Spitzer Space Telescope (Spitzer IRS; Werner et al. 2004; Houck et al. 2004) through GTO programme 40035 (PI: J. Bernard-Salas). The observations were carried out between December 2007 and December 2008 with the staring mode of the IRS using the short-low (SL), short-high (SH) and long-high (LH) modules, each allowing for simultaneous observations from two nod positions, at 1/3 and 2/3 of the way along the observing slit. These produced spectra with resolutions of R  ~ 60–127 in the range 5.2–14.5 μm, and of R  ~ 600 in the range 9.9–36.4 μm.

The sources are listed in Table A.1. These were chosen according to the following two criteria: (a) the sources were located in the direction of the anti-centre (l = 120–240°, b = 0°± 20°); (b) the physicalsizes of the PNe were generally small enough to fit in the widest observing slit of Spitzer IRS (LH, 11.1′′  × 22.3′′), thus minimising the aperture corrections required to account for the different slit sizes.

2.2 Data reduction and extraction

The basic calibration data (bcd) image files obtained from the Spitzer IRS were processed through the Spitzer Science Centre (SSC) pipeline, version S18.18, then reduced and analysed through the Spectroscopic Modelling Analysis and Reduction Tool (SMART; Higdon et al. 2004). Rogue pixels were removed using the IRSCLEAN1 package.

The methods of spectral extraction were based on the diameters of the sources (see Table A.2). For the low-resolution SL module, sources with diameters ≤3′′ were extracted using the Advanced Optimal Extraction (AdOpt) package (Lebouteiller et al. 2010). This method weights each pixel based on their signal-to-noise ratios and is better suited to smaller sources. Those with diameters in the range 3′′ –8′′ were extracted using tapered column extractions, as the SL slit is not wide enough (3.6′′   × 57′′) to detect all the flux in partially extended sources. The four PNe with diameters ≥10′′ were extracted using fixed column extractions as the FWHM of their emission was beyond the point spread function of the objects by factors up to approximately three. For the high-resolution SH and LH modules, full aperture extractions were used in each case; these weight all pixels equally from the aperture, allowing for most of the flux to be obtained from more extended sources if they covered an area larger than the observing slits. The flux values of a PN with a physical size bordering on two of these ranges do not change significantly (≲10%).

Most of the PNe in our sample were chosen such that all of the source flux would fall within the LH module. In fact, 15 of our 23 targets also have a physical diameter of ≤5′′, so in these cases most of the flux would also be detected by the smaller and narrower SH module (4.7′′   × 11.3′′). Some of the source flux might still not be detected by the SL and SH modules, resulting in jumps in the baseline continuum. To account for this, we scale the SL and SH flux values by matching the continua in the overlapping wavelength regions. The scale factors are listed in Table A.2. The full low- and high-resolution spectra of two representative PNe in the sample, M1-16 and M1-8, are shown in Fig. 1.

thumbnail Fig. 1

Full Spitzer IRS spectrum of M1-16 (top) and M1-8 (bottom). The low resolution spectra (SL) are shown on the left, and thehigh resolution spectra (SH and LH) are on the right.

3 Analysis

3.1 Ionic abundances

The ionic abundances obtained from the spectra have been calculated using the equation nionnp=neIionIHβλulλHβαHβAul(nunion)1,\begin{equation*} \frac{n_{ion}}{n_{p}} = n_{e} \frac{I_{ion}}{I_{H\beta}} \frac{\lambda_{ul}}{\lambda_{H\beta}} \frac{\alpha_{H\beta}}{A_{ul}} \Bigg(\frac{n_{u}}{n_{ion}}\Bigg)^{-1},\vspace*{-4pt}\end{equation*}(1)

where np is the proton density, ne is the electron density, Iion is the intensity of the ion, λul is the line wavelength, α is the effective recombination coefficient for , Aul is the Einstein coefficient of spontaneous emission, and nunion is the ratio of the upper level population of the transition to the entire population of the ion.

For the derivation of the neon, sulphur, and argon abundances, we have observed the emission lines for the most populated ions of these elements: Ne+, Ne2+, Ne4+, S2+, S3+, Ar+, Ar2+ and Ar4+. We complemented our IR ionic abundances with those of the missing ionisation states from optical spectra (Ne3+, S+ and Ar3+) in the literature in order to avoid or reduce the need for ICFs. These values were primarily taken from Henry et al. (2010) and references included in Sterling & Dinerstein (2008). We compared IR and optical data, so homogeneity in the slit sizes was assumed. This is a reasonable assumption as Henry et al. (2010) use data from the Apache Point Observatory (APO), for which the slit size is 2′′   ×   360′′, and most of the PNe they have observed from this sample are ≲4′′ in diameter, so most of the flux will have been detected from these sources. We also considered S4+, which is not directly detectable in IR or optical spectra, although it is only expected to contribute to the overall sulphur abundance for PNe with high IP values. Where there are no flux values available for any particular line in a PN, we applied correction factors (see Sect. 3.5).

3.2 Line flux measurements

Flux values for the fine structure ionic emission lines with a ≥3σ detection were determined by applying Gaussian fits through the ISAP line fitting programme (Sturm et al. 1998) in SMART. The raw F(λ) values were calculated for each line in each of the two observing positions (nods), which typically agreed to within 10%. These were then averaged; in the few cases where there were any remaining low-level glitches obstructing a particular line, these were discarded infavour of the flux value from the nod without glitches. The associated uncertainties were propagated from those calculated from the individual flux measurements within each nod, unless the difference in flux between the two nods was greater than the assumed uncertainties. In this latter case, we considered the flux difference to be more representative of the uncertainty. The data were then corrected for extinction using the extinction law from Fluks et al. (1994). Table A.3 shows all values for the extinction corrected intensities, I(λ). Upper limits of 3σ or more were calculated for emission lines of the most important ions of neon, sulphur, and argon when there were no clear detections.

3.3 Hβ intensity

In the wavelength range of Spitzer IRS we observed the emission features of several recombination transitions of atomic hydrogen (H I), the strongest of which were observed at 7.5 μm and 12.4 μm. Both of these emission lines account for at least two transitions; the H I 6–5, 8–6, 11–7 and 17–8 lines are blended around 7.5 μm (H I 6–5 is the strongest of these transitions, contributing 74.43% of the total flux) whereas the H I 7–6 and 11–8 transitions both contribute to the emission line at 12.4 μm (H I 7–6 provides 89.08% of this flux). We applied the Balmer decrement to obtain values of I() from Hummer & Storey (1987), interpolated to account for the electron density and temperature values of our PNe. When both of these IR emission features were observed, their I() values agreed by up to ~25%, and an average was taken.

Our calculated I() values are shown in Table A.4 alongside the F() and I() values from 4861 Å optical line measurements and the extinction coefficients, C(), all taken from literature. In the four cases when neither of the two recombination lines were observed in a spectrum, we applied the Fluks extinction law (Fluks et al. 1994) to these literature F() values. Our values agree with those in literature mostly to within a factor of two (more for K3-90 and M1-7). In these situations we favour our IR values, as the H I lines are measured in the same spectra as the ionic emission that we have derived. Another advantage of using these lines to determine F() is that the extinction corrections are far smaller than those from optical wavelengths; Aλ   <  4.6 from the use of the 4861 Å line, whereas from the IR recombination lines we find that Aλ   <  0.2.

3.4 Electron densities and temperatures

Both ne and Te are needed to determine abundance values; ne is a direct component of their calculation (see Eq. (1)) whereas Te designates the statistical populations of the excited electronic states present within the ion. These are listed in Tables A.5 and A.6 respectively, and the transition probabilities and collision strengths used in calculating these values are shown in Table A.7. These values were taken from TIPbase, part of the IRON project (Hummer et al. 1993).

Infraredlines originate from electronic transitions close to the ground state. Therefore, by analysing the ratio of I(λ) values for two transitions of the same ion, wewere able to obtain ne values that are mostly independent of temperature. Of the line flux ratios available from our spectra, we favoured those of the [S III] 18.7/33.5 μm transitions as both of these lines are easy to measure in high resolution spectra and frequently seen together. While otherline ratios were available in some PNe (e.g. [Ne III], [Ne V]), they were either detected in fewer of the PNe in our sample, or were detected at noisy wavelength regions. For instance, the [Ne III] line at 36.0 μm is found at the upper wavelength region of the LH module, which is highly susceptible to noise above ~35 μm. The 14.3 and 24.3 μm lines of [Ne V] were only observed in 12 of the 23 PNe in the sample, though the associated density values agree well with those of [S III]. Uncertainties averaged ~20% for values of ne > 1000 cm−3, though thisbecomes larger for the few sources where ne < 1000 cm−3. In the four cases where the two [S III] lines were not directly measurable and any other line ratios were either not observed or affected by noise, we used values given in the literature from the [S II] 6716 Å / 6731 Å line intensity ratio. We applied the mean value derived from the [S III] line ratios in our sample of ne = 3700 cm−3 for K3-69 and Y-C 2-5 as these lines were not observed in these PNe and there were no ne values givenin literature. In these cases, while the uncertainty in density is high, the abundances are little affected, with neon and argon showing little change at the density extremes, and sulphur being affected by 20% at most. All density values are shown in Table A.5, with the uncertainties reflecting those of the [S III] 18.7/33.5 μm ratios. We note that for K3-90, Henry et al. (2010) apply the high density limit to estimate ne despite having intensity values for the [S II] 6716 Å and 6731 Å lines. This is due to the two values having high uncertainties. A density of ne = 800 cm−3 would have been calculated with these values, which is almost within the error margins of our IR [S III] line ratio.

Measurements of Te require electronic transitions with large differences in energy. For this study, we relied on temperatures calculated from optical line flux values based on the transition ratios of the [N II], [O II], [O III], [S II] and [S III] lines calculated from the literature. As no Te values could be found in the literature for K3-65 and K3-69, for these two PNe we adopted the average value of Te = 11 900 K. All of these values can be found in Table A.6.

3.5 Elemental abundances

One of the main advantages of analysing spectra at infrared wavelengths is that the main ionisation lines of neon, sulphur, and argon can be observed. From these lines, we have measured the ionic abundances of Ne+, Ne2+, Ne4+, S2+, S3+, Ar+, Ar2+ and Ar4+. We complemented these data with the ionic abundances of S+ and Ar3+ measured by Henry et al. (2010) from optical spectra, hence fewer corrections are required in determining their elemental abundances. We accounted for Ne3+ in sources with observable Ne4+ emission, and we considered S4+ in sources with O3+, which has a greater IP (47.22 eV and 54.94 eV, respectively).

We corrected for these missing ionic abundances with ICFs. ICFs can either be determined empirically (e.g. Surendiranath et al. 2004; Pottasch & Surendiranath 2005), by considering lines with similar IP values (e.g. Peimbert & Costero 1969), or from photoionisation models (e.g. Natta et al. 1980; Kingsburgh & Barlow 1994; Kwitter & Henry 2001; Delgado-Inglada et al. 2014). In many cases, argon and particularly sulphur are highlighted as being complicated to correct for, as the low IPs of higher ionisation states may lead to their greater contributions towards the overall elemental abundances. Many variants of the ICFs for these elements have been given in literature (e.g. Kingsburgh & Barlow 1994; Thuan et al. 1995; Kwitter & Henry 2001) with significant disagreement between some of them (Vermeij & van der Hulst 2002).

To account for the abundances of missing ionisation states, we complemented our IR values with optical abundances derived by Henry et al. (2010) where possible. To correct for S+ and Ar3+ when these values are not available, we calculated the percentage contributions of these ions towards their respective elemental abundances for the PNe in our sample with these ionic abundances and applied the mean values as ICFs. In each case, the minimum and maximum values were taken as the uncertainty limits. We also applied this method to account for Ar+ in three sources for which the 7.0 μm line intensity cannot be measured. For Ne3+ and S4+, we considered the range of contributions of these ions to their respective elemental abundances as given by Bernard-Salas et al. (2008), who calculated these from the analysis of the PN sample of Pottasch & Bernard-Salas (2006) and also from the Galactic PN models of Surendiranath et al. (2004) and Pottasch & Surendiranath (2005). Each of these ranges are shown in Table A.8. We note that we only applied ICFs correcting for these missing ionic states in the PNe for which we observed other ions with greater or similar IP values.

Tables A.9A.11 give the ionic and elemental abundance values for Ne, S and Ar respectively, with the empirical ICFs applied. ICFs can be uncertain, so we also compared the resulting abundances calculated with our empirical ICFs with those calculated using well-established ICFs from the literature, such as those of Kingsburgh & Barlow (1994), Kwitter & Henry (2001) and Delgado-Inglada et al. (2014). In all cases, we applied our ionic abundances. For the equations given in this section, we applied the notation A(X) = ICF(Xm+ + Xn+) × (Xm+ + Xn+)/H, where A(X) is the elemental abundance of X. Table A.12 shows the supplementary data used in the following calculations from optical abundance studies. We note that some ICF equations from other studies use ionic ratios that are particularly sensitive to the electron density and temperature, such as those involving O+ and O2+. Hence, in cases where the [S III] ion ratio does not adequately account for the electron density of an important ion for the ICF calculation, or the uncertainties in Te are large due to the dispersion of values over several studies, there may be additional uncertainty in the abundances calculated using these ICFs.

3.5.1 Neon

At mid-IR wavelengths, lines of Ne+, Ne2+, Ne4+ and Ne5+ can be measured, though Ne3+ is best observed in the optical and near-UV regions, respectively at ~4720 Å and ~2424 Å. Unfortunately, no literature values exist for the abundance of Ne3+ in any of the PNe in the anti-centre sample.

The photoionisation model of Kingsburgh & Barlow (1994), also used by Kwitter & Henry (2001), used the following ICFs: ICF(Ne2++Ne4+) = 1.5,\begin{equation*} \textrm{ICF}(\textrm{Ne}^{2&#x002B;} &#x002B; \textrm{Ne}^{4&#x002B;})~=~1.5,\end{equation*}(2) ICF(Ne2+) = OO2+.\begin{equation*} \textrm{ICF}(\textrm{Ne}^{2&#x002B;})~=~\frac{\textrm{O}}{\textrm{O}^{2&#x002B;}}.\end{equation*}(3)

However, the abundances of sources with weak radiation fields are underestimated due to the disregard of the Ne+ ionic contribution to the total neon abundance (e.g. Tsamis et al. 2013). This problem is also observed in the ionic abundances of several sources in our sample of anti-centre PNe, in which Ne+ sometimes contributes more to the overall elemental abundance than Ne2+. These include K3-66, M1-6 and M1-14 (see Table A.9).

Recently, Delgado-Inglada et al. (2014) produced a newer set of ICF models to account for parameters such as effective temperatures and stellar luminosities. For Ne/H corrections, they apply: ICF(Ne2++Ne4+) = (1.31+12.68ν2.57)0.27\begin{equation*} \textrm{ICF}(\textrm{Ne}^{2&#x002B;} &#x002B; \textrm{Ne}^{4&#x002B;})~=~(1.31 &#x002B; 12.68\nu^{2.57})^{0.27}\end{equation*}(4)

where ν = He2+/(He+ + He2+). However, the requirement for Ne4+ limits the usability of this correction. They also state that the ICFs will overestimate the neon abundances unless 0.4  ≲ ν  ≲ 0.6. The He+ and He2+ abundances given for 12 of the 23 PNe in the anti-centre sample from Henry et al. (2010) give ν values outside this range, ten of these sources having values of ν < 0.2 and 5 of these with ν < 0.005. PNe with very low ν values have small He2+ ionic abundances, which typically indicates low-ionisation sources with little or no Ne4+ emission. However, this is not true for Y-C 2-5, which has a relatively large ν (0.61) but no observable Ne4+ emission in its Spitzer IRS spectrum.

Table A.13 shows a comparison between the neon abundances calculated with both our empirical ICFs and the well-established ICFs of Kingsburgh & Barlow (1994) and Delgado-Inglada et al. (2014). There is good agreement in almost all cases between the two sets of values, though the disregard of Ne+ leads to major underestimates from the Kingsburgh & Barlow (1994) model for M1-6 and M1-14. For Y-C 2-5, the applied ICF of 2.58 leads to a much greater abundance than predicted empirically. However, as the mid-IR spectrum of this PN shows the strong emission line of Ne2+ but not those of Ne+ or Ne4+, it is possible that the Ne3+ for which we are correcting may give a significant contribution to the neon abundance. The fact that we see large amounts of He2+ compared to He+ in this source shows that the radiation field in Y-C 2-5 is greater than 54.4 eV (the IP of He2+), and as Ne3+ ionises at 63.5 eV, a large ICF may be required.

3.5.2 Sulphur

As ionisation models are typically applied to optical spectra, it is normal to only see corrections for S3+ from S+ and S2+ ionic abundance measurements. In these studies, the low IP of S2+ (34.8 eV) has always been taken as an indication that larger ionic states are likely to be present in PNe, and the similar IP of O+ (35.1 eV) is generally considered in obtaining an ICF for sulphur. Dinerstein (1980) carried out an IR spectroscopic survey of 12 PNe and found that the commonly used ICF of O/O+ can overpredict the measured abundances of S3+. However, the presence of even greater ionisation states must also be considered.

Based on models of H II regions from Stasińska (1978), the ICF for sulphur from optical spectra was calculated by Barker (1980) to be ICF(S++S2+)=[1(1O+O)α]1/α\begin{equation*} \textrm{ICF}(\textrm{S}^{&#x002B;} &#x002B; \textrm{S}^{2&#x002B;}) = \bigg[1 - \bigg(1 - \frac{\textrm{O}^{&#x002B;}}{\textrm{O}}\bigg)^{\alpha}\bigg]^{-1/\alpha}\end{equation*}(5)

where α = 3, though subsequent studies argued that α = 2 (French 1981) or 2 ≤ α ≤ 3 (Garnett 1989) better represented the sulphur abundances. Kingsburgh & Barlow (1994) also used this equation, with α = 3.

A different method of determining the ICF for sulphur was calculated by Kwitter & Henry (2001), who considered newer atomic data and incorporated the charge exchange rates into the ICF values. They used the equation ICF(S++S2+)=exp[0.017+0.18β0.11β2+0.072β3]\begin{equation*} \textrm{ICF}(\textrm{S}^{&#x002B;} &#x002B; \textrm{S}^{2&#x002B;}) = \textrm{exp}[-0.017 &#x002B; 0.18\beta - 0.11\beta^{2} &#x002B; 0.072\beta^{3}]\end{equation*}(6)

where β = log(O/O+).

The models of Delgado-Inglada et al. (2014) calculate ICFs for the S/O ratio, and multiply by the O/H abundance: logICF((S++S2+)/O+)=0.020.03ω2.31ω2+2.19ω30.69+2.09ω2.69ω2,\begin{equation*} \log \textrm{ICF}((\textrm{S}^{&#x002B;} &#x002B; \textrm{S}^{2&#x002B;})/\textrm{O}^{&#x002B;}) = \frac{-0.02 - 0.03\omega - 2.31\omega^{2} &#x002B; 2.19\omega^{3}}{0.69 &#x002B; 2.09\omega - 2.69\omega^{2}}, \end{equation*}(7) logICF(O++O2+)=0.08ν+0.006ν20.340.27ν\begin{equation*} \log \textrm{ICF}(\textrm{O}^{&#x002B;} &#x002B; \textrm{O}^{2&#x002B;}) = \frac{0.08\nu &#x002B; 0.006\nu^{2}}{0.34 - 0.27\nu}\end{equation*}(8)

where ν is as defined in Eq. (4) and ω = O2+/(O+ + O2+). In all cases, ω > 0.5 and for K3-67, M2-2 and Y-C 2-5, ω > 0.95.

The sulphur abundances calculated using these ICFs are shown in Table A.14. Again, the values calculated using empirical ICFs compare well with those from photoionisation models. ICFs obtained with Eq. (5) are often large when α = 2, and provide much greater estimates than those of the compared studies. However, the agreement is greatly improved when α = 3. The uncertainties in the ICF calculated from Kwitter & Henry (2001) are larger due to their propagation, but the ICFs themselves are smaller, with most of them having values of ≤1.35. The only exception to this is K3-90, which has a stronger radiation field than the others (S3+/S2+ = 5.4). Despite these uncertainties, the S/H values calculated from Eq. (6) show a good agreement with the other values. The abundances calculated using ICFs from Delgado-Inglada et al. (2014) are similar in that their uncertainties are relatively high, though in almost all cases they show better agreement with the abundances calculated using empirically determined ICFs than those from Kwitter & Henry (2001).

3.5.3 Argon

In optical spectra, the Ar2+, Ar3+ and Ar4+ ions can allbe observed. The most abundant of these ions is thought to be Ar2+, though this largely depends on the radiation field of the source.

Argon abundances calculated by Kingsburgh & Barlow (1994) applied the following ICFs. ICF(Ar2++Ar3++Ar4+)=11(N+/N),\begin{equation*} \textrm{ICF}(\textrm{Ar}^{2&#x002B;} &#x002B; \textrm{Ar}^{3&#x002B;} &#x002B; \textrm{Ar}^{4&#x002B;}) = \frac{1}{1 - (\textrm{N}^{&#x002B;}/\textrm{N})},\end{equation*}(9) ICF(Ar2+)=1.87 ± 0.41.\begin{equation*} \textrm{ICF}(\textrm{Ar}^{2&#x002B;}) = 1.87~\pm~0.41.\end{equation*}(10)

Kwitter & Henry (2001) built upon Eq. (9) by considering the ICF when only the Ar2+ and Ar3+ ionic states could be observed: ICF(Ar2++Ar3+)=11(N+/N)×He++He2+He+.\begin{equation*} \textrm{ICF}(\textrm{Ar}^{2&#x002B;} &#x002B; \textrm{Ar}^{3&#x002B;}) = \frac{1}{1 - (\textrm{N}^{&#x002B;}/\textrm{N})} \times \frac{\textrm{He}^{&#x002B;} &#x002B; \textrm{He}^{2&#x002B;}}{\textrm{He}^{&#x002B;}}.\end{equation*}(11)

Delgado-Inglada et al. (2014) calculated their abundances in terms of Ar/O, before multiplying by the O/H abundance calculated with the ICF from Eq. (8): logICF(Ar2+O++O2+)=0.03ω0.40.3ω0.05(0.5<ω<0.95).\begin{equation*} \log \textrm{ICF}\bigg(\frac{\textrm{Ar}^{2&#x002B;}}{\textrm{O}^{&#x002B;} &#x002B; \textrm{O}^{2&#x002B;}}\bigg) = \frac{0.03\omega}{0.4 - 0.3\omega} - 0.05 \quad (0.5 < \omega < 0.95).\end{equation*}(12)

Table A.15 shows a comparison between the Ar/H abundances using ICFs from all studies. In general, the empirical ICFs agree well with those from each of the three comparative studies. The uncertainties in the abundances calculated from the ICFs of Delgado-Inglada et al. (2014) are large, particularly for PNe with strong ionisation fields (ω ≳ 0.95). This is likely due to the ICF only requiring the ionic abundance of Ar2+.

Overall, the abundances calculated with empirically determined ICFs in all three considered elements compare well with those using ICFs from the literature. The main advantage of using the empirical method is that there is no need for large ICFs due to the number of ionisation states for which we have data. While much larger ICFs are applied in the literature toaccount for missing ionisation states that contribute more to the elemental abundances, the resulting values relate well to those calculated empirically in which some or all of these missing ionisation states have been observed.

3.6 Galactocentric distances

Distances toGalactic PNe are known to be notoriously difficult to measure owing to the variation of bolometric luminosity and effective temperature values within the sources. We considered distances relative to the Galactic centre and exclusively observed PNe towards the anti-centre, which somewhat reduces the errors relative to those of their respective heliocentric distances (Rh).

We converted from Rh to Rg using the following equation: |Rg|=([Rhcos(b)cos(l)R]2+Rh2cos2(b)sin2(l)+Rh2sin2(b))0.5 \begin{multline*} |\overrightarrow{R_{g}}| = ([R_{h} \cos(b) \cos(l) - R_{\odot}]^{2} &#x002B; R_{h}^{2} \cos^{2}(b) \sin^{2}(l) \\ &#x002B; R_{h}^2 \sin^{2}(b))^{0.5}\end{multline*}(13)

where we take R = 8.0 ± 0.5 kpc.

Table A.1 shows the distance values used for each PN. In this paper, we primarily adopt values from Frew et al. (2016) as measured through statistical means, though where possible we prioritised the use of directly determined distances from Giammanco et al. (2011). Table A.16 shows the Rh values of the PNe in the anti-centre sample with distance values given by Giammanco et al. (2011), alongside those from Frew et al. (2016). While there are strong disagreements in most cases, the abundances towards the anti-centre remain lower than elsewhere inthe Galactic disk. The choice of data set will not affect the overall conclusions as each case presents a similar levelof dispersion around the overall gradient, part of which comes from uncertainties in the distance measurements.

The uncertainties quoted by Frew et al. (2016) are estimated to be ~20–30%, which seem small for statistical values, though larger uncertainties will not affect the outcomes of our discussion.

4 Discussion

4.1 Comparison of abundances with literature

In Tables A.9A.11, we have compared the elemental abundances of neon, sulphur, and argon with those available from literature. We accounted for average uncertainties of 17.5% in neon abundances for the missing Ne3+ ionisation state, 25% in sulphur abundances for S+ and S4+, and 30% in argon for Ar3+ and occasionally Ar+. We used thesample of Henry et al. (2010) as our main comparison, as this is a recent study involving 12 of our 23 PNe, though we also considered the abundances shown within the study of Sterling & Dinerstein (2008). These all involved the use of optical spectra.

In most cases, there is good agreement between all sets of abundances. Where we have lower abundance values than in literature, this likely comes from over-estimated ICFs from optical studies; for example, the neon abundance of K3-90 from literature (Table A.9) has an ICF of 23.8 (from Eq. (3)) applied to an uncorrected value relatively similar to ours (Henry et al. 2010), resulting in our elemental abundances disagreeing with theirs by a factor of eight. The neon ICF was high because the oxygen ICF was comparatively large (22.3), which was calculated from the ratio of (He+ + He2+)/He+ abundances (Kwitter & Henry 2001). Kingsburgh & Barlow (1994) give the same oxygen ICF to the power of 2/3, which would have decreased the ICF from 22.3 to 7.9 and the neon abundance to (4.7 ± 2.6) × 10−5. This is narrowly within the uncertainty range of our Ne/H value ((1.8 ± 0.8) × 10−5).

Conversely, our neon abundance of M4-18 is a factor of 20 greater than that given in literature by De Marco & Crowther (1999). Their abundance was calculated using the 12.8 μm [Ne II] line flux from Aitken & Roche (1982) (3.6 × 10−12 erg cm−2 s−1), which compares well with our value, 3.3 × 10−12 erg cm−2 s−1. However, the value of log I() used by De Marco & Crowther (1999) is lower than ours (−11.44 and −11.24, respectively), accounting for a discrepancy of factor 1.6. We have used I() and C() values from Henry et al. (2010), but we were unable to directly compare our neon abundances for this source as these authors were unable to observe the dominant Ne+ ion, and Ne2+ was barely observable, if at all (they gave its ionic abundance with a 90% uncertainty).

The source M1-1 also shows relatively large discrepancies between our abundance and those given in literature for neon and argon. By considering its high ionisation field, from which we can observe Ne5+, we expect that Ne3+ and Ar3+ will contribute significantly to their overall abundances. Aller et al. (1986) estimate the abundances of Ne2+, Ne3+ and Ne4+ in M1-1, obtaining values of a factor of approximately three lower than our values, with Ne+3 contributing 40% of the total of those three ions (taken from measurements of the 2424 Å [Ne IV] line). They also estimate Ar3+ to contribute 31% of the total argon abundance. ICFs of 1.12 and 2.1, for neon and argon respectively, are applied to their overall abundances.Even so, there are still significant discrepancies between our sets of abundance values. The abundances of the anti-centre PNe agree well with those from literature. For the few cases in which there are discrepancies between these values, they are likely to come from the larger ICFs used within the literature.

thumbnail Fig. 2

Abundance gradients of neon, sulphur, and argon in the Milky Way. The dashed lines represent the oxygen abundance gradient from within the Galactic disk with a slope of −0.085 dex/kpc (Pottasch & Bernard-Salas 2006), passing through the solar value at 8.0 kpc (Asplund et al. 2005). The solid lines represent the line of best fit in each plot, with gradients of −0.058 ± 0.021, −0.079 ± 0.012 and −0.062 ± 0.023 dex/kpc respectively.

4.2 The abundance gradient

Figure 2 shows the abundances of neon, sulphur, and argon plotted against Rg for the PNe in the anti-centre sample and for those of Pottasch & Bernard-Salas (2006) from the solar neighbourhood. Both of these samples were analysed and reduced in a similar way.

The abundances are lower than those of the Galactic bulge and the solar neighbourhood, and are consistent with a continuation of the metallicity gradient up to Rg ~ 20 kpc, albeit with a large dispersion within the data. However, when analysing the neon and argon anti-centre data separately from the solar neighbourhood data, there is no clear correlation, with Pearson correlation coefficients (RPCC) of ~ − 0.05 in each case. Hence, we cannot discern with our data whether there is a gradient in the anti-centre (Rg > 10 kpc). For sulphur, the anti-centre data show RPCC = − 0.45 with a corresponding p-value of 0.032, showing that there is a slight negative correlation in the anti-centre that is statistically significant. Together with the solar neighbourhood data, RPCC = − 0.64, − 0.82 and − 0.66 for neon, sulphur, and argon respectively.

Table A.17 compares the radial metallicity gradients from a selection of studies over multiple wavebands, including PNe, H ii regions, young B-type stars and Cepheid variables. Our analysis includes studies of neon, sulphur, and argon gradients. We also considered oxygen; even though its abundance changes over the course of stellar evolution in PNe, Pottasch & Bernard-Salas (2006) suggested that the oxygen gradient of sources that had not undergone hot bottom burning was identical to those of the other three elements from their PN sample. Our metallicity gradient slopes compare well with those from most other studies, though ours have greater uncertainties which arise from the dispersion at greater Rg, primarily due to the uncertainties in the distance measurements. Despite this, the slopes calculated from the PN studies of Maciel & Quireza (1999), Pottasch & Bernard-Salas (2006) and this work are typically steeper than those of other source types. Our metallicity gradients are consistent with a continuation at high Rg, though the slopes suggest an eventual flattening or steepening with distance is possible, particularly for neon and argon. While some of the studies we considered rule out flattening as a possibility, the studies of Andrievsky et al. (2002a,b,c) and Luck et al. (2003), analysing Cepheid variables, find discontinuities in the abundance gradient with Rg. They show that the gradient is seen to be steeper in regions closer to the Galactic centre (Rg ~ 4–6 kpc) and towards the anti-centre (Rg ~ 10–15 kpc) for 25 different elements, including oxygen and sulphur. Several studies that analyse sources beyond ~15 kpc (e.g. Rolleston et al. 2000; Lemasle et al. 2013; Fernández-Martín et al. 2017; this work) show steeper gradients than those found in most other studies. While this appears to agree with the findings of Andrievsky et al. (2002a,b,c), there are relatively few sources at these distances in each named study, hence the effect of the large distance sources on the respective gradients is likely to be minimal. By factoring in the large uncertainties in distance measurements for most sources in these samples, we do not find that the abundance gradients steepen with Rg in the direction of the anti-centre.

Analysis of the time evolution of the Galactic abundance gradient from this PN sample is also difficult, due to the large uncertainties in abundances and Galactocentric distances. Our data are consistent with a continuation of the gradient at large distances, so there is no suggestion that the inner and outer disks of the Milky Way evolved separately (e.g. Stanghellini & Haywood 2010; Kubryk et al. 2015).

thumbnail Fig. 3

Plots comparing the neon, sulphur, and argon abundances. The Galactic disk sample was analysed by Pottasch & Bernard-Salas (2006). The average abundances for the LMC and SMC are also shown, with those of neon and sulphur obtained from Bernard-Salas et al. (2008) and argon from Leisy & Dennefeld (2006).

4.3 α-process elements

In the evolution of low- to intermediate-mass stars, the abundances of elements heavier than carbon are generally not affected, except for those formed during the slow neutron-capture process (known as the s-process) which can occur during the AGB phase (e.g. Lugaro et al. 2012). As a result, the abundances of neon, sulphur, and argon should trace each other. Figure 3 shows that there is good agreement between these abundances, hence proving that they do trace each other well, though it is clear that the plot of sulphur against argon showsa greater dispersion. This could be explained by the need to account for two ions for argon (Ar+ and Ar3+), for which Ar3+ can be the dominantion.

Included on these plots are the abundances of Ne, S and Ar from the Magellanic Clouds. The mean abundance values for neon and sulphur were taken from an IR Spitzer sample of PNe from Bernard-Salas et al. (2008), and those for argon were taken from an optical PN sample from Leisy & Dennefeld (2006). The anti-centre sample shows abundances scattered around the LMC metallicity for each of the three elements, with few reaching values below those of the SMC.

The sulphur anomaly is the term coined by Henry et al. (2004), used to describe the observed sulphur abundances in PNe being lower than those of H II regions at the same metallicity (see also e.g. Henry et al. 2012). It was originally suggested that this could be explained by the lack of measured emission lines of ionisation states of S3+ and above in optical spectra and the need to account for them, particularly as S3+ can be a key stage of ionisation for sulphur. The sulphur anomaly has been seen in multiple galaxies; García-Rojas et al. (2016) observe this anomaly from four H II regions with abundances greater than most of the thirteen PNe in their sample from NGC 6822. Shaw et al. (2012) find the anomaly in the Magellanic Clouds from a combination of IR, optical and UV data, and Shingles & Karakas (2013) find the anomaly in the Milky Way from the PN sample of Pottasch & Bernard-Salas (2010), also from spectra in the same wavebands, compared to the ISM trend of H II regions and blue compact galaxies from the optical sample of Milingo et al. (2010).

In Fig. 4, we compare the sulphur abundances of PNe in the Galactic anti-centre and solar neighbourhood (Pottasch & Bernard-Salas 2006) with the sulphur abundances in two samples of Galactic H II regions, one derived from IR data (Martín-Hernández et al. 2002) and the other from optical data (Fernández-Martín et al. 2017). Both of these samples cover a similar range of Galactocentric distance values to the anti-centre PNe.

The IR H II region data from Martín-Hernández et al. (2002) agree well with the PN abundances, inferring that the sulphur anomaly is not observed from these data. However, their H II region abundances disagree with the interstellar and solar values of sulphur by a factor of approximately two to four. Based on this, they argue that their abundances are underestimated by up to a factor of four, which they ascribe to uncertainties in their ne and Te values, with the lack of S+ abundance values from their IR data from ISO accounting for a further ~15% discrepancy.We note that since the release of the paper of Martín-Hernández et al. (2002), the most widely used solar sulphur abundance value from the literature shows a decrease of ~20% from the value they used (Snow & Witt 1996; Asplund et al. 2009), though both the abundances of the anti-centre PNe and the H II regions remain low in comparison.

Comparing the infrared PN abundances with the optical H II region abundance data of Fernández-Martín et al. (2017), with Galactocentric distances of 11–17 kpc, shows a clear discrepancy between the two sets of data. The PN abundances are lower than those of the H II regions by a factor of approximately two, as shown by the two lines of best fit in Fig. 4. In this case, we clearly observe the sulphur anomaly.

thumbnail Fig. 4

Sulphur abundances of our PNe alongside samples of H II regions from Fernández-Martín et al. (2017) (FM+’17) and Martín-Hernández et al. (2002) (MH+’02), and solar neighbourhood PNe from Pottasch & Bernard-Salas (2006). The solid line is the line of best fit for PN abundances. The dashed line represents that for the H II regions of the FM+’17 sample.

5 Summary

We have presented an infrared spectroscopic study of 23 PNe in the Galactic anti-centre with Rg values of 8–21 kpc using Spitzer IRS to determine the abundances of neon, sulphur, and argon in a region that is a priori assumed to be metal-poor.

We have calculated the abundances in two ways: using empirically calculated ICFs from a combination of IR and optical data, and using more well-established ICFs from the literature. We find that the two methods produce similar results; the empirical ICFs consider a wider range of ionic states and are therefore small in value. We find that the abundances of neon, sulphur, and argon are lower in the anti-centre than those in the solar neighbourhood. The metallicity gradients of these elements seem to continue beyond Rg = 10 kpc despite a large spread of data values. The abundances of the α-process elements trace each other well, though there is a slightly larger dispersion between those of sulphur and argon.

Spitzer IRS has enabled the study of abundances from observations of PNe in the bulge, disk and halo of the Milky Way, as well as in nearby galaxies (primarily the Magellanic Clouds) at infrared wavelengths. With its greater sensitivity, the James Webb Space Telescope (JWST) will be able to continue to obtain spectra for PNe as far as the Local Group of galaxies, enabling us to carry out abundance studies over a wider range of parameter space. In addition, JWST will be able to spatially resolve PNe in the Milky Way, allowing us to investigate how the gas and dust content varies within these nebulae.

Acknowledgements

We thank the anonymous referee for their feedback and comments. This research has made use of the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This work has been supported by a Marie Curie FP7 CIG grant under project number 630861 (FEASTFUL).

Appendix A: Tables

Table A.1

23 PNe observed in the sample.

Table A.2

Aperture corrections applied to each of the line flux values.

Table A.3

Selected line intensity values from the Galactic anti-centre PNe, in units of 10−14 erg cm−2 s−1.

Table A.4

I() values for the sample, in erg cm−2 s−1.

Table A.5

Electron density values of PNe (cm−3).

Table A.6

Electron temperature values of PNe as averages of literature values.

Table A.7

Atomic data for ions shown in Table A.3.

Table A.8

Percentage contributions of ions that have required the use of empirically calculated ICFs towards their respective elemental abundances.

Table A.9

Ionic and total abundances of neon (× 10−5).

Table A.10

Ionic and total abundances of sulphur (× 10−6).

Table A.11

Ionic and total abundances of argon (× 10−7).

Table A.12

Ionic and elemental abundances from optical data in the literature, used in the calculations of neon, sulphur, and argon ICFs in Tables A.13A.15.

Table A.13

Comparison of the neon abundances and ICFs used in this study with those in which ICFs from other sources have been applied (see Sect. 3.5.1).

Table A.14

Comparison of the sulphur abundances and ICFs used in this study with those in which ICFs from other sources have been applied (see Sect. 3.5.2).

Table A.15

Comparison of the argon abundances and ICFs used in this study with those in which ICFs from other sources have been applied (see Sect. 3.5.3).

Table A.16

Comparison of heliocentric distances from Giammanco et al. (2011) and Frew et al. (2016).

Table A.17

Comparison of several abundance gradient studies which use various sources.

References

  1. Acker, A., Raytchev, B., Stenholm, B., & Tylenda, R. 1991, A&AS, 90, 89 [NASA ADS] [Google Scholar]
  2. Aitken, D. K., & Roche, P. F. 1982, MNRAS, 200, 217 [NASA ADS] [CrossRef] [Google Scholar]
  3. Aksaker, N., Yerli, S. K., K"i"z"i"loğlu, Ü., & Atalay, B. 2015, PASA, 32, 3 [NASA ADS] [CrossRef] [Google Scholar]
  4. Aller, L. H., & Czyzak, S. J. 1983, ApJS, 51, 211 [NASA ADS] [CrossRef] [Google Scholar]
  5. Aller, L. H., & Keyes, C. D. 1987, ApJS, 65, 405 [NASA ADS] [CrossRef] [Google Scholar]
  6. Aller, L. H., Keyes, C. D., & Feibelman, W. A. 1986, Proc. Natl. Acad. Sci., 83, 2777 [NASA ADS] [CrossRef] [Google Scholar]
  7. Anders, F., Chiappini, C., Minchev, I., et al. 2017, A&A, 600, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Andrievsky, S. M., Bersier, D., Kovtyukh, V. V., et al. 2002a, A&A, 384, 140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Andrievsky, S. M., Kovtyukh, V. V., Luck, R. E., et al. 2002b, A&A, 381, 32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Andrievsky, S. M., Kovtyukh, V. V., Luck, R. E., et al. 2002c, A&A, 392, 491 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Cosmic Abundances as Records of Stellar Evolution and Nucleosynthesis, eds. T. G. Barnes III, & F. N. Bash, ASP Conf. Ser. (San Fransisco: Astronomical Society of the Pacific), 336, 25 [Google Scholar]
  12. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  13. Barker, T. 1980, ApJ, 240, 99 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bernard-Salas, J. 2003, PhD Thesis, University of Groningen, The Netherlands [Google Scholar]
  15. Bernard-Salas, J., Pottasch, S. R., Gutenkunst, S., Morris, P. W., & Houck, J. R. 2008, ApJ, 672, 274 [NASA ADS] [CrossRef] [Google Scholar]
  16. Böhm-Vitense, E. 1958, Z. Astrophys., 46, 108 [Google Scholar]
  17. Butler, K., & Zeippen, C. J. 1994, A&AS, 108, 1 [NASA ADS] [Google Scholar]
  18. Carrasco, L., Serrano, A., & Costero, R. 1983, Rev. Mex. Astron. Astrofis., 8, 187 [NASA ADS] [Google Scholar]
  19. Chiappini, C., Matteucci, F., & Gratton, R. 1997, ApJ, 477, 765 [NASA ADS] [CrossRef] [Google Scholar]
  20. Chiappini, C., Matteucci, F., & Romano, D. 2001, ApJ, 554, 1044 [NASA ADS] [CrossRef] [Google Scholar]
  21. Condon, J. J., & Kaplan, D. L. 1998, ApJS, 117, 361 [NASA ADS] [CrossRef] [Google Scholar]
  22. Costa, R. D. D., Chiappini, C., Maciel, W. J., & de Freitas Pacheco, J. A. 1996, A&AS, 116, 249 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Costa, R. D. D., Uchida, M. M. M., & Maciel, W. J. 2004, A&A, 423, 199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Cuisinier, F., Acker, A., & Koeppen, J. 1996, A&A, 307, 215 [NASA ADS] [Google Scholar]
  25. de Freitas Pacheco, J. A., Maciel, W. J., Costa, R. D. D., & Barbuy, B. 1991, A&A, 250, 159 [NASA ADS] [Google Scholar]
  26. De Marco, O., & Crowther, P. A. 1999, MNRAS, 306, 931 [NASA ADS] [CrossRef] [Google Scholar]
  27. Delgado-Inglada, G., Morisset, C., & Stasińska, G. 2014, MNRAS, 440, 536 [NASA ADS] [CrossRef] [Google Scholar]
  28. Delgado-Inglada, G., Rodríguez, M., Peimbert, M., Stasińska, G., & Morisset, C. 2015, MNRAS, 449, 1797 [NASA ADS] [CrossRef] [Google Scholar]
  29. Dinerstein, H. L. 1980, ApJ, 237, 486 [NASA ADS] [CrossRef] [Google Scholar]
  30. Esteban, C., García-Rojas, J., & Pérez-Mesa, V. 2015, MNRAS, 452, 1553 [NASA ADS] [CrossRef] [Google Scholar]
  31. Esteban, C., Fang, X., García-Rojas, J., & Toribio San Cipriano, L. 2017, MNRAS, 471, 987 [NASA ADS] [CrossRef] [Google Scholar]
  32. Fernández-Martín, A., Pérez-Montero, E., Vílchez, J. M., & Mampaso, A. 2017, A&A, 597, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Fitzsimmons, A., Dufton, P. L., & Rolleston, W. R. J. 1992, MNRAS, 259, 489 [NASA ADS] [CrossRef] [Google Scholar]
  34. Fluks, M. A., Plez, B., The, P. S., et al. 1994, A&AS, 105, 311 [NASA ADS] [Google Scholar]
  35. French, H. B. 1981, ApJ, 246, 434 [NASA ADS] [CrossRef] [Google Scholar]
  36. Frew, D. J., Bojičić, I. S., & Parker, Q. A. 2013, MNRAS, 431, 2 [NASA ADS] [CrossRef] [Google Scholar]
  37. Frew, D. J., Parker, Q. A., & Bojičić, I. S. 2016, MNRAS, 455, 1459 [NASA ADS] [CrossRef] [Google Scholar]
  38. Friel, E. D. 1995, ARA&A, 33, 381 [NASA ADS] [CrossRef] [Google Scholar]
  39. Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1995, A&AS, 111, 347 [NASA ADS] [Google Scholar]
  40. Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1997, A&AS, 123, 159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1998, A&AS, 131, 499 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. García-Hernández, D. A., Ventura, P., Delgado-Inglada, G., et al. 2016, MNRAS, 458, L118 [NASA ADS] [CrossRef] [Google Scholar]
  43. García-Rojas, J., Peña, M., Flores-Durán, S., & Hernández-Martínez, L. 2016, A&A, 586, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Garnett, D. R. 1989, ApJ, 345, 282 [NASA ADS] [CrossRef] [Google Scholar]
  45. Genovali, K., Lemasle, B., da Silva, R., et al. 2015, A&A, 580, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Giammanco, C., Sale, S. E., Corradi, R. L. M., et al. 2011, A&A, 525, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Goodrich, R. W., & Dahari, O. 1985, ApJ, 289, 342 [NASA ADS] [CrossRef] [Google Scholar]
  48. Griffin, D. C.,& Badnell, N. R. 2000, J. Phys. B At. Mol. Phys., 33, 4389 [NASA ADS] [CrossRef] [Google Scholar]
  49. Griffin, D. C., Mitnik, D. M., & Badnell, N. R. 2001, J. Phys. B At. Mol. Phys., 34, 4401 [NASA ADS] [CrossRef] [Google Scholar]
  50. Henry, R. B. C., Kwitter, K. B., & Balick, B. 2004, AJ, 127, 2284 [NASA ADS] [CrossRef] [Google Scholar]
  51. Henry, R. B. C., Kwitter, K. B., Jaskot, A. E., et al. 2010, ApJ, 724, 748 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  52. Henry, R. B. C., Speck, A., Karakas, A. I., Ferland, G. J., & Maguire, M. 2012, ApJ, 749, 61 [NASA ADS] [CrossRef] [Google Scholar]
  53. Herwig, F., Bloecker, T., Schoenberner, D., & El Eid, M. 1997, A&A, 324, L81 [NASA ADS] [Google Scholar]
  54. Higdon, S. J. U., Devost, D., Higdon, J. L., et al. 2004, PASP, 116, 975 [NASA ADS] [CrossRef] [Google Scholar]
  55. Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 [NASA ADS] [CrossRef] [Google Scholar]
  56. Hummer, D. G., & Storey, P. J. 1987, MNRAS, 224, 801 [NASA ADS] [CrossRef] [Google Scholar]
  57. Hummer, D. G., Berrington, K. A., Eissner, W., et al. 1993, A&A, 279, 298 [NASA ADS] [Google Scholar]
  58. Johnson, C. T., Kingston, A. E., & Dufton, P. L. 1986, MNRAS, 220, 155 [NASA ADS] [CrossRef] [Google Scholar]
  59. Karakas, A. I., & Lattanzio, J. C. 2014, PASA, 31, 30 [Google Scholar]
  60. Kingsburgh, R. L., & Barlow, M. J. 1994, MNRAS, 271, 257 [NASA ADS] [CrossRef] [Google Scholar]
  61. Koeppen, J., Acker, A., & Stenholm, B. 1991, A&A, 248, 197 [NASA ADS] [Google Scholar]
  62. Kubryk, M., Prantzos, N., & Athanassoula, E. 2015, A&A, 580, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Kwitter, K. B., & Henry, R. B. C. 2001, ApJ, 562, 804 [NASA ADS] [CrossRef] [Google Scholar]
  64. Lebouteiller, V., Bernard-Salas, J., Sloan, G. C., & Barry, D. J. 2010, PASP, 122, 231 [NASA ADS] [CrossRef] [Google Scholar]
  65. Leisy, P., & Dennefeld, M. 2006, A&A, 456, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Lemasle, B., François, P., Genovali, K., et al. 2013, A&A, 558, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Luck, R. E., & Lambert, D. L. 2011, AJ, 142, 136 [NASA ADS] [CrossRef] [Google Scholar]
  68. Luck, R. E., Gieren, W. P., Andrievsky, S. M., et al. 2003, A&A, 401, 939 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Lugaro, M., Karakas, A. I., Stancliffe, R. J., & Rijs, C. 2012, ApJ, 747, 2 [NASA ADS] [CrossRef] [Google Scholar]
  70. Maciel, W. J., & Quireza, C. 1999, A&A, 345, 629 [NASA ADS] [Google Scholar]
  71. Maciel, W. J., Costa, R. D. D., & Uchida, M. M. M. 2003, A&A, 397, 667 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Maciel, W. J., Costa, R. D. D., & Cavichia, O. 2015, Rev. Mex. Astron. Astrofis., 51, 165 [NASA ADS] [Google Scholar]
  73. Magrini, L., Coccato, L., Stanghellini, L., Casasola, V., & Galli, D. 2016, A&A, 588, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Magrini, L., Randich, S., Kordopatis, G., et al. 2017, A&A, 603, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Marigo, P. 2001, A&A, 370, 194 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Marigo, P., Bernard-Salas, J., Pottasch, S. R., Tielens, A. G. G. M., & Wesselius, P. R. 2003, A&A, 409, 619 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Martín-Hernández, N. L., Peeters, E., Morisset, C., et al. 2002, A&A, 381, 606 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Mendoza, C. 1983, in Planetary Nebulae, ed. D. R. Flower, IAU Symp., 103, 143 [NASA ADS] [CrossRef] [Google Scholar]
  79. Mendoza, C., & Zeippen, C. J. 1982, MNRAS, 199, 1025 [NASA ADS] [CrossRef] [Google Scholar]
  80. Mendoza, C., & Zeippen, C. J. 1983, MNRAS, 202, 981 [NASA ADS] [CrossRef] [Google Scholar]
  81. Milingo, J. B., Henry, R. B. C., & Kwitter, K. B. 2002, ApJS, 138, 285 [NASA ADS] [CrossRef] [Google Scholar]
  82. Milingo, J. B., Kwitter, K. B., Henry, R. B. C., & Souza, S. P. 2010, ApJ, 711, 619 [NASA ADS] [CrossRef] [Google Scholar]
  83. Minchev, I., Chiappini, C., & Martig, M. 2013, A&A, 558, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Minchev, I., Famaey, B., Quillen, A. C., et al. 2012, A&A, 548, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Minchev, I., Chiappini, C., & Martig, M. 2014, A&A, 572, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Mitnik, D. M., Griffin, D. C., & Badnell, N. R. 2001, J. Phys. B At. Mol. Phys., 34, 4455 [NASA ADS] [CrossRef] [Google Scholar]
  87. Natta, A., Preite-Martinez, A., & Panagia, N. 1980, ApJ, 242, 596 [NASA ADS] [CrossRef] [Google Scholar]
  88. Peimbert, M., & Costero, R. 1969, Boletin de los Observatorios Tonantzintla y Tacubaya, 5, 3 [Google Scholar]
  89. Peimbert, M., Luridiana, V., & Torres-Peimbert, S. 1995, Rev. Mex. Astron. Astrofis., 31, 147 [NASA ADS] [Google Scholar]
  90. Pelan, J., & Berrington, K. A. 1995, A&AS, 110, 209 [NASA ADS] [Google Scholar]
  91. Pereira, C.-B.,& Miranda, L.-F. 2007, A&A, 467, 1249 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Perinotto, M., & Corradi, R. L. M. 1998, A&A, 332, 721 [NASA ADS] [Google Scholar]
  93. Phillips,J. P., & Márquez-Lugo, R. A. 2011, Rev. Mex. Astron. Astrofis., 47, 83 [NASA ADS] [Google Scholar]
  94. Pignatari, M., Herwig, F., Hirschi, R., et al. 2016, ApJS, 225, 24 [Google Scholar]
  95. Pottasch, S. R., & Beintema, D. A. 1999, A&A, 347, 975 [NASA ADS] [Google Scholar]
  96. Pottasch, S. R., & Bernard-Salas, J. 2006, A&A, 457, 189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  97. Pottasch, S. R., & Bernard-Salas, J. 2010, A&A, 517, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Pottasch, S. R., & Bernard-Salas, J. 2015, A&A, 583, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Pottasch, S. R., & Surendiranath, R. 2005, A&A, 432, 139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Rolleston, W. R. J., Smartt, S. J., Dufton, P. L., & Ryans, R. S. I. 2000, A&A, 363, 537 [NASA ADS] [Google Scholar]
  101. Rubin, R. H., Simpson, J. P., Erickson, E. F., & Haas, M. R. 1988, ApJ, 327, 377 [NASA ADS] [CrossRef] [Google Scholar]
  102. Saraph, H. E., & Storey, P. J. 1999, A&AS, 134, 369 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Sellwood, J. A., & Binney, J. J. 2002, MNRAS, 336, 785 [NASA ADS] [CrossRef] [Google Scholar]
  104. Shaver, P. A., McGee, R. X., Newton, L. M., Danks, A. C., & Pottasch, S. R. 1983, MNRAS, 204, 53 [NASA ADS] [CrossRef] [Google Scholar]
  105. Shaw, R. A., Lee, T.-H., Stanghellini, L., et al. 2012, IAU Symp., 283, 502 [NASA ADS] [Google Scholar]
  106. Shingles, L. J., & Karakas, A. I. 2013, MNRAS, 431, 2861 [NASA ADS] [CrossRef] [Google Scholar]
  107. Smartt, S. J. 2000, in The Evolution of the Milky Way: stars versus clusters, eds. F. Matteucci, & F. Giovannelli, Astrophys. Space Sci. Lib. (Dordrecht: Kluwer Academic Publishers), 255, 323 [NASA ADS] [CrossRef] [Google Scholar]
  108. Snow, T. P., & Witt, A. N. 1996, ApJ, 468, L65 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  109. Stanghellini, L., & Haywood, M. 2010, ApJ, 714, 1096 [NASA ADS] [CrossRef] [Google Scholar]
  110. Stasińska, G. 1978, A&A, 66, 257 [Google Scholar]
  111. Sterling, N. C., & Dinerstein, H. L. 2008, ApJS, 174, 158 [NASA ADS] [CrossRef] [Google Scholar]
  112. Sturm, E., Bauer, O. H., Brauer, J., et al. 1998, in Astronomical Data Analysis Software and Systems VII, eds. R. Albrecht, R. N. Hook, & H. A. Bushouse, ASP Conf. Ser., 154, 161 [NASA ADS] [Google Scholar]
  113. Surendiranath, R., Pottasch, S. R., & García-Lario, P. 2004, A&A, 421, 1051 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  114. Tamura, S., & Shaw, R. A. 1987, PASP, 99, 1264 [NASA ADS] [CrossRef] [Google Scholar]
  115. Thuan, T. X., Izotov, Y. I., & Lipovetsky, V. A. 1995, ApJ, 445, 108 [NASA ADS] [CrossRef] [Google Scholar]
  116. Tsamis, Y. G., Flores-Fajardo, N., Henney, W. J., Walsh, J. R., & Mesa-Delgado, A. 2013, MNRAS, 430, 3406 [NASA ADS] [CrossRef] [Google Scholar]
  117. Vera-Ciro, C., D’Onghia, E., Navarro, J., & Abadi, M. 2014, ApJ, 794, 173 [NASA ADS] [CrossRef] [Google Scholar]
  118. Vermeij, R., & van der Hulst, J. M. 2002, A&A, 391, 1081 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  119. Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 [NASA ADS] [CrossRef] [Google Scholar]
  120. Zhang, H. L., Graziani, M., & Pradhan, A. K. 1994, A&A, 283, 319 [NASA ADS] [Google Scholar]

1

Available from the SSC website: http://ssc.spitzer.caltech.edu.

All Tables

Table A.1

23 PNe observed in the sample.

Table A.2

Aperture corrections applied to each of the line flux values.

Table A.3

Selected line intensity values from the Galactic anti-centre PNe, in units of 10−14 erg cm−2 s−1.

Table A.4

I() values for the sample, in erg cm−2 s−1.

Table A.5

Electron density values of PNe (cm−3).

Table A.6

Electron temperature values of PNe as averages of literature values.

Table A.7

Atomic data for ions shown in Table A.3.

Table A.8

Percentage contributions of ions that have required the use of empirically calculated ICFs towards their respective elemental abundances.

Table A.9

Ionic and total abundances of neon (× 10−5).

Table A.10

Ionic and total abundances of sulphur (× 10−6).

Table A.11

Ionic and total abundances of argon (× 10−7).

Table A.12

Ionic and elemental abundances from optical data in the literature, used in the calculations of neon, sulphur, and argon ICFs in Tables A.13A.15.

Table A.13

Comparison of the neon abundances and ICFs used in this study with those in which ICFs from other sources have been applied (see Sect. 3.5.1).

Table A.14

Comparison of the sulphur abundances and ICFs used in this study with those in which ICFs from other sources have been applied (see Sect. 3.5.2).

Table A.15

Comparison of the argon abundances and ICFs used in this study with those in which ICFs from other sources have been applied (see Sect. 3.5.3).

Table A.16

Comparison of heliocentric distances from Giammanco et al. (2011) and Frew et al. (2016).

Table A.17

Comparison of several abundance gradient studies which use various sources.

All Figures

thumbnail Fig. 1

Full Spitzer IRS spectrum of M1-16 (top) and M1-8 (bottom). The low resolution spectra (SL) are shown on the left, and thehigh resolution spectra (SH and LH) are on the right.

In the text
thumbnail Fig. 2

Abundance gradients of neon, sulphur, and argon in the Milky Way. The dashed lines represent the oxygen abundance gradient from within the Galactic disk with a slope of −0.085 dex/kpc (Pottasch & Bernard-Salas 2006), passing through the solar value at 8.0 kpc (Asplund et al. 2005). The solid lines represent the line of best fit in each plot, with gradients of −0.058 ± 0.021, −0.079 ± 0.012 and −0.062 ± 0.023 dex/kpc respectively.

In the text
thumbnail Fig. 3

Plots comparing the neon, sulphur, and argon abundances. The Galactic disk sample was analysed by Pottasch & Bernard-Salas (2006). The average abundances for the LMC and SMC are also shown, with those of neon and sulphur obtained from Bernard-Salas et al. (2008) and argon from Leisy & Dennefeld (2006).

In the text
thumbnail Fig. 4

Sulphur abundances of our PNe alongside samples of H II regions from Fernández-Martín et al. (2017) (FM+’17) and Martín-Hernández et al. (2002) (MH+’02), and solar neighbourhood PNe from Pottasch & Bernard-Salas (2006). The solid line is the line of best fit for PN abundances. The dashed line represents that for the H II regions of the FM+’17 sample.

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.