Stellar population models for young stellar clusters

Stellar population models are a very useful tool to derive the stellar cluster age and luminosity mass from its integrated light. Evolutionary synthesis models depend on the stellar evolutionary tracks and the stellar libraries used to predict the spectral energy distribution of the stellar populations. In this review, I briefly comment on the new models that have incorporated the new evolutionary tracks with rotation and new computations for the evolutionary tracks of the TP-AGB. A more extended summary is also given of the current status of the most recent high-resolution stellar libraries at optical wavelengths and their implementation in evolutionary synthesis models. A comparison of the results obtained fitting the optical spectra of LMC and SMC stellar clusters with different high spectral resolution evolutionary synthesis models is also presented.     

Near-IR age-dating of red supergiant-dominated stellar populations

We present new spectral synthesis models for solar metallicity stellar populations, based on a library of stellar spectra that extends across near-IR wavelengths out to 2.4 µm at a resolution approaching 1000. We show that the spectra of massive star clusters in the starburst galaxy M 82 can be reproduced very well with these models. We compare near-IR spectroscopic ages with optical ages, and discuss the main sources of (systematic) errors that still affect those ages.     

Velocity and abundance precisions for future high‐resolution spectroscopic surveys: A study for 4MOST

In preparation for future, large‐scale, multi‐object, high‐resolution spectroscopic surveys of the Galaxy, we present a series of tests of the precision in radial velocity and chemical abundances that any such project can achieve at a 4 m class telescope. We briefly discuss a number of science cases that aim at studying the chemo‐dynamical history of the major Galactic components (bulge, thin and thick disks, and halo) – either as a follow‐up to the Gaia mission or on their own merits. Based on a large grid of synthetic spectra that cover the full range in stellar parameters of typical survey targets, we devise an optimal wavelength range and argue for a moderately high‐resolution spectrograph. As a result, the kinematic precision is not limited by any of these factors, but will practically only suffer from systematic effects, easily reaching uncertainties <1km s^–1. Under realistic survey conditions (namely, considering stars brighter than r = 16 mag with reasonable exposure times) we prefer an ideal resolving power of R ∼20 000 on average, for an overall wavelength range (with a common two‐arm spectrograph design) of [395;456.5] nm and [587;673] nm. We show for the first time on a general basis that it is possible to measure chemical abundance ratios to better than 0.1 dex for many species (Fe, Mg, Si, Ca, Ti, Na, Al, V, Cr, Mn, Co, Ni, Y, Ba, Nd, Eu) and to an accuracy of about 0.2 dex for other species such as Zr, La, and Sr. While our feasibility study was explicitly carried out for the 4MOST facility, the results can be readily applied to and used for any other conceptual design study for high‐resolution spectrographs. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Isotope spectroscopy

The measurement of isotopic ratios provides a privileged insight both into nucleosynthesis and into the mechanisms operating in stellar envelopes, such as gravitational settling. In this article, we give a few examples of how isotopic ratios can be determined from high‐resolution, high‐quality stellar spectra. We consider examples of the lightest elements, H and He, for which the isotopic shifts are very large and easily measurable, and examples of heavier elements for which the determination of isotopic ratios is more difficult. The presence of ⁶Li in the stellar atmospheres causes a subtle extra depression in the red wing of the ⁷Li 670.7 nm doublet which can only be detected in spectra of the highest quality. But even with the best spectra, the derived ⁶Li abundance can only be as good as the synthetic spectra used for their interpretation. It is now known that 3D non‐LTE modelling of the lithium spectral line profiles is necessary to account properly for the intrinsic line asymmetry, which is produced by convective flows in the atmospheres of cool stars, and can mimic the presence of ⁶Li. We also discuss briefly the case of the carbon isotopic ratio in metal‐poor stars, and provide a new determination of the nickel isotopic ratios in the solar atmosphere. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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).     

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.     

A New Solar Spectrum from 656 to 3088 nm

The solar spectrum is a key parameter for different scientific disciplines such as solar physics, climate research, and atmospheric physics. The SOLar SPECtrometer (SOLSPEC) instrument of the Solar Monitoring Observatory (SOLAR) payload onboard the International Space Station (ISS) has been built to measure the solar spectral irradiance (SSI) from 165 to 3088 nm with high accuracy. To cover the full wavelength range, three double-monochromators with concave gratings are used. We present here a thorough analysis of the data from the third channel/double-monochromator, which covers the spectral range between 656 and 3088 nm. A new reference solar spectrum is therefore obtained in this mainly infrared wavelength range (656 to 3088 nm); it uses an absolute preflight calibration performed with the blackbody of the Physikalisch-Technische Bundesanstalt (PTB). An improved correction of temperature effects is also applied to the measurements using in-flight housekeeping temperature data of the instrument. The new solar spectrum (SOLAR–IR) is in good agreement with the ATmospheric Laboratory for Applications and Science (ATLAS?3) reference solar spectrum from 656 nm to about 1600 nm. However, above 1600 nm, it agrees better with solar reconstruction models than with spacecraft measurements. The new SOLAR/SOLSPEC measurement of solar spectral irradiance at about 1600 nm, corresponding to the minimum opacity of the solar photosphere, is 248.08 ± 4.98 mW?m^?2?nm^?1 (1?\(\sigma\)), which is higher than recent ground-based evaluations.     

Elemental abundances in the Milky Way stellar disk(s), bulge,and halo

We present a review of elemental abundances in the Milky Way stellar disk, bulge, and halo with a focus on data derived from high-resolution stellar spectra. These data are fundamental in disentangling the formation history and subsequent evolution of the Milky Way. Information from such data is still limited and confined to narrowly defined stellar samples. The astrometric Gaia satellite will soon be launched by the European Space Agency. Its final data set will revolutionize information on the motions of a billion stars in the Milky Way. This will be complemented by several ground-based observational campaigns, in particular spectroscopic follow-up to study elemental abundances in the stars in detail. Our review shows the very rich and intriguing picture built from rather small and local samples. The Gaia data deserve to be complemented by data of the same high quality that have been collected for the solar neighborhood.     

A Multi-Wavelength Simultaneous Study of the Composition of the Halley Family Comet 8P/Tuttle

We report on simultaneous optical and infrared observations of the Halley Family comet 8P/Tuttle performed with the ESO Very Large Telescope. Such multi-wavelength and coordinated observations are a good example of what can be done to support space missions. From high resolution optical spectroscopy of the CN (0,0) 388 nm and NH₂ (0,9,0) 610 nm bands using UVES at UT2 we determined ¹²C/¹³C = 90 ± 10 and ¹⁴N/¹⁵N = 150 ± 20 in CN and we derived a nuclear spin temperature of NH₃ of 29 ± 1 K. These values are similar to those found in Oort-Cloud and Jupiter Family comets. From low resolution long slit spectroscopy with FORS1 at UT2 we determined the CN, C₃ and C₂ production rates and the parent and daughter scale lengths up to 5.2 10⁵ km tailward. From high resolution IR spectroscopy with CRIRES at UT1 we measured simultaneously the production rates and mixing ratios of H₂O, HCN, C₂H₂, CH₄, C₂H₆, and CH₃OH.     

Investigation of galaxies in the second byurakan sky survey. I. Spectroscopic observations in a field at 15h30m, +59°

The results of spectroscopic observations of 30 galaxies from the Second Byurakan Spectral Sky Survey in a field of 16 square degrees centered at coordinates a = 15^h30^m, δ = + 53°are presented. The spectra were obtained on the 6-m telescope of the Special Astrophysical Observatory, Russian Academy of Sciences, in the wavelength range of 4000–7400 Å with 15 Å resolution. Emission lines are present in the spectra of most of the galaxies. Redshifts and absolute stellar magnitudes were determined for all the galaxies.     

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     

The nuclear structure and associated electron capture rates on odd-Z nucleus 51V in stellar matter

The Gamow-Teller strength distribution function, B(GT), for the odd Z parent ⁵¹V, N?Z=5, up to 30 MeV of excitation energy in the daughter ⁵¹Ti is calculated in the domain of proton-neutron Quasiparticle Random Phase Approximation (pn-QRPA) theory. The pn-QRPA results are compared against other theoretical calculations, (n, p) and high resolution (d, ²He) reaction experiments. For the case of (d, ²He) reaction the calibration was performed for 0≤E _j ≤5.0 MeV, where the authors stressed that within this excitation energy range the ΔL=0 transition strength can be extracted with high accuracy for ⁵¹V. Within this energy range the current pn-QRPA total B(GT) strength 0.79 is in good agreement with the (d, ²He) experiment’s total strength of 0.9±0.1. The pn-QRPA calculated Gamow-Teller centroid at 4.20 MeV in daughter ⁵¹Ti is also in good agreement with high resolution (d, ²He) experiment which placed the Gamow-Teller centroid at 4.1±0.4 MeV in daughter ⁵¹Ti. The low energy detailed Gamow-Teller structure and Gamow-Teller centroid play a sumptuous role in associated weak decay rates and consequently affect the stellar dynamics. The stellar weak rates are sensitive to the location and structure of these low-lying states in daughter ⁵¹Ti. The calculated electron capture rates on ⁵¹V in stellar matter are also in good agreement with the large scale shell model rates.     

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.     

Theoretical models of planetary nebulae

Theoretical models are calculated for 15 planetary nebulae of medium-to-high excitation, following procedures previously described. Initial stellar energy distributions are adopted from Cassinelli (1971), but are subsequently modified to obtain the best representation of optical spectra for the selected objects. Other adjustable parameters include the stellar radius,R ^(*), the nebular density,N _H, the truncation radius,r _c, for the nebular shell, and the chemical composition. Excitationsensitive ratios are usually well-represented as are the actual observed intensities of spectral lines. Forbidden lines arising from 3p ³ configurations, e.g., those of [SII], [ArIV], and [ClIII] offer difficulties. For this sample of nebulae, the mean abundances seem to agree well with those found in an earlier study where the models were used as interpolation devices (Aller, 1978). Our objective is not to use the models to derive abundances explicitly, but rather to use them to find ionization correction factors. Some cautions and limitations of this procedure are described.     

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.     

Uncertainties of Synthetic Population Models

This paper investigates the uncertainties in the synthetic integrated colours of simple stellar populations, currently the most popular method of estimating the ages of unresolved stellar systems. Three types of uncertainties studied here originate from the stellar models, the population synthesis techniques, and from the stellar spectral libraries. Despite some skepticism, synthetic colours appear to be reliable age indicators as long as they are used for select age ranges. This revised version was published online in September 2006 with corrections to the Cover Date.     

Second‐order error propagation in the Mueller matrix of a spectropolarimeter

Accurate measurements of Stokes IQUV in spectral lines is required for precise reconstruction of stellar magnetic field geometries with Zeeman‐Dopper imaging. Spectral Zeeman features are intrinsically weak and subjected to a number of instrumental uncertainties. The aim of this work is to study the details of the instrumental uncertainties in the Stokes IQUV measurements in spectral lines and ways of their reduction. We make a practical comparison of the polarimetric performances of two high‐resolution échelle spectropolarimeters, namely SOFIN at the NOT, and HARPS at ESO. We show the residual spectra for both instruments to characterize the cross‐talk between the observed Stokes parameters. We employ a self‐calibrating least‐squares fit to eliminate some of the polarization uncertainties to derive the full Stokes vector from stellar spectra (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A high-resolution stellar library for evolutionary population synthesis

We present a library of 1654 high-resolution stellar spectra, with a sampling of 0.3 Å and covering the wavelength range from 3000 to 7000 Å. The library was computed with the latest improvements in stellar atmospheres, incorporating non-local thermodynamic equilibrium (non-LTE) line-blanketed models for hot, massive  ( T _eff≥ 27 500 K)  and line-blanketed models for cool  (3000 ≤ T _eff≤ 4500 K)  stars. The total coverage of the grid is  3000 K ≤ T _eff≤ 55 000 K  and  −0.5 ≤ log g ≤ 5.5  , for four chemical abundance values: twice solar, solar, half solar and 1/10 solar. Evolutionary synthesis models using this library are presented in a companion paper. We tested the general behaviour of the library by calculating and comparing equivalent widths of numerous H and He  i lines, and some of the commonly used metallic indices. We also compared the library with the empirical libraries STELIB and Indo-US. The full set of the synthetic stellar spectra is available from our websites ( http://www.iaa.csic.es/~rosa and http://www.astro.iag.usp.br/~lucimara/library.htm ).     

Evaluating the Uncertainties in the Electron Temperature and Radial Speed Measurements Using White Light Corona Eclipse Observations

We examine the uncertainties in two plasma parameters from their true values in a simulated asymmetric corona. We use the Corona Heliosphere (CORHEL) and Magnetohydrodynamics Around the Sphere (MAS) models in the Community Coordinated Modeling Center (CCMC) to investigate the differences between an assumed symmetric corona and a more realistic, asymmetric one. We were able to predict the electron temperatures and electron bulk flow speeds to within ±?0.5 MK and ±?100 km?s^?1, respectively, over coronal heights up to 5.0 R_⊙ from Sun center. We believe that this technique could be incorporated in next-generation white-light coronagraphs to determine these electron plasma parameters in the low solar corona. We have conducted experiments in the past during total solar eclipses to measure the thermal electron temperature and the electron bulk flow speed in the radial direction in the low solar corona. These measurements were made at different altitudes and latitudes in the low solar corona by measuring the shape of the K-coronal spectra between 350 nm and 450 nm and two brightness ratios through filters centered at 385.0 nm/410.0 nm and 398.7 nm/423.3 nm with a bandwidth of ≈?4 nm. Based on symmetric coronal models used for these measurements, the two measured plasma parameters were expected to represent those values at the points where the lines of sight intersected the plane of the solar limb.     

Automated estimation of stellar fundamental parameters from low resolution spectra: the PLS method

PLS (Partial Least Squares regression) is introduced into an automatic esti-mation of fundamental stellar spectral parameters. It extracts the most correlative spec-tral component to the parameters (Teff, log g and [Fe/H]), and sets up a linear regres-sion function from spectra to the corresponding parameters. Considering the properties of stellar spectra and the PLS algorithm, we present a piecewise PLS regression method for estimation of stellar parameters, which is composed of one PLS model for Teff, and seven PLS models for log g and [Fe/H] estimation. Its performance is investigated by large experiments on flux calibrated spectra and continuum normalized spectra at dif-ferent signal-to-noise ratios (SNRs) and resolutions. The results show that the piecewise PLS method is robust for spectra at the medium resolution of 0.23 nm. For low resolu-tion 0.5 nm and 1 nm spectra, it achieves competitive results at higher SNR. Experiments using ELODIE spectra of 0.23 nm resolution illustrate that our piecewise PLS models trained with MILES spectra are efficient for O ～ G stars: for flux calibrated spectra, the systematic offsets are 3.8%, 0.14dex, and -0.09 dex for Teff, log g and [Fe/H], with error scatters of 5.2%, 0.44 dex and 0.38 dex, respectively; for continuum normalized spectra, the systematic offsets are 3.8%, 0.12dex, and -0.13 dex for Teff, log g and [Fe/H], with error scatters of 5.2%, 0.49 dex and 0.41 dex, respectively. The PLS method is rapid, easy to use and does not rely as strongly on the tightness of a parameter grid of templates to reach high precision as Artificial Neural Networks or minimum distance methods do.