Fitting formulae for the effects of binary interactions on lick indices and colors of stellar populations

Single-star stellar population (ssSSP) models are usually used for spectral stellar population studies. However, more than 50% of stars are in binaries and evolve differently from single stars. This suggests that the effects of binary interactions should be considered when modeling the stellar populations of galaxies and star clusters. Via a rapid spectral stellar population synthesis model, we give detailed studies of the effects of binary interactions on the Lick indices and colors of stellar populations and on the determination of the stellar ages and metallicities of populations. Our results show that binary interactions make stellar populations less luminous and bluer, with larger age-sensitive Lick indices (H beta) and smaller metallicity-sensitive indices (e. g., Mg b, Fe5270, and Fe5335) compared to ssSSPs. It also shows that when ssSSP models are used to determine the ages and metallicities of stellar populations, smaller ages or metallicities are obtained when using two line indices (H beta and [MgFe]) or two colors (e. g., u - R and R - K), respectively. Some relations for linking the stellar population parameters obtained by ssSSPs to those obtained by binary-star stellar populations (bsSSPs) are presented in this work. This can help us get some absolute values for stellar population parameters and is useful for absolute studies. However, it is found that the relative luminosity-weighted stellar ages and metallicities obtained via ssSSPs and bsSSPs are similar. This suggests that ssSSPs can be used for most spectral stellar population studies, except in some special cases.

Establishing the reliability with which stellar population parameters can be measured is vital to extragalactic astronomy. Galactic globular clusters (GCs) provide an excellent medium in which to test the consistency of single stellar population (SSP) models as they should be our best analogue to a homogeneous (single) stellar population. Here we present age, metallicity and α-element abundance measurements for 48 Galactic GCs as determined from integrated spectra using Lick indices and SSP models from Thomas, Maraston & Korn, Lee & Worthey and Vazdekis et al. By comparing our new measurements to independent determinations we are able to assess the ability of these SSPs to derive consistent results – a key requirement before application to heterogeneous stellar populations like galaxies. We find that metallicity determinations are extremely robust, showing good agreement for all models examined here, including a range of enhancement methods. Ages and α-element abundances are accurate for a subset of our models, with the caveat that the range of these parameters in Galactic GCs is limited. We are able to show that the application of published Lick index response functions to models with fixed abundance ratios allows us to measure reasonable α-element abundances from a variety of models. We also examine the age–metallicity and [α/Fe]–metallicity relations predicted by SSP models, and characterize the possible effects of varied model horizontal branch morphology on our overall results.

This article is devoted to the study of the central stellar populations of the brightest cluster galaxies (BCGs). High signal-to-noise ratio, long-slit spectra for a sample of 39 galaxies were fitted against two stellar population models, Pegase.HR (P.HR) and Vazdekis/MILES, to determine the star formation histories of the galaxies using full spectrum fitting, to investigate, in particular, whether a single stellar population (SSP) or composite stellar population model provides the better description. Monte Carlo simulations and χ2 maps were used to check the reliability of the solutions. The ages and [Fe/H] were (i) compared with those derived with the Lick indices and (ii) tested against the internal galaxy properties (the velocity dispersions and absolute K-band magnitudes) and the properties of the host cluster environment (X-ray temperatures, luminosities, offsets and the presence of cooling flows (CFs)), to determine whether any statistically significant correlations existed. The results indicate that 79 per cent of the BCG sample could be represented by an SSP fit, while the remaining 21 per cent likely experienced more than one star formation epoch. The correlations showed that the BCGs hosted in CF clusters are generally found closer to the centre of the clusters, while the BCGs in non-CF clusters are generally found further away. The main results suggested that at least some of the galaxies in the BCG sample had a more complex star formation history than first assumed and that the presence of CFs in the clusters could account for some, but not all, of the star formation activity of BCGs.

In this work, predictions of the spectral energy distribution from populations of single and binary stars are incorporated into a galactic chemical and color evolution model to explore the significance of the effects of the binary interactions on the color evolution of M33. We first constructed a model without binary interactions, and the model is able to reproduce most of the available observational constraints on the distribution of stellar parameters. We then run simulations with the same set of model parameters but with binary interactions considered. By comparing the results for the populations with and without binary interactions, we find that the inclusion of binary interactions makes the surface brightness greater (similar to 0.1 mag arcsec(-2)) in FUV-band but smaller (similar to 0.7 mag arcsec(-2)) in K-band, while it results in the FUV-K color bluer (similar to 0.8 mag). To reproduce the observations, a model that considers the binary interactions should make more gas fall onto the disk in the early time of the galaxy evolution, or increase the total stellar mass, or both.

We first presented an isochrone data base that can be widely used for stellar population synthesis studies and colour-magnitude diagram (CMD) fitting. The data base consists of the isochrones of both single-star and binary-star simple stellar populations (ss-SSPs and bs-SSPs). The ranges for the age and metallicity of populations are 0-15 Gyr and 0.0001-0.03, respectively. All data are available for populations with two widely used initial mass functions (IMFs), that is, Salpeter IMF and Chabrier IMF. The uncertainty caused by the data base (about 0.81 per cent) is designed to be smaller than those caused by the Hurley code and widely used stellar spectra libraries (e.g. BaSeL 3.1) when it is used for stellar population synthesis. Based on the isochrone data base, we then built a rapid stellar population synthesis (RPS) model and calculated the high-resolution (0.3-angstrom) integrated spectral energy distributions, Lick indices and colour indices for bs-SSPs and ss-SSPs. In particular, we calculated the UBVRIJHKLM colours, ugriz colours and some composite colours that consist of magnitudes on different systems. These colours are useful for disentangling the well-known stellar age-metallicity degeneracy according to our previous work. As an example for applying the isochrone data base for CMD fitting, we fitted the CMDs of two star clusters (M67 and NGC 1868) and obtained their distance moduli, colour excesses, stellar metallicities and ages. The results showed that the isochrones of bs-SSPs are closer to those of real star clusters. It suggests that we should take the effects of binary interactions into account in stellar population synthesis. We also discussed on the limitations of the application of the isochrone data base and the results of the RPS model.

The optical integrated spectra of three LMC young stellar clusters (NGC 1984, NGC 1994 and NGC 2011) exhibit concave continua and prominent molecular bands which deviate significantly from the predictions of single stellar population (SSP) models. In order to understand the appearance of these spectra, we create a set of young stellar population (MILES) models, which we make available to the community. We use archival International Ultraviolet Explorer integrated UV spectra to independently constrain the cluster masses and extinction, and rule out strong stochastic effects in the optical spectra. In addition, we also analyze deep colour-magnitude diagrams of the clusters to provide independent age determinations based on isochrone fitting. We explore hypotheses including age-spreads in the clusters, a top-heavy initial mass function, different SSP models and the role of red supergiant stars (RSG). We find that the strong molecular features in the optical spectra can only be reproduced by modeling an increased fraction of about 20 per cent by luminosity of RSG above what is predicted by canonical stellar evolution models. Given the uncertainties in stellar evolution at Myr ages, we cannot presently rule-out the presence of Myr age-spreads in these clusters. Our work combines different wavelengths as well as different approaches (resolved data as well as integrated spectra for the same sample) in order to reveal the complete picture. We show that each approach provides important information but in combination can we better understand the cluster stellar populations.

This paper presents theoretical integrated spectral energy distributions (SEDs) of binary star composite stellar populations (bsCSPs) in early-type galaxies, and how the bsCSP model can be used for spectral studies of galaxies. All bsCSPs are built basing on three adjustable inputs (metallicity, ages of old and young components). The effects of binary interactions and stellar population mixture are taken into account. The results show some UV-upturn SEDs naturally for bsCSPs. The SEDs of bsCSPs are affected obviously by all of three stellar population parameters, and the effects of three parameters are degenerate. This suggests that the effects of metallicity, and the ages of the old (major in stellar mass) and young (minor) components of stellar populations should be taken into account in SED studies of early-type galaxies. The sensitivities of SEDs at different wavelengths to the inputs of a stellar population model are also investigated. It is shown that UV SEDs are sensitive to all of three stellar population parameters, rather than to only stellar age. Special wavelength ranges according to some SED features that are relatively sensitive to stellar metallicity, young-component age, and old-component age of bsCSPs are found by this work. For example, the shapes of SEDs with wavelength ranges of 5110-5250AA, 5250--5310AA, 5310--5350AA, 5830--5970AA, 20950--23550AA are relatively sensitive to the stellar metallicity of bsCSPs. The shapes of SEDs within 965-985AA, 1005--1055AA, 1205--1245AA are sensitive to old-component age, while SED features within the wavelength ranges of 2185--2245AA, 2455--2505AA, 2505--2555AA, 2775--2825AA, 2825--2875AA to young-component age.

We present a photometric and spectroscopic study of stellar populations in the X-ray–luminous cluster of galaxies RX J0142.0+2131 at z = 0.280. This paper analyzes the results of high signal-to-noise ratio spectroscopy, as well as g'-, r'-, and i'-band imaging, using the Gemini Multi-Object Spectrograph on Gemini North. Of 43 spectroscopic targets, we find 30 cluster members over a range in color. Central velocity dispersions and absorption-line strengths for lines in the range 3700 A λrest 5800 A are derived for cluster members and are compared with a low-redshift sample of cluster galaxies and single stellar population (SSP) models. We use a combination of these indicators to estimate luminosity-weighted mean ages, metallicities ([M/H]), and α-element abundance ratios ([α/Fe]). RX J0142.0+2131 is a relatively poor cluster and lacks galaxies with high central velocity dispersions. Although the red sequence and the Faber-Jackson relation are consistent with pure passive evolution of the early-type population with a formation redshift of zform 2, the strengths of the 4000 A break and scaling relations between metal line indices and velocity dispersion reject this model with high significance. By inverting SSP models for the HβG, Mg b, and Fe line indices, we calculate that, at a given velocity dispersion and metallicity, galaxies in RX J0142.0+2131 have luminosity-weighted mean ages 0.14 ± 0.07 dex older than the low-redshift sample. We also find that [α/Fe] in stellar populations in RX J0142.0+2131 is 0.14 ± 0.03 greater than at low redshift. All scaling relations are consistent with these estimated offsets. We speculate that the older luminosity-weighted mean ages and [α/Fe] enhancement can be brought about by a rapidly curtailed burst of star formation in RX J0142.0+2131, such as may be experienced in a cluster-cluster merger. We note that the cluster's velocity dispersion, 1278 ± 134 km s-1, is larger than expected from both its X-ray luminosity and richness. However, the velocity distribution of galaxies in RX J0142.0+2131 is consistent with being drawn from a Gaussian distribution, and no sign of substructure is found. We conclude that stellar populations in RX J0142.0+2131 cannot evolve into stellar populations similar to those seen in our low-redshift sample through passive evolution. This study provides further evidence that a more complex model, possibly involving ongoing or intermittent star formation and galaxy mergers, is required to describe the evolution of cluster galaxies.

Efficiency tests for estimating the gas and stellar population parameters in Type 2 objects

In our evolutionary population synthesis models, the samples of binaries are reproduced by the 'patched' Monte Carlo simulation and the stellar masses, integrated J, H, K, L, L2 and M magnitudes, mass-to-light ratios and broad colours involving infrared bands are presented, for an extensive set of instantaneous-burst binary stellar populations. In addition, the fluctuations in the integrated colours, which have been given by Zhang et al. (2005), are reduced. By comparing the results for binary stellar populations with (Model A) and without (Model B) binary interactions we show that the inclusion of binary interactions makes the stellar mass of a binary stellar population smaller (~3.6-4.5% during the past 15Gyr); magnitudes greater (except U, ~0.18mag at the most); colours smaller (~0.15mag for V-K at the most); the mass-to-light ratios greater (~0.06 for K-band) except those in the U and B passbands at higher metallicities. And, Binary interactions make the V magnitude less sensitive to age, R and I magnitudes more sensitive to metallicity. Given an age, the absolute values of the differences in the stellar mass, magnitudes, mass-to-light ratios (except those in the U and B bands) between Models A and B reach the maximum at Z=0.0001, i.e., the effects of binary interactions on these parameters reach the maximum, while the differences in some colours reach the maximum at Z ~0.01-0.0004. On the contrary, the absolute value of the difference in the stellar mass is minimal at Z=0.03, those in the U,B,V magnitudes and the mass-to-light ratios in the U and B bands reach the minimum at Z ~0.01-0.004.

We perform a direct comparison of two state-of-the art single stellar population (SSP) models that have been used to demonstrate the non-universality of the low-mass end of the Initial Mass Function (IMF) slope. The two public versions of the SSP models are restricted to either solar abundance patterns or solar metallicity, too restrictive if one aims to disentangle elemental enhancements, metallicity changes and IMF variations in massive early-type galaxies (ETGs) with star formation histories different from the solar neighborhood. We define response functions (to metallicity and \alpha-abundance) to extend the parameter space of each set of models. We compare these extended models with a sample of Sloan Digital Sky Survey (SDSS) ETGs spectra with varying velocity dispersions. We measure equivalent widths of optical IMF-sensitive stellar features to examine the effect of the underlying model assumptions and ingredients, such as stellar libraries or isochrones, on the inference of the IMF slope down to ~0.1 solar masses. We demonstrate that the steepening of the low-mass end of the Initial Mass Function (IMF) based on a non-degenerate set of spectroscopic optical indicators is robust against the choice of the stellar population model. Although the models agree in a relative sense (i.e. both imply more bottom-heavy IMFs for more massive systems), we find non-negligible differences on the absolute values of the IMF slope inferred at each velocity dispersion by using the two different models. In particular, we find large inconsistency in the quantitative predictions of IMF slope variations and abundance patterns when sodium lines are used. We investigate the possible reasons for these inconsistencies.

With the goal of assembling a new generation of more realistic single stellar population (SSP) models, we have obtained magnesium abundances for nearly 80% of the stars of the widely employed MILES empirical spectral library. Additional spectroscopic observations of carefully selected stars have recently been obtained to improve the parametric coverage of this library. Here we report on: (i) the framework of Mg abundance determination carried out at mid-resolution, (ii) the newly acquired data, and (iii) the preliminary steps towards modelling stellar populations.

We present the stellar population content of early-type galaxies from the Atlas3D survey. Using spectra integrated within apertures covering up to one effective radius, we apply two methods: one based on measuring line-strength indices and applying single stellar population (SSP) models to derive SSP-equivalent values of stellar age, metallicity, and alpha enhancement; and one based on spectral fitting to derive non-parametric star-formation histories, mass-weighted average values of age, metallicity, and half-mass formation timescales. Using homogeneously derived effective radii and dynamically-determined galaxy masses, we present the distribution of stellar population parameters on the Mass Plane (M_JAM, Sigma_e, R_maj), showing that at fixed mass, compact early-type galaxies are on average older, more metal-rich, and more alpha-enhanced than their larger counterparts. From non-parametric star-formation histories, we find that the duration of star formation is systematically more extended in lower mass objects. Assuming that our sample represents most of the stellar content of today's local Universe, approximately 50% of all stars formed within the first 2 Gyr following the big bang. Most of these stars reside today in the most massive galaxies (>10^10.5 M_sun), which themselves formed 90% of their stars by z~2. The lower-mass objects, in contrast, have formed barely half their stars in this time interval. Stellar population properties are independent of environment over two orders of magnitude in local density, varying only with galaxy mass. In the highest-density regions of our volume (dominated by the Virgo cluster), galaxies are older, alpha-enhanced and have shorter star-formation histories with respect to lower density regions.

The stellar population, metallicity distribution and ionized gas in the elliptical galaxies NGC 6868 and 5903 are investigated in this paper by means of long-slit spectroscopy and stellar population synthesis. Lick indices in both galaxies present a negative gradient indicating an overabundance of Fe, Mg, Na and TiO in the central parts with respect to the external regions. We found that Mg2 correlates both with Fe i λ5270 and Fe i λ5335, suggesting that these elements probably underwent the same enrichment process in NGC 6868. However, only a marginal correlation of Mg2 and Fe i λ5270 occurs in NGC 5903. The lack of correlation between computed galaxy mass and the Mg2 gradient suggests that these elliptical galaxies were formed by merger events. The stellar population synthesis shows the presence of at least two populations with ages of 13 and 5 Gyr old in both galaxies. We have estimated the metallicity of the galaxies using single-aged stellar population (SSP) models. The central region of NGC 6868 (|R| ≲ 0.5 kpc) presents a deficiency of alpha elements with respect to iron and solar metallicity. The external parts present a roughly uniform distribution of [α/Fe] ratios and metallicities ranging from [Z/Z⊙]=− 0.33 and solar. A similar conclusion applies to NGC 5903. Concerning the emitting gas conspicuously detected in NGC 6868, we test three hypotheses as ionizing source: an H ii region, post-AGB (asymptotic giant branch) stars and an active galactic nucleus (AGN). Diagnostic diagrams involving the ratios ([N ii] λ6584)/(Hα), ([O i] λ6300)/(Hα) and ([S ii] λλ6717, 6731)/(Hα) indicate that values measured in the central region of NGC 6868 are typical of LINERs (Low-Ionization Nuclear Emission-Line Regions). Together with the stellar population synthesis, this result suggests that the main source of gas ionization in NGC 6868 is non-thermal, produced by a low-luminosity AGN, probably with some contribution of shocks to explain ionization at distances of ∼3.5 kpc from the nucleus.

We present the integrated J, H, K, L, M and N magnitudes and the colours involving infrared bands, for an extensive set of instantaneous-burst binary stellar populations (BSPs) by using evolutionary population synthesis (EPS). By comparing the results for BSPs WITH and WITHOUT binary interactions we show that the inclusion of binary interactions makes the magnitudes of populations larger (fainter) and the integrated colours smaller (bluer) for τ ≥ 1 Gyr. Also, we compare our model magnitudes and colours with those of Bruzual & Charlot (2003, hereafter BC03) and Maraston (2005, hereafter M05). At last, we compare these model broad colours with Magellanic Clouds globular clusters (GCs) and Milky Way GCs. In (V − R)−[Fe/H] and (V − I)−[Fe/H] diagrams it seems that our models match the observations better than those of BC03 and M05.