Synphot Data User's Guide


A.3 The Castelli and Kurucz 2004 ATLAS9 of Model Atmospheres

The atlas contains 476 stellar atmosphere models for each atmosphere grid. The models were computed for 76 effective temperature values, 11 gravity values (although for a given temperature not all fluxes for all gravity values were calculated), and for 8 different metallicities. Table 1 in Castelli and Kurucz 2003 summarizes the available models for all grids. These LTE models with no convective overshooting computed by Fiorella Castelli, use improved opacity distribution functions (ODFs) and updated abundances upon previously used by Kurucz (1990).

The main improvements as summarized in Castelli & Kurucz 2003, are:

  1. solar abundances from Grevesse & Sauval 1998, in place of the older ones from Andres & Grevesse (1989, GCA,53,197). See Table A.3 of Castelli & Kurucz 2003;
  2. TiO lines from Schwenke 1998 in place of those from Kurucz (1993);
  3. The addition of the H_2 O lines (Partridge & Schwenke 1997) as distributed by Kurucz (1999a, 1999b), and the HI-HI and HI-H+ quasi-molecular absorptions near 1600 A and 1400 A computed according to Allard et al. 1998; plus additional changes to atomic and molecular data.

These models are computed with the same wavelength resolution and a smaller temperature resolution than the Kurucz 1993 models. Most models have the same number of plane parallel layers from log(tau_Ross)=-6.875 to +2.00 in steps of Delta[log(tau_Ross)] = 0.125, and all are computed assuming a mixing-length convection (no overshooting) with 1/Hp=1.25. As for earlier models, a microturbulent velocity of 2 km s-1 is used.

These models are for metallicities [M/H]=0.0, -0.5, -1.0, -1.5, -2.0, -2.5, +0.5, +0.2 and gravity values from log_g= 0.0 to +5.0 in steps of +0.5. The range in effective temperature from 3500 K to 50000 K is covered with an uneven grid (see Table A.3). The model spectra cover the ultraviolet (1000A) to infrared (10 microns) spectral range with non-uniform wavelength spacing (see Table A.4).

Table A.3: Grid of Temperatures for the Castelli-Kurucz Models
Temperature Range Grid Step
3500 - 13000 250
13000 - 50000 1000

Table A.4: Wavelength coverage for the Castelli-Kurucz Models
Wavelength Range Grid Step
microns A
0.09 - 0.29 10
0.29 - 1.00 20
1.00 - 1.60 50
1.60 - 3.20 100
3.20 - 6.4 200
6.4 - 10.0 400

A.3.1 The HST/CDBS Version of the Castelli and Kurucz models

The new atlas is divided in 8 independent subdirectories, according to metallicity (the CDBS tables were created in January 2007). Within each subdirectory the stellar atmosphere models are given in FITS multicolumn table format. Each table consists of 12 columns, the first one containing the wavelength grid and each of the rest containing the spectrum of a star with the same effective temperature but different gravity, ranging from log_g= 0.0 to +5.0. Columns filled with zeros indicate that the model spectrum for that particular metallicity, effective temperature and gravity combination is not covered by the atlas.

The names of the tables are given as ckszz_ttttt.fits where "ck" stand for Castelli & Kurucz; "szz" is the metallicity ([M/H] or log_Z) of the model (zz) with its sign (s); and "ttttt" is the model effective temperature (Teff), using four or five digits depending on the value. For instance, models for an effective temperature of 5000 K with [M/H]= -0.5 and [M/H]= +3.5 are indicated by ttttt= 5000, s= m, zz= 05 and ttttt= 5000, s= p, zz= 35, respectively, and correspond to the files ckm05_5000.fits and ckp35_5000.fits.

Within each individual table file, each column is named "gyy" where "yy" corresponds to 10*log_g. For example, log_g= +0.5 and log_g= +4.0 models are located in columns named g05 and g40, respectively.

In Table A.5, an example of a standard header for the file ckp00_8000.fits which gathers all calculated models for a star of metallicity log_Z= 0.0 and effective temperature Teff= 8000 K is given. Models cover a range of gravities from log_g= +1.0 (g10 in the header) to log_g= +5.0 (g50 in the header). Models for gravities log_g= +0.0 and +0.5 are not available for this particular metallicity and effective temperature combination, and therefore are not listed in the header.Their corresponding columns (g00 and g05) are filled with zeros.

Table A.5: Sample FITS Header for Castelli and Kurucz Model

1 SIMPLE b T Fits standard

2 BITPIX i 16 Bits per pixel

3 NAXIS i 0 Number of axes

4 EXTEND b T File may contain extensions

5 DATE t '2007-01-09T18:55:35' Date FITS file was generated

6 IRAF-TLM t '13:55:47 (09/01/2007)' Time of last modification

7 COMMENT t FITS (Flexible Image Transport System) format is defined in 'Astronomy

8 COMMENT t and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H

9 ORIGIN t 'STScI-STSDAS/TABLES' Tables version 2002-02-22

10 FILENAME t 'ckp00_8000.fits' name of file

11 NEXTEND i 1 number of extensions in file

12 TEFF i 8000

13 LOG_Z d 0.

14 HISTORY t File created by F.R.Boffi

15 HISTORY t ATLAS9 model atmospheres by Castelli and Kurucz (2004).

16 HISTORY t Wavelength is in Angstrom.

17 HISTORY t Fluxes tabulated in units of erg/s/cm^2/A

18 HISTORY t (after converting original units into "flam",

19 HISTORY t as described in README file and the SYNPHOT manual)

20 HISTORY t and are surface fluxes. To transform to observed

21 HISTORY t fluxes multiply by (R/D)^2 where R is the

22 HISTORY t radius of the star and D the distance.

23 HISTORY t Each column in the table represents the

24 HISTORY t spectrum of a star for the same metallicity

25 HISTORY t and effective temperature but different gravity.

26 HISTORY t Gravities range from log_g = +0.0 (g00 in the column

27 HISTORY t header) to log_g = +5.0 (g50 in the column header).

As indicated in the header file, physical fluxes of the spectra are given in FLAM surface flux units, i.e. ergs cm-2 s-1 -1. These flux units differ from those in the Castelli & Kurucz tables by a factor of (3.336 x 10-19 x -2 x (4)-1)-1, i.e. are converted from the original units of ergs cm-2 s-1 Hz-1 steradian-1 to ergs cm-2 s-1 -1 by multiplying the Castelli & Kurucz values by (3.336 x 10-19 x 2 x (4)-1)-1, where lambda is in Angstrom units. To convert to observed flux at Earth, multiply by a factor of (R/D)2 where R is the stellar radius, and D is the distance to Earth.

The files located in each metallicity subdirectory are listed in a file located in the main subdirectory called catalog.fits: for a given temperature a number of models are available for different gravity values.

In the following an excerpt of the file is provided (Table A.6). The first column provides temperature, metallicity, and gravity value, while the second column lists the corresponding file with the gravity column in square brackets. In this case, the two sets of models belong to the same directory, ckm05, for a [M/H]= -0.5, for two temperatures (10000 and 10250 respectively), and for gravities between log_g= +0.0 and log= +5.0. By attaching to the filename a gravity value between square brackets a specific column, with flux values for that specific gravity value, is going to be read and used by synphot. See next paragraph about the use of the catalog with Synphot tasks.

Table A.6: Excerpt from catalog.fits in crgrid$ck04models/
10000,-0.5,0.0 ckm05/ckm05_10000.fits[g00]
10000,-0.5,0.5 ckm05/ckm05_10000.fits[g05]
10000,-0.5,1.0 ckm05/ckm05_10000.fits[g10]
10000,-0.5,1.5 ckm05/ckm05_10000.fits[g15]
10000,-0.5,2.0 ckm05/ckm05_10000.fits[g20]
10000,-0.5,2.5 ckm05/ckm05_10000.fits[g25]
10000,-0.5,3.0 ckm05/ckm05_10000.fits[g30]
10000,-0.5,3.5 ckm05/ckm05_10000.fits[g35]
10000,-0.5,4.0 ckm05/ckm05_10000.fits[g40]
10000,-0.5,4.5 ckm05/ckm05_10000.fits[g45]
10000,-0.5,5.0 ckm05/ckm05_10000.fits[g50]
10250,-0.5,0.0 ckm05/ckm05_10250.fits[g00]
10250,-0.5,0.5 ckm05/ckm05_10250.fits[g05]
10250,-0.5,1.0 ckm05/ckm05_10250.fits[g10]
10250,-0.5,1.5 ckm05/ckm05_10250.fits[g15]
10250,-0.5,2.0 ckm05/ckm05_10250.fits[g20]
10250,-0.5,2.5 ckm05/ckm05_10250.fits[g25]
10250,-0.5,3.0 ckm05/ckm05_10250.fits[g30]
10250,-0.5,3.5 ckm05/ckm05_10250.fits[g35]
10250,-0.5,4.0 ckm05/ckm05_10250.fits[g40]
10250,-0.5,4.5 ckm05/ckm05_10250.fits[g45]
10250,-0.5,5.0 ckm05/ckm05_10250.fits[g50]


A.3.2 Use of the Castelli and Kurucz models with Synphot

Synphot tasks permit the use of spectra selected from one of many columns in a single FITS file. One does this by specifying as the "spectrum" parameter the name of the disk file (as before), and appending the name of the column containing the flux in brackets. Thus, to select any model spectrum characterized by a given metallicity, effective temperature, and gravity, one needs to specify a "spectrum" of the form: crgridck04$m_directory/ckszz_ttttt.fits[gyy], where m_directory is the name of the subdirectory for a given metallicity, szz, ttttt, and gyy are as above. For example, to select the spectrum of a star with a metallicity of +0.1, a temperature of 10,000 K, and log_g of 3.0, the specification would be: crgridck04$ckp01/ckp01_10000.fits[g30].

Please note that the model spectra in the atlas are in surface flux units. Thus, if the number of counts or the calculated absolute flux is needed, the model spectrum must be renormalized appropriately. One can do this in synphot with the "rn" function.

Synphot also allows the use of the cat() and icat() functions to select Castelli & Kurucz spectra. The syntax is "cat(ck04models,t,m,g) (or icat(ck04models,t,m,g)) where "t" is the effective temperature, "m" is the log metallicity [M/H], and "g" is the log_g.

The icat task can give incorrect results without warning if it is used at the edges of the defined grid of models. See above for a discussion of the available ranges for these values, and exercise caution when selecting a value near one of the extremes.

Since the entire atlas is very large, and many cases can be explored by use of the solar metallicity models, only these are made available.

A list of solar metallicity stars of different spectral types and luminosity classes together with their closest Castelli & Kurucz 2004 model spectrum is presented in Table A.7. The physical parameters, Teff and log_g, characterizing each O stars are taken from Martins, Schaerer, & Hiller's (2005) compilation of stellar parameters of Galactic O stars. The physical parameters for later stars are taken from Schmidt-Kaler's (1982) compilation of physical parameters of stars); for these, the U-B and B-V colors of the closest model agree with the characteristic color of each star (see Schmidt-Kaler 1982) to better than 0.06 magnitude.

Table A.7: Suggested Models for Specific Stellar Types
Type Teff log_g Castelli & Kurucz model
O3V 44852 +3.92 ckp00_45000[g45]
O4V 42857 +3.92 ckp00_43000[g45]
O5V 40862 +3.92 ckp00_41000[g45]
O5.5V 39865 +3.92 ckp00_40000[g40]
O6V 38867 +3.92 ckp00_39000[g40]
O6.5V 37870 +3.92 ckp00_38000[g40]
O7V 36872 +3.92 ckp00_37000[g40]
O7.5V 35874 +3.92 ckp00_36000[g40]
O8V 34877 +3.92 ckp00_35000[g40]
O8.5 33879 +3.92 ckp00_34000[g40]
O9V 32882 +3.92 ckp00_33000[g40]
O9.5 31884 +3.92 ckp00_32000[g40]
B0V 30000 +3.90 ckp00_30000[g40]
B1V 25400 +3.90 ckp00_25000[g40]
B3V 18700 +3.94 ckp00_19000[g40]
B5V 15400 +4.04 ckp00_15000[g40]
B8V 11900 +4.04 ckp00_12000[g40]
A0V 9520 +4.14 ckp00_9500[g40]
A1V 9230 +4.10 ckp00_9250[g40]
A3V 8270 +4.20 ckp00_8250[g40]
A5V 8200 +4.29 ckp00_8250[g40]
F0V 7200 +4.34 ckp00_7250[g40]
F2V 6890 +4.34 ckp00_7000[g40]
F5V 6440 +4.34 ckp00_6500[g40]
F8V 6200 +4.40 ckp00_6250[g45]
G0V 6030 +4.39 ckp00_6000[g45]
G2V 5860 +4.40 ckp00_5750[g45]
G5V 5770 +4.49 ckp00_5750[g45]
G8V 5570 +4.50 ckp00_5500[g45]
K0V 5250 +4.49 ckp00_5250[g45]
K2V 4780 +4.5 ckp00_4750[g45]
K4V 4560 +4.5 ckp00_4500[g45]
K5V 4350 +4.54 ckp00_4250[g45]
K7V 4060 +4.5 ckp00_4000[g45]
M0V 3850 +4.59 ckp00_3750[g45]
M2V 3580 +4.64 ckp00_3500[g45]
M4V 3370 +4.80 ckp00_3500[g50]
M5V 3240 +4.94 ckp00_3500[g50]
M6V 3050 +5.00 ckp00_3500[g50]
B0III 29000 +3.34 ckp00_29000[g35]
B5III 15000 +3.49 ckp00_15000[g35]
G0III 5850 +2.94 ckp00_5750[g30]
G5III 5150 +2.54 ckp00_5250[g25]
K0III 4750 +2.14 ckp00_4750[g20]
K5III 3950 +1.74 ckp00_4000[g15]
M0III 3800 +1.34 ckp00_3750[g15]
O5I 40300 +3.34 ckp00_40000[g45]
O6I 39000 +3.24 ckp00_39000[g40]
O8I 34200 +3.24 ckp00_34000[g40]
BOI 26000 +2.84 ckp00_26000[g30]
B5I 13600 +2.44 ckp00_14000[g25]
AOI 9730 +2.14 ckp00_9750[g20]
A5I 8510 +2.04 ckp00_8500[g20]
F0I 7700 +1.74 ckp00_7750[g20]
F5I 6900 +1.44 ckp00_7000[g15]
G0I 5550 +1.34 ckp00_5500[g15]
G5I 4850 +1.14 ckp00_4750[g10]
K0I 4420 +0.94 ckp00_4500[g10]
K5I 3850 +0.34 ckp00_3750[g05]
M0I 3650 +0.14 ckp00_3750[g00]
M2I 3450 -0.06 ckp00_3500[g00]

Space Telescope Science Institute
Voice: (410) 338-1082