Population II Post-AGB Stars. I. Methodology, Initial Results, and the First Luminosity Function

The gist The authors construct a multi-wavelength catalog for Post-asymptotic giant branch (PAGB) stars in 39 Galactic Globular Clusters to define their luminosity function and test their utility as reliable Population II distance indicators.

Methodology

The study constructs a multi-wavelength catalog spanning the ultraviolet (UV) to infrared (IR) for PAGB stars in 39 Galactic GCs In this work, we construct a multi-wavelength catalog, spanning the ultraviolet (UV) to infrared (IR), for PAGB stars in 39 Galactic GCs. They validate cluster membership with Gaia astrometry and identify 21 objects across 16 GCs using ground-based uBVI and space-based UV photometry. They perform Bayesian spectral energy distribution (SED) fitting with stellar atmosphere models to estimate their effective temperatures, luminosities, and radii.

Data Acquisition and Selection

The researchers selected 39 GCs spanning a broad range of metallicities, luminosities, and HB morphology. Almost all clusters were imaged from the UV through the IR using exclusively wide-field instruments. The photometric data sources included The Galaxy Evolution Explorer (GALEX) for far-UV and near-UV photometry, the Ultra-Violet/Optical Telescope (UVOT) for UV filters, and Gaia spacecraft data for cluster membership probabilities and G band magnitudes. They also included Two Micron All Sky Survey (2MASS) J band photometry to constrain the IR tail of the SEDs. Cluster membership was assessed using Gaia Early Data Release 3 (EDR3) parallaxes and proper motions (PMs) via a mixture modeling procedure, initially selecting stars with membership probabilities P > 70%.

PAGB Classification and Characterization

PAGB stars are classified based on their position in the optical MV versus (B − V)0 CMD using an AHB classification scheme from D22. They are categorized into red post-AGB (rPAGB) stars, yellow post-AGB (yPAGB) stars, blue post-AGB (bPAGB) objects, and violet post-AGB (vPAGB) stars. The classification boundaries are defined by tracing a locus of constant bolometric luminosity and using the curve from D22 to separate the AHB region below the PAGB boundary from the HB and AGB stars.

Stellar Parameter Derivation

To characterize these stars, researchers compare multi-wavelength photometry to a grid of ATLAS9 stellar atmosphere models. They first perform a traditional χ2 grid search to identify the best-fitting model and obtain initial parameter estimates. These results are then used to initialize a Bayesian analysis, which interpolates the model grids to refine parameter estimates and infer uncertainties from the posterior probability distributions. The resulting sample consists of 13 non-variable, intermediate-temperature and hot PAGB stars across 11 GCs.

Key Findings

The study found that PAGB stars occupy a narrow luminosity range centered at log(Lbol/L⊙) = 3.25 ± 0.03, with dispersion σ = 0.12 dex, supporting their use as reliable Population II distance indicators. The yPAGB stars have temperatures of Teff ≈ 6500 – 8250 K and luminosities of 1.5 – 1.9 × 103 L⊙, while the bPAGB stars are significantly hotter, with Teff ≈ 12,000 – 18,000 K and L ≈ 1.6 – 2.0×103 L⊙. The PAGB luminosity function has a mean bolometric luminosity of log(Lbol/L⊙) = 3.25 with an observed dispersion of σ = 0.12 dex. Furthermore, the majority of PAGB-hosting GCs have HBR > 0, consistent with previous expectations.

Limitations and Future Work

The analysis demonstrated that broadband photometry alone usually does not uniquely constrain these parameters, as inferred surface gravities and metal abundances are poorly constrained. The study found that the derived temperatures and luminosities from SED fitting are in good agreement with published spectroscopic analyses within the quoted uncertainties. Future work should apply this methodology to the full Galactic GC catalog and expand the sample size, particularly toward the hot end of the PAGB sequence, to improve constraints on stellar evolution models and test their utility as precise extragalactic standard candles. The authors also noted that high-resolution spectroscopic follow-up will be essential for obtaining better constraints on Teff, log g, and [Fe/H]. Additionally, the study suggests that incorporating spatially resolved reddening corrections or maps will be important for improving parameter estimates. The paper also explored the impact of dust extinction laws and filter red-leak effects, finding that they had a minimal impact on derived stellar parameters for PAGB stars.

The study establishes the framework and analysis required for a complete census of PAGB stars Using 39 Galactic GCs we have built a multi-wavelength catalog containing farUV through near-IR photometry assessed cluster membership probabilities using Gaia astrometry to minimize field star contamination defined quantitative PAGB classification boundaries and identified 21 candidate PAGB stars in 16 GCs. fit our broadband photometry to ATLAS9 stellar atmosphere models to infer the stars’ temperatures and luminosities evaluated degeneracies in estimated stellar metal abundances and surface gravities, demonstrating that broadband photometry alone usually does not uniquely constrain these parameters presented this sample’s PAGB luminosity function, which has a mean bolometric luminosity of log(Lbol/L⊙) = 3.25 with an observed dispersion of σ = 0.12 dex examined the correlation between PAGB frequency and HB morphology, finding that the majority of PAGB-hosting GCs have HBR > 0 — consistent with previous expectations compared our results with existing literature and showed that the temperatures and luminosities derived from broadband SED fitting are in good agreement with published spectroscopic analyses

The next step is to apply this methodology to the full Galactic GC catalog, evaluating previously identified PAGB stars and searching for additional candidates Expanding the sample size, particularly toward the hotend of the PAGB sequence, is crucial for improving constraints on stellar evolution models and testing the utility of PAGB stars for extragalactic distance measurements The current sample of 13 non-variable, intermediatetemperature and hot PAGB stars with UV–IR photometry is sufficient to constrain the mean bolometric luminosity of the class to log(Lbol/L⊙) = 3.25 (Mbol = −3.39), with an observed dispersion of 0.12 dex (0.30 mag) and a standard error of the mean of 0.03 dex (0.08 mag) A more complete survey will robustly constrain the intrinsic scatter of the luminosity function and test dependencies on metallicity or HB morphology Finally, a larger and more uniform sample will enable measurements of the frequency of PAGB stars across different cluster environments, the relative populations of vPAGB/bPAGB/yPAGB/rPAGB subclasses, and the evolutionary pathways and timescales of these stars. Ultimately, expanding the sample of PAGB stars is needed to establish their utility as precise extragalactic standard candles As demonstrated by the analysis of the yPAGB in NGC 1904 by T. S¸ahin & D. L. Lambert (2009) and the bPAGB star in NGC 6712 by M19, standard models like ATLAS9 may not fully capture the atmospheric structure and non-LTE effects of these high-luminosity, lowgravity stars High-resolution spectroscopic follow-up of PAGB candidates will provide the best constraints on Teff, log g, and [Fe/H], and enable direct measurements of line profiles, ionization balance, and elemental abundances As discussed in Section 4.3 and explored in Appendix A, assumptions in the reddening parameter RV can significantly impact SED fitting In this study, we adopt the canonical average Milky Way extinction curve (RV = 3.1) However, deviations from this frequently used standard value can introduce systematic errors in the derived stellar properties R. Zhang et al. (2023) built a 2D RV map using color excess ratios of three million stars from Large Sky Area Multi-Object Fiber Spectroscopy Telescope They found an average value of RV = 3.25 ± 0.

1. Bold header: Improved Stellar Parameter Derivation via Bayesian Inference

The improved AI system can derive effective temperatures, surface gravities, metal abundances, and stellar radius R by using a Bayesian framework that interpolates model grids to refine parameter estimates, as described in Section 4.2. This moves beyond traditional grid searches to provide a full posterior distribution for each star.

2. Bold header: Robust Luminosity Function Construction

The system can accurately construct the PAGB luminosity function by adopting the luminosities derived from the Bayesian analysis, which incorporates both measurement uncertainties and prior information, as stated in Section 5. The resulting function shows a mean bolometric luminosity of log(Lbol/L⊙) = 3.25 with an observed dispersion of σ = 0.12 dex.

3. Bold header: Enhanced Extinction Law Modeling

The system can perform more accurate SED fitting by incorporating spatially resolved reddening corrections, as suggested in Section 7.1, moving beyond the canonical average Milky Way extinction curve (RV = 3.1) to account for a varying extinction law across the Galaxy.

4. Bold header: Identification of Evolutionary Pathway Dependencies

The AI can explicitly model and test dependencies between stellar populations by examining the correlation between PAGB frequency and HB morphology, finding that "the majority of PAGB-hosting GCs have HBR > 0 — consistent with previous expectations."

Post-asymptotic giant branch (PAGB) stars share many advantages with the tip of the red giant branch (TRGB) for Population II distance measurements, including minimal internal reddening, easy identification, a solid theoretical foundation, and an expected weak metallicity dependence. At visual wavelengths PAGB stars are the brightest members in old populations, appearing up to about2 mag brighter than the TRGB and about4 mag brighter than RR Lyrae stars. Stellar evolution theory predicts that these objects occupy a narrow luminosity strip across the top of the Hertzsprung-Russell diagram, making them excellent standard candles. Observations in Galactic globular clusters (GCs) support this prediction, but a comprehensive survey is required to define their luminosity function and test for systematic errors. In this work, we construct a multi-wavelength catalog, spanning the ultraviolet (UV) to infrared (IR), for PAGB stars in 39 Galactic GCs. We validate cluster membership with Gaia astrometry, identify 21 objects across 16 GCs using ground-based uBVI and space-based UV photometry, and perform Bayesian spectral energy distribution (SED) fitting with stellar atmosphere models to estimate their effective temperatures, luminosities, and radii. We find that PAGB stars occupy a narrow luminosity range centered at (L bol/L)=3.25 plus or minus0.03, with dispersion σ=0.12 dex, supporting their use as reliable Population II distance indicators. Comparisons with theoretical tracks show excellent agreement with evolutionary models. We confirm that PAGB stars occur preferentially in metal-poor clusters with blue horizontal branches. Our results demonstrate that PAGB stars provide a powerful probe of late stellar evolution and a promising Population II standard candle.