A&A 442, 1127-1134 (2005)
DOI: 10.1051/0004-6361:20042490
U. Munari1 - R. Sordo1,2 - F. Castelli3 - T. Zwitter4
1 - INAF Osservatorio Astronomico di Padova, via dell'Osservatorio 8, 36012 Asiago (VI), Italy
2 -
Dipartimento di Astronomia, Univ. di Padova, vicolo della Specola 5, 35122 Padova, Italy
3 -
INAF - Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, 34131 Trieste, Italy
4 -
University of Ljubljana, Department of Physics, Jadranska 19, 1000 Ljubljana, Slovenia
Received 6 December 2004 / Accepted 19 July 2005
Abstract
We present a complete library of synthetic spectra based on
Kurucz's codes that covers the 2500-10 500 Å wavelength range at
resolving powers
,
11 500 (
GAIA), 8500
(
RAVE), 2000 (
SLOAN) and uniform dispersions of 1 and
10 Å/pix. The library maps the whole HR diagram, exploring 51 288
combinations of atmospheric parameters spanning the ranges: 3500
K,
,
,
[
/Fe] = 0.0,+0.4,
= 1, 2, 4 km s-1, 0
500 km s-1. The spectra are available both
as absolute fluxes as well as continuum normalized. Performance tests and
spectroscopic applications of the library are discussed, including automatic
classification of data from spectroscopic surveys like RAVE, SLOAN,
GAIA. The entire library of synthetic spectra is accessible via the web.
Key words: astronomical data bases: miscellaneous - stars: atmospheres - stars: fundamental parameters - surveys
Large and complete synthetic spectral libraries of normal
stars are of interest for many different applications, including (a)
automatic analysis and classification of large volumes of data, like those
collected by ongoing spectral surveys, (b) derivation of radial velocities
via cross-correlation against best matching templates, (c) calibration of
spectroscopic line/band classification criteria, and (d) calibration of
photometric indexes. Tasks a-c require the synthetic spectra to match
the absorption lines in observed spectra normalized to the local continuum,
generally irrespective of how well the overall spectral energy distribution
is reproduced. The goal of task d is instead to match the energy contained
within photometric bands, irrespective of the performance on the individual
absorption lines. Therefore, synthetic spectral libraries well suited for
one task do not necessarily perform equally well on another. For example,
within the Kurucz's synthetic spectra suite, the inclusion in the
computations of the so-called "predicted lines'' (hereafter PLs, Kurucz 1994)
limits the range of applications. The PLs are calculated absorption lines
due to transitions involving one and sometimes two levels whose locations
are predicted by atomic structure codes. Wavelengths for these lines may be
uncertain by 10 Å, with large uncertainties also affecting their
computed intensity. As a consequence the predicted lines may not correspond
in position and size to the observable counterparts (cf. Bell et al. 1994;
Castelli & Kurucz 2004). For all these reasons the predicted lines are not
intended to be used for the analysis of high resolution spectra, but are an
essential contribution to the total line blanketing in model atmospheres and
spectral energy distribution computations. Very accurate line wavelengths
and shapes are not needed for the statistical computations of the
line opacity in model atmospheres or for the prediction of low
resolution spectra.
Table 1:
Comparison of the principal characteristics of publically
available libraries of synthetic spectra with resolution
better than 5 Å. N: number of combinations of
atmospheric parameters and rotational velocities. new
ODF: use of the new opacity distribution functions.
TiO: inclusion of the TiO line-list. pred. lines:
inclusion of the "predicted lines''. The last three columns
indicate if the spectra are provided for different values of
enhancement of -elements, micro-turbulence velocity
and rotational velocity.
Various libraries of synthetic spectra are becoming available, and they are summarized and compared in Table 1. These libraries widely differ in wavelength interval, spectral resolution, input model atmospheres and range of atmospheric parameters explored, so that their overlap is not significant.
As for other libraries in Table 1, our one is based on Kurucz's atmospheres,
line-lists and computing software. The main goal of our spectral library is
spectroscopic applications of the a and b types above, which imposes
the exclusion of the PLs from computations and therefore does not claim
accurate broad-band photometric performances particularly at lower
temperatures. Other relevant characteristics of our library are the various
resolutions and rotational velocities at which it is provided, its wide
wavelength range, the adoption of improved model atmospheres based on the
new opacity distribution functions (ODFs) by Castelli & Kurucz (2003), the
use of the TiO line list of Schwenke (1998), the inclusion of
-element enhancement and different micro-turbulent velocities.
Table 2: Range of atmospheric parameters explored by our synthetic spectral library.
The whole grid of spectra in our library was computed using the SYNTHE code
by Kurucz (Kurucz & Avrett 1981; Kurucz 1993), running under the VMS operating
system on a Digital Alpha workstation in Asiago. We adopted as input model
atmospheres the ODFNEW models
(http://wwwuser.oat.ts.astro.it/castelli/grids/; Castelli & Kurucz 2003).
They differ from the NOVER models (http://kurucz.harvard.edu/grids.html;
Castelli et al. 1997) in the adoption of new ODFs, replacement of the solar
abundances from Anders & Grevesse (1989) with those from Grevesse & Sauval
(1998), and improvements in the molecular opacities including the adoption
of the molecular line-lists of TiO by Schwenke (1998, as distributed by
Kurucz 1999a) and of H2O by Partridge & Schwenke (1997, as distributed
by Kurucz 1999b). For the combination of atmospheric parameters for which no
ODFNEW models were available at the time of writing, we adopted the
corresponding NOVER input model atmospheres. We plan to update the library
from NOVER to ODFNEW models as more of the latter will become available and
to post the result on the library web page (see below). Both
the NOVER and the ODFNEW models were computed with the overshooting option
for the mixing-length convection switched off, while Kurucz (K) atmospheric
models (also available from http://kurucz.harvard.edu/grids.html) were computed with the overshooting option switched on. Several papers have
demonstrated that for stars with active convection (
K)
the no-overshooting convection treatment provides better agreement with the
observations than the overshooting case does (Castelli et al. 1997; Smalley
& Kupka 1997; Gardiner et al. 1999; Smalley et al. 2002). The
no-overshooting models used by us were computed for the mixing-length
parameter to the scale height of 1.25. This value allows to fit the observed
solar irradiance, while a lower value, like that of 0.5 suggested by Smalley
et al. (2002), does not. On the contrary, 0.5 seems to better fit the
wings of H
,
provided that the position of the continuum is known with
an uncertainty smaller than 1%. In fact, a difference of 1% in the
position of the solar continuum corresponds to the difference between 0.5
and 1.25 of the mixing-length parameter (Castelli et al. 1997). Generally,
it is very difficult to determine the location of the continuum across the wings
of the Balmer lines, especially in Echelle spectra which are notoriously
severely affected by the blaze function.
In computing the synthetic spectra we used the Kurucz atomic and molecular
line-lists (Kurucz 1992). The effects of several molecules were taken into
account, including C2, CN, CO, CH, NH, SiH, SiO, MgH and OH. The TiO
molecular lines (from Schwenke 1998) were included in the computation when
K. The synthetic spectra are calculated and made
available both as fluxes for unit wavelength and as normalized intensities.
The latter have a different meaning from the normalization to the continuum
of the observational spectroscopy, where the location of the
continuum is taken to be represented by the regions away from absorption
lines or bands. In synthetic spectra, a continuum-normalized spectrum is
obtained by dividing the absolute fluxed spectrum by its calculated
continuum. Noticeable differences are limited to the head of the
Paschen and Balmer series of hydrogen, and to the coolest stars dominated by
molecular bands.
The spectra in our library covers from 2500 to 10 500 Å and were
calculated at a resolving power
and then degraded by Gaussian convolution to lower
resolving powers and properly re-sampled to the Nyquist criterion (the FWHM of
the PSF being 2 pixels), to limit the data volume and therefore facilitate
the distribution.
The highest resolving power for our library is
,
which is
typical of moderately high resolution Echelle observations, such as the spectral
database of real stars of Munari & Tomasella (1999) and Marrese et al.
(2003). The library is also provided at the resolving powers of the SLOAN
survey
(York et al. 2000), of the RAVE survey
(Steinmetz 2003, Munari et al. 2005) and of the Gaia space
mission
(Katz et al. 2004, Wilkinson et al. 2005).
The present library is also made available at the uniform
dispersions of 1 and 10 Å/pix (again with the FWHM of the PSF
being 2 pixels). The latter two are intended to be of assistance in
classification of medium and low resolution grating spectra obtained in
single dispersion mode (as with conventional B&C spectrographs or
EFOSC-like imager/spectrographs).
Each spectrum is provided for a range of rotational velocities, as detailed in Table 1, spanning 14 values between 0 and 500 km s-1 for stars hotter than 6000 K and 11 values between 0 and 100 km s-1 for cooler ones.
The explored combinations of parameters in Table 2 give rise to 51 228
spectra, each provided in six different combinations of resolving power and
sampling, and two flux varieties (absolute units and normalized
intensities). The total number of individual spectra contained in this
library is therefore
.
The spectra are
provided as gzipped ASCII files containing only the flux column, the common
wavelength column being given separately only once for each of the six
resolutions. The adopted scheme for the naming of files is outlined in
Table 3.
All spectra in the library can be directly accessed and retrieved through the dedicated web page http://archives.pd.astro.it/2500-10500/. The version of the library at 1 Å/pix is accessible also via ESA's web site http://gaia.esa.int/spectralib/, where browsing facilities based on Virtual Observatory tools are provided. A distribution via DVDs will be possible in special cases (to be arranged directly at munari@pd.astro.it).
Table 3: Naming scheme for the spectra contained in the library.
To evaluate the spectroscopic and photometric performances of our library we have carried out several tests. They support the use of the library in spectroscopic applications (our main goal), while caution should be applied in the photometric use at lower temperatures and/or shorter wavelengths. It is worth noticing that significant discrepancies affect the comparison of available libraries with observational data over blue wavelengths and low temperatures.
A heavily demanding spectroscopic application of our library has been
carried out by Siviero et al. (2004) and Marrese et al. (2005) in
investigation of double-lined solar type eclipsing
binaries, as a source for (i) best matching templates in cross-correlation
determination of radial velocities and (ii) a reference grid in analysis of atmospheric properties of binary components. They used the
version of the library.
![]() |
Figure 2:
The effect of ![]() ![]() |
Open with DEXTER |
The availability of a wide and complete range of rotational velocities for all the spectra in our library always allowed them to select the best template even with binaries characterized by high rotation, a key ingredient to reach absolute precisions better than 1 km s-1 in individual radial velocities even on modest S/N spectra.
The derivation of atmospheric parameters via
fitting to the library
converged independently to the same values as derived geometrically by the
orbital solution. The
and orbital determinations provide the same
gravities within 0.03 dex, the same rotational velocities within
4 km s-1 and the same difference in temperature between the two
components within 42 K (0.6%). The metallicity was indirectly derived by
the orbital solution via matching with theoretical isochrone tracks on the
temperature-luminosity plane. The values derived in this way are confirmed
within 0.06 dex by direct
analysis. Given the significant rotation
exhibited by components of binaries (generally synchronous with orbital
motion) the wide range of rotational velocities covered by our library
proved once again to be essential to obtain accurate
results.
The above results would have been significantly degraded if we had used
a library computed with inclusion of the PLs.
To support this assertion we have re-run exactly the same fits of
Siviero et al. (2004) and Marrese et al. (2005) using an identical grid of
synthetic spectra, this time computed with the inclusion
of PLs. In the case of V570 Per the use of spectra with PLs pushes the
fit away from the geometrical solution by
K and
[M/H] = 0.16, and by
K,
and
[M/H] = 0.11 for V432 Aur
(average values for the components of the given binary). The differences
are always in the sense that inclusion of PLs pushes the
fit away
from the correct one by going toward cooler temperatures and lower
metallicities. This is an expected, systematic trend. In fact, the inclusion
of PLs causes the continuum in synthetic spectra to become fainter (cf.
bottom panel of Fig. 2), causing a reduction in the contrast between strong
absorption lines (those dominating the outcome of the
fit) and the
adjacent continuum. To re-establish the proper contrast, the
fits a
cooler temperature. Also fitting a lower metallicity helps to re-establish
the proper contrast. PLs are mainly a forest of weak, optically
thin lines that respond much faster to a small change in metallicity than
the stronger absorption lines which are much closer to optically thick
conditions than PLs.
We observed also a reduction in the accuracy of the radial velocities, which
is however more difficult to quantify. It appears as a systematic
shift (toward shorter or longer wavelengths) dependent on the actual
wavelength interval and the orbital phase. This is caused by the beating in
the cross-correlation between the pattern of PLs and the observed absorption
lines which is a function of both
and radial velocity separation
between the binary components. The presence of these locally different
systematic shifts becomes evident only when radial velocities from various
wavelength intervals on the same spectrum are inter-compared (like those
provided by individual orders in an Echelle spectrum). The overall effect is
an increase of the scatter of individual measurements, different in amount from
spectrum to spectrum.
This concerns the analysis performed on wavelength regions dominated by weak
PLs, as in the example on the bottom panel of Fig. 2. There are however
limited wavelength intervals where PLs much stronger than usual are
stochastically clustered together. The result of a
fit performed on
any such interval would provide unrealistic results. One such interval is
illustrated in Fig. 3, where the observed spectra of V432 Aur and V570 Per
are taken from Siviero et al. (2004) and Marrese et al. (2005). While the
fit with synthetic spectra computed without PLs provide - for the values of the
atmospheric parameters derived by the geometrical orbital solution - an
almost perfect match to the observed spectra, the same spectra computed with
inclusion of PLs bear no resemblance at all to the observed spectrum. It is
evident that the photometric calibration of a color-index that includes a
narrow band centered on such a wavelength interval would provide unrealistic
results if carried out on synthetic libraries computed with PLs.
![]() |
Figure 3:
The observed spectra (solid thick lines) of two eclipsing
binaries compared with synthetic ones computed without (solid thin
line) and with (dotted line) inclusion of so-called "predicted
lines''. The synthetic spectra are computed for the atmospheric
parameters (
![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
The use of our synthetic spectral
library at resolving power
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We are engaged in the data analysis of the RAVE spectroscopic survey of the
galactic Halo (Steinmetz 2003; Munari et al. 2005) with the 6dF fiber
positioner and spectrograph at the UK-Schmidt at the Anglo-Australian
Observatory. The survey is carried out at high S/N over the 8420-8780 Å interval at
resolving power.
Two low galactic latitude calibration fields rich in HD stars (1607m492 and
1716m425) have been re-observed three times during the survey. Spectral
classification of these stars is taken from the Michigan Project (Houk
1978). To test the performance of our library, we have
fitted the
spectra of these HD stars against our library of synthetic spectra in its
version. The results are
presented in Fig. 4. The right panels illustrate the accuracy of the
achieved fits by over-plotting observed and synthetic spectra for the sample
cases of a Sun analogue and of a A-type star. The
and
values derived by the
fit well match those expected from the G2V
and A7IV classifications following the calibration by Straizys & Kuriliene
(1981). For the whole sample of HD stars observed in the two calibration
fields, the left panels of Fig. 4 compare the
and
values derived by the
fitting with those expected on the basis of
the spectral classification and the Straizys & Kuriliene (1981)
calibration. The results in Fig. 4 show an agreement within 3% for
and better than 0.25 dex in
,
the dispersion being
almost entirely accountable for by the natural width in both
and
of a spectral sub-type and a luminosity class in the MK
classification scheme.
The present library, in the 1 Å/pix version, was used in differential broad-band photometric applications by Munari et al. (2004) while deriving a precise distance to the Pleiades via orbital solution of the eclipsing binary member star HD 23642, whose spectral type is A2V.
Magnitudes in ten different photometric systems were found in the literature for
this star, and they were used to derive temperature, gravity and reddening
of both components of the binary via calibration on our spectral libraries
of the color indexes of each photometric system. The zero point of colors
was fixed to 0.00 for a Vega model characterized by
,
,
[M/H] = -0.5 and
= 0.0 km s-1 downloaded from
Kurucz's web site.
The photometric results on atmospheric and reddening parameters so derived for both components of HD 23642 were confirmed by high resolution spectroscopy and the orbital solution. The successful photometric application of the library was made possible by the proximity in the parameter space of HD 23642 and the calibrator Vega.
To test the broad-band photometric performances of our library we have
investigated the reproduction of observed colors along the Main Sequence
(MS). The results are shown in Figs. 5 and 6, where the
and
of stars along the main sequence were adopted from Straizys &
Kuriliene (1981). We stopped at M2V as the MS faint end, because later types
are characterized by
cooler than the 3500 K limit of our
library (cf. Table 4). We included for comparison all the other
libraries from Table 1 for which data are actually accessible and are
characterized by a wavelength coverage wide enough to cover at least three
consecutive photometric bands. Caution is
necessary when considering the results from Martins et al. (2005) library,
whose red limit is placed at 7000 Å where the V band still has 1%
transmission and reaches 0.0 only at 7400 Å. The effect of missing this
far red wing of the V transmission profile has effects on colors of cooler
stars by making their B-V color marginally bluer.
![]() |
Figure 5: Comparison of synthetic colors along the Main Sequence truncated at M2V. The observed empirical Main Sequence of Fitzgerald (1970) serves as a comparison. The adopted band transmission profiles are those of Buser (1978). |
Open with DEXTER |
Table 4: Comparison between different calibrations for colors and temperatures of late type main sequence stars. a: Bessell (1990); b: Bessell (1995); c: Straizys (1992); d: Bertone et al. (2004) using Kurucz spectra; e: Bertone et al. (2004) using Phoenix spectra; f: Drilling & Landolt (2000); g: Straizys & Kuriliene (1981).
To assist in the interpretation of the result of the photometric tests we have also recomputed the spectra of our library that define the MS, this time with inclusion of PLs (these spectra are available on the library web page). In Figs. 5 and 6 by "NOVER 1221 pt'' and "ODFNEW 1221 pt'' we mean the flux distributions that are predicted directly from the model atmosphere code, not to be confused with the spectra generated by the spectral synthesis code.
![]() |
Figure 6:
The behaviour of photometric colors along
the Main Sequence for individual spectral types. The same data of
Fig. 5 are used. The adopted band transmission profiles for ![]() ![]() |
Open with DEXTER |
Figure 5 illustrates the behaviour on classical 2-color diagrams. Our library performs similarly to the others on the U-B, B-V plane for colors bluer than 0.7 in B-V. For redder colors, the synthetic MS has U-B colors too blue by some hundredths of a magnitude. The inclusion of PLs in the computation brings the synthetic MS in close agreement with the observed one and very close to the ODFNEW 1221 pt flux distributions. The effect of adopting the ODFNEW model atmospheres affects the lowest temperatures where the synthetic spectra (with and without PLs, as well as the 1221 pt flux distributions) fail to match the very red colors observed in real stars, with the synthetic MS displaying a small upward hook.
Figure 6 illustrates the data in a different format and expanded scale. The
difference in the single color between real stars and synthetic spectra are
plotted as function of the spectral type along the MS. This diagram
highlights better the scatter between observed colors and those from
synthetic spectra, particularly in U-B. It shows that, in all colors,
significant departures of our library from observations are limited to
spectral types cooler than K4V (
K).
Having already proved the satisfactory spectroscopic performance of our
library, we refrain from discussing any further its broad-band photometric
performances and that of the other libraries listed in Table 1, which is
clearly beyond the scope of this paper. Figures 5 and 6 should not be
considered as the ultimate comparison of the photometric performances of the
available synthetic libraries, even if we placed all libraries
on the same reference system, using the same UBVband profiles and the same
/
conversion to spectral
type for all of them. Even if different reconstruction of the
UBV
band profiles or
/
conversion would have been used, the spread among different synthetic
libraries and with respect to observed colors would have remained.
Table 4 shows the widely different
calibration that have been
published for the cooler spectral types along the MS considered in Figs. 5
and 6.
Acknowledgements
Support to R.S. by MIUR COFIN2001 program is acknowledged. We would like to thank S.Ansari for managing the library access via the ESA web portal within the framework of GAIA support activities, A. Baruffolo and S. Pastore for library access via the web server of the Astron. Obs. of Padova, and M. Fiorucci for useful support with treatment of RAVE data. The careful scrutiny and stimulating comments by the anonymous referee were deeply appreciated.