Issue |
A&A
Volume 693, January 2025
|
|
---|---|---|
Article Number | A238 | |
Number of page(s) | 11 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202449571 | |
Published online | 21 January 2025 |
Extragalactic globular cluster near-infrared spectroscopy
II. The stellar population synthesis model zero-point problem extends to the near-infrared
1
Instituto de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain
2
Departamento de Astrofísica, Universidad de La Laguna, E-38205 La Laguna, Tenerife, Spain
3
School of Astronomy, Institute for Research in Fundamental Sciences – IPM, 19395-5531 Tehran, Iran
4
Instituto de Física, Universidade Federal do Rio Grande do Sul (UFRGS), Av. Bento Gonçalves, 9500 Porto Alegre, RS, Brazil
5
Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, China
6
Centre for Astrophysics and Supercomputing, Swinburne University, John Street, Hawthorn, VIC 3122, Australia
7
Waterloo Centre for Astrophysics, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
8
Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
9
Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, 05508-900 São Paulo, Brazil
⋆ Corresponding author; eftekhari.eli@gmail.com
Received:
11
February
2024
Accepted:
1
December
2024
Many recent studies have pointed out significant discrepancies between observations and models of stellar populations in the near-infrared (NIR). With current and future observing facilities being focused in this wavelength range, properly assessing and solving these issues is of utmost importance. Here, we present the first application of the extragalactic globular cluster (GC) near-infrared spectroscopy survey, and present evidence that these GCs reveal an age zero-point problem of stellar population synthesis (SPS) models. This problem has already been identified in the optical range for the GCs of the Milky Way. Such an issue arises when derived GC spectroscopic ages appear older than the Universe itself. We extend this discussion for the first time to the NIR, specifically using the Paβ line at 1.28 microns. We focus on the GCs of the nearby Centaurus A galaxy using their NIR spectra. This work broadens our understanding of the age zero-point problem and emphasises the necessity to revisit and refine SPS models, especially in the NIR domain.
Key words: galaxies: star clusters: general / galaxies: stellar content / infrared: galaxies
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
The age zero-point problem presents a significant challenge in stellar population synthesis (SPS) models, notably in spectral indices analysis and the accurate age determination of old stellar populations such as globular clusters (GCs). This issue emerges as a discrepancy where the spectroscopically derived ages of these populations (> 10 Gyr; e.g. Strader et al. 2005; Chies-Santos et al. 2011) appear older than the estimated age of the Universe. This problem was initially identified in the optical spectra of Milky Way (MW) GCs. For example, the study of the optical Balmer line Hγ in the high-metallicity GC 47 Tucanae (47 Tuc) by Gibson et al. (1999) highlighted this issue.
Several hypotheses have been proposed to explain the age zero-point problem. One key area of discussion is the influence of horizontal branch (HB) stars and turn-off (TO) on Balmer lines. These lines are sensitive to the age of the stellar populations and are mainly contributed by the hottest stars, like those at the main sequence turn-off (MSTO) and the HB. In metal-rich regimes, old populations usually have an HB that is predominantly a red clump of cool stars (Teff < 5000 K) and therefore the hotter MSTO stars primarily contribute to the Balmer line. Vazdekis (2001) examined the impact of HB on the Balmer lines in 47 Tuc and found a minimal influence. However, in low-metallicity regimes (< − 1 dex), the HB contains stars with temperatures comparable to the MSTO, leading to an age degeneracy, where very old, low-metallicity populations appear younger from their integrated light. In such metallicity regimes, blue HB stars can complicate the age zero-point problem by predicting stronger Balmer lines for oldest simple stellar populations (SSPs). At intermediate metallicities (−1.0 ≤ [Z/H] ≤ −0.5), where the relative contributions of the HB stars depend on the morphology of this evolutionary branch, Mendel et al. (2007) found that while HB morphology does affect Hβ, its impact on derived ages for most Galactic GCs (GGCs) remains relatively minor, regardless of the HB morphology used. Additionally, atomic diffusion (the movement of elements throughout the layers of a star) alters the surface composition and temperature of stars, especially at the TO point, making them appear cooler and thus older. Incorporating atomic diffusion into modelling isochrones can cause spectroscopic age estimates to move closer to those derived from colour–magnitude diagrams (CMDs) (Vazdekis 2001). Moreover, α-enhanced isochrones affect the temperature and luminosity of stars in the TO point. According to Vazdekis (2001), implementing α-enhanced isochrones in SPS models leads to a closer alignment between spectroscopic age estimates and those derived from CMDs.
The nebular emission fill-in hypothesis also attempts to address the age zero-point problem; it proposes that gas emission might partially fill in Balmer lines, affecting their strengths and consequently the age estimates derived from them. However, the study of the behaviour of higher-order Balmer lines by Vazdekis (2001) reveals inconsistencies in the effects of emission fill-in across different Balmer lines, weakening the support for this hypothesis.
Another aspect to consider in explaining the age zero-point problem is the impact of the systematic uncertainties influencing fundamental stellar atmospheric parameters, such as effective temperature, surface gravity, and iron abundance, on the line-strength predictions of SPS models. The study by Percival & Salaris (2009) suggests that typical systematic offsets of 100 K in Teff, 0.15 dex in [Fe/H], and 0.25 dex in log g affect the Hβ index’s derived ages by up to 20% for both old and intermediate-age populations. The interpretation of ages becomes more challenging when considering the HδF and HγF indices, as their response to systematic shifts in [Fe/H] contrast with those of Hβ. However, the study by Vazdekis et al. (2010) indicates that overcoming the model zero-point problem requires a larger temperature offset of around 150 K, which exceeds current temperature scale uncertainties. These authors used line-strength diagnostic diagrams from their models and alternative models with different temperature scales to determine the ages and metallicities of GGC samples.
All these findings point to the fact that interpreting Balmer lines is complicated, and different factors can influence GC age estimates. As another example, a recent study by Leath et al. (2022) brings new insights into interpreting Balmer lines. Their work, which focuses on integrated spectroscopy of GGCs, points out a split in the Hβo line-strength index at intermediate to high metallicities ([Z/H] > −1), resulting in an apparent ‘upper branch’ and ‘lower branch’ in the Hβo–[MgFe] diagram1. Contrary to previous belief that this was caused by the presence of hot blue straggler stars, leading to an underestimation of spectroscopic ages (Cenarro et al. 2008), their findings suggest that this split is actually due to increased helium abundance.
Interestingly, the age zero-point problem is also observed when using extragalactic GCs. For example, Fig. 8 in Cenarro et al. (2007) shows that most metal-poor GCs in the elliptical galaxy NGC 1407 fall below model predictions for Balmer lines. Likewise, Leath et al. (2022)’s exploration of He abundance effects on the Hβo index in M31 GCs shows a similar pattern (see their Figure 11). Such studies expand our understanding of the age zero-point problem, offering insights into environmental impacts on stellar populations.
An additional aspect of the age zero-point problem that is worth considering is the extension of the analysis to the near-infrared (NIR) spectral range. While past studies have mainly focused on optical wavelengths, recent advancements have shifted attention to the NIR spectral range. This range is often invoked to address optical challenges such as the age–metallicity degeneracy and high sensitivity to dust. However, SPS models in this spectral range are still in their infancy. This is particularly evident in the failure of SPS models to reproduce some NIR absorptions in galaxies and the challenges in interpreting CO absorption bands in massive early-type galaxies (ETGs) (Mármol-Queraltó et al. 2009; Baldwin et al. 2018; Dahmer-Hahn et al. 2018; Riffel et al. 2019; Eftekhari et al. 2022). Aiming to investigate these issues further, we observed the first data set of NIR spectra of extragalactic GCs (Dahmer-Hahn et al. 2024, hereafter Paper I). This sample is composed of 21 GCs in Centaurus A/NGC 5128 (hereafter Cen A). Cen A is among the most massive galaxies in the nearby Universe, with an estimated stellar mass of 1011.21 M⊙ (Fall & Romanowsky 2018). It is classified as a lenticular galaxy (de Vaucouleurs et al. 1991), and its morphological features include a prominent dust lane, extensive lobes, and jets, which are indicative of its active galactic nucleus. The present study, part of the extragalactic GC NIR spectroscopy survey, is the first to investigate the age zero-point problem in the NIR, displaying it in the distinct environment of Cen A, thus providing new insights into NIR SPS modelling and the interpretation of age-sensitive indices.
The paper is organised as follows: Section 2 details the sample and data used in the study. Section 3 describes the stellar population models employed. In Section 4, we describe the spectral analysis, outlining our methodologies for data homogenisation and index measurements. The results of our analysis are presented in Section 5. Finally, we provide a summary of our findings and conclusions in Section 6.
2. Samples and data
Our investigation of the model zero-point problem makes use of different data sets, including NIR and optical spectra. In this section, we outline the selection criteria for GCs in the Cen A galaxy, for which NIR spectra were obtained, and summarise pre-existing optical spectral data. Subsequently, we describe the NIR and optical datasets for MW GCs with which we compare the GCs in Cen A.
2.1. Centaurus A GC data
We selected a sample of 21 GCs in the Cen A galaxy based on the optical catalogue by Beasley et al. (2008), which was designed to cover as many ages and metallicities as possible. As an age proxy, we chose Hβo, and we used the composite index [MgFe]′ (Thomas et al. 2003) as our total metallicity indicator, which was constructed from the combination of three indices: Mgb, Fe 5270, and Fe 5335. This index is insensitive to variations in the [α/Fe] abundance ratio (Vazdekis et al. 2015). The GCs were observed in the NIR using the TripleSpec4 (R = λ/Δλ ≃ 3500) instrument on the SOAR telescope during observing runs in March and April 2020 and 2021. The NIR spectra were reduced following standard procedures. We employed the Spextool IDL pipeline of Vacca et al. (2003), Cushing et al. (2004) for data reduction, incorporating flat-field correction, wavelength calibration, spatial extraction, telluric band correction, and flux calibration. The final spectra were corrected for the Doppler shift of Cen A and MW dust reddening, yielding rest frame spectra covering 0.96–2.46 microns. Further details of the observations and data reduction procedures are described in Paper I.
We also use the 21 optical spectra from the Beasley et al. (2008) survey of Cen A GCs. This survey was conducted using the 2-degree field instrument on the Anglo-Australian Telescope. This provided high-quality spectra over a wide field, which are essential for the metallicity, age, and kinematic analysis of the GC system. The spectra were obtained with a resolution of R ≈ 2000, covering the wavelength range of 3700–5800 Å.
2.2. Milky Way GC data
We incorporated the NIR spectroscopic dataset of MW GCs from Riffel et al. (2011) for comparative purposes. These data were selected from the Bica & Alloin (1986) catalogue, which was constructed in order to provide spectral observational constraints for calibrating evolutionary population models in the NIR. Riffel et al. (2011) obtained cross-dispersed NIR spectra for 12 GCs using the Ohio State Infrared Imager/Spectrometer mounted on the SOAR Telescope, covering the 1.2 micron to 2.35 micron region with a resolving power of R ∼ 1200. Their observational strategy and instrumental configuration were optimised to minimise the impact of seeing and aperture variations, yielding a reliable subset of 12 GCs with S/N ≥ 25 and spanning −1.62 ≤ [Z/H] ≤ −0.04 range. These clusters were originally chosen for their wide ranges of age and metallicity, which is relevant for testing evolutionary models.
We also performed a cross-match of 12 GGCs studied in Riffel et al. (2011) with the optical spectroscopy datasets of Schiavon et al. (2005), Kim et al. (2016), and Usher et al. (2017), which yielded 11 matches. These datasets feature wide-ranging wavelength coverage and relatively high spectral resolution. The library of integrated spectra of MW and Local Group GCs by Usher et al. (2017) covers a wavelength range of 3300 Å–9050 Å with a resolution of R = 6800. The wavelength coverage of integrated spectra of 24 GCs in Kim et al. (2016) is from 4000 to 5400 Å with a spectral resolution of FWHM ∼ 2.0 Å. The library of Schiavon et al. (2005) covers the range ∼3350–6430 Å with a resolution of ∼3.1 Å and integrates the full projected area within cluster core radii.
3. Stellar population models
Stellar population models are required to interpret observations of stellar clusters and galaxies. As each model is constructed using its own assumptions, methods, and distinct ingredients, variations in their outputs can naturally be anticipated. To determine whether the age zero-point problem is prevalent across various sets of models, we employed different SPS models. In the following, we summarise the construction and characteristics of these models.
-
Updated E–MILES (Röck et al. 2015, 2016; Vazdekis et al. 2016): The updated E–MILES models (see section 3 in La Barbera et al. 2024) are now used to make predictions in the NIR that incorporate both the original (Cushing et al. 2004; Rayner et al. 2009) and the Extended IRTF (E-IRTF) stellar library (Villaume et al. 2017). These models mainly improve the predictions in the NIR for metallicities below −0.4 dex. E–MILES wavelength coverage extends from 1680 to 50 000 Å. These models employ theoretical isochrones from the BaSTI (Pietrinferni et al. 2004) and Padova (Girardi et al. 2000) databases converted to the observational plane by employing extensive photometric stellar libraries (see Vazdekis et al. 2016 for a detailed description). The predictions in the NIR are considered reliable for ages above 1 Gyr (although in the optical range, the predictions are considered reliable for much younger populations). The SSP spectral predictions have a resolution of FWHM = 2.5 Å in the optical and R ≈ 2000 in the NIR. The E–MILES models also incorporate a variety of initial mass functions (IMFs), including the universal and revised Kroupa (Kroupa 2001), Chabrier (Chabrier 2001), and unimodal and bimodal IMFs (Vazdekis et al. 1996). Here, we use the updated E–MILES model predictions based on the BaSTI(Padova) isochrones and consider ages ranging from 5 to 13 Gyr (5.01 to 12.59 Gyr) and total metallicities spanning from as low as −2.27 (−2.32) up to the solar value of +0.06 (0.00). We adopt a low-tapered (< 0.6 M⊙) bimodal IMF with a logarithmic slope of 1.3, which almost resembles the standard Kroupa IMF (see Vazdekis et al. 1996).
-
α-enhanced E–MILES (La Barbera et al. 2017) In addition to the scaled-solar models, we also incorporate α-enhanced E-MILES models, which are an updated version of the Na-enhanced models described in La Barbera et al. (2017), to examine the impact of varying α-element abundances. For this study, we used two sets of models: one with [α/Fe] = 0 and another with [α/Fe] = +0.4. Both sets of models were computed using the Padova isochrones. The first set ([α/Fe] = 0) includes stellar populations with isochrones and a stellar library assuming solar-scaled abundances, covering ages from 5 to 13 Gyr and metallicities ranging from [Z/H] = −0.66 to +0.06. The second set ([α/Fe] = +0.4) consists of α-enhanced models with isochrones and a stellar library enhanced by [α/Fe] = +0.4, covering the same age range but with metallicities from [Z/H] = −0.35 to +0.06.
4. Spectral analysis
In this section, we describe the key steps of our analysis. In Sect. 4.1, we show the procedures undertaken to ensure homogeneity and consistency across the different spectral datasets. The aim of this processing step is to prepare the data for subsequent line-strength measurements and analysis. In Sect. 4.2, we present the selected indices, their significance, and the methodology adopted to measure them.
4.1. Data-set homogenisation
4.1.1. Centaurus A GC spectra
The optical spectra of Cen A GCs provided in the air wavelength system were in the observed frame. We therefore used the radial velocities of the GCs taken from Beasley et al. (2008) to bring the original optical spectra to the rest frame. The spectral resolution of these observations varied, with FWHM = 5.3 Å for spectrograph A and 4.0 Å for spectrograph B.
The NIR spectra of Cen A GCs were initially in the vacuum wavelength system, and therefore we transformed them to the air system by applying a conversion formula similar to the one used by SDSS. We then adjusted for the observed frame using the redshift of the Cen A galaxy (z = 0.00183, Graham 1978) and the radial velocities of GCs provided by Beasley et al. (2008), thereby placing these NIR spectra in the rest frame. In this study, our primary focus is on the Paβ1.28 absorption feature. Upon closely comparing the redshifted spectra of the GCs with a reference SSP model, we found that the dip of the Paβ feature in most GCs does not align exactly with the model in terms of wavelength. To resolve this, we applied a second-order wavelength correction, shifting the spectra around the Paβ feature to ensure they were precisely in the rest frame. It is crucial to have the Paβ absorption accurately positioned at the rest frame for line-strength analysis, as Eftekhari et al. (2021) demonstrated that this index is sensitive to wavelength shifts and an uncertainty of more than 92 km s−1 in wavelength calibration can lead to significant variation in the measured index value. Additionally, the NIR spectra of the Cen A GCs include error spectra.
4.1.2. Milky Way GC spectra
The optical spectra of the MW GCs had already been converted to the air wavelength system. The spectra from Schiavon et al. (2005) and Kim et al. (2016) were also provided in the rest frame. However, the spectra from Usher et al. (2017) were not. Therefore, we used the PPXF code (Cappellari & Emsellem 2004; Cappellari 2017) to derive the radial velocities of this sample and shift their spectra to the rest frame. The Usher et al. (2017) and Kim et al. (2016) samples include error spectra as well.
The NIR spectra of MW GCs were originally provided in the air system and in the rest frame. However, similar to the NIR spectra of Cen A GCs, the dip of the Paβ1.28 absorption feature did not exactly match the model. Therefore, we performed a second-order wavelength calibration around this index by applying a spectral shift to align the feature more precisely.
4.1.3. Resolution homogenisation
Before proceeding with line-strength measurements, it was imperative to ensure that all spectra shared a common resolution. To achieve this, we employed a pixel-by-pixel convolution method i.e. convolving each pixel individually to achieve the desired resolution) originally derived from the PPXF code (Cappellari & Emsellem 2004; Cappellari 2017). Specifically, all spectra in the optical data sets were convolved to a uniform resolution of σ = 180 km s−1, aligning them with the lowest resolution among the optical data. Similarly, for the NIR samples, we convolved all spectra to a resolution of σ = 106 km s−1, which corresponds to the lowest resolution among the NIR data. We note that for Cen A GCs data sets, we considered the intrinsic broadening of spectral features for each GC by incorporating the velocity dispersion information provided by Beasley et al. (2008). This homogenisation process is essential, as broadening can affect line-strength measurements.
Throughout these procedures, we carefully processed the spectra to ensure consistency across all data sets. As part of this process, we employed the SPECTRES Python tool (Carnall 2017) to resample the spectra and their associated uncertainties onto an arbitrary wavelength grid.
4.2. Index measurements
In stellar population studies, hydrogen spectral lines play a crucial role in characterising the age of stellar populations. The Balmer line Hβ, which corresponds to the transition from the n = 2 to n = 4 energy level in hydrogen at ∼4861 Å, is highly sensitive to the age of stellar populations. Other Balmer lines in the optical window, such as Hγ and Hδ, are less used for spectroscopic age estimates as they show sensitivity to variations in the abundance of α elements. One notable exception is the Hγσ index definition of Vazdekis & Arimoto (1999) and Vazdekis (2001), which requires extremely high signal-to-noise ratios (Vazdekis et al. 2015), and these are not satisfied by our GC spectra. In our study, we employ the optimised index Hβo, which is derived in a manner that maximises its sensitivity to age while minimising its dependence on metallicity (Cervantes & Vazdekis 2009). Therefore, the Hβo index is able to provide better age estimates for the GCs of Cen A, enabling us to explore the age zero-point problem in extragalactic GCs.
To address the age zero-point problem in the NIR, we rely on the Paschen line centred at 1.28 microns. We use the index definition of Eftekhari et al. (2021) as it is optimised to be maximally sensitive to the age of the stellar populations. Figure 3 in Eftekhari et al. (2021) shows a remarkable drop in the Paβ1.28 line strength with increasing age for populations older than 2 Gyr. Importantly, Paβ1.28 is notably more responsive to changes in age than to variations in metallicity or the IMF. Furthermore, it should be noted that the Paβ1.28 index exhibits a mild dependency on titanium abundance due to potential contamination of the central bandpass by Ti I lines (see Fig. A.1 in E21). Unfortunately, some of the spectra of our sample of Cen A GCs exhibit some artifacts within the central bandpass of the Paβ1.28 original index definition. To minimise the impact of these artifacts, we narrowed the central bandpass to a range of 12 805–12 835 Å. Despite this adjustment, five GCs still displayed artifacts within the central region of the Paβ1.28 index and for that reason, we excluded them from our analysis. A zoom onto the Paβ1.28 index in the spectra of the Cen A GCs is shown in Fig. 1, where the pseudo-continua and central bandpasses are marked in orange and grey, respectively. The spectra in the plot are already de-redshifted to the rest frame and smoothed to a velocity dispersion of 106 km s−1. It is worth noting that our study marks the first analysis of Paβ1.28 absorption in the NIR for GCs, and a further novelty is that we are analysing GCs in an extragalactic context.
![]() |
Fig. 1. Spectra of Cen A GCs in the region of the Paβ1.28 index. All spectra are presented in the rest frame and have been smoothed to a resolution of 106 km s−1. Error bars are shown as vertical lines for all spectra. The black spectrum on top represents the E–MILES SSP of age 11 Gyr, at Z = −0.66, with a Kroupa-like IMF. The central bandpass of the Paβ1.28 index, along with the blue and red pseudo-continua, are highlighted as grey and orange areas, respectively. These bandpasses are based on Eftekhari et al. (2021), with the central bandpass narrowed to 12 805–12 835 Å to avoid spectral artifacts within the index. |
We measured the line strengths of the aforementioned indices in both the observed and model spectra. Line-strength measurements were performed using the INDEXF code, a software tool developed by Cardiel (2010). As we have error spectra for Cen A GCs in the NIR, as well as GGCs from Usher et al. (2017) and Kim et al. (2016) datasets, the INDEXF code estimated line-strength errors using the provided error spectrum. The code employs error-propagation techniques – accounting for the noise characteristics present in the error spectrum – to determine the associated uncertainties in the line-strength measurements. For datasets where error spectra were unavailable, we used the INDEXF program to estimate random errors in the line-strength indices through a simulation-based approach. The program generated synthetic error spectra based on the original data, assuming randomly generated S/Ns (default range: 1−100). Using these synthetic error spectra, INDEXF performed 1000 Monte Carlo simulations (simulsn = 1000), introducing random variations following a Gaussian distribution. The program assumes these errors represent 1σ variations. For each simulation, the line-strength index was measured, the final index value was taken as the mean value of all simulations, and its standard deviation was taken as the uncertainty on the index line-strength value. It is important to note that for the GGCs common to the datasets of Usher et al. (2017), Kim et al. (2016), and Schiavon et al. (2005), we employed the average values of the indices.
5. Results
The model zero-point problem has been explored in the optical spectrum, primarily using MW GCs, but its origin remains a challenge to our understanding. Moving beyond the MW, we have expanded this research to include extragalactic GCs, focusing on those of the Cen A galaxy. First, we measured the optical indices Hβo and [MgFe]′, both expressed in Å, for our sample of Cen A GCs. The results are shown in the upper panel of Fig. 2, which shows the Hβo versus the total metallicity indicator. The blue points represent the measurements derived from the spectra of our sample of GCs in Cen A. The solid grid represents updated E–MILES SSP models with BaSTI isochrones. These models have total metallicities from −2.27 to solar. Each colour on the grid represents a metallicity, as described in the legend. Model ages are shown with grey scale and span from 5 to 13 Gyr. It is evident that the GCs of Cen A have a consistently lower Hβo value compared to the predictions of the E–MILES models with BaSTI isochrones (except for HGHH-34 and HHH86-30). The observed discrepancy between the measurements of the Cen A GCs and our reference model confirms that the model zero-point problem also appears in extragalactic GC systems (e.g. Brodie et al. 2005; Cenarro et al. 2007; Chies-Santos et al. 2012).
![]() |
Fig. 2. Line strengths of two hydrogen line indices, Hβo in the optical (upper panel) and Paβ1.28 in the NIR (lower panel), plotted against the total metallicity indicator [MgFe]′ for two samples of GCs in the Cen A galaxy (blue points) and MW (red squares). The overlaid grids represent predictions from updated E–MILES SSP models with BaSTI and Padova isochrones (solid and dotted lines), as well as E–MILES SSPs with [α/Fe] = 0.0 and +0.4 (dashed and dash-dotted lines). These models span a range of metallicities from Z = −2.27 (−2.32) to solar, as indicated by various colours, and ages ranging from 5 to 13 Gyr (12.59 Gyr) for SSPs with BaSTI (Padova) isochrones. |
To extend our analysis of the model zero-point problem to the NIR, we focus on the Paβ1.28 index. The lower panel in Fig. 2 shows Paβ1.28 line strengths against the total metallicity indicator, [MgFe]′. As in the upper panel, the blue data points represent measurements from the spectra of Cen A GCs. These points are overlaid with the predictions of the E–MILES SSP models with BaSTI isochrones for varying metallicities and ages (solid grid). Similar to the Hβo panel, a majority of the GCs in the Cen A galaxy are situated below the model grid in the NIR. This observation emphasises that the age zero-point problem is not confined to the optical range and reveals limitations within SPS models when extended into the NIR domain as well. Finally, it is remarkable that in this spectral range the model zero-point problem emerges across the entire range of metallicities in our sample of GCs. This observation aligns with previous findings in high- and intermediate-metallicity GCs (e.g. in 47 Tuc GCs, Gibson et al. 1999; Vazdekis 2001; Mendel et al. 2007). Vazdekis (2001) also emphasised the importance of investigating this age gap in metal-poor GCs ([Z/H] ≲ −0.9) and the necessity of enhancing stellar spectral libraries to improve their predictive accuracy for lower metallicity ranges.
It is also noteworthy that the age zero-point problem in both the optical and NIR is not exclusive to the E–MILES model. The XSL model, one of the latest SPS models with high resolution (R ∼ 10 000) and broad wavelength coverage from NUV to NIR, also exhibits the age zero-point problem. A detailed comparison of predictions from three recent and widely used SPS models is presented in Appendix A.
To test the impact of varying the stellar evolutionary prescriptions used in the E-MILES model, we compared grids corresponding to two distinct sets of isochrones: BaSTI and Padova. The grid of E–MILES models based on BaSTI isochrones is shown with solid lines, while the grid based on Padova isochrones is represented by dotted lines in Fig. 2. For both indices, the Padova-based grid consistently lies lower than the BaSTI-based grid. The BaSTI-based model predicts a sharper increase in the Hβo line strength for metallicities of Z < −0.66 compared to the Padova-based model at the oldest ages. However, the predictions for Paβ1.28 show minimal differences between the two isochrones. Similar to the E–MILEs models with BASTI isochrones, the models with Padova isochrones show the age zero-point problem for Hβo and Paβ1.28 indices. A detailed discussion of the impact of isochrones on the predictions of the XSL model for hydrogen lines compared to those of the E–MILES model is presented in Appendix B.
α-enhanced isochrones vary the temperature and luminosity of the stars at the TO point and therefore affect the hydrogen lines. According to Vazdekis (2001), incorporating α-enhanced isochrones into SPS models brings spectroscopic age estimates closer to those derived from CMDs. In order to examine the impact of α-enhancement on the Hβo and Paβ1.28 indices, we compared grids of E–MILES models with two different α-element abundances: [α/Fe] = 0 and +0.4. The grids of these two models are shown in Fig. 2 with dashed and dash-dotted lines, respectively. For the Hβo index, the grid with [α/Fe] = +0.4 predicts only slightly higher line strengths compared to the [α/Fe] = 0 grid. In contrast, the difference is more significant in the case of the Paβ1.28 index. The grid with [α/Fe] = +0.4 predicts markedly higher line strengths for the Paβ1.28 index compared to the [α/Fe] = 0 grid. This suggests that Paβ1.28 is more sensitive to variations in α-element abundances. Despite these differences, both α-enhanced and solar-scaled models continue to exhibit the age zero-point problem in the optical and NIR, and unlike in Vazdekis (2001) for the Hγ line, they do not help to reduce the gap. Therefore, the impact of α-enhancement requires further investigation.
We also include data points for MW GCs (red squares) in Fig. 2 to compare the behaviour of MW GCs with that of extragalactic GCs in the context of the age zero-point problem. Similar to Cen A GCs, most of the MW GCs fall below the E–MILES model grids for both the Hβo and Paβ1.28 indices, suggesting ages that are consistently older than the oldest age of the models and even exceeding the age of the Universe itself.
To further investigate the age zero-point problem in stellar population models, we explored the relationship between the Hβo and Paβ1.28 indices. Figure 3 illustrates this relationship for GCs from the Cen A and MW galaxies. In the plot, Cen A GCs are represented by filled circles, contrasting with the filled squares that represent MW GCs. The [MgFe]′ metallicity values are indicated through colour bars, each assigned a specific colour for the GCs of each galaxy. We overplot the E–MILES model grid with the BaSTI isochrones. We see that the models exhibit a correlation between Hβo and Paβ1.28, with more variation in Hβo than in Paβ1.28. The greater variability in Hβo is expected, given the different nature of the transitions: Hβo originates from the Balmer series involving transitions to the n = 2 level of hydrogen, while Paβ1.28 is part of the Paschen series, involving transitions to the n = 3 level. This results in distinct behaviours with age; the Paβ1.28 strength decreases sharply at approximately 1.5 Gyr, while Hβo shows a more gradual decline over time. We note that contrary to the model prediction, the data reveal a significant spread in Paβ1.28, while Hβo varies within a narrower range. Additionally, no evident distinction is observed between the Cen A and MW GCs in the plot, implying comparable ages for both galactic and extragalactic GCs. Moreover, the data do not show a clear trend with metallicity. A detailed comparison of the Hβo–Paβ1.28 relation among three distict SPS models is presented in Appendix A.
![]() |
Fig. 3. Hβo vs Paβ1.28 line strengths in GCs from the Cen A galaxy and the MW. In this plot, Cen A GCs are represented by filled circles and MW GCs by filled squares. The line strengths are plotted against the [MgFe]′ metallicity indicator, with a bluish colour bar denoting the [MgFe]′ values for Cen A GCs and a reddish colour bar for MW GCs. The overlaid grid shows the predictions of the E–MILES model with BaSTI isochrones. The model grid displays predictions across a range of metallicities for the distinct ages specified in the legend. |
One potential avenue to alleviate the age discrepancies in GCs is the inclusion of atomic diffusion in the isochrones. Vazdekis (2001) demonstrated that atomic diffusion can significantly alter the TO temperatures of low-mass stars at a relatively old age, yielding spectroscopic ages more consistent with those derived from CMDs. Furthermore, the findings of Vazdekis (2001) suggest that the contribution of the HB to the Balmer lines in metal-rich clusters is likely negligible, redirecting the focus towards other factors that might influence the TO temperature of main sequence stars in the models. Further investigation of this area is of paramount importance, but is beyond the scope of this paper; here we aim to show the intrinsic nature of this problem across different spectral ranges.
It is also worth mentioning the study by Schiavon et al. (2002) on the modelling of the luminosity function (LF) at the giant branch level for better reproduction of the integrated spectrum. Their analysis shows that using observed LFs instead of theoretical ones aligns spectroscopic ages with those derived from the CMD, countering the shortcomings of theoretical LFs, which often neglect AGB stars or underestimate bright giants. According to these authors, this miscalculation, which is particularly evident in the case of 47 Tuc, has previously led to exaggerated spectroscopic ages due to the neglect of factors like He diffusion and α-enhancement. A significant open question pertains to the applicability of the red giant excess observed in 47 Tuc – which is a high-metallicity GGC – across GCs of different metallicities and in extragalactic GCs with compositions distinct from those in the MW.
Finally, it is worth considering the impact of incorporating binary stars into SPS modelling. Mass transfer and stellar interactions in binary systems within a stellar population can modify the apparent age of the population by rejuvenating stars or shortening their lifetimes. Models such as BPASS (Stanway & Eldridge 2018), which account for binary evolution, suggest that binaries may lead to younger apparent ages in observed stellar populations. However, we note that BPASS models rely on theoretical spectral libraries, which appear to have difficulty in modelling NIR spectra, particularly for cool stars and molecular bands that dominate the NIR light. The models used in our study are based on empirical spectral libraries, providing reliable predictions in the NIR.
6. Summary and conclusions
This paper presents a pioneering application of the extragalactic globular cluster near-infrared spectroscopy survey, wherein we extend the investigation of the model zero-point problem to the NIR range. We focus on the Paβ line strength at 1.28 microns for GCs in the Cen A galaxy and compare this to values for the GCs in the MW. Our main findings can be summarised as follows.
-
Revealing the age zero-point problem in the NIR: The key conclusion of our study is the recognition of the model zero-point problem in the NIR spectrum by means of the Paβ1.28 index. Previously noted in the optical spectra of MW GCs, this discrepancy has now been identified in the NIR as the majority of the GCs in the Cen A galaxy are observed to appear older than the predictions of the stellar population models with the oldest ages, that is, older than the current age of the Universe. Our results highlight the intrinsic nature of the age zero-point problem across different spectral ranges.
-
The model zero-point problem is seen in extragalactic GCs in the NIR: The inclusion of MW GCs in this study reveals similar trends, showing that the model zero-point problem in the NIR applies to GCs outside the Galaxy.
-
Isochrone selection: Our results indicate that the choice of isochrones cannot be leveraged to solve the age zero-point problem. Both BaSTI and Padova isochrones in the E–MILES model show similar inconsistencies in the optical and NIR ranges.
-
α-enhancement: α-enhanced models still exhibit the age zero-point issue in both optical and NIR ranges. Contrary to some expectations, α-enhancement did not reduce the gap between the oldest models and GCs, highlighting the need for further investigation into its influence on stellar populations predictions.
In summary, our study highlights the challenge of determining GC ages using hydrogen lines. Our study presents insights into the age zero-point problem in SPS models, confirming its reach into the NIR spectrum. This problem is not limited to a specific galactic context but is a broader challenge affecting both galactic and extragalactic GCs. Our analysis shows that this issue is not resolved by varying isochrones or introducing α-enhanced models, as both approaches still exhibit the same discrepancies. With the advent of observing facilities such as the Extremely Large Telescope and the James Webb Space Telescope, which are equipped with enhanced NIR capabilities, we have an unprecedented opportunity to observe extragalactic GCs with greater precision. Their precise data will allow refined calibrations of SPS models, improving our understanding of stellar evolution and the age zero-point problem.
Acknowledgments
We thank the anonymous referee for their insightful comments and suggestions, which have greatly improved the quality and clarity of this manuscript. EE and AV acknowledge support from grants PID2021-123313NA-I00 and PID2022-140869NB-I00 from the Spanish Ministry of Science, Innovation and Universities MCIU. This work has also been supported 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. EE and AV also acknowledge support from grant POKEBOWL PID2021-123313NA-I00 from the MCIN/AEI and the European Regional Development Fund (ERDF). RR acknowledges support from the Fundación Jesús Serra and the Instituto de Astrofísica de Canarias under the Visiting Researcher Programme 2023-2025 agreed between both institutions. RR, also acknowledges support from the ACIISI, Consejería de Economía, Conocimiento y Empleo del Gobierno de Canarias and the European Regional Development Fund (ERDF) under grant with reference ProID2021010079, and the support through the RAVET project by the grant PID2019-107427GB-C32 from the Spanish Ministry of Science, Innovation and Universities MCIU. RR also thanks to Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, Proj. 311223/2020-6, 304927/2017-1 and 400352/2016-8), Fundação de amparo à pesquisa do Rio Grande do Sul (FAPERGS, Proj. 16/2551-0000251-7 and 24/2551-0001282-6), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES, Proj. 0001). LGDH acknowledges support by National Key R&D Program of China No. 2022YFF0503402, and National Natural Science Foundation of China (NSFC) project number E345251001. ACS acknowledges funding from the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), the Rio Grande do Sul Research Foundation (FAPERGS) and the Chinese Academy of Sciences (CAS) President’s International Fellowship Initiative (PIFI) through grants CNPq-11153/2018-6, CNPq-314301/2021-6, FAPERGS/CAPES 19/2551-0000696-9, E085201009.
References
- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Baldwin, C., McDermid, R. M., Kuntschner, H., Maraston, C., & Conroy, C. 2018, MNRAS, 473, 4698 [Google Scholar]
- Beasley, M. A., Bridges, T., Peng, E., et al. 2008, MNRAS, 386, 1443 [NASA ADS] [CrossRef] [Google Scholar]
- Bica, E., & Alloin, D. 1986, A&A, 162, 21 [NASA ADS] [Google Scholar]
- Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Brodie, J. P., Strader, J., Denicoló, G., et al. 2005, AJ, 129, 2643 [NASA ADS] [CrossRef] [Google Scholar]
- Cappellari, M. 2017, MNRAS, 466, 798 [Google Scholar]
- Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 [Google Scholar]
- Cardiel, N. 2010, Astrophysics Source Code Library [record ascl:1010.046] [Google Scholar]
- Carnall, A. C. 2017, ArXiv e-prints [arXiv:1705.05165] [Google Scholar]
- Cenarro, A. J., Beasley, M. A., Strader, J., Brodie, J. P., & Forbes, D. A. 2007, AJ, 134, 391 [NASA ADS] [CrossRef] [Google Scholar]
- Cenarro, A. J., Cervantes, J. L., Beasley, M. A., Marín-Franch, A., & Vazdekis, A. 2008, ApJ, 689, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Cervantes, J. L., & Vazdekis, A. 2009, MNRAS, 392, 691 [NASA ADS] [CrossRef] [Google Scholar]
- Chabrier, G. 2001, ApJ, 554, 1274 [NASA ADS] [CrossRef] [Google Scholar]
- Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525 [Google Scholar]
- Chen, Y., Bressan, A., Girardi, L., et al. 2015, MNRAS, 452, 1068 [Google Scholar]
- Chies-Santos, A. L., Larsen, S. S., Kuntschner, H., et al. 2011, A&A, 525, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Chies-Santos, A. L., Larsen, S. S., & Kissler-Patig, M. 2012, MNRAS, 427, 2349 [NASA ADS] [CrossRef] [Google Scholar]
- Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102 [Google Scholar]
- Conroy, C., Villaume, A., van Dokkum, P. G., & Lind, K. 2018, ApJ, 854, 139 [Google Scholar]
- Cushing, M. C., Vacca, W. D., & Rayner, J. T. 2004, PASP, 116, 362 [NASA ADS] [CrossRef] [Google Scholar]
- Dahmer-Hahn, L. G., Riffel, R., Rodríguez-Ardila, A., et al. 2018, MNRAS, 476, 4459 [Google Scholar]
- Dahmer-Hahn, L. G., Chies-Santos, A. L., Eftekhari, E., et al. 2024, A&A, 690, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., et al. 1991, Third Reference Catalogue of Bright Galaxies (New York: Springer) [Google Scholar]
- Dotter, A. 2016, ApJS, 222, 8 [Google Scholar]
- Eftekhari, E., Vazdekis, A., & La Barbera, F. 2021, MNRAS, 504, 2190 [CrossRef] [Google Scholar]
- Eftekhari, E., La Barbera, F., Vazdekis, A., Allende Prieto, C., & Knowles, A. T. 2022, MNRAS, 512, 378 [NASA ADS] [CrossRef] [Google Scholar]
- Fall, S. M., & Romanowsky, A. J. 2018, ApJ, 868, 133 [Google Scholar]
- Gibson, B. K., Madgwick, D. S., Jones, L. A., Da Costa, G. S., & Norris, J. E. 1999, AJ, 118, 1268 [NASA ADS] [CrossRef] [Google Scholar]
- Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Graham, J. A. 1978, PASP, 90, 237 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, H.-S., Cho, J., Sharples, R. M., et al. 2016, ApJS, 227, 24 [NASA ADS] [CrossRef] [Google Scholar]
- Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
- La Barbera, F., Vazdekis, A., Ferreras, I., et al. 2017, MNRAS, 464, 3597 [Google Scholar]
- La Barbera, F., Vazdekis, A., Pasquali, A., et al. 2024, A&A, 687, A156 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Leath, H. J., Beasley, M. A., Vazdekis, A., Salvador-Rusiñol, N., & Gvozdenko, A. 2022, MNRAS, 512, 548 [NASA ADS] [CrossRef] [Google Scholar]
- Lee, H.-C., & Worthey, G. 2005, ApJS, 160, 176 [NASA ADS] [CrossRef] [Google Scholar]
- Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, ApJ, 804, 64 [Google Scholar]
- Marigo, P., Bressan, A., Nanni, A., Girardi, L., & Pumo, M. L. 2013, MNRAS, 434, 488 [Google Scholar]
- Mármol-Queraltó, E., Cardiel, N., Sánchez-Blázquez, P., et al. 2009, ApJ, 705, L199 [CrossRef] [Google Scholar]
- Mendel, J. T., Proctor, R. N., & Forbes, D. A. 2007, MNRAS, 379, 1618 [NASA ADS] [CrossRef] [Google Scholar]
- Pastorelli, G., Marigo, P., Girardi, L., et al. 2019, MNRAS, 485, 5666 [Google Scholar]
- Pastorelli, G., Marigo, P., Girardi, L., et al. 2020, MNRAS, 498, 3283 [Google Scholar]
- Percival, S. M., & Salaris, M. 2009, ApJ, 703, 1123 [Google Scholar]
- Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 [Google Scholar]
- Prugniel, P., Vauglin, I., & Koleva, M. 2011, A&A, 531, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rayner, J. T., Cushing, M. C., & Vacca, W. D. 2009, ApJS, 185, 289 [Google Scholar]
- Riffel, R., Ruschel-Dutra, D., Pastoriza, M. G., et al. 2011, MNRAS, 410, 2714 [NASA ADS] [CrossRef] [Google Scholar]
- Riffel, R., Rodríguez-Ardila, A., Brotherton, M. S., et al. 2019, MNRAS, 486, 3228 [Google Scholar]
- Röck, B., Vazdekis, A., Peletier, R. F., Knapen, J. H., & Falcón-Barroso, J. 2015, MNRAS, 449, 2853 [CrossRef] [Google Scholar]
- Röck, B., Vazdekis, A., Ricciardelli, E., et al. 2016, A&A, 589, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rosenfield, P., Marigo, P., Girardi, L., et al. 2016, ApJ, 822, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
- Schiavon, R. P., Faber, S. M., Rose, J. A., & Castilho, B. V. 2002, ApJ, 580, 873 [NASA ADS] [CrossRef] [Google Scholar]
- Schiavon, R. P., Rose, J. A., Courteau, S., & MacArthur, L. A. 2005, ApJS, 160, 163 [NASA ADS] [CrossRef] [Google Scholar]
- Sharma, K., Prugniel, P., & Singh, H. P. 2016, A&A, 585, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Stanway, E. R., & Eldridge, J. J. 2018, MNRAS, 479, 75 [NASA ADS] [CrossRef] [Google Scholar]
- Strader, J., Brodie, J. P., Cenarro, A. J., Beasley, M. A., & Forbes, D. A. 2005, AJ, 130, 1315 [NASA ADS] [CrossRef] [Google Scholar]
- Tang, J., Bressan, A., Rosenfield, P., et al. 2014, MNRAS, 445, 4287 [NASA ADS] [CrossRef] [Google Scholar]
- Thomas, D., Maraston, C., & Bender, R. 2003, MNRAS, 339, 897 [NASA ADS] [CrossRef] [Google Scholar]
- Thomas, D., Maraston, C., & Korn, A. 2004, MNRAS, 351, L19 [NASA ADS] [CrossRef] [Google Scholar]
- Usher, C., Pastorello, N., Bellstedt, S., et al. 2017, MNRAS, 468, 3828 [NASA ADS] [CrossRef] [Google Scholar]
- Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389 [NASA ADS] [CrossRef] [Google Scholar]
- Vazdekis, A. 2001, Ap&SS, 276, 839 [NASA ADS] [CrossRef] [Google Scholar]
- Vazdekis, A., & Arimoto, N. 1999, ApJ, 525, 144 [NASA ADS] [CrossRef] [Google Scholar]
- Vazdekis, A., Casuso, E., Peletier, R. F., & Beckman, J. E. 1996, ApJS, 106, 307 [Google Scholar]
- Vazdekis, A., Sánchez-Blázquez, P., Falcón-Barroso, J., et al. 2010, MNRAS, 404, 1639 [NASA ADS] [Google Scholar]
- Vazdekis, A., Coelho, P., Cassisi, S., et al. 2015, MNRAS, 449, 1177 [Google Scholar]
- Vazdekis, A., Koleva, M., Ricciardelli, E., Röck, B., & Falcón-Barroso, J. 2016, MNRAS, 463, 3409 [Google Scholar]
- Verro, K., Trager, S. C., Peletier, R. F., et al. 2022a, A&A, 661, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Verro, K., Trager, S. C., Peletier, R. F., et al. 2022b, A&A, 660, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Villaume, A., Conroy, C., Johnson, B., et al. 2017, ApJS, 230, 23 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Comparison of different SPS models
Discrepancies in GCs age predictions by SSP models in the optical have been observed in the past. Mendel et al. (2007) focused on measuring ages, metallicities, and α-element abundances for GGCs using a variety of stellar population models (namely Thomas et al. 2004; Lee & Worthey 2005; Vazdekis et al. 2010). They found notable discrepancies in the ages predicted by different SPS models based on fitting simultaneously as many spectral indices as possible. Specifically, the average ages derived were 10.74 ± 1.84 Gyr for the Thomas et al. (2004) model, 9.38 ± 1.82 Gyr for the Lee & Worthey (2005) model, and 11.70 ± 3.60 Gyr for the Vazdekis et al. (2010) model. These findings emphasise the need for additional calibration of SSP models, particularly when interpreting age–sensitive indices.
In addition to our examination of the age zero-point problem within the framework of the E–MILES SSP models, we extend our investigation to encompass other recent stellar population models. Below, we summarise the construction and characteristics of these models:
-
Verro et al. (2022a) (hereafter XSL models): The XSL SSP models employ 639 star spectra from the X–shooter Spectral Library (XSL) DR3 (Verro et al. 2022b), which covers a wide range of stellar parameters (effective temperature, surface gravity, and metallicity). These models offer a broad wavelength spectrum from the U to the K bands at moderately high spectral resolution (R∼10000). The XSL SSP models employ two sets of isochrones: the latest PARSEC/COLIBRI (P/C) isochrones (Bressan et al. 2012; Tang et al. 2014; Chen et al. 2014, 2015; Marigo et al. 2013; Rosenfield et al. 2016; Pastorelli et al. 2019, 2020), which provide an advanced treatment of stellar evolutionary stages, particularly the thermally pulsing asymptotic giant branch (TP-AGB) phase, and the Padova isochrones (Girardi et al. 2000), which offer a simpler approach, lacking a third dredge-up and thus excluding the transition to carbon-rich TP-AGB stars. For this study, we selected a representative set of SSPs with P/C (Padova) isochrones, computed with Kroupa IMF, with ages between 5 and 15.8 Gyr (5.01 to 12.6 Gyr), whilst the total metallicity ranges from -1.4 to 0.0 dex (-1.68 to 0.00 dex).
-
Conroy et al. (2018) (hereafter C18): We also incorporated the grids of base SSP models from the C18 SPS model and fitter ALF. These SSPs were constructed using standard isochrone synthesis techniques using the MIST series of isochrones (Choi et al. 2016; Dotter 2016). The abundance scale of the MIST isochrones is set by the Asplund et al. (2009) value of Z⊙ = 0.0142. The Spectral Polynomial Interpolator (SPI; Villaume et al. 2017) was used to assign a stellar spectrum to each isochrone point. SPI was trained on the MILES and E–IRTF stellar libraries, using stellar parameters from (Prugniel et al. 2011) and (Sharma et al. 2016). The training set was further supplemented by empirical spectra from the (Mann et al. 2015) M dwarf library. The C18 models offer options for flexible IMFs, including variable power–law indices and low–mass cutoffs. Specifically, for our research purposes, we utilised the C18 models with stellar ages from 5 to 13.5 Gyr and within the metallicity range of -1.5 to 0.0 dex. In order to represent the MW IMF, we adopted an IMF parametrisation with a dual–slope power–law structure: x1 = 1.3 for the mass range 0.08-0.5 M⊙ and x2 = 2.3 from 0.5 to 1 M⊙. Beyond 1 M⊙, the slope adheres to the Salpeter IMF (Salpeter 1955).
It is noteworthy to mention that the isochrones of all these models are solar–scaled, such that [Z/H]=[Fe/H] for the base SSPs used in this work.
Figure A.1 shows the Hβo and Paβ1.28 line-strengths for GCs in the Cen A galaxy (blue points), as a function of the [MgFe]′ total metallicity indicator. Three distinct model grids are overlaid: E–MILES SSP model with BaSTI isochrone (in orange), XSL SSP model with P/C isochrone (in pink) and C18 with MIST isochrone (in green). The three models have different predictions for the Hβo index in the metal-poor regime, in part due to the differing treatment of the HB phase and mass loss. Most of the Cen A GCs fall below the XSL model grid, similar to the E-MILES grid. However, the oldest SSPs from the C18 model predict the lowest values for Hβo, and, as a result, most of the Cen A GCs fall within the C18 grid. Another interesting aspect of Fig. A.1 lies in the divergent predictions of the E–MILES, XSL, and C18 models in the NIR spectrum, as evidenced by the distinct shapes of their respective grids. The C18 model predicts the lowest line strengths for the Paβ1.28 index compared to other models and encompasses nearly half of the Cen A GC sample. However, a few GCs fall above the C18 models with an age of 5 Gyr, which contrasts with the much older ages derived for these GCs by Beasley et al. (2008) using optical data. Furthermore, an internal inconsistency is observed between the ages determined with these models using the Hβo and Paβ1.28 panels, as the Paβ1.28 panel generally indicates younger ages for the GCs. Similar to the E-MILES models, only a few GCs fall within the XSL grid in the Paβ1.28 panel. Therefore, the predictions of the E–MILES and XSL models show consistency between the optical and NIR spectral ranges. Another interesting point is that, in the metal–poor regime ([MgFe]′ < 2.2 Å), the E-MILES model shows contrasting behavior for the Hβo and Paβ1.28 indices. As metallicity decreases, the strength of Hβo increases, while the strength of Paβ1.28 decreases. This could be attributed to the different dependencies of Hβo and Paβ1.28 on HB stars. In old, metal–poor ([Z/H] < −1 dex) populations, the morphology of the HB includes stars along the blue, intermediate, and red HB. The stars contributing mostly to Hβo in the optical range are the blue and intermediate HB stars, which can be hotter than the TO stars, thus strengthening the Hβo index. Conversely, the stars contributing mostly to Paβ1.28 in the J band are the redder HB stars, which are cooler than the TO stars, thereby weakening the Paβ1.28 index.
Similar to Cen A GCs, most of the MW GCs fall below the E–MILES model grid for both the Hβo and Paβ1.28 indices, suggesting ages that are consistently older than the oldest age of the models and even exceeding the age of the Universe itself. However, the trend differs for the XSL and C18 models. For the Hβo index, only three MW GCs exhibit weaker absorption than the oldest SSPs in the XSL model (one in the case of the C18 model). For the Paβ1.28 index, all reliable MW GC measurements fall below the XSL model predictions, similar to the E–MILES model. In contrast, for the C18 model, only two MW GCs fall below the model grid.
![]() |
Fig. A.1. Hβo and Paβ1.28 line strengths of GCs from the Cen A galaxy (blue points) and the MW (red squares) against the metallicity indicator index, [MgFe]′. Predictions from various SPS models are superimposed: C18 with MIST isochrones in green, XSL with P/C isochrones in pink, and E–MILES with BaSTI isochrones in orange. Different lines within each colour represent predictions for varying ages, as the legends indicate. |
In Fig. A.2, we explored the relationship between the Hβo and Paβ1.28 indices in different SPS models. The metallicity range of each model is as following: BaSTI-based E–MILES models, -1.49< Z< 0.06, P/C-based XSL model, -1.4< Z< 0.0, and MIST-based C18 model, -1.5< Z< 0.0. The E–MILES model shows an increase in the Hβo index for the oldest SSPs with decreasing metallicity, contrasting with the behaviour of the XSL and C18 models. This discrepancy may stem from different approaches to modelling HB morphology in metal–poor SSPs. In older metal–poor populations, the prevalence of hot HB stars leads to increased Balmer line strengths, causing these populations to appear younger. The complex interplay of factors like mass loss, metallicity, and dynamical effects makes detailed HB morphology modelling challenging. The XSL and C18 models likely have employed a more conservative mass loss approach, resulting in cooler evolved stars, unlike the E–MILES model.
![]() |
Fig. A.2. Comparative analysis of Hβo and Paβ1.28 line strengths in GCs from the Cen A galaxy and the MW. In this plot, Cen A GCs are represented by filled circles and MW GCs by filled squares. The line strengths are plotted against the [MgFe]′ metallicity indicator, with a bluish colourbar denoting the [MgFe]′ values for Cen A GCs and a reddish colourbar for MW GCs. Predictive grids from three distinct models are overlaid: C18 with MIST isochrones (green), XSL with P/C isochrones (pink), and E–MILES with BaSTI isochrones (orange). Each model grid displays predictions across a range of metallicity for distinct ages specified in the legend. |
Appendix B: Comparison of models with different isochrones
In Fig. B.1, we compared two different sets of isochrones for the XSL model: PARSEC/COLIBRI (grid with solid lines) and Padova (grid with dotted lines). It is important to mention that the C18 model is based solely on one set of isochrones and is therefore not included in this analysis. The upper panel of Fig. B.1 shows model predictions for the Hβo index versus the [MgFe]′ indicator, and the lower panel shows the same but for the Paβ1.28 index. For both indices, the grid based on the Padova isochrones consistently lies lower than the one based on the P/C isochrones. The oldest SSP with the Padova isochrone predicts a sudden increase in the Hβo line-strength at Z < -0.6, whereas the oldest SSP with the P/C isochrone shows this increase at a lower metallicity of −1.0 dex. Similar to P/C isochrone, the XSL model with Padova isochrone shows the age zero-point problem too.
A particularly notable difference arises in the shape of the Paβ1.28 grid when using the Padova-based XSL model. In the lower panel of Fig. B.1, the Padova grid exhibits a completely different shape compared to the P/C grid, and most of the GCs in our sample fall within the Padova grid. This result contrasts with the behavior of the E–MILES model, where the Padova isochrone does not show significant differences from the BaSTI grid for the Paβ1.28 index (see Fig. 2). The consistency between the E–MILES grids across different isochrones provides more stable results.
![]() |
Fig. B.1. Comparison of different isochrones. Hβo (upper panels) and Paβ1.28 (lower panels) line strengths of GCs from the Cen A galaxy (blue points) and the MW (red squares) are plotted against the metallicity indicator [MgFe]′. XSL models with P/C and Padova isochrones are shown with solid and dotted grids, respectively. |
All Figures
![]() |
Fig. 1. Spectra of Cen A GCs in the region of the Paβ1.28 index. All spectra are presented in the rest frame and have been smoothed to a resolution of 106 km s−1. Error bars are shown as vertical lines for all spectra. The black spectrum on top represents the E–MILES SSP of age 11 Gyr, at Z = −0.66, with a Kroupa-like IMF. The central bandpass of the Paβ1.28 index, along with the blue and red pseudo-continua, are highlighted as grey and orange areas, respectively. These bandpasses are based on Eftekhari et al. (2021), with the central bandpass narrowed to 12 805–12 835 Å to avoid spectral artifacts within the index. |
In the text |
![]() |
Fig. 2. Line strengths of two hydrogen line indices, Hβo in the optical (upper panel) and Paβ1.28 in the NIR (lower panel), plotted against the total metallicity indicator [MgFe]′ for two samples of GCs in the Cen A galaxy (blue points) and MW (red squares). The overlaid grids represent predictions from updated E–MILES SSP models with BaSTI and Padova isochrones (solid and dotted lines), as well as E–MILES SSPs with [α/Fe] = 0.0 and +0.4 (dashed and dash-dotted lines). These models span a range of metallicities from Z = −2.27 (−2.32) to solar, as indicated by various colours, and ages ranging from 5 to 13 Gyr (12.59 Gyr) for SSPs with BaSTI (Padova) isochrones. |
In the text |
![]() |
Fig. 3. Hβo vs Paβ1.28 line strengths in GCs from the Cen A galaxy and the MW. In this plot, Cen A GCs are represented by filled circles and MW GCs by filled squares. The line strengths are plotted against the [MgFe]′ metallicity indicator, with a bluish colour bar denoting the [MgFe]′ values for Cen A GCs and a reddish colour bar for MW GCs. The overlaid grid shows the predictions of the E–MILES model with BaSTI isochrones. The model grid displays predictions across a range of metallicities for the distinct ages specified in the legend. |
In the text |
![]() |
Fig. A.1. Hβo and Paβ1.28 line strengths of GCs from the Cen A galaxy (blue points) and the MW (red squares) against the metallicity indicator index, [MgFe]′. Predictions from various SPS models are superimposed: C18 with MIST isochrones in green, XSL with P/C isochrones in pink, and E–MILES with BaSTI isochrones in orange. Different lines within each colour represent predictions for varying ages, as the legends indicate. |
In the text |
![]() |
Fig. A.2. Comparative analysis of Hβo and Paβ1.28 line strengths in GCs from the Cen A galaxy and the MW. In this plot, Cen A GCs are represented by filled circles and MW GCs by filled squares. The line strengths are plotted against the [MgFe]′ metallicity indicator, with a bluish colourbar denoting the [MgFe]′ values for Cen A GCs and a reddish colourbar for MW GCs. Predictive grids from three distinct models are overlaid: C18 with MIST isochrones (green), XSL with P/C isochrones (pink), and E–MILES with BaSTI isochrones (orange). Each model grid displays predictions across a range of metallicity for distinct ages specified in the legend. |
In the text |
![]() |
Fig. B.1. Comparison of different isochrones. Hβo (upper panels) and Paβ1.28 (lower panels) line strengths of GCs from the Cen A galaxy (blue points) and the MW (red squares) are plotted against the metallicity indicator [MgFe]′. XSL models with P/C and Padova isochrones are shown with solid and dotted grids, respectively. |
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.