Your new experience awaits. Try the new design now and help us make it even better

MINI REVIEW article

Front. Astron. Space Sci., 20 November 2025

Sec. Extragalactic Astronomy

Volume 12 - 2025 | https://doi.org/10.3389/fspas.2025.1673373

This article is part of the Research TopicConnecting Planetary Nebulae with Their Host Galaxy Stellar Population PropertiesView all 3 articles

Planetary nebulae as tracers of accreted stellar populations in massive galaxies in groups and clusters

  • 1Finnish Centre for Astronomy with ESO (FINCA), University of Turku, Turku, Finland
  • 2Tuorla Observatory, Department of Physics and Astronomy, University of Turku, Turku, Finland
  • 3Turku Collegium for Science, Medicine and Technology (TCSMT), University of Turku, Turku, Finland

Planetary nebulae (PNe) are valuable spatial and kinematic tracers of nearby galaxies. In this mini-review, I focus on their role in tracing the halo and intra-cluster/intra-group light assembly in groups and clusters of galaxies within 100 Mpc and, in particular, the link between characteristic PN metrics such as the α-parameter and the PN luminosity function and changes from the underlying in situ to ex-situ stellar populations. These results from nearby groups and clusters are placed into context with integral-field spectroscopic surveys of galaxies, which allow the co-spatial measurement of PN and stellar population properties. I provide an outlook on upcoming instrumentation that will provide new opportunities for the study of extragalactic PN populations. I address the challenges of reconciling observations of extragalactic PN populations with predictions from stellar evolution models and how revised late-stellar-evolution models have alleviated some of the tensions between observations and theory.

1 Introduction

Planetary nebulae (PNe) occur in the very late evolutionary stages of low-mass (initial masses between one and 8M) stars. The PN phase is considered part of the post-asymptotic giant branch (post-AGB) phase. The spectra of PNe are characterised by the absence of a stellar continuum and strong nebular emission lines. It is predicted that up to 15% of the radiation of their central stars is re-emitted in a single line: the [O III] λ5007A line (Dopita et al., 1992). This line can be used for the identification and radial velocity measurement of individual PNe in galaxies at distances as far as 100 Mpc (Gerhard et al., 2005), making them excellent tracers of the spatial distribution and kinematics of galaxies and their haloes.

Extragalactic PN populations are typically characterised with two metrics: the PN luminosity function (PNLF) and the luminosity-specific PN number, α-parameter for short. The PNLF describes the m5007 magnitude distribution, where m5007 is the magnitude corresponding to the flux of the [O III] λ5007A line (Jacoby, 1989). Empirically, the PNLF is near-universal in early- and late-type galaxies, described by the following relation (Ciardullo et al., 1989):

NPNm50071e3Mm5007+μec2m5007μ,(1)

where M denotes the absolute magnitude of the bright cut-off, μ the distance modulus, and c2 the faint-end slope. To first order, the absolute bright cut-off magnitude has been found to be invariant at M=4.5 (Ciardullo, 2022), facilitating the use of the PNLF as a secondary distance indicator by determining the apparent bright cut-off magnitude m and thus the distance modulus μ. While well-established observationally, the invariant bright cut-off cannot be easily reconciled with predictions from stellar evolution models, in particular for early-type galaxies with old stellar populations that were predicted to have much fainter bright cut-off magnitudes than observed (e.g., Marigo et al., 2004; Ciardullo, 2022) until the advent of revised post-AGB stellar evolution models (Miller Bertolami, 2016; Gesicki et al., 2018). The second term of Equation 1 was originally introduced with the motivation of describing the fading brightness of a gas cloud around the PN central star (Henize and Westerlund, 1963), with the slope fixed to c2=0.307 based on observations of M31 and other Local Group galaxies (Ciardullo et al., 1989).

By integrating the PNLF up to a limiting magnitude mlim, one obtains the total number of PNe, NPN, which can be related to the bolometric luminosity Lbol of the underlying stellar population via the α-parameter1 (Jacoby, 1980):

α=NPNLbol=BτPN,(2)

where B is the specific evolutionary flux2 and τPN the visibility lifetime of PNe. Accounting for incompleteness effects, and assuming that the bolometric luminosity Lbol of the host galaxy can be derived from its broadband photometry with a bolometric correction (Buzzoni et al., 2006), the α-parameter can be straightforwardly obtained from PN surveys. For a chemically homogeneous, coeval stellar population in the single stellar population (SSP) framework, the α-parameter can also be related to the visibility lifetime τPN (Buzzoni, 1989; Buzzoni et al., 2006). Assuming that the specific evolutionary flux B does not significantly vary between SSPs with different properties such as metallicity or initial mass function (Renzini and Buzzoni, 1986), the α-parameter can thus be interpreted as a direct proxy for the PN visibility lifetime.

Typical values of α are calculated within a limiting magnitude mlim=m+2.5 and are of the order of 108 to 109L,bol1. If the PN survey is shallower, the α-parameter can be extrapolated:

α2.5=ρmlimαmlimm.(3)

The extrapolation factor ρmlim can be obtained via integration of the PNLF (Equation 1):

ρmlim=mm+2.5Nm5007dm5007mmlimNm5007dm5007.(4)

In this case, the α-parameter can only be used to provide an estimate of the PN visibility lifetime within the limiting magnitude of the survey.

In the past, the primary target of PN surveys for galaxy kinematics have been early-type galaxies (e.g.,Arnaboldi et al., 1994; Hui et al., 1995; Coccato et al., 2009; Cortesi et al., 2013; Pulsoni et al., 2018), for the simple reason that their identification is facilitated by the lack of other strong [O III] emitters (such as supernova remnants and H II regions) in old stellar populations compared to late-type galaxies with younger stellar populations. Historically, surveys for PNe in galaxies made use of the on-off band technique, where the on band image is taken with a narrow-band [O III] filter, and the off band through an adjacent broad-band filter, for example, the V- or g-bands (e.g., Ford et al., 1973; Ford and Jenner, 1975). The detected PNe candidates are then followed up with spectrographs with high multiplexing capabilities to measure their line-of-sight velocity and identify remaining PN mimics such as Lyman-α emitters at high redshift. Counter-dispersed imaging, a variant of slitless spectroscopy, was developed to remove the need for a two-step process, allowing the detection and velocity measurement in one observation with the custom-built Planetary Nebula Spectrograph (PN.S, Douglas and Taylor, 1999; Douglas et al., 2002). However, due to the lack of full spectral information, this technique was predominantly used in early-type galaxies to limit the contamination from other [O III] emitters until a Hα-arm was installed at the PN.S (Aniyan et al., 2018; 2021).

These classical methods excel at detecting PNe in large areas, spanning several hundreds of square-arcseconds, corresponding to multiple effective radii of galaxies in the local Universe, allowing to map the kinematic transition from in situ to ex-situ stellar haloes (e.g., Pulsoni et al., 2018) as well as from haloes to the surrounding intra-cluster and intra-group light (ICL and IGL, respectively) as recently reviewed by Arnaboldi and Gerhard 2022, and references therein. In this review, I want to address how these kinematic transitions coincide with observed changes in the PN population properties, how these PN population properties can be linked to the characteristics of the underlying stellar populations, and how both new observing techniques and numerical stellar evolution models are necessary to contextualise these results.

This review is organised as follows: in Section 2, I review recent observational results concerning the link between PN- and stellar populations based on PN-surveys carried out with classical methods and new insights from integral-field spectroscopic surveys, and provide an outlook on new facilities and opportunities. In Section 3, I address new developments in modelling PN populations in the context of updated post-AGB stellar evolution models and numerical models of galaxy evolution. I conclude this review in Section 4.

2 Characterisation of PN populations from observations

2.1 Classical PN surveys

Surveys for PNe with classical techniques, such as ‘on-off’ imaging and slitless spectroscopy, provided first constraints on the link between stellar-population properties and the α- and PNLF-parameters. Ciardullo et al. (2002) suggested an oxygen-abundance-dependent PNLF bright cut-off M based on a small sample of a dozen galaxies with known Cepheid distances. Buzzoni et al. (2006) investigated the variation of the α-parameter with the underlying stellar population in the context of simple stellar population models, and based on those predicted an increase of the α-parameter with galaxy colour, in agreement with the observations available at the time.

In the past two decades, PNe have increasingly been used as kinematic tracers in early-type galaxies as well as to trace the ICL and IGL. Here, I want to especially focus on the transition from “in situ” to accreted (“ex-situ”) stellar populations, which is not only observationally signalled by radially varying light profiles and kinematics but also by changes in the PN population properties. In particular, several studies have found that the α-parameter increases significantly from the inner to the outer halos of the group- or cluster-dominating galaxy to the surrounding IGL or ICL (e.g., Doherty et al., 2009; Longobardi et al., 2013; 2018; Hartke et al., 2017; 2020). Furthermore, motivated by the observed steeper faint-end slope of the PNLF of Virgo ICL-PNe, Longobardi et al. (2013) proposed the modification of Equation 1 to allow a variable faint-end slope c2, which has since been applied to observations of PNe in nearby groups and clusters (Hartke et al., 2017; 2020) as well as in the nearby M31 (Bhattacharya et al., 2021). While it is difficult to constrain the stellar population properties of the low-surface-brightness haloes, and the surrounding ICL or IGL, pencil-beam surveys suggest that high α-parameters and steep PNLF slopes can be linked with old and metal-poor stellar populations (Williams et al., 2007; Harris et al., 2007; Lee and Jang, 2016), but with the caveat that these measurements are not exactly co-spatial. To interpret the studies of PN populations in groups and clusters of galaxies, it is important to obtain co-spatial measurements of PN and stellar populations, which will be discussed in turn.

2.2 New insights from integral-field spectroscopic surveys

Since the pioneering integral-field spectroscopic studies of Roth et al. (2004) and Sarzi et al. (2011), the Multi Unit Spectroscopic Explorer (MUSE, Bacon et al., 2010) has been transformational (see also the review by Roth et al., 2023) for the study of extragalactic PN populations, especially in late-type galaxies that have been targeted less with classical methods. Furthermore, integral-field spectroscopic observations allow the simultaneous investigation of stellar population properties and discovery and characterisation of PN populations, which are important for putting the results presented in the previous subsection into context.

Kreckel et al. (2017) showed the importance of full spectral information when determining the PNLF distance to NGC 628, as supernova remnants biased previous distance determinations based on narrowband imaging alone (Herrmann et al., 2008) to a shorter distance, despite the MUSE data only covering parts of the disk of the galaxy. With a much higher filling factor, the PHANGS-MUSE survey (Physics at High Angular resolution in Nearby GalaxieS; Emsellem et al., 2022) has been a treasure trove for studying ionized nebulae in nearby star-forming galaxies, with Scheuermann et al. (2022) determining PNLF distances to the 19 galaxies in the sample and can reconcile their derived distances with literature tip of the red-giant branch distances without the need for a metallicity-dependent M. To facilitate the classification of ten thousands of nebulae in the 19 galaxies into H II regions, supernova remnants, and PNe, Congiuet al. (2023) developed a Bayesian algorithm, but noted that the PNe sample is incomplete compared to that of Scheuermann et al. (2022), attributing the differences to different source-detection methods, but finding consistent classifications for the sources that are in both catalogues. Building on the PHANGS-MUSE strategy, Congiu et al. (2025) built the largest MUSE mosaic to date, covering the Sculptor galaxy (NGC 253), and resulting in the detection of 500 PNe. This work demonstrates the importance of accounting for dust in the interstellar medium, which can bias the PNLF distance to larger values.

Different techniques have been developed to detect PNe from integral-field spectroscopic data, starting with the ‘classical’ visual inspection and blinking of on-off images (e.g., Roth et al., 2018). Spriggs et al. (2020). developed a technique based on PSF and pixel-by-pixel spectral fitting to automatically detect PNe from Fornax3D survey data targeting the brightest early-type galaxies in the Fornax cluster (Sarzi et al., 2018). Galán-de Anta et al. (2021) applied this technique to three galaxies from Fornax3D to probe metallicity-dependent variations of the α-parameter in the centres and (inner) haloes of the galaxies, but did not find any evidence for an increase in the α-parameter between metal-rich and metal-poor regions. This differs from the results presented in Section 2.1 on halo and ICL α-parameter variations with stellar populations on significantly larger spatial scales that allowed for probing lower host stellar population metallicities.

Roth et al. (2021) developed the so-called differential emission line filtering (DELF) method with the goal of facilitating photometric measurements precise enough to use PNLF distances to alleviate the Hubble tension. Jacoby et al. (2024) applied this method to a heterogeneous sample of 20 galaxies with MUSE archival data, demonstrating that the method can yield excellent PNLFs and outline the way forward for PNLF distance measurements, both from an observational standpoint, as well as regarding the need for a better understanding and definition of the analytical form of the PNLF.

The aforementioned studies were all carried out with the MUSE integral-field spectrograph, that has a prohibitively small field-of-view (FoV, 1×1 in wide-field mode) to study galaxies in the local Universe in their full extent in a single pointing. The SITELLE instrument (Spectrométre Imageurä Transformeé deFourier pour l’É tude en Long et en Large des raies d’Émission; Grandmont et al., 2012), an optical imaging Fourier transform spectrometer, covers an 11×11 FoV, comparable to that of the PN.S instrument. However, the larger FoV is compromised by a shorter wavelength range that can be sampled in one observation. The ongoing Star formation, Ionized Gas, and Nebular Abundances Legacy Survey (SIGNALS; Rousseau-Nepton et al., 2019) provides observations of 40 late-type galaxies in the local Universe in three filters, SN1 (36503850 Å, covering [O II]λ3727 Å), SN2 (48005200 Å, covering the λλ4959,5007 Å doublet and Hβ) and SN3 (65106850 Å, covering Hα as well as the [Ne II]λλ6548,6583 Å and [S II]λλ6716,6731 Å lines), allowing for the detection and classification of extragalactic PNe in these galaxies, as piloted by Martin et al. (2018) and Vicens-Mouret et al. (2023).

2.3 New facilities and opportunities

The coming years and decades will see the arrival of several telescopes and instruments that may open new and exciting discovery spaces for PN populations outside of the Milky Way. BlueMUSE (Richard et al., 2019) at the Very Large Telescope, will have a similar design to MUSE that has been transformative for the study of PN populations in nearby galaxies, but covering important bluer emission lines (for example, [O II]λ3727 Å, Hγλ4340 Å, [O III]λ4364 Å, and He II λ4686 Å3). The bluer [O III] line is crucial for direct abundance determinations and the He II and Hγ lines provide important constraints for the determination of the excitation class of PNe (e.g., Reid and Parker, 2010; Bhattacharya et al., 2019).

Integral-field units will be complemented by spectroscopic facilities with high multiplexing capabilities and large fields of view, such as the Maunakea Spectroscopic Explorer (Sheinis et al., 2023) and the Wide-field Spectroscopic Telescope (Mainieri et al., 2024) at 10m-class facilities in the northern and southern hemispheres. These will be transformative for the study of individual nebulae, the relation between PNe and their host stellar populations, but also for PNe as tracers of galaxy, halo, and ICL assembly.

While the bulk of instrumentation for the upcoming Extremely Large Telescopes will focus on the infrared wavelength range, several spectrographs will still operate in the optical, covering the important [O III]λλ4959,5007 Å doublet from redshift z=0, facilitating unprecedentedly deep surveys for PNe in the local universe, as well as pushing the redshift boundaries of extragalactic PN populations beyond the Coma Cluster. The first light for the first generation of these instruments, the High Angular Resolution Monolithic Optical and Near-infrared Integral field spectrograph (HARMONI, Thatte et al., 2016) at ESO’s Extremely Large Telescope (ELT) and the Giant Magellan Telescope Multi-object Astronomical and Cosmological Spectrograph (GMACS, Fabricant et al., 2025) is foreseen for the coming decades. Instruments such as the Multi-AO Imaging Camera for Deep Observations (MICADO, Sturm et al., 2024) at the ELT may furthermore complement existing PN surveys with resolved stellar population information for galaxies beyond 10 Mpc.

3 Modelling PN populations and insights from stellar evolution models

The reconciliation of the observed invariance of the absolute magnitude of the PNLF bright cut-off for galaxies of different morphological types with the predictions from stellar evolution models has been a long-standing issue. Furthermore, the use of the PNLF as a distance indicator requires a better theoretical understanding about the PNe that populate the bright end of the PNLF and their origin. To reconcile observed PN properties with stellar evolution models, one needs to rely on synthetic post-AGB evolution tracks that describe the evolution of low-to intermediate-mass stars following the AGB phase. Historically, two model grids were widely used: Vassiliadis and Wood (1994) and Bloecker (1995). While it was possible to reproduce the observed bright cut-off of the PNLF M in galaxies with recent star-formation, the presence of similarly bright PNe in elliptical galaxies with old stellar populations could not be reconciled (e.g., Marigo et al., 2004; Ciardullo, 2012).

Furthermore, based on observations of PNe in the Galactic bulge, Gesicki et al. (2014) found that the Bloecker (1995) tracks4 evolve too slowly to reconcile the predicted with the observed local white dwarf masses. The post-AGB evolutionary tracks of Miller Bertolami (2016) address this issue, including updated descriptions of micro- and macrophysics, resulting in 310× faster post-AGB evolution timescales and overall brighter luminosities. Gesicki et al. (2018) were the first to show that with these faster evolving tracks, even populations with ages between three and 7 Gyrs could produce PNe that reach the bright cut-off M under the assumption that the brightest PNe are optically thick (maximum-nebula hypothesis). However, combining the new post-AGB tracks (Miller Bertolami, 2016) with the PN population modelling prescription of Mendez and Soffner (1997), Valenzuela et al. (2019) showed that under the maximum-nebula hypothesis too few PNe are produced overall, in conflict with observed PNLFs, and that accounting for optically-thin PNe is important.

The majority of the above models (and others in the literature) use recipes based on solar measurements and abundances, or, on larger scales, values derived based on Milky Way properties. The PICS (PNe in cosmological simulations) framework (Valenzuela et al., 2025) overcomes this limitation by modelling PNe for SSPs with different masses, ages, and metallicities. The PICS models furthermore explore the critical role of dust, using the empirical prescription of Jacoby and Ciardullo (2025), as well as the effect of different prescriptions for the initial-to-final mass relation (IFMR) and for the Helium abundances. While the authors reproduce the general trend from previous models that older SSPs produce less luminous PNe, they also demonstrate that metallicity plays an important role: old SSPs with higher metallicities are able to produce brighter PNe. Furthermore, they find the abundance of the bright PNe to be especially sensitive to the IFMR in old stellar populations, with a flatter IFMR (e.g., Cummings et al., 2018) leading to larger core masses and thus brighter PNe, alleviating some of the long-standing tension between observations and models, as also discussed in Jacoby and Ciardullo (2025).

By producing models normalised to the bolometric luminosity, Valenzuela et al. (2025) also determine the α-parameter for a given SSP in a given magnitude range. To facilitate the comparison with the predictions from Buzzoni et al. (2006), α-parameters are calculated eight magnitudes from a fixed M=4.5. While the general trends of α8 with age and metallicity are similar between the two models, the use of Miller Bertolami’s (2016) metallicity-dependent stellar evolution models by Valenzuela et al. (2025) leads to a stronger dependence of α8 on metallicity. To facilitate the comparison with observations, Valenzuela et al. (2025) also provide predictions for α2.5, reducing the need for observers with magnitude-limited data to make assumptions about the shape of the PNLF for extrapolating α, and better taking into account the effect of metallicity on the amount of bright PNe.

4 Conclusions and outlook

In this review, I discussed advances of linking extragalactic PN populations with the underlying stellar population properties both from an observational (Sect. 2) and modelling (Sect. 3) perspective. Especially on large spatial scales (i.e., comparing galaxy centres and extended haloes or the surrounding ICL in massive environments), there is evidence for a significant change of the α-parameter and PNLF shape. However, these measurements are limited by indirect metallicity inferences. To reconcile observations and theory, systematic surveys where PN and stellar population properties are simultaneously measured are needed, covering the entire mass-metallicity relation. As discussed in Section 2.2, integral-field spectrographs open up a new discovery space for PNe and allow for the co-spatial measurement of stellar population properties such as ages, metallicities, and abundances via spectral fitting.

Ongoing surveys, as well as those planned with new instrumentation at 8-, 10- and 30-m-class telescopes, using integral-field spectroscopy and ‘classical’ PN detection techniques, will provide constraints on the variation of the important diagnostics α-parameter and PNLF in stellar populations of in situ and ex-situ origins. For the interpretation of these observations, but also to further constrain theoretical models of late stellar evolution, it is crucial to properly model extragalactic PN populations in cosmological simulations of galaxy evolution, which has finally become possible, as reviewed in Section 3.

In summary, there are promising studies suggesting that, in the coming years, PNe may be elevated to stellar population tracers in low-surface brightness regions such as galaxy haloes and the ICL–where stellar population parameters such as age and metallicity cannot be easily measured directly–in addition to their important role as kinematic tracers. This is fuelled by advances in instrumentation as well as new models of late stellar evolution that will be combined with state-of-the-art numerical simulations of galaxy evolution.

Author contributions

Johanna Hartke: Conceptualization, Writing – original draft, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. JH acknowledges the support from the Turku Collegium for Science, Medicine and Technology (TCSMT) in the form of a starting grant and the financial support from the Visitor and Mobility program of the Finnish Centre for Astronomy with ESO (FINCA).

Acknowledgments

Acknowledgements

JH thanks her colleagues, in particular M. Arnaboldi, S. Bhattacharya, M. Bureau, L. Coccato, E. Congiu, A. Ennis, O. Gerhard, C. Pulsoni, M. Sarzi, S. Penger, C. Spiniello, L. Valenzuela, and N. Yang, for their contributions and discussions leading to this mini-review and also thanks the referee for their constructive comments that improved the manuscript.

Conflict of interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1Not to be confused with the α-abundance of elements.

2The rate of stars that evolve to the post-main sequence stages, normalised per unit light of the stellar population (Renzini and Buzzoni, 1986).

3This line can already be observed with MUSE in the extended spectral configuration, but the majority of extant archival data has been observed in the nominal configuration, starting at 4800 Å.

4The same holds for the tracks of Vassiliadis and Wood (1994), which evolve on even slower timescales.

References

Aniyan, S., Freeman, K. C., Arnaboldi, M., Gerhard, O. E., Coccato, L., Fabricius, M., et al. (2018). Resolving the disc-halo degeneracy - I: a look at NGC 628. Mon. Notices R. Astronomical Soc. 476, 1909–1930. doi:10.1093/mnras/sty310

CrossRef Full Text | Google Scholar

Aniyan, S., Ponomareva, A. A., Freeman, K. C., Arnaboldi, M., Gerhard, O. E., Coccato, L., et al. (2021). Resolving the disc-halo Degeneracy - II: NGC 6946. Mon. Notices R. Astronomical Soc. 500, 3579–3593. doi:10.1093/mnras/staa3106

CrossRef Full Text | Google Scholar

Arnaboldi, M., and Gerhard, O. (2022). Kinematics of the diffuse intragroup and intracluster light in groups and clusters of galaxies in the local universe within 100 Mpc distance. Front. Astronomy Space Sci. 9, 872283. doi:10.3389/fspas.2022.872283

CrossRef Full Text | Google Scholar

Arnaboldi, M., Freeman, K. C., Hui, X., Capaccioli, M., and Ford, H. (1994). The kinematics of the planetary nebulae in the outer regions of NGC 1399. Messenger 76, 40–44.

Google Scholar

Bacon, R., Accardo, M., Adjali, L., Anwand, H., Bauer, S., Biswas, I., et al. (2010). “The MUSE second-generation VLT instrument,”, 7735. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. doi:10.1117/12.856027

CrossRef Full Text | Google Scholar

Bhattacharya, S., Arnaboldi, M., Caldwell, N., Gerhard, O., Blaña, M., McConnachie, A., et al. (2019). The survey of planetary nebulae in andromeda (M 31). II. Age-velocity dispersion relation in the disc from planetary nebulae. Astronomy Astrophysics 631, A56. doi:10.1051/0004-6361/201935898

CrossRef Full Text | Google Scholar

Bhattacharya, S., Arnaboldi, M., Gerhard, O., McConnachie, A., Caldwell, N., Hartke, J., et al. (2021). The survey of planetary nebulae in Andromeda (M 31). III. Constraints from deep planetary nebula luminosity functions on the origin of the inner halo substructures in M 31. Astronomy Astrophysics 647, A130. doi:10.1051/0004-6361/202038366

CrossRef Full Text | Google Scholar

Bloecker, T. (1995). Stellar evolution of low- and intermediate-mass stars. II. Post-AGB evolution. Astronomy Astrophysics 299, 755–299.

Google Scholar

Buzzoni, A. (1989). Evolutionary population synthesis in stellar systems. I. A global approach. Astrophysical J. Suppl. Ser. 71, 817. doi:10.1086/191399

CrossRef Full Text | Google Scholar

Buzzoni, A., Arnaboldi, M., and Corradi, R. L. M. (2006). Planetary nebulae as tracers of galaxy stellar populations. Mon. Notices R. Astronomical Soc. 368, 877–894. doi:10.1111/j.1365-2966.2006.10163.x

CrossRef Full Text | Google Scholar

Ciardullo, R. (2012). The planetary nebula luminosity function at the dawn of gaia. Astrophysics Space Sci. 341, 151–161. doi:10.1007/s10509-012-1061-2

CrossRef Full Text | Google Scholar

Ciardullo, R. (2022). The planetary nebula luminosity function in the era of precision cosmology. Front. Astronomy Space Sci. 9, 896326. doi:10.3389/fspas.2022.896326

CrossRef Full Text | Google Scholar

Ciardullo, R., Jacoby, G. H., Ford, H. C., and Neill, J. D. (1989). Planetary Nebulae as standard candles. II. The calibration in M31 and its companions. Astrophysical J. 339, 53. doi:10.1086/167275

CrossRef Full Text | Google Scholar

Ciardullo, R., Feldmeier, J. J., Jacoby, G. H., de Naray, R. K., Laychak, M. B., and Durrell, P. R. (2002). Planetary Nebulae as standard candles. XII. Connecting the population I and population II distance scales. Astrophysical J. 577, 31–50. doi:10.1086/342180

CrossRef Full Text | Google Scholar

Coccato, L., Gerhard, O., Arnaboldi, M., Das, P., Douglas, N. G., Kuijken, K., et al. (2009). Kinematic properties of early-type galaxy haloes using planetary nebulae. Mon. Notices R. Astronomical Soc. 394, 1249–1283. doi:10.1111/j.1365-2966.2009.14417.x

CrossRef Full Text | Google Scholar

Congiu, E., Blanc, G. A., Belfiore, F., Santoro, F., Scheuermann, F., Kreckel, K., et al. (2023). PHANGS-MUSE: detection and Bayesian classification of ∼40 000 ionised nebulae in nearby spiral galaxies. Astronomy Astrophysics 672, A148. doi:10.1051/0004-6361/202245153

CrossRef Full Text | Google Scholar

Congiu, E., Scheuermann, F., Kreckel, K., Leroy, A., Emsellem, E., Belfiore, F., et al. (2025). The MUSE view of the Sculptor galaxy: survey overview and the luminosity function of planetary nebulae. Astronomy Astrophysics 700, A125. doi:10.1051/0004-6361/202554144

CrossRef Full Text | Google Scholar

Cortesi, A., Arnaboldi, M., Coccato, L., Merrifield, M. R., Gerhard, O., Bamford, S., et al. (2013). The planetary nebula Spectrograph survey of S0 galaxy kinematics. Data and overview. Astronomy Astrophysics 549, A115. doi:10.1051/0004-6361/201220306

CrossRef Full Text | Google Scholar

Cummings, J. D., Kalirai, J. S., Tremblay, P. E., Ramirez-Ruiz, E., and Choi, J. (2018). The White Dwarf initial-Final mass relation for Progenitor Stars from 0.85 to 7.5 M. Astrophysical J. 866, 21. doi:10.3847/1538-4357/aadfd6

CrossRef Full Text | Google Scholar

Doherty, M., Arnaboldi, M., Das, P., Gerhard, O., Aguerri, J. A. L., Ciardullo, R., et al. (2009). The edge of the M 87 halo and the kinematics of the diffuse light in the Virgo cluster core. Astronomy Astrophysics 502, 771–786. doi:10.1051/0004-6361/200811532

CrossRef Full Text | Google Scholar

Dopita, M. A., Jacoby, G. H., and Vassiliadis, E. (1992). A theoretical calibration of the planetary nebular cosmic distance Scale. Astrophysical J. 389, 27. doi:10.1086/171186

CrossRef Full Text | Google Scholar

Douglas, N. G., and Taylor, K. (1999). Galaxy kinematics from counter-dispersed imaging. Mon. Notices R. Astronomical Soc. 307, 190–196. doi:10.1046/j.1365-8711.1999.02614.x

CrossRef Full Text | Google Scholar

Douglas, N. G., Arnaboldi, M., Freeman, K. C., Kuijken, K., Merrifield, M. R., Romanowsky, A. J., et al. (2002). “The planetary nebula spectrograph: the green light for galaxy kinematics,”, 114. Publications of the Astronomical Society of the Pacific, 1234–1251. doi:10.1086/342765

CrossRef Full Text | Google Scholar

Emsellem, E., Schinnerer, E., Santoro, F., Belfiore, F., Pessa, I., McElroy, R., et al. (2022). The PHANGS-MUSE survey. Probing the chemo-dynamical evolution of disc galaxies. Astronomy Astrophysics 659, A191. doi:10.1051/0004-6361/202141727

CrossRef Full Text | Google Scholar

Fabricant, D., Catropa, D., Fata, R., Brown, W., Doherty, P., Durusky, D., et al. (2025). GMACS: a moderate-dispersion optical spectrograph for the giant Magellan Telescope. Publ. Astronomical Soc. Pac. 137, 035002. doi:10.1088/1538-3873/adb0f0

CrossRef Full Text | Google Scholar

Ford, H. C., and Jenner, D. C. (1975). Planetary nebulae in local group galaxies. II. Identifications, positions, number, and production rate of nebulae in NGC 221. Astrophysical J. 202, 365–371. doi:10.1086/153984

CrossRef Full Text | Google Scholar

Ford, H. C., Jenner, D. C., and Epps, H. W. (1973). Planetary Nebulae in Local-Group Galaxies. I. Identifications in NGC 185, NGC 205, and NGC 221. Astrophysical J. 183, L73. doi:10.1086/181255

CrossRef Full Text | Google Scholar

Galán-de Anta, P. M., Sarzi, M., Spriggs, T. W., Nedelchev, B., Pinna, F., Martín-Navarro, I., et al. (2021). The Fornax 3D project: PNe populations and stellar metallicity in edge-on galaxies. Astronomy Astrophysics 652, A109. doi:10.1051/0004-6361/202140834

CrossRef Full Text | Google Scholar

Gerhard, O., Arnaboldi, M., Freeman, K. C., Kashikawa, N., Okamura, S., and Yasuda, N. (2005). Detection of intracluster planetary nebulae in the coma cluster. Astrophysical J. 621, L93–L96. doi:10.1086/429221

CrossRef Full Text | Google Scholar

Gesicki, K., Zijlstra, A. A., Hajduk, M., and Szyszka, C. (2014). Accelerated post-AGB evolution, initial-final mass relations, and the star-formation history of the Galactic bulge. Astronomy Astrophysics 566, A48. doi:10.1051/0004-6361/201118391

CrossRef Full Text | Google Scholar

Gesicki, K., Zijlstra, A. A., and Bertolami, M. M. M. (2018). The mysterious age invariance of the planetary nebula luminosity function bright cut-off. Nat. Astron. 92, 1–5. doi:10.1038/s41550-018-0453-9

CrossRef Full Text | Google Scholar

Grandmont, F., Drissen, L., Mandar, J., Thibault, S., and Baril, M. (2012). “Final design of SITELLE: a wide-field imaging Fourier transform spectrometer for the Canada-France-Hawaii Telescope,”. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. doi:10.1117/12.926782

CrossRef Full Text | Google Scholar

Harris, W. E., Harris, G. L. H., Layden, A. C., and Wehner, E. M. H. (2007). The leo elliptical NGC 3379: a metal-poor halo emerges. Astrophysical J. 666, 903–918. doi:10.1086/520799

CrossRef Full Text | Google Scholar

Hartke, J., Arnaboldi, M., Longobardi, A., Gerhard, O., Freeman, K. C., Okamura, S., et al. (2017). The halo of M 49 and its environment as traced by planetary nebulae populations. Astronomy Astrophysics 603, A104. doi:10.1051/0004-6361/201730463

CrossRef Full Text | Google Scholar

Hartke, J., Arnaboldi, M., Gerhard, O., Coccato, L., Pulsoni, C., Freeman, K. C., et al. (2020). The halo of M 105 and its group environment as traced by planetary nebula populations. I. Wide-field photometric survey of planetary nebulae in the Leo I group. Astronomy Astrophysics 642, A46. doi:10.1051/0004-6361/202038009

CrossRef Full Text | Google Scholar

Henize, K. G., and Westerlund, B. E. (1963). Dimensions of diffuse and planetary nebulae in the Small Magellanic Cloud. Astrophysical J. 137, 747. doi:10.1086/147552

CrossRef Full Text | Google Scholar

Herrmann, K. A., Ciardullo, R., Feldmeier, J. J., and Vinciguerra, M. (2008). Planetary Nebulae in Face-On spiral galaxies. I. Planetary Nebula photometry and distances. Astrophysical J. 683, 630–643. doi:10.1086/589920

CrossRef Full Text | Google Scholar

Hui, X., Ford, H. C., Freeman, K. C., and Dopita, M. A. (1995). The planetary nebula System and dynamics of NGC 5128. III. Kinematics and Halo Mass distributions. Astrophysical J. 449, 592. doi:10.1086/176082

CrossRef Full Text | Google Scholar

Jacoby, G. H. (1980). The luminosity function for planetary nebulae and the number of planetary nebulae in local group galaxies. Astrophysical J. Suppl. Ser. 42, 1–18. doi:10.1086/190642

CrossRef Full Text | Google Scholar

Jacoby, G. H. (1989). Planetary Nebulae as standard candles. I. Evolutionary models. Astrophysical J. 339, 39. doi:10.1086/167274

CrossRef Full Text | Google Scholar

Jacoby, G. H., and Ciardullo, R. (2025). The critical role of dust on the [O III] planetary nebula luminosity function’s bright-end cutoff. Astrophysical J. 983, 129. doi:10.3847/1538-4357/adc0fb

CrossRef Full Text | Google Scholar

Jacoby, G. H., Ciardullo, R., Roth, M. M., Arnaboldi, M., and Weilbacher, P. M. (2024). Toward precision cosmology with improved planetary nebula luminosity function distances using VLT-MUSE. II. A Test sample from archival data. Astrophysical J. Suppl. Ser. 271, 40. doi:10.3847/1538-4365/ad2166

CrossRef Full Text | Google Scholar

Kreckel, K., Groves, B., Bigiel, F., Blanc, G. A., Kruijssen, J. M. D., Hughes, A., et al. (2017). A revised planetary nebula luminosity function distance to NGC 628 using MUSE. Astrophysical J. 834, 174. doi:10.3847/1538-4357/834/2/174

CrossRef Full Text | Google Scholar

Lee, M. G., and Jang, I. S. (2016). Dual stellar halos in the standard elliptical galaxy M105 and Formation of massive early-type galaxies. Astrophysical J. 822, 70. doi:10.3847/0004-637X/822/2/70

CrossRef Full Text | Google Scholar

Longobardi, A., Arnaboldi, M., Gerhard, O., Coccato, L., Okamura, S., and Freeman, K. C. (2013). The planetary nebula population in the halo of M 87. Astronomy Astrophysics 558, A42. doi:10.1051/0004-6361/201321652

CrossRef Full Text | Google Scholar

Longobardi, A., Arnaboldi, M., Gerhard, O., Pulsoni, C., and Söldner-Rembold, I. (2018). Kinematics of the outer halo of M 87 as mapped by planetary nebulae. Astronomy and Astrophysics 620, A111. doi:10.1051/0004-6361/201832729

CrossRef Full Text | Google Scholar

Mainieri, V., Anderson, R. I., Brinchmann, J., Cimatti, A., Ellis, R. S., Hill, V., et al. (2024). The Wide-field Spectroscopic Telescope (WST) Science white paper. doi:10.48550/arXiv.2403.05398

CrossRef Full Text | Google Scholar

Marigo, P., Girardi, L., Weiss, A., Groenewegen, M. A. T., and Chiosi, C. (2004). Evolution of planetary nebulae. II. Population effects on the bright cut-off of the PNLF. Astronomy Astrophysics 423, 995–1015. doi:10.1051/0004-6361:20040234

CrossRef Full Text | Google Scholar

Martin, T. B., Drissen, L., and Melchior, A.-L. (2018). A SITELLE view of M31’s central region - I. Calibrations and radial velocity catalogue of nearly 800 emission-line point-like sources. Mon. Notices R. Astronomical Soc. 473, 4130–4149. doi:10.1093/mnras/stx2513

CrossRef Full Text | Google Scholar

Mendez, R. H., and Soffner, T. (1997). Improved simulations of the planetary nebula luminosity function. Astronomy Astrophysics 321, 898–906. doi:10.48550/arXiv.astro-ph/9611128321

CrossRef Full Text | Google Scholar

Miller Bertolami, M. M. (2016). New models for the evolution of post-asymptotic giant branch stars and central stars of planetary nebulae. Astronomy Astrophysics 588, A25. doi:10.1051/0004-6361/201526577

CrossRef Full Text | Google Scholar

Pulsoni, C., Gerhard, O., Arnaboldi, M., Coccato, L., Longobardi, A., Napolitano, N. R., et al. (2018). The extended Planetary Nebula Spectrograph (ePN.S) early-type galaxy survey: the kinematic diversity of stellar halos and the relation between halo transition scale and stellar mass. Astronomy Astrophysics 618, A94. doi:10.1051/0004-6361/201732473

CrossRef Full Text | Google Scholar

Reid, W. A., and Parker, Q. A. (2010). “An evaluation of the excitation-class parameter for the central stars of planetary nebulae,”, Publ. Astron. Soc. Aust. 187–198. doi:10.1071/as09055

CrossRef Full Text | Google Scholar

Renzini, A., and Buzzoni, A. (1986). Global properties of stellar populations and the spectral evolution of galaxies (Spectral Evolution of Galaxies). Spectr. Evol. Galaxies 122, 195–231. doi:10.1007/978-94-009-4598-2_19

CrossRef Full Text | Google Scholar

Richard, J., Bacon, R., Blaizot, J., Boissier, S., and Boselli, A. (2019). BlueMUSE: project overview and science cases. arXiv e-prints , arXiv:1906.01657doi. doi:10.48550/arXiv.1906.01657

CrossRef Full Text | Google Scholar

Roth, M. M., Becker, T., Kelz, A., and Schmoll, J. (2004). Spectrophotometry of planetary nebulae in the Bulge of M31. Astrophysical J. 603, 531–547. doi:10.1086/381526

CrossRef Full Text | Google Scholar

Roth, M. M., Sandin, C., Kamann, S., Husser, T.-O., Weilbacher, P. M., Monreal-Ibero, A., et al. (2018). MUSE crowded field 3D spectroscopy in NGC 300. I. First results from central fields. Astronomy Astrophysics 618, A3. doi:10.1051/0004-6361/201833007

CrossRef Full Text | Google Scholar

Roth, M. M., Jacoby, G. H., Ciardullo, R., Davis, B. D., Chase, O., and Weilbacher, P. M. (2021). Toward precision cosmology with improved PNLF distances using VLT-MUSEI. Methodology and tests. Astrophysical J. 916, 21. doi:10.3847/1538-4357/ac02ca

CrossRef Full Text | Google Scholar

Roth, M. M., Jacoby, G., Ciardullo, R., Soemitro, A., Weilbacher, P. M., and Arnaboldi, M. (2023). Integral Field spectroscopy: a disruptive innovation for observations of planetary nebulae and the PNLF. arXiv e-prints 19, 50–62. doi:10.1017/s1743921323005768

CrossRef Full Text | Google Scholar

Rousseau-Nepton, L., Martin, R. P., Robert, C., Drissen, L., Amram, P., Prunet, S., et al. (2019). SIGNALS: I. Survey description. Mon. Notices R. Astronomical Soc. 489, 5530–5546. doi:10.1093/mnras/stz2455

CrossRef Full Text | Google Scholar

Sarzi, M., Mamon, G. A., Cappellari, M., Emsellem, E., Bacon, R., Davies, R. L., et al. (2011). The planetary nebulae population in the central regions of M32: the SAURON view. Mon. Notices R. Astronomical Soc. 415, 2832–2843. doi:10.1111/j.1365-2966.2011.18900.x

CrossRef Full Text | Google Scholar

Sarzi, M., Iodice, E., Coccato, L., Corsini, E. M., de Zeeuw, P. T., Falcón-Barroso, J., et al. (2018). Fornax3D project: overall goals, galaxy sample, MUSE data analysis, and initial results. Astronomy Astrophysics 616, A121. doi:10.1051/0004-6361/201833137

CrossRef Full Text | Google Scholar

Scheuermann, F., Kreckel, K., Anand, G. S., Blanc, G. A., Congiu, E., Santoro, F., et al. (2022). Planetary nebula luminosity function distances for 19 galaxies observed by PHANGS-MUSE. Mon. Notices R. Astronomical Soc. 511, 6087–6109. doi:10.1093/mnras/stac110

CrossRef Full Text | Google Scholar

Sheinis, A., Barden, S. C., and Sobeck, J. (2023). The maunakea spectroscopic explorer: thousands of fibers, infinite possibilities. Astron. Nachrichten 344, e20230108. doi:10.1002/asna.20230108

CrossRef Full Text | Google Scholar

Spriggs, T. W., Sarzi, M., Napiwotzki, R., Galán-de Anta, P. M., Viaene, S., Nedelchev, B., et al. (2020). Fornax 3D project: automated detection of planetary nebulae in the centres of early-type galaxies and first results. Astronomy Astrophysics 637, A62. doi:10.1051/0004-6361/201936862

CrossRef Full Text | Google Scholar

Sturm, E., Davies, R., Alves, J., Clénet, Y., Kotilainen, J., Monna, A., et al. (2024). “The MICADO first light imager for the ELT: overview and current status,”. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. doi:10.1117/12.3017752

CrossRef Full Text | Google Scholar

Thatte, N. A., Clarke, F., Bryson, I., Shnetler, H., Tecza, M., Fusco, T., et al. (2016). “The E-ELT first light spectrograph HARMONI: capabilities and modes,” in Ground-based and airborne instrumentation for astronomy VI. doi:10.1117/12.2230629

CrossRef Full Text | Google Scholar

Valenzuela, L. M., Méndez, R. H., and Miller Bertolami, M. M. (2019). Revised simulations of the planetary nebulae luminosity function. Astrophysical J. 887, 65. doi:10.3847/1538-4357/ab4e96

CrossRef Full Text | Google Scholar

Valenzuela, L. M., Miller Bertolami, M. M., Remus, R.-S., and Méndez, R. H. (2025). The PICS project: I. The impact of metallicity and helium abundance on the bright end of the planetary nebula luminosity function. Astronomy Astrophysics 699, A371. doi:10.1051/0004-6361/202553974

CrossRef Full Text | Google Scholar

Vassiliadis, E., and Wood, P. R. (1994). Post–Asymptotic giant branch evolution of Low-to intermediate-mass stars. Astrophysical J. Suppl. Ser. 92, 125. doi:10.1086/191962

CrossRef Full Text | Google Scholar

Vicens-Mouret, S., Drissen, L., Robert, C., Rousseau-Nepton, L., Martin, R. P., and Amram, P. (2023). Planetary nebulae and supernova remnants in NGC 4214 from the SIGNALS survey. Mon. Notices R. Astronomical Soc. 524, 3623–3635. doi:10.1093/mnras/stad2154

CrossRef Full Text | Google Scholar

Williams, B. F., Ciardullo, R., Durrell, P. R., Vinciguerra, M., Feldmeier, J. J., Jacoby, G. H., et al. (2007). The metallicity distribution of intracluster stars in virgo. Astrophysical J. 656, 756–769. doi:10.1086/510149

CrossRef Full Text | Google Scholar

Keywords: Planetary nebulae: general, galaxies: clusters, galaxies: groups, galaxies, elliptical and lenticular, cD, galaxies: halos

Citation: Hartke J (2025) Planetary nebulae as tracers of accreted stellar populations in massive galaxies in groups and clusters. Front. Astron. Space Sci. 12:1673373. doi: 10.3389/fspas.2025.1673373

Received: 25 July 2025; Accepted: 27 October 2025;
Published: 20 November 2025.

Edited by:

Marc Sarzi, Armagh Observatory, United Kingdom

Reviewed by:

Lucas Manuel Valenzuela, Universitäts-Sternwarte, Germany

Copyright © 2025 Hartke. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Johanna Hartke , am9oYW5uYS5oYXJ0a2VAdXR1LmZp

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.