arxiv: v1 [astro-ph.sr] 26 Jul 2016

Transcription

1 Mon. Not. R. Astron. Soc. 000,?? (20) Printed 9 November 20 (MN LaT E X style file v2.2) NLTE carbon abundance determination in selected A- and B-type stars and the interpretation of C i emission lines arxiv:0.02v [astro-ph.sr] 2 Jul 20 S. A. Alexeeva, T. A. Ryabchikova, & L. I. Mashonkina Institute of Astronomy of the Russian Academy of Sciences 8 Pyatnitskaya St. 90, Moscow, Russia Released 20 Xxxxx XX INTRODUCTION ABSTRACT We constructed a comprehensive model atom for Ci Cii using the most up-todate atomic data available and evaluated the non-local thermodynamic equilibrium (NLTE) line formation for Ci and Cii in classical D models representing the atmospheres of A and late B-type stars. Our NLTE calculations predict the emission that appears at effective temperature of 920 to 000 K depending on log g in the Ci 8, 90Å singlet lines and at T eff > 000 K (log g = ) in the Ci 90 9Å, 90 98Å triplet lines. A prerequisite of the emission phenomenon is the overionization-recombination mechanism resulting in a depopulation of the lower levels of Ci to a greater extent than the upper levels. Extra depopulation of the lower levels of the transitions corresponding to the near-infrared lines, is caused by photon loss in the UV lines Ci 29, 90, and Å. We analysed the lines of Ci and Cii in Vega, HD, Sirius, 2 Peg, π Cet, HD 22, and ι Her taking advantage of their observed high-resolution spectra. The C i emission lines were detected in the four hottest stars, and they were well reproduced in our NLTE calculations. For each star, the mean NLTE abundances from lines of the two ionization stages, Ci and Cii, including the Ci emission lines, were found to be consistent. We show that the predicted C i emission phenomenon depends strongly on whether accurate or approximate electron-impact excitation rates are applied. Key words: line: formation, stars: abundances Rapid development of observational techniques over the last decade has resulted in dramatic improvement of the quality of spectral observations in astronomy. Thanks to Echelle spectrographs, we are able to cover wide spectral regions with one exposure, obtaining spectra with a high spectral resolution of up to R = λ/ λ = and high signal-tonoise ratio. Huge amounts of high quality spectra are collected in different archives with the open access for the astronomical community. High quality spectra require immediately an adequate improvement in theoretical methods of spectrum analysis, model atmospheres, line formation scenarios, etc. Fossati et al. (2009) studied three apparently normal A and B sharp-lined stars and found the Ci emission lines at 8, and 90Å in the hottest of them, π Cet. The absorption Ci lines in the 20Å range are rather shallow and allow us to determine only an upper limit for the abundance. In local thermodynamic equilibrium(lte) analysis, π Cet reveals a disparity between Ci and Cii, in line with Roby & Lambert (990). For the cooler star, 2 Peg, Fossati et al. (2009) could not obtain consistent abundances from different Ci lines, namely the abundances obtained from Ci 92, 02, 80 Å were found to be significantly smaller than those from other Ci lines. Fossati et al. (2009) discussed several mechanisms of the Ci emission in π Cet. Emission lines may form, if the star has a chromosphere, however, there is no evidence for chromospheric activity in π Cet. This star was classified as a Herbig Ae/Be star due to small infrared (IR) excess (Malfait et al. 998), suggesting the existence of a circumstellar disc. This is supported by detection of an emission signature in H α (Fossati et al. 2009). However, half-widths of the Ci emission lines do not differ from those of the Ci absorption lines, giving evidence for their common origin in the star s atmosphere. Fossati et al. (2009) proposed that the emission can be caused by the departures from LTE in the Ci line formation. In the literature there are examples of successfully reproducing the observed emission lines by non-local thermodynamic equilibrium (NLTE) calculations,

2 2 S. A. Alexeeva, T. A. Ryabchikova and L. I. Mashonkina for example, Mgi 2 µm in the Sun (Carlsson et al. 992), Mnii22-2Åinthethreelate typebstars (Sigut200), Cii, 2Å in τ Sco (B0V) and Cii 2Å in HR 8 (BV) (Nieva & Przybilla 2008). There are few NLTE studies of Ci/Cii in the early A to late B stars. Rentzsch-Holm(99) computed a negative and small absolute value of less than 0.0 dex NLTE abundance corrections for lines of Cii in A-type stars. Przybilla et al. (200) obtained very small NLTE corrections for the Ci and Cii lines in the visible region using a Vega model atmosphere with an effective temperature T eff = 90 K. The NLTE carbon abundance analysis of 20 sharp-lined early B stars in the effective temperature range of K, including two stars with T eff < 8000 K, was performed by Nieva & Przybilla (202). None of the carbon NLTE papers investigated the NLTE effects for the Ci near-ir lines in the early A to late B stars. This paper aims to understand mechanism(s) of the Ci emission in B-type stars and to treat the method of accurate abundance determination from different lines of C i and Cii based on NLTE line formation. We construct a comprehensive model atom for Ci Cii using the most up-to-date atomic data available so far and analyse lines of Ci and Cii in high-resolution spectra of reference A and B-type stars with well-determined stellar parameters. The paper is organized as follows. Section 2 describes an updated model atom of C i-c iiand discusses departures from LTE for Ci Cii in the model atmospheres of A- B stars and mechanisms driving the Ci emission lines. In Sect., we analyse the Ci near-ir emission lines observed in the four B-type stars and determine the C abundance of the selected A-B stars. We inspect the abundance differences between different lines of a common species and between two ionisation stages, Ci and Cii. Our conclusions are summarized in Sect.. 2 NLTE LINE FORMATION FOR Ci-Cii 2. Model atom and atomic data Alexeeva & Mashonkina (20) developed and described in detail the carbon model atom that represents the term structure of Ci well and only approximately that of Cii. Here, the model atom was updated by extending the Cii term structure and including recent accurate collisional data for Ci. Energy levels of Cii. Energy levels of Cii belong to the doublet terms of the 2s 2 nl (n = 2 0, l = 0 2), 2s 2 nf (n = 0), 2s 2 ng (n = 0), 2s 2 2p nh (n = 0), 2s2p, and 2s2p l (l = 0 ) electronic configurations and the quartet terms of 2s2p 2, 2s2p, 2p, 2s2p l (l = 0 ). We thus include explicitly the Cii energy levels up to 0.0 ev below the ionization threshold. Fine structure splittingwas takenintoaccount everywhere, upton=. All the Ciistates with n = to haveclose together energies, therefore they were combined into superlevels. Energy levels were adopted from the NIST database (Ralchenko et al. 2008). The term diagram for Cii is shown in Fig.. Radiative data. We included 28 allowed boundbound (b b) transitions, whose transition probabilities were taken from the NIST and VALD (Kupka et al. 999) databases, where available, and the Opacity Project (OP) database TOPbase 2 (Cunto et al. 99, Luo et al. 989, Hibbert et al. 99). Their Photo-ionization cross-sections for levels of Cii with n 0, l are from TOPbase, and we adopted the hydrogen-like cross-sections for the higher excitation levels. Collisional data. For all the transitions between levels oftheci2p 2, 2p,2pl (l =0,, 2), and 2ps electronic configurations we use the effective collision strengths computed by Wang et al. (20) on the base of the B-spline R-matrix with pseudo states method. For the transitions involving the higher excitation levels of the 2pp, 2pd, and 2ps electronic configurations we adopted cross-sections calculated by Reid (99) in the close-coupling approximation using the R-matrix method. For the transitions connecting the 0 lowest levels in C ii we use the effective collision strengths from Wilson et al. (200). The remaining transitions in both Ci and Cii were treated using an approximation formula of van Regemorter (92) for the allowed transitions and assuming the effective collision strength Ω ij = for the forbidden transitions. Ionization by electronic collisions was everywhere treated through the Seaton (92) classical path approximation. 2.2 Method of calculations We used the code DETAIL (Butler & Giddings 98) based on the accelerated Λ-iteration method (Rybicki & Hummer 99) for solving the radiative transfer and statistical equilibrium (SE) equations. The opacity package was updated by Przybilla et al. (20). The departure coefficients, b i = n NLTE / n LTE, were then used to calculate the synthetic NLTE line profiles via the binmag code (Kochukhov 200) and synthv-nlte code (Ryabchikova et al. 20). Here, n NLTE and n LTE are the SE and thermal (Saha Boltzmann) number densities, respectively. We used plane-parallel (D), chemically homogeneous model atmospheres calculated with the LLmodels code (Shulyak et al. 200). An exception is Sirius, for which we took the Kurucz s model computed with the atmospheric parameters close to those derived by Hill & Landstreet (99). The carbon abundance was determined from line profile fitting. The uncertainty in fitting the observed profile is less than 0.02 dex for weak lines and 0.0 dex for strong lines. 2. SE of carbon In Fig.2, we show the LTE and NLTE fractions of Ci, Cii, and C iii in the model atmospheres of different temperatures. In each model, C ii dominates the total element abundance throughout the line formation region. The fraction of Ci everywhere does not exceed several thousandth parts, and in the line formation layers, it is smaller in NLTE than in LTE due to ultraviolet (UV) overionization. 2 http//legacy.gsfc.nasa.gov/topbase

3 The interpretation of C i emission lines p * p * s * 2p 2p 98 2 p * 2p p * s * p * d * d * Energy (cm - ) p 2 8 2p 2 2p 2-2 2p 2 0 ns 2 S 2 P np 2 P o 2 nd 2 D 2 D o nf 2 F o ng 2 G nh 2 H o S S o P P o D F o D o Figure. The term diagram for singly-ionized carbon. The dashed lines indicate the seven transitions, where the investigated spectral lines arise. Fig. displays a behaviour of the departure coefficients for selected levels in the model atmospheres with T eff / log g / [Fe/H] = 000 /. / 0 and 00 /.8 / 0. Deep in the atmosphere, the departure coefficients are equal to unity because of a dominance of the collision processes in establishing the statistical equilibrium. In the line formation layers, outward log τ 000 =0., NLTE leads todepleted populations of the Ci levels because of UV overionization. The deviations from LTE are small everywhere for the ground state of Cii. In the model 000/./0, outwards log τ 000 = 0., the Cii high-excitation levels p 2 P, s 2 S, andp 2 P )are overpopulated viatheuvradiative pumping transitions 2p 2 2 D p 2 P (Cii 0Å), p 2 P s 2 S (Cii 98Å), and 2p 2 2 D p 2 P (Cii 2Å). In the model 00/.8/0, these UV lines form in higher atmospheric layers, around log τ 000 =.0, resulting in a depopulation of the upper levels via spontaneous transitions. The NLTE effects for a given spectral line can be understood from analysis of the departure coefficients of the lower and upper levels, b l and b u, at the line formation depths. The line can be strengthened or weakened compared with its LTE strength depending on the b l and b u values. Here, we consider two lines with similar excitation energy of the lower level, E exc, but different gf values. Atomic data for the investigated lines are indicated in Table. In the model 000/./0, the Ci 02Å line (s P p D) is weak and it forms in the layers around logτ 000 = 0, where b l < and b l < b u. NLTE leads to weakening this line and positive NLTE abundance correction of NLTE = logǫ NLTE logǫ LTE = +0. dex. In contrast, Ci 9Å (s P p P) is strong and its core forms around logτ 000 =.2, where b l > b u and the line source function drops below the Planck function resulting in strengthened line and negative NLTE abundance correction of NLTE = Mechanisms driving the C i emission Our calculations predict that the Ci near-ir lines at 8Å (s P p S),90Å(s P p D),90 9Å(s P p P), and 90 98Å (s P p S) may appear as emission lines depending on the atmospheric parameters. The Ci 8 and 90Å singlet lines reveal an emission at the lower effective temperature compared with the triplet lines. The lower log g, the lower T eff is, at which the emission profile is developed for given line. Fig. shows how a variation in T eff and log g affects theoretical profiles of Ci 90, 98, and 9088 Å and at which atmospheric parameters the line absorption changes to emission. For example, the Ci 90Å emission line appears at 920 K < T eff < 900 K in the log g = 2 models and at T eff > 000 K when log g =. In contrast, the triplet line at 98Å comes into emission at T eff > 000 K in the log g = models. Having appeared at specific T eff /log g, the emission is strengthened towards higher temperature, reaches its maximum and gradually disappears. For example, in the models with log g =, the strongest emission in Ci 90Å is predicted for T eff = 000 K, and it almost disappears at T eff = K (Fig. ). Which processes in Ci do drive the emission in the IR lines? How do they depend on atmospheric parameters? In Fig., we present the NET rates, NET = n l (R lu + C lu ) n u(r ul +C ul ), in the model 000/./0 at logτ 000 =., where the Ci 8 and 90Å emission form (Fig.). Here, R lu and C lu are radiative and collisional rates, respectively,

4 S. A. Alexeeva, T. A. Ryabchikova and L. I. Mashonkina log [n(ci / II / III) / N(C)] /.9 / -0. LTE CI LTE CII LTE CIII NLTE CI NLTE CII NLTE CIII τ 000 log [n(ci / II / III) / N(C)] /. / τ /.8 / /.8 / 0.0 log [n(ci / II / III) / N(C)] log [n(ci / II / III) / N(C)] τ τ 000 Figure 2. NLTE and LTE fractions of, I, and II in the model atmospheres of different effective temperature. for the transition from low to upper level, and R ul and C ul for the inverse transition. The high-excitation levels, with E exc >.9 ev, including the upper levels of the investigated transitions, are mainly populated by recombination from C ii that follows an overionization of the Ci ground state (NET = s cm ) and low-excitation levels, including the lower levels of the investigated transitions (NET varies between 0 0 s cm and 0 s cm ). Thus, the overionization-recombination mechanism resulting in a depopulation of the lower levels to a greater extent than the upper levels is a prerequisite of the emission phenomenon. Extra depopulation of the levels s P and s P can be caused by photon loss in the UV transitions 2p 2 S s P (29Å), 2p 2 D s P (90Å), and 2p 2 P s P (Å). An influence of these UV transitions on the SE of Ci in the layers, where the IR lines form, depends on atmospheric parameters. In the model 000/./0, photon loss in Ci 29Å drains population of s P effectively in the layers around logτ 000 =., with NET = 0 0 s cm, while the transitions corresponding to Ci 90Å and Å are in detailed balance because of τ 90 and τ. We performed test calculations by reducing a radiative rate of Ci 29Å by 0%. This resulted in weakened emission, such that the line centre relative flux of Ci 90Å reduced from.0 to.028. With increasing T eff, optical depth of formation of the K 000 K 8000 K 9000 K K K Figure. Changes in the Ci 90Å line profile with temperature increasing. Everywhere, log g =, [C/Fe] = 0, Vsini = 0 kms, and ξ t = kms. The theoretical spectra are convolved with an instrumental profile of R = 000. Ci UV lines shift inwards (Fig. for T eff = 00 K) and photon losses depopulate not only s P, but also s P, resulting in the emission triplet lines at 90 9Å and 90 98Å.

5 The interpretation of C i emission lines Figure. Departure coefficients for selected levels of Ci (top row) and Cii plus the Ciii ground state (bottom row) as a function of logτ 000 in the model atmospheres 000 /. / 0.0 (left-hand column) and 00 /.8 / 0.0 (right-hand column). The wavelengths and vertical marks indicate the location of optical depth unity for the lines under consideration. CARBON LINES IN THE SELECTED STARS. Stellar sample, observations, and atmospheric parameters Our sample includes seven sharp-lined stars, with welldetermined stellar parameters (Table 2). Due to low rotational velocities, with Vsini 0 kms, theyallow spectral analysis to be done at the highest precision, maximizing the chance to identify even weak emission lines (WELs) in their spectra. They are either single stars, or primaries in close binary systems (SB) with much fainter companions, or individual components in a visual binary. For all thestars, thespectrain thevisible spectral range ( Å) were obtained with the Echelle SpectroPolarimetric Device for the Observation of Stars (ESPaDOnS) attached at the. m telescope of the Canada France Hawaii Telescope (CFHT) observatory. Spectra were extracted from the ESPaDOnS archive. The resolving power and signal-to-noise ratio are R = 000 and S/N 00-00, respectively. For Sirius, we also used the UV spectrum that covers the wavelength region between 2 and 8 Å. It was obtained with the Goddard High Resolution Spectrograph on the HubbleSpaceTelescope(HST), and it has R = and S/N of , in general. Characteristics of the observed spectra for individual stars are given in Table. We comment on the programme stars. The star 2 Peg (HD 2099) is known from previous studies as a normal single star with Vsini kms (Sadakane 98, Fossati et al. 2009). The star π Cet (HD 08) is an SB star with Vsini 20 kms. The observed small infrared excess at λ µ makes it a candidate for Herbig Ae/Be star (Malfait et al. 998). Since its spectrum is not visibly contaminated by the companion, it fully serves for our analysis. This star was already used by Smith & Dworetsky (99) as a normal comparison star in the abundance study of chemically peculiar stars. For 2 Peg and π Cet, their effective temperatures and surface gravities were spectroscopically determined from the Balmer lines (Fossati et al. 2009). HD 22 is also a normal star with Vsini kms, showing no chemical peculiarities. Its fundamental parameters were derived by Bailey & Landstreet (20) based on the Geneva and uvbyβ photometry. Effective temperature and surface gravity of Vega

6 S. A. Alexeeva, T. A. Ryabchikova and L. I. Mashonkina Table. Lines of Ci and Cii used in abundance analysis. λ, Å Transition log gf E exc, ev λ, Å Transition log gf E exc, ev Ci Ci p 2 P 2p P s P 2 p P p 2 P 2p P s P 2 p P p 2 P 2p P s P 0 p S p 2 P 2p P s P 2 p S p 2 P 2p P p P d P p 2 D 2 d P s P 2 p D p 2 D 2 d F s P p D p 2 P s P s P 0 p D p 2 P 2 s P s P 2 p D s P p P s P p D s P p P s P 2 p D s P p P Cii. s P 2 p P p 2 P /2 2p 2 2 D / s P 2 p P p 2 2 D /2 2p 2 D / s P p S p 2 2 D /2 2p 2 D / s P p D p 2 2 D /2 2p 2 D / s P p P p 2 2 D /2 2p 2 D / p D s P p 2 P /2 s 2 S / p D d F p 2 P /2 s 2 S / p D d F d 2 D /2 f 2 F / p D 2 s P d 2 D /2 f 2 F / p P d P d 2 D /2 f 2 F / p D d F ps P /2 2pp P / p D d F ps P /2 2pp P / p D 2 d F ps P /2 2pp P / p D s P ps P /2 2pp P / p D s P ps P /2 2pp P / p D 2 s P ps P /2 2pp P / s P p S ps P /2 2pp P / s P p D s 2 S /2 p 2 P / s P p P s 2 S /2 p 2 P / s P 0 p P p 2 P /2 - d 2 D / s P p P p 2 P /2 - d 2 D / s P p P p 2 P /2 - d 2 D / Table 2. Atmospheric parameters of the selected stars and sources of the data. HD Name Sp. t. T eff log g [Fe/H] ξ t Vsini Ref. (K) (CGS) (kms ) (kms ) 08 π Cet B IV E B8V C Sirius AV+DA Cnc AV C ι Her B IV SPB Vega A0Va C Peg B9.V C References: = Fossati et al. (2009); 2 = Bailey & Landstreet (20); = Hill & Landstreet (99); = Fossati et al. (200); = Nieva & Przybilla (202); = Przybilla et al. (2000).

7 The interpretation of C i emission lines K 000 K 00 K K 0000 K 020 K 0.9 log g = log g = K 900 K 90 K 0000 K log g = K 000 K 000 K log g = K 000 K 000 K K 000 K 000 K log g = log g = Figure. Change of the line profile with a variation in T eff and log g for Ci 90Å, 98Å, and 9088Å. Everywhere, [C/Fe] = 0, Vsini = 0 kms, and ξ t = kms. The theoretical spectra are convolved with an instrumental profile of R = 000. (HD 2) were determined by Przybilla et al. (2000) from the Balmer line wings and Mgi and Mgii lines. Vega is a rapidly rotating star seen pole-on. Rapid rotation causes a change of the surface shape from spherical to ellipsoidal one, therefore, the temperature stratification and surface gravity vary from the pole to equator (Hill et al. 200). We ignore the non-spherical effects in our study and analyse Vega s flux spectrum using the average temperature and gravity. Vega was also classified as a mild λ Bootis-type star (Venn & Lambert 990). ι Her (HD 02) is a single bright star, with Vsini kms. Its atmospheric parameters were determined spectroscopically from all available H and He lines and multiple ionisation equilibria, and they were confirmed via the spectral energy distribution and the Hipparcos distance (Nieva & Przybilla 202). HD is ablue Straggler and amember of the Praesepe cluster. It is previously considered as an Ap (Si) star, but appears to have the abundances of a normal A-type star (Fossati et al. 200). This star is a primary component of SB, as is the case for many other Blue Stragglers (Leonard 99). The flux coming from the secondary star is negligible, as previously checked by Burkhart & Coupry (998), so the star was analysed ignoring the presence of the secondary. The physical parameters were taken from Fossati et al. (200). Both Fe i excitation (for the effective temperature) and Fe i/fe ii ionisation (for log g) equilibria were used in the parameter determination.

8 8 S. A. Alexeeva, T. A. Ryabchikova and L. I. Mashonkina Energy (cm - ) E8 -E9 -E8-8E 2E0 -E8 2p -E8-2E9 E -2E8 -E8 -E9 2E0 8 -E9 E0 29 2p E E0 -E9 E -2E9 2p E9-9E9 9 -E8 2E0 -E9 2p 2 e -E9 -E8 -E8 -E8 2p 2p np S ns P o np P nd P o np D nd D o nd F o np S ns P o np P nd P o np D nd D o nd F o S o g o f h Figure. The NET-diagram at logτ 000 =. in the model 000/./0. Wavelengths (Å ) and NET values (s cm ) are quoted for the selected transitions. NETs of negative sign are marked in blue, and positive marked in red (only on-line available). Table. Characteristics of observed spectra. HD V Telescope/ Spectral range t exp Observing run R S/N (mag) spectrograph (Å) (s) year/month / / / / / / / / Notes. V is a visual magnitude from the SIMBAD database. Telescope / spectrograph: = CFHT/ESPaDOnS; 2 = HST/GHRS. The Sirius binary system (HD 89 = HIP 29) is composed of a main sequence AV star and a hot DA white dwarf. It is an astrometric visual binary system at a distance of only 2. pc. Their physical parameters are available with the highest precision and were taken from Hill & Landstreet (99). Sirius A is classified as a hot metallic-line (Am) star..2 Analysis of C i emission lines in B-type stars Emission lines of Ci were detected in the observed spectra of the four stars: 2 Peg, HD 22, π Cet, and ι Her. In 2 Peg, HD 22, and π Cet only Ci 8 and 90Å are in emission, while the hottest star of our sample, ι Her, also has the emission lines at 90 9Å and 90 98Å. Some examples are shown in Fig..

9 The interpretation of C i emission lines Peg 2Peg 2Peg Peg Peg Peg i Her Fe II i Her I i Her i Her.0.02 i Her I i Her Figure. Best fits (solid curves) of the selected Ci and Cii lines in 2 Peg (two top rows) and ι Her (two bottom rows). The observed spectra are shown by bold dots. For each star, our NLTE calculations reproduce well the observed Ci emission lines, using the treated model atom of Ci Cii and adopted atmospheric parameters. The best NLTE fits of the selected emission lines in 2 Peg and ι Her are presented in Fig.. Analysis of stellar emission lines provides an opportunity to test collisional data for Ci. As described in Sect.2., electron-impact excitation data for Ci were computed by Wang et al. (20) and Reid (99). We checked the influence of varying collisional rates on the appearance of the Ci emission lines in the model 000/./0 representing the atmosphere of 2 Peg. Calculations were performed with the three different sets of data, namely (i) 0 transitions from Wang et al. (20) plus 2 transitions from Reid (99) and plus formula of van Regemorter (92) for the remaining transitions; this is our standard model as described in Sect. 2., (ii) transitions from Reid (99) and the formula of van Regemorter (92) and Ω = for the remaining allowed and forbidden transitions, (iii) only the formula of van Regemorter (92) and Ω = were employed. Fig.8 shows that the effect on Ci 90Å is large. With the theoretical approximations, we computed the absorption line profile, and this is not supported by observations of 2 Peg, where Ci 90Å is in emission. Only applying accurate collisional rates allows us to obtain the emission. However, the emission is stronger in calculations with Reid (99) than Wang et al. (20) data. The choice was decided by comparison with Ci 90Å observed in 2 Peg. With the ele-

10 0 S. A. Alexeeva, T. A. Ryabchikova and L. I. Mashonkina WZB+R9+vR2 R9+vR2 vr Figure 8. NLTE line profiles of Ci 90Å in the model 000/./0 from calculations with collisional data of Wang et al. (20) where available, Reid (99), and van Regemorter (92) (solid curve, WZB+R9+vR2), collisional data from Reid (99) for transitions and formula of van Regemorter (92) for the remaining ones (dashed curve, R9+vR2), and using formula of van Regemorter (92) only (dotted curve, vr2). Everywhere, logǫ C = 8.. The theoretical spectra are convolved with an instrumental profile of R = 000. ment abundance determined from the C i absorption lines, a strength of the Ci 90Å emission was reproduced, when using data of Wang et al. (20). They were employed everywhere in further stellar abundance analysis. Table. Carbon NLTE abundances of Sirius λ, Å EW, må NLTE LTE NLTE Ci } 29.0 bl } bl bl bl bl bl bl Mean C i.. σ Cii.0 bl bl Mean C ii.9.9 σ Mean C i+c ii.. σ Determination of stellar carbon abundances In this section, we derive the carbon abundances of the selected stars from Ci and Cii lines using the atomic data from Table. The results for individual stars are as follows. Vega: no emission line of Ci is observed in the Vega s spectrum. For the lines longwards 080 Å we adopted equivalent widths (EW) from Lambert et al. (982). All the Ci lines can be divided into two groups. Lines in the visible spectral range (λ.2å, 0 lines) are weak, with EW < 0 må, while lines in the near-ir spectral range (02.9Å λ 0Å, 0 lines) are strong, with EW > 00 må. NLTE and LTE abundances from individual lines are presented in Table. Everywhere, NLTE leads to strengthened lines, however, the effect is small, with NLTE 0.0 dex by absolute value, for weak lines because they form in deep atmospheric layers. We note a large abundance discrepancy of 0.0 dex between the two weak lines, and Å, that cannot be removed by NLTE. A similar difference was also obtained by Przybilla et al. (200) and Stuerenburg & Holweger (990). For strong lines, NLTE varies mainly between 0.28 and 0. dex. Applying the treated model atom, we can reconcilethecarbonabundancefromthethreestrong(logǫ C = 8.±0.0) and 0 weak (log ǫ C = 8.2±0.) lines, while, in LTE, the IR lines give systematically higher abundances compared with the visible lines. Using EW from Lambert et al. (982) leads to a large line-to-line scatter of determined abundances, however, the mean NLTE abundance is consistent with that from visible and near-ir lines, Ci 908, 9088, 9Å obtained in this study from analysis of observed spectrum of Vega. In Table we show also the NLTE abundances derived by Przybilla et al. (200) and Stuerenburg & Holweger (990). It is worth noting that their abundances were reducedtotheoscillator strengths oftable. OurmeanNLTE abundance derived from all Ci lines, log ǫ C = 8.±0., is consistent with the error bars found in literature. We note, in particular, the Ci 908, 9088, 9Å lines, with large departures from LTE, for which the agreement is perfect. We inspected an impact of varying collisional data on the derived carbon abundances of Vega. In addition to our basic collisional recipe, WZB+R9+vR2, Table presents the NLTE abundances and abundance corrections for R9+vR2 and vr2, as described above. The abundance difference between WZB+R9+vR2 and R9+vR2 is small and nowhere exceeds 0.0 dex. However, applying rough theoretical approximations, i.e. the recipe vr2, leads to underestimated NLTE effects, in particular, for the strong lines. For example, the abundance difference between WZB+R9+vR2 and vr2 amounts to 0.2 dex for Ci 09Å. Sirius: the element abundances were derived from the Ci and Cii lines in a wide spectral region from UV to near- IR (Table ). The NLTE effects were found to be overall small. For lines of Cii and lines of Ci in the UV and visible range, NLTE 0.0 dex. This can be easily understood. Cii is the dominant ionisation stage. The Ci UV lines arise

11 The interpretation of C i emission lines Table. NLTE abundances from C i lines in Vega and NLTE abundance corrections. WZB+R9+vR2 R9+vR2 vr2 P200 SH990 λ, Å NLTE LTE NLTE NLTE / NLTE NLTE / NLTE NLTE NLTE / / / / / / / / / / / / / / / / } / / Mean(vis) σ / / / / / / Mean(IR) σ Mean(vis+IR) σ / / / / / / / / / / / / / / Mean(all) σ Notes: P200: Przybilla et al. (200), SH990: Stuerenburg & Holweger (990). in the transitions between the lowest levels, which are closely coupled together (see the departure coefficients in Fig. ). Both C i visible lines are weak. For C i near-ir lines, the NLTE corrections are negative and do not exceed 0. dex in absolute value. The mean NLTE abundances from lines of the two ionization stages, Ci and Cii, appear to be consistent. Although, we note rather large line-to-line scatter for Ci. Small influence of NLTE on the carbon abundance determination for Sirius makes possible to compare our results with the data obtained by different studies under the LTE assumption. A complete analysis of the literature data was given by Landstreet (20) who analysed the UV spectrum of Sirius. The average abundance, log ǫ C =.9±0.8, deduced by Landstreet (20) (his Table ) agrees well with our determination. For the other five stars of our programme, the abundance results are presented in Table. In all stars but the coolest one, HD, few/all Ci lines are observed in emission. These lines are marked by e in the corresponding NLTE columns. For comparison, we adopted the solar carbon NLTE abundance, log ǫ C = 8., derived by Alexeeva & Mashonkina (20) from the solar Ci atomic lines. HD : there are plentyof Ci lines in the spectrum of this Blue Straggler, while Cii is represented by the only weak line. The meannlte abundancesfrom lines of thetwo ionization stages, Ci and Cii, were found to be consistent, and they are 0. dex higher than the solar value, in line with the overall metallicity of HD (Fossati et al. 200). 2 Peg: this is the coolest star in our sample that shows an emission in the Ci near-ir lines at 8Å and 90Å. As discussed in Sect.2., the NLTE abundance corrections are positive and large for the Ci visible lines, with NLTE up to +0. dex. This explains low carbon abundance obtained from these lines by Fossati et al. (2009) under the LTE assumption. We determined the element abundance from the Ci emission lines, too, and they appear to be consistent with that from the Ci absorption lines. The departures from LTE are minor for Cii. NLTE provides consistent abundances from lines of the two ionisation stages, C iand C ii, including the C iemission lines, in contrast to LTE, where an abundance difference between the individual lines reaches 0. dex and the emission lines cannot be reproduced. The best NLTE fits of the Ci absorption and emission lines in 2 Peg are shown in Fig.. The obtained carbonnlteabundanceof2pegis close tothesolar value. HD 22: wecouldmeasureonlythreelinesofciand four lines of Cii in this star because of heavy blending the near-ir spectrum by the telluric lines. Element abundance

12 2 S. A. Alexeeva, T. A. Ryabchikova and L. I. Mashonkina Table. NLTE abundances and abundance corrections for the program stars. ι Her π Cet HD 22 2 Peg HD λ,å NLTE LTE NLTE NLTE LTE NLTE NLTE LTE NLTE NLTE LTE NLTE NLTE LTE NLTE Ci } e 8.9 e 8. e 8.22 e e e e e 8. e 8.0 e e Mean C i σ Cii Mean C ii σ Mean σ Notes. Emission lines are marked by symbol e.

13 The interpretation of C i emission lines derived from the Ci 90Å emission line is consistent with the NLTE abundances from theciand Ciiabsorption lines. The obtained carbon NLTE abundance of HD 22 is very similar to the solar one. π Cet: this star is hotter by 00 K than HD 22 and slightly faster rotating. Two emission lines, Ci 8 and 90 Å, were detected in its spectrum. The only measurable absorption line of C i is Å. NLTE abundance corrections for the Cii lines are mainly positive and, for each line, NLTE does not exceed 0.0 dex in absolute value. Carbon abundances from the absorption and emission lines of C i agree with each other, but an average abundance from the Ci lines is higher than that derived from the Cii lines, although a difference of 0. dex is still within 2σ. ι Her: this is the hottest star of our sample. No Ci line wasdetectedinthevisiblerange. AlltheCinear-IRlinesappear in emission. NLTE abundance corrections have different sign for different lines of Cii. For most lines they do not exceed 0. dex in absolute value, but they are large and negative, with NLTE = 0.2 dex, for strong Cii 8, 82 Å lines. Our NLTE analysis provides consistent abundances from different groups of lines, i.e. the Ci emission lines and Cii absorption lines. The best NLTE fits of the Ci emission lines and Cii absorption lines are shown in Fig.. The obtained carbon abundance of ι Her, log ǫ C = 8.±0.0, is essentially solar. Similar result, log ǫ C = 8.0±0.0, was obtained by Nieva & Przybilla (202) in their NLTE analysis of the Cii lines in ι Her. Larger dispersion in our abundance determination is mainly caused by the two lines, Cii 8 Å and 82 Å, which give about 0.2 dex higher abundance compared with that from the remaining lines, while no significant abundance discrepancy between different lines was obtained by Nieva & Przybilla (202). CONCLUSIONS The motivation of this work was to solve two problems in stellar astrophysics: the search for explanation of the appearance of emission lines in Ci in the near-ir spectral region and a reliable determination of carbon abundances for ABtype stars. We constructed a comprehensive model atom for Ci Cii using the most up-to-date atomic data and evaluated the NLTE line formation for Ci and Cii in classical D models representing the atmospheres of A and late B-type stars. Our NLTE calculations predict that some lines of Ci in the near IR spectral range may appear as emission lines depending on the atmospheric parameters. The emission appears first in the Ci 8Å and 90Å singlet lines at effective temperature of 920 K to 000 K depending on the value of log g. It is strengthened toward higher T eff, reaches a maximal level at T eff = 000 K (log g = ) and almost disappears at T eff = K. The Ci triplet lines at 90 9Å and 90 98Å come into emission at T eff > 000 K (log g = ). The mechanisms driving the Ci emission can be understood as follows. A prerequisite of the emission phenomenon is the overionization-recombination mechanism resulting in a depopulation of the lower levels of Ci to a greater extent than the upper levels. Extra depopulation of s P and s P, which are the lower levels of the transitions corresponding to the listed near-ir lines, can be caused by photon loss in the UV lines Ci 29, 90, and Å. In the models with T eff of about 0000 K only, Ci 29Å plays a role draining population of the s P singlet level effectively in the layers, where the Ci 8 and 90Å lines form, while these layers are optically thick for radiation at 90 and Å. With increasing T eff, formation depths of all the Ci UV lines shift inwards, resulting in a depopulation of not only s P, but also s P and the lower levels of the emission triplet lines at 90 9 Å and Å. Our theoretical results were confirmed with observations of the reference stars. The stellar sample consists of seven bright and apparently slow-rotating A- and late B- type stars in the solar neighbourhood. Our analysis is based on high S/N and high-resolution spectra with a broad wavelength coverage. The Ci emission lines were measured in the four hottest stars, with T eff 000 K, and theywere well reproduced in our NLTE calculations. For each star, the mean NLTE abundances from lines of the two ionisation stages, Ci and Cii, including the C i emission lines, were found to be consistent. Thus, we settled the dispute on whether the Ci emission in the late-b stars is produced by the circumstellar disc, or vertical stratification / horizontal inhomogeneity of the atmosphere, or the NLTE effects in the atmosphere. The six of our stars reveal highly uniform and closeto-solar carbon abundance. We confirm a significant underabundance of carbon in Sirius, with [C/H] = 0.2. We show an importance of applying accurate atomic data to the statistical equilibrium calculations. In particular, the Ci emission phenomenon turns out to be extremely sensitive to varying electron-impact excitation rates. If the latter are not accounted for properly, the stellar Ci emission lines cannot be reproduced. The results obtained in this study favour the collisional data for C i from predictions of Wang et al. (20). It is worth noting, WELs of Mgii, Siii, Pii, Caii, Crii, Feii, Niii, Cuii, and Hgii were detected in the near-uv and visible spectral regions for several stars not showing any trace of chromosphere (Sigut et al. 2000, Wahlgren & Hubrig 2000, 200, Castelli & Hubrig 200). According to Wahlgren (2008), WELs of some metals are observed in sharp-lined spectra of mid- to late-b type stars. The WELs are detected over a range of element abundance and are found among both chemically-normal and chemically-peculiar stars. The NLTE line formation calculations for these specific ions are highly desirable to understand mechanisms of the observed emission. Acknowledgements. This research is based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/IRFU, at the Canada France Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de l Univers of the Centre National de la Recherche Scientifique(CNRS) of France, and the University of Hawaii. We thank Oleg Zatsarinny for providing us the data on effective collision strengths for the C i transitions published in Wang et al. (20). This work was supported in part by the Russian Foundation for Basic Research (grants and ).

14 S. A. Alexeeva, T. A. Ryabchikova and L. I. Mashonkina REFERENCES Stuerenburg S., Holweger H., 990, A&A, 2, 2 van Regemorter H., 92, ApJ,, 90 Alexeeva S. A., Mashonkina L. I., 20, MNRAS,, 9 Venn K. A., Lambert D. L., 990, ApJ,, 2 Bailey J. D., Landstreet J. D., 20, A&A,, A0 Wahlgren G. M., 2008, Contributions of the Astronomical Burkhart C., Coupry M. F., 998, A&A, 8, 0 Observatory Skalnate Pleso, 8, 29 Butler K., Giddings J., 98, Newsletter on the analysis of Wahlgren G. M., Hubrig S., 2000, A&A, 2, L astronomical spectra, No. 9, University of London Wahlgren G. M., Hubrig S., 200, A&A, 8, 0 Carlsson M., Rutten R. J., Shchukina N. G., 992, A&A, Wang Y., Zatsarinny O., Bartschat K., 20, Phys. Rev. A, 2, 8, 020 Castelli F., Hubrig S., 200, A&A,, 0 Wilson N. J., Bell K. L., Hudson C. E., 200, A&A, 2, Cunto W., Mendoza C., Ochsenbein F., Zeippen C. J., 99, Bulletin d Information du Centre de Donnees Stellaires, 2, 9 Fossati L., Bagnulo S., Monier R., Khan S. A., Kochukhov O., Landstreet J., Wade G., Weiss W., 200, A&A,, 9 Fossati L., Ryabchikova T., Bagnulo S., Alecian E., Grunhut J., Kochukhov O., Wade G., 2009, A&A, 0, 9 Hibbert A., Biemont E., Godefroid M., Vaeck N., 99, A&AS, 99, 9 Hill G., Gulliver A. F., Adelman S. J., 200, ApJ, 2, 20 Hill G. M., Landstreet J. D., 99, A&A, 2, 2 Kochukhov O.,, 200, oleg/binmag.html Kupka F., Piskunov N., Ryabchikova T. A., Stempels H. C., Weiss W. W., 999, A&AS, 8, 9 Lambert D. L., Roby S. W., Bell R. A., 982, ApJ, 2, Landstreet J. D., 20, A&A, 28, A2 Leonard P. J. T., 99, in Milone E. F., Mermilliod J.-C., eds, The Origins, Evolution, and Destinies of Binary Stars in Clusters Vol. 90 of Astronomical Society of the Pacific Conference Series, Blue Stragglers in Star Clusters. p. Luo D., Pradhan A. K., Saraph H. E., Storey P. J., Yu Y., 989, Journal of Physics B Atomic Molecular Physics, 22, 89 Malfait K., Bogaert E., Waelkens C., 998, A&A,, 2 Nieva M. F., Przybilla N., 2008, A&A, 8, 99 Nieva M.-F., Przybilla N., 202, A&A, 9, A Przybilla N., Butler K., Becker S. R., KudritzkiR. P., Venn K. A., 2000, A&A, 9, 08 Przybilla N., Butler K., Kudritzki R. P., 200, A&A, 9, 9 Przybilla N., Nieva M.-F., Butler K., 20, Journal of Physics Conference Series, 28, 020 RalchenkoY.A., KramidaE., ReaderJ., Team N.A., 2008, NIST Atomic Spectra Database (version..). USA Reid R. H. G.,, 99, Electron impact exitation of Rentzsch-Holm I., 99, A&A, 2, 9 Roby S. W., Lambert D. L., 990, ApJS,, Ryabchikova T., Piskunov N., Pakhomov Y., Tsymbal V., Titarenko A., Sitnova T., Alexeeva S., Fossati L., Mashonkina L., 20, MNRAS,, 22 Rybicki G. B., Hummer D. G., 99, A&A, 2, Sadakane K., 98, PASP, 9, 8 Seaton M. J., 92, Proceedings of the Physical Society, 9, 0 Shulyak D., Tsymbal V., Ryabchikova T., Stütz C., Weiss W. W., 200, A&A, 28, 99 Sigut T. A. A., 200, A&A,, L2 Sigut T. A. A., Landstreet J. D., Shorlin S. L. S., 2000, ApJ, 0, L89 Smith K. C., Dworetsky M. M., 99, A&A, 2,