Open Access
Issue
A&A
Volume 691, November 2024
Article Number A311
Number of page(s) 12
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202451174
Published online 22 November 2024

© The Authors 2024

Licence Creative CommonsOpen 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

Thanks to the photometric and spectroscopic campaigns that had been performed over the last twenty years, the picture of the 3D structure, chemistry and kinematics of the Galactic bulge has improved substantially, suggesting a composite formation process. However, no general consensus has been reached yet about the exact mechanisms that generated and shaped the bulge in the Milky Way and analogous galaxies.

Two main formation channels are currently invoked to reproduce the observed chemodynamical properties of the Milky Way’s inner region: on the one hand, its barred and X-shaped morphology (Wegg & Gerhard 2013; Wegg et al. 2015), together with a global cylindrical kinematics (Kunder et al. 2012; Ness et al. 2013; Zoccali et al. 2014; Rojas-Arriagada et al. 2020), would suggest a dynamical secular evolution of the Galactic disk driven by bar and buckling instabilities (Combes et al. 1990; Patsis et al. 2002; Athanassoula 2005; Shen et al. 2010; Fragkoudi et al. 2017; Debattista et al. 2017); on the other hand, the dominant old ages (>10 Gyr, Clarkson et al. 2008, 2011; Valenti et al. 2013; Renzini et al. 2018; Surot et al. 2019) and the chemical pattern of its stars (e.g., Zoccali et al. 2017; Rojas-Arriagada et al. 2017, 2020; Queiroz et al. 2021) requires an intense in-situ star formation occurred in a very short timescale (1–2 Gyr), as predicted by chemical evolution models (see Rojas-Arriagada et al. 2017; Matteucci 2021).

One of the most promising ways to power this second formation path – a first and vigorous assembly phase of the bulge – is the merging of giant star-forming clumps generated by disk instabilities. Indeed, it has been theoretically shown that turbulent gas-rich disks can fragment into massive substructures of the order of 107−109 M (e.g., Behrendt et al. 2016). If long-lived enough (mainly determined by different feedback prescriptions), these massive structures can gravitationally interact, causing a redistribution of angular momentum outward and migration of clumps and gas inward. Numerical simulations demonstrate that this process forms a bulge in a characteristic timescale of ~1 Gyr for a Milky Way-like galaxy (Noguchi 1999; Immeli et al. 2004a,b; Bournaud et al. 2007; Carollo et al. 2007; Elmegreen et al. 2008; Dekel et al. 2009). Clumps developed by hydrodynamical simulations are also characterized by high star formation rate (SFR) densities (Clarke et al. 2019) and provide a natural way of reproducing the chemical trend displayed by bulge field stars (Debattista et al. 2023). Moreover, the clumpy morphology is a common feature of star-forming galaxies at cosmic-noon (see, e.g., Elmegreen et al. 2009; Guo et al. 2015; Shibuya et al. 2016). The majority of the observed clumps substantially contribute to the rest-frame UV light of their host galaxies (up to ~40%), denoting very high SFRs (e.g., Zanella et al. 2015; Messa et al. 2022); however, a fraction of redder substructures is also detected (Guo et al. 2015). The measured stellar mass values range between 105 and 109 M, with sizes from a few to hundred parsecs (Dessauges-Zavadsky et al. 2017; Messa et al. 2022; Vanzella et al. 2022b,a; Claeyssens et al. 2023).

Yet, the actual contribution of high-redshift clumps to the bulge assembly is still a debated open question in both the galactic and extragalactic scientific communities. On the basis of current observations, it is very challenging to robustly evaluate if clumps survive long enough to migrate to the central region of the galaxy (e.g., Cava et al. 2018). On the other hand, recent theoretical studies show that not only the massive clumps are able to migrate to the center, but some of them are also able to survive and to generate stellar clusters (Bournaud 2016; Dekel et al. 2023). It is likely that the primordial clumps were massive enough to retain iron-enriched supernova ejecta, creating a reservoir of gas to fuel, in a self-enrichment fashion, multiple bursts of star formation. This process would give rise to multi-iron and multi-age stellar subpopulations, at odds with what observed in genuine globular clusters (GCs). These in-situ stellar systems would be the empirical witnesses of the hierarchical assembly of the Galactic bulge, and are therefore called Bulge Fossil Fragments (BFFs, Ferraro et al. 2009, 2021). Thus, the detection of BFFs in the bulge would represent a key observational evidence in support to the involvement of clumps in building bulges, which is fundamental for our understanding of the formation of Milky Way-like galaxies in a broad cosmological context (i.e., to guide the formation of these galaxies implemented in current cosmological simulations; to have constraints on the gas fraction, turbulence and kinematics of high-redshift disk galaxies; to understand the feedback process prescriptions; etc.).

To date, only two BFF candidates have been detected in the bulge: Terzan 5 (Ferraro et al. 2009, 2016) and Liller 1 (Ferraro et al. 2021). As suggested by isochrone fits, these massive (few 106 M, Lanzoni et al. 2010; Saracino et al. 2015) stellar systems show a color-magnitude diagram (CMD) simultaneously populated by old (~12 Gyr) and young (up to 1 Gyr) stellar populations. In addition, spectroscopic analyses of these systems demonstrate a tight chemical similarity with bulge field stars. Indeed, their old stars are metal-rich ([Fe/H] ≥ −0.5 dex) and α-enhanced ([α/Fe] ~ +0.3 dex), while young stars are characterized by super-solar metallicity ([Fe/H] ~ +0.3 dex) and solar-scaled [α/Fe] abundance ratio (Origlia et al. 2011, 2013, 2019; Massari et al. 2014; Crociati et al. 2023; Alvarez Garay et al. 2024). Moreover, recent chemical evolutionary models specifically computed for the case of Terzan 5 (Romano et al. 2023) nicely reproduce all the chemical patterns observed in this stellar system in the framework of the self-enrichment evolution of a progenitor with an initial mass of 4 × 107 M. However, alternative formation models have been proposed for the origin of these two stellar systems. In particular, McKenzie & Bekki (2018) and Bastian & Pfeffer (2022) suggested that the young stellar component was originated from the accretion of a field Giant Molecular Clouds (GMCs) by a genuine GC Indeed, according to current numerical simulations (McKenzie & Bekki 2018), these should be extremely rare and fine-tuned events that would give rise to a single and brief star formation burst. Hence, this scenario would predict only two metallicity components and it would be clearly incompatible with a multiple burst and prolonged star formation history Bastian & Pfeffer 2022.

Determining the star formation histories (SFHs) of Terzan 5 and Liller 1 therefore represents a crucial step toward assessing their true origin. In the case of resolved stellar populations, a proper recovering of the SFH is obtained from synthetic CMD fitting methods (e.g., Tosi et al. 1991; Tolstoy & Saha 1996; Dolphin 1997, 2002; Hernandez et al. 1999; Cignoni & Tosi 2010; Aparicio & Hidalgo 2009; Weisz et al. 2012; Ruiz-Lara et al. 2018). These aim to determine key properties of composite stellar populations, such as age, mass, metallicity, and SFR, from the features readable in the CMD. Indeed, this analysis can unlock an amount of crucial information – namely, the duration of the star formation episodes, the presence of intermediate and barely detectable star formation bursts, and the build-up of the stellar mass throughout the lifetime of the systems – that cannot be obtained by simply performing isochrone fitting.

Within this context, in Dalessandro et al. (2022) we derived the SFH of Liller 1. Interestingly, we measured three main and prolonged (Δt > 1 Gyr) star formation events started at ~13,9 and 3 Gyr ago, which generated the 75, 15 and 10% of the system’s mass, respectively. These results are in line with those expected from the BFF scenario, reinforcing the hypothesis that Liller 1 is the surviving remnant of a much more massive primordial structure that experienced self-enrichment and likely contributed to the Galactic bulge formation. Here we present an analogous analysis on Terzan 5. Previous age-dating studies based on isochrone fitting (Ferraro et al. 2016) allowed the detection of an old and a young main sequence turn off (MS-TO), measuring an age of 12 ± 1 Gyr and 4.5 ± 0.5 Gyr, respectively. In fact, as in the case of Liller 1, to properly assess the total number of bursts and their time duration, as well as the overall efficiency of the star formation activity experienced by the system it is necessary to reconstruct the entire SFH. The manuscript is structured as follows: in Section 2 we describe all the steps we performed to obtain the observed CMD of Terzan 5, including the photometric reduction, the correction for the effects of differential reddening and the evaluation of the observational uncertainties; in Section 3 we present the adopted method for reconstructing the SFH and the main outcomes; in Section 4 we summarize and critically evaluate our results.

2 The input photometric catalog

2.1 Observations and data reduction

As demonstrated in Dalessandro et al. (2022), the combination of high spatial resolution infrared (IR) and optical photometric observations is a successful choice for analyzing the SFH of heavily extincted stellar clusters like Liller 1 (Pallanca et al. 2021). Consequently, the data set employed in this study – focused on Liller 1’s sibling, Terzan 5 – is a combination of archival images covering both optical and IR bands.

Specifically, the IR images were acquired with the Multiconjugate Adaptive Optics Demonstrator (MAD, science demonstration proposal, PI: Ferraro), which had been installed for a brief trial period at the ESO/VLT. This data set is composed of 15 exposures in the K filter1 roughly centered on Terzan 5, with a dither pattern resulting in a final field of view (FOV) of 1′ × 1′ (see the black square in Figure 1). The IR images used on this work are a subsample selected on the basis of the seeing quality, i.e. with a stable FWHM in the K-band of 0.1″ across the entire MAD FOV. The integration time of each exposure is texp = 120 s. The optical imagery consists of 10 exposures in the filters V and I acquired with the ACS/WFC on board the Hubble Space Telescope (HST), secured through proposal GO 12933 (PI: Ferraro). All the exposures are texp = 365 s long.

The aforementioned images have been already utilized in previous photometric studies of Terzan 5 (e.g., Massari et al. 2015; Ferraro et al. 2016). However, the reconstruction of the SFH from synthetic CMD fitting methods requires not only large samples of resolved stars, but also high photometric accuracy, the proper assessment of photometric completeness, and accurate corrections for differential reddening. Hence, following the same strategy adopted in the case if Liller 1 (see Ferraro et al. 2021; Pallanca et al. 2021; Dalessandro et al. 2022), we reanalyzed those images exploiting simultaneously the information coming from all the available exposures. This approach has been shown to lead to an increase in the statistics, photometric accuracy and the depth. Moreover, as discussed in Section 2.4, it allowed us to precisely reconstruct all the steps of the photometric reduction, which is necessary to estimate the photometric incompleteness and uncertainties. In tunr, this also brought to an improved estimate of the differential reddening affecting the system (see Sect. 2.3).

The photometric reduction was performed through PSF fitting by using the software DAOPHOTIV (Stetson 1987). First, PSF models were determined via the PSF routine for each exposure, selecting ~200 isolated, bright and not saturated stars homogeneously distributed across the entire FOVs. The analytical PSF functions were complemented by a spatially variable lookup table, following a cubic polynomial spatial variation. We applied the best PSF models to all the sources corresponding to flux peaks at least 3σ above the local background by using the ALLSTAR routine.

Differently from previous photometric analyses of Terzan 5, these first catalogs of instrumental magnitudes and positions were combined together via the cross-matching subroutines DAOMATCH and DAOMASTER (Stetson 1993), generating an input master starlist for ALLFRAME (Stetson 1994) which includes the information coming from both optical and IR images. Specifically, the master list was generated by including the sources detected in more than five images in the K filter, and in more than three images in at least one of the optical ones. The strength of this analysis is that stars detected only in the optical images were then searched and analyzed also in the K images (and vice versa), because the ALLFRAME routine forces a fit at the corresponding positions of stars in the master list. The information provided by the two filters separately is thus maximized, with fruitful results in terms of photometric accuracy and statistics.

The final raw magnitude for each recovered star (in a specific filter) is the weighted mean of multiple magnitude estimates (one for each exposure) once homogenized via DAOMATCH and DAOMASTER, while the photometric error is the corresponding standard deviation. The final catalog includes 30834 objects that have both I and K magnitudes. Therefore, it is limited to the overlapping region between the ACS/WFC and the MAD FOVs (see Figure 1).

We calibrated the optical instrumental magnitudes into the HST VEGAMAG system by applying the recipes and timedependent ACS zero points provided by the STScI web interface2. Aperture corrections were determined with respect to a radius of 8 pixel, selecting a subsample of bright and isolated stars across the FOV, and referring to the value of Encircle Energy tabulated in the HST website. The IR raw magnitudes were calibrated onto the 2MASS photometric system by computing the offset of the common stars between our catalog and the public available VISTA Variables in the Vía Láctea (VVV) survey (Minniti et al. 2010; Surot et al. 2019). Among the stars in common, we retained from the VVV catalog the ones with low photometric error (σK < 0.075), and we adopted as zero-point the median of the difference after a 3σ rejection. The resulting CMDs, both in optical and hybrid (I, IK) set of filters, are shown in Figure 2. The distributions of DAOPHOT photometric parameters χ and sharpness, which are related to the quality of the PSF fitting, are displayed in Figure 3, and they will be used to clean the CMDs from spurious objects or badly fitted stars.

Since the final catalog was built considering one optical exposure as reference frame, we corrected the instrumental coordinates for the ACS/WFC geometric distortions following the prescription of Anderson & Bedin (2010), and subsequently we cross-matched our catalog with the one from Ferraro et al. (2016) to retrieve HST-based high-resolution PMs calculated by Massari et al. (2015) and absolute coordinates (α, δ) based on the astrometric 2MASS catalog.

thumbnail Fig. 1

HST ACS/WFC image of Terzan 5 in the F814W filter. The field of view is 204″ × 204″. The black box corresponds to the VLT-MAD pointing (60″ × 60″). The region used for the SFH reconstruction is beyond the red circle, corresponding to a radius of 15″. North is up, east is to the left.

2.2 Field contamination

The next step is to distinguish genuine member stars of Terzan 5 from Galactic field interlopers. To this purpose, we used the relative PMs determined by Massari et al. (2015). The result of this procedure is shown in Figure 4. The leftmost panel shows the CMD of all the stars observed within the MAD FOV. Their measured PMs are plotted in the vector point diagrams (VPDs) shown in the rightmost panel of the figure. Following Massari et al. (2015), we flagged as members of Terzan 5 all the stars included within a radius of 1.5 mas yr−1 from the origin of the VPD (black dots), while field interlopers are those beyond this radius (red dots). Their respective CMDs are shown in the third and second panels of Figure 4, respectively.

Since the PM distribution of bulge field stars substantially overlap with the population of Terzan 5 in the VPD (see the rightmost panel of Fig. 4), the adopted criterion for the selection of Terzan 5 members (as those with PM < 1.5 mas yr−1) necessarily implies some residual field contamination. To ensure that the recovered SFH will not be affected by bulge or disk interlopers, we thus estimated the fraction of residual field interlopers by using two different methods. Firstly, we used a Monte Carlo approach: we fitted the two PM components (ΔPMx, ΔPMy) observed in the VPD (black and red points in the rightmost panel of Figure 4) with Gaussian functions, and subsequently we run 1000 realizations of such Gaussians, each time counting the number field stars falling within the 1.5 mas yr−1 radial limit. The contamination fraction is computed as the mean over the 1000 realizations. Secondly, we utilized the Python package scikit-learn (Pedregosa et al. 2011) to exploit the Gaussian mixture model (GMM) statistics. We forced the algorithm to fit two Gaussian components to the distribution of points observed in the VPDs and we looked at the weight of the field stars component. In both cases, we obtained a contamination fraction lower than 5%, which is not expected to heavily influence the SFH analysis. The low fraction of field objects is somehow expected, since we are considering a small field of view centered on Terzan 5.

thumbnail Fig. 2

Optical (I, V–I) and hybrid (I, I–K) CMDs of Terzan 5 obtained through the combined HST and MAD observations, corresponding to the FOV delimited by the black box in Figure 1. No quality cuts and differential reddening corrections are applied to these diagrams. It is worth noticing the distortion of the evolutionary sequences along the direction of the reddening vectors, indicated by the black arrows in the upper-left part of the panels. The red line corresponds to a 12 Gyr old isochrone (Bressan et al. 2012; Marigo et al. 2017) with [M/H] = −0.2 dex, E(BV) = 2.38, and μ0 = 13.87 (reference parameters from the database of Valenti et al. 2007), reddenned by adopting RV = 3.1. The dashed black lines show the same isochrone shifted by assuming RV = 2.1, 2.6, 3.6, 4.1 (from left to right). We varied E(BV) and μ0 to fit the optical (I, V–I) CMD. However, despite the remarkable broadening of the evolutionary sequences, these alternative RV values are clearly unable to match the hybrid CMD.

2.3 High spatial resolution differential reddening map

The Galactic bulge is one of the most extincted region of the Milky Way. The large amount of interstellar dust causes the scattering and the absorption of the (blue) light coming from the stars, which appear redder and fainter with respect to their intrinsic colors and magnitudes, depending on the degree of heterogeneity of the medium across the observed FOV. This effect is clearly visible in Figure 2, where all the evolutionary sequences of the CMD result blurred and elongated along the direction of the reddening vector (see the black arrows in Figure 2). Given the high extinction in the direction of Terzan 5 and its location within the Galaxy (E(BV) = 2.38; l, b = 3.81°, 1.67°; Valenti et al. 2007), this problem is exacerbated, and the correction for differential reddening effects becomes a crucial step in the identification and characterization of its subpopulations. Indeed, the study of Ferraro et al. (2016) was based on the differential extinction map presented in Massari et al. (2012), which reached a spatial resolution of 8″ × 8″. This resolution provided an overall good correction for differential reddening effects, but it led the detection of a double MS-TO only after the selection of the least extincted region of the MAD FOV, that corresponds to an area of ~25″ × 25″ in the south-west direction. To improve the accuracy of the correction, we decided to take advantage of the newly performed photometric analysis for building a higher resolution differential reddening map in the direction of Terzan 5.

It has been extensively shown that, when constructing differential reddening maps of stellar clusters, techniques based on the so-called star-by-star method can provide maps with a spatial scale of few arcseconds. Thanks to the improved statistics and photometric accuracy, we could model the differential reddening in the field of view covered by our observations with this novel approach, following the guidelines provided in several recent papers (e.g., Dalessandro et al. 2018; Saracino et al. 2019; Pallanca et al. 2021; Cadelano et al. 2020 and references therein). The star-by-star method estimates the differential reddening of each star in the catalog by computing the shift, along the reddening vector, needed to match the mean ridge line (MRL) of the cluster to the “local” CMD. The latter is defined as the CMD built with the objects spatially close to the targeted star (see below).

The direction and the components of the reddening vector3 are fixed once a precise extinction law (i.e., the dependence of the absorption coefficient A on wavelength λ) is adopted. Indeed, extinction laws are parameterized as (see e.g. Cardelli et al. 1989) Aλ=RV×cλ,RV×E(BV),$\[A_{\lambda}=R_{V} \times c_{\lambda, R_{V}} \times E(B-V),\]$(1)

where E(BV) is the color excess, that is defined as the difference between the observed (BV) and the intrinsic (BV)0 color. The coefficient cλ,RV represents the ratio Aλ/AV, and it depends on the assumed extinction law (i.e., Cardelli et al. 1989, O’Donnell 1994 or Fitzpatrick 1999), the adopted value of the RV parameter, and the effective wavelength of the filter. It is conventionally set equal to 1 in the (Johnson) V band. RV is usually set RV = 3.1, following fundamental studies about diffuse Galactic interstellar medium (Sneden et al. 1978). However, it has long been known that the coefficient RV can vary depending on the observed region of the Milky Way, particularly in the direction of the bulge (see, e.g., Johnson & Borgman 1963; Johnson 1965; Frogel et al. 1995; Popowski 2000; Alonso-García et al. 2011; Nataf et al. 2013; Casagrande & VandenBerg 2014; Alonso-García et al. 2017; Saha et al. 2019; Pallanca et al. 2021; Legnardi et al. 2023). A powerful tool to determine RV in the direction of a stellar system, even in the case of strong stretching of the evolutionary sequences, is to search for the value that allows the simultaneous match of the optical, NIR and hybrid CMDs with an isochrone of appropriate metallicity and age (see, e.g., the case of Liller 1 in Pallanca et al. 2021). In fact, changing the value of RV modifies the value of cλ,RV, making the extinction law more or less steeply dependent on wavelength, especially in the optical regime. Hence, following Pallanca et al. (2021), we computed an isochrone (Bressan et al. 2012; Marigo et al. 2017) with age t = 12 Gyr (Ferraro et al. 2016) and global metallicity [M/H] = −0.2 dex (Valenti et al. 2007, Massari et al. 2014), and we searched for possible combinations of RV, color excess E(BV), and distance modulus μ0 able to make the isochrone simultaneously matching the observed optical and in the hybrid CMDs. We varied RV by steps of 0.5 around the canonical value (RV = 3.1), since Legnardi et al. (2023) showed that, for heavily extincted stellar clusters, only large differences (e.g., ΔRV ~ 0.5) can affect the local reddening correction by using the star-by-star method. As shown in Figure 2, the best solution is obtained by adopting the standard value RV = 3.1 (red solid line), together with the color excess and distance estimated by Valenti et al. (2007) for Terzan 5, namely E(BV) = 2.38 and μ0 = 13.87. Therefore, we adopted the extinction coefficients calculated considering RV = 3.1, the extinction law from O’Donnell (1994), and effective wavelengths of the I (F814W) and K filters tabulated in Rodrigo & Solano (2020)4. Under these assumptions, the extinction coefficients are RI = 1.90 and RK = 0.365.

We then proceeded to estimate the reddening of each star of the catalog, following the procedure extensively explained in Pallanca et al. (2021). The corrections were computed in the PM-cleaned (see Section 2.2 for details) and hybrid (I, I–K) CMD. Briefly, we first determined the MRL that approximates the upper main sequence, subgiant branch and red giant branch of Terzan 5. To do so, we divided the CMD in different magnitude bins adapting the bin size to the evolutionary sequence (from 0.2 to 1 mag). We calculated the mean values of color and magnitude of each bin by considering only a subsample of member high-quality stars, chosen through a 3σ rejection applied to the distributions of χ and sharpness (see colored cells in Figure 3). The next step is to determine the value of differential reddening δE(BV) for all the stars. Specifically, for each star in the catalog, we first built the local CMD by using the N = 30 high-quality closest targets. Subsequently, the MRL is shifted along the direction of the reddening vector in steps of δE(BV), until it fits the local CMD. The value of reddening that minimizes the distance between the MRL and the local CMD is found by computing the residual color ΔIK, defined as: ΔIK=i=1N|IKobs,iIKMRL,i|+wi|IKobs,iIKMRL,i|,$\[\Delta \mathrm{IK}=\sum_{i=1}^{N_{\star}}\left|\mathrm{I} \mathrm{K}_{\mathrm{obs}, i}-\mathrm{IK}_{\mathrm{MRL}, i}\right|+w_{i} \cdot\left|\mathrm{IK}_{\mathrm{obs}, i}-\mathrm{IK}_{\mathrm{MRL}, i}\right|,\]$(2)

where wi is a weight assigned to the stars of the local CMD, in order to valorize stars with small photometric errors or distances from the targeted star. It is defined as wi=1diσi[j=1N(1djσj)]1.$\[w_{i}=\frac{1}{d_{i} \cdot \sigma_{i}}\left[\sum_{j=1}^{N_{\star}}\left(\frac{1}{d_{j} \cdot \sigma_{j}}\right)\right]^{-1}.\]$(3)

In Equations (2) and (3), IKobs,i is the observed color of each star composing the local CMD, IKMRL,i is the color of the MRL at the same magnitude level, d is the projected distance from the selected star, and σ the photometric uncertainty on the color. Finally, the assumed value of δE(BV) is the one that minimizes the normalized residual color (ΔIK/N).

The resulting absolute differential reddening map in the direction of Terzan 5 (assuming the average color excess estimated in Valenti et al. 2007) is shown in Figure 5. This map has a typical spatial resolution of 1″ − 2″ across the observed FOV, except for small regions limited to its very edge where the resolution is equal to ~3″. It is worth noticing that, despite the small area covered by the IR camera (60″ × 60″), the maximum variation of E(BV) across the FOV is ~1 mag, demonstrating that it is necessary to properly model the patchy absorption pattern of the dust clouds on small (arcsecond) scales. Indeed, the differential extinction map from Massari et al. (2012), built with a larger cell size (8″ × 8″), quoted an extinction variation of δE(BV) = 0.67 mag. The benefit from the differential reddening correction can be appreciated from the comparison between the observed CMD (Fig. 2) and the differential reddening corrected (DRC) CMD shown in the left-hand panel of Figure 6: indeed, independently of a slightly smaller color span along the x-axis, all the evolutionary sequences appear substantially less distorted along the reddening vector direction in Fig. 6.

thumbnail Fig. 3

Trends of the photometric quality parameters χ and sharpness as a function of the I magnitude resulting from the combined optical-IR reduction. Darker colors are associated to denser populated cells. Red points are stars excluded by the 3σ clipping procedure.

thumbnail Fig. 4

Cleaning procedure from field interlopers via PMs. Left panel: CMD of the stars included in the MAD footprint with measured PMs from Massari et al. (2015). Middle panel: CMD of the field population as selected from the VPDs, that is stars with PMs not compatible with that of Terzan 5. Right panel: PM-cleaned CMD of Terzan 5 obtained by considering only likely cluster members selected from the VPDs. Rightmost column: VPDs of the measured stars divided in 5 bins of 1.5 mag each, starting at IDRC = 16.5. Black dots represent stars within a radius of 1.5 mas yr−1 and correspond to likely cluster members according to the criterion of Massari et al. (2015).

2.4 Artificial stars experiments

A new photometric analysis has been necessary also for the determination of the incompleteness and photometric errors of the input photometric catalog through artificial star tests. In fact, this method requires to re-perform the same photometric analysis that produced the final catalog once artificial stars are placed onto the images.

Following the same procedure adopted in the case of Liller 1 (Dalessandro et al. 2022), we created a catalog of artificial stars extracting input magnitudes in the I band (Iin) from a flat luminosity function (LF). The extremes of the LF were chosen to match the observed CMD. We added the input artificial stars on the images with the PSF models and Poisson noise resulting from the photometric analysis thanks to the DAOPHOT ADDSTAR package. The stars were placed on the ACS/WFC images following the stellar density profile modeled by Lanzoni et al. (2010), which corresponds to a King model (King 1966) with structural parameters c = 1.49, rcore = 9″, and rt = 4.6′. To avoid overcrowding, we divided the images into cells of 20 × 20 pixel2 (more than 10 times the FWHM measured in the ACS/WFC), where only one artificial star was injected. We assigned to each artificial star a Kin and a Vin magnitude from a randomly extracted color in the range of the observed (IK) and (VI) colors. In this way, we could homogeneously cover the colormagnitude plane of the observed stars. Thanks to the magnitudes homogenization and the geometrical transformation between the different detectors and exposures calculated by DAOMATCH and DAOMASTER, we ensured to place artificial stars in the same location in the V, I and K images, assigning the proper framedependent magnitude. We then performed the same photometric analysis described in Section 2.1 on the images containing artificial stars, and, at the end of the reduction, we cross-matched the output catalog with the one containing only input artificial stars. We iterated this procedure until the final catalog of artificial stars counted ~700 000 objects within the common FOV between the optical and IR data sets (see Figure 1).

In this way, we could retrieve two important measurements: on the one hand, it is possible to count the fraction of stars lost during the photometric procedure, obtaining an estimate of the completeness of the photometric catalog as a function of color and magnitude; on the other hand, for the recovered stars, the difference between the input and the output magnitude provides a more realistic estimate of the photometric errors. Specifically, the completeness was defined as the ratio between the number of recovered stars and that of the injected ones. An artificial star was considered “lost” when it is not recovered in the I or K filter, or when it is recovered more than 0.75 magnitude brighter. Indeed, in this latter case, the artificial star is likely blended with one or more other real stars, creating an “altered” source with more than twice the flux of the input one. This completeness is a sensitive function of the crowding conditions, as shown by the curves displayed in the upper panels of Figure 7, built by dividing the FOV in annuli centered on the cluster center (from Lanzoni et al. 2010). As shown in Figure 7, the completeness significantly drops in the first 15″ from the cluster center. Within this distance bin, the magnitude corresponding to the 50% completeness is I = 21.3 (or K = 16.8), which is ~0.3 mag brighter than the MSTO of the old population (assuming literature parameters, see the red isochrone in Figure 2). On the other hand, the photometric completeness is equal to or better than 80% at the old MS-TO level at larger radii. We concluded that the severe crowding of the innermost 15″ of Terzan 5 irremediably hampers the use of the MS-TO for the determination of its SFH, and therefore we restricted the analysis to larger radial distances. The radial limit of 15″ is marked by a dashed red circle in Figure 1. Considering the region beyond this radial cut in the field of view covered by the observations, we still sample the ~45% of the total mass of Terzan 5, computed from the best-fit King model of Lanzoni et al. (2010).

Photometric errors and blends are assigned to the synthetic stars by using the distributions of (IoutIin) and (KoutKin) displayed in the bottom panels of Figure 7. As expected, the distributions become broader and more asymmetric toward fainter magnitudes, indicating larger photometric uncertainties and blend effects.

thumbnail Fig. 5

Map of the absolute reddening E(BV) in the direction of Terzan 5. Darker colors indicate the more extincted areas within the observed FOV. The white cross marks the gravity center of Terzan 5 (Lanzoni et al. 2010). North is up, while east is on the left.

thumbnail Fig. 6

DRC and PM-cleaned CMDs of Terzan 5 obtained through the photometric reduction and the differential reddening map described in this work. The left-hand panel displays all the member stars measured in the FOV, while the right-hand panel shows the subsample of stars satisfying our quality cuts, specifically a radial distance of r ≥ 15″ from the center and good values of χ and asharpness (i.e., within 3σ of their distributions; see Fig. 3).

thumbnail Fig. 7

Upper panels: completeness curves as a function of the I and K magnitudes in different concentric annuli from the cluster’s center. The colored dots and vertical lines highlight the limit magnitude at which the 50% is reached within the respective radial annulus. Lower panels: photometric errors, i.e. the difference between the input and output magnitudes, in I and K bands recovered from the artificial star tests.

3 Reconstruction of the star formation history

3.1 Brief description of the method

As for the case of Liller 1 (Dalessandro et al. 2022), we reconstructed the SFH exploiting the synthetic CMD fitting technique (e.g., Tosi et al. 1991; Tolstoy & Saha 1996; Dolphin 1997, 2002; Hernandez et al. 1999; Cignoni & Tosi 2010; Aparicio & Hidalgo 2009; Weisz et al. 2012; Ruiz-Lara et al. 2018) through the code SFERA (Star Formation Evolution Recovery Algorithm, Cignoni et al. 2015).

SFERA parameterizes the SFH as a linear combination of j × k star formation events (also called basis functions, BFs), each characterized by a finite duration Δt (centered at the time step tj) and a narrow interval of metallicities Δz (centered at the metallicity step zk). For each synthetic star, we randomly extract age and metallicity from the intervals (tj,tj + Δt) and (zk, zk + Δz), and a stellar mass from the adopted initial mass function (IMF). Magnitudes are assigned by interpolating a library of isochrones. A percentage of synthetic stars is then coupled with a companion star, whose mass is extracted from the same IMF and whose flux is added to the primary one. To overcome the short duration of some evolutionary phases, the BFs are typically populated with a high number of synthetic stars (of the order of 106).

To produce realistic simulations, theoretical photometry is degraded to mimic the observational conditions of the data. This translates into including extinction and distance effects, and convolving the synthetic photometry with incompleteness and photometric errors. These observational uncertainties are quantified by means of extensive artificial stars experiments (see Section 2.4 for details). Finally, BFs’ CMDs and the observational CMD are binned in n bins of color and m bins of magnitude to construct density maps (Hess diagrams). The result is a library of j × k Hess diagrams BFm,n(j, k), whose elements are linearly combined to match the observational counterpart, following the equation below: Nn,m=jka(j,k)×BFm,n(j,k).$\[N_{n, m}=\sum_{j} \sum_{k} a(j, k) \times \mathrm{BF}_{m, n}(j, k).\]$(4)

The coefficients a(j, k), representing the SFR at the time step j and metallicity step k, are found by minimizing a Poissonian likelihood: χP2=2i=1Nbineioi+oi ln oiei,$\[\chi_{P}^{2}=2 \sum_{i=1}^{N_{b i n}} e_{i}-o_{i}+o_{i} ~\ln~ \frac{o_{i}}{e i},\]$(5)

where ei and oi are, respectively, the expected and the measured number of stars contained in the ith cell of the binned CMDs. The uncertainties around the best-fit solution are estimated by bootstrapping the data and re-deriving the solutions.

From the minimization process we obtain the set of weights a(j, k) from which, together with the total mass M(j, k) and duration of each BF, we find the most likely SFH for the data set under analysis. Moreover, the mean metallicity of the “best” model in age step tj can be obtained by adding up all the possible metallicity values covered by the BF zk, weighted by the ratio between their masses and the total mass of the “best” model in the corresponding age range: z(tj)=kzka(j,k)M(j,k)ka(j,k)M(j,k).$\[z\left(t_{j}\right)=\sum_{k} z_{k} \frac{a(j, k) M(j, k)}{\sum_{k} a(j, k) M(j, k)}.\]$(6)

The procedure can be repeated over all the chosen age steps, in order to obtain an estimate of the “true” age-metallicity relation.

To explore the wide parameter space resulting from the choice of the age and metallicity steps, SFERA combines the genetic algorithm (GA) Pikaia6 with a local search routine. On the one hand, the GA allows to explore the parameter space in more points simultaneously, being less sensitive to the initial conditions and local minima; on the other hand, the local search speeds up the convergence time and increases the accuracy of the solution (see Cignoni et al. 2015 for an extensive and detailed description).

3.2 Input parameters

The BF CMDs used for the comparison with the observed one were built from the PARSEC-COLIBRI stellar isochrones (Bressan et al. 2012; Marigo et al. 2017), covering an age range from 2 Myr to 13.4 Gyr, with a logarithmic step of Δ log(t/yr) = 0.005 dex, and a metallicity range from [M/H] = −1.5 to [M/H] = +0.3 dex7, with a step of 0.025 dex. The isochrones were computed in the HST/ACS-WFC and 2MASS photometric systems, since these are the photometric bands used to build the observed CMD.

The SFH is parameterized by 14 contiguous time steps of variable duration: 0.5 Gyr for stars younger than 1 Gyr, 1 Gyr for stars between 1 and 13 Gyr, 1.4 Gyr for older star. The metallicity step is fixed at 0.1 dex. The BF CMDs were generated with a Kroupa (2001) IMF between 0.1 and 300 M. 30% of synthetic stars are coupled with a stellar companion sampled from the same IMF to mimic unresolved binaries.

The distance modulus (μ0) and average reddening E(BV) are initially adopted from the literature (Valenti et al. 2007), then the code is permitted to adjust these parameters to maximize the likelihood. In order to model the residual differential reddening, we added to each synthetic star a Gaussian reddening whose width is selected during the fitting process.

The CMD bin size is chosen as a compromise between the photometric error, resolution of the selected stellar models, and computational time, ranging from 0.1 to 0.05 mag both in magnitude and color. Finally, the routine does not assume any age-metallicity relation a priori.

3.3 SFERA best-fit solution

The result of the differential reddening correction and PM selection procedures discussed in Sects. 2.3 and 2.2 is shown in the left panel of Figure 6. By taking into account the completeness limits and to work with a suitable photometric sample, we run SFERA considering the CMD shown in the right panel of the figure, which includes only stars lying at distances larger than 15″ from the center with optimal values of χ and sharpness (i.e., stars not excluded by the 3σ clipping procedure, see Figure 3). We also limited the analysis to stars brighter than I = 22, which corresponds to a completeness larger than 50% in the considered area. Following this selection, the final sample of stars amounts to 8104.

SFERA found a best-fit extinction corrected distance modulus μ0 = 13.75 and average foreground reddening E(BV) = 2.5. The distance modulus agrees, within the uncertainties, with the value quoted in Valenti et al. (2007), while the color excess is ~0.1 mag higher, being only marginally in agreement with the E(BV) calculated in the same work8. However, it is worth recalling the remarkable variation of E(BV) across the FOV of the observations (see Section 2.3, and also Massari et al. 2012).

The red line in Figure 8 shows the best-fit SFH. The shaded areas represent the statistical uncertainties, corresponding to the 5th and 95th percentiles of the distribution of all the synthetic CMDs produced by SFERA. The SFH of Terzan 5 exhibits a main burst of star formation (SF) at t ⩾ 12 Gyr, that confirms the age of the old population inferred by the previous isochrone fitting analysis (Ferraro et al. 2016). Then, for younger ages, we can appreciate a prolonged star formation composed of several minor peaks, with the most relevant ones occurring around 8–7 Gyr, and between 4–2 Gyr ago.

The upper plot of Figure 9 displays the observed data (lefthand panel) and the Hess diagram for the synthetic CMD, corresponding to the best-fit SFH (right-hand panel). From a visual inspection, the best-fit solution nicely reproduces both the extension and the color distribution of the blue plume, the location of the old MS-TO, and the complexity of the RC region. To compare the observed and synthetic CMDs in a more quantitative way, in the lower panel of Figure 9 we display the observed (gray histogram) and best-fit solution (red lines) color distributions in four different magnitude bins. Overall, the model color function does a fairly good job in reproducing both the width and the mean color distributions in all magnitude bins.

The cumulative stellar mass function (CMF) associated to this solution is reported in Figure 10. It shows that the first burst of SF builts ~65% of the observed stellar mass. Afterward, the CMF of Terzan 5 gradually grows reaching the ~90% at t < 8 Gyr. The missing 10% was then built-up during the latest SF episode (2–4 Gyr ago). The fraction of mass built during the young SF episode barely reconciles with the estimate from Lanzoni et al. (2010), who concluded, based on the RC star counts, that the young population accounts for ~38% of the total mass. However, it is important to note that, due to observational limitations, we excluded from the analysis the innermost 15″ from the cluster center, where the metal-richer and younger population is preferentially concentrated (Ferraro et al. 2009; Lanzoni et al. 2010). Hence, the fraction of mass ascribable to the young stellar populations is inevitably underestimated in this analysis.

Finally, in Figure 11 we show the metallicity distribution function (MDF) recovered from the best-fit solution. The distribution turns out to be clearly bimodal with two major components: a subsolar peak at −0.7 < [M/H] < −0.5 dex, and a supersolar one at [M/H] ~ +0.25 dex. Qualitatively, the two metallicity peaks correspond to the oldest and youngest star formation episodes, but a one-to-one association is not expected because intermediate bursts may introduce some spreads, instead of further isolated peaks. The derived distribution appears to nicely reproduce the overall shape of the one obtained by Origlia et al. (2013); Massari et al. (2014) by using high-resolution spectroscopy of about 200 stars, including also the minor metal poor sub-component discussed in Origlia et al. (2013)9.

thumbnail Fig. 8

SFH of Terzan 5 associated to the SFERA best-fit solution. The red shaded areas are the uncertainties associated to the model computed from a bootstrap method, and they mark the 5th and 95th confidence level. To obtain this result, an average foreground extinction of E(BV) = 2.5 and distance modulus μ0 = 13.75 were adopted.

thumbnail Fig. 9

Upper panels: Hess diagrams of the observed DRC (left) and best-fit synthetic (right) CMD. Lower panels: color functions of the CMDs in four different magnitude bins; the data and relative Poissonian uncertainty are displayed in gray, while the models are the superimposed red line.

thumbnail Fig. 10

Cumulative stellar mass function for the recovered best-fit SFH.

thumbnail Fig. 11

Metallicity distribution function as derived from the best-fit solution.

4 Summary and conclusions

The photometric and spectroscopic observations recently gathered for two massive and metal-rich stellar clusters in the Galactic bulge - Terzan 5 and Liller 1 - disclosed their intriguing peculiarity of hosting multiple stellar populations with remarkable metallicity (Δ[Fe/H] > 0.5 dex) and age (Δt > 4 Gyr) spreads (Ferraro et al. 2009, 2016, 2021; Origlia et al. 2011, 2013; Massari et al. 2014; Dalessandro et al. 2022; Crociati et al. 2023). At present, these are the only in-situ (see Massari et al. 2019; Callingham et al. 2022) stellar systems showing an old and subsolar population coexisting with a young and supersolar one, and their origin is still debated.

The first hypothesis links these systems to one of the mechanisms that possibly contributed to the formation of the Milky Way’s bulge (Ferraro et al. 2009, 2021). Specifically, massive and highly star-forming complexes, customarily dubbed clumps, are commonly observed in high-redshift, still bulgeless disk galaxies (z ≥ 1.5; see e.g., Guo et al. 2015). An increasing number of observational (Guo et al. 2018; Ambachew et al. 2022) and theoretical (e.g., Elmegreen et al. 2008; Bournaud 2016; Garver et al. 2023) studies are suggesting that these disk structures can migrate to the center of the host galaxy and contribute to the assembly of its bulge. A small fraction of them is also predicted to survive this coalescence process and be nowadays observable as cluster-like stellar systems (Dekel et al. 2023). Therefore, Terzan 5 and Liller 1 can represent the present-day remnants of two such massive clumps. Given their high progenitor mass, these systems may have undergone a self-enrichment process. An alternative hypothesis, firstly advanced by McKenzie & Bekki (2018, see also Bastian & Pfeffer 2022) through numerical simulations, involves the gravitational capture of field GMCs by old massive GCs. Under highly fine-tuned conditions (regarding the orbits, relative velocities, masses and dimensions of the systems, the number of GMCs supposed to be present in the bulge at the time of the merging, etc.), a strongly bound GC-GMC collision would fuel subsequent a burst of star formation. However, this is an intrinsically rare event, which very unlikely can happen more than once in a cluster life.

An efficient way to disentangle between these two formation scenarios is the study of the SFH of these systems, to verify whether they experienced just one or more star formation bursts, or even a prolonged activity. Indeed, the reconstructed SFH of Liller 1 (Dalessandro et al. 2022) shows evidence of multiple and long-lasting star formation bursts. Motivated by those results, here we conducted an analogous investigation on Terzan 5. To retrieve its SFH, we employed the code SFERA comparing the observed CMD of Terzan 5 with a set of synthetic ones through the CMD fitting technique. The observed CMD used in this comparison has taken advantage of the powerful combination of deep optical (HST) and IR Adaptive Optics (VLT/MAD) observations. Indeed, the hybrid (I, I–K) CMD turns out to be the ideal tool for obtaining a clear view of the stellar evolutionary sequences, especially in regions of strong reddening as the Galactic bulge. In addition, such a data set allowed the accurate determination of relative PMs for the removal of Galactic field interlopers, and the construction of a high spatial resolution (below 1.5″) reddening map, used to correct for the effects of differential extinction.

The best-fit solution shows a main and narrow burst of SF around 13 Gyr ago, combined with a decreasing and prolonged SF activity followed by another significant SF episode starting 4 Gyr ago and extending for ~2 Gyr. A possible additional burst occurred ~7–8 Gyr ago. The overall time extension of the main bursts is remarkable, and it would point to a constant supply of gas to fuel a SF activity with a duration of more than 1 Gyr. Beside the details, the cumulative mass function derived from this study clearly shows a gradual build-up of the stellar mass in this system, which suggests a continuous, low rate SF activity characterized by a few main bursts.

The overall structure of the reconstructed SFH of Terzan 5 appears astonishingly similar to that derived for Liller 1 (Dalessandro et al. 2022), showing an uninterrupted star formation activity with at least three relevant star formation events in both systems. The first, major episode appears to be quite prolonged in time (1–2 Gyr), at odds with what observed for genuine GCs. This is followed by two additional bursts emerging from an underlying continuous star formation activity with very low intensity. The most recent bust occurs at different lookback times in the two systems (1–2 Gyr ago in Liller 1, and 4–5 Gyr ago in Terzan 5). It has been possibly triggered by violent interactions with other Bulge substructures (as the bar) and thus likely depends on the details of the orbits of each stellar system. Intriguingly, instead, the intermediate-age burst occurred at the same cosmic age in both the system (8–9 Gyr ago). This curious synchrony of the two events needs to be further investigated since it might trace a major event in the past of the Galactic Bulge. In both cases, the reconstructed SFH is hardly reconcilable with the predicted short timescales (<10 Myr, McKenzie & Bekki 2018) needed to a GMC to form stars after a bound collision with a GC. The results presented in this paper thus further support a scenario where the complexity of Terzan 5 and Liller 1 is due to self-enrichment processes in massive structures that contributed to form the bulge of our Galaxy.

Acknowledgements

This work is part of the project Cosmic-Lab (“Globular Clusters as Cosmic Laboratories”) at the Physics and Astronomy Department “A. Righi” of the Bologna University (http://www.cosmic-lab.eu/Cosmic-Lab/Home.html). C.C. acknowledges the ESO Studentship Programme for funding the period spent at the ESO Headquarters in Garching bei München (Germany) where part of this work was carried out. E.V. acknowledges the Excellence Cluster ORIGINS funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC209439078331.

References

  1. Alonso-García, J., Mateo, M., Sen, B., Banerjee, M., & von Braun, K. 2011, AJ, 141, 146 [CrossRef] [Google Scholar]
  2. Alonso-García, J., Minniti, D., Catelan, M., et al. 2017, ApJ, 849, L13 [Google Scholar]
  3. Alvarez Garay, D. A., Fanelli, C., Origlia, L., et al. 2024, A&A, 686, A198 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Ambachew, L., Fisher, D. B., Glazebrook, K., et al. 2022, MNRAS, 512, 3079 [NASA ADS] [CrossRef] [Google Scholar]
  5. Anderson, J., & Bedin, L. R. 2010, PASP, 122, 1035 [CrossRef] [Google Scholar]
  6. Aparicio, A., & Hidalgo, S. L. 2009, AJ, 138, 558 [NASA ADS] [CrossRef] [Google Scholar]
  7. Athanassoula, E. 2005, MNRAS, 358, 1477 [NASA ADS] [CrossRef] [Google Scholar]
  8. Barbary, K. 2016, https://doi.org/10.5281/zenodo. 804967 [Google Scholar]
  9. Bastian, N., & Pfeffer, J. 2022, MNRAS, 509, 614 [Google Scholar]
  10. Behrendt, M., Burkert, A., & Schartmann, M. 2016, ApJ, 819, L2 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bournaud, F. 2016, Astrophys. Space Sci. Lib., 418, 355 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bournaud, F., Elmegreen, B. G., & Elmegreen, D. M. 2007, ApJ, 670, 237 [Google Scholar]
  13. Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
  14. Cadelano, M., Saracino, S., Dalessandro, E., et al. 2020, ApJ, 895, 54 [NASA ADS] [CrossRef] [Google Scholar]
  15. Callingham, T. M., Cautun, M., Deason, A. J., et al. 2022, MNRAS, 513, 4107 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
  17. Carollo, C. M., Scarlata, C., Stiavelli, M., Wyse, R. F. G., & Mayer, L. 2007, ApJ, 658, 960 [NASA ADS] [CrossRef] [Google Scholar]
  18. Casagrande, L., & VandenBerg, D. A. 2014, MNRAS, 444, 392 [Google Scholar]
  19. Cava, A., Schaerer, D., Richard, J., et al. 2018, Nat. Astron., 2, 76 [Google Scholar]
  20. Cignoni, M., & Tosi, M. 2010, Adv. Astron., 2010, 158568 [CrossRef] [Google Scholar]
  21. Cignoni, M., Sabbi, E., van der Marel, R. P., et al. 2015, ApJ, 811, 76 [Google Scholar]
  22. Claeyssens, A., Adamo, A., Richard, J., et al. 2023, MNRAS, 520, 2180 [NASA ADS] [CrossRef] [Google Scholar]
  23. Clarke, A. J., Debattista, V. P., Nidever, D. L., et al. 2019, MNRAS, 484, 3476 [NASA ADS] [CrossRef] [Google Scholar]
  24. Clarkson, W., Sahu, K., Anderson, J., et al. 2008, ApJ, 684, 1110 [NASA ADS] [CrossRef] [Google Scholar]
  25. Clarkson, W. I., Sahu, K. C., Anderson, J., et al. 2011, ApJ, 735, 37 [NASA ADS] [CrossRef] [Google Scholar]
  26. Combes, F., Debbasch, F., Friedli, D., & Pfenniger, D. 1990, A&A, 233, 82 [NASA ADS] [Google Scholar]
  27. Crociati, C., Valenti, E., Ferraro, F. R., et al. 2023, ApJ, 951, 17 [NASA ADS] [CrossRef] [Google Scholar]
  28. Dalessandro, E., Lardo, C., Cadelano, M., et al. 2018, A&A, 618, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Dalessandro, E., Crociati, C., Cignoni, M., et al. 2022, ApJ, 940, 170 [NASA ADS] [CrossRef] [Google Scholar]
  30. Debattista, V. P., Ness, M., Gonzalez, O. A., et al. 2017, MNRAS, 469, 1587 [Google Scholar]
  31. Debattista, V. P., Liddicott, D. J., Gonzalez, O. A., et al. 2023, ApJ, 946, 118 [NASA ADS] [CrossRef] [Google Scholar]
  32. Dekel, A., Sari, R., & Ceverino, D. 2009, ApJ, 703, 785 [Google Scholar]
  33. Dekel, A., Tziperman, O., Sarkar, K. C., et al. 2023, MNRAS, 521, 4299 [NASA ADS] [CrossRef] [Google Scholar]
  34. Dessauges-Zavadsky, M., Schaerer, D., Cava, A., Mayer, L., & Tamburello, V. 2017, ApJ, 836, L22 [Google Scholar]
  35. Dolphin, A. 1997, New A, 2, 397 [NASA ADS] [CrossRef] [Google Scholar]
  36. Dolphin, A. E. 2002, MNRAS, 332, 91 [NASA ADS] [CrossRef] [Google Scholar]
  37. Elmegreen, B. G., Bournaud, F., & Elmegreen, D. M. 2008, ApJ, 688, 67 [NASA ADS] [CrossRef] [Google Scholar]
  38. Elmegreen, D. M., Elmegreen, B. G., Marcus, M. T., et al. 2009, ApJ, 701, 306 [CrossRef] [Google Scholar]
  39. Ferraro, F. R., Dalessandro, E., Mucciarelli, A., et al. 2009, Nature, 462, 483 [Google Scholar]
  40. Ferraro, F. R., Massari, D., Dalessandro, E., et al. 2016, ApJ, 828, 75 [NASA ADS] [CrossRef] [Google Scholar]
  41. Ferraro, F. R., Pallanca, C., Lanzoni, B., et al. 2021, Nat. Astron., 5, 311 [Google Scholar]
  42. Fitzpatrick, E. L. 1999, PASP, 111, 63 [Google Scholar]
  43. Fragkoudi, F., Di Matteo, P., Haywood, M., et al. 2017, A&A, 606, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Frogel, J. A., Kuchinski, L. E., & Tiede, G. P. 1995, AJ, 109, 1154 [NASA ADS] [CrossRef] [Google Scholar]
  45. Garver, B. R., Nidever, D. L., Debattista, V. P., Beraldo e Silva, L., & Khachaturyants, T. 2023, ApJ, 953, 128 [NASA ADS] [CrossRef] [Google Scholar]
  46. Guo, Y., Ferguson, H. C., Bell, E. F., et al. 2015, ApJ, 800, 39 [NASA ADS] [CrossRef] [Google Scholar]
  47. Guo, Y., Rafelski, M., Bell, E. F., et al. 2018, ApJ, 853, 108 [NASA ADS] [CrossRef] [Google Scholar]
  48. Hernandez, X., Valls-Gabaud, D., & Gilmore, G. 1999, MNRAS, 304, 705 [NASA ADS] [CrossRef] [Google Scholar]
  49. Immeli, A., Samland, M., Gerhard, O., & Westera, P. 2004a, A&A, 413, 547 [CrossRef] [EDP Sciences] [Google Scholar]
  50. Immeli, A., Samland, M., Westera, P., & Gerhard, O. 2004b, ApJ, 611, 20 [CrossRef] [Google Scholar]
  51. Johnson, H. L. 1965, ApJ, 141, 923 [NASA ADS] [CrossRef] [Google Scholar]
  52. Johnson, H. L., & Borgman, J. 1963, Bull. Astron. Inst. Netherlands, 17, 115 [NASA ADS] [Google Scholar]
  53. King, I. R. 1966, AJ, 71, 64 [Google Scholar]
  54. Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
  55. Kunder, A., Koch, A., Rich, R. M., et al. 2012, AJ, 143, 57 [NASA ADS] [CrossRef] [Google Scholar]
  56. Lanzoni, B., Ferraro, F. R., Dalessandro, E., et al. 2010, ApJ, 717, 653 [Google Scholar]
  57. Legnardi, M. V., Milone, A. P., Cordoni, G., et al. 2023, MNRAS, 522, 367 [CrossRef] [Google Scholar]
  58. Marigo, P., Girardi, L., Bressan, A., et al. 2017, ApJ, 835, 77 [Google Scholar]
  59. Massari, D., Mucciarelli, A., Dalessandro, E., et al. 2012, ApJ, 755, L32 [NASA ADS] [CrossRef] [Google Scholar]
  60. Massari, D., Mucciarelli, A., Ferraro, F. R., et al. 2014, ApJ, 795, 22 [NASA ADS] [CrossRef] [Google Scholar]
  61. Massari, D., Dalessandro, E., Ferraro, F. R., et al. 2015, ApJ, 810, 69 [CrossRef] [Google Scholar]
  62. Massari, D., Koppelman, H. H., & Helmi, A. 2019, A&A, 630, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Matteucci, F. 2021, A&A Rev., 29, 5 [NASA ADS] [CrossRef] [Google Scholar]
  64. McKenzie, M., & Bekki, K. 2018, MNRAS, 479, 3126 [CrossRef] [Google Scholar]
  65. Messa, M., Dessauges-Zavadsky, M., Richard, J., et al. 2022, MNRAS, 516, 2420 [Google Scholar]
  66. Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, New A, 15, 433 [Google Scholar]
  67. Nataf, D. M., Gould, A., Fouqué, P., et al. 2013, ApJ, 769, 88 [Google Scholar]
  68. Ness, M., Freeman, K., Athanassoula, E., et al. 2013, MNRAS, 430, 836 [NASA ADS] [CrossRef] [Google Scholar]
  69. Noguchi, M. 1999, ApJ, 514, 77 [NASA ADS] [CrossRef] [Google Scholar]
  70. O’Donnell, J. E. 1994, ApJ, 422, 158 [Google Scholar]
  71. Origlia, L., Rich, R. M., Ferraro, F. R., et al. 2011, ApJ, 726, L20 [CrossRef] [Google Scholar]
  72. Origlia, L., Massari, D., Rich, R. M., et al. 2013, ApJ, 779, L5 [NASA ADS] [CrossRef] [Google Scholar]
  73. Origlia, L., Mucciarelli, A., Fiorentino, G., et al. 2019, ApJ, 871, 114 [NASA ADS] [CrossRef] [Google Scholar]
  74. Pallanca, C., Ferraro, F. R., Lanzoni, B., et al. 2021, ApJ, 917, 92 [NASA ADS] [CrossRef] [Google Scholar]
  75. Patsis, P. A., Skokos, C., & Athanassoula, E. 2002, MNRAS, 337, 578 [NASA ADS] [CrossRef] [Google Scholar]
  76. Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, J. Mach. Learn. Res., 12, 2825 [Google Scholar]
  77. Popowski, P. 2000, ApJ, 528, L9 [NASA ADS] [CrossRef] [Google Scholar]
  78. Queiroz, A. B. A., Chiappini, C., Perez-Villegas, A., et al. 2021, A&A, 656, A156 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Renzini, A., Gennaro, M., Zoccali, M., et al. 2018, ApJ, 863, 16 [NASA ADS] [CrossRef] [Google Scholar]
  80. Rodrigo, C., & Solano, E. 2020, in XIV.0 Scientific Meeting (virtual) of the Spanish Astronomical Society, 182 [Google Scholar]
  81. Rojas-Arriagada, A., Recio-Blanco, A., de Laverny, P., et al. 2017, A&A, 601, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Rojas-Arriagada, A., Zasowski, G., Schultheis, M., et al. 2020, MNRAS, 499, 1037 [Google Scholar]
  83. Romano, D., Ferraro, F. R., Origlia, L., et al. 2023, ApJ, 951, 85 [NASA ADS] [CrossRef] [Google Scholar]
  84. Ruiz-Lara, T., Gallart, C., Beasley, M., et al. 2018, A&A, 617, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Saha, A., Vivas, A. K., Olszewski, E. W., et al. 2019, ApJ, 874, 30 [NASA ADS] [CrossRef] [Google Scholar]
  86. Saracino, S., Dalessandro, E., Ferraro, F. R., et al. 2015, ApJ, 806, 152 [NASA ADS] [CrossRef] [Google Scholar]
  87. Saracino, S., Dalessandro, E., Ferraro, F. R., et al. 2019, ApJ, 874, 86 [Google Scholar]
  88. Shen, J., Rich, R. M., Kormendy, J., et al. 2010, ApJ, 720, L72 [NASA ADS] [CrossRef] [Google Scholar]
  89. Shibuya, T., Ouchi, M., Kubo, M., & Harikane, Y. 2016, ApJ, 821, 72 [NASA ADS] [CrossRef] [Google Scholar]
  90. Sneden, C., Gehrz, R. D., Hackwell, J. A., York, D. G., & Snow, T. P. 1978, ApJ, 223, 168 [CrossRef] [Google Scholar]
  91. Stetson, P. B. 1987, PASP, 99, 191 [Google Scholar]
  92. Stetson, P. B. 1993, IAU Colloq., 136, 291 [NASA ADS] [Google Scholar]
  93. Stetson, P. B. 1994, PASP, 106, 250 [Google Scholar]
  94. Surot, F., Valenti, E., Hidalgo, S. L., et al. 2019, A&A, 623, A168 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Tolstoy, E., & Saha, A. 1996, ApJ, 462, 672 [Google Scholar]
  96. Tosi, M., Greggio, L., Marconi, G., & Focardi, P. 1991, AJ, 102, 951 [NASA ADS] [CrossRef] [Google Scholar]
  97. Valenti, E., Ferraro, F. R., & Origlia, L. 2007, AJ, 133, 1287 [Google Scholar]
  98. Valenti, E., Zoccali, M., Renzini, A., et al. 2013, A&A, 559, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Vanzella, E., Castellano, M., Bergamini, P., et al. 2022a, A&A, 659, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Vanzella, E., Castellano, M., Bergamini, P., et al. 2022b, ApJ, 940, L53 [NASA ADS] [CrossRef] [Google Scholar]
  101. Wegg, C., & Gerhard, O. 2013, MNRAS, 435, 1874 [Google Scholar]
  102. Wegg, C., Gerhard, O., & Portail, M. 2015, MNRAS, 450, 4050 [NASA ADS] [CrossRef] [Google Scholar]
  103. Weisz, D. R., Zucker, D. B., Dolphin, A. E., et al. 2012, ApJ, 748, 88 [CrossRef] [Google Scholar]
  104. Zanella, A., Daddi, E., Le Floc’h, E., et al. 2015, Nature, 521, 54 [Google Scholar]
  105. Zoccali, M., Gonzalez, O. A., Vasquez, S., et al. 2014, A&A, 562, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  106. Zoccali, M., Vasquez, S., Gonzalez, O. A., et al. 2017, A&A, 599, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

1

Throughout the manuscript, we refer to the filters Ks, F606W and F814W with the labels K, V and I, respectively.

3

Specifically for the hybrid CMD, the components of the reddening vector are Vx = AIAK; Vy = AI.

4

λe f f, F814W = 7973.39 Å ; λe f f,K = 21590 Å. We emphasize that also Massari et al. (2012) adopted RV = 3.1 in the Johnson V-band, while the value of 2.83 quoted in that paper refers to the F606W filter of the HST ACS/WFC.

5

We calculated the extinction coefficients by means of the Python package extinction (Barbary 2016) adopting the O’Donnell (1994) extinction law.

6

This is a public available routine developed at the High Altitude Observatory.

7

We assume [M/H] = log (Z/Z), with Z = 0.0152.

8

In this respect, it is important to specify that, in the case of Valenti et al. (2007), the mean color excess associated to Terzan 5 was determined from the pure NIR plane with no differential reddening corrections applied to the CMD.

9

Small differences in the absolute value of the metallicity of the subsolar components are possibly due to a variety of factors, such as the adoption of solar-scaled isochrones, systematic offsets between isochrones and observations, and a possible metallicity-reddening degeneracy, which can be non-negligible in the extreme reddening conditions of Terzan 5.

All Figures

thumbnail Fig. 1

HST ACS/WFC image of Terzan 5 in the F814W filter. The field of view is 204″ × 204″. The black box corresponds to the VLT-MAD pointing (60″ × 60″). The region used for the SFH reconstruction is beyond the red circle, corresponding to a radius of 15″. North is up, east is to the left.

In the text
thumbnail Fig. 2

Optical (I, V–I) and hybrid (I, I–K) CMDs of Terzan 5 obtained through the combined HST and MAD observations, corresponding to the FOV delimited by the black box in Figure 1. No quality cuts and differential reddening corrections are applied to these diagrams. It is worth noticing the distortion of the evolutionary sequences along the direction of the reddening vectors, indicated by the black arrows in the upper-left part of the panels. The red line corresponds to a 12 Gyr old isochrone (Bressan et al. 2012; Marigo et al. 2017) with [M/H] = −0.2 dex, E(BV) = 2.38, and μ0 = 13.87 (reference parameters from the database of Valenti et al. 2007), reddenned by adopting RV = 3.1. The dashed black lines show the same isochrone shifted by assuming RV = 2.1, 2.6, 3.6, 4.1 (from left to right). We varied E(BV) and μ0 to fit the optical (I, V–I) CMD. However, despite the remarkable broadening of the evolutionary sequences, these alternative RV values are clearly unable to match the hybrid CMD.

In the text
thumbnail Fig. 3

Trends of the photometric quality parameters χ and sharpness as a function of the I magnitude resulting from the combined optical-IR reduction. Darker colors are associated to denser populated cells. Red points are stars excluded by the 3σ clipping procedure.

In the text
thumbnail Fig. 4

Cleaning procedure from field interlopers via PMs. Left panel: CMD of the stars included in the MAD footprint with measured PMs from Massari et al. (2015). Middle panel: CMD of the field population as selected from the VPDs, that is stars with PMs not compatible with that of Terzan 5. Right panel: PM-cleaned CMD of Terzan 5 obtained by considering only likely cluster members selected from the VPDs. Rightmost column: VPDs of the measured stars divided in 5 bins of 1.5 mag each, starting at IDRC = 16.5. Black dots represent stars within a radius of 1.5 mas yr−1 and correspond to likely cluster members according to the criterion of Massari et al. (2015).

In the text
thumbnail Fig. 5

Map of the absolute reddening E(BV) in the direction of Terzan 5. Darker colors indicate the more extincted areas within the observed FOV. The white cross marks the gravity center of Terzan 5 (Lanzoni et al. 2010). North is up, while east is on the left.

In the text
thumbnail Fig. 6

DRC and PM-cleaned CMDs of Terzan 5 obtained through the photometric reduction and the differential reddening map described in this work. The left-hand panel displays all the member stars measured in the FOV, while the right-hand panel shows the subsample of stars satisfying our quality cuts, specifically a radial distance of r ≥ 15″ from the center and good values of χ and asharpness (i.e., within 3σ of their distributions; see Fig. 3).

In the text
thumbnail Fig. 7

Upper panels: completeness curves as a function of the I and K magnitudes in different concentric annuli from the cluster’s center. The colored dots and vertical lines highlight the limit magnitude at which the 50% is reached within the respective radial annulus. Lower panels: photometric errors, i.e. the difference between the input and output magnitudes, in I and K bands recovered from the artificial star tests.

In the text
thumbnail Fig. 8

SFH of Terzan 5 associated to the SFERA best-fit solution. The red shaded areas are the uncertainties associated to the model computed from a bootstrap method, and they mark the 5th and 95th confidence level. To obtain this result, an average foreground extinction of E(BV) = 2.5 and distance modulus μ0 = 13.75 were adopted.

In the text
thumbnail Fig. 9

Upper panels: Hess diagrams of the observed DRC (left) and best-fit synthetic (right) CMD. Lower panels: color functions of the CMDs in four different magnitude bins; the data and relative Poissonian uncertainty are displayed in gray, while the models are the superimposed red line.

In the text
thumbnail Fig. 10

Cumulative stellar mass function for the recovered best-fit SFH.

In the text
thumbnail Fig. 11

Metallicity distribution function as derived from the best-fit solution.

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.