ELODIE: THE STELLAR LIBRARY

Philippe Prugniel and Caroline Soubiran

 

HOMEDownloadsDescriptionUsageLinks

              CONTENT  REDUCTION  FITS KEYWORDS  

The ELODIE.3.1 library is a subset of  The ELODIE archive distributing the original data, in instrumental flux, of all the observations obtained at the OHP 1.93m telescope with the corresponding spectrograph, since 1994.



      

ELODIE.3.1 contains 1962 spectra of 1388 stars with a large coverage in Teff, log g, and [Fe/H]. Fig. 1 presents the distribution of the stars in the HR diagram, with three metallicity bins differentiated with symbols of different colors. Fig. 2 presents the distribution in the Teff - [Fe/H] plane for dwarfs and giants. Stellar parameters and useful data have been searched in the literature or determined for all stars. They are described here, and can be retreived in this table (README).
 
 

   Figure 1





The log file of the modifications in the different versions the ELODIE library is here.


Since the first version of the library, the reduction procedure described in 2001A&A...369.1048P has been improved in several aspects but the general philosophy remains the same. The basic steps in the reduction are to

(1) correct the different orders in the spectra for the blaze effect and connect the orders together

(2) mask the telluric lines and spikes due to cosmic rays : in the new version, the mask of the spikes due to cosmic rays has been suppressed because it was found to alter real features. Because of it, in the previous version of the library Lick indices measured on our spectra presented a slope compared to measurements on Jones spectra, the strongest features were found weaker than real (thanks to G. Worthey who stressed our attention on this important problem).

(3) make the flux calibration

(4) build an interpolator which in turn is used to generate the grid that feeds the population synthesis program (PEGASE.HR).

The strengthening of the spectra have been improved: artifacts near Balmer's lines were corrected. The wavelength interval has been extended by adding 5 orders in the blue. It was 410-680 nm, it is now 400-680. The correction for telluric lines has been improved. We corrected some small bugs in the FITS keywords (in particular an unfortunate bad rounding of CRVAL1).

Flux calibration:
The main improvement in the flux calibration is due to the increase of the number of external and internal comparisons. We made 2 observing runs to have more intercomparisons. We have also improved the color equations. Finally, we have improved the determination of the Galactic extinction. It is determined from either Klare and Neckel, Hakkila et al. or Chen et al. after a consistency test between these various 3D model. Correction of the galactic extinction is primordial for population
synthesis and quite critical because there exists a correlation in the library between the spectral type and the extinction (both hot
stars and red giants are often in the disk and have high extinction).

Interpolator:
The main goal of project was to make a stellar library for synthesis of stellar population. This require to provide a spectrum at any desired point of the parameter space, and actually the spectral systhesis tools we are using, PEGASE.HR, uses a grid of spectra with a defined mesh.

The interpolator produces a spectrum at given Teff, log g and [Fe/H] by evaluating polynomials in a way quite similar to the Lick fitting functions (Worthey 1994). The polynomials include with up to 21 terms and up to the power 4 in each of the parameters with a selection of cross terms. Each wavelength element is represented with an independent polynomial but the functional form is the same for the whole spectrum. As well as for the Lick fitting functions the parameter space is divided in a few overlapping regions where different polynomials are fitted. The overlap between these regions is used to give a smooth transition between the solutions obtained in each region. In the overlap the fluxes are a linear interpolation between the polynomials of each region. The terms choosen, and the validity range of the polynomials were chosen empirically by trial and errors.

The quality of the interpolator can be checked by comparing the interpolated fluxes with the observed fluxes for the spectra in the library. So, at each iteration we searched these residuals for discrepancies and systematic effects and we modified the form of the polynomials and the validity regions in order to improve the model. The solution adopted in the present version reflect the intrinsic variation of the lines and the actual distribution of the stars in the space of parameters.

Because the library contains a large number of F-K stars, we have also weighted the spectra in order to re-equilibrate the distributions. The polynomials were fitted by least-square with a kappa-sigma rejection of outliers (un-corrected bad pixels in some spectra). In the new version we have much more hot stars (>10000K) and cool stars (<4500K) than in previous version and the interpolator is therefore considerably improved with respect to the previous version.

Absolute and internal interpolators: The interpolator was first computed using the stars with known atmospheric parameters, giving each stars a weight function of the quality of these determination. This version is called the "absolute interpolator". Inverting the interpolator for a spectrum corresponding to a star with unknown parameters allows to determine these parameters. For the stars used to build the interpolator, the comparison between the parameters resulting for inversion and those from our catalogue allow to detect errors on the input parameters of some stars (or errors on the identification of some stars at the telescope). In some cases these errors could be corrected and the interpolator were re-computed. This comparison also allowed to assess the weaknesses of the interpolator and drived its improvement (see above) and finally the statistics on this comparison reflect the errors on the parameters in our catalogue plus the errors due to the processing and modelling.

We then determined the parameters for all the stars using the absolute interpolator and used them to compute a new interpolator: the internal interpolator. This internal interpolator may suffer from some biases introduced by the absolute interpolator, but it has smaller ramdom errors because the errors on the input atmospheric parameters are decreased. It is also more stable because it is computed on more spectra, which is particularly sensitive in the less populated regions of the space of parameters. We used the internal interpolator to compute the internal determination of the atmospheric parameters listed in table 1.

Physical flux interpolator, correction for Galactic extinction: The interpolator was yet computed using the fluxes normalized to the continuum, because it is better suited to the determination of atmospheric parameters. The reason is that our flux calibration is not absolutely precise and since the information on the atmospheric parameters is essentially contained in the lines, we wanted to avoid to bias the atmospheric parameters with errors in the flux calibration. But to make the population synthesis we need to generate spectra calibrated in fluxes as emitted by the stars (ie. above the earth atmosphere and corrected for interstellar extinction). The next step was therefore to correct galactic extinction. We used for that the Schild (1977) extinction law (the choice of the extinction law is not critical in this region of wavelengths) scaled with the E(B-V) excess listed in our catalogue. We build an interpolator for source flux spectra using the internal determination of the atmospheric parameters. Again it was possible to compare the results of the interpolation with all the input spectra. This allowed to detect spectra with discrepant SED, either due to flux calibration, to erroneous value of the Galactic extinction. The most discrepant case were closely examined, searching the literature for possible explanations or hints for errors on the parameters and checking the consistency of the estimates of the Galactic extinction and the various photometric measurements for the objects. Errors in the catalogues were detected and corrected and in other cases we substituted the catalogued value of the Galactic extinction with a value estimated from the photometry (and atmospheric parameters).

A final test on the quality of the source flux internal interpolator was to compare the colour of the stars in different regions of the HR diagram with the colors measured on the spectra (Table xxx). The result of the comparison is quite accpetable but we still see a systematic effect that we believe is due to the color equations used to convert between the colors measured on Elodie spectra, the Johnson colors and the Tycho colors. We would like to solve this discrepancy in the near future.

Consistency and external comparisons: Measurements of Lick indices and determination of the atmospheric parameters are described in the article and the results are presented in.


Description of "calibration" files
These files, whose archive number starts with "99", are primarly intended to be used by the specialized reduction programs and by the pipeline processing. Their description may not be accurate enough for a standalone usage, and for such applications the users are invited to contact us for additional help.

LL_ELODIE/99001 Blaze function used as reference This is a "s2d" file produced by the ELODIE reduction program (TACOS). "s2d" means: spectra extracted by TACOS, each line is one order. The Blaze function is obtained from a Tungsten lamp spectrum. Each order is smoothed using a degree 17 polynomial and is normalized to its maximum value.

LL_ELODIE/99002 Correction to the TACOS blaze (LL_ELODIE/99001) Dividing a spectrum by a TACOS blaze produces a stair case result (consecutive orders do not connect because of the independent nomalizations). In order to connect the orders and to reduce the instrumental response to a "typical" stellar SED, we computed this correction using 3 metal-deficient stars observed during the same observing run. The original "s2d" spectra are divided by LL_ELODIE/99001 and the continuum of each order is modelled with a degree 2 polynomial.

LL_ELODIE/99003 Pseudo-continuum mask It is computed after the whole archive has been reduced. In each continuum-normalized spectrum, a simple threshold is used to diagnostic if a given wavelength point belong to the pseudo-continuum. We choosed 0.97 as threshold. The corresonding wavelength point is given the value 1 or 0 if it belongs respectively to the continuum or not. LL_ELODIE/99003 is the average of this mask over the whole archive. In some sense it indicates the probability for each wavelength point to be a continuum point. It is used as a weight for iterrating the computation of LL_ELODIE/99002 and of the whole archive. (determining the continuum is required to fit and correct variations of the blaze, see the text of the article).

LL_ELODIE/99004 Primary flux calibration. The first axis is the wavelength direction The first line is the instrumental response, The second the airmass dependence and the third the "haze" dependence (see text). LL_ELODIE/99004 is obtained after a comparison between the instrumental fluxes and templates from the datasets L1985BURN L1998GLUS and L1997ALEK

LL_ELODIE/99005 Smoothed flux calibration file Each line of LL_ELODIE/99004 is modeled by a degree 17 polynomial in order to discard the wavy residuals due to an approximate modeling of the spectral resolution of the templates

LL_ELODIE/99006 Secondary flux calibration. LL_ELODIE/99005 with correction after pairwyse comparisons

LL_ELODIE/99007 Corrected flux calibration file The spectra calibrated in physical flux with LL_ELODIE/99005 are compared with spectra from L1984JACO L1996SERO L96111KP1. The details of these comparison are given in LL_ELODIE/99007. The mean comparison is smoothed and used to compute a correction to the instrumental response. LL_ELODIE/99006 serves to compute the final flux calibration relations associated to each spectrum in the archive.

LL_ELODIE/99008 Comparison with other templates. The first axis is the wavelength direction. The 9 first lines contains respectively the comparisons with (1) L96111KP1 and (2) L96111KP2 (Jone's library, respectively red and blue bands), (3) L1996SERO (OHP observations), (4) L1996SERO (CFH observations), (5) L1984JACO (Jacoby's library) (6) L1987KIEH (Kiehling 1987) (7) L1994DANK (Danks \& Dennefeld 1994) (8) L1983GUNN (Gunn & Stryker 1983 ) and (9) L1984NSOA (Solar flux atlas). The 10th line combines the first 9 lines. And the 11th is the final correction obtained by smoothing the 10th.

LL_ELODIE/99010 Absolute TGM model The first axis is the wavelength direction. The second axis contains the 20 parameters of the TGM model fitted to the spectra for which the atmosphric parameters are known. The functions of the models are: 1 Cst 2 T 3 Z 4 G 5 T**2 6 T**3 7 T**4 8 T*Z 9 T*G 10 G**2 11 Z**2 12 T**2*G 13 T**2*Z 14 G**3 15 Z**3 16 T*G**2 17 T*Z**2 18 10**T 19 10**(T*2) 20 G*Z If Teff is the effective temperature, log(G) the decimal logarithm of the surface gravity, T and G are defined as: T=log(Teff)-3.7 and G=log(G)-4.44. Z is [Fe/H]. Each wavelength point is fittedt to the 20 functions above and each line in 99010 are the coefficients of these functions. This is a preliminary model and the final version will be described in the paper presenting TGMET.

LL_ELODIE/99011 Internal TGM model Same format as LL_ELODIE/99009 but the coefficient are fitted on the whole archive using the atmospheric parameters obtained from a fit to LL_ELODIE/99010. 


Description of FITS keywords specific to these datasets

TEFF 	[K]	Teff from literature
Q_TEFF 		Quality flag for Teff 
Q_TEFF=1:Poor determination -> Q_Teff=4: Excellent (see Sect. 2)
Q_TEFF=0:means that we have estimated Teff from the B-V colour (Tycho2 catalogue),assuming the empirical colour-temperature relation for a main sequence star and neglecting the interstellar extinction
Q_Teff=-1: Internal determination of the effective temperature
LOGG 	[log cm/s**2)] 	log g from literature

Q_LOGG Quality flag for LOGG
Q_LOGG=1 means that the determination of LOGG is taken from the literature.
Q_LOGG=0 means that LOGG was converted from the V absolute magnitude from Hiparcos and Teff, using a bolometric correction valid for a main sequence start and an empirical mass-to-light relation.
Q_LOGG=-1 Internal determination of the surface gravity

FE_H [Sun] [Fe/H] from literature

Q_FE_H Quality flag for [Fe/H]
Q_Fe/H=1: Poor determination -> Q_Fe/H=4: Excellent
Q_Fe/H=-1: Internal determination of [Fe/H]

Mv [mag] Mv from Hipparcos Q_MV Quality flag for Mv: 4=excellent

S/N S/N-per-pixel of original spectrum at 555 nm

VR [km/s] heliocentric radial velocity

FWHM [km/s] FWHM line broadening

Q_VR comment on Vr determination

SPTYPE spectral type from INCA

VMAG [mag] V-magnitude

E_VMAG [mag] error on V-magnitude

B-V [mag] B-V Johnson

R_PHOT Source of photometry: TYCHO2, INCA
webmaster