Chemical composition of the young massive cluster NGC 1569-B

We present a detailed chemical abundance analysis of the young massive cluster (YMC) NGC 1569-B. The host galaxy, NGC~1569, is a dwarf irregular starburst galaxy at a distance of 3.36$\pm$0.20 Mpc. We determined the abundance ratios from the analysis of an optical integrated-light spectrum of NGC 1569-B, obtained with the HIRES echelle spectrograph on the Keck I telescope. We considered different red-to-blue supergiant ratios, namely: the ratio obtained from a theoretical isochrone, the ratio obtained from a resolved colour-magnitude diagram of the YMC, and the ratio that minimises the $\chi^2$ when comparing our model spectra with the observations. We adopted the latter ratio for our resulting chemical abundances. The derived iron abundance is sub-solar with [Fe/H] = $-0.74\pm0.05$. In relation to the scaled solar composition, we find enhanced $\alpha$-element abundances, $\text{[/Fe]}=+0.25\pm$0.11, with a particularly high Ti abundance of +0.49$\pm$0.05. Other super-solar elements include $\text{[Cr/Fe]}=+0.50\pm$0.11, $\text{[Sc/Fe]}=+0.78\pm$0.20, and $\text{[Ba/Fe]}=+1.28\pm$0.14, while other Fe-peak elements are close to scaled solar abundances: ($\text{[Mn/Fe]}=-0.22\pm$0.12 and $\text{[Ni/Fe]}=+0.13\pm$0.11). The composition of NGC 1569-B resembles the stellar populations of the YMC NGC 1705-1, located in a blue compact dwarf galaxy. The two YMCs agree with regard to $\alpha$-elements and the majority of the Fe-peak elements, except for Sc and Ba, which are extremely super-solar in NGC~1569-B -- and higher than in any YMC studied so far. The blue part of the optical spectrum of a young population is still a very challenging wavelength region to analyse using IL spectroscopic studies. This is due to the uncertain contribution to the light from blue supergiant stars, which can be difficult to disentangle from turn-off stars, even when resolved photometry is available.


Introduction
The element abundances of star clusters can reveal information on their star formation (SF) history and the evolution of their host galaxies. Indeed, obtaining accurate abundances for star clusters is a crucial goal for both galactic and extragalactic targets, since these parameters facilitate studies of the evolution of both stars and galaxies.
In recent decades, the availability of 8-10 m class telescopes has made it possible to study the spectra coming from stars outside of the Milky Way (MW), including those located in neighbouring galaxies such as the Small and Large Magellanic Clouds (SMC and LMC) as well as nearby dwarf galaxies (Hill 1999). It has been found that for old stellar populations with low iron abundances ([Fe/H]), the α-abundances of stars in dwarf galaxies exhibit similar patterns to those characterising the MW, with enhanced abundances (relative to scaled-solar composition). However, as [Fe/H] increases, the [α/Fe] ratio in dwarfs reaches lower values than in the MW (Tolstoy et al. 2009).
Beyond the Local Group, detailed information on the chemical composition of stellar populations is much more scarce, especially with regard to detailed abundance patterns. Based on observations of H ii regions, abundances of a limited set of elements (such as oxygen and nitrogen) can be inferred in galaxies with ongoing or recent SF (Shields 1990), and the spectroscopy of individual supergiant (SG) stars has been used to measure overall metallicities (Larsen et al. 2008;Davies et al. 2013). However, differences in the elements typically included in these types of analyses have made it difficult to compare these results directly with the data for old stellar populations in nearby galaxies.
The chemical composition of young massive clusters (YMCs) provides valuable information on stellar population. These objects allow us to gain insight into SF under more extreme conditions than those encountered in nearby star-forming regions, and they are of particular interest, particularly since an ongoing debate in the literature has questioned whether these objects might be the progenitors of the old globular clusters (GCs) that we see today (Lardo et al. 2017;Cabrera-Ziri et al. 2020). Additionally, it is interesting to compare the chemical abundance patterns that we see in our Local Group with systems outside of it. Thus, it is important to study extragalactic clusters to explore their chemical composition, and YMCs can serve as ideal targets thanks to their intrinsic high luminosities.
The high light-to-mass ratios of young clusters, such as YMCs, allow us to study their abundances using integrated-light (IL) spectroscopy. The IL spectra for the majority of star clusters are broadened by only up to 10 km s −1 which implies that relatively high spectral resolutions (λ/∆λ > ∼ 10,000) are accessible; thus, detailed chemical abundances can be extracted from their IL spectra. In general, this is done by modelling the IL spectra using information from the Hertzsprung-Russell diagram (HRD), combined with atmospheric models and synthetic stellar spectra (Davies et al. 2010;Larsen et al. 2012;Sakari et al. 2021). One method that uses this technique was introduced and described in Larsen et al. (2012), and it has been used to derive abundances for iron and α-elements, as well as light-and n-capture elements for old GCs in Local Group dwarf and spiral galaxies (Larsen et al. 2014(Larsen et al. , 2022. Extragalactic YMCs studied with this method include the clusters NGC 1313-379 (in the barred spiral NGC 1313) and NGC 1705-1, in a blue compact dwarf (Hernandez et al. 2017). The derived metallicities were: It is important to increase the number of YMCs for which the chemical composition has been measured in detail. In Larsen et al. (2006), iron and α abundances were derived for a young super-star cluster SSC in the nearby spiral galaxy NGC 6946 using near-infrared (NIR) observations from the Near Infrared Spectrograph (NIRSPEC) instrument on the Keck II telescope. These authors found a super-solar [α/Fe] and sub-solar metallicity of [Fe/H] = −0.45±0.08. Other studies of YMCs include Lardo et al. (2017) and Cabrera-Ziri et al. (2016) where YMCs were studied to search for GC-like abundance patterns.
Such patterns can be used to infer star-to-star variations of some light elements (e.g. Na-O anticorrelation), which are linked to the presence of 'multiple populations' (MPs) in GCs (Bastian & Lardo 2018). In particular, it is interesting to try and identify the first signs of the MPs within the YMCs that would lead to the patterns seen in the GCs. However, no clear indications of the abundance patterns associated with MPs have yet been found in any clusters younger than about 2 Gyr (Martocchia et al. 2018(Martocchia et al. , 2019. In this study, we derive the detailed chemical abundances of a YMC NGC 1569-B. It is located in a dwarf irregular galaxy (dIrr) NGC 1569, which has two bright YMCs: NGC 1569-A (a binary cluster) and NGC 1569-B. The masses were estimated to be ≈ 10 6 M for the binary cluster and (4.4 ± 1.1) × 10 5 M for NGC 1569-B (Ho & Filippenko 1996;Sternberg 1998;Larsen et al. 2008). Both of these clusters are more massive than any of the young clusters in the MW or LMC (Anders et al. 2004). The cluster investigated in this work is the older and less bright one out of the two YMCs. Figure 1 shows a Hubble Space Telescope (HST) image of the target with the slit used for our spectroscopic observations. NGC 1569-B is particularly interesting given its mass, since such clusters are rare in the Milky Way and even in the Local Group. The massive clusters can provide clues about the MPs seen in old GCs. Table 1 lists the basic parameters of NGC 1569. The stellar population in this dIrr is dominated by stars that are younger than a few tens of Myrs, while some old stars are also present (Greggio et al. 1998). Anders et al. (2004) confirmed a major burst 25 Myr ago in NGC 1569. During this time, a large number of clusters (≈150) were formed, including the massive clusters NGC 1569-A and NGC 1569-B. Evidence of older episodes of SF from 1.5 Gyr to 150 Myr ago has been presented by Vallenari  Notes. D -the distance to the galaxy; d -maximum optical size; M gastotal atomic gas mass (H I + He); M tot -total galaxy mass.
We wanted to see what effect the different ratios of red-toblue supergiants, N RS G /N BS G , would have on the derived chemical abundances. This ratio in NGC 1569-B has been previously discussed in the literature. For example, in Larsen et al. (2008), the derived ratio was N RS G /N BS G = 2.50±0.66. However, Larsen et al. (2011) used resolved photometry and identified 69 red SGs (RSGs) and 45 blue SGs (BSGs) giving N RS G /N BS G = 1.53±0.09. The difference between the ratios in the Larsen et al.
Article number, page 2 of 18 A. Gvozdenko et al.: Chemical composition of the young massive cluster NGC 1569-B (2008) and Larsen et al. (2011) studies is mainly due to the way the BSGs were selected, which illustrates that it is challenging to separate BSGs and the bright main sequence (MS) stars into two distinguishable groups.
Here, we study the chemical composition of YMC NGC 1569-B, using a high-resolution optical IL spectrum and applying the same method and software created by Larsen et al. (2012). We derive abundance ratios for a larger set of elements than those previously measured, including several α-elements (Mg, Si, Ca, Ti), Fe-peak elements (Sc, Cr, Mn, Ni) and the neutron capture element: Ba. Section 2 presents the observations we used, which are described together with the data reduction pipeline. The main stages of the analysis method are explained in Section 3, including the investigation of an important uncertainty for this study: the red-to-blue supergiant ratio (Section 3.2). The results are given in Section 4. We discuss and compare our findings to other YMCs, the MW stars, and LMC stars in Section 5. We offer a discussion of the supergiant ratios in Section 5.5. In Section 6, we present our main conclusions.

IL spectrum
NGC 1569-B was observed with the high-resolution echelle spectrometer (HIRES) on the Keck I 10-m telescope on the 1st of November 2018. The exposure time was 2 × 1800 s with the C5 decker, a 7 × 1 . 148 slit. This gives a resolving power of R = λ/∆λ = 37500. The right ascension and declination are 04:30:49.01 and +64:50:52.7 (J2000), respectively. These data were reduced using the highly automated MAKEE package (written by T. Barlow) that was specifically designed for HIRES data. The pipeline performs bias subtraction and flat fielding, determines the location of the echelle orders on the images and extracts the spectra. The spectra of Th-Ar calibration lamps mounted in HIRES were used to calibrate the wavelength. Afterwards, the reduced two spectra were used to produce the final IL spectrum and an error spectrum. The final spectrum consists of a sum of two observations covering the ≈3900-8150 Å wavelength range; whereas the error spectrum is a difference between the two observations, providing an uncertainty for every pixel in the spectrum. The error spectrum was convolved with a boxcar kernel, using the running mean with window of N=100 pixels, this helped to remove the anomalies from the error spectrum, such as telluric absorption (TA), cosmic rays, and so on. The signal-to-noise ratio (S/N) increases from the blue to the red end of the spectrum, with typical values of ∼50 at 4000 Å and ∼70 at 6000 Å.

HST photometry
The photometry data were obtained with the Advanced Camera for Surveys (ACS) onboard HST, for the programme GTO-9300 (P.I.: H. Ford). This provided us with F555W and F814W bands, which can be approximated to V and I. The exposure time was 390 s (3 × 130) for both filters. These data were previously used in Larsen et al. (2008) and the photometry is taken directly from it. As described in this paper, the photometry was foreground extinction-corrected and is not insignificant in the blue-optical (A V =1. 74 Israel (1988)). The available NGC1569-B photometry was used to find the best-fit isochrone, indicating an age of ∼25 Myr. The isochrone was used to perform the photometry test to see the contribution of light from red and blue supergiants (see Section 3.2).

Models and stellar parameters
The analysis of high-resolution IL data for star clusters is described in detail in Larsen et al. (2012Larsen et al. ( , 2014; Hernandez et al. (2017). Here, we give a brief overview of the steps particularly relevant to the present study.
The first step in the analysis is to split the cluster HRD into bins (up to 100) of different stellar evolutionary stages. This can be done using a combination of the isochrone and available photometry or the isochrones alone. Here, only the isochrone was used to create the bins because the photometry of the NGC 1569-B is not sufficiently deep to properly sample the CMD, preventing us from performing an accurate cluster and field separation. Each bin is described by mass, effective temperature (T eff ), logarithm of surface gravity (log(g)), radius of the star, weight of the bin, logarithm of microturbulent velocity (log(ν t )), and the atmospheric model used to produce the synthetic spectra. The two different atmospheric models used are: ATLAS9 for stars with T eff > 4000 K (Kurucz 1970) and for the cool stars (T eff < 4000 K) MARCS models (Gustafsson et al. 2008). The weight of each bin is calculated using the IMF and the mass of stars in the bins. We use the power law for the IMF, dN/dM ∝ M −α , with the exponent α = 2.35 from Salpeter (1955). The microturbulent velocities were allocated depending on the stellar effective temperature (similarly to Lyubimkov et al. (2004); Hernandez et al. (2017)): for stars with T eff < 6000 K the microturbulent velocities were ν t =2.0 km s −1 ; for stars with 6000 < T eff < 22000 K -ν t =4.0 km s −1 ; and stars with T eff > 22000 K had ν t =8.0 km s −1 .
Using the Kurucz and MARCS models, the model atmospheres and synthetic spectra are computed for each bin. The codes used for the spectral synthesis are SYNTHE (Kurucz & Furenlid 1979) and Turbospectrum (Plez 2012). Then these spectra are scaled to the stellar luminosity with respect to each bin. Afterwards, these are all combined into a single integrated synthetic spectrum. The abundances used in the modelling are then iteratively adjusted until the best fit to the observed spectra is obtained.
The method we use in this paper is based on the spectra of Arcturus (Hinkle et al. 2000) and the Sun (2005 version of Kurucz et al. (1984)) to define the spectral windows used for the fitting procedures. A crucial step is to ensure that each window contains lines with a range of equivalent widths, while also including the strong lines (equivalent width > 100 mÅ). The full and updated list of spectral windows can be found in Larsen et al. (2022), whereas the list used for this study is presented in B.1.
Since the method uses full spectral synthesis, blends are automatically accounted for and it is possible to access a wide variety of elements: α-abundances (Mg, Si, Ca, Ti), Fe-peak elements (Cr, Mn, Fe, Ni, and Sc), heavy elements (Cu, Zn, Zr, Ba, and Eu), and lighter elements such as Na.
The blending of the lines is an important detail and a main challenge for the IL spectrum analysis. The Larsen et al. (2012) method was previously applied to high-resolution spectra of GCs, however, young stellar clusters such as NGC 1569-B tend to have a higher degree of blending compared to metal-poor, old GCs due to the higher metallicity (although GCs up to nearly solar metallicity were included in Larsen et al. 2022). Therefore, to account for this, the analysed windows had to be carefully arranged to ensure that the element lines do not overlap with others produced by different types of elements. The spectral windows fitted are a combination of the previously used set of wavelength  (Girardi et al. 2000). Red empty circles show the theoretical isochrone bins used for the HRD file (Girardi et al. 2000).
regions (Larsen et al. 2022) and a new set of lines for the blue part of the spectrum (Table 2).

Ratio of red-to-blue supergiants
One of the major uncertainties that we aim to tackle in this analysis is the impact of the ratio of red-to-blue supergiant stars. RSG stars are cool stars which are the brightest stars in NIR light for young stellar populations. These stars evolve from the MS as their hydrogen is depleted and heavy elements come closer to the stellar surface. The following step of the stellar evolution is the BSG star phase, which are hot stars and the most luminous in the blue optical part of the spectrum (Pagel 2009). Given the very different effective temperatures of blue and red SG stars, they affect different parts of the spectrum. This effect is shown in Figures 3 and 4, which provide a comparison between the best-fit model of a stellar population with only BSG (shown in blue) and with only RSG (shown in red). It can be seen that at shorter wavelengths, the blue model spectrum tends to provide a better match overall to the observations; while at longer wavelengths, the RSGs with their more numerous features start dominating the integrated light. The wavelength range of the analysed IL spectrum is ∼ 4000 − 8000 Å and the main contribution in the blue up to about 5500 Å is from BSGs; while the red part of the spectrum is more dominated by the light from RSGs. Since the light from RSGs dominates the spectrum at redder wavelengths and has almost no effect at the shorter wavelengths (see Section 5.5.2), by increasing N RS G /N BS G the IL spectrum will become redder. As expected, the BSGs have the opposite effect on the spectrum, and decreasing the ratio will mainly affect the bluer part of the spectrum. These wavelength ranges indicate that the optimal contribution between the blue and the red SG stars has to be achieved.
The RSGs and BSGs bins were selected using the surface gravity and effective temperature from an isochrone: log(g) ≤ 2.5 and T eff ≤ 5000 K for RSGs, while log(g) ≤ 2.5 and T eff ≥ 5000 K for BSGs. The number of blue and red SGs were then calculated using the weight of each bin based on the IMF (Salpeter 1955). The original ratio within the isochrone used  for 25 Myrs and Z=0.008 from the 2008 version of the Padua isochrones has N RS G /N BS G = 1.24 (Girardi et al. 2000Marigo et al. 2008). However, another way to find the optimal red-to-blue SG ratio is to use statistics on the existing spectroscopic data. It is possible to combine the reduced χ 2 values of multiple elements best fits to calculate the total χ 2 value characterising the runs with different red-to-blue SGs ratios. Figure 5 shows how χ 2 varies for ratios between 0.5 and 2.0, namely, higher ratios between 2.0 and 4.0 fail to fit multiple Fe lines; hence, those ratios were not considered in the search for the lowest χ 2 . A visible trend is that the higher ratios of RSGs to BSGs are favoured by the analysis. The ratio with the lowest χ 2 value and therefore the most statistically accurate is N RS G /N BS G =1.90.

Chemical abundances
We derived the abundances for a number of individual elements from the NGC 1569-B IL spectrum, starting with the elements with the highest number of lines. The first step was to fit for the broadening and overall scaling of the solar reference abundances, [Z/H], modelling the whole wavelength range for every order. Once derived, the broadening can be fixed for all the subsequent runs of this YMC. The first individual element to fit is Fe since it has the largest number of absorption lines, we subsequently fit for Ti and so on. This procedure allows us to obtain the individual element abundances values for NGC 1569-B (Table B.1). Table 3 lists all the values of the averaged abundances derived for three different ratios of red to blue SGs: 1.24 -the original ratio from the isochrone, 1.53 -the observed ratio from Larsen et al. (2011), and 1.90 -based on the best χ 2 value from iteratively fitting the spectrum ( Figure 5). Section 5.5.2 shows that the isochrones reproduce the observed fraction of light from the RSGs quite well. However, for this particular spectral data the χ 2 argument favours the higher ratio. For some elements, there is some disagreement between  the blue and the red part of the spectrum (see the bottom plot in Figure 6 or the numerical values for Ca and Cr in the

Uncertainties
Each element had a number of spectral windows fitted. Table B.1 provides the full list of wavelength ranges, the derived values, and their uncertainties from iSPy3. Those values were then combined to produce the final value for single elements. The values and uncertainties listed in the Table 3 are the weighted average of the values derived for the spectral windows for every individual element: where the weights are defined as and where σ 0 = 0.05 dex is used to account for nonrandom uncertainties on the derived individual window abundances. Then the uncertainties were calculated using the following formula: A different approach is to take the standard errors (S X ) on the mean of the weighted standard deviations (S D X,w ):  Reddy et al. (2003Reddy et al. ( , 2006  The comparison between the random and systematic uncertainties are listed in Table 4. In particular, σ X is always lower than S X , except for the errors of Cu (this element had a single spectral window). This indicates that the errors calculated using the weighted average tend to underestimate the uncertainties, even if the floor value is used (σ 0 ). Therefore, weighted standard deviation average errors are used as a more accurate estimate of the true uncertainties.

Discussion
Observations of field stars and star clusters indicate that NGC 1569 experienced a major recent burst of SF about 25 Myr ago (Greggio et al. 1998;Anders et al. 2004). According to the chemical evolution models from Romano et al. (2006), following such bursty SF the cluster would be enriched with α elements, which is consistent with the enhanced α-element abundances that we measure ([α/Fe] = +0.25 ± 0.11). Currently, we have no information about abundance variations between the field star populations in NGC 1569 and its YMCs, NGC 1569-A, and NGC 1569-B. Therefore, it would be interesting to search for signatures of a bursty SFH also in other clusters of this dwarf galaxy.
Analyses of the blue part of the spectrum are not a common method for studying young stellar populations such as NGC 1569-B (usually studied using the IR band), due to challenges inherent in modelling BSG stars. This aspect can be clearly seen in Figure 6, where the scatter in individual elements in blue region is significantly larger than that seen at redder wavelengths. However, the blue region of the IL spectrum not only provides more lines for abundance analysis, but it also helps with determining and constraining the blue-to-red SGs ratios. Thus, while it is not perhaps the most efficient wavelength range for abundance analysis, it still gives us additional constraints on the properties of the BSG and MS stars and remains to be fully explored.  Figure 7 shows the α-abundance versus the iron abundance. These abundance ratios indicate a rapid and bursty SF, where chemical enrichment is dominated by SNII, which have a massive progenitor star and a short life (McWilliam 1997).
Individual α-elements are consistent with other YMCs previously studied, such as NGC 1705-1 whose abundance ratios [Mg/Fe], [Ca/Fe], and [Ti/Fe] resemble those in NGC 1569-B very closely (see Figure 8), as expected in the bursty models from Romano et al. (2006). For [Si/Fe] Figure 8 suggests the similarity with Fornax, LMC, and MW thick disk, however, these stars are older and currently there are no YMCs with a similar [Fe/H]. Therefore, it is more difficult to conclude whether this super-solar value of [Si/Fe] is expected, nevertheless this α element is close to the total α-abundance.
In addition to atomic Ti i and Ti ii lines, we also fit the TiO1 band, which is sensitive to T eff and the IMF. The derived Ti abundances using only the TiO1 molecular band for three different ratios are: [Ti/Fe] 1.24 = +0.31±0.04, [Ti/Fe] 1.53 = +0.34±0.04, and [Ti/Fe] 1.90 = +0.33±0.04, which behaves similarly to the overall α-abundance ratio. This behaviour of the TiO band compared to the total Ti abundance ratio is expected because as we go redder the RSGs will increasingly dominate, making the analysis less sensitive to the ratio of RSG to BSG stars. TiO1 is slightly more sensitive than TiO2 to [Ti/Fe] (La Barbera et al. 2013). This method only uses the IMF from Salpeter (1955), hence, the above stated values might change if a different IMF is used. However, given that we haven't explored other IMF choices, it is currently unclear how strong the constraint on the IMF from the TiO bands is for NGC 1569-B.

Fe-peak element abundances
The Fe-peak element abundances are shown in (Figure 9). These elements are the heaviest elements produced through thermonuclear reactions in SNIa and massive stars (Kobayashi et al. 2020;Minelli et al. 2021 (2017), offering two possible explanations: either there is a different contribution from massive stars and SNIa in these clusters versus the contribution in the MW and LMC or due to the inclusion of a 7700-7800 Å bin, which increased the final value. In particular, if the bin was excluded from the calculation, the [Ni/Fe] values for NGC 1705-1 and NGC 1313-379 would become +0.25 and +0.31, respectively -instead of +0.74±0.49 and +0.46±0.14 (as shown at the bottom plot of Figure 9).
The main outlier is the abundance ratio of Sc (produced in both SNIa and SNII), which seems extremely enhanced. We considered that this might be caused by an issue in measurement. However, despite changing the fitting combinations of the spline function used for continuum matching of the observed and model spectra (e.g. the power and the number of knots), [Sc/Fe] remained high at all times. We see a large range of Sc values across the various spectral window regions which might indicate that our Sc measurements are quite uncertain. For example, if we would exclude the spectral windows at 5522.2-5531.0 Å and 5638.0-5690.0 Å, the value of [Sc/Fe] would become 0.17±0.08, which places it very close to the stars of NGC 1705-1 with a Sc abundance ratio of 0.192±0.052. The analysis of NGC1705-1 had only a single window fitted between 6222.0-6244.0 Å, therefore, there is no overlap in the two analyses for Sc.
Neither of the previously studied SSCs and YMCs have shown such high abundance ratios for this element (see the top plot in Figure 9). It is true that Sc is often a tracer of α-abundance and while we do find that [α/Fe] is enhanced, this does not fully explain the high value of [Sc/Fe]. If the derived enhancement of Sc is real, this would imply that the contribution from SNIa and SNII in NGC 1569-B differs from the contributions in the MW, LMC, and other studied galaxies.

Heavy elements abundances
Heavy elements are those with atomic number Z > 30 (e.g. Ba and Eu) and they are neutron-capture elements. In this study, we derived a very high abundance ratio of Ba (the fourth plot in Figure 9). Ba is an s-process element, which is produced mainly in low-mass asymptotic giant branch (AGB) stars (1-3 M ) (Busso et al. 1999). However, it is highly unlikely that AGB stars would be present in a young population of 15-25 Myrs. Unlike the case of Sc, all of the fitted spectral windows yield consistently high Ba values, which suggests that it is a more reliable measurement. One possible explanation for such a high value of [Ba/Fe] could be pollution of the ISM in NGC 1569 coming from AGB stars belonging to earlier bursts of SF.
The Ba abundance is often different from the MW disk stars. For example, in Van der Swaelmen et al. (2013) [BaII/FeI] increases along with an increase of metallicity in the LMC bar, which indicates that LMC chemical enrichment was slower than the one of the MW. Other high values of Ba were obtained for Fornax, which is dominated by slow neutron capture (s-process) at iron-abundance of ≈ −1; this shows the strong role of the (metal-poor) AGB in its evolution (Letarte 2007). A minor contribution of Ba comes from massive stars (Busso et al. 1999). Conceivably, at low Z, the Ba peak nuclei do not receive the main contribution from the AGB stars, but from the short-lived massive stars instead. These stars produce heavy r-process nuclei that dominate the contribution (Truran 1981;Cowan et al. 1996;Busso et al. 1999).
Article number, page 10 of 18 A. Gvozdenko et al.: Chemical composition of the young massive cluster NGC 1569-B 5.5. Red-to-blue supergiants

Constraints on trends with metallicity from field stars in different environments
The red-to-blue SGs ratio greatly influences the IL spectra of young stellar populations. However, it is extremely sensitive to details of the stellar evolution models such as the mass loss, convection, and mixing processes, and this has posed a problem in stellar astrophysics for some time (Langer & Maeder 1995). This ratio was first shown to vary significantly among different galaxies by van den Bergh (1968). Eggenberger et al. (2002) spectroscopically confirmed the increase of this ratio when the metallicity decreases, additionally providing a normalised relation. A study of this ratio in stellar clusters (in the MW, the LMC and the SMC) was performed by Meylan & Maeder (1982), who reach a similar conclusion, namely, that the ratio of red-to-blue decreases steeply with increasing metallicity. These authors found a difference in the ratio between the MW and the SMC of approximately an order of magnitude. These results were revisited and corrected for incompleteness by Humphreys & McElroy (1984) who found that the centre of the MW has a red-to-blue SGs ratio that is ten times lower than of the SMC. In summary, the red-to-blue SGs ratio is a decreasing function of metallicity, as confirmed by studies in the MW, the Magellanic Clouds, and the Triangulum galaxy, M33 (Robertson 1973(Robertson , 1974Hagen & van den Bergh 1974;Ivanov 1998).
A study of NGC 1569-B using the CMD has shown a significant difference between the observed and the simulated (with single stellar populations (SSP) and multiple-bursts) red-to-blue SG ratios, the latter being three times smaller for NGC 1569-B (Larsen et al. 2011). The number of observationally detected BSGs is affected by the magnitude cuts and the lack of a clear Blue Hertzsprung Gap (BHG) (see below). This makes it challenging to clearly distinguish between BSGs and stars on the upper part of the main sequence. Eldridge et al. (2008) compared the observations with their models and found that the best results are achieved by including binary stars; otherwise, either there are too many RSGs, too few BSG, or a combination of the two issues. However, this alone is not enough to explain the observed ratios, and the challenge of the effect of enhanced mass loss and rotation still stands.
In the recent past, the study of red and blue SGs became one of the fundamental ways to explore the extragalactic young stellar populations. This tool allows one to obtain accurate abundance measurements outside of the Local Group (Bresolin et al. 2006;Bergemann et al. 2013;Gazak et al. 2014). The advantage of working with the RSGs is that the IL spectrum of a stellar cluster can be modelled with a single red SG star, simply because the NIR continuum of young stellar populations is almost entirely dominated by RSG stars (Davies et al. 2013). This approach to modelling the IL as a single RSG was used in Larsen et al. (2008) to study NGC 1569-B in the NIR.

Fraction of light versus wavelength (NGC 1569-B photometry versus isochrone)
The reason why the red-to-blue SG ratio is crucial for NGC 1569-B is that it is too young to have AGB or RGB stars. Therefore, most of the light of its spectrum comes from the RSGs and BSGs. Moreover, the impact of BSGs becomes more important at optical wavelengths, compared to previous studies in the IR, and must be properly taken into account. A key difficulty in implementing red and blue SGs in this analysis is the disagreement between the stellar evolution tracks and the observations. Evolutionary tracks tend to underestimate the extent of blue loops (He burning stage). The isochrones suggest a gap in the CMD between the H and He burning stellar stages -the BHG (Mengel et al. 1979). However, this gap is not always clear in the CMD, which results in larger uncertainties of the observed N RS G /N BS G , rather than when the uncertainties are predicted by Poissonian errors.
In our spectral analysis, the ratio of red-to-blue SGs can be varied within the input file containing the stellar population of the cluster, HRD file. We altered the ratio while keeping the total sum of the SG stars constant. Keeping the total number of SGs constant while varying the N RSG /N BSG ratio is motivated by the fact that the total amount of time a star spends as a SG (blue or red) is determined mainly by the total luminosity and amount of fuel available. Since the (bolometric) luminosity does not vary significantly between the BSG and RSG phases, we then expect that the total number of SGs is roughly independent of the ratio.
To test how well the isochrone reproduces the observed number of RSGs, we calculated the total integrated magnitudes in F555W (V-band) and F814W (I-band) for the observed CMD of NGC 1569-B and for the isochrone used. For NGC 1569-B, this was measured on the HST image within the radial range 15 < r < 50 pixels, where photometry for individual RSG stars was also available.
The total integrated magnitudes were compared with those of RSGs only, measured in the magnitude range of 21.0 < M V < 23.6 and for colours 2.0 < M V − M I < 3.0. The measured contribution of RSGs was shown to be in full agreement with that predicted by the isochrone in the V-band (22% and 21%); whereas for the I-band, it varied slightly: 48% for the isochrone and 53% for the photometry. Therefore, for NGC 1569-B, the isochrones reproduce the observed contributions of RSGs to the V-and Iband luminosities quite accurately. However, the fact that stars other than RSGs contribute significantly to the light shows that full SSP modelling is indeed necessary.

Abundance trends for varying red-to-blue supergiant ratios
The abundance trends seen with the increase of the N RS G /N BS G ratio act in an expected manner. Particularly there is an apparent decrease in metallicity as N RS G /N BS G increases. This occurs because in general, the spectral lines become stronger as a greater number of RSGs dominate the spectrum. To compensate, the abundances have to decrease to match the observed spectrum. However, α-abundance ratios relative to Fe stays approximately constant since the abundances of two of the α elements, [Mg/Fe] and [Si/Fe], increase as the ratio increases, while two other element abundances, [Ca/Fe] and [Ti/Fe], decrease. The rest of the element abundance ratios relative to Fe measured tend to decrease with the increase of the N RS G /N BS G . The ratio of red-to-blue SGs within the isochrone uses an SSP which does not always describe the population of this YMC accurately. For our case, the isochrone ratio was not the optimal one since it lead to significantly poorer solutions in the spectral fit although it reproduced the observed fraction of light from RSGs quite well. We also explored N RS G /N BS G derived from the photometry in Larsen et al. (2011). However, this photometry was not of ideal quality, hence, the task of resolving and sifting the members of the cluster was challenging. The final ratio used was derived from the minimum χ 2 method using the spectroscopic data available. We decided to use the set of abundance ra-Article number, page 11 of 18 A&A proofs: manuscript no. output tios derived with this approach (with N RS G /N BS G =1.90), because it was determined statistically for the data used in this study.
The availability of the photometry can help to constrain some of the stellar types, for instance, BSGs or RSGs. However, not all targets have observed photometry or it is simply not possible to obtain the photometric data. For these cases, the solution with a total χ 2 of varying the ratio of the most luminous stellar types (such as BSGs and RSGs for the YMCs) might help in determining the favoured and optimal solution. Additionally, a comparison among the blue and the red sides of spectrum might provide a useful route to determine the ratio in the integrated spectrum.

Multiple populations
The detailed chemical composition of a YMC can indicate the presence of MP, for example, Al variations or a strong [Na/Fe] enhancement (Lardo et al. 2017;Bastian et al. 2020). While neither of these elements were measured for NGC 1569-B, another sign for MP could be a lower [Mg/Fe] abundance when compared to [Ca/Fe] and [Ti/Fe] (Shetrone 1996;Kraft et al. 1997;Gratton et al. 2004;Colucci et al. 2009;Larsen et al. 2014

Conclusions
We performed a high-resolution IL spectrum analysis of the YMC NGC 1569-B. This is the first detailed abundance study of this YMC using the full optical wavelength range. The CMD of this cluster is well fit with an isochrone of 25 Myr. From our analysis of the IL spectrum, we find that NGC 1569-B is slightly metal-poor at [Fe/H] = −0.77 ± 0.05. The cluster displays α and the majority of the Fe-peak (Cr, Mn and Ni) element abundance ratios typical for MW thick disc or LMC bar stars (Reddy et al. 2006;Van der Swaelmen et al. 2013). We find that the cluster contains super-solar [Ba/Fe] abundances and it is unclear whether the contribution from massive stars can account for the strikingly high derived value of [Ba/Fe]. Another unexpectedly high abundance ratio is [Sc/Fe], one of the Fe-peak elements, usually produced in the SNII and SNIa. Although [α/Fe] is enhanced, and Sc is a tracer of α-abundance, this does not fully explain the derived abundances. A closer look at the fitting procedure also did not resolve the issue. Thus, we are presently unable to explain the extremely high value of [Sc/Fe].
Overall, this YMC resembles the element abundance ratios of another YMC that was the subject of previous studies, namely, NGC 1705-1, which is in agreement with the similar bursty star formations in both of the galaxies modelled by Romano et al. (2006). We note the importance of extending the number of studied YMCs, as, currently, the chemical composition of only four extragalactic YMC have been studied in detail. NGC 1569-B is particularly interesting due to its mass since clusters of this nature are rare in the MW and even in the Local Group. Such massive clusters may provide insights into the multiple populations seen in GCs. The biggest challenge for the analysis of YMCs in the blue spectral region is the proper treatment of the relative contributions of RSGs and BSGs. These stars have a major impact on the IL spectra of young clusters, hence, they must be modelled with care. While it remains challenging, the study of the optical IL of YMCs offers potentially important insights into the evolution of this stage of clusters and their host galaxies.