A&A 366, 241-253 (2001)
P. A. Zaal1 - A. de Koter1 - L. B. F. M. Waters1,2
1 - Astronomical Institute Anton Pannekoek, University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands
2 - Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, Celestijnenlaan 200 B, 3001 Heverlee, Belgium
Received 4 November 1999 / Accepted 6 November 2000
We present observations of the Br, Br and Pf lines of 16 dwarf and (sub)giant stars in the spectral range O9-B3. The observations were done using the Short Wavelength Spectrometer on board the Infrared Space Observatory, and have a signal-to-noise of 20 to >150 and a resolving power varying from 1400 to 2100. We compare the equivalent widths of these lines with predictions using non-LTE model atmospheres to investigate to what extent these infrared lines can be used to derive effective temperatures. We find that Pf is a sensitive diagnostic for the range of spectral types investigated, and Br for types O9-B2, yielding agreement with optical results to within 1-4 kK or one-three spectral sub-types. We find evidence for a gradient in the turbulent velocity, increasing from 5 kms-1 for the atmospheric region in which Br is formed to 15 kms-1 for the regime where Pf originates. When this gradient in turbulent velocity is taken into account, the accuracy of the spectral type calibration is improved to 1 kK or one spectral sub-type. The gravity dependence of the strengths of the investigated infrared lines is relatively weak, and could not be used to constrain luminosity class. This failure is in part a result of the modest S/N and resolution and in part a result of a cancelation of gravity effects in the line core and line wing. Our line predictions show that He I 2.058 is relatively unsuited for spectral classification of O9-B3 stars. Hu, however, is expected to be an even better diagnostic as are Pf and Br. This line may be observed with the VLT Imager and Spectrometer for mid InfraRed when it is installed on the Very Large Telescope.
Key words: line: formation - radiative transfer - stars: atmospheres - stars: early-type - stars: fundamental parameters - infrared: stars
In addition to the above mentioned broad perspective, the infrared offers an unique opportunity to study an important part of the atmosphere of these early-type stars, namely that of the outer photosphere. Due to the increasing importance of free-free processes towards longer wavelengths, the prominent infrared lines are formed in different regions of the outer photosphere. This allows one to study possible gradients in atmospheric properties, such as in the amount of (micro)turbulence, in this regime.
In this paper, we will study the behaviour of
a number of IR lines observed in 16 stars using the Short
Wavelength Spectrometer (SWS, de Graauw et al. 1996) on board the
Infrared Space Observatory (ISO, Kessler et al. 1996). The SWS covers
the wavelength range from 2.4 to 45 m. The principal goal is to
investigate to what extent the Br
4.05 (wavelength given in
7.46 and Pf
4.65 lines may be used for spectral type calibration. This requires one
to establish which processes are relevant and, therefore, which physics
needs to be included in the model calculations providing this
In O9 through B0.5 stars, Br and Pf are typically observed to be in emission. A grid of non-LTE plane-parallel hydrostatic models calculated in the spectral analysis study of the B0 V star Sco (Zaal et al. 1999; hereafter Paper I), shows that such profiles are expected and result from non-LTE effects in the line forming region. In B-type stars, already small departures from LTE are amplified by stimulated emission (Sigut & Lester 1996) facilitating infrared lines to go into emission (Paper I). Murdoch et al. (1994) also sought to explain the Br and Br line emission in 10 Lac (O9 V) and Sco by means of this non-LTE effect.
In this study, we will compare the grid of models from Paper I with the ISO spectra of our programme stars, in this way deriving effective temperatures based on infrared diagnostics alone. Although the models in this grid are relatively standard, i.e. they account for hydrogen and helium only and assume plane-parallel layers and hydrostatic equilibrium, they already show very promising results in reproducing the line strengths of Br, Br and Pf. We further investigate these results by addressing the presence of turbulent velocity fields; the possible effect of the onset of a stellar wind, and the effects of metal line blanketing on the structure of the atmosphere.
The paper is organized in the following way: the observations are described in Sect. 2. In Sect. 3, we compare observed and predicted line strengths using standard non-LTE model atmospheres. The influence of turbulence, a stellar wind and line blanketing is investigated in Sect. 4. In Sect. 5, we compare the diagnostic value of Br and Pf lines with other infrared diagnostics, especially with classifications using H- or K-band spectroscopy. Finally, Sect. 6 gives a summary of the main results.
|Figure 1: Normalized Br spectra ordered according to spectral type. Indicated are the H I and He I lines detected in at least one of the spectra. The thick horizontal lines at the outer edges of the spectral region represent the normalization level. The errorbar on the right indicates the noise per bin in the central part of the spectrum. The edges of the spectra show higher noise values, as these parts are not scanned by all detectors|
|Open with DEXTER|
|Figure 2: Normalized Br spectra ordered according to spectral type. Indicated are the H I and He I lines detected in at least one of the spectra. The errorbars are defined as in Fig. 1|
|Open with DEXTER|
|Figure 3: Normalized Pf spectra ordered according to spectral type. Indicated are the H I lines detected in at least one of the spectra. The errorbars are defined as in Fig. 1|
|Open with DEXTER|
The spectra were reduced using the SWS Interactive Analysis (IA3) programs. For the data reduction we made use of the algorithms developed by Valentijn & Thi (2000). A short summary of the data reduction method is given in Paper I. We optimized the signal-to-noise (S/N) ratio of the spectra using the interactive tool "LOW_SIGNAL'', as described by Lahuis et al. (1998). The SWS spectra suffered from fringing. The spectra were de-fringed using the fourier-filtering method in cases where the fringe was clearly present as a sharp peak in the fourier domain. Unfortunately, the Pf line observations were of insufficient S/N to be used in our analysis, and we discard these data in the remainder of this study. Br is observed for all programme stars. For the other lines we note that the dataset is not homogeneous, i.e. not all lines have been observed for all stars. For instance, the Pf profile in the late O-type stars has not been observed as these stars are relatively more distant and consequently have too low fluxes to obtain sufficient S/N. Figures 1-3 show Br, Br and Pf spectra respectively, arranged in order of progressively later spectral type. The observed line parameters for both the H I and the strongest He I lines are given in Table 2. The errorbar given in the figures represents the mean of the noise per bin in the central part of the spectrum. As the spectra are obtained in scanning mode the outer edges of the spectra show significantly lower S/N ratios compared to the central part of the scan. Only at Pf, several of the continua show fluctuations larger than the given noise level. We attribute these fluctuations to fringes which, because of insufficient S/N, could not be filtered out.
Using the ISO Spectral Analysis Package (ISAP), we measured the line strength and width for all detected hydrogen and helium lines. The results are presented in Table 2. The listed profile parameters follow from multiple Gaussian fits to all identified line profiles within the observed spectral range. We distinguished the line absorption from the line emission in case both are clearly present. The derived line fluxes and equivalent widths for Br are lower limits, as the continuum is difficult to determine for this broad absorption line, especially for the late O- and earliest B-type stars (see Fig. 1).
All three hydrogen lines show a decrease in line core emission going
from early to late spectral types, where at some point the core reverts
to absorption. The line emission is strongest for Pf
for which emission is detected down to spectral type B1 and B0
respectively. For Br
weak central emission is only found in the late
O-type stars and in Sco. The line core emission is not fully
resolved and one should realize that in high resolution spectra core
emission is typically present in Br
down to spectral type
B3. For example, Zaal et al. (1997) showed that high resolution UKIRT
spectra of two slowly rotating stars, Her (B3IV)
and Oph (B2IV) still show weak core emission.
|Pf||Not suff. S/N||20|
|Br||Not suff. S/N||30|
|Pf||Not suff. S/N|
|Br||Not suff. S/N||20|
|Br||a||He I 1s6f-1s4d (S=3)a||26192.1|
|26258.7||b||He I 1s6f-1s4d (S=1)b||26205.6|
|c||He I 1s6g-1s4f (S=1,3)c||26240.9|
|d||He I 1s5f-1s4d (S=3)d||40377.3|
|Br||e||He I 1s5f-1s4d (S=1)e||40409.3|
|40522.6||f||He I 1s5g-1s4f (S=1)f||40490.0|
|g||He I 1s5g-1s4f (S=1,3)g||40490.3|
As mentioned, the observed spectra also show the presence of several He I lines. In some cases these lines blend with the observed H I profile and affect the shape of these profiles. The most prominent one is He I 4.049 situated at -240 kms-1 relative to Br. It shows weak emission for the later spectral type stars: Cen, Peg and Pav (and possibly Dra), causing a strong asymmetry in the Br absorption profile (see Fig. 2). Zaal et al. (1997) also found He I 4.049 emission for Her (B3IV) and Oph (B2IV). This suggests that it shows significant emission for stars of spectral type B1-B3. As in the case of Her the measured velocity of the He I 4.049 line is small (i.e. 5 kms-1 relative to the center of mass of the star), its origin is most likely photospheric. He I 4.049 is also suspected to be in emission for Cir, 10 Lac, Sco, and Car. We conclude that the line asymmetry in Br for stars with spectral type B1-B3 can qualitatively be explained by line emission in the He I 4.049 line. We also suspect that the weak asymmetry in the Br line profile of Peg and Pav can be attributed to weak emission in He I 1s6g-1s4f at -203 kms-1 relative to Br. However, the blue-shifted peak-velocity of the Br absorption in Pav (-65 kms-1, see Fig. 1) can not be explained by this He I blend and may be related to its binary nature.
Other weak He I lines in the profiles of Br
only show weak
absorption, possibly with the exception of He I 1s5f-1s4d (at -1075
kms-1 relative to Br)
which may be in emission in case of Sco.
Although it is hard to define a
proper continuum in the region of Br,
the observations indicate that
the blue absorption wing of this line is more enhanced for the earliest
spectral types in our sample. We attribute this asymmetry to the
presence of weak He I absorption lines. For late B-type stars, such as
Dra, the absorption in these
He I lines is much weaker, resulting in more symmetric Br
(see Fig. 1). Within the spectral range of Pf
only detect the H I 8-6 and H I 11-7 line blends in the spectrum of
Cru. In both stars these two lines are, like
|Figure 4: The observed equivalent width of Br (square), Br (circle) and Pf (triangle) profiles versus . The EWs derived from non-blanketed model calculations are given for three different values: 4.3 (solid), 4.0 (dotted) and 3.7 (dashed). Both observed and predicted equivalent widths were measured over a wavelength of kms-1 relative to the rest wavelength of each line. We included also four Br EWs (small symbol size) obtained with UKIRT (Zaal et al. 1997). The adopted effective temperatures are derived from optical analyses (references are given in Table 5)|
|Open with DEXTER|
The presence of He I absorption and/or emission can not explain all
asymmetries. For the two O stars,
Cir and 10 Lac, and
for Sco the spectra show additional asymmetric features. In case of
10 Lac, for instance, the central reversal in Br
is blue-shifted to
-54 kms-1. Also the peak velocity of Br
emission of these three
stars is significantly blue-shifted. The Br
spectrum of 10
Cir even show additional asymmetric features.
In both stars a redshifted absorption feature is present at
respectively +200 and +300 kms-1 (see Fig. 2).
Zaal et al. (1997) found the same asymmetries in Br
for 10 Lac and
Sco. In case of 10 Lac they also found blue-shifted line core
emission, similar to that found here for Br.
These asymmetries can not
be explained by contamination of He I lines and must be intrinsic to
the velocity field in the line forming region. The fact that these
asymmetries occur only for stars with the earliest spectral types in our
sample may suggest that the H I lines are affected by a stellar wind.
However, in Paper I we demonstrated that for Sco the line
displacements of the observed H I lines do not seem to fit to the
scenario of an accelerating outflow. If these lines are affected by a
stellar wind one would expect to find a correlation between depth of
formation and line displacement, which is not present for Sco.
The observed line displacements for 10 Lac point to the
same conclusion. Possibly, the observed asymmetries in the H I infrared lines of these early-type stars are due to (macro) turbulent
motions in the outer photospheric region.
|= 0 kms-1|
|= 25 kms-1|
A comparison of observed and predicted values demonstrates that the
EWs agree reasonably well. One can clearly see that each individual
line is temperature sensitive only in a limited range. For Pf
is strongest, showing a decrement of more than 17
Å going from 34 to 20 kK. For Br
the dependence on
also significant, showing a decrement of about 8 Å going from 34 to
26 kK. For Br,
however, the dependence on temperature is rather weak.
Note that between 17 and 27 kK the observed Pf
strengths seem to be systematically lower than predicted by our
models. We will return to this discrepancy below. In general, the
sensitivity of the line strengths to gravity appears to be weak. Even
when the equivalent width is determined over a larger wavelength
interval, i.e. including the more gravity sensitive line wings, we
found that the
sensitivity to this parameter remains weak. For temperatures larger
than 34 kK, the EWs of Br
do show some dependence on
although in view of the errorbars on the observed EWs (see
Fig. 4) this sensitivity is too weak to be useful. The
same holds for the EW of Pf
We now apply the model results in Fig. 4 to derive effective temperatures from the infrared. Table 5 presents these derived temperatures and compares them with and values from optical analyses. In those regimes where the observed EW changes rapidly with temperature, the diagnostic can be used to derive fairly accurate values. The range in spectral types for which this can be done is largest for Pf. This line may be used for spectral types between O9 and B2, i.e. for stars with temperatures between 33 20 kK. For Br the interval is somewhat smaller, extending from O9 to B1. The Br line shows the weakest dependence on temperature. The temperature interval in which it can be applied is similar to that for Br. As a result of a rather poor continuum determination, yielding relatively large errors in equivalent width, we could only apply Br for the B0 V star Sco. The line yields upper limits for several other programme stars.
How accurate are the derived temperatures? The absolute flux level of Pf is about a factor of 3 to 4 lower compared to the flux at the Br wavelength. Also, the SWS sensitivity at Pf (observed in band 2) is significantly lower compared to the sensitivity at Br. These differences imply that the error in the line strength determination of Pf will be much larger than the error estimate for Br. This cancels out the larger predicted sensitivity of Pf. In practice the reliability of both lines as a diagnostic is about equal in the range in which they both may be applied and amounts to an intrinsic error 1-2 kK.
As noted above, a discrepancy exists between observed and predicted Pf
strength in the range
kK (corresponding to
spectral types B0.5-B2). Here the predicted EW is systematically
higher by about 1-2 Å, which relates to an infrared temperature
about 1-2 kK higher than the temperatures derived from the optical.
We will return to this systematic shift in Sect. 4, where we will
show that it is connected with the neglect of turbulence in the set of
models shown in Fig. 4.
|Figure 5: The equivalent width for a) Hu and b) He I 1.083 m as derived from the non-LTE models for three different values: 3.7 (dashed), 4.0 (dotted) and 4.3 (solid)|
|Open with DEXTER|
The above implies that the EW of Br may be used directly for spectral classification. We find that the line allows for an identification to within about one subclass in the spectral range O9-B0.5. The Pf line yields an accuracy of one to three subtypes within the range B0.5-B2. However, this last error may be narrowed down to about 1 to 1.5 subclasses when turbulence is included (see Sect. 4). The equivalent widths of Br as obtained with [CGS4] UKIRT (Zaal et al. 1997) have errorbars about two times smaller than those obtained with SWS (see Fig. 4). If one can secure spectra of similar quality, i.e. spectra with resolution R 10000 and signal-to-noise S/N >50, the Br line may even be suited to derive the spectral type to within half a subclass for stars within the spectral range O9-B1.
We may also use our models to predict to what extent other infrared lines, not observable with ISO, may be used as a temperature diagnostic. First we discuss Hu at 12.4 m. This particular line may be observed using VISIR (VLT Image and Spectrometer for mid InfraRed), scheduled to be installed in 2001 on the Very Large Telescope (VLT). With this instrument it should be possible to obtain Hu spectra with sufficient S/N for most of our program stars. Figure 5a shows the predicted EW of this line as a function of temperature. Compared to Pf, Hu shows an even larger sensitivity to extending down to temperatures as low as 17 kK (spectral type B3). We conclude that this line is expected to be an excellent infrared diagnostic for O9 to B3 type stars. Next, we investigate to what extent He I infrared lines may be used as temperature diagnostics. Figure 5b shows the equivalent width behaviour for He I 1.083 m. The gravity sensitivity of this line is like that of the hydrogen lines rather weak. Unfortunately also the dependence of this He I line is weak, rendering the line almost useless for determination unless very high quality spectra can be secured ( ). A second prominent He I line in the infrared is He I 2.058. The strength of this line is, however, strongly coupled to the strength of the UV continuum through the resonance transition at 584 Å (Morris et al. 1996). This makes this line unsuited for spectral classification. We conclude that helium infrared lines are not expected to yield a sensitive temperature calibration method for stars of spectral type O9-B3, unless very high S/N data is available for He I 1.083 m.
The predicted equivalent widths discussed in the previous section are based on non-LTE hydrostatic H&He model calculations. In this section we will test to what extent some of the assumptions in these models will affect the temperature determination. First, we will compare the observed line strengths with LTE calculations to demonstrate that non-LTE is a crucial physical ingredient for modeling the observed EWs. Secondly, we will show that when accounting for line blanketing, one does not expect to affect the predicted line strength significantly. Thirdly, we describe the effect of turbulence. By including a turbulent velocity of 15 kms-1 in the model calculations we are able to explain the mismatch between the observed Pf EW and those predicted. Finally, we will discuss whether a stellar wind may be expected to have an effect on the H I and He I infrared lines.
Figure 6a shows the observed EWs together with model predictions from LTE calculations. Only one gravity ( = 4.0) was used. The predicted values for Pf are strikingly different from the non-LTE results for all of the investigated range. Large discrepancies also occur for Br at kK. At kK the LTE assumption seems to be reasonable for both Br and Br. However, close inspection of the predicted profiles indicates that at these temperatures non-LTE effects are also important. The agreement between the LTE and non-LTE models at kK is merely coincidental: in the case of LTE the line center is deeper, whereas the line wings are weaker compared to the non-LTE profiles. These differences happen to cancel out approximately, yielding the same value for the equivalent width. In general we conclude that LTE calculations fail to predict the strengths and shapes of the investigated infrared hydrogen lines for stars with spectral type earlier than B3. The same conclusion is reached for the strongest infrared neutral helium lines.
The term "line blanketing'' signifies the effect that millions of spectral lines may have on the structure of a stellar atmosphere. In Paper I, we demonstrated that for Sco a model accounting for line blanketing in general yielded better infrared hydrogen profile fits. Typically, the profiles show a slightly enhanced core emission and a moderately enhanced wing absorption compared to H&He predictions. Essentially, these effects are a direct result of differences in the temperature structure, the temperature gradient in the line blanketed models being somewhat steeper. In Sco the differences in line shape are most pronounced for Br and Br. Compared to the non-blanketed case the equivalent width of Br is increased by about 1 Å (i.e. overall stronger line absorption). An effect of this magnitude may explain the systematic offset between observed and predicted strength of this line at kK, but should not occur at lower temperatures when the line fits the non-blanketed models used in Fig. 4 well. The overall strength of Br and Pf was not found to be much different between the blanketed and non-blanketed case.
Although one should not blindly extrapolate the results for Sco to the entire temperature regime investigated in this paper, the test case does seem to indicate that the effect of blanketing on the equivalent width is relatively modest. This also seems to be validated by the success of the H&He models in reproducing the observed infrared line strengths. To some extent this may be explained by the fact that the EW is calculated using bounds of kms-1. For this interval the "extra'' absorption in the wings of lines of models incorporating blanketing about cancels the "extra'' emission in the line center.
Unfortunately, it is not easy to extend our blanketing calculations such that they cover the temperature range of all late-O and early-B stars. The computational costs of such a grid are presently too large.
In the models the effects of turbulence have been neglected. Turbulence is expected to enhance non-LTE effects in the strongest IR hydrogen and helium lines. The presence of strong turbulence is likely the reason why the observations of Sco differ from our model predictions (see Paper I). To illustrate this: in Sco a broad emission is seen in Br reaching velocities of 150 kms-1.
Figure 6b shows models similar to those given in
Fig. 4, now using a
of 25 kms-1.
Table 4 provides a listing of equivalent width values
plotted in Fig. 6. For all temperatures we find that the
EW of Br
decrease (i.e. stronger line emission) by
respectively 2 and 3.5 Å. The introduction of
the overall fit worse for Br
whereas for Pf
15 kms-1 will remove the overall discrepancy between
observed and predicted equivalent width using H&He models. For the
lines which show line emission, i.e. those
with spectral type B0.5-B2, the line is expected to be purely
photospheric as possible wind contamination can likely be ruled out for
such spectral types (see Sect. 2.1). As the overall line formation of
occurs further out in the photosphere compared to Br,
increases outward in the atmosphere of dwarfs of
spectral type B0.5-B2. A detailed quantification of the turbulent
velocity gradient requires a detailed analysis of line contribution
functions, based on models with much higher spatial resolution than the
ones used in this study (which have 50 radial depth points). This is
beyond the scope of the present paper. In any case,
this interesting result demonstrates the
potential of H I infrared lines in studying the outer photospheric
structure of early-type stars.
|Figure 6: a) Comparison between observed equivalent widths for Br (solid), Br (dot-dashed) and Pf (dashed) and predicted values assuming LTE and = 0 kms-1, and b) assuming non-LTE and = 25 kms-1. For the LTE models we only used = 4.0. In the case of non-LTE models the same gravity values (and linestyles) are used as in Fig. 4|
|Open with DEXTER|
We will now address the possible presence of a stellar wind and its effect on the infrared profiles. The H&He model calculations assume a hydrostatic atmosphere. In Paper I we have demonstrated that in case of Sco a wind is not expected to affect the H I and He I line profiles. This conclusion is based on both an observational and a theoretical argument.
First, the observed line displacements for the different hydrogen lines do not seem to correlate with the rest wavelengths ( 's) of these lines in the sense that a larger displacement is expected for lines with larger . This is expected as due to the proportionality of free-free processes, the infrared continuum forming layers shift outward in the atmosphere. If a measurable velocity gradient would already be present in these layers a larger displacement should have been seen in, for instance, Pf compared to Br. Second, after linking the density structure to a velocity structure using the mass continuity equation, inspection of the Sco model showed that the cores of the lines discussed in this study all formed in highly subsonic layers.
We expect these arguments also to apply to the low mass-loss rates of the B-type dwarfs and (sub)giants investigated in this paper. (For B-supergiants the argument no longer holds as these stars suffer a much larger mass loss. See e.g. Najarro et al. 1998). As discussed in Sect. 5.2.1, the earliest spectral types in our sample show evidence for asymmetries in the observed hydrogen profiles which could not be attributed to He I blends. We can not fully exclude that these asymmetries are the result of a stellar wind, however, at the present moment it seems more reasonable to attribute them to macro turbulent velocities in the outer photosphere, as already suggested for Sco (Paper I).
Earlier work discussing infrared spectroscopy for spectral classification of OB stars was published by Hanson et al. (1996, 1998) and Blum et al. (1997). These groups focus on the classification of supergiants, although their samples also include some (sub)dwarfs.
Hanson et al. (1996) used Br, He I 2.058, 2.116, He II 2.190 and some metal line transitions in the K-band (2.04-2.2 m) for the classification of 180 optically visible, luminous OB stars. For the late O- and B-type dwarfs they could derive only a coarse classification as only Br and He I 2.116 could be used. In Sect. 3 we have demonstrated that Br is not ideal for spectral classification. As the spectral dependence for Br is even weaker (see Paper I), this explains why Hanson et al. (1996) are only able to provide a coarse spectral classification, with "coarse'' meaning the stars may be grouped in "early'', "mid'' and "late'' type. First attempts using the H-band for spectral classification (Blum et al. 1997; Hanson et al. 1998) yield similar results. Hanson et al. (1998) examined the spectral characteristics of 34 normal OB stars using the He I 1.700, He II 1.693 and the Br 11 1.7681 m lines in the H-band. For late O- and B-type dwarfs again only a coarse spectral classification could be derived. Compared to the K-band it is also more difficult to find a luminosity effect in these lines, although Hanson et al. (1998) appear to reveal a weak luminosity dependence using the ratio of Br 11 to He I 1.700 m.
In comparison with these spectral classifications using the K- and H-band, spectral classification using the Br and Pf line appears much more promising as these lines show a significantly stronger dependence on stellar temperature for the late O- and early B-type dwarfs allowing for a classification within one or two subclasses. The accuracy of the spectral classification is expected to improve when Hu becomes observable with the installation of VISIR in 2001 on the VLT. Contamination by excess flux at longer wavelengths ( m) due to warm circumstellar dust is not anticipated to be an important problem, as dust is not expected to be present in dwarfs or subgiants. However, dust contamination is expected for supergiants. Therefore, the use of Pf or Hu for spectral classification of supergiants may not be feasible, unless the dust contamination is weak or can be modeled.
With the observed line profiles we were not able to put constraints on
the value of .
Gravity causes only minor changes in the derived
equivalent widths for the observed H I infrared lines. As one also
has to deal with the uncertainty in microturbulent velocity, we expect
it to be difficult to determine
using the EW of H I and/or
He I lines. Alternatively one could derive the gravity by fitting
the Stark broadened wings of hydrogen lines like, Br
However, this requires high quality spectra with well determined
continua and also a more complex modeling accounting for the effects
of line blanketing.
|Figure 7: Comparison between the equivalent width of the observed He I 2.058 m line as given by Hanson et al. (1996) and those predicted in Paper I. The observed EWs are represented by dots. The large (medium, small) dots represent supergiants (giants, dwarfs). The model calculations are shown for four different values of : 4.3, 3.85, 3.4 and 3.0 (increasing line thickness). The dashed lines represent similar model calculations using a of 25 kms-1; these calculations have not been performed for the lowest gravity value|
|Open with DEXTER|
The K-band atlas of Hanson et al. (1996) includes the He I 2.058 line. These data allows for a comparison of observed and predicted equivalent widths. Figure 7 shows this comparison for the normal stars in their sample, i.e. excluding the emission line stars. In order to link the spectral type adopted by Hanson et al. (1996) to a temperature value we used the spectral class calibrations given by Vacca et al. (1996) for the O-type stars and those from Meyer et al. (1998) for the B-type stars. To allow for a comparison of the more evolved stars, we extended our gravity calculations to include = 3.4 and 3.0.
Figure 7 shows that for the O-type stars the predicted equivalent widths are in fair agreement with the observations. The supergiants (large dots) show the largest equivalent widths, whereas the main sequence stars (small dots) have an . The observed equivalent widths of the supergiants may be reproduced when including low gravity models. Note however that also an increase in the microturbulent velocity (the dashed lines are for = 25 kms-1) helps bring the models in range of these observations. The large scatter in the observed EWs for the O supergiants indicates that the correlation between He I 2.058 strength and luminosity is not pronounced. In the case of the B1-B2 supergiants the observed He I strength also shows a large scatter. For these later type supergiants the hydrostatic H&He models completely fail to predict the observed strong He I 2.058 line emission. The reason for this large discrepancy is most likely due to the presence of a strong stellar wind. However, one should be very prudent in applying this line. The resonance 1s2p-1s1s transition at 584 Å shares the same upper level with the He I 2.058 transition. The population of this upper level is dominated through the resonance channel. This implies that the He I 2.058 line is strongly sensitive to variations of the flux in the far-UV part of the spectrum (Morris et al. 1996). The line will thus respond to even minor density enhancements near the surface of the star. Such density variations are expected to be present in the onset of the stellar wind of early-type stars (Owocki et al. 1988). This may render He I 2.058 unsuited for spectral classification, especially for the more luminous supergiants. For the dwarfs the He I 2.058 line is weak. In Paper I, we observed the He I 2.058 profile of Sco which showed a weak red-shifted line emission. The line shape nor strength could be explained using a hydrostatic model. This serves to demonstrate the complexity of the He I 2.058 line formation, not only in supergiants but also in case of dwarfs.
We have demonstrated the potential of determining the stellar temperature from the equivalent width of Pf and Br. These lines are sensitive to temperature for dwarfs and giants within the spectral range O9-B2 (O9-B1 in case of Br). The Br sensitivity to is less strong compared to Pf. The quality of our observations allow one to derive a spectral type within one or two subclasses. Using high resolution and high S/N spectra obtained with UKIRT, the classification may be improved to within one subclass, using Br. The EW of Br depends only weakly on over the entire investigated spectral range. Therefore this line as well as higher Brackett series lines have only very limited potential for spectral classification.
Compared to the optical spectral classification, the effective temperatures estimated from Br and Pf lines have errorbars about a factor of two to four larger. As the dependence is weak for the observed infrared lines, we presently need to rely on the optical for a determination of the luminosity class. In the future one may be able to determine gravity from fitting Stark broadened lines like Br and Br. However, this requires high quality spectra with well determined continua and also a more realistic modeling, accounting for line blanketing.
The Br, Br and Pf line strengths as determined from models (Paper I) show a promising agreement with the programme stars observed with ISO. This indicates that the hydrostatic non-LTE H&He models are valid and allow for a meaningful prediction of temperature based on the H I and He I line strengths. This may seem surprising as these models are rather basic, i.e. they neglect line blanketing and/or effects of a wind. The most important ingredient in the model is to consider departures from LTE, as such departures strongly influence the equivalent width as well as the line shape.
Line blanketing seems only modestly important in determining the strengths of these infrared lines. This may in part be due to the adopted line bounds of kms-1 used to measure the equivalent width. These bounds are such that they result in a cancelation of additional wind absorption - compared to non-blanketed predictions - by additional core emission.
Microturbulence competes in importance with line blanketing when discussing line strengths. The inclusion of these small scale motions leads to an increased line emission. Essentially, turbulence tends to desaturate the lines and may enhance non-LTE effects (Sigut & Lester 1996). In order to obtain a better overall fit between observed and predicted Pf strength we needed to include turbulent velocities of magnitude kms-1. For Br an optimal fit was found for < 5 kms-1. This points to a gradient in the turbulent velocity in the photospheres of late-O and early-B dwarfs and subgiants, in the sense that increases with decreasing density.
In comparison with the lines observed in the H- and K-band spectroscopy (Hanson et al. 1996, 1998), Br and Pf show a larger sensitivity to effective temperature for late O- and early B-type dwarfs. The accuracy of the spectral classification is expected to improve when using the equivalent width of Hu. This line may be observed with VISIR when installed in 2001 on the VLT.
We thank Conny Aerts (KUL, Leuven, Belgium) for providing the effective temperature of the B-type stars. The data presented was analysed with the support of the Dutch ISO Data Analysis Centre (DIDAC) at the Space Research Organization Netherlands (SRON) in Groningen. This work was supported by NWO "Pionier'' grant 600-78-333 and by Spinoza grant 08-0 to E. P. J. van den Heuvel.