---
author:
- 'P. Westera'
- 'T. Lejeune'
- 'R. Buser'
- 'F. Cuisinier'
- 'G. Bruzual A.'
date: 'Received &lt;date&gt;; accepted &lt;date&gt;'
subtitle: 'III. Metallicity calibration'
title: 'A standard stellar library for evolutionary synthesis.'
---

Introduction
============

As present grids of theoretical spectral energy distributions (SEDs) suffer from intrinsic inhomogeneities and incompleteness and show large systematic discrepancies with empirical calibrations due to unavailable molecular opacity [see @paperi hereafter Paper I], we have undertaken the construction of a comprehensive combined library of realistic stellar flux distributions. Empirical ${T_{\rm eff}}$-colour relations in $UBVRIJHKL$ photometry are used to adjust the spectra using an algorithm developed by Cuisinier et al. [see @buser Paper I].\
The current state of the art is the following: The semi-empirical BaSeL (Basel Stellar Library) 2.2 SED library [@paperii hereafter Paper II] has been widely used successfully in different areas [@ortolani; @kauffmann; @westera; @gonzalez; @origlia; @leitherer; @bruz; @brocato; @maraston; @barmby; @liu; @oblak; @lotz; @molla; @nikolaev; @kong; @yoon; @marleau; @vacca; @lastennet; @kotilainen; @fricke]. However, as a result of being calibrated from solar metallicity data only, BaSeL 2.2 still has its weaknesses at low metallicities (${\rm [Fe/H]} < -1$), especially in the ultraviolet ($U-B$) and the infrared ($V-K$, $J-H$, $H-K$, $J-K$, $K-L$). In these colours and at these lower metallicities, synthetical globular cluster CMDs appear too blue using the BaSeL 2.2 semi-empirical library. Furthermore, the transition between dwarfs and giants produces a discontinuity in some colours (for more details, see Papers I and II as well as T. Lejeune’s thesis (1997)).\
The purpose of the present paper is to remove the weaknesses of the BaSeL 2.2 semi-empirical library by extending the colour-calibration to low metallicities, and to create a library that reproduces empirical colour-temperature relations and globular cluster CMDs (using existing grids of isochrones) at all metallicities.\
The outline of the paper is the following: In section 2, we discuss the compilation and properties of the calibration data, i. e. globular cluster CMDs and empirical colour-${T_{\rm eff}}$ relations. In section 3, we briefly describe the calibration algorithm and the changes made relative to the previous algorithm (described in Paper I). As a result, we present the BaSeL 3.1 ”WLBC 99” SED library, which is able to reproduce empirical colour-temperature relations in all ($UBVRIJHKL$) colours, and is therefore projected to be a powerful tool for studies of individual stars. Unfortunately, the library doesn’t provide any improvements in the representation of globular cluster CMDs, and it even proved impossible to provide a library that satisfies both requirements at the same time. As a pragmatic solution, we also produce an application-oriented library, the BaSeL ”Padova 2000” library, which, if used along with the Padova 2000 isochrones, successfully reproduces globular cluster CMDs at all levels of metallicity; this version of the library is presented and discussed in section 4. The conclusions are summarised in section 5, where an outlook on future work and on a first application can also be found.\
A more detailed description of this work is given in @diss.

Compilation and discussion of calibration data {#compilation}
==============================================

![image](MS1686f1.eps){width="\textwidth"}

![image](MS1686f2.eps){width="\textwidth"}

  cluster, [\[Fe/H\]]{}, $(m-M)_{0}$, $E_{B-V}$   colour bands   source                  range (in $m_{V}$)
  ----------------------------------------------- -------------- ----------------------- -----------------------------------
  47 Tuc (NGC 104), –0.70, 13.14, 0.04            $B-V$          @wysocka                $\stackrel{_<}{_\sim} 19^{m}$
                                                  $B-V$          @hesser                 $\sim 24^{m}$ - $\sim 19^{m}$
                                                  $V-I$          @barbuy                 $\stackrel{_<}{_\sim} 13.5^{m}$
                                                  $V-I$          @wysocka                $\sim 21^{m}$ - $\sim 13.5^{m}$
                                                  $U-B$          @gcpd$^{a}$             $\sim 16^{m}$ - $\sim 12^{m}$
                                                  $V-R^{b}$      @cathey                 $\sim 16^{m}$ - $\sim 12^{m}$
                                                  $VJK$          @montegriffo            $\stackrel{_<}{_\sim} 17^{m}$
                                                  $H-K$          @frogel                 $\stackrel{_<}{_\sim} 14.5^{m}$
  M 5 (NGC 5904), –1.11, 14.32, 0.03              $B-V$          @sandquist              $\stackrel{_<}{_\sim} 22^{m}$
                                                  $B-V$          @richer                 $\sim 24^{m}$ - $\sim 22^{m}$
                                                  $V-I$          @sandquist              $\stackrel{_<}{_\sim} 22^{m}$
                                                  $U-B$          @drissen,               $\stackrel{_<}{_\sim} 19.5^{m}$
                                                                 @vonbraun               
                                                  $V-R$          @vonbraun               $\stackrel{_<}{_\sim} 17.5^{m}$
                                                  $J-K$          @persson                $\stackrel{_<}{_\sim} 17.5^{m}$
  M 3 (NGC 5272), –1.34, 15.01, 0.01              $B-V$          @ferraro                $\stackrel{_<}{_\sim} 22^{m}$
                                                  $V-I$          @ferraro                $\stackrel{_<}{_\sim} 16^{m}$
                                                  $V-I$          @bolte                  $\sim 23^{m}$ - $\sim 16^{m}$
                                                  $V-I$          @marconi                $\sim 27^{m}$ - $\sim 23^{m}$
                                                  $U-B$          @gcpd$^{c}$             $\stackrel{_<}{_\sim} 21.5^{m}$
                                                  $V-R^{b}$      Arribas and             $\stackrel{_<}{_\sim} 14.5^{m}$
                                                                 Martínez-Roger (1987)   
                                                  $V-K$          @frogel                 $\stackrel{_<}{_\sim} 15^{m}$
                                                  $J-K$          @kim                    $\stackrel{_<}{_\sim} 15^{m}$
  NGC 6397, –1.82, 11.75, 0.18                    $B-V$          @kaluzny                $\stackrel{_<}{_\sim} 18.5^{m}$
                                                  $B-V$          @alvarado               $\sim 21.5^{m}$ - $\sim 18.5^{m}$
                                                  $U-B$          @alvarado               $\sim 20^{m}$ - $\sim 16.5^{m}$
                                                  $U-B$          @gcpd$^{d}$             $\stackrel{_<}{_\sim} 16.5^{m}$
                                                  $V-I$          @alvarado               $\stackrel{_<}{_\sim} 12.5^{m}$
                                                  $V-I$          @alvarado               $\sim 16.5^{m}$ - $\sim 12.5^{m}$
                                                  $V-I$          @king                   $\sim 27^{m}$ - $\sim 16.5^{m}$
                                                  $R-I$          @alcaino                $\stackrel{_<}{_\sim} 16.25^{m}$
  M 92 (NGC 6341), –2.16, 14.53, 0.02             $B-V$          @sandage                $\stackrel{_<}{_\sim} 21.5^{m}$
                                                  $B-V$          @stetson                $\sim 24^{m}$ - $\sim 21.5^{m}$
                                                  $V-I$          @vonbraun               $\stackrel{_<}{_\sim} 14^{m}$
                                                  $V-I$          @bolte                  $\sim 20^{m}$ - $\sim 14^{m}$
                                                  $V-I$          @piotto                 $\sim 28^{m}$ - $\sim 20^{m}$
                                                  $U-B$          @gcpd$^{e}$             $\stackrel{_<}{_\sim} 20.5^{m}$
                                                  $V-K$          @frogel                 $\stackrel{_<}{_\sim} 15^{m}$
                                                  $JHK^{f}$      @cohen                  $\stackrel{_<}{_\sim} 15^{m}$

[$^{a}$@alcaino [@evans; @norris; @lee; @demarque; @hartwick; @cannon; @menzies; @eggen]\
$^{b}$transformed to Johnson-Cousins using @bessellb\
$^{c}$@johnson [@sandageb; @sandage]\
$^{d}$@woolley [@newell; @cannon; @alcain; @vandenbergh]\
$^{e}$@walker [@sandageb; @sandage; @eggen]\
$^{f}$transformed to Johnson using @besselld]{}

In order to carry out the colour-calibration in a metallicity-dependent way, we collected empirical photometry from the literature of the best-studied Galactic globular clusters, spanning wide ranges in metallicity and evolutionary stages. The clusters used were 47 Tuc, M5, M3, NGC6397 and M92, with the respective ${\rm [Fe/H]}$-values of $-0.70$, $-1.11$, $-1.34$, $-1.82$, and $-2.16$ [@carretta]. To reach the best possible coverage of the multi-colour ($UBVRIJHKL$) - magnitude - diagram of each cluster, the photometry was combined from numerous sources summarised in table \[clustasources\]. Where necessary, the photometry was brought onto the Johnson- Cousins system using colour-colour relations by Bessell and collaborators [@bessella; @bessellb; @bessellc; @besselld]. Where no fiducial lines were given, they were drawn by eye. In some colours (mainly $B-V$ and $V-I$), different sources were found for the same part of the CMD of the same cluster. The fiducials derived from or given by such sources usually didn’t differ from each other by more than a few hundredths of a magnitude.\
In order to put these colour-magnitude diagrams onto the intrinsic system, the $E_{B-V}$ and $(m$-$M)_{V}$ values from the Harris online catalog of globular cluster parameters [@harris] were used. The adopted distance scale is based on the luminosity level of the horizontal branch (HB), where the absolute magnitude of the HB was derived from an empirical, metallicity-dependent relation ($M_{V}(HB) = 0.15 {\rm [Fe/H]} + 0.80$). Harris estimates that the absolute uncertainty of the predicted $M_{V}(HB)$ is of the order of 0.1-$0.2^{m}$, in the non-extreme cases of the used clusters more likely $0.1^{m}$. For the extinction coefficients, he gives an uncertainty of $\Delta E_{B-V} = 0.1 E_{B-V}$. This shouldn’t pose any problems, because we used clusters with extinction coefficients of only a few hundredths of a magnitude (only NGC 6397 has an $E_{B-V}$ of 0.18). The result are the fiducials shown in figs. \[CMDsa\] and \[CMDsb\] (only the main sequence, the subgiant branch and the red giant branch are shown because the other parts of the CMD weren’t used for the calibration).\
Two properties of the empirical CMDs are very striking. First, there are large gaps in the data. Measurements down to a few magnitudes below the main sequence turn-off exist only in $BVI$. In the far infrared, hardly any data at all are available. In order to create a sufficient photometric basis for population synthesis, more work clearly needs to be done on the observational side, in order to improve results based on either observational or synthetic spectra.\
The second striking point about the CMDs shown in figs. \[CMDsa\] and \[CMDsb\] is the fact that their behaviour with increasing metallicity isn’t as systematic or monotonous as one would expect. The colours of the RGB at a given ($V$) magnitude can differ up to $\sim 0.1^{m}$ from any trend with ${\rm [Fe/H]}$. The scatter around trends with metallicity is not of the same amplitude for all colours. While in $U-B$, $B-V$, and $R-I$ it is the most pronounced, $V-I$ behaves much more systematically with ${\rm [Fe/H]}$, as is also seen in homogeneous databases like the one from @saviane[^1]. Part of the scatter (a few hundredths of a mag, see above) can be explained by the heterogeneity of the database which however cannot account for the entire scatter in those colours where it is observed. Obviously, there are parameters apart from age and metallicity that govern the appearance of the CMD. It is therefore more than amazing how in the literature, one often finds perfect agreement between synthetic isochrones and empirical CMDs for whole sets of globular clusters at the same time.\
From these CMDs, combined multi-colour ($UBVRIJHKL$) - ${\rm [Fe/H]}$ - ${T_{\rm eff}}$ - $\log g$ relations were synthesised, using the ${T_{\rm eff}}$ - ($V-K$) relation from the BaSeL 2.2 library. This relation incorporates @ridgway complemented with the differential properties in ${T_{\rm eff}}$ for solar metallicity of the original (uncalibrated) grid, as the latter is the only relation which covers the entire range in temperature. This resulting relation was used for all metallicities, because the $V-K$ colour is expected to be metallicity insensitive [@vonbraun; @martinez]. Finally, $\log g$ values were added to the calibration files, using empirical ${T_{\rm eff}}$ - $\log g$ relations for red giants from Cohen, Frogel and Persson [@cohen; @frogel; @persson ${\rm [Fe/H]}$-dependent] and for dwarfs from @angelov [${\rm [Fe/H]}$-independent].

  Luminosity class     ${\rm [Fe/H]}$ Colours          temperature range \[K\]
  ------------------ ---------------- ---------------- -------------------------
  giants                         –0.5 $UBVRIJHK^{a}$   3350 - 6000
                                 –1.0 $UBVRIJHK$       3750 - 6000
                                 –1.5 $UBVRIJHK$       3750 - 6000
                                 –2.0 $UBVRIJHK$       4000 - 6000
                                  all $K-L$            -
  dwarfs                          all $UBVI^{b}$       4500 - 10000
                                  all $RJHKL^{c}$      -

[$^{a}UBVRIJHK$: $U-B$, $B-V$, $V-I$, $R-I$, $V-K$, $J-H$, $H-K$, $J-K$\
$^{b}UBVI$: $U-B$, $B-V$, $V-I$\
$^{c}RJHKL$: $R-I$, $V-K$, $J-H$, $H-K$, $J-K$, $K-L$]{}

We could only produce colour-temperature relations in the ranges specified in Table \[ranges\]. For giants, the ranges include the entire RGB (except in $K-L$). Outside these ranges, we used the differential properties with regard to ${T_{\rm eff}}$ of the semi-empirical (BaSeL 2.2) grid, as there was no other information available (neither empirical nor theoretical) about the behaviour of these colours for these temperatures. In $K-L$ for all stars, and for dwarfs in $R-I$, $V-K$, $J-H$, $H-K$, and $J-K$, where no observed data were available, but which are needed to complete the calibration files, we had no choice but to adopt the ${T_{\rm eff}}$ - $K-L$ relation directly from the BaSeL 2.2 library.\
This set (i. e. metallicity-dependent colour - temperature - $\log g$) was complemented with the solar relations synthesised by Lejeune et al. to calibrate the BaSeL 2.2 library (see Paper II).

![image](MS1686f3.eps){width="\textwidth"}

![image](MS1686f4.eps){width="\textwidth"}

The so-derived calibration files (two for each metallicity (${\rm [Fe/H]} = 0, -0.5, -1, -1.5, -2$), one for (red) giants and one for dwarfs) are available by public ftp from the university of Basel.

Metallicity-dependent library calibration, the BaSeL 3.1 ”WLBC 99” library
==========================================================================

These calibration files were then used to colour-calibrate the hybrid library of theoretical SEDs (consisting of the Kurucz-, the Allard and Hauschildt-, and the Scholz-models [@kurucz; @allard; @scholz]) using the (modified) algorithm developed by Cuisinier et al. (1997, see Paper I). For each temperature, it calculates a correction function which must be multiplied with the theoretical spectrum corresponding to the parameters given in the calibration file for that temperature, such that the corrected spectra reproduce the colours given in the calibration file. This correction function being smooth, it leaves line strength indices untouched and only changes the continuum. This correction function is then used to correct all model spectra of this temperature, whereby the original differential properties with $\log g$ and ${\rm [Fe/H]}$ are usually conserved, and hence, continue to be well represented by the models. In the end, all SEDs are rescaled to make the bolometric fluxes match the temperatures (for a more precise description, see Paper I). Thus, the new feature of the calibration algorithm is to waive this differential approach in favour of calculating a separate correction function each temperature and for each level of ${\rm [Fe/H]}$, the differential properties with $\log g$ of the models are preserved, however. On top of that, the calibration programs were modified to include a smooth transition between the giant-calibrated and the dwarf-calibrated range.\
In this way, we created a metallicity-calibrated library of synthetic stellar SEDs, designed to reproduce observed colour-colour - and colour-temperature relations, which (for historical reasons) will from now on be called the BaSeL 3.1 ”WLBC 99” library.\
As a first test of the new SED library, it was used to derive spectra and colours from the Yale and the Padova isochrones for the grid metallicities and for age 0 and 14 Gyr and also for the observed cluster metallicities at age 14 Gyr. As can be seen from figs. \[WLBC99coltemprelsgi\] and \[WLBC99coltemprelsdw\] - where the input (calibration) relations (solid) are shown against the output (model) relations (dashed) and the BaSeL 2.2 semi-empirical relations (dotted) - the used colour - temperature relations are well-reproduced. This mainly confirms the functionality of the calibration algorithm, but also suggests that the colour - temperature relations of the models can be trusted. Because these relations have been calibrated in a metallicity-dependent way, the present results show improvements over earlier results obtained from the semi-empirical (BaSeL 2.2) library calibration in $U-B$, $J-H$, $H-K$ and $J-K$, especially at low temperatures. They also appear smoother than their predecessors. Below the lower temperature limits of the calibration relations specified in table \[ranges\], they still show the same discontinuities as the BaSeL 2.2 models, due to lack of data, but for population synthesis, this shouldn’t pose major problems, as stars of these low temperatures either don’t show up at all (giants), or contribute only negligibly to the integrated light of a population (dwarfs).

![Empirical $(B-V)$ - $M_{V}$ CMD of the globular clusters 47 Tuc and M92 (solid). Overlaid are the CMDs created by combining the Padova 10 Gyr isochrone for ${\rm [Fe/H]} = -0.70$ and the 16 Gyr isochrone for ${\rm [Fe/H]} = -2.16$ with the BaSeL 3.1 ”WLBC 99” library (dashed) resp. with the BaSeL 2.2 semi-empirical library (dotted).[]{data-label="WLBC99isochronesBV"}](MS1686f5.eps){width="\columnwidth"}

![Empirical $(V-K)$ - $M_{V}$ CMD of the globular clusters 47 Tuc and M92 (solid). Overlaid are the CMDs created by combining the Padova 10 Gyr isochrone for ${\rm [Fe/H]} = -0.70$ and the 16 Gyr isochrone for ${\rm [Fe/H]} = -2.16$ with the BaSeL 3.1 ”WLBC 99” library (dashed) resp. with the BaSeL 2.2 semi-empirical library (dotted).[]{data-label="WLBC99isochronesVK"}](MS1686f6.eps){width="\columnwidth"}

The second test unfortunately yields much less satisfying results. $M_{V} - 
   UBVRIJHK$ CMDs were produced using Yale and Padova isochrones for the observed metallicities of the calibrating globulars and typical cluster ages (eg. 10, 12, 12, 14 and 16 Gyr for 47 Tuc, M 5, M 3, NGC 6397, and M 92 respectively), in order to reproduce their CMDs shown in figs. \[CMDsa\] and \[CMDsb\]. Contrary to expectations, these isochrones do not reproduce the cluster CMDs better than the BaSeL 2.2 semi-empirical models. The RGBs come out too steep and in certain colours ($U-B$, $B-V$, $V-I$) too red for high metallicities, and too blue in all colours for low values of ${\rm [Fe/H]}$ (see figs. \[WLBC99isochronesBV\] and \[WLBC99isochronesVK\]).\
Obviously, it is impossible at the present epoch to establish a unique calibration of $UBVRIJHKL$ colours in terms of stellar metallicity that reproduces both empirical colour-temperature relations of stars and the CMDs of stellar populations (using theoretical isochrones). Whether these discrepancies are due to the used colour - temperature relations, or to the isochrones remains to be investigated, but recent indications point towards shortcomings in the isochrones, especially in the convection treatment [@salarisb].\
Figs. \[WLBC99coltemprelsgi\], \[WLBC99coltemprelsdw\], \[WLBC99isochronesBV\], and \[WLBC99isochronesVK\] are only shown for Padova isochrones, but the results also hold for Yale isochrones, and the corresponding figures look very similar. The SED library is available by public ftp from the university of Basel as the BaSeL 3.1 ”WLBC 99” version. Because of its good agreement for high metallicities, the (unmodified) ${\rm [Fe/H]}=+0.5$ file from the BaSeL 2.2 library has been added to extend the metallicity range.

Pragmatic solution: application-oriented library calibration, the BaSeL 3.1 ”Padova 2000” library
=================================================================================================

As obviously both the colour-temperature relations of stars and the CMDs of stellar populations (using existing theoretical isochrones) cannot be matched at the same time, we decided to also produce a library that is able to reproduce the CMDs of globulars well, because such a library should at least yield reliable integrated colours for stellar populations and hopefully reliable integrated spectra. Another motivation was the fact that according to some authors [@brocato], colour-temperature relations, in particular for cool stars, aren’t that well-known anyway. These authors even consider them a free parameter.\
As a theoretical isochrone library, which, combined with our to-be-produced SED library, should reproduce the CMDs of our calibration globular clusters, we chose the Padova 2000 isochrones [@padova], because it is widely in use nowadays.

We produced this library in an iterative process that is based on the fact (observed by us), that the $M_{V}$ magnitude, derived for a set of parameters ${\rm [Fe/H]}$, ${T_{\rm eff}}$, and $\log g$ from a stellar library, is practically independent of the choice of the spectral library. This is easily explained by the fact that $M_{bol}$ depends only on the stellar parameters (thus is independent of the choice of library), so the difference in $M_{V}$ stems only from the difference in the bolometric correction $BC(V)$ derived from the different libraries, which is negligible on the scale of absolute magnitudes. This library independence of $M_{V}$ can now be used to assign the colours of one’s choice to a certain $M_{V}$ value, by assigning them to the ${T_{\rm eff}}$ and $\log g$ values that will reproduce this $M_{V}$ value in the to-be-calibrated library, which one knows already from an existing library. This way one can in principle shape a synthetical colour-magnitude diagram in the way one desires[^2].\
The iterative process employed was the following:

1.  New calibration files were created from the ${T_{\rm eff}}$-$\log g$ relations of the Padova 2000 isochrones for 0 Gyr for dwarfs (the age at which the main sequence is still complete) and 14 Gyr for giants (a typical globular cluster age), and colour-temperature relations derived via the colour-$M_{V}$ relations from the globular cluster fiducials presented in Section \[compilation\] and the $M_{V}$-${T_{\rm eff}}$ relations derived with the before-mentioned Padova 2000 isochrones from the previous library (in the first step, the ”previous library” was the ”WLBC 99” library, and in the following steps, it was the library produced in the previous step). For solar metallicity, the calibration file was kept the same as for the ”WLBC 99” library (thus also the same as for the BaSeL 2.2 library), with the small difference, that the ${T_{\rm eff}}$-$\log g$ relation was taken from the Padova 2000 isochrones for consistency.

2.  These calibration files were used to calibrate a new library.

3.  Steps 1 and 2 were repeated until there was no significant difference anymore between subsequent libraries.

After four iteration steps this was the case. \[The calibration files used for the final iteration are also available by ftp from Basel university, although they are probably of little use to the user\].

![image](MS1686f7.eps){width="\textwidth"}

![image](MS1686f8.eps){width="\textwidth"}

As can be seen from figs. \[PAD2000isochronesa\] and \[PAD2000isochronesb\], the above procedure indeed provides the expected improvements. They show the globular cluster CMDs of figs. \[CMDsa\] and \[CMDsb\] compared to the Padova 2000 isochrones of the same (cluster) metallicities and the above-mentioned ages evaluated with the new BaSeL 3.1 ”Padova 2000” library.\
The (metallicity-dependent) shapes and locations of the RGB and the main sequence are well-reproduced for the entire range from ${\rm [Fe/H]}$ = $-2.16$ to $-0.70$, apart from the RGB tip of 47 Tuc, of which the ”bend-down” in $V$ proved virtually impossible to reproduce (for solar ${\rm [Fe/H]}$, the quality of the library has already been confirmed extensively, as in its present form, it is almost identical with the @paperii BaSeL 2.2 library, which is widely used). At a given metallicity, the CMDs of individual clusters, however, can still differ by as much as 0.1$^{m}$ in colour, due to the above-mentioned large intrinsic cluster-by-cluster scatter in the observed CMDs. From the generally good overall agreement, it can be concluded that this new library should be useful for population synthesis, if combined with Padova 2000 isochrones or tracks. It can be retrieved by ftp from Basel university as the BaSeL 3.1 ”Padova 2000” library. Here too, the unmodified ${\rm [Fe/H]}=+0.5$ file from the BaSeL 2.2 library was added.

![image](MS1686f9.eps){width="\textwidth"}

The good agreement in the CMDs comes at the price of rather unusual colour temperature - relations. Fig. \[PAD2000coltemprels\] shows the empirical ${T_{\rm eff}}$ - $V-K$ relation from @ridgway and empirical ${T_{\rm eff}}$ - $\log g$ relations for red giants from Cohen, Frogel and Persson [@cohen; @frogel; @persson ${\rm [Fe/H]}$-dependent] and dwarfs from @angelov [${\rm [Fe/H]}$-independent] for three metallicities from the grid (solid lines) versus the BaSeL 3.1 ”Padova 2000” (dashed) and the BaSeL 2.2 semi-empirical (dotted) relations for the same parameters. For giants, the disagreement reaches up to 500 K or 1$^{m}$ in $V-K$ (for colours with a smaller baseline, this error can be scaled down accordingly); for the lowest metallicities and temperatures (${\rm [Fe/H]} \leq -1.5$ and ${T_{\rm eff}} \leq 4000$ K) it can even reach 800 K, but fortunately, these stars don’t show up in the Padova 2000 isochrones. For dwarfs, these problems aren’t as bad as for giants, but it is clear that temperatures derived from the BaSeL 3.1 ”Padova 2000” models should be treated with scepticism. On the positive side, most of the discontinuities have disappeared in the colour temperature - relations.

![image](MS1686f10.eps){width="\textwidth"}

$$

[lrrrrr]{} $template$ & ${\rm [Fe/H]}$ & U-B & B-V & V-I & R-I\
$G177$ & 0.25 & - & - & 1.235 & 0.653\
$G1-M31cl$ & 0.00 & - & 0.993 & 1.238 & 0.663\
$G2b-M31cl$ & -0.50 & 0.296 & 0.869 & 1.117 & 0.585\
$G3r$ & -1.00 & 0.093 & 0.739 & 1.044 & 0.549\
$G4b$ & -1.50 & 0.041 & 0.676 & 1.032 & 0.551\
$G5$ & -2.00 & -0.001 & 0.639 & 0.948 & 0.507\

$$\end{center}
   \end{table}
%
%__________________________________________________ Padova2000 colours
   \begin{table*}
   \begin{center}
      \caption[]{Magnitudes and colours for single stellar populations
                    sythesised using the BaSeL 3.1 ''Padova 2000'' SED library.}
      \label{padova2000colours}$$

[lrrrrrrrrrrrrr]{} $age$ & ${\rm [Fe/H]}$ & M\_[bol]{} & M\_[V]{} & U-B & B-V & V-I & V-K & R-I & J-H & H-K & J-K & K-L & K-M\
12.000 & 0.25 & 6.582 & 7.491 & 0.776 & 1.045 & 1.269 & 3.346 & 0.644 & 0.698 & 0.229 & 0.927 & 0.178 & 1.094\
14.000 & 0.00 & 6.616 & 7.407 & 0.599 & 0.978 & 1.224 & 3.130 & 0.634 & 0.691 & 0.200 & 0.892 & 0.164 & 1.020\
14.000 & -0.50 & 6.453 & 7.020 & 0.320 & 0.856 & 1.124 & 2.672 & 0.587 & 0.621 & 0.164 & 0.785 & 0.129 & 0.839\
18.000 & -1.00 & 6.523 & 6.976 & 0.155 & 0.771 & 1.051 & 2.366 & 0.532 & 0.555 & 0.135 & 0.691 & 0.120 & 0.735\
20.000 & -1.50 & 6.541 & 6.958 & 0.053 & 0.717 & 1.011 & 2.201 & 0.507 & 0.519 & 0.128 & 0.647 & 0.120 & 0.724\
20.000 & -2.00 & 6.563 & 6.978 & -0.009 & 0.684 & 0.980 & 2.135 & 0.498 & 0.499 & 0.121 & 0.620 & 0.116 & 0.731\

$$\end{center}
   \end{table*}
%
%__________________________________________________ Padova2000 residuals
   \begin{table}
   \begin{center}
      \caption[]{Differences in colours between
                    single stellar populations sythesised using 
                    the BaSeL 3.1 ''Padova 2000'' SED library and
                    template globular cluster integrated energy distributions
                    from Bica et al., in the sense BaSeL 3.1 $-$ template.}
      \label{residuals}$$

[lrrrrr]{} $template$ & ${\rm [Fe/H]}$ & U-B & B-V & V-I & R-I\
$G177$ & 0.25 & - & - & 0.034 & -0.009\
$G1-M31cl$ & 0.00 & - & -0.015 & -0.014 & -0.029\
$G2b-M31cl$ & -0.50 & 0.024 & -0.013 & 0.007 & 0.002\
$G3r$ & -1.00 & 0.062 & 0.032 & 0.007 & -0.017\
$G4b$ & -1.50 & 0.012 & 0.041 & -0.021 & -0.044\
$G5$ & -2.00 & -0.008 & 0.045 & 0.032 & -0.009\

$$

In a last test, the effect of the metallicity-calibration on the integrated spectra of synthetic single-burst stellar populations was investigated. Synthetic populations were created in order to compare their integrated spectra with empirical template globular cluster integrated energy distributions from Bica et al. (1996, private communication). The synthetic populations were created using the GISSEL (Galaxy Isochrone Synthesis Spectral Evolution Library) software by G. Bruzual A. and S. Charlot [2000, for a short description, see @charlot; @bruzual] for the template metallicities (${\rm [Fe/H]} = 0.25$, $0.00$, $-0.50$, $-1.00$, $-1.50$, $-2.00$), and the same ages as Lejeune used in his thesis (1997) for comparison of the BaSeL 2.2 semi-empirical library with the same template spectra (12, 14, 14, 18, 20 Gyr). Both the BaSeL 2.2 semi-empirical - and the BaSeL 3.1 ”Padova 2000” SED libraries were used, along with a Salpeter IMF (mass range: $0.1 - 50 M_{\odot}$) and the Padova 2000 tracks (except for ${\rm [Fe/H]} = -2.0$, where the Padova 1995 isochrone [@fagotto; @girardi] was used, because the corresponding Padova 2000 isochrone wasn’t implemented in the evolutionary synthesis program). The resulting integrated spectra are compared in fig. \[spectra\], where the solid lines stand for the empirical template spectra, the dotted lines for the BaSeL 2.2 semi-empirical library and the dashed lines for the metallicity-calibrated library. Residuals in the sense BaSeL 2.2 $-$ empirical spectra (dotted), resp. BaSeL 3.1 $-$ empirical spectra (solid) are shown in the same figure. To guide the eye, the zero-line is also shown in solid. Synthetic ($UBVRIJHKLM$) colours and magnitudes from the empirical - and the synthetical (BaSeL 2.2 ”Padova 2000”) spectra are given in tables \[bicacolours\], and \[padova2000colours\], and the differences between these colours in table \[residuals\] (due to the limited wavelength range of the Bica et al. spectra, only $UBVRI$ colours can be given for the empirical spectra). As expected, the differences between the two libraries are negligible for high ${\rm [Fe/H]}$. They are hardly visible in fig. \[spectra\], and in the colours, they are of the order of $0.001^{m}$. For low metallicities (${\rm [Fe/H]}=-1.65$), the metallicity-calibration seems to add flux in the red ($R$, $I$) and to reduce the flux in the visible a bit, but clearly, the dramatic differences in stellar colour-temperature relations between the two libraries have only a minor effect on the spectral properties of entire (synthetic) populations. The main improvements can be seen in the visible (the spectral shape in this region is matched perfectly for all metallicities), whereas in the ultraviolet one can still see some strong deviations. However, the differences between empirical integrated colours and the ones of the synthetic population only amount to a few hundredths of a magnitude (see table \[residuals\]), confirming that the BaSeL 3.1 ”Padova 2000” SED library should be a useful tool for evolutionary synthesis, if used in combination with the Padova 2000 tracks/isochrones.\

Conclusions and future work
===========================

1.  We have constructed two comprehensive libraries of theoretical stellar energy distributions combining the Kurucz -, the Allard and Hauschildt -, and the Scholz-models [@kurucz; @allard; @scholz]. They provide synthetic stellar spectra with useful resolution on a homogeneous wavelength grid, from 9.1 nm to 160 $\mu$m, over large ranges of fundamental parameters. In addition to the preceding version, the BaSeL 2.2 semi-empirical library (see Papers I and II), the present libraries have been ($UBVRIJHKL$) colour-calibrated independently at all levels of metallicity, using Galactic globular cluster photometric data, to overcome the weaknesses of the BaSeL 2.2 library at low ${\rm [Fe/H]}$.

2.  It proved impossible to calibrate the library in such a way that both the empirical colour-temperature relations (taken from @ridgway) and the CMDs of the calibrating clusters (using Yale or Padova 2000 isochrones) could be matched at the same time. Whether this is due to the isochrones or to irrealistic colour-temperature relations still remains to be investigated. Of great help would also be a more complete (colour) data basis. For this reason, we created two metallicity-calibrated libraries, the BaSeL 3.1 ”WLBC 99” library, which is able to reproduce empirical colour-temperature relations and should therefore be a useful tool for investigations on a stellar level, and the BaSeL 3.1 ”Padova 2000” library which, in combination with the Padova 2000 isochrones, successfully reproduces the CMDs (i. e. the location and shape of the different branches) of Galactic globular clusters, and which should therefore be a particularly suitable tool for spectral evolutionary synthesis studies of stellar systems and other synthetic photometry applications, especially if used in combination with the Padova 2000 isochrones/tracks.

3.  As the calibration algorithm is very flexible, new libraries can easily be produced as soon as better data and/or SED models are available. The present results, therefore, do not necessarily provide the final version of the BaSeL models.

4.  As a first application, we are using the ”Padova 2000 tracks” library for the analysis of the colours of chemo-evolutionary dynamical models of spiral galaxies [@diss; @westerasamland].

This work was supported by the Swiss National Science Foundation.

Alcaino, G. 1977, A&AS, 29, 397

Alcaino, G., Liller, W. 1985, A&A, 146, 389

Alcaino, G., Liller, W. 1986, AJ, 91, 87

Alcaino, G., Liller, W., Alvarado, et al. 1997, AJ, 114, 1067

Allard, F., & Hauschildt, P. H. 1995, ApJ, 445, 433

Alonso, A., Arribas, S., & Martínez-Roger, C. 1999, A&AS, 140, 261

Angelov, T. 1996, BOBeo, 153, 19

Arribas, S., & Martínez-Roger, C. 1987, A&A, 178, 106

Barmby, P., & Huchra, J. P. 2000, ApJ, 531L, 29

van den Bergh, S. 1988, AJ, 95, 106

Bessell, M. S. 1979, PASP, 91, 589

Bessell, M. S. 1983, PASP, 95, 480

Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134

Bessell, M. S., & Weiss, E. W. 1987, PASP, 99, 642

Bica, E., Ortolani, S., & Barbuy, B. 1994, A&A, 106, 161

Bica, E., Alloin, D., Bonatto, C., et al. 1996 (private communication to T. Lejeune)

von Braun, K., Chiboucas, K., Minske, J. C., Salgado, J. F., & Worthey, G. 1998, PASP, 110, 810

Brocato, E., Castellani, V., Poli, F. M., & Raimondo, G. 2000, A&AS, 146, 91

Bruzual A., G. 1999, in: Galaxy Interactions at Low and High Redshift, J. E. Barnes & D. B. Sanders (eds), IAU Symp., 186, Dordrecht: Kluwer, p.459

Bruzual A., G., & Charlot, S. 1993, ApJ, 405, 538

Bruzual A., G., & Charlot, S. 2000, Galaxy isochrone spectral synthesis evolution library (private communication)

Bruzual A., G., Barbuy, B., Ortolani, S. et al. 1997, AJ, 114, 1531

Buser, R., & Kurucz, R. L. 1992, A&A, 264, 557

Cannon, R. D. 1974, MNRAS, 167, 551

Carretta, E., & Gratton, R. G. 1997, A&AS, 121, 95

Cathey L. R. 1974, AJ, 79 1370

Charlot, S., & Bruzual A., G. 1991, ApJ, 367, 126

Cohen, J. G., Frogel, J. A., & Persson, S. E. 1978, ApJ, 222, 165

daCosta, G. S., & Armandroff, T. E. 1990, AJ, 100, 162

Demarque, P., & McClure, R. D. 1977, ApJ, 213, 716

Drissen, L., & Shara, M. M. 1998, AJ, 115, 725

Eggen, O. J. 1972, ApJ, 172, 639

Evans, T. L. 1983, SAAOC, 7, 86

Fagotto F., Bressan A., Bertelli G., & Chiosi C., 1994 A&AS 105, 39

Ferraro, F. R., Carretta, E., Corsi, C. E., et al. 1997, A&A, 320, 757

Fricke, K. J., Izotov, Y. I., Papaderos, P., Guseva, N. G., & Thuan, T. X. 2001, AJ, 121, 169

Frogel, J. A., Persson, S. E., Matthews, K., Aaronson, M. 1978, ApJ, 220, 75

Frogel, J. A., Persson, S. E., & Cohen, J. G. 1981, ApJ, 246, 842

Frogel, J. A. , Persson, S. E., & Cohen, J. G. 1983, ApJ, 53, 713

Girardi L., Bressan A., Chiosi C., Bertelli G., & Nasi E., 1996, A&AS 117, 113

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 (”Padova 2000” isochrones)

González Delgado, R. M., García-Vargas, M. L., Goldader, J. et al. 1999, ApJ, 513, 707

Grebel, E. K., & Roberts, W. J. 1995, A&AS, 109, 293

Harris, W. E. 1996, AJ, 112, 1487

Hesser, J. E., & Hartwick, F. D. A. 1977, ApJS, 33, 361

Hesser, J. E., Harris, W. E., VandenBerg, D. A. et al. 1987, PASP, 99, 739

Johnson, H. L., & Sandage, A. R. 1956, ApJ, 124, 379

Johnson, J. A., & Bolte, M. 1998, AJ, 115, 693

Johnson, K. E., Leitherer, C., Vacca, W. D., & Conti, P. S. 2000, AJ, 120, 1273

Kaluzny, J. 1997, A&AS, 122, 1

Kaluzny, J., Wysocka, A., Stanek, K. Z., & Krzeminsky, W. 1998, Acta Astron., 48, 439

Kauffmann, G., & Charlot, S. 1998, MNRAS, 294, 705

King, I. R., Anderson, J., Cool, A. M., & Piotto, G. 1998, ApJ, 492L, 37

Kong, X., Zhou, X., Chen, J. et al. 2000, AJ, 119, 2745

Kotilainen, J. K., Reunanen, J., Laine, S., & Ryder, S. D. 2001, A&A, 366, 439

Kurucz, R. L. 1995 (private communication to R. Buser)

Lastennet, E., Lejeune, T., Westera, P., & Buser, R. 1999, A&A, 341, 857

Lastennet, E., Valls-Gabaud, D., & Oblak, E. 2000, in: Birth and Evolution of Binary Stars, B. Reipurth & H. Zinnecker (eds), IAU Symp., 200, Dordrecht: Kluwer, p. 164

Lastennet, E., Lignières, F., Buser, R., et al. 2001, A&A, 365, 535

Lee, S. W. 1977, A&AS, 27, 381

Lee, S-G., Lee, M. G., & Kim, E. 1996, Journal of the Korean Astron. Soc., 29, 171

Lee, H., Yoon, S, & Lee, Y. 2000, AJ, 120, 998

Leitherer, C., Schaerer, D., Goldader, J. D. et al. 1999, ApJS, 123, 3

Lejeune, T., Cuisinier, F., & Buser, R. 1997, A&AS, 125, 229 (Paper I)

Lejeune, T., Cuisinier, F., & Buser, R. 1998, A&AS, 130, 65 (Paper II)

Lejeune, T. 1997, Contribution à la synthèse spectrale évolutive d’amas et de galaxies par calibration de bibliothèques stellaires, PhD thesis, Univ. Strasbourg & Univ. Basel, 334 pp.

Liu, M. C., Charlot, S., & Graham, J. R. 2000, ApJ, 543, 644

Lotz, J. M., Ferguson, H. C., & Bohlin, R. C. 2000, ApJ, 532, 830

Maraston, C., & Thomas, D. 2000, ApJ, 541, 126

Marconi, G., Buonanno, R., Carretta, E. et al. 1998, MNRAS, 293, 479

Marleau, F. R., Graham, J. R., Liu, M. C., & Charlot, S. 2000, AJ, 120, 1779

Menzies, J. 1973, MNRAS, 163, 323

Mermilliod, J.-C. 1998, The General Catalogue of Photometric Data, http://obswww.unige.ch/gcpd/

Mollá, M., & García-Vargas, M. L. 2000, A&A, 359, 18

Montegriffo P., Ferraro F. R., Fusi Pecci F., & Origlia L. 1995, MNRAS, 276, 739

Newell, E. B., Rodgers, A. W., & Searle, L. 1969, ApJ, 156, 597

Nikolaev, S., & Weinberg, M. D. 2000, ApJ, 542, 804

Norris, J., & Freeman, K. C. 1982, ApJ, 254, 143

Origlia, L., Goldader, J. D., Leitherer, C., Schaerer, D., & Oliva, E. 1999, ApJ, 514, 96

Piotto, G., Cool, A. M., & King, I. R. 1997, AJ, 113, 1345

Richer, H. B., & Fahlman, G. G. 1987, ApJ, 316, 189

Ridgway, S. T., Joyce, R. R., White, N. M., & Wing, R. F. 1980, ApJ, 235, 126

Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618

Salaris, M. 2001, Theoretical uncertainties in Red Giant Branch stellar models, in: Observed HR diagrams and stellar evolution: the interplay between observational constraints and theory, T. Lejeune, & J. Fernandes (eds), ASP Conference Series (submitted)

Salaris, M., Chieffi, A., & Straniero O. 1993, ApJ, 414, 580

Sandage, A. 1969, ApJ, 157, 515

Sandage, A. 1970, ApJ, 162, 841

Sandage, A., & Walker, M. F. 1966, ApJ, 143, 313

Sandquist, E. L., Bolte, M., Stetson, P. B., & Hesser, J. E. 1996, ApJ, 470, 910

Saviane, I., Rosenberg, A., Piotto, G., & Aparicio, A.. 2000, A&A, 355, 966

Scholz, M. 1997 (private communication to R. Buser)

Stetson, P. B., & Harris, W. E. 1988, AJ, 96, 909

Westera, P. 2001, The BaSeL 3.1 models: Metallicity calibration of a theoretical stellar spectral library and its application to chemo-dynamical galaxy models, PhD thesis, Univ. of Basel, 378 pp.

Westera, P., Samland, M., Buser, R., & Gerhard, O. E., 2001, A&A, Colour Evolution of Disk Galaxy Models. (submitted)

Woolley, R., Alexander, J. B., Mather, L., Epps, E. 1961, Royal Observatory Bulletins, 43, 303

[^1]: In the infrared colours, there were barely enough data to derive a trend, but certainly not enough to determine the amplitude of a possible scatter around it.

[^2]: Note that this procedure implies giving up consistency with empirically established ${T_{\rm eff}}$-colour relations.
