arxiv: v1 [astro-ph.sr] 9 Apr 2014

Transcription

1 Astronomy& Astrophysics manuscript no. 7 c ESO 8 June 7, 8 The virtual observatory service TheoSSA: Establishing a database of synthetic stellar flux standards II. NLTE spectral analysis of the OB-type subdwarf Feige, T. Rauch, A. Rudkowski, D. Kampka, K. Werner, J. W. Kruk, and S. Moehler arxiv:44.446v [astro-ph.sr] 9 Apr 4 Institute for Astronomy and Astrophysics, Kepler Center for Astro and Particle Physics, Eberhard Karls University, Sand, 776 Tübingen, Germany rauch@astro.uni-tuebingen.de NASA Goddard Space Flight Center, Greenbelt, MD 77, USA European Southern Observatory, Karl-Schwarzschild-Str., Garching, Germany Received 6 February 4/ Accepted 9 April 4 ABSTRACT Context. In the framework of the Virtual Observatory (VO), the German Astrophysical Virtual Observatory (GAVO) developed the registered service TheoSSA (Theoretical Stellar Spectra Access). It provides easy access to stellar spectral energy distributions (SEDs) and is intended to ingest SEDs calculated by any model-atmosphere code, generally for all effective temperature, surface gravities, and elemental compositions. We will establish a database of SEDs of flux standards that are easily accessible via TheoSSA s web interface. Aims. The OB-type subdwarf Feige is a standard star for flux calibration. State-of-the-art non-local thermodynamic equilibrium (NLTE) stellar-atmosphere models that consider opacities of species up to trans-iron elements will be used to provide a reliable synthetic spectrum to compare with observations. Methods. In case of Feige, we demonstrate that the model reproduces not only its overall continuum shape from the far-ultraviolet (FUV) to the optical wavelength range but also the numerous metal lines exhibited in its FUV spectrum. Results. We present a state-of-the-art spectral analysis of Feige. We determined T eff = 47 5± K, logg=6.±., and the abundances of He, N, P, S, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, and Ge. Ti, V, Mn, Co, Zn, and Ge were identified for the first time in this star. Upper abundance limits were derived for C, O, Si, Ca, and Sc. Conclusions. The TheoSSA database of theoretical SEDs of stellar flux standards guarantees that the flux calibration of astronomical data and cross-calibration between different instruments can be based on models and SEDs calculated with state-of-the-art modelatmosphere codes. Key words. Standards Stars: abundances Stars: atmospheres Stars: individual: Feige Stars: subdwarfs Virtual observatory tools. Introduction Feige is a bright (m V =.845±., Kharchenko & Roeser 9), subluminous OB-star (type sdob, Heber et al. 984a; type sdo D,Vennes et al. ). It is widely used as a spectrophotometric standard star (e.g. Oke 99; Turnshek et al. 99; Bohlin et al. 99). Since Feige will be used as a reference star for the flux calibration of X-SHOOTER (Vernet et al. ) observations from Å to 5 Å (Moehler et al. 4), we decided to reanalyze its spectrum with state-of-the-art model-atmosphere techniques. An early spectral analysis with approximate LTE, lineblanketed hydrogen model atmospheres yielded an effective tem- Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. instruments/xshooter.html local thermodynamic equilibrium perature T eff = 9 K and a surface gravity log(g/cm/s )= 6.5 (Newell 97). Kudritzki (976) showed that both, the consideration of deviations from LTE (NLTE ) as well as of opacities of elements heavier than H, have a significant influence on the determination of T eff and logg in an analysis of optical spectra (Table ). Heber et al. (984a) extended the analysis of Feige to the ultraviolet (UV) wavelength range (IUE 4 observations, 5 Å < λ < Å) in addition to highresolution optical spectra (4 Å< λ < 5 Å) and derived T eff = 4 +5 K, logg=5.±., and He/H=.+.. (by number) using H+He (with subsequent C+N+Si line-formation calculations) NLTE models. With the FUSE 5 mission, the interstellar deuterium and oxygen column densities toward Feige were measured. Friedman et al. () used optical spectra and estimated the atmospheric parameters by comparison with a grid of synthetic NLTE model-atmosphere spectra (using TLUSTY to compute non-local thermodynamic equilibrium 4 International Ultraviolet Explorer 5 Far Ultraviolet Spectroscopic Explorer Article number, page of 9

2 the stellar atmosphere model and SYNSPEC to generate the SED, Hubeny & Lanz 995, just TLUSTY here after), that considered H and He. They achieved T eff = 4 ± K, logg=5.95±.5, and He/H=.±.5. With the higher log g (in agreement with Kudritzki 976), their spectroscopic distance of d=88±4 pc agreed with the Hipparcos 6 parallax distance of d= pc. In the following, we describe our analysis in detail. In Sect., we give some remarks on the observations. Then, we introduce our models and the considered atomic data (Sect. ) and start with a preliminary analysis (Sect. 4) of the optical spectrum based on H+He models followed by a highly sophisticated analysis with metal-line blanketed models (Sect. 5). We summarize our results and conclude in Sect. 6.. Observations Our main optical spectrum is a median of 9 X-SHOOTER observations, taken between Oct 6 and July 5 with a 5 slit (the seeing was below during the observations) and an exposure time of s each. The achieved resolving power was R=λ/ λ 48. All spectra were extracted with ESO s standard pipeline-reduction software (with the actual version at the time of the respective observation). Heliocentric correction and correction to an airmass = were applied. In addition, we used optical HST/STIS 7 spectra (ObsIds O48 and O48 co-added) from the archive for the determination of the interstellar reddening. Our far-ultraviolet spectrum consists of two observations of Feige that were performed by FUSE, both in June and both through the LWRS spectrograph aperture. The dataset IDs were M88 and P44, with exposure times of 6. ks and.8 ks, respectively. Alignment of the four FUSE telescope channels was excellent throughout both observations, with RMS exposure-to-exposure variations in flux under.5 % in all channels. The processing of individual exposures to produce a combined spectrum spanning Å was the same as that described for G9 BB in Rauch et al. () and won t be repeated here. The signal to noise per. Å pixel in the continuum for the combined spectrum is typically 8: shortward of Å and : longward of Å. Approximately 7 % of the exposure time was obtained during orbital night. Comparison of the spectra obtained during day and night portions of the orbit found a discernible difference only in the cores of Ly β and Ly γ. Only the night data were used for these spectral regions. Weak airglow emission was still present at Ly β during orbital night, but the affected pixels had no impact on the analysis of the stellar spectrum. Additional ultraviolet spectra were retrieved from MAST. We used all available low-resolution IUE spectra (SWP77, SWP9, SWP888, SWP89, SWP89, SWP89, LWP9, LWP94, LWP95, LWP55, LWP56, LWP57, LWP58, and LWR785 co-added) and an HST/STIS spectrum (ObsId OBIE, exposure time 74. s, start time -- 8::54 UT, grating G4M, 9 Å < λ < 46 Å, aperture 5.5, resolution=. Å) Hubble Space Telescope/Space Telescope Imaging Spectrograph. Model atmospheres and atomic data For our model-atmosphere calculations, we use the Tübingen NLTE model-atmosphere package 8 (Werner et al. ; Rauch & Deetjen ), that assumes a plane-parallel geometry and considers opacities of elements from H to Ni (Rauch 997, ). The models are in hydrostatic and radiative equilibrium. TMAP was successfully used for many spectral analyses of hot, compact stars (e.g. Rauch et al. 7; Wassermann et al. ; Klepp & Rauch ; Ziegler et al. ; Rauch et al. ). The model atoms used in our model-atmosphere calculations were either retrieved from the Tübingen model-atom database 9 or compiled via the registered Virtual Observatory (VO) tool TIRO that uses Kurucz s atomic data and line lists (Kurucz 99, 9,, and priv. comm.). Table shows the statistics of our model atoms. 4. Preliminary analysis For a preliminary analysis, or verification of basic previous results, we employ the registered Virtual Observatory (VO) service TheoSSA and the related registered VO tool TMAW, to download pre-calculated synthetic spectral energy distributions (SEDs) or to calculate individual SEDs, respectively (cf. Rauch et al. ). Figure shows a comparison of SEDs with model parameters of Heber et al. (984a), Friedman et al. (), and of this work with the observed optical spectrum. The Balmer decrement is a sensitive indicator of logg (e.g. Rauch et al. 998), and we derive logg=5.9±.. At this logg, the Hei/Heii ionization equilibrium, i.e. the measured equivalent-width ratio of Hei and Heii lines, is well reproduced by our model at T eff = 46 5± K. The He line strengths are matched at a photospheric He abundance of 8± % by mass. Although the theoretical H and He line profiles agree well with the observation, the central depressions are not matched perfectly. This may be a hint of a weak Balmer-line problem (cf. Napiwotzki & Rauch 994; Rauch ) exists because metal opacities are neglected. For the same reason, our synthetic H+He SEDs are not suitable for an analysis of the Hi Lyman lines in FUV 4 spectrum. Fully metal-line blanketed model-atmospheres are mandatory for this purpose (Sect. 5). However, from our derived T eff, we are well in the parameter range where no deviation between Hi Lyman- and Balmer-line analysis is expected (Good et al. 4, their Fig. 4). We adopt our derived T eff and logg values, that also reproduce well the HST/STIS observation (Fig. ), for our further analysis and will verify them with our final, fully metal-line blanketed model. Within error limits, our preliminary values agree well with those of Friedman et al. (). Only T eff =± from their χ fit appears to be too optimistic. It is worthwhile to note, that the result of Kudritzki (976) (Table, his lower He abundance model) is relatively close to our result. 8 TMAP, 9 TMAD, Tübingen iron-opacity interface Theoretical Stellar Spectra Access, Tübingen Model-Atmosphere WWW Interface 4 far-ultraviolet Article number, page of 9

3 T. Rauch et al.: TheoSSA: Establishing a database of synthetic stellar flux standards. II. Feige relative flux He I.7.6 4/5. H:He = 89: (Heber) 4/5.95 H:He = 96: 4 (Friedman) 465/5.9 H:He = 9: 8 (our work) 4 45 Fig.. Comparison of three synthetic spectra with our optical observation of Feige. T eff, logg, and the H:He ratio by mass are indicated. f λ / erg/cm /s/a o Fig.. Comparison of the optical HST/STIS observation with our final model SED. The synthetic spectrum is convolved with a Gaussian (FWHM=5 Å) to match the resolution of the observation. The error bar indicates the visual brightness (m V =.847±.). 5. Line identification and detailed analysis Friedman et al. () identified photospheric lines from N iii - v, Svi, Criv -v, Feiii -iv, and Niiv in the FUSE observation. Their SED calculation (SYNSPEC, Hubeny & Lanz 995) included all elements from H to Zn, all with solar abundances but He (. times solar), C ( 4 6 times solar), Si ( 7 times solar), Cr ( times solar). We decided to include H, He, C, N, O, Si, P, S, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, and Ge in our calculations. Figure shows the ionization fractions of these elements. For the irongroup elements (here Ca - Ni), the dominant ionization stages areiv -v. All SEDs that were calculated for this analysis are available via TheoSSA. For the line identification in the FUSE wavelength range (95 88 Å), it was necessary to determine the stellar continuum flux precisely. We started with a measurement of the Fig. 4. FUSE observation of Feige (gray) compared with two synthetic spectra calculated from our final model (thin, blue in the online version: with Kurucz s POS lines; thick, red: with Kurucz s LIN lines). All spectra are convolved with a Gaussian (FWHM= Å) for clarity. interstellar neutral hydrogen column density. From Ly α in the STIS spectrum and the higher members of the Lyman series in the FUSE spectrum, we determined n Hi =.8±.8 cm in agreement with Friedman et al. (, n Hi =.4±.5 cm determined from the high-resolution IUE spectrum SWP57). To measure the interstellar reddening, we normalized our synthetic spectrum to the MASS H brightness because the interstellar reddening is negligible there and adjusted E B V to match the IUE, STIS, and FUSE flux levels. Our result is E B V =.7±.7. This very small value is in agreement with the absence of the 75 Å bump in the IUE LWP spectra (Sect. ). The Galactic reddening law of Liszt (4a,b, valid for.5 < E B V <.75 and b ), N Hi /E B V = 8. cm mag, predicts. E B V. in agreement with our value. The comparison of our models to the FUSE observations shows that we can reproduce well the observed flux level (Fig. 4), if we include all the lines from Kurucz s LIN lists (Sect. ). These include laboratory-measured lines with good wavelengths as well as theoretical lines 5. The lines with good wavelengths are presented in Kurucz s POS lists. Unfortunately, the ratio of LIN to POS lines is about and thus, most line wavelengths are uncertain. Moreover, the continuum flux of the POS-line spectrum appears artificially high compared to the LIN spectrum due to the neglected line opacity (Fig. 5). The line-identification process is easy (the comparison of two SEDs calculated from our final model where the oscillator strengths of one individual atom/ion was artificially reduced for one SED). It enabled us to unambiguously identify hundreds of lines of N, O, P, S, Ti, V, Cr, Mn, Fe, Mn, Ni, Zn, and Ge. Our metal-line blanketed models have a different atmospheric structure compared to the H-He models that were used in the preliminary determination of T eff = 46 5 K and logg=5.9 (Sect. 4). Figure 6 shows the typical surface-cooling (log m < -.5) and backwarming effects (log m >.5), that are an impact of the additionally considered metal opacities. A detailed evaluation of the optical spectrum shows that slightly higher T eff = 47 5 K and logg=6. values are necessary to reproduce the He i/ He ii ionization equilibrium and the observed H i, He i, and He ii line profiles best. The abundance analysis follows a standard procedure. Identified lines are reproduced by an abundance adjustment of the respective species. For elements with no lines identified, we in- 5 Kurucz s LIN lists are used in our model-atmosphere calculations. Article number, page of 9

4 Si IV P V Si IV Si IV.5. T / kk 5 He I Fe II is Fig. 5. Same as Fig. 4, for a section of the FUSE observation (top panel: POS, bottom panel: LIN). The POS lines are identified (in black, Fe and Ni lines in blue and green for clarity, respectively) at top of the top panel. Lines of interstellar origin are marked in blue (with subscript is ) at the bottom of the top panel. The synthetic spectra a convolved with a Gaussian (FWHM=.6 Å) to match FUSE s resolution. Fe II is log (m / g/cm ) Fig. 6. Temperature structure of our H+He model (thin, blue line: T eff = 46 5 K, logg=5.9), a metal-line blanketed model with the same T eff and logg (dashed, red), and of our final model (thick, red: T eff = 47 5 K, logg=6.). The formation depths of the lines cores of optical lines of H i (H α is the most outside), He i, and He ii are shown. [X] H He C N O Mg Si P S CaScTi V Cr MnFeCoNi Zn Ge AA Dor 4/5.46 Feige 465/5.9 EC 48-55/ atomic number creased the abundances until their line-detection limit. The optical spectrum was used to further constrain the upper limit because lines of lower ionization stages, that are not observed, appear there at too-high abundances in the synthetic spectrum. Table summarizes the lines that were used and the derived abundances. To identify ISM 6 absorption lines in the FUSE observation (Friedman et al., cf.) and to judge the contamination of photospheric lines, we follow our standard procedure and model the stellar spectrum simultaneously with the ISM line absorption (e.g., Rauch et al. ). We modeled the latter with the program OWENS (Hébrard et al. ; Hébrard & Moos ), that considers different clouds with individual radial and turbulent velocities, temperatures, column densities and chemical compositions. Lines are represented by Voigt profiles. The best fit is determined via aχ method. Our ISM model includes lines of H (J= 9), Hi, Di, Cii-iii, Ni-ii, Oi, Sii-ii, Pii, Ari, and Feii. Our results for Di and Oi are consistent with those of Friedman et al. (). The complete FUSE observation is compared (including line identifications) with our final model in Fig interstellar medium Fig. 8. Comparison of the determined photospheric abundances (arrows indicate upper limits) of the three OB-type subdwarfs AA Dor (Klepp & Rauch ), EC 48 Rauch et al. (), and Feige. Their T eff and logg are shown in the legend. 6. Results and conclusions We performed a comprehensive spectral analysis of Feige, based on observations from the FUV to the optical wavelength range. We determined T eff = 47 5± K and logg=6.±.. The ionization equilibria of Hei/Heii, Niii / Niv/Nv, Piv/Pv, Siv/Sv/Svi, Tiiv/Tiv, Viv/Vv, Criv / Crv/Crvi, Mnv/Mnvi, Fev/Fevi, Cov/Covi, and Niv / Ni vi are well reproduced with these values. The photospheric abundances were determined based on the FUSE and optical observations (Table ). Figure 8 shows a comparison of the photospheric abundances patterns of three hot O(B)-type subdwarfs. While the intermediate-mass metals are solar or subsolar in all these stars, the iron-group elements but Fe have strongly supersolar values. An exception is Fe in AA Dor and EC 48, that appears to be solar. Neither this Fe peculiarity nor the extremely low C and Si abundances in Feige can be explained. The position of Feige in the T eff logg plane shows that it is located directly on the He main sequence (Fig. 9). Feige belongs, like AA Dor or EC 48, to the hottest post- Article number, page 4 of 9

5 T. Rauch et al.: TheoSSA: Establishing a database of synthetic stellar flux standards. II. Feige Table. Strategic lines and determined element abundances (mass fraction, error±. dex). [X]=log (abundance/ solar abundance) of species X (solar values from Asplund et al. 9). element analyzed lines abundance [X] H optical Hi lines 9.. He optical Hei -ii lines C C IIIλλ Å <.4 7 < 4.6 C IVλλ 8.4,.49, Å N λλ 8.97, 8., 84.5, 84.57Å λλ 998.6, 4.58, 479., 45.9, , 464.4, Å λλ 9.99, 9.5, 9.6, 9., 9.68, 94.8Å N Vλλ 8.8, 4.8Å O O IIIλ 5.78 Å <.78 6 <.5 O IV λλ 9., 9.6, 9.7, 9.4Å Si Si IIIλ. Å <.96 7 <. Si IVλ.49 Å P P IV λλ 5.56, 8.9,.5,.5,., 5.5Å P Vλλ 7.98, 8.Å S S IV λλ 6.66, 7.97, 7.5, 98.9, 99.48Å λλ 9.9,., 8.67, 8.78Å Iλλ 9.8, 944.5Å Ca optical Caiv lines < <.6 Sc Sc IVλ 9.4 Å <.8 4 <.8 Sc Vλλ 99.4, 944.4Å Ti optical Ti iv lines, Ti IV λ 8.6 Å Ti V λλ 5.8, 6.5Å V V IVλ.5 Å V V λλ 978.6, 4.74, 57.58Å Cr many Criv -vi lines in the FUV, e.g. λλ 4.46, 65.6, 7., 96.64, 6.5Å λλ., 5.4, 4.55, 45.4, 6.65 Å Iλ 957. Å.9.6 Mn λλ 4.4, 4.65, 48.6, 49.4, 55.98, 6.49, 7.6 Å.9.5 λλ 8.9,.58Å Fe many Fev -vi lines in the FUV, e.g..8.8 λλ.87, 5.,.6Å Iλλ.9, 67.7Å Co many Co v - vi lines in the FUV, e.g λλ 79.59, 8.9, 84.6Å λλ.7, 4.77, 5., 69.55, 75.6 Å Ni many Niv -vi lines in the FUV, e.g..8.5 λλ 4., 78.9Å λλ.9, 57.55, 59., 78.7Å Zn Zn Vλλ 6.84,., 58.76Å Ge Ge Vλλ 6.67, 69., 7.66, 6.95, 65.6Å EHB 7 stars. From a comparison to post-ehb tracks (Dorman et al. 99), we can extrapolate a stellar mass of M=.469±. M. With R= GM/g (G is the gravitational constant), we calculated the stellar radius of R= R. We determined the distance of Feige following the flux calibration of Heber et al. (984b) forλ eff = 5454 Å, d spec = 7. 4 H ν M.4 m V logg pc, () with m Vo = m V.75c, c=.475e B V, and the Eddington flux H ν = 7.4±.7 4 erg/cm /s/hz atλ eff of our final model atmosphere. We used E B V =.7±.7 (Sect. 4), 7 extended horizontal branch M=.469±+. M, and m V =.847±. (Kharchenko & Roeser 9) and derived a distance of d spec = pc and a height below the Galactic plane of z= pc. This distance is about a factor of three larger than the new Hipparcos parallaxmeasurement reduction (van Leeuwen 7, HIP595, π = 9.76±.44 mas) of d parallax = pc. Interestingly, the older Hipparcos measurement published by Perryman et al. (997,π=5.59±.4 mas, d parallax = pc) deviates from this new value by a factor of almost two and would be in agreement with our spectroscopic distance within error limits. The discrepancy between spectroscopic and parallax distances is a significant problem. log g cannot be higher by about.5 to achieve a distance agreement, because the spectral lines in the models appear too broad and too shallow. This apparently Article number, page 5 of 9

6 log (g / cm/s ) EC AA Dor Feige He main sequence 5 4 T eff / kk terminal-age EHB zero-age EHB The established database of spectrophotometric standard stars in TheoSSA was extended by the OB-type subdwarf Feige. The successfully launched GAIA 8 mission will provide accurate parallax measurements for spectrophotometric standard stars. This will strengthen the importance of a VOcompliant database like TheoSSA that provides easy access to the best synthetic spectra calculated for these stars. Acknowledgements. TR is supported by the German Aerospace Center (DLR, grant 5 OR ). The GAVO project at Tübingen has been supported by the Federal Ministry of Education and Research (BMBF, grants 5 AC 6 VTB, 5 AC VTB). This work used the profile-fitting procedure OWENS developed by M. Lemoine and the FUSE French Team. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. This research has made use of NASA s Astrophysics Data System. Some of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS Support for MAST for non-hst data is provided by the NASA Office of Space Science via grant NNX9AF8G and by other grants and contracts. The TIRO service ( astro-uni-tuebingen.de/~tiro) used to calculate opacities for this paper was constructed as part of the activities of the German Astrophysical Virtual Observatory. The TMAW service ( used to calculate theoretical spectra for this paper was constructed as part of the activities of the German Astrophysical Virtual Observatory. Fig. 9. Location of Feige in the T eff logg plane compared to sdbs and sdobs from Edelmann (), Rauch et al. (, EC 48 ), and Klepp & Rauch (, AA Dor). Post-EHB tracks from Dorman et al. (99, Y HB =.88, labeled with the respective stellar masses in M ) are also shown. Their start and kink points are used to illustrate the location of the zero-age and terminal age EHB for this He composition. The He main sequence is taken from Paczyński (97). is not a problem of our TMAP code, because Friedman et al. (, d spec = 88±4 pc) used the TLUSTY code and encountered the same problem. Similar discrepancies are reported by Rauch et al. (7, for LSV+46 with TMAP: d spec = pc vs. d parallax= pc) and by Latour et al. (, for BD+8 4 with TLUSTY, T eff = 8 K, logg=6., and an assumed M=.5 M : d spec = 57 pc (no error estimate given) vs. d parallax = 9 + pc). Latour et al. () mentioned that a relatively high logg value and/or a low mass may be the solution and since they regard the HIPPARCOS measurement as fully reliable and their TLUSTY results reasonably reliable, the mass of BD+8 4 must be much less than the canonical post-ehb mass of about.5 M. For their d parallax /d spec =.59, the mass has to be about.7 M. In case of Feige, with d parallax /d spec =., the mass has to be about. M. In both cases, the mass can be higher, if logg is higher. Thus, since logg is also the main error source in the spectroscopic distance (Eq. ), one might speculate about the applied broadening theory for lines that are used to determine logg. For the relevant Hi and Heii lines (linear Stark effect), TMAP as well as TLUSTY use the same data of Tremblay & Bergeron (9) and Schöning & Butler (989), respectively. However, all the narrow metal lines (e.g. of the iron-group element) in the UV, that are broadened by the quadratic Stark effect, cannot be reproduced at a much higher log g. To summarize, the distance discrepancy is as yet unexplained. The analysis of the FUV spectrum has shown that the lack of reliably measured wavelengths of lines of the iron-group elements (Ca - Ni) and of elements heavier than Ni hampers the line-identification. highly desirable. Efforts in this field in the near future are References Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 9, ARA&A, 47, 48 Bohlin, R. C., Harris, A. W., Holm, A. V., & Gry, C. 99, ApJS, 7, 4 Dorman, B., Rood, R. T., & O Connell, R. W. 99, ApJ, 49, 596 Edelmann, H., Dissertation, University Erlangen-Nuremberg, Germany Friedman, S. D., Howk, J. C., Chayer, P., et al., ApJS, 4, 7 Good, S. A., Barstow, M. A., Holberg, J. B., et al. 4, MNRAS, 55, Heber, U., Hamann, W.-R., Hunger, K., et al. 984a, A&A, 6, Heber, U., Hunger, K., Jonas, G., & Kudritzki, R. P. 984b, A&A,, 9 Hébrard, G., Friedman, S. D., Kruk, J. W., et al., Planet. Space Sci., 5, 69 Hébrard, G. & Moos, H. W., ApJ, 599, 97 Hubeny, I. & Lanz, T. 995, ApJ, 49, 875 Kharchenko, N. V. & Roeser, S. 9, VizieR Online Data Catalog, 8, Klepp, S. & Rauch, T., A&A, 5, L7 Kudritzki, R.-P. 976, A&A, 5, Kurucz, R. L. 99, in NATO ASIC Proc. 4: Stellar Atmospheres - Beyond Classical Models, ed. L. Crivellari, I. Hubeny, & D. G. Hummer, 44 Kurucz, R. L. 9, in American Institute of Physics Conference Series, Vol. 7, American Institute of Physics Conference Series, ed. I. Hubeny, J. M. Stone, K. MacGregor, & K. Werner, 4 Kurucz, R. L., Canadian Journal of Physics, 89, 47 Latour, M., Fontaine, G., Chayer, P., & Brassard, P., ApJ, 77, 84 Liszt, H. 4a, ApJ, 78, 7 Liszt, H. 4b, ApJ, 78, Moehler, S., Modigliani, A., Freudling, W., et al. 4, A&A, in press Napiwotzki, R. & Rauch, T. 994, A&A, 85, 6 Newell, E. B. 97, ApJS, 6, 7 Oke, J. B. 99, AJ, 99, 6 Paczyński, B. 97, Acta Astron.,, Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 997, A&A,, L49 Rauch, T. 997, A&A,, 7 Rauch, T., A&A, 56, 665 Rauch, T., A&A, 4, 79 Rauch, T. & Deetjen, J. L., in Astronomical Society of the Pacific Conference Series, Vol. 88, Stellar Atmosphere Modeling, ed. I. Hubeny, D. Mihalas, & K. Werner, Rauch, T., Dreizler, S., & Wolff, B. 998, A&A, 8, 65 Rauch, T., Werner, K., Bohlin, R., & Kruk, J. W., A&A, 56, A6 Rauch, T., Werner, K., & Kruk, J. W., Ap&SS, 9, Rauch, T., Ziegler, M., Werner, K., et al. 7, A&A, 47, 7 Schöning, T. & Butler, K. 989, A&AS, 78, 5 Tremblay, P.-E. & Bergeron, P. 9, ApJ, 696, 755 Turnshek, D. A., Bohlin, R. C., Williamson, II, R. L., et al. 99, AJ, 99, 4 van Leeuwen, F. 7, A&A, 474, 65 Vennes, S., Kawka, A., & Németh, P., MNRAS, 4, 95 Vernet, J., Dekker, H., D Odorico, S., et al., A&A, 56, A5 Wassermann, D., Werner, K., Rauch, T., & Kruk, J. W., A&A, 54, A9 Werner, K., Deetjen, J. L., Dreizler, S., et al., in Astronomical Society of the Pacific Conference Series, Vol. 88, Stellar Atmosphere Modeling, ed. I. Hubeny, D. Mihalas, & K. Werner, Ziegler, M., Rauch, T., Werner, K., Köppen, J., & Kruk, J. W., A&A, 548, A9 8 Article number, page 6 of 9

7 Table. T eff and logg of Feige determined by Kudritzki (976). He/H gives his models abundance ratio by number. LTE NLTE T eff / K log (g/cm/s ) T eff / K log (g/cm/s ) A&A 7, Online Material p 7 He/H

8 A&A 7, Online Material p 8 Table. Statistics the atoms used in our calculations. In the case of iron-group elements (Ca Ni), the super lines include the sample lines (Kurucz s LIN lines, cf. Rauch & Deetjen ). levels levels ion NLTE LTE lines ion NLTE LTE super lines sample lines 5 5 Ca II II III He I IV 7 9 II 9 V III VI C II 6 7 VII III 54 Sc II IV III V IV N II 5 8 V III 4 9 VI IV 6 78 VII V Ti II VI III O II 6 6 IV 7 7 III 54 8 V IV VI V VII VI V II Si III III IV IV V 5 59 V 7 6 VI VI P IV VII V 8 7 Cr II VI III S III 5 IV IV 7 8 V V VI VI 5 5 VII VII Mn II Zn III III IV 75 IV V V VI VI Ge III 5 VII IV 8 8 Fe II V III VI IV V VI VII Co II III IV V VI VII Ni II III IV V VI VII total

9 A&A 7, Online Material p 9 5 log ionization fractions log (n e / cm - ) T / kk T n e H He C Si P S Ca Sc Cr Mn Fe Co Ni log ionization fractions - N Ti Zn - - O V Ge - log (m / g/cm ) I II III IV V VI VII Fig.. Temperature and density structure and ionizations fractions of our final model with T eff = 47 5 K and logg=6.. m is the column mass, measured from the outer boundary of our model atmosphere.

10 N II is N II P IV Si V O IV O IV O IV O IV I Fig. 7. FUSE spectrum of Feige (gray) compared with synthetic spectra calculated from our final model (red: Kurucz s LIN lines and interstellar line absorption, blue: Kurucz s POS lines). The locations of the strongest stellar (marked in black, Fe and Ni which have the most lines are marked in blue and green for clarity, respectively, for Ca - Ni, only Kurucz s POS lines are identified) and interstellar lines (marked in blue, with subscript is ) are indicated at the top of the panels. The small, green identification marks at the bottom of the panels indicate the locations of interstellar H lines. A&A 7, Online Material p

11 A&A 7, Online Material p I P IV N I is N I is N I is I I P II is V V N I is N I is " Fig. 7. Continued.

12 Fig. 7. Continued O IV Fe II is O V Si I is O V O V O V P V P V I I C III is C III V V O III O III O III O III O IV A&A 7, Online Material p

13 A&A 7, Online Material p I I I Ge V O V O V O V Si II is P IV O I is P IV Fig. 7. Continued.

14 A&A 7, Online Material p 4 P IV P IV O VI P IV O IV O III O III O III O IV P IV C II is C II is O VI Co IV Ar I is Fig. 7. Continued.

15 A&A 7, Online Material p 5 S IV Si IV Si IV Si IV Ge V Ge V S IV S IV N II is I Fig. 7. Continued.

16 Fig. 7. Continued Ti V V IV Si III Si III V IV V IV Si III Si III Si III Ti V Zn V Ge V I P V C IV P IV S IV S IV I A&A 7, Online Material p 6

17 Fig. 7. Continued Ge V V V V V Fe II is Zn V Fe II is Si IV Ge V Fe II is P V Si IV Si IV V IV Ti V Ni IV N II is N II is A&A 7, Online Material p 7

18 Fig. 7. Continued Ge V I I I Ti V C IV C IV C IV I Fe IV V IV Ti V Sc V I Ti V Ti V V V Zn V V V I Ge V Ti V A&A 7, Online Material p 8

19 A&A 7, Online Material p 9 Ti IV I Fig. 7. Continued.