Issue |
A&A
Volume 652, August 2021
|
|
---|---|---|
Article Number | A109 | |
Number of page(s) | 7 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202140834 | |
Published online | 19 August 2021 |
The Fornax 3D project: PNe populations and stellar metallicity in edge-on galaxies
1
Armagh Observatory and Planetarium, College Hill, Armagh BT61 9DG, UK
2
Astrophysics Research centre, School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 INN, UK
e-mail: pgalandeanta01@qub.ac.uk
3
Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield AL10 9AB, UK
4
Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany
5
Instituto de Astrofísica de Canarias, Vía Láctea s/n, 38205 La Laguna, Tenerife, Spain
6
Depto. Astrofísica, Universidad de La Laguna, Calle Astrofísico Francisco Sánchez s/n, 38206 La Laguna, Tenerife, Spain
7
European Southern Observatory, Karl-Schwarzschild-Straße 2, 85748 Garching bei München, Germany
8
Dipartimento di Fisica e Astronomia ‘G. Galilei’, Università di Padova, vicolo dell’Osservatorio 3, 35122 Padova, Italy
9
INAF–Osservatorio Astronomico di Padova, vicolo dell’Osservatorio 5, 35122 Padova, Italy
10
Sterrewacht Leiden, Leiden University, Postbus 9513, 2300 RA Leiden, The Netherlands
11
Max-Planck-Institut fuer extraterrestrische Physik, Giessenbachstrasse, 85741 Garching bei Muenchen, Germany
12
INAF – Astronomical Observatory of Capodimonte, Salita Moiariello 16, 80131 Naples, Italy
13
Department of Physics and Astronomy, Macquarie University, Sydney, NSW 2109, Australia
14
ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Sydney, Australia
15
Instituto de Astronomía y Ciencias Planetarias, Universidad de Atacama, Avenida Copayapu 485, Copiapó, Chile
16
Department of Astrophysics, University of Vienna, Türkenschanzstrasse 17, 1180 Vienna, Austria
17
Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281, 9000 Gent, Belgium
18
Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, PR China
Received:
19
March
2021
Accepted:
16
June
2021
Context. Extragalactic planetary nebulae (PNe) are useful distance indicators and are often used to trace the dark-matter content in external galaxies. At the same time, PNe can also be used as probes of their host galaxy stellar populations and to help understand the later stages of stellar evolution. Previous works have indicated that a specific number of PNe per stellar luminosity can vary across different galaxies and as a function of stellar-population properties, for instance increasing with decreasing stellar metallicity.
Aims. In this study we further explore the importance of stellar metallicity in driving the properties of the PNe population in early-type galaxies, using three edge-on galaxies in the Fornax cluster offering a clear view into their predominantly metal-rich and metal-poor regions near the equatorial plane or both below and above it, respectively.
Methods. Using very large telescope-multi unit spectroscopic explorer (VLT-MUSE) integral-field observations and dedicated PNe detection procedures, we constructed the PNe luminosity function and computed the luminosity-specific number of PNe α in both in- and off-plane regions of our edge-on systems.
Results. Comparing these α values with metallicity measurements also based on the same MUSE data, we find no evidence for an increase in the specific abundance of PNe when transitioning between metal-rich and metal-poor regions.
Conclusions. Our analysis highlights the importance of ensuring spatial consistency to avoid misleading results when investigating the link between PNe and their parent stellar populations, and suggest that in passively evolving systems variations in the specific number of PNe may pertain to rather extreme metallicity regimes found either in the innermost or outermost regions of galaxies.
Key words: planetary nebulae: general / galaxies: abundances / galaxies: elliptical and lenticular, cD / techniques: imaging spectroscopy / methods: observational
© ESO 2021
1. Introduction
Planetary nebulae (PNe) can be detected in external galaxies out to great distances thanks to their intense [O III] λ5007 nebular emission, emitting as much as 103 ∼ 104 L⊙ (e.g., O’dell 1963; Paczynski 1970; Rose & Smith 1970; Gesicki et al. 2018). This makes PNe useful tracers for the kinematics of their parent stellar populations in the very outskirts of galaxies allowing one to constrain their dark matter content (e.g., Romanowsky et al. 2003; Douglas et al. 2007; Coccato et al. 2009; Kafle et al. 2018). The PNe populations of external galaxies have also been found to have similar PN luminosity functions (PNLF) with a common bright-end cutoff that in turn makes PNe reliable cosmic distance indicators (Ciardullo et al. 1989; Ciardullo 2012).
As tracers of their parent stellar populations, extragalactic PNe can also help to understand the late stage of stellar evolution in galactic environments different from our Milky Way. In fact, both the detailed shape of the PNLF and the fact that the PNe population in different kinds of galaxies seems to share the same cutoff are still not fully understood. For instance, the actual presence of bright PNe in old, passively evolving systems is a long-standing problem (Marigo et al. 2004; Ciardullo 2005) that only recent theoretical development has started to address (e.g., Gesicki et al. 2018; Valenzuela et al. 2019). Likewise, both the origin of the PNLF bright cut-off (e.g., Ciardullo 2012) and its slope at fainter regimes are not fully understood (e.g., in star-forming regions, as found in the Large Magellanic Cloud by Reid & Parker 2010 or in the outskirts of early-type galaxies, as shown by Hartke et al. 2020).
At a more basic level, the abundance of PNe should also vary across different kinds of galaxies, as proposed by Buzzoni et al. (2006, hereafter B06), and in relation to variations to their overall far-UV flux (the so-called UV-upturn, Burstein et al. 1988). Furthermore, B06 showed the first piece of observational evidence for an anti-correlation between the specific number of PNe per luminosity (α) and stellar metallicity, finding at the same time that the presence of a UV-upturn also corresponds to less abundant PNe populations (consistent with previous theoretical work by Greggio & Renzini 1990 as well). It is important to note, however, that the analysis of B06 suffers from a spatial inconsistency in their comparison of PNe and stellar population measurements. Indeed, whereas their α values come from previous halo PNe studies (Ciardullo et al. 1989; Jacoby et al. 1989; Hui et al. 1993), their literature measurements for the Mgb index values (used to trace metallicity) and the UV-upturn pertain to the innermost regions of their sample galaxies.
Nowadays, deep imaging facilities allow us to trace radial variations in α out to the most metal-poor halo regions of galaxies (e.g., Bhattacharya et al. 2019; Hartke et al. 2020). Furthermore, the advent of integral-field spectroscopy can both facilitate the measurement of the stellar metallicity in faint galaxy outskirts (by collecting large aperture spectra, e.g., Weijmans et al. 2009) and detection of PNe deeply embedded in the central bright regions of galaxies (thanks to a detailed spectral modelling, e.g., Sarzi et al. 2011; Pastorello et al. 2013; Spriggs et al. 2020), where they are otherwise inaccessible to narrow-band imaging or counter-dispersed slit-less spectroscopy (e.g., Gerhard et al. 2005; Douglas et al. 2007; Ventimiglia et al. 2011).
In this respect, edge-on galaxies can be considered as ideal laboratories to test ideas on the connection between PNe and their parent stellar populations. Indeed, the particular inclination of these galaxies makes it possible to directly compare PNe populations in metal-rich disc regions and their counterparts in the more metal-poor bulge or halo, thus allowing one to start redressing the spatial inconsistencies in the B06 analysis.
Using MUSE integral-field spectroscopic data, we explore the PNe populations for the three edge-on galaxies FCC 153, FCC 170, and FCC 177 in this paper. These targets were observed during the magnitude-limited Fornax3D survey for bright galaxies within the virial radius of the Fornax cluster (Sarzi et al. 2018). We place our findings in the context of the stellar population measurements previously reported by Pinna et al. (2019a,b, hereafter P19). These edge-on galaxies host a range of stellar populations that can be divided into a predominantly metal-rich bulge/thin-disc and metal-poor off-plane populations. In this respect, PNe can be used to trace the kinematics of these parent stellar populations and check whether the abundance of hosted PNe is dependent on different star formation signatures such as the metallicity.
This paper is organised as follows. In Sect. 2, we present a brief summary of the sample galaxies as well as the properties of the observations. In Sect. 3, we explain the methodology and procedures to identify and confirm the detected PNe, as well as to estimate the PNLF, light-weighted metallicity, and luminosity specific PNe number. In Sect. 4, we present and discuss the results of both the PNLF for the three galaxies and the specific PNe number per galaxy as a function of the metallicity of different in-plane and off-plane regions. Lastly, in Sect. 5 we give our conclusions.
2. Observations and data reduction
The MUSE data for FCC 153, FCC 170, and FCC 177 (IC 1963, NGC 1381, and NGC 1380A, respectively) were obtained in the wide-field mode, which ensures a spatial sampling of 0.2″ × 0.2″ on a 1′×1′ field-of-view (FoV). The wavelength range of MUSE cubes ranges between 4650 Å and 9300 Å with a spectral sampling of 1.25 Å pixel−1 and an average instrumental spectral resolution of FWHMint = 2.5 Å (Sarzi et al. 2018). In the particular case of our three edge-on galaxies, the MUSE data are comprised of a central pointing and offset pointing further covering the outer disc and halo regions, and which from now on we refer to as the halo pointing. Different integration times of 1 h and 1h30min for a central and halo pointing, respectively, allow one to reach the same limiting surface brightness of μB = 25 mag arcsec−2. Central pointings were required to be observed under good seeing conditions (FWHM < 0.8″), whereas halo pointings had less stringent image quality constraints (FWHM < 1.5″). This is accounted in our PNe analysis (see Sect. 3.2).
The data reduction was performed using the MUSE pipeline (Weilbacher et al. 2012) applying the ESOREFLEX environment (Freudling et al. 2013) as described in Sarzi et al. (2018) and Iodice et al. (2019). This includes key steps such as sky-subtraction, telluric correction, and both relative and absolute flux calibration. Normally, the single pointings would be aligned through reference stars and further combined to produce final MUSE mosaics, but here this last step was skipped since our PNe analysis will need to be separately applied to our individual, finally reduced central and halo pointings due the aforementioned differences between their imaging quality.
Maps for the stellar kinematics of these objects are presented both in Pinna et al. (2019a,b) and in Iodice et al. (2019), with the former works additionally showing maps for the stellar age and metallicity. These studies also took notice of the absence of diffuse ionised-gas emission in these galaxies, although the presence of PNe in these and other passively evolving objects had already been mentioned in Iodice et al. (2019). This is thanks to a careful spaxel-by-spaxel separation of the stellar continuum and nebular emission using the GandALF code (e.g., Sarzi et al. 2006) within the novel GIST pipeline of Bittner et al. (2019), as also described in Sarzi et al. (2018).
3. Methodology
In order to characterise the PNe populations of our edge-on galaxies, we followed the method of Spriggs et al. (2020). This consists of five separate steps: (i) the identification of PNe candidate using [O III] λ5007 signal-to-noise maps obtained from our dedicated spaxel-by-spaxel GandALF fits; (ii) a dedicated 3D-fitting of such PNe candidates for their kinematics and total [O III] λ5007 flux while imposing a fixed spatial profile for their [O III] λ5007 emission according to the (pre-determined) spatial Moffat point-spread function (PSF)1; (iii) isolating and removing PNe interlopers based on the comparison with the host-galaxy stellar kinematics and unresolved HII-regions or supernovae remnant PNe impostors using line diagnostics; (iv) the construction of the PNLF and of our PNe detection incompleteness function (Sect. 3.1); and (v) finally, the estimation of the total luminosity-specific number of PNe within a given magnitude limit (Sect. 3.2). This last step requires either adopting or deriving the galaxy distance using the PNLF, and it involves re-scaling the incompleteness-corrected model for the PNLF to match the observed number of PNe.
The first three steps of this analysis deliver a catalogue for the PNe contained within the FoV of our MUSE observations. Separate PNe catalogues for each pointing are needed in the first place since the central and halo data were obtained under rather different seeing conditions (Table 1), which has an impact on the shape and extent of the incompleteness functions that we then folded in the PNLF modelling. Our detected PNe are in agreement with those in Spriggs et al. (2021), after that, the only difference is that here we separately applied the last two steps of our PNe analysis to the predominantly metal-rich and metal-poor regions near and off the equatorial plane of our edge-on galaxies, respectively. We obtained separate PNLF and luminosity-specific PNe numbers in order to explore variations in the PNe populations of regions with significantly different parent stellar populations. For this, we applied the stellar population differentiation adopted in Pinna et al. (2019b,a), which is shown by horizontal lines in the top panels of Figs. 1–3. Finally, using two complementary approaches, we proceeded to estimate the stellar metallicity in these regions (Sect. 3.3) in order to place our PNe population results in the context of the B06 relation, which we also re-evaluated using the stellar metallicity measurements from the literature (Sect. 3.4).
![]() |
Fig. 1. PNLF of FCC 153 for in-plane and off-plane regions. Top panel: map of the [O III] λ5007 line. The blue circles represent the detected PNe in central pointing whereas the green ones show the detected objects in the halo pointing. The black dashed lines indicate the separation between in-plane and off-plane regions. Bottom panels: PNLF of both components: in-plane (left panel) and off-plane (right panel). The blue histograms as well as the blue solid and dashed lines represent the PNLF for PNe in the central pointing, whereas the green histograms along with the green solid and dashed lines indicate the PNLF for the PNe in the halo pointing. The red contours are the isophotes of the flux white image. |
![]() |
Fig. 2. Same as Fig. 1, but for FCC 170. The black circular area corresponds to the excluded template mismatch region. |
Point-spread function across the MUSE pointings of the sample galaxies, as characterised using a Moffat function fit to either a foreground star or a selection of bright PNe sources.
3.1. Observed PNLF, completeness corrections, and best-matched PNLF models
Following the procedure outlined in the previous section, we obtained separate catalogues for PNe populations encompassed by our central and halo MUSE pointings. The location of our detected PNe is shown in the top panels of Figs. 1–3 for each galaxy, with blue and green circles for PNe identified in the central and halo pointing, respectively. We then further divided the central and halo PNe into two sub-catalogues for regions close to the equatorial plane and those lying above and below it, and we constructed corresponding PNe luminosity functions as shown in the lower panels of Figs. 1–3. To first order, after accounting for the incompleteness of our observations, adjusting for distance and applying an appropriate re-scaling, our observed PNLF should compare well with the standard form of the PNLF function introduced by Ciardullo et al. (1989). This is given by
![]() |
Fig. 3. Same as Fig. 1, but for FCC 177. The black elliptical area corresponds to the excluded foreground galaxy. |
where is the characteristic bright-end cut-off (according to the latest calibration of Ciardullo 2012) that makes the PNLF a useful distance indicator.
Following Spriggs et al. (2020), we defined the PNe detection completeness at a given apparent magnitude m5007 as the fraction of galaxy stellar light within the MUSE FoV where a PNe of that particular magnitude can be detected. For this, we applied our PNe detection criteria on a spaxel-by-spaxel basis by first computing the peak [O III] λ5007 flux for PNe of apparent magnitude m5007 (given our PSF model) and subsequently checking if the corresponding peak A[O III] spectral amplitude (given the MUSE spectral resolution and the typical PNe instrinsic σ of ∼40 km s−1) exceeds three times the local residual-noise level from our previous spaxel-by-spaxel GandALF spectral fitting.
This procedure is equivalent to randomly populating the MUSE FoV with simulated PNe of magnitude m5007, while considering that larger numbers would be expected in brighter regions, and then simply evaluating the completeness as the fraction of PNe that would be detected (see also Sarzi et al. 2011; Pastorello et al. 2013). During this procedure, we also excluded regions affected by either background or foreground sources (e.g., in FCC 177) and where we have evidence of significant systematic effects in our spectral fitting (i.e., template-mismatch during the spaxel-by-spaxel GandALF fit) that would lead to spurious PNe detection or biased [OIII] flux measurements (e.g., in the central regions of FCC 170). The outlined procedure was applied to the in and off-plane regions separately in both central and halo pointings to construct the PNLFs shown in Figs. 1–3.
With these completeness-correction functions at hand, we first considered the best distance to our sample objects from Spriggs et al. (in prep.), using all the data from the central pointing, where the seeing conditions are less prominent. In short, this procedure adjusts the distance modulus until the observed cumulative PNLF distribution is optimally matched by the corresponding distance-shifted and completeness-corrected cumulative Ciardullo PNLF, through a minimisation of the Kolgorov-Smirnov statistics. Once this best distance was obtained, we computed appropriately shifted and completeness-corrected model PNLF for each of the in- and off-plane regions in the central and halo pointing, further rescaling them until their integrated value matched the total number of PNe observed in these regions. Such adjusted model PNLFs are also shown in the lower panels of Figs. 1–3, and they are based on best-fitting distances of 21.7, 19.6, and 17.8 Mpc for FCC153, FCC170, and FCC177, respectively (Spriggs et al., in prep.). As detailed in that paper, these distances are in good agreement with other measurements, such as those based on surface brightness fluctuations as in Blakeslee et al. (2009).
3.2. Total and luminosity-specific PNe numbers
Once the Ciardullo et al. (1989) PNLF had been matched to the observed PNLF in the in- and off-plane regions of our edge-on galaxies, we could estimate the total number of PNe from the bright cut-off down to a chosen magnitude limit by simply integrating the shifted and re-scaled model PNLF as it is. Following Buzzoni et al. (2006), in this way we estimated the true number N2.5 of PNe 2.5 magnitude from the PNLF bright cut-off. To this extrapolated total, we associated errors obtained by appropriately rescaling the Poisson uncertainty on the actual number of PNe that we observed within a given region, computed following the prescription of Gehrels (1986).
From this number, we can estimate the so-called luminosity specific number of PNe
where Lbol is the bolometric luminosity for the region where we have been looking for PNe. In turn, this is given by
where Mλ and Mλ, ⊙ are the absolute magnitude for the region of interest within the galaxy and the Sun within a chosen broad-band filter, while BCλ and BCλ, ⊙ are the corresponding bolometric magnitude corrections from the given broad-band filter. See also Torres (2010, Eq. (9)) and Hartke et al. (2020, Eq. (14)).
In the case of our MUSE observations, we chose the SDSS g-band to compute the observed mg and absolute Mg magnitude of the in- and off-plane regions in our edge-on galaxies. As for deriving the bolometric correction, we followed the procedure described by Spriggs et al. (2020), whereby BCg was derived using the optimal combination of EMILES stellar templates (Vazdekis et al. 2016) that best matches the MUSE integrated spectrum observed within the region of interest (i.e., our in- and off-plane region in both a central and halo pointing). Adopting a bolometric correction for the Sun in the g-band of BCg, ⊙ = −1.78 mag and a g-band absolute magnitude of Mg, ⊙ = 5.23 mag (Willmer 2018), as well as using Eq. (3), we estimated the bolometric luminosity of our in- and of-plane regions, which are reported in Table 2. Errors in Lbol are dominated by Mλ errors from distance uncertainties.
Results of the analysis for the MUSE central and halo pointings of the sample galaxies.
3.3. Stellar metallicity measurements from integrated MUSE spectra
In order to place our PNe findings in relation to the properties of their parent stellar population and in particular to check on any dependence of α2.5 with stellar metallicity, we proceeded to measure the latter in both the in- and off-plane regions in and above the disc of our edge-on galaxies, both in MUSE central and halo pointings. For this, we integrated our MUSE spectra in these regions (as defined by P19; Sect. 3) over a minimum signal-to-noise of 1, obtaining four different spectra for each galaxy, and we derived the stellar metallicity according to the following two approaches: (i) based on full spectral fitting using the pPXF code (Cappellari & Emsellem 2004; Cappellari 2017) in the 4750–5500 Å wavelength range as done in P19, using a set of single-age stellar population spectral models from the MILES library (Falcón-Barroso et al. 2011) to estimate an average value for the stellar metallicity; (ii) spectral fitting focused on the wavelength regions around a specific set of absorption-line features particularly sensitive to either age (e.g., Hβ), metallicity (e.g., Fe 5270 and Fe 5335), or α-element abundance (e.g., Mg 5177) as done in Martín-Navarro et al. (2019), using also the single-age stellar population models from the MILES library to derive a single, best-fitting value for the metallicity [Z/H]SSP.
These two methods are complementary. The spectral-fitting approach of P19 allows one to account for an extended star-formation and metal-enrichment history through a combination of single-age models of varying age and metallicity instead of adopting only the values corresponding to the best matching single-age stellar population model as in Martín-Navarro et al. (2019). The latter approach, however, also accounts for more subtle variations in α-element abundances (and in fact, also in the slope at the low-mass-end of the stellar initial-mass function), to which a full-spectral fitting can be rather insensitive 2.
3.4. Revising the Buzzoni et al. relation
The evidence for the anti-correlation between the specific number of PNe α2.5 and stellar metallicity of B06 was based on literature measurements for the strength of the Mg2 absorption-line index rather than an actual metallicity measurement. To quantify the α2.5 variations that they would have found with actual metallicity measurements, we took all the objects in the B06 sample that were also observed over the course of the ATLAS3D survey (Cappellari et al. 2011). Using the measurements of McDermid et al. (2015), we revisited the α2.5-metallicity relation of B06. As in the case of our two adopted approaches to estimate the stellar metallicity, McDermid et al. (2015) provide average (light-weighted) values based on full-spectral fitting [Z/H]SFH and single-age values from the best stellar-population model [Z/H]SSP that best matched the observed strength of a set of absorption line-strength indices, while accounting for [α/Fe] variations.
Figure 4 shows that such a revised α2.5 – metallicity relation remains in place only when considering [Z/H]SFH values, indicating at best a factor ∼6 variation in the specific number of PNe per 0.4 dex variation in metallicity. No visible trend is present when considering [Z/H]SSP values. This suggests that either [α/Fe] may also play a part in the original trend with the Mg2 line-strength index, which also subsists when looking at the ATLAS3D MgbValues, or caution is needed when considering the widely different spatial scales probed by B06.
![]() |
Fig. 4. PNe specific number vs stellar metallicity for Buzzoni et al. (2006) sample galaxies covered also by the Atlas3D integral-field spectroscopic survey (see McDermid et al. 2015), for two different star formation history measurements: full spectral fitting from a single-age stellar population (left panel) and full-index-fitting based on a focused set of absorption lines (right panel). |
4. Results and discussion
We give our final estimates for the specific number of PNe 2.5 magnitudes from the PNLF bright cut-off and stellar metallicity in the regions near or above and below the equatorial plane of our sample galaxies summarised in Table 2. Figure 5 finally shows how these quantities fare against each other across such different regions and as observed both in our central and halo MUSE pointings.
![]() |
Fig. 5. Luminosity specific PNe number α as a function of metallicity for the best-fitting single-age stellar population to selected spectral regions (top panel) following Martín-Navarro et al. (2019) and for the light-weighted metallicity resulting from a full spectral fitting (bottom panel) following P19. Circles and triangles show values from in-plane and off-plane regions, respectively whereas filled and open symbols denote whether these were found within our central or halo MUSE pointings. The colour of the symbols indicate measurements from each galaxy, with dashed lines showing linear fits to all values for each object. |
We observe no significant difference in α2.5 between in- and off-plane regions, irrespective of the methodology used to estimate the stellar metallicity, finding at most a factor ∼2.6 increase in α2.5 between the off-plane measurements of FCC 177 where there is a 0.33 dex difference. This lack of a trend appears to contradict the results of B06. Indeed, even though the metallicity estimates for the predominantly metal-poor off-plane regions can span a wide range of values, on average they are ∼0.4 dex lower than their counterparts closer to the equatorial plane and therefore ought to show up to ∼6 times higher values for α2.5 according to our most optimistic revision of the B06 analysis (Sect. 3.4). Despite the large uncertainties in the α2.5 values for such fainter off-plane regions, there is simply no room to accommodate a factor 6 (i.e., ∼0.8 dex) increase in log(α2.5).
This result is at odds with the recent work of Hartke et al. (2020), who find a clear increase in the specific PNe number in the outskirts of the closer Leo Group early-type galaxy NGC 3379 (M105), based on deep narrow-band on-off imaging extending out to nearly 23 effective radii (Re). This study reveals, however, that such an enhancement in α values begins from 8Re onwards, deep in the outer stellar halo or even in the intra-group light medium and at a metallicity regime most likely below the one found in the metal-poor extra-planar regions of our edge-on galaxies. On the other hand, at smaller radii, Hartke et al. (2020) find no α variation despite the presence of large stellar population gradients (e.g., within 3–4Re, Weijmans et al. 2009) and this is consistent with our findings. Changing perspective, the lack of a trend between α and metallicity also contrasts with the decrease in the number of bright post-AGB stars and the increase in the far-UV flux towards the nuclear metal-rich regions of M31 found by Rosenfield et al. (2012), and this corresponds to tentative evidence for a decrease in the number of bright PNe (Pastorello et al. 2013). Although edge-on systems are not particularly suited to explore the innermost and crowded PNe populations of galaxies, future MUSE studies assisted by adaptive optics have the potential to shed more light in this respect.
5. Conclusions
Using MUSE observations, we have explored the PNe populations in three edge-on galaxies in the Fornax cluster, offering an unique benchmark for testing the presence of systematic differences in the PNe content between predominantly metal-rich regions near the equatorial plane and generally metal-poor off-plane populations. We have found no significant evidence of a change in the specific number of PNe between the metal-poor and metal-rich regions of our edge-on galaxies, despite us probing a range in stellar metallicity values (0.5 dex) similar to that of previous investigations which were, however, comparing the properties of PNe and their parent stellar population on rather different spatial scales. Presently, our results further suggest that variations in the specific number of PNe in passively evolving systems may pertain to extreme metallicity regimes and to either the innermost or outermost regions of galaxies.
The PSF evaluation assumes a Moffat (1969) profile and is achieved using foreground stars by simultaneously applying our 3D-fitting procedure to the [O III] λ5007 emission of a few bright PNe. The radial extent and kurtosis of the Moffat profile is defined by its parameters α and β, which relate to the .
Acknowledgments
This work was supported by Science and Technology Facilities Council [grant number ST/R504786/1]. We really appreciate the computational facilities provided by University of Hertfordshire and Swinburne University to provide us with the necessary storage and analysis of the MUSE data, necessary to carry out our analysis. F.P., I. M-N, and J. F-B acknowledge support through the RAVET project by the grant PID2019-107427GB-C32 from the Spanish Ministry of Science, Innovation and Universities (MCIU), and through the IAC project TRACES which is partially supported through the state budget and the regional budget of the Consejería de Economía, Industria, Comercio y Conocimiento of the Canary Islands Autonomous Community. E.M.C. acknowledges support by Padua University grants DOR1885254/18, DOR1935272/19, and DOR2013080/20 and by MIUR grant PRIN 2017 20173ML3WW001.GvdV acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No 724857 (Consolidator Grant ArcheoDyn).
References
- Bhattacharya, S., Arnaboldi, M., Hartke, J., et al. 2019, A&A, 624, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bittner, A., Falcón-Barroso, J., Nedelchev, B., et al. 2019, A&A, 628, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Blakeslee, J. P., Jordán, A., Mei, S., et al. 2009, ApJ, 694, 556 [NASA ADS] [CrossRef] [Google Scholar]
- Burstein, D., Bertola, F., Buson, L. M., Faber, S. M., & Lauer, T. R. 1988, ApJ, 328, 440 [NASA ADS] [CrossRef] [Google Scholar]
- Buzzoni, A., Arnaboldi, M., & Corradi, R. L. M. 2006, MNRAS, 368, 877 [NASA ADS] [CrossRef] [Google Scholar]
- Cappellari, M. 2017, MNRAS, 466, 798 [NASA ADS] [CrossRef] [Google Scholar]
- Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 [NASA ADS] [CrossRef] [Google Scholar]
- Cappellari, M., Emsellem, E., Krajnović, D., et al. 2011, MNRAS, 416, 1680 [NASA ADS] [CrossRef] [Google Scholar]
- Ciardullo, R. 2005, in Planetary Nebulae as Astronomical Tools, eds. R. Szczerba, G. Stasińska, & S. K. Gorny, AIP Conf. Ser., 804, 277 [NASA ADS] [CrossRef] [Google Scholar]
- Ciardullo, R. 2012, Ap&SS, 341, 151 [NASA ADS] [CrossRef] [Google Scholar]
- Ciardullo, R., Jacoby, G. H., Ford, H. C., & Neill, J. D. 1989, ApJ, 339, 53 [NASA ADS] [CrossRef] [Google Scholar]
- Coccato, L., Gerhard, O., Arnaboldi, M., et al. 2009, MNRAS, 394, 1249 [NASA ADS] [CrossRef] [Google Scholar]
- Douglas, N. G., Napolitano, N., Romanowsky, A. J., et al. 2007, ApJ, 664, 257 [NASA ADS] [CrossRef] [Google Scholar]
- Falcón-Barroso, J., Sánchez-Blázquez, P., Vazdekis, A., et al. 2011, A&A, 532, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gehrels, N. 1986, ApJ, 303, 336 [NASA ADS] [CrossRef] [Google Scholar]
- Gerhard, O., Arnaboldi, M., Freeman, K. C., et al. 2005, ApJ, 621, L93 [NASA ADS] [CrossRef] [Google Scholar]
- Gesicki, K., Zijlstra, A. A., & Miller Bertolami, M. M. 2018, Nat. Astron., 2, 580 [NASA ADS] [CrossRef] [Google Scholar]
- Greggio, L., & Renzini, A. 1990, in Windows on Galaxies, eds. G. Fabbiano, J. S. Gallagher, & A. Renzini, Astrophys. Space Sci. Lib., 160, 63 [NASA ADS] [CrossRef] [Google Scholar]
- Hartke, J., Arnaboldi, M., Gerhard, O., et al. 2020, A&A, 642, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hui, X., Ford, H. C., Ciardullo, R., & Jacoby, G. H. 1993, ApJ, 414, 463 [NASA ADS] [CrossRef] [Google Scholar]
- Iodice, E., Sarzi, M., Bittner, A., et al. 2019, A&A, 627, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jacoby, G. H., Ciardullo, R., Ford, H. C., & Booth, J. 1989, ApJ, 344, 704 [NASA ADS] [CrossRef] [Google Scholar]
- Kafle, P. R., Sharma, S., Lewis, G. F., Robotham, A. S., & Driver, S. P. 2018, MNRAS, 475, 4043 [NASA ADS] [CrossRef] [Google Scholar]
- Marigo, P., Girardi, L., Weiss, A., Groenewegen, M. A. T., & Chiosi, C. 2004, A&A, 423, 995 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martín-Navarro, I., Lyubenova, M., van de Ven, G., et al. 2019, A&A, 626, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- McDermid, R. M., Alatalo, K., Blitz, L., et al. 2015, MNRAS, 448, 3484 [NASA ADS] [CrossRef] [Google Scholar]
- Moffat, A. 1969, A&A, 3, 455 [Google Scholar]
- O’dell, C. 1963, ApJ, 138, 67 [CrossRef] [Google Scholar]
- Paczynski, B. 1970, Acta Astron., 20, 47 [Google Scholar]
- Pastorello, N., Sarzi, M., Cappellari, M., et al. 2013, MNRAS, 430, 1219 [NASA ADS] [CrossRef] [Google Scholar]
- Pinna, F., Falcón-Barroso, J., Martig, M., et al. 2019a, A&A, 623, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pinna, F., Falcón-Barroso, J., Martig, M., et al. 2019b, A&A, 625, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Reid, W. A., & Parker, Q. A. 2010, MNRAS, 405, 1349 [NASA ADS] [Google Scholar]
- Romanowsky, A. J., Douglas, N. G., Kuijken, K., et al. 2003, in American Astronomical Society Meeting Abstracts, BAAS, 35, 1312 [Google Scholar]
- Rose, W. K., & Smith, R. L. 1970, ApJ, 159, 903 [NASA ADS] [CrossRef] [Google Scholar]
- Rosenfield, P., Johnson, L. C., Girardi, L., et al. 2012, ApJ, 755, 131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sarzi, M., Falcón-Barroso, J., Davies, R. L., et al. 2006, MNRAS, 366, 1151 [Google Scholar]
- Sarzi, M., Mamon, G. A., Cappellari, M., et al. 2011, MNRAS, 415, 2832 [NASA ADS] [CrossRef] [Google Scholar]
- Sarzi, M., Iodice, E., Coccato, L., et al. 2018, A&A, 616, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Spriggs, T. W., Sarzi, M., Napiwotzki, R., et al. 2020, A&A, 637, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Torres, G. 2010, AJ, 140, 1158 [NASA ADS] [CrossRef] [Google Scholar]
- Valenzuela, L. M., Méndez, R. H., & Miller Bertolami, M. M. 2019, ApJ, 887, 65 [NASA ADS] [CrossRef] [Google Scholar]
- Vazdekis, A., Koleva, M., Ricciardelli, E., Röck, B., & Falcón-Barroso, J. 2016, MNRAS, 463, 3409 [Google Scholar]
- Ventimiglia, G., Arnaboldi, M., & Gerhard, O. 2011, A&A, 528, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Weijmans, A.-M., Cappellari, M., Bacon, R., et al. 2009, MNRAS, 398, 561 [NASA ADS] [CrossRef] [Google Scholar]
- Weilbacher, P. M., Streicher, O., Urrutia, T., et al. 2012, in Software and Cyberinfrastructure for Astronomy II, eds. N. M. Radziwill, G. Chiozzi, SPIE Conf. Ser., 8451, 84510B [CrossRef] [Google Scholar]
- Willmer, C. N. 2018, ApJS, 236, 47 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Point-spread function across the MUSE pointings of the sample galaxies, as characterised using a Moffat function fit to either a foreground star or a selection of bright PNe sources.
Results of the analysis for the MUSE central and halo pointings of the sample galaxies.
All Figures
![]() |
Fig. 1. PNLF of FCC 153 for in-plane and off-plane regions. Top panel: map of the [O III] λ5007 line. The blue circles represent the detected PNe in central pointing whereas the green ones show the detected objects in the halo pointing. The black dashed lines indicate the separation between in-plane and off-plane regions. Bottom panels: PNLF of both components: in-plane (left panel) and off-plane (right panel). The blue histograms as well as the blue solid and dashed lines represent the PNLF for PNe in the central pointing, whereas the green histograms along with the green solid and dashed lines indicate the PNLF for the PNe in the halo pointing. The red contours are the isophotes of the flux white image. |
In the text |
![]() |
Fig. 2. Same as Fig. 1, but for FCC 170. The black circular area corresponds to the excluded template mismatch region. |
In the text |
![]() |
Fig. 3. Same as Fig. 1, but for FCC 177. The black elliptical area corresponds to the excluded foreground galaxy. |
In the text |
![]() |
Fig. 4. PNe specific number vs stellar metallicity for Buzzoni et al. (2006) sample galaxies covered also by the Atlas3D integral-field spectroscopic survey (see McDermid et al. 2015), for two different star formation history measurements: full spectral fitting from a single-age stellar population (left panel) and full-index-fitting based on a focused set of absorption lines (right panel). |
In the text |
![]() |
Fig. 5. Luminosity specific PNe number α as a function of metallicity for the best-fitting single-age stellar population to selected spectral regions (top panel) following Martín-Navarro et al. (2019) and for the light-weighted metallicity resulting from a full spectral fitting (bottom panel) following P19. Circles and triangles show values from in-plane and off-plane regions, respectively whereas filled and open symbols denote whether these were found within our central or halo MUSE pointings. The colour of the symbols indicate measurements from each galaxy, with dashed lines showing linear fits to all values for each object. |
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.