No Arabic abstract
Fitting whole spectra at intermediate spectral resolution (R = 1000 -- 3000), to derive physical properties of stellar populations, appears as an optimized alternative to methods based on spectrophotometric indices: it uses all the redundant information contained in the signal. This paper addresses the validation of the method and it investigates the quality of the population models together with the reliability of the fitting procedures. We are using two algorithms: STECKMAP, a non-parametric regularized program and NBURSTS a parametric non-linear minimization. We compare three spectral synthesis models for single stellar populations: Pegase-HR, Galaxev (BC03) and Vazdekis/Miles, and we analyse spectra of Galactic clusters whose populations are known from studies of color-magnitude diagrams (CMD) and spectroscopy of individual stars. We find that: (1) The quality of the models critically depends on the stellar library they use. Pegase-HR and Vazdekis/Miles are consistent, while the comparison between Pegase-HR and BC03 shows some systematics reflecting the limitations of the stellar library (STELIB) used to generate the latter models; (2) The two fitting programs are consistent; (3) For globular clusters and M67 spectra, the method restitutes metallicities in agreement with spectroscopy of stars within 0.14 dex; (4) The spectroscopic ages are very sensitive to the presence of a blue horizontal branch (BHB) or of blue stragglers. A BHB morphology results in a young SSP-equivalent age. Fitting a free amount of blue stars in addition to the SSP model to mimic the BHB improves and stabilizes the fit and restores ages in agreement with CMDs studies. This method is potentially able to disentangle age or BHB effects in extragalactic clusters.
We investigate the precision of the ages and metallicities of 21,000 mock simple stellar populations (SSPs) determined through full-spectrum fitting. The mock SSPs cover an age range of 6.8 $<$ log (age/yr) $<$ 10.2, for three wavelength ranges in the optical regime, using both Padova and MIST isochrone models. Random noise is added to the model spectra to achieve S/N ratios between 10 to 100 per wavelength pixel. We find that for S/N $geq$ 50, this technique can yield ages of SSPs to an overall precision of $Delta,mbox{log(age/yr)} sim 0.1$ for ages in the ranges 7.0 $leq$ log (age/yr) $leq$ 8.3 and 8.9 $leq$ log (age/yr) $leq$ 9.4. For the age ranges of 8.3 $leq$ log (age/yr) $leq$ 8.9 and log (age/yr) $geq$ 9.5, which have significant flux contributions from asymptotic giant branch (AGB) and red giant branch (RGB) stars, respectively, the age uncertainty rises to about $pm 0.3$ dex. The precision of age and metallicity estimation using this method depends significantly on the S/N and the wavelength range used in the fitting. We quantify the systematic differences in age predicted by the MIST and Padova isochrone models, due to their different assumptions about stellar physics in various important (i.e., luminous) phases of stellar evolution, which needs to be taken in consideration when comparing ages of star clusters obtained using these popular models. Knowing the strengths and limitations of this technique is crucial in interpreting the results obtained for real star clusters and for deciding the optimal instrument setup before performing the observations.
Dwarf galaxies are generally faint. To derive their age and metallicity distributions, it is critical to optimize the use of any collected photon. Koleva et al., using full spectrum fitting, have found strong population gradients in some dwarf elliptical galaxies. Here, we show that the population profiles derived with this method are consistent and more precise than those obtained with spectrophotometric indices. This allows studying fainter objects in less telescope time.
Stellar evolution theory has been extraordinarily successful at explaining the different phases under which stars form, evolve and die. While the strongest constraints have traditionally come from binary stars, the advent of asteroseismology is bringing unique measures in well-characterised stars. For stellar populations in general, however, only photometric measures are usually available, and the comparison with the predictions of stellar evolution theory have mostly been qualitative. For instance, the geometrical shapes of isochrones have been used to infer ages of coeval populations, but without any proper statistical basis. In this chapter we provide a pedagogical review on a Bayesian formalism to make quantitative inferences on the properties of single, binary and small ensembles of stars, including unresolved populations. As an example, we show how stellar evolution theory can be used in a rigorous way as a prior information to measure the ages of stars between the ZAMS and the Helium flash, and their uncertainties, using photometric data only.
In this work we present chemical abundances of heavy elements (Z$>$28) for a homogeneous sample of 1059 stars from HARPS planet search program. We also derive ages using parallaxes from Hipparcos and Gaia DR1 to compare the results. We study the [X/Fe] ratios for different populations and compare them with models of Galactic chemical evolution. We find that thick disk stars are chemically disjunt for Zn and Eu. Moreover, the high-alpha metal-rich population presents an interesting behaviour, with clear overabundances of Cu and Zn and lower abundances of Y and Ba with respect to thin disk stars. Several abundance ratios present a significant correlation with age for chemically separated thin disk stars (regardless of their metallicity) but thick disk stars do not present that behaviour. Moreover, at supersolar metallicities the trends with age tend to be weaker for several elements.
We present a new age-dating technique that combines gyrochronology with isochrone fitting to infer ages for FGKM main-sequence and subgiant field stars. Gyrochronology and isochrone fitting are each capable of providing relatively precise ages for field stars in certain areas of the Hertzsprung-Russell diagram: gyrochronology works optimally for cool main-sequence stars, and isochrone fitting can provide precise ages for stars near the main-sequence turnoff. Combined, these two age-dating techniques can provide precise and accurate ages for a broader range of stellar masses and evolutionary stages than either method used in isolation. We demonstrate that the position of a star on the Hertzsprung- Russell or color-magnitude diagram can be combined with its rotation period to infer a precise age via both isochrone fitting and gyrochronology simultaneously. We show that incorporating rotation periods with 5% uncertainties into stellar evolution models improves age precision for FGK stars on the main sequence, and can, on average, provide age estimates up to three times more precise than isochrone fitting alone. In addition, we provide a new gyrochronology relation, calibrated to the Praesepe cluster and the Sun, that includes a variance model to capture the rotational behavior of stars whose rotation periods do not lengthen with the square-root of time, and parts of the Hertzsprung-Russell diagram where gyrochronology has not been calibrated. This publication is accompanied by an open source Python package, stardate, for inferring the ages of main-sequence and subgiant FGKM stars from rotation periods, spectroscopic parameters and/or apparent magnitudes and parallaxes.