Synthetic stellar and SSP libraries as templates for Gaia simulations

ESA’s Gaia mission will collect low resolution spectroscopy in the optical range for ～10⁹ objects. Complete and up-to-date libraries of synthetic stellar spectra are needed to built algorithms aimed to automatically derive the classification and the parametrization of this huge amount of data. In addition, libraries of stellar spectra are one of the main ingredients of stellar population synthesis models, aiming to derive the properties of unresolved stellar populations from their integrated light. We present (a) the newly computed libraries of synthetic spectra built by the Gaia community, covering the whole optical range (300–1100 nm) at medium-high resolution of (0.3 nm) for stars spanning the most different types, from M to O, from A-peculiar to Emission lines to White Dwarfs, and (b) the implementation of those libraries in our SSP code (Tantalo in The Initial Mass Function 50 Years Later, 327:235 2005), exploring different stellar evolution models.     

Empirical calibration of the near-infrared Ca ii triplet – I. The stellar library and index definition

A new stellar library at the near-IR spectral region developed for the empirical calibration of the Ca  ii triplet and stellar population synthesis modelling is presented. The library covers the range λλ 8348–9020 at 1.5-Å (FWHM) spectral resolution, and consists of 706 stars spanning a wide range in atmospheric parameters. We have defined a new set of near-IR indices, CaT*, CaT and PaT, which mostly overcome the limitations of previous definitions, the former being specially suited for the measurement of the Ca  ii triplet strength corrected for the contamination from Paschen lines. We also present a comparative study of the new and the previous Ca indices, as well as the corresponding transformations between the different systems. A thorough analysis of the sources of index errors and the procedure to calculate them is given. Finally, index and error measurements for the whole stellar library are provided together with the final spectra.     

Atmospheric parameters of the B-supergiant HD 198478 from the UV spectra

We determined atmospheric parameters of the Galactic early B-supergiant HD 198478 (55 Cyg) from the UV silicon lines taken from the high-resolution 1150–1980 Å IUE spectra. TLUSTY numerical code was used to model the stellar atmosphere and to determine the temperature and surface gravity assuming a non-LTE plane parallel hydrostatic stellar atmosphere with microturbulence. The synthesized spectra were broadened by the IUE instrumental profile, rotational and macroturbulent velocity with ROTIN numerical code. The silicon 1264 Å, 1309 Å, 1312 Å, 1417 Å and 1294–1303 Å multiplet lines of different stages of ionization (Si II and Si III) and Balmer Hδ 4101 Å line were modeled, leading to the temperature, surface gravity, rotational and macroturbulent velocity values. Our results have shown that the line broadening cannot be explained by rotational velocity only, but additional macroturbulent velocity component should be taken into account. HD 198478 shows a significant degeneracy in velocity, which means that the individual contributions of the macroturbulence and rotation in the total velocity broadening cannot be distinguished. Adequate fit of TLUSTY models to the observed non-resonant silicon lines suggests that the non-LTE plane-parallel hydrostatic stellar model without wind contribution can be used to explain such lines. We have obtained similar results using the HST STIS spectra in the same procedure, showing that the IUE spectra, despite their lacking quality compared to the STIS spectra, are reliable enough in determination of the B supergiants’ photospheric parameters.     

The formation and evolution of massive stellar clusters in IC4662

We present a multiwavelength study of the formation of massive stellar clusters, their emergence from cocoons of gas and dust, and their feedback on surrounding matter. Using data that span from radio to optical wavelengths, including Spitzer and Hubble Space Telescope ACS observations, we examine the population of young star clusters in the central starburst region of the irregular Wolf–Rayet galaxy IC4662. We model the radio-to-infrared (IR) spectral energy distributions of embedded clusters to determine the properties of their Hii regions and dust cocoons (sizes, masses, densities, temperatures), and use near-IR and optical data with mid-IR spectroscopy to constrain the properties of the embedded clusters themselves (mass, age, extinction, excitation, abundance). The two massive star-formation regions in IC4662 are excited by stellar populations with ages of ～4 Myr and masses of ～3×10⁵ M_⊙ (assuming a Kroupa initial mass function). They have high excitation and subsolar abundances, and they may actually be comprised of several massive clusters rather than the single monolithic massive compact objects known as ‘super star clusters’ (SSCs). Mid-IR spectra reveal that these clusters have very high extinction values, A _V ～20–25 mag, and that the dust in IC4662 is well mixed with the emitting gas, not in a foreground screen.     

The stellar populations of starburst galaxies through near-infrared spectroscopy

We study the central (inner few hundred parsecs) stellar populations of four starburst galaxies (NGC 34, 1614, 3310 and 7714) in the near-infrared (NIR), from 0.8 to 2.4 μm, by fitting combinations of stellar population models of various ages and metallicities. The NIR spectra of these galaxies feature many absorption lines. For the first time, we fit simultaneously as much as 15 absorption features in the NIR. The observed spectra are best explained by stellar populations containing a sizable amount (20–56 per cent by mass) of ∼1-Gyr-old stellar population with thermally pulsing asymptotic giant branch stars. We found that the metallicity of the stars which dominates the light is solar. Metallicities substantially different from solar give a worse fit. Though the ages and metallicities we estimate using the NIR spectroscopy are in agreement with values from the literature based on the ultraviolet/optical, we find older ages and a larger age spread. This may be due to the fact that the optical is mostly sensitive to the last episode of star formation, while the NIR better maintains the record of previous stellar generations. Another interesting result is that the reddening estimated from the whole NIR spectrum is considerably lower than that based on emission lines. Finally, we find a good agreement of the free emission-line spectrum with photoionization models, using as input spectral energy distribution the synthetic composite template we derived as best fit.     

PopStar I: evolutionary synthesis model description

We present new evolutionary synthesis models for simple stellar populations for a wide range of ages and metallicities. The models are based on the Padova isochrones. The core of the spectral library is provided by the medium resolution Lejeune et al. atmosphere models. These spectra are complemented by Non Local Thermodynamic Equilibrium (NLTE) atmosphere models for hot stars that have an important impact on the stellar cluster's ionizing spectra: O, B and WR stellar spectra at the early ages, and spectra of post asymptotic giant branch stars and planetary nebulae, at intermediate and old ages. At young ages, our models compare well with other existing models, but we find that the inclusion of the nebular continuum, not considered in several other models, significantly reddens the integrated colours of very young stellar populations. This is consistent with the results of spectral synthesis codes particularly devised for the study of starburst galaxies. At intermediate and old ages, the agreement with the literature model is good and, in particular, we reproduce the observed colours of star clusters in Large Magellanic Cloud well. Given the ability to produce good integrated spectra from the far-ultraviolet to the infrared at any age, we consider that our models are particularly suited for the study of high-redshift galaxies. These models are available on the web site http://www.fractal-es.com/SEDmod.htm and also through the Virtual Observatory Tools on the PopStar server.     

Near-infrared template spectra of normal galaxies: k-corrections, galaxy models and stellar populations

We have observed 28 local galaxies in the wavelength range between 1 and 2.4 μm in order to define template spectra of the normal galaxies along the Hubble sequence. Five galaxies per morphological type were observed in most cases, and the resulting rms spread of the normalized spectra of each class, including both intrinsic differences and observational uncertainties, is about 1 per cent in K , 2 per cent in H and 3 per cent in J . Many absorption features can be accurately measured. The target galaxies and the spectroscopic aperture (7×53 arcsec²) were chosen to be similar to those used by Kinney et al. to define template UV and optical spectra. The two data sets are matched in order to build representative spectra between 0.1 and 2.4 μm. The continuum shape of the optical spectra and the relative normalization of the near-IR ones were set to fit the average effective colours of the galaxies of the various Hubble classes. The resulting spectra are used to compute the k -corrections of the normal galaxies in the near-IR bands, and to check the predictions of various spectral synthesis models: while the shape of the continuum is generally well predicted, large discrepancies are found in the absorption lines. Among the other possible applications, here we also show how these spectra can be used to place constraints on the dominant stellar population in local galaxies. Spectra and k -corrections are publicly available and can be downloaded from the web site http://www.arcetri.astro.it/~filippo/spectra .     

A comparison of optical and near-infrared colours of Magellanic Cloud star clusters with predictions of simple stellar population models

We present integrated JHK _S Two-Micron All-Sky Survey photometry and a compilation of integrated-light optical photoelectric measurements for 84 star clusters in the Magellanic Clouds. These clusters range in age from ≈200 Myr to >10 Gyr, and have [Fe/H] values from −2.2 to −0.1 dex. We find a spread in the intrinsic colours of clusters with similar ages and metallicities, at least some of which is due to stochastic fluctuations in the number of bright stars residing in low-mass clusters. We use 54 clusters with the most-reliable age and metallicity estimates as test particles to evaluate the performance of four widely used simple stellar population models in the optical/near-infrared (near-IR) colour–colour space. All models reproduce the reddening-corrected colours of the old (≥10 Gyr) globular clusters quite well, but model performance varies at younger ages. In order to account for the effects of stochastic fluctuations in individual clusters, we provide composite   B − V , B − J , V − J , V − K _S  and   J − K _S  colours for Magellanic Cloud clusters in several different age intervals. The accumulated masses for most composite clusters are higher than that needed to keep luminosity variations due to stochastic fluctuations below the 10 per cent level. The colours of the composite clusters are clearly distinct in optical–near-IR colour–colour space for the following intervals of age: >10 Gyr, 2–9 Gyr, 1–2 Gyr, and 200 Myr−1 Gyr. This suggests that a combination of optical plus near-IR colours can be used to differentiate clusters of different age and metallicity.     

The ultraviolet morphology of evolved populations

In this paper I present a summary of the recent investigations we have developed at the Stellar Atmospheres and Populations Research Group (GrAPEs-for its designation in Spanish) at INAOE and collaborators in Italy. These investigations have aimed at providing updated stellar tools for the analysis of the UV spectra of a variety of stellar aggregates, mainly evolved ones. The sequence of material here presented roughly corresponds to the steps we have identified as mandatory to properly establish the applicability of synthetic populations in the analyses of observational data of globular clusters and more complex aged aggregates. The sequence is composed of four main stages, namely, (a) the creation of a theoretical stellar data base that we have called UVBLUE, (b) the comparison of such data base with observational stellar data, (c) the calculation of a set of synthetic spectral energy distributions (SEDs) of simple stellar populations (SSPs) and their validation through a comparison with observations of a sample of galactic globular clusters (GGCs), (d) construction of models for dating local ellipticals and distant red-envelope galaxies. Most of the work relies on the analysis of absorption line spectroscopic indices. The global results are more than satisfactory in the sense that theoretical indices closely follow the overall trends with chemical composition depicted by their empirical counterparts (stars and GGCs).     

Population synthesis of H ii galaxies

We study the stellar population of galaxies with active star formation, determining ages of the stellar components by means of spectral population synthesis of their absorption spectra. The data consist of optical spectra of 185 nearby ( z 0.075) emission-line galaxies . They are mostly H  ii galaxies, but we also include some starbursts and Seyfert 2s, for comparison purposes. They were grouped into 19 high signal-to-noise ratio template spectra, according to their continuum distribution, absorption- and emission-line characteristics. The templates were then synthesized with a star cluster spectral base.
The synthesis results indicate that H  ii galaxies are typically age-composite stellar systems, presenting important contributions from generations up to as old as 500 Myr. We detect a significant contribution of populations with ages older than 1 Gyr in two groups of H  ii galaxies. The age distributions of stellar populations among starbursts can vary considerably despite similarities in the emission-line spectra. In the case of Seyfert 2 groups we obtain important contributions from the old population, consistent with a bulge.
From the diversity of star formation histories, we conclude that typical H  ii galaxies in the local Universe are not systems presently forming their first stellar generation.     

Evolutionary stellar population synthesis at high spectral resolution: optical wavelengths

We present the single stellar population (SSP) synthesis results of our new synthetic stellar atmosphere models library with a spectral sampling of 0.3 Å, covering the wavelength range from 3000 to 7000 Å for a wide range of metallicities (twice solar, solar, half solar and 1/10 solar). The stellar library is composed of 1650 spectra computed with the latest improvements in stellar atmospheres. In particular, it incorporates non-local thermodynamic equilibrium (LTE) line-blanketed models for hot  ( T _eff≥ 27 500 K)  , and LTE line-blanketed models (Phoenix) for cool  (3000 ≤ T _eff≤ 4500 K)  stars. Because of the high spectral resolution of this library, evolutionary synthesis models can be used to predict the strength of numerous weak absorption lines and the evolution of the profiles of the strongest lines over a wide range of ages. The SSP results have been calculated for ages from 1 Myr to 17 Gyr using the stellar evolutionary tracks provided by the Geneva and Padova groups. For young stellar populations, our results have a very detailed coverage of high-temperature stars with similar results for the Padova and Geneva isochrones. For intermediate and old stellar populations, our results, once degraded to a lower resolution, are similar to the ones obtained by other groups (limitations imposed by the stellar evolutionary physics notwidthstanding). The limitations and advantages of our models for the analysis of integrated populations are described. The full set of the stellar library and the evolutionary models are available for retrieval at the websites http://www.iaa.csic.es/~rosa and http://www.iaa.csic.es/~mcs/sed@ , or on request from the first two authors.     

Age and metallicity of compact stellar systems in Virgo and Fornax

From high signal-to-noise ratio GMOS-N and AAOmega spectra, we have measured line-strength indices in the Lick/IDS system for several luminous and spatially dispersed compact stellar systems (CSSs) located in the Virgo and Fornax galaxy clusters. We estimate their [α/Fe] ratios, ages and metallicities using a simple (single-burst) stellar population model. We confirm that the Virgo core region luminous CSSs (     ) contain ancient stellar populations with subsolar total metallicities, suggesting that they comprise the bright tail of M87's GC distribution. The two Virgo intracluster globular clusters have ages and metallicities consistent with the cluster core CSSs. Two Fornax luminous CSSs also have ancient stellar populations but are at the upper end of the Virgo CSS metallicity range, while the third (UCD3) appears to be relatively young, metal-rich and with a core + halo radial profile. Our results suggest that Fornax may contain an extra population of luminous CSSs formed more recently than the ancient GC-like systems found in both clusters.     

Spatially resolved kinematics and stellar populations of brightest cluster and group galaxies

We present an examination of the kinematics and stellar populations of a sample of three brightest group galaxies (BGGs) and three brightest cluster galaxies (BCGs) in X-ray groups and clusters. We have obtained high signal-to-noise ratio Gemini/Gemini South Multi-Object Spectrograph (GMOS) long-slit spectra of these galaxies and use Lick indices to determine ages, metallicities and α-element abundance ratios out to at least their effective radii. We find that the BGGs and BCGs have very uniform masses, central ages and central metallicities. Examining the radial dependence of their stellar populations, we find no significant velocity dispersion, age, or α-enhancement gradients. However, we find a wide range of metallicity gradients, suggesting a variety of formation mechanisms. The range of metallicity gradients observed is surprising, given the homogeneous environment these galaxies probe and their uniform central stellar populations. However, our results are inconsistent with any single model of galaxy formation and emphasize the need for more theoretical understanding of both the origins of metallicity gradients and galaxy formation itself. We postulate two possible physical causes for the different formation mechanisms.     

Near-IR spectral evolution of dusty starburst galaxies

We propose a multicomponent analysis of starburst galaxies, based on a model that takes into account the young and evolved stellar components and the gas emission, with their respective extinction, in the frame of a coherent dust distribution pattern. Near-IR signatures are preferentially investigated, in order to penetrate as deep as possible into the dusty starburst cores. We computed the 1.4-2.5 μm spectra of synthetic stellar populations evolving through strong, short timescale bursts of star formation (continuum and lines, R ? 500). The evolution model is specifically sensitive to cool stellar populations (AGB and red supergiant stars). It takes advantage of the stellar library of Lançon & Rocca-Volmerange (1992) [A&ASS, 96, 593], observed with the same instrument (FTS/CFHT) as the analysed galaxy sample, so that the instrumental effects are minimised. The main near-IR observable constraints are the molecular signatures of CO and H₂O and the slope of the continuum, observed over a range exceptionally broad for spectroscopic data. The H - K colour determined from the spectra measures the intrinsic stellar energy distribution but also differential extinction, which is further constrained by optical emission line ratios. Other observational constraints are the near-IR emission lines (Brγ, He I 2.06 μm, [Fe II] 1.64 μm, H₂ 2.12 μm) and the far-IR luminosity. The coherence of the results relies on the interpretation in terms of stellar populations from which all observable properties are derived, so that the link between the various wavelength ranges is secured. The luminosity L_K is used for the absolute calibration.We apply this approach to the typical spectrum of the core of NGC 1614. Consistent solutions for the starburst characteristics (star-formation rate, IMF, burst age, morphology) are found and the role of each observational constraint in deriving satisfactory models is extensively discussed. The acceptable contamination of the K band light by the underlying population amounts ≥ 15% even through a 5 arcsec aperture. The model leads to a limit on the direct absorption of Lyman continuum photons by dust situated inside the ionised areas, which in turn, with standard gas-to-dust ratios, translates into small characteristic sizes for the individual coexisting H II regions of the massive starburst area (clusters containing ∼ 10² ionising stars). We show that room is left for IMFs extending to 120 M_⊙, rather than truncated at ∼ 60 M_⊙ as most conservative studies conclude. High internal velocity dispersions (≥ 20 km s⁻¹) are then needed for the H II regions. An original feature of this work is to base the analysis of near-infrared spectral galaxy observations on a large wavelength range, using models constructed with spectral stellar data observed with the same instrument. However a broader use of this spectral evolution model on other spectral or photometric data samples is possible if the spectral resolution of the model is adapted to observations or if colours are derived from the energy distributions.Catherine J. Cesarsky     

<Emphasis Type="SmallCaps">pacce</Emphasis>: Perl algorithm to compute continuum and equivalent widths

We present Perl Algorithm to Compute continuum and Equivalent Widths (pacce). We describe the methods used in the computations and the requirements for its usage. We compare the measurements made with pacce and “manual” ones made using iraf splot task. These tests show that for synthetic simple stellar population (SSP) models the equivalent widths strengths are very similar (differences ?0.2 Å) for both measurements. In real stellar spectra, the correlation between both values is still very good, but with differences of up to 0.5 Å. pacce is also able to determine mean continuum and continuum at line center values, which are helpful in stellar population studies. In addition, it is also able to compute the uncertainties in the equivalent widths using photon statistics. The code is made available for the community through the web at http://www.if.ufrgs.br/~riffel/software.html.     

Evolutionary history of the elliptical galaxy NGC 1052

We have obtained Keck spectra for 16 globular clusters (GCs) associated with the merger remnant elliptical NGC 1052, as well as a long-slit spectrum of the galaxy. We derive ages, metallicities and abundance ratios from simple stellar population models using the recently published methods of Proctor & Sansom , applied to extragalactic GCs for the first time. A number of GCs indicate the presence of strong blue horizontal branches that are not fully accounted for in the current stellar population models. We find all of the GCs to be ∼13 Gyr old according to simple stellar populations, with a large range of metallicities. From the galaxy spectrum we find NGC 1052 to have a luminosity-weighted central age of ∼2 Gyr and metallicity of  [Fe/H]∼+0.6  . No strong gradients in either age or metallicity were found to the maximum radius measured  (0.3  r _e≃ 1 kpc)  . However, we do find a strong radial gradient in α-element abundance, which reaches a very high central value. The young central starburst age is consistent with the age inferred from the H  i tidal tails and infalling gas of ∼1 Gyr. Thus, although NGC 1052 shows substantial evidence for a recent merger and an associated starburst, it appears that the merger did not induce the formation of new GCs, perhaps suggesting that little recent star formation occurred. This interpretation is consistent with 'frosting' models for early-type galaxy formation.     

<Emphasis Type="Italic">GALEX</Emphasis> and star formation

Wide-field far-UV (FUV, 1344–1786 Å) and near-UV (NUV, 1771–2831 Å) imaging from GALEX provides a deep, comprehensive view of the young stellar populations in hundreds of nearby galaxies, shedding new light on the process of star formation (SF) in different environments, and on the interplay between dust and SF. GALEX’s FUV-NUV color is extremely sensitive to stellar populations of ages up to a few hundred Myrs, unambiguously probing their presence and enabling age-dating and stellar mass estimate, together with the characterization of interstellar dust extinction. The deep sensitivity, combined with the wide field-of-view, made possible in particular the discovery and characterization of star formation in extremely low-density, diffuse gas environments such as outer galaxy disks, tidal tails, low-surface-brightness galaxies (LSB) and dwarf Irregular galaxies, and of rejuvenation episodes in early-type galaxies. Such results provide several missing links for interpreting galaxy classes in an evolutionary context, extend our knowledge of the star-formation process to previously unexplored conditions, constrain models of galaxy disk formation, and clarify the mutual role of dust and star formation. We review a variety of star-forming environments studied by GALEX, and provide some model analysis tools useful for interpretation of GALEX measurements, and potentially as basic science planning tools for next-generation UV instruments.     

Optical and near-infrared colours as a discriminant of the age and metallicity of stellar populations

We present a comprehensive analysis of the ability of current stellar population models to reproduce the optical ( ugriz ) and near-infrared ( JHK ) colours of a small sample of well-studied nearby elliptical and S0 galaxies. We find broad agreement between the ages and metallicities derived using different population models, although different models show different systematic deviations from the measured broad-band fluxes. Although it is possible to constrain simple stellar population models to a well-defined area in age–metallicity space, there is a clear degeneracy between these parameters even with such a full range of precise colours. The precision to which age and metallicity can be determined independently, using only broad-band photometry with realistic errors, is  Δ[Fe/H]≃ 0.18  and  Δlog Age ≃ 0.25  . To constrain the populations and therefore the star formation history further, it will be necessary to combine broad-band optical–IR photometry with either spectral line indices, or else photometry at wavelengths outside this range.     

Comparisons between stellar models and reliability of the theoretical models

The high quality of the asteroseismic data provided by space missions such as CoRoT (Michel et al. in The CoRoT Mission, ESA Spec. Publ. vol. 1306, p. 39, 2006) or expected from new operating missions such as Kepler (Christensen-Dalsgaard et al. in Commun. Asteroseismol. 150:350, 2007) requires the capacity of stellar evolution codes to provide accurate models whose numerical precision is better than the expected observational errors (i.e. below 0.1 μHz on the frequencies in the case of CoRoT). We present a review of some thorough comparisons of stellar models produced by different evolution codes, involved in the CoRoT/ESTA activities (Monteiro in Evolution and Seismic Tools for Stellar Astrophysics, 2009). We examine the numerical aspects of the computations as well as the effects of different implementations of the same physics on the global quantities, physical structure and oscillations properties of the stellar models. We also discuss a few aspects of the input physics.