And Log G Determinations: Barry Smalley
And Log G Determinations: Barry Smalley
And Log G Determinations: Barry Smalley
8, 130
Memorie della
c SAIt 2005
Supplementi
example, the HIPPARCOS parallax catalogue ibration with uses Vega as a zero-point must
(ESA 1997), or the membership of a cluster assume a value for log g. However, detailed
with a known distance, provided that distance model atmosphere analyses give a value of
has not been obtained using model-dependent log g= 3.95 ± 0.05 (Castelli & Kurucz 1994).
methods. An interesting discussion on the accuracy
Single star mass determination is exceed- of the visible and near-infrared absolute flux
ingly difficult, with microlensing the only calibrations is given by Mégessier (1995).
known direct method (Alcock et al. 2001; These uncertainties place a limit on our cur-
Jiang et al. 2004). This relies on chance align- rent direct determinations of stellar fundamen-
ments and is considerably less accurate than tal parameters.
that possible with eclipsing binary systems.
3. Indirect Methods
2.4. Accuracy of Direct Measurements The direct determination of T eff and log g is not
2.4.1. Sun possible for most stars. Hence, we have to use
indirect methods. In this section we discuss the
Our nearly stellar companion, the Sun, has the use of various techniques used to determine the
most accurately known stellar parameters. The atmospheric parameters.
measured total solar flux at the earth, the Solar When determining T eff and log g, using
Constant, is f = 1367 ± 4 W m−2 (Mendoza model-dependent techniques, we must not ne-
2005). Variations due to the Solar Cycle and glect metallicity ([M/H]). An incorrect metal-
rotation, contribute 0.1% and 0.2%, respec- licity can have a significant effect on perceived
tively (Zahid et al. 2004). This equates to ± values of these parameters.
4 K in the Solar effective temperature. A value
of T eff = 5777 ± 10 K is obtained from the
Solar Constant and the measured Solar radius, 3.1. Photometric Grid Calibrations
including calibration uncertainties. The Solar There have been many photometric systems
surface gravity is exceedingly well known; developed to describe the shape of stellar
log g = 4.4374 ± 0.0005 (Gray 1992). flux distributions via magnitude (colour) dif-
ferences. Since they use wide band passes ob-
2.4.2. Vega servations can be obtained in a fraction of the
time required by spectrophotometry and can be
The bright star Vega is our primary stellar flux extended to much fainter magnitudes. The use
calibrator (Hayes & Latham 1975; Bohlin & of standardized filter sets allows for the quan-
Gilliland 2004). The measured total flux at titative analysis of stars over a wide magnitude
the earth is f⊕ = 29.83 ± 1.20 × 10−9 W m−2 range.
(Alonso et al. 1994), which is an uncertainty Theoretical photometric indices from
of some 4%. There have been reports that Vega ATLAS flux calculations are normalized using
may be variable (Fernie 1981; Vasil’yev et al. the observed colours and known atmospheric
1989), but these have not been substantiated, parameters of Vega. Vega was originally
and may well be spurious. Nevertheless, this is chosen because it is the primary flux standard
something that ought to be investigated. Using with the highest accuracy spectrophotometry.
the interferometric angular diameter of Ciardi An alternative, semi-empirical approach, is
et al. (2001), θ = 3.223 ± 0.008 mas, we ob- to adjust the theoretical photometry to min-
tain T eff = 9640 ± 100 K. Most of the uncer- imize discrepancy with observations of stars
tainty (∼95K) is due to the uncertainties in the with known parameters. Moon & Dworetsky
measured fluxes, while the error in the angular (1985) used stars with fundamental values
diameter only contributes ∼10K. to shift the grids in order to reduce the dis-
Since Vega is a single star, there is no direct crepancy between the observed and predicted
fundamental log g measurement. Thus any cal- colours. In contrast, Lester, Gray & Kurucz
Smalley: T eff and log g determinations 133
(1986) treated the raw model colours in the calculations (Blackwell & Lynas-Gray 1994;
same manner as raw stellar photometry. The Mégessier 1994):
model colours were placed on the standard
system using the usual relations of photomet- f⊕ F∗ σT eff 4
ric transformation. However, both these ap- ≡ =
fλ0 F λ0 φ(T eff , log g, λ0 )
proached have the potential to mask physical
problems with models. The method requires a complete flux dis-
Overall, photometry can give very good tribution in order to obtain the total integrated
first estimates of atmospheric parameters. In ( f⊕ ) stellar flux. In practice, however, all of
the absence of any other suitable observations, the flux is not observable, especially in the far-
the values obtainable from photometry are of ultraviolet. But, this is only a serious problem
sufficient accuracy for most purposes, with in the hottest stars, where model atmospheres
typical uncertainties of ±200 K and ±0.2 dex can be used to insert the missing flux, in or-
in T eff and log g, respectively. der to obtain the total integrated flux. Accurate
infrared fluxes are, of course, essential for this
method to produce reliable results.
3.2. T eff –colour Relationships The method is sensitive to the presence of
Effective temperatures can be estimated from any cooler companion stars. The effect of the
photometric colour indices. Empirical calibra- companion is to lower the T eff derived for the
tions are based on stars with known tem- primary. A modified method was proposed and
peratures, often obtained using the IRFM. discussed by Smalley (1993). This method re-
There are many examples in the literature, lies on the relative radii of the two components
for example, (Alonso et al. 1996; Houdashelt in the binary system. The effect of allowing for
et al. 2000; Sekiguchi & Fukugita 2000; the companion can be dramatic; the T eff de-
VandenBerg & Clem 2003; Clem et al. 2004; termined for the primary can be increased by
200 K or more.
Ramírez & Meléndez 2005b).
A very useful by-product of the IRFM is
Particularly useful are V − K calibrations, that it also gives the angular diameter (θ) of the
since this index is much less sensitive to metal- star.
licity than B − V (Alonso et al. 1996; Kinman Given good spectrophotometry, the IRFM
& Castelli 2002; Ramírez & Meléndez should give estimates of T eff , which are clos-
2005b). However, this index is more sensitive est to the ‘true’ fundamental value. In fact it
to the presence of a cool companion. has been used as the basis of other calibrations
Often, there are several steps involved in (e.g. Ramírez & Meléndez 2005a). Typically
obtaining the calibrations. The uncertainties we can obtain temperatures to an accuracy of
and final error on the parameters obtained are 1∼2% (Blackwell et al. 1990). The IRFM re-
not always immediately obvious. sults for Vega have an uncertainty of ∼150K.
Uncertainties in absolute calibration of IR
photometry are important. For example, for
3.3. InfraRed Flux Method 2MASS an error of ∼50K, for a T eff of 6500K,
The InfraRed Flux Method (IRFM), devel- arises from the uncertainty in the absolute cal-
oped by Blackwell & Shallis (1977) and ibration alone.
Blackwell, Petford & Shallis (1980), can be
used to determine T eff . The method relies on 3.4. Flux Fitting
the fact that the stellar surface flux at an in-
frared wavelength (λ0 ) is relatively insensitive The emergent flux distribution of a star is re-
to temperature. The method is almost model lated to its atmospheric parameters. We can
independent (hence near fundamental), with use spectrophotometry to determine values for
only the infrared flux at the stellar surface, these parameters, by fitting model atmosphere
φ(T eff , log g, λ0 ), requiring the use of model fluxes to the observations. Figure 2 shows the
134 Smalley: T eff and log g determinations
sensitivity of the flux distribution to the var- agree. This gives a line in a T eff – log g di-
ious atmospheric parameters. However, inter- agram.
stellar reddening must be allowed for, since it Excitation Potential
can have a significant effect on the observed Abundances from the same element and
flux distribution and derived parameters. ionization stage should agree for all exci-
The currently available optical flux dis- tation potentials
tributions need are in need of revision. This Microturbulence
is something that will be done by ASTRA The same abundance of an element should
(Adelman et al. 2005). be obtained irrespective of the equivalent
width of the lines. This is the technique
used to obtain the mictroturbulence param-
3.5. Balmer Profiles eter (ξturb ). See Magain (1984) for a dis-
The Balmer lines provide an excellent T eff di- cussion of the systematic errors in micro-
agnostic for stars cooler than about 8000 K due turbulence determinations. Typically, can
to their virtually nil gravity dependence (Gray expect to get ξturb to no better than ± 0.1
1992; Heiter et al. 2002). By fitting these the- km s−1 .
oretical profiles to observations, we can deter-
mine T eff . For stars hotter than 8000 K, how-
ever, the profiles are sensitive to both temper- Using these techniques it is possible to get
ature and gravity. For these stars, the Balmer a self-consistent determination of a star’s at-
lines can be used to obtain values of log g, pro- mospheric parameters.
vided that the T eff can be determined from a
different method.
3.8. Global Spectral Fitting
3.6. Spectral Line Ratios
An alternative to a detailed analysis of individ-
Spectral lines are sensitive to temperature vari- ual spectral line measurements, is to use the
ations within the line-forming regions. Line whole of the observed stellar spectrum and find
strength ratios can be used as temperature di- the best-fitting synthetic spectrum. The normal
agnostics, similar to their use in spectral clas- procedure is to take a large multi-dimensional
sification. Gray & Johanson (1991) used line grid of synthetic spectra computed with var-
depth ratios to determine stellar effective tem- ious combinations of T eff , log g, ξturb , [M/H]
peratures with a precision of ±10 K. While and locate the best-fitting solution by least
this method can yield very precise relative squares techniques.
temperatures, the absolute calibration on to
The benefit of this method is that it can be
the T eff scale is much less well determined
automated for vast quantities of stellar obser-
(Gray 1994). This method is ideal for inves-
vations and that it can be used for spectra that
tigating stellar temperature variations (Gray &
are severely blended due to low resolution or
Livingston 1997).
rapid rotation.
Naturally, the final parameters are model
3.7. Metal Line Diagnostics dependent and only as good as the quality of
In a detailed spectral analysis, the equivalent the model atmospheres used. The internal fit-
width of many lines are often measured. These ting error only gives a measure of the precision
can be used to determine the atmospheric pa- of the result and is thus a lower limit uncer-
rameters via metal line diagnostics. tainty of the parameters on the absolute scale.
Determination of the accuracy of the parame-
Ionization Balance ters requires the assessment of the results of fit-
The abundances obtained from differing ting, using the exact same methods, to spectra
ionization stages of the same element must of fundamental stars.
136 Smalley: T eff and log g determinations
Fig. 2. The sensitivity of flux distributions to T eff , log g and [M/H]. The base model (T eff = 7500, log g =
4.0, [M/H] = 0.0) is represented by a solid line. The dotted and dashed lines indicate models with one of
the parameters adjusted, as indicated.
Smalley: T eff and log g determinations 137
Fig. 3. Sensitivity of Balmer profiles to [M/H] and v sin i. The synthetic spectra (solid line) have been
calculated with T eff = 7500 and log g = 4.0 and a simulated resolution of around 0.4Å. The true shape of
the Hγ profile (with no metal lines) is shown as the dotted line.
138 Smalley: T eff and log g determinations
Fig. 4. A T eff –log g diagram for 63 Tau. The results from four methods are shown as follows: the filled
square is from the Moon & Dworetsky (1985) uvbyβ girds, the filled circle is from spectrophotometric flux
fitting, the dashed line that from fitting Hβ profiles and the dotted line the IRFM result. Photometry and
Balmer lines agree very well, but are significantly hotter than the results from Spectrophotometry and the
IRFM. The solid arrows indicate the effect of using [M/H] = +0.5 models. Now spectrophotometry is in
good agreement with photometry and Balmer line, but the IRFM is still significantly lower. However, by
introducing a cool companion (5000 K) the IRFM can be brought into agreement with the other methods
(dotted arrow). The solid line is the Hyades isochrone, based on the evolutionary calculations of Schaller et
al. (1992). (Adapted from Smalley 1996)
Moon T.T., Dworetsky M.M., 1985, MNRAS, Smalley B., 1996, ASP Conf, 108, 43
217, 305 Smalley B., Dworetsky M.M., 1995, A&A,
Mozurkewich D., Armstrong J.T., Hindsley 293, 446
R.B., Quirrenbach A., Hummel C.A. et al., Smalley B., Kupka F., 1997, A&A, 328, 349
2003, AJ, 126, 2502 Smalley B., Gardiner R.B., Kupka F., Bessell
Pasinetti-Fracassini L.E., Pastori L., Covino S., M.S., 2002, A&A, 395, 601
Pozzi A., 2001, A&A, 367, 521 [II/224] Southworth J., Maxted P.F.L., Smalley B.,
Patience J., Ghez A.M., Reid I.N., Weinberger 2005, A&A, 429, 645
A.J., Matthews, K., 1998, AJ, 115, 1972 Steffen M., 1985, A&AS, 59, 403
Perevozkina E.L., Svechnikov M.A., 1999 Thompson G.I., Nandy K., Jamar C., Monfils
[V/118] A., Houziaux L., Carnochan D.J., Wilson R.,
Popper D.M., 1980, ARA&A, 18, 115 1978, The Science Research Council, U.K.
Ramírez I., Meléndez J., 2005a, ApJ, 626, [II/59B]
4446 VandenBerg D.A, Clem J.L., 2003, AJ, 126,
Ramírez I., Meléndez J., 2005b, ApJ, 626, 465 778
Richichi A., Percheron I., Khristoforova M., Vasil’Yev I.A., Merezhin V.P., Nalimov V.N.,
A&A, 431, 773 Novosyolov V.A., 1989, IBVS, 3308, 1
Schaller G., Schaerer D., Meynet G., Maeder Verhoelst T., Bordé P.J., Perrin G., Decin L.,
A., 1992, A&AS, 96, 269 Eriksson K., et al., 2005, A&A, 435, 289
Sekiguchi M., Fukugita M., 2000, AJ, 120, Zahid H.J., Hudson H.S., Fröhlich C., 2004,
1072 Solar Physics, 222, 1
Smalley B., 1993, MNRAS, 265, 1035